Smart resolution replica exchange: an efficient algorithm for exploring complex energy landscapes.

A coarse-grained representation of a condensed phase system can significantly reduce the number of system degrees of freedom, making coarse-grained simulations very computationally efficient. Moreover, coarse graining can smoothen the free energy landscape of the system. Thus coarse-grained dynamics is usually faster than its fully atomistic counterpart. In this work, the smart resolution replica exchange method is introduced that incorporates the information from coarse-grained simulations into atomistic simulations in order to accelerate the sampling of rough, complex atomistic energy landscapes. Within this methodology, interactions between particles are defined by a potential energy that interpolates between a fully atomistic potential and a fully coarse-grained effective potential according to a parameter lambda. Instead of exchanging the configurations from neighboring resolutions directly, as has been done in the resolution replica exchange methods [E. Lyman et al., Phys. Rev. Lett. 96, 028105 (2006); M. Christen and W. F. v. Gunsteren, J. Chem. Phys. 124, 154106 (2006)], the configuration described at the coarser resolution is first relaxed before an exchange is attempted, similar to the smart walking method [R. Zhou and B. J. Berne, J. Chem. Phys. 107, 9185 (1997)]. This approach greatly increases the acceptance ratio of exchange and only two replicas, one at the atomistic level and one at the coarse-grained level, are usually required (although more can be implemented if desired). This new method can approximately obtain the correct canonical sampling if the exchange interval is sufficiently large to allow the system to explore the local energy landscape. The method is demonstrated for a two-dimensional model system, where the ideal population distribution can be recovered, and also for an alanine polypeptide (Ala(15)) model with explicit water, where its native structure, an alpha helix, is obtained from the extended structure within 1 ns.

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

[2]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[3]  Gregory A Voth,et al.  Multiscale coarse graining of liquid-state systems. , 2005, The Journal of chemical physics.

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

[5]  M. Klein,et al.  Computer simulation studies of biomembranes using a coarse grain model , 2002 .

[6]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

[7]  D. L. Freeman,et al.  Reducing Quasi-Ergodic Behavior in Monte Carlo Simulations by J-Walking: Applications to Atomic Clusters , 1990 .

[8]  Bruce J. Berne,et al.  Multicanonical jump walking: A method for efficiently sampling rough energy landscapes , 1999 .

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

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

[11]  M. Karplus,et al.  Effective energy functions for protein structure prediction. , 2000, Current opinion in structural biology.

[12]  B. Berne,et al.  A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems , 1995 .

[13]  D. Thirumalai,et al.  Anisotropic coarse-grained statistical potentials improve the ability to identify nativelike protein structures , 2003, physics/0302009.

[14]  Janna K. Maranas,et al.  Speed up of dynamic observables in coarse-grained molecular-dynamics simulations of unentangled polymers. , 2005, The Journal of chemical physics.

[15]  M. Hoare,et al.  Statistical mechanics and morphology of very small atomic clusters , 1976 .

[16]  Marek Cieplak,et al.  Molecular dynamics of folding of secondary structures in Go-type models of proteins , 1999, cond-mat/9911488.

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

[18]  N. Go,et al.  Noninteracting local‐structure model of folding and unfolding transition in globular proteins. I. Formulation , 1981, Biopolymers.

[19]  C. Hall,et al.  α‐Helix formation: Discontinuous molecular dynamics on an intermediate‐resolution protein model , 2001, Proteins.

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

[21]  Wilfred F van Gunsteren,et al.  Multigraining: an algorithm for simultaneous fine-grained and coarse-grained simulation of molecular systems. , 2006, The Journal of chemical physics.

[22]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[23]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[24]  Scott Brown,et al.  Cool walking: A new Markov chain Monte Carlo sampling method , 2003, J. Comput. Chem..

[25]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[26]  A. Liwo,et al.  A united‐residue force field for off‐lattice protein‐structure simulations. II. Parameterization of short‐range interactions and determination of weights of energy terms by Z‐score optimization , 1997 .

[27]  Bruce J. Berne,et al.  Quantum path minimization: An efficient method for global optimization , 2003 .

[28]  Qiang Shi,et al.  Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel. , 2006, The journal of physical chemistry. B.

[29]  Edward Lyman,et al.  Resolution exchange simulation. , 2005, Physical review letters.

[30]  Gregory A Voth,et al.  Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching. , 2004, The Journal of chemical physics.

[31]  Edward Lyman,et al.  Resolution Exchange Simulation with Incremental Coarsening. , 2006, Journal of chemical theory and computation.

[32]  B. Berne,et al.  Smart walking: A new method for Boltzmann sampling of protein conformations , 1997 .

[33]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[34]  A. Voter,et al.  Smart Darting Monte Carlo , 2001 .

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

[36]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[37]  Sergei Izvekov,et al.  Multiscale Coarse-Graining of Mixed Phospholipid/Cholesterol Bilayers. , 2006, Journal of chemical theory and computation.

[38]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[39]  R. Luo,et al.  Overcoming entropic barrier with coupled sampling at dual resolutions. , 2005, The Journal of chemical physics.

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

[41]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

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

[43]  Michael L. Klein,et al.  Coarse grain models and the computer simulation of soft materials , 2004 .

[44]  T. Head-Gordon,et al.  Minimalist models for protein folding and design. , 2003, Current opinion in structural biology.

[45]  Gregory A Voth,et al.  Multiscale coarse-graining of ionic liquids. , 2006, The journal of physical chemistry. B.

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

[47]  J. Schofield,et al.  Extended state-space Monte Carlo methods. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  D Thirumalai,et al.  Development of novel statistical potentials for protein fold recognition. , 2004, Current opinion in structural biology.

[49]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .