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

Insights into Early Steps of Decanoic Acid Self-Assemblies under Prebiotic Temperatures Using Molecular Dynamics Simulations

The origin of life possibly required processes in confined systems that facilitated simple chemical reactions and other more complex reactions impossible to achieve under the condition of infinite dilution. In this context, the self-assembly of micelles or vesicles derived from prebiotic amphiphilic molecules is a cornerstone in the chemical evolution pathway. A prime example of these building blocks is decanoic acid, a short-chain fatty acid capable of self-assembling under ambient conditions. This study explored a simplified system made of decanoic acids under temperatures ranging from 0 °C to 110 °C to replicate prebiotic conditions. The study revealed the first point of aggregation of decanoic acid into vesicles and examined the insertion of a prebiotic-like peptide in a primitive bilayer. The information gathered from this research provides critical insights into molecule interactions with primitive membranes, allowing us to understand the first nanometric compartments needed to trigger further reactions that were essential for the origin of life.

Peptide-Based Supramolecular Systems Chemistry

Peptide-based supramolecular systems chemistry seeks to mimic the ability of life forms to use conserved sets of building blocks and chemical reactions to achieve a bewildering array of functions. Building on the design principles for short peptide-based nanomaterials with properties, such as self-assembly, recognition, catalysis, and actuation, are increasingly available. Peptide-based supramolecular systems chemistry is starting to address the far greater challenge of systems-level design to access complex functions that emerge when multiple reactions and interactions are coordinated and integrated. We discuss key features relevant to systems-level design, including regulating supramolecular order and disorder, development of active and adaptive systems by considering kinetic and thermodynamic design aspects and combinatorial dynamic covalent and noncovalent interactions. Finally, we discuss how structural and dynamic design concepts, including preorganization and induced fit, are critical to the ability to develop adaptive materials with adaptive and tunable photonic, electronic, and catalytic properties. Finally, we highlight examples where multiple features are combined, resulting in chemical systems and materials that display adaptive properties that cannot be achieved without this level of integration.

An Isosymmetric High-Pressure Phase Transition in α-Glycylglycine: A Combined Experimental and Theoretical Study

We investigated the effects of hydrostatic pressure on α-glycylglycine (α-digly) using a combined experimental and theoretical approach. The results of powder X-ray diffraction show a change in compressibility of the axes above 6.7 GPa, but also indicate that the structure remains in the same monoclinic space group, suggesting an isosymmetric phase transition. A noticeable change in the Raman spectra between 6 and 7.5 GPa further supports the observed phase transition. First-principles-based calculations combined with the crystal structure prediction code USPEX predict a number of possible polymorphs at high pressure. An orthorhombic structure with a bent peptide backbone is the lowest enthalpy polymorph above 6.4 GPa; however, it is not consistent with experimental observations. A second monoclinic structure isosymmetric to α-digly, α'-digly, is predicted to become more stable above 11.4 GPa. The partial atomic charges in α'-digly differ from α-digly, and the molecule is bent, possibly indicating different reactivity of α'-digly. The similarity in the lattice parameters predicted from calculations and the axial changes observed experimentally support that the α'-digly phase is likely observed at high pressure. A possible explanation for the isosymmetric phase transition is discussed in terms of relaxing strained hydrogen bonding interactions. Such combined experimental and modeling efforts provide atomic-level insight into how pressure-driven conformational changes alter hydrogen-bonding networks in complicated molecular crystals.

Thermodynamic Impact of Mineral Surfaces on Amino Acid Polymerization: Aspartate Dimerization on Goethite

This article presents a thermodynamic predictive scheme for amino acid polymerization in the presence of minerals as a function of various environmental parameters (pH, ionic strength, amino acid concentration, and the solid/water ratio) using l-aspartate (Asp) and goethite as a model combination. This prediction is enabled by the combination of the surface adsorption constants of amino acid and its polymer, determined from the extended triple layer model characterization of the corresponding experimental results, with the thermodynamic data of these organic compounds in water reported in the literature. Calculations for the Asp-goethite system showed that the goethite surface drastically shifts the Asp monomer-dipeptide equilibrium toward the dipeptide side; when the dimerization of 0.1 mM Asp was considered in the presence of 10 m² L^-1 of goethite, an Asp dipeptide concentration around 10⁵ times larger was computed to be thermodynamically attainable compared with that in the absence of goethite at acidic pH (4-5) and low ionic strength (0.1 mM NaCl). Under this condition, the dipeptide-to-monomer molecular ratio in the adsorbed state reached 20%. In contrast, no significant enhancement by goethite was predicted at alkaline pH (>8), where the electrostatic interactions of the goethite surface with Asp and Asp dipeptide are weak. Thus, mineral surfaces should have had a significant impact on the thermodynamics of prebiotic peptide bond formation on the early Earth, although the influences likely depended largely on the environmental conditions. Future experimental studies for various amino acid-mineral interactions using our proposed methodology will provide a quantitative constraint on favorable geochemical settings for the chemical evolution on Earth. This approach can also offer important clues for future exploration of extraterrestrial life.

Geochemistry and the Origin of Life: From Extraterrestrial Processes, Chemical Evolution on Earth, Fossilized Life's Records, to Natures of the Extant Life

In 2001, the first author (S.N.) led the publication of a book entitled "Geochemistry and the origin of life" in collaboration with Dr. Andre Brack aiming to figure out geo- and astro-chemical processes essential for the emergence of life. Since then, a great number of research progress has been achieved in the relevant topics from our group and others, ranging from the extraterrestrial inputs of life's building blocks, the chemical evolution on Earth with the aid of mineral catalysts, to the fossilized records of ancient microorganisms. Here, in addition to summarizing these findings for the origin and early evolution of life, we propose a new hypothesis for the generation and co-evolution of photosynthesis with the redox and photochemical conditions on the Earth's surface. Besides these bottom-up approaches, we introduce an experimental study on the role of water molecules in the life's function, focusing on the transition from live, dormant, and dead states through dehydration/hydration. Further spectroscopic studies on the hydrogen bonding behaviors of water molecules in living cells will provide important clues to solve the complex nature of life.

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.

Generating Converged Accurate Free Energy Surfaces for Chemical Reactions with a Force-Matched Semiempirical Model

We demonstrate the capability of creating robust density functional tight binding (DFTB) models for chemical reactivity in prebiotic mixtures through force matching to short time scale quantum free energy estimates. Molecular dynamics using density functional theory (DFT) is a highly accurate approach to generate free energy surfaces for chemical reactions, but the extreme computational cost often limits the time scales and range of thermodynamic states that can feasibly be studied. In contrast, DFTB is a semiempirical quantum method that affords up to a thousandfold reduction in cost and can recover DFT-level accuracy. Here, we show that a force-matched DFTB model for aqueous glycine condensation reactions yields free energy surfaces that are consistent with experimental observations of reaction energetics. Convergence analysis reveals that multiple nanoseconds of combined trajectory are needed to reach a steady-fluctuating free energy estimate for glycine condensation. Predictive accuracy of force-matched DFTB is demonstrated by direct comparison to DFT, with the two approaches yielding surfaces with large regions that differ by only a few kcal mol^-1.

Prebiotic Synthesis of Glycine from Ethanolamine in Simulated Archean Alkaline Hydrothermal Vents

Submarine hydrothermal vents are generally considered as the likely habitats for the origin and evolution of early life on Earth. In recent years, a novel hydrothermal system in Archean subseafloor has been proposed. In this model, highly alkaline and high temperature hydrothermal fluids were generated in basalt-hosted hydrothermal vents, where H2 and CO2 could be abundantly provided. These extreme conditions could have played an irreplaceable role in the early evolution of life. Nevertheless, sufficient information has not yet been obtained for the abiotic synthesis of amino acids, which are indispensable components of life, at high temperature and alkaline condition. This study aims to propose a new method for the synthesis of glycine in simulated Archean submarine alkaline vent systems. We investigated the formation of glycine from ethanolamine under conditions of high temperature (80-160 °C) and highly alkaline solutions (pH = 9.70). Experiments were performed in an anaerobic environment under mild pressure (0.1-8.0 MPa) at the same time. The results suggested that the formation of glycine from ethanolamine occurred rapidly and efficiently in the presence of metal powders, and was favored by high temperatures and high pressures. The experiment provides a new pathway for prebiotic glycine formation and points out the phenomenal influence of high-temperature alkaline hydrothermal vents in origin of life in the early ocean.

Origin of the Reductive Tricarboxylic Acid (rTCA) Cycle-Type CO2 Fixation: A Perspective

The reductive tricarboxylic acid (rTCA) cycle is among the most plausible candidates for the first autotrophic metabolism in the earliest life. Extant enzymes fixing CO₂ in this cycle contain cofactors at the catalytic centers, but it is unlikely that the protein/cofactor system emerged at once in a prebiotic process. Here, we discuss the feasibility of non-enzymatic cofactor-assisted drive of the rTCA reactions in the primitive Earth environments, particularly focusing on the acetyl-CoA conversion to pyruvate. Based on the energetic and mechanistic aspects of this reaction, we propose that the deep-sea hydrothermal vent environments with active electricity generation in the presence of various sulfide catalysts are a promising setting for it to progress. Our view supports the theory of an autotrophic origin of life from primordial carbon assimilation within a sulfide-rich hydrothermal vent.

Glycine Polymerization on Oxide Minerals

It has long been suggested that mineral surfaces played an important role in peptide bond formation on the primitive Earth. However, it remains unclear which mineral species was key to the prebiotic processes. This is because great discrepancies exist among the reported catalytic efficiencies of minerals for amino acid polymerizations, owing to mutually different experimental conditions. This study examined polymerization of glycine (Gly) on nine oxide minerals (amorphous silica, quartz, α-alumina and γ-alumina, anatase, rutile, hematite, magnetite, and forsterite) using identical preparation, heating, and analytical procedures. Results showed that a rutile surface is the most effective site for Gly polymerization in terms of both amounts and lengths of Gly polymers synthesized. The catalytic efficiency decreased as rutile > anatase > γ-alumina > forsterite > α- alumina > magnetite > hematite > quartz > amorphous silica. Based on reported molecular-level information for adsorption of Gly on these minerals, polymerization activation was inferred to have arisen from deprotonation of the NH3+ group of adsorbed Gly to the nucleophilic NH2 group, and from withdrawal of electron density from the carboxyl carbon to the surface metal ions. The orientation of adsorbed Gly on minerals is also a factor influencing the Gly reactivity. The examination of Gly-mineral interactions under identical experimental conditions has enabled the direct comparison of various minerals' catalytic efficiencies and has made discussion of polymerization mechanisms and their relative influences possible Further systematic investigations using the approach reported herein (which are expected to be fruitful) combined with future microscopic surface analyses will elucidate the role of minerals in the process of abiotic peptide bond formation.

Dissolved Divalent Metal and pH Effects on Amino Acid Polymerization: A Thermodynamic Evaluation

Polymerization of amino acids is a fundamentally important step for the chemical evolution of life. Nevertheless, its response to changing environmental conditions has not yet been well understood because of the lack of reliable quantitative information. For thermodynamics, detailed prediction over diverse combinations of temperature and pH has been made only for a few amino acid-peptide systems. This study used recently reported thermodynamic dataset for the polymerization of the simplest amino acid "glycine (Gly)" to its short peptides (di-glycine and tri-glycine) to examine chemical and structural characteristics of amino acids and peptides that control the temperature and pH dependence of polymerization. Results showed that the dependency is strongly controlled by the intramolecular distance between the amino and carboxyl groups in an amino acid structure, although the side-chain group role is minor. The polymerization behavior of Gly reported earlier in the literature is therefore expected to be a typical feature for those of α-amino acids. Equilibrium calculations were conducted to examine effects of dissolved metals as a function of pH on the monomer-polymer equilibria of Gly. Results showed that metals shift the equilibria toward the monomer side, particularly at neutral and alkaline pH. Metals that form weak interaction with Gly (e.g., Mg2+) have no noticeable influence on the polymerization, although strong interaction engenders significant decrease of the equilibrium concentrations of Gly peptides. Considering chemical and structural characteristics of Gly and Gly peptides that control their interactions with metals, it can be expected that similar responses to the addition of metals are applicable in the polymerization of neutral α-amino acids. Neutral and alkaline aqueous environments with dissolved metals having high affinity with neutral α-amino acids (e.g., Cu2+) are therefore not beneficial places for peptide bond formation on the primitive Earth.

Predicting Thermodynamic Behaviors of Non-Protein Amino Acids as a Function of Temperature and pH

Why does life use α-amino acids exclusively as building blocks of proteins? To address that fundamental question from an energetic perspective, this study estimated the standard molal thermodynamic data for three non-α-amino acids (β-alanine, γ-aminobutyric acid, and ε-aminocaproic acid) and α-amino-n-butyric acid in their zwitterionic, negative, and positive ionization states based on the corresponding experimental measurements reported in the literature. Temperature dependences of their heat capacities were described based on the revised Helgeson-Kirkham-Flowers (HKF) equations of state. The obtained dataset was then used to calculate the standard molal Gibbs energies (∆G (o)) of the non-α-amino acids as a function of temperature and pH. Comparison of their ∆G (o) values with those of α-amino acids having the same molecular formula showed that the non-α-amino acids have similar ∆G (o) values to the corresponding α-amino acids in physiologically relevant conditions (neutral pH, <100 °C). In acidic and alkaline pH, the non-α-amino acids are thermodynamically more stable than the corresponding α-ones over a broad temperature range. These results suggest that the energetic cost of synthesis is not an important selection pressure to incorporate α-amino acids into biological systems.

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.

Polymerization of Building Blocks of Life on Europa and Other Icy Moons

The outer Solar System may provide a potential habitat for extraterrestrial life. Remote sensing data from the Galileo spacecraft suggest that the jovian icy moons--Europa, Ganymede, and possibly Callisto--may harbor liquid water oceans underneath their icy crusts. Although compositional information required for the discussion of habitability is limited because of significantly restricted observation data, organic molecules are ubiquitous in the Universe. Recently, in situ spacecraft measurements and experiments suggest that amino acids can be formed abiotically on interstellar ices and comets. These amino acids could be continuously delivered by meteorite or comet impacts to icy moons. Here, we show that polymerization of organic monomers, in particular amino acids and nucleotides, could proceed spontaneously in the cold environment of icy moons, in particular the jovian icy moon Europa as a typical example, based on thermodynamic calculations, though kinetics of formation are not addressed. Observed surface temperature on Europa is 120 and 80 K in the equatorial region and polar region, respectively. At such low temperatures, Gibbs energies of polymerization become negative, and the estimated thermal structure of the icy crust should contain a shallow region (i.e., at a depth of only a few kilometers) favorable for polymerization. Investigation of the possibility of organic monomer polymerization on icy moons could provide good constraints on the origin and early evolution of extraterrestrial life.