How Efficient Is Replica Exchange Molecular Dynamics? An Analytic Approach.

Replica exchange molecular dynamics (REMD) has become a standard technique for accelerating relaxation in biosimulations. Despite its widespread use, questions remain about its efficiency compared with conventional, constant temperature molecular dynamics (MD). An analytic approach is taken to describe the relative efficiency of REMD with respect to MD. This is applied to several simple two-state models and to several real proteins-protein L and the B domain of protein A-to predict the relative efficiency of REMD with respect to MD in actual applications. In agreement with others, we find the following: as long as there is a positive activation energy for folding, REMD is more efficient than MD; the effectiveness of REMD is strongly dependent on the activation enthalpy; and the efficiency of REMD for actual proteins is a strong function of the maximum temperature. Choosing the maximum temperature too high can result in REMD becoming significantly less efficient than conventional MD. A good rule of thumb appears to be to choose the maximum temperature of the REMD simulation slightly above the temperature at which the enthalpy for folding vanishes. Additionally, we find that the number of replicas in REMD, while important for simulations shorter than one or two relaxation times, has a minimal effect on the asymptotic efficiency of the method.

[1]  V. Hornak,et al.  Folding Cooperativity in a Three-stranded β-Sheet Model , 2005 .

[2]  Jane Clarke,et al.  What contributions to protein side-chain dynamics are probed by NMR experiments? A molecular dynamics simulation analysis. , 2005, Journal of molecular biology.

[3]  C. Brooks,et al.  Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. , 2006, Journal of the American Chemical Society.

[4]  Yuguang Mu,et al.  Folding, misfolding, and amyloid protofibril formation of WW domain FBP28. , 2006, Biophysical journal.

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

[6]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[7]  R. Levy,et al.  Protein folding pathways from replica exchange simulations and a kinetic network model. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

[9]  A. Mark,et al.  Convergence and sampling efficiency in replica exchange simulations of peptide folding in explicit solvent. , 2007, The Journal of chemical physics.

[10]  R. Levy,et al.  Simulating replica exchange simulations of protein folding with a kinetic network model , 2007, Proceedings of the National Academy of Sciences.

[11]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[12]  Daniel J. Sindhikara,et al.  Exchange frequency in replica exchange molecular dynamics. , 2008, The Journal of chemical physics.

[13]  A. Voter.  Parallel replica method for dynamics of infrequent events , 1998 .

[14]  Yuko Okamoto,et al.  Structures of a peptide fragment of β2‐microglobulin studied by replica‐exchange molecular dynamics simulations – Towards the understanding of the mechanism of amyloid formation , 2005, FEBS letters.

[15]  OU Zhaoyang,et al.  半屈曲性高分子電解質のLangevin動力学 ロッド-トロイド-グロビュール-コイル構造及び対イオン分布 , 2005 .

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

[17]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[18]  C. Brooks,et al.  Generation of native-like protein structures from limited NMR data, modern force fields and advanced conformational sampling , 2005, Journal of biomolecular NMR.

[19]  K. Sanbonmatsu,et al.  Exploring the energy landscape of a beta hairpin in explicit solvent. , 2001, Proteins.

[20]  M. Troyer,et al.  Optimized parallel tempering simulations of proteins. , 2006, The Journal of chemical physics.

[21]  Joan-Emma Shea,et al.  Effects of frustration, confinement, and surface interactions on the dimerization of an off-lattice beta-barrel protein. , 2005, The Journal of chemical physics.

[22]  H. Nymeyer,et al.  Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water. , 2007, Journal of structural biology.

[23]  K. Hellingwerf,et al.  Predicting the signaling state of photoactive yellow protein. , 2005, Biophysical journal.

[24]  J. Onuchic,et al.  Folding a protein in a computer: An atomic description of the folding/unfolding of protein A , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[25]  Bernd A. Berg.  Error Analysis for Markov Chain Data , 2004 .

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

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

[28]  S. Jang,et al.  Free energy surfaces of miniproteins with a ββα motif: Replica exchange molecular dynamics simulation with an implicit solvation model , 2005 .

[29]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Seokmin Shin,et al.  Amyloid beta-peptide oligomerization in silico: dimer and trimer. , 2006, The journal of physical chemistry. B.

[31]  Andrzej Kolinski,et al.  A minimal proteinlike lattice model: an alpha-helix motif. , 2005, The Journal of chemical physics.

[32]  R. Zhou.  Exploring the protein folding free energy landscape: coupling replica exchange method with P3ME/RESPA algorithm. , 2004, Journal of molecular graphics & modelling.

[33]  C. Geyer,et al.  Annealing Markov chain Monte Carlo with applications to ancestral inference , 1995 .

[34]  Cezary Czaplewski,et al.  Protein structure prediction with the UNRES force‐field using Replica‐Exchange Monte Carlo‐with‐Minimization; Comparison with MCM, CSA, and CFMC , 2005, J. Comput. Chem..

[35]  Chun Wu,et al.  Convergence of replica exchange molecular dynamics. , 2005, The Journal of chemical physics.

[36]  K. Kuczera,et al.  Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations. , 2004, Biophysical journal.

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

[38]  V. Hornak,et al.  Folding cooperativity in a three-stranded beta-sheet model. , 2005, Journal of molecular biology.

[39]  Edward Lyman,et al.  A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation. , 2006, Journal of chemical theory and computation.

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

[41]  A. Baumketner,et al.  Aggregation of polyalanine in a hydrophobic environment. , 2006, The Journal of chemical physics.

[42]  S Gnanakaran,et al.  Simulations of the pressure and temperature unfolding of an alpha-helical peptide. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[43]  D Baker,et al.  Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[44]  P. Nguyen,et al.  Structure and energy landscape of a photoswitchable peptide: A replica exchange molecular dynamics study , 2005, Proteins.

[45]  A. Garcia,et al.  Simulation of the pressure and temperature folding/unfolding equilibrium of a small RNA hairpin. , 2008, Journal of the American Chemical Society.

[46]  V. Pande,et al.  Multiplexed-replica exchange molecular dynamics method for protein folding simulation. , 2003, Biophysical journal.

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

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

[49]  I. Thorpe,et al.  Evaluating the conformational entropy of macromolecules using an energy decomposition approach. , 2006, The Journal of chemical physics.

[50]  Yuko Okamoto,et al.  Molecular mechanism for stabilizing a short helical peptide studied by generalized-ensemble simulations with explicit solvent. , 2005, Biophysical journal.

[51]  Lucjan Piela,et al.  Theoretical model of prion propagation: A misfolded protein induces misfolding , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[52]  Jian Zhang,et al.  Folding mechanism of β‐hairpins studied by replica exchange molecular simulations , 2005, Proteins.

[53]  H. Godwin,et al.  High metal concentrations are required for self-association of synaptotagmin II. , 2004, Biophysical journal.

[54]  Michael R. Shirts,et al.  Mathematical analysis of coupled parallel simulations. , 2001, Physical review letters.

[55]  C. Brooks,et al.  Interfacial folding and membrane insertion of designed peptides studied by molecular dynamics simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Bernd A. Berg,et al.  Markov chain monte carlo simulations and their statistical analysis: with web-based fortran code , 2004 .

[57]  Sheena E Radford,et al.  Microsecond folding dynamics of the F13W G29A mutant of the B domain of staphylococcal protein A by laser-induced temperature jump. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[58]  H. Nymeyer,et al.  Folding is not required for bilayer insertion: Replica exchange simulations of an α‐helical peptide with an explicit lipid bilayer , 2004, Proteins.

[59]  Emilio Gallicchio,et al.  Conformational equilibrium of cytochrome P450 BM-3 complexed with N-palmitoylglycine: a replica exchange molecular dynamics study. , 2006, Journal of the American Chemical Society.

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

[61]  Soonmin Jang,et al.  Amyloid β-Peptide Oligomerization in Silico: Dimer and Trimer , 2006 .

[62]  Ken A. Dill,et al.  Folding Very Short Peptides Using Molecular Dynamics , 2006, PLoS Comput. Biol..