Determination of effective potentials for the stretching of C(α) ⋯ C(α) virtual bonds in polypeptide chains for coarse-grained simulations of proteins from ab initio energy surfaces of N-methylacetamide and N-acetylpyrrolidine.

The potentials of mean force (PMF's) for the deformation of the C(α) ⋯ C(α) virtual bonds in polypeptide chains were determined from the diabatic energy surfaces of N-methylacetamide (modeling regular peptide groups) and N-acetylpyrrolidine (modeling the peptide groups preceding proline), calculated at the Møller-Plesset (MP2) ab initio level of theory with the 6-31G(d,p) basis set. The energy surfaces were expressed in the C(α) ⋯ C(α) virtual-bond length (d) and the H-N-C(α) ⋯ C' improper dihedral angle (α) that describes the pyramidicity of the amide nitrogen, or in the C(α)-C'(O)-N-C(α) dihedral angle (ω) and the angle α. For each grid point, the potential energy was minimized with respect to all remaining degrees of freedom. The PMF's obtained from the (d, α) energy surfaces produced realistic free-energy barriers to the trans-cis transition (10 kcal/mol and 13 kcal/mol for the regular and proline peptide groups, respectively, compared to 12.6 - 13.9 kcal/mol and 17.3 - 19.6 kcal/mol determined experimentally for glycylglycine and N-acylprolines, respectively), while those obtained from the (ω, α) energy maps produced either low-quality PMF curves when direct Boltzmann summation was implemented to compute the PMF's or too-flat curves with too-low free-energy barriers to the trans-cis transition if harmonic extrapolation was used to estimate the contributions to the partition function. An analytical bimodal logarithmic-Gaussian expression was fitted to the PMF's, and the potentials were implemented in the UNRES force field. Test Langevin-dynamics simulations were carried out for the Gly-Gly and Gly-Pro dipeptides, which showed a 10(6)-fold increase of the simulated rate of the trans-cis isomerization with respect to that measured experimentally; effectively the same result was obtained with the analytical Kramers theory of reaction rate applied to the UNRES representation of the peptide groups. Application of Kramers' theory to compute the rate constants from the all-atom ab initio energy surfaces of the model compounds studied resulted in isomerization rates close to the experimental values, which demonstrates that the increase of the isomerization rate in UNRES simulations results solely from averaging out the secondary degrees of freedom.

[1]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[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]  A. Liwo,et al.  Ab initio simulations of protein-folding pathways by molecular dynamics with the united-residue model of polypeptide chains. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[4]  K. Orrell,et al.  Evaluation of the thermodynamic parameters for the interconversion of the conformational isomers of two N-acyl-prolines by nuclear magnetic resonance spectroscopy , 1971 .

[5]  R. Huber,et al.  Accurate Bond and Angle Parameters for X-ray Protein Structure Refinement , 1991 .

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

[7]  B. Zwolinski,et al.  High-precision viscosity of supercooled water and analysis of the extended range temperature coefficient , 1971 .

[8]  Y. Kang,et al.  Cis-trans isomerization and puckering of proline residue. , 2004, Biophysical chemistry.

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

[10]  D. Thirumalai,et al.  Protein folding kinetics: timescales, pathways and energy landscapes in terms of sequence-dependent properties. , 1996, Folding & design.

[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]  G. Voth.  Coarse-Graining of Condensed Phase and Biomolecular Systems , 2008 .

[13]  A. Liwo,et al.  A method for optimizing potential-energy functions by a hierarchical design of the potential-energy landscape: Application to the UNRES force field , 2002, Proceedings of the National Academy of Sciences of the United States of America.

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

[15]  A. Liwo,et al.  Molecular dynamics with the united-residue model of polypeptide chains. II. Langevin and Berendsen-bath dynamics and tests on model alpha-helical systems. , 2005, The journal of physical chemistry. B.

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

[17]  R. Kubo.  GENERALIZED CUMULANT EXPANSION METHOD , 1962 .

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

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

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

[21]  Adam Liwo,et al.  Implementation of molecular dynamics and its extensions with the coarse-grained UNRES force field on massively parallel systems; towards millisecond-scale simulations of protein structure, dynamics, and thermodynamics. , 2010, Journal of chemical theory and computation.

[22]  H. Halvorson,et al.  Consideration of the Possibility that the slow step in protein denaturation reactions is due to cis-trans isomerism of proline residues. , 1975, Biochemistry.

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

[24]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[25]  J. Skolnick,et al.  Discretized model of proteins. I. Monte Carlo study of cooperativity in homopolypeptides , 1992 .

[26]  A. Liwo,et al.  A united‐residue force field for off‐lattice protein‐structure simulations. II. Parameterization of short‐range interactions and determination of weights of energy terms by Z‐score optimization , 1997 .

[27]  H. Scheraga,et al.  Proline cis-trans isomerization and protein folding. , 2002, Biochemistry.

[28]  Adam Liwo,et al.  Relation between free energy landscapes of proteins and dynamics. , 2010, Journal of chemical theory and computation.

[29]  H. Kramers.  Brownian motion in a field of force and the diffusion model of chemical reactions , 1940 .

[30]  Michael R. Shirts,et al.  Atomistic protein folding simulations on the submillisecond time scale using worldwide distributed computing. , 2003, Biopolymers.

[31]  P. Gennes.  Scaling Concepts in Polymer Physics , 1979 .

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

[33]  Vijay S. Pande,et al.  Accelerating molecular dynamic simulation on graphics processing units , 2009, J. Comput. Chem..

[34]  S. Asher,et al.  UV resonance Raman studies of cis-to-trans isomerization of glycyglycine derivatives , 1999 .

[35]  A. Liwo,et al.  Molecular dynamics with the united-residue model of polypeptide chains. I. Lagrange equations of motion and tests of numerical stability in the microcanonical mode. , 2005, The journal of physical chemistry. B.

[36]  J. P. Grossman,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2008, CACM.

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

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

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

[40]  Young Kee Kang,et al.  Imide Cis−Trans Isomerization of N-Acetyl-N‘-methylprolineamide and Solvent Effects , 1999 .

[41]  Marek Cieplak,et al.  Thermal folding and mechanical unfolding pathways of protein secondary structures , 2001, Proteins.

[42]  R. Jakob,et al.  Elimination of a cis-proline-containing loop and turn optimization stabilizes a protein and accelerates its folding. , 2010, Journal of molecular biology.

[43]  R. Raines,et al.  Thermodynamic Origin of Prolyl Peptide Bond Isomers. , 1993, Tetrahedron letters.