Ultrafast Excited-State Dynamics of Cytosine Aza-Derivative and Analogues

Ribose and Deoxyribose Group Alter Excited-State Dynamics of 5-Azacytosine in Solution

5-Azacytosine (5-AC) is one of the best interesting noncanonical nucleobases due to its functionalization and structural imitation of natural bases. 5-AC can be used as the scaffold of two important chemotherapeutic medicines, 5-azacytidine and 2'-deoxy-5-azacytidine. Furthermore, increased sensitivity to UV leads to the photochemical effects of 5-AC also attracted attention. Yet, no study has been reported to explore the effect of glycosyl groups on the photophysical and photochemical properties of 5-AC, which can help to reveal the photostability of related actual clinic drugs. In this study, the excited-state dynamics of 5-azacytidine and 2'-deoxy-5-azacytidine are studied by femtosecond transient absorption and quantum-chemical calculations while revisiting that of 5-AC with a wider probe spectral range. It is shown that glycosyl substitution on the N1 position leads to ultrafast excited-state relaxation within several picoseconds in both nucleosides, which is distinct compared with the 17 ps lifetime seen in 5-AC. It is proposed that these changes are due to altering the energy level of the dark nπ* state. Moreover, our results suggest that it should be cautioned to simply replace sugar groups with methyl groups when doing a theoretical calculation study on nucleobases and their derivatives.

Quantum mechanics/molecular mechanics studies on the excited-state decay mechanisms of cytidine aza-analogues: 5-azacytidine and 2'-deoxy-5-azacytidine in aqueous solution

The excited state properties and deactivation pathways of two DNA methylation inhibitors, i.e., 5-azacytidine (5ACyd) and 2'-deoxy-5-azacytidine (5AdCyd) in aqueous solution, are comprehensively explored with the QM(CASPT2//CASSCF)/MM protocol. We systematically map the feasible decay mechanisms based on the obtained excited-state decay paths involving all the identified minimum-energy structures, conical intersections, and crossing points driving the different internal conversion (IC) and intersystem crossing (ISC) routes in and between the 1ππ*, 1nπ*, 3ππ*, 3nπ*, and S0 states. Unlike the 1nπ* state below the 1ππ* state in 5ACyd, deoxyribose group substitution at the N1 position leads to the 1ππ* state becoming the S1 state in 5AdCyd. In 5ACyd and 5AdCyd, the initially populated 1ππ* state mainly deactivates to the S0 state through the direct 1ππ* → S0 IC or mediated by the 1nπ* state. The former nearly barrierless IC channel of 1ππ* → S0 occurs ultrafast via the nearby low-lying 1ππ*/S0 conical intersection. In the latter IC channel of 1ππ* → 1nπ* → S0, the initially photoexcited 1ππ* state first approaches the nearby S2/S1 conical section 1ππ*/1nπ* and then undergoes efficient IC to the 1nπ* state, followed by the further IC to the initial S0 state via the S1/S0 conical intersection 1nπ*/S0. The 1nπ*/S0 conical intersection is estimated to be located 6.0 and 4.9 kcal mol-1 above the 1nπ* state minimum in 5ACyd and 5AdCyd, respectively, at the QM(CASPT2)/MM level. In addition to the efficient singlet-mediated IC channels, the minor ISC routes would populate 1ππ* to T1(ππ*) through 1ππ* → T1 or 1ππ* → 1nπ* → T1. Relatively, the 1ππ* → 1nπ* → T1 route benefits from the spin-orbit coupling (SOC) of 1nπ*/3ππ* of 8.7 cm-1 in 5ACyd and 10.2 cm-1 in 5AdCyd, respectively. Subsequently, the T1 system will approach the nearby T1/S0 crossing point 3ππ*/S0 driving it back to the S0 state. Given the 3ππ*/S0 crossing point located above the T1 minimum and the small T1/S0 SOC, i.e., 8.4 kcal mol-1 and 2.1 cm-1 in 5ACyd and 6.8 kcal mol-1 and 1.9 cm-1 in 5AdCyd, respectively, the slow T1 → S0 would trap the system in the T1 state for a while. The present work could contribute to understanding the mechanistic photophysics and photochemistry of similar aza-nucleosides and their derivatives.

Vertical Excitation Energies and Lifetimes of the Two Lowest Singlet Excited States of Cytosine, 5-Aza-cytosine, and the Triazine Family: Quantum Mechanics-Molecular Mechanics Studies

A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ₁=0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S₁ → S₀ de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s-triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH₂ ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.

Quantum mechanics/molecular mechanics studies on mechanistic photophysics of cytosine aza-analogues: 2,4-diamino-1,3,5-triazine and 2-amino-1,3,5-triazine in aqueous solution

The excited-state properties and photophysics of cytosine aza-analogues, i.e., 2,4-diamino-1,3,5-triazine (2,4-DT) and 2-amino-1,3,5-triazine (2-AT) in solution have been systematically explored using the QM(MS-CASPT2//CASSCF)/MM approach. The excited-state nonradiative relaxation mechanisms for the initially photoexcited S₁(ππ*) state decay back to the S₀ state are proposed in terms of the present computed minima, surface crossings (conical intersections and singlet-triplet crossings), and excited-state decay paths in the S₁, S₂, T₁, T₂, and S₀ states. Upon photoexcitation to the bright S₁(ππ*) state, 2,4-DT quickly relaxes to its S₁ minimum and then overcomes a small energy barrier of 5.1 kcal mol^-1 to approach a S₁/S₀ conical intersection, where the S₁ system hops to the S₀ state through S₁ → S₀ internal conversion (IC). In addition, at the S₁ minimum, the system could partially undergo intersystem crossing (ISC) to the T₁ state, followed by further ISC to the S₀ state via the T₁/S₀ crossing point. In the T₁ state, an energy barrier of 7.9 kcal mol^-1 will trap 2,4-DT for a while. In parallel, for 2-AT, the system first relaxes to the S₁ minimum and then S₁ → S₀ IC or S₁ → T₁ → S₀ ISCs take place to the S₀ state by surmounting a large barrier of 15.3 kcal mol^-1 or 11.9 kcal mol^-1, respectively, which heavily suppress electronic transition to the S₀ state. Different from 2,4-DT, upon photoexcitation in the Franck-Condon region, 2-AT can quickly evolve in an essentially barrierless manner to nearby S₂/S₁ conical intersection, where the S₂ and T₁ states can be populated. Once it hops to the S₂ state, the system will overcome a relatively small barrier (6.6 kcal mol^-1vs. 15.3 kcal mol^-1) through IC to the S₀ state. Similarly, an energy barrier of 11.9 kcal mol^-1 heavily suppresses the T₁ state transformation to the S₀ state. The present work manifests that the amination/deamination of the triazine rings can affect some degree of different vertical and adiabatic excitation energies and nonradiative decay pathways in solution. It not only rationalizes excited-state decay dynamics of 2,4-DT and 2-AT in aqueous solution but could also provide insights into the understanding of the photophysics of aza-nucleobases.

Excited State Dynamics of Methylated Guanosine Derivatives Revealed by Femtosecond Time-resolved Spectroscopy

Methylated DNA/RNA nucleobases are important epigenetic marks in living species and play an important role for targeted therapies. Moreover, they could bring significant changes to the photo-stability of nucleic acid, leading these sites become mutational hotspots for disease such as skin cancer. While a number of studies have demonstrated the relationship between excited state dynamics and the biological function of methylated cytosine in DNA, investigations aimed at unraveling the excited state dynamics of methylated guanosine in RNA have been largely overlooked. In this work, influence of methylation on the excited state dynamics of guanosine is studied by using femtosecond time-resolved spectroscopy. Our results suggest that the effect of methyl substitution on the photophysical properties of guanosine is position sensitive. N1-methylguanosine shows very similar excited state dynamics as that in guanosine, while almost one order of magnitude longer lifetime of the La state is observed in N2, N2-dimethylguanosine. Notably, N7-methylation can lead to a new minimum on the La state, which shows a two orders of magnitude longer excited state lifetime compared with guanosine. These findings not only help understanding excited state dynamics of methylated guanosines, but also lay the foundation for further studying DNA/RNA strands incorporated with these bases.

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning's cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.

Solvent-Dependent Stabilization of a Charge Transfer State is the Key to Ultrafast Triplet State Formation in an Epigenetic DNA Nucleoside

2'-Deoxy-5-formylcytidine (5fdCyd), a naturally occurring nucleoside found in mammalian DNA and mitochondrial RNA, exhibits important epigenetic functionality in biological processes. Because it efficiently generates triplet excited states, it is an endogenous photosensitizer capable of damaging DNA, but the intersystem crossing (ISC) mechanism responsible for ultrafast triplet state generation is poorly understood. In this study, time-resolved mid-IR spectroscopy and quantum mechanical calculations reveal the distinct ultrafast ISC mechanisms of 5fdCyd in water versus acetonitrile. Our experiment indicates that in water, ISC to triplet states occurs within 1 ps after 285 nm excitation. PCM-TD-DFT computations suggest that this ultrafast ISC is mediated by a singlet state with significant cytosine-to-formyl charge-transfer (CT) character. In contrast, ISC in acetonitrile proceeds via a dark ¹ nπ* state with a lifetime of ∼3 ps. CT-induced ISC is not favored in acetonitrile because reaching the minimum of the gateway CT state is hampered by intramolecular hydrogen bonding, which enforces planarity between the aldehyde group and the aromatic group. Our study provides a comprehensive picture of the non-radiative decay of 5fdCyd in solution and new insights into the factors governing ISC in biomolecules. We propose that the intramolecular CT state observed here is a key to the excited-state dynamics of epigenetic nucleosides with modified exocyclic functional groups, paving the way to study their effects in DNA strands.

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.

Excited State Decay Pathways of 2'-Deoxy-5-methylcytidine and Deoxycytidine Revisited in Solution: A Comprehensive Kinetic Study by Femtosecond Transient Absorption

Methylated cytosine is proved to have an important role as an epigenetic signal in gene regulation and is often referred to "the fifth base of DNA". A comprehensive understanding of the electronic excited state relaxation in cytosine and its methylated derivatives is crucial for revealing UV-induced photodamage to the biological genome. Because of the existence of multiple closely lying "bright" and "dark" excited states, the decay pathways in these DNA nucleosides are the most complex and the least understood so far. In this study, femtosecond transient absorption with different excitation wavelengths (240-296 nm) was used to study the relaxation of excited electronic states of 5-methylcytosine (5mC) and 2'-deoxy-5-methylcytidine (5mdCyd) in phosphate buffered aqueous solution and in acetonitrile solution. Two distinct nonradiative decay channels were directly observed. The first one is a several picosecond internal conversion channel that involves two bright ππ* states (ππ*₂ and ππ*₁) when ππ*₂ state is initially populated. The second channel contains the lower energy ππ*₁ state and a so far experimental unidentified long-lived state which exhibits a several nanosecond lifetime. The long-lived state can only be accessed by the initially excited ππ*₁ state. Inspired by this new discovery in 5mC and 5mdCyd, we revisited the decay of excited state of 2'-deoxycytidine (dCyd), revealing very similar decay pathways. Additionally, a well-known dark n_Oπ* state (carbonyl lone pair) with ∼30 ps lifetime is present in both decay channels in dCyd. With our detailed experimental results, we successfully reconcile the long history debate of cytosine excited state relaxation mechanism by pointing out that the reason for the complex dynamics under traditional 266 nm excitation is mixed signals from the above-mentioned two distinct decay pathways. Our findings lead to a dramatically different and new picture of electronic energy relaxation in 5mdCyd/dCyd and could help to understand photostability as well as UV-induced photodamage of these nucleotides and related DNAs.

How nature covers its bases

The response of nucleobases to UV radiation depends on structure in subtle ways, as revealed by gas-phase experiments.