Relaxation Mechanisms of 5-Azacytosine

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

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

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

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

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

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

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

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

On the photorelease of nitric oxide by nitrobenzene derivatives: A CASPT2//CASSCF model

Nitroaromatic compounds can photorelease nitric oxide after UV absorption. The efficiency of the photoreaction depends on the molecular structure, and two features have been pointed out as particularly important for the yield of the process: the presence of methyl groups at the ortho position with respect to the nitro group and the degree of conjugation of the molecule. In this paper, we provide a theoretical characterization at the CASPT2//CASSCF (complete active space second-order perturbation theory//complete active space self-consistent field) level of theory of the photorelease of NO for four molecules derived from nitrobenzene through the addition of ortho methyl groups and/or the elongation of the conjugation. Our previously described mechanism obtained for the photorelease of NO in nitrobenzene has been adopted as a model for the process. According to this model, the process proceeds through a reactive singlet-triplet crossing (STC) region that the system can reach from the triplet 3(πOπ*) minimum. The energy barrier that must be surmounted in order to populate the reactive STC can be associated with the efficiency of the photoreaction. Here, the obtained results display clear differences in the efficiency of the photoreaction in the studied systems and can be correlated with experimental results. Thus, the model proves its ability to highlight the differences in the photoreaction efficiency for the nitroaromatic compounds studied here.

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

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

Building a Functionalizable, Potent Chemiluminescent Agent: A Rational Design Study on 6,8-Substituted Luminol Derivatives

Luminol is a prominent chemiluminescent (CL) agent, finding applications across numerous fields, including forensics, immunoassays, and imaging. Different substitution patterns on the aromatic ring can enhance or decrease its CL efficiency. We herein report a systematic study on the synthesis and photophysics of all possible 6,8-disubstituted luminol derivatives bearing H, Ph, and/or Me substituents. Their CL responses are monitored at three pH values (8, 10, and 12), thus revealing the architecture with the optimum CL efficiency. The most efficient pattern is used for the synthesis of a strongly CL luminol derivative, bearing a functional group for further, straightforward derivatization. This adduct exhibits an unprecedented increase in chemiluminescence efficiency at pH = 12, pH = 10, and especially at pH = 8 (closer to the biologically relevant conditions) compared to luminol. Complementary work on the fluorescence of the emissive species as well as quantum chemistry computations are employed for the rationalization of the observed results.

On the Intrinsically Low Quantum Yields of Pyrimidine DNA Photodamages: Evaluating the Reactivity of the Corresponding Minimum Energy Crossing Points

The low quantum yield of photoformation of cyclobutane pyrimidine dimers and pyrimidine-pyrimidone (6-4) adducts in DNA bases is usually associated with the presence of more favorable nonreactive decay paths and with the unlikeliness of exciting the system in a favorable conformation. Here, we prove that the ability of the reactive conical intersection to bring the system either back to the absorbing conformation or to the photoproduct must be considered as a fundamental factor in the low quantum yields of the mentioned photodamage. In support of the proposed model, the one order of magnitude difference in the quantum yield of formation of the cyclobutane thymine dimer with respect to the thymine-thymine (6-4) adduct is rationalized here by comparing the reactive ability of the seam of intersections leading respectively to the cyclobutane thymine dimer and the oxetane precursor of the thymine-thymine (6-4) adduct at the CASPT2 level of theory.

Photochemical Relaxation Pathways of 9 H-8-Azaguanine and 8 H-8-Azaguanine

The photochemical reaction path approach, the MS-CASPT2 quantum-chemical method, and double-ζ basis sets (cc-pVDZ) were used to study 9 H-8-azaguanine and 8 H-8-azaguanine relaxation pathways. Several potential energy hypersurfaces were characterized by means of equilibrium geometries, surface crossings (conical intersections and singlet-triplet intersystem crossings), minimum energy paths, and linear interpolation in internal coordinates. The 9 H-8-azaguanine main photochemical event begins with the direct population of the ¹(ππ* L_a) state, which evolves toward a conical intersection with the ground state after surmounting a small energy barrier, explaining why it is nonfluorescent. For 8 H-8-azaguanine, two relaxation mechanisms are possible, depending on the excitation energy. If the S₁ ¹(ππ*) state is initially populated (lower-energy region), the system evolves barrierless to the S₁ ¹(ππ*)_min region, from where the excess energy is released. If the ¹(ππ* L_a) state is populated (higher-energy radiation), the system will encounter conical intersections with the S₂ ¹(n_Oπ*) and S₁ ¹(ππ*) states before evolving to the ¹(ππ* L_a)_min region, from where a conical intersection with the ground state is accessible, favoring radiationless deactivation to the ground state. However, because a fraction of the population can be transferred from ¹(ππ* L_a) to the S₁ ¹(ππ*) state, emission from the S₁ ¹(ππ*)_min region is also expected, although weaker than it would be if the S₁ ¹(ππ*) state were populated directly. Irrespective of the excitation energy, the emissive state is the same and a single fluorescence band is observed, with the strongest emission occurring upon excitation in the lower-energy region, as observed experimentally. Therefore, our computational study corroborates experimental results, attributing the emission of the neutral form of 8-azaguanine in solution to the presence of the minor 8 H-8-azaguanine tautomer, while the 9 H-8-azaguanine major tautomer is nonfluorescent.

Similar chemical structures, dissimilar triplet quantum yields: a CASPT2 model rationalizing the trend of triplet quantum yields in nitroaromatic systems

S₁/S₀ accessibility strongly influences the triplet quantum yields of nitronaphthalenes.

Ab initio molecular dynamics relaxation and intersystem crossing mechanisms of 5-azacytosine

Nonadiabatic SHARC dynamics simulations reveal the molecular deformations involved in the photodeactivation pathways of 5-azacytosine.