Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G*·C* DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson-Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices.Communicated by Ramaswamy H. Sarma.

Atomistic mechanisms of the tautomerization of the G·C base pairs through the proton transfer: quantum-chemical survey

This study is devoted to the investigation of the G·C*^t_O2(WC)↔G*_NH3·C*^t(WC), G·C*_O2(WC)↔G*_NH3·C*(WC) and G*·C*_O2(WC)↔G*_NH3·C(w_WC)_↓ tautomerization reactions occurring through the proton transfer, obtained at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory in gas phase under normal conditions ('WC' means base pair in Watson-Crick configuration, T=298.15 K). These reactions lead to the formation of the G*_NH3·C*^t(WC), G*_NH3·C*(WC) and G*_NH3·C(w_WC)_↓ base pairs by the participation of the G*_NH3 base with NH₃ group. Gibbs free energies of activation for these reactions are 6.43, 11.00 and 1.63 kcal·mol^-1, respectively. All of these tautomerization reactions are dipole active. Finally, we believe that these non-dissociative processes, which are tightly connected with the tautomeric transformations of the G·C base pairs, play an outstanding role in supporting of the spatial structure of the DNA and RNA molecules with various functional purposes.

Kinetic Study of Transition Mutations from G-C to A-T Base Pairs in Watson-Crick DNA Base Pairs: Double Proton Transfers

According to the Löwdin model [ Rev. Mod. Phys. 1963, 35, 724-732], the Watson-Crick guanine-cytosine (G-C) base pair is tautomerized (G*-C*) with a small probability and then replication of G*-C* produces G*-thymine (T) and adenine (A)-C* base pairs. On the basis of this model and our previous work [ J. Phys. Chem. B 2020, 124, 1715-1722], we first calculated the intrinsic reaction coordinates from G*-T to G-T* using density functional theory and evaluated the probability of G*-T tautomerization to G-T* by double proton transfer (DPT) on the basis of the transition state theory. Similarly, we calculated the probability of A-C* tautomerization to A*-C by DPT. Then, according to these probabilities, we calculated the probability of transition mutations from G-C to A-T after 2 replications. The calculated probability was 1.31 × 10^-8, a value consistent with the mutation rate previously reported by Drake et al. [ Proc. Natl. Acad. Sci. U.S.A. 1991, 88, 7160-7164]. Our results suggest that DPT is one cause of the G-C → A-T transition. To investigate differences in the optical properties between G*-T and G-T* and between A-C* and A*-C, we also evaluated the infrared absorption spectra and Raman intensities for these base pairs.

The Research of G-Motif Construction and Chirality in Deoxyguanosine Monophosphate Nucleotide Complexes

A detailed understanding of the mismatched base-pairing interactions in DNA will help reveal genetic diseases and provide a theoretical basis for the development of targeted drugs. Here, we utilized mononucleotide fragment to simulate mismatch DNA interactions in a local hydrophobic microenvironment. The bipyridyl-type bridging ligands were employed as a mild stabilizer to stabilize the GG mismatch containing complexes, allowing mismatch to be visualized based on X-ray crystallography. Five single crystals of 2′-deoxyguanosine–5′–monophosphate (dGMP) metal complexes were designed and obtained via the process of self-assembly. Crystallographic studies clearly reveal the details of the supramolecular interaction between mononucleotides and guest intercalators. A novel guanine–guanine base mismatch pattern with unusual (high anti)–(high anti) type of arrangement around the glycosidic angle conformations was successfully constructed. The solution state ¹H–NMR, ESI–MS spectrum studies, and UV titration experiments emphasize the robustness of this g–motif in solution. Additionally, we combined the methods of single-crystal and solution-, solid-state CD spectrum together to discuss the chirality of the complexes. The complexes containing the g–motif structure, which reduces the energy of the system, following the solid-state CD signals, generally move in the long-wave direction. These results provided a new mismatched base pairing, that is g–motif. The interaction mode and full characterizations of g–motif will contribute to the study of the mismatched DNA interaction.

Structural dynamics of double-stranded DNA with epigenome modification

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson-Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson-Crick base pair to allow base-flipping in living cells.

Novel mechanisms of the conformational transformations of the biologically important G·C nucleobase pairs in Watson–Crick, Hoogsteen and wobble configurations via the mutual rotations of the bases around the intermolecular H-bonds: a QM/QTAIM study

It was established conformational transformations of the G·C nucleobase pairs, occurring via the mutual rotation of the G and C bases around the intermolecular H-bonds.

A Quantum-Mechanical Looking Behind the Scene of the Classic G·C Nucleobase Pairs Tautomerization

For the first time, at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory, a comprehensive quantum-mechanical investigation of the physico-chemical mechanism of the tautomeric wobblization of the four biologically-important G·C nucleobase pairs by the participation of the monomers in rare, in particular mutagenic, tautomeric forms (marked with an asterisk) was provided. These novel tautomeric transformations (wobblization or shifting of the bases within the pair) are intrinsically inherent properties of the G·C nucleobase pairs. In this study, we have obtained intriguing results, lying far beyond the existing representations. Thus, it was shown that Löwdin's G*·C*(WC) base pair does not tautomerize according to the wobblization mechanism. Tautomeric wobblization of the G*·C*(rWC) (relative Gibbs free energy ΔG = 0.00/relative electronic energy ΔE = 0.00 kcal·mol-1) ("r"-means the configuration of the base pair in reverse position; "WC"-the classic Watson-Crick configuration) and G*t·C*(H) (ΔG = -0.19/ΔE = 0.29 kcal·mol-1) ("H"-Hoogsteen configuration;"t" denotes the O6H hydroxyl group in the trans position) base pairs are preceded by the stages of the base pairs tautomerization by the single proton transfer (SPT). It was established that the G*t·C*(rH) (ΔG = 2.21/ΔE = 2.81 kcal·mol-1) base pair can be wobbled through two different pathways via the traditional one-stage mechanism through the TSs, which are tight G+·C- ion pairs, stabilized by the participation of only two intermolecular H-bonds. It was found out that the G·C base pair is most likely incorporated into the DNA/RNA double helix with parallel strands in the G*·C*(rWC), G·C*(rwwc), and G*·C(rwwc) ("w"-wobble configuration of the pair) tautomeric forms, which are in rapid tautomeric equilibrium with each other. It was proven that the G*·C*(rWC) nucleobase pair is also in rapid tautomeric equilibrium with the eight tautomeric forms of the so-called Levitt base pair. It was revealed that a few cases of tautomerization via the DPT of the nucleobase pairs by the participation of the C8H group of the guanine had occurred. The biological role of the obtained results was also made apparent.

Can modified DNA base pairs with chalcogen bonding expand the genetic alphabet? A combined quantum chemical and molecular dynamics simulation study

A comprehensive (DFT and MD) computational study is presented with the goal to design and analyze model chalcogen-bonded modified nucleobase pairs that replace one (i.e., A_XY:T, G:C_XY, G_XY:C) or two (G_XY:C_X'Y', X/X' = S, Se and Y/Y' = F, Cl, Br) Watson-Crick (WC) hydrogen bonds of the canonical A:T or G:C pair with chalcogen bond(s). DFT calculations on 18 base pair combinations that replace one WC hydrogen bond with a chalcogen bond reveal that the bases favorably interact in the gas phase (binding strengths up to -140 kJ mol^-1) and water (up to -85 kJ mol^-1). Although the remaining hydrogen bond(s) exhibits similar characteristics to those in the canonical base pairs, the structural features of the (Y-XO) chalcogen bond(s) change significantly with the identity of X and Y. The 36 doubly-substituted (G_XY:C_X'Y') base pairs have structural deviations from canonical G:C similar to those of the singly-substituted modifications (G:C_XY or G_XY:C). Furthermore, despite the replacement of two strong hydrogen bonds with chalcogen bonds, some G_XY:C_X'Y' pairs possess comparable binding energies (up to -132 kJ mol^-1 in the gas phase and up to -92 kJ mol^-1 in water) to the most stable G:C_XY or G_XY:C pairs, as well as canonical G:C. More importantly, G:C-modified pairs containing X = Se (high polarizability) and Y = F (high electronegativity) are the most stable, with comparable or slightly larger (by up to 13 kJ mol^-1) binding energies than G:C. Further characterization of the chalcogen bonding in all modified base pairs (AIM, NBO and NCI analyses) reveals that the differences in the binding energies of modified base pairs are mainly dictated by the differences in the strengths of their chalcogen bonds. Finally, MD simulations on DNA oligonucleotides containing the most stable chalcogen-bonded base pair from each of the four classifications (A_XY:T, G:C_XY, G_XY:C and G_XY:C_X'Y') reveal that the singly-modified G:C pairs best retain the local helical structure and pairing stability to a greater extent than the modified A:T pair. Overall, our study identifies two (G:C_SeF and G_SeF:C) promising pairs that retain chalcogen bonding in DNA and should be synthesized and further explored in terms of their potential to expand the genetic alphabet.

Structural Patterns and Stabilities of Hydrogen-Bonded Pairs Involving Ribonucleotide Bases and Arginine, Glutamic Acid, or Glutamine Residues of Proteins from Quantum Mechanical Calculations

Ribonucleotide:protein interactions play crucial roles in a number of biological processes. Unlike the RNA:protein interface where van der Waals contacts are prevalent, the recognition of a single ribonucleotide such as ATP by a protein occurs predominantly through hydrogen-bonding interactions. As a first step toward understanding the role of hydrogen bonding in ribonucleotide:protein recognition, the present work employs density functional theory to provide a detailed quantum-mechanical analysis of the structural and energetic characteristics of 18 unique hydrogen-bonded pairs involving the nucleobase/nucleoside moiety of four canonical ribonucleotides and the side chains of three polar amino-acid residues (arginine, glutamine, and glutamic acid) of proteins. In addition, we model five new pairs that are till now not observed in crystallographically identified ribonucleotide:protein complexes but may be identified in complexes crystallized in the future. We critically examine the characteristics of each pair in its ribonucleotide:protein crystal structure occurrence and (gas phase and water phase) optimized intrinsic structure. We further evaluated the interaction energy of each pair and characterized the associated hydrogen bonds using a number of quantum mechanics-based relationships including natural bond orbital analysis, quantum theory atoms in molecules analysis, Iogansen relationships, Nikolaienko-Bulavin-Hovorun relationships, and noncovalent interaction-reduced density gradient analysis. Our analyses reveal rich variability in hydrogen bonds in the crystallographic as well as intrinsic structure of each pair, which includes conventional O/N-H···N/O and C-H···O hydrogen bonds as well as donor/acceptor-bifurcated hydrogen bonds. Further, we identify five combinations of nucleobase and amino acid moieties; each of which exhibits at least two alternate (i.e., multimodal) structures that interact through the same nucleobase edge. In fact, one such pair exhibits four multimodal structures; one of which possesses unconventional "amino-acceptor" hydrogen bonding with comparable (-9.4 kcal mol^-1) strength to the corresponding conventional (i.e., amino:donor) structure (-9.2 kcal mol^-1). This points to the importance of amino-acceptor hydrogen bonds in RNA:protein interactions and suggests that such interactions must be considered in the future while studying the dynamics in the context of molecular recognition. Overall, our study provides preliminary insights into the intrinsic features of ribonucleotide:amino acid interactions, which may help frame a clearer picture of the molecular basis of RNA:protein recognition and further appreciate the role of such contacts in biology.