Energetics of amino acid synthesis in hydrothermal ecosystems

Solid-State Single-Molecule Sensing with the Electronic Life-Detection Instrument for Enceladus/Europa (ELIE)

Growing evidence of the potential habitability of Ocean Worlds across our solar system is motivating the advancement of technologies capable of detecting life as we know it-sharing a common ancestry or physicochemical origin with life on Earth-or don't know it, representing a distinct emergence of life different than our one known example. Here, we propose the Electronic Life-detection Instrument for Enceladus/Europa (ELIE), a solid-state single-molecule instrument payload that aims to search for life based on the detection of amino acids and informational polymers (IPs) at the parts per billion to trillion level. As a first proof-of-principle in a laboratory environment, we demonstrate the single-molecule detection of the amino acid L-proline at a 10 μM concentration in a compact system. Based on ELIE's solid-state quantum electronic tunneling sensing mechanism, we further propose the quantum property of the HOMO-LUMO gap (energy difference between a molecule's highest energy-occupied molecular orbital and lowest energy-unoccupied molecular orbital) as a novel metric to assess amino acid complexity. Finally, we assess the potential of ELIE to discriminate between abiotically and biotically derived α-amino acid abundance distributions to reduce the false positive risk for life detection. Nanogap technology can also be applied to the detection of nucleobases and short sequences of IPs such as, but not limited to, RNA and DNA. Future missions may utilize ELIE to target preserved biosignatures on the surface of Mars, extant life in its deep subsurface, or life or its biosignatures in a plume, surface, or subsurface of ice moons such as Enceladus or Europa. One-Sentence Summary: A solid-state nanogap can determine the abundance distribution of amino acids, detect nucleic acids, and shows potential for detecting life as we know it and life as we don't know it.

Evaluating the abiotic synthesis potential and the stability of building blocks of life beneath an impact-induced steam atmosphere

A prerequisite for prebiotic chemistry is the accumulation of critical building blocks of life. Some studies argue that more frequent impact events on the primitive Earth could have induced a more reducing steam atmosphere and thus favor widespread and more efficient synthesis of life building blocks. However, elevated temperature is also proposed to threaten the stability of organics and whether life building blocks could accumulate to appreciable levels in the reducing yet hot surface seawater beneath the steam atmosphere is still poorly examined. Here, we used a thermodynamic tool to examine the synthesis affinity of various life building blocks using inorganic gasses as reactants at elevated temperatures and corresponding steam pressures relevant with the steam-seawater interface. Our calculations show that although the synthesis affinity of all life building blocks decreases when temperature increases, many organics, including methane, methanol, and carboxylic acids, have positive synthesis affinity over a wide range of temperatures, implying that these species were favorable to form (>10-6 molal) in the surface seawater. However, cyanide and formaldehyde have overall negative affinities, suggesting that these critical compounds would tend to undergo hydrolysis in the surface seawaters. Most of the 18 investigated amino acids have positive affinities at temperature <220°C and their synthesis affinity increases under more alkaline conditions. Sugars, ribose, and nucleobases have overall negative synthesis affinities at the investigated range of temperatures. Synthesis affinities are shown to be sensitive to the hydrogen fugacity. Higher hydrogen fugacity (in equilibrium with FQI or IW) favors the synthesis and accumulation of nearly all the investigated compounds, except for HCN and its derivate products. In summary, our results suggest that reducing conditions induced by primitive impacts could indeed favor the synthesis/accumulation of some life building blocks, but some critical species, particularly HCN and nucleosides, were still unfavorable to accumulate to appreciable levels. Our results can provide helpful guidance for future efforts to search for or understand the stability of biomolecules on other planets like Mars and icy moons. We advocate examining craters formed by more reducing impactors to look for the preservation of prebiotic materials.

Heterotrophic Growth Dominates in the Most Extremotolerant Extremophile Cultures

Due to their ability to withstand "extreme" conditions, Earth's extremophilic organisms can constrain habitability windows for other planetary systems. However, there are many other considerations to microbial growth requirements beyond environmental extremes, such as nutrient availability. Here, we conduct a literature review of the most extremotolerant extremophiles in culture, since working with cultured organisms allows environmental and nutrient variables to be constrained with a high level of specificity. We generated a database that includes the isolation environment, carbon source(s) used, and growth preferences across temperature, pressure, salinity, and pH extremes. We found that the "most extreme" conditions were primarily sustained by heterotrophs, except for hyperthermophiles. These results highlight the importance of considering organic carbon availability when using extremophiles for habitability constraints. We also interrogated polyextreme potential across temperature, pressure, salinity, and pH conditions. Our findings suggest that the investigation of growth tolerance rather than growth optimum may reveal wider habitability parameters. Overall, these results highlight the potential polyextremes, environments, nutrient requirements, and additional analyses that could improve the application of cultured investigations to astrobiology questions.

Using thermodynamics to obtain geochemical information from genomes

Thermodynamic characterization of the relative stabilities of chemical compounds is a pillar of conceptual models in various fields of geosciences. Analogous models applied to genomes can yield new information about the relationship between genomes and their geochemical environments. In this perspective article, we present a chemical and thermodynamic analysis of prokaryotic lineages that have been the target of previous phylogenomic studies of evolutionary adaptation to varying redox conditions. The thermodynamic model development begins by quantifying the effects of hydrogen activity (aH₂ ) and temperature on the relative stabilities of organic compounds with different carbon oxidation state. When applied to proteins instead of metabolites, the same techniques can be used to identify combinations of aH₂ and temperature at which reference proteomes for Class I or Class II methanogens are relatively stable. The calculated aH₂ values are compatible with reported measurements for habitats of methanogens ranging from highly reducing submarine hydrothermal systems to less reducing environments including methanogenic sediments. In contrast to the transition between the two classes of methanogenic archaea, that between basal and terrestrial groups of Thaumarchaeota (denoting the origin of ammonia-oxidizing archaea) occurs at a less-reducing redox boundary. These examples reveal the consequences of energy minimization driving evolution and show how geochemical calculations involving biomolecules can be used to quantify and better understand the coevolution of the geosphere and biosphere.

Reconstruction and Analysis of Thermodynamically Constrained Models Reveal Metabolic Responses of a Deep-Sea Bacterium to Temperature Perturbations

Microbial acclimation to different temperature conditions can involve broad changes in cell composition and metabolic efficiency. A systems-level view of these metabolic responses in nonmesophilic organisms, however, is currently missing. In this study, thermodynamically constrained genome-scale models were applied to simulate the metabolic responses of a deep-sea psychrophilic bacterium, Shewanella psychrophila WP2, under suboptimal (4°C), optimal (15°C), and supraoptimal (20°C) growth temperatures. The models were calibrated with experimentally determined growth rates of WP2. Gibbs free energy change of reactions (Δ_rG'), metabolic fluxes, and metabolite concentrations were predicted using random simulations to characterize temperature-dependent changes in the metabolism. The modeling revealed the highest metabolic efficiency at the optimal temperature, and it suggested distinct patterns of ATP production and consumption that could lead to lower metabolic efficiency under suboptimal or supraoptimal temperatures. The modeling also predicted rearrangement of fluxes through multiple metabolic pathways, including the glycolysis pathway, Entner-Doudoroff pathway, tricarboxylic acid (TCA) cycle, and electron transport system, and these predictions were corroborated through comparisons to WP2 transcriptomes. Furthermore, predictions of metabolite concentrations revealed the potential conservation of reducing equivalents and ATP in the suboptimal temperature, consistent with experimental observations from other psychrophiles. Taken together, the WP2 models provided mechanistic insights into the metabolism of a psychrophile in response to different temperatures. IMPORTANCE Metabolic flexibility is a central component of any organism's ability to survive and adapt to changes in environmental conditions. This study represents the first application of thermodynamically constrained genome-scale models in simulating the metabolic responses of a deep-sea psychrophilic bacterium to various temperatures. The models predicted differences in metabolic efficiency that were attributed to changes in metabolic pathway utilization and metabolite concentration during growth under optimal and nonoptimal temperatures. Experimental growth measurements were used for model calibration, and temperature-dependent transcriptomic changes corroborated the model-predicted rearrangement of metabolic fluxes. Overall, this study highlights the utility of modeling approaches in studying the temperature-driven metabolic responses of an extremophilic organism.

Metabolic Strategies Shared by Basement Residents of the Lost City Hydrothermal Field

Alkaline fluids venting from chimneys of the Lost City hydrothermal field flow from a potentially vast microbial habitat within the seafloor where energy and organic molecules are released by chemical reactions within rocks uplifted from Earth's mantle. In this study, we investigated hydrothermal fluids venting from Lost City chimneys as windows into subseafloor environments where the products of geochemical reactions, such as molecular hydrogen (H₂), formate, and methane, may be the only available sources of energy for biological activity. Our deep sequencing of metagenomes and metatranscriptomes from these hydrothermal fluids revealed a few key species of archaea and bacteria that are likely to play critical roles in the subseafloor microbial ecosystem. We identified a population of Thermodesulfovibrionales (belonging to phylum Nitrospirota) as a prevalent sulfate-reducing bacterium that may be responsible for much of the consumption of H₂ and sulfate in Lost City fluids. Metagenome-assembled genomes (MAGs) classified as Methanosarcinaceae and Candidatus Bipolaricaulota were also recovered from venting fluids and represent potential methanogenic and acetogenic members of the subseafloor ecosystem. These genomes share novel hydrogenases and formate dehydrogenase-like sequences that may be unique to hydrothermal environments where H₂ and formate are much more abundant than carbon dioxide. The results of this study include multiple examples of metabolic strategies that appear to be advantageous in hydrothermal and subsurface alkaline environments where energy and carbon are provided by geochemical reactions. IMPORTANCE The Lost City hydrothermal field is an iconic example of a microbial ecosystem fueled by energy and carbon from Earth's mantle. Uplift of mantle rocks into the seafloor can trigger a process known as serpentinization that releases molecular hydrogen (H₂) and creates unusual environmental conditions where simple organic carbon molecules are more stable than dissolved inorganic carbon. This study provides an initial glimpse into the kinds of microbes that live deep within the seafloor where serpentinization takes place, by sampling hydrothermal fluids exiting from the Lost City chimneys. The metabolic strategies that these microbes appear to be using are also shared by microbes that inhabit other sites of serpentinization, including continental subsurface environments and natural springs. Therefore, the results of this study contribute to a broader, interdisciplinary effort to understand the general principles and mechanisms by which serpentinization-associated processes can support life on Earth and perhaps other worlds.

Energy at Origins: Favorable Thermodynamics of Biosynthetic Reactions in the Last Universal Common Ancestor (LUCA)

Though all theories for the origin of life require a source of energy to promote primordial chemical reactions, the nature of energy that drove the emergence of metabolism at origins is still debated. We reasoned that evidence for the nature of energy at origins should be preserved in the biochemical reactions of life itself, whereby changes in free energy, ΔG, which determine whether a reaction can go forward or not, should help specify the source. By calculating values of ΔG across the conserved and universal core of 402 individual reactions that synthesize amino acids, nucleotides and cofactors from H2, CO2, NH3, H2S and phosphate in modern cells, we find that 95-97% of these reactions are exergonic (ΔG ≤ 0 kJ⋅mol-1) at pH 7-10 and 80-100°C under nonequilibrium conditions with H2 replacing biochemical reductants. While 23% of the core's reactions involve ATP hydrolysis, 77% are ATP-independent, thermodynamically driven by ΔG of reactions involving carbon bonds. We identified 174 reactions that are exergonic by -20 to -300 kJ⋅mol-1 at pH 9 and 80°C and that fall into ten reaction types: six pterin dependent alkyl or acyl transfers, ten S-adenosylmethionine dependent alkyl transfers, four acyl phosphate hydrolyses, 14 thioester hydrolyses, 30 decarboxylations, 35 ring closure reactions, 31 aromatic ring formations, and 44 carbon reductions by reduced nicotinamide, flavins, ferredoxin, or formate. The 402 reactions of the biosynthetic core trace to the last universal common ancestor (LUCA), and reveal that synthesis of LUCA's chemical constituents required no external energy inputs such as electric discharge, UV-light or phosphide minerals. The biosynthetic reactions of LUCA uncover a natural thermodynamic tendency of metabolism to unfold from energy released by reactions of H2, CO2, NH3, H2S, and phosphate.

Factors in Protobiomonomer Selection for the Origin of the Standard Genetic Code

Natural selection of specific protobiomonomers during abiogenic development of the prototype genetic code is hindered by the diversity of structural, spatial, and rotational isomers that have identical elemental composition and molecular mass (M), but can vary significantly in their physicochemical characteristics, such as the melting temperature T_m, the T_m:M ratio, and the solubility in water, due to different positions of atoms in the molecule. These parameters differ between cis- and trans-isomers of dicarboxylic acids, spatial monosaccharide isomers, and structural isomers of α-, β-, and γ-amino acids. The stable planar heterocyclic molecules of the major nucleobases comprise four (C, H, N, O) or three (C, H, N) elements and contain a single -C=C bond and two nitrogen atoms in each heterocycle involved in C-N and C=N bonds. They exist as isomeric resonance hybrids of single and double bonds and as a mixture of tautomer forms due to the presence of -C=O and/or -NH₂ side groups. They are thermostable, insoluble in water, and exhibit solid-state stability, which is of central importance for DNA molecules as carriers of genetic information. In M-T_m diagrams, proteinogenic amino acids and the corresponding codons are distributed fairly regularly relative to the distinct clusters of purine and pyrimidine bases, reflecting the correspondence between codons and amino acids that was established in different periods of genetic code development. The body of data on the evolution of the genetic code system indicates that the elemental composition and molecular structure of protobiomonomers, and their M, T_m, photostability, and aqueous solubility determined their selection in the emergence of the standard genetic code.

Plausible Emergence and Self Assembly of a Primitive Phospholipid from Reduced Phosphorus on the Primordial Earth

How life arose on the primitive Earth is one of the biggest questions in science. Biomolecular emergence scenarios have proliferated in the literature but accounting for the ubiquity of oxidized (+ 5) phosphate (PO₄^3-) in extant biochemistries has been challenging due to the dearth of phosphate and molecular oxygen on the primordial Earth. A compelling body of work suggests that exogenous schreibersite ((Fe,Ni)₃P) was delivered to Earth via meteorite impacts during the Heavy Bombardment (ca. 4.1-3.8 Gya) and there converted to reduced P oxyanions (e.g., phosphite (HPO₃^2-) and hypophosphite (H₂PO₂^-)) and phosphonates. Inspired by this idea, we review the relevant literature to deduce a plausible reduced phospholipid analog of modern phosphatidylcholines that could have emerged in a primordial hydrothermal setting. A shallow alkaline lacustrine basin underlain by active hydrothermal fissures and meteoritic schreibersite-, clay-, and metal-enriched sediments is envisioned. The water column is laden with known and putative primordial hydrothermal reagents. Small system dimensions and thermal- and UV-driven evaporation further concentrate chemical precursors. We hypothesize that a reduced phospholipid arises from Fischer-Tropsch-type (FTT) production of a C8 alkanoic acid, which condenses with an organophosphinate (derived from schreibersite corrosion to hypophosphite with subsequent methylation/oxidation), to yield a reduced protophospholipid. This then condenses with an α-amino nitrile (derived from Strecker-type reactions) to form the polar head. Preliminary modeling results indicate that reduced phospholipids do not aggregate rapidly; however, single layer micelles are stable up to aggregates with approximately 100 molecules.

Carbon Oxidation State in Microbial Polar Lipids Suggests Adaptation to Hot Spring Temperature and Redox Gradients

The influence of oxidation-reduction (redox) potential on the expression of biomolecules is a topic of ongoing exploration in geobiology. In this study, we investigate the novel possibility that structures and compositions of lipids produced by microbial communities are sensitive to environmental redox conditions. We extracted lipids from microbial biomass collected along the thermal and redox gradients of four alkaline hot springs in Yellowstone National Park (YNP) and investigated patterns in the average oxidation state of carbon (Z_C), a metric calculated from the chemical formulae of lipid structures. Carbon in intact polar lipids (IPLs) and their alkyl chains becomes more oxidized (higher Z_C) with increasing distance from each of the four hot spring sources. This coincides with decreased water temperature and increased concentrations of oxidized inorganic solutes, such as dissolved oxygen, sulfate, and nitrate. Carbon in IPLs is most reduced (lowest Z_C) in the hot, reduced conditions upstream, with abundance-weighted Z_C values between −1.68 and −1.56. These values increase gradually downstream to around −1.36 to −1.33 in microbial communities living between 29.0 and 38.1°C. This near-linear increase in Z_C can be attributed to a shift from ether-linked to ester-linked alkyl chains, a decrease in average aliphatic carbons per chain (nC), an increase in average degree of unsaturation per chain (nUnsat), and increased cyclization in tetraether lipids. The Z_C of lipid headgroups and backbones did not change significantly downstream. Expression of lipids with relatively reduced carbon under reduced conditions and oxidized lipids under oxidized conditions may indicate microbial adaptation across environmental gradients in temperature and electron donor/acceptor supply.

Possible mechanisms of CO2 reduction by H2 via prebiotic vectorial electrochemistry

Methanogens are putatively ancestral autotrophs that reduce CO2 with H2 to form biomass using a membrane-bound, proton-motive Fe(Ni)S protein called the energy-converting hydrogenase (Ech). At the origin of life, geologically sustained H+ gradients across inorganic barriers containing Fe(Ni)S minerals could theoretically have driven CO2 reduction by H2 through vectorial chemistry in a similar way to Ech. pH modulation of the redox potentials of H2, CO2 and Fe(Ni)S minerals could in principle enable an otherwise endergonic reaction. Here, we analyse whether vectorial electrochemistry can facilitate the reduction of CO2 by H2 under alkaline hydrothermal conditions using a microfluidic reactor. We present pilot data showing that steep pH gradients of approximately 5 pH units can be sustained over greater than 5 h across Fe(Ni)S barriers, with H+-flux across the barrier about two million-fold faster than OH--flux. This high flux produces a calculated 3-pH unit-gradient (equating to 180 mV) across single approximately 25-nm Fe(Ni)S nanocrystals, which is close to that required to reduce CO2. However, the poor solubility of H2 at atmospheric pressure limits CO2 reduction by H2, explaining why organic synthesis has so far proved elusive in our reactor. Higher H2 concentration will be needed in future to facilitate CO2 reduction through prebiotic vectorial electrochemistry.

The fitness of chemotrophs increases when their catabolic by-products are consumed by other species

Chemotrophic microorganisms synthesise biomass by utilising energy obtained from a set of chemical reactions that convert resources to by-products, forming catabolic interactions. The amount of energy obtained per catabolic reaction decreases with the abundance of the by-product named as the 'abundant resource premium'. Consider two species, Species 1 and 2, Species 1 obtains energy from a reaction that converts resource A to by-product B. Species 2 then utilises B as its resource, extracting energy from a reaction that converts B to C. Thus, the presence of Species 2 reduces the abundance of B, which improves the fitness of Species 1 by increasing the energy acquisition per reaction of A to B. We discuss the population dynamic implication of this effect and its importance in expanding a realised niche, boosting material flow through the ecosystem and providing mutualistic interactions among species linked by the material flow. Introducing thermodynamics into population ecology could offer us fundamental ecological insights into understanding the ecology of chemotrophic microorganisms dominating the subsurface realm.

Deciphering Biosignatures in Planetary Contexts

Microbial life permeates Earth's critical zone and has likely inhabited nearly all our planet's surface and near subsurface since before the beginning of the sedimentary rock record. Given the vast time that Earth has been teeming with life, do astrobiologists truly understand what geological features untouched by biological processes would look like? In the search for extraterrestrial life in the Universe, it is critical to determine what constitutes a biosignature across multiple scales, and how this compares with "abiosignatures" formed by nonliving processes. Developing standards for abiotic and biotic characteristics would provide quantitative metrics for comparison across different data types and observational time frames. The evidence for life detection falls into three categories of biosignatures: (1) substances, such as elemental abundances, isotopes, molecules, allotropes, enantiomers, minerals, and their associated properties; (2) objects that are physical features such as mats, fossils including trace-fossils and microbialites (stromatolites), and concretions; and (3) patterns, such as physical three-dimensional or conceptual n-dimensional relationships of physical or chemical phenomena, including patterns of intermolecular abundances of organic homologues, and patterns of stable isotopic abundances between and within compounds. Five key challenges that warrant future exploration by the astrobiology community include the following: (1) examining phenomena at the "right" spatial scales because biosignatures may elude us if not examined with the appropriate instrumentation or modeling approach at that specific scale; (2) identifying the precise context across multiple spatial and temporal scales to understand how tangible biosignatures may or may not be preserved; (3) increasing capability to mine big data sets to reveal relationships, for example, how Earth's mineral diversity may have evolved in conjunction with life; (4) leveraging cyberinfrastructure for data management of biosignature types, characteristics, and classifications; and (5) using three-dimensional to n-D representations of biotic and abiotic models overlain on multiple overlapping spatial and temporal relationships to provide new insights.

Nonequilibrium abundances for the building blocks of life

The difficulty of obtaining appreciable quantities of biologically important molecules in thermodynamic equilibrium has long been identified as an obstacle to life's emergence, and determining the specific nonequilibrium conditions that might have given rise to life is challenging. To address these issues, we investigate how the concentrations of life's building blocks change as a function of the distance from equilibrium on average, in two example settings: (i) the synthesis of heavy amino acids and (ii) their polymerization into peptides. We find that relative concentrations of the heaviest amino acids can be boosted by four orders of magnitude, and concentrations of the longest peptide chains can be increased by hundreds of orders of magnitude. The average nonequilibrium distribution does not depend on the details of how the system was driven from equilibrium, indicating that environments might not have to be fine-tuned to support life.

The Cell and the Sum of Its Parts: Patterns of Complexity in Biosignatures as Revealed by Deep UV Raman Spectroscopy

The next NASA-led Mars mission (Mars 2020) will carry a suite of instrumentation dedicated to investigating Martian history and the in situ detection of potential biosignatures. SHERLOC, a deep UV Raman/Fluorescence spectrometer has the ability to detect and map the distribution of many organic compounds, including the aromatic molecules that are fundamental building blocks of life on Earth, at concentrations down to 1 ppm. The mere presence of organic compounds is not a biosignature: there is widespread distribution of reduced organic molecules in the Solar System. Life utilizes a select few of these molecules creating conspicuous enrichments of specific molecules that deviate from the distribution expected from purely abiotic processes. The detection of far from equilibrium concentrations of a specific subset of organic molecules, such as those uniquely enriched by biological processes, would comprise a universal biosignature independent of specific terrestrial biochemistry. The detectability and suitability of a small subset of organic molecules to adequately describe a living system is explored using the bacterium Escherichia coli as a model organism. The DUV Raman spectra of E. coli cells are dominated by the vibrational modes of the nucleobases adenine, guanine, cytosine, and thymine, and the aromatic amino acids tyrosine, tryptophan, and phenylalanine. We demonstrate that not only does the deep ultraviolet (DUV) Raman spectrum of E. coli reflect a distinct concentration of specific organic molecules, but that a sufficient molecular complexity is required to deconvolute the cellular spectrum. Furthermore, a linear combination of the DUV resonant compounds is insufficient to fully describe the cellular spectrum. The residual in the cellular spectrum indicates that DUV Raman spectroscopy enables differentiating between the presence of biomolecules and the complex uniquely biological organization and arrangements of these molecules in living systems. This study demonstrates the ability of DUV Raman spectroscopy to interrogate a complex biological system represented in a living cell, and differentiate between organic detection and a series of Raman features that derive from the molecular complexity inherent to life constituting a biosignature.

Abiotic synthesis of amino acids in the recesses of the oceanic lithosphere

Abiotic hydrocarbons and carboxylic acids are known to be formed on Earth, notably during the hydrothermal alteration of mantle rocks. Although the abiotic formation of amino acids has been predicted both from experimental studies and thermodynamic calculations, its occurrence has not been demonstrated in terrestrial settings. Here, using a multimodal approach that combines high-resolution imaging techniques, we obtain evidence for the occurrence of aromatic amino acids formed abiotically and subsequently preserved at depth beneath the Atlantis Massif (Mid-Atlantic Ridge). These aromatic amino acids may have been formed through Friedel-Crafts reactions catalysed by an iron-rich saponite clay during a late alteration stage of the massif serpentinites. Demonstrating the potential of fluid-rock interactions in the oceanic lithosphere to generate amino acids abiotically gives credence to the hydrothermal theory for the origin of life, and may shed light on ancient metabolisms and the functioning of the present-day deep biosphere.

Thermodynamic favorability and pathway yield as evolutionary tradeoffs in biosynthetic pathway choice

The structure of the metabolic network contains myriad organism-specific variations across the tree of life, but the selection basis for pathway choices in different organisms is not well understood. Here, we examined the metabolic capabilities with respect to cofactor use and pathway thermodynamics of all sequenced organisms in the Kyoto Encyclopedia of Genes and Genomes Database. We found that (i) many biomass precursors have alternate synthesis routes that vary substantially in thermodynamic favorability and energy cost, creating tradeoffs that may be subject to selection pressure; (ii) alternative pathways in amino acid synthesis are characteristically distinguished by the use of biosynthetically unnecessary acyl-CoA cleavage; (iii) distinct choices preferring thermodynamic-favorable or cofactor-use-efficient pathways exist widely among organisms; (iv) cofactor-use-efficient pathways tend to have a greater yield advantage under anaerobic conditions specifically; and (v) lysine biosynthesis in particular exhibits temperature-dependent thermodynamics and corresponding differential pathway choice by thermophiles. These findings present a view on the evolution of metabolic network structure that highlights a key role of pathway thermodynamics and cofactor use in determining organism pathway choices.

Organic synthesis on Mars by electrochemical reduction of CO2

The sources and nature of organic carbon on Mars have been a subject of intense research. Steele et al. (2012) showed that 10 martian meteorites contain macromolecular carbon phases contained within pyroxene- and olivine-hosted melt inclusions. Here, we show that martian meteorites Tissint, Nakhla, and NWA 1950 have an inventory of organic carbon species associated with fluid-mineral reactions that are remarkably consistent with those detected by the Mars Science Laboratory (MSL) mission. We advance the hypothesis that interactions among spinel-group minerals, sulfides, and a brine enable the electrochemical reduction of aqueous CO₂ to organic molecules. Although documented here in martian samples, a similar process likely occurs wherever igneous rocks containing spinel-group minerals and/or sulfides encounter brines.

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.

Green Rust: The Simple Organizing 'Seed' of All Life?

Korenaga and coworkers presented evidence to suggest that the Earth's mantle was dry and water filled the ocean to twice its present volume 4.3 billion years ago. Carbon dioxide was constantly exhaled during the mafic to ultramafic volcanic activity associated with magmatic plumes that produced the thick, dense, and relatively stable oceanic crust. In that setting, two distinct and major types of sub-marine hydrothermal vents were active: ~400 °C acidic springs, whose effluents bore vast quantities of iron into the ocean, and ~120 °C, highly alkaline, and reduced vents exhaling from the cooler, serpentinizing crust some distance from the heads of the plumes. When encountering the alkaline effluents, the iron from the plume head vents precipitated out, forming mounds likely surrounded by voluminous exhalative deposits similar to the banded iron formations known from the Archean. These mounds and the surrounding sediments, comprised micro or nano-crysts of the variable valence FeII/FeIII oxyhydroxide known as green rust. The precipitation of green rust, along with subsidiary iron sulfides and minor concentrations of nickel, cobalt, and molybdenum in the environment at the alkaline springs, may have established both the key bio-syntonic disequilibria and the means to properly make use of them-the elements needed to effect the essential inanimate-to-animate transitions that launched life. Specifically, in the submarine alkaline vent model for the emergence of life, it is first suggested that the redox-flexible green rust micro- and nano-crysts spontaneously precipitated to form barriers to the complete mixing of carbonic ocean and alkaline hydrothermal fluids. These barriers created and maintained steep ionic disequilibria. Second, the hydrous interlayers of green rust acted as engines that were powered by those ionic disequilibria and drove essential endergonic reactions. There, aided by sulfides and trace elements acting as catalytic promoters and electron transfer agents, nitrate could be reduced to ammonia and carbon dioxide to formate, while methane may have been oxidized to methyl and formyl groups. Acetate and higher carboxylic acids could then have been produced from these C1 molecules and aminated to amino acids, and thence oligomerized to offer peptide nests to phosphate and iron sulfides, and secreted to form primitive amyloid-bounded structures, leading conceivably to protocells.

Temperature-Dependent Estimation of Gibbs Energies Using an Updated Group-Contribution Method

Reaction-equilibrium constants determine the metabolite concentrations necessary to drive flux through metabolic pathways. Group-contribution methods offer a way to estimate reaction-equilibrium constants at wide coverage across the metabolic network. Here, we present an updated group-contribution method with 1) additional curated thermodynamic data used in fitting and 2) capabilities to calculate equilibrium constants as a function of temperature. We first collected and curated aqueous thermodynamic data, including reaction-equilibrium constants, enthalpies of reaction, Gibbs free energies of formation, enthalpies of formation, entropy changes of formation of compounds, and proton- and metal-ion-binding constants. Next, we formulated the calculation of equilibrium constants as a function of temperature and calculated the standard entropy change of formation (ΔfS∘) using a model based on molecular properties. The median absolute error in estimating ΔfS∘ was 0.013 kJ/K/mol. We also estimated magnesium binding constants for 618 compounds using a linear regression model validated against measured data. We demonstrate the improved performance of the current method (8.17 kJ/mol in median absolute residual) over the current state-of-the-art method (11.47 kJ/mol) in estimating the 185 new reactions added in this work. The efforts here fill in gaps for thermodynamic calculations under various conditions, specifically different temperatures and metal-ion concentrations. These, to our knowledge, new capabilities empower the study of thermodynamic driving forces underlying the metabolic function of organisms living under diverse conditions.

Acetyl Phosphate as a Primordial Energy Currency at the Origin of Life

Metabolism is primed through the formation of thioesters via acetyl CoA and the phosphorylation of substrates by ATP. Prebiotic equivalents such as methyl thioacetate and acetyl phosphate have been proposed to catalyse analogous reactions at the origin of life, but their propensity to hydrolyse challenges this view. Here we show that acetyl phosphate (AcP) can be synthesised in water within minutes from thioacetate (but not methyl thioacetate) under ambient conditions. AcP is stable over hours, depending on temperature, pH and cation content, giving it an ideal poise between stability and reactivity. We show that AcP can phosphorylate nucleotide precursors such as ribose to ribose-5-phosphate and adenosine to adenosine monophosphate, at modest (~2%) yield in water, and at a range of pH. AcP can also phosphorylate ADP to ATP in water over several hours at 50 °C. But AcP did not promote polymerization of either glycine or AMP. The amino group of glycine was preferentially acetylated by AcP, especially at alkaline pH, hindering the formation of polypeptides. AMP formed small stacks of up to 7 monomers, but these did not polymerise in the presence of AcP in aqueous solution. We conclude that AcP can phosphorylate biologically meaningful substrates in a manner analogous to ATP, promoting the origins of metabolism, but is unlikely to have driven polymerization of macromolecules such as polypeptides or RNA in free solution. This is consistent with the idea that a period of monomer (cofactor) catalysis preceded the emergence of polymeric enzymes or ribozymes at the origin of life.

Chance, necessity and the origins of life: a physical sciences perspective

Earth's 4.5-billion-year history has witnessed a complex sequence of high-probability chemical and physical processes, as well as 'frozen accidents'. Most models of life's origins similarly invoke a sequence of chemical reactions and molecular self-assemblies in which both necessity and chance play important roles. Recent research adds two important insights into this discussion. First, in the context of chemical reactions, chance versus necessity is an inherently false dichotomy-a range of probabilities exists for many natural events. Second, given the combinatorial richness of early Earth's chemical and physical environments, events in molecular evolution that are unlikely at limited laboratory scales of space and time may, nevertheless, be inevitable on an Earth-like planet at time scales of a billion years.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Bacterial Diversity and Biogeochemistry of Two Marine Shallow-Water Hydrothermal Systems off Dominica (Lesser Antilles)

Shallow-water hydrothermal systems represent extreme environments with unique biogeochemistry and high biological productivity, at which autotrophic microorganisms use both light and chemical energy for the production of biomass. Microbial communities of these ecosystems are metabolically diverse and possess the capacity to transform a large range of chemical compounds. Yet, little is known about their diversity or factors shaping their structure or how they compare to coastal sediments not impacted by hydrothermalism. To this end, we have used automated ribosomal intergenic spacer analysis (ARISA) and high-throughput Illumina sequencing combined with porewater geochemical analysis to investigate microbial communities along geochemical gradients in two shallow-water hydrothermal systems off the island of Dominica (Lesser Antilles). At both sites, venting of hydrothermal fluids substantially altered the porewater geochemistry by enriching it with silica, iron and dissolved inorganic carbon, resulting in island-like habitats with distinct biogeochemistry. The magnitude of fluid flow and difference in sediment grain size, which impedes mixing of the fluids with seawater, were correlated with the observed differences in the porewater geochemistry between the two sites. Concomitantly, individual sites harbored microbial communities with a significantly different community structure. These differences could be statistically linked to variations in the porewater geochemistry and the hydrothermal fluids. The two shallow-water hydrothermal systems of Dominica harbored bacterial communities with high taxonomical and metabolic diversity, predominated by heterotrophic microorganisms associated with the Gammaproteobacterial genera Pseudomonas and Pseudoalteromonas, indicating the importance of heterotrophic processes. Overall, this study shows that shallow-water hydrothermal systems contribute substantially to the biogeochemical heterogeneity and bacterial diversity of coastal sediments.

Origins of life: a problem for physics, a key issues review

The origins of life stands among the great open scientific questions of our time. While a number of proposals exist for possible starting points in the pathway from non-living to living matter, these have so far not achieved states of complexity that are anywhere near that of even the simplest living systems. A key challenge is identifying the properties of living matter that might distinguish living and non-living physical systems such that we might build new life in the lab. This review is geared towards covering major viewpoints on the origin of life for those new to the origin of life field, with a forward look towards considering what it might take for a physical theory that universally explains the phenomenon of life to arise from the seemingly disconnected array of ideas proposed thus far. The hope is that a theory akin to our other theories in fundamental physics might one day emerge to explain the phenomenon of life, and in turn finally permit solving its origins.

Relative Importance of Chemoautotrophy for Primary Production in a Light Exposed Marine Shallow Hydrothermal System

The unique geochemistry of marine shallow-water hydrothermal systems promotes the establishment of diverse microbial communities with a range of metabolic pathways. In contrast to deep-sea vents, shallow-water vents not only support chemosynthesis, but also phototrophic primary production due to the availability of light. However, comprehensive studies targeting the predominant biogeochemical processes are rare, and consequently a holistic understanding of the functioning of these ecosystems is currently lacking. To this end, we combined stable isotope probing of lipid biomarkers with an analysis of the bacterial communities to investigate if chemoautotrophy, in parallel to photoautotrophy, plays an important role in autotrophic carbon fixation and to identify the key players. The study was carried out at a marine shallow-water hydrothermal system located at 5 m water depth off Dominica Island (Lesser Antilles), characterized by up to 55°C warm hydrothermal fluids that contain high amounts of dissolved Fe²⁺. Analysis of the bacterial diversity revealed Anaerolineae of the Chloroflexi as the most abundant bacterial class. Furthermore, the presence of key players involved in iron cycling generally known from deep-sea hydrothermal vents (e.g., Zetaproteobacteria and Geothermobacter), supported the importance of iron-driven redox processes in this hydrothermal system. Uptake of ¹³C-bicarbonate into bacterial fatty acids under light and dark conditions revealed active photo- and chemoautotrophic communities, with chemoautotrophy accounting for up to 65% of the observed autotrophic carbon fixation. Relatively increased ¹³C-incorporation in the dark allowed the classification of aiC_15:0, C_15:0, and iC_16:0 as potential lipid biomarkers for bacterial chemoautotrophy in this ecosystem. Highest total ¹³C-incorporation into fatty acids took place at the sediment surface, but chemosynthesis was found to be active down to 8 cm sediment depth. In conclusion, this study highlights the relative importance of chemoautotrophy compared to photoautotrophy in a shallow-water hydrothermal system, emphasizing chemosynthesis as a prominent process for biomass production in marine coastal environments influenced by hydrothermalism.

"Hot" acetogenesis

Thermophilic microorganisms as well as acetogenic bacteria are both considered ancient. Interestingly, only a few species of bacteria, all belonging to the family Thermoanaerobacteraceae, are described to conserve energy from acetate formation with hydrogen as electron donor and carbon dioxide as electron acceptor. This review reflects the metabolic differences between Moorella spp., Thermoanaerobacter kivui and Thermacetogenium phaeum, with focus on the biochemistry of autotrophic growth and energy conservation. The potential of these thermophilic acetogens for biotechnological applications is discussed briefly.

Proteomic indicators of oxidation and hydration state in colorectal cancer

New integrative approaches are needed to harness the potential of rapidly growing datasets of protein expression and microbial community composition in colorectal cancer. Chemical and thermodynamic models offer theoretical tools to describe populations of biomacromolecules and their relative potential for formation in different microenvironmental conditions. The average oxidation state of carbon (Z_C) can be calculated as an elemental ratio from the chemical formulas of proteins, and water demand per residue (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\overline{n}}_{{\mathrm{H}}_{2}\mathrm{O}}$\end{document}n¯H2O) is computed by writing the overall formation reactions of proteins from basis species. Using results reported in proteomic studies of clinical samples, many datasets exhibit higher mean Z_C or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\overline{n}}_{{\mathrm{H}}_{2}\mathrm{O}}$\end{document}n¯H2O of proteins in carcinoma or adenoma compared to normal tissue. In contrast, average protein compositions in bacterial genomes often have lower Z_C for bacteria enriched in fecal samples from cancer patients compared to healthy donors. In thermodynamic calculations, the potential for formation of the cancer-related proteins is energetically favored by changes in the chemical activity of H₂O and fugacity of O₂ that reflect the compositional differences. The compositional analysis suggests that a systematic change in chemical composition is an essential feature of cancer proteomes, and the thermodynamic descriptions show that the observed proteomic transformations in host tissue could be promoted by relatively high microenvironmental oxidation and hydration states.

The energetics of anabolism in natural settings

The environmental conditions that describe an ecosystem define the amount of energy available to the resident organisms and the amount of energy required to build biomass. Here, we quantify the amount of energy required to make biomass as a function of temperature, pressure, redox state, the sources of C, N and S, cell mass and the time that an organism requires to double or replace its biomass. Specifically, these energetics are calculated from 0 to 125 °C, 0.1 to 500 MPa and -0.38 to +0.86 V using CO2, acetate or CH4 for C, NO3(-) or NH4(+) for N and SO4(2-) or HS(-) for S. The amounts of energy associated with synthesizing the biomolecules that make up a cell, which varies over 39 kJ (g cell)(-1), are then used to compute energy-based yield coefficients for a vast range of environmental conditions. Taken together, environmental variables and the range of cell sizes leads to a ~4 orders of magnitude difference between the number of microbial cells that can be made from a Joule of Gibbs energy under the most (5.06 × 10(11) cells J(-1)) and least (5.21 × 10(7) cells J(-1)) ideal conditions. When doubling/replacement time is taken into account, the range of anabolism energies can expand even further.

One step beyond a ribosome: The ancient anaerobic core

Life arose in a world without oxygen and the first organisms were anaerobes. Here we investigate the gene repertoire of the prokaryote common ancestor, estimating which genes it contained and to which lineages of modern prokaryotes it was most similar in terms of gene content. Using a phylogenetic approach we found that among trees for all 8779 protein families shared between 134 archaea and 1847 bacterial genomes, only 1045 have sequences from at least two bacterial and two archaeal groups and retain the ancestral archaeal–bacterial split. Among those, the genes shared by anaerobes were identified as candidate genes for the prokaryote common ancestor, which lived in anaerobic environments. We find that these anaerobic prokaryote common ancestor genes are today most frequently distributed among methanogens and clostridia, strict anaerobes that live from low free energy changes near the thermodynamic limit of life. The anaerobic families encompass genes for bifunctional acetyl-CoA-synthase/CO-dehydrogenase, heterodisulfide reductase subunits C and A, ferredoxins, and several subunits of the Mrp-antiporter/hydrogenase family, in addition to numerous S-adenosyl methionine (SAM) dependent methyltransferases. The data indicate a major role for methyl groups in the metabolism of the prokaryote common ancestor. The data furthermore indicate that the prokaryote ancestor possessed a rotor stator ATP synthase, but lacked cytochromes and quinones as well as identifiable redox-dependent ion pumping complexes. The prokaryote ancestor did possess, however, an Mrp-type H⁺/Na⁺ antiporter complex, capable of transducing geochemical pH gradients into biologically more stable Na⁺-gradients. The findings implicate a hydrothermal, autotrophic, and methyl-dependent origin of life. This article is part of a Special Issue entitled ‘EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2–6, 2016’, edited by Prof. Paolo Bernardi.

•

Life arose without oxygen, the universal ancestor (Luca) was an anaerobe.

•

We used phylogenetic and physiological criteria to identify genes present in Luca.

•

An ancient core of 65 metabolic genes shed light on Luca's anaerobic lifestyle.

•

Ancient core genes are most widespread among modern methanogens and clostridia.

•

The data implicate a major role for methyl groups in Luca's anaerobic metabolism.

Early Microbial Evolution: The Age of Anaerobes

In this article, the term "early microbial evolution" refers to the phase of biological history from the emergence of life to the diversification of the first microbial lineages. In the modern era (since we knew about archaea), three debates have emerged on the subject that deserve discussion: (1) thermophilic origins versus mesophilic origins, (2) autotrophic origins versus heterotrophic origins, and (3) how do eukaryotes figure into early evolution. Here, we revisit those debates from the standpoint of newer data. We also consider the perhaps more pressing issue that molecular phylogenies need to recover anaerobic lineages at the base of prokaryotic trees, because O2 is a product of biological evolution; hence, the first microbes had to be anaerobes. If molecular phylogenies do not recover anaerobes basal, something is wrong. Among the anaerobes, hydrogen-dependent autotrophs--acetogens and methanogens--look like good candidates for the ancestral state of physiology in the bacteria and archaea, respectively. New trees tend to indicate that eukaryote cytosolic ribosomes branch within their archaeal homologs, not as sisters to them and, furthermore tend to root archaea within the methanogens. These are major changes in the tree of life, and open up new avenues of thought. Geochemical methane synthesis occurs as a spontaneous, abiotic exergonic reaction at hydrothermal vents. The overall similarity between that reaction and biological methanogenesis fits well with the concept of a methanogenic root for archaea and an autotrophic origin of microbial physiology.

Energetics of Amino Acid Synthesis in Alkaline Hydrothermal Environments

Alkaline hydrothermal systems have received considerable attention as candidates for the origin and evolution of life on the primitive Earth. Nevertheless, sufficient information has not yet been obtained for the thermodynamic properties of amino acids, which are necessary components for life, at high temperatures and alkaline pH. These properties were estimated using experimental high-temperature volume and heat capacity data reported in the literature for several amino acids, together with correlation algorithms and the revised Helgeson-Kirkham-Flowers (HKF) equations of state. This approach enabled determination of a complete set of the standard molal thermodynamic data and the revised HKF parameters for the 20 protein amino acids in their zwitterionic and ionization states. The obtained dataset was then used to evaluate the energetics of amino acid syntheses from simple inorganic precursors (CO2, H2, NH3 and H2S) in a simulated alkaline hydrothermal system on the Hadean Earth. Results show that mixing between CO2-rich seawater and the H2-rich hydrothermal fluid can produce energetically favorable conditions for amino acid syntheses, particularly in the lower-temperature region of such systems. Together with data related to the pH and temperature dependences of the energetics of amino acid polymerizations presented in earlier reports, these results suggest the following. Hadean alkaline hydrothermal settings, where steep pH and temperature gradients may have existed between cool, slightly acidic Hadean ocean water and hot, alkaline hydrothermal fluids at the vent-ocean interface, may be energetically the most suitable environment for the synthesis and polymerization of amino acids.

Oxidative Self-Polymerization of Dopamine in an Acidic Environment

A weak alkaline condition (pH > 8) is a general requirement for oxidative self-polymerization of dopamine. Here, we first demonstrated the generation of polydopamine in an acidic environment via a hydrothermal method. The pH scope of self-polymerization of dopamine is extended to pH ∼ 1 in a hydrothermal process. Polydopamine generated via a hydrothermal method shows similar chemical features and radical scavenging activity with that generated in a basic environment.

Average oxidation state of carbon in proteins

The formal oxidation state of carbon atoms in organic molecules depends on the covalent structure. In proteins, the average oxidation state of carbon (Z(C)) can be calculated as an elemental ratio from the chemical formula. To investigate oxidation-reduction (redox) patterns, groups of proteins from different subcellular locations and phylogenetic groups were selected for comparison. Extracellular proteins of yeast have a relatively high oxidation state of carbon, corresponding with oxidizing conditions outside of the cell. However, an inverse relationship between Z(C) and redox potential occurs between the endoplasmic reticulum and cytoplasm. This trend provides support for the hypothesis that protein transport and turnover are ultimately coupled to the maintenance of different glutathione redox potentials in subcellular compartments. There are broad changes in Z(C) in whole-genome protein compositions in microbes from different environments, and in Rubisco homologues, lower Z(C) tends to occur in organisms with higher optimal growth temperature. Energetic costs calculated from thermodynamic models are consistent with the notion that thermophilic organisms exhibit molecular adaptation to not only high temperature but also the reducing nature of many hydrothermal fluids. Further characterization of the material requirements of protein metabolism in terms of the chemical conditions of cells and environments may help to reveal other linkages among biochemical processes with implications for changes on evolutionary time scales.

Enhanced Indirect Photochemical Transformation of Histidine and Histamine through Association with Chromophoric Dissolved Organic Matter

Photochemical transformations greatly affect the stability and fate of amino acids (AAs) in sunlit aquatic ecosystems. Whereas the direct phototransformation of dissolved AAs is well investigated, their indirect photolysis in the presence of chromophoric dissolved organic matter (CDOM) is poorly understood. In aquatic systems, CDOM may act both as sorbent for AAs and as photosensitizer, creating microenvironments with high concentrations of photochemically produced reactive intermediates, such as singlet oxygen (1O2). This study provides a systematic investigation of the indirect photochemical transformation of histidine (His) and histamine by 1O2 in solutions containing CDOM as a function of solution pH. Both His and histamine showed pH-dependent enhanced phototransformation in the CDOM systems as compared to systems in which model, low-molecular-weight 1O2 sensitizers were used. Enhanced reactivity resulted from sorption of His and histamine to CDOM and thus exposure to elevated 1O2 concentrations in the CDOM microenvironment. The extent of reactivity enhancement depended on solution pH via its effects on the protonation state of His, histamine, and CDOM. Sorption-enhanced reactivity was independently supported by depressed rate enhancements in the presence of a cosorbate that competitively displaced His and histamine from CDOM. Incorporating sorption and photochemical transformation processes into a reaction rate prediction model improved the description of the abiotic photochemical transformation rates of His in the presence of CDOM.

The production of methane, hydrogen, and organic compounds in ultramafic-hosted hydrothermal vents of the Mid-Atlantic Ridge

Both hydrogen and methane are consistently discharged in large quantities in hydrothermal fluids issued from ultramafic-hosted hydrothermal fields discovered along the Mid-Atlantic Ridge. Considering the vast number of these fields discovered or inferred, hydrothermal fluxes represent a significant input of H2 and CH4 to the ocean. Although there are lines of evidence of their abiogenic formation from stable C and H isotope results, laboratory experiments, and thermodynamic data, neither their origin nor the reaction pathways generating these gases have been fully constrained yet. Organic compounds detected in the fluids may also be derived from abiotic reactions. Although thermodynamics are favorable and extensive experimental work has been done on Fischer-Tropsch-type reactions, for instance, nothing is clear yet about their origin and formation mechanism from actual data. Since chemolithotrophic microbial communities commonly colonize hydrothermal vents, biogenic and thermogenic processes are likely to contribute to the production of H2, CH4, and other organic compounds. There seems to be a consensus toward a mixed origin (both sources and processes) that is consistent with the ambiguous nature of the isotopic data. But the question that remains is, to what proportions? More systematic experiments as well as integrated geochemical approaches are needed to disentangle hydrothermal geochemistry. This understanding is of prime importance considering the implications of hydrothermal H2, CH4, and organic compounds for the ocean global budget, global cycles, and the origin of life.

An Interpretation of the Ancestral Codon from Miller's Amino Acids and Nucleotide Correlations in Modern Coding Sequences

Purine bias, which is usually referred to as an “ancestral codon”, is known to result in short-range correlations between nucleotides in coding sequences, and it is common in all species. We demonstrate that RWY is a more appropriate pattern than the classical RNY, and purine bias (Rrr) is the product of a network of nucleotide compensations induced by functional constraints on the physicochemical properties of proteins. Through deductions from universal correlation properties, we also demonstrate that amino acids from Miller’s spark discharge experiment are compatible with functional primeval proteins at the dawn of living cell radiation on earth. These amino acids match the hydropathy and secondary structures of modern proteins.

Autocatalytic sets in E. coli metabolism

Background

A central unsolved problem in early evolution concerns self-organization towards higher complexity in chemical reaction networks. In theory, autocatalytic sets have useful properties to help model such transitions. Autocatalytic sets are chemical reaction systems in which molecules belonging to the set catalyze the synthesis of other members of the set. Given an external supply of starting molecules – the food set – and the conditions that (i) all reactions are catalyzed by at least one molecule, and (ii) each molecule can be constructed from the food set by a sequence of reactions, the system becomes a reflexively autocatalytic food-generated network (RAF set). Autocatalytic networks and RAFs have been studied extensively as mathematical models for understanding the properties and parameters that influence self-organizational tendencies. However, despite their appeal, the relevance of RAFs for real biochemical networks that exist in nature has, so far, remained virtually unexplored.

Results

Here we investigate the best-studied metabolic network, that of Escherichia coli, for the existence of RAFs. We find that the largest RAF encompasses almost the entire E. coli cytosolic reaction network. We systematically study its structure by considering the impact of removing catalysts or reactions. We show that, without biological knowledge, finding the minimum food set that maintains a given RAF is NP-complete. We apply a randomized algorithm to find (approximately) smallest subsets of the food set that suffice to sustain the original RAF.

Conclusions

The existence of RAF sets within a microbial metabolic network indicates that RAFs capture properties germane to biological organization at the level of single cells. Moreover, the interdependency between the different metabolic modules, especially concerning cofactor biosynthesis, points to the important role of spontaneous (non-enzymatic) reactions in the context of early evolution.

Graphical Abstract

Extraordinarily adaptive properties of the genetically encoded amino acids

Using novel advances in computational chemistry, we demonstrate that the set of 20 genetically encoded amino acids, used nearly universally to construct all coded terrestrial proteins, has been highly influenced by natural selection. We defined an adaptive set of amino acids as one whose members thoroughly cover relevant physico-chemical properties, or “chemistry space.” Using this metric, we compared the encoded amino acid alphabet to random sets of amino acids. These random sets were drawn from a computationally generated compound library containing 1913 alternative amino acids that lie within the molecular weight range of the encoded amino acids. Sets that cover chemistry space better than the genetically encoded alphabet are extremely rare and energetically costly. Further analysis of more adaptive sets reveals common features and anomalies, and we explore their implications for synthetic biology. We present these computations as evidence that the set of 20 amino acids found within the standard genetic code is the result of considerable natural selection. The amino acids used for constructing coded proteins may represent a largely global optimum, such that any aqueous biochemistry would use a very similar set.

Novel type II and monomeric NAD+ specific isocitrate dehydrogenases: phylogenetic affinity, enzymatic characterization, and evolutionary implication

NAD(+) use is an ancestral trait of isocitrate dehydrogenase (IDH), and the NADP(+) phenotype arose through evolution as an ancient adaptation event. However, no NAD(+)-specific IDHs have been found among type II IDHs and monomeric IDHs. In this study, novel type II homodimeric NAD-IDHs from Ostreococcus lucimarinus CCE9901 IDH (OlIDH) and Micromonas sp. RCC299 (MiIDH), and novel monomeric NAD-IDHs from Campylobacter sp. FOBRC14 IDH (CaIDH) and Campylobacter curvus (CcIDH) were reported for the first time. The homodimeric OlIDH and monomeric CaIDH were determined by size exclusion chromatography and MALDI-TOF/TOF mass spectrometry. All the four IDHs were demonstrated to be NAD(+)-specific, since OlIDH, MiIDH, CaIDH and CcIDH displayed 99-fold, 224-fold, 61-fold and 37-fold preferences for NAD(+) over NADP(+), respectively. The putative coenzyme discriminating amino acids (Asp326/Met327 in OlIDH, Leu584/Asp595 in CaIDH) were evaluated, and the coenzyme specificities of the two mutants, OlIDH R(326)H(327) and CaIDH H(584)R(595), were completely reversed from NAD(+) to NADP(+). The detailed biochemical properties, including optimal reaction pH and temperature, thermostability, and metal ion effects, of OlIDH and CaIDH were further investigated. The evolutionary connections among OlIDH, CaIDH, and all the other forms of IDHs were described and discussed thoroughly.

Miller experiments in atomistic computer simulations

The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.

The drive to life on wet and icy worlds

This paper presents a reformulation of the submarine alkaline hydrothermal theory for the emergence of life in response to recent experimental findings. The theory views life, like other self-organizing systems in the Universe, as an inevitable outcome of particular disequilibria. In this case, the disequilibria were two: (1) in redox potential, between hydrogen plus methane with the circuit-completing electron acceptors such as nitrite, nitrate, ferric iron, and carbon dioxide, and (2) in pH gradient between an acidulous external ocean and an alkaline hydrothermal fluid. Both CO2 and CH4 were equally the ultimate sources of organic carbon, and the metal sulfides and oxyhydroxides acted as protoenzymatic catalysts. The realization, now 50 years old, that membrane-spanning gradients, rather than organic intermediates, play a vital role in life's operations calls into question the idea of "prebiotic chemistry." It informs our own suggestion that experimentation should look to the kind of nanoengines that must have been the precursors to molecular motors-such as pyrophosphate synthetase and the like driven by these gradients-that make life work. It is these putative free energy or disequilibria converters, presumably constructed from minerals comprising the earliest inorganic membranes, that, as obstacles to vectorial ionic flows, present themselves as the candidates for future experiments. Key Words: Methanotrophy-Origin of life. Astrobiology 14, 308-343. The fixation of inorganic carbon into organic material (autotrophy) is a prerequisite for life and sets the starting point of biological evolution. (Fuchs, 2011 ) Further significant progress with the tightly membrane-bound H(+)-PPase family should lead to an increased insight into basic requirements for the biological transport of protons through membranes and its coupling to phosphorylation. (Baltscheffsky et al., 1999 ).

Thermodynamic prediction of glycine polymerization as a function of temperature and pH consistent with experimentally obtained results

Prediction of the thermodynamic behaviors of biomolecules at high temperature and pressure is fundamental to understanding the role of hydrothermal systems in the origin and evolution of life on the primitive Earth. However, available thermodynamic dataset for amino acids, essential components for life, cannot represent experimentally observed polymerization behaviors of amino acids accurately under hydrothermal conditions. This report presents the thermodynamic data and the revised HKF parameters for the simplest amino acid "Gly" and its polymers (GlyGly, GlyGlyGly and DKP) based on experimental thermodynamic data from the literature. Values for the ionization states of Gly (Gly(+) and Gly(-)) and Gly peptides (GlyGly(+), GlyGly(-), GlyGlyGly(+), and GlyGlyGly(-)) were also retrieved from reported experimental data by combining group additivity algorithms. The obtained dataset enables prediction of the polymerization behavior of Gly as a function of temperature and pH, consistent with experimentally obtained results in the literature. The revised thermodynamic data for zwitterionic Gly, GlyGly, and DKP were also used to estimate the energetics of amino acid polymerization into proteins. Results show that the Gibbs energy necessary to synthesize a mole of peptide bond is more than 10 kJ mol(-1) less than previously estimated over widely various temperatures (e.g., 28.3 kJ mol(-1) → 17.1 kJ mol(-1) at 25 °C and 1 bar). Protein synthesis under abiotic conditions might therefore be more feasible than earlier studies have shown.