Theoretical study on the excitation energies of six tautomers of guanine: evidence for the assignment of the rare tautomers

The CASPT2//CASSCF method with the 6-31G basis set and an active space up to (16,12) was used to calculate the excitation energies for six tautomers of guanine. Our calculations provide further support on the recent reassignment of the near-UV resonant two-photon ionization (R2PI) spectrum, in which two rare tautomers of the 7H-oxo-imino form were proposed to replace the previously assigned 7H/9H-oxo-amino tautomers. The adiabatic excitation energies of the 7H-oxo-imino tautomers are calculated to be 0.3-0.5 eV higher than those of the 7H/9H-oxo-amino tautomers. Our calculations also indicate that the missing most stable tautomers (7H/9H-oxo-amino tautomers) in the R2PI experiment is possibly due to the existence of an ultrafast nonradiative deactivation process in the excited-state of these two tautomers.     

An ab initio study of excited states of guanine in the gas phase and aqueous media: Electronic transitions and mechanism of spectral oscillations

Ground state geometries of the four tautomeric forms keto‐N9H, keto‐N7H, enol‐N9H, and enol‐N7H of guanine were optimized in the gas phase at the RHF level using a mixed basis set consisting of the 4‐31G basis set for all the atoms except the nitrogen atom of the amino group for which the 6‐311+G* basis set was used. These calculations were also extended to hydrogen‐bonded complexes of three water molecules with each of the keto‐N9H (G9‐3W) and keto‐N7H (G7‐3W) forms of guanine. Relative stabilities of the four above‐mentioned tautomers of guanine as well as those of G9‐3W and G7‐3W complexes in the ground state in the gas phase were studied employing the MP2 correlation correction. In aqueous solution, relative stabilities of these systems were studied using the MP2 correlation correction and polarized continuum model (PCM) or the isodensity surface polarized continuum model (IPCM) of the self‐consistent reaction field (SCRF) theory. Geometry optimization in the gas phase at the RHF level using the 6‐31+G* basis set for all atoms and the solvation calculations in water at the MP2 level using the same basis set were also carried out for the nonplanar keto‐N9H and keto‐N7H forms of guanine. Thus, it is shown that among the different tautomers of guanine, the keto‐N7H form is most stable in the gas phase, while the keto‐N9H form is most stable in aqueous solution. It appears that both the keto‐N9H and keto‐N7H forms of guanine would be present in the ground state, particularly near the aqueous solution–air interface. Vertical excitation and excited state geometry optimization calculations were performed using configuration interaction involving single electron excitation (CIS). It is found that the absorption spectrum of guanine would arise mainly due to its keto‐N9H form but the keto‐N7H form of the same would also make some contribution to it. The enol‐N9H and enol‐N7H forms of the molecule are not expected to occur in appreciable abundance in the gas phase or aqueous media. The normal fluorescence spectrum of guanine in aqueous solution with a peak near 332 nm seems to originate from the lowest singlet excited state of the keto‐N7H form of the molecule while the fluorescence of oxygen‐rich aqueous solutions of guanine with a peak near 450 nm appears to originate from the lowest singlet excited state of the keto‐N9H form of the molecule. The origin of the slow damped spectral oscillation observed in the absorption spectrum of guanine has been explained. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 826–846, 2000     

Valence and dipole-bound anions of the most stable tautomers of guanine

Anionic states of guanine, which is the only nucleic acid base of which the anions have not yet been studied in either photoelectron spectroscopic (PES) or Rydberg electron transfer (RET) experiments, have been characterized for the four most stable tautomers of neutral guanine using a broad spectrum of electronic structure methods from the density functional theory, with the B3LYP exchange-correlation functional, to the coupled-cluster method, with single, double, and perturbative triple excitations. Both valence and dipole-bound anionic states were addressed. We identified some of the difficulties facing future PES or RET experiments on the anion of guanine. Even if guanine is successfully transferred to the gas phase without thermal decomposition, it is critical to have the canonical amino-oxo (G) and both amino-hydroxy (GH and GHN7H) tautomers in the beam, not only the most stable, a noncanonical, amino-oxo tautomer (GN7H), as the latter does not support an adiabatically bound anionic state. We also suggested a scheme for enrichment of gas-phase guanine with the canonical tautomer, which is not the most stable in the gas phase, but which is of main interest due to its biological relevance. The tautomers G, GN7H, and GHN7H support vertically bound valence anionic states with the CCSD(T) value of vertical detachment energy of +0.58, +0.21, and +0.39 eV, respectively. These anionic states are, however, adiabatically unbound and thus metastable. The vertical electronic stability of these valence anionic states is accompanied by serious "buckling" of the molecular skeleton. The G and GHN7H tautomers support dipole-bound states with the CCSD(T) values of adiabatic electron affinity of 65 and 36 meV, respectively. A contribution from higher-than-second-order correlation terms represents, respectively, 48 and 68% of the total vertical electron detachment energy determined at the CCSD(T) level.     

Excited-state ab initio calculations and multidimensional Franck-Condon simulations on guanine

The guanine enol and keto N7H and N9H tautomers have been optimized at the CASSCF/cc-pVDZ levels of theory. Except for the enol N7H tautomer, CASSCF predicts distorted nonplanar S1 state geometries. Among the vibronic simulations carried out with the optimized structures only the enol N7H tautomer qualitatively mirrors the appearance of the experimental R2PI spectrum. Refined symmetry-adapted cluster configuration interaction (SACCI) geometries of the enol N7H tautomer produce simulations in good agreement with experiment and support the assignment of the first vibronic band and associated vibronic features of the R2PI spectrum to this tautomer. The sharp spectral features and the fact that Franck-Condon simulations based on the harmonic approximation allow for a faithful reproduction of the spectral signature associated with the enol N7H tautomer indicate that within the simulated energy window the S1 potential energy surface of this isomer is fairly harmonic and free from conical intersections involved in the S1 state lifetime-shortening relaxation processes of other DNA bases and possibly the remaining tautomers of guanine.     

Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.     

Quantum chemical investigation of the electronic spectra of the keto, enol, and keto-imine tautomers of cytosine

The low-lying excited singlet states of the keto, enol, and keto-imine tautomers of cytosine have been investigated employing a combined density functional/multireference configuration interaction (DFT/MRCI) method. Unconstrained geometry optimizations have yielded out-of-plain distorted structures of the pi --> pi and n --> pi excited states of all cytosine forms. For the keto tautomer, the DFT/MRCI adiabatic excitation energy of the pi --> pi state (4.06 eV including zero-point vibrational energy corrections) supports the resonant two-photon ionization (R2PI) spectrum (Nir et al. Phys. Chem. Chem. Phys. 2002, 5, 4780). On its S1 potential energy surface, a conical intersection between the 1pipi state and the electronic ground state has been identified. The barrier height of the reaction along a constrained minimum energy path amounts to merely 0.2 eV above the origin and explains the break-off of the R2PI spectrum. The 1pipi minimum of the enol tautomer is found at considerably higher excitation energies (4.50 eV). Because of significant geometry shifts with respect to the ground state, long vibrational progressions are expected, in accord with experimental observations. For the keto-imine tautomer, a crossing of the 1pipi potential energy surface with the ground-state surface has been found, too. Its n --> pi minimum (3.27 eV) is located well below the conical intersection between the pi --> pi and S0 states, but it will be difficult to observe because of its small transition moment. The identified conical intersections of the pi --> pi excited states of the keto cytosine tautomers are made responsible for the ultrafast decay to the electronic ground states and thus may explain their subpicoseconds lifetimes.     

Molecular beam resonant two-photon ionization study of caffeine and its hydrated clusters

We investigated electronically excited states of caffeine and its 1:1 complex with water by using resonant two-photon ionization (R2PI) and UV-UV hole-burning techniques. Strong vibronic coupling between a pair of close-lying pi-pi* and n-pi* transitions is proposed to be responsible for the broad spectral feature observed. By comparing the experimental results with those of theoretical calculations, both the O-bonded and N-bonded forms were suggested to be initially produced for the 1:1 complex. The electronic transitions of the O-bonded complex were blueshifted in the R2PI spectrum. For the N-bonded complex, the excited state undergoes an ultrafast decay process, followed by dissociation on a repulsive potential energy surface, which gives rise to a characteristically anomalous cluster distribution in nanosecond experiments.     

On the accessibility to conical intersections in purines: hypoxanthine and its singly protonated and deprotonated forms

The dynamics following electronic excitation of hypoxanthine and its nucleoside inosine were studied by femtosecond fluorescence up-conversion. Our objective was to explore variants of the purinic DNA bases in order to determine the molecular parameters that increase or reduce the accessibility to ground state conical intersections. From experiments in water and methanol solution we conclude that both dominant neutral tautomers of hypoxanthine exhibit ultrashort excited state lifetimes (τ < 0.2 ps), which are significantly shorter than in the related nucleobase guanine. This points to a more accessible conical intersection for the fluorescent state upon removal of the amino group, present in guanine but absent in hypoxanthine. The excited state dynamics of singly protonated hypoxanthine were also studied, showing biexponential decays with a 1.1 ps component (5%) besides a sub-0.2 ps ultrafast component. On the other hand, the S(1) lifetimes of the singly deprotonated forms of hypoxanthine and inosine show drastic differences, where the latter remains ultrafast but the singly deprotonated hypoxanthine shows a much longer lifetime of 19 ps. This significant variation is related to the different deprotonation sites in hypoxanthine versus inosine, which gives rise to significantly different resonance structures. In our study we also include multireference perturbation theory (MRMP2) excited state calculations in order to determine the nature of the initial electronic excitation in our experiments and clarify the ordering of the states in the singlet manifold at the ground state geometry. In addition, we performed multireference configuration interaction calculations (MR-CIS) that identify the presence of low-lying conical intersections for both prominent neutral tautomers of hypoxanthine. In both cases, the surface crossings occur at geometries reached by out of plane opposite motions of C2 and N3. The study of this simpler purine gives several insights into how small structural modifications, including amino substitution and protonation site and state, determine the accessibility to conical intersections in this kind of heterocycles.     

Structure and singly occupied molecular orbital analysis of anionic tautomers of guanine

Recently, we reported the discovery of adiabatically bound anions of guanine that might be involved in the processes of DNA damage by low-energy electrons and in charge transfer through DNA. These anions correspond to some tautomers that have been ignored thus far. They were identified using a hybrid quantum mechanical-combinatorial approach in which an energy-based screening was performed on the library of 499 tautomers with their relative energies calculated with quantum chemistry methods. In the current study, we analyze the adiabatically bound anions of guanine in two aspects: (1) the geometries and excess electron distributions are analyzed and compared with anions of the most stable neutrals to identify the sources of stability; (2) the chemical space of guanine tautomers is explored to verify if these new tautomers are contained in a particular subspace of the tautomeric space. The first task involves the development of novel approaches-the quantum chemical data like electron density, orbital, and information on its bonding/antibonding character are coded into holograms and analyzed using chemoinformatics techniques. The second task is completed using substructure analysis and clustering techniques performed on molecules represented by 2D fingerprints. The major conclusion is that the high stability of adiabatically bound anions originates from the bonding character of the pi orbital occupied by the excess electron. This compensates for the antibonding character that usually causes significant buckling of the ring. Also, the excess electron is more homogenously distributed over both rings than in the case of anions of the most stable neutral species. In terms of 2D substructure, the most stable anionic tautomers generally have additional hydrogen atoms at C8 and/or C2 and they do not have hydrogen atoms attached to C4, C5, and C6. They also form an "island of stability" in the tautomeric space of guanine.     

Na+, Mg2+, and Zn2+ binding to all tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine: a correlated ab initio quantum chemical study

Interactions of adenine, cytosine, guanine, and thymine with Na(+), Mg(2+), and Zn(2+) cations were studied using an approximate resolution of identity correlated second-order MP2 (RI-MP2) method with the TZVPP ([5s3p2d1f/3s2p1d]) basis set. All existing tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine were considered. Cations bind mostly in a bidentate manner, and stabilization energies of these complexes are larger than those in the case when cations bind in a unidentate manner. The cation...Y (Y equal to N or O) distances for divalent metals are shorter than those for Na(+) and for Zn(2+) are mostly shorter than the Mg(2+)...Y distance. The intermolecular distances between the cation and the base for complexes containing adenine and cytosine are systematically shorter than those for complexes containing guanine and thymine. Only for cytosine the canonical keto/amino tautomer structure with ions represents the global minimum. For guanine, the metalated canonical form is again the most stable, but its stabilization energy is within less than 5% of the stabilization energies of the two other rare tautomers, which indicates that the canonical form and these two rare tautomers could coexist. The canonical structures of adenine and thymine in the presence of ions are considerably less stable (by more than 10%) than the complexes of the rare tautomers. It can be concluded that the interaction of Na(+), Mg(2+), and Zn(2+) cations with cytosine in the gas phase will not induce the change of the canonical form to the rare tautomeric form. In the case of isolated guanine, the equilibrium of the canonical form with rare tautomers can be found. For isolated adenine and thymine the presence of rare tautomers is highly probable.     

生色团连接的苯骈三氮唑衍生物的激发态分子内质子转移

用从头算和密度泛函理论研究了对硝基二苯乙烯作为生色团连接的2-(2-羟基-苯基)-苯骈三氮唑的衍生物2-羟基-5-[对硝基-二苯乙烯基-氧亚甲基]-苯基-(2H-苯骈三氮唑)(C1)和4′-硝基-3,4-二[2-羟基-(2H-苯骈三氮唑)-苄氧基]-二苯乙烯(C2)发生激发态分子内质子转移(ESIPT)的可能性.系统研究了C1和C2发生ESIPT的互变异构体的基态与激发态的性质变化,包括相关的键长、键角等结构参数,Mulliken电荷和偶极矩,前线轨道以及势能曲线.计算结果表明,对于C1来讲,酮式(keto)的基态(K)不存在稳定结构,因此发生基态分子内质子转移(GSIPT)可能性很小.酮式的激发态(K*)的氢键强度要远强于烯醇式(enol)的激发态(E*)的氢键强度.分子在光致激发后,质子供体所带负电荷减小而质子受体所带负电荷增加.在K*,HOMO→LUMO的电子跃迁导致电子密度从"酚环"向质子化杂环转移.E*→K*跃迁只需要克服较小的能垒(约41 kJ.mol-1).计算结果表明C1发生ESIPT的可能性很大.C2由于具有高能量,其具有基态的单质子转移特征的异构体EK(同时含烯醇E与酮K结构)、具有基态的双质子转移特征的异构体2K(含有双酮结构),以及具有双酮结构特征的激发态2K*均无法获得它们的稳定结构,因此,基态分子内单或双质子转移和激发态分子内双重质子转移发生的可能性极小.然而,由于双烯醇式的激发态(2E*)和EK的激发态(EK*)存在稳定结构,且2E*→EK*跃迁具有低能垒,因此C2有可能发生激发态分子内单重质子转移.本文进一步计算了两个分子的紫外-可见吸收光谱与荧光发射光谱,获得了具有较大斯托克位移的ESIPT的荧光发射峰.     

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

Altogether eight keto and enol tautomers of guanine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in bulk water. The structures of isolated, as well as mono- and dihydrated tautomers were determined by means of the RI-MP2 method using the extended TZVPP (5s3p2d1f/3s2p1d) basis set. The relative energies of isolated tautomers included the correction to higher correlation energy terms evaluated at the CCSD(T)/aug-cc-pVDZ level. The relative enthalpies at 0 K and relative free energies at 298 K were based on the above-mentioned relative energies and zero-point vibration energies, temperature-dependent enthalpy terms and entropies evaluated at the MP2/6-31G level. The keto form having hydrogen atom at N7 is the global minimum while the canonical form having hydrogen atom at N9 represents the first local minimum at all theoretical levels in vacuo and in the presence of 1 and 2 water molecules. All three unusual rare tautomers having hydrogens at N3 and N7, at N3 and N9, and also at N9 and N7 are systematically considerably less stable and can be hardly detected in the gas phase. The theoretical predictions fully agree with existing theoretical as well as experimental results. The effect of bulk solvent on the relative stability of guanine tautomers was studied by self-consistent reaction field and molecular dynamics free energy calculations using the thermodynamic integration method. Bulk solvent, surprisingly, strongly favored these three rare tautomers over all remaining low-energy tautomers and probably only these forms can exist in water phase. The global minimum (tautomer with hydrogens at N3 and N7) is by 13 kcal/mol more stable than the canonical form (3rd local minimum). Addition of one or two water molecules does not change the relative stability order of isolated guanine tautomers but the respective trend clearly supports the surprising stabilization of three rare forms.     

The guanine tautomer puzzle: quantum chemical investigation of ground and excited states

Combined density functional and multireference configuration interaction methods have been employed to explore the ground and low-lying electronically excited states of the most important tautomeric and rotameric forms of guanine with the purpose of resolving the conflicting assignments of IR-UV bands found in the literature. The calculations predict sharp 1(pi-->pi*) origin transitions for the RN1 rotamer of the 7H-amino-hydroxy species and the RN7 rotamer of the 9H-amino-hydroxy species. The other 9H-amino-hydroxy rotamer, RN1, undergoes ultrafast nonradiative decay and is thus missing in the UV spectra. Because of its very small Franck-Condon factor and the presence of a conical intersection close by, it appears questionable, whether the 1(pi-->pi*) origin transition of 9H-amino-oxo-guanine can be observed experimentally. Vibrational overlap is more favorable for the 1(pi-->pi*) origin transition of the 7H- amino-oxo form, but also this tautomer is predicted to undergo ultrafast nonradiative decay of the 1(pi-->pi*) population. The good agreement of calculated IR frequencies of the amino-oxo species with recent IR spectra in He droplets and their mismatch with peaks observed in IR-UV spectra indicate that none of the bands stem from 7H- or 9H-amino-oxo guanine. Instead, our results suggest that these bands originate from 7H-imino-oxo guanine tautomers. In the excited-state dynamics of the biologically relevant 9H-amino-oxo tautomer, a diffuse charge transfer state is predicted to play a significant role.     

Combined experimental/theoretical investigation of the He+ICl interactions. I. Rovibronic spectrum of He...ICl complexes in the ICl B-X, 3-0 region

Transitions of two different stereoisomers of the He...ICl(X,v" = 0) weakly bound complex, one with a T-shaped orientation and another that is most likely linear, have been observed in laser-induced fluorescence experiments performed in the ICl B-X region. Here we present experimental and theoretical results aimed at confirming the previous assignments and at gaining additional insights into the He+ICl interactions. High resolution action spectra were recorded in the same region to identify those features that could be attributed to transitions of the He...I35Cl(X,v" = 0) isomers and not to higher-order complexes, Hen...I35Cl, where n > or = 2, or I37Cl containing species. Calculations of the rovibronic spectra of the He...I35Cl complexes in the ICl B-X, 2-0 and 3-0 regions were performed using an ab initio potential energy surface for the He+ICl(X,v" = 0) ground state and two different pairwise additive potentials for the He+ICl(B,v' = 2,3) excited states. The rotation-vibration energies and wave functions for the He cdots, three dots, centered I35Cl complexes were obtained for all bound states with total angular momentum J < 10 using both of these potentials. Electronic spectra were generated using these results, assuming that the transition moment lies along the ICl bond and is not perturbed by the presence of the helium atom. The calculations qualitatively reproduce the He cdots, three dots, centered I35Cl action spectrum and strongly support the previous assignments. The calculations also indicate that some of the spectral congestion observed near the linear band may be attributed to transitions of the linear isomer to multiple intermolecular levels in the excited state. Coriolis coupling strongly mixes He cdots, three dots, centered ICl(B,v') states with rotational excitation, making simulations and assignments of the linear band observed in the experimental spectrum difficult.     

Tautomerism of acridin-9-amines substituted at the exocyclic nitrogen atom: spectroscopic investigations and theoretical studies

In the ground electronic state, the first two (1 and 2) of the compounds investigated--9-(methoxyamino)acridine (1), 9-hydrazinoacridine (2), N-(2-chloroethyl)acridin-9-amine (3) and N-(5-methylpyridin-2-yl)acridin-9-amine (4)--exist principally in the imino tautomeric form, while the other two (3 and 4) can coexist as amino and imino tautomers in solvents of various polarities and abilities to enter into specific interactions. These features of the molecules are reflected in the experimental absorption spectra and the predicted thermodynamic and spectral data. The predicted thermodynamic characteristics suggest that in the S1 state, the imino tautomers of 1, 2 and 4 and the amino tautomer of 3 are more stable than their tautomeric counterparts. However, the predicted rates of intersystem crossing suggest that the imino tautomers of 1-3 and the amino tautomer of 4 lose excitation energy very rapidly, so that only their counterpart tautomers in fact emit radiation. This explains why 1 and 2 do not fluoresce and why the amino form of 3 and the imino form of 4 are the emitters. 3 and 4 thus represent acridin-9-amines whose imino forms are preferred in the ground state, but whose respective amino and imino forms are preferred in the excited state. Because 3 and 4 are capable of tautomeric transformations in the ground and excited states, and also enter into specific interactions with solvents, these compounds could be potent spectral indicators or probes of environmental properties.     

Tautomers and electronic states of jet-cooled 2-aminopurine investigated by double resonance spectroscopy and theory

We present resonant two-photon ionization (R2PI), IR-UV, and UV-UV double resonance spectra of jet-cooled 2-aminopurine (2AP) as well as Fourier transform infrared (FTIR) gas phase spectra. 2AP is a fluorescing isomer of the nucleobase adenine. The results show that there is only one tautomer of 2AP which absorbs in the wavelength range 32,300-34,500 cm(-1). The comparison with the calculated IR spectra of 9H- and 7H-2AP points to 9H-2AP as the dominating tautomer in the gas phase but the spectra are too similar to allow an unambiguous assignment to the respective tautomer. Hence, we determined vertical and adiabatic excitation energies of both tautomers employing combined density functional theory and multi-reference configuration interaction techniques. For the 0-0 band of the first 1pipi* transition of 9H-2AP we obtain a theoretical value of 32,328 cm(-1), in excellent agreement with the band origin of our R2PI spectrum at 32,371 cm(-1). The first singlet pipi* transition of the less stable 7H-2AP tautomer is predicted to be red-shifted by about 1700 cm(-1) with respect to the corresponding transition in 9H-2AP. From the absence of experimental bands in the energy region between 30,300 and 32,350 cm(-1) we conclude that 7H-2AP is not present to an appreciable extent in the molecular beam. Our calculations yield nearly equal energies for the 1npi* and 1pipi* minima of isolated 2AP, similar to the situation in adenine. The hitherto existing argument that the energetic order of states is responsible for the different spectroscopic properties of these isomers therefore does not hold. Rather, vibronic levels close to the origin of the 1pipi* transition cannot access the conical intersection between the 1pipi* and S(0) states along a puckering coordinate of the six-membered ring, in contrast to the situation in electronically excited 9H-adenine. As a consequence, a rich vibrational structure can be observed in the R2PI spectrum of 2AP whereas the spectrum of 9H-adenine breaks off at low energies.     

Nonradiative decay mechanisms of the biologically relevant tautomer of guanine

We present the excited-state potential energy profiles of the biologically relevant 9H-keto-amino tautomer of guanine with respect to the radiationless decay via the out-of-plane deformation of the six-membered ring as well as the dissociation of NH bonds. The CASPT2//CASSCF method is employed for the reaction-path calculations. The reaction path for the out-of-plane deformation in the (1)pi pi* state leads in a barrierless way to a conical intersection with the electronic ground state. For the NH dissociation via the (1)pi sigma* state, the 9H-keto-amino tautomer is shown to have lower energy barriers than the 7H tautomers which we have studied recently. These two radiationless decay mechanisms explain the unexpected missing of the biologically relevant form in the resonant two-photon ionization spectrum of guanine in a supersonic jet. It is suggested that these ultrafast deactivation processes provide the biologically relevant tautomer of guanine with a high degree of photostability.     

The electronic spectrum of protonated adenine: theory and experiment

In this work we present the results of a combined experimental and theoretical study concerned with the question how a proton changes the electronic spectrum and dynamics of adenine. In the experimental part, isolated adenine ions have been formed by electro-spray ionisation, stored, mass-selected and cooled in a Paul trap and dissociated by resonant photoexcitation with ns UV laser pulses. The S(0)-S1 spectrum of protonated adenine recorded by fragment ion detection lies in a similar energy range as the first pipi* transition of neutral 9H-adenine. It shows a flat onset with a broad substructure, indicating a large S(0)-S1 geometry shift and an ultra-short lifetime. In the theoretical part, relative energies of the ground and the excited states of the most important tautomers have been calculated by means of a combined density functional theory and multi-reference configuration interaction approach. Protonation at the nitrogen in position 1 of the neutral 9H-adenine tautomer yields the most stable protonated adenine species, 1H-9H-A+. The 3H-7H-A+ and the 3H-9H-A+ tautomers, formed by protonation of 7H- and 9H-adenine in 3-position, are higher in energy by 162 cm(-1) and 688 cm(-1), respectively. Other tautomers lie at considerably higher energies. Calculated vertical absorption spectra are reported for all investigated tautomers whereas geometry optimisations of excited states have been carried out only for the most interesting ones. The S1 state energies and geometries are found to depend on the protonation site. The theoretical data match best with the experimental onset of the spectrum for the 1H-9H-A+ tautomer although we cannot definitely exclude contributions to the experimental spectrum from the 3H-7H-A+ tautomer at higher energies. The vertical S(0)--> S1 excitation energy is similar to the one in neutral 9H-adenine. As for the neutral adenine, we find a conical intersection of the S1 of protonated adenine with the ground state in an out-of-plane coordinate but at lower energies and accessible without barrier.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

Imino tautomers of gas-phase guanine from mid-infrared laser spectroscopy

We reinvestigated the assignment of the three major guanine conformers detected via resonance enhanced two-photon ionization (R2PI) in supersonic expansions and present IR/UV double resonance spectra in the spectral region between 1500 and 1800 cm(-1). Comparison with B3LYP/TZVPP and RI-MP2/cc-pVQZ calculations shows that both conformers B and C are 7H-keto tautomers with an imine group in the 2-position. They differ only in the local conformation of the imine group but are otherwise identical. Conformer A is an amino-enol form with the OH group in the trans position.