w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations

Free energy governs the equilibrium extent of many biological processes. High barriers separating free energy minima often limit the sampling in molecular dynamics (MD) simulations, leading to inaccurate free energies. Here, we demonstrate enhanced sampling and improved free energy calculations, relative to conventional MD, using windowed accelerated MD within a Hamiltonian replica exchange framework (w-REXAMD). We show that for a case in which multiple conformations are separated by large free energy barriers, w-REXAMD is a useful enhanced sampling technique, without any necessary reweighting.

[1]  J. Onuchic,et al.  Funnels, pathways, and the energy landscape of protein folding: A synthesis , 1994, Proteins.

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

[3]  G. Bouvignies,et al.  Exploring multiple timescale motions in protein GB3 using accelerated molecular dynamics and NMR spectroscopy. , 2007, Journal of the American Chemical Society.

[4]  Wei Yang,et al.  Simulated scaling method for localized enhanced sampling and simultaneous "alchemical" free energy simulations: a general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations. , 2007, The Journal of chemical physics.

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

[6]  Donald Hamelberg,et al.  Estimating kinetic rates from accelerated molecular dynamics simulations: alanine dipeptide in explicit solvent as a case study. , 2007, The Journal of chemical physics.

[7]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[8]  J Andrew McCammon,et al.  Conformational selection in G-proteins: lessons from Ras and Rho. , 2010, Biophysical journal.

[9]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[10]  Ray Luo,et al.  Implicit nonpolar solvent models. , 2007, The journal of physical chemistry. B.

[11]  Michelle Zhou,et al.  Large-Scale Conformational Changes of Trypanosoma cruzi Proline Racemase Predicted by Accelerated Molecular Dynamics Simulation , 2011, PLoS Comput. Biol..

[12]  J. Andrew McCammon,et al.  Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration , 2008, Journal of chemical theory and computation.

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

[14]  César Augusto F. de Oliveira,et al.  Coupling Accelerated Molecular Dynamics Methods with Thermodynamic Integration Simulations , 2008, Journal of chemical theory and computation.

[15]  David L. Mobley,et al.  Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations , 2010, J. Comput. Aided Mol. Des..

[16]  Alan E. Mark,et al.  Estimating the Relative Free Energy of Different Molecular States with Respect to a Single Reference State , 1996 .

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

[18]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[19]  Levi C. T. Pierce,et al.  On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations , 2011, Journal of chemical theory and computation.

[20]  Wei Yang,et al.  A divide-and-conquer strategy to improve diffusion sampling in generalized ensemble simulations. , 2008, The Journal of chemical physics.

[21]  Levi C. T. Pierce,et al.  Accelerating chemical reactions: exploring reactive free-energy surfaces using accelerated ab initio molecular dynamics. , 2011, The Journal of chemical physics.

[22]  Levi C. T. Pierce,et al.  Routine Access to Millisecond Time Scale Events with Accelerated Molecular Dynamics , 2012, Journal of chemical theory and computation.

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

[24]  Maxim V Fedorov,et al.  Toward a universal model to calculate the solvation thermodynamics of druglike molecules: the importance of new experimental databases. , 2011, Molecular pharmaceutics.

[25]  William Sinko,et al.  Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations , 2011, Journal of chemical theory and computation.

[26]  David Chandler,et al.  Barrier crossings:. classical theory of rare but important events , 1998 .

[27]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[28]  J Andrew McCammon,et al.  Statistical mechanics and molecular dynamics in evaluating thermodynamic properties of biomolecular recognition , 2011, Quarterly Reviews of Biophysics.

[29]  Rommie E. Amaro,et al.  Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy , 2010, Proteins.

[30]  Donald Hamelberg,et al.  A statistical analysis of the precision of reweighting-based simulations. , 2008, The Journal of chemical physics.

[31]  Donald Hamelberg,et al.  Relating kinetic rates and local energetic roughness by accelerated molecular-dynamics simulations. , 2005, The Journal of chemical physics.