Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes.

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.

[1]  G. Schulz,et al.  Adenylate kinase motions during catalysis: an energetic counterweight balancing substrate binding. , 1996, Structure.

[2]  G. Schulz,et al.  Structure of a mutant adenylate kinase ligated with an ATP-analogue showing domain closure over ATP. , 1996, Journal of molecular biology.

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

[4]  G. Schulz,et al.  The structure of bovine mitochondrial adenylate kinase: Comparison with isoenzymes in other compartments , 1996, Protein science : a publication of the Protein Society.

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

[6]  G. Schulz,et al.  Structure of the complex between adenylate kinase from Escherichia coli and the inhibitor Ap5A refined at 1.9 A resolution. A model for a catalytic transition state. , 1992, Journal of molecular biology.

[7]  Herbert S. Rosenkranz,et al.  Studies on deoxyribonucleic acid after exposure to tritium gas , 1959 .

[8]  Bernard R Brooks,et al.  Normal-modes-based prediction of protein conformational changes guided by distance constraints. , 2005, Biophysical journal.

[9]  Haiyan Liu,et al.  Efficiently explore the energy landscape of proteins in molecular dynamics simulations by amplifying collective motions , 2003 .

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

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

[12]  V. Pande,et al.  Can conformational change be described by only a few normal modes? , 2006, Biophysical journal.

[13]  A. Atilgan,et al.  Vibrational Dynamics of Folded Proteins: Significance of Slow and Fast Motions in Relation to Function and Stability , 1998 .

[14]  Martin Karplus,et al.  Simulation of conformational transitions by the restricted perturbation-targeted molecular dynamics method. , 2005, The Journal of chemical physics.

[15]  Dror Tobi,et al.  Allosteric changes in protein structure computed by a simple mechanical model: hemoglobin T<-->R2 transition. , 2003, Journal of molecular biology.

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

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

[18]  T. Haliloglu,et al.  Conformational dynamics of chymotrypsin inhibitor 2 by coarse‐grained simulations , 1999, Proteins.

[19]  D Perahia,et al.  Motions in hemoglobin studied by normal mode analysis and energy minimization: evidence for the existence of tertiary T-like, quaternary R-like intermediate structures. , 1996, Journal of molecular biology.

[20]  B. Vallone,et al.  The structures of deoxy human haemoglobin and the mutant Hb Tyrα42His at 120 K , 2000 .

[21]  I. Bahar,et al.  Coarse-grained normal mode analysis in structural biology. , 2005, Current opinion in structural biology.

[22]  R L Jernigan,et al.  Vibrational dynamics of transfer RNAs: comparison of the free and synthetase-bound forms. , 1998, Journal of molecular biology.

[23]  P. Rogers,et al.  A third quaternary structure of human hemoglobin A at 1.7-A resolution. , 1992, The Journal of biological chemistry.

[24]  T. Haliloglu.  Coarse-grained simulations of the conformational dynamics of proteins , 1999 .

[25]  Osamu Miyashita,et al.  Conformational transitions of adenylate kinase: switching by cracking. , 2007, Journal of molecular biology.

[26]  Ivet Bahar,et al.  DYNAMICS OF PROTEINS AND BIOMOLECULAR COMPLEXES: INFERRING FUNCTIONAL MOTIONS FROM STRUCTURE , 1999 .

[27]  C. Brooks,et al.  Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism , 2007, Proceedings of the National Academy of Sciences.

[28]  Mark A. Wilson,et al.  Intrinsic motions along an enzymatic reaction trajectory , 2007, Nature.

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

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

[31]  Robert L Jernigan,et al.  Rigid-cluster models of conformational transitions in macromolecular machines and assemblies. , 2005, Biophysical journal.

[32]  Jianpeng Ma,et al.  Folding of small helical proteins assisted by small-angle X-ray scattering profiles. , 2005, Structure.

[33]  M. Karplus,et al.  The allosteric mechanism of yeast chorismate mutase: a dynamic analysis. , 2006, Journal of molecular biology.

[34]  I. Bahar,et al.  Coarse‐grained simulations of conformational dynamics of proteins: Application to apomyoglobin , 1998, Proteins.

[35]  G. Phillips,et al.  Crystal structure of ADP/AMP complex of Escherichia coli adenylate kinase , 2005, Proteins.

[36]  Robert I Cukier,et al.  Molecular dynamics of apo-adenylate kinase: a principal component analysis. , 2006, The journal of physical chemistry. B.

[37]  J. G. Lewis,et al.  A Shifted Block Lanczos Algorithm for Solving Sparse Symmetric Generalized Eigenproblems , 1994, SIAM J. Matrix Anal. Appl..

[38]  M. Gerstein,et al.  The morph server: a standardized system for analyzing and visualizing macromolecular motions in a database framework. , 2000, Nucleic acids research.

[39]  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.

[40]  I. Bahar,et al.  Gaussian Dynamics of Folded Proteins , 1997 .

[41]  Bernard R Brooks,et al.  Modeling protein conformational changes by iterative fitting of distance constraints using reoriented normal modes. , 2006, Biophysical journal.

[42]  L. Mouawad,et al.  New insights into the allosteric mechanism of human hemoglobin from molecular dynamics simulations. , 2002, Biophysical journal.

[43]  R L Jernigan,et al.  Short‐range conformational energies, secondary structure propensities, and recognition of correct sequence‐structure matches , 1997, Proteins.

[44]  G. Rose,et al.  The T-to-R transformation in hemoglobin: a reevaluation. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[46]  Ruth Nussinov,et al.  HingeProt: Automated prediction of hinges in protein structures , 2008, Proteins.

[47]  Martin Karplus,et al.  Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase. , 2005, Journal of molecular biology.

[48]  Martin Karplus,et al.  Minimum free energy pathways and free energy profiles for conformational transitions based on atomistic molecular dynamics simulations. , 2007, The Journal of chemical physics.

[49]  O. Marques.  BLZPACK: Description and User's Guide , 1995 .

[50]  G. Phillips,et al.  The closed conformation of a highly flexible protein: The structure of E. coli adenylate kinase with bound AMP and AMPPNP , 1994, Proteins.

[51]  B. Shaanan,et al.  Structure of human oxyhaemoglobin at 2.1 A resolution. , 1983, Journal of molecular biology.

[52]  G. Schulz,et al.  Crystal structures of two mutants of adenylate kinase from Escherichia coli that modify the Gly‐loop , 1993, Proteins.

[53]  Guang Song,et al.  How well can we understand large-scale protein motions using normal modes of elastic network models? , 2007, Biophysical journal.

[54]  B. Vallone,et al.  The structures of deoxy human haemoglobin and the mutant Hb Tyralpha42His at 120 K. , 2000, Acta crystallographica. Section D, Biological crystallography.

[55]  R. Jernigan,et al.  Inter-residue potentials in globular proteins and the dominance of highly specific hydrophilic interactions at close separation. , 1997, Journal of molecular biology.

[56]  G. Chirikjian,et al.  Efficient generation of feasible pathways for protein conformational transitions. , 2002, Biophysical journal.