Calculation of the standard molal thermodynamic properties of aqueous biomolecules at elevated temperatures and pressures II. Unfolded proteins

The minimum energy required to build a cell

Understanding the energy requirements for cell synthesis accurately and comprehensively has been a longstanding challenge. We introduce a computational model that estimates the minimum energy necessary to build any cell from its constituent parts. This method combines omics and internal cell compositions from various sources to calculate the Gibbs Free Energy of biosynthesis independently of specific metabolic pathways. Our public tool, Synercell, can be used with other models for minumum species-specific energy estimations in any well-sequenced species. The energy for synthesising the genome, transcriptome, proteome, and lipid bilayer of four cell types: Escherichia coli, Saccharomyces cerevisiae, an average mammalian cell and JCVI-syn3A were estimated. Their modelled minimum synthesis energies at 298 K were 9.54 × 10 - 11 J/cell, 4.99 × 10 - 9 J/cell, 3.71 × 10 - 7 J/cell and 3.69 × 10 - 12 respectively. Gram-for-gram synthesis of lipid bilayers requires the most energy, followed by the proteome, genome, and transcriptome. The average per gram cost of biomass synthesis is in the 300s of J/g for all four cells. Implications for the generalisability of cell construction and applications to biogeosciences, cellular biology, biotechnology, and astrobiology are discussed.

Reaction pathways and kinetics for tetra-alanine in hot, compressed liquid water

pH alters the prevalence of the dissociation states of peptides, through acid–base equilibrium, and influences reaction selectivity under hydrothermal conditions.

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.

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.

The energetics of organic synthesis inside and outside the cell

Thermodynamic modelling of organic synthesis has largely been focused on deep-sea hydrothermal systems. When seawater mixes with hydrothermal fluids, redox gradients are established that serve as potential energy sources for the formation of organic compounds and biomolecules from inorganic starting materials. This energetic drive, which varies substantially depending on the type of host rock, is present and available both for abiotic (outside the cell) and biotic (inside the cell) processes. Here, we review and interpret a library of theoretical studies that target organic synthesis energetics. The biogeochemical scenarios evaluated include those in present-day hydrothermal systems and in putative early Earth environments. It is consistently and repeatedly shown in these studies that the formation of relatively simple organic compounds and biomolecules can be energy-yielding (exergonic) at conditions that occur in hydrothermal systems. Expanding on our ability to calculate biomass synthesis energetics, we also present here a new approach for estimating the energetics of polymerization reactions, specifically those associated with polypeptide formation from the requisite amino acids.

Calculation of the relative chemical stabilities of proteins as a function of temperature and redox chemistry in a hot spring

Uncovering the chemical and physical links between natural environments and microbial communities is becoming increasingly amenable owing to geochemical observations and metagenomic sequencing. At the hot spring known as Bison Pool in Yellowstone National Park, the cooling of the water in the outflow channel is associated with an increase in oxidation potential estimated from multiple field-based measurements. Representative groups of proteins whose sequences were derived from metagenomic data also exhibit an increase in average oxidation state of carbon in the protein molecules with distance from the hot-spring source. The energetic requirements of reactions to form selected proteins used in the model were computed using amino-acid group additivity for the standard molal thermodynamic properties of the proteins, and the relative chemical stabilities of the proteins were investigated by varying temperature, pH and oxidation state, expressed as activity of dissolved hydrogen. The relative stabilities of the proteins were found to track the locations of the sampling sites when the calculations included a function for hydrogen activity that increases with temperature and is higher, or more reducing, than values consistent with measurements of dissolved oxygen, sulfide and oxidation-reduction potential in the field. These findings imply that spatial patterns in the amino acid compositions of proteins can be linked, through energetics of overall chemical reactions representing the formation of the proteins, to the environmental conditions at this hot spring, even if microbial cells maintain considerably different internal conditions. Further applications of the thermodynamic calculations are possible for other natural microbial ecosystems.

Peptide synthesis in early Earth hydrothermal systems

We report here results from experiments and thermodynamic calculations that demonstrate a rapid, temperature-enhanced synthesis of oligopeptides from the condensation of aqueous glycine. Experiments were conducted in custom-made hydrothermal reactors, and organic compounds were characterized with ultraviolet-visible procedures. A comparison of peptide yields at 260 degrees C with those obtained at more moderate temperatures (160 degrees C) gives evidence of a significant (13 kJ . mol(-1)) exergonic shift. In contrast to previous hydrothermal studies, we demonstrate that peptide synthesis is favored in hydrothermal fluids and that rates of peptide hydrolysis are controlled by the stability of the parent amino acid, with a critical dependence on reactor surface composition. From our study, we predict that rapid recycling of product peptides from cool into near-supercritical fluids in mid-ocean ridge hydrothermal systems will enhance peptide chain elongation. It is anticipated that the abundant hydrothermal systems on early Earth could have provided a substantial source of biomolecules required for the origin of life.

Calculation of the relative metastabilities of proteins using the CHNOSZ software package

BACKGROUND

Proteins of various compositions are required by organisms inhabiting different environments. The energetic demands for protein formation are a function of the compositions of proteins as well as geochemical variables including temperature, pressure, oxygen fugacity and pH. The purpose of this study was to explore the dependence of metastable equilibrium states of protein systems on changes in the geochemical variables.

RESULTS

A software package called CHNOSZ implementing the revised Helgeson-Kirkham-Flowers (HKF) equations of state and group additivity for ionized unfolded aqueous proteins was developed. The program can be used to calculate standard molal Gibbs energies and other thermodynamic properties of reactions and to make chemical speciation and predominance diagrams that represent the metastable equilibrium distributions of proteins. The approach takes account of the chemical affinities of reactions in open systems characterized by the chemical potentials of basis species. The thermodynamic database included with the package permits application of the software to mineral and other inorganic systems as well as systems of proteins or other biomolecules.

CONCLUSION

Metastable equilibrium activity diagrams were generated for model cell-surface proteins from archaea and bacteria adapted to growth in environments that differ in temperature and chemical conditions. The predicted metastable equilibrium distributions of the proteins can be compared with the optimal growth temperatures of the organisms and with geochemical variables. The results suggest that a thermodynamic assessment of protein metastability may be useful for integrating bio- and geochemical observations.

Group additivity schemes for the calculation of the partial molar heat capacities and volumes of unfolded proteins in aqueous solution

A critical review is given of the present state of group additivity schemes for the calculation of partial molar volumes and heat capacities of unfolded proteins. The comparison between the experimental values and the predictions based on the different models shows clearly that only the peptide-based additivity scheme represents properly both the absolute values and the temperature dependence of these thermodynamic quantities.

Group additivity calculations of the thermodynamic properties of unfolded proteins in aqueous solution: a critical comparison of peptide-based and HKF models

A recent paper in this journal [Amend and Helgeson, Biophys. Chem. 84 (2000) 105] presented a new group additivity model to calculate various thermodynamic properties of unfolded proteins in aqueous solution. The parameters given for the revised Helgeson-Kirkham-Flowers (HKF) equations of state for all the constituent groups of unfolded proteins can be used, in principle, to calculate the partial molar heat capacity, C(o)p.2, and volume, V2(0), at infinite dilution of any polypeptide. Calculations of the values of C(o)p.2 and V2(0) for several polypeptides have been carried out to test the predictive utility of the HKF group additivity model. The results obtained are in very poor agreement with experimental data, and also with results calculated using a peptide-based group additivity model. A critical assessment of these two additivity models is presented.