Simulated tempering distributed replica sampling: A practical guide to enhanced conformational sampling

Simulated tempering distributed replica sampling (STDR) is a generalized-ensemble method designed specifically for simulations of large molecular systems on shared and heterogeneous computing platforms [Rauscher, Neale and Pomes (2009) J. Chem. Theor. Comput. 5, 2640]. The STDR algorithm consists of an alternation of two steps: (1) a short molecular dynamics (MD) simulation; and (2) a stochastic temperature jump. Repeating these steps thousands of times results in a random walk in temperature, which allows the system to overcome energetic barriers, thereby enhancing conformational sampling. The aim of the present paper is to provide a practical guide to applying STDR to complex biomolecular systems. We discuss the details of our STDR implementation, which is a highly-parallel algorithm designed to maximize computational efficiency while simultaneously minimizing network communication and data storage requirements. Using a 35-residue disordered peptide in explicit water as a test system, we characterize the efficiency of the STDR algorithm with respect to both diffusion in temperature space and statistical convergence of structural properties. Importantly, we show that STDR provides a dramatic enhancement of conformational sampling compared to a canonical MD simulation.

[1]  Daniel M Zuckerman,et al.  Ensemble-based convergence analysis of biomolecular trajectories. , 2006, Biophysical journal.

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

[3]  Patrik Edén,et al.  Accounting for one-channel depletion improves missing value imputation in 2-dye microarray data , 2008, BMC Genomics.

[4]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[5]  J. S. Sodhi,et al.  Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. , 2004, Journal of molecular biology.

[6]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

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

[8]  Fred W. Keeley,et al.  Structural disorder and dynamics of elastin. , 2010, Biochemistry and cell biology = Biochimie et biologie cellulaire.

[9]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[10]  Sarah Rauscher,et al.  Molecular simulations of protein disorder. , 2010, Biochemistry and cell biology = Biochimie et biologie cellulaire.

[11]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[12]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

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

[14]  Chris Neale,et al.  Simulated Tempering Distributed Replica Sampling, Virtual Replica Exchange, and Other Generalized-Ensemble Methods for Conformational Sampling. , 2009, Journal of chemical theory and computation.

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

[16]  K. Plaxco,et al.  Using protein folding rates to test protein folding theories. , 2003, Annual review of biochemistry.

[17]  Régis Pomès,et al.  Distributed Replica Sampling. , 2006, Journal of chemical theory and computation.

[18]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .

[19]  Régis Pomès,et al.  Calculation of absolute protein-ligand binding free energy using distributed replica sampling. , 2008, The Journal of chemical physics.

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

[21]  Sarah Rauscher,et al.  Proline and glycine control protein self-organization into elastomeric or amyloid fibrils. , 2006, Structure.

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

[23]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[24]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

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

[26]  Christophe Chipot,et al.  Comprar Free Energy Calculations · Theory and Applications in Chemistry and Biology | Chipot, Christophe | 9783540736172 | Springer , 2007 .

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

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

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