Solvent effect on the singlet excited-state lifetimes of nucleic acid bases: A computational study of 5-fluorouracil and uracil in acetonitrile and water

The first comprehensive quantum mechanical study of solvent effects on the behavior of the two lowest energy excited states of uracil derivatives is presented. The absorption and emission spectra of uracil and 5-fluorouracil in acetonitrile and aqueous solution have been computed at the time-dependent density-functional theory level, using the polarizable continuum model (PCM) to take into account bulk solvent effects. The computed spectra and the solvent shifts provided by our method are close to their experimental counterpart. The S0/S1 conical intersection, located in the presence of hydrogen-bonded solvent molecules by CASSCF (8/8) calculations, indicates that the mechanism of ground-state recovery, involving out-of-plane motion of the 5 substituent, does not depend on the nature of the solvent. Extensive explorations of the excited-state surfaces in the Franck-Condon (FC) region show that solvent can modulate the accessibility of an additional decay channel, involving a dark n/pi* excited state. This finding provides the first unifying explanation for the experimental trend of 5-fluorouracil excited-state lifetime in different solvents. The microscopic mechanisms underlying solvent effects on the excited-state behavior of nucleobases are discussed.     

Solvent effect on the singlet excited-state dynamics of 5-fluorouracil in acetonitrile as compared with water

The excited-state dynamics of 5-fluorouracil in acetonitrile has been investigated by femtosecond fluorescence upconversion spectroscopy in combination with quantum chemistry TD-DFT calculations ((PCM/TD-PBE0). Experimentally, it was found that when going from water to acetonitrile solution the fluorescence decay of 5FU becomes much faster. The calculations show that this is related to the opening of an additional decay channel in acetonitrile solution since the dark n/pi* excited state becomes near degenerate with the bright pi/pi* state, forming a conical intersection close to the Franck-Condon region. In both solvents, a S1-S0 conical intersection, governed by the out-of-plane motion of the fluorine atom, is active, allowing an ultrafast internal conversion to the ground state.     

Ultrafast deactivation of an excited cytosine-guanine base pair in DNA

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.     

The decay from the dark npi* excited state in uracil: an integrated CASPT2/CASSCF and PCM/TD-DFT study in the gas phase and in water

By integrating the results of MS-CASPT2/CASSCF and TD-PBE0 calculations, we propose a mechanism for the decay of the excited dark state in pyrimidine, fully consistent with all the available experimental results. An effective conical intersection (CI-npi) exists between the spectroscopic pi/pi* excited state (Spi) and a dark n/pi* state (Sn), and a fraction of the population decays to the minimum of Sn (Sn-min). The conical intersection between Sn and the ground-state is not involved in the decay mechanism, because of its high energy gap with respect to Sn-min. On the other hand, especially in hydrogen bonding solvents, the energy gap between Sn-min and CI-npi is rather small. After thermalization in Sn-min, the system can thus recross CI-npi and then quickly proceed on the Spi barrierless path toward the conical intersection with the ground state.     

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.     

Quenching of tryptophan (1)(pi,pi*) fluorescence induced by intramolecular hydrogen abstraction via an aborted decarboxylation mechanism

CASSCF computations show that the hydrogen-transfer-induced fluorescence quenching of the (1)(pi,pi*) excited state of zwitterionic tryptophan occurs in three steps: (1) formation of an intramolecular excited-state complex, (2) hydrogen transfer from the amino acid side chain to the indole chromophore, and (3) radiationless decay through a conical intersection, where the reaction path bifurcates to a photodecarboxylation and a phototautomerization route. We present a general model for fluorescence quenching by hydrogen donors, where the radiationless decay occurs at a conical intersection (real state crossing). At the intersection, the reaction responsible for the quenching is aborted, because the reaction path bifurcates and can proceed forward to the products or backward to the reactants. The position of the intersection along the quenching coordinate depends on the nature of the states and, in turn, affects the formation of photoproducts during the quenching. For a (1)(n,pi*) model system reported earlier (Sinicropi, A.; Pogni, R.; Basosi, R.; Robb, M. A.; Gramlich, G.; Nau, W. M.; Olivucci, M. Angew. Chem., Int. Ed. 2001, 40, 4185-4189), the ground and the excited state of the chromophore are hydrogen acceptors, and the excited-state hydrogen transfer is nonadiabatic and leads directly to the intersection point. There, the hydrogen transfer is aborted, and the reaction can return to the reactant pair or proceed further to the hydrogen-transfer products. In the tryptophan case, the ground state is not a hydrogen acceptor, and the excited-state hydrogen transfer is an adiabatic, sequential proton and electron transfer. The decay to the ground state occurs along a second reaction coordinate associated with decarboxylation of the amino acid side chain and the corresponding aborted conical intersection. The results show that, for (1)(pi,pi*) states, the hydrogen transfer alone is not sufficient to induce the quenching, and explain why fluorescence quenching induced by hydrogen donors is less general for (1)(pi,pi*) than for (1)(n,pi*) states.     

Theoretical study on the singlet excited state of pterin and its deactivation pathway

The excited-state properties and related photophysical processes of the acidic and basic forms of pterin have been investigated by the density functional theory and ab initio methodologies. The solvent effects on the low-lying states have been estimated by the polarized continuum model and combined QM/MM calculations. Calculations reveal that the observed two strong absorptions arise from the strong pi --> pi* transitions to 1(pipi*L(a)) and 1(pipi*L(b)) in the acidic and basic forms of pterin. The first 1(pipi*L(a)) excited state is exclusively responsible for the experimental emission band. The vertical 1(n(N)pi*) state with a small oscillator strength, slightly higher in energy than the 1(pipi*L(a)) state, is less accessible by the direct electronic transition. The 1(n(N)pi*) state may be involved in the photophysical process of the excited pterin via the 1(pipi*L(a)/n(N)pi*) conical intersection. The radiationless decay of the excited PT to the ground state experiences a barrier of 13.8 kcal/mol for the acidic form to reach the (S(1)/S(0)) conical intersection. Such internal conversion can be enhanced with the increase in excitation energy, which will reduce the fluorescence intensity as observed experimentally.     

Energetics of cytosine singlet excited-state decay paths--a difficult case for CASSCF and CASPT2

Three deactivation paths for singlet excited cytosine are calculated at the CASPT2//CASSCF (complete active space second-order perturbation//complete active space self-consistent field) level of theory, using extended active spaces that allow for a reliable characterization of the paths and their energies. The lowest energy path, with a barrier of approximately 0.1 eV, corresponds to torsion of the C5-C6 bond, and the decay takes place at a conical intersection analogous to the one found for ethylene and its derivatives. There is a further path with a low energy barrier of approximately 0.2 eV associated with the (n(N),pi*) state which could also be populated with a low energy excitation. The path associated with a conical intersection between the ground and (n(O),pi*) states is significantly higher in energy (> 1 eV). The presence of minima on the potential energy surface for the (n,pi*) states that could contribute to the biexponential decay found in the gas phase was investigated, but could not be established unequivocally.     

Excited-state relaxation of protonated adenine.

The excited state dynamics of protonated adenine in the gas phase were investigated by femtosecond pump-probe transient mass spectroscopy. Adenine was protonated in an electrospray ionization source and transferred to a Paul trap. Two femtosecond laser pulses at 266 nm and 800 nm excited the lowest electronic pipi* state and probed the excited-state dynamics by monitoring ion fragment formation. The measured excited state decay is monoexponential with a lifetime shorter than 161 fs. This agrees with a theoretical prediction of very fast internal conversion via a conical intersection with the ground state.     

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.     

Unified model for the ultrafast decay of pyrimidine nucleobases

Ultrafast decay processes detected after absorption of UV radiation in gas-phase pyrimidine nucleobases uracil, thymine, and cytosine are ascribed to the barrierless character of the pathway along the low-lying 1(pipi*) hypersurface connecting the Franck-Condon region with an out-of-plane distorted ethene-like conical intersection with the ground state. Longer lifetime decays and low quantum yield emission are on the other hand related to the presence of a 1(pipi*) state planar minimum on the S1 surface and the barriers to access other conical intersections. A unified model for the three systems is established on the basis of accurate multiconfigurational CASPT2 calculations, whereas the effect of the different levels of theory on the results is carefully analyzed.     

Dynamics of photoinduced processes in adenine and thymine base pairs

The excited-state dynamics of adenine and thymine dimers and the adenine-thymine base pair were investigated by femtosecond pump-probe ionization spectroscopy with excitation wavelengths of 250-272 nm. The base pairs showed a characteristic ultrafast decay of the initially excited pi pi* state to an n pi* state (lifetime tau(pi pi*) approximately 100 fs) followed by a slower decay of the latter with tau(n pi*) approximately 0.9 ps for (adenine)2, tau(n pi*) = 6-9 ps for (thymine)2, and tau(n pi*) approximately 2.4 ps for the adenine-thymine base pair. In the adenine dimer, a competing decay of the pi pi* state via the pi sigma* state greatly suppressed the n pi* state signals. Similarities of the excited-state decay parameters in the isolated bases and the base pairs suggest an intramonomer relaxation mechanism in the base pairs.     

Ultrafast decay of electronically excited singlet cytosine via a pi,pi* to n(O),pi* state switch

Singlet fluorescence lifetimes of adenosine, cytidine, guanosine, and thymidine, determined by femtosecond pump-probe spectroscopy (Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2000, 122, 9348. Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2001, 123, 10370), show that the excited states produced by 263 nm light in these nucleosides decay in the subpicosecond range (290-720 fs). Ultrafast radiationless decay to the ground state greatly reduces the probability of photochemical damage. In this work we present a theoretical study of isolated cytosine, the chromophore of cytidine. The experimental lifetime of 720 fs indicates that there must be an ultrafast decay channel for this species. We have documented the possible decay channels and approximate energetics, using a valence-bond derived analysis to rationalize the structural details of the paths. The mechanism favored by our calculations and the experimental data involves (1) a two-mode decay coordinate composed of initial bond length inversion followed by internal vibrational energy redistribution (IVR) to populate a carbon pyramidalization mode, (2) a state switch between the pi,pi* and nO,pi* (excitation from oxygen lone pair) excited states, and (3) decay to the ground state through a conical intersention. A second decay path through the nN,pi* state (excitation from the nitrogen lone pair), with a higher barrier, involves out-of-plane bending of the amino substituent.     

Cyclooctatetraene computational photo- and thermal chemistry: a reactivity model for conjugated hydrocarbons

We use ab initio CASSCF and CASPT2 computations to construct the composite multistate relaxation path relevant to cycloocta-1,3,5,7-tetraene singlet photochemistry. The results show that an efficient population of the dark excited state (S(1)) takes place after ultrafast decay from the spectroscopic excited state (S(2)). A planar D(8)(h)-symmetric minimum represents the collecting point on S(1). Nonadiabatic transitions to S(0) appear to be controlled by two different tetraradical-type conical intersections, which are directly accessible from the S(1) minimum following specific excited-state reaction paths. The higher-energy conical intersection belongs to the same type of intersections previously documented in linear and cyclic conjugated hydrocarbons and features a triangular -(CH)(3)- kink. This point mediates both cis --> trans photoisomerization and cyclopropanation reactions. The lowest energy conical intersection has a boat-shaped structure. This intersection accounts for production of semibullvalene or for double-bond shifting. The mapping of both photochemical and thermal reaction paths (including also Cope rearrangements, valence isomerizations, ring inversions, and double-bond shifting) has allowed us to draw a comprehensive reactivity scheme for cyclooctatetraene, which rationalizes the experimental observations and documents the complex network of photochemical and thermal reaction path interconnections. The factors controlling the selection and accessibility of a number of conjugated hydrocarbon prototype conical intersections and ground-state relaxation channels are discussed.     

CASCSF study on the photochemical transposition reactions of pyrazines

Several reaction pathways for the photochemical transformations of methyl-substituted pyrazine in its first excited state 1(pi --> pi*) have been determined using the CASSCF (six-orbital/six-electron active space) and MP2-CAS methods with the 6-311G(d) basis set. Our model investigations suggest that conical intersections play a crucial role in the photoisomerization of pyrazines. Moreover, the present theoretical findings indicate that all of the photoisomerizations of pyrazines adopt the same reaction path as follows: pyrazine --> Franck-Condon region --> conical intersection --> pyrimidine. That is, although an excited-state pyrazine molecule can initiate a phototransposition process easily, this process can be completed on the ground-state potential energy surface after passage through a conical intersection where a fast, radiationless decay is possible. The existence of these nonadiabatic reaction pathways is consistent with the available experimental observations of the photochemistry and photophysics of pyrazine and its methyl derivatives. In the present work, we propose a simple p-pi orbital topology model, which can be used as a diagnostic tool to predict the location of the conical intersections, as well as the geometries of the phototransposition products of various heterocycles.     

A comparative study on the photochemistry of two bipyridyl derivatives: [2,2'-bipyridyl]-3,3'-diamine and [2,2'-bipyridyl]-3,3'-diol.

The two isoelectronic bipyridyl derivatives, [2,2'-bipyridyl]-3,3'-diamine and [2,2'-bipyridyl]-3,3'-diol, are experimentally known to undergo very different excited-state double-proton-transfer processes, which result in fluorescence quantum yields that differ by four orders of magnitude. Herein, density functional theory (DFT), time-dependent DFT (TDDFT), and complete active space self-consistent field (CASSCF) calculations are used to study the double-proton-transfer processes in the ground and first singlet pi-->pi* excited state. The quantum-chemistry calculations indicate 1) the existence of only one energy minimum in the ground electronic state corresponding to reactants (thus avoiding the possibility of a fast fluorescent relaxation process from the photoproducts region), 2) an endoergic process of the complete double proton transfer, and 3) the presence of a conical intersection in the excited intermediate region of [2,2'-bipyridyl]-3,3'-diamine. These facts explain the very low fluorescence quantum yield in [2,2'-bipyridyl]-3,3'-diamine compared to [2,2'-bipyridyl]-3,3'-diol.     

Singlet excited-state lifetimes of cytosine derivatives measured by femtosecond transient absorption

Lifetimes of the lowest excited singlet (S1) electronic states of various derivatives of the pyrimidine nucleobase cytosine (Cyt) were measured by the femtosecond transient absorption technique. The bases were excited in room-temperature aqueous solution at 265 nm using approximately 200 fs pump pulses from a titanium-sapphire laser system. The decay of excited-state absorption (ESA) at visible probe wavelengths was used to determine the S1 lifetimes of a variety of modified Cyt compounds at different pH values by global fitting. Identical lifetimes were observed for Cyt and cytidine (Cyd) within experimental uncertainty, but ESA by the ribonucleoside was considerably stronger, suggesting that the ribose group increases the oscillator strength of the S1 --> SN transition. The S1 lifetime of the important minor base 5-methylcytosine (m5Cyt) is 7.2 +/- 0.4 ps at pH 6.8. The same lifetime was measured for the ribonucleoside 5-methylcytidine, but sugar substitution again increased the strength of the ESA signal. Protonation of Cyd and m5Cyt at low pH led to a modest decrease in their S1 lifetimes. On the other hand, deprotonation of Cyt and m5Cyt significantly increased the lifetime of their respective S1 states. These trends support the intermediacy of the n,pi* state localized on the carbonyl oxygen in the nonradiative decay mechanism of Cyt. Longer S1 lifetimes were observed for 5-fluorocytosine and N4-acetylcytosine. Collectively, these results illustrate the great potential of femtosecond laser spectroscopy for investigating excited-state dynamics in DNA and DNA components.     

Nonradiative deactivation mechanisms of uracil, thymine, and 5-fluorouracil: a comparative ab initio study

The mechanisms of the ultrafast nonradiative deactivation of uracil and its substituted derivatives thymine (5-methyluracil) and 5-fluorouracil after absorption of UV light are explored and compared by means of ab initio multistate (MS) CASPT2 calculations. The MS-CASPT2 method is applied for the calculation of potential energy profiles, especially for the geometry optimization in the electronically excited state, with the aim of an accurate prediction of deactivation pathways. The resulting energy curves of each molecule exhibit that the conical intersection between the (1)ππ* and ground states is accessible via small energy barriers from the minimum in the (1)ππ* state as well as from that in the (1)nπ* state. The barrier of 5-fluorouracil in the (1)ππ* state is calculated to be definitely higher than those of uracil and thymine, which is consistent with experiments and suggests that the elongation of the excited-state lifetime of uracil by fluorine substitution is significantly contributed from intrinsic electronic effect of the molecule. However, no evidence of the experimentally observed longer excited-state lifetime of thymine than uracil is found in the presently calculated MS-CASPT2 potential energy curves in the (1)ππ* and (1)nπ* states, implying nonnegligible contribution of other factors such as solvation effect and substituent mass to the photoinduced dynamics of uracil derivatives.     

The different photoisomerization efficiency of azobenzene in the lowest n pi* and pi pi* singlets: the role of a phantom state

Azobenzene E<==>Z photoisomerization, following excitation to the bright S(pi pi*) state, is investigated by means of ab initio CASSCF optimizations and perturbative CASPT2 corrections. Specifically, by elucidating the S(pi pi*) deactivation paths, we explain the mechanism responsible for azobenzene photoisomerization, the lower isomerization quantum yields observed for the S(pi pi*) excitation than for the S1(n pi*) excitation in the isolated molecule, and the recovery of the Kasha rule observed in sterically hindered azobenzenes. We find that a doubly excited state is a photoreaction intermediate that plays a very important role in the decay of the bright S(pi pi*). We show that this doubly excited state, which is immediately populated by molecules excited to S(pi pi*), drives the photoisomerization along the torsion path and also induces a fast internal conversion to the S1(n pi*) at a variety of geometries, thus shaping (all the most important features of) the S(pi pi*) decay pathway and photoreactivity. We reach this conclusion by determining the critical structures, the minimum energy paths originating on the bright S(pi pi*) state and on other relevant excited states including S1(n pi*), and by characterizing the conical intersection seams that are important in deciding the photochemical outcome. The model is consistent with the most recent time-resolved spectroscopic and photochemical data.     

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.