RNA Environment Is Responsible for Decreased Photostability of Uracil

Theoretical insight into photodeactivation mechanisms of adenine-uracil and adenine-thymine nucleobase pairs

In this work, several plausible intra- and intermolecular photoinduced processes of the Watson-Crick base pairs of adenine with uracil (A-U) or thymine (A-T) according to the results of spin component scaling variant of algebraic diagrammatic construction up to the second order [SCS-ADC(2)] calculations are discussed. Although widely explored, these systems lack complete characterization of possible intramolecular relaxation channels perturbed by intermolecular interactions. In particular, we address the still open debate on photodeactivation via purine-ring puckering at the C2 or C6-atom position of adenine. We also show that the presence of low-lying, long-lived 1nπ* states can be a significant factor in hindering relaxation via an electron-driven proton transfer process, as the population of these states can lead to an efficient intersystem crossing to a triplet manifold, the estimated rate of which is 1.6 × 1010 s-1 which exceeds the corresponding internal conversion to the ground state by an order of magnitude. Additionally, the SCS variant of the ADC(2) method is shown to provide a more balanced description of valence and charge-transfer excited states.

Following the Nonadiabatic Ultrafast Dynamics of Uracil via Simulated X-ray Absorption Spectra

The nucleobase uracil exhibits high photostability due to ultrafast relaxation processes mediated by conical intersections (CoIns), where the interplay between nuclear and electron dynamics becomes crucial. In our previous study, we observed seemingly long-lived traces of electronic coherence for the relaxation process through the S2/S1 CoIn by applying our ansatz for coupled nuclear and electron dynamics in molecules (NEMol). In this work, we theoretically investigate how time-dependent transient X-ray absorption spectroscopy can be used to observe this ultrafast dynamics. Therefore, we calculated X-ray absorption spectra (XAS) for the oxygen K-edge, using a multireference protocol in combination with NEMol dynamics. Thus, we have access to both the transient XAS based on the nuclear wavepacket dynamics and the modulation of the signals caused by the electronic coherence induced by the excitation process and the presence of a CoIn seam. In both cases, we were able to qualitatively predict its influence on the resulting XAS.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

On the propensity of formation of cyclobutane dimers in face-to-face and face-to-back uracil stacks in solution

UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.

How the Interplay among Conformational Disorder, Solvation, Local, and Charge-Transfer Excitations Affects the Absorption Spectrum and Photoinduced Dynamics of Perylene Diimide Dimers: A Molecular Dynamics/Quantum Vibronic Approach

In this contribution we present a mixed quantum-classical dynamical approach for the computation of vibronic absorption spectra of molecular aggregates and their nonadiabatic dynamics, taking into account the coupling between local excitations (LE) and charge-transfer (CT) states. The approach is based on an adiabatic (Ad) separation between the soft degrees of freedom (DoFs) of the system and the stiff vibrations, which are described by the quantum dynamics (QD) of wave packets (WPs) moving on the coupled potential energy surfaces (PESs) of the LE and CT states. These PESs are described with a linear vibronic coupling (LVC) Hamiltonian, parameterized by an overlap-based diabatization on the grounds of time-dependent density functional theory computations. The WPs time evolution is computed with the multiconfiguration time-dependent Hartree method, using effective modes defined through a hierarchical representation of the LVC Hamiltonian. The soft DoFs are sampled with classical molecular dynamics (MD), and the coupling between the slow and fast DoFs is included by recomputing the key parameters of the LVC Hamiltonians, specifically for each MD configuration. This method, named Ad-MD|gLVC, is applied to a perylene diimide (PDI) dimer in acetonitrile and water solutions, and it is shown to accurately reproduce the change in the vibronic features of the absorption spectrum upon aggregation. Moreover, the microscopic insight offered by the MD trajectories allows for a detailed understanding of the role played by the fluctuation of the aggregate structure on the shape of the vibronic spectrum and on the population of LE and CT states. The nonadiabatic QD predicts an extremely fast (∼50 fs) energy transfer between the two LEs. CT states have only a moderate effect on the absorption spectrum, despite the fact that after photoexcitation they are shown to acquire a fast and non-negligible population, highlighting their relevance in dictating the charge separation and transport in PDI-based optical devices.

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

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

Conical Intersection Passages of Molecules Probed by X-ray Diffraction and Stimulated Raman Spectroscopy

Conical intersections (CoIns) play an important role in ultrafast relaxation channels. Their monitoring remains a formidable experimental challenge. We theoretically compare the probing of the S₂ → S₁ CoIn passage in 4-thiouracil by monitoring its vibronic coherences, using off-resonant X-ray-stimulated Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). The quantum nuclear wavepacket (WP) dynamics provides an accurate picture of the photoinduced dynamics. Upon photoexcitation, the WP oscillates among the Franck-Condon point, the S₂ minimum, and the CoIn with a 70 fs period. A vibronic coherence first emerges at 20 fs and can be observed until the S₂ state is fully depopulated. The distribution of the vibronic frequencies involved in the coherence is recorded by the TRUECARS spectrogram. The TRXD signal provides spatial images of electron densities associated with the CoIn. In combination, the two signals provide a complementary picture of the nonadiabatic passage, which helps in the study of the underlying photophysics in thiobases.

The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes

This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.

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

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

Waveform control of molecular dynamics close to a conical intersection

Conical intersections are ubiquitous in chemical systems but, nevertheless, extraordinary points on the molecular potential energy landscape. They provide ultra-fast radiationless relaxation channels, their topography influences the product branching, and they equalize the timescales of the electron and nuclear dynamics. These properties reveal optical control possibilities in the few femtosecond regime. In this theoretical study, we aim to explore control options that rely on the carrier envelope phase of a few-cycle IR pulse. The laser interaction creates an electronic superposition just before the wave packet reaches the conical intersection. The imprinted phase information is varied by the carrier envelope phase to influence the branching ratio after the conical intersection. We test and analyze this scenario in detail for a model system and show to what extent it is possible to transfer this type of control to a realistic system like uracil.

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.

Photoprotecting Uracil by Coupling with Lossy Nanocavities

We analyze how the photorelaxation dynamics of a molecule can be controlled by modifying its electromagnetic environment using a nanocavity mode. In particular, we consider the photorelaxation of the RNA nucleobase uracil, which is the natural mechanism to prevent photodamage. In our theoretical work, we identify the operative conditions in which strong coupling with the cavity mode can open an efficient photoprotective channel, resulting in a relaxation dynamics twice as fast as the natural one. We rely on a state-of-the-art chemically detailed molecular model and a non-Hermitian Hamiltonian propagation approach to perform full-quantum simulations of the system dissipative dynamics. By focusing on the photon decay, our analysis unveils the active role played by cavity-induced dissipative processes in modifying chemical reaction rates, in the context of molecular polaritonics. Remarkably, we find that the photorelaxation efficiency is maximized when an optimal trade-off between light-matter coupling strength and photon decay rate is satisfied. This result is in contrast with the common intuition that increasing the quality factor of nanocavities and plasmonic devices improves their performance. Finally, we use a detailed model of a metal nanoparticle to show that the speedup of the uracil relaxation could be observed via coupling with a nanosphere pseudomode, without requiring the implementation of complex nanophotonic structures.

Experimental and theoretical studies on thymine photodimerization mediated by oxidatively generated DNA lesions and epigenetic intermediates

Combined spectroscopic and computational studies reveal that, in spite of their structural similarities, 5-formyluracil and 5-formylcytosine photosensitize cyclobutane thymine dimers through two different types of mechanisms.

UV photostability of three 2-aminoazoles with key roles in prebiotic chemistry on the early earth

Three related molecules in the 2-aminoazole family are potentially important for prebiotic chemistry: 2-aminooxazole, 2-aminoimidazole, and 2-aminothiazole, which can provide critical functions as an intermediate in nucleotide synthesis, a nucleotide activating agent, and a selective agent, respectively. Here, we examine the wavelength-dependent photodegradation of these three molecules under mid-range UV light (210-290 nm). We then assess the implications of the observed degradation rates for the proposed prebiotic roles of these compounds. We find that all three 2-aminoazoles degrade under UV light, with half lives ranging from ≈7-100 hours under a solar-like spectrum. 2-Aminooxazole is the least photostable, while 2-aminoimidazole is the most photostable. The relative photostabilities are consistent with the order in which these molecules would be used prebiotically: AO is used first to build nucleotides and AI is used last to activate them.