Challenges in Simulating Light-Induced Processes in DNA

The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study

Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.

Short-Range Charge Transfer in DNA Base Triplets: Real-Time Tracking of Coherent Fluctuation Electron Transfer

The short-range charge transfer of DNA base triplets has wide application prospects in bioelectronic devices for identifying DNA bases and clinical diagnostics, and the key to its development is to understand the mechanisms of short-range electron dynamics. However, tracing how electrons are transferred during the short-range charge transfer of DNA base triplets remains a great challenge. Here, by means of ab initio molecular dynamics and Ehrenfest dynamics, the nuclear-electron interaction in the thymine-adenine-thymine (TAT) charge transfer process is successfully simulated. The results show that the electron transfer of TAT has an oscillating phenomenon with a period of 10 fs. The charge density difference proves that the charge transfer proportion is as high as 59.817% at 50 fs. The peak position of the hydrogen bond fluctuates regularly between -0.040 and -0.056. The time-dependent Marcus-Levich-Jortner theory proves that the vibrational coupling between nucleus and electron induces coherent electron transfer in TAT. This work provides a real-time demonstration of the short-range coherent electron transfer of DNA base triplets and establishes a theoretical basis for the design and development of novel biological probe molecules.

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.

The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics

The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Excited-state dynamics of CH2I2 and CH2IBr studied with UV-pump VUV-probe momentum-resolved photoion spectroscopy

We perform time-resolved ionization spectroscopy measurements of the excited state dynamics of CH₂I₂ and CH₂IBr following photoexcitation in the deep UV. The fragment ions produced by ionization with a vacuum-ultraviolet probe pulse are measured with velocity map imaging, and the momentum resolved yields are compared with trajectory surface hopping calculations of the measurement observable. Together with recent time-resolved photoelectron spectroscopy measurements of the same dynamics, these results provide a detailed picture of the coupled electronic and nuclear dynamics involved. Our measurements highlight the non-adiabatic coupling between electronic states, which leads to notable differences in the dissociation dynamics for the two molecules.

Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities

We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.

Transient X-ray Absorption Spectral Fingerprints of the S1 Dark State in Uracil

Low-lying dark nπ* states play an important role in many photophysical and photochemical processes of organic chromophores. Transient X-ray absorption spectroscopy (TXAS) provides a powerful technique for probing the dynamics of valence states by exciting the electrons into high-lying core excited states. We employ multiconfigurational self-consistent field calculations to investigate the TXAS of uracil along its nonradiative photodecay pathways. An open issue is whether dark nπ* state S₁ (n is the lone pair localized on an oxygen atom) is accessible when bright ππ* state S₂ is selectively excited. Vertical core excitations were calculated along the potential energy surfaces of the three lowest states, S₀-S₂, interpolated between two minima and two minimum-energy conical intersections. Computed TXAS data from the C, N, and O K edges show distinct spectral fingerprints of the dark state in all spectral regimes. At the O 1s edge, the nπ* state has a very strong absorption at 526-527 eV, while at the C (N) 1s edge, by contrast, there is almost zero (very weak) absorption at 279-282 eV (397-398 eV). All K-edge spectra can be used to sensitively detect the dark states. Our proposed O 1s feature has already been observed in a recent TXAS experiment with thymine. Natural transition orbital analysis is used to interpret all dominant features of the three lowest-valence states along the reaction coordinate and reveal some important valence fine-structure information from the core excitation.

Electronic Relaxation Dynamics of UV-Photoexcited 2-Aminopurine-Thymine Base Pairs in Watson-Crick and Hoogsteen Conformations

The fluorescent analogue 2-aminopurine (2AP) of the canonical nucleobase adenine (6-aminopurine) base-pairs with thymine (T) without disrupting the helical structure of DNA. It therefore finds frequent use in molecular biology for probing DNA and RNA structures and conformational dynamics. However, detailed understanding of the processes responsible for fluorescence quenching remains largely elusive on a fundamental level. Although attempts have been made to ascribe decreased excited-state lifetimes to intrastrand charge-transfer and stacking interactions, possible influences from dynamic interstrand H-bonding have been widely ignored. Here, we investigate the electronic relaxation of UV-excited 2AP·T in Watson-Crick (WC) and Hoogsteen (HS) conformations. Although the WC conformation features slowed-down, monomer-like electronic relaxation in τ ∼ 1.6 ns toward ground-state recovery and triplet formation, the dynamics associated with 2AP·T in the HS motif exhibit faster deactivation in τ ∼ 70 ps. As recent research has revealed abundant transient interstrand H-bonding in the Hoogsteen motif for duplex DNA, the established model for dynamic fluorescence quenching may need to be revised in the light of our results. The underlying supramolecular photophysical mechanisms are discussed in terms of a proposed excited-state double-proton transfer as an efficient deactivation channel for recovery of the HS species in the electronic ground state.

Exciton Localization on Ru-Based Photosensitizers Induced by Binding to Lipid Membranes

The characterization of electronic properties of metal complexes embedded in membrane environments is of paramount importance to develop efficient photosensitizers in optogenetic applications. Molecular dynamics and QM/MM simulations together with quantitative wave function analysis reveal a directional electronic redistribution of the exciton formed upon excitation of [Ru(bpy)₂(bpy-C17)]²⁺ when going from water to a lipid bilayer, despite the fact that the media influence neither the metal-to-ligand charge-transfer character nor the excitation energy of the absorption spectra. When the photosensitizer is embedded into the DOPC lipid membrane, exciton population is mainly located in the bypyridyl sites proximal to the positively charged surface of the bilayer due to electrostatic interactions. This behavior shows that the electronic structure of metal complexes can be controlled through the binding to external species, underscoring the crucial role of the environment in directing the electronic flow upon excitation and thus helping rational tuning of optogenetic agents.

Photorelaxation and Photorepair Processes in Nucleic and Amino Acid Derivatives

Understanding the fundamental interaction between electromagnetic radiation and matter is essential for a large number of phenomena, with significance to civilization.[...].

Machine learning molecular dynamics for the simulation of infrared spectra

Machine learning has emerged as an invaluable tool in many research areas. In the present work, we harness this power to predict highly accurate molecular infrared spectra with unprecedented computational efficiency. To account for vibrational anharmonic and dynamical effects - typically neglected by conventional quantum chemistry approaches - we base our machine learning strategy on ab initio molecular dynamics simulations. While these simulations are usually extremely time consuming even for small molecules, we overcome these limitations by leveraging the power of a variety of machine learning techniques, not only accelerating simulations by several orders of magnitude, but also greatly extending the size of systems that can be treated. To this end, we develop a molecular dipole moment model based on environment dependent neural network charges and combine it with the neural network potential approach of Behler and Parrinello. Contrary to the prevalent big data philosophy, we are able to obtain very accurate machine learning models for the prediction of infrared spectra based on only a few hundreds of electronic structure reference points. This is made possible through the use of molecular forces during neural network potential training and the introduction of a fully automated sampling scheme. We demonstrate the power of our machine learning approach by applying it to model the infrared spectra of a methanol molecule, n-alkanes containing up to 200 atoms and the protonated alanine tripeptide, which at the same time represents the first application of machine learning techniques to simulate the dynamics of a peptide. In all of these case studies we find an excellent agreement between the infrared spectra predicted via machine learning models and the respective theoretical and experimental spectra.

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

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