Protocol for Computational Enzymatic Reactivity Based on Geometry Optimisation

Abstract

Enzymes play a biologically essential role in performing and controlling an important share of the chemical processes occurring in life. However, despite their critical role in nature, attaining a clear understanding of the way an enzyme acts is still cumbersome. Computational enzymology is playing an increasingly important role in this field of research. It allows the elucidation of a complete and detailed mechanism of an enzymatic reaction, including the characterization of reaction intermediates and transition states from both structural and energetic points of view, which is something that no other single experiment can produce alone. In this review, we present a general computational strategy to study enzymatic mechanisms based on adiabatic mapping and free geometry optimization. These methods allow chemical reactions to be studied with high theoretical levels, and allow a more exhaustive exploration of the chemical reactional space than other available methods, albeit being limited to the extent that they explore the enzyme conformational space. Special attention is given to the choice of the theoretical levels, as well as describing the model systems that are currently used to study enzymatic reactions. With this, we aim to provide a good introduction for non-specialised users in this field of research. We also provide a selection of hand-picked examples from our own work that illustrate the power of computational enzymology to study catalytic mechanisms. Some of these studies constitute pioneering work in the field that were later validated by experimental means.

Conformational diversity induces nanosecond-timescale chemical disorder in the HIV-1 protease reaction pathway

“…Finally, it is important to discuss the use of the X-ray conformation (the most abundant conformation) for the determination of the chemical mechanism and barrier height. A very large body of studies 19,[57][58][59][60][61][62][63] shows that a correct X-ray conformation almost always provides barrier heights in agreement with experiments, and chemical mechanisms that are widely accepted to also be correct. This is particularly true for Xray structures co-crystallized with a substrate/transition state analogue.…”

Solving the Catalytic Mechanism of Tryptophan Synthase: an Emergent Drug Target in the Treatment of Tuberculosis

“…The ZPE, thermal, and entropic energy corrections were calculated at 310.15 K and 1.0 bar. A more detailed overview about the methodology that was followed in this study can be found in the literature …”

“…A more detailed overview about the methodology that was followed in this study can be found in the literature. [44] The final energies of minima in all steps were calculated using a more complete basis set, 6-311 + + G(2df,pd) in the HL layer. The activation and reaction free energies were calculated through the difference between the Gibbs free energy of TS and reactant, or product and reactant, respectively.…”

molUP: A VMD plugin to handle QM and ONIOM calculations using the gaussian software

“…The user can also save his/her own template and include it in the dropdown menu to have it accessible when future input files will be prepared. For example, in ONIOM calculations, the user can easily interchange between mechanical and electrostatic embedding type of calculations, which only differ in a single keyword (for more details see Supporting Information) …”