Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification

While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step towards quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region, and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol \textit{O}-methyltransferase, and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.

The enzyme mechanism of patchoulol synthase

Different mechanisms for the cyclisation of farnesyl pyrophosphate to patchoulol by the patchoulol synthase are discussed in the literature. They are based on isotopic labelling experiments, but the results from these experiments are contradictory. The present work reports on a reinvestigation of patchoulol biosynthesis by isotopic labelling experiments and computational chemistry. The results are in favour of a pathway through the neutral intermediates germacrene A and α-bulnesene that are both reactivated by protonation for further cyclisation steps, while previously discussed intra- and intermolecular hydrogen transfers are not supported. Furthermore, the isolation of the new natural product (2S,3S,7S,10R)-guaia-1,11-dien-10-ol from patchouli oil is reported.