Dynamics of photoinduced processes in adenine and thymine base pairs

Correlated rotational alignment spectroscopy: a new tool for high-resolution spectroscopy and the analysis of heterogeneous samples

Correlated rotational alignment spectroscopy correlates observables of ultrafast gas-phase spectroscopy with high-resolution, broad-band rotational Raman spectra. This article reviews the measurement principle of CRASY, existing implementations for mass-correlated measurements, and the potential for future developments. New spectroscopic capabilities are discussed in detail: signals for individual sample components can be separated even in highly heterogeneous samples. Isotopologue rotational spectra can be observed at natural isotope abundance. Fragmentation channels are readily assigned in molecular and cluster mass spectra. And finally, rotational Raman spectra can be measured with sub-MHz resolution, an improvement of several orders-of-magnitude as compared to preceding experiments.

Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT

Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

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.

Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects

We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.

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.

A Plausible Mechanism of Uracil Photohydration Involves an Unusual Intermediate

It is well-known that photolysis of pyrimidine nucleobases, such as uracil, in an aqueous environment results in the formation of hydrate as one of the main products. Although several hypotheses regarding photohydration have been proposed in the past, e.g., the zwitterionic and "hot" ground-state mechanisms, its detailed mechanism remains elusive. Here, theoretical nonadiabatic simulations of the uracil photodynamics reveal the formation of a highly energetic but kinetically stable intermediate that features a half-chair puckered pyrimidine ring and a strongly twisted intracyclic double bond. The existence and the kinetic stability of the intermediate are confirmed by a variety of computational chemistry methods. According to the simulations, the unusual intermediate is mainly formed almost immediately (∼50-200 fs) upon photoabsorption and survives long enough to engage in a hydration reaction with a neighboring water. A plausible mechanism of uracil photohydration is proposed on the basis of the modeling of nucleophilic insertion of water into the twisted double bond of the intermediate.

Excited-state proton transfer in the 2-aminopyridine dimer: A surface hopping study

We present a theoretical investigation of the excited-state intermolecular proton transfer process in the 2-aminopyridine dimer. Prior experimental and theoretical studies on this doubly hydrogen bonded system have attributed an...

Adenine ultrafast photorelaxation via electron-driven proton transfer

Photorelaxation of adenine in water was reported to be ultrafast (within 180 fs) primarily due to radiationless relaxation. However, in the last two decades, several experimental and theoretical investigations on photoexcitation of adenine have revealed diverse types of decay mechanisms. Using time-dependent density functional excited-state nonadiabatic dynamics simulations we show that it is the water to adenine electron-driven proton transfer (EDPT) barrierless reaction responsible for the ultrafast component of the adenine relaxation, which, however, occurred only in the case of the 7H isomer of adenine with five water molecules. This result reveals a known reaction pathway, however not found in previous simulations, with inference for the ultrafast relaxation mechanisms of adenine reported in experiments. The 9H isomer of adenine with six water molecules relaxing in a water cluster followed the previously known structural distortion (C2) decay pathway. The observations of the adenine EDPT reaction with water provide the origin of the experimental ultrafast adenine decay component and present a possible method to tackle future computational challenges in molecular-level biological processes.

Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S₁-trapping, S₂-trapping, and S₁&S₂-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S₂-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S₁-trapping mechanism with two lifetimes, τ₁ = 30 ± 1 fs and τ₂ = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S₂ (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S₂→S₁ internal conversion and resides around the minimum on the S₁ (n_Oπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S₂→S₁ internal conversion is mediated by fast bond length alternation distortion, the subsequent S₁→S₀ occurs through several conical intersections, involving a slow puckering motion.

Excitation energy transport in DNA modelled by multi-chromophoric field-induced surface hopping

Revealing the extended excited state lifetime due to excitation energy transport in DNA by multi-chromophoric field-induced surface-hopping (McFISH).

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.

Internal conversion and intersystem crossing pathways in UV excited, isolated uracils and their implications in prebiotic chemistry

The photodynamic properties of molecules determine their ability to survive in harsh radiation environments. As such, the photostability of heterocyclic aromatic compounds to electromagnetic radiation is expected to have been one of the selection pressures influencing the prebiotic chemistry on early Earth. In the present study, the gas-phase photodynamics of uracil, 5-methyluracil (thymine) and 2-thiouracil-three heterocyclic compounds thought to be present during this era-are assessed in the context of their recently proposed intersystem crossing pathways that compete with internal conversion to the ground state. Specifically, time-resolved photoelectron spectroscopy measurements evidence femtosecond to picosecond timescales for relaxation of the singlet (1)ππ* and (1)nπ* states as well as for intersystem crossing to the triplet manifold. Trapping in the excited triplet state and intersystem crossing back to the ground state are investigated as potential factors contributing to the susceptibility of these molecules to ultraviolet photodamage.

New Insights into the State Trapping of UV-Excited Thymine

After UV excitation, gas phase thymine returns to a ground state in 5 to 7 ps, showing multiple time constants. There is no consensus on the assignment of these processes, with a dispute between models claiming that thymine is trapped either in the first (S₁) or in the second (S₂) excited states. In the present study, a nonadiabatic dynamics simulation of thymine is performed on the basis of ADC(2) surfaces, to understand the role of dynamic electron correlation on the deactivation pathways. The results show that trapping in S₂ is strongly reduced in comparison to previous simulations considering only non-dynamic electron correlation on CASSCF surfaces. The reason for the difference is traced back to the energetic cost for formation of a CO π bond in S₂.

Sunlight as an energetic driver in the synthesis of molecules necessary for life

This review considers how photochemistry and sunlight-driven reactions can abiotically generate prebiotic molecules necessary for the evolution of life.

Computational modeling of photoexcitation in DNA single and double strands

The photoexcitation of DNA strands triggers extremely complex photoinduced processes, which cannot be understood solely on the basis of the behavior of the nucleobase building blocks. Decisive factors in DNA oligomers and polymers include collective electronic effects, excitonic coupling, hydrogen-bonding interactions, local steric hindrance, charge transfer, and environmental and solvent effects. This chapter surveys recent theoretical and computational efforts to model real-world excited-state DNA strands using a variety of established and emerging theoretical methods. One central issue is the role of localized vs delocalized excitations and the extent to which they determine the nature and the temporal evolution of the initial photoexcitation in DNA strands.

Early events of DNA photodamage

Ultraviolet (UV) radiation is a leading external hazard to the integrity of DNA. Exposure to UV radiation triggers a cascade of chemical reactions, and many molecular products (photolesions) have been isolated that are potentially dangerous for the cellular system. The early steps that take place after UV absorption by DNA have been studied by ultrafast spectroscopy. The review focuses on the evolution of excited electronic states, the formation of photolesions, and processes suppressing their formation. Emphasis is placed on lesions involving two thymine bases, such as the cyclobutane pyrimidine dimer, the (6-4) lesion, and its Dewar valence isomer.

Solvent effects on the ultrafast nonradiative deactivation mechanisms of thymine in aqueous solution: excited-state QM/MM molecular dynamics simulations

On-the-fly excited-state quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations of thymine in aqueous solution are performed to investigate the role of solvent water molecules on the nonradiative deactivation process. The complete active space second-order perturbation theory (CASPT2) method is employed for a thymine molecule as the QM part in order to provide a reliable description of the excited-state potential energies. It is found that, in addition to the previously reported deactivation pathway involving the twisting of the C-C double bond in the pyrimidine ring, another efficient deactivation pathway leading to conical intersections that accompanies the out-of-plane displacement of the carbonyl group is observed in aqueous solution. Decay through this pathway is not observed in the gas phase simulations, and our analysis indicates that the hydrogen bonds with solvent water molecules play a key role in stabilizing the potential energies of thymine in this additional decay pathway.

Photoinduced processes in nucleic acids

Photoinduced processes in nucleic acids are phenomena of fundamental interest in diverse fields, from prebiotic studies, through medical research on carcinogenesis, to the development of bioorganic photodevices. In this contribution we survey many aspects of the research across the boundaries. Starting from a historical background, where the main milestones are identified, we review the main findings of the physical-chemical research of photoinduced processes on several types of nucleic-acid fragments, from monomers to duplexes. We also discuss a number of different issues which are still under debate.

Recent advances in experimental techniques to probe fast excited-state dynamics in biological molecules in the gas phase: dynamics in nucleotides, amino acids and beyond

In many chemical reactions, an activation barrier must be overcome before a chemical transformation can occur. As such, understanding the behaviour of molecules in energetically excited states is critical to understanding the chemical changes that these molecules undergo. Among the most prominent reactions for mankind to understand are chemical changes that occur in our own biological molecules. A notable example is the focus towards understanding the interaction of DNA with ultraviolet radiation and the subsequent chemical changes. However, the interaction of radiation with large biological structures is highly complex, and thus the photochemistry of these systems as a whole is poorly understood. Studying the gas-phase spectroscopy and ultrafast dynamics of the building blocks of these more complex biomolecules offers the tantalizing prospect of providing a scientifically intuitive bottom-up approach, beginning with the study of the subunits of large polymeric biomolecules and monitoring the evolution in photochemistry as the complexity of the molecules is increased. While highly attractive, one of the main challenges of this approach is in transferring large, and in many cases, thermally labile molecules into vacuum. This review discusses the recent advances in cutting-edge experimental methodologies, emerging as excellent candidates for progressing this bottom-up approach.

Photoreaction channels of the guanine-cytosine base pair explored by long-range corrected TDDFT calculations

Photoinduced processes in the Watson-Crick guanine-cytosine base pair are comprehensively studied by means of long-range corrected (LC) TDDFT calculations of potential energy profiles using the LC-BLYP and CAM-B3LYP functionals. The ab initio CC2 method and the conventional TDDFT method with the B3LYP functional are also employed to assess the reliability of the LC-TDDFT method. The present approach allows us to compare the potential energy profiles at the same computational level for excited-state reactions of the base pair, including single and double proton transfer between the bases and nonradiative decay via ring puckering in each base. In particular, long-range correction to the TDDFT method is critical for a qualitatively correct description of the proton transfer reactions. The calculated energy profiles exhibit low barriers for out-of-plane deformation of the guanine moiety in the locally-excited state, which is expected to lead to a conical intersection with the ground state, as well as for single proton transfer from guanine to cytosine with the well-known electron-driven proton transfer mechanism. Thus the present results suggest that both processes can compete in hydrogen-bonded base pairs and play a significant role in the mechanism of photostability.

Excited-state energies and electronic couplings of DNA base dimers

The singlet excited electronic states of two pi-stacked thymine molecules and their splittings due to electronic coupling have been investigated with a variety of computational methods. Focus has been given on the effect of intermolecular distance on these energies and couplings. Single-reference methods, CIS, CIS(2), EOM-CCSD, TDDFT, and the multireference method CASSCF, have been used, and their performance has been compared. It is found that the excited-state energies are very sensitive to the applied method but the couplings are not as sensitive. Inclusion of diffuse functions in the basis set also affects the excitation energies significantly but not the couplings. TDDFT is inadequate in describing the states and their coupling, while CIS(2) gives results very similar to EOM-CCSD. Excited states of cytosine and adenine pi-stacked dimers were also obtained and compared with those of thymine dimers to gain a more general picture of excited states in pi-stacked DNA base dimers. The coupling is very sensitive to the relative position and orientation of the bases, indicating great variation in the degree of delocalization of the excited states between stacked bases in natural DNA as it fluctuates.

Photoinduced nonadiabatic dynamics of pyrimidine nucleobases: on-the-fly surface-hopping study with semiempirical methods

The photoinduced relaxation dynamics of pyrimidine nucleobases (uracil, thymine, and cytosine) was studied using the surface-hopping approach at the semiempirical OM2/MRCI level of theory. The relevant potential energy surfaces were characterized by performing geometry optimizations of the energy minima of the lowest electronic states and of the most important conical intersections and by computing excitation energies at each configuration. Surface-hopping molecular dynamics simulations were performed to describe the nonadiabatic dynamics after excitation into the optically active state. In each of the molecules, the two lowest excited singlet states are involved in the dynamics, and there are competing relaxation paths. The dynamics is dominated by a two-step relaxation mechanism in uracil and thymine, while the direct decay to the ground state is most important in cytosine. For all three molecules, the simulations yield ultrafast S(2)-S(1) deexcitation within 50 fs and internal conversion to the ground state in less than 1 ps, consistent with recent experimental results from time-resolved photoelectron spectroscopy.