Grand canonical Monte Carlo simulations of water in protein environments.

The grand canonical simulation algorithm is considered as a general methodology to sample the configuration of water molecules confined within protein environments. First, the probability distribution of the number of water molecules and their configuration in a region of interest for biochemical simulations, such as the active site of a protein, is derived by considering a finite subvolume in open equilibrium with a large system serving as a bulk reservoir. It is shown that the influence of the bulk reservoir can be represented as a many-body potential of mean force acting on the atoms located inside the subvolume. The grand canonical Monte Carlo (GCMC) algorithm, augmented by a number of technical advances to increase the acceptance of insertion attempts, is implemented, and tested for simple systems. In particular, the method is illustrated in the case of a pure water box with periodic boundary conditions. In addition, finite spherical systems of pure water and containing a dialanine peptide, are simulated with GCMC while the influence of the surrounding infinite bulk is incorporated using the generalized solvent boundary potential [W. Im, S. Berneche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)]. As a last illustration of water confined in the interior of a protein, the hydration of the central cavity of the KcsA potassium channel is simulated.

[1]  A. Warshel,et al.  Energetics of ion permeation through membrane channels. Solvation of Na+ by gramicidin A. , 1989, Biophysical journal.

[2]  L. Pratt.  Molecular theory of hydrophobic effects: "She is too mean to have her name repeated.". , 2001, Annual review of physical chemistry.

[3]  W. Im,et al.  Generalized solvent boundary potential for computer simulations , 2001 .

[4]  M. Mezei,et al.  Grand canonical ensemble Monte Carlo simulation of a lipid bilayer using extension biased rotations , 1999 .

[5]  A. Panagiotopoulos,et al.  Application of excluded volume map sampling to phase equilibrium calculations in the Gibbs ensemble , 1990 .

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

[7]  Gerhard Hummer,et al.  New perspectives on hydrophobic effects , 2000 .

[8]  Mihaly Mezei,et al.  Use of the Grand Canonical Ensemble in Potential of Mean Force Calculations , 1996 .

[9]  Wei Yang,et al.  Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. , 2004, The Journal of chemical physics.

[10]  L. Scriven,et al.  Efficient molecular simulation of chemical potentials , 1989 .

[11]  R. MacKinnon,et al.  The cavity and pore helices in the KcsA K+ channel: electrostatic stabilization of monovalent cations. , 1999, Science.

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

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

[14]  R. MacKinnon,et al.  Chemistry of ion coordination and hydration revealed by a K+ channel–Fab complex at 2.0 Å resolution , 2001, Nature.

[15]  G. N. Patey,et al.  A CONFIGURATION BIAS MONTE CARLO METHOD FOR WATER , 1995 .

[16]  D. Beglov,et al.  Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations , 1994 .

[17]  Roderick MacKinnon,et al.  Energetic optimization of ion conduction rate by the K+ selectivity filter , 2001, Nature.

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

[19]  J. Mccammon,et al.  Efficient trajectory simulation methods for diffusional barrier crossing processes , 1980 .

[20]  A. Thompson,et al.  Direct molecular simulation of gradient-driven diffusion of large molecules using constant pressure , 1999 .

[21]  Jonathan W. Essex,et al.  An empirical boundary potential for water droplet simulations , 1995, J. Comput. Chem..

[22]  P. Pohl,et al.  Massively parallel molecular dynamics simulation of gas permeation across porous silica membranes , 1999 .

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

[24]  B. Roux,et al.  Energetics of ion conduction through the K + channel , 2022 .

[25]  J A McCammon,et al.  Theory of biomolecular recognition. , 1998, Current opinion in structural biology.

[26]  Benoît Roux,et al.  Calculating excess chemical potentials using dynamic simulations in the fourth dimension , 1999 .

[27]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

[28]  D. Frenkel,et al.  Unexpected length dependence of the solubility of chain molecules , 1992 .

[29]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[30]  William L. Jorgensen,et al.  Free energy of TIP4P water and the free energies of hydration of CH4 and Cl- from statistical perturbation theory , 1989 .

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

[32]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[33]  A. Warshel,et al.  Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor. , 1985, Journal of molecular biology.

[34]  M. Mezei.  Excess free energy of different water models computed by Monte Carlo methods , 1982 .

[35]  J. Hermans,et al.  Excess free energy of liquids from molecular dynamics simulations. Application to water models. , 1988, Journal of the American Chemical Society.

[36]  Robert M Stroud,et al.  Glycerol facilitator GlpF and the associated aquaporin family of channels. , 2003, Current opinion in structural biology.

[37]  Mihaly Mezei,et al.  Grand Canonical Monte Carlo Simulation of Water Positions in Crystal Hydrates , 1994 .

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

[39]  Benoît Roux,et al.  Electrostatic free energy calculations using the generalized solvent boundary potential method , 2002 .

[40]  Barry Honig.  Theory and simulation , 1995 .

[41]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[42]  D. Chandler,et al.  Theory of the hydrophobic effect , 1977 .

[43]  Franci Merzel,et al.  Is the first hydration shell of lysozyme of higher density than bulk water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  B. Roux,et al.  Molecular dynamics of the KcsA K(+) channel in a bilayer membrane. , 2000, Biophysical journal.

[45]  Mihaly Mezei,et al.  A cavity-biased (T, V, μ) Monte Carlo method for the computer simulation of fluids , 1980 .

[46]  D. J. Adams,et al.  Grand canonical ensemble Monte Carlo for a Lennard-Jones fluid , 1975 .

[47]  I. Shrivastava,et al.  Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer. , 2000, Biophysical journal.

[48]  B. Roux,et al.  Statistical mechanical equilibrium theory of selective ion channels. , 1999, Biophysical journal.

[49]  W. Im,et al.  A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. , 2000, Biophysical journal.

[50]  M Mezei,et al.  Grand canonical ensemble Monte Carlo simulation of the dCpG/proflavine crystal hydrate. , 1996, Biophysical journal.

[51]  M. Sansom,et al.  Water in channel-like cavities: structure and dynamics. , 1996, Biophysical journal.

[52]  G. Phillips,et al.  Structure and dynamics of the water around myoglobin , 1995, Protein science : a publication of the Protein Society.