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.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations

Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field

Improved Peptide and Protein Torsional Energetics with the OPLS-AA Force Field

Efficient Generation of Torsional Energy Profiles by Multifidelity Gaussian Processes for Hindered-Rotor Corrections

Automated Force Field Developer and Optimizer Platform: Torsion Reparameterization

A Rotational/Roto-translational Constraint Method for Condensed Matter

Optimal Solution to the Torsional Coefficient Fitting Problem in Force Field Parametrization

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Molecular Dynamics of Steroidal Rotors Probed by Theoretical, Spectroscopic and Dielectric Methods

pmx: Automated protein structure and topology generation for alchemical perturbations

GROMACS: Fast, flexible, and free

MoViES: Molecular Vibrations Evaluation Server for Analysis of Fluctuational Dynamics of Proteins and Nucleic Acids

Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

Coupled Two-Dimensional Main-Chain Torsional Potential for Protein Dynamics II: Performance and Validation.

A Random Rotor Molecule: Vibrational Analysis and Molecular Dynamics Simulations.

Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles

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

Identification of Torsional Modes in Complex Molecules Using Redundant Internal Coordinates: the Multistructural Method with Torsional Anharmonicity with a Coupled Torsional Potential and Delocalized Torsions.

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field