On the use of mass scaling for stable and efficient simulated tempering with molecular dynamics

Simulated tempering (ST) is a generalized‐ensemble algorithm that employs trajectories exploring a range of temperatures to effectively sample rugged energy landscapes. When implemented using the molecular dynamics method, ST can require the use of short time steps for ensuring the stability of trajectories at high temperatures. To address this shortcoming, a mass‐scaling ST (MSST) method is presented in which the particle mass is scaled in proportion to the temperature. Mass scaling in the MSST method leads to velocity distributions that are independent of temperature and eliminates the need for velocity scaling after the accepted temperature updates that are required in conventional ST simulations. The homogeneity in time scales with changing temperature improves the stability of simulations and allows for the use of longer time steps at high temperatures. As a result, the MSST is found to be more efficient than the standard ST method, particularly for cases in which a large temperature range is employed. © 2016 Wiley Periodicals, Inc.

[1]  Joan-Emma Shea,et al.  Computational studies of protein aggregation: methods and applications. , 2015, Annual review of physical chemistry.

[2]  C. Brooks Computer simulation of liquids , 1989 .

[3]  Bernd A. Berg Markov Chain Monte Carlo Simulations and Their Statistical Analysis , 2004 .

[4]  Charles H. Bennett,et al.  Mass tensor molecular dynamics , 1975 .

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

[6]  H. C. Andersen Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[7]  Vijay S. Pande,et al.  Accelerating molecular dynamic simulation on graphics processing units , 2009, J. Comput. Chem..

[8]  Yuko Okamoto,et al.  The Generalized-Ensemble Approach for Protein Folding Simulations , 1999 .

[9]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

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

[11]  Markus Deserno,et al.  Membrane composition-mediated protein-protein interactions , 2008, Biointerphases.

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

[13]  Takuya Takahashi,et al.  Mass-scaling replica-exchange molecular dynamics optimizes computational resources with simpler algorithm. , 2014, The Journal of chemical physics.

[14]  J. Straub,et al.  Statistical temperature molecular dynamics: application to coarse-grained beta-barrel-forming protein models. , 2007, The Journal of chemical physics.

[15]  W. Janke,et al.  Crossover scaling in the two-dimensional three-state Potts model , 2013, 1307.3865.

[16]  Yuko Okamoto,et al.  Simulated tempering and magnetizing: application of two-dimensional simulated tempering to the two-dimensional Ising model and its crossover. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[18]  Phuong H Nguyen Replica exchange simulation method using temperature and solvent viscosity. , 2010, The Journal of chemical physics.

[19]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[20]  John E Straub,et al.  Replica exchange statistical temperature Monte Carlo. , 2009, The Journal of chemical physics.

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

[22]  John E Straub,et al.  Exploring the solid-liquid phase change of an adapted Dzugutov model using generalized replica exchange method. , 2012, The journal of physical chemistry. B.

[23]  A. Roitberg,et al.  Long-Time-Step Molecular Dynamics through Hydrogen Mass Repartitioning. , 2015, Journal of chemical theory and computation.

[24]  Y. Okamoto Institute for Molecular Science,et al.  Replica-exchange simulated tempering method for simulations of frustrated systems , 2000 .

[25]  Diwakar Shukla,et al.  OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. , 2013, Journal of chemical theory and computation.

[26]  Yuko Okamoto,et al.  Analysis of helix-helix interactions of bacteriorhodopsin by replica-exchange simulations. , 2009, Biophysical journal.

[27]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[28]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[29]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

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

[31]  Yuko Okamoto,et al.  Replica-Exchange Molecular Dynamics Simulations for Various Constant Temperature Algorithms , 2010 .

[32]  Yuko Okamoto,et al.  Communication: Simulated tempering with fast on-the-fly weight determination. , 2013, The Journal of chemical physics.

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

[34]  G M Maggiora,et al.  Mass‐weighted molecular dynamics simulation of cyclic polypeptides , 1991, Biopolymers.

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

[36]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[37]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[38]  J. Straub,et al.  Generalized simulated tempering for exploring strong phase transitions. , 2010, The Journal of chemical physics.

[39]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[40]  Raymond D. Mountain,et al.  Quantative measure of efficiency of Monte Carlo simulations , 1994 .

[41]  Markus Deserno,et al.  Mesoscopic membrane physics: concepts, simulations, and selected applications. , 2009, Macromolecular rapid communications.

[42]  Sanghyun Park,et al.  Comparison of the serial and parallel algorithms of generalized ensemble simulations: an analytical approach. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  Vijay S Pande,et al.  Convergence of folding free energy landscapes via application of enhanced sampling methods in a distributed computing environment. , 2008, The Journal of chemical physics.

[44]  Hisashi Okumura,et al.  Liquid–vapor coexistence curves of several interatomic model potentials , 2000 .

[45]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[46]  V. Pande,et al.  Choosing weights for simulated tempering. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  J. Straub,et al.  Replica exchange statistical temperature molecular dynamics algorithm. , 2012, The journal of physical chemistry. B.

[48]  Yuko Okamoto,et al.  Multidimensional generalized-ensemble algorithms for complex systems. , 2009, The Journal of chemical physics.

[49]  C. Zhang,et al.  Comparison of sampling efficiency between simulated tempering and replica exchange. , 2008, The Journal of chemical physics.

[50]  John E Straub,et al.  Generalized replica exchange method. , 2010, The Journal of chemical physics.

[51]  Mark E Tuckerman,et al.  Enhanced conformational sampling of peptides via reduced side-chain and solvent masses. , 2010, The journal of physical chemistry. B.

[52]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[53]  Z. Tan,et al.  Theory of binless multi-state free energy estimation with applications to protein-ligand binding. , 2012, The Journal of chemical physics.

[54]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[55]  Eiji Tsuchida,et al.  Ab initio mass tensor molecular dynamics. , 2010, The Journal of chemical physics.

[56]  Y. Okamoto,et al.  Molecular dynamics, Langevin, and hybrid Monte Carlo simulations in multicanonical ensemble , 1996, physics/9710018.

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

[58]  S. Nosé A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[59]  Yuko Okamoto,et al.  From multidimensional replica-exchange method to multidimensional multicanonical algorithm and simulated tempering. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[60]  B. Mao,et al.  Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide. , 1991, Biophysical journal.

[61]  M. Parrinello,et al.  Accurate sampling using Langevin dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[62]  M G E da Luz,et al.  Comparing parallel- and simulated-tempering-enhanced sampling algorithms at phase-transition regimes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[64]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[65]  Vijay S. Pande,et al.  Efficient nonbonded interactions for molecular dynamics on a graphics processing unit , 2010, J. Comput. Chem..

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

[67]  J. Andrew McCammon,et al.  Mass and step length optimization for the calculation of equilibrium properties by molecular dynamics simulation , 1990 .

[68]  D. Landau,et al.  Determining the density of states for classical statistical models: a random walk algorithm to produce a flat histogram. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[69]  Yingkai Zhang,et al.  Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations , 2009, J. Comput. Chem..

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

[71]  Ayori Mitsutake,et al.  Simulated-tempering replica-exchange method for the multidimensional version. , 2009, The Journal of chemical physics.

[72]  Kurt Kremer,et al.  Tunable generic model for fluid bilayer membranes. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[73]  A. Kidera,et al.  Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of Peptides , 1997 .

[74]  Michael R. Shirts,et al.  Replica exchange and expanded ensemble simulations as Gibbs sampling: simple improvements for enhanced mixing. , 2011, The Journal of chemical physics.

[75]  Yuko Okamoto,et al.  Application of simulated tempering and magnetizing to a two-dimensional Potts model , 2012, 1212.3084.

[76]  Robert H. Swendsen,et al.  Replica Monte Carlo Simulation (Revisited) , 2005 .

[77]  Tetsuya Morishita,et al.  Decomposition-order effects of time integrator on ensemble averages for the Nosé-Hoover thermostat. , 2013, The Journal of chemical physics.

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

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