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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

On the particular importance of vibrational contributions to the static electrical properties of model linear molecules under spatial confinement.

The influence of the spatial confinement on the electronic and vibrational contributions to longitudinal electric-dipole properties of model linear molecules including HCN, HCCH and CO2 is discussed. The effect of confinement is represented by two-dimensional harmonic oscillator potential of cylindrical symmetry, which mimics the key features of various types of trapping environments like, for instance, nanotubes or quantum well wires. Our results indicate that in general both (electronic and vibrational) contributions to (hyper)polarizabilities diminish upon spatial confinement. However, since the electronic term is particularly affected, the relative importance of vibrational contributions is larger for confined species. This effect increases also with the degree of anharmonicity of vibrational motion.

Molecular mechanism of diaminomaleonitrile to diaminofumaronitrile photoisomerization: an intermediate step in the prebiotic formation of purine nucleobases.

The photoinduced isomerization of diaminomaleonitrile (DAMN) to diaminofumaronitrile (DAFN) was suggested to play a key role in the prebiotically plausible formation of purine nucleobases and nucleotides. In this work we analyze two competitive photoisomerization mechanisms on the basis of state-of-the-art quantum-chemical calculations. Even though it was suggested that this process might occur on the triplet potential-energy surface, our results indicate that the singlet reaction channel should not be disregarded either. In fact, the peaked topography of the S1 /S0 conical intersection suggests that the deexcitation should most likely occur on a sub-picosecond timescale and the singlet photoisomerization mechanism might effectively compete even with a very efficient intersystem crossing. Such a scenario is further supported by the relatively small spin-orbit coupling of the S1 and T2 states in the Franck-Condon region, which does not indicate a very effective triplet bypass for this photoreaction. Therefore, we conclude that the triplet reaction channel in DAMN might not be as prominent as was previously thought.

Solvent effects on the photochemistry of 4-aminoimidazole-5-carbonitrile, a prebiotically plausible precursor of purines.

4-Aminoimidazole-5-carbonitrile (AICN) was suggested as a prebiotically plausible precursor of purine nucleobases and nucleotides. Although it can be formed in a sequence of photoreactions, AICN is immune to further irradiation with UV-light. We present state-of-the-art multi-reference quantum-chemical calculations of potential energy surface cuts and conical intersection optimizations to explain the molecular mechanisms underlying the photostability of this compound. We have identified the N-H bond stretching and ring-puckering mechanisms that should be responsible for the photochemistry of AICN in the gas phase. We have further considered the photochemistry of AICN-water clusters, while including up to six explicit water molecules. The calculations reveal charge transfer to solvent followed by formation of an H3O(+) cation, both of which occur on the (1)πσ* hypersurface. Interestingly, a second proton transfer to an adjacent water molecule leads to a (1)πσ*/S0 conical intersection. We suggest that this electron-driven proton relay might be characteristic of low-lying (1)πσ* states in chromophore-water clusters. Owing to its nature, this mechanism might also be responsible for the photostability of analogous organic molecules in bulk water.

On the performance of long-range-corrected density functional theory and reduced-size polarized LPol-n basis sets in computations of electric dipole (hyper)polarizabilities of π-conjugated molecules.

Static longitudinal electric dipole (hyper)polarizabilities are calculated for six medium-sized π-conjugated organic molecules using recently developed LPol-n basis set family to assess their performance. Dunning's correlation-consistent basis sets of triple-ζ quality combined with MP2 method and supported by CCSD(T)/aug-cc-pVDZ results are used to obtain the reference values of analyzed properties. The same reference is used to analyze (hyper)polarizabilities predicted by selected exchange-correlation functionals, particularly those asymptotically corrected.

Resonance-assisted hydrogen bonds revisited. Resonance stabilization vs. charge delocalization.

The origins of stabilization in the short strong hydrogen bonds commonly referred to as "resonance-assisted" (RAHB) have been revisited using the modern valence-bond theory, the hybrid variational-perturbational interaction energy decomposition scheme and atoms in molecules (AIM) analysis. Dimers of carboxylic acids and amides have been chosen as the model structures for intermolecular RAHBs, while for the intramolecular case malondialdehyde and its substituted derivatives have been selected. The estimated (negligible) resonance stabilization energies and relative magnitudes of interaction energy components indicate that the origin of stabilization in the studied complexes is charge-delocalization. Although in the case of intramolecular RAHBs the resonance effects are much more pronounced, still they are a relatively minor contribution to the total stabilization energy. In fact, the estimated resonance stabilization energies diminish with an increasing strength of the hydrogen bond (as indicated by AIM and structural descriptors).

Theoretical studies of the mechanism of 2-aminooxazole formation under prebiotically plausible conditions.

2-Aminooxazole is generally considered to play a central role in the origin of informational polymers. In the current contribution we use density functional calculations to investigate the detailed mechanism of 2-aminooxazole formation from the prebiotic soup according to the scenario suggested by M. W. Powner, B. Gerland and J. D. Sutherland, Nature, 2009, 459, 239-242. Parallel to the phosphate-catalyzed reaction pathway we also describe its water-assisted variant. Our calculations show that phosphate-catalysis is indispensable not only in the cyclization and the subsequent water-elimination steps, as previously suggested, but also in the very first reaction step leading to the formation of the carbinolamine intermediate. In addition, we suggest concurrent reaction channels for the cyclization and water-elimination reaction steps, both involving catalytic phosphate ions.