Application of inhomogeneous fluid solvation theory to model the distribution and thermodynamics of water molecules around biomolecules.

The structures of biomolecules and the strengths of association between them depend critically on interactions with water molecules. Thus, understanding these interactions is a prerequisite for understanding the structure and function of all biomolecules. Inhomogeneous fluid solvation theory provides a framework to derive thermodynamic properties of individual water molecules from a statistical mechanical analysis. In this work, two biomolecules are analysed to probe the distribution and thermodynamics of surrounding water molecules. The great majority of hydration sites are predicted to contribute favourably to the total free energy with respect to bulk water, though hydration sites close to non-polar regions of the solute do not contribute significantly. Analysis of a biomolecule with a positively and negatively charged functional group predicts that a charged species perturbs the free energy of water molecules to a distance of approximately 6.0 Å. Interestingly, short simulations are found to provide converged predictions if samples are taken with sufficient frequency, a finding that has the potential to significantly reduce the required computational cost of such analysis. In addition, the predicted thermodynamic properties of hydration sites with the potential for direct hydrogen bonding interactions are found to disagree significantly for two different water models. This study provides important information on how inhomogeneous fluid solvation theory can be employed to understand the structures and intermolecular interactions of biomolecules.

[1]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[2]  C. Chakravarty,et al.  Comparison of Tetrahedral Order, Liquid State Anomalies, and Hydration Behavior of mTIP3P and TIP4P Water Models. , 2011, Journal of chemical theory and computation.

[3]  Matsuda,et al.  Physical nature of higher-order mutual information: intrinsic correlations and frustration , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  J. Andrew McCammon,et al.  How Can Hydrophobic Association Be Enthalpy Driven? , 2010, Journal of chemical theory and computation.

[5]  Mike C. Payne,et al.  Thermodynamic Properties of Water Molecules at a Protein–Protein Interaction Surface , 2011, Journal of chemical theory and computation.

[6]  J. Zielkiewicz.  Structural properties of water: comparison of the SPC, SPCE, TIP4P, and TIP5P models of water. , 2005, The Journal of chemical physics.

[7]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

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

[9]  R. Best,et al.  Free‐energy landscape of the GB1 hairpin in all‐atom explicit solvent simulations with different force fields: Similarities and differences , 2011, Proteins.

[10]  Jeremy C. Smith,et al.  Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study. , 1996, Biophysical journal.

[11]  Themis Lazaridis,et al.  Solvent Reorganization Energy and Entropy in Hydrophobic Hydration , 2000 .

[12]  Franz Saija,et al.  Entropy from Correlations in TIP4P Water. , 2010, Journal of chemical theory and computation.

[13]  Woody Sherman,et al.  High‐energy water sites determine peptide binding affinity and specificity of PDZ domains , 2009, Protein science : a publication of the Protein Society.

[14]  R. Levy,et al.  Enthalpy−Entropy and Cavity Decomposition of Alkane Hydration Free Energies: Numerical Results and Implications for Theories of Hydrophobic Solvation , 2000 .

[15]  Gerhard Hummer,et al.  Molecular binding: Under water's influence. , 2010, Nature chemistry.

[16]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

[17]  C. Vega,et al.  A general purpose model for the condensed phases of water: TIP4P/2005. , 2005, The Journal of chemical physics.

[18]  Alan K. Soper,et al.  The radial distribution functions of water and ice from 220 to 673 K and at pressures up to 400 MPa , 2000 .

[19]  T. Lazaridis.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 1. Theory , 1998 .

[20]  Robert Abel,et al.  Motifs for molecular recognition exploiting hydrophobic enclosure in protein–ligand binding , 2007, Proceedings of the National Academy of Sciences.

[21]  S. Rick.  A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums. , 2004, The Journal of chemical physics.

[22]  Themis Lazaridis,et al.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 2. Applications to Simple Fluids , 1998 .

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

[24]  Robert J Woods,et al.  On the Role of Water Models in Quantifying the Binding Free Energy of Highly Conserved Water Molecules in Proteins: The Case of Concanavalin A. , 2011, Journal of chemical theory and computation.

[25]  Themis Lazaridis,et al.  Thermodynamics of buried water clusters at a protein-ligand binding interface. , 2006, The journal of physical chemistry. B.

[26]  W. Sherman,et al.  Understanding Kinase Selectivity Through Energetic Analysis of Binding Site Waters , 2010, ChemMedChem.

[27]  R. Best,et al.  Protein simulations with an optimized water model: cooperative helix formation and temperature-induced unfolded state collapse. , 2010, The journal of physical chemistry. B.

[28]  Alexander D. MacKerell,et al.  A polarizable model of water for molecular dynamics simulations of biomolecules , 2006 .

[29]  D. Hamelberg,et al.  A dry ligand-binding cavity in a solvated protein , 2008, Proceedings of the National Academy of Sciences.

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

[31]  W. Sherman,et al.  Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization , 2011, Proteins.

[32]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[33]  M. Otyepka,et al.  Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact. , 2010, Journal of chemical theory and computation.

[34]  Themis Lazaridis,et al.  Orientational correlations and entropy in liquid water , 1996 .

[35]  James A. Krumhansl,et al.  Triplet Correlation in Liquid Argon by Monte Carlo Method: Low Densities , 1972 .

[36]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[37]  B. Hess,et al.  Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models. , 2006, The journal of physical chemistry. B.

[38]  Themis Lazaridis,et al.  Thermodynamic contributions of the ordered water molecule in HIV-1 protease. , 2003, Journal of the American Chemical Society.

[39]  Fumio Hirata,et al.  Water molecules in a protein cavity detected by a statistical-mechanical theory. , 2005, Journal of the American Chemical Society.

[40]  David J. Huggins,et al.  Benchmarking the thermodynamic analysis of water molecules around a model beta sheet , 2012, J. Comput. Chem..

[41]  Jeremy C. Smith,et al.  Dual function of the hydration layer around an antifreeze protein revealed by atomistic molecular dynamics simulations. , 2008, Journal of the American Chemical Society.

[42]  S. Hovmöller,et al.  Conformations of amino acids in proteins. , 2002, Acta crystallographica. Section D, Biological crystallography.

[43]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[44]  A. Soper,et al.  Structural studies on the hydration of L-glutamic acid in solution. , 2006, The journal of physical chemistry. B.

[45]  R. Baron,et al.  Water in Cavity−Ligand Recognition , 2010, Journal of the American Chemical Society.

[46]  B. Berne,et al.  Dewetting transitions in protein cavities , 2010, Proteins.

[47]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[48]  Olexandr Isayev,et al.  Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk. , 2011, The journal of physical chemistry. A.

[49]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[50]  Michael K Gilson,et al.  Extraction of configurational entropy from molecular simulations via an expansion approximation. , 2007, The Journal of chemical physics.

[51]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[52]  David J Huggins,et al.  Correlations in liquid water for the TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP models. , 2012, The Journal of chemical physics.

[53]  Robert Abel,et al.  Thermodynamic properties of liquid water: an application of a nonparametric approach to computing the entropy of a neat fluid. , 2009, Journal of chemical theory and computation.

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

[55]  W. L. Jorgensen,et al.  Energetics of displacing water molecules from protein binding sites: consequences for ligand optimization. , 2009, Journal of the American Chemical Society.

[56]  Michael K Gilson,et al.  Grid inhomogeneous solvation theory: hydration structure and thermodynamics of the miniature receptor cucurbit[7]uril. , 2012, The Journal of chemical physics.

[57]  J. Zielkiewicz,et al.  Structural properties of hydration shell around various conformations of simple polypeptides. , 2010, The journal of physical chemistry. B.