Comparison of the non-radiative decay mechanisms of 4-pyrimidinone and uracil: an ab initio study

Experimental and theoretical study of the cytosine tautomerism through excited states

CONTEXT

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

METHODS

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

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Non-radiative deactivation of excited cytosine: probing of different DFT functionals and basis sets in solvents with different polarity

The vertical excitation and emission energies of the cytosine oxo-amino form were calculated at the TD DFT level of theory with several functionals, basis sets, and solvents with different polarity (PCM). They were compared with the experimental UV absorption and fluorescence maxima, revealing that the minimal deviation of the vertical excitation energies from the of UV absorption maxima can be achieved when the hybrid functional B3LYP is applied within the calculations. Regular correlations like curves with saturation between the vertical excitation/emission energies and the dielectric constants were registered. The relaxation of the ¹ππ^* excited state through an ethylene-like conical intersection S₀/S₁ should occur through decrease of the decay rate (commented qualitatively here) with the rise of the solvent polarity.

Disclosing the Role of C4-Oxo Substitution in the Photochemistry of DNA and RNA Pyrimidine Monomers: Formation of Photoproducts from the Vibrationally Excited Ground State

Oxo and amino substituted purines and pyrimidines have been suggested as protonucleobases participating in ancient pre-RNA forms. Considering electromagnetic radiation as a key environmental selection pressure on early Earth, the investigation of the photophysics of modified nucleobases is crucial to determine their viability as nucleobases' ancestors and to understand the factors that rule the photostability of natural nucleobases. In this Letter, we combine femtosecond transient absorption spectroscopy and quantum mechanical simulations to reveal the photochemistry of 4-pyrimidinone, a close relative of uracil. Irradiation of 4-pyrimidinone with ultraviolet radiation populates the S₁(ππ*) state, which decays to the vibrationally excited ground state in a few hundred femtoseconds. Analysis of the postirradiated sample in water reveals the formation of a 6-hydroxy-5H-photohydrate and 3-(N-(iminomethyl)imino)propanoic acid as the primary photoproducts. 3-(N-(Iminomethyl)imino)propanoic acid originates from the hydrolysis of an unstable ketene species generated from the C4-N3 photofragmentation of the pyrimidine core.

A molecular insight into the photophysics of barbituric acid, a candidate for canonical nucleobases' ancestor

This work investigates the photophysics of barbituric acid at different pH conditions using ab initio methods. Our calculations ascribe the most intense bands at ca. 260 nm at neutral pH and 210 nm at acidic pH conditions in the absorption spectra of this chromophore to the lowest lying ππ* transitions. Consistently with the ultrashort excited state lifetimes experimentally registered, the potential energy landscapes of both the neutral and deprotonated forms of barbituric acid combined with the interpretation of their transient absorption spectra suggest the deactivation of these systems along the singlet manifold. Compared to uracil, its closest natural nucleobase, barbituric acid presents a red shifted absorption spectrum, due to the lowering by more than 0.5 eV of the lowest-energy ππ* excited state, and a much more complex topography of the S₁ potential energy surface, with several energetically accessible local minima. This fact, however, does not affect the excited state lifetimes, which for barbituric acid were experimentally registered in the sub-ps time scale.

Photoinduced phenomena in water solution of melamine explaining the photostability of the compound

Two tautomers of melamine (triamino and imino-diamino) were studied at the BLYP/aug-cc-pVDZ theoretical level. It was found that the two tautomers are bridged with the ¹πσ^* excited-state reaction path. The high photostability of melamine in water solution was explained with the mechanism of ring deformation which occurs along the ¹ππ^* excited-state reaction path. The two mechanisms are investigated at the TD BLYP level of theory using the linear interpolation in internal coordinates approach.

Theoretical insights on the excited-state-deactivation mechanisms of protonated thymine and cytosine

Ab initio and surface-hopping nonadiabatic dynamics simulation methods were employed to investigate relaxation mechanisms in protonated thymine (TH⁺) and cytosine (CH⁺). A few conical intersections were located between ¹ππ* and S₀ states for each system with the CASSCF (8,8) theoretical model and relevant contributions to the deactivation mechanism of titled systems were addressed by the determination of potential energy profiles at the CASPT2 (12,10) theoretical level. It was revealed that the relaxation of the ¹ππ* state of the most stable conformer of both systems to the ground state is mostly governed by the accessible S₁/S₀ conical intersection resulting from the barrier-free out-of-plane deformation. Interestingly, it was exhibited that the ring puckering coordinate driven from the C₆ position of the heterocycle ring in TH⁺ and CH⁺ plays the most prominent role in the deactivation mechanism of considered systems. Our ab initio results are also supported by excited-state nonadiabatic dynamics simulations based on ADC(2), describing the ultrashort S₁ lifetime of TH⁺/CH⁺ by analyzing trajectories leading excited systems to the ground. It was confirmed that the excited-state population mostly relaxes to the ground via the ring puckering coordinate from the C₆ moiety. Overall, the theoretical results of this study shed light on the deactivation mechanism of protonated DNA bases.

New Insights about the Photostability of DNA/RNA Bases: Triplet nπ* State Leads to Effective Intersystem Crossing in Pyrimidinones

The high photostability of DNA/RNA nucleobases is attributed to the effective internal conversions of their bright ¹ππ* states to the ground state through conical intersections. Intersystem crossing (ISC) from singlet to triplet excited states is a minor decay pathway in nucleobases and it is observed with ∼1-2% quantum yields (QYs) in pyrimidine bases. Presumably, ISC in pyrimidines takes place from the dark singlet ¹nπ* state to the lowest triplet ³ππ* state. However, recent studies showed that ISC from the initial populated bright ¹ππ* state to higher energy triplet ³nπ* states indeed occurs in the subpicosecond timescale. Such a mechanism is still poorly understood since direct observation of this pathway is challenging. Herein, excited state dynamics of three pyrimidinones, which share the same skeleton with pyrimidine bases, is investigated in different solvents. Compared to canonical pyrimidine bases, removing the oxygen atom at the C4 position revokes the low-lying dark ¹nπ* state in pyrimidinones, resulting in direct ISC from the S₁ (¹ππ*) state to triplet T₃ (³nπ*) state with much higher QYs. Meanwhile, hydrogen bonding between the carbonyl group in pyrimidinones and protic solvents can accelerate vibrational cooling of the hot S₁ (¹ππ*) state, leading to higher fluorescence QYs and smaller ISC rate constants. These results not only evidence the hypothesis of the direct ¹ππ* → ³nπ* ISC mechanism, but also contribute to a better understanding of triplet formation in pyrimidines.

Influence of the N atom position on the excited state photodynamics of protonated azaindole

We present a study of the photofragmentation of three protonated azaindole molecules - 7-azaindole, 6-azaindole, and 5-azaindole - consisting of fused pyrrole-pyridine bicyclic aromatic systems, in which the pyridinic (protonated) nitrogen heteroatom is located at the 7, 6, and 5 positions, respectively. Photofragmentation electronic spectra of the isolated aforementioned azaindolinium cations reveal that their photodynamics extends over timescales covering nine orders of magnitude and provide evidence about the resultant fragmentation pathways. Moreover, we show how the position of the heteroatom in the aromatic skeleton influences the excited state energetics, fragmentation pathways, and fragmentation timescales. Computed ab initio adiabatic transition energies are used to assist the assignation of the spectra, while geometry optimisation in the excited electronic states as well as ab initio calculations along the potential surfaces demonstrate the role of ππ*/πσ* coupling and/or large geometry changes in the dynamics of these species. Evidence supporting the formation of Dewar valence isomers as intermediates involved in sub-picosecond relaxation processes is discussed.

Excited State Deactivation Mechanism in Protonated Uracil: New Insights from Theoretical Studies

We have conducted here a theoretical exploration, discussing the distinct excited state lifetimes reported experimentally for the two lowest lying protonated isomers of uracil. In this regard, the first-principal computational levels as well as the nonadiabatic surface hopping dynamics have been employed. It has been revealed that relaxation of the ¹ππ* state of enol-enol form (EE⁺) to the ground is barrier-free via out-of-plane coordinates, resulting in an ultrashort S₁ lifetime of this species. For the second most stable isomer (EK⁺), however, a significant barrier predicted in the CASPT2 S₁ potential energy profile along the twisting coordinate has been proposed to explain the relevant long lifetime reported experimentally.

Photoprotection: extending lessons learned from studying natural sunscreens to the design of artificial sunscreen constituents

Evolution has ensured that plants and animals have developed effective protection mechanisms against the potentially harmful effects of incident ultraviolet radiation (UVR). Tanning is one such mechanism in humans, but tanning only occurs post-exposure to UVR. Hence, there is ever growing use of commercial sunscreens to pre-empt overexposure to UVR. Key requirements for any chemical filter molecule used in such a photoprotective capacity include a large absorption cross-section in the UV-A and UV-B spectral regions and the availability of one or more mechanisms whereby the absorbed photon energy can be dissipated without loss of the molecular integrity of the chemical filter. Here we summarise recent experimental (mostly ultrafast pump-probe spectroscopy studies) and computational progress towards unravelling various excited state decay mechanisms that afford the necessary photostability in chemical filters found in nature and those used in commercial sunscreens. We also outline ways in which a better understanding of the photophysics and photochemistry of sunscreen molecules selected by nature could aid the design of new and improved commercial sunscreen formulations.

Role of Solvent Polarity and Hydrogen-Bonding on Excited-State Fluorescence of 3-[(E)-{4-[Dimethylamino]benzylidene}amino]-2-naphthoic Acid (DMAMN): Isomerization vs Rotomerization

The present experimental and theoretical study on a new chromophore DMAMN of the type push-π-pull (push = dimethylaniline, π = imine, pull = 2-naphthoic acid), allows understanding of the mechanism by which the molecular conformational undergoes isomerization/rotomerization following electronic excitation. The steady-state fluorescence spectra of this compound, carried out in solvents of different polarities and proticities, showed significant changes in both the shape and peak positions. The wavelength and intensity change depend on the polarity and hydrogen-bonding environment. In highly polar solvents, the emission is weak and red-shifted compared to that for cyclohexane, but it is more red-shifted in moderate aprotic polar solvents. In hydroxyl solvents, a new weak low-energy emission band appears at ∼525 nm, attributed to the intermolecularly H-bonded open conformer. On the basis of the generated potential energy landscapes of the ground state and low-lying excited state in the gas phase and solution, we found that selective photon absorption, brings this molecule to a "bright" state, from which N═C isomerization Z → E, takes place. This isomerization in gas-phase and low-polarity solvents leads to two minima with a barrier, whereas in highly polar-protic media, it gives one minimum on the S₁ surface with low ΔE_S1/T1 (0.17 eV), facilitating deactivation via ISC.

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

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

Hydrogen detachment driven by a repulsive 1πσ* state – an electron localization function study of 3-amino-1,2,4-triazole

Electron localization function analysis reveals the details of a charge induced hydrogen detachment mechanism of 3-amino-1,2,4-triazole, identified recently as responsible for phototautomerization of the molecule.

Non-radiative processes in protonated diazines, pyrimidine bases and an aromatic azine

The excited state lifetimes of DNA bases are often very short due to very efficient non-radiative processes assigned to the ππ*-nπ* coupling. A set of protonated aromatic diazine molecules (pyridazine, pyrimidine and pyrazine C4H5N2(+)) and protonated pyrimidine DNA bases (cytosine, uracil and thymine), as well as the protonated pyridine (C5H6N(+)), have been investigated. For all these molecules except one tautomer of protonated uracil (enol-keto), electronic spectroscopy exhibits vibrational line broadening. Excited state geometry optimization at the CC2 level has been conducted to find out whether the excited state lifetimes measured from line broadening can be correlated to the calculated ordering of the ππ* and nπ* states and the ππ*-nπ* energy gap. The short lifetimes, observed when one nitrogen atom of the ring is not protonated, can be rationalized by relaxation of the ππ* state to the nπ* state or directly to the electronic ground state through ring puckering.

Electronic and structural elements that regulate the excited-state dynamics in purine nucleobase derivatives

The excited-state dynamics of the purine free base and 9-methylpurine are investigated using experimental and theoretical methods. Femtosecond broadband transient absorption experiments reveal that excitation of these purine derivatives in aqueous solution at 266 nm results primarily in ultrafast conversion of the S2(ππ*) state to the vibrationally excited (1)nπ* state. Following vibrational and conformational relaxation, the (1)nπ* state acts as a doorway state in the efficient population of the triplet manifold with an intersystem crossing lifetime of hundreds of picoseconds. Experiments show an almost 2-fold increase in the intersystem crossing rate on going from polar aprotic to nonpolar solvents, suggesting that a solvent-dependent energy barrier must be surmounted to access the singlet-to-triplet crossing region. Ab initio static and surface-hopping dynamics simulations lend strong support to the proposed relaxation mechanism. Collectively, the experimental and computational results demonstrate that the accessibility of the nπ* states and the topology of the potential energy surfaces in the vicinity of conical intersections are key elements in controlling the excited-state dynamics of the purine derivatives. From a structural perspective, it is shown that the purine chromophore is not responsible for the ultrafast internal conversion in the adenine and guanine monomers. Instead, C6 functionalization plays an important role in regulating the rates of radiative and nonradiative relaxation. C6 functionalization inhibits access to the (1)nπ* state while simultaneously facilitating access to the (1)ππ*(La)/S0 conical intersection, such that population of the (1)nπ* state cannot compete with the relaxation pathways to the ground state involving ring puckering at the C2 position.

Excited states of proton-bound DNA/RNA base homodimers: pyrimidines

We are presenting the electronic photofragment spectra of the protonated pyrimidine DNA base homodimers. Only the thymine dimer exhibits a well structured vibrational progression, while the protonated monomer shows broad vibrational bands. This shows that proton bonding can block some nonradiative processes present in the monomer.

Effect of Ag(+) on the Excited-State Properties of a Gas-Phase (Cytosine)2Ag(+) Complex: Electronic Transition and Estimated Lifetime

Recently, DNA molecules have received great attention because of their potential applications in material science. One interesting example is the production of highly fluorescent and tunable DNA-Agn clusters with cytosine (C)-rich DNA strands. Here, we report the UV photofragmentation spectra of gas-phase cytosine···Ag(+)···cytosine (C2Ag(+)) and cytosine···H(+)···cytosine (C2H(+)) complexes together with theoretical calculations. In both cases, the excitation energy does not differ significantly from that of isolated cytosine or protonated cytosine, indicating that the excitation takes place on the DNA base. However, the excited-state lifetime of the C2H(+) (τ = 85 fs), estimated from the bandwidth of the spectrum, is at least 2 orders of magnitude shorter than that of the C2Ag(+) (τ > 5000 fs). The increased excited-state lifetime upon silver complexation is quite unexpected, and it clearly opens the question about what factors are controlling the nonradiative decay in pyrimidine DNA bases. This is an important result for the expanding field of metal-mediated base pairing and may also be important to the photophysical properties of DNA-templated fluorescent silver clusters.

Solvent influence on the excited states of the oxo form of barbituric acid and the mechanisms of the out-of-plane non-radiative elongation of the N-H bond: a comparative theoretical and experimental study

The solvent influence on the excited states, emission and absorption energies of the oxo for of barbituric acid was studied with experimental (UV and fluorescence spectra) and theoretical methods. The excited-state reaction paths of the out-pf-plane elongation of the N-H bond of the oxo form of barbituric acid were also investigated (TD DFT level) to the conical intersections mediating internal conversions to the ground state. The (1)nσ* excited state was found to be the driven electronic state. We found that the increase of the polarity of the solvent reduces the (1)nσ* excited state decay rate through a non-planar conical intersection and increases its energy. Thus, solvents with higher polarity disfavor the non-radiative decay through conical intersections.

Probing the Conical Intersection Dynamics of the RNA Base Uracil by UV-Pump Stimulated-Raman-Probe Signals; Ab Initio Simulations

Nonadiabatic electron and nuclear dynamics of photoexcited molecules involving conical intersections is of fundamental importance in many reactions such as the self-protection mechanism of DNA and RNA bases against UV irradiation. Nonlinear vibrational spectroscopy can provide an ultrafast sensitive probe for these processes. We employ a simulation protocol that combines nonadiabatic on-the-fly molecular dynamics with a mode-tracking algorithm for the simulation of femtosecond stimulated Raman spectroscopy (SRS) signals of the high frequency C–H- and N–H-stretch vibrations of the photoexcited RNA base uracil. The simulations rely on a microscopically derived expression that takes into account the path integral of the excited state evolution and the pulse shapes. Analysis of the joint time/frequency resolution of the technique reveals a matter chirp contribution that limits the inherent temporal resolution. Characteristic signatures of relaxation dynamics mediated in the vicinity of conical intersection are predicted. The C–H and N–H spectator modes provide high sensitivity to their local environment and act as local probes with submolecular and high temporal resolution.

Monitoring Non-Adiabatic Dynamics of the RNA Base Uracil by UV-Pump-IR-Probe Spectroscopy

Resolving the excited state dynamics of DNA- and RNA- nucleobases has attracted considerably attention. UV irradiation of the isolated nucleobases leads to the population of an electronic excited state which is quenched by internal conversion mediated by conical intersections on an ultrafast timescale. We present non-adiabatic on-the-fly molecular dynamics simulations of the UV-pump-IR-probe signal of the pyrimidine nucleobase uracil using a novel semiclassical protocol which takes into account the path integral over the excited state vibrational dynamics and properly describes the joint temporal and spectral resolution of the technique. Simulations of vibrational motions of carbonyl fingerprint modes in the electronically excited states reveal clear signatures of different relaxation pathways on a timescale of hundreds of femtoseconds which arise from an ultrafast branching in the excited state. We show that the inherent temporal and spectral resolution of the technique is not purely instrumental but also depends on the vibrational fluctuation timescale.

ONIOM approach for non-adiabatic on-the-fly molecular dynamics demonstrated for the backbone controlled Dewar valence isomerization

Non-adiabatic on-the-fly molecular dynamics (NA-O-MD) simulations require the electronic wavefunction, energy gradients, and derivative coupling vectors in every timestep. Thus, they are commonly restricted to the excited state dynamics of molecules with up to ≈20 atoms. We discuss an approximation that combines the ONIOM(QM:QM) method with NA-O-MD simulations to allow calculations for larger molecules. As a proof of principle we present the excited state dynamics of a (6-4)-lesion containing dinucleotide (63 atoms), and especially the importance to include the confinement effects of the DNA backbone. The method is able to include electron correlation on a high level of theory and offers an attractive alternative to QM:MM approaches for moderate sized systems with unknown force fields.