Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations

An approach is developed to obtain statistical properties similar to those of an infinite bulk system from computer simulations of a finite cluster. A rigorous theoretical formulation is given for the solvent boundary potential which takes the influence of the surrounding bulk into account. The solvent boundary potential is the configuration‐dependent solvation free energy of an effective cluster composed of an arbitrary solute and a finite number of explicit solvent molecules embedded inside a hard sphere of variable radius; the hard sphere does not act directly on the solute or the explicit solvent molecules, and its radius varies according to the instantaneous configurations. The formulation follows from an exact separation of the multidimensional configurational Boltzmann integral in terms of the solvent molecules nearest to the solute and the remaining bulk solvent molecules. An approximation to the solvent boundary potential is constructed for simulations of bulk water at constant pressure, includin...

[1]  J. Kirkwood,et al.  Statistical Mechanics of Liquid Solutions. , 1936 .

[2]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[3]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[4]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

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

[6]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[7]  J. Murrell,et al.  Molecular dynamics and chemical reactivity , 1977 .

[8]  John Burgess,et al.  Metal Ions in Solution , 1978 .

[9]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[10]  William L. Jorgensen,et al.  Monte Carlo simulation of n‐butane in water. Conformational evidence for the hydrophobic effect , 1982 .

[11]  J. Mccammon,et al.  Molecular dynamics with stochastic boundary conditions , 1982 .

[12]  M. Karplus,et al.  Deformable stochastic boundaries in molecular dynamics , 1983 .

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

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

[15]  J. Andrew McCammon,et al.  The structure of liquid water at an extended hydrophobic surface , 1984 .

[16]  M. Karplus,et al.  Stochastic boundary conditions for molecular dynamics simulations of ST2 water , 1984 .

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

[18]  A. Warshel,et al.  Polarization constraints in molecular dynamics simulation of aqueous solutions: The surface constraint all atom solvent (SCAAS) model , 1985 .

[19]  Alan C. Belch,et al.  Molecular dynamics simulations of tips2 water restricted by a spherical hydrophobic boundary , 1985 .

[20]  R. Zauhar,et al.  A new method for computing the macromolecular electric potential. , 1985, Journal of molecular biology.

[21]  T. Straatsma,et al.  Free energy of hydrophobic hydration: A molecular dynamics study of noble gases in water , 1986 .

[22]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[23]  B. Montgomery Pettitt,et al.  Alkali halides in water: Ion–solvent correlations and ion–ion potentials of mean force at infinite dilution , 1986 .

[24]  William L. Jorgensen,et al.  Use of statistical perturbation theory for computing solvent effects on molecular conformation: butane in water , 1987 .

[25]  J. Rullmann,et al.  Analysis of discrete and continuum dielectric models : Application to the calculation of protonation energies in solution , 1987 .

[26]  J. Valleau,et al.  Water‐like particles at surfaces. I. The uncharged, unpolarized surface , 1987 .

[27]  Harold A. Scheraga,et al.  Free energies of hydration of solute molecules. 1. Improvement of the hydration shell model by exact computations of overlapping volumes , 1987 .

[28]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[29]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[30]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

[31]  V. P. Chuprina,et al.  DNA B to D transition can be explained in terms of hydration economy of the minor groove atoms. , 1989, Journal of molecular biology.

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

[33]  Charles L. Brooks,et al.  The thermodynamics of solvophobic effects: A molecular‐dynamics study of n‐butane in carbon tetrachloride and water , 1990 .

[34]  C Sander,et al.  Excluded volume approximation to protein-solvent interaction. The solvent contact model. , 1990, Biophysical journal.

[35]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[36]  B. Montgomery Pettitt,et al.  Molecular dynamics with a variable number of molecules , 1991 .

[37]  J Tirado-Rives,et al.  Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water. , 1991, Biochemistry.

[38]  C. Brooks,et al.  Thermodynamics and mechanism of alpha helix initiation in alanine and valine peptides. , 1991, Biochemistry.

[39]  Anna Tempczyk,et al.  Electrostatic contributions to solvation energies: comparison of free energy perturbation and continuum calculations , 1991 .

[40]  S J Wodak,et al.  Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[41]  Kim A. Sharp,et al.  Incorporating solvent and ion screening into molecular dynamics using the finite‐difference Poisson–Boltzmann method , 1991 .

[42]  Charles L. Brooks,et al.  Conformational equilibrium in the alanine dipeptide in the gas phase and aqueous solution : a comparison of theoretical results , 1992 .

[43]  Joel S. Bader,et al.  Computer simulation study of the mean forces between ferrous and ferric ions in water , 1992 .

[44]  H J Berendsen,et al.  Molecular dynamics simulation of the stability of a 22‐residue α‐helix in water and 30% trifluoroethanol , 1993, Biopolymers.