Excited state dynamics of DNA bases

Theoretical insight into photodeactivation mechanisms of adenine-uracil and adenine-thymine nucleobase pairs

In this work, several plausible intra- and intermolecular photoinduced processes of the Watson-Crick base pairs of adenine with uracil (A-U) or thymine (A-T) according to the results of spin component scaling variant of algebraic diagrammatic construction up to the second order [SCS-ADC(2)] calculations are discussed. Although widely explored, these systems lack complete characterization of possible intramolecular relaxation channels perturbed by intermolecular interactions. In particular, we address the still open debate on photodeactivation via purine-ring puckering at the C2 or C6-atom position of adenine. We also show that the presence of low-lying, long-lived 1nπ* states can be a significant factor in hindering relaxation via an electron-driven proton transfer process, as the population of these states can lead to an efficient intersystem crossing to a triplet manifold, the estimated rate of which is 1.6 × 1010 s-1 which exceeds the corresponding internal conversion to the ground state by an order of magnitude. Additionally, the SCS variant of the ADC(2) method is shown to provide a more balanced description of valence and charge-transfer excited states.

UV wavelength-dependent photoionization quantum yields for the dark 1nπ* state of aqueous thymidine

Despite the important role of the dark 1nπ* state in the photostability of thymidine in aqueous solution, no detailed ultraviolet (UV) wavelength-dependent investigation of the 1nπ* quantum yield (QY) in aqueous thymidine has been experimentally performed. Here, we investigate the wavelength-dependent photoemission spectra of aqueous thymidine from 266.7 to 240 nm using liquid-microjet photoelectron spectroscopy. Two observed ionization channels are assigned to resonant ionizations from 1ππ* to the cationic ground state D0 (π-1) and 1nπ* to the cationic excited state D1 (n-1). The weak 1nπ* → D1 ionization channel appears due to ultrafast 1ππ* → 1nπ* internal conversion within the pulse duration of ∼180 fs. The obtained 1nπ* quantum yields exhibit a strong wavelength dependence, ranging from 0 to 0.27 ± 0.01, suggesting a hitherto uncharacterized 1nπ* feature. The corresponding vertical ionization energies (VIEs) of D0 and D1 of aqueous thymidine are experimentally determined to be 8.47 ± 0.12 eV and 9.22 ± 0.29 eV, respectively. Our UV wavelength-dependent QYs might indicate that different structural critical points to connect the multidimensional 1ππ*/1nπ* conical intersection seam onto the multidimensional potential energy surface of the 1ππ* state might exist and determine the relaxation processes of aqueous thymidine upon UV excitation.

The Pigment World: Life's Origins as Photon-Dissipating Pigments

Many of the fundamental molecules of life share extraordinary pigment-like optical properties in the long-wavelength UV-C spectral region. These include strong photon absorption and rapid (sub-pico-second) dissipation of the induced electronic excitation energy into heat through peaked conical intersections. These properties have been attributed to a "natural selection" of molecules resistant to the dangerous UV-C light incident on Earth's surface during the Archean. In contrast, the "thermodynamic dissipation theory for the origin of life" argues that, far from being detrimental, UV-C light was, in fact, the thermodynamic potential driving the dissipative structuring of life at its origin. The optical properties were thus the thermodynamic "design goals" of microscopic dissipative structuring of organic UV-C pigments, today known as the "fundamental molecules of life", from common precursors under this light. This "UV-C Pigment World" evolved towards greater solar photon dissipation through more complex dissipative structuring pathways, eventually producing visible pigments to dissipate less energetic, but higher intensity, visible photons up to wavelengths of the "red edge". The propagation and dispersal of organic pigments, catalyzed by animals, and their coupling with abiotic dissipative processes, such as the water cycle, culminated in the apex photon dissipative structure, today's biosphere.

Ultrafast excitonic dynamics in DNA: Bridging correlated quantum dynamics and sequence dependence

After photoexcitation of DNA, the excited electron (in the LUMO) and the remaining hole (in the HOMO) localized on the same DNA base form a bound pair, called the Frenkel exciton, due to their mutual Coulomb interaction. In this study, we demonstrate that a tight-binding (TB) approach, using TB parameters for electrons and holes available in the literature, allows us to correlate relaxation properties, average charge separation, and dipole moments to a large ensemble of double-stranded DNA sequences (all 16384 possible sequences with 14 nucleobases). This way, we are able to identify a relatively small subensemble of sequences responsible for long-lived excited states, high average charge separation, and high dipole moment. Further analysis shows that these sequences are particularly T rich. By systematically screening the impact of electron-hole interaction (Coulomb forces), we verify that these correlations are relatively robust against finite-size variations of the interaction parameter, not directly accessible experimentally. This methodology combines simulation methods from quantum physics and physical chemistry with statistical analysis known from genetics and epigenetics, thus representing a powerful bridge to combine information from both fields.

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.

Galván I, Rivalta I, Aquilante F, Garavelli M, Lindh R, Segarra-Martí J. Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

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.

The Non-Equilibrium Thermodynamics of Natural Selection: From Molecules to the Biosphere

Evolutionary theory suggests that the origin, persistence, and evolution of biology is driven by the "natural selection" of characteristics improving the differential reproductive success of the organism in the given environment. The theory, however, lacks physical foundation, and, therefore, at best, can only be considered a heuristic narrative, of some utility for assimilating the biological and paleontological data at the level of the organism. On deeper analysis, it becomes apparent that this narrative is plagued with problems and paradoxes. Alternatively, non-equilibrium thermodynamic theory, derived from physical law, provides a physical foundation for describing material interaction with its environment at all scales. Here we describe a "natural thermodynamic selection" of characteristics of structures (or processes), based stochastically on increases in the global rate of dissipation of the prevailing solar spectrum. Different mechanisms of thermodynamic selection are delineated for the different biotic-abiotic levels, from the molecular level at the origin of life, up to the level of the present biosphere with non-linear coupling of biotic and abiotic processes. At the levels of the organism and the biosphere, the non-equilibrium thermodynamic description of evolution resembles, respectively, the Darwinian and Gaia descriptions, although the underlying mechanisms and the objective function of selection are fundamentally very different.

Photoionization Dynamics and Proton Transfer within the Adenine-Thymine Nucleobase Pair

Studying the stability of hydrogen-bonded nucleobase pairs, at the heart of the genetic code, is of utmost importance for an in-depth understanding of basic mechanisms of life and biomolecular evolution. We present here a VUV single photon ionization dynamic study of the nucleobase pair adenine-thymine (AT), revealing its ionization and dissociative ionization thresholds via double imaging electron/ion coincidence spectroscopy. The experimental data, consisting of cluster mass-resolved threshold photoelectron spectra and photon energy-dependent ion kinetic energy release distributions, allow the unambiguous distinction of the dissociation of AT into protonated adenine AH⁺ and a dehydrogenated thymine radical T(-H) from dissociative ionization processes of other nucleobase clusters. Comparison to high-level ab initio calculations indicates that our experimental observations can be explained by a single hydrogen-bonded conformer present in our molecular beam and allows the estimation of an upper limit of the barrier of the proton transfer in the ionized AT pair.

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.

Aggregation of nucleobases and metabolites: Adenine-theobromine trimers

The selection of cytosine, guanine, thymine, and adenine as components of the information biopolymers was a complex process influenced by several factors. Among them, the intermolecular interactions may have played a determinant role. Thus, a deep understanding of the intermolecular interactions between nucleobases and other prebiotic molecules may help understand the first instants of chemical evolution. Following this hypothesis, we present here a combined spectroscopic and computational study of theobromine₂-adenine and thebromine-adenine₂ trimers. While adenine is a nucleobase, theobromine was probably part of the prebiotic chemistry. The trimers were formed in jets and probed by a combination of UV and IR spectroscopic techniques. The spectra were interpreted in light of the predictions obtained using density-functional methods. The results suggest the existence of a subtle balance between formation of hydrogen bonds and π-π interactions. Thus, while theobromine₂-adenine tends to form complex in stacked structures, theobromine-adenine₂ prefers formation of planar structures, maximizing the interaction by hydrogen bonds. The small energy difference between planar and stacked structures highlights the importance of accurately modeling the dispersion forces in the functionals to produce reliable predictions.

Photochemistry of 2-thiooxazole: a plausible prebiotic precursor to RNA nucleotides

Potentially prebiotic chemical reactions leading to RNA nucleotides involve periods of UV irradiation, which are necessary to promote selectivity and destroy biologially irrelevant side products. Nevertheless, UV light has only been applied to promote specific stages of prebiotic reactions and its effect on complete prebiotic reaction sequences has not been extensively studied. Here, we report on an experimental and computational investigation of the photostability of 2-thiooxazole (2-TO), a potential precursor of pyrimidine and 8-oxopurine nucleotides on early Earth. Our UV-irradiation experiments resulted in rapid decomposition of 2-TO into unidentified small molecule photoproducts. We further clarify the underlying photochemistry by means of accurate ab initio calculations and surface hopping molecular dynamics simulations. Overall, the computational results show efficient rupture of the aromatic ring upon the photoexcitation of 2-TO via breaking of the C-O bond. Consequently, the initial stage of the divergent prebiotic synthesis of pyrimidine and 8-oxopurine nucleotides would require periodic shielding from UV light either with sun screening chromophores or through a planetary scenario that would protect 2-TO until it is transformed into a more stable intermediate compound, e.g. oxazolidinone thione.

Environmentally sensitive fluorescent nucleoside analogues as probes for nucleic acid - protein interactions: molecular design and biosensing applications

Fluorescent nucleoside analogues (FNAs) are indispensable in studying the interactions of nucleic acids with nucleic acid-binding proteins. By replacing one of the poorly emissive natural nucleosides, FNAs enable real-time optical monitoring of the binding interactions in solutions, under physiologically relevant conditions, with high sensitivity. Besides that, FNAs are widely used to probe conformational dynamics of biomolecular complexes using time-resolved fluorescence methods. Because of that, FNAs are tools of high utility for fundamental biological research, with potential applications in molecular diagnostics and drug discovery. Here I review the structural and physical factors that can be used for the conversion of the molecular binding events into a detectable fluorescence output. Typical environmentally sensitive FNAs, their properties and applications, and future challenges in the field are discussed.

Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, namely, equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the ππ* states (one for pyrimidine, two for purines) associated with the lowest-energy absorption band, and for the close-lying nπ* state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effects by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the nπ* state. Finally, ESA spectra have been computed from the minima of the lowest energy ππ* state, and found in line with the available experimental transient absorption spectra of the nucleosides in solution, providing further validation of our computational approach.

Molecular Identification of the Transient Species Mediating the Deactivation Dynamics of Solvated Guanosine and Deazaguanosine

Small structural alterations of the purine/pyrimidine core have been related to important photophysical changes, such as the loss of photostability. Similarly to canonical nucleobases, solute-solvent interactions can lead to a change in the excited state lifetimes and/or to the interplay of different states in the photophysics of these modified nucleobases. To shed light on both effects, we here report a complete picture of the absorption spectra and excited state deactivation of deoxyguanosine and its closely related derivative, deoxydeazaguanosine, in water and methanol through the mapping of the excited state potential energy surfaces and molecular dynamics simulations at the TD-DFT level of theory. We show that the N by CH exchange in the imidazole ring of deoxyguanosine translates into a small red-shift of the bright states and slightly faster dynamics. In contrast, changing solvent from water to methanol implies the opposite, i.e., that the deactivation of both systems to the ground state is significantly hindered.

Theoretical Study of Light-Induced Crosslinking Reaction Between Pyrimidine DNA Bases and Aromatic Amino Acids

Low-lying electronic excited states and their relaxation pathways as well as energetics of the crosslinking reaction between uracil as a model system for pyrimidine-type building blocks of DNA and RNA and benzene as a model system for aromatic groups of tyrosine (Tyr) and phenylalanine (Phe) amino acids have been studied in the framework of density functional theory. The equilibrium geometries of the ground and electronic excited states as well as the crossing points between the potential energy surfaces of the uracil–benzene complex were computed. Based on these results, different relaxation pathways of the electronic excited states that lead to either back to the initial geometry configuration or the dimerization between the six-membered rings of the uracil–benzene complex have been identified, and the energetic conditions for their occurrence are discussed. It can be concluded that the DNA–protein crosslinking reaction can be induced by the external electromagnetic field via the dimerization reaction between the six-membered rings of the uracil–benzene pair at the electronic excited-state level of the complex. In the case of the uracil–phenol complex, the configuration of the cyclic adduct (dimerized) conformation is less likely to be formed.

Revealing the Structure and Noncovalent Interactions of Isolated Molecules by Laser-Desorption/Ionization-Loss Stimulated Raman Spectroscopy and Quantum Calculations

The structural and dynamical characteristics of isolated molecules are essential, yet obtaining this information is difficult. We demonstrate laser-desorption jet-cooling/ionization-loss stimulated Raman spectroscopy to obtain Raman spectral signatures of nonvolatile molecules in the gas phase. The vibrational features of a test substance, the most abundant conformer of tryptamine, are compared and found to match those resulting from the scaled harmonic Raman spectrum obtained by density functional theory calculations. The vibrational signatures serve to identify the most prominent gauche conformer and evaluate its predicted electronic structure. These findings, together with noncovalent interaction (NCI) analysis, provide new insights into electron densities and reduced density gradients, assessing the hydrogen bonds (N-H···π and C-H···H-C) and interplay between attractive and repulsive NCIs affecting the structure. This approach accesses vibrational signatures of isolated nonvolatile molecules by tabletop lasers at uniform resolution and in a broad frequency range, promising great benefit to future studies.

Structural Reconstruction of Optically Invisible State in a Single Molecule via Scanning Tunneling Microscope

Molecular dark states, participating in various energy- and electron-transfer processes, are typically beyond direct optical-spectroscopic measurements because of the forbidden transition dictated by the selection rule. In this work, we demonstrate a direct profile of the dark-state transition density of a single molecule on the subnanometer scale by using a scanning tunneling microscope. Our method allows one to resolve the four-lobe configuration in a 1 nm region for the example molecule. The current proposal will bring about a new methodology to study the single-molecule properties in electro-optical devices and light-assisted biological processes.

A theoretical study on the excited-state deactivation paths for the A-5FU dimer

The photostability of DNA plays a key role in the normal function of organisms. A-5FU is a base pair derivative of the A-T dimer where the methyl group is replaced by a F atom. Here, accurate static TDDFT calculations and non-adiabatic dynamic simulations are used to systematically investigate the excited-state decay paths of the A-5FU dimer related to the proton transfer and the out-of-plane twisting deformation motion of A and 5FU in the 1ππ* and 1nπ* states. CC2 is used to check the accuracy of the current TDDFT calculations. Our results show that the deformation of the C[double bond, length as m-dash]C or C[double bond, length as m-dash]N double bond in A and 5FU provides an efficient pathway for the depopulation of the lowest excited states, which can compete with the excited-state proton transfer paths in the dimer. This finding indicates that monomer-like decay paths could be important for the photostability of weakly hydrogen-bonded DNA base pairs and provide a new insight into the excited-state decay paths in base pairs and their analogues.

Unravelling the full relaxation dynamics of superexcited helium nanodroplets

The relaxation dynamics of superexcited superfluid He nanodroplets is thoroughly investigated by means of extreme-ultraviolet (XUV) femtosecond electron and ion spectroscopy complemented by time-dependent density functional theory (TDDFT). Three main paths leading to the emission of electrons and ions are identified: droplet autoionization, pump-probe photoionization, and autoionization induced by re-excitation of droplets relaxing into levels below the droplet ionization threshold. The most abundant product ions are He2+, generated by droplet autoionization and by photoionization of droplet-bound excited He atoms. He+ appear with some pump-probe delay as a result of the ejection He atoms in their lowest excited states from the droplets. The state-resolved time-dependent photoelectron spectra reveal that intermediate excited states of the droplets are populated in the course of the relaxation, terminating in the lowest-lying metastable singlet and triplet He atomic states. The slightly faster relaxation of the triplet state compared to the singlet state is in agreement with the simulation showing faster formation of a bubble around a He atom in the triplet 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.

Ultrafast Dynamics of the Isolated Adenosine-5'-triphosphate Dianion Probed by Time-Resolved Photoelectron Imaging

The excited state dynamics of the doubly deprotonated dianion of adenosine-5'-triphosphate, [ATP-H₂]^2-, has been spectroscopically explored by time-resolved photoelectron spectroscopy following excitation at 4.66 eV. Time-resolved photoelectron spectra show that two competing processes occur for the initially populated ¹ππ* state. The first is rapid electron emission by tunneling through a repulsive Coulomb barrier as the ¹ππ* state is a resonance. The second is nuclear motion on the ¹ππ* state surface leading to an intermediate that no longer tunnels and subsequently decays by internal conversion to the ground electronic state. The spectral signatures of the features are similar to those observed for other adenine-derivatives, suggesting that this nucleobase is quite insensitive to the nearby negative charges localized on the phosphates, except of course for the appearance of the additional electron tunneling channel, which is open in the dianion.

Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

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.

Nonadiabatic Absorption Spectra and Ultrafast Dynamics of DNA and RNA Photoexcited Nucleobases

We have recently proposed a protocol for Quantum Dynamics (QD) calculations, which is based on a parameterisation of Linear Vibronic Coupling (LVC) Hamiltonians with Time Dependent (TD) Density Functional Theory (TD-DFT), and exploits the latest developments in multiconfigurational TD-Hartree methods for an effective wave packet propagation. In this contribution we explore the potentialities of this approach to compute nonadiabatic vibronic spectra and ultrafast dynamics, by applying it to the five nucleobases present in DNA and RNA. For all of them we computed the absorption spectra and the dynamics of ultrafast internal conversion (100 fs timescale), fully coupling the first 2-3 bright states and all the close by dark states, for a total of 6-9 states, and including all the normal coordinates. We adopted two different functionals, CAM-B3LYP and PBE0, and tested the effect of the basis set. Computed spectra are in good agreement with the available experimental data, remarkably improving over pure electronic computations, but also with respect to vibronic spectra obtained neglecting inter-state couplings. Our QD simulations indicate an effective population transfer from the lowest energy bright excited states to the close-lying dark excited states for uracil, thymine and adenine. Dynamics from higher-energy states show an ultrafast depopulation toward the more stable ones. The proposed protocol is sufficiently general and automatic to promise to become useful for widespread applications.

Long living excited state of protonated adenosine unveiled by ultrafast fluorescence spectroscopy and density functional theoretical study

Adenosine (Ado) possesses ultrafast nonradiative dynamics accounting for its remarkably high photostability. The deactivation dynamics of Ado after protonation in an aqueous solution remains an elusive issue. Herein we report an investigation of the excited state dynamics of protonated Ado (AdoH⁺) performed using ultrafast time-resolved fluorescence spectroscopy combined with density functional theoretical calculation. The result obtained from comparison of conformers with protonation at different sites revealed that the syn-conformer with protonation occurring at the N3 position (syn-N3) is the predominant form of AdoH⁺ in the ground state, similar to that of Ado. In contrast, the fluorescence of AdoH⁺ with maximum intensity at 385 nm, significantly red-shifted from that of Ado, displaying decay dynamics composed of an ultrafast component with the lifetime of ∼0.5 ps and a slower one of ∼2.9 ns. The former is because of the decay of the syn-N3 conformer, similar to that reported for AdoH⁺ under the gas phase condition. The latter is due to the syn-N1 conformer formed via ultrafast proton transfer of the syn-N3. The excited state of syn-N1 has a peculiar nonplanar conformation over the purine molecule, which is responsible for the substantial Stokes shift showed in the fluorescence spectrum and correlates with a large energy barrier for nonradiative decay likely involving a reversed proton transfer. This study demonstrates the importance of protonation and solvent environment in altering dramatically the excited states of Ado, providing insight for better understanding nonradiative dynamics of both the monomeric bases and the oligomeric or polymeric DNAs.

The effect of hydrogen bonds on the ultrafast relaxation dynamics of a BODIPY dimer

The influence of hydrogen bonds (H-bonds) in the structure, dynamics, and functionality of biological and artificial complex systems is the subject of intense investigation. In this broad context, particular attention has recently been focused on the ultrafast H-bond dependent dynamical properties in the electronic excited state because of their potentially dramatic consequences on the mechanism, dynamics, and efficiency of photochemical reactions and photophysical processes of crucial importance for life and technology. Excited-state H-bond dynamics generally occur on ultrafast time scales of hundreds of femtoseconds or less, making the characterization of associated mechanisms particularly challenging with conventional time-resolved techniques. Here, 2D electronic spectroscopy is exploited to shed light on this still largely unexplored dynamic mechanism. An H-bonded molecular dimer prepared by self-assembly of two boron-dipyrromethene dyes has been specifically designed and synthesized for this aim. The obtained results confirm that upon formation of H-bonds and the dimer, a new ultrafast relaxation channel is activated in the ultrafast dynamics, mediated by the vibrational motions of the hydrogen donor and acceptor groups. This relaxation channel also involves, beyond intra-molecular relaxations, an inter-molecular transfer process. This is particularly significant considering the long distance between the centers of mass of the two molecules. These findings suggest that the design of H-bonded structures is a particularly powerful tool to drive the ultrafast dynamics in complex materials.

Ultrafast excited-state dynamics of promising nucleobase ancestor 2,4,6-triaminopyrimidine

The ultrafast excited-state dynamics of 2,4,6-triaminopyrimidine - thought to be a promising candidate for a proto-RNA nucleobase - have been investigated via static multireference quantum-chemical calculations and mixed-quantum-classical/trajectory surface-hopping dynamics with a focus on the lowest-lying electronic states of the singlet manifold and with a view towards understanding the UV(C)/UV(B) photostability of the molecule. Ultrafast internal conversion channels have been identified that connect the lowest-lying ππ* electronically-excited state of 2,4,6-triaminopyrimidine with the ground electronic state, and non-radiative decay has been observed to take place on the picosecond timescale via a ππ* out-of-plane NH2 ("oop-NH2") minimum-energy crossing point. The short excited-state lifetime is competitive with the excited-state lifetimes of the canonical pyrimidine nucleobases, affirming the promise of 2,4,6-triaminopyrimidine as an ancestor. Evidence for energy-dependent excited-state dynamics is presented, and the open question of intersystem crossing is discussed speculatively.

The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine

The non-equilibrium thermodynamics and the photochemical reaction mechanisms are described which may have been involved in the dissipative structuring, proliferation and complexation of the fundamental molecules of life from simpler and more common precursors under the UVC photon flux prevalent at the Earth’s surface at the origin of life. Dissipative structuring of the fundamental molecules is evidenced by their strong and broad wavelength absorption bands in the UVC and rapid radiationless deexcitation. Proliferation arises from the auto- and cross-catalytic nature of the intermediate products. Inherent non-linearity gives rise to numerous stationary states permitting the system to evolve, on amplification of a fluctuation, towards concentration profiles providing generally greater photon dissipation through a thermodynamic selection of dissipative efficacy. An example is given of photochemical dissipative abiogenesis of adenine from the precursor HCN in water solvent within a fatty acid vesicle floating on a hot ocean surface and driven far from equilibrium by the incident UVC light. The kinetic equations for the photochemical reactions with diffusion are resolved under different environmental conditions and the results analyzed within the framework of non-linear Classical Irreversible Thermodynamic theory.

Spectroscopic identification of fragment ions of DNA/RNA building blocks: the case of pyrimidine

Pyrimidine (Pym, 1,3-diazine, 1,3-diazabenzene) is an important N-heterocyclic building block of nucleobases. Understanding the structures of its fragment and precursor ions provides insight into its prebiotic and abiotic synthetic route. The long-standing controversial debate about the structures of the primary fragment ions of the Pym+ cation (C4H4N2+, m/z 80) resulting from loss of HCN, C3H3N+ (m/z 53), is closed herein with the aid of a combined approach utilizing infrared photodissociation (IRPD) spectroscopy in the CH and NH stretch ranges (νCH/NH) and density functional theory (DFT) calculations. IRPD spectra of cold Ar/N2-tagged fragment ions reveal that the C3H3N+ population is dominated by cis-/trans-HCCHNCH+ ions (∼90%) along with a minor contribution of the most stable H2CCCNH+ and cis-/trans-HCCHCNH+ isomers (∼10%). We also spectroscopically confirm that the secondary fragment resulting from further loss of HCN, C2H2+ (m/z 26), is the acetylene cation (HCCH+). The spectroscopic characterization of the identified C3H3N+ isomers and their hydrogen-bonded dimers with Ar and N2 provides insight into the acidity of their CH and NH groups. Finally, the vibrational properties of Pym+ in the 3 μm range are probed by IRPD of Pym+-(N2)1-2 clusters, which shows a high π-binding affinity of Pym+ toward a nonpolar hydrophobic ligand. Its νCH spectrum confirms the different acidity of the three nonequivalent CH groups.

Observation of Enhanced Dissociative Photochemistry in the Non-Native Nucleobase 2-Thiouracil

We present the first study to measure the dissociative photochemistry of 2-thiouracil (2-TU), an important nucleobase analogue with applications in molecular biology and pharmacology. Laser photodissociation spectroscopy is applied to the deprotonated and protonated forms of 2-TU, which are produced in the gas-phase using electrospray ionization mass spectrometry. Our results show that the deprotonated form of 2-thiouracil ([2-TU-H]⁻) decays predominantly by electron ejection and hence concomitant production of the [2-TU-H]· free-radical species, following photoexcitation across the UVA-UVC region. Thiocyanate (SCN⁻) and a m/z 93 fragment ion are also observed as photodecay products of [2-TU-H]⁻ but at very low intensities. Photoexcitation of protonated 2-thiouracil ([2-TU·H]⁺) across the same UVA-UVC spectral region produces the m/z 96 cationic fragment as the major photofragment. This ion corresponds to ejection of an HS· radical from the precursor ion and is determined to be a product of direct excited state decay. Fragment ions associated with decay of the hot ground state (i.e., the ions we would expect to observe if 2-thiouracil was behaving like UV-dissipating uracil) are observed as much more minor products. This behaviour is consistent with enhanced intersystem crossing to triplet excited states compared to internal conversion back to the ground state. These are the first experiments to probe the effect of protonation/deprotonation on thionucleobase photochemistry, and hence explore the effect of pH at a molecular level on their photophysical properties.

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.

Excited State Deactivation Mechanism in Protonated Uracil: New Insights from Theoretical Studies

We have conducted here a theoretical exploration, discussing the distinct excited state lifetimes reported experimentally for the two lowest lying protonated isomers of uracil. In this regard, the first-principal computational levels as well as the nonadiabatic surface hopping dynamics have been employed. It has been revealed that relaxation of the ¹ππ* state of enol-enol form (EE⁺) to the ground is barrier-free via out-of-plane coordinates, resulting in an ultrashort S₁ lifetime of this species. For the second most stable isomer (EK⁺), however, a significant barrier predicted in the CASPT2 S₁ potential energy profile along the twisting coordinate has been proposed to explain the relevant long lifetime reported experimentally.

Relaxation Dynamics of Hydrated Thymine, Thymidine, and Thymidine Monophosphate Probed by Liquid Jet Time-Resolved Photoelectron Spectroscopy

The relaxation dynamics of thymine and its derivatives thymidine and thymidine monophosphate are studied using time-resolved photoelectron spectroscopy applied to a water microjet. Two absorption bands are studied; the first is a bright ππ* state which is populated using tunable-ultraviolet light in the range 4.74-5.17 eV and probed using a 6.20 eV probe pulse. By reversing the order of these pulses, a band containing multiple ππ* states is populated by the 6.20 eV pulse and the lower energy pulse serves as the probe. The lower lying ππ* state is found to decay in ∼400 fs in both thymine and thymidine independent of pump photon energy, while thymidine monophosphate decays vary from 670 to 840 fs with some pump energy dependence. The application of a computational quantum mechanical/molecular mechanical scheme at the XMS-CASPT2//CASSCF/AMBER level of theory suggests that conformational differences existing between thymidine and thymidine monophosphate in solution account for this difference. The higher lying ππ* band is found to decay in ∼600 fs in all three cases, but it is only able to be characterized when the 5.17 eV probe pulse is used. Notably, no long-lived signal from an nπ* state can be identified in either experiment on any of the three molecules.

UV photostability of three 2-aminoazoles with key roles in prebiotic chemistry on the early earth

Three related molecules in the 2-aminoazole family are potentially important for prebiotic chemistry: 2-aminooxazole, 2-aminoimidazole, and 2-aminothiazole, which can provide critical functions as an intermediate in nucleotide synthesis, a nucleotide activating agent, and a selective agent, respectively. Here, we examine the wavelength-dependent photodegradation of these three molecules under mid-range UV light (210-290 nm). We then assess the implications of the observed degradation rates for the proposed prebiotic roles of these compounds. We find that all three 2-aminoazoles degrade under UV light, with half lives ranging from ≈7-100 hours under a solar-like spectrum. 2-Aminooxazole is the least photostable, while 2-aminoimidazole is the most photostable. The relative photostabilities are consistent with the order in which these molecules would be used prebiotically: AO is used first to build nucleotides and AI is used last to activate them.

Comprehensive model for X-ray-induced damage in protein crystallography

Acquisition of X-ray crystallographic data is always accompanied by structural degradation owing to the absorption of energy. The application of high-fluency X-ray sources to large biomolecules has increased the importance of finding ways to curtail the onset of X-ray-induced damage. A significant effort has been under way with the aim of identifying strategies for protecting protein structure. A comprehensive model is presented that has the potential to explain, both qualitatively and quantitatively, the structural changes induced in crystalline protein at ∼100 K. The first step is to consider the qualitative question: what are the radiation-induced intermediates and expected end products? The aim of this paper is to assist in optimizing these strategies through a fundamental understanding of radiation physics and chemistry, with additional insight provided by theoretical calculations performed on the many schemes presented.

Excited-state dynamics of mononucleotides and DNA strands in a deep eutectic solvent

The photophysics of several mono- and oligonucleotides were investigated in a deep eutectic solvent for the first time. The solvent glyceline, prepared as a 1 : 2 mole ratio mixture of choline chloride and glycerol, was used to study excited-state deactivation in a non-aqueous solvent by the use of steady-state and time-resolved spectroscopy. DNA strands in glyceline retain the secondary structures that are present in aqueous solution to some degree, thus enabling a study of the effects of solvent properties on the excited states of stacked bases and stacked base pairs. The excited-state lifetime of the mononucleotide 5'-AMP in glyceline is 630 fs, or twice as long as in aqueous solution. Even slower relaxation is seen for 5'-TMP in glyceline, and a possible triplet state with a lifetime greater than 3 ns is observed. Circular dichroism spectra show that the single strand (dA)18 and the duplex d(AT)9·d(AT)9 adopt similar structures in glyceline and in aqueous solution. Despite having similar conformations in both solvents, femtosecond transient absorption experiments reveal striking changes in the dynamics. Excited-state decay and vibrational cooling generally take place more slowly in glyceline than in water. Additionally, the fraction of long-lived excited states in both oligonucleotide systems is lower in glyceline than in aqueous solution. For a DNA duplex, water is suggested to favor decay pathways involving intrastrand charge separation, while the deep eutectic solvent favors interstrand deactivation channels involving neutral species. Slower solvation dynamics in the viscous deep eutectic solvent may also play a role. These results demonstrate that the dynamics of excitations in stacked bases and stacked base pairs depend not only on conformation, but are also highly sensitive to the solvent.

Initial excited-state structural dynamics of 2′-deoxyadenosine

Purine nucleobases (adenosine and guanosine) are prone to spontaneous breaking of the nucleosidic bond to form abasic sites in both DNA and RNA. However, the purines also undergo photochemical reactions, including oxidation and cycloaddition, to form damage sites, albeit at lower photochemical quantum yields than the pyrimidines. In this study, we use ultraviolet resonance Raman spectroscopy to measure the initial excited-state structural dynamics in the nucleoside, 2′-deoxyadenosine. The resonance Raman-derived initial excited-state structural dynamics throughout the 260 nm La excited electronic state of adenine are found to be smaller in the nucleoside than in the previously reported 9-methyladenine nucleobase derivative, consistent with what is found for the pyrimidines thymine and uracil. Interestingly, resonance-enhanced vibrational modes in this electronic state also contain internal coordinates localized on the sugar, which may represent a different energy dissipation mechanism than in the pyrimidine nucleosides. The results will be discussed in terms of the initial excited-state photophysics and photochemistry of DNA and RNA.

State-Dependent Fragmentation of Protonated Uracil and Uridine

Action spectroscopy using photon excitation in the VUV range (photon energy 4.5-9 eV) was performed on protonated uracil (UraH⁺) and uridine (UrdH⁺). The precursor ions with m/ z 113 and m/ z 245, respectively, were produced by an electrospray ionization source and accumulated inside a quadrupole ion trap mass spectrometer. After irradiation with tunable synchrotron radiation, product ion mass spectra were obtained. Fragment yields as a function of excitation energy show several maxima that can be attributed to the photoexcitation into different electronic states. For uracil, vertically excited states were calculated using the equation-of-motion coupled cluster approach and compared to the observed maxima. This allows to establish correlations between electronic states and the resulting fragment masses and can thus help to disentangle the complex de-excitation and fragmentation pathways of nucleic acid building blocks. Photofragmentation of the nucleoside uridine shows a significantly lower variety of fragments, indicating stabilization of the nucleobase by the attached sugar.

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.

Imaging the Renner-Teller effect using laser-induced electron diffraction

Structural information on electronically excited neutral molecules can be indirectly retrieved, largely through pump-probe and rotational spectroscopy measurements with the aid of calculations. Here, we demonstrate the direct structural retrieval of neutral carbonyl disulfide (CS₂) in the [Formula: see text] excited electronic state using laser-induced electron diffraction (LIED). We unambiguously identify the ultrafast symmetric stretching and bending of the field-dressed neutral CS₂ molecule with combined picometer and attosecond resolution using intrapulse pump-probe excitation and measurement. We invoke the Renner-Teller effect to populate the [Formula: see text] excited state in neutral CS₂, leading to bending and stretching of the molecule. Our results demonstrate the sensitivity of LIED in retrieving the geometric structure of CS₂, which is known to appear as a two-center scatterer.

Probing photodissociation dynamics using ring polymer molecular dynamics

The performance of the ring polymer molecular dynamics (RPMD) approach to simulate typical photodissociation processes is assessed. The correct description of photodissociation requires the calculation of correlation functions or expectation values associated with non-equilibrium initial conditions, which was shown to be possible with RPMD very recently [J. Chem. Phys. 145, 204118 (2016)]. This approach is combined with treatment of the nonadiabatic dynamics employing the ring polymer surface hopping approach (RPSH), which is based on Tully's fewest switches surface hopping (FSSH) approach. The performance is tested using one-dimensional photodissociation models. It is found that RPSH with non-equilibrium initial conditions can well reproduce the time-dependent dissociation probability, and adiabatic and diabatic populations for cases where the crossing point is below and above the Franck-Condon point, respectively, while standard FSSH fails to reproduce the exact quantum dynamics in the latter case. Thus, it is shown that RPSH is an efficient and accurate alternative to standard FSSH, which is one of the most widely employed approaches to study photochemistry.