DNA Polymerase λ Active Site Favors a Mutagenic Mispair between the Enol Form of Deoxyguanosine Triphosphate Substrate and the Keto Form of Thymidine Template: A Free Energy Perturbation Study

Correct Nucleotide Selection Is Confined at the Binding Site of Polymerase Enzymes

DNA polymerases (Pols) add incoming nucleotides (deoxyribonucleoside triphosphate (dNTPs)) to growing DNA strands, a crucial step for DNA synthesis. The insertion of correct (vs incorrect) nucleotides relates to Pols' fidelity, which defines Pols' ability to faithfully replicate DNA strands in a template-dependent manner. We and others have demonstrated that reactant alignment and correct base pairing at the Pols catalytic site are crucial structural features to fidelity. Here, we first used equilibrium molecular simulations to demonstrate that the local dynamics at the protein-DNA interface in the proximity of the catalytic site is different when correct vs incorrect dNTPs are bound to polymerase β (Pol β). Formation and dynamic stability of specific interatomic interactions around the incoming nucleotide influence the overall binding site architecture. This explains why certain Pols' mutants can affect the local catalytic environment and influence the selection of correct vs incorrect nucleotides. In particular, this is here demonstrated by analyzing the interaction network formed by the residue R283, whose mutant R283A has an experimentally measured lower capacity of differentiating correct (G:dCTP) vs incorrect (G:dATP) base pairing in Pol β. We also used alchemical free-energy calculations to quantify the G:dCTP →G:dATP transformation in Pol β wild-type and mutant R283A. These results correlate well with the experimental trend, thus corroborating our mechanistic insights. Sequence and structural comparisons with other Pols from the same family suggest that these findings may also be valid in similar enzymes.

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.

Rare Tautomers of Artificially Expanded Genetic Letters and their Effects on the Base pair Stabilities

To expand the existing genetic letters, it is necessary to design robust nucleotides that can function naturally in living cells. Therefore, it is desirable to examine the roles of recently proposed second-generation artificially expanded genetic letters in producing stable duplex DNA. Here, a reliable dispersion-corrected density functional theory method is used to understand the electronic structures and properties of different rare tautomers of proposed expanded genetic letters and their effects on the base pair stabilities in the duplex DNA. It is found that the rare tautomers are not only stable in the aqueous medium but can also base pair with natural bases to produce stable mispairs. Except for J and V, all the artificial genetic letters are found to produce mispairs that are about 1-7 kcal/mol more stable than their complementary counterparts. They are also appreciably more stable than the naturally occurring G:C, A:T, and G:T pairs. The higher base pair stabilities are found to be mainly because of the polarity of monomers and attractive electrostatic interactions.

Constraints on error rate revealed by computational study of G•U tautomerization in translation

In translation, G•U mismatch in codon-anticodon decoding is an error hotspot likely due to transition of G•U from wobble (wb) to Watson-Crick (WC) geometry, which is governed by keto/enol tautomerization (wb-WC reaction). Yet, effects of the ribosome on the wb-WC reaction and its implications for decoding mechanism remain unclear. Employing quantum-mechanical/molecular-mechanical umbrella sampling simulations using models of the ribosomal decoding site (A site) we determined that the wb-WC reaction is endoergic in the open, but weakly exoergic in the closed A-site state. We extended the classical ‘induced-fit’ model of initial selection by incorporating wb-WC reaction parameters in open and closed states. For predicted parameters, the non-equilibrium exoergic wb-WC reaction is kinetically limited by the decoding rates. The model explains early observations of the WC geometry of G•U from equilibrium structural studies and reveals discrimination capacity for the working ribosome operating at non-equilibrium conditions. The equilibration of the exoergic wb-WC reaction counteracts the equilibration of the open-closed transition of the A site, constraining the decoding accuracy and potentially explaining the persistence of the G•U as an error hotspot. Our results unify structural and mechanistic views of codon-anticodon decoding and generalize the ‘induced-fit’ model for flexible substrates.

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

Energetically unfavorable Watson–Crick (WC)-like tautomeric forms of nucleobases are known to introduce spontaneous mutations, and contribute to replication, transcription, and translation errors. Recent NMR relaxation dispersion techniques were able to show that wobble (w) G•U mispair exists in equilibrium with the short-lived, low-population WC-like enolic tautomers. Presently, we have investigated the wG•U → WC-like enolic reaction pathway using various theoretical methods: quantum mechanics (QM), molecular dynamics (MD), and combined quantum mechanics/molecular mechanics (QM/MM). The previous studies on QM gas phase calculations were inconsistent with experimental data. We have also explored the environmental effects on the reaction energies by adding explicit water. While the QM-profile clearly becomes endoergic in the presence of water, the QM/MM-profile remains consistently endoergic in the presence and absence of water. Hence, by including microsolvation and QM/MM calculations, the experimental data can be explained. For the G•Uenol→ Genol•U pathway, the latter appears to be energetically more favorable throughout all computational models. This study can be considered as a benchmark of various computational models of wG•U to WC-like tautomerization pathways with and without the environmental effects, and may contribute on further studies of other mispairs as well.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Unexpected Routes of the Mutagenic Tautomerization of the T Nucleobase in the Classical A·T DNA Base Pairs: A QM/QTAIM Comprehensive View

In this paper using quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) in the continuum with ε = 1, we have theoretically demonstrated for the first time that revealed recently highly-energetic conformers of the classical A·T DNA base pairs - Watson-Crick [A·T(wWC)], reverse Watson-Crick [A·T(wrWC)], Hoogsteen [A·T(wH)] and reverse Hoogsteen [A·T(wrH)] - act as intermediates of the intrapair mutagenic tautomerization of the T nucleobase owing to the novel tautomerisation pathways: A·T(wWC)↔A·T*(w⊥ WC); A·T(wrWC)↔A·T O 2 * (w⊥ rWC); A·T(wH)↔A·T*(w⊥ H); A·T(wrH)↔A·T O 2 * (w⊥ rH). All of them occur via the transition states as tight ion pairs (A+, protonated by the N6H2 amino group)·(T-, deprotonated by the N3H group) with quasi-orthogonal geometry, which are stabilized by the participation of the strong (A)N6+H···O4-/O2-(T) and (A)N6+H···N3-(T) H-bonds. Established tautomerizations proceed through a two-step mechanism of the protons moving in the opposite directions along the intermolecular H-bonds. Initially, proton moves from the N3H imino group of T to the N6H2 amino group of A and then subsequently from the protonated N6+H3 amino group of A to the O4/O2 oxygen atom of T, leading to the products - A·T*(w⊥ WC), A·T O 2 * (w⊥ rWC), A·T*(w⊥ H), and A·T O 2 * (w⊥ rH), which are substantially non-planar, conformationally-labile complexes. These mispairs are stabilized by the participation of the (A)N6H/N6H'···N3(T) and (T)O2H/O4H···N6(A) H-bonds, for which the pyramidalized amino group of A is their donor and acceptor. The Gibbs free energy of activation of these mutagenic tautomerizations lies in the range of 27.8-29.8 kcal·mol-1 at T = 298.15 K in the continuum with ε = 1.

Non-dissociative structural transitions of the Watson-Crick and reverse Watson-Crick А·Т DNA base pairs into the Hoogsteen and reverse Hoogsteen forms

In this study it was theoretically shown that discovered by us recently (Brovarets' et al., Frontiers in Chemistry, 2018, 6:8; doi: 10.3389/fchem.2018.00008) high-energetical, significantly non-planar (symmetry C₁), short-lived wobbled conformers of the classical Watson-Crick А·Т(WC), reverse Watson-Crick А·Т(rWC), Hoogsteen А·Т(Н) and reverse Hoogsteen А·Т(rН) DNA base pairs are the intermediates of their pairwise А∙Т(WC)/А∙Т(rWC) ↔ А∙Т(H)/А∙Т(rH) conformational transformations. These transitions do not require for their realization the energy-consumable anisotropic rotation of the amino group of A around the exocyclic C6-N6 bond. They are controlled by the non-planar transition states with quasi-orthogonal geometry (symmetry C₁) joined by the single intermolecular (Т)N3H···N6(А) H-bond (~4 kcal∙mol^-1). The Gibbs free energies of activation for these non-dissociative, dipole-active conformational transitions consist 7.33 and 7.81 kcal∙mol^-1, accordingly. Quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) have been performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory in the continuum with ε = 4 under normal conditions.

Novel pathway for mutagenic tautomerization of classical А∙Т DNA base pairs via sequential proton transfer through quasi-orthogonal transition states: A QM/QTAIM investigation

In this paper we have theoretically predicted a novel pathway for the mutagenic tautomerization of the classical A∙T DNA base pairs in the free state, the Watson-Crick A·Т(WC), reverse Watson-Crick A·Т(rWC), Hoogsteen A·Т(H) and reverse Hoogsteen A·Т(rH) pairs, via sequential proton transfer accompanied by a significant change in the mutual orientation of the bases. Quantum-mechanical (QM) calculations were performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level in vacuum phase, along with Bader's quantum theory of Atoms in Molecules (QTAIM). These processes involve transition states (TSs) with quasi-orthogonal structures (symmetry C1), which are highly polar, tight ion pairs (A-, N6H2-deprotonated)∙(T+, O4/O2-protonated). Gibbs free energies of activation for the A∙T(WC) / A∙T(rWC) ↔ A*∙Т(rwWC) / A*∙Т(wWC) tautomeric transitions (~43.5 kcal∙mol-1) are lower than for the A∙T(H) / A∙T(rH) ↔ A*N7∙Т(rwH) / A*N7∙Т(wH) tautomerisations (~53.0 kcal∙mol-1) (rare tautomers are marked by an asterisk; w-wobble configured tautomerisation products). The (T)N3+H⋯N1-(A), (T)O4+H⋯N1-(A) / (T)N3+H⋯N1-(A) and (T)O2+H⋯N1-(A) H-bonds are found in the transition states TSA-·T+A·T(WC)↔A*·T(rwWC) / TSA-·T+A·T(rWC)↔A*·T(wWC). However, in the transition state TSA-·T+A·Т(H)↔A*N7·T(rwH) / TSA-·T+A·Т(rH)↔A*N7·T(wH), the (T)N3+H⋯N7-(A), (T)O4+H⋯N7-(A) / (T)N3+H⋯N7-(A) and (T)O2+H⋯N7-(A) H-bonds are supplemented by the attractive (T)O4+/O2+⋯N6-(A) van der Waals contacts. It was demonstrated that the products of the tautomerization of the classical A∙T DNA base pairs-A*∙Т(rwWC), A*N7∙Т(rwH) and A*N7∙Т(wH) (symmetry Cs)-further transform via double proton transfer into the energetically favorable wobble A∙T*(rwWC), A∙T*(rwH) and A∙T*O2(wH) base mispairs (symmetry Cs).

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Dynamic basis for dG•dT misincorporation via tautomerization and ionization

Tautomeric and anionic Watson-Crick-like mismatches have important roles in replication and translation errors through mechanisms that are not fully understood. Here, using NMR relaxation dispersion, we resolve a sequence-dependent kinetic network connecting G•T/U wobbles with three distinct Watson-Crick mismatches: two rapidly exchanging tautomeric species (Genol•T/UG•Tenol/Uenol; population less than 0.4%) and one anionic species (G•T-/U-; population around 0.001% at neutral pH). The sequence-dependent tautomerization or ionization step was inserted into a minimal kinetic mechanism for correct incorporation during replication after the initial binding of the nucleotide, leading to accurate predictions of the probability of dG•dT misincorporation across different polymerases and pH conditions and for a chemically modified nucleotide, and providing mechanisms for sequence-dependent misincorporation. Our results indicate that the energetic penalty for tautomerization and/or ionization accounts for an approximately 10-2 to 10-3-fold discrimination against misincorporation, which proceeds primarily via tautomeric dGenol•dT and dG•dTenol, with contributions from anionic dG•dT- dominant at pH 8.4 and above or for some mutagenic nucleotides.

Computational Study of Uracil Tautomeric Forms in the Ribosome: The Case of Uracil and 5-Oxyacetic Acid Uracil in the First Anticodon Position of tRNA

Tautomerism is important in many biomolecular interactions, not least in RNA biology. Crystallographic studies show the possible presence of minor tautomer forms of transfer-RNA (tRNA) anticodon bases in the ribosome. The hydrogen positions are not resolved in the X-ray studies, and we have used ab initio calculations and molecular dynamics simulations to understand if and how the minor enol form of uracil (U), or the modified uracil 5-oxyacetic acid (cmo⁵U), can be accommodated in the tRNA-messenger-RNA interactions in the ribosome decoding center. Ab initio calculations on isolated bases show that the modification affects the keto-enol equilibrium of the uracil base only slightly; the keto form is dominant (>99.99%) in both U and cmo⁵U. Other factors such as interactions with the surrounding nucleotides or ions would be required to shift the equilibrium toward the enol tautomer. Classical molecular simulations show a better agreement with the X-ray structures for the enol form, but free energy calculations indicate that the most stable form is the keto. In the ribosome, the enol tautomers of U and cmo⁵U pair with a guanine forming two hydrogen bonds, which do not involve the enol group. The oxyacetic acid modification has a minor effect on the keto-enol equilibrium.

The A·T(rWC)/A·T(H)/A·T(rH) ↔ A·T*(rwWC)/A·T*(wH)/A·T*(rwH) mutagenic tautomerization via sequential proton transfer: a QM/QTAIM study

In this study for the first time we have revealed by QM and QTAIM calculations at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory the novel routes of the mutagenic tautomerization of three biologically important A·T DNA base pairs – reverse Watson–Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) – followed by their rebuilding into the wobble (w) A·T*(rw_WC), A·T*(w_H) and A·T*(rw_H) base mispairs by the participation of the mutagenic tautomers of the DNA bases (denoted by asterisk) and vice versa, thus complementing the physico-chemical property of the canonical A·T(WC) Watson–Crick DNA base pair reported earlier (Brovarets' et al., RSC Adv., 2015, 5, 99594–99605). These non-dissociative tautomeric transformations in the classical A·T(rWC), A·T(H) and A·T(rH) DNA base pairs proceed similarly to the canonical A·T(WC) DNA base pair via the intrapair sequential proton transfer with shifting towards major or minor grooves of DNA followed by further double proton transfer along the intermolecular H-bonds and are controlled by the plane symmetric and highly stable transition states – tight ion pairs formed by the A⁺ nucleobase, protonated by the N1/N7 nitrogen atoms, and T⁻ nucleobase, deprotonated by the N3H imino group. Comparison of the estimated populations of the tautomerised states (10⁻²¹ to 10⁻¹⁴) with similar characteristics for the canonical A·T(WC) DNA base pair (10⁻⁸ to 10⁻⁷) leads authors to the conclusion, that only a base pair with WC architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development. Among all four classical DNA base pairs, only A·T(WC) DNA base pair can ensure the proper rate of the spontaneous point errors of replication in DNA.

We discovered tautomeric wobbling of the classical A·T DNA base pairs. This data evidence, that only a base pair with Watson–Crick architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development.