Strategies to model the near‐solute solvent molecular density/polarization

The solvent molecular distribution significantly affects the behavior of the solute molecules and is thus important in studying many biological phenomena. It can be described by the solvent molecular density distribution, g, and the solvent electric dipole distribution, p. The g and p can be computed directly by counting the number of solvent molecules/dipoles in a microscopic volume centered at r during a simulation or indirectly from the mean force F and electrostatic field E acting on the solvent molecule at r, respectively. However, it is not clear how the g and p derived from simulations depend on the solvent molecular center or the solute charge and if the gF and pE computed from the mean force and electric field acting on the solvent molecule, respectively, could reproduce the corresponding g and p obtained by direct counting. Hence, we have computed g, p, gF, and pE using different water centers from simulations of a solute atom of varying charge solvated in TIP3P water. The results show that gF and pE can reproduce the g and p obtained using a given count center. This implies that rather than solving the coordinates of each water molecule by MD simulations, the distribution of water molecules could be indirectly obtained from analytical formulas for the mean force F and electrostatic field E acting on the solvent molecule at r. Furthermore, the dependence of the g and p distributions on the solute charge revealed provides an estimate of the change in g and p surrounding a biomolecule upon a change in its conformation. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009

[1]  M. Nakasako.  Water-protein interactions from high-resolution protein crystallography. , 2004, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[2]  B. Sung,et al.  Monte Carlo simulations and integral equation theory for the structure of telechelic polymers , 2003 .

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

[4]  D. Frenkel,et al.  Simulation of liquids and solids : molecular dynamics and Monte Carlo methods in statistical mechanics , 1987 .

[5]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[6]  F. Lado,et al.  An integral equation study of a simple point charge model of water , 1999 .

[7]  Nathan A. Baker,et al.  Solvation forces on biomolecular structures: A comparison of explicit solvent and Poisson–Boltzmann models , 2004, J. Comput. Chem..

[8]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. II. Solvation structure from the three-dimensional reference interaction site model approach, and comparison with simulations , 2000 .

[9]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

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

[11]  M. Taylor,et al.  A Born–Green–Yvon integral equation theory for self-interacting lattice polymers , 1998 .

[12]  Fumio Hirata,et al.  Solvation thermodynamics of protein studied by the 3D-RISM theory , 2004 .

[13]  P. Ball.  Water as an active constituent in cell biology. , 2008, Chemical reviews.

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

[15]  C. Lim,et al.  Nonconvergence of the solute potential in an infinite solvent and its implications in continuum models , 2002 .

[16]  Ruhong Zhou,et al.  Comment on: Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? The Poisson-Boltzmann equation by Zhou, R. et al. , 2004 .

[17]  K. Binder,et al.  Molecular Dynamics Simulations , 2003, cond-mat/0308148.

[18]  Benoît Roux,et al.  Grand canonical Monte Carlo simulations of water in protein environments. , 2004, The Journal of chemical physics.

[19]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[20]  D. C. Clary,et al.  The Water Dipole Moment in Water Clusters , 1997, Science.

[21]  R. Zhou.  Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

[22]  B. Roux,et al.  Solvation Free Energy of Polar and Nonpolar Molecules in Water: An Extended Interaction Site Integral Equation Theory in Three Dimensions , 2000 .

[23]  Michel Daune,et al.  Molecular Biophysics: Structures in Motion , 1999 .

[24]  H. Modarress,et al.  Solving integral equations for binary and ternary systems , 2002 .

[25]  Gerhard Hummer,et al.  Molecular Theories and Simulation of Ions and Polar Molecules in Water , 1998 .

[26]  Greg L. Hura,et al.  Determining the role of hydration forces in protein folding , 1999 .

[27]  B. Berne,et al.  Reply to Comment on “Electrostatic Potentials and Free Energies of Solvation of Polar and Charged Molecules” , 1998 .

[28]  S. Noskov,et al.  Free energy decomposition of protein-protein interactions. , 2001, Biophysical journal.

[29]  S. Ten-no,et al.  On the connection between the reference interaction site model integral equation theory and the partial wave expansion of the molecular Ornstein–Zernike equation , 1999 .

[30]  Gerhard Hummer,et al.  Electrostatic Potentials and Free Energies of Solvation of Polar and Charged Molecules , 1997, physics/9704020.

[31]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

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