Quantum Chemical Investigations on the Nonradiative Deactivation Pathways of Cytosine Derivatives

Quantum mechanics/molecular mechanics studies on mechanistic photophysics of epigenetic C5-halogenated DNA nucleosides: 2'-deoxy-5-chlorocytidine and 2'-deoxy-5-bromocytidine in aqueous solution

In this work, we have employed the high-level QM(CASPT2//CASSCF)/MM method to study the photophysical mechanisms of two important metabolized DNA/RNA nucleoside byproducts, i.e., 2'-deoxy-5-chlorocytidine (5CldCyd) and 2'-deoxy-5-bromocytidine (5BrdCyd), in aqueous solution. On the basis of our optimized minimum-energy structures, conical intersections, and crossing points, as well as the computed associated excited-state relaxation pathways involving the different internal conversion (IC) and intersystem crossing (ISC) processes in and between the S1, T1, T2, and S0 states, we have suggested the feasible excited-state relaxation mechanisms of these two important epigenetic halogenated DNA nucleosides. The initially populated spectroscopic bright 1ππ* state in the Franck-Condon (FC) region is the S1 state both for 5CldCyd and 5BrdCyd under 295 nm irradiation. The excited S1 state first evolves into its minimum S1-MIN and rapidly undergoes efficient IC to the S0 state via the nearby low-lying S1/S0 conical intersection. The corresponding energy barrier of the S1 → S0 IC path in 5CldCyd is estimated to be 4.6 kcal mol-1 at the QM(CASPT2)/MM level, while it is found to be an almost barrierless process in 5BrdCyd. In addition to this very efficient IC, the S1 state can partially slowly undergo ISC to transfer to the T1 state. Because the small spin-orbit couplings (SOCs) of S1/T1 and S1/T2 are estimated to be less than 5.0 cm-1 at the QM(CASPT2)/MM level, the ISC involved T1 formation is not so efficient. The resulting T1 state from the minor S1 → T1 and S1 → T2 → T1 ISCs will first relax to its minimum T1-MIN and continue to approach the nearby accessible T1/S0 crossing point, followed by further T1 → S0 ISC to the S0 state. Relatively, the T1 → S0 ISC of 5BrdCyd is significantly enhanced by a large T1/S0 SOC of 32.9 cm-1 at the T1/S0 crossing point. The present work rationalizes the excited-state dynamics of 5CldCyd and 5BrdCyd in aqueous solution and could provide mechanistic insights into understanding the photophysics of similar halogenated DNA nucleosides and their derivatives.

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

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.

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.

Analytic gradients for multiconfiguration pair-density functional theory with density fitting: Development and application to geometry optimization in the ground and excited states

Density fitting reduces the computational cost of both energy and gradient calculations by avoiding the computation and manipulation of four-index electron repulsion integrals. With this algorithm, one can efficiently optimize the geometries of large systems with an accurate multireference treatment. Here, we present the derivation of multiconfiguration pair-density functional theory for energies and analytic gradients with density fitting. Six systems are studied, and the results are compared to those obtained with no approximation to the electron repulsion integrals and to the results obtained by complete active space second-order perturbation theory. With the new approach, there is an increase in the speed of computation with a negligible loss in accuracy. Smaller grid sizes have also been used to reduce the computational cost of multiconfiguration pair-density functional theory with little effect on the optimized geometries and gradient values.

Quantum dynamics of the ππ*/nπ* decay of the epigenetic nucleobase 1,5-dimethyl-cytosine in the gas phase

We study the ultrafast dynamics of 1,5-dimethyl-cytosine, a model for 5-methyl-cytidine, after photoexcitation to the first two bright ππ* states, focusing on the possible population transfer to dark nπ* states. To that end we propagate the initial wave packets on the coupled potential energy surfaces of the seven lowest energy excited states modelled with a diabatic linear vibronic coupling (LVC) model, considering all the vibrational coordinates. Time-evolution is computed by the multilayer version of the multiconfigurational time dependent Hartree (ML-MCTDH) method. The LVC Hamiltonian is parametrized with time-dependent density functional theory (TD-DFT) calculations adopting PBE0 and CAM-B3LYP functionals, which provide a different energy gap between the lowest energy nπ* states and the spectroscopic ππ* state. Population of the lowest ππ* flows to a dark nπ* state which involves a lone pair (LP) of the carbonyl oxygen (nOπ*), but the extent of such transfer is much larger according to PBE0 than to CAM-B3LYP. Photoexcitation to the second bright state gives rise to much richer dynamics with an ultrafast (50 fs) complete decay to the lowest ππ*, to nOπ* and to another nπ* in which the excited electron comes from the LP of the ring nitrogen. We perform a detailed analysis of the vibronic dynamics both in terms of normal modes and valence coordinates (bond lengths and angles). The comparison with the analogous dynamics in 1-methyl-cytosine, a model for cytidine, provides insights into the effect of methylation at carbon 5 on the electronic and nuclear dynamics.

Analytic gradients for state-averaged multiconfiguration pair-density functional theory

Analytic gradients are important for efficient calculations of stationary points on potential energy surfaces, for interpreting spectroscopic observations, and for efficient direct dynamics simulations. For excited electronic states, as are involved in UV-Vis spectroscopy and photochemistry, analytic gradients are readily available and often affordable for calculations using a state-averaged complete active space self-consistent-field (SA-CASSCF) wave function. However, in most cases, a post-SA-CASSCF step is necessary for quantitative accuracy, and such calculations are often too expensive if carried out by perturbation theory or configuration interaction. In this work, we present the analytic gradients for multiconfiguration pair-density functional theory based on SA-CASSCF wave functions, which is a more affordable alternative. A test set of molecules has been studied with this method, and the stationary geometries and energetics are compared to values in the literature as obtained by other methods. Excited-state geometries computed with state-averaged pair-density functional theory have similar accuracy to those from complete active space perturbation theory at the second-order.

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.

Detection of DNA Modifications by Sequence-Specific Transcription Factors

The establishment, detection, and alteration or elimination of epigenetic DNA modifications are essential to controlling gene expression ranging from bacteria to mammals. The DNA methylations occurring at cytosine and adenine are carried out by SAM-dependent methyltransferases. Successive oxidations of 5-methylcytosine (5mC) by Tet dioxygenases generate 5-hydroxymethyl (5hmC), 5-formyl (5fC), and 5-carboxyl (5caC) derivatives; thus, DNA elements with multiple methylation sites can have a wide range of modification states. In contrast, oxidation of N6-methyladenine by homologs of Escherichia coli AlkB removes the methyl group directly. Both Tet and AlkB enzymes are 2-oxoglutarate- and Fe(II)-dependent dioxygenases. DNA-binding proteins decode the modification status of specific genomic regions. This article centers on two families of sequence-specific transcription factors: bZIP (basic leucine-zipper) proteins, exemplified by the AP-1 and CEBPβ recognition of 5mC; and bHLH (basic helix-loop-helix) proteins, exemplified by MAX and TCF4 recognition of 5caC. We discuss the impact of template strand DNA modification on the activities of DNA and RNA polymerases, and the varied tendencies of modifications to alter base pairing and their interactions with DNA repair enzymes.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Detecting and interpreting DNA methylation marks

The generation, alteration, recognition, and erasure of epigenetic modifications of DNA are fundamental to controlling gene expression in mammals. These covalent DNA modifications include cytosine methylation by AdoMet-dependent methyltransferases and 5-methylcytosine oxidation by Fe(II)-dependent and α-ketoglutarate-dependent dioxygenases. Sequence-specific transcription factors are responsible for interpreting the modification status of specific regions of chromatin. This review focuses on recent developments in characterizing the functional and structural links between the modification status of two DNA bases: 5-methylcytosine and 5-methyluracil (thymine).

Planarizing cytosine: The S1 state structure, vibrations, and nonradiative dynamics of jet-cooled 5,6-trimethylenecytosine

We measure the S₀ → S₁ spectrum and time-resolved S₁ state nonradiative dynamics of the "clamped" cytosine derivative 5,6-trimethylenecytosine (TMCyt) in a supersonic jet, using two-color resonant two-photon ionization (R2PI), UV/UV holeburning, and ns time-resolved pump/delayed ionization. The experiments are complemented with spin-component scaled second-order approximate coupled cluster (SCS-CC2), time-dependent density functional theory, and multi-state second-order perturbation-theory (MS-CASPT2) ab initio calculations. While the R2PI spectrum of cytosine breaks off ∼500 cm^-1 above its 0₀⁰ band, that of TMCyt extends up to +4400 cm^-1 higher, with over a hundred resolved vibronic bands. Thus, clamping the cytosine C⁵-C⁶ bond allows us to explore the S₁ state vibrations and S₀ → S₁ geometry changes in detail. The TMCyt S₁ state out-of-plane vibrations ν₁^', ν₃^', and ν₅^' lie below 420 cm^-1, and the in-plane ν₁₁^', ν₁₂^', and ν₂₃^' vibrational fundamentals appear at 450, 470, and 944 cm^-1. S₀  →  S₁ vibronic simulations based on SCS-CC2 calculations agree well with experiment if the calculated ν₁^', ν₃^', and ν₅^' frequencies are reduced by a factor of 2-3. MS-CASPT2 calculations predict that the ethylene-type S₁ ⇝ S₀ conical intersection (CI) increases from +366 cm^-1 in cytosine to >6000 cm^-1 in TMCyt, explaining the long lifetime and extended S₀ → S₁ spectrum. The lowest-energy S₁ ⇝ S₀ CI of TMCyt is the "amino out-of-plane" (OP_X) intersection, calculated at +4190 cm^-1. The experimental S₁ ⇝ S₀ internal conversion rate constant at the S₁(v^'=0) level is k_IC=0.98-2.2⋅10⁸ s^-1, which is ∼10 times smaller than in 1-methylcytosine and cytosine. The S₁(v^'=0) level relaxes into the T₁(³ππ*) state by intersystem crossing with k_ISC=0.41-1.6⋅10⁸ s^-1. The T₁ state energy is measured to lie 24 580±560 cm^-1 above the S₀ state. The S₁(v^'=0) lifetime is τ=2.9 ns, resulting in an estimated fluorescence quantum yield of Φ_fl=24%. Intense two-color R2PI spectra of the TMCyt amino-enol tautomers appear above 36 000 cm^-1. A sharp S₁ ionization threshold is observed for amino-keto TMCyt, yielding an adiabatic ionization energy of 8.114±0.002 eV.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

Remarkable effects of solvent and substitution on the photo-dynamics of cytosine: a femtosecond broadband time-resolved fluorescence and transient absorption study

Cytosine (Cyt) among all the nucleic acid bases features the most complex and least understood nonradiative deactivation, a process that is crucially important for its photostability. Herein, the excited state dynamics of Cyt and a series of its N1- and C5-derivatives, including the full set of Cyt nucleosides and nucleotides in DNA and RNA and the nucleosides of 5-methyl cytosine, 5-methylcytidine and 2'-deoxy-5-methylcytidine, have been investigated in water and in methanol employing femtosecond broadband time-resolved fluorescence coupled with fs transient absorption spectroscopy. The results reveal remarkable state-specific effects of the substitution and solvent in tuning distinctively the timescales and pathways of the nonradiative decays. For Cyt and the N1-derivatives, the nonradiative deactivations occur in a common two-state process through three channels, two from the light-absorbing ππ* state with respectively the sub-picosecond (∼0.2 ps) and the picosecond (∼1.5 ps) time constant, and the third is due to an optically dark nπ* state with the lifetime ranging from several to hundreds of picoseconds depending on solvents and substitutions. Compared to Cyt, the presence of the ribose or deoxyribose moiety at the N1 position of N1-derivatives facilitates the formation of the nπ* at the sub-picosecond timescale and at the same time increases its lifetime by ∼4-6 times in both water and methanol. In sharp contrast, the existence of the methyl group at the C5 position of the C5-derivatives eliminates completely the sub-picosecond ππ* channel and the channel due to the nπ*, but on the other hand slows down the decay of the ππ* state which after relaxation exhibits a single time constant of ∼4.1 to ∼7.6 ps depending on solvents. Varying the solvent from water to methanol accelerates only slightly the decay of the ππ* state in all the compounds; while for Cyt and its N1-derivatives, this change of solvent also retards strongly the nπ* channel, prolongs its lifetime from such as ∼7.7 ps in water to ∼52 ps in methanol for Cyt and from ∼30 ps in water to ∼186 ps in methanol for deoxycytidine. The spectral signatures we obtained for the ππ* and the nπ* states allow unambiguous evidence for clarifying uncertainties in the excited states of Cyt and the derivatives. The results provide a unifying experimental characterization at an improved level of detail about the photophysics of Cyt and its analogues under biologically relevant conditions and may help in understanding the photostability as well as photo-damages of the bases and related DNAs.

Gas-Phase Cytosine and Cytosine-N1-Derivatives Have 0.1-1 ns Lifetimes Near the S1 State Minimum

Ultraviolet radiative damage to DNA is inefficient because of the ultrafast S1 ⇝ S0 internal conversion of its nucleobases. Using picosecond pump-ionization delay measurements, we find that the S1((1)ππ*) state vibrationless lifetime of gas-phase keto-amino cytosine (Cyt) is τ = 730 ps or ∼ 700 times longer than that measured by femtosecond pump-probe ionization at higher vibrational excess energy, Eexc. N1-Alkylation increases the S1 lifetime up to τ = 1030 ps for N1-ethyl-Cyt but decreases it to 100 ps for N1-isopropyl-Cyt. Increasing the vibrational energy to Eexc = 300-550 cm(-1) decreases the lifetimes to 20-30 ps. The nonradiative dynamics of S1 cytosine is not solely a property of the amino-pyrimidinone chromophore but is strongly influenced by the N1-substituent. Correlated excited-state calculations predict that the gap between the S2((1)nOπ*) and S1((1)ππ*) states decreases along the series of N1-derivatives, thereby influencing the S1 state lifetime.