Modern views of ancient metabolic networks

Primitive purine biosynthesis connects ancient geochemistry to modern metabolism

An unresolved question in the origin and evolution of life is whether a continuous path from geochemical precursors to the majority of molecules in the biosphere can be reconstructed from modern-day biochemistry. Here we identified a feasible path by simulating the evolution of biosphere-scale metabolism, using only known biochemical reactions and models of primitive coenzymes. We find that purine synthesis constitutes a bottleneck for metabolic expansion, which can be alleviated by non-autocatalytic phosphoryl coupling agents. Early phases of the expansion are enriched with enzymes that are metal dependent and structurally symmetric, supporting models of early biochemical evolution. This expansion trajectory suggests distinct hypotheses regarding the tempo, mode and timing of metabolic pathway evolution, including a late appearance of methane metabolisms and oxygenic photosynthesis consistent with the geochemical record. The concordance between biological and geological analyses suggests that this trajectory provides a plausible evolutionary history for the vast majority of core biochemistry.

Carbonyl Migration in Uronates Affords a Potential Prebiotic Pathway for Pentose Production

Carbohydrate biosynthesis is fundamental to modern terrestrial biochemistry, but how this collection of metabolic pathways originated remains an open question. Prebiotic sugar synthesis has focused primarily on the formose reaction and Kiliani-Fischer homologation; however, how they can transition to extant biochemical pathways has not been studied. Herein, a nonenzymatic pathway for pentose production with similar chemical transformations as those of the pentose phosphate pathway is demonstrated. Starting from a C6 aldonate, namely, gluconate, nonselective chemical oxidation yields a mixture of 2-oxo-, 4-oxo-, 5-oxo-, and 6-oxo-uronate regioisomers. Regardless at which carbinol the oxidation takes place, carbonyl migration enables β-decarboxylation to yield pentoses. In comparison, the pentose phosphate pathway selectively oxidizes 6-phosphogluconate to afford the 3-oxo-uronate derivative, which undergoes facile subsequent β-decarboxylation and carbonyl migration to afford ribose 5-phosphate. The similarities between these two pathways and the potential implications for prebiotic chemistry and protometabolism are discussed.

Free energy and inference in living systems

Organisms are non-equilibrium, stationary systems self-organized via spontaneous symmetry breaking and undergoing metabolic cycles with broken detailed balance in the environment. The thermodynamic free-energy (FE) principle describes an organism's homeostasis as the regulation of biochemical work constrained by the physical FE cost. By contrast, recent research in neuroscience and theoretical biology explains a higher organism's homeostasis and allostasis as Bayesian inference facilitated by the informational FE. As an integrated approach to living systems, this study presents an FE minimization theory overarching the essential features of both the thermodynamic and neuroscientific FE principles. Our results reveal that the perception and action of animals result from active inference entailed by FE minimization in the brain, and the brain operates as a Schrödinger's machine conducting the neural mechanics of minimizing sensory uncertainty. A parsimonious model suggests that the Bayesian brain develops the optimal trajectories in neural manifolds and induces a dynamic bifurcation between neural attractors in the process of active inference.

The Special and General Mechanism of Cyanobacterial Harmful Algal Blooms

Cyanobacterial harmful algal blooms (CyanoHABs) are longstanding aquatic hazards worldwide, of which the mechanism is not yet fully understood, i.e., the process in which cyanobacteria establish dominance over coexisting algae in the same eutrophic waters. The dominance of CyanoHABs represents a deviation from their low abundance under conventional evolution in the oligotrophic state, which has been the case since the origin of cyanobacteria on early Earth. To piece together a comprehensive mechanism of CyanoHABs, we revisit the origin and adaptive radiation of cyanobacteria in oligotrophic Earth, demonstrating ubiquitous adaptive radiation enabled by corresponding biological functions under various oligotrophic conditions. Next, we summarize the biological functions (ecophysiology) which drive CyanoHABs and ecological evidence to synthesize a working mechanism at the population level (the special mechanism) for CyanoHABs: CyanoHABs are the consequence of the synergistic interaction between superior cyanobacterial ecophysiology and elevated nutrients. Interestingly, these biological functions are not a result of positive selection by water eutrophication, but an adaptation to a longstanding oligotrophic state as all the genes in cyanobacteria are under strong negative selection. Last, to address the relative dominance of cyanobacteria over coexisting algae, we postulate a "general" mechanism of CyanoHABs at the community level from an energy and matter perspective: cyanobacteria are simpler life forms and thus have lower per capita nutrient demand for growth than coexisting eukaryotic algae. We prove this by comparing cyanobacteria and eukaryotic algae in cell size and structure, genome size, size of genome-scale metabolic networks, cell content, and finally the golden standard-field studies with nutrient supplementation in the same waters. To sum up, the comprehensive mechanism of CyanoHABs comprises a necessary condition, which is the general mechanism, and a sufficient condition, which is the special mechanism. One prominent prediction based on this tentative comprehensive mechanism is that eukaryotic algal blooms will coexist with or replace CyanoHABs if eutrophication continues and goes over the threshold nutrient levels for eukaryotic algae. This two-fold comprehensive mechanism awaits further theoretic and experimental testing and provides an important guide to control blooms of all algal species.

Hydrogen Drives Part of the Reverse Krebs Cycle under Metal or Meteorite Catalysis

Hydrogen (H₂ ) is a geological source of reducing electrons that is thought to have powered the metabolism of the last universal common ancestor to all extant life, and that is still metabolized by various modern organisms. It has been suggested that H₂ drove a geochemical analogue of some or all of the reverse Krebs cycle at the emergence of the metabolic network, catalyzed by metals, but this has yet to be demonstrated experimentally. Herein, we show that three consecutive steps of the reverse Krebs cycle, converting oxaloacetate into succinate, can be driven without enzymes and in one-pot by H₂ as the reducing agent under mild conditions compatible with biological chemistry. Low catalytic amounts of nickel (10-20 mol %) or platinum group metals (0.1-1 mol %) or even small amounts of ground meteorites were found to promote the reductive chemistry at temperatures between 5 and 60 °C and over a wide pH range, including pH 7. These results lend additional support to the hypothesis that geologically produced hydrogen and metal catalysts could have initiated early metabolic networks.

Building the uracil skeleton in primitive ponds at the origins of life: carbamoylation of aspartic acid

A large set of nucleobases and amino acids is found in meteorites, implying that several chemical reservoirs are present in the solar system. The "geochemical continuity" hypothesis explores how protometabolic paths developed from so-called "bricks" in an enzyme-free prebiotic world and how they affected the origins of life. In the living cell, the second step of synthesizing uridine and cytidine RNA monomers is a carbamoyl transfer from a carbamoyl donor to aspartic acid. Here we compare two enzyme-free scenarios: aqueous and mineral surface scenarios in a thermal range up to 250 °C. Both processes could have happened in ponds under open atmosphere on the primeval Earth. Carbamoylation of aspartic acid with cyanate in aqueous solutions at 25 °C gives high N-carbamoyl aspartic acid yields within 16 h. It is important to stress that, while various molecules could be efficient carbamoylating agents according to thermodynamics, kinetics plays a determining role in selecting prebiotically possible pathways.

The limits of metabolic heredity in protocells

The universal core of metabolism could have emerged from thermodynamically favoured prebiotic pathways at the origin of life. Starting with H 2 and CO 2 , the synthesis of amino acids and mixed fatty acids, which self-assemble into protocells, is favoured under warm anoxic conditions. Here, we address whether it is possible for protocells to evolve greater metabolic complexity, through positive feedbacks involving nucleotide catalysis. Using mathematical simulations to model metabolic heredity in protocells, based on branch points in protometabolic flux, we show that nucleotide catalysis can indeed promote protocell growth. This outcome only occurs when nucleotides directly catalyse CO 2 fixation. Strong nucleotide catalysis of other pathways (e.g. fatty acids and amino acids) generally unbalances metabolism and slows down protocell growth, and when there is competition between catalytic functions cell growth collapses. Autocatalysis of nucleotide synthesis can promote growth but only if nucleotides also catalyse CO 2 fixation; autocatalysis alone leads to the accumulation of nucleotides at the expense of CO 2 fixation and protocell growth rate. Our findings offer a new framework for the emergence of greater metabolic complexity, in which nucleotides catalyse broad-spectrum processes such as CO 2 fixation, hydrogenation and phosphorylation important to the emergence of genetic heredity at the origin of life.

Protein cost minimization promotes the emergence of coenzyme redundancy

Metabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes [e.g., NAD(H) and NADP(H)] associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling and structural and sequence analysis indicates that coenzyme redundancy may not be essential for metabolism but could rather constitute an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate and acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules [e.g., NAD(H) and NADP(H)] that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here, we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy and that coenzymes typically thought to exist in a single pool (e.g., coenzyme A [CoA]) may exist in more than one form (e.g., dephospho-CoA).

Computational Analysis of a Prebiotic Amino Acid Synthesis with Reference to Extant Codon-Amino Acid Relationships

Novel density functional theory calculations are presented regarding a mechanism for prebiotic amino acid synthesis from alpha-keto acids that was suggested to happen via catalysis by dinucleotide species. Our results were analysed with comparison to the original hypothesis (Copley et al., PNAS, 2005, 102, 4442–4447). It was shown that the keto acid–dinucleotide hypothesis for possible prebiotic amino acid synthesis was plausible based on an initial computational analysis, and details of the structures for the intermediates and transition states showed that there was wide scope for interactions between the keto acid and dinucleotide moieties that could affect the free energy profiles and lead to the required proto-metabolic selectivity.

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

The Semi-Enzymatic Origin of Metabolic Pathways: Inferring a Very Early Stage of the Evolution of Life

The early evolution of life is a period with many important events with a lot of big and open questions. One of them is the evolution of metabolic pathways, which means the origin and assembly of enzymes that act together. The retrograde hypothesis was the first attempt to explain the origin and evolutionary history of metabolic pathways; Norman Horowitz developed this first significant hypothesis. This idea was followed by relevant proposals developed by Sam Granick, who proposed the "forward direction hypothesis," and then the successful idea of "Patchwork" assembly proposed independently by Martynas Yčas and Roy Jensen. Since then, a few new hypotheses were proposed; one of the most influential was made by Antonio Lazcano and Stanley Miller in the Journal Molecular Evolution, the "semi-enzymatic origin" of metabolic pathways. This article was cited more than 160 times, including in most papers published about the early evolution of metabolism, placing it as influential work in the field. The ideas proposed in this work and their effects on studying the origin and early evolution of life are analyzed.

A thermodynamic atlas of carbon redox chemical space

Carbon redox chemistry plays a fundamental role in biology. However, the thermodynamic and physicochemical principles underlying the rise of metabolites involved in redox biochemistry remain poorly understood. Our work introduces the theory and techniques that allow us to quantify and understand the global energy landscape of carbon redox biochemistry. We analyze the space of all possible oxidation states of linear-chain molecules with two to five carbon atoms and generate a detailed atlas of the thermodynamic stability of metabolites in comparison to nonbiological molecules. Although the emergence of life required the underlying chemistry to bootstrap itself out of equilibrium, a quantitative understanding of the environment-dependent thermodynamic landscape of prebiotic molecules will be extremely valuable for future origins of life models.

Redox biochemistry plays a key role in the transduction of chemical energy in living systems. However, the compounds observed in metabolic redox reactions are a minuscule fraction of chemical space. It is not clear whether compounds that ended up being selected as metabolites display specific properties that distinguish them from nonbiological compounds. Here, we introduce a systematic approach for comparing the chemical space of all possible redox states of linear-chain carbon molecules to the corresponding metabolites that appear in biology. Using cheminformatics and quantum chemistry, we analyze the physicochemical and thermodynamic properties of the biological and nonbiological compounds. We find that, among all compounds, aldose sugars have the highest possible number of redox connections to other molecules. Metabolites are enriched in carboxylic acid functional groups and depleted of ketones and aldehydes and have higher solubility than nonbiological compounds. Upon constructing the energy landscape for the full chemical space as a function of pH and electron-donor potential, we find that metabolites tend to have lower Gibbs energies than nonbiological molecules. Finally, we generate Pourbaix phase diagrams that serve as a thermodynamic atlas to indicate which compounds are energy minima in redox chemical space across a set of pH values and electron-donor potentials. While escape from thermodynamic equilibrium toward kinetically driven states is a hallmark of life and its origin, we envision that a deeper quantitative understanding of the environment-dependent thermodynamic landscape of putative prebiotic molecules will provide a crucial reference for future origins-of-life models.

Future climates: Markov blankets and active inference in the biosphere

We formalize the Gaia hypothesis about the Earth climate system using advances in theoretical biology based on the minimization of variational free energy. This amounts to the claim that non-equilibrium steady-state dynamics-that underwrite our climate-depend on the Earth system possessing a Markov blanket. Our formalization rests on how the metabolic rates of the biosphere (understood as Markov blanket's internal states) change with respect to solar radiation at the Earth's surface (i.e. external states), through the changes in greenhouse and albedo effects (i.e. active states) and ocean-driven global temperature changes (i.e. sensory states). Describing the interaction between the metabolic rates and solar radiation as climatic states-in a Markov blanket-amounts to describing the dynamics of the internal states as actively inferring external states. This underwrites climatic non-equilibrium steady-state through free energy minimization and thus a form of planetary autopoiesis.

An ecological framework for the analysis of prebiotic chemical reaction networks

It is becoming widely accepted that very early in life's origin, even before the emergence of genetic encoding, reaction networks of diverse small chemicals might have manifested key properties of life, namely self-propagation and adaptive evolution. To explore this possibility, we formalize the dynamics of chemical reaction networks within the framework of chemical ecosystem ecology. To capture the idea that life-like chemical systems are maintained out of equilibrium by fluxes of energy-rich food chemicals, we model chemical ecosystems in well-mixed compartments that are subject to constant dilution by a solution with a fixed concentration of input chemicals. Modelling all chemical reactions as fully reversible, we show that seeding an autocatalytic cycle with tiny amounts of one or more of its member chemicals results in logistic growth of all member chemicals in the cycle. This finding justifies drawing an instructive analogy between an autocatalytic cycle and a biological species. We extend this finding to show that pairs of autocatalytic cycles can exhibit competitive, predator-prey, or mutualistic associations just like biological species. Furthermore, when there is stochasticity in the environment, particularly in the seeding of autocatalytic cycles, chemical ecosystems can show complex dynamics that can resemble evolution. The evolutionary character is especially clear when the network architecture results in ecological precedence, which makes a system's trajectory historically contingent on the order in which cycles are seeded. For all its simplicity, the framework developed here helps explain the onset of adaptive evolution in prebiotic chemical reaction networks, and can shed light on the origin of key biological attributes such as thermodynamic irreversibility and genetic encoding.

First arrived, first served: competition between codons for codon-amino acid stereochemical interactions determined early genetic code assignments

Stereochemical nucleotide-amino acid interactions, in the form of noncovalent nucleotide-amino acid interactions, potentially produced the genetic code's codon-amino acid assignments. Empirical estimates of single nucleotide-amino acid affinities on surfaces and in solution are used to test whether trinucleotide-amino acid affinities determined genetic code assignments pending the principle "first arrived, first served": presumed early amino acids have greater codon-amino acid affinities than ulterior ones. Here, these single nucleotide affinities are used to approximate all 64 × 20 trinucleotide-amino acid affinities. Analyses show that (1) on surfaces, genetic code codon-amino acid assignments tend to match high affinities for the amino acids that integrated earliest the genetic code (according to Wong's metabolic coevolution hypothesis between nucleotides and amino acids) and (2) in solution, the same principle holds for the anticodon-amino acid assignments. Affinity analyses match best genetic code assignments when assuming that trinucleotides competed for amino acids, rather than amino acids for trinucleotides. Codon-amino acid affinities stick better to genetic code assignments than anticodon-amino acid affinities. Presumably, two independent coding systems, on surfaces and in solution, converged, and formed the current translation system. Proto-translation on surfaces by direct codon-amino acid interactions without tRNA-like adaptors coadapted with a system emerging in solution by proto-tRNA anticodon-amino acid interactions. These systems assigned identical or similar cognates to codons on surfaces and to anticodons in solution. Results indicate that a prebiotic metabolism predated genetic code self-organization.

Small protein folds at the root of an ancient metabolic network

Life on Earth is driven by electron transfer reactions catalyzed by a suite of enzymes that comprise the superfamily of oxidoreductases (Enzyme Classification EC1). Most modern oxidoreductases are complex in their structure and chemistry and must have evolved from a small set of ancient folds. Ancient oxidoreductases from the Archean Eon between ca. 3.5 and 2.5 billion years ago have been long extinct, making it challenging to retrace evolution by sequence-based phylogeny or ancestral sequence reconstruction. However, three-dimensional topologies of proteins change more slowly than sequences. Using comparative structure and sequence profile-profile alignments, we quantify the similarity between proximal cofactor-binding folds and show that they are derived from a common ancestor. We discovered that two recurring folds were central to the origin of metabolism: ferredoxin and Rossmann-like folds. In turn, these two folds likely shared a common ancestor that, through duplication, recruitment, and diversification, evolved to facilitate electron transfer and catalysis at a very early stage in the origin of metabolism.

On Computing Structural and Behavioral Complexities of Threshold Boolean Networks : Application to Biological Networks

Various threshold Boolean networks (TBNs), a formalism used to model different types of biological networks (genes notably), can produce similar dynamics, i.e. share same behaviors. Among them, some are complex (according to Kolmogorov complexity), others not. By computing both structural and behavioral complexities, we show that most TBNs are structurally complex, even those having simple behaviors. For this purpose, we developed a new method to compute the structural complexity of a TBN based on estimates of the sizes of equivalence classes of the threshold Boolean functions composing the TBN.

The Future of Origin of Life Research: Bridging Decades-Old Divisions

Research on the origin of life is highly heterogeneous. After a peculiar historical development, it still includes strongly opposed views which potentially hinder progress. In the 1st Interdisciplinary Origin of Life Meeting, early-career researchers gathered to explore the commonalities between theories and approaches, critical divergence points, and expectations for the future. We find that even though classical approaches and theories-e.g. bottom-up and top-down, RNA world vs. metabolism-first-have been prevalent in origin of life research, they are ceasing to be mutually exclusive and they can and should feed integrating approaches. Here we focus on pressing questions and recent developments that bridge the classical disciplines and approaches, and highlight expectations for future endeavours in origin of life research.

Environmental boundary conditions for the origin of life converge to an organo-sulfur metabolism

It has been suggested that a deep memory of early life is hidden in the architecture of metabolic networks, whose reactions could have been catalyzed by small molecules or minerals before genetically encoded enzymes. A major challenge in unravelling these early steps is assessing the plausibility of a connected, thermodynamically consistent proto-metabolism under different geochemical conditions, which are still surrounded by high uncertainty. Here we combine network-based algorithms with physico-chemical constraints on chemical reaction networks to systematically show how different combinations of parameters (temperature, pH, redox potential and availability of molecular precursors) could have affected the evolution of a proto-metabolism. Our analysis of possible trajectories indicates that a subset of boundary conditions converges to an organo-sulfur-based proto-metabolic network fuelled by a thioester- and redox-driven variant of the reductive tricarboxylic acid cycle that is capable of producing lipids and keto acids. Surprisingly, environmental sources of fixed nitrogen and low-potential electron donors are not necessary for the earliest phases of biochemical evolution. We use one of these networks to build a steady-state dynamical metabolic model of a protocell, and find that different combinations of carbon sources and electron donors can support the continuous production of a minimal ancient 'biomass' composed of putative early biopolymers and fatty acids.

Twenty Years of "Lipid World": A Fertile Partnership with David Deamer

"The Lipid World" was published in 2001, stemming from a highly effective collaboration with David Deamer during a sabbatical year 20 years ago at the Weizmann Institute of Science in Israel. The present review paper highlights the benefits of this scientific interaction and assesses the impact of the lipid world paper on the present understanding of the possible roles of amphiphiles and their assemblies in the origin of life. The lipid world is defined as a putative stage in the progression towards life's origin, during which diverse amphiphiles or other spontaneously aggregating small molecules could have concurrently played multiple key roles, including compartment formation, the appearance of mutually catalytic networks, molecular information processing, and the rise of collective self-reproduction and compositional inheritance. This review brings back into a broader perspective some key points originally made in the lipid world paper, stressing the distinction between the widely accepted role of lipids in forming compartments and their expanded capacities as delineated above. In the light of recent advancements, we discussed the topical relevance of the lipid worldview as an alternative to broadly accepted scenarios, and the need for further experimental and computer-based validation of the feasibility and implications of the individual attributes of this point of view. Finally, we point to possible avenues for exploring transition paths from small molecule-based noncovalent structures to more complex biopolymer-containing proto-cellular systems.

Genomic Evidence for Simultaneous Optimization of Transcription and Translation through Codon Variants in the pmoCAB Operon of Type Ia Methanotrophs

Microbial methane oxidation plays a fundamental role in the biogeochemical cycle of Earth’s system. Recent reports have provided evidence for the acquisition of methane monooxygenases by horizontal gene transfer in methane-oxidizing bacteria from different environments, but how evolution has shaped the coding sequences to execute methanotrophy efficiently remains unexplored. In this work, we provide genomic evidence that among the different types of methanotrophs, type Ia methanotrophs possess a unique coding sequence of the pmoCAB operon that is under positive selection for optimal resource allocation and efficient synthesis of transcripts and proteins. This adaptive trait possibly enables type Ia methanotrophs to respond robustly to fluctuating methane availability and explains their global prevalence.

Understanding the interplay between genotype and phenotype is a fundamental goal of functional genomics. Methane oxidation is a microbial phenotype with global-scale significance as part of the carbon biogeochemical cycle and a sink for greenhouse gas. Microorganisms that oxidize methane (methanotrophs) are taxonomically diverse and widespread around the globe. In methanotrophic bacteria, enzymes in the methane oxidation metabolic module (KEGG module M00174, conversion of methane to formaldehyde) are encoded in four operons (pmoCAB, mmoXYZBCD, mxaFI, and xoxF). Recent reports have suggested that methanotrophs in Proteobacteria acquired methane monooxygenases through horizontal gene transfer. Here, we used a genomic meta-analysis to infer the transcriptional and translational advantages of coding sequences from the methane oxidation metabolic modules of different types of methanotrophs. By analyzing isolate and metagenome-assembled genomes from phylogenetically and geographically diverse sources, we detected an anomalous nucleotide composition bias in the coding sequences of particulate methane monooxygenase genes (pmoCAB) from type Ia methanotrophs. We found that this nucleotide bias increases the level of codon bias by decreasing the GC content in the third base of codons, a strategy that contrasts with that of other coding sequences in the module. Further codon usage analyses uncovered that codon variants of the type Ia pmoCAB coding sequences deviate from the genomic signature to match ribosomal protein-coding sequences. Subsequently, computation of transcription and translation metrics revealed that the pmoCAB coding sequences of type Ia methanotrophs optimize the usage of codon variants to maximize translation efficiency and accuracy, while minimizing the synthesis cost of transcripts and proteins.

IMPORTANCE Microbial methane oxidation plays a fundamental role in the biogeochemical cycle of Earth’s system. Recent reports have provided evidence for the acquisition of methane monooxygenases by horizontal gene transfer in methane-oxidizing bacteria from different environments, but how evolution has shaped the coding sequences to execute methanotrophy efficiently remains unexplored. In this work, we provide genomic evidence that among the different types of methanotrophs, type Ia methanotrophs possess a unique coding sequence of the pmoCAB operon that is under positive selection for optimal resource allocation and efficient synthesis of transcripts and proteins. This adaptive trait possibly enables type Ia methanotrophs to respond robustly to fluctuating methane availability and explains their global prevalence.

Synthesis and breakdown of universal metabolic precursors promoted by iron

Life builds its molecules from carbon dioxide (CO₂) and breaks them back down again through the intermediacy of just five metabolites, which are the universal hubs of biochemistry¹. However, it is unclear how core biological metabolism began and why it uses the intermediates, reactions and pathways that it does. Here we describe a purely chemical reaction network promoted by ferrous iron, in which aqueous pyruvate and glyoxylate-two products of abiotic CO₂ reduction^2-4-build up 9 of the 11 intermediates of the biological Krebs (or tricarboxylic acid) cycle, including all 5 universal metabolic precursors. The intermediates simultaneously break down to CO₂ in a life-like regime that resembles biological anabolism and catabolism⁵. Adding hydroxylamine^6-8 and metallic iron into the system produces four biological amino acids in a manner that parallels biosynthesis. The observed network overlaps substantially with the Krebs and glyoxylate cycles^9,10, and may represent a prebiotic precursor to these core metabolic pathways.

Autocatalytic Networks at the Basis of Life's Origin and Organization

Life is more than the sum of its constituent molecules. Living systems depend on a particular chemical organization, i.e., the ways in which their constituent molecules interact and cooperate with each other through catalyzed chemical reactions. Several abstract models of minimal life, based on this idea of chemical organization and also in the context of the origin of life, were developed independently in the 1960s and 1970s. These models include hypercycles, chemotons, autopoietic systems, (M,R)-systems, and autocatalytic sets. We briefly compare these various models, and then focus more specifically on the concept of autocatalytic sets and their mathematical formalization, RAF theory. We argue that autocatalytic sets are a necessary (although not sufficient) condition for life-like behavior. We then elaborate on the suggestion that simple inorganic molecules like metals and minerals may have been the earliest catalysts in the formation of prebiotic autocatalytic sets, and how RAF theory may also be applied to systems beyond chemistry, such as ecology, economics, and cognition.

Serpentinization: Connecting Geochemistry, Ancient Metabolism and Industrial Hydrogenation

Rock⁻water⁻carbon interactions germane to serpentinization in hydrothermal vents have occurred for over 4 billion years, ever since there was liquid water on Earth. Serpentinization converts iron(II) containing minerals and water to magnetite (Fe₃O₄) plus H₂. The hydrogen can generate native metals such as awaruite (Ni₃Fe), a common serpentinization product. Awaruite catalyzes the synthesis of methane from H₂ and CO₂ under hydrothermal conditions. Native iron and nickel catalyze the synthesis of formate, methanol, acetate, and pyruvate-intermediates of the acetyl-CoA pathway, the most ancient pathway of CO₂ fixation. Carbon monoxide dehydrogenase (CODH) is central to the pathway and employs Ni⁰ in its catalytic mechanism. CODH has been conserved during 4 billion years of evolution as a relic of the natural CO₂-reducing catalyst at the onset of biochemistry. The carbide-containing active site of nitrogenase-the only enzyme on Earth that reduces N₂-is probably also a relic, a biological reconstruction of the naturally occurring inorganic catalyst that generated primordial organic nitrogen. Serpentinization generates Fe₃O₄ and H₂, the catalyst and reductant for industrial CO₂ hydrogenation and for N₂ reduction via the Haber⁻Bosch process. In both industrial processes, an Fe₃O₄ catalyst is matured via H₂-dependent reduction to generate Fe₅C₂ and Fe₂N respectively. Whether serpentinization entails similar catalyst maturation is not known. We suggest that at the onset of life, essential reactions leading to reduced carbon and reduced nitrogen occurred with catalysts that were synthesized during the serpentinization process, connecting the chemistry of life and Earth to industrial chemistry in unexpected ways.

An appeal to magic? The discovery of a non-enzymatic metabolism and its role in the origins of life

Until recently, prebiotic precursors to metabolic pathways were not known. In parallel, chemistry achieved the synthesis of amino acids and nucleotides only in reaction sequences that do not resemble metabolic pathways, and by using condition step changes, incompatible with enzyme evolution. As a consequence, it was frequently assumed that the topological organisation of the metabolic pathway has formed in a Darwinian process. The situation changed with the discovery of a non-enzymatic glycolysis and pentose phosphate pathway. The suite of metabolism-like reactions is promoted by a metal cation, (Fe(II)), abundant in Archean sediment, and requires no condition step changes. Knowledge about metabolism-like reaction topologies has accumulated since, and supports non-enzymatic origins of gluconeogenesis, the S-adenosylmethionine pathway, the Krebs cycle, as well as CO2 fixation. It now feels that it is only a question of time until essential parts of metabolism can be replicated non-enzymatically. Here, I review the 'accidents' that led to the discovery of the non-enzymatic glycolysis, and on the example of a chemical network based on hydrogen cyanide, I provide reasoning why metabolism-like non-enzymatic reaction topologies may have been missed for a long time. Finally, I discuss that, on the basis of non-enzymatic metabolism-like networks, one can elaborate stepwise scenarios for the origin of metabolic pathways, a situation that increasingly renders the origins of metabolism a tangible problem.

The last universal common ancestor between ancient Earth chemistry and the onset of genetics

All known life forms trace back to a last universal common ancestor (LUCA) that witnessed the onset of Darwinian evolution. One can ask questions about LUCA in various ways, the most common way being to look for traits that are common to all cells, like ribosomes or the genetic code. With the availability of genomes, we can, however, also ask what genes are ancient by virtue of their phylogeny rather than by virtue of being universal. That approach, undertaken recently, leads to a different view of LUCA than we have had in the past, one that fits well with the harsh geochemical setting of early Earth and resembles the biology of prokaryotes that today inhabit the Earth's crust.

Systems protobiology: origin of life in lipid catalytic networks

Life is that which replicates and evolves, but there is no consensus on how life emerged. We advocate a systems protobiology view, whereby the first replicators were assemblies of spontaneously accreting, heterogeneous and mostly non-canonical amphiphiles. This view is substantiated by rigorous chemical kinetics simulations of the graded autocatalysis replication domain (GARD) model, based on the notion that the replication or reproduction of compositional information predated that of sequence information. GARD reveals the emergence of privileged non-equilibrium assemblies (composomes), which portray catalysis-based homeostatic (concentration-preserving) growth. Such a process, along with occasional assembly fission, embodies cell-like reproduction. GARD pre-RNA evolution is evidenced in the selection of different composomes within a sparse fitness landscape, in response to environmental chemical changes. These observations refute claims that GARD assemblies (or other mutually catalytic networks in the metabolism first scenario) cannot evolve. Composomes represent both a genotype and a selectable phenotype, anteceding present-day biology in which the two are mostly separated. Detailed GARD analyses show attractor-like transitions from random assemblies to self-organized composomes, with negative entropy change, thus establishing composomes as dissipative systems-hallmarks of life. We show a preliminary new version of our model, metabolic GARD (M-GARD), in which lipid covalent modifications are orchestrated by non-enzymatic lipid catalysts, themselves compositionally reproduced. M-GARD fills the gap of the lack of true metabolism in basic GARD, and is rewardingly supported by a published experimental instance of a lipid-based mutually catalytic network. Anticipating near-future far-reaching progress of molecular dynamics, M-GARD is slated to quantitatively depict elaborate protocells, with orchestrated reproduction of both lipid bilayer and lumenal content. Finally, a GARD analysis in a whole-planet context offers the potential for estimating the probability of life's emergence. The invigorated GARD scrutiny presented in this review enhances the validity of autocatalytic sets as a bona fide early evolution scenario and provides essential infrastructure for a paradigm shift towards a systems protobiology view of life's origin.