Multiscale modeling of macromolecular conformational changes combining concepts from rigidity and elastic network theory

The development of a two‐step approach for multiscale modeling of macromolecular conformational changes is based on recent developments in rigidity and elastic network theory. In the first step, static properties of the macromolecule are determined by decomposing the molecule into rigid clusters by using the graph‐theoretical approach FIRST and an all‐atom representation of the protein. In this way, rigid clusters are not limited to consist of residues adjacent in sequence or secondary structure elements as in previous studies. Furthermore, flexible links between rigid clusters are identified and can be modeled as such subsequently. In the second step, dynamical properties of the molecule are revealed by the rotations‐translations of blocks approach (RTB) using an elastic network model representation of the coarse‐grained protein. In this step, only rigid body motions are allowed for rigid clusters, whereas links between them are treated as fully flexible. The approach was tested on a data set of 10 proteins that showed conformational changes on ligand binding. For efficiency, coarse‐graining the protein results in a remarkable reduction of memory requirements and computational times by factors of 9 and 27 on average and up to 25 and 125, respectively. For accuracy, directions and magnitudes of motions predicted by our approach agree well with experimentally determined ones, despite embracing in extreme cases >50% of the protein into one rigid cluster. In fact, the results of our method are in general comparable with when no or a uniform coarse‐graining is applied; and the results are superior if the movement is dominated by loop or fragment motions. This finding indicates that explicitly distinguishing between flexible and rigid regions is advantageous when using a simplified protein representation in the second step. Finally, motions of atoms in rigid clusters are also well predicted by our approach, which points to the need to consider mobile protein regions in addition to flexible ones when modeling correlated motions. Proteins 2006. © 2006 Wiley‐Liss, Inc.

[1]  J. Maxwell,et al.  The Scientific Papers of James Clerk Maxwell: On the Calculation of the Equilibrium and Stiffness of Frames , 1864 .

[2]  WM. CHAPPELL,et al.  Molecular Vibrations , 1879, Nature.

[3]  G. Laman On graphs and rigidity of plane skeletal structures , 1970 .

[4]  M. Karplus,et al.  Harmonic dynamics of proteins: normal modes and fluctuations in bovine pancreatic trypsin inhibitor. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[5]  N. Go,et al.  Dynamics of a small globular protein in terms of low-frequency vibrational modes. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Y. Sanejouand,et al.  A new approach for determining low‐frequency normal modes in macromolecules , 1994 .

[7]  D. Case Normal mode analysis of protein dynamics , 1994 .

[8]  R. Brüschweiler Collective protein dynamics and nuclear spin relaxation , 1995 .

[9]  J. Kyte,et al.  Structure in Protein Chemistry , 1995 .

[10]  C. Mandl,et al.  Oligomeric rearrangement of tick-borne encephalitis virus envelope proteins induced by an acidic pH , 1995, Journal of virology.

[11]  Y. Sanejouand,et al.  Hinge‐bending motion in citrate synthase arising from normal mode calculations , 1995, Proteins.

[12]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[13]  Jacobs,et al.  Generic rigidity percolation: The pebble game. , 1995, Physical review letters.

[14]  L. Kay,et al.  Correlation between dynamics and high affinity binding in an SH2 domain interaction. , 1996, Biochemistry.

[15]  Tirion,et al.  Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. , 1996, Physical review letters.

[16]  W. Junge,et al.  Intersubunit rotation in active F-ATPase , 1996, Nature.

[17]  A. Mildvan,et al.  15N NMR relaxation studies of free and inhibitor-bound 4-oxalocrotonate tautomerase: backbone dynamics and entropy changes of an enzyme upon inhibitor binding. , 1996, Biochemistry.

[18]  S. Fesik,et al.  Backbone dynamics of the C-terminal domain of Escherichia coli topoisomerase I in the absence and presence of single-stranded DNA. , 1996, Biochemistry.

[19]  Jacobs,et al.  Generic rigidity percolation in two dimensions. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[20]  E. Olejniczak,et al.  Changes in the NMR-derived motional parameters of the insulin receptor substrate 1 phosphotyrosine binding domain upon binding to an interleukin 4 receptor phosphopeptide. , 1997, Biochemistry.

[21]  A. Horwich,et al.  The crystal structure of the asymmetric GroEL–GroES–(ADP)7 chaperonin complex , 1997, Nature.

[22]  M. Hodsdon,et al.  Ligand binding alters the backbone mobility of intestinal fatty acid-binding protein as monitored by 15N NMR relaxation and 1H exchange. , 1997, Biochemistry.

[23]  A. Atilgan,et al.  Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. , 1997, Folding & design.

[24]  L. Kay,et al.  Correlation between binding and dynamics at SH2 domain interfaces , 1998, Nature Structural Biology.

[25]  K. Hinsen Analysis of domain motions by approximate normal mode calculations , 1998, Proteins.

[26]  Milos V. Novotny,et al.  Increased protein backbone conformational entropy upon hydrophobic ligand binding , 1999, Nature Structural Biology.

[27]  Julie D. Forman-Kay,et al.  The 'dynamics' in the thermodynamics of binding , 1999, Nature Structural Biology.

[28]  R. Nussinov,et al.  Folding and binding cascades: shifts in energy landscapes. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[29]  R. Nussinov,et al.  Folding funnels, binding funnels, and protein function , 1999, Protein science : a publication of the Protein Society.

[30]  Joachim Frank,et al.  A ratchet-like inter-subunit reorganization of the ribosome during translocation , 2000, Nature.

[31]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[32]  Y. Sanejouand,et al.  Building‐block approach for determining low‐frequency normal modes of macromolecules , 2000, Proteins.

[33]  Andrew L. Lee,et al.  Redistribution and loss of side chain entropy upon formation of a calmodulin–peptide complex , 2000, Nature Structural Biology.

[34]  Y. Sanejouand,et al.  Conformational change of proteins arising from normal mode calculations. , 2001, Protein engineering.

[35]  D. Jacobs,et al.  Protein flexibility predictions using graph theory , 2001, Proteins.

[36]  R E Oswald,et al.  An increase in side chain entropy facilitates effector binding: NMR characterization of the side chain methyl group dynamics in Cdc42Hs. , 2001, Biochemistry.

[37]  R. Jernigan,et al.  Anisotropy of fluctuation dynamics of proteins with an elastic network model. , 2001, Biophysical journal.

[38]  D. Jacobs,et al.  Protein flexibility and dynamics using constraint theory. , 2001, Journal of molecular graphics & modelling.

[39]  M. Delarue,et al.  Simplified normal mode analysis of conformational transitions in DNA-dependent polymerases: the elastic network model. , 2002, Journal of molecular biology.

[40]  H. Carlson Protein flexibility and drug design: how to hit a moving target. , 2002, Current opinion in chemical biology.

[41]  A. Rader,et al.  Identifying protein folding cores from the evolution of flexible regions during unfolding. , 2002, Journal of molecular graphics & modelling.

[42]  Florence Tama,et al.  The mechanism and pathway of pH induced swelling in cowpea chlorotic mottle virus. , 2002, Journal of molecular biology.

[43]  Robert L. Jernigan,et al.  Dynamics of large proteins through hierarchical levels of coarse‐grained structures , 2002, J. Comput. Chem..

[44]  Leslie A Kuhn,et al.  Protein unfolding: Rigidity lost , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[45]  M. Thorpe,et al.  Rigidity theory and applications , 2002 .

[46]  Leslie A. Kuhn,et al.  Flexible and Rigid Regions in Proteins , 2002 .

[47]  J. Mccammon,et al.  Changes in flexibility upon binding: Application of the self-consistent pair contact probability method to protein-protein interactions , 2002 .

[48]  W. Wriggers,et al.  Exploring global distortions of biological macromolecules and assemblies from low-resolution structural information and elastic network theory. , 2002, Journal of molecular biology.

[49]  Guohui Li,et al.  A coarse-grained normal mode approach for macromolecules: an efficient implementation and application to Ca(2+)-ATPase. , 2002, Biophysical journal.

[50]  D. Ming,et al.  How to describe protein motion without amino acid sequence and atomic coordinates , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[51]  Donald J. Jacobs,et al.  GENERIC RIGIDITY OF NETWORK GLASSES , 2002 .

[52]  Mark Gerstein,et al.  MolMovDB: analysis and visualization of conformational change and structural flexibility , 2003, Nucleic Acids Res..

[53]  S. Teague Implications of protein flexibility for drug discovery , 2003, Nature Reviews Drug Discovery.

[54]  Haiyan Liu,et al.  Molecular dynamics simulations of peptides and proteins with amplified collective motions. , 2003, Biophysical journal.

[55]  Robert L Jernigan,et al.  Functional motions can be extracted from on‐lattice construction of protein structures , 2003, Proteins.

[56]  Keith Brew,et al.  Increased backbone mobility in beta-barrel enhances entropy gain driving binding of N-TIMP-1 to MMP-3. , 2003, Journal of molecular biology.

[57]  Sebastian Doniach,et al.  A comparative study of motor-protein motions by using a simple elastic-network model , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[58]  Leslie A Kuhn,et al.  Modeling correlated main‐chain motions in proteins for flexible molecular recognition , 2004, Proteins.

[59]  Ming Lei,et al.  Sampling protein conformations and pathways , 2004, J. Comput. Chem..

[60]  Brian K Shoichet,et al.  Testing a flexible-receptor docking algorithm in a model binding site. , 2004, Journal of molecular biology.

[61]  Holger Gohlke,et al.  Change in protein flexibility upon complex formation: Analysis of Ras‐Raf using molecular dynamics and a molecular framework approach , 2004, Proteins.

[62]  Martin Zacharias,et al.  Rapid protein–ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: Binding of FK506 to FKBP , 2004, Proteins.

[63]  Gregory S Chirikjian,et al.  Normal mode analysis of proteins: a comparison of rigid cluster modes with C(alpha) coarse graining. , 2004, Journal of molecular graphics & modelling.

[64]  Gregory S. Chirikjian,et al.  Normal mode analysis of proteins: a comparison of rigid cluster modes with Cα coarse graining , 2004 .

[65]  M. Gerstein,et al.  Conformational changes associated with protein-protein interactions. , 2004, Current opinion in structural biology.

[66]  Qiang Cui,et al.  Analysis of functional motions in Brownian molecular machines with an efficient block normal mode approach: myosin-II and Ca2+ -ATPase. , 2004, Biophysical journal.

[67]  Valentina Tozzini,et al.  Coarse-grained models for proteins. , 2005, Current opinion in structural biology.

[68]  Osamu Miyashita,et al.  Simple energy landscape model for the kinetics of functional transitions in proteins. , 2005, The journal of physical chemistry. B.

[69]  Jianpeng Ma,et al.  Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. , 2005, Structure.