Electronic and structural elements that regulate the excited-state dynamics in purine nucleobase derivatives

When Tautomers Matter: UV-Vis Absorption Spectra of Hypoxanthine in Aqueous Solution from Fully Atomistic Simulations

The UV-Vis spectrum of the solvated purine derivative Hypoxanthine (HYX) is investigated using the Quantum Mechanics/Fluctuating Charges (QM/FQ) multiscale approach combined with a sampling of configurations through atomistic Molecular Dynamics (MD) simulations. Keto 1H7H and 1H9H tautomeric forms of HYX are the most stable in aqueous solution and form different stable complexes with the surrounding water molecules, ultimately affecting the electronic absorption spectra. The final simulated spectrum resulting from the combination of the individual spectra of tautomers agrees very well with most of the characteristics in the measured spectrum. The importance of considering the effect of the solute tautomers and, in parallel, the contribution of the different solvent arrangements around the solute when modeling spectral properties, is highlighted. In addition, the high quality of the computed spectra leads to suggesting an alternative way for acquiring tautomeric populations from combined computational/experimental spectra.

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

Electronic relaxation mechanism of 9-methyl-2,6-diaminopurine and 2,6-diaminopurine-2'-deoxyribose in solution

Prolonged ultraviolet exposure results in the formation of cyclobutane pyrimidine dimers (CPDs) in RNA. Consequently, prebiotic photolesion repair mechanisms should have played an important role in the maintenance of the structural integrity of primitive nucleic acids. 2,6-Diaminopurine is a prebiotic nucleobase that repairs CPDs with high efficiency when incorporated into polymers. We investigate the electronic deactivation pathways of 2,6-diaminopurine-2'-deoxyribose and 9-methyl-2,6-diaminopurine in acetonitrile and aqueous solution to shed light on the photophysical and excited state properties of the 2,6-diaminopurine chromophore. Evidence is presented that both are photostable compounds exhibiting similar deactivation mechanisms upon the population of the S1 (ππ* La ) state at 290 nm. The mechanism involves deactivation through the C2- and C6-reaction coordinates and >99% of the excited state population decays through nonradiative pathways involving two conical intersections with the ground state. The radiative and nonradiative lifetimes are longer in aqueous solution compared to acetonitrile. While τ1 is similar in both derivatives, τ2 is ca. 1.5-fold longer in 2,6-diaminopurine-2'-deoxyribose due to a more efficient trapping in the S1 (ππ* La ) minimum. Therefore, 2,6-diaminopurine could have accumulated in significant quantities during prebiotic times to be incorporated into non-canonical RNA and play a significant role in its photoprotection.

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.

Dual Functionality of 6-Methylthioguanine: Synergistic Effects Enhancing the Photolability of DNA Nucleobases

A small chemical modification of the nucleobase structure can significantly enhance the photoactivity of DNA, which may incur DNA damage, thus holding promising applications in photochemotherapy treatment of cancers or pathogens. However, single substitution confers only limited phototoxicity to DNA. Herein, we combine femtosecond and nanosecond time-resolved spectroscopy with high-level ab initio calculations to disentangle the excited-state dynamics of 6-methylthioguanine (me6-TG) under variable wavelength UVA excitation (310-330 nm). We find that double substitution of nucleobases (thionation and methylation) boosts the photoactivity by introducing more reactive channels. Intriguingly, 1nNπ*, rather than 1nSπ*, acts as the doorway state engendering the formation of the long-lived reactive triplet state in me6-TG. The 1nNπ* induces a low spin-orbit coupling of 8.3 cm-1, which increases the intersystem crossing (ISC) time (2.91 ± 0.14 ns). Despite the slowed ISC, the triplet quantum yield (ΦT) still accounts for a large fraction (0.6 ± 0.1), consistent with the potential energy surface that favors excited-state bifurcation to 1nNπ*min (3.36 ± 0.15 ps) rather than 1ππ*min (5.05 ± 0.26 ps), such that the subsequent ISC to triplet via 1nNπ*min constitutes the main relaxation pathway in me6-TG. Although this ΦT is inferior to its single-substituted predecessor 6-thioguanine (6-TG, 0.8 ± 0.2), the effect of thionation in synergy with methylation opens a unique C-S bond cleavage pathway through crossing to a repulsive 1πσ* state, generating thiyl radicals as highly reactive intermediates that may invoke biological damage. This photodissociation channel is extremely difficult for conventional nucleobases. These findings demonstrate the synergistic effects of double functionality substitution in modulating excited-state dynamics and enhancing the photolabile character of DNA nucleobases, providing inspirations for the rational design of advanced photodynamic and photochemotherapy approaches.

Primary Photophysics of Nicotinamide Chromophores in Their Oxidized and Reduced Forms

Nicotinamide adenine dinucleotide (NADH) is an important enzyme cofactor with emissive properties that allow it to be used in fluorescence microscopies to study cell metabolism. Its oxidized form NAD+, on the other hand, is considered to produce negligible fluorescence. In this contribution, we describe the photophysics of the isolated nicotinamidic system in both its reduced and oxidized states. This was achieved through the study of model molecules that do not carry the adenine nucleotide since its absorbance would overlap with the absorption spectrum of the nicotinamidic chromophores. We studied three model molecules: nicotinamide (niacinamide, an oxidized form without nitrogen substitution), the oxidized chromophore 1-benzyl-3-carbamoyl-pyridinium bromide (NBzOx), and its reduced form 1-benzyl-1,4-dihydronicotinamide (NBz). For a full understanding of the dynamics, we performed both femtosecond-resolved emission and transient absorption experiments. The oxidized systems, nicotinamide and NBzOx, have similar photophysics, where the originally excited bright state decays on an ultrafast timescale of less than 400 fs. The depopulation of this state is followed by excited-state positive absorption signals, which evolve in two timescales: the first one is from 1 to a few picoseconds and is followed by a second decaying component of 480 ps for nicotinamide in water and of 80-90 ps for nicotinamide in methanol and NBzOx in aqueous solution. The long decay times are assigned as the S1 lifetimes populated from the original higher-lying bright singlet, where this state is nonemissive but can be detected by transient absorption. While for NBzOx in aqueous solution and for nicotinamide in methanol, the S1 signal decays to the solvent-only level, for the aqueous solutions of nicotinamide, a small transient absorption signal remains after the 480 ps decay. This residual signal was assigned to a small population of triplet states formed during the slower S1 decay for nicotinamide in water. The experimental results were complemented by XMS-CASPT2 calculations, which reveal that in the oxidized forms, the rapid evolution of the initial π-π* state is due to a direct crossing with lower-energy dark n-π* singlet states. This coincides with the experimental observation of long-lived nonemissive states (80 to 480 ps depending on the system). On the other hand, the reduced model compound NBz has a long-lived emissive π-π* S1 state, which decays with a 510 ps time constant, similarly to the parent compound NADH. This is consistent with the XMS-CASPT2 calculations, which show that for the reduced chromophore, the dark states lie at higher energies than the bright π-π* S1 state.

Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring

Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure-dynamics-optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.

Rationalizing the Unexpected Sensitivity in Excited State Lifetimes of Adenine to Tautomerization by Nonadiabatic Molecular Dynamics

The remarkable photostability of canonical nucleobases makes them ideal building blocks for DNA and RNA. Even minor structural changes are expected to lead to drastic alteration of their subpicosecond excited state lifetimes. However, it is interesting to note that while the 9H- and 7H-amino tautomers of adenine possess drastically different lifetimes, 9H- and 7H-keto guanine possess similar excited state lifetimes. With an aim to explain this unexpected difference in sensitivity of lifetimes to tautomerization, we have investigated the excited state relaxation mechanism of UV-excited adenine and guanine tautomers using surface hopping based nonadiabatic molecular dynamics. We find that internal conversion in both guanine tautomers is almost barrierless while both adenine tautomers encounter significant barriers before they can deactivate. Moreover, the major deactivation channel (C2-puckering) in 9H-amino adenine is overall more efficient than the one (C6-puckering) in the 7H-amino form. We trace this difference to the frequent rotation of the amino group which disrupts its conjugation with the heterocyclic ring thereby reducing the strength of nonadiabatic coupling and, hence, delaying internal conversion.

Excited state dynamics of 2'-deoxyisoguanosine and isoguanosine in aqueous solution

Photostability is thought to be an inherent property of nucleobases required to survive the extreme ultraviolet radiation conditions of the prebiotic era. Previous studies have shown that absorption of ultraviolet radiation by the canonical nucleosides results in ultrafast internal conversion to the ground state, demonstrating that these nucleosides efficiently dissipate the excess electronic energy to the environment. In recent years, studies on the photophysical and photochemical properties of nucleobase derivatives have revealed that chemical substitution influences the electronic relaxation pathways of purine and pyrimidine nucleobases. It has been suggested that amino or carbonyl substitution at the C6 position could increase the photostability of the purine derivatives more than the substitution at the C2 position. This investigation aims to elucidate the excited state dynamics of 2'-deoxyisoguanosine (dIsoGuo) and isoguanosine (IsoGuo) in aqueous solution at pH 7.4 and 1.4, which contain an amino group at the C6 position and a carbonyl group at the C2 position of the purine chromophore. The study of these derivatives is performed using absorption and emission spectroscopies, broadband transient absorption spectroscopy, and density functional and time-dependent density functional levels of theory. It is shown that the primary relaxation mechanism of dIsoGuo and IsoGuo involves nonradiative decay pathways, where the population decays from the S₁(ππ*) state through internal conversion to the ground state via two relaxation pathways with lifetimes of hundreds of femtoseconds and less than 2 ps, making these purine nucleosides photostable in aqueous solution.

Disclosing the Role of C4-Oxo Substitution in the Photochemistry of DNA and RNA Pyrimidine Monomers: Formation of Photoproducts from the Vibrationally Excited Ground State

Oxo and amino substituted purines and pyrimidines have been suggested as protonucleobases participating in ancient pre-RNA forms. Considering electromagnetic radiation as a key environmental selection pressure on early Earth, the investigation of the photophysics of modified nucleobases is crucial to determine their viability as nucleobases' ancestors and to understand the factors that rule the photostability of natural nucleobases. In this Letter, we combine femtosecond transient absorption spectroscopy and quantum mechanical simulations to reveal the photochemistry of 4-pyrimidinone, a close relative of uracil. Irradiation of 4-pyrimidinone with ultraviolet radiation populates the S₁(ππ*) state, which decays to the vibrationally excited ground state in a few hundred femtoseconds. Analysis of the postirradiated sample in water reveals the formation of a 6-hydroxy-5H-photohydrate and 3-(N-(iminomethyl)imino)propanoic acid as the primary photoproducts. 3-(N-(Iminomethyl)imino)propanoic acid originates from the hydrolysis of an unstable ketene species generated from the C4-N3 photofragmentation of the pyrimidine core.

Photostability of 2,6-diaminopurine and its 2'-deoxyriboside investigated by femtosecond transient absorption spectroscopy

Ultraviolet radiation (UVR) from the sun is essential for the prebiotic syntheses of nucleotides, but it can also induce photolesions such as the cyclobutane pyrimidine dimers (CPDs) to RNA or DNA oligonucleotide in prebiotic Earth. 2,6-Diaminopurine (26DAP) has been proposed to repair CPDs in high yield under prebiotic conditions and be a key component in enhancing the photostability of higher-order prebiotic DNA structures. However, its electronic relaxation pathways have not been studied, which is necessary to know whether 26DAP could have survived the intense UV fluxes of the prebiotic Earth. We investigate the electronic relaxation mechanism of both 26DAP and its 2'-deoxyribonucleoside (26DAP-d) in aqueous solution using steady-state and femtosecond transient absorption measurements that are complemented with electronic-structure calculations. The results demonstrate that both purine derivatives are significantly photostable to UVR. It is shown that upon excitation at 287 nm, the lowest energy ¹ππ* state is initially populated. The population then branches following two relaxation coordinates in the ¹ππ* potential energy surface, which are identified as the C2- and C6-relaxation coordinates. The population following the C6-coordinate internally converts to the ground state nonradiatively through a nearly barrierless conical intersection within 0.7 ps in 26DAP or within 1.1 ps in 26DAP-d. The population that follows the C2-relaxation coordinate decays back to the ground state by a combination of nonradiative internal conversion via a conical intersection and fluorescence emission from the ¹ππ* minimum in 43 ps and 1.8 ns for the N9 and N7 tautomers of 26DAP, respectively, or in 70 ps for 26DAP-d. Fluorescence quantum yields of 0.037 and 0.008 are determined for 26DAP and 26DAP-d, respectively. Collectively, it is demonstrated that most of the excited state population in 26DAP and 26DAP-d decays back to the ground state via both nonradiative and radiative relaxation pathways. This result lends support to the idea that 26DAP could have accumulated in large enough quantities during the prebiotic era to participate in the formation of prebiotic RNA or DNA oligomers and act as a key component in the protection of the prebiotic genetic alphabet.

Photophysical and Electrical Properties of Highly Luminescent 2/6-Triazolyl-Substituted Push-Pull Purines

New push-pull N(9)-alkylated 6-piperidino-2-triazolylpurine and 2-piperidino-6-triazolylpurine derivatives are synthesized, and their optical and optoelectronic properties are comprehensively characterized with experimental and computational methods. The compounds possess intense violet or blue fluorescence with fluorescence quantum yields of up to 91% in solution and 40% in host-free films. Depending on their structural composition, the compounds have ionization energy in the range of 5.25-6.04 eV, electron affinity of 2.18-3.15 eV, and triplet energy of 2.52-2.95 eV. Due to the presence of hole-transporting purine and electron-transporting triazole fragments, compounds exhibit bipolar charge-transportation ability. Despite the favorable emissive properties of the studied push-pull purines, their electroluminescence in thin films is quenched owing to large current densities that are present even at a moderate driving voltage. This marks application directions related to a predominantly charge-transportation functionality as the most suitable for this compound class.

A molecular insight into the photophysics of barbituric acid, a candidate for canonical nucleobases' ancestor

This work investigates the photophysics of barbituric acid at different pH conditions using ab initio methods. Our calculations ascribe the most intense bands at ca. 260 nm at neutral pH and 210 nm at acidic pH conditions in the absorption spectra of this chromophore to the lowest lying ππ* transitions. Consistently with the ultrashort excited state lifetimes experimentally registered, the potential energy landscapes of both the neutral and deprotonated forms of barbituric acid combined with the interpretation of their transient absorption spectra suggest the deactivation of these systems along the singlet manifold. Compared to uracil, its closest natural nucleobase, barbituric acid presents a red shifted absorption spectrum, due to the lowering by more than 0.5 eV of the lowest-energy ππ* excited state, and a much more complex topography of the S₁ potential energy surface, with several energetically accessible local minima. This fact, however, does not affect the excited state lifetimes, which for barbituric acid were experimentally registered in the sub-ps time scale.

Implementation of the interacting quantum atom energy decomposition using the CASPT2 method

We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.

Excited state dynamics of 7-deazaguanosine and guanosine 5'-monophosphate

Minor structural modifications to the DNA and RNA nucleobases have a significant effect on their excited state dynamics and electronic relaxation pathways. In this study, the excited state dynamics of 7-deazaguanosine and guanosine 5'-monophosphate are investigated in aqueous solution and in a mixture of methanol and water using femtosecond broadband transient absorption spectroscopy following excitation at 267 nm. The transient spectra are collected using photon densities that ensure no parasitic multiphoton-induced signal from solvated electrons. The data can be fit satisfactorily using a two- or three-component kinetic model. By analyzing the results from steady-state, time-resolved, computational calculations, and the methanol-water mixture, the following general relaxation mechanism is proposed for both molecules, L_b → L_a → ¹πσ^*(ICT) → S₀, where the ¹πσ^*(ICT) stands for an intramolecular charge transfer excited singlet state with significant πσ^* character. In general, longer lifetimes for internal conversion are obtained for 7-deazaguanosine compared to guanosine 5'-monophosphate. Internal conversion of the ¹πσ^*(ICT) state to the ground state occurs on a similar time scale of a few picoseconds in both molecules. Collectively, the results demonstrate that substitution of a single nitrogen atom for a methine (C-H) group at position seven of the guanine moiety stabilizes the ¹ππ^* L_b and L_a states and alters the topology of their potential energy surfaces in such a way that the relaxation dynamics in 7-deazaguanosine are slowed down compared to those in guanosine 5'-monophosphate but not for the internal conversion of ¹πσ^*(ICT) state to the ground state.

Fundamental photophysics of isomorphic and expanded fluorescent nucleoside analogues

Fluorescent nucleoside analogues (FNAs) are structurally diverse mimics of the natural essentially non-fluorescent nucleosides which have found numerous applications in probing the structure and dynamics of nucleic acids as well as their interactions with various biomolecules. In order to minimize disturbance in the labelled nucleic acid sequences, the FNA chromophoric groups should resemble the natural nucleobases in size and hydrogen-bonding patterns. Isomorphic and expanded FNAs are the two groups that best meet the criteria of non-perturbing fluorescent labels for DNA and RNA. Significant progress has been made over the past decades in understanding the fundamental photophysics that governs the spectroscopic and environmentally sensitive properties of these FNAs. Herein, we review recent advances in the spectroscopic and computational studies of selected isomorphic and expanded FNAs. We also show how this information can be used as a rational basis to design new FNAs, select appropriate sequences for optimal spectroscopic response and interpret fluorescence data in FNA applications.

Modeling Spin-Crossover Dynamics

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.

Theoretical insights on the excited-state-deactivation mechanisms of protonated thymine and cytosine

Ab initio and surface-hopping nonadiabatic dynamics simulation methods were employed to investigate relaxation mechanisms in protonated thymine (TH⁺) and cytosine (CH⁺). A few conical intersections were located between ¹ππ* and S₀ states for each system with the CASSCF (8,8) theoretical model and relevant contributions to the deactivation mechanism of titled systems were addressed by the determination of potential energy profiles at the CASPT2 (12,10) theoretical level. It was revealed that the relaxation of the ¹ππ* state of the most stable conformer of both systems to the ground state is mostly governed by the accessible S₁/S₀ conical intersection resulting from the barrier-free out-of-plane deformation. Interestingly, it was exhibited that the ring puckering coordinate driven from the C₆ position of the heterocycle ring in TH⁺ and CH⁺ plays the most prominent role in the deactivation mechanism of considered systems. Our ab initio results are also supported by excited-state nonadiabatic dynamics simulations based on ADC(2), describing the ultrashort S₁ lifetime of TH⁺/CH⁺ by analyzing trajectories leading excited systems to the ground. It was confirmed that the excited-state population mostly relaxes to the ground via the ring puckering coordinate from the C₆ moiety. Overall, the theoretical results of this study shed light on the deactivation mechanism of protonated DNA bases.

New Insights about the Photostability of DNA/RNA Bases: Triplet nπ* State Leads to Effective Intersystem Crossing in Pyrimidinones

The high photostability of DNA/RNA nucleobases is attributed to the effective internal conversions of their bright ¹ππ* states to the ground state through conical intersections. Intersystem crossing (ISC) from singlet to triplet excited states is a minor decay pathway in nucleobases and it is observed with ∼1-2% quantum yields (QYs) in pyrimidine bases. Presumably, ISC in pyrimidines takes place from the dark singlet ¹nπ* state to the lowest triplet ³ππ* state. However, recent studies showed that ISC from the initial populated bright ¹ππ* state to higher energy triplet ³nπ* states indeed occurs in the subpicosecond timescale. Such a mechanism is still poorly understood since direct observation of this pathway is challenging. Herein, excited state dynamics of three pyrimidinones, which share the same skeleton with pyrimidine bases, is investigated in different solvents. Compared to canonical pyrimidine bases, removing the oxygen atom at the C4 position revokes the low-lying dark ¹nπ* state in pyrimidinones, resulting in direct ISC from the S₁ (¹ππ*) state to triplet T₃ (³nπ*) state with much higher QYs. Meanwhile, hydrogen bonding between the carbonyl group in pyrimidinones and protic solvents can accelerate vibrational cooling of the hot S₁ (¹ππ*) state, leading to higher fluorescence QYs and smaller ISC rate constants. These results not only evidence the hypothesis of the direct ¹ππ* → ³nπ* ISC mechanism, but also contribute to a better understanding of triplet formation in pyrimidines.

Photoexcitation dynamics in perylene diimide dimers

We utilize first-principles theory to investigate photo-induced excited-state dynamics of functionalized perylene diimide. This class of materials is highly suitable for solar energy conversion because of the strong optical absorbance, efficient energy transfer, and chemical tunability. We couple time-dependent density functional theory to a recently developed time-resolved non-adiabatic dynamics approach based on a semi-empirical description. By studying the monomer and dimer, we focus on the role stacking plays on the time-scales associated with excited-state non-radiative relaxation from a high excitonic state to the lowest energy exciton. We predict that the time-scale for energy conversion in the dimer is significantly faster than that in the monomer when equivalent excited states are accounted for. Additionally, for the dimer, the decay from the second to the nearly degenerate lowest energy excited-state involves two time-scales: a rapid decay on the order of ∼10 fs followed by a slower decay of ∼100 fs. Analysis of the spatial localization of the electronic transition density during the internal conversion process points out the existence of localized states on individual monomers, indicating that the strength of thermal fluctuations exceeds electronic couplings between the states such that the exciton hops between localized states throughout the simulation.

On the Origin of the Photostability of DNA and RNA Monomers: Excited State Relaxation Mechanism of the Pyrimidine Chromophore

Today's genetic composition is the result of continual refinement processes on primordial heterocycles present in prebiotic Earth and at least partially regulated by ultraviolet radiation. Femtosecond transient absorption spectroscopy and state-of-the-art ab initio calculations are combined to unravel the electronic relaxation mechanism of pyrimidine, the common chromophore of the nucleobases. The excitation of pyrimidine at 268 nm populates the S₁(nπ*) state directly. A fraction of the population intersystem crosses to the triplet manifold within 7.8 ps, partially decaying within 1.5 ns, while another fraction recovers the ground state in >3 ns. The pyrimidine chromophore is not responsible for the photostability of the nucleobases. Instead, C2 and C4 amino and/or carbonyl functionalization is essential for shaping the topography of pyrimidine's potential energy surfaces and results in accessible conical intersections between the initially populated electronic excited state and the ground state.

Radiation Induced One-Electron Oxidation of 2-Thiouracil in Aqueous Solutions

Oxidative damage to 2-thiouracil (2-TU) by hydroxyl (•OH) and azide (●N3) radicals produces various primary reactive intermediates. Their optical absorption spectra and kinetic characteristics were studied by pulse radiolysis with UV-vis spectrophotometric and conductivity detection and by time-dependent density functional theory (TD-DFT) method. The transient absorption spectra recorded in the reactions of •OH with 2-TU depend on the concentration of 2-TU, however, only slightly on pH. At low concentrations, they are characterized by a broad absorption band with a weakly pronounced maxima located at λ = 325, 340 and 385 nm, whereas for high concentrations, they are dominated by an absorption band with λmax ≈ 425 nm. Based on calculations using TD-DFT method, the transient absorption spectra at low concentration of 2-TU were assigned to the ●OH-adducts to the double bond at C5 and C6 carbon atoms (3●, 4●) and 2c-3e bonded ●OH adduct to sulfur atom (1…●OH) and at high concentration of 2-TU also to the dimeric 2c-3e S-S-bonded radical in neutral form (2●). The dimeric radical (2●) is formed in the reaction of thiyl-type radical (6●) with 2-TU and both radicals are in an equilibrium with Keq = 4.2 × 103 M-1. Similar equilibrium (with Keq = 4.3 × 103 M-1) was found for pH above the pKa of 2-TU which involves admittedly the same radical (6●) but with the dimeric 2c-3e S-S bonded radical in anionic form (2●-). In turn, ●N3-induced oxidation of 2-TU occurs via radical cation with maximum spin location on the sulfur atom which subsequently undergoes deprotonation at N1 atom leading again to thiyl-type radical (6●). This radical is a direct precursor of dimeric radical (2●).

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

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.

pH Controlled Intersystem Crossing and Singlet Oxygen Generation of 8-Azaadenine in Aqueous Solution

Azabases are intriguing DNA and RNA analogues and have been used as effective antiviral and anticancer medicines. However, photosensitivity of these drugs has also been reported. Here, pH-controlled intersystem crossing (ISC) process of 9H 8-azaadenine (8-AA) in aqueous solution is reported. Broadband transient absorption measurements reveal that the hydrogen atom at N9 position can greatly affect ISC of 8-AA and ISC is more favorable when 8-AA is in its neutral form in aqueous solution. The initial excited ππ* (S₂ ) state evolves through ultrafast internal conversion (IC) (4.2 ps) to the lower-lying nπ* state (S₁ ), which further stands as a door way state for ISC with a time constant of 160 ps. The triplet state has a lifetime of 6.1 μs. On the other hand, deprotonation at N9 position promotes the IC from the ππ* (S₂ ) state to the ground state (S₀ ) and the lifetime of the S₂ state is determined to be 10 ps. The experimental results are further supported by time-dependent density functional theory (TDDFT) calculations. Singlet oxygen generation yield is measured to be 13.8 % for the neutral 8-AA while the deprotonated one exhibit much lower yield (<2 %), implying that this compound could be a potential pH-sensitized photodynamic therapy agent.

Tracking the origin of photostability in purine nucleobases: the photophysics of 2-oxopurine

Molding purine PES through functionalization: whilst purine C₂-substitution maintains the features of the spectroscopic PES of the heterocycle, C₆-functionalization reshapes its topography leading to photostable systems.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Photochemical and Photodynamical Properties of Sulfur-Substituted Nucleic Acid Bases

Sulfur-substituted nucleobases (a.k.a., thiobases) are among the world's leading prescriptions for chemotherapy and immunosuppression. Long-term treatment with azathioprine, 6-mercaptopurine and 6-thioguanine has been correlated with the photoinduced formation of carcinomas. Establishing an in-depth understanding of the photochemical properties of these prodrugs may provide a route to overcoming these carcinogenic side effects, or, alternatively, a basis for developing effective compounds for targeted phototherapy. In this review, a broad examination is undertaken, surveying the basic photochemical properties and excited-state dynamics of sulfur-substituted analogs of the canonical DNA and RNA nucleobases. A molecular-level understanding of how sulfur substitution so remarkably perturbs the photochemical properties of the nucleobases is presented by combining experimental results with quantum-chemical calculations. Structure-property relationships demonstrate the impact of site-specific sulfur substitution on the photochemical properties, particularly on the population of the reactive triplet state. The value of fundamental photochemical investigations for driving the development of ultraviolet-A chemotherapeutics is showcased. The most promising photodynamic agents identified thus far have been investigated in various carcinoma cell lines and shown to decrease cell proliferation upon exposure to ultraviolet-A radiation. Overarching principles have been elucidated for the impact that sulfur substitution of the carbonyl oxygen has on the photochemical properties of the nucleobases.

Photochemical relaxation pathways of S6-methylthioinosine and O6-methylguanosine in solution

S6-Methylthioinosine and O6-methylguanosine are byproducts resulting from the enzymatic reactions of sulfur-substituted prodrugs in cells and from the interaction of alkylating agents with cellular DNA, respectively. Their photochemistry has not been investigated, and it is currently unknown whether light absorption by these byproducts may pose any threat to the cell. In this contribution, their photoinduced processes upon absorption of UVB radiation are reported using broadband transient absorption spectroscopy. Plausible electronic relaxation mechanisms are proposed for both biological molecules, which are supported by steady-state absorption and emission measurements, and by singlet and triplet vertical excitation energies performed on a large subset of ground-state optimized conformational isomers in solution. The results are compared to the body of knowledge gathered in the scientific literature about the light-induced processes in the sulfur-substituted and canonical purine monomers. In particular, it is shown that S6-methylation decreases the rate to populate the lowest-energy triplet state and blueshifts the ground-state absorption spectrum compared to those for the sulfur-substituted prodrugs and for the 6-thioguanosine metabolite. Similarly, O6-methylation decreases the rate of internal conversion to the ground state observed in the guanine monomers by more than 10-fold in acetonitrile and 40-fold in aqueous solution, while it redshifts the ground-state absorption spectrum. Collectively, this investigation provides relevant new insights about the relationship between structural modifications of the purine chromophore and the electronic relaxation mechanisms in this important group of biological molecules.

Unusually Short-Lived Solvent-Dependent Excited State in a Half-Sandwich Ru(II) Complex Induced by Low-Lying 3MC States

A ruthenium complex with a half-sandwich geometry ([(p-cymene)Ru(Cl)(curcuminoid)]) was synthesized, characterized, and investigated regarding its ultrafast photophysics. These photophysical investigations of the complex revealed a weak and short-lived emission from the initially populated ¹MLCT state and solvent-dependent photoinduced dynamics, where the secondarily populated ³MC state is stabilized by nonpolar solvents. Overall the decay of the ³dd-MC state to the ground state is completed within picoseconds. This short excited-state lifetime is in stark contrast to the typically observed long-lived ³MLCT states with lifetimes of nanoseconds or microseconds in unstrained, octahedral ruthenium complexes but is in good agreement with the findings for distorted octahedral complexes. This is pointing to the half-sandwich geometry as a new and easy approach to study these otherwise often concealed dd states.

Efficient intersystem crossing in 2-aminopurine riboside probed by femtosecond time-resolved transient vibrational absorption spectroscopy

The photophysical dynamics of 2-aminopurine, a fluorescent analogue of the canonical nucleobase adenine, has been studied by femtosecond transient vibrational absorption spectroscopy.

An ab initio multiple cloning approach for the simulation of photoinduced dynamics in conjugated molecules

Multidimensional wave function: a superposition of Gaussian coherent states guided by Ehrenfest trajectories suited to clone and swap their electronic amplitudes.

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.

Excited-State Dynamics in O6-Methylguanosine: Impact of O6-Methylation on the Relaxation Mechanism of Guanine Monomers

Absorption of ultraviolet radiation by DNA bases results in ultrafast internal conversion to the ground state, which minimizes photodamage. However, exogenous and endogenous alkylating agents present in the cellular environment can methylate the nucleobases in DNA. In particular, methylation of guanosine at the O⁶ position in DNA leads to the formation of the O⁶-methylguanosine adduct, which may alter the photostability of DNA. This contribution demonstrates that O⁶-methylation of guanosine red shifts its ground-state absorption spectrum and slows down the rate of internal conversion to the ground state by ∼40-fold in aqueous solution. The 40-fold decrease in the rate of excited-state decay increases the probability of photodamage within cellular DNA. It is proposed that the longer decay lifetime corresponds to relaxation of the excited-state population in O⁶-methylguanosine along a C6-puckered reaction coordinate in the ¹ππ*(L_a) potential energy surface that runs parallel to an ultrafast internal conversion pathway along a C2-puckered coordinate.

Trajectory Surface-Hopping Dynamics Including Intersystem Crossing in [Ru(bpy)3]2

Surface-hopping dynamics coupled to linear response TDDFT and explicit nonadiabatic and spin-orbit couplings have been used to model the ultrafast intersystem crossing (ISC) dynamics in [Ru(bpy)₃]²⁺. Simulations using an ensemble of trajectories starting from the singlet metal-to-ligand charge transfer (¹MLCT) band show that the manifold of ³MLCT triplet states is first populated from high-lying singlet states within 26 ± 3 fs. ISC competes with an intricate internal conversion relaxation process within the singlet manifold to the lowest singlet state. Normal-mode analysis and principal component analysis, combined with further dynamical simulations where the nuclei are frozen, unequivocally demonstrate that it is not only the high density of states and the large spin-orbit couplings of the system that promote ISC. Instead, geometrical relaxation involving the nitrogen atoms is required to allow for state mixing and efficient triplet population transfer.

Sub-50 fs excited state dynamics of 6-chloroguanine upon deep ultraviolet excitation

The photophysical properties of natural nucleobases and their respective nucleotides are ascribed to the sub-picosecond lifetime of their first singlet states in the UV-B region (260-350 nm). Electronic transitions of the ππ* type, which are stronger than those in the UV-B region, lie at the red edge of the UV-C range (100-260 nm) in all isolated nucleobases. The lowest energetic excited states in the UV-B region of nucleobases have been investigated using a plethora of experimental and theoretical methods in gas and solution phases. The sub-picosecond lifetime of these molecules is not a general attribute of all nucleobases but specific to the five primary nucleobases and a few xanthine and methylated derivatives. To determine the overall UV photostability, we aim to understand the effect of more energetic photons lying in the UV-C region on nucleobases. To determine the UV-C initiated photophysics of a nucleobase system, we chose a halogen substituted purine, 6-chloroguanine (6-ClG), that we had investigated previously using resonance Raman spectroscopy. We have performed quantitative measurements of the resonance Raman cross-section across the Bb absorption band (210-230 nm) and constructed the Raman excitation profiles. We modeled the excitation profiles using Lee and Heller's time-dependent theory of resonance Raman intensities to extract the initial excited state dynamics of 6-ClG within 30-50 fs after photoexcitation. We found that imidazole and pyrimidine rings of 6-ClG undergo expansion and contraction, respectively, following photoexcitation to the Bb state. The amount of distortions of the excited state structure from that of the ground state structure is reflected by the total internal reorganization energy that is determined at 112 cm(-1). The contribution of the inertial component of the solvent response towards the total reorganization energy was obtained at 1220 cm(-1). In addition, our simulation also yields an instantaneous response of the first solvation shell within an ultrafast timescale of less than 30 fs following photoexcitation.

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.

Distortion dependent intersystem crossing: A femtosecond time-resolved photoelectron spectroscopy study of benzene, toluene, and p-xylene

The competition between ultrafast intersystem crossing and internal conversion in benzene, toluene, and p-xylene is investigated with time-resolved photoelectron spectroscopy and quantum chemical calculations. By exciting to S2 out-of-plane symmetry breaking, distortions are activated at early times whereupon spin-forbidden intersystem crossing becomes (partly) allowed. Natural bond orbital analysis suggests that the pinnacle carbon atoms distorting from the aromatic plane change hybridization between the planar Franck-Condon geometry and the deformed (boat-shaped) S2 equilibrium geometry. The effect is observed to increase in the presence of methyl-groups on the pinnacle carbon-atoms, where largest extents of σ and π orbital-mixing are observed. This is fully consistent with the time-resolved spectroscopy data: Toluene and p-xylene show evidence for ultrafast triplet formation competing with internal conversion, while benzene appears to only decay via internal conversion within the singlet manifold. For toluene and p-xylene, internal conversion to S1 and intersystem crossing to T3 occur within the time-resolution of our instrument. The receiver triplet state (T3) is found to undergo internal conversion in the triplet manifold within ≈100-150 fs (toluene) or ≈180-200 fs (p-xylene) as demonstrated by matching rise and decay components of upper and lower triplet states. Overall, the effect of methylation is found to both increase the intersystem crossing probability and direct the molecular axis of the excited state dynamics.

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.

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.

Ultrafast Electronic Deactivation Dynamics of Xanthosine Monophosphate

Ultrafast energy dissipation is a crucial factor for the photostability of DNA and RNA, but even some of the key electronic deactivation pathways in monomeric nucleic acid building stones are still controversial. Here, we report on the excited-state dynamics of the rare nucleotide xanthosine monophosphate as a function of deprotonation state (XMP vs. XMP−) and excitation wavelength (λpump= 278–243 nm) by femtosecond time-resolved fluorescence and absorption spectroscopy. We show that the predominating relaxation channel leads to a return of the photo-excited molecules to the electronic ground state in τ∼1 ps. The mechanism likely involves an out-of-plane deformation of the five-membered ring, different from the main electronic deactivation pathways in the canonical purine bases adenine and guanine. The results are discussed in terms of the structural and electronic differences of XMP compared to the canonical nucleotides.

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.

Ultrafast structural dynamics of photoexcited adenine

Ultraviolet Resonance Raman (UVRR) spectroscopy derives distinct electronic properties of adenine in the L_a (260 nm) and B_b (210 nm) excited states.

Electronic relaxation pathways of the biologically relevant pterin chromophore

Femtosecond-to-microsecond transient absorption spectroscopy is used to report the ultrafast relaxation mechanism of 2-amino-1H-pteridin-4-one (pterin) for the first time.