Abstract:Hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations on aqueous solvated uridine show that TD‐DFT and CASPT2 reference methods predict very similar topologies of the potential energy surface of the lowest ππ* excited state. The discrepant behavior, reported in previous publications, between the two methods in describing DNA in polar environments is addressable to the different relative energetics of the ππ* and nπ* states. We demonstrate that by adopting an appropriate scheme for modeling the environment, the considerably cheaper TD‐DFT method could be positioned as the standout protocol to model photoinduced dynamics in DNA monomers and their aggregates. 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 by nπ* 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.

QM/MM-MD analysis of dissociation of Ag+ and H+ mediated cytidines

Processes of Dna Condensation Induced by Multivalent Cations: Approximate Annealing Experiments and Molecular Dynamics Simulations

Quantum mechanics/molecular mechanics studies on the excited-state decay mechanisms of cytidine aza-analogues: 5-azacytidine and 2ˊ-deoxy-5-azacytidine in aqueous solution

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions.

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

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes.

Revealing the Excited-state Dynamics of Cytidine and the Role of Excited-state Proton Transfer Process

Elucidating the stabilization mechanism of a K + -depleted c-MYC DNA G-quadruplex in hydrophobic imidazolium-based ionic liquids using spectroscopy coupled with molecular dynamics simulations

Unravelling Mg$^{2+}$-RNA binding with atomistic molecular dynamics

Semi-empirical and Linear-Scaling DFT Methods to Characterize duplex DNA and G-quadruplexes in Presence of Interacting Small Molecules

QM/MM(ABEEM) Study on the Ligand Substitution Processes of Ruthenium(III) Complex NAMI‐A

Negative Electron Transfer Collision-Induced Dissociation of G-Quadruplexes: Uncovering the Guanine Radical Anion Loss Pathway

Quantum mechanical calculations related to ionization and charge transfer in DNA

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis.

Excited States in Single-Stranded and i-Motif DNA with Silver Ions

Excess-electron Capture and Energy Transfer to Bulk Water for Aqueous DNA Nucleotide

Is 1-methylcytosine a faithful model for ultrafast deactivation dynamics of cytosine nucleosides in solution?

Differential Mg2+ Deposition on DNA Holliday Junctions Dictates the Rate and Stability of Conformational Exchange

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

On the supramolecular interactions into a pH- and Metal-Actuated Molecular Shuttle: some insights from QTAIM modeling