Combining Elastic Network Analysis and Molecular Dynamics Simulations by Hamiltonian Replica Exchange.

Coarse-grained elastic network models (ENM) of proteins can be used efficiently to explore the global mobility of a protein around a reference structure. A new Hamiltonian-replica exchange molecular dynamics (H-RexMD) method has been designed that effectively combines information extracted from an ENM analysis with atomic-resolution MD simulations. The ENM analysis is used to construct a distance-dependent penalty (flooding or biasing) potential that can drive the structure away from its current conformation in directions compatible with the ENM model. Various levels of the penalty or biasing potential are added to the force field description of the MD simulation along the replica coordinate. One replica runs at the original force field. By focusing the penalty potential on the relevant soft degrees of freedom the method avoids the rapid increase of the replica number with increasing system size to cover a desired temperature range in conventional (temperature) RexMD simulations. The application to domain motions in lysozyme of bacteriophage T4 and to peptide folding indicates significantly improved conformational sampling compared to conventional MD simulations.

[1]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[2]  B. Berne,et al.  The free energy landscape for beta hairpin folding in explicit water. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Bert L de Groot,et al.  Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange. , 2007, Biophysical journal.

[4]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

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

[6]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[7]  R. Zhou.  Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

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

[9]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[10]  B. Matthews,et al.  Protein flexibility and adaptability seen in 25 crystal forms of T4 lysozyme. , 1995, Journal of molecular biology.

[11]  Yuko Okamoto,et al.  Generalized-ensemble algorithms: enhanced sampling techniques for Monte Carlo and molecular dynamics simulations. , 2003, Journal of molecular graphics & modelling.

[12]  Chin-Kun Hu,et al.  Free energy landscape and folding mechanism of a β‐hairpin in explicit water: A replica exchange molecular dynamics study , 2005, Proteins.

[13]  Asim Okur,et al.  Improved Efficiency of Replica Exchange Simulations through Use of a Hybrid Explicit/Implicit Solvation Model. , 2006, Journal of chemical theory and computation.

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

[15]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[16]  Martin Zacharias,et al.  Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing‐potential , 2006, Proteins.

[17]  B. Matthews,et al.  Structure of a thermostable disulfide-bridge mutant of phage T4 lysozyme shows that an engineered cross-link in a flexible region does not increase the rigidity of the folded protein. , 1990, Biochemistry.

[18]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[19]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[20]  L. Kay,et al.  What is the average conformation of bacteriophage T4 lysozyme in solution? A domain orientation study using dipolar couplings measured by solution NMR. , 2001, Journal of molecular biology.

[21]  H. Nymeyer,et al.  Simulation of the folding equilibrium of α-helical peptides: A comparison of the generalized Born approximation with explicit solvent , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[22]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[23]  F. Rao,et al.  Replica exchange molecular dynamics simulations of reversible folding , 2003 .

[24]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[25]  M. Seibert,et al.  Reproducible polypeptide folding and structure prediction using molecular dynamics simulations. , 2005, Journal of molecular biology.

[26]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  B. Matthews,et al.  Structure of bacteriophage T4 lysozyme refined at 1.7 A resolution. , 1987, Journal of molecular biology.

[28]  Ruhong Zhou,et al.  Replica exchange with solute tempering: efficiency in large scale systems. , 2007, The journal of physical chemistry. B.

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

[30]  Zhongwei Zhu,et al.  Using novel variable transformations to enhance conformational sampling in molecular dynamics. , 2002, Physical review letters.

[31]  Andrea Amadei,et al.  A molecular dynamics study of the 41‐56 β‐hairpin from B1 domain of protein G , 1999, Protein science : a publication of the Protein Society.

[32]  M. Zacharias,et al.  Accounting for global protein deformability during protein-protein and protein-ligand docking. , 2005, Biochimica et biophysica acta.

[33]  Replica-exchange molecular dynamics simulation of small peptide in water and in ethanol , 2005 .

[34]  Gregory D. Hawkins,et al.  Pairwise solute descreening of solute charges from a dielectric medium , 1995 .

[35]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[36]  Martin Zacharias,et al.  Energy minimization in low‐frequency normal modes to efficiently allow for global flexibility during systematic protein–protein docking , 2008, Proteins.

[37]  C. Predescu,et al.  On the efficiency of exchange in parallel tempering monte carlo simulations. , 2004, The journal of physical chemistry. B.

[38]  M. Zacharias,et al.  Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations. , 2007, Biophysical journal.

[39]  Structural and Dynamic Properties of the CAGQW Peptide in Water: A Molecular Dynamics Simulation Study Using Different Force Fields , 2004 .

[40]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[41]  X. Daura,et al.  Reversible peptide folding in solution by molecular dynamics simulation. , 1998, Journal of molecular biology.

[42]  S Gnanakaran,et al.  Peptide folding simulations. , 2003, Current opinion in structural biology.

[43]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  I. Tavernelli,et al.  A Novel Hamiltonian Replica Exchange MD Protocol to Enhance Protein Conformational Space Sampling. , 2006, Journal of chemical theory and computation.

[45]  D. Bruce,et al.  Replica Exchange Molecular Dynamics Modeling of Foldamers , 2004 .

[46]  I. Bahar,et al.  Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[47]  K. Sanbonmatsu,et al.  Exploring the energy landscape of a β hairpin in explicit solvent , 2001 .

[48]  H. Berendsen,et al.  Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data , 1998, Proteins.

[49]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[50]  A. Roitberg,et al.  All-atom structure prediction and folding simulations of a stable protein. , 2002, Journal of the American Chemical Society.

[51]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[52]  Seokmin Shin,et al.  Replica-exchange method using the generalized effective potential. , 2003, Physical review letters.

[53]  Julia M. Goodfellow,et al.  Molecular dynamics study , 1997 .

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

[55]  V S Pande,et al.  Molecular dynamics simulations of unfolding and refolding of a beta-hairpin fragment of protein G. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[56]  V. Hornak,et al.  Modified replica exchange simulation methods for local structure refinement. , 2005, The journal of physical chemistry. B.

[57]  Yuko Okamoto,et al.  Replica-exchange Monte Carlo method for the isobaric isothermal ensemble , 2001 .

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