Exploration of Biochemical Reactivity with a QM/MM Growing String Method

03 October 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

In this work, we have implemented the single-ended growing string method within our in-house QM/MM package, QoMMMa. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations – the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, and more than 88% of the pathways generated were deemed to be of production quality. Compared to using adiabatic mapping for similar problems, we found that our technique is slightly less efficient, but readily produces better pathways and requires less user intervention. Given that we generated a relatively large number of pathways to fully validate our technique, we have also investigated the underlying distribution of QM/MM barrier heights. We show that this distribution is not strictly Gaussian and discuss the implications for studying reactivity with static QM/MM approaches.

Keywords

QM/MM
Minimum Energy Path
Reaction Pathway

Supplementary materials

Title
Description
Actions
Title
Supporting Information
Description
Supporting information document containing any additional methodology, tables, figures, and forcefield parameters.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.