The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study

Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.

The photophysics of protonated cytidine and hemiprotonated cytidine base pair: A computational study

We here study the effect that a lowering of the pH has on the excited state processes of cytidine and a cytidine/cytidine pair in solution, by integrating time-dependent density functional theory and CASSCF/CASPT2 calculations, and including solvent by a mixed discrete/continuum model. Our calculations reproduce the effect of protonation at N3 on the steady-state infrared and absorption spectra of a protonated cytidine (CH+ ), and predict that an easily accessible non-radiative deactivation route exists for the spectroscopic state, explaining its sub-ps lifetime. Indeed, an extremely small energy barrier separates the minimum of the lowest energy bright state from a crossing region with the ground electronic state, reached by out-of-plane motion of the hydrogen substituents of the CC double bond, the so-called ethylenic conical intersection typical of cytidine and other pyrimidine bases. This deactivation route is operative for the two bases forming an hemiprotonated cytidine base pair, [CH·C]+ , the building blocks of I-motif secondary structures, whereas interbase processes play a minor role. N3 protonation disfavors instead the nπ* transitions, associated with the long-living components of cytidine photoactivated dynamics.

Vibrational relaxation by methylated xanthines in solution: Insights from 2D IR spectroscopy and calculations

Two-dimensional infrared (2D IR) spectroscopy, infrared pump-infrared probe spectroscopy, and density functional theory calculations were used to study vibrational relaxation by ring and carbonyl stretching modes in a series of methylated xanthine derivatives in acetonitrile and deuterium oxide (heavy water). Isotropic signals from the excited symmetric and asymmetric carbonyl stretch modes decay biexponentially in both solvents. Coherent energy transfer between the symmetric and asymmetric carbonyl stretching modes gives rise to a quantum beat in the time-dependent anisotropy signals. The damping time of the coherent oscillation agrees with the fast decay component of the carbonyl bleach recovery signals, indicating that this time constant reflects intramolecular vibrational redistribution (IVR) to other solute modes. Despite their similar frequencies, the excited ring modes decay monoexponentially with a time constant that matches the slow decay component of the carbonyl modes. The slow decay times, which are faster in heavy water than in acetonitrile, approximately match the ones observed in previous UV pump-IR probe measurements on the same compounds. The slow component is assigned to intermolecular energy transfer to solvent bath modes from low-frequency solute modes, which are populated by IVR and are anharmonically coupled to the carbonyl and ring stretch modes. 2D IR measurements indicate that the carbonyl stretching modes are weakly coupled to the delocalized ring modes, resulting in slow exchange that cannot explain the common solvent-dependence. IVR is suggested to occur at different rates for the carbonyl vs ring modes due to differences in mode-specific couplings and not to differences in the density of accessible states.

Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study

Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)₅ sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.

Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform

We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.

Nucleic Acids as a Playground for the Computational Study of the Photophysics and Photochemistry of Multichromophore Assemblies

ConspectusThe interaction between light and multichromophoric assemblies (MCAs) is the primary event of many fundamental processes, from photosynthesis to organic photovoltaics, and it triggers dynamical processes that share remarkable similarities at the molecular scale: light absorption, energy and charge transfer, internal conversions, emission, and so on. Those events often involve many chromophores and different excited electronic states that are coupled on an ultrafast time scale. This Account aims to discuss some of the chemical physical effects ruling these processes, a fundamental step toward their control, based on our experience on nucleic acids.In the last 15 years, we have, indeed, studied the photophysics and photochemistry of DNA and its components. By combining different quantum mechanical methods, we investigated the molecular processes responsible for the damage of the genetic code or, on the contrary, those preventing it by dissipating the excess energy deposited in the system by UV absorption. Independently of its fundamental biological role, DNA, with its fluctuating closely stacked bases stabilized by weak nonbonding interactions, can be considered a prototypical MCA. Therefore, it allows one to tackle within a single system many of the conceptual and methodological challenges involved in the study of photoinduced processes in MCA.In this Account, by using the outcome of our studies on oligonucleotides as a guideline, we thus highlight the most critical modellistic issues to be faced when studying, either experimentally or computationally, the interaction between UV light and DNA and, at the same time, bring out their general relevance for the study of MCAs.We first discuss the rich photoactivated dynamics of nucleobases (the chromophores), highlighting the main effects modulating the interplay between their excited states and how the latter can affect the photoactivated dynamics of the polynucleotides, either providing effective monomer-like nonradiative decay routes or triggering reactive processes (e.g., triplet generation).We then tackle the reaction paths involving multiple bases, showing that in the DNA duplex the most important ones involve two stacked bases, forming a neutral excimer or a charge transfer (CT) state, which exhibit different spectral signatures and photochemical reactivity. In particular, we analyze the factors affecting the dynamic equilibrium between the excimer and CT, such as the fluctuations of the backbone or the rearrangement of the solvent.Next, we highlight the importance of the effects not directly connected to the chromophores, such as the flexibility of the backbone or the solvent effect. The former, affecting the stacking geometry of the bases, can determine the preference between different deactivation paths. The latter is particularly influential for CT states, making very important an accurate treatment of dynamical solvation effects, involving both the solvent bulk and specific solute-solvent interactions.In the last section, we describe the main methodological challenges related to the study of polynucleotide excited states and stress the benefits derived by the integration of complementary approaches, both computational and experimental. Only exploiting different point of views, in our opinion, it is possible to shed light on the complex phenomena triggered by light absorption in DNA, as in every MCA.

Thienoguanosine brightness in DNA duplexes is governed by the localization of its ππ* excitation in the lowest energy absorption band

Thienoguanosine (thG) is an isomorphic fluorescent guanosine (G) surrogate, which almost perfectly mimics the natural G in DNA duplexes and may therefore be used to sensitively investigate for example protein-induced local conformational changes. To fully exploit the information given by the probe, we carefully re-investigated the thG spectroscopic properties in 12-bp duplexes, when the Set and Ring Associated (SRA) domain of UHRF1 flips its 5' flanking methylcytosine (mC). The SRA-induced flipping of mC was found to strongly increase the fluorescence intensity of thG, but this increase was much larger when thG was flanked in 3' by a C residue as compared to an A residue. Surprisingly, the quantum yield and fluorescence lifetime values of thG were nearly constant, regardless of the presence of SRA and the nature of the 3' flanking residue, suggesting that the differences in fluorescence intensities might be related to changes in absorption properties. We evidenced that thG lowest energy absorption band in the duplexes can be deconvoluted into two bands peaking at ~350 nm and ~310 nm, respectively red-shifted and blue-shifted, compared to the spectrum of thG monomer. Using quantum mechanical calculations, we attributed the former to a nearly pure * excitation localized on thG and the latter to excited states with charge transfer character. The amplitude of thG red-shifted band strongly increased when its 3' flanking C residue was replaced by an A residue in the free duplex, or when its 5' flanking mC residue was flipped by SRA. As only the species associated with the red-shifted band were found to be emissive, the highly unusual finding of this work is that the brightness of thG in free duplexes as well as its changes on SRA-induced mC flipping almost entirely depend on the relative population and/or absorption coefficient of the red-shifted absorbing species.

The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model

In this contribution we present a quantum dynamical study of the photoexcited hydrogen bonded base pair adenine–thymine (AT) in a Watson–Crick arrangement. To that end, we parametrize Linear Vibronic Coupling (LVC) models with Time-Dependent Density Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization scheme (FrD) we have developed to define diabatic states on the basis of individual chromophores in a multichromophoric system. Wavepacket propagations were run with the multilayer extension of the Multiconfiguration Time-Dependent Hartree method. We considered excitations to the three lowest bright states, a ππ* state of thymine and two ππ* states (L_a and L_b) of adenine, and we found that on the 100 fs time scale the main decay pathways involve intramonomer population transfers toward nπ* states of the same nucleobase. In AT this transfer is less effective than in the isolated nucleobases, because hydrogen bonding destabilizes the nπ* states. The population transfer to the A → T charge transfer state is negligible, making the ultrafast (femtosecond) decay through the proton coupled electron transfer mechanism unlikely, in line with experimental results in apolar solvents. The excitation energy transfer is also very small. We carefully compare the predictions of LVC Hamiltonians obtained with different sets of diabatic states, defined so to match either local states of the two separated monomers or the base pair adiabatic states in the Franck–Condon region. To that end we also extend the flexibility of the FrD-LVC approach, introducing a new strategy to define fragments diabatic states that account for the effect of the rest of the multichromohoric system through a Molecular Mechanics potential.

Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge-Transfer State in the Photoexcited Guanine-Cytosine Pair

We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a linear vibronic coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intrachromophore decays, like those at conical intersections, and interchromophore processes, like exciton localization/delocalization and the involvement of charge-transfer (CT) states. We used FrD-LVC parametrized with time-dependent density functional theory (TD-DFT) calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a guanine-cytosine (GC) hydrogen-bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all of the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G → C CT state, which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ∼50% of that on C decay to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show that it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.

Photoinduced DNA Lesions in Dormant Bacteria: The Peculiar Route Leading to Spore Photoproducts Characterized by Multiscale Molecular Dynamics*

Some bacterial species enter a dormant state in the form of spores to resist to unfavorable external conditions. Spores are resistant to a wide series of stress agents, including UV radiation, and can last for tens to hundreds of years. Due to the suspension of biological functions, such as DNA repair, they accumulate DNA damage upon exposure to UV radiation. Differently from active organisms, the most common DNA photoproducts in spores are not cyclobutane pyrimidine dimers, but rather the so-called spore photoproducts. This noncanonical photochemistry results from the dry state of DNA and its binding to small, acid-soluble proteins that drastically modify the structure and photoreactivity of the nucleic acid. Herein, multiscale molecular dynamics simulations, including extended classical molecular dynamics and quantum mechanics/molecular mechanics based dynamics, are used to elucidate the coupling of electronic and structural factors that lead to this photochemical outcome. In particular, the well-described impact of the peculiar DNA environment found in spores on the favored formation of the spore photoproduct, given the small free energy barrier found for this path, is rationalized. Meanwhile, the specific organization of spore DNA precludes the photochemical path that leads to cyclobutane pyrimidine dimer formation.

What Makes Thienoguanosine an Outstanding Fluorescent DNA Probe?

Thienoguanosine (^thG) is an isomorphic guanosine (G) surrogate that almost perfectly mimics G in nucleic acids. To exploit its full potential and lay the foundation for future applications, 20 DNA duplexes, where the bases facing and neighboring ^thG were systematically varied, were thoroughly studied using fluorescence spectroscopy, molecular dynamics simulations, and mixed quantum mechanical/molecular mechanics calculations, yielding a comprehensive understanding of its photophysics in DNA. In matched duplexes, ^thG's hypochromism was larger for flanking G/C residues but its fluorescence quantum yield (QY) and lifetime values were almost independent of the flanking bases. This was attributed to high duplex stability, which maintains a steady orientation and distance between nucleobases, so that a similar charge transfer (CT) mechanism governs the photophysics of ^thG independently of its flanking nucleobases. ^thG can therefore replace any G residue in matched duplexes, while always maintaining similar photophysical features. In contrast, the local destabilization induced by a mismatch or an abasic site restores a strong dependence of ^thG's QY and lifetime values on its environmental context, depending on the CT route efficiency and solvent exposure of ^thG. Due to this exquisite sensitivity, ^thG appears ideal for monitoring local structural changes and single nucleotide polymorphism. Moreover, ^thG's dominant fluorescence lifetime in DNA is unusually long (9-29 ns), facilitating its selective measurement in complex media using a lifetime-based or a time-gated detection scheme. Taken together, our data highlight ^thG as an outstanding emissive substitute for G with good QY, long fluorescence lifetimes, and exquisite sensitivity to local structural changes.

Dynamics of the excited-state hydrogen transfer in a (dG)·(dC) homopolymer: intrinsic photostability of DNA

The intrinsic photostability of nucleic acids is intimately related to evolution of life, while its understanding at the molecular and electronic levels remains a challenge for modern science. Among the different decay pathways proposed in the last two decades, the excited-state hydrogen transfer between guanine-cytosine base pairs has been identified as an efficient non-reactive channel to dissipate the excess of energy provided by light absorption. The present work studies the dynamics of such phenomena taking place in a (dG)·(dC) B-DNA homopolymer in water solution using state-of-the-art molecular modelling and simulation methods. A dynamic effect that boosts the photostability of the inter-strand hydrogen atom transfers, inherent to the Watson-Crick base pairing, is unveiled and ascribed to the energy released during the proton transfer step. Our results also reveal a novel mechanism of DNA decay named four proton transfer (FPT), in which two protons of two adjacent G-C base pairs are transferred to form a biradical zwitterionic intermediate. Decay of the latter intermediate to the ground state triggers the transfer of the protons back to the guanine molecules recovering the Watson-Crick structure of the tetramer. This FPT process is activated by the close interaction of a nearby Na+ counterion with the oxygen atoms of the guanine nucleobases and hence represents a photostable channel operative in natural nucleic acids.

Radicals generated in alternating guanine-cytosine duplexes by direct absorption of low-energy UV radiation

Recent studies have evidenced that oxidatively damaged DNA, which potentially leads to carcinogenic mutations and aging, may result from the direct absorption of low-energy photons (>250 nm). Herein, the primary species, i.e., ejected electrons and base radicals associated with such damage in duplexes with an alternating guanine-cytosine sequence are quantified by nanosecond transient absorption spectroscopy. The one-photon ionization quantum yield at 266 nm is 1.2 × 10-3, which is similar to those reported previously for adenine-thymine duplexes. This means that the simple presence of guanine, the nucleobase with the lowest ionization potential, does not affect photo-ionization. The transient species detected after 3 μs are identified as deprotonated guanine radicals, which decay with a half-time of 2.5 ms. Spectral assignment is made with the help of quantum chemistry calculations (TD-DFT), which for the first time, provide reference absorption spectra for guanine radicals in duplexes. In addition, our computed spectra predict the changes in transient absorption expected for hole localization as well as deprotonation (to cytosine and bulk water) and hydration of the radical cation.