Combining a polarizable force‐field and a coarse‐grained polarizable solvent model. II. Accounting for hydrophobic effects

A revised and improved version of our efficient polarizable force‐field/coarse grained solvent combined approach (Masella, Borgis, and Cuniasse, J. Comput. Chem. 2008, 29, 1707) is described. The polarizable pseudo‐particle solvent model represents the macroscopic solvent polarization by induced dipoles placed on mobile pseudo‐particles. In this study, we propose a new formulation of the energy term handling the nonelectrostatic interactions among the pseudo‐particles. This term is now able to reproduce the energetic and structural response of liquid water due to the presence of a hydrophobic spherical cavity. Accordingly, the parameters of the energy term handling the nonpolar solute/solvent interactions have been refined to reproduce the free‐solvation energy of small solutes, based on a standard thermodynamic integration scheme. The reliability of this new approach has been checked for the properties of solvated methane and of the solvated methane dimer, as well as by performing 10 × 20 ns molecular dynamics (MD) trajectories for three solvated proteins. A long‐time stability of the protein structures along the trajectories is observed. Moreover, our method still provides a measure of the protein solvation thermodynamic at the same accuracy as standard Poisson–Boltzman continuum methods. These results show the relevance of our approach and its applicability to massively coupled MD schemes to accurately and intensively explore solvated macromolecule potential energy surfaces. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011

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

[2]  Daniel Borgis,et al.  Electrostatics on particles: Phenomenological and orientational density functional theory approach , 2002 .

[3]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[4]  Piotr Cieplak,et al.  Polarization effects in molecular mechanical force fields , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[5]  A. Warshel,et al.  Electrostatic contributions to protein–protein binding affinities: Application to Rap/Raf interaction , 1998, Proteins.

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

[7]  Arieh Warshel,et al.  Langevin Dipoles Model for ab Initio Calculations of Chemical Processes in Solution: Parametrization and Application to Hydration Free Energies of Neutral and Ionic Solutes and Conformational Analysis in Aqueous Solution , 1997 .

[8]  Amish J. Patel,et al.  An improved coarse-grained model of solvation and the hydrophobic effect. , 2010, The Journal of chemical physics.

[9]  A. Bonvin,et al.  NMR Relaxation and Internal Dynamics of Ubiquitin from a 0.2 μs MD Simulation. , 2005, Journal of chemical theory and computation.

[10]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[11]  T. Darden,et al.  Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions , 2000 .

[12]  Richard A. Friesner,et al.  Accurate ab Initio Quantum Chemical Determination of the Relative Energetics of Peptide Conformations and Assessment of Empirical Force Fields , 1997 .

[13]  I. Pagonabarraga,et al.  Dissipative particle dynamics for interacting systems , 2001, cond-mat/0105075.

[14]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Sergio A Hassan,et al.  Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: Calculation of structural and dynamic properties , 2005, Proteins.

[16]  R. Jernigan,et al.  Solvent effect on binding thermodynamics of biopolymers , 1990, Biopolymers.

[17]  Daniel Borgis,et al.  Combining a polarizable force‐field and a coarse‐grained polarizable solvent model: Application to long dynamics simulations of bovine pancreatic trypsin inhibitor , 2008, J. Comput. Chem..

[18]  Bernard R. Brooks,et al.  Solvent-Induced Forces between Two Hydrophilic Groups , 1994 .

[19]  D. Chandler,et al.  Scaling of Hydrophobic Solvation Free Energies , 2001 .

[20]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

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

[22]  S. F. Lienin,et al.  Anisotropic Intramolecular Backbone Dynamics of Ubiquitin Characterized by NMR Relaxation and MD Computer Simulation , 1998 .

[23]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[24]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[25]  Valeria Molinero,et al.  Water modeled as an intermediate element between carbon and silicon. , 2009, The journal of physical chemistry. B.

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

[27]  J. Pitera,et al.  Dielectric properties of proteins from simulation: the effects of solvent, ligands, pH, and temperature. , 2001, Biophysical journal.

[28]  L. Pratt,et al.  Role of attractive methane-water interactions in the potential of mean force between methane molecules in water. , 2008, The Journal of chemical physics.

[29]  Wilfred F. van Gunsteren,et al.  SWARM-MD: Searching Conformational Space by Cooperative Molecular Dynamics , 1998 .

[30]  B. Rupp,et al.  Structure of bovine pancreatic trypsin inhibitor at 125 K definition of carboxyl-terminal residues Gly57 and Ala58. , 1996, Acta crystallographica. Section D, Biological crystallography.

[31]  Y. Jan,et al.  Probing Protein Electrostatics with a Synthetic Fluorescent Amino Acid , 2002, Science.

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

[33]  D. Chandler,et al.  Hydrophobicity at Small and Large Length Scales , 1999 .

[34]  Wei Wang,et al.  Fast evaluation of polarizable forces. , 2005, The Journal of chemical physics.

[35]  Dongping Zhong,et al.  Protein hydration dynamics and molecular mechanism of coupled water-protein fluctuations. , 2009, Journal of the American Chemical Society.

[36]  N. Wingreen,et al.  Hydrophobic interaction and hydrogen-bond network for a methane pair in liquid water , 2007, Proceedings of the National Academy of Sciences.

[37]  G L Gilliland,et al.  Two crystal structures of the B1 immunoglobulin-binding domain of streptococcal protein G and comparison with NMR. , 1994, Biochemistry.

[38]  A. Ben-Naim.  Solvent effects on protein association and protein folding , 1990, Biopolymers.

[39]  Carolyn A. Koh,et al.  Water ordering around methane during hydrate formation , 2000 .

[40]  Daniel Borgis,et al.  A semi‐implicit solvent model for the simulation of peptides and proteins , 2004, J. Comput. Chem..

[41]  Johann Fischer,et al.  Molecular dynamics simulation of the liquid–vapor interface: The Lennard-Jones fluid , 1997 .

[42]  M. Yáñez,et al.  Cooperative (nonpairwise) effects in water trimers: An ab initio molecular orbital study , 1992 .

[43]  M. Nakasako,et al.  Probability distributions of hydration water molecules around polar protein atoms obtained by a database analysis. , 2009, The journal of physical chemistry. B.

[44]  P. Cuniasse,et al.  A many-body model to study proteins. I. Applications to MLnm+ complexes, Mm+=Li+, Na+, K+, Mg2+, Ca2+, and Zn2+, L=H2O, CH3OH, HCONH2, n=1–6, and to small hydrogen bonded systems , 2003 .

[45]  S. Hassan,et al.  A critical analysis of continuum electrostatics: The screened Coulomb potential–implicit solvent model and the study of the alanine dipeptide and discrimination of misfolded structures of proteins , 2002, Proteins.

[46]  A. Ben-Naim.  Solvent-induced forces in protein folding , 1990 .

[47]  A. Ben-Naim.  Strong forces between hydrophilic macromolecules: Implications in biological systems , 1990 .

[48]  R. Winter,et al.  Molecular Dynamics Simulations of Staphylococcal Nuclease: Properties of Water at the Protein Surface , 2004 .

[49]  B. Schimmelpfennig,et al.  Quantum chemical and molecular dynamics study of the coordination of Th(IV) in aqueous solvent. , 2010, The journal of physical chemistry. B.

[50]  Sergio A Hassan,et al.  Liquid-structure forces and electrostatic modulation of biomolecular interactions in solution. , 2007, The journal of physical chemistry. B.

[51]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[52]  Claudio N. Cavasotto,et al.  Quantum mechanical dynamics of charge transfer in ubiquitin in aqueous solution. , 2009, Chemphyschem : a European journal of chemical physics and physical chemistry.

[53]  J. Feigon,et al.  Structural determinants for the binding of ubiquitin‐like domains to the proteasome , 2003, The EMBO journal.

[54]  Peter G. Wolynes,et al.  Classical solvent dynamics and electron transfer. 1. Continuum theory , 1983 .

[55]  M Karplus,et al.  Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme. , 1989, Journal of molecular biology.

[56]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[57]  J. Flament,et al.  A pairwise and two many-body models for water: Influence of nonpairwise effects upon the stability and geometry of (H2O)n cyclic (n=3–6) and cagelike (n=6–20) clusters , 1997 .

[58]  Daniel Borgis,et al.  Particle-Based Implicit Solvent Model for Biosimulations:  Application to Proteins and Nucleic Acids Hydration. , 2006, Journal of chemical theory and computation.

[59]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[60]  Teresa Head-Gordon,et al.  Hydrophobic potential of mean force as a solvation function for protein structure prediction. , 2007, Structure.

[61]  Sergio A. Hassan,et al.  Intermolecular Potentials of Mean Force of Amino Acid Side Chain Interactions in Aqueous Medium , 2004 .

[62]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[63]  S. Hassan,et al.  A General Treatment of Solvent Effects Based on Screened Coulomb Potentials , 2000 .

[64]  Sotiris S. Xantheas,et al.  Ab initio studies of cyclic water clusters (H2O)n, n=1–6. II. Analysis of many‐body interactions , 1994 .

[65]  Richard A Friesner,et al.  Structure and dynamics of the solvation of bovine pancreatic trypsin inhibitor in explicit water: a comparative study of the effects of solvent and protein polarizability. , 2005, The journal of physical chemistry. B.

[66]  D. Chandler.  Interfaces and the driving force of hydrophobic assembly , 2005, Nature.

[67]  Kazuyuki Akasaka,et al.  Simulation of Hydrated BPTI at High Pressure: Changes in Hydrogen Bonding and Its Relation with NMR Experiments , 2001 .

[68]  Alan K. Soper,et al.  Site–site pair correlation functions of water from 25 to 400 °C: Revised analysis of new and old diffraction data , 1997 .

[69]  Michel Masella,et al.  The multiple time step r-RESPA procedure and polarizable potentials based on induced dipole moments , 2006 .

[70]  J Andrew McCammon,et al.  Computing accurate potentials of mean force in electrolyte solutions with the generalized gradient-augmented harmonic Fourier beads method. , 2008, The Journal of chemical physics.

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