Strikingly Different Effects of Hydrogen Bonding on the Photodynamics of Individual Nucleobases in DNA: Comparison of Guanine and Cytosine

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

Point-of-Care Determination of Acetaminophen Levels with Multi-Hydrogen Bond Manipulated Single-Molecule Recognition (eMuHSiR)

This work aims to face the challenge of monitoring small molecule drugs accurately and rapidly for point-of-care (POC) diagnosis in current clinical settings. Overdose of acetaminophen (AP), a commonly used over the counter (OTC) analgesic drug, has been determined to be a major cause of acute liver failure in the US and the UK. However, there is no rapid and accurate detection method available for this drug in the emergency room. The present study examined an AP sensing strategy that relies on a previously unexplored strong interaction between AP and the arginine (Arg) molecule. It was found that as many as 4 hydrogen bonds can be formed between one Arg molecule and one AP molecule. By taking advantages of this structural selectivity and high tenability of hydrogen bonds, Arg, immobilized on a graphene surface via electrostatic interactions, was utilized to structurally capture AP. Interestingly, bonded AP still remained the perfect electrochemical activities. The extent of Arg-AP bonds was quantified using a newly designed electrochemical (EC) sensor. To verify the feasibility of this novel assay, based on multihydrogen bond manipulated single-molecule recognition (eMuHSiR), both pharmaceutical and serum sample were examined. In commercial tablet measurement, no significant difference was seen between the results of eMuHSiR and other standard methods. For measuring AP concentration in the mice blood, the substances in serum, such as sugars and fats, would not bring any interference to the eMuHSiR in a wide concentration range. This eMuHSiR method opens the way for future development of small molecule detection for the POC testing.

An ab initio multiple cloning approach for the simulation of photoinduced dynamics in conjugated molecules

Multidimensional wave function: a superposition of Gaussian coherent states guided by Ehrenfest trajectories suited to clone and swap their electronic amplitudes.

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.

Photophysical deactivation pathways in adenine oligonucleotides

In this work we study deactivation processes in adenine oligomers after absorption of UV radiation using Quantum Mechanics combined with Molecular Mechanics (QM/MM). Correlated electronic structure methods appropriate for describing the excited states are used to describe a π-stacked dimer of adenine bases incorporated into (dA)20(dT)20. The results of these calculations reveal three different types of excited state minima which play a role in deactivation processes. Within this set of minima there are minima where the excited state is localized on one adenine (monomer-like) as well as minima where the excited state is delocalized on two adenines, forming different types of excimers and bonded excimers of varying but inter-related character. The proximity of their energies reveals that the minima can decay into one another along a flat potential energy surface dependent on the interbase separation. Additionally, analysis of the emissive energies and other physical properties, including theoretical anisotropy calculations, and comparison with fluorescence experiments, provides evidence that excimers play an important role in long-lived signals in adenine oligonucleotides while the subpicosecond decay is attributed to monomer-like minima. The necessity for a close approach of the nucleobases reveals that the deactivation mechanism is tied to macro-molecular motion.

A novel method developed for estimating mineralization efficiencies and its application in PC and PEC degradations of large molecule biological compounds with unknown chemical formula

A new method to estimate the photocatalytic (PC) and photoelectrocatalytic (PEC) mineralization efficiencies of large molecule biological compounds with unknown chemical formula in water was firstly developed and experimentally validated. The method employed chemical oxidation under the standard dichromate chemical oxygen demand (COD) conditions to obtain QCOD values of model compounds with unknown chemical formula. The measured QCOD values were used as the reference to replace QCOD values of model compounds for calculation of the mineralization efficiencies (in %) by assuming the obtained QCOD values are the measure of the theoretical charge required for the complete mineralization of organic pollutants. Total organic carbon (TOC) was also employed as a reference to confirm the mineralization capacity of dichromate chemical oxidation. The developed method was applied to determine the degradation extent of model compounds, such as bovine serum albumin (BSA), lecithin and bacterial DNA, by PC and PEC. Incomplete PC mineralization of all large molecule biological compounds was observed, especially for BSA. But the introduction of electrochemical technique into a PC oxidation process could profoundly improve the mineralization efficiencies of model compounds. PEC mineralization efficiencies of bacterial DNA was the highest, while that of lecithin was the lowest. Overall, PEC degradation method was found to be much effective than PC method for all large molecule biological compounds investigated, with PEC/PC mineralization ratios followed an order of BSA > lecithin > DNA.

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.

Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations

Sequential vs. concerted S₁ relaxation pathways.

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.

Role of excitonic coupling and charge-transfer states in the absorption and CD spectra of adenine-based oligonucleotides investigated through QM/MM simulations

In this work, we study the photophysical properties of an adenine-based oligonucleotide using an ensemble of about 200 configurations obtained from molecular dynamics simulations. Specifically, a QM/MM approach is used to obtain the excited-state energies and properties of (dA)20(dT)20 with a dimer of π-stacked adenine bases included in the quantum region. The absorption and circular dichroism spectra are computed and analyzed using the algebraic diagrammatic construction through second order level of theory method (ADC(2)) combined with classical mechanics. We find that the experimentally observed red-shifted shoulder in the absorption spectrum is due to excitonic interactions, while charge-transfer states are present within the absorption band at the higher-energy end of the spectrum. More importantly, low-energy states with charge-transfer mixing exist, which could lead to excimers and bonded excimers. These observations suggest that mixing between charge-transfer and excitonic states plays an important role in the photophysics of oligonucleotides. They also highlight the importance of taking into account the conformational flexibility of the oligonucleotide when investigating photophysical properties.

Electronic excitation and structural relaxation of the adenine dinucleotide in gas phase and solution

The excited states and potential surfaces of the adenine dinucleotide are analyzed in gas phase and in solution using a correlated ab initio methodology in a QM/MM framework. In agreement with previous studies, a rather flat S1 surface with a number of minima of different character is found. Specifically, our results suggest that exciplexes with remarkably short intermolecular separation down to ~2.0 Å are formed. A detailed analysis shows that due to strong orbital interactions their character differs significantly from any states present in the Franck-Condon region. The lowest S1 energy minimum is a ππ* exciplex with only a small amount of charge transfer. It possesses appreciable oscillator strength with a polarization almost perpendicular to the planes of the two adenine molecules.

Synchronous concerted multiple-body photodissociation of oxalyl chloride explored by ab initio-based dynamics simulations

Photo-induced multiple body dissociation is of fundamental interest in chemistry and physics. A description of the mechanism associated with n-body (n ≥ 3) photodissociation has proven to be an intriguing and yet challenging issue in the field of chemical dynamics. Oxalyl chloride, (ClCO)2, is the sole molecule reported up to date that can undergo four-body dissociation following absorption of a single UV photon, with a rich history of mechanistic debate. In the present work, the combined electronic structure calculations and dynamics simulations have been performed at the advanced level, which provides convincing evidence for resolving the mechanistic debate. More importantly, synchronous and asynchronous concertedness were explored for the first time for the (ClCO)2 photodissociation, which is based on the simulated time constants for the C-C and C-Cl bond fissions. Upon photoexcitation of (ClCO)2 to the S1 state, the adiabatic C-C or C-Cl fission takes place with little possibility. The four-body dissociation to 2Cl((2)P) and 2CO((1)Σ) was determined to a dominant channel with its branch of ∼0.7, while the three-body dissociation to ClCO((2)A(')) + CO((1)Σ) + Cl((2)P) was predicted to play a minor role in the (ClCO)2 photodissociation at 193 nm. Both the four-body and three-body dissociations are non-adiabatic processes, which proceed in a synchronous concerted way as a result of the S1 → S0 internal conversion. There is a little possibility for two-body dissociation to occur in the S0 and S1 states.

A QM/MM MD insight into photodynamics of hypoxanthine: distinct nonadiabatic decay behaviors between keto-N7H and keto-N9H tautomers in aqueous solution

QM/MM MD simulations reveal different H-bonding networks around DNA base analogues, and the π-electron H-bond in the solvated keto-N7H may facilitate its S₁ → S₀ nonadiabatic decay.

Ultrafast electronic deactivation dynamics of the inosine dimer--a model case for H-bonded purine bases

The structural properties and ultrafast electronic deactivation dynamics of the inosine dimer in CHCl3 have been investigated by two-dimensional (1)H NMR and static FTIR spectroscopy and by femtosecond time-resolved transient absorption spectroscopy, respectively. The (1)H NMR and IR spectra show the formation of a well-defined, symmetric dimer with an association equilibrium constant of KI·I = 690 ± 100 M(-1). The excited-state dynamics after photoexcitation at λpump = 260 nm monitored by ultrafast absorption spectroscopy show great similarity with those of the monomer inosine in an aqueous solution and are governed by a decay time of τ = 90 ± 10 fs, which is one of the shortest electronic lifetimes of all nucleobases and nucleobase dimers studied so far. On the basis of these observations, the inosine dimer is expected to follow a similar relaxation pathway as the monomer, involving an out-of-plane deformation of the six-membered ring. The importance of the C(2) position for the electronic deactivation of hypoxanthine and guanine is discussed. The obtained well-determined structure and straightforward dynamics qualify the inosine dimer as an excellent reference case for more complicated systems such as the G·G dimer and the G·C and A·T Watson-Crick pairs.

Nonadiabatic photodynamics of a retinal model in polar and nonpolar environment

The nonadiabatic photodynamics of the all-trans-2,4-pentadiene-iminium cation (protonated Schiff base 3, PSB3) and the all-trans-3-methyl-2,4-pentadiene-iminium cation (MePSB3) were investigated in the gas phase and in polar (aqueous) and nonpolar (n-hexane) solutions by means of surface hopping using a multireference configuration-interaction (MRCI) quantum mechanical/molecular mechanics (QM/MM) level. Spectra, lifetimes for radiationless deactivation to the ground state, and structural and electronic parameters are compared. A strong influence of the polar solvent on the location of the crossing seam, in particular in the bond length alternation (BLA) coordinate, is found. Additionally, inclusion of the polar solvent changes the orientation of the intersection cone from sloped in the gas phase to peaked, thus enhancing considerably its efficiency for deactivation of the molecular system to the ground state. These factors cause, especially for MePSB3, a substantial decrease in the lifetime of the excited state despite the steric inhibition by the solvent.