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

Effect of A-DNA and B-DNA Conformation on the Interplay between Local Excitations and Charge-Transfer States in the Ultrafast Decay of Guanine-Cytosine Stacked Dimers: A Quantum Dynamical Investigation

We here simulate in the gas phase the population dynamics of guanine/cytosine (GC) and cytosine/guanine (CG) stacked dimers in B-DNA and A-DNA arrangement, following excitation in the lowest-energy band, and considering the four lowest-energy ππ* bright excited states, the three lowest-energy nπ* states, and the G → C charge-transfer (CT) state. We resort to a generalized Linear Vibronic Coupling (LVC) model parametrized with time-dependent density functional theory (TD-DFT) computations, exploiting a fragment-based diabatization and we run nonadiabatic quantum dynamical simulations with the multilayer version of the Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. G → C CT results in a major decay process for GC in B-DNA but less in A-DNA arrangement, where also the population transfer to the lowest-energy excited state localized on C is an important intermonomer process. In CG arrangements, mostly intramonomeric decays take place. We simulate the dynamics of several other GC structures whose arrangement is intermediate between B-DNA and A-DNA, obtaining further insights on the effect that the sequence and, especially, the stacking geometry have on the population transfer to the G → C CT.

Simulations of photoinduced processes with the exact factorization: state of the art and perspectives

This perspective offers an overview of the applications of the exact factorization of the electron-nuclear wavefunction to the domain of theoretical photochemistry, where the aim is to gain insights into the ultrafast dynamics of molecular systems via simulations of their excited-state dynamics beyond the Born-Oppenheimer approximation. The exact factorization offers an alternative viewpoint to the Born-Huang representation for the interpretation of dynamical processes involving the electronic ground and excited states as well as their coupling through the nuclear motion. Therefore, the formalism has been used to derive algorithms for quantum molecular-dynamics simulations where the nuclear motion is treated using trajectories and the electrons are treated quantum mechanically. These algorithms have the characteristic features of being based on coupled and on auxiliary trajectories, and have shown excellent performance in describing a variety of excited-state processes, as this perspective illustrates. We conclude with a discussion on the authors' point of view on the future of the exact factorization.

Excitonic Approach for Nonadiabatic Dynamics: Extending Beyond the Frenkel Exciton Model

We report the formulation and implementation of an extended Frenkel exciton model (EFEM) designed for simulating the dynamics of multichromophoric systems, taking into account the possible presence of interchromophore charge transfer states, as well as other states in which two chromophores are simultaneously excited. Our approach involves constructing a Hamiltonian based on calculations performed on monomers and selected dimers within the multichromophoric aggregate. Nonadiabatic molecular dynamics is addressed using a surface hopping approach, while the electronic wave functions and energies required for constructing the EFEM are computed utilizing the semiempirical floating occupation molecular orbitals-configuration interaction (FOMO-CI) electronic structure method. To validate our approach, we simulate the singlet fission process in a trimer of 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene (ThBF) molecules, embedded in their crystal environment, comparing the results of the EFEM to the standard "supermolecule" approach.

Simulation of the non-adiabatic dynamics of an enone-Lewis acid complex in an explicit solvent

Unlocking the full potential of Lewis acid catalysis for photochemical transformations requires a comprehensive understanding of the ultrafast dynamics of substrate-Lewis acid complexes. In a previous article [Peschel et al., Angew. Chem. Int. Ed., 2021, 60, 10155], time-resolved spectroscopy supported by static calculations revealed that the Lewis acid remains attached during the relaxation of the model complex cyclohexenone-BF3. In contrast to the experimental observation, surface-hopping dynamics in the gas phase predicted ultrafast heterolytic dissociation. We attributed the discrepancy to missing solvent interactions. Thus, in this work, we present an interface between the SHARC and FermiONs++ program packages, which enables us to investigate the ultrafast dynamics of cyclohexenone-BF3 in an explicit solvent environment. Our simulations demonstrate that the solvent prevents the dissociation of the complex, leading to an intriguing dissociation-reassociation mechanism. Comparing the dynamics with and without triplet states highlights their role in the relaxation process and shows that the Lewis acid inhibits intersystem crossing. These findings provide a clear picture of the relaxation process, which may aid in designing future Lewis acid catalysts for photochemical applications. They underscore that an explicit solvent model is required to describe relaxation processes in weakly bound states, as energy transfer to the solvent is crucial for the system to reach its minimum geometries.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain onπ π * state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained byn π * trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of theπ π * state leading to a wrongπ π * / n π * energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Unraveling Excited State Dynamics of a Single-Stranded DNA-Assembled Conjugated Polyelectrolyte

Conformational templating of conjugated polyelectrolytes with single-stranded DNAs (ssDNAs) has the prospect of tailoring excited state dynamics for specific optoelectronic applications. We use ultrafast time-resolved infrared spectroscopy to study the photophysics of a cationic polythiophene assembled with different ssDNAs, inducing distinct conformations (flexible disordered structures vs more rigid complexes with increased backbone planarity). Intrachain polarons are always produced upon selective excitation of the polymer, the extent being dependent on backbone torsional order. Polaron formation and decay were monitored through evolution of IR-active vibrational modes that interfere with mid-IR polaron electronic absorption giving rise to Fano-antiresonances. Selective UV excitation of ssDNAs revealed that stacking interactions between thiophene rings and nucleic acid bases can promote the formation of an intermolecular charge transfer complex. The findings inform designers of functional conjugated polymers by identifying that involvement of the scaffold in the photophysics needs to be considered when developing such structures for optoelectronic applications.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

The Ubiquity of High-Energy Nanosecond Fluorescence in DNA Duplexes

During the past few years, several studies reported that a significant part of the intrinsic fluorescence of DNA duplexes decays with surprisingly long lifetimes (1-3 ns) at wavelengths shorter than the ππ* emission of their monomeric constituents. This high-energy nanosecond emission (HENE), hardly discernible in the steady-state fluorescence spectra of most duplexes, was investigated by time-correlated single-photon counting. The ubiquity of HENE contrasts with the paradigm that the longest-lived excited states correspond to low-energy excimers/exciplexes. Interestingly, the latter were found to decay faster than the HENE. So far, the excited states responsible for HENE remain elusive. In order to foster future studies for their characterization, this Perspective presents a critical summary of the experimental observations and the first theoretical approaches. Moreover, some new directions for further work are outlined. Finally, the obvious need for computations of the fluorescence anisotropy considering the dynamic conformational landscape of duplexes is stressed.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

Excimer evolution hampers symmetry-broken charge-separated states

Achieving long-lived symmetry-broken charge-separated states in chromophoric assemblies is quintessential for enhanced performance of artificial photosynthetic mimics. However, the occurrence of energy trap states hinders exciton and charge transport across photovoltaic devices, diminishing power conversion efficiency. Herein, we demonstrate unprecedented excimer formation in the relaxed excited-state geometry of bichromophoric systems impeding the lifetime of symmetry-broken charge-separated states. Core-annulated perylenediimide dimers (SC-SPDI2 and SC-NPDI2) prefer a near-orthogonal arrangement in the ground state and a π-stacked foldamer structure in the excited state. The prospect of an excimer-like state in the foldameric arrangement of SC-SPDI2 and SC-NPDI2 has been rationalized by fragment-based excited state analysis and temperature-dependent photoluminescence measurements. Effective electronic coupling matrix elements in the Franck-Condon geometry of SC-SPDI2 and SC-NPDI2 facilitate solvation-assisted ultrafast symmetry-breaking charge-separation (SB-CS) in a high dielectric environment, in contrast to unrelaxed excimer formation (Ex*) in a low dielectric environment. Subsequently, the SB-CS state dissociates into an undesired relaxed excimer state (Ex) due to configuration mixing of a Frenkel exciton (FE) and charge-separated state in the foldamer structure, downgrading the efficacy of the charge-separated state. The decay rate constant of the FE to SB-CS (k FE→SB-CS) in polar solvents is 8-17 fold faster than that of direct Ex* formation (k FE→Ex*) in non-polar solvent (k FE→SB-CS≫k FE→Ex*), characterized by femtosecond transient absorption (fsTA) spectroscopy. The present investigation establishes the impact of detrimental excimer formation on the persistence of the SB-CS state in chromophoric dimers and offers the requisite of conformational rigidity as one of the potential design principles for developing advanced molecular photovoltaics.

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.

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.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Colour-tuneable solid-state fluorescence of crystalline powders formed from push–pull substituted 2,5-diphenyl-stilbenes

Polycrystalline powders of push–pull substituted stilbenes with various acceptors emit from blue to infrared. Exciton localization on a monomer (in J-like packing) or a stacked dimer (for H-aggregates) avoid exciton migration to the quenching sites.

Deprotonation Dynamics of Guanine Radical Cations†

This review is dedicated to guanine radical cations (G⁺ )^· that are precursors to oxidatively generated damage to DNA. (G⁺ )^· are unstable in neutral aqueous solution and tend to lose a proton. The deprotonation process has been studied by time-resolved absorption experiments in which (G⁺ )^· radicals are produced either by an electron abstraction reaction, using an external oxidant, or by low-energy/low-intensity photoionization of DNA. Both the position of the released proton and the dynamics of the process depend on the secondary DNA structure. While deprotonation in duplex DNA leads to (G-H1)^· radicals, in guanine quadruplexes the (G-H2)^· analogs are observed. Deprotonation in monomeric guanosine proceeds with a time constant of ˜60 ns; in genomic DNA, it is completed within 2 µs; and in guanine quadruplexes, it spans from at least 30 ns to over 50 µs. Such a deprotonation dynamics in four-stranded structures, extended over more than three decades of times, is correlated with the anisotropic structure of DNA and the mobility of its hydration shell. In this case, commonly used second-order reaction models are inappropriate for its description.

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

Nucleobases thin films deposited on nanostructured transparent conductive electrodes for optoelectronic applications

Environmentally-friendly bio-organic materials have become the centre of recent developments in organic electronics, while a suitable interfacial modification is a prerequisite for future applications. In the context of researches on low cost and biodegradable resource for optoelectronics applications, the influence of a 2D nanostructured transparent conductive electrode on the morphological, structural, optical and electrical properties of nucleobases (adenine, guanine, cytosine, thymine and uracil) thin films obtained by thermal evaporation was analysed. The 2D array of nanostructures has been developed in a polymeric layer on glass substrate using a high throughput and low cost technique, UV-Nanoimprint Lithography. The indium tin oxide electrode was grown on both nanostructured and flat substrate and the properties of the heterostructures built on these two types of electrodes were analysed by comparison. We report that the organic-electrode interface modification by nano-patterning affects both the optical (transmission and emission) properties by multiple reflections on the walls of nanostructures and the electrical properties by the effect on the organic/electrode contact area and charge carrier pathway through electrodes. These results encourage the potential application of the nucleobases thin films deposited on nanostructured conductive electrode in green optoelectronic devices.

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.

Excited-State Properties and Relaxation Pathways of Selenium-Substituted Guanine Nucleobase in Aqueous Solution and DNA Duplex

The excited-state properties and relaxation mechanisms after light irradiation of 6-selenoguanine (6SeG) in water and in DNA have been investigated using a quantum mechanics/molecular mechanics (QM/MM) approach with the multistate complete active space second-order perturbation theory (MS-CASPT2) method. In both environments, the S₁ ¹(n_Seπ₅^*) and S₂ ¹(π_Seπ₅^*) states are predicted to be the spectroscopically dark and bright states, respectively. Two triplet states, T₁ ³(π_Seπ₅^*) and T₂ ³(n_Seπ₅^*), are found energetically below the S₂ state. Extending the QM region to include the 6SeG-Cyt base pair slightly stabilizes the S₂ state and destabilizes the S₁, due to hydrogen-bonding interactions, but it does not affect the order of the states. The optimized minima, conical intersections, and singlet-triplet crossings are very similar in water and in DNA, so that the same general mechanism is found. Additionally, for each excited state geometry optimization in DNA, three kind of structures ("up", "down", and "central") are optimized which differ from each other by the orientation of the C═Se group with respect to the surrounding guanine and thymine nucleobases. After irradiation to the S₂ state, 6SeG evolves to the S₂ minimum, near to a S₂/S₁ conical intersection that allows for internal conversion to the S₁ state. Linear interpolation in internal coordinates indicate that the "central" orientation is less favorable since extra energy is needed to surmount the high barrier in order to reach the S₂/S₁ conical intersection. From the S₁ state, 6SeG can further decay to the T₁ ³(π_Seπ₅^*) state via intersystem crossing, where it will be trapped due to the existence of a sizable energy barrier between the T₁ minimum and the T₁/S₀ crossing point. Although this general S₂ → T₁ mechanism takes place in both media, the presence of DNA induces a steeper S₂ potential energy surface, that it is expected to accelerate the S₂ → S₁ internal conversion.

Stacking Effects on Anthraquinone/DNA Charge-Transfer Electronically Excited States

The design of more efficient photosensitizers is a matter of great importance in the field of cancer treatment by means of photodynamic therapy. One of the main processes involved in the activation of apoptosis in cancer cells is the oxidative stress on DNA once a photosensitizer is excited by light. As a consequence, it is very relevant to investigate in detail the binding modes of the chromophore with DNA, and the nature of the electronically excited states that participate in the induction of DNA damage, for example, charge-transfer states. In this work, we investigate the electronic structure of the anthraquinone photosensitizer intercalated into a double-stranded poly(dG-dC) decamer model of DNA. First, the different geometric configurations are analyzed by means of classical molecular dynamics simulations. Then, the excited states for the most relevant poses of anthraquinone inside the binding pocket are computed by an electrostatic-embedding quantum mechanics/molecular mechanics approach, where anthraquinone and one of the nearby guanine residues are described quantum mechanically to take into account intermolecular charge-transfer states. The excited states are characterized as monomer, exciton, excimer, and charge-transfer states based on the analysis of the transition density matrix, and each of these contributions to the total density of states and absorption spectrum is discussed in terms of the stacking interactions. These results are relevant as they represent the footing for future studies on the reactivity of anthraquinone derivatives with DNA and give insights on possible geometrical configurations that potentially favor the oxidative stress of DNA.

Enhanced Rigidity Changes Ultraviolet Absorption: Effect of a Merocyanine Binder on G-Quadruplex Photophysics

The urge to discover selective fluorescent binders to G-quadruplexes (G4s) for rapid diagnosis must be linked to understand the effect that those have on the DNA photophysics. Herein, we report on the electronic excited states of a bound merocyanine dye to c-Myc G4 using extensive multiscale quantum mechanics/molecular mechanics calculations. We find that the absorption spectra of c-Myc G4, both without and with the intercalated dye, are mainly composed of exciton states and mixed local/charge-transfer states. The presence of merocyanine hardly affects the energy range of the guanine absorption or the number of guanines excited. However, it triggers a substantial amount (16%) of detrimental pure charge-transfer states involving oxidized guanines. We identify the rigidity introduced by the probe in G4, reducing the overlap among guanines, as the one responsible for the changes in the exciton and charge-transfer states, ultimately leading to a redshift of the absorption maximum.