User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse




Journal Article

On the faithfulness of molecular mechanics representations of proteins towards quantum-mechanical energy surfaces


König,  Gerhard
Research Department Thiel, Max-Planck-Institut für Kohlenforschung, Max Planck Society;
Laboratory of Physical Chemistry, ETH Zürich;

External Ressource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available

König, G., & Riniker, S. (2020). On the faithfulness of molecular mechanics representations of proteins towards quantum-mechanical energy surfaces. Interface Focus, 10(6): 20190121. doi:10.1098/rsfs.2019.0121.

Cite as: http://hdl.handle.net/21.11116/0000-0007-6F9D-D
Force fields based on molecular mechanics (MM) are the main computational tool to study the relationship between protein structure and function at the molecular level. To validate the quality of such force fields, high-level quantum-mechanical (QM) data are employed to test their capability to reproduce the features of all major conformational substates of a series of blocked amino acids. The phase-space overlap between MM and QM is quantified in terms of the average structural reorganization energies over all energy minima. Here, the structural reorganization energy is the MM potential-energy difference between the structure of the respective QM energy minimum and the structure of the closest MM energy minimum. Thus, it serves as a measure for the relative probability of visiting the QM minimum during an MM simulation. We evaluate variants of the AMBER, CHARMM, GROMOS and OPLS biomolecular force fields. In addition, the two blocked amino acids alanine and serine are used to demonstrate the dependence of the measured agreement on the QM method, the phase, and the conformational preferences. Blocked serine serves as an example to discuss possible improvements of the force fields, such as including polarization with Drude particles, or using tailored force fields. The results show that none of the evaluated force fields satisfactorily reproduces all energy minima. By decomposing the average structural reorganization energies in terms of individual energy terms, we can further assess the individual weaknesses of the parametrization strategies of each force field. The dominant problem for most force fields appears to be the van der Waals parameters, followed to a lesser degree by dihedral and bonded terms. Our results show that performing a simple QM energy optimization from an MM-optimized structure can be a first test of the validity of a force field for a particular target molecule.