Tautomerism of the DNA Base Guanine and Its Methylated Derivatives as Studied by Gas-Phase Infrared and Ultraviolet Spectroscopy

Experimental and theoretical study of the cytosine tautomerism through excited states

CONTEXT

The irradiation of water solution of cytosine with UV light (λmax = 254 nm) shows oxo-hydroxy tautomerism with a rate constant of 6.297 × 10-3 min-1. The order of the reaction implies a tautomeric conversion. After removing the UV light source, we observed a dark reaction with a rate constant of 1.473 × 10-3 min-1 which leads to a restoration of the initial tautomer as before the irradiation. The mechanism of oxo-hydroxy tautomerism of cytosine in water solution was studied in the excited state. It was found that the transformations occur along the 1πσ* excited-state reaction paths which link the Franck-Condon geometries of the tautomers and the conical intersections S0/S1 connected with the H-detachment processes of the corresponding bonds. Furthermore, we established that the conical intersections S0/S1 are also mutually accessible along the 1πσ* excited-state reaction paths.

METHODS

The ground-state equilibrium geometries were optimized at the B3LYP/aug-cc-pVDZ level of theory in water environment according to PCM as well as at the CC2/aug-cc-pVDZ level in the gas phase. The TD B3LYP and CC2 methods were applied for the study of the excited states. The tautomerization mechanisms were studied with the use of the linear interpolation in internal coordinates approach using the optimized geometries of tautomers minima and conical intersections S0/S1 at the CASSCF(6,6)/6-31G* level. All calculations were performed with the GAUSSIAN 16 commercial software.

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.

The Structure of DNA Fragments: Quantum-Chemical Modelling

In this review, we analyze and systematize our computational studies of the nucleic acid duplex formations and thermodynamic stability under the different factors of investigation. The proposed structural models of mini-helix contains N nucleobase pairs (N = 3-5); QM structural data suggest that the helical conformations of mini-helix adopt geometrical parameters comparable to those of natural A- and B-DNA forms under specific conditions as micro hydration and charge compensation. The gas-phase models adopt non regular conformations between the helical form and a ladder form.. The natural helical shape of DNA mini-helix is stabilized by the presence of counterions or by explicit micro-hydration of the major and minor groves. The presence of aqueous solution is shown as a minor factor for the helical shape formation. The studies are performed at the level of density functional theory.

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.

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

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

Cryogenic IR and UV spectroscopy of isomer-selected cytosine radical cation

The UV photodissociation of cryogenic-cooled isomer-selected cytosine–silver complex leads to the production of cytosine radical cation without isomerization.

Structural Characterization of 6-Thioguanosine and Its Monohydrate in the Gas Phase

Detailed structural analysis of 6-thioguanosine (6TGs) in relation to its tautomerization and sugar conformation is performed in the gas phase using UV and IR spectroscopy combined with ab initio calculations. We have observed a thiol tautomer of 6TGs with its sugar moiety in the syn conformation that is stabilized by a strong intramolecular H-bonding between O5'H of the sugar and N3 atom of the guanine moiety. This observation is consistent with previous results for guanosine (Gs) in which the corresponding enol form is solely detected. We have also identified a monohydrate of 6TGs consisting of a thiol tautomer with the water linking guanine moiety and sugar OH group. It is demonstrated that hydration behavior of 6TGs is significantly different from that of Gs as a result of a weaker H-bonding ability of the thiol group.

Guanine Tautomerism in Ionic Complexes with Ag+ Investigated by IRMPD Spectroscopy and Mass Spectrometry

In this paper, we present the IRMPD spectra of three ionic complexes between guanine (G) and silver (Ag⁺): [GAg-H₂O]⁺, [GAgG]⁺ produced in the electrospray ionization source of the mass spectrometer, and [GAg]⁺ produced by collision induced dissociation of the [GAgG]⁺ complex. On the basis of the comparison of theoretically calculated IR spectra, we show that there are two isomers of each complex containing two different keto-amino (KA) tautomers of G (GKA(1,9) and GKA(1,7)). The observed isomers are the most stable structures in aqueous solution, and their experimentally estimated relative populations are in better agreement with the calculated relative populations in solution than in the gas phase, both at 298 K. We concluded that these observations suggest that GKA(1,9) and GKA(1,7) coexist in solution according to previous theoretical reports (Colominas, C.; et al. J. Am. Chem. Soc. 1996, 118, 6811). We were unable to find any evidence of the presence of the GEA(9), GKA(3,7), GKA(3,9), or GKA(7,9), whose relative stabilities in solution are strongly dependent on the theoretical method used to account for the solvent effect (Hanus, M.; et al. J. Am. Chem. Soc. 2003, 125, 7678).

Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds

In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm^-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH₂) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.

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.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Exploration of the theobromine-water dimer: comparison with DNA microhydration

Understanding the molecular basis of the appearance of life on Earth is an exciting research field. Many factors may have influenced the election of the molecules used by living beings and evolution may have modified those original compounds. In an attempt to understand the role played by intermolecular interactions in the election of CGAT as the alphabet of life, we present here a thorough experimental and computational study on the interaction of theobromine with water. Theobromine is a xanthine derivative, structurally related to the nucleobases, and also present in many living beings. The experimental results demonstrate that the most stable isomer of theobromine-water was formed and detected in supersonic expansions. This isomer very well resembles the structure of the dimers between nucleobases and water, offering similar values of binding energy. A comparison between the results obtained for theobromine-water with those reported in the literature for monohydrates of nucleobases is also offered.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

Tautomerism of Guanine Analogues

Tautomerism of nucleic acid (NA) bases is a crucial factor for the maintenance and translation of genetic information in organisms. Only canonical tautomers of NA bases can form hydrogen-bonded complexes with their natural counterparts. On the other hand, rare tautomers of nucleobases have been proposed to be involved in processes catalysed by NA enzymes. Isocytosine, which can be considered as a structural fragment of guanine, is known to have two stable tautomers both in solution and solid states. The tautomer equilibrium of isocytosine contrasts with the remarkable stability of the canonical tautomer of guanine. This paper investigates the factors contributing to the stability of the canonical tautomer of guanine by a combination of NMR experiments and theoretical calculations. The electronic effects of substituents on the stability of the rare tautomers of isocytosine and guanine derivatives are studied by density functional theory (DFT) calculations. Selected derivatives are studied by variable-temperature NMR spectroscopy. Rare tautomers can be stabilised in solution by intermolecular hydrogen-bonding interactions with suitable partners. These intermolecular interactions give rise to characteristic signals in proton NMR spectra, which make it possible to undoubtedly confirm the presence of a rare tautomer.

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.

Simulated NIR spectra as sensitive markers of the structure and interactions in nucleobases

Near-infrared (near-IR; NIR) spectroscopy is continuously advancing in biophysical and biochemical fields of investigation. For instance, recent progresses in NIR hyperspectral imaging of biological systems may be noted. However, interpretation of NIR bands for biological samples is difficult and creates a considerable barrier in exploring the full potential of NIR spectroscopy in bioscience. For this reason, we carried out a systematic study of NIR spectra of adenine, cytosine, guanine, and thymine in polycrystalline state. Interpretation of NIR spectra of these nucleobases was supported by anharmonic vibrational analysis using Deperturbed Vibrational Second-Order Perturbation Theory (DVPT2). A number of molecular models of nucleobases was applied to study the effect of the inter-molecular interactions on the NIR spectra. The accuracy of simulated NIR spectra appears to depend on the intra-layer interactions; in contrast, the inter-layer interactions are less influential. The best results were achieved by combining the simulated spectra of monomers and dimers. It is of particular note that in-plane deformation bands are far more populated than out-of-plane ones and the importance of ring modes is relatively small. This trend is in contrast to that observed in mid-IR region. As shown, the local, short-range chemical neighborhood of nucleobase molecules influence their NIR spectra more considerably. This suggests that NIR spectra are more sensitive probe of the nucleobase pairing than mid-IR ones. The obtained results allow, for the first time, to construct a frequency correlation table for NIR spectra of purines and pyrimidines.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

UV-promoted radical formation, and near-IR-induced and spontaneous conformational isomerization in monomeric 9-methylguanine isolated in low-temperature Ar matrices

Three low-energy isomers of 9-methylguanine, the amino-oxo (AO) form and two amino-hydroxy (AH1 and AH2) conformers, were trapped from the gas phase into low-temperature argon matrices. The AH1 and AH2 isomers, differing in the orientation of the OH group, were found to transform into each other upon excitation with near-IR light. The population of the AO form of the compound was not changed upon any near-IR irradiation of the matrix samples. Using monochromatic near-IR light, generated by a frequency-tunable laser source, it was possible to selectively induce the AH1 → AH2 or AH2 → AH1 conversion. Photoreversibility of this conformational transformation was then demonstrated. Exposure of matrix-isolated monomers of 9-methylguanine to broadband near-IR light also led to conformational conversions within the amino-hydroxy tautomeric form; the final stage of this process was always the same photostationary state independent of the initial ratio of AH1 and AH2 populations. Spontaneous conformational conversion, transforming the higher-energy AH2 form into the lower-energy AH1 isomer, was observed for matrix-isolated monomers of 9-methylguanine kept in the dark. The mechanism of this process must rely on quantum tunneling of the light hydrogen atom. Irradiation of matrix-isolated 9-methylguanine with UV laser light at λ = 288 or 285 nm led to a substantial consumption of the two AH forms, while the amount of AO isomer remained unchanged. On the other hand, a decrease in the population of the AO isomer occurred upon excitations at shorter wavelengths, λ = 280 or 275 nm. The spectral changes observed after UV-irradiation suggest the generation (and stabilization in the matrix) of a radical species, resulting from the photocleavage of the O-H or N1-H bonds, in the AH or AO isomer, respectively.

Vibrational investigation of nucleobases by means of divide and conquer semiclassical dynamics

In this work, we report a computational study of the vibrational features of four different nucleobases employing the divide-and-conquer semiclassical initial value representation molecular dynamics method. Calculations are performed on uracil, cytosine, thymine, and adenine. Results show that the overall accuracy with respect to experiments is within 20 wavenumbers, regardless of the dimensionality of the nucleobase. Vibrational estimates are accurate even in the complex case of cytosine, where two relevant conformers are taken into account. These results are promising in the perspective of future studies on more complex systems, such as nucleotides or nucleobase pairs.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

Two-dimensional self-assembled nanostructures of nucleobases and their related derivatives on Au(111)

The construction of two-dimensional (2D) self-assembled nanostructures has been one of the considerably interesting areas of on-surface chemistry in the past few decades, and has benefited from the rapid development and improvement of scanning probe microscopy techniques. In this research field, many attempts have been made in the controllable fabrication of well-ordered and multifunctional surface nanostructures, which attracted interest because of the prospect for artificial design of functional molecular nanodevices. DNA and RNA are considered to be programmable self-assembly systems and it is possible to use their base sequences to encode instructions for assembly in a predetermined fashion at the nanometer scale. As important constituents of nucleic acids, nucleobases, with intrinsic functional groups for hydrogen bonding, coordination bonding, and electrostatic interactions, can be employed as a potential system for the versatile construction of various biomolecular nanostructures, which may be used to structure the self-assembly of DNA-based artificial molecular constructions and play an important role in novel biosensors based on surface functionalization. In this article, we will review the recent progress of on-surface self-assembly of nucleobases and their derivatives together with different reactants (e.g., metals, halogens, salts and water), and as a result, various 2D surface nanostructures are summarized.

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

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

A 'bottom up', ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs

The availability of non-radiative decay mechanisms by which photoexcited molecules can revert to their ground electronic state, without experiencing potentially deleterious chemical transformation, is fundamental to molecular photostability. This Perspective Article combines results of new ab initio electronic structure calculations and prior experimental data in an effort to systematise trends in the non-radiative decay following UV excitation of selected families of heterocyclic molecules. We start with the prototypical uni- and bicyclic molecules phenol and indole, and explore the structural and photophysical consequences of incorporating progressively more nitrogen atoms within the respective ring structures en route to the DNA bases thymine, cytosine, adenine and guanine. For each of the latter, we identify low energy non-radiative decay pathways via conical intersections with the ground state potential energy surface accessed by out-of-plane ring deformations. This is followed by summary descriptions and illustrations of selected rival (electron driven H atom transfer) non-radiative excited state decay processes that demand consideration once the nucleobases are merely components in larger biomolecular systems like nucleosides, and both individual and stacked base-pairs.

How nature covers its bases

The response of nucleobases to UV radiation depends on structure in subtle ways, as revealed by gas-phase experiments.

TDDFT Study of the Optical Excitation of Nucleic Acid Bases-C60 Complexes

The potential of C₆₀ as a nucleic acid base (NAB) optical sensor is theoretically explored. We investigate the adsorption of four NABs, namely, adenine, cytosine, guanine, and thymine, on C₆₀ in the gas phase. For the optimal NAB@C₆₀ adsorption configurations, obtained using a dispersion-corrected density functional, we calculate the vis-near-ultraviolet optical response using time-dependent density functional theory. While the isolated C₆₀ and NAB molecules do not exhibit visible optical excitation, we find that C₆₀/NAB conjugation gives rise to distinct spectral features in the visible range. These results suggest that C₆₀ conjugation can be applied for photodetection of individual NABs.

Excited-State Dynamics in O6-Methylguanosine: Impact of O6-Methylation on the Relaxation Mechanism of Guanine Monomers

Absorption of ultraviolet radiation by DNA bases results in ultrafast internal conversion to the ground state, which minimizes photodamage. However, exogenous and endogenous alkylating agents present in the cellular environment can methylate the nucleobases in DNA. In particular, methylation of guanosine at the O⁶ position in DNA leads to the formation of the O⁶-methylguanosine adduct, which may alter the photostability of DNA. This contribution demonstrates that O⁶-methylation of guanosine red shifts its ground-state absorption spectrum and slows down the rate of internal conversion to the ground state by ∼40-fold in aqueous solution. The 40-fold decrease in the rate of excited-state decay increases the probability of photodamage within cellular DNA. It is proposed that the longer decay lifetime corresponds to relaxation of the excited-state population in O⁶-methylguanosine along a C6-puckered reaction coordinate in the ¹ππ*(L_a) potential energy surface that runs parallel to an ultrafast internal conversion pathway along a C2-puckered coordinate.

Site-dependent effects of methylation on the electronic spectra of jet-cooled methylated xanthine compounds

We obtained the electronic spectra of various methylated xanthine compounds including caffeine in a supersonic jet by resonant two-photon ionization spectroscopy. The methyl group in the tested methylated xanthine compounds has a distinct, site-dependent effect on the electronic spectrum. Methylation at the N3 position causes a significant red shift of the ππ* state, whereas methylation at the N1 position has only minimal effects on the electronic spectrum. The notably broad spectra of theobromine and caffeine result from methyl substitution at the N7 position, which causes a large displacement between the potential energy surfaces of the S₀ and S₁ states, and a strong vibronic coupling. We also investigated the internal rotation of the methyl group and its effect on the electronic spectrum of the methylated xanthine compounds. We found that the barrier height for the torsional motion in the ground state is significantly affected by a carbonyl or methyl group that lies close to the methyl group of interest. In contrast, the torsional barrier in the excited state is governed by the hyperconjugation interaction in the lowest unoccupied molecular orbital. The agreement between the experimental and simulated spectra of torsional vibronic bands suggested that the low frequency torsional vibrations arising from the tunneling splitting and the coupling between the torsional and molecular motions give theobromine and theophylline the multiplet nature of their origin bands. This study provides a new level of understanding for the methyl substitution effects on the electronically excited states of xanthine compounds, which may very well be applicable to many other methyl substituted biomolecules including DNAs and proteins.

Sugar-peptidic bond interactions: spectroscopic characterization of a model system

Sugars are small carbohydrates which play numerous roles in living organisms such as storage of energy or as structural components. Modifications of specific sites within the glycan chain can modulate a carbohydrate's overall biological function as it happens with nucleic acids and proteins. Hence, identifying discrete carbohydrate modifications and understanding their biological effects is essential. A study of such processes requires of a deep knowledge of the interaction mechanism at the molecular level. Here, we use a combination of laser spectroscopy in jets and quantum mechanical calculations to characterize the interaction between phenyl-β-d-glucopyranoside and N-methylacetamide as a model to understand the interaction between a sugar and a peptide bond. The most stable structure of the molecular aggregate shows that the main interaction between the peptide fragment and the sugar proceeds via a C[double bond, length as m-dash]OH-O2 hydrogen bond. A second conformer was also found, in which the peptide establishes a C[double bond, length as m-dash]OH-O6 hydrogen bond with the hydroxymethyl substituent of the sugar unit. All the conformers present an additional interaction point with the aromatic ring. This particular preference of the peptide for the hydroxyl close to the aromatic ring could explain why glycogenin uses tyrosine in order to convert glucose into glycogen by exposing the O4H hydroxyl group for the other glucoses for the polymerization to take place.

Gas-Phase Conformations and N-Glycosidic Bond Stabilities of Sodium Cationized 2'-Deoxyguanosine and Guanosine: Sodium Cations Preferentially Bind to the Guanine Residue

2'-Deoxyguanosine (dGuo) and guanosine (Guo) are fundamental building blocks of DNA and RNA nucleic acids. In order to understand the effects of sodium cationization on the gas-phase conformations and stabilities of dGuo and Guo, infrared multiple photon dissociation (IRMPD) action spectroscopy experiments and complementary electronic structure calculations are performed. The measured IRMPD spectra of [dGuo+Na]⁺ and [Guo+Na]⁺ are compared to calculated IR spectra predicted for the stable low-energy structures computed for these species to determine the most favorable sodium cation binding sites, identify the structures populated in the experiments, and elucidate the influence of the 2'-hydroxyl substituent on the structures and IRMPD spectral features. These results are compared with those from a previous IRMPD study of the protonated guanine nucleosides to elucidate the differences between sodium cationization and protonation on structure. Energy-resolved collision-induced dissociation (ER-CID) experiments and survival yield analyses of protonated and sodium cationized dGuo and Guo are performed to compare the effects of these cations toward activating the N-glycosidic bonds of these nucleosides. For both [dGuo+Na]⁺ and [Guo+Na]⁺, the gas-phase structures populated in the experiments are found to involve bidentate binding of the sodium cation to the O6 and N7 atoms of guanine, forming a 5-membered chelation ring, with guanine found in both anti and syn orientations and C2'-endo (²T₃ or ³T₂) puckering of the sugar. The ER-CID results, IRMPD yields and the computed C1'-N9 bond lengths indicate that sodium cationization activates the N-glycosidic bond less effectively than protonation for both dGuo and Guo. The 2'-hydroxyl substituent of Guo is found to impact the preferred structures very little except that it enables a 2'OH···3'OH hydrogen bond to be formed, and stabilizes the N-glycosidic bond relative to that of dGuo in both the sodium cationized and protonated complexes.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated 2'-deoxyguanosine and guanosine

Experimental and theoretical investigations suggest that hydrolysis of N-glycosidic bonds generally involves a concerted SN2 or a stepwise SN1 mechanism. While theoretical investigations have provided estimates for the intrinsic activation energies associated with N-glycosidic bond cleavage reactions, experimental measurements to validate the theoretical studies remain elusive. Here we report experimental investigations for N-glycosidic bond cleavage of the protonated guanine nucleosides, [dGuo+H](+) and [Guo+H](+), using threshold collision-induced dissociation (TCID) techniques. Two major dissociation pathways involving N-glycosidic bond cleavage, resulting in production of protonated guanine or the elimination of neutral guanine are observed in competition for both [dGuo+H](+) and [Guo+H](+). The detailed mechanistic pathways for the N-glycosidic bond cleavage reactions observed are mapped via electronic structure calculations. Excellent agreement between the measured and B3LYP calculated activation energies and reaction enthalpies for N-glycosidic bond cleavage of [dGuo+H](+) and [Guo+H](+) in the gas phase is found indicating that these dissociation pathways involve stepwise E1 mechanisms in analogy to the SN1 mechanisms that occur in the condensed phase. In contrast, MP2 is found to significantly overestimate the activation energies and slightly overestimate the reaction enthalpies. The 2'-hydroxyl substituent is found to stabilize the N-glycosidic bond such that [Guo+H](+) requires ∼25 kJ mol(-1) more than [dGuo+H](+) to activate the glycosidic bond.

Capturing Transient Endoperoxide in the Singlet Oxygen Oxidation of Guanine

The chemistry of singlet O2 toward the guanine base of DNA is highly relevant to DNA lesion, mutation, cell death, and pathological conditions. This oxidative damage is initiated by the formation of a transient endoperoxide through the Diels-Alder cycloaddition of singlet O2 to the guanine imidazole ring. However, no endoperoxide formation was directly detected in native guanine or guanosine, even at -100 °C. Herein, gas-phase ion-molecule scattering mass spectrometry was utilized to capture unstable endoperoxides in the collisions of hydrated guanine ions (protonated or deprotonated) with singlet O2 at ambient temperature. Corroborated by results from potential energy surface exploration, kinetic modeling, and dynamics simulations, various aspects of endoperoxide formation and transformation (including its dependence on guanine ionization and hydration states, as well as on collision energy) were determined. This work has pieced together reaction mechanisms, kinetics, and dynamics data concerning the early stage of singlet O2 induced guanine oxidation, which is missing from conventional condensed-phase studies.

Hydrogen-bonding interactions of uric acid complexes with water/melamine by mid-infrared spectroscopy

Hydrogen (H)-bonding interactions of uric acid (UA) with water have been investigated via IR-UV double resonance measurements in the mid-IR region. Comparison of the present results with those obtained previously in the near-IR region enables us to examine microscopic hydration effects that are specific to the H-bonding acceptor sites of UA. It is shown that hydration of the C8O site promotes the mode coupling of this stretch with the C2O stretch. The occurrence of this coupling is manifested in the IR intensity pattern, in which the transition associated with the in-phase contribution C8O + C2O is significantly suppressed, whereas the corresponding out-of-phase contribution gives rise to a strong peak. We also measured the mid-IR spectra of the 1 : 1 complex formed between UA and melamine (MEL) and carried out a structural analysis using the spectroscopic signature of the H-bonding derived from the result of the monohydrated cluster. It is shown that the complex possesses a triple H-bonding structure with the C2O acceptor site of UA H-bonded to MEL. Furthermore, the IR-depleted UV spectroscopy technique was employed in order to ascertain whether other structural isomers are present in the probe UV spectra.

Gas-Phase IR Spectroscopy of Nucleobases

IR spectroscopy of nucleobases in the gas phase reflects simultaneous advances in both experimental and computational techniques. Important properties, such as excited state dynamics, depend in subtle ways on structure variations, which can be followed by their infrared signatures. Isomer specific spectroscopy is a particularly powerful tool for studying the effects of nucleobase tautomeric form and base pair hydrogen-bonding patterns.

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.

Direct observation of multiple tautomers of oxythiamine and their recognition by the thiamine pyrophosphate riboswitch

Structural diversification of canonical nucleic acid bases and nucleotide analogues by tautomerism has been proposed to be a powerful on/off switching mechanism allowing regulation of many biological processes mediated by RNA enzymes and aptamers. Despite the suspected biological importance of tautomerism, attempts to observe minor tautomeric forms in nucleic acid or hybrid nucleic acid-ligand complexes have met with challenges due to the lack of sensitive methods. Here, a combination of spectroscopic, biochemical, and computational tools probed tautomerism in the context of an RNA aptamer-ligand complex; studies involved a model ligand, oxythiamine pyrophosphate (OxyTPP), bound to the thiamine pyrophosphate (TPP) riboswitch (an RNA aptamer) as well as its unbound nonphosphorylated form, oxythiamine (OxyT). OxyTPP, similarly to canonical heteroaromatic nucleic acid bases, has a pyrimidine ring that forms hydrogen bonding interactions with the riboswitch. Tautomerism was established using two-dimensional infrared (2D IR) spectroscopy, variable temperature FTIR and NMR spectroscopies, binding isotope effects (BIEs), and computational methods. All three possible tautomers of OxyT, including the minor enol tautomer, were directly identified, and their distributions were quantitated. In the bound form, BIE data suggested that OxyTPP existed as a 4'-keto tautomer that was likely protonated at the N1'-position. These results also provide a mechanistic framework for understanding the activation of riboswitch in response to deamination of the active form of vitamin B1 (or TPP). The combination of methods reported here revealing the fine details of tautomerism can be applied to other systems where the importance of tautomerism is suspected.

Increasing occurrences and functional roles for high energy purine-pyrimidine base-pairs in nucleic acids

There are a growing number of studies reporting the observation of purine-pyrimidine base-pairs that are seldom observed in unmodified nucleic acids because they entail the loss of energetically favorable interactions or require energetically costly base ionization or tautomerization. These high energy purine-pyrimidine base-pairs include G•C(+) and A•T Hoogsteen base-pairs, which entail ∼180° rotation of the purine base in a Watson-Crick base-pair, protonation of cytosine N3, and constriction of the C1'-C1' distance by ∼2.5Å. Other high energy pure-pyrimidine base-pairs include G•T, G•U, and A•C mispairs that adopt Watson-Crick like geometry through either base ionization or tautomerization. Although difficult to detect and characterize using biophysical methods, high energy purine-pyrimidine base-pairs appear to be more common than once thought. They further expand the structural and functional diversity of canonical and non-canonical nucleic acid base-pairs.

Coupled-cluster and density functional theory studies of the electronic 0–0 transitions of the DNA bases

The 0–0 transitions of the electronic excitation spectra of the lowest tautomers of the four nucleotide (DNA) bases have been studied using linear-response approximate coupled-cluster singles and doubles (CC2) calculations.

Time-resolved photoelectron imaging of the isolated deprotonated nucleotides

Time-resolved photoelectron spectroscopy of deprotonated nucleotides provides new insights into their relaxation dynamics.

Spectroscopic signatures and structural motifs in isolated and hydrated caffeine: a computational study

Blue and red shift of the VEE of the S₁ (¹ππ*) state of caffeine after forming complex with water at isolated and conjugated carbonyl site.