Insights Into the Origin of Life: Did It Begin from HCN and H2O?

Discovery of New Synthetic Routes of Amino Acids in Prebiotic Chemistry

The origin of life on Earth remains one of the most perplexing challenges in biochemistry. While numerous bottom-up experiments under prebiotic conditions have provided valuable insights into the spontaneous chemical genesis of life, there remains a significant gap in the theoretical understanding of the complex reaction processes involved. In this study, we propose a novel approach using a roto-translationally invariant potential (RTIP) formulated with pristine Cartesian coordinates to facilitate the simulation of chemical reactions. By employing RTIP pathway sampling to explore the reactivity of primitive molecules, we identified several low-energy reaction mechanisms, such as two-hydrogen-transfer hydrogenation and HCOOH-catalyzed hydration and amination. This led to the construction of a comprehensive reaction network, illustrating the synthesis pathways for glycine, serine, and alanine. Further thermodynamic analysis highlights the pivotal role of formaldimine as a key precursor in amino acid synthesis, owing to its more favorable reactivity in coupling reactions compared to the traditionally recognized hydrogen cyanide. Our study demonstrates that the RTIP methodology, coupled with a divide-and-conquer strategy, provides new insights into the simulation of complex reaction processes, offering promising applications for advancing organic design and synthesis.

Space Molecules in Superacidic Media─The Synthesis and Structure of Monoprotonated Cyanoacetylene and Diprotonated Dicyanoacetylene

The molecules cyanoacetylene and dicyanoacetylene detected in the atmosphere of Saturn's largest moon Titan were investigated in superacidic media XF/MF5 (X = H, D; M = As, Sb), XF/GeF4, and XF/BF3. Cyanoacetylene is obtained as a monoprotonated salt, while only the diprotonated salt of dicyanoacetylene has been isolated. The salts were characterized by vibrational spectroscopy, NMR spectroscopy, and single-crystal X-ray diffraction. The protonations lead to a significant shortening of the C≡N bond. The salts of [C3HXN][SbF6] (X = H, D) show very rare coupled oscillations in vibrational spectroscopy. The experimental results are discussed together with quantum chemical calculations at the M062X/aug-cc-pVTZ level of theory. The -M and -I effects of the C≡N moiety decrease by protonation of the nitrile moiety.

Unravelling the prebiotic origins of the simplest α-ketoacids in cometary ices: a computational investigation

We have employed the ab initio nanoreactor (AINR) and DFT calculations to explore how the soft impact of comets entering early earth's dense atmosphere could induce chemical reactions in trapped interstellar ice components, leading to the origin of glyoxylic and pyruvic acids the simplest α-ketoacids essential for prebiotic metabolic cycles.

Mimicking lightning-induced electrochemistry on the early Earth

To test the hypothesis that an abiotic Earth and its inert atmosphere could form chemically reactive carbon- and nitrogen-containing compounds, we designed a plasma electrochemical setup to mimic lightning-induced electrochemistry under steady-state conditions of the early Earth. Air-gap electrochemical reactions at air-water-ground interfaces lead to remarkable yields, with up to 40 moles of carbon dioxide being reduced into carbon monoxide and formic acid, and 3 moles of gaseous nitrogen being fixed into nitrate, nitrite, and ammonium ions, per mole of transmitted electrons. Interfaces enable reactants (e.g., minerals) that may have been on land, in lakes, and in oceans to participate in radical and redox reactions, leading to higher yields compared to gas-phase-only reactions. Cloud-to-ground lightning strikes could have generated high concentrations of reactive molecules locally, establishing diverse feedstocks for early life to emerge and survive globally.

Automated Exploration of Reaction Networks and Mechanisms Based on Metadynamics Nanoreactor Simulations

We developed an automated approach to construct a complex reaction network and explore the reaction mechanisms for numerous reactant molecules by integrating several theoretical approaches. Nanoreactor-type molecular dynamics was used to generate possible chemical reactions, in which the metadynamics was used to overcome the reaction barriers, and the semiempirical GFN2-xTB method was used to reduce the computational cost. Reaction events were identified from trajectories using the hidden Markov model based on the evolution of the molecular connectivity. This provided the starting points for further transition-state searches at the electronic structure levels of density functional theory to obtain the reaction mechanism. Finally, the entire reaction network containing multiple pathways was built. The feasibility and efficiency of the automated construction of the reaction network were investigated using the HCHO and NH3 biomolecular reaction and the reaction network for a multispecies system comprising dozens of HCN and H2O molecules. The results indicated that the proposed approach provides a valuable and effective tool for the automated exploration of the reaction networks.

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

Thermally Accessible Prebiotic Pathways for Forming Ribonucleic Acid and Protein Precursors from Aqueous Hydrogen Cyanide

The search for prebiotic chemical pathways to biologically relevant molecules is a long-standing puzzle that has generated a menagerie of competing hypotheses with limited experimental prospects for falsification. However, the advent of computational network exploration methodologies has created the opportunity to compare the kinetic plausibility of various channels and even propose new pathways. Here, the space of organic molecules that can be formed within four polar or pericyclic reactions from water and hydrogen cyanide (HCN), two established prebiotic candidates for generating biological precursors, was comprehensively explored with a state-of-the-art exploration algorithm. A surprisingly diverse reactivity landscape was revealed within just a few steps of these simple molecules. Reaction pathways to several biologically relevant molecules were discovered involving lower activation energies and fewer reaction steps compared with recently proposed alternatives. Accounting for water-catalyzed reactions qualitatively affects the interpretation of the network kinetics. The case-study also highlights omissions of simpler and lower barrier reaction pathways to certain products by other algorithms that qualitatively affect the interpretation of HCN reactivity.

The transition state region in nonsynchronous concerted reactions

The critical and vanishing points of the reaction force F(ξ) = -dV(ξ)/dξ yield five important coordinates (ξ_R, ξ_R^* , ξ_TS, ξ_P^* , ξ_P) along the intrinsic reaction coordinate (IRC) for a given concerted reaction or reaction step. These points partition the IRC into three well-defined regions, reactants (ξ_R→ξ_R^* ), transition state (ξ_R^* →ξ_P^* ), and products (ξ_P^* →ξ_P), with traditional roles of mostly structural changes associated with the reactants and products regions and mostly electronic activity associated with the transition state (TS) region. Following the evolution of chemical bonding along the IRC using formal descriptors of synchronicity, reaction electron flux, Wiberg bond orders, and their derivatives (or, more precisely, the intensity of the electron activity) unambiguously indicates that for nonsynchronous reactions, electron activity transcends the TS region and takes place well into the reactants and products regions. Under these circumstances, an extension of the TS region toward the reactants and products regions may occur.

Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.

Characterization of non-covalent contacts in mono- and di-halo substituted acetaldehydes: probing the substitution effects of electron donating and withdrawing groups

In the current work, a systematic evaluation of the different types of non-covalent interactions (NCIs) in acetaldehyde dimers, including dimers of mono-halo (XCH₂CHO)₂, di-halo (X₂CHCHO)₂ and tri-halo substituted (X₃CCHO)₂ acetaldehydes via the associated stabilization energy of these dimers has been performed. Furthermore, a topological analysis of the electron density based on the quantum theory of atoms in molecules (QTAIM) and non-covalent interaction reduced density gradient (NCI-RDG) isosurfaces has also been performed to evaluate the nature of these NCIs. The geometrical and electronic characteristics have been evaluated via the presence of different electron-donating groups (EDGs) and electron-withdrawing groups (EWGs) or substituents in dimers of these molecules, namely, XCH(Y)CHO and X₂C(Y)CHO (wherein X = -F, -Cl, and -Br and Y = -SO₃H, -CN, -NO₂, -NH₂, -CH₃, -OCH₃, and -SMe₃). The C-H⋯O, C-H⋯X, X⋯X, X⋯O and C⋯O tetrel bonded contacts have been recognized to play an important role in the stabilization of the formed dimers. This study also establishes the fact that the overall stability of the dimeric assemblies is governed by the contributions from the mutual and complex interplay of a variety of interactions in the investigated dimers. Hence considerations based on strong H-bond donor-acceptor characteristics hold relevance for simple systems only, but slight alteration in the electronic environment can affect the overall stabilization energies of the system being investigated and the nature of the interactions that contribute towards the same.

Multivariate Analysis Applied to Microwave-Driven Cyanide Polymerization: A Statistical View of a Complex System

For the first time, chemometrics was applied to the recently reported microwave-driven cyanide polymerization. Fast, easy, robust, low-cost, and green-solvent processes are characteristic of these types of reactions. These economic and environmental benefits, originally inspired by the constraints imposed by plausible prebiotic synthetic conditions, have taken advantage of the development of a new generation of HCN-derived multifunctional materials. HCN-derived polymers present tunable properties by temperature and reaction time. However, the apparently random behavior observed in the evolution of cyanide polymerizations, assisted by microwave radiation over time at different temperatures, leads us to study this highly complex system using multivariate analytical tools to have a proper view of the system. Two components are sufficient to explain between 84 and 98% of the total variance in the data in all principal component analyses. In addition, two components explain more than 91% of the total variance in the data in the case of principal component analysis for categorical data. These consistent statistical results indicate that microwave-driven polymerization is a more robust process than conventional thermal syntheses but also that plausible prebiotic chemistry in alkaline subaerial environments could be more complex than in the aerial part of these systems, presenting a clear example of the "messy chemistry" approach of interest in the research about the origins of life. In addition, the methodology discussed herein could be useful for the data analysis of extraterrestrial samples and for the design of soft materials, in a feedback view between prebiotic chemistry and materials science.

Algorithmic Explorations of Unimolecular and Bimolecular Reaction Spaces

Algorithmic reaction exploration based on transition state searches has already made inroads into many niche applications, but its potential as a general-purpose tool is still largely unrealized. Computational cost and the absence of benchmark problems involving larger molecules remain obstacles to further progress. Here an ultra-low cost exploration algorithm is implemented and used to explore the reactivity of unimolecular and bimolecular reactants, comprising a total of 581 reactions involving 51 distinct reactants. The algorithm discovers all established reaction pathways, where such comparisons are possible, while also revealing a much richer reactivity landscape, including lower barrier reaction pathways and a strong dependence of reaction conformation in the apparent barriers of the reported reactions. The diversity of these benchmarks illustrate that reaction exploration algorithms are approaching general-purpose capability.

Fully Automated Generation of Prebiotically Relevant Reaction Networks from Optimized Nanoreactor Simulations

The nanoreactor approach first introduced by the group of Martı́nez [Wang et al. Nat. Chem. 2014, 6, 1044-1048] has recently attracted much attention because of its ability to accelerate the discovery of reaction pathways. Here, we provide a comprehensive study of various simulation parameters and present an alternative implementation for the reactivity-enhancing spherical constraint function, as well as for the detection of reaction events. In this context, a fully automated postsimulation evaluation procedure based on RDKit and NetworkX analysis is introduced. The chemical and physical robustness of the procedure is examined by investigating the reactivity of selected homogeneous systems. The optimized procedure is applied at the GFN2-xTB level of theory to a system composed of HCN molecules and argon atoms, acting as a buffer, yielding prebiotically plausible primary and secondary precursors for the synthesis of RNA. Furthermore, the formose reaction network is explored leading to numerous sugar precursors. The discovered compounds reflect experimental findings; however, new synthetic routes and a large collection of exotic, highly reactive molecules are observed, highlighting the predictive power of the nanoreactor approach for unraveling the reactive manifold.

Formation and evolution of C-C, C-O, C[double bond, length as m-dash]O and C-N bonds in chemical reactions of prebiotic interest

A series of prebiotic chemical reactions yielding the precursor building blocks of amino acids, proteins and carbohydrates, starting solely from HCN and water is studied here. We closely follow the formation and evolution of the pivotal C-C, C-O, C[double bond, length as m-dash]O, and C-N bonds, which dictate the chemistry of the molecules of life. In many cases, formation of these bonds is set in motion by proton transfers in which individual water molecules act as catalysts so that water atoms end up in the products. Our results indicate that the prebiotic formation of carbon dioxide, formaldehyde, formic acid, formaldimine, glycolaldehyde, glycine, glycolonitrile, and oxazole derivatives, among others, are best described as highly nonsynchronous concerted single step processes. Nonetheless, for all reactions involving double proton transfer, the formation and breaking of O-H bonds around a particular O atom occur in a synchronous fashion, apparently independently from other primitive processes. For the most part, the first process to initiate seems to be the double proton transfer in the reactions where they are present, then bond breaking/formation around the reactive carbon in the carbonyl group and finally rupture of the C-N bonds in the appropriate cases, which are the most reluctant to break. Remarkably, within the limitations of our non-dynamical computational model, the wide ranges of temperature and pressure in which these reactions occur, downplay the problematic determination of the exact constraints on the early Earth.

Prebiotic synthesis of noncanonical nucleobases under plausible alkaline hydrothermal conditions

Herein, the potential of alkaline hydrothermal environments for the synthesis of possible ancestral pre-RNA nucleobases using cyanide as a primary source of carbon and nitrogen is described. Water cyanide polymerizations were assisted by microwave radiation to obtain high temperature and a relatively high pressure (MWR, 180 °C, 15 bar) and were also carried out using a conventional thermal system (CTS, 80 °C, 1 bar) to simulate subaerial and aerial hydrothermal conditions, respectively, on the early Earth. For these syntheses, the initial concentration of cyanide and the diffusion effects were studied. In addition, it is well known that hydrolysis conditions are directly related to the amount and diversity of organic molecules released from cyanide polymers. Thus, as a first step, we studied the effect of several hydrolysis procedures, generally used in prebiotic chemistry, on some of the potential pre-RNA nucleobases of interest, together with some of their isomers and/or deamination products, also presumably formed in these complex reactions. The results show that the alkaline hydrothermal scenarios with a relatively constant pH are good geological scenarios for the generation of noncanonical nucleobases using cyanide as a prebiotic precursor.

O-Acetyl Migration within the Sialic Acid Side Chain: A Mechanistic Study Using the Ab Initio Nanoreactor

Many disease-causing viruses target sialic acids on the surface of host cells. Some viruses bind preferentially to sialic acids with O-acetyl modification at the hydroxyl group of C7, C8, or C9 on the glycerol-like side chain. Studies of proteins binding to sialosides containing O-acetylated sialic acids are crucial in understanding the related diseases but experimentally difficult due to the lability of the ester group. We recently showed that O-acetyl migration among hydroxyl groups of C7, C8, and C9 in sialic acids occurs in all directions in a pH-dependent manner. In the current study, we elucidate a full mechanistic pathway for the migration of O-acetyl among C7, C8, and C9. We used an ab initio nanoreactor to explore potential reaction pathways and density functional theory, pK_a calculations, and umbrella sampling to investigate elementary steps of interest. We found that when a base is present, migration is easy in any direction and involves three key steps: deprotonation of the hydroxyl group, cyclization between the two carbons, and the migration of the O-acetyl group. This dynamic equilibrium may play a defensive role against pathogens that evolve to gain entry to the cell by binding selectively to one acetylation state.

Dissipative Photochemical Abiogenesis of the Purines

We have proposed that the abiogenesis of life around the beginning of the Archean may have been an example of "spontaneous" microscopic dissipative structuring of UV-C pigments under the prevailing surface ultraviolet solar spectrum. The thermodynamic function of these Archean pigments (the "fundamental molecules of life"), as for the visible pigments of today, was to dissipate the incident solar light into heat. We have previously described the non-equilibrium thermodynamics and the photochemical mechanisms which may have been involved in the dissipative structuring of the purines adenine and hypoxanthine from the common precursor molecules of hydrogen cyanide and water under this UV light. In this article, we extend our analysis to include the production of the other two important purines, guanine and xanthine. The photochemical reactions are presumed to occur within a fatty acid vesicle floating on a hot (∼80 ∘C) neutral pH ocean surface exposed to the prevailing UV-C light. Reaction-diffusion equations are resolved under different environmental conditions. Significant amounts of adenine (∼10-5 M) and guanine (∼10-6 M) are obtained within 60 Archean days, starting from realistic concentrations of the precursors hydrogen cyanide and cyanogen (∼10-5 M).

Horizontal gene transfer provides insights into the deep evolutionary history and biology of Trichinella

Evolution involves temporal changes in the characteristics of a species that are subsequently propagated or rejected through natural selection. In the case of parasites, host switching also plays a prominent role in the evolutionary process. These changes are rooted in genetic variation and gene flow where genes may be deleted, mutated (sequence), duplicated, rearranged and/or translocated and then transmitted through vertical gene transfer. However, the introduction of new genes is not driven only by Mendelian inheritance and mutation but also by the introduction of DNA from outside a lineage in the form of horizontal gene transfer between donor and recipient organisms. Once introduced and integrated into the biology of the recipient, vertical inheritance then perpetuates the newly acquired genetic factor, where further functionality may involve co-option of what has become a pre-existing physiological capacity. Upon sequencing the Trichinella spiralis (Clade I) genome, a cyanate hydratase (cyanase) gene was identified that is common among bacteria, fungi, and plants, but rarely observed among other eukaryotes. The sequence of the Trichinella cyanase gene clusters with those derived from the Kingdom Plantae in contrast to the genes found in some Clade III and IV nematodes that cluster with cyanases of bacterial origin. Phylogenetic analyses suggest that the Trichinella cyanase was acquired during the Devonian period and independently from those of other nematodes. These data may help inform us of the deep evolutionary history and ecological connectivity of early ancestors within the lineage of contemporary Trichinella. Further, in many extant organisms, cyanate detoxification has been largely superseded by energy requirements for metabolism. Thus, deciphering the function of Trichinella cyanase may provide new avenues for treatment and control.

Protocells: Milestones and Recent Advances

The origin of life is still one of humankind's great mysteries. At the transition between nonliving and living matter, protocells, initially featureless aggregates of abiotic matter, gain the structure and functions necessary to fulfill the criteria of life. Research addressing protocells as a central element in this transition is diverse and increasingly interdisciplinary. The authors review current protocell concepts and research directions, address milestones, challenges and existing hypotheses in the context of conditions on the early Earth, and provide a concise overview of current protocell research methods.

Identifying the wide diversity of extraterrestrial purine and pyrimidine nucleobases in carbonaceous meteorites

The lack of pyrimidine diversity in meteorites remains a mystery since prebiotic chemical models and laboratory experiments have predicted that these compounds can also be produced from chemical precursors found in meteorites. Here we report the detection of nucleobases in three carbonaceous meteorites using state-of-the-art analytical techniques optimized for small-scale quantification of nucleobases down to the range of parts per trillion (ppt). In addition to previously detected purine nucleobases in meteorites such as guanine and adenine, we identify various pyrimidine nucleobases such as cytosine, uracil, and thymine, and their structural isomers such as isocytosine, imidazole-4-carboxylic acid, and 6-methyluracil, respectively. Given the similarity in the molecular distribution of pyrimidines in meteorites and those in photon-processed interstellar ice analogues, some of these derivatives could have been generated by photochemical reactions prevailing in the interstellar medium and later incorporated into asteroids during solar system formation. This study demonstrates that a diversity of meteoritic nucleobases could serve as building blocks of DNA and RNA on the early Earth.

CARNOT: a Fragment-Based Direct Molecular Dynamics and Virtual-Reality Simulation Package for Reactive Systems

Traditionally, the study of reaction mechanisms of complex reaction systems such as combustion has been performed on an individual basis by optimizations of transition structure and minimum energy path or by reaction dynamics trajectory calculations for one elementary reaction at a time. It is effective, but time-consuming, whereas important and unexpected processes could have been missed. In this article, we present a direct molecular dynamics (DMD) approach and a virtual-reality simulation program, CARNOT, in which plausible chemical reactions are simulated simultaneously at finite temperature and pressure conditions. A key concept of the present ab initio molecular dynamics method is to partition a large, chemically reactive system into molecular fragments that can be adjusted on the fly of a DMD simulation. The theory represents an extension of the explicit polarization method to reactive events, called ReX-Pol. We propose a highest-and-lowest adapted-spin approximation to define the local spins of individual fragments, rather than treating the entire system by a delocalized wave function. Consequently, the present ab initio DMD can be applied to reactive systems consisting of an arbitrarily varying number of closed and open-shell fragments such as free radicals, zwitterions, and separate ions found in combustion and other reactions. A graph-data structure algorithm was incorporated in CARNOT for the analysis of reaction networks, suitable for reaction mechanism reduction. Employing the PW91 density functional theory and the 6-31+G(d) basis set, the capabilities of the CARNOT program were illustrated by a combustion reaction, consisting of 28 650 atoms, and by reaction network analysis that revealed a range of mechanistic and dynamical events. The method may be useful for applications to other types of complex reactions.

Automated Exploration of Prebiotic Chemical Reaction Space: Progress and Perspectives

Prebiotic chemistry often involves the study of complex systems of chemical reactions that form large networks with a large number of diverse species. Such complex systems may have given rise to emergent phenomena that ultimately led to the origin of life on Earth. The environmental conditions and processes involved in this emergence may not be fully recapitulable, making it difficult for experimentalists to study prebiotic systems in laboratory simulations. Computational chemistry offers efficient ways to study such chemical systems and identify the ones most likely to display complex properties associated with life. Here, we review tools and techniques for modelling prebiotic chemical reaction networks and outline possible ways to identify self-replicating features that are central to many origin-of-life models.

Direct evidence for a radiation-induced synthesis of acetonitrile and isoacetonitrile from a 1 : 1 CH4HCN complex at cryogenic temperatures: is it a missing link between inorganic and prebiotic astrochemistry?

Nitriles are important constituents of extraterrestrial media. Nitriles are supposed to play a crucial role in prebiotic chemistry occurring in the interstellar medium. In this work, we have investigated the low-temperature radiation-induced transformations of a 1 : 1 CH₄HCN complex as a plausible precursor of the simplest nitriles using the matrix isolation approach with FTIR spectroscopic detection. The parent complexes isolated in a noble gas (Ng) matrix were obtained by deposition of the CH₄/HCN/Ng gaseous mixture and characterized by comparison of experimental complexation-induced shifts of the HCN fundamentals with the results of the ab initio calculations. It was found that the X-ray irradiation of low-temperature matrices containing the isolated 1 : 1 CH₄HCN complex resulted in the formation of acetonitrile (CH₃CN) and isoacetonitrile (CH₃NC) and it appears to be the first experimental evidence for the formation of C₂ nitriles (acetonitrile and isoacetonitrile) from such a "building block". Additionally, a 1 : 1 CH₄HNC complex was tentatively assigned to the irradiated Ar and Kr matrices. It is demonstrated that the matrix has a strong effect on the CH₃CN/CH₃NC yield ratio, which dramatically increases in the row Ar < Kr < Xe. Also, the efficiency of the radiation-induced formation of the CH₄HNC complex was shown to decrease from Ar to Kr. It is believed that the proposed pathway for acetonitrile formation may be a significant step in the radiation-induced evolution leading to complex organic molecules and biomolecules under astrochemical conditions. Furthermore, the obtained results provide a prominent example of the impact of very weak intermolecular interactions on the radiation-induced transformations in cold media.

Self-catalytic mechanism of prebiotic reactions: from formamide to pterins and guanine

Reaction pathway of prebiotic reactions for formation of the pteridines: pterin, xanthopterine, isoxanthopterine and leucopterine, as well as the purine nucleobase guanine from pure formamide are presented. In these reactions, formamide or its tautomer, formimidic acid, play the role of proton-carrying catalyst. All required raw materials, such as hydrogen cyanide, ammonia, water, formic acid, urea, 2-aminomalononitrile, glyoxal, glyoxylic acid and oxalic acid needed in the self-catalyzed reactions are obtained by partial decomposition of formamide. We show that the prebiotic formation of nucleobases and pterins is closely linked and they probably coexisted at the beginning of chemical evolution.

Data-Driven Prediction of Formation Mechanisms of Lithium Ethylene Monocarbonate with an Automated Reaction Network

Interfacial reactions are notoriously difficult to characterize, and robust prediction of the chemical evolution and associated functionality of the resulting surface film is one of the grand challenges of materials chemistry. The solid-electrolyte interphase (SEI), critical to Li-ion batteries (LIBs), exemplifies such a surface film, and despite decades of work, considerable controversy remains regarding the major components of the SEI as well as their formation mechanisms. Here we use a reaction network to investigate whether lithium ethylene monocarbonate (LEMC) or lithium ethylene dicarbonate (LEDC) is the major organic component of the LIB SEI. Our data-driven, automated methodology is based on a systematic generation of relevant species using a general fragmentation/recombination procedure which provides the basis for a vast thermodynamic reaction landscape, calculated with density functional theory. The shortest pathfinding algorithms are employed to explore the reaction landscape and obtain previously proposed formation mechanisms of LEMC as well as several new reaction pathways and intermediates. For example, we identify two novel LEMC formation mechanisms: one which involves LiH generation and another that involves breaking the (CH₂)O-C(═O)OLi bond in LEDC. Most importantly, we find that all identified paths, which are also kinetically favorable under the explored conditions, require water as a reactant. This condition severely limits the amount of LEMC that can form, as compared with LEDC, a conclusion that has direct impact on the SEI formation in Li-ion energy storage systems. Finally, the data-driven framework presented here is generally applicable to any electrochemical system and expected to improve our understanding of surface passivation.

A Lizardite-HCN Interaction Leading the Increasing of Molecular Complexity in an Alkaline Hydrothermal Scenario: Implications for Origin of Life Studies

Hydrogen cyanide, HCN, is considered a fundamental molecule in chemical evolution. The named HCN polymers have been suggested as precursors of important bioorganics. Some novel researches have focused on the role of mineral surfaces in the hydrolysis and/or polymerization of cyanide species, but until now, their role has been unclear. Understanding the role of minerals in chemical evolution processes is crucial because minerals undoubtedly interacted with the organic molecules formed on the early Earth by different process. Therefore, we simulated the probable interactions between HCN and a serpentinite-hosted alkaline hydrothermal system. We studied the effect of serpentinite during the thermolysis of HCN at basic conditions (i.e., HCN 0.15 M, 50 h, 100 °C, pH > 10). The HCN-derived thermal polymer and supernatant formed after treatment were analyzed by several complementary analytical techniques. The results obtained suggest that: (I) the mineral surfaces can act as mediators in the mechanisms of organic molecule production such as the polymerization of HCN; (II) the thermal and physicochemical properties of the HCN polymer produced are affected by the presence of the mineral surface; and (III) serpentinite seems to inhibit the formation of bioorganic molecules compared with the control (without mineral).

Insights Into Chemical Reactions at the Beginning of the Universe: From HeH+ to H3. Front Chem 2021;9:679750. [PMID: 34222195 PMCID: PMC8249737 DOI: 10.3389/fchem.2021.679750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

At the dawn of the Universe, the ions of the light elements produced in the Big Bang nucleosynthesis recombined with each other. In our present study, we have tried to mimic the conditions in the early Universe to show how the recombination process would have led to the formation of the first ever formed diatomic species of the Universe: HeH+, as well as the subsequent processes that would have led to the formation of the simplest triatomic species: H3 +. We have also studied some special cases: higher positive charge with fewer number of hydrogen atoms in a dense atmosphere, and the formation of unusual and interesting linear, dicationic He chains beginning from light elements He and H in a positively charged atmosphere. For all the simulations, the ab initio nanoreactor (AINR) dynamics method has been employed.

An XPS study of HCN-derived films on pyrite surfaces: a prebiotic chemistry standpoint towards the development of protective coatings

Traditionally, the effect of mineral surfaces on increasing molecular complexity has been considered a major issue in studies about the origin of life. In contrast, herein, the effects of organic films derived from cyanide over an important prebiotic mineral, pyrite, are considered. An XPS spectroscopy study was carried out to understand the surface chemistry of the HCN-derived polymer/pyrite system. As a result, the simulation of a plausible prebiotic alkaline hydrothermal environment led to the identification of an NH4CN-based film with protective corrosion properties that immediately prevented the oxidation of the highly reactive pyrite surface. In addition, the effect of coating with antioxidant properties was preserved over a relatively long time, and the polymeric film was very stable under ambient conditions. These results increase the great potential of HCN polymers for development as a cheap and easily produced new class of multifunctional polymeric materials that also show promising and attractive insights into prebiotic chemistry.

The non-adiabatic nanoreactor: towards the automated discovery of photochemistry

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S₀/S₁ and S₁/S₂ conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

Characterization of the Simplest Thiolimine: The Higher Energy Tautomer of Thioformamide

As sulfur‐containing organic molecules thioamides and their isomers are conceivable intermediates in prebiotic chemistry, for example, in the formation of amino acids and thiazoles and resemble viable candidates for detection in interstellar media. Here, we report the characterization of parent thioformamide in the solid state via single‐crystal X‐ray diffraction and its photochemical interconversion to its hitherto unreported higher energy tautomer thiolimine in inert argon and dinitrogen matrices. Upon photogeneration, four conformers of thiolimine form, whose ratio depends on the employed wavelength. One of these conformers interconverts due to quantum mechanical tunneling with a half‐life of 30–45 min in both matrix materials at 3 and 20 K. A spontaneous reverse reaction from thiolimine to thioformamide is not observed. To support our experimental findings, we explored the potential energy surface of the system at the AE‐CCSD(T)/aug‐cc‐pCVTZ level of theory and computed tunneling half‐lives with the CVT/SCT approach applying DFT methods.

Formation of HCN+ in collisions of N+ and N2+ with a self-assembled propanethiol surface on gold

Collisions of N⁺ and N₂⁺ with C3 hydrocarbons, represented by a self assembled monolayer of propanethiol on a polcrystalline gold surface, were investigated by experiments over the incident energy range between 5 eV and 100 eV. For N⁺, formation of HCN⁺ is observed at incident energies of projectile ions as low as 20 eV. In the case of N₂⁺ projectile ions, the yield of HCN⁺ increased above zero only at incident energies of about 50 eV. This collision energy in the laboratory frame corresponds to an activation energy of about 3 eV to 3.5 eV. In the case of N⁺ projectile ions, the yield of HCN⁺ was large for most of the incident energy range, but decreased to zero at incident energies below 20 eV. This may indicate a very small energy threshold for the surface reaction between N⁺ and C3 hydrocarbons of a few tenths of an eV. Such a threshold for the formation of HCN⁺ may exist also for collisions of N⁺ with an adsorbed mixture of hydrocarbon molecules.

A Comprehensive Review of HCN-Derived Polymers

HCN-derived polymers are a heterogeneous group of complex substances synthesized from pure HCN; from its salts; from its oligomers, specifically its trimer and tetramer, amino-nalono-nitrile (AMN) and diamino-maleo-nitrile (DAMN), respectively; or from its hydrolysis products, such as formamide, under a wide range of experimental conditions. The characteristics and properties of HCN-derived polymers depend directly on the synthetic conditions used for their production and, by extension, their potential applications. These puzzling systems have been known mainly in the fields of prebiotic chemistry and in studies on the origins of life and astrobiology since the first prebiotic production of adenine by Oró in the early years of the 1960s. However, the first reference regarding their possible role in prebiotic chemistry was mentioned in the 19th century by Pflüger. Currently, HCN-derived polymers are considered keys in the formation of the first and primeval protometabolic and informational systems, and they may be among the most readily formed organic macromolecules in the solar system. In addition, HCN-derived polymers have attracted a growing interest in materials science due to their potential biomedical applications as coatings and adhesives; they have also been proposed as valuable models for multifunctional materials with emergent properties such as semi-conductivity, ferroelectricity, catalysis and photocatalysis, and heterogeneous organo-synthesis. However, the real structures and the formation pathways of these fascinating substances have not yet been fully elucidated; several models based on either computational approaches or spectroscopic and analytical techniques have endeavored to shed light on their complete nature. In this review, a comprehensive perspective of HCN-derived polymers is presented, taking into account all the aspects indicated above.

The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine

The non-equilibrium thermodynamics and the photochemical reaction mechanisms are described which may have been involved in the dissipative structuring, proliferation and complexation of the fundamental molecules of life from simpler and more common precursors under the UVC photon flux prevalent at the Earth’s surface at the origin of life. Dissipative structuring of the fundamental molecules is evidenced by their strong and broad wavelength absorption bands in the UVC and rapid radiationless deexcitation. Proliferation arises from the auto- and cross-catalytic nature of the intermediate products. Inherent non-linearity gives rise to numerous stationary states permitting the system to evolve, on amplification of a fluctuation, towards concentration profiles providing generally greater photon dissipation through a thermodynamic selection of dissipative efficacy. An example is given of photochemical dissipative abiogenesis of adenine from the precursor HCN in water solvent within a fatty acid vesicle floating on a hot ocean surface and driven far from equilibrium by the incident UVC light. The kinetic equations for the photochemical reactions with diffusion are resolved under different environmental conditions and the results analyzed within the framework of non-linear Classical Irreversible Thermodynamic theory.

Plausible Emergence of Biochemistry in Enceladus Based on Chemobrionics

Saturn's satellite Enceladus is proposed to have a soda-type subsurface ocean with temperature able to support life and an iron ore-based core. Here, it was demonstrated that ocean chemistry related to Enceladus can support the development of Fe-based hydrothermal vents, one of the places suggested to be the cradle of life. The Fe-based chemical gardens were characterized with Fourier-transform (FT)IR spectroscopy and XRD. The developed chemobrionic structures catalyzed the condensation polymerization of simple organic prebiotic molecules to kerogens. Further, they could passively catalyze the condensation of the prebiotic molecule formamide to larger polymers, suggesting that elementary biochemical precursors could have emerged in Enceladus.

A dual perspective on the microwave-assisted synthesis of HCN polymers towards the chemical evolution and design of functional materials

In this paper, the first study on NH4CN polymerization induced by microwave radiation is described, where a singular kinetic behaviour, especially when this reaction is conducted in the absence of air, is found. As a result, a complex conjugated N-heterocyclic polymer system is obtained, whose properties are very different, and even improved according to morphological features, characterized by their X-ray diffraction patterns and scanning electron microscopy analysis, with respect to those produced under conventional thermal treatment. In addition, a wide variety of relevant bioorganics have been identified, such as amino acids, nucleobases, co-factors, etc., from the synthetized NH4CN polymers. These particular families of polymers are of high interest in the fields of astrobiology and prebiotic chemistry and, more recently, in the development of smart multifunctional materials. From an astrobiological perspective, microwave-driven syntheses may simulate hydrothermal environments, which are considered ideal niches for increasing organic molecular complexity, and eventually as scenarios for an origin of life. From an industrial point of view and for potential applications, a microwave irradiation process leads to a notable decrease in the reaction times, and tune the properties of these new series macromolecular systems. The characteristics found for these materials encourage the development of further systematic research on this alternative HCN polymerization.

AFM Images of Viroid-Sized Rings That Self-Assemble from Mononucleotides through Wet-Dry Cycling: Implications for the Origin of Life

It is possible that early life relied on RNA polymers that served as ribozyme-like catalysts and for storing genetic information. The source of such polymers is uncertain, but previous investigations reported that wet–dry cycles simulating prebiotic hot springs provide sufficient energy to drive condensation reactions of mononucleotides to form oligomers. The aim of the study reported here was to visualize the products by atomic force microscopy. In addition to globular oligomers, ring-like structures ranging from 10–200 nm in diameter, with an average around 30–40 nm, were abundant, particularly when nucleotides capable of base pairing were present. The thickness of the rings was consistent with single stranded products, but some had thicknesses indicating base pair stacking. Others had more complex structures in the form of short polymer attachments and pairing of rings. These observations suggest the possibility that base-pairing may promote polymerization during wet–dry cycling followed by solvation of the rings. We conclude that RNA-like rings and structures could have been synthesized non-enzymatically on the prebiotic Earth, with sizes sufficient to fold into ribozymes and genetic molecules required for life to begin.

Synthetic connectivity, emergence, and self-regeneration in the network of prebiotic chemistry

The challenge of prebiotic chemistry is to trace the syntheses of life’s key building blocks from a handful of primordial substrates. Here we report a forward-synthesis algorithm that generates a full network of prebiotic chemical reactions accessible from these substrates under generally accepted conditions. This network contains both reported and previously unidentified routes to biotic targets, as well as plausible syntheses of abiotic molecules. It also exhibits three forms of nontrivial chemical emergence, as the molecules within the network can act as catalysts of downstream reaction types; form functional chemical systems, including self-regenerating cycles; and produce surfactants relevant to primitive forms of biological compartmentalization. To support these claims, computer-predicted, prebiotic syntheses of several biotic molecules as well as a multistep, self-regenerative cycle of iminodiacetic acid were validated by experiment.

Structure, Stability, and Spectroscopic Properties of Small Acetonitrile Cation Clusters

Ionization and fragmentation pathways induced by ionizing agents are key to understanding the formation of complex molecules in astrophysical environments. Acetonitrile (CH₃CN), the simplest organic nitrile, is an important molecule present in the interstellar medium. In this work, DFT and MP2 calculations were performed in order to obtain the low energy structures of the most relevant cations formed from electron-stimulated ion desorption of CH₃CN ices. Selected reaction pathways and spectroscopic properties were also calculated. Our results indicate that the most stable acetonitrile cation structure is CH₂CNH⁺ and that hydrogenation can occur successively without isomerization steps until its complete saturation. Moreover, the stability of distinct cluster families formed from the interaction of acetonitrile with small fragments, such as CH_n⁺, C₂H_n⁺, and CH_nCNH⁺, is discussed in terms of their respective binding energies. Some of these molecular clusters are stabilized by hydrogen bonds, leading to species whose infrared features are characterized by a strong redshift of the N-H stretching mode. Finally, the rotational spectra of CH₃CN and protonated acetonitrile, CH₃CNH⁺, were simulated using distinct computational protocols based on DFT, MP2, and CCSD(T) considering centrifugal distortion, vibrational-rotational coupling, and vibrational anharmonicity corrections. By adopting an empirical scaling procedure for calculating spectroscopic parameters, we were able to estimate the rotational frequencies of CH₃CNH⁺ with an expected average error below 1 MHz for J values up to 10.

Characterization of HCN-Derived Thermal Polymer: Implications for Chemical Evolution

Hydrogen cyanide (HCN)-derived polymers have been recognized as sources of relevant organic molecules in prebiotic chemistry and material sciences. However, there are considerable gaps in the knowledge regarding the polymeric nature, the physicochemical properties, and the chemical pathways along polymer synthesis. HCN might have played an important role in prebiotic hydrothermal environments; however, only few experiments use cyanide species considering hydrothermal conditions. In this work, we synthesized an HCN-derived thermal polymer simulating an alkaline hydrothermal environment (i.e., HCN (l) 0.15 M, 50 h, 100 °C, pH approximately 10) and characterized its chemical structure, thermal behavior, and the hydrolysis effect. Elemental analysis and infrared spectroscopy suggest an important oxidation degree. The thermal behavior indicates that the polymer is more stable compared to other HCN-derived polymers. The mass spectrometric thermal analysis showed the gradual release of several volatile compounds along different thermal steps. The results suggest a complicate macrostructure formed by amide and hydroxyl groups, which are joined to the main reticular chain with conjugated bonds (C=O, N=O, –O–C=N). The hydrolysis treatment showed the pH conditions for the releasing of organics. The study of the synthesis of HCN-derived thermal polymers under feasible primitive hydrothermal conditions is relevant for considering hydrothermal vents as niches of chemical evolution on early Earth.

Chemical evolution: from formamide to nucleobases and amino acids without the presence of catalyst

Abiotic synthesis of nucleobases and amino acids is of critical importance as it sheds light on potential prebiotic chemical reactions. During thermal decomposition of formamide in vacuum conditions, purine, cytosine, adenine, hypoxanthine, uracil, pterin, urea, urocanic acid, glycine, alanine and norvaline were detected. The compounds were obtained without catalyst by heating at 100-180 °C or microwave heating of formamide. Reaction network of self-catalyzed chemical reactions is suggested, showing how from only one parent molecule, nucleobases, urea and the amino acid glycine can be produced. The reaction pathways are theoretically determined using SCS-MP2 calculations.Communicated by Ramaswamy H. Sarma.

Layered Double Hydroxides in Bioinspired Nanotechnology

Layered Double Hydroxides (LDHs) are a relevant class of inorganic lamellar nanomaterials that have attracted significant interest in life science-related applications, due to their highly controllable synthesis and high biocompatibility. Under a general point of view, this class of materials might have played an important role for the origin of life on planet Earth, given their ability to adsorb and concentrate life-relevant molecules in sea environments. It has been speculated that the organic–mineral interactions could have permitted to organize the adsorbed molecules, leading to an increase in their local concentration and finally to the emergence of life. Inspired by nature, material scientists, engineers and chemists have started to leverage the ability of LDHs to absorb and concentrate molecules and biomolecules within life-like compartments, allowing to realize highly-efficient bioinspired platforms, usable for bioanalysis, therapeutics, sensors and bioremediation. This review aims at summarizing the latest evolution of LDHs in this research field under an unprecedented perspective, finally providing possible challenges and directions for future research.

Theoretical Investigations on the Possibility of Prebiotic HCN Formation via O-Addition Reactions

Until now, reactions between methane photolysis products (CH₃^•, CH₂) and active N atom or reactive NO radical are proposed as routes of HCN formation in the prebiotic Earth. Scientists think that the reducing atmosphere of primitive Earth was made of H₂, He, N₂, NO, CH₄, H₂O, CO₂, etc., and there was no molecular oxygen. However, it has been evident from experiments that the vacuum ultraviolet (VUV) photolysis of CO₂ can produce atomic oxygen. Therefore, it can be presumed that atomic oxygen was likely present in early Earth's atmosphere. Was there any impact of atomic oxygen in production of early atmospheric HCN for the emergence of life? To hunt for the answer, we have employed computational methods to study the mechanism and kinetics of CH₃NO + O(¹D) and CH₂NO^• + O(³P) addition reactions. Current study suggests that the addition of O(¹D) into nitrosomethane (CH₃NO) and the addition of O(³P) into nitrosomethylene radical (CH₂NO^•) can efficiently produce HCN through an effectively barrierless pathway. At STP, Bartis-Widom phenomenological loss rate coefficients of O(¹D) and O(³P) are obtained as 2.47 × 10^-12 and 4.67 × 10^-11 cm³ molecule^-1 s^-1, respectively. We propose that addition reactions of atomic oxygen with CH₃NO and CH₂NO^• might act as a potential source for early atmospheric HCN.

CH3NO as a potential intermediate for early atmospheric HCN: a quantum chemical insight

Hydrogen cyanide (HCN) has played a central role in the production of several biological molecules under prebiotic conditions on primitive Earth. Previously, K. J. Zahnle (J. Geophys. Res.: Atmos., 1986, 91, 2819) and Tian et al. (Earth Planet. Sci. Lett., 2011, 308, 417) emphasized that HCN production in the early Earth's CH₄-rich atmosphere could have been possible through the reaction between active nitrogen atoms (N) and methane photolysis products. Here, we have proposed alternative pathways for the formation of early atmospheric HCN via the decomposition of CH₃NO as an intermediate. In the early Earth's O₂-free atmosphere, CH₃˙ could preferentially attach to NO, which was generated via early atmospheric volcanism or lightning and photochemical processes. We have quantum chemically explored both unimolecular and bimolecular decomposition pathways of CH₃NO via the assistance of another CH₃NO molecule and via H₂O, NH₃, HCl, HCOOH, HNO₃ and H₂SO₄ catalysis. Both energetic and kinetic analyses reveal that H₂SO₄ is more efficient in this regard than other atmospheric species. Overall, it has been suggested that the proposed bimolecular decomposition pathways might have been alternative pathways for the formation of HCN under certain conditions on prebiotic Earth, while the unimolecular decomposition of CH₃NO could lead to the formation of HCN in the high temperature volcanic environment on early Earth.