Energetic and Dynamic Aspects of the Affinity Maturation Process: Characterizing Improved Variants from the Bevacizumab Antibody with Molecular Simulations

Antibody affinity maturation is one of the fundamental processes of immune defense against invading pathogens. From the biological point of view, the clonal selection hypothesis represents the most accepted mechanism to explain how mutations increasing the affinity for target antigens are introduced and selected in antibody molecules. However, understanding at the molecular level how protein modifications, such as point mutation, can modify and modulate the affinity of an antibody for its antigen is still a major open issue in molecular biology. In this paper, we address various aspects of this problem by analyzing and comparing atomistic simulations of 17 variants of the bevacizumab antibody, all directed against the common target protein VEGF-A. In particular, we examine MD-based descriptors of the internal energetics and dynamics of mutated antibodies and their possible correlations with experimentally determined affinities for the antigens. Our results show that affinity improvement is correlated with a variation of the internal stabilization energy of the antibody molecule when bound to the antigen, compensated by the variation in the interaction energy between the antigen and the antibody, paralleled by an overall modulation of internal coordination within the antibody molecular structure. A possible model of the mechanism of rigidification and of the main residues involved is proposed. Overall, our results can help in understanding the molecular determinants of antigen recognition and have implications in the rational design of new antibodies with optimized affinities.

[1]  G. Colombo,et al.  Investigating dynamic and energetic determinants of protein nucleic acid recognition: analysis of the zinc finger zif268-DNA complexes , 2010, BMC Structural Biology.

[2]  G. Tiana,et al.  Similar folds with different stabilization mechanisms: the cases of prion and doppel proteins , 2006, BMC Structural Biology.

[3]  M. Oda,et al.  Protein Dynamics and the Diversity of an Antibody Response* , 2012, The Journal of Biological Chemistry.

[4]  Giorgio Colombo,et al.  Identification of domains in protein structures from the analysis of intramolecular interactions. , 2012, The journal of physical chemistry. B.

[5]  C. Milstein,et al.  Memory in the B-cell compartment: antibody affinity maturation. , 2000, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[6]  Leonard,et al.  Humanization of an anti-vascular endothelial growth factor monoclonal antibody for the therapy of solid tumors and other disorders. , 1997, Cancer research.

[7]  Matthias J. Feige,et al.  An unfolded CH1 domain controls the assembly and secretion of IgG antibodies. , 2009, Molecular cell.

[8]  Ian F. Thorpe,et al.  Antibody evolution constrains conformational heterogeneity by tailoring protein dynamics , 2006, Proceedings of the National Academy of Sciences.

[9]  Giorgio Colombo,et al.  Modeling Signal Propagation Mechanisms and Ligand-Based Conformational Dynamics of the Hsp90 Molecular Chaperone Full-Length Dimer , 2009, PLoS Comput. Biol..

[10]  I. Roterman,et al.  Analysis of correlated domain motions in IgG light chain reveals possible mechanisms of immunological signal transduction , 2005, Proteins.

[11]  H. Nakamura,et al.  Junctional amino acids determine the maturation pathway of an antibody. , 1999, Immunity.

[12]  P. Rousseeuw.  Silhouettes: a graphical aid to the interpretation and validation of cluster analysis , 1987 .

[13]  C. Robinson,et al.  Regional and segmental flexibility of antibodies in interaction with antigens of different size , 2006, The FEBS journal.

[14]  Andrew T. Fenley,et al.  Entropy–enthalpy transduction caused by conformational shifts can obscure the forces driving protein–ligand binding , 2012, Proceedings of the National Academy of Sciences.

[15]  G. Colombo,et al.  HMGB1-carbenoxolone interactions: dynamics insights from combined nuclear magnetic resonance and molecular dynamics. , 2011, Chemistry, an Asian journal.

[16]  Bruce Tidor,et al.  Computational design of antibody-affinity improvement beyond in vivo maturation , 2007, Nature Biotechnology.

[17]  Giorgio Colombo,et al.  Determinants of protein stability and folding: Comparative analysis of beta‐lactoglobulins and liver basic fatty acid binding protein , 2005, Proteins.

[18]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[19]  Ricardo A Broglia,et al.  Understanding the determinants of stability and folding of small globular proteins from their energetics , 2003, Protein science : a publication of the Protein Society.

[20]  Giorgio Colombo,et al.  Investigating allostery in molecular recognition: insights from a computational study of multiple antibody-antigen complexes. , 2013, The journal of physical chemistry. B.

[21]  A. D. de Vos,et al.  Selection and analysis of an optimized anti-VEGF antibody: crystal structure of an affinity-matured Fab in complex with antigen. , 1999, Journal of molecular biology.

[22]  Paolo Valerio,et al.  Vascular endothelial growth factor (VEGF) as a target of bevacizumab in cancer: from the biology to the clinic. , 2006, Current medicinal chemistry.

[23]  J. Schlitter.  Estimation of absolute and relative entropies of macromolecules using the covariance matrix , 1993 .

[24]  Peter J. Rousseeuw,et al.  Finding Groups in Data: An Introduction to Cluster Analysis , 1990 .

[25]  M. Buckle,et al.  Can immunoglobulin C(H)1 constant region domain modulate antigen binding affinity of antibodies? , 1996, The Journal of clinical investigation.

[26]  F. Peale,et al.  Mice expressing a humanized form of VEGF-A may provide insights into the safety and efficacy of anti-VEGF antibodies , 2007, Proceedings of the National Academy of Sciences.

[27]  Bernhardt L. Trout,et al.  Design of therapeutic proteins with enhanced stability , 2009, Proceedings of the National Academy of Sciences.

[28]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[29]  A. Sali,et al.  Modeling of loops in protein structures , 2000, Protein science : a publication of the Protein Society.

[30]  Charles Eigenbrot,et al.  Crystal Structure at 1.7 Å Resolution of VEGF in Complex with Domain 2 of the Flt-1 Receptor , 1997, Cell.

[31]  Kenneth J. Hillan,et al.  Discovery and development of bevacizumab, an anti-VEGF antibody for treating cancer , 2004, Nature Reviews Drug Discovery.

[32]  Levi C. T. Pierce,et al.  Routine Access to Millisecond Time Scale Events with Accelerated Molecular Dynamics , 2012, Journal of chemical theory and computation.

[33]  Gunnar Jeschke,et al.  DEER distance measurements on proteins. , 2012, Annual review of physical chemistry.

[34]  A. Casadevall,et al.  The Immunoglobulin Heavy Chain Constant Region Affects Kinetic and Thermodynamic Parameters of Antibody Variable Region Interactions with Antigen* , 2007, Journal of Biological Chemistry.

[35]  Giorgio Colombo,et al.  Relationship between energy distribution and fold stability: Insights from molecular dynamics simulations of native and mutant proteins , 2008, Proteins.

[36]  L. Cavacini,et al.  Structure and function of immunoglobulins. , 2010, The Journal of allergy and clinical immunology.

[37]  A. Casadevall,et al.  The immunoglobulin constant region contributes to affinity and specificity. , 2008, Trends in immunology.

[38]  Luke N Robinson,et al.  Redesign of a cross-reactive antibody to dengue virus with broad-spectrum activity and increased in vivo potency , 2013, Proceedings of the National Academy of Sciences.

[39]  Johannes Buchner,et al.  How antibodies fold. , 2010, Trends in biochemical sciences.

[40]  M. Oda,et al.  Exploring the energy landscape of antibody-antigen complexes: protein dynamics, flexibility, and molecular recognition. , 2008, Biochemistry.

[41]  R C Stevens,et al.  Structural insights into the evolution of an antibody combining site. , 1997, Science.

[42]  N. Ferrara,et al.  Pituitary follicular cells secrete a novel heparin-binding growth factor specific for vascular endothelial cells. 1989. , 2012, Biochemical and biophysical research communications.

[43]  I. Roterman,et al.  The Indirect Generation of Long‐distance Structural Changes in Antibodies upon their Binding to Antigen , 2006, Chemical biology & drug design.

[44]  L. Lopalco,et al.  Isotype modulates epitope specificity, affinity, and antiviral activities of anti–HIV-1 human broadly neutralizing 2F5 antibody , 2012, Proceedings of the National Academy of Sciences.

[45]  M. Parrinello,et al.  Funnel metadynamics as accurate binding free-energy method , 2013, Proceedings of the National Academy of Sciences.

[46]  Ian F. Thorpe,et al.  Molecular evolution of affinity and flexibility in the immune system , 2007, Proceedings of the National Academy of Sciences.

[47]  K. Rajewsky.  Clonal selection and learning in the antibody system , 1996, Nature.

[48]  Napoleone Ferrara,et al.  Vascular endothelial growth factor: basic science and clinical progress. , 2004, Endocrine reviews.

[49]  N. Shimba,et al.  Role of the domain-domain interaction in the construction of the antigen combining site. A comparative study by 1H-15N shift correlation NMR spectroscopy of the Fv and Fab fragments of anti-dansyl mouse monoclonal antibody. , 1994, Journal of molecular biology.

[50]  D. Benjamin,et al.  Reduction in the amide hydrogen exchange rates of an anti-lysozyme Fv fragment due to formation of the Fv-lysozyme complex. , 1997, Journal of molecular biology.

[51]  Andrew C. R. Martin,et al.  Analysis and improvements to Kabat and structurally correct numbering of antibody variable domains. , 2008, Molecular immunology.

[52]  A. D. de Vos,et al.  VEGF and the Fab fragment of a humanized neutralizing antibody: crystal structure of the complex at 2.4 A resolution and mutational analysis of the interface. , 1998, Structure.

[53]  Giorgio Colombo,et al.  Corresponding Functional Dynamics across the Hsp90 Chaperone Family: Insights from a Multiscale Analysis of MD Simulations , 2012, PLoS Comput. Biol..

[54]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .