Photoinduced water-chromophore electron transfer causes formation of guanosine photodamage.

UV-induced photolysis of aqueous guanine nucleosides produces 8-oxo-guanine and Fapy-guanine, which can induce various types of cellular malfunction. The mechanistic rationale underlying photodestructive processes of guanine nucleosides is still largely obscure. Here, we employ accurate quantum chemical calculations and demonstrate that an excited-state non-bonding interaction of guanosine and a water molecule facilitates the electron-driven proton transfer process from water to the chromophore fragment. This subsequently allows for the formation of a crucial intermediate, namely guanosine photohydrate. Further (photo)chemical reactions of this intermediate lead to the known products of guanine photodamage.

Light-Induced Modulation of Chiral Functions in G-Quadruplex-Photochrome Systems.

The design of artificially engineered chiral structures has received much attention, but the implementation of dynamic functions to modulate the chiroptical response of the systems is less explored. Here, we present a light-responsive G-quadruplex (G4)-based assembly in which chirality enrichment is induced, tuned, and fueled by molecular switches. In particular, the mirror-image dependence on photoactivated azo molecules, undergoing trans-to-cis isomerization, shows chiral recognition effects on the inherent flexibility and conformational diversity of DNA G4s having distinct handedness (right- and left-handed). Through a detailed experimental and computational analysis, we bring compelling evidence on the binding mode of the photochromes on G4s, and we rationalize the origin of the chirality effect that is associated with the complexation event.

Ribose Alters the Photochemical Properties of the Nucleobase in Thionated Nucleosides.

Substitution of exocyclic oxygen with sulfur was shown to substantially influence the properties of RNA/DNA bases, which are crucial for prebiotic chemistry and photodynamic therapies. Upon UV irradiation, thionucleobases were shown to efficiently populate triplet excited states and can be involved in characteristic photochemistry or generation of singlet oxygen. Here, we show that the photochemistry of a thionucleobase can be considerably modified in a nucleoside, that is, by the presence of ribose. Our transient absorption spectroscopy experiments demonstrate that thiocytosine exhibits 5 times longer excited-state lifetime and different excited-state absorption features than thiocytidine. On the basis of accurate quantum chemical simulations, we assign these differences to the dominant population of a shorter-lived triplet nπ* state in the nucleoside and longer-lived triplet ππ* states in the nucleobase. This explains the distinctive photoanomerziation of thiocytidine and indicates that the nucleoside will be a less efficient phototherapeutic agent with regard to singlet oxygen generation.

An effective potential for Frenkel excitons.

Excitation energy transfer (EET) is a ubiquitous process in life and materials sciences. Here, a new and computationally efficient method of evaluating the electronic EET couplings between interacting chromophores is introduced that is valid in a wide range of intermolecular distances. The proposed approach is based on the effective elimination of electron repulsion integrals from the excitonic Hamiltonian matrix elements via the density-fitting approach and distributed multipole approximation. The excitonic Hamiltonian represented in a basis including charge transfer (CT) states is re-cast in terms of the effective one-electron potential functions (EOPs) and adapted into the effective fragment parameter (EFP) framework. Calculations for model systems indicate that the speedup of at least three orders of magnitude, as compared to the state-of-the-art methods, can be achieved while maintaining the accuracy of the EET couplings even at short intermolecular distances.

Ultrafast nonradiative deactivation of photoexcited 8-oxo-hypoxanthine: a nonadiabatic molecular dynamics study.

In the scientific endeavor to understand the chemical origins of life, the photochemistry of the smallest life building blocks, nucleobases, has been a constant object of focus and intense research. Here, we report the results of the first theoretical study on the photo-properties of an 8-oxo-hypoxanthine molecule, the chromophore of 8-oxo-inosine, which is relevant to the recently proposed, prebiotically plausible synthetic routes to the formation of purine- and pyrimidine-nucleotides. With ab initio and semi-empirical OM2/MRCI quantum-chemistry calculations, we predict a strong photostability of the 8-oxo-hypoxanthine system and see the origin of this effect in ultrafast nonradiative relaxation through puckering of the 6-membered heterocyclic ring.

Ab initio effective one-electron potential operators: Applications for charge-transfer energy in effective fragment potentials.

The concept of effective one-electron potentials (EOPs) has proven to be extremely useful in efficient description of electronic structure of chemical systems, especially extended molecular aggregates such as interacting molecules in condensed phases. Here, a general method for EOP-based elimination of electron repulsion integrals is presented, that is tuned toward the fragment-based calculation methodologies such as the second generation of the effective fragment potentials (EFP2) method. Two general types of the EOP operator matrix elements are distinguished and treated either via the distributed multipole expansion or the extended density fitting (DF) schemes developed in this work. The EOP technique is then applied to reduce the high computational costs of the effective fragment charge-transfer (CT) terms being the bottleneck of EFP2 potentials. The alternative EOP-based CT energy model is proposed, derived within the framework of intermolecular perturbation theory with Hartree-Fock noninteracting reference wavefunctions, compatible with the original EFP2 formulation. It is found that the computational cost of the EFP2 total interaction energy calculation can be reduced by up to 38 times when using the EOP-based formulation of CT energy, as compared to the original EFP2 scheme, without compromising the accuracy for a wide range of weakly interacting neutral and ionic molecular fragments. The proposed model can thus be used routinely within the EFP2 framework.

Selective prebiotic formation of RNA pyrimidine and DNA purine nucleosides.

The nature of the first genetic polymer is the subject of major debate1. Although the 'RNA world' theory suggests that RNA was the first replicable information carrier of the prebiotic era-that is, prior to the dawn of life2,3-other evidence implies that life may have started with a heterogeneous nucleic acid genetic system that included both RNA and DNA4. Such a theory streamlines the eventual 'genetic takeover' of homogeneous DNA from RNA as the principal information-storage molecule, but requires a selective abiotic synthesis of both RNA and DNA building blocks in the same local primordial geochemical scenario. Here we demonstrate a high-yielding, completely stereo-, regio- and furanosyl-selective prebiotic synthesis of the purine deoxyribonucleosides: deoxyadenosine and deoxyinosine. Our synthesis uses key intermediates in the prebiotic synthesis of the canonical pyrimidine ribonucleosides (cytidine and uridine), and we show that, once generated, the pyrimidines persist throughout the synthesis of the purine deoxyribonucleosides, leading to a mixture of deoxyadenosine, deoxyinosine, cytidine and uridine. These results support the notion that purine deoxyribonucleosides and pyrimidine ribonucleosides may have coexisted before the emergence of life5.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. II. Halogen-bonded systems.

Following our study on hydrogen-bonded (HB) complexes [Phys. Chem. Chem. Phys., 2018, 20, 19841], the physical nature of interaction-induced (non)linear optical properties of another important class of molecular complexes, namely halogen-bonded (XB) systems, was analyzed in this study. The excess electronic and nuclear relaxation (hyper)polarizabilities of nine representative XB complexes covering a wide range of halogen-bond strengths were computed. The partitioning of the excess properties into individual interaction-energy components (electrostatic, exchange, induction, dispersion) was performed by using the variational-perturbational energy decomposition scheme at the MP2/aug-cc-pVTZ level of theory and further supported by calculations with the SCS-MP2 method. In the case of the electronic interaction-induced properties, the physical composition of Δαel and Δγel was found to be very similar for the two types of bonding, despite the different nature of the binding. For Δßel, the XB complexes exhibit a more systematic interplay of interaction-energy contributions compared to the HB systems studied in the previous work. Our analysis revealed that the patterns of interaction-energy contributions to the interaction-induced nuclear-relaxation contributions to the linear polarizability and the first hyperpolarizability are very similar. For both properties the exchange repulsion term is canceled out by the electrostatic and delocalization terms. The physical composition of these contributions is analogous to those observed for the HB complexes.

UV-induced hydrogen transfer in DNA base pairs promoted by dark nπ* states.

Dark nπ* states were shown to have substantial contribution to the destructive photochemistry of pyrimidine nucleobases. Based on quantum-chemical calculations, we demonstrate that the characteristic hydrogen bonding pattern of the GC base pair could facilitate the formation of a wobble excited-state charge-transfer complex. This entails a barrierless electron-driven proton transfer (EDPT) process which enables damageless photodeactivation of the base pair. These photostabilizing properties are retained even when guanine is exchanged to hypoxanthine. The inaccessibility of this process in the AT base pair sheds further light on the reasons why cytosine is less susceptible to the formation of photodimers in double-stranded DNA.

Photostability of oxazoline RNA-precursors in UV-rich prebiotic environments.

Pentose aminooxazolines and oxazolidinone thiones are considered as the key precursors which could have enabled the formation of RNA nucleotides under the conditions of early Earth. UV-irradiation experiments and quantum-chemical calculations demonstrate that these compounds are remarkably photostable and could accumulate over long periods of time in UV-rich prebiotic environments to undergo stereoisomeric purification.

Electron-driven proton transfer enables nonradiative photodeactivation in microhydrated 2-aminoimidazole.

2-Aminoimidazole (2-AIM) was proposed as a plausible nucleotide activating group in a nonenzymatic copying and polymerization of short RNA sequences under prebiotically plausible conditions. One of the key selection factors controlling the lifespan and importance of organic molecules on early Earth was ultraviolet radiation from the young Sun. Therefore, to assess the suitability of 2-AIM for prebiotic chemistry, we performed non-adiabatic molecular dynamics simulations and static explorations of potential energy surfaces of the photoexcited 2-AIM-(H2O)5 model system by means of the algebraic diagrammatic construction method to the second order [ADC(2)]. Our quantum mechanical simulations demonstrate that 1πσ* excited states play a crucial role in the radiationless deactivation of the UV-excited 2-AIM-(H2O)5 system. More precisely, electron-driven proton transfer (EDPT) along water wires is the only photorelaxation pathway leading to the formation of 1πσ*/S0 conical intersections. The availability of this mechanism and the lack of destructive photochemistry indicate that microhydrated 2-AIM is characterized by substantial photostability and resistance to prolonged UV irradiation.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. I. Hydrogen-bonded systems.

Understanding the effects of different fundamental intermolecular interactions on nonlinear optical properties is crucial for proposing efficient strategies to obtain new materials with tailored properties. In this study, we computed the electronic and vibrational (hyper)polarizabilities of ten hydrogen-bonded molecular complexes employing the MP2, CCSD and CCSD(T) methods combined with the aug-cc-pVTZ basis set. The vibrational contributions to hyperpolarizabilities included nuclear-relaxation anharmonic corrections. The effect of intermolecular interactions was analyzed in terms of excess properties, which are defined as the difference between a property of the complex and the net properties of the noninteracting subsystems. Considering systems covering a wide range of hydrogen bond strengths, the electronic and vibrational excess (hyper)polarizabilities were decomposed into different interaction energy contributions (electrostatic, exchange, induction and dispersion). This systematic study, the very first of this kind, revealed that the physical origins of the electronic and vibrational excess properties are completely different. In the case of vibrational contributions, the decomposition pattern is very similar for the polarizability and first and second hyperpolarizabilities. The exchange contributions to excess vibrational properties are the largest and they have an opposite sign to the electrostatic, induction and dispersion terms. On the other hand, no general patterns can be established for the electronic excess properties.

Water-chromophore electron transfer determines the photochemistry of cytosine and cytidine.

Many of the UV-induced phenomena observed experimentally for aqueous cytidine were lacking the mechanistic interpretation for decades. These processes include the substantial population of the puzzling long-lived dark state, photohydration, cytidine to uridine conversion and oxazolidinone formation. Here, we present quantum-chemical simulations of excited-state spectra and potential energy surfaces of N1-methylcytosine clustered with two water molecules using the second-order approximate coupled cluster (CC2), complete active space with second-order perturbation theory (CASPT2), and multireference configuration interaction with single and double excitation (MR-CISD) methods. We argue that the assignment of the long-lived dark state to a singlet nπ* excitation involving water-chromophore electron transfer might serve as an explanation for the numerous experimental observations. While our simulated spectra for the state are in excellent agreement with experimentally acquired data, the electron-driven proton transfer process occurring on the surface may initiate the subsequent damage in the vibrationally hot ground state of the chromophore.

A prebiotically plausible synthesis of pyrimidine ß-ribonucleosides and their phosphate derivatives involving photoanomerization.

Previous research has identified ribose aminooxazoline as a potential intermediate in the prebiotic synthesis of the pyrimidine nucleotides with remarkable properties. It crystallizes spontaneously from reaction mixtures, with an enhanced enantiomeric excess if initially enantioenriched, which suggests that reservoirs of this compound might have accumulated on the early Earth in an optically pure form. Ribose aminooxazoline can be converted efficiently into α-ribocytidine by way of 2,2'-anhydroribocytidine, although anomerization to ß-ribocytidine by ultraviolet irradiation is extremely inefficient. Our previous work demonstrated the synthesis of pyrimidine ß-ribonucleotides, but at the cost of ignoring ribose aminooxazoline, using arabinose aminooxazoline instead. Here we describe a long-sought route through ribose aminooxazoline to the pyrimidine ß-ribonucleosides and their phosphate derivatives that involves an extraordinarily efficient photoanomerization of α-2-thioribocytidine. In addition to the canonical nucleosides, our synthesis accesses ß-2-thioribouridine, a modified nucleoside found in transfer RNA that enables both faster and more-accurate nucleic acid template-copying chemistry.

Photorelaxation of imidazole and adenine via electron-driven proton transfer along H2O wires.

Photochemically created πσ* states were classified among the most prominent factors determining the ultrafast radiationless deactivation and photostability of many biomolecular building blocks. In the past two decades, the gas phase photochemistry of πσ* excitations was extensively investigated and was attributed to N-H and O-H bond fission processes. However, complete understanding of the complex photorelaxation pathways of πσ* states in the aqueous environment was very challenging, owing to the direct participation of solvent molecules in the excited-state deactivation. Here, we present non-adiabatic molecular dynamics simulations and potential energy surface calculations of the photoexcited imidazole-(H2O)5 cluster using the algebraic diagrammatic construction method to the second-order [ADC(2)]. We show that electron driven proton transfer (EDPT) along a wire of at least two water molecules may lead to the formation of a πσ*/S0 state crossing, similarly to what we suggested for 2-aminooxazole. We expand on our previous findings by direct comparison of the imidazole-(H2O)5 cluster to non-adiabatic molecular dynamics simulations of imidazole in the gas phase, which reveal that the presence of water molecules extends the overall excited-state lifetime of the chromophore. To embed the results in a biological context, we provide calculations of potential energy surface cuts for the analogous photorelaxation mechanism present in adenine, which contains an imidazole ring in its structure.

On the physical origins of interaction-induced vibrational (hyper)polarizabilities.

This paper presents the results of a pioneering exploration of the physical origins of vibrational contributions to the interaction-induced electric properties of molecular complexes. In order to analyze the excess nuclear relaxation (hyper)polarizabilities, a new scheme was proposed which relies on the computationally efficient Bishop-Hasan-Kirtman method for determining the nuclear relaxation contributions to electric properties. The extension presented herein is general and can be used with any interaction-energy partitioning method. As an example, in this study we employed the variational-perturbational interaction-energy decomposition scheme (at the MP2/aug-cc-pVQZ level) and the extended transition state method by employing three exchange-correlation functionals (BLYP, LC-BLYP, and LC-BLYP-dDsC) to study the excess properties of the HCN dimer. It was observed that the first-order electrostatic contribution to the excess nuclear relaxation polarizability cancels with the negative exchange repulsion term out to a large extent, resulting in a positive value of Δα(nr) due to the contributions from the delocalization and the dispersion terms. In the case of the excess nuclear relaxation first hyperpolarizability, the pattern of interaction contributions is very similar to that for Δα(nr), both in terms of their sign as well as relative magnitude. Finally, our results show that the LC-BLYP and LC-BLYP-dDsC functionals, which yield smaller values of the orbital relaxation term than BLYP, are more successful in predicting excess properties.

Ultrafast excited-state dynamics of isocytosine.

The alternative nucleobase isocytosine has long been considered as a plausible component of hypothetical primordial informational polymers. To examine this hypothesis we investigated the excited-state dynamics of the two most abundant forms of isocytosine in the gas phase (keto and enol). Our surface-hopping nonadiabatic molecular dynamics simulations employing the algebraic diagrammatic construction to the second order [ADC(2)] method for the electronic structure calculations suggest that both tautomers undergo efficient radiationless deactivation to the electronic ground state with time constants which amount to τketo = 182 fs and τenol = 533 fs. The dominant photorelaxation pathways correspond to ring-puckering (ππ* surface) and C[double bond, length as m-dash]O stretching/N-H tilting (nπ* surface) for the enol and keto forms respectively. Based on these findings, we infer that isocytosine is a relatively photostable compound in the gas phase and in these terms resembles biologically relevant nucleobases. The estimated S1 [radiolysis arrow - arrow with voltage kink] T1 intersystem crossing rate constant of 8.02 × 10(10) s(-1) suggests that triplet states might also play an important role in the overall excited-state dynamics of the keto tautomer. The reliability of ADC(2)-based surface-hopping molecular dynamics simulations was tested against multireference quantum-chemical calculations and the potential limitations of the employed ADC(2) approach are briefly discussed.

Prebiotic synthesis of nucleic acids and their building blocks at the atomic level - merging models and mechanisms from advanced computations and experiments.

The origin of life on Earth is one of the most fascinating questions of contemporary science. Extensive research in the past decades furnished diverse experimental proposals for the emergence of first informational polymers that could form the basis of the early terrestrial life. Side by side with the experiments, the fast development of modern computational chemistry methods during the last 20 years facilitated the use of in silico modelling tools to complement the experiments. Modern computations can provide unique atomic-level insights into the structural and electronic aspects as well as the energetics of key prebiotic chemical reactions. Many of these insights are not directly obtainable from the experimental techniques and the computations are thus becoming indispensable for proper interpretation of many experiments and for qualified predictions. This review illustrates the synergy between experiment and theory in the origin of life research focusing on the prebiotic synthesis of various nucleic acid building blocks and on the self-assembly of nucleotides leading to the first functional oligonucleotides.

Distributed Multipolar Expansion Approach to Calculation of Excitation Energy Transfer Couplings.

We propose a new approach for estimating the electrostatic part of the excitation energy transfer (EET) coupling between electronically excited chromophores based on the transition density-derived cumulative atomic multipole moments (TrCAMM). In this approach, the transition potential of a chromophore is expressed in terms of truncated distributed multipolar expansion and analytical formulas for the TrCAMMs are derived. The accuracy and computational feasibility of the proposed approach is tested against the exact Coulombic couplings, and various multipole expansion truncation schemes are analyzed. The results of preliminary calculations show that the TrCAMM approach is capable of reproducing the exact Coulombic EET couplings accurately and efficiently and is superior to other widely used schemes: the transition charges from electrostatic potential (TrESP) and the transition density cube (TDC) method.