Photophysics of Methyl-Substituted Uracils and Thymines and Their Water Complexes in the Gas Phase

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.

A Plausible Mechanism of Uracil Photohydration Involves an Unusual Intermediate

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

Time Resolved Photoelectron Spectroscopy as a Test of Electronic Structure and Nonadiabatic Dynamics

We compare different levels of theory for simulating excited state molecular dynamics and use time-resolved photoelectron spectroscopy measurements to benchmark the theory. We perform trajectory surface hopping simulations for uracil excited to the first bright state (ππ*) using three different levels of theory (CASSCF, MRCIS, and XMS-CASPT2) in order to understand the role of dynamical correlation in determining the excited state dynamics, with a focus on the coupling between different electronic states and internal conversion back to the ground state. These dynamics calculations are used to simulate the time-resolved photoelectron spectra. The comparison of the calculated and measured spectra allows us to draw conclusions regarding the relative insights and quantitative accuracy of the calculations at the three different levels of theory, demonstrating that detailed quantitative comparisons of time-resolved photoelectron spectra can be used to benchmark methodology.

Effect of dynamic correlation on the ultrafast relaxation of uracil in the gas phase

High level multi-reference non-adiabatic dynamics simulations reveal that uracil’s photoexcited S₂ state decays very quickly without any significant trapping.

Time-resolved dynamics in iodide-uracil-water clusters upon excitation of the nucleobase

The dynamics of iodide-uracil-water (I^-·U·H₂O) clusters following π-π^* excitation of the nucleobase are probed using time-resolved photoelectron spectroscopy. Photoexcitation of this cluster at 4.77 eV results in electron transfer from the iodide moiety to the uracil, creating a valence-bound anion within the cross correlation of the pump and probe laser pulses. This species can decay by a number of channels, including autodetachment and dissociation to I^- or larger anion fragments. Comparison of the energetics of the photoexcited cluster and its decay dynamics with those of the bare iodide-uracil (I^-·U) complex provides a sensitive probe of the effects of microhydration on these species.

Time-resolved radiation chemistry: femtosecond photoelectron spectroscopy of electron attachment and photodissociation dynamics in iodide-nucleobase clusters

Iodide-nucleobase (I-·N) clusters studied by time-resolved photoelectron spectroscopy (TRPES) are an opportune model system for examining radiative damage of DNA induced by low-energy electrons. By initiating charge transfer from iodide to the nucleobase and following the dynamics of the resulting transient negative ions (TNIs) with femtosecond time resolution, TRPES provides a novel window into the chemistry triggered by the attachment of low-energy electrons to nucleobases. In this Perspective, we examine and compare the dynamics of electron attachment, autodetachment, and photodissociation in a variety of I-·N clusters, including iodide-uracil (I-·U), iodide-thymine (I-·T), iodide-uracil-water (I-·U·H2O), and iodide-adenine (I-·A), to develop a more unified representation of our understanding of nucleobase TNIs. The experiments probe whether dipole-bound or valence-bound TNIs are formed initially and the subsequent time evolution of these species. We also provide an outlook for forthcoming applications of TRPES to larger iodide-containing complexes to enable the further investigation of microhydration dynamics in nucleobases, as well as electron attachment and photodissociation in more complex nucleic acid constituents.

Dynamics of electron attachment and photodissociation in iodide-uracil-water clusters via time-resolved photoelectron imaging

The dynamics of low energy electron attachment to monohydrated uracil are investigated using time-resolved photoelectron imaging to excite and probe iodide-uracil-water (I^-·U·H₂O) clusters. Upon photoexcitation of I^-·U·H₂O at 4.38 eV, near the measured cluster vertical detachment energy of 4.40 eV ± 0.05 eV, formation of both the dipole bound (DB) anion and valence bound (VB) anion of I^-·U·H₂O is observed and characterized using a probe photon energy of 1.58 eV. The measured binding energies for both anions are larger than those of the non-hydrated iodide-uracil (I^-·U) counterparts, indicating that the presence of water stabilizes the transient negative ions. The VB anion exhibits a somewhat delayed 400 fs rise when compared to I^-·U, suggesting that partial conversion of the DB anion to form the VB anion at early times is promoted by the water molecule. At a higher probe photon energy, 3.14 eV, I^- re-formation is measured to be the major photodissociation channel. This product exhibits a bi-exponential rise; it is likely that the fast component arises from DB anion decay by internal conversion to the anion ground state followed by dissociation to I^-, and the slow component arises from internal conversion of the VB anion.

Direct observation of vibrational energy dispersal via methyl torsions

Explicit evidence for the role of methyl rotor levels in promoting energy dispersal is reported. A set of coupled zero-order vibration/vibration-torsion (vibtor) levels in the S₁ state of para-fluorotoluene (pFT) are investigated. Two-dimensional laser-induced fluorescence (2D-LIF) and two-dimensional zero-kinetic-energy (2D-ZEKE) spectra are reported, and the assignment of the main features in both sets of spectra reveals that the methyl torsion is instrumental in providing a route for coupling between vibrational levels of different symmetry classes. We find that there is very localized, and selective, dissipation of energy via doorway states, and that, in addition to an increase in the density of states, a critical role of the methyl group is a relaxation of symmetry constraints compared to direct vibrational coupling.

Optical absorption and magnetic circular dichroism spectra of thiouracils: a quantum mechanical study in solution

The excited electronic states of thiouracils, the analogues of uracil where the carbonyl oxygens are substituted by sulphur atoms, have been investigated by computing the magnetic circular dichroism (MCD) and one-photon absorption (OPA) spectra at the TD-DFT level of theory.

Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Intramolecular vibrational redistribution in the non-radiative excited state decay of uracil in the gas phase: an ab initio molecular dynamics study

We report a study of intramolecular vibrational distribution (IVR) occurring in the electronic ground state of uracil (S0) in the gas phase, following photoexcitation in the lowest energy bright excited state (Sπ) and decay through the ethylene-like Sπ/S0 Conical Intersection (CI-0π). To this aim we have performed 20 independent ab initio molecular dynamics simulations starting from CI-0π (ten of them with 1 eV kinetic energy randomly distributed over the different molecular degrees of freedom) and 10 starting from the ground state minimum (Franck-Condon, FC, point), with the excess kinetic energy equal to the energy gap between CI-0π and the FC point. The simulations, exploiting PBE0/6-31G(d) calculations, were performed over an overall period of 10 ps. A thorough statistical analysis of the variation of the geometrical parameters of uracil during the simulation time and of the distribution of the kinetic energy among the different vibrational degrees of freedom provides a consistent picture of the IVR process. In the first 0-200 fs the structural dynamics involve mainly the recovery of the average planarity. In the 200-600 fs time range, a substantial activation of CO and NH degrees of freedom is observed. After 500-600 fs most of the geometrical parameters reach average values similar to those found after 10 ps, though the system cannot be considered to be in equilibrium yet.

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

The excited-state dynamics of the purine free base and 9-methylpurine are investigated using experimental and theoretical methods. Femtosecond broadband transient absorption experiments reveal that excitation of these purine derivatives in aqueous solution at 266 nm results primarily in ultrafast conversion of the S2(ππ*) state to the vibrationally excited (1)nπ* state. Following vibrational and conformational relaxation, the (1)nπ* state acts as a doorway state in the efficient population of the triplet manifold with an intersystem crossing lifetime of hundreds of picoseconds. Experiments show an almost 2-fold increase in the intersystem crossing rate on going from polar aprotic to nonpolar solvents, suggesting that a solvent-dependent energy barrier must be surmounted to access the singlet-to-triplet crossing region. Ab initio static and surface-hopping dynamics simulations lend strong support to the proposed relaxation mechanism. Collectively, the experimental and computational results demonstrate that the accessibility of the nπ* states and the topology of the potential energy surfaces in the vicinity of conical intersections are key elements in controlling the excited-state dynamics of the purine derivatives. From a structural perspective, it is shown that the purine chromophore is not responsible for the ultrafast internal conversion in the adenine and guanine monomers. Instead, C6 functionalization plays an important role in regulating the rates of radiative and nonradiative relaxation. C6 functionalization inhibits access to the (1)nπ* state while simultaneously facilitating access to the (1)ππ*(La)/S0 conical intersection, such that population of the (1)nπ* state cannot compete with the relaxation pathways to the ground state involving ring puckering at the C2 position.

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

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

Specific photodynamics in thymine clusters: the role of hydrogen bonding

A photoionization detected IR study of thymine and 1-methylthymine monohydrates and of their homodimers was carried out to shed some light on the structure of the thymine clusters whose complex photodynamics has recently been the subject of great interest. Under supersonic jet conditions, thymine forms doubly H-bonded cyclic clusters with water or another base preferentially via its N1-H group and the adjacent carbonyl group. This hydrate is of no biological relevance since the N1-H group is the sugar binding site in thymidine. On the other hand, 1-methylthymine forms the donor H-bonds only via the N3-H group. Hence, properties of the N1-H and the N3-H bound clusters of thymine can be studied using thymine and 1-methylthymine molecules, respectively. No biologically relevant conformations of the dimers and hydrates of thymine, contrary to those of 1-methylthymine, are observed under supersonic jet conditions. Thymine homodimer, which extensively fragments upon UV ionization by formation of a protonated monomer, exhibits two N1-H···O═C2 hydrogen bonds. The photodynamics of hydrated thymines is found to be extremely sensitive to the hydration site: ranging from an ultrafast relaxation in less than 100 fs up to formation of a dark state with the lifetime on the microsecond time scale.

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

To explore the excited state dynamics of pyrimidine derivatives, we performed a combined experimental and theoretical study. We present resonant two-photon ionization (R2PI) and IR-UV double resonance spectra of 2,4-diaminopyrimidine and 2,6-diaminopurine seeded in a supersonic jet by laser desorption. For 2,4-diaminopyrimidine (S(0)-->S(1) 34,459 cm(-1)), we observed only the diamino tautomer with an excited state lifetime bracketed between experimental limits of 10 ps and 1 ns. For 2,6-diaminopurine, we observed two tautomers, the 9H- (S(0)-->S(1) 34,881 cm(-1)) and 7H- (S(0)-->S(1) 32,215 cm(-1)) diamino forms, with excited state lifetimes of 6.3±0.4 ns and 8.7±0.8 ns, respectively. We investigated the nature of the excited state of 2,4-diaminopyrimidine by means of multi-reference ab initio methods. The calculations of stationary points in the ground and excited states, minima on the S(0)/S(1) crossing seam and connecting reaction paths show that several paths with negligible barriers exist, allowing ultrafast radiationless deactivation if excited at energies slightly higher than the band origin. The sub-nanosecond lifetime found experimentally is in good agreement with this finding.

An ab initio and AIM investigation into the hydration of 2-thioxanthine

BACKGROUND

Hydration is a universal phenomenon in nature. The interactions between biomolecules and water of hydration play a pivotal role in molecular biology. 2-Thioxanthine (2TX), a thio-modified nucleic acid base, is of significant interest as a DNA inhibitor yet its interactions with hydration water have not been investigated either computationally or experimentally. Here in, we reported an ab initio study of the hydration of 2TX, revealing water can form seven hydrated complexes.

RESULTS

Hydrogen-bond (H-bond) interactions in 1:1 complexes of 2TX with water are studied at the MP2/6-311G(d, p) and B3LYP/6-311G(d, p) levels. Seven 2TX...H2O hydrogen bonded complexes have been theoretically identified and reported for the first time. The proton affinities (PAs) of the O, S, and N atoms and deprotonantion enthalpies (DPEs) of different N-H bonds in 2TX are calculated, factors surrounding why the seven complexes have different hydrogen bond energies are discussed. The theoretical infrared and NMR spectra of hydrated 2TX complexes are reported to probe the characteristics of the proposed H-bonds. An improper blue-shifting H-bond with a shortened C-H bond was found in one case. NBO and AIM analysis were carried out to explain the formation of improper blue-shifting H-bonds, and the H-bonding characteristics are discussed.

CONCLUSION

2TX can interact with water by five different H-bonding regimes, N-H...O, O-H...N, O-H...O, O-H...S and C-H...O, all of which are medium strength hydrogen bonds. The most stable H-bond complex has a closed structure with two hydrogen bonds (N(7)-H...O and O-H...O), whereas the least stable one has an open structure with one H-bond. The interaction energies of the studied complexes are correlated to the PA and DPE involved in H-bond formation. After formation of H-bonds, the calculated IR and NMR spectra of the 2TX-water complexes change greatly, which serves to identify the hydration of 2TX.