Hydrogen bonding regulates the monomeric nonradiative decay of adenine in DNA strands

New insights into the photocyclization reaction of a popular diarylethene switch: a nonadiabatic molecular dynamics study

Diarylethene (DAE) molecular switches have continued to attract the attention of researchers for over 20 years. Their remarkable photophysical properties endow them with countless applications in photonics and molecular technologies. However, despite extensive experimental and theoretical research, the mechanism of DAE photoswitching is not yet fully rationalized. In this work, we investigate the ring closure dynamics of a popular DAE switch, 1,2-bis(3-methyl-5-phenyl-2 thienyl)perfluorocyclopentene (PT), using nonadiabatic molecular dynamics (NAMD) simulations. Employing the fewest switches surface hopping protocol, along with the semi-empirical multireference ODM2/MRCI-SD method, we investigate possible reaction pathways for this photoprocess, as well as their timescales and resulting photoproducts. Furthermore, using a dynamic configuration-space sampling procedure, we elucidate the role of triplet states in the photocyclization of PT, supporting available experimental data for the closely related DMPT molecule, which indicate an ultrafast intersystem crossing (ISC) transition competing with the singlet-driven photoswitching reaction. Our findings not only corroborate experimental studies on DAE switches, but also provide new mechanistic insights into the potential use in the rational design of DAE switches tailored for specific technological applications.

Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube

Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E₁₁ excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E₁₁ excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out^-in⁺〉 and |out⁺in^-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E₂₂ state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.

Photo-cycloreversion mechanism in diarylethenes revisited: A multireference quantum-chemical study at the ODM2/MRCI level

Photoswitchable diarylethenes (DAEs), over years of intense fundamental and applied research, have been established among the most commonly chosen molecular photoswitches, often employed as controlling units in molecular devices and smart materials. At the same time, providing reliable explanation for their photophysical behavior, especially the mechanism of the photo-cycloreversion transformation, turned out to be a highly challenging task. Herein, we investigate this mechanism in detail by means of multireference semi-empirical quantum chemistry calculations, allowing, for the first time, for a balanced treatment of the static and dynamic correlation effects, both playing a crucial role in DAE photochemistry. In the course of our study, we find the second singlet excited state of double electronic-excitation character to be the key to understanding the nature of the photo-cycloreversion transformation in DAE molecular photoswitches.

"On-The-Fly" Non-Adiabatic Dynamics Simulations on Photoinduced Ring-Closing Reaction of a Nucleoside-Based Diarylethene Photoswitch

Nucleoside-based diarylethenes are emerging as an especial class of photochromic compounds that have potential applications in regulating biological systems using noninvasive light with high spatio-temporal resolution. However, relevant microscopic photochromic mechanisms at atomic level of these novel diarylethenes remain to be explored. Herein, we have employed static electronic structure calculations (MS-CASPT2//M06-2X, MS-CASPT2//SA-CASSCF) in combination with non-adiabatic dynamics simulations to explore the related photoinduced ring-closing reaction of a typical nucleoside-based diarylethene photoswitch, namely, PS-IV. Upon excitation with UV light, the open form PS-IV can be excited to a spectroscopically bright S₁ state. After that, the molecule relaxes to the conical intersection region within 150 fs according to the barrierless relaxed scan of the C1–C6 bond, which is followed by an immediate deactivation to the ground state. The conical intersection structure is very similar to the ground state transition state structure which connects the open and closed forms of PS-IV, and therefore plays a crucial role in the photochromism of PS-IV. Besides, after analyzing the hopping structures, we conclude that the ring closing reaction cannot complete in the S₁ state alone since all the C1–C6 distances of the hopping structures are larger than 2.00 Å. Once hopping to the ground state, the molecules either return to the original open form of PS-IV or produce the closed form of PS-IV within 100 fs, and the ring closing quantum yield is estimated to be 56%. Our present work not only elucidates the ultrafast photoinduced pericyclic reaction of the nucleoside-based diarylethene PS-IV, but can also be helpful for the future design of novel nucleoside-based diarylethenes with better performance.

The influence of base pair tautomerism on single point mutations in aqueous DNA

The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

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).

Nonadiabatic dynamics simulations on internal conversion and intersystem crossing processes in gold(i) compounds

The position at which the second gold(i)-phosphine group is attached was experimentally found to play a noticeable role in intersystem crossing rates of gold(i) naphthalene derivatives. However, the physical origin is ambiguous. Herein we have employed generalized trajectory-based surface-hopping dynamics simulations to simulate the excited-state relaxation dynamics of these gold(i) naphthalene compounds including both the intersystem crossing process from the initially populated first excited singlet states S₁ to triplet manifolds and internal conversion processes within these triplet states. Our predicted intersystem crossing rates are consistent with experiments very well. On the basis of the present results, we have found that (1) ultrafast and subpicosecond intersystem crossing processes are mainly caused by small energy gaps and large spin-orbit couplings between S₁ and T_n; (2) adding the second gold(i)-phosphine group does not increase spin-orbit couplings between S₁ and T_n but decrease their values remarkably, which implies that heavy-atom effects are state-specific, not state-universal; (3) the position at which the second gold(i)-phosphine group is attached has a remarkable influence on the electronic structures of S₁ and T_n and their relative energies, which affect energy gaps and spin-orbit couplings between S₁ and T_n and eventually modulate intersystem crossing rates from S₁ to T_n. These new insights are very useful for the design of gold-containing compounds with excellent photoluminescence properties. Finally, this work also exemplifies that different isomers of a compound could have distinct excited-state relaxation dynamics.

Ultrafast action chemistry in slow motion: atomistic description of the excitation and fluorescence processes in an archetypal fluorescent protein

We report quantum mechanical/molecular mechanical non-adiabatic molecular dynamics simulations on the electronically excited state of green fluorescent protein mutant S65T/H148D. We examine the driving force of the ultrafast (τ < 50 fs) excited-state proton transfer unleashed by absorption in the A band at 415 nm and propose an atomistic description of the two dynamical regimes experimentally observed [Stoner Ma et al., J. Am. Chem. Soc., 2008, 130, 1227]. These regimes are explained in terms of two sets of successive dynamical events: first the proton transfers quickly from the chromophore to the acceptor Asp148. Thereafter, on a slower time scale, there are geometrical changes in the cavity of the chromophore that involve the distance between the chromophore and Asp148, the planarity of the excited-state chromophore, and the distance between the chromophore and Tyr145. We find two different non-radiative relaxation channels that are operative for structures in the reactant region and that can explain the mismatch between the decay of the emission of A* and the rise of the emission of I*, as well as the temperature dependence of the non-radiative decay rate.

Ultrafast intersystem crossing for nitrophenols: ab initio nonadiabatic molecular dynamics simulation

Ultrafast intersystem crossing mechanisms for two p- and m-nitrophenol groups (PNP and MNP) have been investigated using ab initio nonadiabatic molecular dynamics simulations at the 6SA-CASSCF level of theory.

Excited-State Decay Paths in Tetraphenylethene Derivatives

The photophysical properties of tetraphenylethene (TPE) compounds may differ widely depending on the substitution pattern, for example, with regard to the fluorescence quantum yield ϕ_f and the propensity to exhibit aggregation-induced emission (AIE). We report combined electronic structure calculations and nonadiabatic dynamics simulations to study the excited-state decay mechanisms of two TPE derivatives with four methyl substituents, either in the meta position (TPE-4mM, ϕ_f = 0.1%) or in the ortho position (TPE-4oM, ϕ_f = 64.3%). In both cases, two excited-state decay pathways may be relevant, namely, photoisomerization around the central ethylenic double bond and photocyclization involving two adjacent phenyl rings. In TPE-4mM, the barrierless S₁ cyclization is favored; it is responsible for the ultralow fluorescence quantum yield observed experimentally. In TPE-4oM, both the S₁ photocyclization and photoisomerization paths are blocked by non-negligible barriers, and fluorescence is thus feasible. Nonadiabatic dynamics simulations with more than 1000 surface hopping trajectories show ultrafast cyclization upon photoexcitation of TPE-4mM, whereas TPE-4oM remains unreactive during the 1 ps simulations. We discuss the chances for spectroscopic detection of the postulated cyclic photoproduct of TPE-4mM and the relevance of our findings for the AIE process.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Tuning excited-state-intramolecular-proton-transfer (ESIPT) process and emission by cocrystal formation: a combined experimental and theoretical study

The formation of two-component molecular cocrystals can lead to the tunable excited state intramolecular proton transfer (ESIPT) process and emission, as first confirmed by both experimental and computational studies.

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.

Life in the light: nucleic acid photoproperties as a legacy of chemical evolution

Photophysical investigations of the canonical nucleobases that make up DNA and RNA during the past 15 years have revealed that excited states formed by the absorption of UV radiation decay with subpicosecond lifetimes (i.e., <10(-12) s). Ultrashort lifetimes are a general property of absorbing sunscreen molecules, suggesting that the nucleobases are molecular survivors of a harsh UV environment. Encoding the genome using photostable building blocks is an elegant solution to the threat of photochemical damage. Ultrafast excited-state deactivation strongly supports the hypothesis that UV radiation played a major role in shaping molecular inventories on the early Earth before the emergence of life and the subsequent development of a protective ozone shield. Here, we review the general physical and chemical principles that underlie the photostability, or "UV hardiness", of modern nucleic acids and discuss the possible implications of these findings for prebiotic chemical evolution. In RNA and DNA strands, much longer-lived excited states are observed, which at first glance appear to increase the risk of photochemistry. It is proposed that the dramatically different photoproperties that emerge from assemblies of photostable building blocks may explain the transition from a world of molecular survival to a world in which energy-rich excited electronic states were eventually tamed for biological purposes such as energy transduction, signaling, and repair of the genetic machinery.

Multiple Decay Mechanisms and 2D-UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine-Uracil Monophosphate

The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D-UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine (1) La state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine (1) Lb and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter-base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long-lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the (1) Lb , S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D-UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm(-1) in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D-UV spectroscopy to disentangle the photophysics of multichromophoric systems.

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.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters

Semiempirical orthogonalization-corrected methods (OM1, OM2, and OM3) go beyond the standard MNDO model by explicitly including additional interactions into the Fock matrix in an approximate manner (Pauli repulsion, penetration effects, and core–valence interactions), which yields systematic improvements both for ground-state and excited-state properties. In this Article, we describe the underlying theoretical formalism of the OMx methods and their implementation in full detail, and we report all relevant OMx parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. For a standard set of mostly organic molecules commonly used in semiempirical method development, the OMx results are found to be superior to those from standard MNDO-type methods. Parametrized Grimme-type dispersion corrections can be added to OM2 and OM3 energies to provide a realistic treatment of noncovalent interaction energies, as demonstrated for the complexes in the S22 and S66×8 test sets.

Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations

Sequential vs. concerted S₁ relaxation pathways.

Generalized trajectory surface-hopping method for internal conversion and intersystem crossing

Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

Stabilization of Mixed Frenkel-Charge Transfer Excitons Extended Across Both Strands of Guanine-Cytosine DNA Duplexes

The photoreactive pathways that may lead to DNA damage depend crucially upon the nature of the excited electronic states. The study of alternating guanine-cytosine duplexes by fluorescence spectroscopy and quantum mechanical calculations identifies a novel type of excited states that can be populated following UVB excitation. These states, denoted High-energy Emitting Long-lived Mixed (HELM) states, extend across both strands and arise from mixing between cytosine Frenkel excitons and guanine-to-cytosine charge transfer states. They emit at energies higher than ππ* states localized on single bases, survive for several nanoseconds, are sensitive to the ionic strength of the solution, and are strongly affected by the structural transition from the B form to the Z form. Their impact on the formation of lesions of the genetic code needs to be assessed.

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.

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.

Photoswitching of salicylidene methylamine: a theoretical photodynamics study

Photoswitching of simple photochromic molecules attracts substantial attention because of its possible role in future photon-driven molecular electronics. Here we model the full photoswitching cycle of a minimal photochromic Schiff base–salicylidene methylamine (SMA). We perform semiempirical nonadiabatic on-the-fly photodynamics simulations at the OM2/MRCI level and thoroughly analyze the structural time evolution and switching efficiency of the system. We also identify and examine in detail the crucial steps in the SMA photochemistry ruled by excited-state intramolecular proton transfer. The results place the investigated model aromatic Schiff base among the promising candidates for novel photoswitching molecular materials. Our study also shows the potential of the semiempirical multireference photodynamics simulations as a tool for early stage molecular photodevice design.

Slow deactivation channels in UV-photoexcited adenine DNA

The molecular mechanism for removing the excess energy in DNA bases is responsible for the high photostability of DNA and is thus the subject of intense theoretical/computational investigation. To understand why the excited state decay of the stacked bases is significantly longer than that of the monomers, we carried out electronic structure calculations on an adenine monomer and an aqueous (dA)5 oligonucleotide employing the CASPT2//CASSCF and CASPT2//CASSCF/AMBER levels of theory. The newly-found bright excited state pair Sstack1((1)ππ*) and Sstack2((1)ππ*) of d(A)5, originated from base stacking, is of intra-base charge transfer nature and occurs in different stacked bases with charge transfer along opposite directions. Two slow deactivation channels of d(A)5 were proposed as a result of the sizable barriers along the relaxation paths starting from the FC point of the Sstack1((1)ππ*) state. The SN1P((1)nπ*) state of d(A)5 serves as an intermediate state in one relaxation channel, to which a nonadiabatic decay from the Sstack1((1)ππ*) state occurs in an energy degeneracy region. A relatively high barrier in this state is found and attributed to the steric hindrance of the DNA environment due to the large NH2 group twisting, which gives a weak and red-shifted fluorescence. Another direct relaxation channel, induced by the C2-H2 bond twisting motion, is found to go through a conical intersection between the Sstack1((1)ππ*) and the ground state. The barrier found here enables fluorescence from the Sstack1((1)ππ*) state and may explain the bright state emission observed in the fluorescence upconversion measurements. The inter-molecular SCT((1)ππ*) state may be involved in the slow relaxation process of the photoexcited adenine oligomers through efficient internal conversion to the intra-base Sstack1((1)ππ*) state.

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

In this work, we study the photophysical properties of an adenine-based oligonucleotide using an ensemble of about 200 configurations obtained from molecular dynamics simulations. Specifically, a QM/MM approach is used to obtain the excited-state energies and properties of (dA)20(dT)20 with a dimer of π-stacked adenine bases included in the quantum region. The absorption and circular dichroism spectra are computed and analyzed using the algebraic diagrammatic construction through second order level of theory method (ADC(2)) combined with classical mechanics. We find that the experimentally observed red-shifted shoulder in the absorption spectrum is due to excitonic interactions, while charge-transfer states are present within the absorption band at the higher-energy end of the spectrum. More importantly, low-energy states with charge-transfer mixing exist, which could lead to excimers and bonded excimers. These observations suggest that mixing between charge-transfer and excitonic states plays an important role in the photophysics of oligonucleotides. They also highlight the importance of taking into account the conformational flexibility of the oligonucleotide when investigating photophysical properties.