Canonical Tautomer Research Articles

CanSAR chemistry registration and standardization pipeline

BackgroundIntegration of medicinal chemistry data from numerous public resources is an increasingly important part of academic drug discovery and translational research because it can bring a wealth of important knowledge related to compounds in one place. However, different data sources can report the same or related compounds in various forms (e.g., tautomers, racemates, etc.), thus highlighting the need of organising related compounds in hierarchies that alert the user on important bioactivity data that may be relevant. To generate these compound hierarchies, we have developed and implemented canSARchem, a new compound registration and standardization pipeline as part of the canSAR public knowledgebase. canSARchem builds on previously developed ChEMBL and PubChem pipelines and is developed using KNIME. We describe the pipeline which we make publicly available, and we provide examples on the strengths and limitations of the use of hierarchies for bioactivity data exploration. Finally, we identify canonicalization enrichment in FDA-approved drugs, illustrating the benefits of our approach.ResultsWe created a chemical registration and standardization pipeline in KNIME and made it freely available to the research community. The pipeline consists of five steps to register the compounds and create the compounds’ hierarchy: 1. Structure checker, 2. Standardization, 3. Generation of canonical tautomers and representative structures, 4. Salt strip, and 5. Generation of abstract structure to generate the compound hierarchy. Unlike ChEMBL’s RDKit pipeline, we carry out compound canonicalization ahead of getting the parent structure, similar to PubChem’s OpenEye pipeline. canSARchem has a lower rejection rate compared to both PubChem and ChEMBL. We use our pipeline to assess the impact of grouping the compounds in hierarchies for bioactivity data exploration. We find that FDA-approved drugs show statistically significant sensitivity to canonicalization compared to the majority of bioactive compounds which demonstrates the importance of this step.ConclusionsWe use canSARchem to standardize all the compounds uploaded in canSAR (> 3 million) enabling efficient data integration and the rapid identification of alternative compound forms with useful bioactivity data. Comparison with PubChem and ChEMBL pipelines evidenced comparable performances in compound standardization, but only PubChem and canSAR canonicalize tautomers and canSAR has a slightly lower rejection rate. Our results highlight the importance of compound hierarchies for bioactivity data exploration. We make canSARchem available under a Creative Commons Attribution-ShareAlike 4.0 International License (CC BY-SA 4.0) at https://gitlab.icr.ac.uk/cansar-public/compound-registration-pipeline.

Molecular structure of 5-fluorouracil from gas-phase electron diffraction data and quantum-chemical calculations

The accurate equilibrium molecular structure of a canonical tautomer of 5-fluorouracil was determined by electron diffraction taking into account anharmonic vibrational corrections calculated with ab initio force field. This structure was applied to benchmark the results of CCSD(T) computations used for the analysis of substitution effects in uracil derivatives.

Ground and Excited States of Gas-Phase DNA Nucleobase Cation-Radicals. A UV-vis Photodisociation Action Spectroscopy and Computational Study of Adenine and 9-Methyladenine.

Cation radicals of adenine (A•+) and 9-methyladenine (MA•+) were generated in the gas phase by collision-induced intramolecular electron transfer in copper-terpyridine-nucleobase ternary complexes and characterized by collision-induced dissociation (CID) mass spectra and UV-vis photodissociation action spectroscopy in the 210-700 nm wavelength region. The action spectra of both A•+ and MA•+ displayed characteristic absorption bands in the near-UV and visible regions. Another tautomer of A•+ was generated as a minor product by multistep CID of protonated 9-(2-bromoethyl)adenine. Structure analysis by density functional theory and coupled-clusters ab initio calculations pointed to the canonical 9-H-tautomer Ad1•+ as the global energy minimum of adenine cation radicals. The canonical tautomer MA1•+ was also calculated to be a low-energy structure among methyladenine cation radicals. However, two new noncanonical tautomers were found to be energetically comparable to MA1•+. Vibronic absorption spectra were calculated for several tautomers of A•+ and MA•+ and benchmarked on equation-of-motion coupled-clusters excited-state calculations. Analysis of the vibronic absorption spectra of A•+ tautomers pointed to the canonical tautomer Ad1•+ as providing the best match with the action spectrum. Likewise, the canonical tautomer MA1•+ was the unequivocal best match for the MA•+ ion generated in the gas phase. According to potential-energy mapping, MA1•+ was separated from energetically favorable noncanonical cation radicals by a high-energy barrier that was calculated to be above the dissociation threshold for loss of a methyl hydrogen atom, thus preventing isomerization. Structures and energetics of all four DNA nucleobase cation radicals are compared and discussed.

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.

Proton Transfer in Guanine-Cytosine Base Pairs in B-DNA.

A double proton transfer reaction in a guanine-cytosine (GC) base pair has been proposed as a possible mechanism for rare tautomer (G*C*) formation and thus a source of spontaneous mutations. We analyze this system with free energy calculations based on extensive Quantum Mechanics/Molecular Mechanics simulations to properly consider the influence of the DNA biomolecular environment. We find that, although the G*C* rare tautomer is metastable in the gas phase, it is completely unstable in the conditions found in cells. Thus, our calculations show that a double proton reaction cannot be the source of spontaneous point mutations. We have also analyzed the intrabase H transfer reactions in guanine. Our results show that the DNA environment gives rise to a large free energy difference between the rare and canonical tautomers. These results show the key role of the DNA biological environment for the stability of the genetic code.

The Role of Proton Transfer on Mutations.

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid–base characteristics of the nucleobases in Watson–Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs–DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Dinuclear Metal-Mediated Guanine–Uracil Base Pairs: Theoretical Studies of GUM22+ (M = Cu, Ag, and Au) Ions

Dinuclear metal-mediated hetero base pairs with the d10–d10 closed-shell interactions have significant stability. It is interesting to identify whether coinage metal-mediated Wobble base pairs are also stable. Geometric and electronic structures of the lowest-lying isomers of GUM22+ (G = guanine, U = uracil, M = Cu, Ag, and Au) cluster ions were investigated with density functional theory. In the lowest-lying isomers of these dinuclear metal-mediated base pairs, the 2-oxo-4-hydroxy-trans-N1H isomer of uracil is derived from the canonical tautomer of uracil by the hydrogen atom transfer. M22+ cations remain as an unbroken unit and interact with the G···U ligand through two sets of closely linear N···M···O units, while the reciprocal hydrogen bonds between the Wobble base pair (G–U) are entirely substituted by the M–N or M–O interactions in these complexes. The atoms in molecules and EDA–NOCV calculations really reveal that the σ interactions in GUM22+ cations are the paramount term of ΔEOrb. The obtained instantaneous interaction energies ΔEint and bond dissociation energies of the metal–ligand interactions give the trend of the bond strength as Cu > Au > Ag.

Ammoniated Complexes of Uracil and Transition Metal Ions: Structures of [M(Ura-H)(Ura)(NH3)]+ by IRMPD Spectroscopy and Computational Methods (M = Fe, Co, Ni, Cu, Zn, Cd).

The structures of deprotonated d-block metal dication bound uracil dimers, solvated by a single ammonia molecule, were explored in the gas phase using infrared multiple photon dissociation (IRMPD) spectroscopy in a Fourier transform ion cyclotron resonance-mass spectrometer. The IRMPD spectra were then compared with computed IR spectra for various isomers. Calculations were performed using B3LYP with the 6-31+G(d,p) basis set for all atoms, with the exception of Cd, for which the LANL2DZ basis set with relativistic core potentials was used. The calculations were then repeated using the def2-TZVPP basis set on all atoms and were compared to the first set of calculations. The lowest-energy structures are those in which one uracil is deprotonated at the N3 position and, aside from the Cu complex, the intact uracil is a tautomer in which the N3 hydrogen is at the O4 carbonyl oxygen. The metal displays a tetradentate interaction to the uracil moieties, with the exception of Cu, which is tridentate, and the ammonia molecule is bound directly to the metal center. In the Cu complex, a square planar geometry is observed about the metal center, consistent with Jahn-Teller distortions commonly observed in Cu(II) complexes, and the intact uracil assumes its canonical tautomer. All other metal cation complexes are five-coordinate, square pyramidal complexes, with the intact uracil adopting a tautomer in which the N3 hydrogen is on O4. The IRMPD spectroscopic data are consistent with the computed infrared spectra for the lowest-energy structures in all cases.

Gas-phase quasi-degeneracy of zwitterionic and canonical tautomers of glycine and proline induced by the presence of the MAlF4 (M = Li, Na, K) salts

ABSTRACTThe ionic MAlF4 salts (M = Li, Na, K) are postulated to act as effective stabilising agents with respect to glycine and proline zwitterionic tautomers. On the basis of ab-initio MP2/6-311++G(d,p) theoretical calculations, it is demonstrated that strong electrostatic interactions between the polar MAlF4 moiety and either Gly and Pro zwitterion should allow to suppress the intra-molecular proton transfer in the amino acid and thus render the zwitterion tautomer stable in the gas phase. Among the systems proposed, the K+AlF4– turned out to be the most efficient stabilising compound as its presence is expected to lead to the energetic quasi-degeneracy of the canonical and zwitterionic tautomers of glycine and proline.

Structure, stability, energy barrier and ionization energies of chemically modified DNA-bases: Quantum chemical calculations on 37 favored and rare tautomeric forms of tetraphosphoadenine

Density Functional Theory (BMK and B3LYP) and ab initio complete basis set (CBS-QB3) methods are used to study the stability order and tautomerization processes of all possible 37 tetraphosphoadenine tautomers in the gas phase for the first time. Adiabatic ionization energy (IE) and electron affinity (EA) of all tetraphosphoadenine tautomers were computed at B3LYP/6–311+G(d,p) level of theory at 0K. Positive calculated EA values suggest that the tetraphosphoadenine molecule is an attractive candidate for electron capture. It has been found that all tautomers of tetraphosphoadenine exists within the energy range of 0–150kJmol−1. Similar to adenine, canonical tautomers PA7 and PA9 of tetraphosphoadenine are shown to be the most stable tautomeric forms as predicted by all of the above computational methods. The transition states for the tetraphosphadenine tautomerization processes were investigated for selected systems at B3LYP/6–311+G(d,p) level of theory in the gas phase. Energy barrier for the proton transfer processes of tetraphosphadenine, is up to 97kJmol−1. Dipole moment, full atomic natural charges and bond orders have been computed and discussed.

Gas-phase interactions between lead(II) ions and cytosine: tandem mass spectrometry and infrared multiple-photon dissociation spectroscopy study.

Gas-phase interactions between Pb(2+) ions and cytosine (C) were studied by combining tandem mass spectrometry, infrared multiple photon dissociation spectroscopy, and density functional theory (DFT) calculations. Both singly and doubly charged complexes were generated by electrospray. The [Pb(C)-H](+) complex was extensively studied, and this study shows that two structures, involving the interaction of the metal with the deprotonated canonical keto-amino tautomer of cytosine, are generated in the gas phase; the prominent structure is the bidentate form involving both the N1 and O2 electronegative centers. The DFT study also points out a significant charge transfer from the nucleobase to the low-lying p orbitals of the metal and a strong polarization of the base upon complexation. The various potential energy surfaces explored to account for the fragmentation observed are consistent with the high abundance of the [PbNH2](+) fragment ion.

Infrared multiple photon dissociation action spectroscopy of sodium cationized halouracils: Effects of sodium cationization and halogenation on gas-phase conformation

The gas-phase structures of sodium cationized complexes of 5- and 6-halo-substituted uracils are examined via infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical electronic structure calculations. The halouracils examined in this investigation include: 5-flourouracil, 5-chlorouracil, 5-bromouracil, 5-iodouracil, and 6-chlorouracil. Experimental IRMPD action spectra of the sodium cationized halouracil complexes are measured using a 4.7T Fourier transform ion cyclotron resonance mass spectrometer coupled to the FELIX free electron laser (FEL). Irradiation of the mass selected sodium cationized halouracil complexes by the FEL was carried out over the range of frequencies extending from 950 to 1900cm−1. Theoretical linear IR spectra predicted for the stable low-energy conformations of the sodium cationized halouracils, calculated at B3LYP/6-31G(d) level of theory, are compared with the measured IRMPD action spectra to identify the structures accessed in the experiments. Relative stabilities of the low-energy conformations are determined from single-point energy calculations performed at the B3LYP/6-311+G(2d,2p) level of theory. The evolution of IRMPD spectral features as a function of the size (F, Cl, Br, and I) and position (5 versus 6) of the halogen substituent are examined to elucidate the effects of the halogen substituent and noncovalent interactions with sodium cations on the structure of the nucleobase. Present results are compared with results from energy-resolved collision-induced dissociation and IRMPD action spectroscopy studies previously reported for the protonated and sodium cationized forms of uracil, and halo-, methyl-, and thioketo-substituted uracils. The present results suggest that only a single conformer is accessed for all of the 5-halouracil complexes, whereas multiple conformers are accessed for the Na+(6ClU) complex. In all cases, the experimental IRMPD action spectra confirm that the sodium cation binds to the O4 carbonyl oxygen atom of the canonical diketo tautomer in the ground-state conformers, and gains additional stabilization via chelation interactions with the halogen substituent in the complexes to the 5-halouracils as predicted by theory.

Does the G·G*syn DNA mismatch containing canonical and rare tautomers of the guanine tautomerise through the DPT? A QM/QTAIM microstructural study

We have established that the asynchronous concerted double proton transfer (DPT), moving with a time gap and without stable intermediates, is the underlying mechanism for the tautomerisation of the G·G*syn DNA base mispair (C1 symmetry), formed by the keto and enol tautomers of the guanine in the anti- and syn-configurations, into the G*·G*syn base mispair (C1), formed by the enol and imino tautomers of the G base, using quantum-mechanical calculations and Bader's quantum theory of atoms in molecules. By constructing the sweeps of the geometric, electron-topological, energetic, polar and natural bond orbital properties along the intrinsic reaction coordinate of the G·G*syn↔G*·G*syn DPT tautomerisation, the nine key points, that are critical for the atomistic understanding of the tautomerisation reaction, were set and comprehensively analysed. It was found that the G·G*syn mismatch possesses pairing scheme with the formation of the O6···HO6 (7.01) and N1H···N7 (6.77) H-bonds, whereas the G*·G*syn mismatch – of the O6H···O6 (10.68) and N1···HN7 (9.59 kcal mol−1) H-bonds. Our results highlight that these H-bonds are significantly cooperative and mutually reinforce each other in both mismatches. The deformation energy necessary to apply for the G·G*syn base mispair to acquire the Watson–Crick sizes has been calculated. We have shown that the thermodynamically stable G*·G*syn base mispair is dynamically unstable structure with a lifetime of 4.1 × 10−15 s and any of its six low-lying intermolecular vibrations can develop during this period of time. These data exclude the possibility to change the tautomeric status of the bases under the dissociation of the G·G*syn mispair into the monomers during DNA replication. Finally, it has been made an attempt to draw from the physico-chemical properties of all four incorrect purine–purine DNA base pairs a general conclusion, which claims the role of the transversions in spontaneous point mutagenesis.

Excited states of protonated DNA/RNA bases.

The very fast relaxation of the excited states to the ground state in DNA/RNA bases is a necessary process to ensure the photostability of DNA and its rate is highly sensitive to the tautomeric form of the bases. Protonation of the bases plays a crucial role in many biochemical and mutagenic processes and it can result in alternative tautomeric structures, thus making important the knowledge of the properties of protonated DNA/RNA bases. We report here the photofragmentation spectra of the five protonated DNA/RNA bases. In most of the cases, the spectra exhibit well resolved vibrational structures, with broad bands associated with very short excited state lifetimes. The similarity between the electronic properties, e.g. excitation energy and very short excited state lifetimes for the canonical tautomers of protonated and neutral DNA bases, suggests that the former could also play an important role in the photostability mechanism of DNA.

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

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

Performance of M06, M06-2X, and M06-HF Density Functionals for Conformationally Flexible Anionic Clusters: M06 Functionals Perform Better than B3LYP for a Model System with Dispersion and Ionic Hydrogen-Bonding Interactions

We present a comparative assessment of the performance of the M06 suite of density functionals (M06, M06-2X, and M06-HF) against an MP2 benchmark for calculating the relative energies and geometric structures of the Cl(-)·arginine and Br(-)·arginine halide ion-amino acid clusters. Additional results are presented for the popular B3LYP density functional. The Cl(-)·arginine and Br(-)·arginine complexes are important prototypes for the phenomenon of anion-induced zwitterion formation. Results are presented for the canonical (noncharge separated) and zwitterionic (charge separated) tautomers of the clusters, as well as the numerous conformational isomers of the clusters. We find that all of the M06 functions perform well in terms of predicting the general trends in the conformer relative energies and identifying the global minimum conformer. This is in contrast to the B3LYP functional, which performed significantly less well for the canonical tautomers of the clusters where dispersion interactions contribute more significantly to the conformer energetics. We find that the M06 functional gave the lowest mean unsigned error for the relative energies of the canonical conformers (2.10 and 2.36 kJ/mol for Br(-)·arginine and Cl(-)·arginine), while M06-2X gave the lowest mean unsigned error for the zwitterionic conformers (0.85 and 1.23 kJ/mol for Br(-)·arginine and Cl(-)·arginine), thus providing insight into the types of physical systems where each of these functionals should perform best.

Complexation of anions to gas-phase amino acids: Conformation is critical in determining if the global minimum is canonical or zwitterionic

We report the first MP2 computational study of the Cl−·Arginine and Br−·Arginine complexes, with the clusters being systematically classified into conformational families prior to optimization. The canonical tautomers in particular adopt numerous conformeric structures with close relative energies. Conformer energy differences are of similar magnitude to the canonical–zwitterionic tautomer differences, indicating that it is essential to obtain reliable conformer energies to determine whether the global minimum structure is zwitterionic or canonical for these prototype anionic systems. The zwitterionic and canonical minima are found to be very close in energy for Br−·Arginine, whereas the canonical form is significantly lower for Cl−·Arginine.

Discovery of Most Stable Structures of Neutral and Anionic Phenylalanine through Automated Scanning of Tautomeric and Conformational Spaces.

We have developed a software tool for combinatorial generation of tautomers and conformers of small molecules. We have demonstrated it by performing a systematic search for the most stable structures of neutral and anionic phenylalanine (Phe) using electronic structure methods. For the neutral canonical tautomer we found out that the conformers with and without the intramolecular (O)H···NH2 hydrogen bond are similarly stable, within the error bars of our method. A unique IR signature of the conformer without the hydrogen bond has been identified. We also considered anions of Phe, both valence type and dipole-bound. We have found out that tautomers resulting from proton transfer from the carboxylic OH to the phenyl ring do support valence anions that are vertically strongly bound, with electron vertical detachment energies (VDE) in a range of 3.2-3.5 eV. The most stable conformer of these valence anions remains adiabatically unbound with respect to the canonical neutral by only 2.17 kcal/mol at the CCSD(T)/aug-cc-pVDZ level. On the basis of our past experience with valence anions of nucleic acid bases, we suggest that the valence anions of Phe identified in this report can be observed experimentally. The most stable conformer of canonical Phe is characterized by an adiabatic electron affinity of 53 meV (a dipole-bound state).

On relation between prototropy and electron delocalization for neutral and redox adenine – DFT studies

Electron delocalization, which plays a principal role for the tautomeric adenine system, was studied for all possible major, minor, and rare tautomers of neutral, oxidized, and reduced adenine. The recently extended geometry-based HOMED (harmonic oscillator model of electron delocalization) index, which measures any type of resonance conjugation (π–π, n–π, and σ–π), was applied to the geometries of adenine isomers optimized at the DFT(B3LYP)/6-311+G(d,p) level. For neutral adenine, variations of the HOMED indices estimated for the whole tautomeric systems (eleven bonds) are almost parallel to those of the relative Gibbs energies (ΔG), which measure the thermodynamic stabilities of individual isomers. Some subtle effects, being a consequence of intramolecular interactions of neighboring groups and variability of electron delocalization for the imidazole and pyrimidine fragments, perturb this relation. One-electron oxidation has slight effect on the general HOMED/ΔG relation. The favored canonical NHNH tautomer is well delocalized for both neutral and oxidized adenine. One-electron reduction dramatically changes electron delocalization and thermodynamic stability of individual tautomers that two subfamilies can be distinguished on the HOMED vs ΔG plot, one for the NHNH tautomers and the other one for the NHCH isomers. The reduced canonical NHNH tautomer (although well delocalized) is not the favored one. The non-aromatic NHCH tautomer predominates in the tautomeric mixture, indicating that aromaticity is not the main factor that dictates the tautomeric preference for reduced adenine.

Photochemistry and photophysics of the amino and imino tautomers of 1-methylcytosine: tautomerisation as a side product of the radiationless decay

The photochemistry and photophysics of the amino-oxo (AO) and imino-oxo (IO) tautomers of 1-methylcytosine are investigated with ab initio calculations, using the CASPT2//CASSCF approach. Our aim is to rationalize recent experimental results which show that the AO tautomer has an excited state lifetime of about 1 ps, similar to cytosine (J.-W. Ho, H.-C. Yen, W.-K. Chou, C.-N. Weng, L.-H. Cheng, H.-Q. Shi, S.-H. Lai and P.-Y. Cheng, J. Phys. Chem. A, 2011, 115, 8406-8418), and that irradiation with wavelengths shorter than 308 nm induces the AO→IO tautomerization (I. Reva, M. J. Nowak, L. Lapinski and R. Fausto, J. Phys. Chem. B, 2012, 116, 5703-5710). For the canonical AO tautomer we find two analogous decay mechanisms to those described previously for cytosine, involving a conical intersection of ethylenic type and one where the amino group is bent out of plane. Decay through these intersections provides an unreactive return path to the AO ground state species. More importantly, we have identified new decay paths that lead from the two intersections to the trans-IO tautomer, without a barrier. These paths provide the possibility of forming the IO tautomer, presumably in small yields, as a side product of the radiationless decay. Thus, we have established for the first time computationally the mechanism of the UV-induced tautomerization, which is compatible with the well-established decay mechanism for cytosine. For the IO tautomer, we also find a mechanism for the excited state interconversion of the cis and trans forms and for efficient radiationless decay through a conical intersection where the imino group is perpendicular to the ring and bent out of the plane. These results are likely to be valid also for cytosine.