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.

Relative energetics of acetyl-histidine protomers with and without Zn2+ and a benchmark of energy methods

Multireference protonation energetics of a dimeric model of nitrogenase iron-sulfur clusters

Upgrading and Validation of the AMBER Force Field for Histidine and Cysteine Zinc(II)-Binding Residues in Sites with Four Protein Ligands

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Cysteine and histidine pKa prediction: alchemical free energy calculations and force field reparameterization

Comparison of different approaches to derive classical bonded force-field parameters for a transition metal cofactor: a case study for non-heme iron site of ectoine synthase

Upgraded AMBER Force Field for Zinc-Binding Residues and Ligands for Predicting Structural Properties and Binding Affinities in Zinc-Proteins

Reliable Dimerization Energies for Modeling of Supramolecular Junctions

Comparison of Empirical Zn2+ Models in Protein–DNA Complexes

Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo

Accurate Prediction of Noncovalent Interaction Energies with the Effective Fragment Potential Method: Comparison of Energy Components to Symmetry-Adapted Perturbation Theory for the S22 Test Set

Vibrational analysis of amino acids and short peptides in hydrated media. VII. Energy landscapes, energetic and geometrical features of L-histidine with protonated and neutral side chains

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins

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

Proton affinity Revisited: Benchmarking computational approaches for accurate predictions

Systematic Derivation of AMBER Force Field Parameters Applicable to Zinc-Containing Systems.

A hierarchical method for generating low‐energy conformers of a protein‐ligand complex

Benchmark ab initio characterization of the multi-channel Cl + CH3X [X = F, Cl, Br, I] reactive potential energy surfaces

ChemInform Abstract: Redox Potentials and Acidity Constants from Density Functional Theory Based Molecular Dynamics

The multiple conformational charge states of zinc(II) coordination by 2His-2Cys oligopeptide investigated by ion mobility-mass spectrometry, density functional theory and theoretical collision cross sections

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations