Correlated ab initio study of nucleic acid bases and their tautomers in the gas phase, in a microhydrated environment and in aqueous solution. guanine: surprising stabilization of rare tautomers in aqueous solution

IR and UV Spectroscopy of Gas-Phase Monohydrated Protonated Guanine

We use UV and infrared photodissociation spectroscopy to study monohydrated protonated guanine in a dual cryogenic ion trap spectrometer. The monohydrated complexes are formed through helium-mediated collisions between bare electrosprayed protonated guanine and low-pressure water vapor in a clustering trap maintained at 180 K, before being transferred to a quadrupole ion trap at 10 K. The spectrum of the monohydrated complex exhibits sharp vibronic transitions at the band origin and becomes broader and higher in intensity further in blue, which is very similar to protonated guanine but with a notable blue shift of ∼1850 cm-1 (∼0.23 eV). The UV hole-burning experiments showed that the vibronic bands recorded in the region of the band origin belong to a single conformer under our experimental conditions. The IR photodissociation spectrum in the 3000-3600 cm-1 range, with the aid of theoretical calculations (SCS-CC2/aug-cc-pVDZ), allowed us to assign the structure to the lowest energy N7-O conformer.

A rotational spectroscopy study of microsolvation effects on intramolecular proton transfer in trifluoroacetylacetone-(H2O)1-3

Trifluoroacetylacetone (TFAA) has two enol forms, which can switch to each other via proton transfer. While much attention has been paid to their conformational preferences, the influence of microsolvation on regulating the proton position remains unexplored. Herein, we report the rotational spectra of trifluoroacetylacetone-(water)n (n = 1-3) investigated by chirped pulse Fourier transform microwave spectroscopy in the 2-8 GHz frequency range. Two conformers were identified for both TFAA-H2O and TFAA-(H2O)2, while only one conformer was characterized for TFAA-(H2O)3. The results indicate that water binding on the CH3 side stabilizes the enolF form, whereas water binding on the CF3 side stabilizes the enolH form. The enolF form predominates over the enolH form in these hydrated complexes, which contrasts with the fact that only enolH exists in isolated TFAA. EnolH becomes preferred only when water inserts itself into the intramolecular hydrogen bond. Instanton theory calculations reveal that the proton transfer reaction is dominated by quantum tunneling at low temperatures, leading to the stable existence of only one enol form in each configuration of the hydrated clusters.

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

Water-Controlled Keto-Enol Tautomerization of a Prebiotic Nucleobase

Barbituric acid is believed to be a proto-RNA nucleobase that was used for biological information transfer on prebiotic earth before DNA and RNA in their present forms evolved. Nucleobases have various tautomeric forms and the relative stability of these forms is critical to their biological function. It has been shown that barbituric acid has a tri-keto form in the gas phase and an enol form in the solid state. However, its dominant tautomeric form in aqueous medium that is most relevant for biology has been investigated only to a limited extent and the findings are inconclusive. We have used multiple approaches, namely, molecular dynamics, quantum chemistry, NMR, and IR spectroscopy to determine the most stable tautomer of barbituric acid in solution. We find a delicate balance in the stability of the two tautomers, tri-keto and enol, which is tipped toward the enol as the extent of solvation by water increases.

Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches

The determination of the redox properties of nucleobases is of paramount importance to get insight into the charge-transfer processes in which they are involved, such as those occurring in DNA-inspired biosensors. Although many theoretical and experimental studies have been conducted, the value of the one-electron oxidation potentials of nucleobases is not well-defined. Moreover, the most appropriate theoretical protocol to model the redox properties has not been established yet. In this work, we have implemented and evaluated different static and dynamic approaches to compute the one-electron oxidation potentials of solvated nucleobases. In the static framework, two thermodynamic cycles have been tested to assess their accuracy against the direct determination of oxidation potentials from the adiabatic ionization energies. Then, the introduction of vibrational sampling, the effect of implicit and explicit solvation models, and the application of the Marcus theory have been analyzed through dynamic methods. The results revealed that the static direct determination provides more accurate results than thermodynamic cycles. Moreover, the effect of sampling has not shown to be relevant, and the results are improved within the dynamic framework when the Marcus theory is applied, especially in explicit solvent, with respect to the direct approach. Finally, the presence of different tautomers in water does not affect significantly the one-electron oxidation potentials.

Noncanonical Isomers of Nucleoside Cation Radicals: An Ab Initio Study of the Dark Matter of DNA Ionization

Cation radicals of DNA nucleosides, 2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and 2'-deoxythymidine, can exist in standard canonical forms or as noncanonical isomers in which the charge is introduced by protonation of the nucleobase, whereas the radical predominantly resides in the deoxyribose moiety. Density functional theory as well as correlated ab initio calculations with coupled clusters (CCSD(T)) that were extrapolated to the complete basis set limit showed that noncanonical nucleoside ion isomers were thermodynamically more stable than their canonical forms in both the gas phase and as water-solvated ions. This indicated the possibility of exothermic conversion of canonical to noncanonical forms. The noncanonical isomers were calculated to have very low adiabatic ion-electron recombination energies (RE_ad) for the lowest-energy isomers 2'-deoxy-(N-3H)adenos-1'-yl (4.74 eV), 2'-deoxy-(N-7H)guanos-1'-yl (4.66 eV), 2'-deoxy-(N-3H)cytid-1'-yl (5.12 eV), and 2'-deoxy-5-methylene-(O-2H)uridine (5.24 eV). These were substantially lower than the RE_ad value calculated for the canonical 2'-deoxyadenosine, 2'-deoxy guanosine, 2'-deoxy cytidine, and 2'-deoxy thymidine cation radicals, which were 7.82, 7.46, 8.14, and 8.20 eV, respectively, for the lowest-energy ion conformers of each type. Charge and spin distributions in noncovalent cation-radical dA⊂dT and dG⊂dC nucleoside pairs and dAT, dCT, dTC, and dGC dinucleotides were analyzed to elucidate the electronic structure of the cation radicals. Born-Oppenheimer molecular dynamics trajectory calculations of the dinucleotides and nucleoside pairs indicated rapid exothermic proton transfer from noncanonical T^+· to A in both dAT^+· and dA⊂dT^+·, leading to charge and radical separation. Noncanonical T^+· in dCT^+· and dTC^+· initiated rapid proton transfer to cytosine, whereas the canonical dCT^+· dinucleotide ion retained the cation radical structure without isomerization. No spontaneous proton transfer was found in dGC^+· and dG⊂dC^+· containing canonical neutral and noncanonical ionized deoxycytidine.

The Effects of Potential and pH on the Adsorption of Guanine on the Au(111) Electrode

The adsorption of nucleobase at a gold electrode has been a model system to study the interaction between biomolecule and metal, which is relevant to the development of sensors and molecular electronics. The current study has employed in situ scanning tunneling microscopy (STM) and voltammetry to investigate the adsorption configuration and spatial structure of guanine (G) on a well-defined Au(111) electrode in perchloric acid (HClO₄) and neutral phosphate buffer solution (PBS) containing 50 μM G. Potential control had a profound effect on the adsorption of the G molecule on the Au(111) electrode. No adsorption of G was observed at a potential more negative than 0 V in HClO₄ and -0.2 V (versus Ag/AgCl) in PBS; shifting potential positively triggered a rapid adsorption of G to yield a well-ordered G array. Different spatial structures of G admolecules were imaged with STM in HClO₄ and PBS, suggesting that ions in the electrolyte were important in this adsorption event. Shifting potential positively caused a more compact G adlayer with molecules adopting a tilted orientation. Meanwhile, G molecules continued to deposit on the Au(111) electrode leading to a multilayer G film. These processes were reversible to the potential modulation. G admolecules on the Au(111) electrode could be irreversibly oxidized in 0.1 M PBS, which resulted in a prominent peak at 0.74 V in the voltammogram. This oxidation process could be used to analyze the G molecule in a sample.

Investigation of interactions of doxorubicin with purine nucleobases by molecular modeling

Doxorubicin, an anthracycline antibiotic with anti-tumor activity, is produced by the bacterium Streptomyces peucetius. The interactions between doxorubicin and genetic material and the details of the intercalation with DNA have been controversial issues. Thus, the interactions of doxorubicin with purine nucleobases were studied by quantum mechanical methods. Initially, conformer analyses of doxorubicin were performed with Spartan 08 software and 319 different conformers from 422 initial structures for doxorubicin were obtained. Geometry optimizations and frequency analyses were performed for each structure using density functional theory (DFT) at B3LYP/6-31G** level using Gaussian 09 software. The most stable 20 conformers of doxorubicin and tautomers of purine nucleobases were optimized again with ɷB97XD/6-31G** level and their interactions were also analyzed at the same level. The Discovery Studio 3.5 Visualizer was used to draw the initial and optimized structures of investigated geometries. The noncovalent interactions (NCIs) were visualized by calculating reduced density gradient (RDG) with Multiwfn program. The color-filled isosurfaces and RDG scatter maps of most stable interaction geometries were plotted by Visual Molecular Dynamics (VMD) software and Gnuplot 5.3 software, respectively. This study showed that adenine, guanine, and hypoxanthine nucleobases interact with doxorubicin by forming strong hydrogen bonds and π-π interactions. Considering the normal cellular conditions, the effect of solvent (water) on the interaction geometries were also analyzed and when compared to gas phase it was determined that the movements of the molecules were restricted and there was a minimal change between initial and optimized structures in the aqueous phase.

Photophysical properties of fluorescent nucleobase P-analogues expected to monitor DNA replication

In this work, we computationally design a series of fluorescent purine analogues based on the 2-amino-8-(1'-β-D-2'-deoxyribofuranosyl)-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one (P) to monitor the DNA replication process with merely a minimal perturbation to the natural structure of nucleic acid. The P-modified fluorescent probes present red-shifted absorption spectra and enhanced photoluminescence due to the additional π-conjugation resulting from the fluorophore modification and the ring-expansion. Efficient fluorescence quenching of P-analogues occurs upon pairing with the complementary 6-amino-5-nitro-3-(1'-β-D-2'-deoxyribofuranosyl)-2(1H)-pyridone (Z) due to the nonradiative relaxation from the low-lying dark excited state to the ground state of Z moiety. Especially, the P3 and the P7, which have high fluorescence intensity in both gas and liquid phases, are proposed as the sensors for studying conformational switching in the presence and absence of a complementary sequence. Also examined are the influences of hydration and the linking to deoxyribose on absorption and emission processes. Besides, the potential phosphorescence emission of these modified base pairs is taken into account by constructing the relaxed potential energy curves of S₀, T₁ and S₁ states.

Self-catalytic mechanism of prebiotic reactions: from formamide to pterins and guanine

Reaction pathway of prebiotic reactions for formation of the pteridines: pterin, xanthopterine, isoxanthopterine and leucopterine, as well as the purine nucleobase guanine from pure formamide are presented. In these reactions, formamide or its tautomer, formimidic acid, play the role of proton-carrying catalyst. All required raw materials, such as hydrogen cyanide, ammonia, water, formic acid, urea, 2-aminomalononitrile, glyoxal, glyoxylic acid and oxalic acid needed in the self-catalyzed reactions are obtained by partial decomposition of formamide. We show that the prebiotic formation of nucleobases and pterins is closely linked and they probably coexisted at the beginning of chemical evolution.

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).

Energetic Tricyclic Polynitropyrazole and Its Salts: Proton-Locking Effect of Guanidium Cations

An axisymmetric polynitro-pyrazole molecule, 3,5-di(3,5-dinitropyrazol-4-yl)]-4-nitro-1H-pyrazole (5), and its salts (6-12) were prepared and fully characterized. These compounds not only show promising energetic properties but also show a unique tautomeric switch via combining different cations with the axisymmetric compound (5). Its salts (6-9) remain axisymmetric when the cations are potassium, ammonium, or amino-1,2,4-triazolium. However, when the cations are guanidiums, the salts (10-12) dramatically become asymmetric owing to the fixed proton. The introduction of guanidium cations breaks the tautomeric equilibrium by blocking the prototropic transformations and results in the switch-off effect to tautomerism. The structural constraints of ¹H NMR and ¹³C NMR spectra provide strong evidence for the unusual structural constraint phenomenon. These stabilized asymmetric tautomers are very important from the point of molecular recognition, and this research may promote further developments in synthetic and isolation methodologies for novel bioactive pyrazole-based compounds.

Tackling a Curious Case: Generation of Charge-Tagged Guanosine Radicals by Gas-Phase Electron Transfer and Their Characterization by UV-vis Photodissociation Action Spectroscopy and Theory

We report the generation of gas-phase riboguanosine radicals that were tagged at ribose with a fixed-charge 6-(trimethylammonium)hexane-1-aminocarbonyl group. The radical generation relied on electron transfer from fluoranthene anion to noncovalent dibenzocrown-ether dication complexes which formed nucleoside cation radicals upon one-electron reduction and crown-ether ligand loss. The cation radicals were characterized by collision-induced dissociation (CID), photodissociation (UVPD), and UV-vis action spectroscopy. Identification of charge-tagged guanosine radicals was challenging because of spontaneous dissociations by loss of a hydrogen atom and guanine that occurred upon storing the ions in the ion trap without further excitation. The loss of H proceeded from an exchangeable position on N-7 in guanine that was established by deuterium labeling and was the lowest energy dissociation of the guanosine radicals according to transition-state energy calculations. Rate constant measurements revealed an inverse isotope effect on the loss of either hydrogen or deuterium with rate constants k_H = 0.25-0.26 s^-1 and k_D = 0.39-0.54 s^-1. We used time-dependent density functional theory calculations, including thermal vibronic effects, to predict the absorption spectra of several protomeric radical isomers. The calculated spectra of low-energy N-7-H guanine-radical tautomers closely matched the action spectra. Transition-state-theory calculations of the rate constants for the loss of H-7 and guanine agreed with the experimental rate constants for a narrow range of ion effective temperatures. Our calculations suggest that the observed inverse isotope effect does not arise from the isotope-dependent differences in the transition-state energies. Instead, it may be caused by the dynamics of post-transition-state complexes preceding the product separation.

A theoretical study on the role of stability of cytosine and its tautomers in DNA (deoxyribonucleic acid), and investigation of interactions of Na+, K+, Mg2+, Ca2+, Zn2+ metal ions and OH radical with cytosine tautomers

In the present study, 21 cytosine tautomers were investigated so that some tautomers were reported for the first time in the gas phase and aqueous solution. C₃ tautomer was the most stable tautomer in gas phase but C₁ was the most stable structure in aqueous solution. The potential energy surface of all trajectories was determined for 21 tautomers and 22 transition states. Also, interactions of cytosine tautomers with Na⁺, K⁺, Mg²⁺, Ca²⁺ and Zn²⁺ metal ions were studied in gas phase and aqueous solution. Three types of interactions among metal ions and (N1 and O10), (N3 and O10) and (N3 and N9) of cytosine tautomers were investigated. The study of interaction energies of all complexes showed the stability of complexes in which interactions among Mg²⁺ and Zn²⁺ with tautomers were stronger than interactions among Ca²⁺, Na⁺ and K⁺ with tautomers, respectively. Some interactions of metal ions with cytosine tautomers made the most stable tautomers. So, the stability of rare tutomeric forms had a significant effect on stabilization of anomalous DNA (deoxyribonucleic acid) double helix and spontaneous mutations. Also, one of the most important causes of mutations in DNA (deoxyribonucleic acid) was the reaction of OH radical with nucleotide bases. So, interactions of OH radical with cytosine and its tautomers were investigated in gas phase and aqueous solution.Communicated by Ramaswamy H. Sarma.

C5-Substituted 2-Selenouridines Ensure Efficient Base Pairing with Guanosine; Consequences for Reading the NNG-3' Synonymous mRNA Codons

5-Substituted 2-selenouridines (R5Se2U) are post-transcriptional modifications present in the first anticodon position of transfer RNA. Their functional role in the regulation of gene expression is elusive. Here, we present efficient syntheses of 5-methylaminomethyl-2-selenouridine (1, mnm5Se2U), 5-carboxymethylaminomethyl-2-selenouridine (2, cmnm5Se2U), and Se2U (3) alongside the crystal structure of the latter nucleoside. By using pH-dependent potentiometric titration, pKa values for the N3H groups of 1–3 were assessed to be significantly lower compared to their 2-thio- and 2-oxo-congeners. At physiological conditions (pH 7.4), Se2-uridines 1 and 2 preferentially adopted the zwitterionic form (ZI, ca. 90%), with the positive charge located at the amino alkyl side chain and the negative charge at the Se2-N3-O4 edge. As shown by density functional theory (DFT) calculations, this ZI form efficiently bound to guanine, forming the so-called “new wobble base pair”, which was accepted by the ribosome architecture. These data suggest that the tRNA anticodons with wobble R5Se2Us may preferentially read the 5′-NNG-3′ synonymous codons, unlike their 2-thio- and 2-oxo-precursors, which preferentially read the 5′-NNA-3′ codons. Thus, the interplay between the levels of U-, S2U- and Se2U-tRNA may have a dominant role in the epitranscriptomic regulation of gene expression via reading of the synonymous 3′-A- and 3′-G-ending codons.

Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet

Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N₁-H and N₃-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.

Why Purine Nucleoside Phosphorylase Ribosylates 2,6-Diamino-8-azapurine in Noncanonical Positions? A Molecular Modeling Study

Protein nucleoside phosphorylase (PNP) is an enzyme that catalyzes a reversible conversion process (ribosylation and phosphorolysis) between nucleobases (purines) and their nucleosides. Experimental studies showed that calf PNP ribosylates purine analogues in specific positions: 2,6-diamino-8-azapurine in position 7 or 8 and 8-azaguanine in position 9 of the triazole ring. The reason for this phenomenon can be a result of different expositions of purine substrates to the channel leading to the binding site. This hypothesis was verified by the application of molecular modeling techniques to two complexes of purine analogues 2,6-diamino-azapurine, calf PNP (pdb-code: 1LVU), and 8-azaguanine, calf PNP (pdb-code: 2AI1). The results obtained with a combination of quantum chemistry, docking, and molecular dynamics methods showed qualitative validity of our hypothesis. Binding free energies of protein-ligand systems showed that most probable binding poses expose N8 nitrogen for 2,6-diamino-8-azapurine and N9 nitrogen for 8-azaguanine into the binding channel and ruled out the exposition of N9 for 2,6-diamino-8-azapurine and N7 for 8-azaguanine, partially in agreement with the experimental data. The other important result obtained in this study is a significantly higher population of the protonated form of crucial residue Glu-201 present in the binding pocket, compared to the standard protonation of free glutamic acid in solution. This result combined with populations of tautomeric forms of both investigated systems strongly suggests that 2,6-diamino-8-azapurine and 8-azaguanine are recognized by proteins with deprotonated and protonated Glu-201 residues, respectively. A comparison of computed binding poses of the investigated ligands to the inhibitors present in crystal structures suggests that the modification of the (S)-PMPDAP inhibitor, in which a 2-(phosphonomethoxy)propyl chain is attached at position 8 instead of position 9, might increase its binding affinity.

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.

Excited-State Dynamics of Thienoguanosine, an Isomorphic Highly Fluorescent Analogue of Guanosine

Thienoguanosine (^th G) is an isomorphic analogue of guanosine with promising potentialities as fluorescent DNA label. As a free probe in protic solvents, ^th G exists in two tautomeric forms, identified as the H1, being the only one observed in nonprotic solvents, and H3 keto-amino tautomers. We herein investigate the photophysics of ^th G in solvents of different polarity, from water to dioxane, by combining time-resolved fluorescence with PCM/TD-DFT and CASSCF calculations. Fluorescence lifetimes of 14.5-20.5 and 7-13 ns were observed for the H1 and H3 tautomers, respectively, in the tested solvents. In methanol and ethanol, an additional fluorescent decay lifetime (≈3 ns) at the blue emission side (λ≈430 nm) as well as a 0.5 ns component with negative amplitude at the red edge of the spectrum, typical of an excited-state reaction, were observed. Our computational analysis explains the solvent effects observed on the tautomeric equilibrium. The main radiative and nonradiative deactivation routes have been mapped by PCM/TD-DFT calculations in solution and CASSCF in the gas phase. The most easily accessible conical intersection, involving an out-of plane motion of the sulfur atom in the five-membered ring of ^th G, is separated by a sizeable energy barrier (≥0.4 eV) from the minimum of the spectroscopic state, which explains the large experimental fluorescence quantum yield.

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).

Nickel Adatoms Induced Tautomeric Dehydrogenation of Thymine Molecules on Au(111)

Tautomerization of nucleobases may induce base mismatches resulting in the abnormal disturbance of gene replication and expression, which has therefore attracted widespread interests in many disciplines. Metal atoms participating in a variety of important biological processes are found to be able to affect the nucleobase tautomerization as evidenced by many theoretical and spectroscopic studies. To get the real-space evidence and to unravel the underlying mechanism for the metal-induced tautomerization, especially from the keto form to the enol one, the interplay of high-resolution scanning tunneling microscopy imaging/manipulation and density functional theory (DFT) calculations has been employed. We present a process showing the Ni adatom-induced keto-enol tautomeric dehydrogenation of thymine molecules on Au(111). The key to making such a process feasible is the Ni atoms which greatly lower the energy barrier for the tautomerization from keto to enol form, which is rationalized by extensive DFT-based transition-state search calculations.

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.

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.

A general spin-complete spin-flip configuration interaction method

A new spin-complete spin-flip configuration interaction (SF-CI) method was developed using the ORMAS-CI algorithm, along with a perturbative correction for capturing dynamic and non-dynamic correlation.

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.

Tautomeric equilibria of isoguanine and related purine analogs

Nucleobase pairs in DNA match hydrogen-bond donor and acceptor groups on the nucleobases. However, these can adopt more than one tautomeric form, and can consequently pair with nucleobases other than their canonical complements, possibly a source of natural mutation. These issues are now being re-visited by synthetic biologists increasing the number of replicable pairs in DNA by exploiting unnatural hydrogen bonding patterns, where tautomerism can also create mutation. Here, we combine spectroscopic measurements on methylated analogs of isoguanine tautomers and tautomeric mixtures with statistical analyses to a set of isoguanine analogs, the complement of isocytosine, the 5th and 6th "letters" in DNA.

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.

Reaction electronic flux and its role in DNA intramolecular proton transfers

Proton transfer reactions present a key step in many biological and chemical processes. Here, we focused on the electronic changes in the proton transfer reactions of the four DNA bases. In combination with the previous structural analysis the reaction electronic flux together with local descriptors as the Hirshfeld-I charges allow us to identify chemical events and rationalize the underlying reaction mechanism. Our results show that imine-enamine in adenine and citosyne, and keto-enol tautomerizations in thymine and guanine have different reaction mechanisms. The former involve net structural rearrangements driven by favoured electrostatic interactions between the proton and the acceptor atom whereas the keto-enol tautomerizations require electronic changes reflected in the reaction electronic flux and changes in the NBO bond orders which favour the proton transfer reaction.

Molecular thermodynamics of metabolism: quantum thermochemical calculations for key metabolites

The present work is the first of a series of papers aiming at a coherent and unified development of the thermodynamics of metabolism and the rationalization of feasibility analysis of metabolic pathways. The focus in this part is on high-level quantum chemical calculations of the thermochemical quantities of relatively heavy metabolites such as amino acids/oligopeptides, nucleosides, saccharides and their derivatives in the ideal gas state. The results of this study will be combined with the corresponding hydration/solvation results in subsequent parts of this work in order to derive the desired thermochemical quantities in aqueous solutions. The above metabolites exist in a vast conformational/isomerization space including rotational conformers, tautomers or anomers exhibiting often multiple or cooperative intramolecular hydrogen bonding. We examine the challenges posed by these features for the reliable estimation of thermochemical quantities. We discuss conformer search, conformer distribution and averaging processes. We further consider neutral metabolites as well as protonated and deprotonated metabolites. In addition to the traditional presentation of gas-phase acidities, basicities and proton affinities, we also examine heats and free energies of ionic species. We obtain simple linear relations between the thermochemical quantities of ions and the formation quantities of their neutral counterparts. Furthermore, we compare our calculations with reliable experimental measurements and predictive calculations from the literature, when available. Finally, we discuss the next steps and perspectives for this work.

Real-Space Evidence of Rare Guanine Tautomer Induced by Water

Water is vital for life as a solvent. Specifically, it has been well established that DNA molecules are hydrated in vivo, and water has been found to be responsible for the presence of some noncanonical DNA base tautomers. Theoretical investigations have shown that the existence of water could significantly influence the relative stability of different DNA base tautomers, reduce the energy barrier of tautomeric conversions, and thus promote the formation of some rare base tautomers. In this work, we report the real-space experimental evidence of rare base tautomers. From the high-resolution scanning tunneling microscopy imaging, we surprisingly find the formation of the rare guanine tautomer, i.e., G/(3H,7H) form, on the Au(111) surface by delicately introducing water into the system. The key to the formation of this rare tautomer is proposed to be the "water bridge" that largely reduces the energy barriers of intramolecular proton-transfer processes as revealed by extensive density functional theory calculations. The real-space experimental evidence and the proposed mechanism make a step forward toward the fundamental understanding of water-assisted base tautomerization processes.

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.

Nanohydration of uracil: emergence of three-dimensional structures and proton-induced charge transfer

Stepwise hydration of uracil proceeds three dimensionally above three molecules and qualitatively changes the response to proton damage.

Gas-Phase Conformations and Energetics of Protonated 2'-Deoxyguanosine and Guanosine: IRMPD Action Spectroscopy and Theoretical Studies

The gas-phase structures of protonated 2'-deoxyguanosine, [dGuo+H](+), and its RNA analogue protonated guanosine, [Guo+H](+), are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical electronic structure calculations. IRMPD action spectra are measured over the range extending from ∼550 to 1900 cm(-1) using the FELIX free electron laser and from ∼2800 to 3800 cm(-1) using an optical parametric oscillator/amplifier (OPO/OPA) laser system. The measured IRMPD spectra of [dGuo+H](+) and [Guo+H](+) are compared to each other and to B3LYP/6-311+G(d,p) linear IR spectra predicted for the stable low-energy conformations computed for these species to determine the most favorable site of protonation, identify the structures accessed in the experiments, and elucidate the influence of the 2'-hydroxyl substituent on the structures and the IRMPD spectral features. Theoretical energetics and the measured IRMPD spectra find that N7 protonation is preferred for both [dGuo+H](+) and [Guo+H](+), whereas O6 and N3 protonated conformers are found to be much less stable. The 2'-hydroxyl substituent does not exert a significant influence on the structures and relative stabilities of the stable low-energy conformations of [dGuo+H](+) versus [Guo+H](+) but does provide additional opportunities for hydrogen bonding such that more low-energy structures are found for [Guo+H](+). [dGuo+H](+) and [Guo+H](+) share very parallel IRMPD spectral features in the FELIX and OPO regions, whereas the effect of the 2'-hydroxyl substituent is primarily seen in the relative intensities of the measured IR bands. The measured OPO/OPA spectral signatures, primarily reflecting the IR features associated with the O-H and N-H stretches, provide complementary information to that of the FELIX region and enable the conformers that arise from different protonation sites to be more readily distinguished. Insight gained from this and parallel studies of other DNA and RNA nucleosides and nucleotides should help better elucidate the effects of the local environment on the overall structures of DNA and RNA.