Exploring Hamiltonian dielectric solvent molecular dynamics

Abstract Hamiltonian dielectric solvent (HADES) is a recent method [7] , [25] , which enables Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric continua. Sample simulations of an α -helical decapeptide with and without explicit solvent demonstrate the high efficiency of HADES-MD. Addressing the folding of this peptide by replica exchange MD we study the properties of HADES by comparing melting curves, secondary structure motifs and salt bridges with explicit solvent results. Despite the unoptimized ad hoc parametrization of HADES, calculated reaction field energies correlate well with numerical grid solutions of the dielectric Poisson equation.

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

[2]  J. Mccammon,et al.  Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation , 2000 .

[3]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[4]  P. Tavan,et al.  Electrostatics of proteins in dielectric solvent continua. I. Newton's third law marries qE forces. , 2007, The Journal of chemical physics.

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

[6]  W. V. van Gunsteren,et al.  A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations , 2001 .

[7]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[8]  Paul Tavan,et al.  Optimizing the Accuracy and Efficiency of Fast Hierarchical Multipole Expansions for MD Simulations. , 2012, Journal of chemical theory and computation.

[9]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[10]  Massimo Marchi,et al.  A dielectric continuum model of solvation for complex solutes , 2005, Comput. Phys. Commun..

[11]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[12]  P. Tavan,et al.  Continuum description of solvent dielectrics in molecular-dynamics simulations of proteins , 2003 .

[13]  Tomasz Grycuk,et al.  Deficiency of the Coulomb-field approximation in the generalized Born model: An improved formula for Born radii evaluation , 2003 .

[14]  Paul Tavan,et al.  Simulated Solute Tempering. , 2009, Journal of chemical theory and computation.

[15]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[16]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[17]  Eric Vanden-Eijnden,et al.  Second-order integrators for Langevin equations with holonomic constraints , 2006 .

[18]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[19]  Guo-Wei Wei,et al.  Multiscale molecular dynamics using the matched interface and boundary method , 2011, J. Comput. Phys..

[20]  V. Hornak,et al.  Investigation of Salt Bridge Stability in a Generalized Born Solvent Model. , 2006, Journal of chemical theory and computation.

[21]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[22]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[23]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[24]  David A. Case,et al.  Effective Born radii in the generalized Born approximation: The importance of being perfect , 2002, J. Comput. Chem..

[25]  M. Marchi,et al.  A dielectric continuum molecular dynamics method , 2001 .

[26]  James C. Phillips,et al.  Parallel Generalized Born Implicit Solvent Calculations with NAMD. , 2011, Journal of chemical theory and computation.

[27]  Graham R. Fleming,et al.  Femtosecond solvation dynamics of water , 1994, Nature.

[28]  V. Barone,et al.  Combining the Fluctuating Charge Method, Non-Periodic Boundary Conditions and Meta-Dynamics: Aqua Ions as case studies. , 2014, Journal of chemical theory and computation.

[29]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[30]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

[31]  Andrew T. Fenley,et al.  An analytical approach to computing biomolecular electrostatic potential. I. Derivation and analysis. , 2008, The Journal of chemical physics.

[32]  Paul Tavan,et al.  Electrostatics of proteins in dielectric solvent continua. II. Hamiltonian reaction field dynamics. , 2014, The Journal of chemical physics.

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

[34]  W. Im,et al.  Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation , 1998 .

[35]  Paul Tavan,et al.  Electrostatics of proteins in dielectric solvent continua. I. An accurate and efficient reaction field description. , 2014, The Journal of chemical physics.

[36]  P. Tavan,et al.  Coupling density functional theory to polarizable force fields for efficient and accurate Hamiltonian molecular dynamics simulations. , 2013, The Journal of chemical physics.

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

[38]  Paul Tavan,et al.  A fast multipole method combined with a reaction field for long-range electrostatics in molecular dynamics simulations: The effects of truncation on the properties of water , 2003 .

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

[40]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[41]  P. Tavan,et al.  Including the Dispersion Attraction into Structure-Adapted Fast Multipole Expansions for MD Simulations. , 2014, Journal of chemical theory and computation.

[42]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Johannes Buchner,et al.  Protein folding handbook , 2005 .

[44]  H. Nymeyer,et al.  Simulation of the folding equilibrium of α-helical peptides: A comparison of the generalized Born approximation with explicit solvent , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Rudolf Reichold.  Rechnergestützte Beschreibung der Struktur und Dynamik von Peptiden und ihren Bausteinen , 2009 .

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

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

[48]  P. Tavan,et al.  Electrostatics of proteins in dielectric solvent continua. II. First applications in molecular dynamics simulations. , 2007, The Journal of chemical physics.

[49]  Paul Tavan,et al.  Optimal temperature ladders in replica exchange simulations , 2009 .