Resolving Ultrafast Photoinduced Deactivations in Water-Solvated Pyrimidine Nucleosides

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain onπ π * state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained byn π * trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of theπ π * state leading to a wrongπ π * / n π * energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light-matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. We present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes-Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

Photochemical Stability of 5-Methylcytidine Relative to Cytidine: Photophysical Insight for mRNA Therapeutic Applications

5-Methylcytidine (5mCyd) has recently been investigated with renewed interest in its utilization in mRNA therapeutics. However, its photostability following exposure to electromagnetic radiation has been overlooked. This Letter compares the photostability and excited-state dynamics of 5mCyd with those of the canonical RNA nucleoside, cytidine (Cyd), using steady-state and femtosecond transient absorption spectroscopy under physiologic conditions. 5mCyd is shown to have a 5-fold higher fluorescence yield and a 5-fold longer 1ππ* excited-state decay lifetime. Importantly, however, the excited-state population in 5mCyd decays primarily by internal conversion, with a photodegradation rate 3 times smaller than that in Cyd. In Cyd, the population of a 1nπ* state with a lifetime of ca. 45 ps is implicated in the formation 6-hydroxycytidine and other photoproducts.

Revealing the excited-state dynamics of cytidine and the role of excited-state proton transfer process

The high photostability of DNAs and RNAs is inextricably related to the photochemical and photophysical properties of their building blocks, nucleobases and nucleosides, which can dissipate the absorbed UV light energy in a harmless manner. The deactivation mechanism of the nucleosides, especially the decay pathways of cytidine (Cyd), has been a matter of intense debate. In the current study, we employ high-level electronic structure calculations combined with excited state non-adiabatic dynamic simulations to provide a clear picture of the excited state deactivation of Cyd in both gas phase and aqueous solution. In both environments, a barrierless decay path driven by the ring-puckering motion and a relaxation channel with a small energy barrier driven by the elongation motion of CO bond are assigned to <200 fs and sub-picosecond decay time component, respectively. The presence of ribose group has a subtle effect on the dynamic behavior of Cyd in gas phase as the ribose-to-base hydrogen/proton transfer process is energetically inaccessible with a sizable energy barrier of about 1.4 eV. However, this energy barrier is significantly reduced in water, especially when an explicit water molecule is present. Therefore, we argue that the long-lived decay channel found in aqueous solution could be assigned to the Cyd-water intermolecular hydrogen/proton transfer process. The present study postulates a novel scenario toward deep understanding the intrinsic photostability of DNAs and RNAs and provides solid evidence to disclose the long history debate of cytidine excited-state decay mechanism, especially for the assignment of experimentally observed time components.

Galván I, Rivalta I, Aquilante F, Garavelli M, Lindh R, Segarra-Martí J. Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

Direct Observation of Excitation Wavelength-Dependent Ultrafast Intersystem Crossing in Cytosine Nucleoside Solution

A triplet excited state can lead to different DNA photolesions, especially in cytosine and its nucleoside/nucleotide as they are hotspots for DNA mutations. However, the triplet state generation mechanism is in controversy, and experimental evidence of ultrafast intersystem crossing (ISC) has not been registered in these molecules. In this work, ultrafast ISC is directly observed in 2'-deoxycytidine (dCyd) solution by using femtosecond transient absorption spectroscopy. Surprisingly, we demonstrate that ISC in dCyd is sensitive to the excitation wavelength, and a spin-vibronic ISC mechanism is proposed. This finding is the last piece of the dCyd excited-state deactivation mechanism puzzle and sets the base for further investigation of triplet state-involved photophysics and photochemistry in dCyd-containing DNA.

Role of the Hydrogen Bond on the Internal Conversion of Photoexcited Adenosine

Experiments and theory have revealed that hydrogen bonds modify the excited-state lifetimes of nucleosides compared to nucleobases. Nevertheless, how these bonds impact the internal conversion is still unsettled. This work simulates the non-adiabatic dynamics of adenosine conformers in the gas phase with and without hydrogen bonds between the sugar and adenine moieties. The isomer containing the hydrogen bond (syn) exhibits a significantly shorter excited-state lifetime than the isomer without it (anti). However, internal conversion through electron-driven proton transfer between sugar and adenine plays only a minor (although non-negligible) role in the photophysics of adenosine. Either with or without hydrogen bonds, photodeactivation preferentially occurs following the ring-puckering pathways. The role of the hydrogen bond is to avoid the sugar rotation relative to adenine, shortening the distance to the ring-puckering internal conversion.

Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, namely, equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the ππ* states (one for pyrimidine, two for purines) associated with the lowest-energy absorption band, and for the close-lying nπ* state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effects by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the nπ* state. Finally, ESA spectra have been computed from the minima of the lowest energy ππ* state, and found in line with the available experimental transient absorption spectra of the nucleosides in solution, providing further validation of our computational approach.

On the Photostability of Cyanuric Acid and Its Candidature as a Prebiotic Nucleobase

Cyanuric acid is a triazine derivative that has been identified from reactions performed under prebiotic conditions and has been proposed as a prospective precursor of ancestral RNA. For cyanuric acid to have played a key role during the prebiotic era, it would have needed to survive the harsh electromagnetic radiation conditions reaching the Earth’s surface during prebiotic times (≥200 nm). Therefore, the photostability of cyanuric acid would have been crucial for its accumulation during the prebiotic era. To evaluate the putative photostability of cyanuric acid in water, in this contribution, we employed density functional theory (DFT) and its time-dependent variant (TD-DFT) including implicit and explicit solvent effects. The calculations predict that cyanuric acid has an absorption maximum at ca. 160 nm (7.73 eV), with the lowest-energy absorption band extending to ca. 200 nm in an aqueous solution and exhibiting negligible absorption at longer wavelengths. Excitation of cyanuric acid at 160 nm or longer wavelengths leads to the population of S_5,6 singlet states, which have ππ* character and large oscillator strengths (0.8). The population reaching the S_5,6 states is expected to internally convert to the S_1,2 states in an ultrafast time scale. The S_1,2 states, which have nπ* character, are predicted to access a conical intersection with the ground state in a nearly barrierless fashion (ca. ≤ 0.13 eV), thus efficiently returning the population to the ground state. Furthermore, based on calculated spin–orbit coupling elements of ca. 6 to 8 cm⁻¹, the calculations predict that intersystem crossing to the triplet manifold should play a minor role in the electronic relaxation of cyanuric acid. We have also calculated the vertical ionization energy of cyanuric acid at 8.2 eV, which predicts that direct one-photon ionization of cyanuric acid should occur at ca. 150 nm. Collectively, the quantum-chemical calculations predict that cyanuric acid would have been highly photostable under the solar radiation conditions reaching the Earth’s surface during the prebiotic era in an aqueous solution. Of relevance to the chemical origin of life and RNA-first theories, these observations lend support to the idea that cyanuric acid could have accumulated in large quantities during the prebiotic era and thus strengthens its candidature as a relevant prebiotic nucleobase.

Molecular Identification of the Transient Species Mediating the Deactivation Dynamics of Solvated Guanosine and Deazaguanosine

Small structural alterations of the purine/pyrimidine core have been related to important photophysical changes, such as the loss of photostability. Similarly to canonical nucleobases, solute-solvent interactions can lead to a change in the excited state lifetimes and/or to the interplay of different states in the photophysics of these modified nucleobases. To shed light on both effects, we here report a complete picture of the absorption spectra and excited state deactivation of deoxyguanosine and its closely related derivative, deoxydeazaguanosine, in water and methanol through the mapping of the excited state potential energy surfaces and molecular dynamics simulations at the TD-DFT level of theory. We show that the N by CH exchange in the imidazole ring of deoxyguanosine translates into a small red-shift of the bright states and slightly faster dynamics. In contrast, changing solvent from water to methanol implies the opposite, i.e., that the deactivation of both systems to the ground state is significantly hindered.

Increased Photostability of the Integral mRNA Vaccine Component N1 -Methylpseudouridine Compared to Uridine

N1 -Methylation of pseudouridine (m1 ψ) replaces uridine (Urd) in several therapeutics, including the Moderna and BioNTech-Pfizer COVID-19 vaccines. Importantly, however, it is currently unknown if exposure to electromagnetic radiation can affect the chemical integrity and intrinsic stability of m1 ψ. In this study, the photochemistry of m1 ψ is compared to that of uridine by using photoirradiation at 267 nm, steady-state spectroscopy, and quantum-chemical calculations. Furthermore, femtosecond transient absorption measurements are collected to delineate the electronic relaxation mechanisms for both nucleosides under physiologically relevant conditions. It is shown that m1 ψ exhibits a 12-fold longer 1 ππ* decay lifetime than uridine and a 5-fold higher fluorescence yield. Notably, however, the experimental results also demonstrate that most of the excited state population in both molecules decays back to the ground state in an ultrafast time scale and that m1 ψ is 6.7-fold more photostable than Urd following irradiation at 267 nm.

Tracking excited state decay mechanisms of pyrimidine nucleosides in real time

DNA owes its remarkable photostability to its building blocks-the nucleosides-that efficiently dissipate the energy acquired upon ultraviolet light absorption. The mechanism occurring on a sub-picosecond time scale has been a matter of intense debate. Here we combine sub-30-fs transient absorption spectroscopy experiments with broad spectral coverage and state-of-the-art mixed quantum-classical dynamics with spectral signal simulations to resolve the early steps of the deactivation mechanisms of uridine (Urd) and 5-methyluridine (5mUrd) in aqueous solution. We track the wave packet motion from the Franck-Condon region to the conical intersections (CIs) with the ground state and observe spectral signatures of excited-state vibrational modes. 5mUrd exhibits an order of magnitude longer lifetime with respect to Urd due to the solvent reorganization needed to facilitate bulky methyl group motions leading to the CI. This activates potentially lesion-inducing dynamics such as ring opening. Involvement of the 1nπ* state is found to be negligible.

Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form

Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.

Simulation and Analysis of the Transient Absorption Spectrum of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Acetonitrile

4-(N,N-Dimethylamino)benzonitrile (DMABN) is a well-known model compound for dual fluorescence-in sufficiently polar solvents, it exhibits two distinct fluorescence emission bands. The interpretation of its transient absorption (TA) spectrum in the visible range is the subject of a long-standing controversy. In the present study, we resolve this issue by calculating the TA spectrum on the basis of nonadiabatic molecular dynamics simulations. An unambiguous assignment of spectral signals to specific excited-state structures is achieved by breaking down the calculated spectrum into contributions from twisted and nontwisted molecular geometries. In particular, the much-discussed excited-state absorption band near 1.7 eV (ca. 700 nm) is attributed to the near-planar locally excited (LE) minimum on the S1 state. On the technical side, our study demonstrates that the second-order approximate coupled cluster singles and doubles (CC2) method can be used successfully to calculate the TA spectra of moderately large organic molecules, provided that the system in question does not approach a crossing between the lowest excited state and the singlet ground state within the time frame of the simulation.

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.

Broadband cavity-enhanced ultrafast spectroscopy

Broadband ultrafast optical spectroscopy methods, such as transient absorption spectroscopy and 2D spectroscopy, are widely used to study molecular dynamics. However, these techniques are typically restricted to optically thick samples, such as solids and liquid solutions. In this article we discuss a cavity-enhanced ultrafast transient absorption spectrometer covering almost the entire visible range with a detection limit of ΔOD < 1 × 10-9, extending broadband all-optical ultrafast spectroscopy techniques to dilute beams of gas-phase molecules and clusters. We describe the technical innovations behind the spectrometer and present transient absorption data on two archetypical molecular systems for excited-state intramolecular proton transfer, 1'-hydroxy-2'-acetonapthone and salicylideneaniline, under jet-cooled and Ar cluster conditions.

Excited-State Absorption of Uracil in the Gas Phase: Mapping the Main Decay Paths by Different Electronic Structure Methods

We present a computational study of the one-photon and excited-state absorption (ESA) from the two lowest energy excited states of uracil in the gas phase: an nπ* dark state (1n) and the lowest energy bright ππ* state (1π). The predictions of six different linear response electronic structure methods, namely, TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, ADC(2), ADC(2)-x, and ADC(3) are critically compared. In general, the spectral shapes predicted by TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, and ADC(3) are fairly similar, though the quality of TD-CAM-B3LYP slightly deteriorates in the high-energy region. By computing the spectra at some key structures on different potential energy surfaces (PES), that is, the Franck-Condon point, the 1n minimum, and structures representative of different regions of the 1π PES, we obtain important insights into the shift of the ESA spectra, following the motion of the wavepacket on the excited-state PES. Though 1π has larger ESA than 1n, some spectral regions are dominated by these latter signals. Aside from its methodological interest, we thus obtain interesting indications to interpret transient absorption spectra to disentangle the photoactivated dynamics of nucleobases.

Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study

The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (π_Hπ_L^*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (n_Oπ_L^*). The population of the second ππ* state decays almost completely, in <100 fs, not only to π_Hπ_L^* and to n_Oπ_L^* states but also to another nπ_L^* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.

Selenium substitution effects on excited-state properties and photophysics of uracil: a MS-CASPT2 study

The photophysics of selenium-substituted nucleobases has attracted recent experimental attention because they could serve as potential photosensitizers in photodynamic therapy. Herein, we present a comprehensive MS-CASPT2 study on the spectroscopic and excited-state properties, and photophysics of 2-selenouracil (2SeU), 4-selenouracil (4SeU), and 2,4-selenouracil (24SeU). Relevant minima, conical intersections, crossing points, and excited-state relaxation paths in the lowest five electronic states (i.e., S₀, S₁, S₂, T₂, and T₁) are explored. On the basis of these results, their photophysical mechanisms are proposed. Upon photoirradiation to the bright S₂ state, 2SeU quickly relaxes to its S₂ minimum and then moves in an essentially barrierless way to a nearby S₂/S₁ conical intersection near which the S₁ state is populated. Next, the S₁ system arrives at an S₁/T₂/T₁ intersection where a large S₁/T₁ spin-orbit coupling of 430.8 cm^-1 makes the T₁ state populated. In this state, a barrier of 6.8 kcal mol^-1 will trap 2SeU for a while. In parallel, for 4SeU or 24SeU, the system first relaxes to the S₂ minimum and then overcomes a small barrier to approach an S₂/S₁ conical intersection. Once hopping to the S₁ state, there exists an extended region with very close S₁, T₂, and T₁ energies. Similarly, a large S₁/T₁ spin-orbit coupling of 426.8 cm^-1 drives the S₁→ T₁ intersystem crossing process thereby making the T₁ state populated. Similarly, an energy barrier heavily suppresses electronic transition to the S₀ state. The present work manifests that different selenium substitutions on uracil can lead to a certain extent of different vertical and adiabatic excitation energies, excited-state properties, and relaxation pathways. These insights could help understand the photophysics of selenium-substituted nucleobases.

Locating Cytosine Conical Intersections by Laser Experiments and Ab Initio Calculations

The decay mechanism of S₀ → S₁ excited cytosine (Cyt) and the effect of substitution are studied combining jet-cooled spectroscopy (nanosecond resonant two-photon ionization (R2PI) and picosecond lifetime measurements) with CASPT2//CASSCF computations for eight derivatives. For Cyt and five derivatives substituted at N1, C5, and C6, rapid internal conversion sets in at 250-1200 cm^-1 above the 0₀⁰ bands. The break-off in the spectra correlates with the calculated barriers toward the "C5-C6 twist" conical intersection, which unambiguously establishes the decay mechanism at low S₁ state vibrational energies. The barriers increase with substituents that stabilize the charge shifts at C4, C5, and C6 following (¹ππ*) excitation. The R2PI spectra of the clamped derivatives 5,6-trimethyleneCyt (TMCyt) and 1-methyl-TMCyt (1M-TMCyt), which decay along an N3 out-of-plane coordinate, extend up to +3500 and +4500 cm^-1.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Modeling multidimensional spectral lineshapes from first principles: application to water-solvated adenine

We theoretically describe spectral lineshape from first principles, providing insight into solvent–solute interactions in terms of static and dynamic disorder and how these shape experimental signals in linear and non-linear optical spectroscopies.

A multi-scale time-resolved study of photoactivated dynamics in 5-benzyl uracil, a model for DNA/protein interactions

We combine fluorescence up-conversion and time correlated single photon counting experiments to investigate the 5-benzyl uracil excited state dynamics in methanol from 100 fs up to several ns. This molecule has been proposed as a model for DNA/protein interactions. Our results show emission bands at about 310 and 350 nm that exhibit bi-exponential sub-ps decays. Calculations, including solvent effects by a mixed discrete-continuum model, indicate that the Franck Condon region is characterized by significant coupling between the excited states of the benzyl and the uracil moieties, mirrored by the short-lived emission at 310 nm. Two main ground state recovery pathways are identified, both contributing to the 350 nm emission. The first 'photophysical' decay path involves a ππ* excited state localized on the uracil and is connected to the ground electronic state by an easily accessible crossing with S0, accounting for the short lifetime component. Simulations indicate that a possible second pathway is characterized by exciplex formation, with partial benzene → uracil charge transfer character, that may lead instead to photocyclization. The relevance of our results is discussed in view of the photoactivated dynamics of DNA/protein complexes, with implications on their interaction mechanisms.

Relaxation Dynamics of Hydrated Thymine, Thymidine, and Thymidine Monophosphate Probed by Liquid Jet Time-Resolved Photoelectron Spectroscopy

The relaxation dynamics of thymine and its derivatives thymidine and thymidine monophosphate are studied using time-resolved photoelectron spectroscopy applied to a water microjet. Two absorption bands are studied; the first is a bright ππ* state which is populated using tunable-ultraviolet light in the range 4.74-5.17 eV and probed using a 6.20 eV probe pulse. By reversing the order of these pulses, a band containing multiple ππ* states is populated by the 6.20 eV pulse and the lower energy pulse serves as the probe. The lower lying ππ* state is found to decay in ∼400 fs in both thymine and thymidine independent of pump photon energy, while thymidine monophosphate decays vary from 670 to 840 fs with some pump energy dependence. The application of a computational quantum mechanical/molecular mechanical scheme at the XMS-CASPT2//CASSCF/AMBER level of theory suggests that conformational differences existing between thymidine and thymidine monophosphate in solution account for this difference. The higher lying ππ* band is found to decay in ∼600 fs in all three cases, but it is only able to be characterized when the 5.17 eV probe pulse is used. Notably, no long-lived signal from an nπ* state can be identified in either experiment on any of the three molecules.

Luminescent and DFT Study of Keto-Enol Tautomers of 5-Fluorouracil and Its Derivatives in Aqueous Solutions

Fluorescence (FL) spectra were recorded and identified for the following 5-fluorouracil (FU) tautomers: 5-fluoropyrimidine-2,4(1H,3H)-dione (A), 5-fluoro-2-hydroxypyrimidine-4(3H)-one (B), and 5-fluoro-4-hydroxypyrimidine-2(1H)-one (D), as well as the corresponding tautomers A and D of 1-(tetrahydrofuranyl-2)-5-fluorouracil and 1-methyl-5-fluorouracil, tautomers A and B of 3-methyl-5-fluorouracil, and the diketo tautomer of 1,3-dimethyl-5-fluorouracil. It was shown that the FL of rare tautomers of FU derivatives occurs by the excitation of uracil homoassociates followed by intramolecular proton transfer (IPT), formation of a pair of rare tautomers, and radiative deactivation of one of them. The FL quantum yields φ_i were estimated. In the case of rare tautomers, the φ_i value includes the probability of IPT in the homoassociate, and density functional theory calculations in the TPSSTPSS/6-311+G(d,p) approximation were used to determine the relative energies of tautomers, associated FU forms, and their spectral properties that agree with the observed luminescence pattern.

Excited-State Dynamics in the RNA Nucleotide Uridine 5'-Monophosphate Investigated Using Femtosecond Broadband Transient Absorption Spectroscopy

Damage to RNA from ultraviolet radiation induces chemical modifications to the nucleobases. Unraveling the excited states involved in these reactions is essential; however, investigations aimed at understanding the electronic-energy relaxation pathways of the RNA nucleotide uridine 5'-monophosphate (UMP) have not received enough attention. In this Letter, the excited-state dynamics of UMP is investigated in aqueous solution. Excitation at 267 nm results in a trifurcation event that leads to the simultaneous population of the vibrationally excited ground state, a long-lived ¹nπ* state, and a receiver triplet state within 200 fs. The receiver state internally converts to the long-lived ³ππ* state in an ultrafast time scale. The results elucidate the electronic relaxation pathways and clarify earlier transient absorption experiments performed for uracil derivatives in solution. This mechanistic information is important because long-lived nπ* and ππ* excited states of both singlet and triplet multiplicities are thought to lead to the formation of harmful photoproducts.

Conical intersection dynamics of pyrimidine nucleosides tracked with sub-20-fs UV pulses

By combining transient absorption spectroscopy with sub-20-fs UV pulses and ab initio numerical simulations we follow the ultrafast dynamics in pyrimidine nucleosides and visualize the passage through conical intersections presiding excited state deactivation.

UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level

The most relevant 'dark' electronic excited states in DNA/RNA pyrimidine nucleosides are mapped in water employing hybrid MS-CASPT2/MM optimisations with explicit solvation and including the sugar. Conical intersections (CIs) between initially accessed bright ¹ππ* and the lowest energy dark ¹nπ* excited states, involving the lone pair localised on the oxygen and/or nitrogen atoms are characterised. They are found in the vicinities of the Franck-Condon (FC) region and are shown to facilitate non-adiabatic population transfer. The excited state population of the ¹n_Oπ* state, localised in the carbonyl moiety on all pyrimidine nucleosides, is predicted to rapidly evolve to its minimum, displaying non-negligible potential energy barriers along its non-radiative decay, and accounting for the ps signal registered in pump-probe experiments as well as for an efficient population of the triplet state. Cytidine displays an additional ¹n_Nπ* state localised in the N3 atom and that leads to its excited state minimum displaying large potential energy barriers in the pathway connecting to the CI with the ground state. Sugar-to-base hydrogen/proton transfer processes are assessed in solution for the first time, displaying a sizable barrier along its decay and thus being competitive with other slow decay channels in the ps and ns timescales. A unified deactivation scheme for the long-lived channels of pyrimidine nucleosides is delivered, where the ¹n_Oπ* state is found to mediate the long-lived decay in the singlet manifold and act as the doorway for triplet population and thus accounting for the recorded phosphorescence and, more generally, for the transient/photoelectron spectral signals registered up to the ns timescale.

Manipulating azobenzene photoisomerization through strong light-molecule coupling

The formation of hybrid light–molecule states (polaritons) offers a new strategy to manipulate the photochemistry of molecules. To fully exploit its potential, one needs to build a toolbox of polaritonic phenomenologies that supplement those of standard photochemistry. By means of a state-of-the-art computational photochemistry approach extended to the strong-coupling regime, here we disclose various mechanisms peculiar of polaritonic chemistry: coherent population oscillations between polaritons, quenching by trapping in dead-end polaritonic states and the alteration of the photochemical reaction pathway and quantum yields. We focus on azobenzene photoisomerization, that encompasses the essential features of complex photochemical reactions such as the presence of conical intersections and reaction coordinates involving multiple internal modes. In the strong coupling regime, a polaritonic conical intersection arises and we characterize its role in the photochemical process. Our chemically detailed simulations provide a framework to rationalize how the strong coupling impacts the photochemistry of realistic molecules.

Manipulation of the photochemistry of molecules is traditionally achieved through synthetic chemical modifications. Here the authors use computational photochemistry to show how to control azobenzene photoisomerization through hybrid light–molecule states (polaritons).

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.

Decay Pathways of Thymine Revisited

The decay of electronically excited states of thymine (Thy) and thymidine 5'-monophosphate (TMP) was studied by time-resolved UV/vis and IR spectroscopy. In addition to the well-established ultrafast internal conversion to the ground state, a so far unidentified UV-induced species is observed. In D₂O, this species decays with a time constant of 300 ps for thymine and of 1 ns for TMP. The species coexists with the lowest triplet state and is formed with a comparably high quantum yield of about 10% independent of the solvent. The experimentally determined spectral signatures are discussed in the light of quantum chemical calculations of the singlet and triplet excited states of thymine.

Evolution of the Excitonic State of DNA Stacked Thymines: Intrabase ππ* → S0 Decay Paths Account for Ultrafast (Subpicosecond) and Longer (>100 ps) Deactivations

Monomer-like ring puckering decay paths for two stacked quantum mechanical thymines inside a solvated DNA duplex described at the molecular mechanics level are mapped using a hybrid CASPT2//CASSCF/MM protocol that accounts for steric, electronic and electrostatic interactions within the nucleobases native environment. Asymmetric stacking between nucleobases open ups different intrabase ππ* decay paths accounting for distinctive excited state lifetimes, spanning the subps to subns time window.

UV-Light-Induced Vibrational Coherences: The Key to Understand Kasha Rule Violation in trans-Azobenzene

We combine sub-20 fs transient absorption spectroscopy with state-of-the-art computations to study the ultrafast photoinduced dynamics of trans-azobenzene (AB). We are able to resolve the lifetime of the ππ* state, whose decay within ca. 50 fs is correlated to the buildup of the nπ* population and to the emergence of coherences in the dynamics, to date unobserved. Nonlinear spectroscopy simulations call for the CNN in-plane bendings as the active modes in the subps photoinduced coherent dynamics out of the ππ* state. Radiative to kinetic energy transfer into these modes drives the system to a high-energy planar nπ*/ground state conical intersection, inaccessible upon direct excitation of the nπ* state, that triggers an ultrafast (0.45 ps) nonproductive decay of the nπ* state and is thus responsible for the observed Kasha rule violation in UV excited trans-AB. On the other hand, cis-AB is built only after intramolecular vibrational energy redistribution and population of the NN torsional mode.

Two-dimensional electronic spectroscopy as a tool for tracking molecular conformations in DNA/RNA aggregates

A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows characterising ground state conformations of flexible nucleobase aggregates that play a crucial role in nucleic acid photochemistry.

Excited-State Dynamics of Isocytosine: A Hybrid Case of Canonical Nucleobase Photodynamics

We present resonant two-photon ionization (R2PI) spectra of isocytosine (isoC) and pump-probe results on two of its tautomers. IsoC is one of a handful of alternative bases that have been proposed in scenarios of prebiotic chemistry. It is structurally similar to both cytosine (C) and guanine (G). We compare the excited-state dynamics with the Watson-Crick (WC) C and G tautomeric forms. These results suggest that the excited-state dynamics of WC form of G may primarily depend on the heterocyclic substructure of the pyrimidine moiety, which is chemically identical to isoC. For WC isoC we find a single excited-state decay with a rate of ∼10¹⁰ s^-1, while the enol form has multiple decay rates, the fastest of which is 7 times slower than for WC isoC. The excited-state dynamics of isoC exhibits striking similarities with that of G, more so than with the photodynamics of C.

Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments

Interest in understanding the photophysics and photochemistry of thiated nucleobases has been awakened because of their possible involvement in primordial RNA or their potential use as photosensitizers in medicinal chemistry. The interpretation of the photodynamics of these systems, conditioned by their intricate potential energy surfaces, requires the powerful interplay between experimental measurements and state of the art molecular simulations. In this review, we provide an overview on the photophysics of natural nucleobases’ thioanalogs, which covers the last 30 years and both experimental and computational contributions. For all the canonical nucleobase’s thioanalogs, we have compiled the main steady state absorption and emission features and their interpretation in terms of theoretical calculations. Then, we revise the main topographical features, including stationary points and interstate crossings, of their potential energy surfaces based on quantum mechanical calculations and we conclude, by combining the outcome of different spectroscopic techniques and molecular dynamics simulations, with the mechanism by which these nucleobase analogs populate their triplet excited states, which are at the origin of their photosensitizing properties.

Water–chromophore electron transfer determines the photochemistry of cytosine and cytidine

Irradiation of aqueous cytidine may result in a water–chromophore electron transfer process which explains numerous experimental observations, including photodamage.