Testing of the GROMOS Force-Field Parameter Set 54A8: Structural Properties of Electrolyte Solutions, Lipid Bilayers, and Proteins

The GROMOS 54A8 force field [Reif et al. J. Chem. Theory Comput.2012, 8, 3705–3723] is the first of its kind to contain nonbonded parameters for charged amino acid side chains that are derived in a rigorously thermodynamic fashion, namely a calibration against single-ion hydration free energies. Considering charged moieties in solution, the most decisive signature of the GROMOS 54A8 force field in comparison to its predecessor 54A7 can probably be found in the thermodynamic equilibrium between salt-bridged ion pair formation and hydration. Possible shifts in this equilibrium might crucially affect the properties of electrolyte solutions or/and the stability of (bio)molecules. It is therefore important to investigate the consequences of the altered description of charged oligoatomic species in the GROMOS 54A8 force field. The present study focuses on examining the ability of the GROMOS 54A8 force field to accurately model the structural properties of electrolyte solutions, lipid bilayers, and proteins. It is found that (i) aqueous electrolytes involving oligoatomic species (sodium acetate, methylammonium chloride, guanidinium chloride) reproduce experimental salt activity derivatives for concentrations up to 1.0 m (1.0-molal) very well, and good agreement between simulated and experimental data is also reached for sodium acetate and methylammonium chloride at 2.0 m concentration, while not even qualitative agreement is found for sodium chloride throughout the whole range of examined concentrations, indicating a failure of the GROMOS 54A7 and 54A8 force-field parameter sets to correctly account for the balance between ion–ion and ion–water binding propensities of sodium and chloride ions; (ii) the GROMOS 54A8 force field reproduces the liquid crystalline-like phase of a hydrated DPPC bilayer at a pressure of 1 bar and a temperature of 323 K, the area per lipid being in agreement with experimental data, whereas other structural properties (volume per lipid, bilayer thickness) appear underestimated; (iii) the secondary structure of a range of different proteins simulated with the GROMOS 54A8 force field at pH 7 is maintained and compatible with experimental NMR data, while, as also observed for the GROMOS 54A7 force field, α-helices are slightly overstabilized with respect to 310-helices; (iv) with the GROMOS 54A8 force field, the side chains of arginine, lysine, aspartate, and glutamate residues appear slightly more hydrated and present a slight excess of oppositely-charged solution components in their vicinity, whereas salt-bridge formation properties between charged residues at the protein surface, as assessed by probability distributions of interionic distances, are largely equivalent in the GROMOS 54A7 and 54A8 force-field parameter sets.

[1]  M. Shimizu.  [Electrolyte solutions]. , 2019, [Kango] Japanese journal of nursing.

[2]  Alexander D. MacKerell,et al.  CHARMM Additive All-Atom Force Field for Acyclic Polyalcohols, Acyclic Carbohydrates and Inositol. , 2009, Journal of chemical theory and computation.

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

[4]  B. Zagrovic,et al.  Mechanism and thermodynamics of binding of the polypyrimidine tract binding protein to RNA. , 2007, Biochemistry.

[5]  John C. Slater,et al.  The Van Der Waals Forces in Gases , 1931 .

[6]  Peter J. Artymiuk,et al.  The Structures of the Monoclinic and Orthorhombic Forms of Hen Egg-White Lysozyme at 6 Angstroms Resolution , 1981 .

[7]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[8]  O. Bonner.  Osmotic and activity coefficients of methyl-substituted ammonium chlorides , 1981 .

[9]  J. Grindlay,et al.  The Electrostatic Potential in a Finite Ionic Crystal , 1972 .

[10]  P. Hünenberger,et al.  A fast-Fourier transform method to solve continuum-electrostatics problems with truncated electrostatic interactions: Algorithm and application to ionic solvation and ion–ion interaction , 2003 .

[11]  R. Diamond,et al.  Treatment of NOE constraints involving equivalent or nonstereoassigned protons in calculations of biomacromolecular structures , 1996, Journal of biomolecular NMR.

[12]  G. Mie.  Zur kinetischen Theorie der einatomigen Körper , 1903 .

[13]  Paul E. Smith,et al.  Molecular Association in Solution: A Kirkwood−Buff Analysis of Sodium Chloride, Ammonium Sulfate, Guanidinium Chloride, Urea, and 2,2,2-Trifluoroethanol in Water , 2002 .

[14]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[15]  Markus Christen,et al.  Architecture, implementation and parallelisation of the GROMOS software for biomolecular simulation , 2012, Comput. Phys. Commun..

[16]  Philippe H. Hünenberger,et al.  Lattice-sum methods for computing electrostatic interactions in molecular simulations , 1999 .

[17]  Wilfred F van Gunsteren,et al.  GROMOS++ Software for the Analysis of Biomolecular Simulation Trajectories. , 2011, Journal of chemical theory and computation.

[18]  Wilfred F. van Gunsteren,et al.  Empirical Classical Force Fields for Molecular Systems , 1999 .

[19]  Arieh Ben-Naim,et al.  Statistical Thermodynamics for Chemists and Biochemists , 1992, Springer US.

[20]  Philippe H. Hünenberger,et al.  Optimal charge-shaping functions for the particle–particle—particle–mesh (P3M) method for computing electrostatic interactions in molecular simulations , 2000 .

[21]  P. Kast,et al.  1.6 A crystal structure of the secreted chorismate mutase from Mycobacterium tuberculosis: novel fold topology revealed. , 2006, Journal of molecular biology.

[22]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[23]  Janet E. Jones.  On the determination of molecular fields. —II. From the equation of state of a gas , 1924 .

[24]  Wilfred F. van Gunsteren,et al.  A new force field for simulating phosphatidylcholine bilayers , 2010, J. Comput. Chem..

[25]  P. Hünenberger,et al.  Influence of cut-off truncation and artificial periodicity of electrostatic interactions in molecular simulations of solvated ions: A continuum electrostatics study , 2003 .

[26]  C. Dobson,et al.  A refined solution structure of hen lysozyme determined using residual dipolar coupling data , 2001, Protein science : a publication of the Protein Society.

[27]  K Wüthrich,et al.  Pseudo-structures for the 20 common amino acids for use in studies of protein conformations by measurements of intramolecular proton-proton distance constraints with nuclear magnetic resonance. , 1983, Journal of molecular biology.

[28]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 3. The van der Waals' potentials and crystal data for aliphatic and aromatic hydrocarbons , 1989 .

[29]  Jean-Marc Simon,et al.  Thermodynamics of a small system in a μT reservoir , 2011 .

[30]  Roberto D. Lins,et al.  A new GROMOS force field for hexopyranose‐based carbohydrates , 2005, J. Comput. Chem..

[31]  E. A. Guggenheim.  L. The specific thermodynamic properties of aqueous solutions of strong electrolytes , 1935 .

[32]  S. Engelsen,et al.  THE CONSISTENT FORCE FIELD. II: AN OPTIMIZED SET OF POTENTIAL ENERGY FUNCTIONS FOR THE ALKANES , 1994 .

[33]  Wilfred F. van Gunsteren,et al.  Lattice‐sum methods for calculating electrostatic interactions in molecular simulations , 1995 .

[34]  Jeremy Pencer,et al.  Lipid bilayer structure determined by the simultaneous analysis of neutron and X-ray scattering data. , 2008, Biophysical journal.

[35]  B. Hess,et al.  Cation specific binding with protein surface charges , 2009, Proceedings of the National Academy of Sciences.

[36]  K. D. Collins,et al.  Ions in water: characterizing the forces that control chemical processes and biological structure. , 2007, Biophysical chemistry.

[37]  Arieh Ben-Naim,et al.  Molecular Theory of Solutions , 2006 .

[38]  H. Höltje,et al.  Consistent valence force‐field parameterization of bond lengths and angles with quantum chemical ab initio methods applied to some heterocyclic dopamine D3‐receptor agonists , 1998 .

[39]  R. Nussinov,et al.  Salt bridge stability in monomeric proteins. , 1999, Journal of molecular biology.

[40]  Gerhard Hummer,et al.  Simulation and Theory of Electrostatic Interactions in Solution: Computational Chemistry, Biophysics and Aqueous Solutions, Santa Fe, New Mexico, U. S. A., 23-25 June 1999 , 1999 .

[41]  P. Cremer,et al.  Interactions between macromolecules and ions: The Hofmeister series. , 2006, Current opinion in chemical biology.

[42]  R. L. Baldwin,et al.  How Hofmeister ion interactions affect protein stability. , 1996, Biophysical journal.

[43]  J. Missimer,et al.  Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1 , 2010, Journal of biomolecular NMR.

[44]  T. Halgren,et al.  Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rules , 1996 .

[45]  Maria M. Reif,et al.  Single-ion solvation : experimental and theoretical approaches to elusive thermodynamic quantities , 2011 .

[46]  Pramod C. Nair,et al.  An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. , 2011, Journal of chemical theory and computation.

[47]  Karl Nicholas Kirschner,et al.  GLYCAM06: A generalizable biomolecular force field. Carbohydrates , 2008, J. Comput. Chem..

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

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

[50]  R. Suter,et al.  X-ray structure determination of fully hydrated L alpha phase dipalmitoylphosphatidylcholine bilayers. , 1996, Biophysical journal.

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

[52]  J. Nagle,et al.  Structure of lipid bilayers. , 2000, Biochimica et biophysica acta.

[53]  André Bardow,et al.  Calculating thermodynamic properties from fluctuations at small scales. , 2011, The journal of physical chemistry. B.

[54]  H. Berendsen,et al.  Simulations of Proteins in Water a , 1986, Annals of the New York Academy of Sciences.

[55]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[56]  G. Guella,et al.  Puckering free energy of pyranoses: A NMR and metadynamics-umbrella sampling investigation. , 2010, The Journal of chemical physics.

[57]  Ming-Jing Hwang,et al.  Derivation of Class II Force Fields. 2. Derivation and Characterization of a Class II Force Field, CFF93, for the Alkyl Functional Group and Alkane Molecules , 1994 .

[58]  R. Bates,et al.  Osmotic coefficients and activity coefficients of guanidinium chloride in concentrated aqueous solutions at 25.degree.C , 1977 .

[59]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[60]  M. Billeter,et al.  Calibration of the angular dependence of the amide proton-Cα proton coupling constants, 3JHNα, in a globular protein: Use of 3JHNα for identification of helical secondary structure , 1984 .

[61]  W. V. van Gunsteren,et al.  New Interaction Parameters for Oxygen Compounds in the GROMOS Force Field: Improved Pure-Liquid and Solvation Properties for Alcohols, Ethers, Aldehydes, Ketones, Carboxylic Acids, and Esters. , 2011, Journal of chemical theory and computation.

[62]  Arnold T. Hagler,et al.  On the use of quantum energy surfaces in the derivation of molecular force fields , 1994 .

[63]  Chris Oostenbrink,et al.  An improved nucleic acid parameter set for the GROMOS force field , 2005, J. Comput. Chem..

[64]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

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

[66]  Wilfred F van Gunsteren,et al.  Biomolecular modeling: Goals, problems, perspectives. , 2006, Angewandte Chemie.

[67]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[68]  U Heinemann,et al.  Crystal structure of CspA, the major cold shock protein of Escherichia coli. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[69]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[70]  Xavier Daura,et al.  Molecular dynamics simulations of a reversibly folding beta-heptapeptide in methanol: influence of the treatment of long-range electrostatic interactions. , 2009, The journal of physical chemistry. B.

[71]  Werner Kunz,et al.  Hofmeister series and specific interactions of charged headgroups with aqueous ions. , 2009, Advances in colloid and interface science.

[72]  X. Daura,et al.  Derivation of an improved simple point charge model for liquid water: SPC/A and SPC/L , 2002 .

[73]  Maria M. Reif,et al.  New Interaction Parameters for Charged Amino Acid Side Chains in the GROMOS Force Field. , 2012, Journal of chemical theory and computation.

[74]  Andreas P. Eichenberger,et al.  Definition and testing of the GROMOS force-field versions 54A7 and 54B7 , 2011, European Biophysics Journal.

[75]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

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

[77]  J A McCammon,et al.  Molecular dynamics simulations of the hyperthermophilic protein sac7d from Sulfolobus acidocaldarius: contribution of salt bridges to thermostability. , 1999, Journal of molecular biology.

[78]  S. Edmondson,et al.  Solution structure of the DNA-binding protein Sac7d from the hyperthermophile Sulfolobus acidocaldarius. , 1995, Biochemistry.

[79]  H. Berendsen.  Simulating the Physical World , 2004 .

[80]  Y. Yamagata,et al.  Contribution of salt bridges near the surface of a protein to the conformational stability. , 2000, Biochemistry.

[81]  Wolfgang Jahnke,et al.  Molecular basis of coiled-coil formation , 2007, Proceedings of the National Academy of Sciences.

[82]  H. Friedman,et al.  Theory of mixed electrolyte solutions and application to a model for aqueous lithium chloride-cesium chloride , 1970 .

[83]  J. A. Barker,et al.  Monte Carlo studies of the dielectric properties of water-like models , 1973 .

[84]  J. Nagle,et al.  Closer look at structure of fully hydrated fluid phase DPPC bilayers. , 2006, Biophysical journal.

[85]  S. Weerasinghe,et al.  A Kirkwood–Buff derived force field for sodium chloride in water , 2003 .

[86]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[87]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[88]  Tim N. Heinz,et al.  Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations , 2001 .

[89]  Berk Hess,et al.  Osmotic coefficients of atomistic NaCl (aq) force fields. , 2006, The Journal of chemical physics.

[90]  D. Black,et al.  Molecular basis of RNA recognition by the human alternative splicing factor Fox‐1 , 2006, The EMBO journal.

[91]  Philippe H. Hünenberger,et al.  A reoptimized GROMOS force field for hexopyranose‐based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers , 2011, J. Comput. Chem..

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

[93]  Janet E. Jones.  On the Determination of Molecular Fields. I. From the Variation of the Viscosity of a Gas with Temperature , 1924 .

[94]  Siewert J. Marrink,et al.  Methodological issues in lipid bilayer simulations , 2003 .

[95]  M. Karplus.  Contact Electron‐Spin Coupling of Nuclear Magnetic Moments , 1959 .

[96]  Victor H Rusu,et al.  GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates. , 2012, Journal of chemical theory and computation.

[97]  S L Mayo,et al.  Contribution of surface salt bridges to protein stability. , 2000, Biochemistry.

[98]  Lennart Nilsson,et al.  Advances in biomolecular simulations: methodology and recent applications , 2003, Quarterly Reviews of Biophysics.

[99]  L. J. Lis,et al.  Interactions between neutral phospholipid bilayer membranes. , 1982, Biophysical journal.

[100]  K. Dill.  Dominant forces in protein folding. , 1990, Biochemistry.

[101]  O. Redlich,et al.  On the Thermodynamics of Solutions. IV. The Determination of Liquid--Vapor Equilibria by Measuring the Total Pressure , 1949 .

[102]  E. Jakobsson,et al.  Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane. , 1995, Biophysical journal.

[103]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Solutions. I , 1951 .