The effect of tautomeric constant on the specificity of nucleotide incorporation during DNA replication: support for the rare tautomer hypothesis of substitution mutagenesis

How to measure the probability of uracil chain tautomers as nucleotide bases in RNA?

The current research focuses on exploring tautomerism in uracil. 47 tautomers were found that varied in significance in RNA and stability. To discover these molecules, diverse potential energy levels were explored, and corresponding transition states were found in these pathways. But the imperative thing that was taken note in this investigation is that for the first time, a method was detailed for the probability of forming distinctive molecules relative to each other. In this method, the conversion of uracil and its tautomers, which together turn into 47 molecules, was composed as a Markov chain. Then, the transition matrix was explained using its support, whose components are the probability of creating molecules from each step. At last, by multiplying this matrix by n times, the probability of forming different molecules was obtained. Moreover, by solving this matrix at different times, it is conceivable to appear which molecules can be converted to uracil sooner. It was appeared that a few tautomers act as transitory absorption point or temporary terminal states and other molecules, to begin with convert to these molecules before turning into uracil.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Insights into the substrate discrimination mechanisms of methyl-CpG-binding domain 4

G:T mismatches, the major mispairs generated during DNA metabolism, are repaired in part by mismatch-specific DNA glycosylases such as methyl-CpG-binding domain 4 (MBD4) and thymine DNA glycosylase (TDG). Mismatch-specific DNA glycosylases must discriminate the mismatches against million-fold excess correct base pairs. MBD4 efficiently removes thymine opposite guanine but not opposite adenine. Previous studies have revealed that the substrate thymine is flipped out and enters the catalytic site of the enzyme, while the estranged guanine is stabilized by Arg468 of MBD4. To gain further insights into the mismatch discrimination mechanism of MBD4, we assessed the glycosylase activity of MBD4 toward various base pairs. In addition, we determined a crystal structure of MBD4 bound to T:O6-methylguanine-containing DNA, which suggests the O6 and N2 of purine and the O4 of pyrimidine are required to be a substrate for MBD4. To understand the role of the Arg468 finger in catalysis, we evaluated the glycosylase activity of MBD4 mutants, which revealed the guanidinium moiety of Arg468 may play an important role in catalysis. D560N/R468K MBD4 bound to T:G mismatched DNA shows that the side chain amine moiety of the Lys stabilizes the flipped-out thymine by a water-mediated phosphate pinching, while the backbone carbonyl oxygen of the Lys engages in hydrogen bonds with N2 of the estranged guanine. Comparison of various DNA glycosylase structures implies the guanidinium and amine moieties of Arg and Lys, respectively, may involve in discriminating between substrate mismatches and nonsubstrate base pairs.

Investigation of tautomerism of 1,3,5-triazine derivative, stability, and acidity of its tautomers from density functional theory

Recent studies have identified N2,N4-bis(4-fluorophenethyl)-N6-(3-(dimethylamino)propyl)-1,3,5-triazine-2,4,6-triamine (1TZ(7,8,9)) as a potent, pure antagonist that inhibits thermosensory transient receptor potential vanilloid 1 channel (TRPV1) channel activity. This study provides theoretical data on the stability and acidity of the tautomers of this molecule. We show that this triazine can exist as three predominant tautomers (2TZ(5,7,8), 4TZ(3,7,9), 7TZ(1,8,9)). In the aqueous phase, equilibrium constants calculations show that only the tautomeric equilibria between 1TZ(7,8,9) and the three most stable triazines can be present which suggests that these three tautomeric equilibria would be the basis of 1TZ(7,8,9)'s biological activity.

Probing conformational transitions towards mutagenic Watson-Crick-like G·T mismatches using off-resonance sugar carbon R1ρ relaxation dispersion

NMR off-resonance R_1ρ relaxation dispersion measurements on base carbon and nitrogen nuclei have revealed that wobble G·T/U mismatches in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-abundance, and mutagenic Watson-Crick-like conformations. As Watson-Crick-like G·T mismatches have base pairing geometries similar to Watson-Crick base pairs, we hypothesized that they would mimic Watson-Crick base pairs with respect to the sugar-backbone conformation as well. Using off-resonance R_1ρ measurements targeting the sugar C3' and C4' nuclei, a structure survey, and molecular dynamics simulations, we show that wobble G·T mismatches adopt sugar-backbone conformations that deviate from the canonical Watson-Crick conformation and that transitions toward tautomeric and anionic Watson-Crick-like G·T mismatches restore the canonical Watson-Crick sugar-backbone. These measurements also reveal kinetic isotope effects for tautomerization in D₂O versus H₂O, which provide experimental evidence in support of a transition state involving proton transfer. The results provide additional evidence in support of mutagenic Watson-Crick-like G·T mismatches, help rule out alternative inverted wobble conformations in the case of anionic G·T^-, and also establish sugar carbons as new non-exchangeable probes of this exchange process.

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.

Pol μ dGTP mismatch insertion opposite T coupled with ligation reveals promutagenic DNA repair intermediate

Incorporation of mismatched nucleotides during DNA replication or repair leads to transition or transversion mutations and is considered as a predominant source of base substitution mutagenesis in cancer cells. Watson-Crick like dG:dT base pairing is considered to be an important source of genome instability. Here we show that DNA polymerase (pol) μ insertion of 7,8-dihydro-8′-oxo-dGTP (8-oxodGTP) or deoxyguanosine triphosphate (dGTP) into a model double-strand break DNA repair substrate with template base T results in efficient ligation by DNA ligase. These results indicate that pol μ-mediated dGTP mismatch insertion opposite template base T coupled with ligation could be a feature of mutation prone nonhomologous end joining during double-strand break repair.

Incorporation of mismatched nucleotides during DNA replication or repair can lead to mutagenesis. Here the authors reveal that DNA ligase can ligate NHEJ intermediates following incorporation of 8-oxodGTP or dGTP opposite T by DNA Polymerase mu (Pol mu) in vitro, which suggests that Pol mu could cause promutagenic mismatches during DSB repair.

The influence of Cu+ binding to hypoxanthine on stabilization of mismatches involving hypoxanthine and DNA bases: a DFT study

In the present work, the influence of Cu⁺ binding to N3- and N7-positions of hypoxanthine on energetic, geometrical and topological properties of hypoxanthine-guanine, hypoxanthine-adenine, hypoxanthine-cytosine, hypoxanthine-thymine and hypoxanthine-hypoxanthine mismatches is theoretically investigated. The calculations, in gas phase, are performed at B3LYP/6-311++G(3df,3pd) level of theory. Unlike the other mispairs, Cu⁺ binding to N3-position of hypoxanthine causes the proton transfer process from enol form of hypoxanthine to imino forms of adenine and cytosine. This process also occurs in all mismatches having enol form of hypoxanthine when Cu⁺ binds to N7-position of hypoxanthine. The mismatches are stabilized by hydrogen bonds. The influence of Cu⁺ on hydrogen bonds is also examined by atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Communicated by Ramaswamy H. Sarma.

The intrinsic basicity of the phosphate backbone exceeds that of uracil and thymine residues: protonation of the phosphate moiety is preferred over the nucleobase for pdThd and pUrd

The gas-phase conformations of the protonated forms of thymidine-5'-monophosphate and uridine-5'-monophosphate, [pdThd+H]⁺ and [pUrd+H]⁺, are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and electronic structure calculations. The IRMPD action spectra of [pdThd+H]⁺ and [pUrd+H]⁺ are measured over the IR fingerprint and hydrogen-stretching regions using the FELIX free electron laser and an OPO/OPA laser system. Low-energy conformations of [pdThd+H]⁺ and [pUrd+H]⁺ and their relative stabilities are computed at the MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparisons of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers indicate that the dominant conformers of [pdThd+H]⁺ and [pUrd+H]⁺ populated in the experiments are protonated at the phosphate oxo oxygen atom, with a syn nucleobase orientation that is stabilized by strong P[double bond, length as m-dash]OH⁺O2 and P-OHO4' hydrogen-bonding interactions, and C2'-endo sugar puckering. Minor abundance of conformers protonated at the O2 carbonyl of the nucleobase residue may also contribute for [pdThd+H]⁺, but do not appear to be important for [pUrd+H]⁺. Comparisons to previous IRMPD spectroscopy investigations of the protonated forms of thymidine and uridine, [dThd+H]⁺ and [Urd+H]⁺, and the deprotonated forms of pdThd and pUrd, [pdThd-H]^- and [pUrd-H]^-, provide insight into the effects of the phosphate moiety and protonation on the conformational features of the nucleobase and sugar moieties. Most interestingly, the thymine and uracil nucleobases remain in their canonical forms for [pdThd+H]⁺ and [pUrd+H]⁺, unlike [dThd+H]⁺ and [Urd+H]⁺, where protonation occurs on the nucleobases and induces tautomerization of the thymine and uracil residues.

Sources of spontaneous mutagenesis in bacteria

Mutations in an organism's genome can arise spontaneously, that is, in the absence of exogenous stress and prior to selection. Mutations are often neutral or deleterious to individual fitness but can also provide genetic diversity driving evolution. Mutagenesis in bacteria contributes to the already serious and growing problem of antibiotic resistance. However, the negative impacts of spontaneous mutagenesis on human health are not limited to bacterial antibiotic resistance. Spontaneous mutations also underlie tumorigenesis and evolution of drug resistance. To better understand the causes of genetic change and how they may be manipulated in order to curb antibiotic resistance or the development of cancer, we must acquire a mechanistic understanding of the major sources of mutagenesis. Bacterial systems are particularly well-suited to studying mutagenesis because of their fast growth rate and the panoply of available experimental tools, but efforts to understand mutagenic mechanisms can be complicated by the experimental system employed. Here, we review our current understanding of mutagenic mechanisms in bacteria and describe the methods used to study mutagenesis in bacterial systems.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

Direct NMR Evidence that Transient Tautomeric and Anionic States in dG·dT Form Watson-Crick-like Base Pairs

The replicative and translational machinery utilizes the unique geometry of canonical G·C and A·T/U Watson-Crick base pairs to discriminate against DNA and RNA mismatches in order to ensure high fidelity replication, transcription, and translation. There is growing evidence that spontaneous errors occur when mismatches adopt a Watson-Crick-like geometry through tautomerization and/or ionization of the bases. Studies employing NMR relaxation dispersion recently showed that wobble dG·dT and rG·rU mismatches in DNA and RNA duplexes transiently form tautomeric and anionic species with probabilities (≈0.01-0.40%) that are in concordance with replicative and translational errors. Although computational studies indicate that these exceptionally short-lived and low-abundance species form Watson-Crick-like base pairs, their conformation could not be directly deduced from the experimental data, and alternative pairing geometries could not be ruled out. Here, we report direct NMR evidence that the transient tautomeric and anionic species form hydrogen-bonded Watson-Crick-like base pairs. A guanine-to-inosine substitution, which selectively knocks out a Watson-Crick-type (G)N2H2···O2(T) hydrogen bond, significantly destabilized the transient tautomeric and anionic species, as assessed by lack of any detectable chemical exchange by imino nitrogen rotating frame spin relaxation (R1ρ) experiments. An 15N R1ρ NMR experiment targeting the amino nitrogen of guanine (dG-N2) provides direct evidence for Watson-Crick (G)N2H2···O2(T) hydrogen bonding in the transient tautomeric state. The strategy presented in this work can be generally applied to examine hydrogen-bonding patterns in nucleic acid transient states including in other tautomeric and anionic species that are postulated to play roles in replication and translational errors.

A QM/QTAIM research under the magnifying glass of the DPT tautomerisation of the wobble mispairs involving 2-aminopurine

In this study, a comprehensive survey of the changes of the physico-chemical parameters at each point of the IRC for the biologically important T·2AP*(w) ↔ T*·2AP(w) and G·2AP*(w) ↔ G*·2AP(w) DPT tautomerisation reactions involved in the point mutations (transitions and transversions) induced by 2-aminopurine (2AP) in DNA is provided.

Whether 2-aminopurine induces incorporation errors at the DNA replication? A quantum-mechanical answer on the actual biological issue

In this paper, we consider the mutagenic properties of the 2-aminopurine (2AP), which has intrigued molecular biologists, biophysicists and physical chemists for a long time and been widely studied by both experimentalists and theorists. We have shown for the first time using QM calculations, that 2AP very effectively produces incorporation errors binding with cytosine (C) into the wobble (w) C·2AP(w) mispair, which is supported by the N4H⋯N1 and N2H⋯N3 H-bonds and is tautomerized into the Watson-Crick (WC)-like base mispair C*·2AP(WC) (asterisk denotes the mutagenic tautomer of the base), that quite easily in the process of the thermal fluctuations acquires enzymatically competent conformation. 2AP less effectively produces transversions forming the wobble mispair with A base - A·2AP(w), stabilized by the participation of the N6H⋯N1 and N2H⋯N1 H-bonds, followed by further tautomerization A·2AP(w) → A*·2AP(WC) and subsequent conformational transition A*·2AP(WC) → A*·2AP_syn thus acquiring enzymatically competent structure. In this case, incorporation errors occur only in those case, when 2AP belongs to the incoming nucleotide. Thus, answering the question posed in the title of the article, we affirm for certain that 2AP induces incorporation errors at the DNA replication. Obtained results are consistent well with numerous experimental data.

Different Divalent Cations Alter the Kinetics and Fidelity of DNA Polymerases

Divalent metal ions are essential components of DNA polymerases both for catalysis of the nucleotidyl transfer reaction and for base excision. They occupy two sites, A and B, for DNA synthesis. Recently, a third metal ion was shown to be essential for phosphoryl transfer reaction. The metal ion in the A site is coordinated by the carboxylate of two highly conserved acidic residues, water molecules, and the 3'-hydroxyl group of the primer so that the A metal is in an octahedral complex. Its catalytic function is to lower the pK_a of the hydroxyl group, making it a highly effective nucleophile that can attack the α phosphorous atom of the incoming dNTP. The metal ion in the B site is coordinated by the same two carboxylates that are affixed to the A metal ion as well as the non-bridging oxygen atoms of the incoming dNTP. The carboxyl oxygen of an adjacent peptide bond serves as the sixth ligand that completes the octahedral coordination geometry of the B metal ion. Similarly, two metal ions are required for proofreading; one helps to lower the pK_a of the attacking water molecule, and the other helps to stabilize the transition state for nucleotide excision. The role of different divalent cations are discussed in relation to these two activities as well as their influence on base selectivity and misincorporation by DNA polymerases. Some, but not all, of the effects of these different metal ions can be rationalized based on their intrinsic properties, which are tabulated in this review.

In silico phylogenomics using complete genomes: a case study on the evolution of hominoids

The increasing availability of complete genome data is facilitating the acquisition of phylogenomic data sets, but the process of obtaining orthologous sequences from other genomes and assembling multiple sequence alignments remains piecemeal and arduous. We designed software that performs these tasks and outputs anonymous loci (AL) or anchored enrichment/ultraconserved element loci (AE/UCE) data sets in ready-to-analyze formats. We demonstrate our program by applying it to the hominoids. Starting with human, chimpanzee, gorilla, and orangutan genomes, our software generated an exhaustive data set of 292 ALs (∼1 kb each) in ∼3 h. Not only did analyses of our AL data set validate the program by yielding a portrait of hominoid evolution in agreement with previous studies, but the accuracy and precision of our estimated ancestral effective population sizes and speciation times represent improvements. We also used our program with a published set of 512 vertebrate-wide AE "probe" sequences to generate data sets consisting of 171 and 242 independent loci (∼1 kb each) in 11 and 13 min, respectively. The former data set consisted of flanking sequences 500 bp from adjacent AEs, while the latter contained sequences bordering AEs. Although our AE data sets produced the expected hominoid species tree, coalescent-based estimates of ancestral population sizes and speciation times based on these data were considerably lower than estimates from our AL data set and previous studies. Accordingly, we suggest that loci subjected to direct or indirect selection may not be appropriate for coalescent-based methods. Complete in silico approaches, combined with the burgeoning genome databases, will accelerate the pace of phylogenomics.

Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)-quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute-solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the 'bare' finite solute-solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in 'bulk' solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar - 11.74 kcal mol^-1 for guanosine and 11.16 kcal mol^-1 for acyclovir, and the respective rate constants (k = 1.5 × 10¹ s^-1, guanosine and k = 4.09 × 10¹ s^-1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.

A novel conception for spontaneous transversions caused by homo-pyrimidine DNA mismatches: a QM/QTAIM highlight

We have firstly shown that the T·T(w) and C·C(w) DNA mismatches with wobble (w) geometry stay in slow tautomeric equilibrium with short T·T*(WC) and C·C*(WC) Watson-Crick (WC) mispairs. These non-dissociative tautomeric rearrangements are controlled by the plane-symmetric, highly stable, highly polar and zwitterionic transition states. The obtained results allow us to understand in what way the T·T(w) and C·C(w) mismatches acquire enzymatically competent T·T*(WC) and C·C*(WC) conformations directly in the hydrophobic recognition pocket of a high-fidelity DNA-polymerase, thereby producing thermodynamically non-equilibrium spontaneous transversions. The simplest numerical estimation of the frequency ratio of the TT to CC spontaneous transversions satisfactorily agrees with experimental data.

2,4-Dihydroxy and O2 Protonated Tautomers of dThd and Thd Coexist in the Gas Phase: Methylation Alters Protonation Preferences versus dUrd and Urd

The gas-phase structures of protonated thymidine, [dThd + H](+), and its modified form, protonated 5-methyluridine, [Thd + H](+), are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy combined with electronic structure calculations. IRMPD action spectra are measured over the ranges extending from ~600 to 1900 cm(-1) and ~2800 to 3800 cm(-1) using the FELIX free electron laser and an optical parametric oscillator/amplifier (OPO/OPA) laser system, respectively. Comparisons between the B3LYP/6-311+G(d,p) linear IR spectra calculated for the stable low-energy conformers and the measured IRMPD spectra are used to determine the most favorable tautomeric conformations of [dThd + H](+) and [Thd + H](+) and to identify those populated in the experiments. Both B3LYP and MP2 levels of theory predict a minor 2,4-dihydroxy tautomer as the ground-state conformer of [dThd + H](+) and [Thd + H](+) indicating that the 2'-hydroxyl substituent of Thd does not exert a significant impact on the structural features. [dThd + H](+) and [Thd + H](+) share parallel IRMPD spectral profiles and yields in both the FELIX and OPO regions. Comparisons between the measured IRMPD and calculated IR spectra suggest that minor 2,4-dihydroxy tautomers and O2 protonated conformers of [dThd + H](+) and [Thd + H](+) are populated in the experiments. Comparison of this work to our previous IRMPD spectroscopy study of protonated 2'-deoxyuridine and uridine suggests that the 5-methyl substituent alters the preferences of O2 versus O4 protonation.

Polymerase specific error rates and profiles identified by single molecule sequencing

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated adenine nucleosides: N3 protonation induces base rotation and enhances N-glycosidic bond stability

N3 protonation induces base rotation and stabilizes the syn orientation of the adenine nucleobase of [dAdo+H]⁺ and [Ado+H]⁺via formation of a strong intramolecular N3H⁺⋯O5′ hydrogen-bonding interaction, which in turn influences the mechanisms and energetics for N-glycosidic bond cleavage.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).

Proton tunneling in the A∙T Watson-Crick DNA base pair: myth or reality?

The results and conclusions reached by Godbeer et al. in their recent work, that proton tunneling in the A∙T(WC) Watson-Crick (WC) DNA base pair occurs according to the Löwdin's (L) model, but with a small (~10(-9)) probability were critically analyzed. Here, it was shown that this finding overestimates the possibility of the proton tunneling at the A∙T(WC)↔A*∙T*(L) tautomerization, because this process cannot be implemented as a chemical reaction. Furthermore, it was outlined those biologically important nucleobase mispairs (A∙A*↔A*∙A, G∙G*↔G*∙G, T∙T*↔T*∙T, C∙C*↔C*∙C, H∙H*↔H*∙H (H - hypoxanthine)) - the players in the field of the spontaneous point mutagenesis - where the tunneling of protons is expected and for which the application of the model proposed by Godbeer et al. can be productive.

The effects of tautomerization and protonation on the adenine-cytosine mismatches: a density functional theory study

In the present work, we demonstrate the results of a theoretical study concerned with the question how tautomerization and protonation of adenine affect the various properties of adenine-cytosine mismatches. The calculations, in gas phase and in water, are performed at B3LYP/6-311++G(d,p) level. In gas phase, it is observed that any tautomeric form of investigated mismatches is more stabilized when adenine is protonated. As for the neutral mismatches, the mismatches containing amino form of cytosine and imino form of protonated adenine are more stable. The role of aromaticity on the stability of tautomeric forms of mismatches is investigated by NICS(1)ZZ index. The stability of mispairs decreases by going from gas phase to water. It can be explained using dipole moment parameter. The influence of hydrogen bonds on the stability of mismatches is examined by atoms in molecules and natural bond orbital analyses. In addition to geometrical parameters and binding energies, the study of the topological properties of electron charge density aids in better understanding of these mispairs.

Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes

Rare tautomeric and anionic nucleobases are believed to have fundamental biological roles, but their prevalence and functional importance has remained elusive because they exist transiently, in low abundance, and involve subtle movements of protons that are difficult to visualize. Using NMR relaxation dispersion, we show here that wobble dG•dT and rG•rU mispairs in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-populated Watson-Crick-like mispairs that are stabilized by rare enolic or anionic bases. These mispairs can evade Watson-Crick fidelity checkpoints and form with probabilities (10(-3) to 10(-5)) that strongly imply a universal role in replication and translation errors. Our results indicate that rare tautomeric and anionic bases are widespread in nucleic acids, expanding their structural and functional complexity beyond that attainable with canonical bases.

Mispairs with Watson-Crick base-pair geometry observed in ternary complexes of an RB69 DNA polymerase variant

Recent structures of DNA polymerase complexes with dGMPCPP/dT and dCTP/dA mispairs at the insertion site have shown that they adopt Watson-Crick geometry in the presence of Mn(2+) indicating that the tautomeric or ionization state of the base has changed. To see whether the tautomeric or ionization state of base-pair could be affected by its microenvironment, we determined 10 structures of an RB69 DNA polymerase quadruple mutant with dG/dT or dT/dG mispairs at position n-1 to n-5 of the Primer/Template duplex. Different shapes of the mispairs, including Watson-Crick geometry, have been observed, strongly suggesting that the local environment of base-pairs plays an important role in their tautomeric or ionization states.

New structural hypostases of the A·T and G·C Watson–Crick DNA base pairs caused by their mutagenic tautomerisation in a wobble manner: a QM/QTAIM prediction

Our investigation reveals the hitherto unknown ability of the canonical Watson–Crick DNA base pairs to switch into wobble mismatches with mutagenic tautomers, clarifying the nature of genome instability.

O6-Carboxymethylguanine in DNA forms a sequence context-dependent wobble base-pair structure with thymine

N-Nitrosation of glycine and its derivatives generates potent alkylating agents that can lead to the formation ofO6-carboxymethylguanine (O6-CMG) in DNA.O6-CMG has been identified in DNA derived from human colon tissue and its occurrence has been linked to diets high in red and processed meats, implying an association with the induction of colorectal cancer. By analogy toO6-methylguanine,O6-CMG is expected to be mutagenic, inducing G-to-A mutations that may be the molecular basis of increased cancer risk. Previously, the crystal structure of the DNA dodecamer d(CGCG[O6-CMG]ATTCGCG) has been reported, in whichO6-CMG forms a Watson–Crick-type pair with thymine similar to the canonical A:T pair. In order to further investigate the versatility ofO6-CMG in base-pair formation, the structure of the DNA dodecamer d(CGC[O6-CMG]AATTTGCG) containingO6-CMG at a different position has been determined by X-ray crystallography using four crystal forms obtained under conditions containing different solvent ions (Sr2+, Ba2+, Mg2+, K+or Na+) with and without Hoechst 33258. The most striking finding is that the pairing modes ofO6-CMG with T are quite different from those previously reported. In the present dodecamer, the T bases are displaced (wobbled) into the major groove to form a hydrogen bond between the thymine N3N—H and the carboxyl group ofO6-CMG. In addition, a water molecule is bridged through two hydrogen bonds between the thymine O2atom and the 2-amino group ofO6-CMG to stabilize the pairing. These interaction modes commonly occur in the four crystal forms, regardless of the differences in crystallization conditions. The previous and the present results show thatO6-CMG can form a base pair with T in two alternative modes: the Watson–Crick type and a high-wobble type, the nature of which may depend on the DNA-sequence context.

Recognition of Watson-Crick base pairs: constraints and limits due to geometric selection and tautomerism

The natural bases of nucleic acids have a strong preference for one tautomer form, guaranteeing fidelity in their hydrogen bonding potential. However, base pairs observed in recent crystal structures of polymerases and ribosomes are best explained by an alternative base tautomer, leading to the formation of base pairs with Watson-Crick-like geometries. These observations set limits to geometric selection in molecular recognition of complementary Watson-Crick pairs for fidelity in replication and translation processes.

DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments

We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

Theoretical study of the pre- and post-translational effects of adenine and thymine tautomers and methyl derivatives

The study of pre-translational effects (ionization, tautomerization) and post-translational effects (methylation) of adenine and thymine has only recently been the focus of some studies. These effects can potentially help regulate gene expression as well as potentially disrupt normal gene function. Because of this wide array of roles, greater insight into these effects in deoxyribonucleic acids (DNA) are paramount. There has been considerable research of each phenomenon (tautomerization, methylation and ionization) individually. In this work, we attempt to shed light upon the pre-translational effects and post translational effects of adenine and thymine by investigating the electron affinities (EAs) and ionization potentials (IPs) of the major and minor tautomers and their methyl derivatives. We performed all calculations using the density functional theory (DFT) B3LYP functional accompanied with 6-311G(d,p), 6-311+G(d,p) and 6-311++G(df,pd) basis sets. Our results reveal that the thymine tautomer has a higher EA and IP than the adenine tautomers. The higher EA suggests that an electron that attaches to the AT base pair would predominately attach to the thymine instead of adenine. The higher IP would suggest that an electron that is removed from the AT base pair would be predominately removed from the adenine within the base pair. Understanding how tautomerization, ionization and methylation differences change effects, discourages, or promotes one another is lacking. In this work, we begin the steps of integrating these effects with one another, to gain a greater understanding of molecular changes in DNA bases.

New structural insights into the decoding mechanism: translation infidelity via a G·U pair with Watson-Crick geometry

Pioneer crystallographic studies of the isolated 30S ribosomal subunit provided the first structural insights into the decoding process. Recently, new crystallographic data on full 70S ribosomes with mRNA and tRNAs have shown that the formation of the tight decoding centre is ensured by conformational rearrangement of the 30S subunit (domain closure), which is identical for cognate or near-cognate tRNA. When a G·U forms at the first or second codon-anticodon positions (near-cognate tRNA), the ribosomal decoding centre forces the adoption of Watson-Crick G·C-like geometry rather than that of the expected Watson-Crick wobble pair. Energy expenditure for rarely occuring tautomeric base required for Watson-Crick G·C-like G·U pair or the repulsion energy due to steric clash within the mismatched base pair could constitute the only cause for efficient rejection of a near-cognate tRNA. Our data suggest that "geometrical mimicry" can explain how wrong aminoacyl-tRNAs with G·U pairs in the codon-anticodon helix forming base pairs with Watson-Crick geometry in the decoding center can be incorporated into the polypeptide chain.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β(μ)) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree-Fock, Møller-Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning's correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β(μ) (9H)/β(μ) (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β(μ) values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.