Effect of substituents on the excited-state dynamics of the modified DNA bases 2,4-diaminopyrimidine and 2,6-diaminopurine

The copious photochemistry of 2,6-diaminopurine: Luminescence, triplet population, and ground state recovery

9H- and 7H-2,6-Diaminopurine (26DAP) photoinduced events in vacuum were studied at the MS-CASPT2/cc-pVDZ level of theory. The S1 1 (ππ* La ) state is initially populated evolving barrierless towards its minimum energy structure, from where two photochemical events can take place in both tautomers. The first is the return of the electronic population to the ground state via the C6 conical intersection (CI-C6). The second involves an internal conversion to the ground through the C2 conical intersection (CI-C2). According to our geodesic interpolated paths connecting the critical structures, the second route is less favorable in both tautomers, due to the presence of high energy barriers. Our calculations suggest a competition between fluorescence and ultrafast relaxation to the electronic ground state via internal conversion process. Based on our calculated potential energy surfaces and experimental excited state lifetimes from the literature, we can infer that the 7H- must have a greater fluorescence yield than the 9H-tautomer. We also explored the triplet state population mechanisms on the 7H-26DAP to understand their long-lived components observed experimentally.

The tautomer-specific excited state dynamics of 2,6-diaminopurine using resonance-enhanced multiphoton ionization and quantum chemical calculations

2,6-Diaminopurine (2,6-dAP) is an alternative nucleobase that potentially played a role in prebiotic chemistry. We studied its excited state dynamics in the gas phase by REMPI, IR-UV hole burning, and ps pump-probe spectroscopy and performed quantum chemical calculations at the SCS-ADC(2) level of theory to interpret the experimental results. We found the 9H tautomer to have a small barrier to ultrafast relaxation via puckering of its 6-membered ring. The 7H tautomer has a larger barrier to reach a conical intersection and also has a sizable triplet yield. These results are discussed relative to other purines, for which 9H tautomerization appears to be more photostable than 7H and homosubstituted purines appear to be less photostable than heterosubstituted or singly substituted purines.

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

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.

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.

Exploring the Influence of Intermolecular Interactions in Prebiotic Chemistry Using Laser Spectroscopy and Calculations

One of the most fascinating questions in chemistry is why nature chose CGAT as the alphabet of life. Very likely, such selection was the result of multiple factors and a long period of refinement. Here, we explore how the intermolecular interactions influenced such process, by characterizing the formation of dimers between adenine, theobromine and 4-aminopyrimidine. Using a combination of mass-resolved excitation spectroscopy and DFT calculations, we determined the structure of adenine-theobromine and 4-aminopyrimidine-theobromine dimers. The binding energy of these dimers is very close to the canonical adenine-thymine nucleobases. Likewise, the dimers are able to adopt Watson-Crick conformations. These findings seem to indicate that there were many options available to build the first versions of the informational polymers, which also had to compete with other molecules, such as 4-aminopyrimidine, which does not have a valid attaching point for a saccharide. For some reason, nature did not select the most strongly-bonded partners or if it did, such proto-bases were later replaced by the nowadays canonical CGAT.

Photodynamics of alternative DNA base isoguanine

Pump–probe experiments and quantum-chemical simulations of UV-excited isoguanine elucidate its tautomer dependent photochemical properties.

How nature covers its bases

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

Excited State Dynamics of 6-Thioguanine

Here we present the excited state dynamics of jet-cooled 6-thioguanine (6-TG), using resonance-enhanced multiphoton ionization (REMPI), IR-UV double resonance spectroscopy, and pump-probe spectroscopy in the nanosecond and picosecond time domains. We report data on two thiol tautomers, which appear to have different excited state dynamics. These decay to a dark state, possibly a triplet state, with rates depending on tautomer form and on excitation wavelength, with the fastest rate on the order of 10¹⁰ s^-1. We also compare 6-TG with 9-enolguanine, for which we observed decay to a dark state with a 2 orders of magnitude smaller rate. At increased excitation energy (∼+500 cm^-1) an additional pathway appears for the predominant thiol tautomer. Moreover, the excited state dynamics for 6-TG thiols is different from that recently predicted for thiones.

Expanding a fluorescent RNA alphabet: synthesis, photophysics and utility of isothiazole-derived purine nucleoside surrogates

A series of emissive ribonucleoside purine mimics, all comprised of an isothiazolo[4,3-d]pyrimidine core, was prepared using a divergent pathway involving a key Thorpe-Ziegler cyclization. In addition to an adenosine and a guanosine mimic, analogues of the noncanonical xanthosine, isoguanosine, and 2-aminoadenosine were also synthesized and found to be emissive. Isothiazolo 2-aminoadenosine, an adenosine surrogate, was found to be particularly emissive and effectively deaminated by adenosine deaminase. Competitive studies with adenosine deaminase with each analogue in combination with native adenosine showed preference for the native substrate while still deaminating the isothiazolo analogues.

Competition between stacked and hydrogen bonded structures of cytosine aggregates

The four bases of DNA constitute what is known as the “alphabet of life”.

Investigation of the mechanisms of photo-induced formation of cyclobutane dimers of cytosine and 2,4-diaminopyrimidine

The mechanisms of the formation of cyclobutane dimers (CBD) of cytosine and 2,4-diaminopyrimidine were studied at the CC2 theoretical level and cc-pVDZ basis functions. Four orientations of the two monomers are explored: cys-syn, cis-anti, trans-syn, and trans-anti. The research revealed that in all cases the cyclobutane structures are formed along the (1)ππ* excited-state reaction paths of the stacked aggregates. We localized the S1/S0 conical intersections mediating those transformations. The results obtained agree well with the previously reported investigations on the cis-syn cyclodimer formations of other pyrimidines.

Intersystem Crossing Pathways in the Noncanonical Nucleobase 2-Thiouracil: A Time-Dependent Picture

The deactivation mechanism after ultraviolet irradiation of 2-thiouracil has been investigated using nonadiabatic dynamics simulations at the MS-CASPT2 level of theory. It is found that after excitation the S2 quickly relaxes to S1, and from there intersystem crossing takes place to both T2 and T1 with a time constant of 400 fs and a triplet yield above 80%, in very good agreement with recent femtosecond experiments in solution. Both indirect S1 → T2 → T1 and direct S1 → T1 pathways contribute to intersystem crossing, with the former being predominant. The results contribute to the understanding of how some noncanonical nucleobases respond to harmful ultraviolet light, which could be relevant for prospective photochemotherapeutic applications.

Protonation effect on the electronic properties of 2-pyridone monomer, dimer and its water clusters: a theoretical study

The CC2 (second order approximate coupled cluster method) has been applied to investigate protonation effect on electronic transition energies of 2-pyridone (2PY), 2-pyridone dimer, and micro-solvated 2-pyridone (0-2 water molecules). The PE profiles of protonated 2-pyridone (2PYH(+)) as well as monohydrated 2PYH(+) at the different electronic states have been investigated. The (1)πσ∗ state in protonated species (2PYH(+)) is a barrier free and dissociative state along the O-H stretching coordinate. In this reaction coordinate, the lowest lying (1)πσ∗ predissociates the bound S1((1)ππ∗) state, connecting the latter to a conical intersection with the S0 state. These conical intersections lead the (1)ππ∗ state to proceed as predissociative state and finally direct the excited system to the ground state. Furthermore, in presence of water molecule, the (1)πσ∗ state still remains dissociative but the conical intersection between (1)πσ∗ and ground state disappears. In addition, according to the CC2 calculation results, it has been predicted that protonation significantly blue shifts the S1-S0 electronic transition of monomer, dimer, and microhydrated 2-pyridone.

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.

On the Origin of the Canonical Nucleobases: An Assessment of Selection Pressures across Chemical and Early Biological Evolution

The native bases of RNA and DNA are prominent examples of the narrow selection of organic molecules upon which life is based. How did nature "decide" upon these specific heterocycles? Evidence suggests that many types of heterocycles could have been present on the early Earth. It is therefore likely that the contemporary composition of nucleobases is a result of multiple selection pressures that operated during early chemical and biological evolution. The persistence of the fittest heterocycles in the prebiotic environment towards, for example, hydrolytic and photochemical assaults, may have given some nucleobases a selective advantage for incorporation into the first informational polymers. The prebiotic formation of polymeric nucleic acids employing the native bases remains, however, a challenging problem to reconcile. Hypotheses have proposed that the emerging RNA world may have included many types of nucleobases. This is supported by the extensive utilization of non-canonical nucleobases in extant RNA and the resemblance of many of the modified bases to heterocycles generated in simulated prebiotic chemistry experiments. Selection pressures in the RNA world could have therefore narrowed the composition of the nucleic acid bases. Two such selection pressures may have been related to genetic fidelity and duplex stability. Considering these possible selection criteria, the native bases along with other related heterocycles seem to exhibit a certain level of fitness. We end by discussing the strength of the N-glycosidic bond as a potential fitness parameter in the early DNA world, which may have played a part in the refinement of the alphabetic bases.

Relaxation mechanisms of UV-photoexcited DNA and RNA nucleobases

A comprehensive effort in photodynamical ab initio simulations of the ultrafast deactivation pathways for all five nucleobases adenine, guanine, cytosine, thymine, and uracil is reported. These simulations are based on a complete nonadiabatic surface-hopping approach using extended multiconfigurational wave functions. Even though all five nucleobases share the basic internal conversion mechanisms, the calculations show a distinct grouping into purine and pyrimidine bases as concerns the complexity of the photodynamics. The purine bases adenine and guanine represent the most simple photodeactivation mechanism with the dynamics leading along a diabatic ππ* path directly and without barrier to the conical intersection seam with the ground state. In the case of the pyrimidine bases, the dynamics starts off in much flatter regions of the ππ* energy surface due to coupling of several states. This fact prohibits a clear formation of a single reaction path. Thus, the photodynamics of the pyrimidine bases is much richer and includes also nπ* states with varying importance, depending on the actual nucleobase considered. Trapping in local minima may occur and, therefore, the deactivation time to the ground state is also much longer in these cases. Implications of these findings are discussed (i) for identifying structural possibilities where singlet/triplet transitions can occur because of sufficient retention time during the singlet dynamics and (ii) concerning the flexibility of finding other deactivation pathways in substituted pyrimidines serving as candidates for alternative nucleobases.