Enhanced Conformational Sampling Using Replica Exchange with Collective-Variable Tempering

The computational study of conformational transitions in RNA and proteins with atomistic molecular dynamics often requires suitable enhanced sampling techniques. We here introduce a novel method where concurrent metadynamics are integrated in a Hamiltonian replica-exchange scheme. The ladder of replicas is built with different strengths of the bias potential exploiting the tunability of well-tempered metadynamics. Using this method, free-energy barriers of individual collective variables are significantly reduced compared with simple force-field scaling. The introduced methodology is flexible and allows adaptive bias potentials to be self-consistently constructed for a large number of simple collective variables, such as distances and dihedral angles. The method is tested on alanine dipeptide and applied to the difficult problem of conformational sampling in a tetranucleotide.

[1]  Guanghong Wei,et al.  Effect of the disulfide bond on the monomeric structure of human amylin studied by combined Hamiltonian and temperature replica exchange molecular dynamics simulations. , 2010, The journal of physical chemistry. B.

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

[3]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[4]  W. L. Jorgensen Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water , 2002 .

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

[6]  Daniel R. Roe,et al.  Evaluation of Enhanced Sampling Provided by Accelerated Molecular Dynamics with Hamiltonian Replica Exchange Methods , 2014, The journal of physical chemistry. B.

[7]  Rafael Brüschweiler,et al.  In silico relationship between configurational entropy and soft degrees of freedom in proteins and peptides. , 2009, Physical review letters.

[8]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[9]  B. Berne,et al.  Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). , 2011, The journal of physical chemistry. B.

[10]  Michele Parrinello,et al.  Variational approach to enhanced sampling and free energy calculations. , 2014, Physical review letters.

[11]  Ruhong Zhou,et al.  Hydrophobic aided replica exchange: an efficient algorithm for protein folding in explicit solvent. , 2006, The journal of physical chemistry. B.

[12]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[13]  Massimiliano Bonomi,et al.  Efficient Simulation of Explicitly Solvated Proteins in the Well-Tempered Ensemble. , 2012, Journal of chemical theory and computation.

[14]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[15]  R. Nussinov,et al.  The role of dynamic conformational ensembles in biomolecular recognition. , 2009, Nature chemical biology.

[16]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[17]  R. Broglia,et al.  Exploring the protein G helix free‐energy surface by solute tempering metadynamics , 2007, Proteins.

[18]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[19]  Paul Tavan,et al.  Efficiency reduction and pseudo-convergence in replica exchange sampling of peptide folding–unfolding equilibria , 2008 .

[20]  Niel M. Henriksen,et al.  Reliable oligonucleotide conformational ensemble generation in explicit solvent for force field assessment using reservoir replica exchange molecular dynamics simulations. , 2013, The journal of physical chemistry. B.

[21]  M. Zacharias,et al.  Enhanced conformational sampling of nucleic acids by a new Hamiltonian replica exchange molecular dynamics approach. , 2009, The Journal of chemical physics.

[22]  M. Zacharias,et al.  Influence of 8-oxoguanosine on the fine structure of DNA studied with biasing-potential replica exchange simulations. , 2013, Biophysical journal.

[23]  Giovanni Bussi,et al.  Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration , 2013, Entropy.

[24]  Angel E García,et al.  High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations , 2013, Proceedings of the National Academy of Sciences.

[25]  A. Laio,et al.  A bias-exchange approach to protein folding. , 2007, The journal of physical chemistry. B.

[26]  T. P. Straatsma,et al.  Treatment of rotational isomeric states. III. The use of biasing potentials , 1994 .

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

[28]  A. Voter Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

[29]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

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

[31]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

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

[33]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[34]  J. P. Grossman,et al.  Biomolecular simulation: a computational microscope for molecular biology. , 2012, Annual review of biophysics.

[35]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

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

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

[38]  I. Tavernelli,et al.  A Novel Hamiltonian Replica Exchange MD Protocol to Enhance Protein Conformational Space Sampling. , 2006, Journal of chemical theory and computation.

[39]  J. Šponer,et al.  Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles , 2011, Journal of chemical theory and computation.

[40]  M. Karplus,et al.  A hierarchy of timescales in protein dynamics is linked to enzyme catalysis , 2007, Nature.

[41]  Christophe Chipot,et al.  Enhanced Sampling of Multidimensional Free-Energy Landscapes Using Adaptive Biasing Forces , 2010, SIAM Journal on Applied Mathematics.

[42]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[43]  M. Parrinello,et al.  Metadynamics with Adaptive Gaussians. , 2012, Journal of chemical theory and computation.

[44]  Giovanni Bussi,et al.  Accurate multiple time step in biased molecular simulations. , 2014, Journal of chemical theory and computation.

[45]  Martin Zacharias,et al.  Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations , 2009, Proteins.

[46]  Sushil Kumar Mishra,et al.  Enhanced conformational sampling of carbohydrates by Hamiltonian replica-exchange simulation. , 2014, Glycobiology.

[47]  Andrew E. Torda,et al.  Local elevation: A method for improving the searching properties of molecular dynamics simulation , 1994, J. Comput. Aided Mol. Des..

[48]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[49]  Hugh Nymeyer,et al.  How Efficient Is Replica Exchange Molecular Dynamics? An Analytic Approach. , 2008, Journal of chemical theory and computation.

[50]  Martin Zacharias,et al.  Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing‐potential , 2006, Proteins.

[51]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[52]  C. Brooks,et al.  Hierarchy of RNA functional dynamics. , 2014, Annual review of biochemistry.

[53]  Michele Parrinello,et al.  Well-tempered metadynamics converges asymptotically. , 2014, Physical review letters.

[54]  D. Turner,et al.  Benchmarking AMBER Force Fields for RNA: Comparisons to NMR Spectra for Single-Stranded r(GACC) Are Improved by Revised χ Torsions , 2011, The journal of physical chemistry. B.

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

[56]  J. VandeVondele,et al.  Canonical Adiabatic Free Energy Sampling (CAFES): A Novel Method for the Exploration of Free Energy Surfaces , 2002 .

[57]  Adrian E. Roitberg,et al.  Multidimensional Replica Exchange Molecular Dynamics Yields a Converged Ensemble of an RNA Tetranucleotide , 2013, Journal of chemical theory and computation.

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

[59]  Massimiliano Bonomi,et al.  PLUMED 2: New feathers for an old bird , 2013, Comput. Phys. Commun..

[60]  Stephen Neidle,et al.  Principles of nucleic acid structure , 2007 .

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

[62]  Giovanni Bussi,et al.  Hamiltonian replica exchange in GROMACS: a flexible implementation , 2013, 1307.5144.

[63]  Amedeo Caflisch,et al.  Calculation of conformational transitions and barriers in solvated systems: Application to the alanine dipeptide in water , 1999 .

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

[65]  M. Tuckerman,et al.  An Adiabatic Molecular Dynamics Method for the Calculation of Free Energy Profiles , 2002 .

[66]  E. Vanden-Eijnden,et al.  Single-sweep methods for free energy calculations. , 2007, The Journal of chemical physics.