Tautomerism provides a molecular explanation for the mutagenic properties of the anti-HIV nucleoside 5-aza-5,6-dihydro-2'-deoxycytidine

Conformation-dependent lesion bypass of bulky arylamine-dG adducts generated from 2-nitrofluorene in epigenetic sequence contexts

Sequence context influences structural characteristics and repair of DNA adducts, but there is limited information on how epigenetic modulation affects conformational heterogeneity and bypass of DNA lesions. Lesions derived from the environmental pollutant 2-nitrofluorene have been extensively studied as chemical carcinogenesis models; they adopt a sequence-dependent mix of two significant conformers: major groove binding (B) and base-displaced stacked (S). We report a conformation-dependent bypass of the N-(2'-deoxyguanosin-8-yl)-7-fluoro-2-aminofluorene (dG-FAF) lesion in epigenetic sequence contexts (d[5'-CTTCTC#G*NCCTCATTC-3'], where C# is C or 5-methylcytosine (5mC), G* is G or G-FAF, and N is A, T, C or G). FAF-modified sequences with a 3' flanking pyrimidine were better bypassed when the 5' base was 5mC, whereas sequences with a 3' purine exhibited the opposite effect. The conformational basis behind these variations differed; for -CG*C- and -CG*T-, bypass appeared to be inversely correlated with population of the duplex-destabilizing S conformer. On the other hand, the connection between conformation and a decrease in bypass for flanking purines in the 5mC sequences relative to C was more complex. It could be related to the emergence of a disruptive non-S/B conformation. The present work provides novel conformational insight into how 5mC influences the bypass efficiency of bulky DNA damage.

Tautomerism and Rotamerism of Favipiravir and Halogenated Analogues in Solution and in the Solid State

Favipiravir is an important selective antiviral against RNA-based viruses, and currently, it is being repurposed as a potential drug for the treatment of COVID-19. This type of chemical system presents different carboxamide-rotameric and hydroxyl-tautomeric states, which could be essential for interpreting its selective antiviral activity. Herein, the tautomeric 3-hydroxypyrazine/3-pyrazinone pair of favipiravir and its 6-substituted analogues, 6-Cl, 6-Br, 6-I, and 6-H, were fully investigated in solution and in the solid state through ultraviolet-visible, 1H nuclear magnetic resonance, infrared spectroscopy, and X-ray diffraction techniques. Also, a study of the gas phase was performed using density functional theory calculations. In general, the keto-enol balance in these 3-hydroxy-2-pyrazinecarboxamides is finely modulated by external and internal electrical variations via changes in solvent polarity or by replacement of substituents at position 6. The enol tautomer was prevalent in an apolar environment, whereas an increase in the level of the keto tautomer was favored by an increase in solvent polarity and, even moreso, with a strong hydrogen-donor solvent. Keto tautomerization was favored either in solution or in the solid state with a decrease in 6-substituent electronegativity as follows: H ≫ I ≈ Br > Cl ≥ F. Specific rotameric states based on carboxamide, "cisoide" and "transoide", were identified for the enol and keto tautomer, respectively; their rotamerism is dependent on the tautomerism and not the aggregation state.

Voltage-driven control of single-molecule keto-enol equilibrium in a two-terminal junction system

Keto-enol tautomerism, describing an equilibrium involving two tautomers with distinctive structures, provides a promising platform for modulating nanoscale charge transport. However, such equilibria are generally dominated by the keto form, while a high isomerization barrier limits the transformation to the enol form, suggesting a considerable challenge to control the tautomerism. Here, we achieve single-molecule control of a keto-enol equilibrium at room temperature by using a strategy that combines redox control and electric field modulation. Based on the control of charge injection in the single-molecule junction, we could access charged potential energy surfaces with opposite thermodynamic driving forces, i.e., exhibiting a preference for the conducting enol form, while the isomerization barrier is also significantly reduced. Thus, we could selectively obtain desired and stable tautomers, which leads to significant modulation of the single-molecule conductance. This work highlights the concept of single-molecule control of chemical reactions on more than one potential energy surface.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Hyperpolarizing DNA Nucleobases via NMR Signal Amplification by Reversible Exchange

The present work investigates the potential for enhancing the NMR signals of DNA nucleobases by parahydrogen-based hyperpolarization. Signal amplification by reversible exchange (SABRE) and SABRE in Shield Enables Alignment Transfer to Heteronuclei (SABRE-SHEATH) of selected DNA nucleobases is demonstrated with the enhancement (ε) of ¹H, ¹⁵N, and/or ¹³C spins in 3-methyladenine, cytosine, and 6-O-guanine. Solutions of the standard SABRE homogenous catalyst Ir(1,5-cyclooctadeine)(1,3-bis(2,4,6-trimethylphenyl)imidazolium)Cl ("IrIMes") and a given nucleobase in deuterated ethanol/water solutions yielded low ¹H ε values (≤10), likely reflecting weak catalyst binding. However, we achieved natural-abundance enhancement of ¹⁵N signals for 3-methyladenine of ~3300 and ~1900 for the imidazole ring nitrogen atoms. ¹H and ¹⁵N 3-methyladenine studies revealed that methylation of adenine affords preferential binding of the imidazole ring over the pyrimidine ring. Interestingly, signal enhancements (ε~240) of both ¹⁵N atoms for doubly labelled cytosine reveal the preferential binding of specific tautomer(s), thus giving insight into the matching of polarization-transfer and tautomerization time scales. ¹³C enhancements of up to nearly 50-fold were also obtained for this cytosine isotopomer. These efforts may enable the future investigation of processes underlying cellular function and/or dysfunction, including how DNA nucleobase tautomerization influences mismatching in base-pairing.

Unravelling the molecular mechanism of mutagenic factors impacting human health

Environmental mutagens are chemical and physical substances in the environment that has a potential to induce a wide range of mutations and generate multiple physiological, biochemical, and genetic modifications in humans. Most mutagens are having genotoxic effects on the following generation through germ cells. The influence of germinal mutations on health will be determined by their frequency, nature, and the mechanisms that keep a specific mutation in the population. Early prenatal lethal mutations have less public health consequences than genetic illnesses linked with long-term medical and social difficulties. Physical and chemical mutagens are common mutagens found in the environment. These two environmental mutagens have been associated with multiple neurological disorders and carcinogenesis in humans. Thus in this study, we aim to unravel the molecular mechanism of physical mutagens (UV rays, X-rays, gamma rays), chemical mutagens (dimethyl sulfate (DMS), bisphenol A (BPA), polycyclic aromatic hydrocarbons (PAHs), 5-chlorocytosine (5ClC)), and several heavy metals (Ar, Pb, Al, Hg, Cd, Cr) implicated in DNA damage, carcinogenesis, chromosomal abnormalities, and oxidative stress which leads to multiple disorders and impacting human health. Biological tests for mutagen detection are crucial; therefore, we also discuss several approaches (Ames test and Mutatox test) to estimate mutagenic factors in the environment. The potential risks of environmental mutagens impacting humans require a deeper basic knowledge of human genetics as well as ongoing research on humans, animals, and their tissues and fluids.

Lethal Mutagenesis of RNA Viruses and Approved Drugs with Antiviral Mutagenic Activity

In RNA viruses, a small increase in their mutation rates can be sufficient to exceed their threshold of viability. Lethal mutagenesis is a therapeutic strategy based on the use of mutagens, driving viral populations to extinction. Extinction catastrophe can be experimentally induced by promutagenic nucleosides in cell culture models. The loss of HIV infectivity has been observed after passage in 5-hydroxydeoxycytidine or 5,6-dihydro-5-aza-2′-deoxycytidine while producing a two-fold increase in the viral mutation frequency. Among approved nucleoside analogs, experiments with polioviruses and other RNA viruses suggested that ribavirin can be mutagenic, although its mechanism of action is not clear. Favipiravir and molnupiravir exert an antiviral effect through lethal mutagenesis. Both drugs are broad-spectrum antiviral agents active against RNA viruses. Favipiravir incorporates into viral RNA, affecting the G→A and C→U transition rates. Molnupiravir (a prodrug of β-d-N⁴-hydroxycytidine) has been recently approved for the treatment of SARS-CoV-2 infection. Its triphosphate derivative can be incorporated into viral RNA and extended by the coronavirus RNA polymerase. Incorrect base pairing and inefficient extension by the polymerase promote mutagenesis by increasing the G→A and C→U transition frequencies. Despite having remarkable antiviral action and resilience to drug resistance, carcinogenic risks and genotoxicity are important concerns limiting their extended use in antiviral therapy.

Tautomerization and Isomerization in Quantitative NMR: A Case Study with 4-Deoxynivalenol (DON)

The regulated mycotoxin 4-deoxynivalenol (DON) has a heterocyclic structure that is readily amenable to tautomerization and conformational isomerization in solution. An analysis of DON in solution by NMR revealed the presence of hemiacetal tautomer(s) and putative conformational isomers, which maintain the intact enone functional group. The extent and type of tautomerization/isomerization vary according to the NMR solvent used and produce different signal patterns in the NMR spectra. Thus, the same proton produces multiple signals depending on which isomer/tautomer it belongs to. To maintain the accuracy of quantitative NMR (qNMR) measurements, it was essential to conclusively identify all signals belonging to the same proton to avoid underestimating its integral value. A strategy to overcome the complications of DON tautomerization and isomerization in solution during qNMR is reported. Of all proton atoms on the DON carbo-skeleton, H-10 produced clearly defined signals centered at 6.6 ppm for suspected conformational isomers and at 5.5 ppm for hemiacetal tautomers. To determine the purity of DON by quantitative proton NMR, the collective integrals of all isomeric and tautomeric signals belonging to H-10 provided the most accurate value. The purity of DON obtained with this protocol is highly accurate and suitable for the value assignment of certified reference materials (CRMs).

Surface-Induced Keto-Enol Tautomerization of DNA Base Molecules and Consequent [4 + 2]-like Cycloaddition on Si(111)7×7

Adsorption and film growth of deoxyribonucleic acid (DNA) base molecules (cytosine, guanine, thymine, and adenine) on Si(111)7×7 have been studied by combining X-ray photoelectron spectroscopy (XPS) with ab initio calculations based on the density functional theory (DFT). Multiple tautomeric forms and keto-enol tautomerization are revealed by the O 1s, N 1s, and C 1s XPS spectra of the O-containing DNA bases: cytosine, guanine, and thymine. While the carbonyl group-containing keto tautomer is more stable in a thick film and in powder, the hydroxyl group-containing enol tautomer is found at the interface. The keto-enol tautomerization, as induced by the reactive Si(111)7×7 surface, leads to the formation of a conjugated aromatic six-membered ring with a delocalized π electron system and to the consequent [4 + 2]-like cycloaddition between the enol tautomer and the 7×7 surface. The DFT calculation suggests that the enol tautomer exhibits a kinetic advantage over the keto one for the [4 + 2]-like cycloaddition. Among the several plausible pathways for the cycloaddition provided by the enol tautomer, the experimentally determined one involves a ring N and ring C atom (a polar pair), rather than two ring C atoms (a nonpolar pair), to better match the polar Si adatom-restatom pair of the 7×7 surface. Furthermore, the reacted ring C atom does not have any attached terminal functional group (e.g., -NH₂ and -OH). Further deposition leads to continuous film growth in the keto tautomeric form for cytosine and guanine. For the only O-free DNA base molecule, adenine, dative bonding N → Si, rather than the [4 + 2]-like cycloaddition, is observed on the 7×7 surface. Of the four DNA base molecules, adenine is also the only one with its aromaticity maintained when adsorbed on the Si(111)7×7 surface. A reactive surface like the 7×7 surface could therefore provide a new control to trigger tautomerization that is often associated with genetic mutation.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

OUP accepted manuscript

BlueScreen HC is a mammalian cell-based assay for measuring the genotoxicity and cytotoxicity of chemical compounds and mixtures. The BlueScreen HC assay has been utilized at the Research Institute for Fragrance Materials in a safety assessment program as a screening tool to prioritize fragrance materials for higher-tier testing, as supporting evidence when using a read-across approach, and as evidence to adjust the threshold of toxicological concern. Predictive values for the BlueScreen HC assay were evaluated based on the ability of the assay to predict the outcome of in vitro and in vivo mutagenicity and chromosomal damage genotoxicity assays. A set of 371 fragrance materials was assessed in the BlueScreen HC assay along with existing or newly generated in vitro and in vivo genotoxicity data. Based on a weight-of-evidence approach, the majority of materials in the data set were deemed negative and concluded not to have the potential to be genotoxic, while only a small proportion of materials were determined to show genotoxic effects in these assays. Analysis of the data set showed a combination of high positive agreement but low negative agreement between BlueScreen HC results, in vitro regulatory genotoxicity assays, and higher-tier test results. The BlueScreen HC assay did not generate any false negatives, thereby providing robustness when utilizing it as a high-throughput screening tool to evaluate the large inventory of fragrance materials. From the perspective of protecting public health, it is desirable to have no or minimal false negatives, as a false-negative result may incorrectly indicate the lack of a genotoxicity hazard. However, the assay did have a high percentage of false-positive results, resulting in poor positive predictivity of the in vitro genotoxicity test battery outcome. Overall, the assay generated 100% negative predictivity and 3.9% positive predictivity. In addition to the data set of 371 fragrance materials, 30 natural complex substances were evaluated for BlueScreen HC, Ames, and in vitro micronucleus assay, and a good correlation in all three assays was observed. Overall, while a positive result may have to be further investigated, these findings suggest that the BlueScreen HC assay can be a valuable screening tool to detect the genotoxic potential of fragrance materials and mixtures.

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.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

How Hydrogen Bonding Amplifies Isomeric Differences in Pyridones toward Strong Changes in Acidity and Tautomerism

Steric hindrance of hydration and hydrogen bond enhancement by localized charges have been identified as key factors for the massive chemical differences between the hydroxypyridine/pyridone isomers in aqueous solution. While all isomers occur mainly in the hydroxypyridine form in the gas phase, they differ by more than 3 orders of magnitude both in their acidity and tautomeric equilibrium constants upon hydration. By monitoring the electronic and solvation structures as a function of the protonation state and the O^– substitution position on the pyridine ring, the amplification of the isomeric differences in aqueous solution has been investigated. Near-edge X-ray absorption fine structure (NEXAFS) measurements at the N K-edge served as the probe of the chemical state. The combination of molecular dynamics simulations, complete active space self-consistent field (CASSCF), and time-dependent density functional theory (TD-DFT) spectral calculations contributes to unraveling the principles of tautomerism and acidity in multiple biochemical systems based on tautomerism.

Unnatural bases for recognition of noncoding nucleic acid interfaces

The notion of using synthetic heterocycles instead of the native bases to interface with DNA and RNA has been explored for nearly 60 years. Unnatural bases compatible with the DNA/RNA coding interface have the potential to expand the genetic code and co-opt the machinery of biology to access new macromolecular function; accordingly, this body of research is core to synthetic biology. While much of the literature on artificial bases focuses on code expansion, there is a significant and growing effort on docking synthetic heterocycles to noncoding nucleic acid interfaces; this approach seeks to illuminate major processes of nucleic acids, including regulation of transcription, translation, transport, and transcript lifetimes. These major avenues of research at the coding and noncoding interfaces have in common fundamental principles in molecular recognition. Herein, we provide an overview of foundational literature in biophysics of base recognition and unnatural bases in coding to provide context for the developing area of targeting noncoding nucleic acid interfaces with synthetic bases, with a focus on systems developed through iterative design and biophysical study.

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.

Tautomer Standardization in Chemical Databases: Deriving Business Rules from Quantum Chemistry

Databases of small, potentially bioactive molecules are ubiquitous across the industry and academia. Designed such that each unique compound should appear only once, the multiplicity of ways in which many compounds can be represented means that these databases require methods for standardizing the representation of chemistry. This is commonly achieved through the use of "Chemistry Business Rules", sets of predefined rules that describe the "house style" of the database in question. At Syngenta, the historical approach to the design of chemistry business rules has been to focus on consistency of representation, with chemical relevance given secondary consideration. In this work, we overturn that convention. Through the use of quantum chemistry calculations, we define a set of chemistry business rules for tautomer standardization that reproduces gas-phase energetic preferences. We go on to show that, compared to our historic approach, this method yields tautomers that are in better agreement with those observed experimentally in condensed phases and that are better suited for use in predictive models.

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.

A dielectric and spectrophotometric study of the tautomerization of 2-hydroxypyridine and 2-mercaptopyridine in water

The basic ionization (pk₁) and acidic ionization (pk₂) constants and equilibrium constant (K_T) of 2HPy and 2MPy were determined. The pk₁(s) of their N- and X-methyl derivatives (X = O, S) were also determined. The equilibrium constant of 2MPy is approximately 60 times greater than its oxygen analog, 2HPy. The micro-ionization constants of the functional groups, –NH (pk_A and pk_C) and –XH (pk_B and pk_D), were determined to provide further insights into the ionization equilibria of these N-heteroaromatic XH compounds (2HPy and 2MPy). The relaxation time of water (τ) in aqueous solutions of 2HPy and 2MPy are collectively used with the K_T values to determine the forward (k_f) and backward (k_b) rate constants of tautomerization. Subsequently, the k_f and k_b are used to provide the rationale for the K_T and τ values.

The basic ionization (pk₁) and acidic ionization (pk₂) constants and equilibrium constant (K_T) of 2HPy and 2MPy were determined.

Role of different tautomers in the base-pairing abilities of some of the vital antiviral drugs used against COVID-19

Base-pair mutations induced by different tautomers of anti-viral drugs are the main reasons for their anti-viral activities.

Biological Evaluation of DNA Biomarkers in a Chemically Defined and Site-Specific Manner

As described elsewhere in this Special Issue on biomarkers, much progress has been made in the detection of modified DNA within organisms at endogenous and exogenous levels of exposure to chemical species, including putative carcinogens and chemotherapeutic agents. Advances in the detection of damaged or unnatural bases have been able to provide correlations to support or refute hypotheses between the level of exposure to oxidative, alkylative, and other stresses, and the resulting DNA damage (lesion formation). However, such stresses can form a plethora of modified nucleobases, and it is therefore difficult to determine the individual contribution of a particular modification to alter a cell's genetic fate, as measured in the form of toxicity by stalled replication past the damage, by subsequent mutation, and by lesion repair. Chemical incorporation of a modification at a specific site within a vector (site-specific mutagenesis) has been a useful tool to deconvolute what types of damage quantified in biologically relevant systems may lead to toxicity and/or mutagenicity, thereby allowing researchers to focus on the most relevant biomarkers that may impact human health. Here, we will review a sampling of the DNA modifications that have been studied by shuttle vector techniques.

An unprecedented intramolecular to intermolecular mechanistic switch in 1,1-diaminoazines leading to differential product formation during the I2-induced tandem oxidative transformation

Intermolecular cyclization of 1,1-diaminoazine.

Inhibition of the prototropic tautomerism in chrysazine by p-sulfonatocalixarene hosts

This study reveals the unusual inhibition of excited-state prototropic tautomerism of Chrysazine by p-sulfonatocalix[4,6]arene hosts.

New strategies against drug resistance to herpes simplex virus

Herpes simplex virus (HSV), a member of the Herpesviridae family, is a significant human pathogen that results in mucocutaneous lesions in the oral cavity or genital infections. Acyclovir (ACV) and related nucleoside analogues can successfully treat HSV infections, but the emergence of drug resistance to ACV has created a barrier for the treatment of HSV infections, especially in immunocompromised patients. There is an urgent need to explore new and effective tactics to circumvent drug resistance to HSV. This review summarises the current strategies in the development of new targets (the DNA helicase/primase (H/P) complex), new types of molecules (nature products) and new antiviral mechanisms (lethal mutagenesis of Janus-type nucleosides) to fight the drug resistance of HSV.

Scaling law characterizing the dynamics of the transition of HIV-1 to error catastrophe

Increasing the mutation rate, μ of viruses above a threshold, μ(c) has been predicted to trigger a catastrophic loss of viral genetic information and is being explored as a novel intervention strategy. Here, we examine the dynamics of this transition using stochastic simulations mimicking within-host HIV-1 evolution. We find a scaling law governing the characteristic time of the transition: τ ≈ 0.6/(μ - μ(c)). The law is robust to variations in underlying evolutionary forces and presents guidelines for treatment of HIV-1 infection with mutagens. We estimate that many years of treatment would be required before HIV-1 can suffer an error catastrophe.

Structural Characterization of Single-Stranded DNA Monolayers Using Two-Dimensional Sum Frequency Generation Spectroscopy

DNA-covered materials are important in technological applications such as biosensors and microarrays, but obtaining structural information on surface-bound biomolecules is experimentally challenging. In this paper, we structurally characterize single-stranded DNA monolayers of poly(thymine) from 10 to 25 bases in length with an emerging surface technique called two-dimensional sum frequency generation (2D SFG) spectroscopy. These experiments are carried out by adding a mid-IR pulse shaper to a femtosecond broad-band SFG spectrometer. Cross peaks and 2D line shapes in the 2D SFG spectra provide information about structure and dynamics. Because the 2D SFG spectra are heterodyne detected, the monolayer spectra can be directly compared to 2D infrared (2D IR) spectra of poly(thymine) in solution, which aids interpretation. We simulate the 2D SFG spectra using DFT calculations and an excitonic Hamiltonian that relates the molecular geometry to the vibrational coupling. Intrabase cross peaks help define the orientation of the bases and interbase cross peaks, created by coupling between bases, and resolves features not observed in 1D SFG spectra that constrain the relative geometries of stacked bases. We present a structure for the poly(T) oligomer that is consistent with the 2D SFG data. These experiments provide insight into the DNA monolayer structure and set precedent for studying complex biomolecules on surfaces with 2D SFG spectroscopy.

Intrinsic mutagenic properties of 5-chlorocytosine: A mechanistic connection between chronic inflammation and cancer

During chronic inflammation, neutrophil-secreted hypochlorous acid can damage nearby cells inducing the genomic accumulation of 5-chlorocytosine (5ClC), a known inflammation biomarker. Although 5ClC has been shown to promote epigenetic changes, it has been unknown heretofore if 5ClC directly perpetrates a mutagenic outcome within the cell. The present work shows that 5ClC is intrinsically mutagenic, both in vitro and, at a level of a single molecule per cell, in vivo. Using biochemical and genetic approaches, we have quantified the mutagenic and toxic properties of 5ClC, showing that this lesion caused C→T transitions at frequencies ranging from 3-9% depending on the polymerase traversing the lesion. X-ray crystallographic studies provided a molecular basis for the mutagenicity of 5ClC; a snapshot of human polymerase β replicating across a primed 5ClC-containing template uncovered 5ClC engaged in a nascent base pair with an incoming dATP analog. Accommodation of the chlorine substituent in the template major groove enabled a unique interaction between 5ClC and the incoming dATP, which would facilitate mutagenic lesion bypass. The type of mutation induced by 5ClC, the C→T transition, has been previously shown to occur in substantial amounts both in tissues under inflammatory stress and in the genomes of many inflammation-associated cancers. In fact, many sequence-specific mutational signatures uncovered in sequenced cancer genomes feature C→T mutations. Therefore, the mutagenic ability of 5ClC documented in the present study may constitute a direct functional link between chronic inflammation and the genetic changes that enable and promote malignant transformation.

How many tautomerization pathways connect Watson-Crick-like G*·T DNA base mispair and wobble mismatches?

In this study, we have theoretically demonstrated the intrinsic ability of the wobble G·T(w)/G*·T*(w)/G·T(w1)/G·T(w2) and Watson-Crick-like G*·T(WC) DNA base mispairs to interconvert into each other via the DPT tautomerization. We have established that among all these transitions, only one single G·T(w) ↔ G*·T(WC) pathway is eligible from a biological perspective. It involves short-lived intermediate - the G·T*(WC) base mispair - and is governed by the planar, highly stable, and zwitterionic [Formula: see text] transition state stabilized by the participation of the unique pattern of the five intermolecular O6(+)H⋯O4(-), O6(+)H⋯N3(-), N1(+)H⋯N3(-), N1(+)H⋯O2(-), and N2(+)H⋯O2(-) H-bonds. This non-dissociative G·T(w) ↔ G*·T(WC) tautomerization occurs without opening of the pair: Bases within mispair remain connected by 14 different patterns of the specific intermolecular interactions that successively change each other along the IRC. Novel kinetically controlled mechanism of the thermodynamically non-equilibrium spontaneous point GT/TG incorporation errors has been suggested. The mutagenic effect of the analogues of the nucleotide bases, in particular 5-bromouracil, can be attributed to the decreasing of the barrier of the acquisition by the wobble pair containing these compounds of the enzymatically competent Watson-Crick's geometry via the intrapair mutagenic tautomerization directly in the essentially hydrophobic recognition pocket of the replication DNA-polymerase machinery. Proposed approaches are able to explain experimental data, namely growth of the rate of the spontaneous point incorporation errors during DNA biosynthesis with increasing temperature.

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.

Two-dimensional IR spectroscopy of the anti-HIV agent KP1212 reveals protonated and neutral tautomers that influence pH-dependent mutagenicity

Antiviral drugs designed to accelerate viral mutation rates can drive a viral population to extinction in a process called lethal mutagenesis. One such molecule is 5,6-dihydro-5-aza-2'-deoxycytidine (KP1212), a selective mutagen that induces A-to-G and G-to-A mutations in the genome of replicating HIV. The mutagenic property of KP1212 was hypothesized to originate from its amino-imino tautomerism, which would explain its ability to base pair with either G or A. To test the multiple tautomer hypothesis, we used 2D IR spectroscopy, which offers subpicosecond time resolution and structural sensitivity to distinguish among rapidly interconverting tautomers. We identified several KP1212 tautomers and found that >60% of neutral KP1212 is present in the enol-imino form. The abundant proportion of this traditionally rare tautomer offers a compelling structure-based mechanism for pairing with adenine. Additionally, the pKa of KP1212 was measured to be 7.0, meaning a substantial population of KP1212 is protonated at physiological pH. Furthermore, the mutagenicity of KP1212 was found to increase dramatically at pH <7, suggesting a significant biological role for the protonated KP1212 molecules. Overall, our data reveal that the bimodal mutagenic properties of KP1212 result from its unique shape shifting ability that utilizes both tautomerization and protonation.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

Evaluation of anti-HIV-1 mutagenic nucleoside analogues

Because of their high mutation rates, RNA viruses and retroviruses replicate close to the threshold of viability. Their existence as quasi-species has pioneered the concept of "lethal mutagenesis" that prompted us to synthesize pyrimidine nucleoside analogues with antiviral activity in cell culture consistent with an accumulation of deleterious mutations in the HIV-1 genome. However, testing all potentially mutagenic compounds in cell-based assays is tedious and costly. Here, we describe two simple in vitro biophysical/biochemical assays that allow prediction of the mutagenic potential of deoxyribonucleoside analogues. The first assay compares the thermal stabilities of matched and mismatched base pairs in DNA duplexes containing or not the nucleoside analogues as follows. A promising candidate should display a small destabilization of the matched base pair compared with the natural nucleoside and the smallest gap possible between the stabilities of the matched and mismatched base pairs. From this assay, we predicted that two of our compounds, 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine, should be mutagenic. The second in vitro reverse transcription assay assesses DNA synthesis opposite nucleoside analogues inserted into a template strand and subsequent extension of the newly synthesized base pairs. Once again, only 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine are predicted to be efficient mutagens. The predictive potential of our fast and easy first line screens was confirmed by detailed analysis of the mutation spectrum induced by the compounds in cell culture because only compounds 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine were found to increase the mutation frequency by 3.1- and 3.4-fold, respectively.