Dimerization of Amino Acid Side Chains: Lessons from the Comparison of Different Force Fields.

The interactions between amino acid side chains govern protein secondary, tertiary, and quaternary structure formation. For molecular modeling approaches to be able to realistically describe these phenomena, the underlying force fields have to represent these interactions as accurately as possible. Here, we compare the side chain-side chain interactions for a number of commonly used force fields, namely the all-atom OPLS, the united-atom GROMOS, and the coarse-grain MARTINI force field. We do so by calculating the dimerization free energies between selected pairs of side chains and structural characterization of their binding modes. To mimic both polar and nonpolar environments, the simulations are performed in water, n-octanol, and decane. In general, reasonable correlations are found between all three force fields, with deviations on the order of 1 kT in aqueous solvent. In apolar solvent, however, significantly larger differences are found, especially for charged amino acid pairs between the OPLS and GROMOS force fields, and for polar interactions in the MARTINI force field in comparison to the higher resolution models. Interestingly, even in cases where the dimerization free energies are similar, the binding mode may differ substantially between the force fields. This was found to be especially the case for aromatic residues. In addition to the inter-force-field comparison, we compared the various force fields to a knowledge-based potential. The two independent approaches show good correlation in aqueous solvent with an exception of aromatic residues for which the interaction strength is lower in the knowledge-based potentials.

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

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

[3]  M. Betancourt,et al.  Another look at the conditions for the extraction of protein knowledge‐based potentials , 2009, Proteins.

[4]  H. Scheraga,et al.  Medium- and long-range interaction parameters between amino acids for predicting three-dimensional structures of proteins. , 1976, Macromolecules.

[5]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[6]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[7]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[8]  Miguel Jorge,et al.  1-Octanol/Water Partition Coefficients of n-Alkanes from Molecular Simulations of Absolute Solvation Free Energies. , 2009, Journal of chemical theory and computation.

[9]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

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

[11]  M. Sippl.  Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins. , 1990, Journal of molecular biology.

[12]  C. Dobson,et al.  Protein Dynamics: Moore's Law in Molecular Biology , 2011, Current Biology.

[13]  F. Diederich,et al.  Interactions with aromatic rings in chemical and biological recognition. , 2003, Angewandte Chemie.

[14]  Durba Sengupta,et al.  Polarizable Water Model for the Coarse-Grained MARTINI Force Field , 2010, PLoS Comput. Biol..

[15]  M. Karplus,et al.  Dynamics of folded proteins , 1977, Nature.

[16]  Marcos R Betancourt,et al.  Pairwise energies for polypeptide coarse-grained models derived from atomic force fields. , 2009, The Journal of chemical physics.

[17]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[18]  Cecilia Clementi,et al.  Coarse-grained models of protein folding: toy models or predictive tools? , 2008, Current opinion in structural biology.

[19]  Adam Liwo,et al.  Protein-folding dynamics: overview of molecular simulation techniques. , 2007, Annual review of physical chemistry.

[20]  Helmut Grubmüller,et al.  Determining equilibrium constants for dimerization reactions from molecular dynamics simulations , 2011, J. Comput. Chem..

[21]  S. Wold,et al.  New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. , 1998, Journal of medicinal chemistry.

[22]  R. Jernigan,et al.  Estimation of effective interresidue contact energies from protein crystal structures: quasi-chemical approximation , 1985 .

[23]  Justin L MacCallum,et al.  Structures of neat and hydrated 1-octanol from computer simulations. , 2002, Journal of the American Chemical Society.

[24]  Dusan Bratko,et al.  Protein aggregation in silico. , 2007, Trends in biotechnology.

[25]  N. Ben-Tal,et al.  Free Energy of Amide Hydrogen Bond Formation in Vacuum, in Water, and in Liquid Alkane Solution , 1997 .

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

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

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

[29]  Vincent A Voelz,et al.  Taming the complexity of protein folding. , 2011, Current opinion in structural biology.

[30]  R. Jernigan,et al.  Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. , 1996, Journal of molecular biology.

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

[32]  Klaus Schulten,et al.  Challenges in protein-folding simulations , 2010 .

[33]  Lewis Stiller,et al.  Computation of the mean residence time of water in the hydration shells of biomolecules , 1993, J. Comput. Chem..

[34]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

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

[36]  D. Tieleman,et al.  Using the Wimley-White Hydrophobicity Scale as a Direct Quantitative Test of Force Fields: The MARTINI Coarse-Grained Model. , 2011, Journal of chemical theory and computation.

[37]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[38]  Matthias Heinig,et al.  New Insights into the Genetic Control of Gene Expression using a Bayesian Multi-tissue Approach , 2010, PLoS Comput. Biol..

[39]  R. Jernigan,et al.  Self‐consistent estimation of inter‐residue protein contact energies based on an equilibrium mixture approximation of residues , 1999, Proteins.

[40]  K. Dill,et al.  The protein folding problem. , 1993, Annual review of biophysics.

[41]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[42]  Thomas Huber,et al.  G protein-coupled receptors self-assemble in dynamics simulations of model bilayers. , 2007, Journal of the American Chemical Society.