Role of excitonic coupling and charge-transfer states in the absorption and CD spectra of adenine-based oligonucleotides investigated through QM/MM simulations

Factors Affecting the Population of Excited Charge Transfer States in Adenine/Guanine Dinucleotides: A Joint Computational and Transient Absorption Study

There is compelling evidence that the absorption of low-energy UV radiation directly by DNA in solution generates guanine radicals with quantum yields that are strongly dependent on the secondary structure. Key players in this unexpected phenomenon are the photo-induced charge transfer (CT) states, in which an electric charge has been transferred from one nucleobase to another. The present work examines the factors affecting the population of these states during electronic relaxation. It focuses on two dinucleotides with opposite orientation: 5'-dApdG-3' (AG) and 5'-dGpdA-3' (GA). Quantum chemistry calculations determine their ground state geometry and the associated Franck-Condon states, map their relaxation pathways leading to excited state minima, and compute their absorption spectra. It has been shown that the most stable conformer is anti-syn for AG and anti-anti for GA. The ground state geometry governs both the excited states populated upon UV photon absorption and the type of excited state minima reached during their relaxation. Their fingerprints are detected in the transient absorption spectra recorded with excitation at 266 nm and a time resolution of 30 fs. Our measurements reveal that in the large majority of dinucleotides, chromophore coupling is already operative in the ground state and that the charge transfer process occurs within ~120 fs. The competition among various relaxation pathways affects the quantum yields of the CT state formation in each dinucleotide, which are estimated to be 0.18 and 0.32 for AG and GA, respectively.

On the Use of the Intrinsic DNA Fluorescence for Monitoring Its Damage: A Contribution from Fundamental Studies

The assessment of DNA damage by means of appropriate fluorescent probes is widely spread. In the specific case of UV-induced damage, it has been suggested to use the emission of dimeric photoproducts as an internal indicator for the efficacy of spermicidal lamps. However, in the light of fundamental studies on the UV-induced processes, outlined in this review, this is not straightforward. It is by now well established that, in addition to photodimers formed via an electronic excited state, photoionization also takes place with comparable or higher quantum yields, depending on the irradiation wavelength. Among the multitude of final lesions, some have been fully characterized, but others remain unknown; some of them may emit, while others go undetected upon monitoring fluorescence, the result being strongly dependent on both the irradiation and the excitation wavelength. In contrast, the fluorescence of undamaged nucleobases associated with emission from ππ* states, localized or excitonic, appearing at wavelengths shorter than 330 nm is worthy of being explored to this end. Despite its low quantum yield, it is readily detected nowadays. Its intensity decreases due to the disappearance of the reacting nucleobases and the loss of exciton coherence provoked by the presence of lesions, independently of their type. Thus, it could potentially provide valuable information about the DNA damage induced, not only by UV radiation but also by other sanitizing or therapeutic agents.

Multistate Energy Decomposition Analysis of Molecular Excited States

A multistate energy decomposition analysis (MS-EDA) method is described to dissect the energy components in molecular complexes in excited states. In MS-EDA, the total binding energy of an excimer or an exciplex is partitioned into a ground-state term, called local interaction energy, and excited-state contributions that include exciton excitation energy, superexchange stabilization, and orbital and configuration-state delocalization. An important feature of MS-EDA is that key intermediate states associated with different energy terms can be variationally optimized, providing quantitative insights into widely used physical concepts such as exciton delocalization and superexchange charge-transfer effects in excited states. By introducing structure-weighted adiabatic excitation energy as the minimum photoexcitation energy needed to produce an excited-state complex, the binding energy of an exciplex and excimer can be defined. On the basis of the nature of intermolecular forces through MS-EDA analysis, it was found that molecular complexes in the excited states can be classified into three main categories, including (1) encounter excited-state complex, (2) charge-transfer exciplex, and (3) intimate excimer or exciplex. The illustrative examples in this Perspective highlight the interplay of local excitation polarization, exciton resonance, and superexchange effects in molecular excited states. It is hoped that MS-EDA can be a useful tool for understanding photochemical and photobiological processes.

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.

High-Energy Long-Lived Emitting Mixed Excitons in Homopolymeric Adenine-Thymine DNA Duplexes

The publication deals with polymeric pA●pT and oligomeric A₂₀●T₂₀ DNA duplexes whose fluorescence is studied by time-correlated single photon counting. It is shown that their emission on the nanosecond timescale is largely dominated by high-energy components peaking at a wavelength shorter than 305 nm. Because of their anisotropy (0.02) and their sensitivity to base stacking, modulated by the duplex size and the ionic strength of the solution, these components are attributed to mixed ππ*/charge transfer excitons. As high-energy long-lived excited states may be responsible for photochemical reactions, their identification via theoretical studies is an important challenge.

Specific chiroptical sensing of cysteine via ultrasound-assisted formation of disulfide bonds in aqueous solution

Cysteine (Cys) can serve as a biomarker to indicate diseases or disorders, and its chiral sensing has attracted increasing attention. Herein, we established an ultrasound-facilitated chiral sensing method for Cys using 4-chloro-7-nitro-1,2,3-benzoxadiazole (NBD-Cl) and electronic circular dichroism (ECD) spectroscopy. The formation of chiral disulfide bonds induced degenerate exciton coupling between two NBD chromophores, resulting in intense Cotton effects (CEs) of the sensing product. The anisotropy factor (g) was linearly correlated with the enantiomeric excess of Cys across the visible region (400-500 nm), and other natural amino acids or biothiols did not interfere with the detection. This ultrasound-promoted efficient and specific chiral sensing method of Cys has potential for application in the diagnosis of related diseases.

Dual Time-Scale Proton Transfer and High-Energy, Long-Lived Excitons Unveiled by Broadband Ultrafast Time-Resolved Fluorescence in Adenine-Uracil RNA Duplexes

In contrast to the immense amount of research on electronically excited DNA, surprisingly little has been done about the excited states of RNA. Herein, we demonstrate an ultrafast broadband time-resolved fluorescence and fluorescence anisotropy study to probe directly the intrinsic fluorescence and overall dynamics of the fluorescence from a homopolymeric adenine·uracil RNA duplex adopting the A-form structure. The results unveiled complex deactivation through distinctive multichannels mediated by states of varied energy, a character of charge transfer, and a lifetime from sub-picosecond to nanoseconds. In particular, we observed an unprecedented kinetic isotopic effect and participation of unusual proton transfer from states in two discrete energies and time domains. We also identified a high-energy nanosecond emission that we attributed to its fluorescence anisotropy to long-lived weakly emissive excitons not reported in DNA. These distinguishing features originate from the stacking, pairing, and local hydration environment specific to the A-form conformation of the adenine·uracil double helix.

ECD exciton chirality method today: a modern tool for determining absolute configurations

The application of the exciton chirality method (ECM) to interpret electronic circular dichroism (ECD) spectra is a well‐established and still popular approach to assign the absolute configuration (AC) of natural products, chiral organic compounds, and organometallic species. The method applies to compounds containing at least two chromophores with electric dipole allowed transitions (e.g., π–π* transitions). The exciton chirality rule correlates the sign of an exciton couplet (two ECD bands with opposite sign and similar intensity) with the overall molecular stereochemistry, including the AC. A correct application of the ECM requires three main prerequisites: (a) the knowledge of the molecular conformation, (b) the knowledge of the directions of the electric transition moments (TDMs), and (c) the assumption that the exciton coupling mechanism must be the major source of the observed ECD signals. All these prerequisites can be easily verified by means of quantum‐mechanical (QM) calculations. In the present review, we shortly introduce the general principles that underpin the use of the ECM for configurational assignments and survey its applications, both classic ones and some reported in the recent literature. Based on these examples, we will stress the advantages of the ECM but also the key requisites for its correct application. Additionally, we will discuss the dependence of the couplet sign on geometrical parameters (angles α,β,γ between TDMs), which can be helpful for discerning the sign of exciton chirality in ambiguous situations. Finally, we will present a molecular orbital (MO) description of the exciton coupling phenomenon.

Fundamentals of the Intrinsic DNA Fluorescence

The intrinsic fluorescence of nucleic acids is extremely weak compared to that of the fluorescent labels used to probe their structural and functional behavior. Thus, for technical reasons, the investigation of the intrinsic DNA fluorescence was limited for a long time. But with the improvement in spectroscopic techniques, the situation started to change around the turn of the century. During the past two decades, various factors modulating the static and dynamic properties of the DNA fluorescence have been determined; it was shown that, under certain conditions, quantum yields may be up 100 times higher than what was known so far. The ensemble of these studies opened up new paths for the development of label-free DNA fluorescence for biochemical applications. In parallel, these studies have shed new light on the primary processes leading to photoreactions that damage DNA when it absorbs UV radiation.We have been studying a variety of DNA systems, ranging from the monomeric nucleobases to double-stranded and four-stranded structures using fluorescence spectroscopy. The specificity of our work resides in the quantitative association of the steady-state fluorescence spectra with time-resolved data recorded from the femtosecond to the nanosecond timescales, made possible by the development of specific methodologies.Among others, our fluorescence studies provide information on the energy and the polarization of electronic transitions. These are valuable indicators for the evolution of electronic excitations in complex systems, where the electronic coupling between chromophores plays a key role. Highlighting collective effects that originate from electronic interactions in DNA multimers is the objective of the present Account.In contrast to the monomeric chromophores, whose fluorescence decays within a few picoseconds, that of DNA multimers persists on the nanosecond timescale. Even if long-lived states represent only a small fraction of electronic excitations, they may be crucial to the DNA photoreactivity because the probability to reach reactive conformations increases over time, owing to the incessant structural dynamics of nucleic acids.Our femtosecond studies have revealed that an ultrafast excitation energy transfer takes place among the nucleobases within duplexes and G-quadruplexes. Such an ultrafast process is possible when collective states are populated directly upon photon absorption. At much longer times, we discovered an unexpected long-lived high-energy emission stemming from what was coined "HELM excitons". These collective states, whose emission increases with the duplex size, could be responsible for the delayed fluorescence of ππ* states observed for genomic DNA.Most studies dealing with excited-state relaxation in DNA were carried out with excitation in the absorption band peaking at around 260 nm. We went beyond this and also performed the first time-resolved study with excitation in the UVA spectral range, where a very weak absorption tail is present. The resulting fluorescence decays are much slower and the fluorescence quantum yields are much higher than for UVC excitation. We showed that the base pairing of DNA strands enhances the UVA fluorescence and, in parallel, increases the photoreactivity because it modifies the nature of the involved collective excited states.

Benchmarking Quantum Mechanical Methods for the Description of Charge-Transfer States in π-Stacked Nucleobases

Charge-transfer (CT) states are of special interest in photochemical research because they can facilitate chemical reactions through the rearrangement of electrons and subsequently chemical bonds in a molecular system. Of particular importance to this research is the transfer of electrons between π-stacked nucleobases in DNA because they play an important role in its photophysics and photochemistry. Computational methods are paramount for the study of CT states because of the current inability of experimental methods to easily detect such states. However, many ab-initio wavefunction-based and density functional theory (DFT) methods fail to accurately describe these CT states. Here, we benchmark how 40 different quantum mechanical methods describe the excited states of a guanine-thymine π-stacked nucleobase dimer system, both in 5'-TG-3' and 5'-GT-3' conformations. We find that the distance between the nucleobases plays a major role in the energy of the CT state and in the difference of the dipole moments between the CT and ground state. There is a larger range of values (and errors) for the energies of CT states compared to those of states localized on one nucleobase. Wavefunction-based methods have similar errors for the CT and localized valence states, while DFT methods are very sensitive to the amount of Hartree-Fock exchange. Long-range-corrected functionals with a careful balance of the Hartree-Fock exchange included can predict very accurate CT states and a balanced description with the localized states.

Guanine Radicals Induced in DNA by Low-Energy Photoionization

Guanine (G) radicals are precursors to DNA oxidative damage, correlated with carcinogenesis and aging. During the past few years, we demonstrated clearly an intriguing effect: G radicals can be generated upon direct absorption of UV radiation with energy significantly lower than the G ionization potential. Using nanosecond transient absorption spectroscopy, we studied the primary species, ejected electrons and guanine radicals, which result from photoionization of various DNA systems in aqueous solution.The DNA propensity to undergo electron detachment at low photon energies greatly depends on its secondary structure. Undetected for monomers or unstacked oligomers, this propensity may be 1 order of magnitude higher for G-quadruplexes than for duplexes. The experimental results suggest nonvertical processes, associated with the relaxation of electronic excited states. Theoretical studies are required to validate the mechanism and determine the factors that come into play. Such a mechanism, which may be operative over a broad excitation wavelength range, explains the occurrence of oxidative damage observed upon UVB and UVA irradiation.Quantification of G radical populations and their time evolution questions some widespread views. It appears that G radicals may be generated with the same probability as pyrimidine dimers, which are considered to be the major lesions induced upon absorption of low-energy UV radiation by DNA. As most radical cations undergo deprotonation, the vast majority of the final reaction products is expected to stem from long-lived deprotonated radicals. Consequently, when G radical cations are involved, the widely used oxidation marker 8-oxodG is not representative of the oxidative damage.Beyond the biological consequences, photogeneration of electron holes in G-quadruplexes may inspire applications in nanoelectronics; although four-stranded structures are currently studied as molecular wires, their behavior as photoconductors has not been explored so far.In the present Account, after highlighting some key experimental issues, we first describe the photoionization process, and then, we focus on radicals. We use as show-cases new results obtained for genomic DNA and Oxytricha G-quadruplexes. Generation and reaction dynamics of G radicals in these systems provide a representative picture of the phenomena reported previously for duplexes and G-quadruplexes, respectively.

Chirality of a rhodamine heterodimer linked to a DNA scaffold: an experimental and computational study

The chiroptical properties of multi-chromophoric systems are governed by the intermolecular arrangement of the monomeric units. We report on a computational and experimental study of the linear optical properties and supramolecular structure of a rhodamine heterodimer assembled on a DNA scaffold. The experimental absorption and circular dichroism (CD) profiles confirm the dimer formation. Computationally, starting from low-cost DFT/TDDFT simulations of the bare dimer we attribute the measured -/+ CD sign sequence of the S1/S2 bands to a specific chiral conformation of the heterodimer. In the monomers, as typical for rhodamine dyes, the electric transition dipole of the lowest π-π* transition is parallel to the long axis of the xanthene planes. We show that in the heterodimer the sign sequence of the two CD bands is related to the orientation of these long axes. To account explicitly for environment effects, we use molecular dynamics (MD) simulations for characterizing the supramolecular structure of the two optical isomers tethered on DNA. Average absorption and CD-profiles were modeled using ab initio TDDFT calculations at the geometries sampled along a few nanosecond MD run. The absorption profiles computed for both optical isomers are in good agreement with the experimental absorption spectrum and do not allow one to discriminate between them. The computed averaged CD profiles provide the orientation of monomers in the enantiomer that is dominant under the experimental conditions.

Azobenzene as a photoregulator covalently attached to RNA: a quantum mechanics/molecular mechanics-surface hopping dynamics study

Azobenzene covalently attached to RNA undergoes trans-to-cis photo-switching on a time scale of ∼15 picoseconds – 30 times slower than in vacuo.

The photoregulation of nucleic acids by azobenzene photoswitches has recently attracted considerable interest in the context of emerging biotechnological applications. To understand the mechanism of photoinduced isomerisation and conformational control in these complex biological environments, we employ a Quantum Mechanics/Molecular Mechanics (QM/MM) approach in conjunction with nonadiabatic Surface Hopping (SH) dynamics. Two representative RNA–azobenzene complexes are investigated, both of which contain the azobenzene chromophore covalently attached to an RNA double strand via a β-deoxyribose linker. Due to the pronounced constraints of the local RNA environment, it is found that trans-to-cis isomerization is slowed down to a time scale of ∼10–15 picoseconds, in contrast to 500 femtoseconds in vacuo, with a quantum yield reduced by a factor of two. By contrast, cis-to-trans isomerization remains in a sub-picosecond regime. A volume-conserving isomerization mechanism is found, similarly to the pedal-like mechanism previously identified for azobenzene in solution phase. Strikingly, the chiral RNA environment induces opposite right-handed and left-handed helicities of the ground-state cis-azobenzene chromophore in the two RNA–azobenzene complexes, along with an almost completely chirality conserving photochemical pathway for these helical enantiomers.

Sequential electron transfer governs the UV-induced self-repair of DNA photolesions

Cyclobutane pyrimidine dimers (CpDs) are among the most common DNA lesions occurring due to the interaction with ultraviolet light. While photolyases have been well known as external factors repairing CpDs, the intrinsic self-repairing capabilities of the GAT 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T DNA sequence were discovered only recently and are still largely obscure. Here, we elucidate the mechanistic details of this self-repair process by means of MD simulations and QM/MM computations involving the algebraic diagrammatic construction to the second order [ADC(2)] method. We show that local UV-excitation of guanine may be followed by up to three subsequent electron transfers, which may eventually enable efficient CpD ring opening when the negative charge resides on the T 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T dimer. Consequently, the molecular mechanism of GAT 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T self-repair can be envisaged as sequential electron transfer (SET) occurring downhill along the slope of the S₁ potential energy surface. Even though the general features of the SET mechanism are retained in both of the studied stacked conformers, our optimizations of different S₁/S₀ state crossings revealed minor differences which could influence their self-repair efficiencies. We expect that such assessment of the availability and efficiency of the SET process in other DNA oligomers could hint towards other sequences exhibiting similar photochemical properties. Such explorations will be particularly fascinating in the context of the origins of biomolecules on Earth, owing to the lack of external repairing factors in the Archean age.

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.

Photophysical deactivation pathways in adenine oligonucleotides

In this work we study deactivation processes in adenine oligomers after absorption of UV radiation using Quantum Mechanics combined with Molecular Mechanics (QM/MM). Correlated electronic structure methods appropriate for describing the excited states are used to describe a π-stacked dimer of adenine bases incorporated into (dA)20(dT)20. The results of these calculations reveal three different types of excited state minima which play a role in deactivation processes. Within this set of minima there are minima where the excited state is localized on one adenine (monomer-like) as well as minima where the excited state is delocalized on two adenines, forming different types of excimers and bonded excimers of varying but inter-related character. The proximity of their energies reveals that the minima can decay into one another along a flat potential energy surface dependent on the interbase separation. Additionally, analysis of the emissive energies and other physical properties, including theoretical anisotropy calculations, and comparison with fluorescence experiments, provides evidence that excimers play an important role in long-lived signals in adenine oligonucleotides while the subpicosecond decay is attributed to monomer-like minima. The necessity for a close approach of the nucleobases reveals that the deactivation mechanism is tied to macro-molecular motion.

Characterization of the excited states of DNA building blocks: a coupled cluster computational study

DNA building blocks consisting of up to four nucleobases are investigated using the EOM-CCSD and CC2-LR methods in two B-DNA-like arrangements of a poly-adenine:poly-thymine (poly-A:poly-T) system. Excitation energies and oscillator strengths are presented and the characteristics of the excited states are discussed. Excited states of single-stranded poly-A systems are highly delocalized, especially the spectroscopically bright states, where delocalization over up to four fragments can be observed. In the case of poly-T systems, the states are somewhat less delocalized, extending to maximally about three fragments. A single A:T Watson-Crick pair has highly localized states, while delocalization over base pairs can be observed for some excited states of the (A)2:(T)2 system, but intrastrand delocalization is more pronounced in this case, as well. As for the characteristics of the simulated UV absorption spectra, a significant decrease of intensity can be observed in the case of single strands with increasing chain length; this is due to the stacking interactions and is in accordance with previous results. On the other hand, the breaking of H-bonds between the two strands does not alter the spectral intensity considerably, it only causes a redshift of the absorption band, thus it is unable to explain the experimentally observed DNA hyperchromism on its own, and stacking interactions need to be considered for the description of this effect as well.

Excimers and Exciplexes in Photoinitiated Processes of Oligonucleotides

A lot has been learned about the physical and chemical transformations that originate from the absorption of light by DNA, and computational chemistry has played a critical role in revealing the mechanisms of how these transformations occur. Nucleic acids consist of chromophores interacting via π stacking and hydrogen bonding. The fate of these systems after they absorb light is determined by the interplay and competition between pathways involving one chromophore or interacting chromophores. This Perspective highlights the role of π stacking in photophysical and photochemical processes in oligonucleotides and reveals the importance of excimers and exciplexes. Special types of excimers/exciplexes, characterized as bonded excimers/exciplexes, are also found to be important.

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.

QM/MM studies reveal pathways leading to the quenching of the formation of thymine dimer photoproduct by flanking bases

It is known that the formation of the photochemical product of thymine-thymine cyclobutane pyrimidine dimer (TT-CPD) formed upon UV excitation in DNA is significantly affected by the nature of the flanking bases, and that the oxidation potential of the flanking base correlates with the quenching of TT-CPD formation. However, the electronic details of this correlation have remained controversial. The quenching of thymine dimer formation exerted by flanking bases was suggested to be driven by both conformational and electronic effects. In the present study, we examine both of these effects using umbrella sampling and a quantum mechanical/molecular mechanical (QM/MM) approach for selected model systems. Our results demonstrate that a charge transfer (CT) state between the flanking base and the adjacent thymine base can provide a decay pathway for the population to escape from dimer formation, which eventually leads to the formation of an exciplex. The QM/MM vertical excitation energies also reveal that the oxidation potential of flanking bases correlates with the energy level of the CT state, thereby determining whether the CT state intersects with the state that can lead to dimer formation. The consistency between these results and experimentally obtained dimer formation rates implies that the quenching of dimer formation is mainly attributed to the decay pathway via the CT state. The present results further underline the importance of the electronic effects in quenching.

UV-induced DNA Damage: The Role of Electronic Excited States

The knowledge of the fundamental processes induced by the direct absorption of UV radiation by DNA allows extrapolating conclusions drawn from in vitro studies to the in-vivo DNA photoreactivity. In this respect, the characterization of the DNA electronic excited states plays a key role. For a long time, the mechanisms of DNA lesion formation were discussed in terms of generic "singlet" and "triplet" excited state reactivity. However, since the beginning of the 21(st) century, both experimental and theoretical studies revealed the existence of "collective" excited states, i.e. excited states delocalized over at least two bases. Two limiting cases are distinguished: Frenkel excitons (delocalized ππ* states) and charge-transfer states in which positive and negative charges are located on different bases. The importance of collective excited states in photon absorption (in particular in the UVA spectral domain), the redistribution of the excitation energy within DNA, and the formation of dimeric pyrimidine photoproducts is discussed. The dependence of the behavior of the collective excited states on conformational motions of the nucleic acids is highlighted.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

Probing deactivation pathways of DNA nucleobases by two-dimensional electronic spectroscopy: first principles simulations

The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.

Excited state evolution of DNA stacked adenines resolved at the CASPT2//CASSCF/Amber level: from the bright to the excimer state and back

L_a and excimer state population exchange, along the common puckering decay coordinate, explains the longest DNA lifetime component.