Ultrafast Deactivation of an Excited Cytosine−Guanine Base Pair in DNA

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.

XMECP: Reaching State-of-the-Art MECP Optimization in Multiscale Complex Systems

The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.

Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Nonadiabatic Coupling in Trajectory Surface Hopping: How Approximations Impact Excited-State Reaction Dynamics

Photochemical reactions are widely modeled using the popular trajectory surface hopping (TSH) method, an affordable mixed quantum-classical approximation to the full quantum dynamics of the system. TSH is able to account for nonadiabatic effects using an ensemble of trajectories, which are propagated on a single potential energy surface at a time and which can hop from one electronic state to another. The occurrences and locations of these hops are typically determined using the nonadiabatic coupling between electronic states, which can be assessed in a number of ways. In this work, we benchmark the impact of some approximations to the coupling term on the TSH dynamics for several typical isomerization and ring-opening reactions. We have identified that two of the schemes tested, the popular local diabatization scheme and a scheme based on biorthonormal wave function overlap implemented in the OpenMOLCAS code as part of this work, reproduce at a much reduced cost the dynamics obtained using the explicitly calculated nonadiabatic coupling vectors. The other two schemes tested can give different results, and in some cases, even entirely incorrect dynamics. Of these two, the scheme based on configuration interaction vectors gives unpredictable failures, while the other scheme based on the Baeck-An approximation systematically overestimates hopping to the ground state as compared to the reference approaches.

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

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

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

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

Excited-state double proton transfer of 1,8-dihydroxy-2-naphthaldehyde: A MS-CASPT2//CASSCF study

Excited-state double proton transfer (ESDPT) is a controversial issue which has long been plagued with theoretical and experimental communities. Herein, we took 1,8-dihydroxy-2-naphthaldehyde (DHNA) as a prototype and used combined complete active space self-consistent field (CASSCF) and multi-state complete active-space second-order perturbation (MS-CASPT2) methods to investigate ES-DPT and excited-state deactivation pathways of DHNA. Three different tautomer minima of S1-ENOL, S1-KETO-1, and S1-KETO-2 and two crucial conical intersections of S1S0-KETO-1 and S1S0-KETO-2 in.and between the S0 and S1 states were obtained. S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands. In addition, two-dimensional potential energy surfaces (2D-PESs) and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism. Specifically, the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless, whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca. 6.0 kcal/mol. The linear interpolated internal coordinate path connecting S1-KETO-1 (S1-KETO-2) and S1S0-KETO-1 (S1S0-KETO-2) is uphill with a barrier of ca. 12.0 kcal/mol, which will trap DHNA in the S1 state while therefore enabling dual-emission bands. On the other hand, the S1/S0 conical intersections would also prompt the S1 system to decay to the S0 state, which could be to certain extent suppressed by locking the rotation of the C5−C8−C9−O10 dihedral angle. These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.

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

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

A theoretical study on the excited-state deactivation paths for the A-5FU dimer

The photostability of DNA plays a key role in the normal function of organisms. A-5FU is a base pair derivative of the A-T dimer where the methyl group is replaced by a F atom. Here, accurate static TDDFT calculations and non-adiabatic dynamic simulations are used to systematically investigate the excited-state decay paths of the A-5FU dimer related to the proton transfer and the out-of-plane twisting deformation motion of A and 5FU in the 1ππ* and 1nπ* states. CC2 is used to check the accuracy of the current TDDFT calculations. Our results show that the deformation of the C[double bond, length as m-dash]C or C[double bond, length as m-dash]N double bond in A and 5FU provides an efficient pathway for the depopulation of the lowest excited states, which can compete with the excited-state proton transfer paths in the dimer. This finding indicates that monomer-like decay paths could be important for the photostability of weakly hydrogen-bonded DNA base pairs and provide a new insight into the excited-state decay paths in base pairs and their analogues.

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

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

Intermolecular conical intersections in molecular aggregates

Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.

Proton Transfer and Nitro Rotation Tuned Photoisomerization of Artificial Base Pair-ZP

Recently, the successful incorporation of artificial base pairs in genetics has made a significant progress in synthetic biology. The present work reports the proton transfer and photoisomerization of unnatural base pair ZP, which is synthesized from the pyrimidine analog 6-amino-5-nitro-3-(1-β-D-2'-deoxyribo-furanosyl)-2 (1H)-pyridone (Z) and paired with its Watson-Crick complement, the purine analog 2-amino-8-(1'-β-D-2'- deoxyribofuranosyl)-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one (P). To explain the mechanism of proton transfer process, we constructed the relaxed potential energy surfaces (PESs) linking the different tautomers in both gas phase and solution. Our results show that the double proton transfer in the gas phase occurs in a concerted way both in S₀ and S₁ states, while the stepwise mechanism becomes more favorable in solution. The solvent effect can promote the single proton transfer, which undergoes a lower energy barrier in S₁ state due to the strengthened hydrogen bond. In contrast to the excited state ultrafast deactivation process of the natural bases, there is no conical intersection between S₀ and S₁ states along the proton transfer coordinate to activate the decay mechanism in ZP. Of particular relevance to the photophysical properties, charge-transfer character is obviously related to the nitro rotation in S₁ state. We characterized the molecular vibration effect on the electronic properties, which reveals the electronic excitation can be tuned by the rotation-induced structural distortion accompanied with the electron localization on nitro group.

Multiscale Conformational Sampling Reveals Excited-State Locality in DNA Self-Repair Mechanism

Ultraviolet (UV) irradiation is known to be responsible for DNA damage. However, experimental studies in DNA oligonucleotides have shown that UV light can also induce sequence-specific self-repair. Following charge transfer from a guanine adenine sequence adjacent to a cyclobutane pyrimidine dimer (CPD), the covalent bond between the two thymines could be cleaved, recovering the intact base sequence. Mechanistic details promoting the self-repair remained unclear, however. In our theoretical study, we investigated whether optical excitation could directly lead to a charge-transfer state, thereby initiating the repair, or whether the initial excited state remains localized on a single nucleobase. We performed conformational sampling of 200 geometries of the damaged DNA double strand solvated in water and used a hybrid quantum and molecular mechanics approach to compute excited states at the complete active space perturbation level of theory. Analysis of the conformational data set clearly revealed that the excited-state properties are uniformly distributed across the fluctuations of the nucleotide in its natural environment. From the electronic wavefunction, we learned that the electronic transitions remained predominantly local on either adenine or guanine, and no direct charge transfer occurred in the experimentally accessed energy range. The investigated base sequence is not only specific to the CPD repair mechanism but ubiquitously occurs in nucleic acids. Our results therefore give a very general insight into the charge locality of UV-excited DNA, a property that is regarded to have determining relevance in the structural consequences following absorption of UV photons.

Electron-driven proton transfer relieves excited-state antiaromaticity in photoexcited DNA base pairs

The Watson–Crick A·T and G·C base pairs are not only electronically complementary, but also photochemically complementary. Upon UV irradiation, DNA base pairs undergo efficient excited-state deactivation through electron driven proton transfer (EDPT), also known as proton-coupled electron transfer (PCET), at a rate too fast for other reactions to take place. Why this process occurs so efficiently is typically reasoned based on the oxidation and reduction potentials of the bases in their electronic ground states. Here, we show that the occurrence of EDPT can be traced to a reversal in the aromatic/antiaromatic character of the base upon photoexcitation. The Watson–Crick A·T and G·C base pairs are aromatic in the ground state, but the purines become highly antiaromatic and reactive in the first ¹ππ* state, and transferring an electron and a proton to the pyrimidine relieves this excited-state antiaromaticity. Even though proton transfer proceeds along the coordinate of breaking a N–H σ-bond, the chromophore is the π-system of the base, and EDPT is driven by the strive to alleviate antiaromaticity in the π-system of the photoexcited base. The presence and absence of alternative excited-state EDPT routes in base pairs also can be explained by sudden changes in their aromatic and antiaromatic character upon photoexcitation.

The Watson–Crick A·T and G·C base pairs are not only electronically complementary, but also photochemically complementary.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

UV-induced hydrogen transfer in DNA base pairs promoted by dark nπ* states

Formation of an excited-state complex enables ultrafast photorelaxation of dark nπ* states in GC and HC base pairs.

Two-dimensional electronic-vibrational spectroscopic study of conical intersection dynamics: an experimental and electronic structure study

The relaxation from the lowest singlet excited state of the triphenylmethane dyes, crystal violet and malachite green, is studied via two-dimensional electronic-vibrational (2DEV) spectroscopy.

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

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

Nonadiabatic dynamics simulation of photoisomerization mechanism of the second stablest isomer of N-salicilydenemethylfurylamine

The photoisomerization processes of the second stablest isomer in the aromatic Schiff base, N-salicilydenemethylfurylamine, in the gas phase have been studied by static electronic structure calculations and surface-hopping dynamics simulations based on the Zhu-Nakamura theory. Various stable structures are obtained in the optimization because of different orientations of methyl-furyl part with respect to the salicylaldimine part and different orientations of hydroxy group with respect to the benzene ring. Upon photoexcitation into the first excited state, bond isomerization in the salicylaldimine part is completely suppressed until the strong excited-state hydrogen bond is broken. The decay pathway involves two excited-state minima, one in cis-enol form and the other in cis-keto form. After the excited-state proton transfer, twists of bonds lead to a conical intersection between the ground and excited states. After internal conversion around a conical intersection, the molecule is stabilized in cis- or trans-keto form. If the reverse hydrogen transfer process occurs in the ground state, the molecule will finally end up in the cis-enol region. The cis-keto and trans-keto isomers are observed as photoproducts. According to our full-dimensional nonadiabatic dynamics simulations, we find the excited-state intramolecular proton transfer and torsions of three single bonds in the chain to be responsible for photoisomerization of the second stablest isomer of N-salicilydenemethylfurylamine.

A 'bottom up', ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs

The availability of non-radiative decay mechanisms by which photoexcited molecules can revert to their ground electronic state, without experiencing potentially deleterious chemical transformation, is fundamental to molecular photostability. This Perspective Article combines results of new ab initio electronic structure calculations and prior experimental data in an effort to systematise trends in the non-radiative decay following UV excitation of selected families of heterocyclic molecules. We start with the prototypical uni- and bicyclic molecules phenol and indole, and explore the structural and photophysical consequences of incorporating progressively more nitrogen atoms within the respective ring structures en route to the DNA bases thymine, cytosine, adenine and guanine. For each of the latter, we identify low energy non-radiative decay pathways via conical intersections with the ground state potential energy surface accessed by out-of-plane ring deformations. This is followed by summary descriptions and illustrations of selected rival (electron driven H atom transfer) non-radiative excited state decay processes that demand consideration once the nucleobases are merely components in larger biomolecular systems like nucleosides, and both individual and stacked base-pairs.

Proton diffusion dynamics along a diol as a proton-conducting wire in a photo-amphiprotic model system

We investigated the dynamics of excited-state proton transfer (ESPT) of photo-amphiprotic 7-hydroxyquinoline (7HQ) in the presence of a hydrogen (H)-bond bridging diol in a polar aprotic medium. The formation of 1 : 1 H-bonded complexes of 7HQ with various diols of different alkane chain lengths was revealed using steady-state electronic spectroscopy. With femtosecond-resolved fluorescence spectroscopy, cyclic H-bonded 1 : 1 complexes were found to undergo facile ESPT from the acidic enol to the basic imine group of 7HQ via the H-bond bridge. Through quantum chemical calculations, we found that the proton-transfer rate of the well-configured H-bonded complex correlated with the intramolecular H-bond length of a H-bond wiring diol molecule. Noncyclic, singly H-bonded 7HQ with a diol molecule was observed to undergo ESPT once another diol molecule diffuses to the noncyclic complex and accomplishes the formation of a reactive cyclic H-bonded 7HQ-(diol)₂ complex, which was evidenced by the observation that the overall proton-transfer rate constant decreases when a longer-chain diol was used as the bridging wire part. The kinetic isotope effect on the proton relay was investigated to confirm that the nature of the activation barrier for the proton diffusion along the wire is isotope-sensitive proton tunnelling, while for the non-cyclic configuration, the isotope-insensitive H-bond bridge formation is a prerequisite for ESPT.

Photoactivated proton coupled electron transfer in DNA: insights from quantum mechanical calculations

The energetics of the two main proton coupled electron transfer processes that could occur in DNA are determined by means of time dependent-DFT calculations, using the M052X functional and the polarizable continuum model to include solvent effect.

Is UV-Induced Electron-Driven Proton Transfer Active in a Chemically Modified A·T DNA Base Pair?

Transient electronic and vibrational absorption spectroscopies have been used to investigate whether UV-induced electron-driven proton transfer (EDPT) mechanisms are active in a chemically modified adenine-thymine (A·T) DNA base pair. To enhance the fraction of biologically relevant Watson-Crick (WC) hydrogen-bonding motifs and eliminate undesired Hoogsteen structures, a chemically modified derivative of A was synthesized, 8-(tert-butyl)-9-ethyladenine (8tBA). Equimolar solutions of 8tBA and silyl-protected T nucleosides in chloroform yield a mixture of WC pairs, reverse WC pairs, and residual monomers. Unlike previous transient absorption studies of WC guanine-cytosine (G·C) pairs, no clear spectroscopic or kinetic evidence was identified for the participation of EDPT in the excited-state relaxation dynamics of 8tBA·T pairs, although ultrafast (sub-100 fs) EDPT cannot be discounted. Monomer-like dynamics are proposed to dominate in 8tBA·T.

Role of Electron-Driven Proton-Transfer Processes in the Ultrafast Deactivation of Photoexcited Anionic 8-oxoGuanine-Adenine and 8-oxoGuanine-Cytosine Base Pairs

It has been reported that 8-oxo-7,8-dihydro-guanosine (8-oxo-G), which is the main product of oxidative damage of DNA, can repair cyclobutane pyrimidine dimer (CPD) lesions when incorporated into DNA or RNA strands in proximity to such lesions. It has therefore been suggested that the 8-oxo-G nucleoside may have been a primordial precursor of present-day flavins in DNA or RNA repair. Because the electron transfer leading to the splitting of a thymine-thymine pair in a CPD lesion occurs in the photoexcited state, a reasonably long excited-state lifetime of 8-oxo-G is required. The neutral (protonated) form of 8-oxo-G exhibits a very short (sub-picosecond) intrinsic excited-state lifetime which is unfavorable for repair. It has therefore been argued that the anionic (deprotonated) form of 8-oxo-G, which exhibits a much longer excited-state lifetime, is more likely to be a suitable cofactor for DNA repair. Herein, we have investigated the exited-state quenching mechanisms in the hydrogen-bonded complexes of deprotonated 8-oxo-G^- with adenine (A) and cytosine (C) using ab initio wave-function-based electronic-structure calculations. The calculated reaction paths and potential-energy profiles reveal the existence of barrierless electron-driven inter-base proton-transfer reactions which lead to low-lying S₁/S₀ conical intersections. The latter can promote ultrafast excited-state deactivation of the anionic base pairs. While the isolated deprotonated 8-oxo-G^- nucleoside may have been an efficient primordial repair cofactor, the excited states of the 8-oxo-G^--A and 8-oxo-G^--C base pairs are likely too short-lived to be efficient electron-transfer repair agents.

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

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

Solvent effects on the excited-state double proton transfer mechanism in the 7-azaindole dimer: a TDDFT study with the polarizable continuum model

Solvent effects on the excited-state double proton transfer (ESDPT) mechanism in the 7-azaindole (7AI) dimer were investigated using the time-dependent density functional theory (TDDFT) method.

Ultrafast structural dynamics of photoexcited adenine

Ultraviolet Resonance Raman (UVRR) spectroscopy derives distinct electronic properties of adenine in the L_a (260 nm) and B_b (210 nm) excited states.

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.

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.

Accidental degeneracy in the spiropyran radical cation: charge transfer between two orthogonal rings inducing ultra-efficient reactivity

We unravel an original photoswitching mechanism in spiropyran radical cation.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

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.

Electronic excitations in Guanine quadruplexes

Guanine rich DNA strands, such as those encountered at the extremities of human chromosomes, have the ability to form four-stranded structures (G-quadruplexes) whose building blocks are guanine tetrads. G-quadruplex structures are intensively studied in respect of their biological role, as targets for anticancer therapy and, more recently, of their potential applications in the field of molecular electronics. Here we focus on their electronic excited states which are compared to those of non-interacting mono-nucleotides and those of single and double stranded structures. Particular emphasis is given to excited state relaxation processes studied by time-resolved fluorescence spectroscopy from femtosecond to nanosecond time scales. They include ultrafast energy transfer and trapping of ππ* excitations by charge transfer states. The effect of various structural parameters, such as the nature of the metal cations located in the central cavity of G-quadruplexes, the number of tetrads or the conformation of the constitutive single strands, are examined.

Non-radiative deactivation in phenol–pyridine complex: theoretical study

A comparative theoretical study of the nonradiative decay mechanisms to the ground state of phenol-pyridine complex under Cs symmetry (Pl and Pe conformers) and without symmetry constraints has been performed by ab initio calculations.