Combining ab initio multiconfigurational and Free Energy Gradient methods to study the π–π* excited state structure and properties of uracil in water

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

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

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins

Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.

Resolving Ultrafast Photoinitiated Dynamics of the Hachimoji 5-Aza-7-Deazaguanine Nucleobase: Impact of Synthetically Expanding the Genetic Alphabet

The guanine derivative, 5-aza-7-deazaguanine (^5N7C G) has recently been proposed as one of four unnatural bases, termed Hachimoji (8-letter) to expand the genetic code. We apply steady-state and time-resolved spectroscopy to investigate its electronic relaxation mechanism and probe the effect of atom substitution on the relaxation mechanism in polar protic and polar aprotic solvents. Mapping of the excited state potential energy surfaces is performed, from which the critical points are optimized by using the state-of-art Extended Multi-State Complete Active Space Second-Order Perturbation Theory. It is demonstrated that excitation to the lowest energy ¹ ππ* state of ^5N7C G results in complex dynamics leading to ca. 10 to 30-fold slower relaxation (depending on solvent) compared to guanine. A significant conformational change occurs at the S₁ minimum, resulting in a 10-fold greater fluorescence quantum yield compared to guanine. The fluorescence quantum yield and S₁ decay lifetime increase going from water to acetonitrile to propanol. The solvent-dependent results are supported by the quantum chemical calculations showing an increase in the energy barrier between the S₁ minimum and the S₁ /S₀ conical intersection going from water to propanol. The longer lifetimes might make ^5N7C G more photochemical active to adjacent nucleobases than guanine or other nucleobases within DNA.

Free-Energy Landscape of the SN2 Reaction CH3Br + Cl- → CH3Cl + Br- in Different Liquid Environments

This work describes in detail the reaction path of the well-known S_N2 reaction CH₃Br + Cl^- → CH₃Cl + Br^-, whose reaction rate has a huge variation with the solvent in the gas phase and in protic and aprotic liquid environments. We employed the ASEC-FEG method to optimize for minima (reactants and products) and saddle points (transition states) in the in-solution free-energy hypersurface. The method takes atomistic details of the solvent into account. A polarizable continuum model (PCM) has also been employed for comparison. The most perceptive structural changes are noted in aqueous solution by using the ASEC-FEG approach. The activation energies in all solvents, estimated by means of free-energy perturbation calculations, are in good agreement with the experimental data. The total solute-solvent hydrogen bonds play an important role in the increased barrier height observed in water and are therefore crucial to explain the huge decrease in the kinetic constant. It is also found that the hydration shell around the ions breaks itself spontaneously to accommodate the molecule, thus forming minimum energy complexes.

Quantifying Fluctuations of Average Solvent Environments for Embedding Calculations

The viability and effectiveness of replacing an ensemble of embedded solute calculations by a single calculation using an average description of the solvent environment are evaluated. This work explores the fluctuations of the average description of the system obtained in two ways: from calculations on an ensemble of geometries and from an average environment constructed from the same ensemble. To this end, classical molecular dynamics simulations of a rigid acetone solute in SPCE water are performed in order to generate an ensemble of solvent environments. From this ensemble of solvent configurations, a number of different approaches for constructing an average solvent environment are employed. We perform a thorough numerical analysis of the fluctuations of the electrostatic potential experienced by the solute, as well as the resulting fluctuations of the solute's electronic density, measured through its dipole moment and fitted atomic point charges. At the same time, we inspect the accuracy of the methods used to construct average environments. Finally, the proposed method for generating the embedding potential from an average environment density is applied to estimate the solvatochromic shift of the first excitation of acetone. In order to account for quantum confinement effects, which may be important in certain cases, the fluctuations in the shift due to the interaction with the solvent are evaluated using frozen-density-embedding theory. Our results demonstrate that, for normally distributed environments, the constructed average environment is a reasonably good representation of a fluctuating molecular solvent environment. We then provide guidance for future comparisons between these theoretical treatments of solute/solvent systems to experimental measurements.

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning's cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.

Free Energy Computation for an Isomerizing Chromophore in a Molecular Cavity via the Average Solvent Electrostatic Configuration Model: Applications in Rhodopsin and Rhodopsin-Mimicking Systems

We present a novel technique for computing the free energy differences between two chromophore "isomers" hosted in a molecular environment (a generalized solvent). Such an environment may range from a relatively rigid protein cavity to a flexible solvent environment. The technique is characterized by the application of the previously reported "average electrostatic solvent configuration" method, and it is based on the idea of using the free energy perturbation theory along with a chromophore annihilation procedure in thermodynamic cycle calculations. The method is benchmarked by computing the ground-state room-temperature relative stabilities between (i) the cis and trans isomers of prototypal animal and microbial rhodopsins and (ii) the analogue isomers of a rhodopsin-like light-driven molecular switch in methanol. Furthermore, we show that the same technology can be used to estimate the activation free energy for the thermal isomerization of systems i-ii by replacing one isomer with a transition state. The results show that the computed relative stability and isomerization barrier magnitudes for the selected systems are in line with the available experimental observation in spite of their widely diverse complexity.

Two-dimensional self-assembled nanostructures of nucleobases and their related derivatives on Au(111)

The construction of two-dimensional (2D) self-assembled nanostructures has been one of the considerably interesting areas of on-surface chemistry in the past few decades, and has benefited from the rapid development and improvement of scanning probe microscopy techniques. In this research field, many attempts have been made in the controllable fabrication of well-ordered and multifunctional surface nanostructures, which attracted interest because of the prospect for artificial design of functional molecular nanodevices. DNA and RNA are considered to be programmable self-assembly systems and it is possible to use their base sequences to encode instructions for assembly in a predetermined fashion at the nanometer scale. As important constituents of nucleic acids, nucleobases, with intrinsic functional groups for hydrogen bonding, coordination bonding, and electrostatic interactions, can be employed as a potential system for the versatile construction of various biomolecular nanostructures, which may be used to structure the self-assembly of DNA-based artificial molecular constructions and play an important role in novel biosensors based on surface functionalization. In this article, we will review the recent progress of on-surface self-assembly of nucleobases and their derivatives together with different reactants (e.g., metals, halogens, salts and water), and as a result, various 2D surface nanostructures are summarized.

Engineering the vibrational coherence of vision into a synthetic molecular device

The light-induced double-bond isomerization of the visual pigment rhodopsin operates a molecular-level optomechanical energy transduction, which triggers a crucial protein structure change. In fact, rhodopsin isomerization occurs according to a unique, ultrafast mechanism that preserves mode-specific vibrational coherence all the way from the reactant excited state to the primary photoproduct ground state. The engineering of such an energy-funnelling function in synthetic compounds would pave the way towards biomimetic molecular machines capable of achieving optimum light-to-mechanical energy conversion. Here we use resonance and off-resonance vibrational coherence spectroscopy to demonstrate that a rhodopsin-like isomerization operates in a biomimetic molecular switch in solution. Furthermore, by using quantum chemical simulations, we show why the observed coherent nuclear motion critically depends on minor chemical modifications capable to induce specific geometric and electronic effects. This finding provides a strategy for engineering vibrationally coherent motions in other synthetic systems.

The ultrafast, vibrationally coherent photoisomerization of rhodopsin is a model of efficient photomechanical energy conversion at the molecular scale. Here, the authors demonstrate a similar photoreaction in synthetic compounds, unraveling the underlying mechanism and discussing its implications.

Zwitterionization of glycine in water environment: Stabilization mechanism and NMR spectral signatures

At physiological conditions, myriads of biomolecules (e.g., amino acids, peptides, and proteins) exist predominantly in the zwitterionic structural form and their biological functions will result in these conditions. However these geometrical structures are inaccessible energetically in the gas phase, and at this point, stabilization of amino-acids in physiological conditions is still under debate. In this paper, the electronic properties of a glycine molecule in the liquid environment were studied by performing a relaxation of the glycine geometry in liquid water using the free energy gradient method combined with a sequential quantum mechanics/molecular mechanics approach. A series of Monte Carlo Metropolis simulations of the glycine molecule embedded in liquid water, followed by only a quantum mechanical calculation in each of them were carried out. Both the local and global liquid environments were emphasized to obtain nuclear magnetic resonance (NMR) parameters for the glycine molecule in liquid water. The results of the equilibrium structure in solution and the systematic study of the hydrogen bonds were used to discard the direct proton transfer from the carboxyl group to the ammonium group of the glycine molecule in water solution. The calculations of the Density Functional Theory (DFT) were performed to study the polarization of the solvent in the parameters of nuclear magnetic resonance of the glycine molecule in liquid water. DFT calculations predicted isotropic chemical changes on the H, C, N, and O atoms of glycine in liquid water solution which agree with the available experimental data.

Mixed Quantum/Classical Method for Nonadiabatic Quantum Dynamics in Explicit Solvent Models: The ππ*/nπ* Decay of Thymine in Water as a Test Case

We present a novel mixed quantum classical dynamical method to include solvent effects on internal conversion (IC) processes. All the solute degrees of freedom are represented by a wavepacket moving according to nonadiabatic quantum dynamics, while the motion of an explicit solvent model is described by an ensemble of classical trajectories. The mutual coupling of the solute and solvent dynamics is included within a mean-field framework and the quantum and classical equations of motions are solved simultaneously. As a test case we apply our method to the ultrafast ππ* → nπ* decay of thymine in water. Solvent dynamical response modifies IC yield already on the 50 fs time scale. This effect is due to water librational motions that stabilize the most populated state. Pure static disorder, that is, the existence of different solvent configurations when photoexcitation takes place, also has a remarkable impact on the dynamics.

Electronic structure and absorption spectra of fluorescent nucleoside analogues

This work presents a systematic investigation of the electronic and conformational properties of five new fluorescent nucleobases belonging to the alphabet based on the isothiazole[4,3-d]pyrimidine molecule, very recently synthesized. This is of particular importance in the characterization of the main electronic aspects of these fluorescent nucleosides. The solvent effects of 1,4-dioxane and water were included combining the Sequential Monte Carlo/CASPT2 and the Free Energy Gradient (FEG) methods. For comparison, the Polarizable Continuum method was also used. The geometries of all compounds were optimized in solvent with the largest effects observed in water using the average solvent electrostatic configuration (ASEC) and the FEG approaches. Statistical analysis of the solute-solvent hydrogen bonds is performed and their effect on the absorption spectra analyzed. The dipole moments were calculated and the value obtained from the ASEC-FEG method in water follows the same trend as the natural canonical bases (adenine → uracil → guanine → cytosine). The theoretical results for the absorption spectra obtained from CASPT2(18,13) calculations using the geometries obtained with the ASEC-FEG procedure are in very good agreement with the experimental data. A detailed elucidation of the main aspects of the absorption spectra of the five new fluorescent nucleoside analogues is successfully attempted.

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.