Simple physics-based analytical formulas for the potentials of mean force of the interaction of amino-acid side chains in water. VI. Oppositely charged side chains.

The two-site coarse-grained model for the interactions of charged side chains, to be used with our coarse-grained UNRES force field for protein simulations proposed in the accompanying paper, has been extended to pairs of oppositely charged side chains. The potentials of mean force of four pairs of molecules modeling charged amino-acid side chains, i.e., propionate-n-pentylamine cation (for aspartic acid-lysine), butyrate-n-pentylamine cation (for glutamic acid-lysine), propionate-1-butylguanidine (for aspartic acid-arginine), and butyrate-1-butylguanidine (for glutamic acid-arginine) pairs were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation, and the analytical expression was fitted to the potentials of mean force. Compared to pairs of like-charged side chains discussed in the accompanying paper, an average quadrupole-quadrupole interaction term had to be introduced to reproduce the Coulombic interactions, and a multistate model of charge distribution had to be introduced to fit the potentials of mean force of all oppositely charged pairs well. The model reproduces all salt-bridge minima and, consequently, is likely to improve the performance of the UNRES force field.

[1]  Adam Liwo,et al.  beta-hairpin-forming peptides; models of early stages of protein folding. , 2010, Biophysical chemistry.

[2]  A. Liwo,et al.  Parametrization of Backbone−Electrostatic and Multibody Contributions to the UNRES Force Field for Protein-Structure Prediction from Ab Initio Energy Surfaces of Model Systems† , 2004 .

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

[4]  A. Ben-Naim.  STATISTICAL POTENTIALS EXTRACTED FROM PROTEIN STRUCTURES : ARE THESE MEANINGFUL POTENTIALS? , 1997 .

[5]  A. Liwo,et al.  Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. 1. Approximate expression for the free energy of hydrophobic association based on a Gaussian-overlap model. , 2007, The journal of physical chemistry. B.

[6]  Harold A. Scheraga,et al.  A combined iterative and boundary-element approach for solution of the nonlinear Poisson-Boltzmann equation , 1992 .

[7]  Haipeng Gong,et al.  Influence of nonlinear electrostatics on transfer energies between liquid phases: Charge burial is far less expensive than Born model , 2008, Proceedings of the National Academy of Sciences.

[8]  A. Liwo,et al.  Cumulant-based expressions for the multibody terms for the correlation between local and electrostatic interactions in the united-residue force field , 2001 .

[9]  A. Liwo,et al.  Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. 2. Tests with simple spherical systems. , 2007, The journal of physical chemistry. B.

[10]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[11]  Adam Liwo,et al.  An improved functional form for the temperature scaling factors of the components of the mesoscopic UNRES force field for simulations of protein structure and dynamics. , 2009, The journal of physical chemistry. B.

[12]  Energetic and entropic contributions to the interactions between like‐charged groups in cationic peptides: A molecular dynamics simulation study , 2002, Protein science : a publication of the Protein Society.

[13]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

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

[15]  Christophe Chipot,et al.  Modeling ion-ion interaction in proteins: A molecular dynamics free energy calculation of the guanidinium-acetate association , 2000 .

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

[17]  A. Liwo,et al.  Modification and optimization of the united-residue (UNRES) potential energy function for canonical simulations. I. Temperature dependence of the effective energy function and tests of the optimization method with single training proteins. , 2007, The journal of physical chemistry. B.

[18]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. II. Solvation structure from the three-dimensional reference interaction site model approach, and comparison with simulations , 2000 .

[19]  Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. 3. Calculation and parameterization of the potentials of mean force of pairs of identical hydrophobic side chains. , 2007, The journal of physical chemistry. B.

[20]  H. Scheraga,et al.  Influence of water structure and of hydrophobic interactions on the strength of side‐chain hydrogen bonds in proteins , 1963 .

[21]  K. Dill,et al.  Statistical potentials extracted from protein structures: how accurate are they? , 1996, Journal of molecular biology.

[22]  S. Rackovsky,et al.  Calculation of protein backbone geometry from α‐carbon coordinates based on peptide‐group dipole alignment , 1993, Protein science : a publication of the Protein Society.

[23]  Adam Liwo,et al.  Determination of side‐chain‐rotamer and side‐chain and backbone virtual‐bond‐stretching potentials of mean force from AM1 energy surfaces of terminally‐blocked amino‐acid residues, for coarse‐grained simulations of protein structure and folding. I. The method , 2010, J. Comput. Chem..

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

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

[26]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[27]  P. Kollman,et al.  The effect of water models on the interaction of the sodium–chloride ion pair in water: Molecular dynamics simulations , 1990 .

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

[29]  A. Liwo,et al.  A united‐residue force field for off‐lattice protein‐structure simulations. I. Functional forms and parameters of long‐range side‐chain interaction potentials from protein crystal data , 1997 .

[30]  W. L. Jorgensen,et al.  Energetics and hydration of the constituent ion pairs of tetramethylammonium chloride , 1989 .

[31]  M. B. Pinto,et al.  Optimized δ expansion for relativistic nuclear models , 1997, nucl-th/9709049.

[32]  Adam Liwo,et al.  Exploring the parameter space of the coarse‐grained UNRES force field by random search: Selecting a transferable medium‐resolution force field , 2009, J. Comput. Chem..

[33]  Adam Liwo,et al.  Theoretical calculations of homoconjugation equilibrium constants in systems modeling acid–base interactions in side chains of biomolecules using the potential of mean force , 2005, J. Comput. Chem..

[34]  Jay W Ponder,et al.  Polarizable Atomic Multipole Solutes in a Generalized Kirkwood Continuum. , 2007, Journal of chemical theory and computation.

[35]  A. Liwo,et al.  Optimization of the UNRES Force Field by Hierarchical Design of the Potential-Energy Landscape. 3. Use of Many Proteins in Optimization , 2004 .

[36]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[37]  Andrzej Kolinski,et al.  Multiscale Approaches to Protein Modeling , 2011 .

[38]  K. Higashitani,et al.  A Reexamination of Mean Force Potentials for the Methane Pair and the Constituent Ion Pairs of NaCl in Water , 2003 .

[39]  C. Chipot,et al.  Effect of Polarizability on the Potential of Mean Force of Two Cations. The Guanidinium−Guanidinium Ion Pair in Water , 1997 .

[40]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[41]  Adam Liwo,et al.  Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. IV. Pairs of different hydrophobic side chains. , 2008, The journal of physical chemistry. B.

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

[43]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

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

[45]  H. Scheraga,et al.  Contribution of unusual Arginine-Arginine short-range interactions to stabilization and recognition in proteins , 1994, Journal of protein chemistry.

[46]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[47]  Adam Liwo,et al.  Molecular Simulation Study of the Potentials of Mean Force for the Interactions between Models of Like-Charged and between Charged and Nonpolar Amino Acid Side Chains in Water , 2003 .

[48]  B. Pettitt,et al.  A theoretical study of like ion pairs in solution , 1990 .

[49]  B. Berne.  Modification of the overlap potential to mimic a linear site-site potential , 1981 .

[50]  William L. Jorgensen,et al.  Monte Carlo simulations of liquid acetonitrile with a three-site model , 1988 .

[51]  T. Lazaridis,et al.  Potentials of mean force between ionizable amino acid side chains in water. , 2003, Journal of the American Chemical Society.

[52]  O. Karim.  Potential of mean force for an aqueous chloride ion pair: Simulation with a polarizable model , 1992 .

[53]  A. Liwo,et al.  United‐residue force field for off‐lattice protein‐structure simulations: III. Origin of backbone hydrogen‐bonding cooperativity in united‐residue potentials , 1998 .

[54]  Benoît Roux,et al.  Solvation thermodynamics: An approach from analytic temperature derivatives , 1990 .

[55]  On the correlation between like ion pairs in water , 1992 .

[56]  Bogdan Lesyng,et al.  Generalized Born model: Analysis, refinement, and applications to proteins , 2004 .

[57]  Elvira Guàrdia,et al.  Na+–Na+ and Cl−–Cl− ion pairs in water: Mean force potentials by constrained molecular dynamics , 1991 .

[58]  Mihaly Mezei,et al.  The Potentials of Mean Force of Sodium Chloride and Sodium Dimethylphosphate in Water : An Application of Adaptive Umbrella Sampling , 1995 .

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

[60]  Adam Liwo,et al.  Determination of virtual-bond-angle potentials of mean force for coarse-grained simulations of protein structure and folding from ab initio energy surfaces of terminally-blocked glycine, alanine, and proline , 2007 .

[61]  B. Montgomery Pettitt,et al.  Alkali halides in water: Ion–solvent correlations and ion–ion potentials of mean force at infinite dilution , 1986 .

[62]  Adam Liwo,et al.  Determination of the Potentials of Mean Force for Rotation about Cα−Cα Virtual Bonds in Polypeptides from the ab Initio Energy Surfaces of Terminally Blocked Glycine, Alanine, and Proline† , 2003 .

[63]  Elvira Guàrdia,et al.  Potential of mean force by constrained molecular dynamics: A sodium chloride ion-pair in water , 1991 .

[64]  A. Liwo,et al.  Theoretical calculations of heteroconjugation equilibrium constants in systems modeling acid-base interactions in side chains of biomolecules using the potential of mean force , 2004 .