Electronic excitations and spectra in single-stranded DNA

Exciton decay mechanism in DNA single strands: back-electron transfer and ultrafast base motions

The photochemistry of DNA systems is characterized by the ultraviolet (UV) absorption of π-stacked nucleobases, resulting in exciton states delocalized over several bases. As their relaxation sensitively depends on local stacking conformations, disentangling the ensuing electronic and structural dynamics has remained an experimental challenge, despite their fundamental role in protecting the genome from potentially harmful UV radiation. Here we use transient absorption and transient absorption anisotropy spectroscopy with broadband femtosecond deep-UV pulses (250–360 nm) to resolve the exciton dynamics of UV-excited adenosine single strands under physiological conditions. Due to the exceptional deep-UV bandwidth and polarization sensitivity of our experimental approach, we simultaneously resolve the population dynamics, charge-transfer (CT) character and conformational changes encoded in the UV transition dipoles of the π-stacked nucleotides. Whilst UV excitation forms fully charge-separated CT excitons in less than 0.3 ps, we find that most decay back to the ground state via a back-electron transfer. Based on the anisotropy measurements, we propose that this mechanism is accompanied by a structural relaxation of the photoexcited base-stack, involving an inter-base rotation of the nucleotides. Our results finally complete the exciton relaxation mechanism for adenosine single strands and offer a direct view into the coupling of electronic and structural dynamics in aggregated photochemical systems.

Despite its key role in DNA photochemistry, the decay mechanism of excitons in stacked bases has remained difficult to resolve. Ultrafast polarization spectroscopy now reveals a back-electron transfer and ultrafast base motions in adenosine strands.

Excimer Intermediates en Route to Long-Lived Charge-Transfer States in Single-Stranded Adenine DNA as Revealed by Nonadiabatic Dynamics

The ultrafast time evolution of a single-stranded adenine DNA is studied using a hybrid multiscale quantum mechanics/molecular mechanics (QM/MM) scheme coupled to nonadiabatic surface hopping dynamics. As a model, we use (dA)₂₀ where a stacked adenine tetramer is treated quantum chemically. The dynamical simulations combined with on-the-fly quantitative wave function analysis evidence the nature of the long-lived electronically excited states formed upon absorption of UV light. After a rapid decrease of the initially excited excitons, relaxation to monomer-like states and excimers occurs within 100 fs. The former monomeric states then relax into additional excimer states en route to forming stabilized charge-transfer states on a longer timescale of hundreds of femtoseconds. The different electronic-state characters is reflected on the spatial separation between the adenines: excimers and charge-transfer states show a much smaller spatial separation than the monomer-like states and the initially formed excitons.

Excitation energy transport in DNA modelled by multi-chromophoric field-induced surface hopping

Revealing the extended excited state lifetime due to excitation energy transport in DNA by multi-chromophoric field-induced surface-hopping (McFISH).

Double-buffered, heterogeneous CPU + GPU integral digestion algorithm for single-excitation calculations involving a large number of excited states

The most widely used quantum-chemical models for excited states are single-excitation theories, a category that includes configuration interaction with single substitutions, time-dependent density functional theory, and also a recently developed ab initio exciton model. When a large number of excited states are desired, these calculations incur a significant bottleneck in the "digestion" step in which two-electron integrals are contracted with density or density-like matrices. We present an implementation that moves this step onto graphical processing units (GPUs), and introduce a double-buffer scheme that minimizes latency by computing integrals on the central processing units (CPUs) concurrently with their digestion on the GPUs. An automatic code generation scheme simplifies the implementation of high-performance GPU kernels. For the exciton model, which requires separate excited-state calculations on each electronically coupled chromophore, the heterogeneous implementation described here results in speedups of 2-6× versus a CPU-only implementation. For traditional time-dependent density functional theory calculations, we obtain speedups of up to 5× when a large number of excited states is computed. © 2018 Wiley Periodicals, Inc.

Electronic delocalization, charge transfer and hypochromism in the UV absorption spectrum of polyadenine unravelled by multiscale computations and quantitative wavefunction analysis

The characterization of the electronically excited states of DNA strands populated upon solar UV light absorption is essential to unveil light-induced DNA damage and repair processes. We report a comprehensive analysis of the electronic properties of the UV spectrum of single-stranded polyadenine based on theoretical calculations that include excitations over eight nucleobases of the DNA strand and environmental effects by a multiscale quantum mechanics/molecular mechanics scheme, conformational sampling by molecular dynamics, and a meaningful interpretation of the electronic structure by quantitative wavefunction analysis. We show that electronic excitations are extended mainly over two nucleobases with additional important contributions of monomer-like excitations and excitons delocalized over three monomers. Half of the spectral intensity derives from locally excited and Frenkel exciton states, while states with partial charge-transfer character account for the other half and pure charge-transfer states represent only a minor contribution. The hypochromism observed when going from the isolated monomer to the strand occurs independently from delocalization and charge transfer and is instead explained by long-range environmental perturbations of the monomer states.

Solvent effects on the excited state characteristics of adenine–thymine base pairs

A systematic analysis of the excited state characteristics of the DNA base pair adenine–thymine in stacked and Watson–Crick hydrogen bonded configurations has been carried out in this study.

Simulation of Electronic Circular Dichroism of Nucleic Acids: From the Structure to the Spectrum

We present a quantum mechanical (QM) simulation of the electronic circular dichroism (ECD) of nucleic acids (NAs). The simulation combines classical molecular dynamics, to obtain the structure and its temperature-dependent fluctuations, with a QM excitonic model to determine the ECD. The excitonic model takes into account environmental effects through a polarizable embedding and uses a refined approach to calculate the electronic couplings in terms of full transition densities. Three NAs with either similar conformations but different base sequences or similar base sequences but different conformations have been investigated and the results were compared with experimental observations; a good agreement was seen in all cases. A detailed analysis of the nature of the ECD bands in terms of their excitonic composition was also carried out. Finally, a comparison between the QM and the DeVoe models clearly revealed the importance of including fluctuations of the excitonic parameters and of accurately determining the electronic couplings. This study demonstrates the feasibility of the ab initio simulation of the ECD spectra of NAs, that is, without the need of experimental structural or electronic data.

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.

Computational modeling of photoexcitation in DNA single and double strands

The photoexcitation of DNA strands triggers extremely complex photoinduced processes, which cannot be understood solely on the basis of the behavior of the nucleobase building blocks. Decisive factors in DNA oligomers and polymers include collective electronic effects, excitonic coupling, hydrogen-bonding interactions, local steric hindrance, charge transfer, and environmental and solvent effects. This chapter surveys recent theoretical and computational efforts to model real-world excited-state DNA strands using a variety of established and emerging theoretical methods. One central issue is the role of localized vs delocalized excitations and the extent to which they determine the nature and the temporal evolution of the initial photoexcitation in DNA strands.

Influence of base stacking geometry on the nature of excited states in G-quadruplexes: a time-dependent DFT study

G-quadruplexes are four-stranded structures of nucleic acids that are formed from the association of guanine nucleobases into cyclical arrangements known as tetrads. G-quadruplexes are involved in a host of biological processes and are of interest in nanomaterial applications. However, not much is known about their electronic properties. In this paper, we analyze electronic excited states of G-quadruplexes using a combination of time-dependent DFT calculations and molecular dynamics simulations. We systematically consider experimentally observed arrangements of stacked guanine tetrads. The effects of structural features on exciton delocalization and photoinduced charge separation are explored using a quantitative analysis of the transition electron density. It is shown that collective coherent excitations shared between two guanine nucleobases dominate in the absorption spectrum of stacked G-tetrads. These excitations may also include a significant contribution of charge transfer states. Large variation in exciton localization is also observed between different structures with a general propensity toward localization between two bases. We reveal large differences in how charge separation occurs within different nucleobase arrangements, with some geometries favoring separation within a single tetrad and others favoring separation between tetrads. We also investigate the effects of the coordinating K(+) ion located in the central cavity of G-quadruplexes on the relative excited state properties of such systems. Our results demonstrate how the nature of excited states in G-quadruplexes depends on the nucleobase stacking geometry resulting from the mutual arrangement of guanine tetrads.

Electronic coupling between photo-excited stacked bases in DNA and RNA strands with emphasis on the bright states initially populated

In biology the interplay between multiple light-absorbers gives rise to complex quantum effects such as superposition states that are of extreme importance for life, both for harvesting solar energy and likely protecting nucleic acids from radiation damage. Still the characteristics of these states and their quantum dynamics are a much debated issue. While the electronic properties of single bases are fairly well understood, the situation for strands is complicated by the fact that stacked bases electronically couple when photoexcited. These newly arising states are denoted as exciton states and are simply linear combinations of localised wavefunctions that involve N - 1 ground-state bases and one base in its excited state (cf. the Frenkel exciton model). There is disagreement over the number of bases, N, that coherently couple, i.e., the spatial extent of the exciton, and how electronic deexcitation back to the ground state occurs. The importance of dark charge-transfer states has been inferred both from time-resolved fluorescence and transient absorption experiments. These states were suggested to be responsible for long deexcitation times but it is unclear whether 'long' is tens of picoseconds or nanoseconds. In this review paper, we focus on the bright states initially populated and discuss their nature based on information obtained from systematic absorption and circular dichroism experiments on single strands of different lengths. Our results from the last five years are compared with those from other groups, and are discussed in the context of successive deexcitation schemes. Pieces to the puzzle have come from different experiments and theory but a complete description has yet to emerge. As such the story about DNA/RNA photophysical decay mechanisms resembles the tale about the blind men and the elephant where all see the beast in different, correct but incomplete ways.

Dynamical treatment of charge transfer through duplex nucleic acids containing modified adenines

We address the issue of whether chemical alterations of nucleobases are an effective tool to modulate charge transfer through DNA molecules. Our investigation uses a multilevel computational approach based on classical molecular dynamics and quantum chemistry. We find yet another piece of evidence that structural fluctuations are a key factor to determine the electronic structure of double-stranded DNA. We argue that the electronic structure and charge transfer ability of flexible polymers is the result of a complex intertwining of various structural, dynamical and chemical factors. Chemical intuition may be used to design molecular wires, but this is not the sole component in the complex charge transfer mechanism through DNA.

Solvent dependency of the UV-Vis spectrum of indenoisoquinolines: role of keto-oxygens as polarity interaction probes

Indenoisoquinolines are the most promising non-campthotecins topoisomerase IB inhibitors. We present an integrated experimental/computational investigation of the UV-Vis spectra of the IQNs parental compound (NSC314622) and two of its derivatives (NSC724998 and NSC725776) currently undergoing Phase I clinical trials. In all the three compounds a similar dependence of the relative absorption intensities at 270 nm and 290 nm on solvent polarity is found. The keto-oxygens in positions 5 and 11 of the molecular scaffold of the molecule are the principal chromophores involved in this dependence. Protic interactions on these sites are also found to give rise to absorptions at wavelength <250 nm observed in water solution, due to the stabilization of highly polarized tautomers of the molecule. These results suggest that the keto-oxygens are important polarizable groups that can act as useful interactors with the molecular receptor, providing at the same time an useful fingerprint for the monitoring of the drug binding to topoisomerase IB.

On the nature of DNA hyperchromic effect

A combined theoretical-experimental study of the hyperchromic effect as occurring in the denaturation of a double stranded polyA-polyT is presented. Our theoretical/computational procedure allows us to reproduce the essential features of the experimental spectra and to characterize those molecular interactions responsible for the changes in the UV absorbance. We found that although excitonic intrastrand interactions strongly affect the absorbance, they are almost fully maintained in the single-stranded DNA. Our data indicate that hyperchromic effect originates from the higher delocalization of the excitonic states in the denaturated DNA with respect to the double-stranded conformation.

Coarse-grained model DNA: structure, sequences, stems, circles, hairpins

A coarse-grained model for DNA that is intended to function realistically at the level of individual bases is reported. The model is composed of residues with up to eight coarse-grained beads each, which is sufficient for DNA-like base stacking and base-base recognition by hydrogen bonding. The beads interact by means of short-ranged pair potentials and a simple implicit solvent model. Movement is simulated by Brownian dynamics without hydrodynamic coupling. The main stabilizing forces are base stacking and hydrogen bonding, as modified by the effects of solvation. Complementary double-stranded chains of such residues form stable double helices over long runs (~10 μs) at or near room temperature, with structural parameters close to those of B-form DNA. Most mismatched chains or mismatched regions within a complementary molecule melt and become disordered. Long-range fluctuations and elastic properties, as measured by bending and twisting persistence lengths, are close to experimental values. Single-stranded chains are flexible, with transient stretches of free bases in equilibrium with globules stabilized by intrastrand stacking and hydrogen bonding. Model DNAs in covalently closed loops form right- or left-handed supercoils, depending on the sign of overtwist or undertwist. Short stem-loop structures melt at elevated temperatures and reanneal when the temperature is carefully lowered. Overall, most qualitative properties of real DNA arise naturally in the model from local interactions at the base-pair level.

Density functional theory studies of the extent of hole delocalization in one-electron oxidized adenine and guanine base stacks

This study investigates the extent of hole delocalization in one-electron oxidized adenine (A) and guanine (G) stacks and shows that new IR vibrational bands are predicted that are characteristic of hole delocalization within A-stacks. The geometries of A-stacks (A(i); i = 2-8) and G-stacks (GG and GGG) in their neutral and one-electron oxidized states were optimized with the bases in a B-DNA conformation using the M06-2X/6-31G* method. The highest occupied molecular orbital (HOMO) is localized on a single adenine in A-stacks and on a single guanine in GG and GGG stacks located at the 5'-site of the stack. On one-electron oxidation (removal of an electron from the HOMO of the neutral A- and G-stacks) a "hole" is created. Mulliken charge analysis shows that these "holes" are delocalized over two to three adenine bases in the A-stack. The calculated spin density distribution of A(i)(•+) (i = 2-8) also showed delocalization of the hole predominantly on two adenine bases, with some delocalization on a neighboring base. For GG and GGG radical cations, the hole was found to be localized on a single G in the stack. The calculated HFCCs of GG and GGG are in good agreement with the experiment. Further, from the vibrational frequency analysis, it was found that IR spectra of neutral and the corresponding one-electron oxidized adenine stacks are quite different. The IR spectra of A(2)(•+) has intense IR peaks between 900 and 1500 cm(-1) that are not present in the neutral A(2) stack. The presence of A(2)(•+) in the adenine stack has a characteristic intense peak at ~1100 cm(-1). Thus, IR and Raman spectroscopy has potential for monitoring the extent of hole delocalization in A stacks.

Photodynamics of the adenine model 4-aminopyrimidine embedded within double strand of DNA

On-the-fly surface hopping nonadiabatic photodynamical simulations using hybrid quantum mechanical/molecular mechanical approach of 4-aminopyrimidine were performed to model the relaxation mechanism of adenine within DNA double strand. The surrounding bases do not affect the overall ring-puckering relaxation mechanisms significantly, however, interesting hydrogen-bond dynamics is observed. First, formation of intra-strand hydrogen bonds is found. It is shown that this effect speeds up the decay process. In addition, the Watson–Crick structure is altered by breaking one of the inter-strand hydrogen bonds also leading to a decrease of the life time.

Effect of oligonucleotide length on the assembly of DNA materials: molecular dynamics simulations of layer-by-layer DNA films

DNA strand length has been found to be an important factor in many DNA-based nanoscale systems. Here, we apply molecular dynamics simulations in a synergistic effort with layer-by-layer experimental data to understand the effect of DNA strand length on the assembly of DNA films. The results indicate that short (less than 10 bases) and long (more than 30 bases) single-stranded DNAs do not exhibit optimal film growth, and this can be associated with the limited accessibility of the bases on the surface due to formation of self-protected interactions that prevent efficient hybridization. Interestingly, the presence of a duplex attached to a single strand significantly alters the persistence length of the polyT strands. Our study suggests that restrained polyT, compared to labile suspensions of free polyT, are more capable of hybridization and hence DNA-based assembly.

Vacuum-ultraviolet circular dichroism spectroscopy of DNA: a valuable tool to elucidate topology and electronic coupling in DNA

Circular dichroism (CD) is a powerful technique to obtain information on electronic transitions and has been used extensively for studies on DNA. Most experiments are done in the UV region but new information is often revealed from extending the wavelength region down into the vacuum ultraviolet (VUV) region. Such experiments are most easily carried out with synchrotron radiation (SR) light sources that provide large photon fluxes. Here we provide a summary of the SRCD data taken on different DNA strands with emphasis on results from our own laboratory within the last five years.(1-3) Signal intensities in the VUV are often significantly larger than those in the UV, and the electronic coupling between bases may increase with excitation energy. CD spectroscopy is particularly useful for investigating the extent of electronic coupling within a strand, i.e., the degree of delocalisation of the excited-state electronic wavefunction. The spatial extent of the wavefunction may be limited to just one base or it extends over two or more bases in a stack or between bases on different strands.(4,5) The actual character of the electronically excited state is linked to base composition and sequence as well as DNA folding motif (A-, B-, Z-DNA, triplexes, quadruplexes, etc.). The latter depends on experimental conditions such as solution acidity, temperature, ionic strength, and solvent.

Vibronic model for the quantum dynamical study of the competition between bright and charge-transfer excited states in single-strand polynucleotides: the adenine dimer case

A simple vibronic model aimed at investigating the interplay between bright excitonic states and dark charge-transfer (CT) states in stacked adenine (Ade) nucleobases is presented. Two orbitals (the HOMO and the LUMO) for each Ade site have been included in the electronic Hamiltonian, whose parameters have been fitted to reproduce the main features of the absorption spectra of two stacked 9-methyladenine (9Me-A) molecules, computed in aqueous solution at the PCM/TD-PBE0 level. Three modes for each adenine unit have been included in the Hamiltonian, to describe the main structural changes among the different excited state minima of the adenine stacked dimer, as described at the TD-DFT level. The developed vibronic Hamiltonian (four electronic states and six nuclear coordinates) has been adopted to perform quantum dynamical calculations of a photoexcited Ade stacked dimer, utilizing the multiconfigurational time-dependent Hartree method. The obtained results indicate that the transfer between the bright excitonic state and the CT state is fast and effective.

Fingerprints of bonding motifs in DNA duplexes of adenine and thymine revealed from circular dichroism: synchrotron radiation experiments and TDDFT calculations

Synchrotron radiation circular dichroism (SRCD) spectra were recorded for a family of 12 DNA duplexes that all contain nine adenines (A) and nine thymines (T) in each strand but in different combinations. The total number of AT Watson-Crick (WC) base pairs is constant (18), but the number of cross-strand (CS) hydrogen bonds between A and T varies between 0 and 16, the maximum possible. Eleven of the duplexes have one or more A tracts, and one duplex has T tracts. The signals due to hybridization were found from subtraction of spectra of single strands from spectra of the duplexes. The residual spectrum of the T-tract duplex T(9)A(9):A(9)T(9) (5'-3':3'-5') significantly differs from that of the A-tract duplex A(9)T(9):T(9)A(9), but only below 210 nm, which suggests that the signal in this region depends on the superhelicity of the duplex. A principal component analysis of all residual spectra reveals that spectra of A-tract duplexes can be obtained to a good approximation as a linear combination of just two basis spectra. The first component is assigned to the spectrum of 18 WC and 8 CS pairs, whereas the second component is that of 8 CS pairs. This interpretation is supported by separate experiments on duplexes of varying lengths but with similar arrangements of the A and T's and by experiments on two other duplex families of 14 and 30 base pairs. The best correlation is obtained by the assumption that cross-strand interactions occur as long as there are two adenine neighbors in a strand. Our data indicate that a circular dichroism spectrum of a duplex containing only A and T can simply be inferred from the number of WC base pairs and the number of CS interactions, and we provide reference spectra for these two interactions. Finally, time dependent density functional theory calculations of the circular dichroism spectra for an isolated WC base pair and two different CS base pairs (between adenine N-6 amine and thymine O-4 or between adenine C-2-H and thymine O-2) were performed to provide some additional support for the interpretation of the experimental spectra. We find large differences between the two calculated CS spectra. However, there is a reasonable qualitative agreement between the calculated WC and the C-2-H...O-2 CS spectra and those deduced from the experimental data.