Implicit Solvation Parameters Derived from Explicit Water Forces in Large-Scale Molecular Dynamics Simulations

Implicit solvation is a mean force approach to model solvent forces acting on a solute molecule. It is frequently used in molecular simulations to reduce the computational cost of solvent treatment. In the first instance, the free energy of solvation and the associated solvent–solute forces can be approximated by a function of the solvent-accessible surface area (SASA) of the solute and differentiated by an atom–specific solvation parameter σiSASA. A procedure for the determination of values for the σiSASA parameters through matching of explicit and implicit solvation forces is proposed. Using the results of Molecular Dynamics simulations of 188 topologically diverse protein structures in water and in implicit solvent, values for the σiSASA parameters for atom types i of the standard amino acids in the GROMOS force field have been determined. A simplified representation based on groups of atom types σgSASA was obtained via partitioning of the atom–type σiSASA distributions by dynamic programming. Three groups of atom types with well separated parameter ranges were obtained, and their performance in implicit versus explicit simulations was assessed. The solvent forces are available at http://mathbio.nimr.mrc.ac.uk/wiki/Solvent_Forces.

[1]  Thomas Simonson,et al.  Testing the Coulomb/Accessible Surface Area solvent model for protein stability, ligand binding, and protein design , 2008, BMC Bioinformatics.

[2]  J. Viñals,et al.  Collapse transition of a hydrophobic self-avoiding random walk in a coarse-grained model solvent. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  Charles L Brooks,et al.  Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions. , 2008, Physical chemistry chemical physics : PCCP.

[4]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[5]  L. C. Martin,et al.  Using information theory to search for co-evolving residues in proteins , 2005, Bioinform..

[6]  C. Schiffer,et al.  The role of protein-solvent interactions in protein unfolding. , 1996, Current opinion in biotechnology.

[7]  Alessandro Pandini,et al.  MinSet: a general approach to derive maximally representative database subsets by using fragment dictionaries and its application to the SCOP database , 2007, Bioinform..

[8]  M. Malim,et al.  Rationalisation of the Differences between APOBEC3G Structures from Crystallography and NMR Studies by Molecular Dynamics Simulations , 2010, PloS one.

[9]  H. Scheraga,et al.  Identifying native‐like protein structures with scoring functions based on all‐atom ECEPP force fields, implicit solvent models and structure relaxation , 2009, Proteins.

[10]  F. Fraternali,et al.  Toward the understanding of MNEI sweetness from hydration map surfaces. , 2006, Biophysical journal.

[11]  M. Karplus,et al.  Effective energy functions for protein structure prediction. , 2000, Current opinion in structural biology.

[12]  Nathan A. Baker,et al.  Assessing implicit models for nonpolar mean solvation forces: the importance of dispersion and volume terms. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

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

[14]  Andreas Vitalis,et al.  ABSINTH: A new continuum solvation model for simulations of polypeptides in aqueous solutions , 2009, J. Comput. Chem..

[15]  Ronald M. Levy,et al.  AGBNP: An analytic implicit solvent model suitable for molecular dynamics simulations and high‐resolution modeling , 2004, J. Comput. Chem..

[16]  Alessandro Pandini,et al.  Structural alphabets derived from attractors in conformational space , 2010, BMC Bioinformatics.

[17]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[18]  Andrew J. Martin,et al.  The ups and downs of protein topology; rapid comparison of protein structure. , 2000, Protein engineering.

[19]  W. V. van Gunsteren,et al.  An efficient mean solvation force model for use in molecular dynamics simulations of proteins in aqueous solution. , 1996, Journal of molecular biology.

[20]  Jakob P Ulmschneider,et al.  A generalized born implicit-membrane representation compared to experimental insertion free energies. , 2007, Biophysical journal.

[21]  W. C. Still,et al.  A rapid approximation to the solvent accessible surface areas of atoms , 1988 .

[22]  F. Fraternali,et al.  Leap‐dynamics: efficient sampling of conformational space of proteins and peptides in solution , 2000, FEBS letters.

[23]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[24]  Arthur M Lesk,et al.  Contact patterns between helices and strands of sheet define protein folding patterns , 2007, Proteins.

[25]  Wilfred F van Gunsteren,et al.  A refined, efficient mean solvation force model that includes the interior volume contribution. , 2011, The journal of physical chemistry. B.

[26]  Peijuan Zhu,et al.  Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth‐permittivity finite difference Poisson–Boltzmann method , 2004, J. Comput. Chem..

[27]  William N. Venables,et al.  Modern Applied Statistics with S , 2010 .

[28]  Wang Lu,et al.  ON THE APPROXIMATION OF SOLVENT EFFECTS ON THE CONFORMATION AND DYNAMICS OF CYCLOSPORIN A BY STOCHASTIC DYNAMICS SIMULATION TECHNIQUES , 1988 .

[29]  Seishi Shimizu,et al.  Anti‐cooperativity and cooperativity in hydrophobic interactions: Three‐body free energy landscapes and comparison with implicit‐solvent potential functions for proteins , 2002, Proteins.

[30]  A. Pandini,et al.  Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics , 2012, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[31]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[32]  X. Daura,et al.  Parametrization of aliphatic CHn united atoms of GROMOS96 force field , 1998 .

[33]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

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

[35]  J. Apostolakis,et al.  Evaluation of a fast implicit solvent model for molecular dynamics simulations , 2002, Proteins.

[36]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

[37]  Chandra S Verma,et al.  Prion and water: tight and dynamical hydration sites have a key role in structural stability. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[38]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[39]  Patrice Koehl,et al.  The ASTRAL compendium for protein structure and sequence analysis , 2000, Nucleic Acids Res..

[40]  Themis Lazaridis,et al.  Implicit solvent simulations of peptide interactions with anionic lipid membranes , 2004, Proteins.

[41]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[42]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[43]  Franca Fraternali,et al.  Parameter optimized surfaces (POPS): analysis of key interactions and conformational changes in the ribosome , 2002, Nucleic Acids Res..

[44]  Franca Fraternali,et al.  Thermal unfolding simulations of apo‐calmodulin using leap‐dynamics , 2003, Proteins.

[45]  Emilio Gallicchio,et al.  On the nonpolar hydration free energy of proteins: surface area and continuum solvent models for the solute-solvent interaction energy. , 2003, Journal of the American Chemical Society.