Multiscale Modeling of Wobble to Watson-Crick-Like Guanine-Uracil Tautomerization Pathways in RNA

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Quantum Tunnelling Effects in the Guanine-Thymine Wobble Misincorporation via Tautomerism

The misincorporation of a noncomplementary DNA base in the polymerase active site is a critical source of replication errors that can lead to genetic mutations. In this work, we model the mechanism of wobble mispairing and the subsequent rate of misincorporation errors by coupling first-principles quantum chemistry calculations to an open quantum systems master equation. This methodology allows us to accurately calculate the proton transfer between bases, allowing the misincorporation and formation of mutagenic tautomeric forms of DNA bases. Our calculated rates of genetic error formation are in excellent agreement with experimental observations in DNA. Furthermore, our quantum mechanics/molecular mechanics model predicts the existence of a short-lived "tunnelling-ready" configuration along the wobble reaction pathway in the polymerase active site, dramatically increasing the rate of proton transfer by a hundredfold, demonstrating that quantum tunnelling plays a critical role in determining the transcription error frequency of the polymerase.