A convective replica‐exchange method for sampling new energy basins

Replica‐exchange is a powerful simulation method for sampling the basins of a rugged energy landscape. The replica‐exchange method's sampling is efficient because it allows replicas to perform round trips in temperature space, thereby visiting both low and high temperatures in the same simulation. However, replicas have a diffusive walk in temperature space, and the round trip rate decreases significantly with the system size. These drawbacks make convergence of the simulation even more difficult than it already is when bigger systems are tackled. Here, we present a simple modification of the exchange method. In this method, one of the replicas steadily raises or lowers its temperature. We tested the convective replica‐exchange method on three systems of varying complexity: the alanine dipeptide in implicit solvent, the GB1 β‐hairpin in explicit solvent and the Aβ25–35 homotrimer in a coarse grained representation. For the highly frustrated Aβ25–35 homotrimer, the proposed “convective” replica‐exchange method is twice as fast as the standard method. It discovered 24 out of 27 free‐energy basins in less than 500 ns. It also prevented the formation of groups of replicas that usually form on either side of an exchange bottleneck, leading to a more efficient sampling of new energy basins than in the standard method. © 2012 Wiley Periodicals, Inc.

[1]  Ulrich H E Hansmann,et al.  Efficient sampling of protein structures by model hopping. , 2005, Physical review letters.

[2]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[3]  Ulrich H E Hansmann,et al.  Folding proteins by first-passage-times-optimized replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  David A. Kofke,et al.  ARTICLES On the acceptance probability of replica-exchange Monte Carlo trials , 2002 .

[5]  N. Mousseau,et al.  Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent. , 2009, The journal of physical chemistry. B.

[6]  M. Mattson,et al.  Aβ25–35 induces rapid lysis of red blood cells: contrast with Aβ1–42 and examination of underlying mechanisms , 1997, Brain Research.

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

[8]  Wolfhard Janke,et al.  Make Life Simple: Unleash the Full Power of the Parallel Tempering Algorithm , 2008 .

[9]  Ulrich H E Hansmann,et al.  Generalized ensemble and tempering simulations: a unified view. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  J. D. de Pablo,et al.  Optimal allocation of replicas in parallel tempering simulations. , 2005, The Journal of chemical physics.

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

[12]  Yannick G. Spill,et al.  Impact of Thermostats on Folding and Aggregation Properties of Peptides Using the Optimized Potential for Efficient Structure Prediction Coarse-Grained Model. , 2011, Journal of chemical theory and computation.

[13]  Ulrich H E Hansmann,et al.  Microcanonical replica exchange molecular dynamics simulation of proteins. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  P. Derreumaux,et al.  A coarse‐grained protein force field for folding and structure prediction , 2007, Proteins.

[15]  Ulrich H E Hansmann,et al.  Velocity scaling for optimizing replica exchange molecular dynamics. , 2011, The Journal of chemical physics.

[16]  R. Tycko,et al.  Abeta40-Lactam(D23/K28) models a conformation highly favorable for nucleation of amyloid. , 2005, Biochemistry.

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

[18]  D. Huse,et al.  Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[20]  G. Forloni,et al.  Apoptosis mediated neurotoxicity induced by chronic application of beta amyloid fragment 25-35. , 1993, Neuroreport.

[21]  D. van der Spoel,et al.  A temperature predictor for parallel tempering simulations. , 2008, Physical chemistry chemical physics : PCCP.

[22]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[23]  Normand Mousseau,et al.  Coarse-grained protein molecular dynamics simulations. , 2007, The Journal of chemical physics.

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

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

[26]  Werner Krauth,et al.  Event-chain Monte Carlo algorithms for hard-sphere systems. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

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

[29]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[30]  R. Best,et al.  Free‐energy landscape of the GB1 hairpin in all‐atom explicit solvent simulations with different force fields: Similarities and differences , 2011, Proteins.

[31]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[32]  D. Kofke,et al.  Selection of temperature intervals for parallel-tempering simulations. , 2005, The Journal of chemical physics.

[33]  Matthias Troyer,et al.  Feedback-optimized parallel tempering Monte Carlo , 2006, cond-mat/0602085.

[34]  C. Tsallis.  Possible generalization of Boltzmann-Gibbs statistics , 1988 .

[35]  D. Kofke.  Erratum: "On the acceptance probability of replica-exchange Monte Carlo trials" [J. Chem. Phys. 117, 6911 (2002)] , 2004 .

[36]  Michael Nilges,et al.  Bayesian inference applied to macromolecular structure determination. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  P. Derreumaux,et al.  Computer simulations aimed at structure prediction of supersecondary motifs in proteins , 2001, Proteins.

[38]  S. Whittington,et al.  Monte carlo study of the interacting self-avoiding walk model in three dimensions , 1996 .

[39]  A. Schug,et al.  All‐atom folding of the three‐helix HIV accessory protein with an adaptive parallel tempering method , 2004, Proteins.

[40]  Michael W. Deem,et al.  Strict detailed balance is unnecessary in Monte Carlo simulation , 1999 .

[41]  U. Hansmann.  Temperature random walk sampling of protein configurations , 2010 .

[42]  Radford M. Neal,et al.  ANALYSIS OF A NONREVERSIBLE MARKOV CHAIN SAMPLER , 2000 .

[43]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[44]  M. Gallagher,et al.  A specific amyloid-β protein assembly in the brain impairs memory , 2006, Nature.

[45]  Jianyin Shao,et al.  Clustering Molecular Dynamics Trajectories: 1. Characterizing the Performance of Different Clustering Algorithms. , 2007, Journal of chemical theory and computation.

[46]  H. Takayama,et al.  APPLICATION OF AN EXTENDED ENSEMBLE METHOD TO SPIN GLASSES , 1996 .

[47]  A. Gronenborn,et al.  A novel, highly stable fold of the immunoglobulin binding domain of streptococcal protein G. , 1993, Science.