In silicoSupport for Eschenmoser’s Glyoxylate Scenario

Representing Catalytic Mechanisms with Rule Composition

An "imaginary transition structure" overlays the molecular graphs of the educt and product sides of an elementary chemical reaction in a single graph to highlight the changes in bond structure. We generalize this idea to reactions with complex mechanisms in a formally rigorous approach based on composing arrow-pushing steps represented as graph-transformation rules to construct an overall composite rule and a derived transition structure. This transition structure retains information about transient bond changes that are invisible at the overall level and can be constructed automatically from an existing database of detailed enzymatic mechanisms. We use the construction to (i) illuminate the distribution of catalytic action across enzymes and substrates and (ii) to search in a large database for reactions of known or unknown mechanisms that are compatible with the mechanism captured by the constructed composite rule.

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.

Prebiotic Reaction Networks in Water

A prevailing strategy in origins of life studies is to explore how chemistry constrained by hypothetical prebiotic conditions could have led to molecules and system level processes proposed to be important for life's beginnings. This strategy has yielded model prebiotic reaction networks that elucidate pathways by which relevant compounds can be generated, in some cases, autocatalytically. These prebiotic reaction networks provide a rich platform for further understanding and development of emergent "life-like" behaviours. In this review, recent advances in experimental and analytical procedures associated with classical prebiotic reaction networks, like formose and Miller-Urey, as well as more recent ones are highlighted. Instead of polymeric networks, i.e., those based on nucleic acids or peptides, the focus is on small molecules. The future of prebiotic chemistry lies in better understanding the genuine complexity that can result from reaction networks and the construction of a centralised database of reactions useful for predicting potential network evolution is emphasised.

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.

Chemical Transformation Motifs-Modelling Pathways as Integer Hyperflows

We present an elaborate framework for formally modelling pathways in chemical reaction networks on a mechanistic level. Networks are modelled mathematically as directed multi-hypergraphs, with vertices corresponding to molecules and hyperedges to reactions. Pathways are modelled as integer hyperflows and we expand the network model by detailed routing constraints. In contrast to the more traditional approaches like Flux Balance Analysis or Elementary Mode analysis we insist on integer-valued flows. While this choice makes it necessary to solve possibly hard integer linear programs, it has the advantage that more detailed mechanistic questions can be formulated. It is thus possible to query networks for general transformation motifs, and to automatically enumerate optimal and near-optimal pathways. Similarities and differences between our work and traditional approaches in metabolic network analysis are discussed in detail. To demonstrate the applicability of the mathematical framework to real-life problems we first explore the design space of possible non-oxidative glycolysis pathways and show that recent manually designed pathways can be further optimized. We then use a model of sugar chemistry to investigate pathways in the autocatalytic formose process. A graph transformation-based approach is used to automatically generate the reaction networks of interest.

Cyanamide as a prebiotic phosphate activating agent - catalysis by simple 2-oxoacid salts

Cyanamide is a prebiotically plausible compound that has previously been invoked as a phosphate activating agent. However, its reactions with phosphate monoesters are very slow and tend to be low yielding. We now report a fast and efficient phosphate activation reaction using cyanamide in the presence of glyoxylate or pyruvate. These simple 2-oxoacid salts are shown to function as catalysts and in an optimised system, adenosine-3'-phosphate was converted to adenosine-2',3'-cyclic phosphate in 95% yield.

Exploring astrobiology using in silico molecular structure generation

The origin of life is typically understood as a transition from inanimate or disorganized matter to self-organized, 'animate' matter. This transition probably took place largely in the context of organic compounds, and most approaches, to date, have focused on using the organic chemical composition of modern organisms as the main guide for understanding this process. However, it has gradually come to be appreciated that biochemistry, as we know it, occupies a minute volume of the possible organic 'chemical space'. As the majority of abiotic syntheses appear to make a large set of compounds not found in biochemistry, as well as an incomplete subset of those that are, it is possible that life began with a significantly different set of components. Chemical graph-based structure generation methods allow for exhaustive in silico enumeration of different compound types and different types of 'chemical spaces' beyond those used by biochemistry, which can be explored to help understand the types of compounds biology uses, as well as to understand the nature of abiotic synthesis, and potentially design novel types of living systems.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Computational exploration of the chemical structure space of possible reverse tricarboxylic acid cycle constituents

The reverse tricarboxylic acid (rTCA) cycle has been explored from various standpoints as an idealized primordial metabolic cycle. Its simplicity and apparent ubiquity in diverse organisms across the tree of life have been used to argue for its antiquity and its optimality. In 2000 it was proposed that chemoinformatics approaches support some of these views. Specifically, defined queries of the Beilstein database showed that the molecules of the rTCA are heavily represented in such compound databases. We explore here the chemical structure “space,” e.g. the set of organic compounds which possesses some minimal set of defining characteristics, of the rTCA cycle’s intermediates using an exhaustive structure generation method. The rTCA’s chemical space as defined by the original criteria and explored by our method is some six to seven times larger than originally considered. Acknowledging that each assumption in what is a defining criterion making the rTCA cycle special limits possible generative outcomes, there are many unrealized compounds which fulfill these criteria. That these compounds are unrealized could be due to evolutionary frozen accidents or optimization, though this optimization may also be for systems-level reasons, e.g., the way the pathway and its elements interface with other aspects of metabolism.

Chemical Transformations in Proto-Cytoplasmic Media. Phosphorus Coupling in the Silica Hydrogel Phase

It has been proposed that prebiotic chemical studies on the emergence of primitive life would be most relevant when performed in a hydrogel, rather than an aqueous, environment. In this paper we describe the ambient temperature coupling of phosphorus oxyacids [Pi] mediated by Fe(II) under aerobic conditions within a silica hydrogel (SHG) environment. We have chosen to examine SHGs as they have considerable geological precedence as key phases in silicification en route to rock formation. Following a description of the preparation and characterization studies on our SHG formulations, coupling experiments between Pi species are described across multiple permutations of (i) Pi compound; (ii) gel formulation; (iii) metal salt additive; and (iv) pH-modifying agent. The results suggest that successful Pi coupling, indicated by observation of pyrophosphate [PPi(V)] via ³¹P-NMR spectroscopy, takes place when the following components are present: (i) a mixture of mixture of Pi(III) and Pi(V) or pure PPi(III– V); (ii) Fe(II); (iii) acetic or formic acid (not hydrochloric acid); (iv) aerobic conditions or the presence of H₂O₂ as an oxidant; and (v) the presence of a gel system. On the basis of these, and aqueous control reactions, we suggest mechanistic possibilities.