Incorporating solvent and ion screening into molecular dynamics using the finite‐difference Poisson–Boltzmann method

The finite difference method for solving the Poisson–Boltzmann equation is used to calculate the reaction field acting on a macromolecular solute due to the surrounding water and ions. Comparisons with analytical test cases indicate that the solvation forces can be calculated rapidly and accurately with this method. These forces act to move charged solute atoms towards the solvent where they are better solvated, and to screen interactions between charges. A way of combining such calculations with conventional molecular dynamics force fields is proposed which requires little modification of existing molecular dynamics programs. Simulations on the alanine dipeptide show that solvent forces affect the conformational dynamics by reducing the preference for internal H‐bonding forms, increasing the R‐alpha helix preference and reducing transition barriers. These solvent effects are similar to previous explicit solvent simulations, but require little more computation than vacuum simulations. The method should scale up with little increase in computational cost to larger molecules such as proteins and nucleic acids.

[1]  G L Seibel,et al.  A molecular dynamics simulation of double-helical B-DNA including counterions and water. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[2]  T. L. Hill,et al.  Influence of Electrolyte on effective Dielectric constants in Enzymes, proteins and other molecules , 1956 .

[3]  H. Berendsen,et al.  Inclusion of reaction fields in molecular dynamics. Application to liquid water , 1978 .

[4]  M Levitt,et al.  Computer simulation of DNA double-helix dynamics. , 1983, Cold Spring Harbor symposia on quantitative biology.

[5]  H. Scheraga,et al.  Prediction of the native conformation of a polypeptide by a statistical‐mechanical procedure. I. Backbone structure of enkephalin , 1985, Biopolymers.

[6]  Martin Neumann,et al.  Dielectric relaxation in water. Computer simulations with the TIP4P potential , 1986 .

[7]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[8]  J. Mccammon,et al.  Intramolecular flexibility in phenylalanine transfer RNA , 1981, Nature.

[9]  B. Hingerty,et al.  Carcinogen–base stacking and base–base stacking in dCpdG modified by (+) and (−) anti‐BPDE , 1985, Biopolymers.

[10]  Kurt Wüthrich,et al.  Secondary structure of the α-amylase polypeptide inhibitor Tendamistat from Streptomyces tendae determined in solution by 1H nuclear magnetic resonance , 1985 .

[11]  P. Kollman,et al.  Water–water and water–ion potential functions including terms for many body effects , 1985 .

[12]  Michiel Sprik,et al.  A polarizable model for water using distributed charge sites , 1988 .

[13]  Barry Honig,et al.  Focusing of electric fields in the active site of Cu‐Zn superoxide dismutase: Effects of ionic strength and amino‐acid modification , 1986, Proteins.

[14]  M. Karplus,et al.  Active site dynamics in protein molecules: A stochastic boundary molecular‐dynamics approach , 1985, Biopolymers.

[15]  B Jayaram,et al.  The electrostatic potential of B‐DNA , 1989, Biopolymers.

[16]  S. Harvey.  Treatment of electrostatic effects in macromolecular modeling , 1989, Proteins.

[17]  Benoît Roux,et al.  Molecular basis for the Born model of ion solvation , 1990 .

[18]  J. Warwicker,et al.  Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. , 1982, Journal of molecular biology.

[19]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[20]  E. Clementi,et al.  Molecular dynamics simulation of liquid water with three-body forces included , 1986 .

[21]  C. Tanford,et al.  Theory of Protein Titration Curves. I. General Equations for Impenetrable Spheres , 1957 .

[22]  R. Watts.  Electric polarisation of water: Monte Carlo studies , 1981 .

[23]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[24]  Barry Honig,et al.  Reevaluation of the Born model of ion hydration , 1985 .

[25]  B. Honig,et al.  On the calculation of electrostatic interactions in proteins. , 1985, Journal of molecular biology.

[26]  R. H. Ritchie,et al.  Dielectric effects in biopolymers: The theory of ionic saturation revisited , 1985 .

[27]  M Karplus,et al.  Picosecond dynamics of tyrosine side chains in proteins. , 1979, Biochemistry.

[28]  Bernard Pettitt,et al.  The potential of mean force surface for the alanine dipeptide in aqueous solution: a theoretical approach , 1985 .

[29]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[30]  J. Moult,et al.  Computer simulation of the solvent structure around biological macromolecules , 1978, Nature.

[31]  Barry Honig,et al.  Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation , 1990 .

[32]  J. E. Quinn,et al.  Cooperative effects in simulated water , 1979, Nature.

[33]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[34]  J. A. McCammon,et al.  Solving the finite difference linearized Poisson‐Boltzmann equation: A comparison of relaxation and conjugate gradient methods , 1989 .

[35]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[36]  Bhyravabhotla Jayaram,et al.  Free energy calculations of ion hydration: an analysis of the Born model in terms of microscopic simulations , 1989 .

[37]  Martin Karplus,et al.  SOLVATION. A MOLECULAR DYNAMICS STUDY OF A DIPEPTIDE IN WATER. , 1979 .

[38]  Alexander A. Rashin,et al.  A simple method for the calculation of hydration enthalpies of polar molecules with arbitrary shapes , 1987 .

[39]  D. Ritson,et al.  The Dielectric Properties of Water and Heavy Water , 1948 .

[40]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[41]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[42]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[43]  A. Warshel,et al.  How do serine proteases really work? , 1989, Biochemistry.

[44]  H. Berendsen,et al.  A consistent empirical potential for water–protein interactions , 1984 .

[45]  Mihaly Mezei,et al.  Monte Carlo determination of the free energy and internal energy of hydration for the Ala dipeptide at 25.degree.C , 1985 .

[46]  J Hermans,et al.  Microfolding: Conformational probability map for the alanine dipeptide in water from molecular dynamics simulations , 1988, Proteins.

[47]  Harold L. Friedman,et al.  Image approximation to the reaction field , 1975 .

[48]  William L. Jorgensen,et al.  Energy component analysis for dilute aqueous solutions of lithium(1+), sodium(1+), fluoride(1-), and chloride(1-) ions , 1984 .

[49]  H. Berendsen,et al.  Computer simulation of the dynamics of hydrated protein crystals and its comparison with x-ray data. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[50]  William L. Jorgensen,et al.  Cis-trans energy difference for the peptide bond in the gas phase and in aqueous solution , 1988 .