Abstract: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.

Potential Energy Surfaces Sampled in Cremer–Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

Accurate Potential Energy Surfaces for the Excited State of CF2 Molecule

On the Potential Energy Surface of the Pyrene Dimer

Data Quality, Data Sampling and Data Fitting: A Tutorial Guide for Constructing Full-Dimensional Accurate Potential Energy Surfaces (PESs) of Molecules and Reactions

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Theoretical assessment of vibrationally resolved C1s X-ray photoelectron spectra of simple cyclic molecules

Comparison of curvilinear coordinates within vibrational structure calculations based on automatically generated potential energy surfaces

DFT-Based Permutationally Invariant Polynomial Potentials Capture the Twists and Turns of C14H30

Advanced Potential Energy Surfaces for Molecular Simulation

Ab initio potential energy surfaces calculation via restricted Hartree–Fock for molecular dynamics simulation: a comprehensive review

Stationary Points on Potential Energy Surface of Cyclic C3H3 with Coupled-Cluster Approaches and Density Functional Theory.

Computational Design of a Tetrapericyclic Cycloaddition and the Nature of Potential Energy Surfaces with Multiple Bifurcations.

Explicitly correlated MRCI-F12 potential energy surfaces for methane fit with several permutation invariant schemes and full-dimensional vibrational calculations

DFT-Based Permutationally Invariant Polynomial Potentials Capture the Twists and Turns of C 14 H 30

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

Optimizing The Structures Of Minimum And Transition State On The Free Energy Surface

FESTA: A Polygon-Based Approach for Extracting Relevant Structures from Free Energy Surfaces Obtained in Molecular Simulations

The Potential Surface in the Ground Electronic State of HCP with the Isomerization Process: the Validity of Calculating Potential Surface with DFT Methods

Theoretical Study of the Potential-Energy Surface of C2NP

Combined Self-Consistent-Field and Spin-Flip Tamm-Dancoff Density Functional Approach to Potential Energy Surfaces for Photochemistry

Anharmonic Terms of the Potential Energy Surface: A Group Theoretical Approach