The unique features of glycolytic pathways in Archaea

Chimeric origins and dynamic evolution of central carbon metabolism in eukaryotes

The origin of eukaryotes was a key event in the history of life. Current leading hypotheses propose that a symbiosis between an asgardarchaeal host cell and an alphaproteobacterial endosymbiont represented a crucial step in eukaryotic origin and that metabolic cross-feeding between the partners provided the basis for their subsequent evolutionary integration. A major unanswered question is whether the metabolism of modern eukaryotes bears any vestige of this ancestral syntrophy. Here we systematically analyse the evolutionary origins of the eukaryotic gene repertoires mediating central carbon metabolism. Our phylogenetic and sequence analyses reveal that this gene repertoire is chimeric, with ancestral contributions from Asgardarchaeota and Alphaproteobacteria operating predominantly in glycolysis and the tricarboxylic acid cycle, respectively. Our analyses also reveal the extent to which this ancestral metabolic interplay has been remodelled via gene loss, transfer and subcellular retargeting in the >2 billion years since the origin of eukaryotic cells, and we identify genetic contributions from other prokaryotic sources in addition to the asgardarchaeal host and alphaproteobacterial endosymbiont. Our work demonstrates that, in contrast to previous assumptions, modern eukaryotic metabolism preserves information about the nature of the original asgardarchaeal-alphaproteobacterial interactions and supports syntrophy scenarios for the origin of the eukaryotic cell.

Sugar-sensing swodkoreceptors and swodkocrine signaling

Sugars are one of the major metabolites and are essential for nucleic acid synthesis and energy production. In addition, sugars can act as signaling molecules. To study sugar signaling at the systemic level, there is an urgent need to systematically identify sugar-sensing proteins and nucleic acids. I propose the terms "swodkoreceptor" and "swodkocrine signaling," derived from the Polish word "słodki" meaning "sweet," to comprise all sugar-sensing proteins and signaling events, respectively, regardless of their cellular location and signaling domains. This proposal is intended to facilitate the inclusion of proteins such as the Escherichia coli LacI repressor as an allolactose receptor, human glucokinase regulatory protein (GCKR) as a fructose receptor, and other sugar-binding based allosterically regulated enzymes and transcription factors as sugar-sensing receptors. In addition, enzyme-interacting proteins whose interaction state is regulated by sugar binding have also been proposed as sugar receptors. The systemic study of protein- and nucleic-acid-based swodkoreceptors may help to identify organelle-specific swodkoreceptors and to also address receptor duality. The study of intra- and inter-organism swodkocrine signaling and its crosstalk with gasocrine signaling may help to understand the etiology of diseases due to dysregulation in sugar homeostasis and signaling.

Formate addition enhanced hydrogen production by Thermococcus paralvinellae when grown on brewery wastewater

The hyperthermophilic archaeon Thermococcus paralvinellae produces H2 when grown on carbohydrates or protein with increased H2 production when cultures are grown on formate. This study examined the use of brewery wastewater as a feedstock for H2 production, the addition of formate to enhance H2 production, and the activities of hydrogenases and formate hydrogenlyase under varying growth conditions as markers of performance. T. paralvinellae was grown at 80°C on maltose only (a model brewery wastewater), formate only, and maltose plus formate media as well as brewery wastewater with and without the addition of formate. Growth rates were higher on formate only medium than on maltose only and brewery wastewater only media. H2 yield per cell was higher in all media containing formate relative to those without formate. Hydrogenase and formate hydrogenlyase specific activities were not affected by the presence of formate and were largely consistent across all growth conditions. Growth rates were consistent in media containing 0.05 to 2.5% (wt/vol) maltose only, but total H2 production doubled from medium containing 0.05% maltose to 0.5% maltose and remained unchanged at higher maltose concentrations. Cells grown in a 2 L N2 flushed batch bioreactor at 80°C on brewery wastewater with and without formate showed no difference in growth rates but the amount of H2 in the headspace was six times higher when formate was present. However, the amount of H2 produced by cells grown on brewery wastewater plus formate peaked in mid-logarithmic growth phase and then decreased to amounts produced by cells without formate addition by late logarithmic growth phase. When the bioreactor was run as a chemostat, the addition of formate to brewery wastewater led to a 12-fold increase in the amount of H2 present in the headspace that was sustained over time relative to growth without formate addition. Therefore, T. paralvinellae grows on brewery wastewater as its sole source of organic carbon and produces biohydrogen at a steady rate in a pilot-scale bioreactor, and H2 production is enhanced by formate addition.

Expanding the Enzymatic Toolbox for Carboligation: Increasing the Diversity of the 'Split' Transketolase Sequence Space

Transketolases (TKs) are thiamine diphosphate (ThDP)-dependent enzymes that catalyze the transfer of two-carbon units in a stereoselective manner, making them valuable biocatalysts for sustainable processes. Most known TKs are about 650 amino acids long; however, a second type found in Archaea and many Bacteria consists of two proteins, each of about 300 amino acids. Exploring the unique features and differences of split TKs may help in assessing their potential use in biocatalysis and for uncovering new reactivities. Additionally, it could provide valuable information on how their structure relates to their function, especially compared to full-length TKs. In this study, we significantly expanded the known repertoire of split TKs approximately 14-fold to the best of our knowledge, by identifying and providing accessions of nearly 500 putative split-TK subunit pairs. Moreover, we doubled the number of experimentally produced and tested split TKs by cloning, purifying, and testing ten candidates retrieved from genomes and in-house metagenomes. Interestingly, pQR2809 and pQR2812, derived from hyperthermophilic organisms, showed enhanced thermostability compared to other TK examples in the literature, maintaining partial activity after heating at 90 °C or 100 °C for 1 hour, respectively.

Comparative genomics of the highly halophilic Haloferacaceae

The Haloferacaceae are a family of extremely halophilic archaea with many species producing enzymes and products beneficial for industrial biotechnology. They are, however, relatively under-characterised with regards to genetics and gene products. This study aims to use existing sequence data to highlight genetic diversity, create pangenomes for three genera, and provide secondary metabolite and pathway analysis. This will establish current knowledge and identify key gaps in research. We show that the Haloferacaceae have significant genetic diversity between genera, with numerous gene gain and loss events in key genera. It also found that the model genus, Haloferax, has relatively low identified secondary metabolites compared to other genera within the family. Additionally, this study has identified potential biotechnology targets for heterologous expression in model organisms.

Navigating the archaeal frontier: insights and projections from bioinformatic pipelines

Archaea continues to be one of the least investigated domains of life, and in recent years, the advent of metagenomics has led to the discovery of many new lineages at the phylum level. For the majority, only automatic genomic annotations can provide information regarding their metabolic potential and role in the environment. Here, genomic data from 2,978 archaeal genomes was used to perform automatic annotations using bioinformatics tools, alongside synteny analysis. These automatic classifications were done to assess how good these different tools perform in relation to archaeal data. Our study revealed that even with lowered cutoffs, several functional models do not capture the recently discovered archaeal diversity. Moreover, our investigation revealed that a significant portion of archaeal genomes, approximately 42%, remain uncharacterized. In comparison, within 3,235 bacterial genomes, a diverse range of unclassified proteins is obtained, with well-studied organisms like Escherichia coli having a substantially lower proportion of uncharacterized regions, ranging from <5 to 25%, and less studied lineages being comparable to archaea with the range of 35-40% of unclassified regions. Leveraging this analysis, we were able to identify metabolic protein markers, thereby providing insights into the metabolism of the archaea in our dataset. Our findings underscore a substantial gap between automatic classification tools and the comprehensive mapping of archaeal metabolism. Despite advances in computational approaches, a significant portion of archaeal genomes remains unexplored, highlighting the need for extensive experimental validation in this domain, as well as more refined annotation methods. This study contributes to a better understanding of archaeal metabolism and underscores the importance of further research in elucidating the functional potential of archaeal genomes.

Isolation of a methyl-reducing methanogen outside the Euryarchaeota

Methanogenic archaea are main contributors to methane emissions, and have a crucial role in carbon cycling and global warming. Until recently, methanogens were confined to Euryarchaeota, but metagenomic studies revealed the presence of genes encoding the methyl coenzyme M reductase complex in other archaeal clades1-4, thereby opening up the premise that methanogenesis is taxonomically more widespread. Nevertheless, laboratory cultivation of these non-euryarchaeal methanogens was lacking to corroborate their potential methanogenic ability and physiology. Here we report the isolation of a thermophilic archaeon LWZ-6 from an oil field. This archaeon belongs to the class Methanosuratincolia (originally affiliated with 'Candidatus Verstraetearchaeota') in the phylum Thermoproteota. Methanosuratincola petrocarbonis LWZ-6 is a strict hydrogen-dependent methylotrophic methanogen. Although previous metagenomic studies speculated on the fermentative potential of Methanosuratincolia members, strain LWZ-6 does not ferment sugars, peptides or amino acids. Its energy metabolism is linked only to methanogenesis, with methanol and monomethylamine as electron acceptors and hydrogen as an electron donor. Comparative (meta)genome analysis confirmed that hydrogen-dependent methylotrophic methanogenesis is a widespread trait among Methanosuratincolia. Our findings confirm that the diversity of methanogens expands beyond the classical Euryarchaeota and imply the importance of hydrogen-dependent methylotrophic methanogenesis in global methane emissions and carbon cycle.

Kinetic characterization and phylogenetic analysis of human ADP-dependent glucokinase reveal new insights into its regulatory properties

Although ADP-dependent sugar kinases were first described in archaea, at present, the presence of an ADP-dependent glucokinase (ADP-GK) in mammals is well documented. This enzyme is mainly expressed in hematopoietic lineages and tumor tissues, although its role has remained elusive. Here, we report a detailed kinetic characterization of the human ADP-dependent glucokinase (hADP-GK), addressing the influence of a putative signal peptide for endoplasmic reticulum (ER) destination by characterizing a truncated form. The truncated form revealed no significant impact on the kinetic parameters, showing only a slight increase in the Vmax value, higher metal promiscuity, and the same nucleotide specificity as the full-length enzyme. hADP-GK presents an ordered sequential kinetic mechanism in which MgADP is the first substrate to bind and AMP is the last product released, being the same mechanism described for archaeal ADP-dependent sugar kinases, in agreement with the protein topology. Substrate inhibition by glucose was observed due to sugar binding to nonproductive species. Although Mg²⁺ is an essential component for kinase activity, it also behaves as a partial mixed-type inhibitor for hADP-GK, mainly by decreasing the MgADP affinity. Regarding its distribution, phylogenetic analysis shows that ADP-GK´s are present in a wide diversity of eukaryotic organisms although it is not ubiquitous. Eukaryotic ADP-GKs sequences cluster into two main groups, showing differences in the highly conserved sugar-binding motif reported for archaeal enzymes [NX(N)XD] where a cysteine residue is found instead of asparagine in a significant number of enzymes. Site directed mutagenesis of the cysteine residue by asparagine produces a 6-fold decrease in V_max, suggesting a role for this residue in the catalytic process, probably by facilitating the proper orientation of the substrate to be phosphorylated.

Discovery of the Streamlined Haloarchaeon Halorutilus salinus, Comprising a New Order Widespread in Hypersaline Environments across the World

The class Halobacteria is one of the most diverse groups within the Euryarchaeota phylum, whose members are ubiquitously distributed in hypersaline environments, where they often constitute the major population. Here, we report the discovery and isolation of a new halophilic archaeon, strain F3-133^T exhibiting ≤86.3% 16S rRNA gene identity to any previously cultivated archaeon, and, thus, representing a new order. Analysis of available 16S rRNA gene amplicon and metagenomic data sets showed that the new isolate represents an abundant group in intermediate-to-high salinity ecosystems and is widely distributed across the world. The isolate presents a streamlined genome, which probably accounts for its ecological success in nature and its fastidious growth in culture. The predominant osmoprotection mechanism appears to be the typical salt-in strategy used by other haloarchaea. Furthermore, the genome contains the complete gene set for nucleotide monophosphate degradation pathway through archaeal RuBisCO, being within the first halophilic archaea representatives reported to code this enzyme. Genomic comparisons with previously described representatives of the phylum Euryarchaeota were consistent with the 16S rRNA gene data in supporting that our isolate represents a novel order within the class Halobacteria for which we propose the names Halorutilales ord. nov., Halorutilaceae fam. nov., Halorutilus gen. nov. and Halorutilus salinus sp. nov. IMPORTANCE The discovery of the new halophilic archaeon, Halorutilus salinus, representing a novel order, family, genus, and species within the class Halobacteria and phylum Euryarchaeota clearly enables insights into the microbial dark matter, expanding the current taxonomical knowledge of this group of archaea. The in-depth comparative genomic analysis performed on this new taxon revealed one of the first known examples of an Halobacteria representative coding the archaeal RuBisCO gene and with a streamlined genome, being ecologically successful in nature and explaining its previous non-isolation. Altogether, this research brings light into the understanding of the physiology of the Halobacteria class members, their ecological distribution, and capacity to thrive in hypersaline environments.

The Fate of Duplicated Enzymes in Prokaryotes: The Case of Isomerases

The isomerases are a unique enzymatic class of enzymes that carry out a great diversity of chemical reactions at the intramolecular level. This class comprises about 300 members, most of which are involved in carbohydrate and terpenoid/polyketide metabolism. Along with oxidoreductases and translocases, isomerases are one of the classes with the highest ratio of paralogous enzymes. Due to its relatively small number of members, it is plausible to explore it in greater detail to identify specific cases of gene duplication. Here, we present an analysis at the level of individual isomerases and identify different members that seem to be involved in duplication events in prokaryotes. As was suggested in a previous study, there is no homogeneous distribution of paralogs, but rather they accumulate into a few subcategories, some of which differ between Archaea and Bacteria. As expected, the metabolic processes with more paralogous isomerases have to do with carbohydrate metabolism but also with RNA modification (a particular case involving an rRNA-modifying isomerase is thoroughly discussed and analyzed in detail). Overall, our findings suggest that the most common fate for paralogous enzymes is the retention of the original enzymatic function, either associated with a dosage effect or with differential expression in response to changing environments, followed by subfunctionalization and, to a much lesser degree, neofunctionalization, which is consistent with what has been reported elsewhere.

Functional analyses of a highly thermostable hexokinase from Pyrobaculum calidifontis

The gene encoding a repressor open reading frame sugar kinase (ROK) family protein from hyperthermophilic crenarchaeon Pyrobaculum calidifontis, Pcal-HK, was cloned and expressed in Escherichia coli. The recombinant protein was produced in soluble and highly active form. Purified Pcal-HK was highly thermostable and existed in a monomeric form in solution. The enzyme was specific to ATP as phosphoryl donor but showed broad specificity to phosphoryl acceptors. It catalyzed the phosphorylation of a number of hexoses, including glucose, glucosamine, N-acetyl glucosamine, fructose and mannose, at nearly the same rate and similar affinity. The enzyme was metal ion dependent exhibiting highest activity at 90-95 °C and pH 8.5. Mg²⁺ was most effective metal ion, which could be partially replaced by Mn²⁺, Ni²⁺ or Zn²⁺. Kinetic parameters were determined at 90 °C and the enzyme showed almost similar catalytic efficiency (k_cat/K_m) towards the above mentioned hexoses. To the best of our knowledge, Pcal-HK is the most active thermostable ROK family hexokinase characterized to date which catalyzes the phosphorylation of various hexoses with nearly similar affinity.

A non-carboxylating pentose bisphosphate pathway in halophilic archaea

Bacteria and Eucarya utilize the non-oxidative pentose phosphate pathway to direct the ribose moieties of nucleosides to central carbon metabolism. Many archaea do not possess this pathway, and instead, Thermococcales utilize a pentose bisphosphate pathway involving ribose-1,5-bisphosphate (R15P) isomerase and ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (Rubisco). Intriguingly, multiple genomes from halophilic archaea seem only to harbor R15P isomerase, and do not harbor Rubisco. In this study, we identify a previously unrecognized nucleoside degradation pathway in halophilic archaea, composed of guanosine phosphorylase, ATP-dependent ribose-1-phosphate kinase, R15P isomerase, RuBP phosphatase, ribulose-1-phosphate aldolase, and glycolaldehyde reductase. The pathway converts the ribose moiety of guanosine to dihydroxyacetone phosphate and ethylene glycol. Although the metabolic route from guanosine to RuBP via R15P is similar to that of the pentose bisphosphate pathway in Thermococcales, the downstream route does not utilize Rubisco and is unique to halophilic archaea.

Genomic Insights Into New Species of the Genus Halomicroarcula Reveals Potential for New Osmoadaptative Strategies in Halophilic Archaea

Metagenomic studies on prokaryotic diversity of hypersaline soils from the Odiel saltmarshes, South-west Spain, revealed a high proportion of genomic sequences not related to previously cultivated taxa, that might be related to haloarchaea with a high environmental and nutritional flexibility. In this study, we used a culturomics approach in order to isolate new haloarchaeal microorganisms from these hypersaline soils. Four haloarchaeal strains, designated strains F24A^T, F28, F27^T, and F13^T, phylogenetically related to the genus Halomicroarcula, were isolated and characterized in detail. The phylogenomic tree based on the 100 orthologous single-copy genes present in the genomes of these four strains as well as those of the type strains of the species Halomicroarcula pellucida CECT 7537^T, Halomicroarcula salina JCM 18369^T and Halomicroarcula limicola JCM 18640^T, that were determined in this study, revealed that these four new isolates clustered on three groups, with strains F24A^T and F28 within a single cluster, and altogether with the species of Halomicroarcula. Additionally, Orthologous Average Nucleotide Identity (OrthoANI), digital DNA-DNA hybridization (dDDH) and Average Amino-acid Identity (AAI) values, likewise phenotypic characteristics, including their polar lipids profiles, permitted to determine that they represent three new species, for which we propose the names Halomicroarcula rubra sp. nov. (type strain F13^T), Halomicroarcula nitratireducens sp. nov. (type strain F27^T) and Halomicroarcula salinisoli sp. nov. (type strain F24A^T). An in deep comparative genomic analysis of species of the genus Halomicroarcula, including their metabolism, their capability to biosynthesize secondary metabolites and their osmoregulatory adaptation mechanisms was carried out. Although they use a salt-in strategy, the identification of the complete pathways for the biosynthesis of the compatible solutes trehalose and glycine betaine, not identified before in any other haloarchaea, might suggest alternative osmoadaptation strategies for this group. This alternative osmoregulatory mechanism would allow this group of haloarchaea to be versatile and eco-physiologically successful in hypersaline environments and would justify the capability of the species of this genus to grow not only on environments with high salt concentrations [up to 30% (w/v) salts], but also under intermediate to low salinities.

The Piezo-Hyperthermophilic Archaeon Thermococcus piezophilus Regulates Its Energy Efficiency System to Cope With Large Hydrostatic Pressure Variations

Deep-sea ecosystems share a common physical parameter, namely high hydrostatic pressure (HHP). Some of the microorganisms isolated at great depths have a high physiological plasticity to face pressure variations. The adaptive strategies by which deep-sea microorganisms cope with HHP variations remain to be elucidated, especially considering the extent of their biotopes on Earth. Herein, we investigated the gene expression patterns of Thermococcus piezophilus, a piezohyperthermophilic archaeon isolated from the deepest hydrothermal vent known to date, under sub-optimal, optimal and supra-optimal pressures (0.1, 50, and 90 MPa, respectively). At stressful pressures [sub-optimal (0.1 MPa) and supra-optimal (90 MPa) conditions], no classical stress response was observed. Instead, we observed an unexpected transcriptional modulation of more than a hundred gene clusters, under the putative control of the master transcriptional regulator SurR, some of which are described as being involved in energy metabolism. This suggests a fine-tuning effect of HHP on the SurR regulon. Pressure could act on gene regulation, in addition to modulating their expression.

Enzymes of an alternative pathway of glucose metabolism in obligate methanotrophs

Aerobic methanotrophic bacteria utilize methane as a growth substrate but are unable to grow on any sugars. In this study we have shown that two obligate methanotrophs, Methylotuvimicrobium alcaliphilum 20Z and Methylobacter luteus IMV-B-3098, possess functional glucose dehydrogenase (GDH) and gluconate kinase (GntK). The recombinant GDHs from both methanotrophs were homotetrameric and strongly specific for glucose preferring NAD⁺ over NADP⁺. GDH from Mtm. alcaliphilum was most active at pH 10 (V_max = 95 U/mg protein) and demonstrated very high K_m for glucose (91.8 ± 3.8 mM). GDH from Mb. luteus was most active at pH 8.5 (V_max = 43 U/mg protein) and had lower K_m for glucose (16 ± 0.6 mM). The cells of two Mtm. alcaliphilum double mutants with deletions either of the genes encoding GDH and glucokinase (gdh^─/glk^─) or of the genes encoding gluconate kinase and glucokinase (gntk^─/glk^─) had the lower glycogen level and the higher contents of intracellular glucose and trehalose compared to the wild type strain. The gntk^─/glk^─ knockout mutant additionally accumulated gluconic acid. These data, along with bioinformatics analysis, demonstrate that glycogen derived free glucose can enter the Entner–Doudoroff pathway or the pentose phosphate cycle in methanotrophs, bypassing glycolysis via the gluconate shunt.

ADP-dependent glucose/glucosamine kinase from Thermococcus kodakarensis: cloning and characterization

The genome sequence of Thermococcus kodakarensis contains an open reading frame, TK1110, annotated as ADP-dependent glucokinase. The encoding gene was expressed in Escherichia coli and the gene product, TK-GLK, was produced in soluble and active form. The recombinant enzyme was extremely thermostable. Thermostability was increased significantly in the presence of ammonium sulfate. ADP was the preferred co-factor for TK-GLK, which could be replaced with CDP but with a 60% activity. TK-GLK was a metal ion-dependent enzyme which exhibited glucokinase, glucosamine kinase and glucose 6-phosphatase activities. It catalyzed the phosphorylation of both glucose and glucosamine with nearly the same rate and affinity. The apparent K_m values for glucose and glucosamine were 0.48 ± 0.03 and 0.47 ± 0.09 mM, respectively. The catalytic efficiency (k_cat/K_m) values against these two substrates were 6.2 × 10⁵ ± 0.25 and 5.8 × 10⁵ ± 0.75 M^-1 s^-1. The apparent K_m value for dephosphorylation of glucose 6-phosphate was ~14-fold higher than that of glucose phosphorylation. Similarly, catalytic efficiency (k_cat/K_m) for phosphatase reaction was ~19-fold lower than that for the kinase reaction. To the best of our knowledge, this is the first report that describes the reversible nature of a euryarchaeal ADP-dependent glucokinase.

Nonenzymatic Metabolic Reactions and Life's Origins

Prebiotic chemistry aims to explain how the biochemistry of life as we know it came to be. Most efforts in this area have focused on provisioning compounds of importance to life by multistep synthetic routes that do not resemble biochemistry. However, gaining insight into why core metabolism uses the molecules, reactions, pathways, and overall organization that it does requires us to consider molecules not only as synthetic end goals. Equally important are the dynamic processes that build them up and break them down. This perspective has led many researchers to the hypothesis that the first stage of the origin of life began with the onset of a primitive nonenzymatic version of metabolism, initially catalyzed by naturally occurring minerals and metal ions. This view of life's origins has come to be known as "metabolism first". Continuity with modern metabolism would require a primitive version of metabolism to build and break down ketoacids, sugars, amino acids, and ribonucleotides in much the same way as the pathways that do it today. This review discusses metabolic pathways of relevance to the origin of life in a manner accessible to chemists, and summarizes experiments suggesting several pathways might have their roots in prebiotic chemistry. Finally, key remaining milestones for the protometabolic hypothesis are highlighted.

Alternative carbohydrate pathways - enzymes, functions and engineering

Metabolic engineering is crucial in the development of production strains for platform chemicals, pharmaceuticals and biomaterials from renewable resources. The central carbon metabolism (CCM) of heterotrophs plays an essential role in the conversion of biomass to the cellular building blocks required for growth. Yet, engineering the CCM ultimately aims toward a maximization of flux toward products of interest. The most abundant dissimilative carbohydrate pathways amongst prokaryotes (and eukaryotes) are the Embden-Meyerhof-Parnas (EMP) and the Entner-Doudoroff (ED) pathways, which build the basics for heterotrophic metabolic chassis strains. Although the EMP is regarded as the textbook example of a carbohydrate pathway owing to its central role in production strains like Escherichia coli, Saccharomyces cerevisiae and Bacillus subtilis, it is either modified, complemented or even replaced by alternative carbohydrate pathways in different organisms. The ED pathway also plays key roles in biotechnological relevant bacteria, like Zymomonas mobilis and Pseudomonas putida, and its importance was recently discovered in photoautotrophs and marine microorganisms. In contrast to the EMP, the ED pathway and its variations are not evolutionary optimized for high ATP production and it differs in key principles such as protein cost, energetics and thermodynamics, which can be exploited in the construction of unique metabolic designs. Single ED pathway enzymes and complete ED pathway modules have been used to rewire carbon metabolisms in production strains and for the construction of cell-free enzymatic pathways. This review focuses on the differences of the ED and EMP pathways including their variations and discusses the use of alternative pathway strategies for in vivo and cell-free metabolic engineering.

Culturomics-based genomics sheds light on the ecology of the new haloarchaeal genus Halosegnis

The development of culture-independent techniques has revolutionized our understanding of microbial ecology, especially through the illustration of the vast gap between the environmentally abundant microbial diversity and that accessible through cultivation. However, culture-based approaches are not only crucial for understanding the evolutionary, metabolic and ecological milieu of microbial diversity but also for the development of novel biotechnological applications. In this study, we used a culturomics-based approach in order to isolate novel microbial taxa from hypersaline environments (i.e. Isla Cristina and Isla Bacuta salterns in Huelva, Spain). We managed to obtain axenic cultures of four haloarchaeal strains that belong to a new haloarchaeal genus and to obtain their genomic sequences. The phylogenomic and phylogenetic analyses (together with AAI, ANI and digital DDH indices) showed that the isolates constitute two new species, for which we propose the names Halosegnis longus sp. nov. and Halosegnis rubeus sp. nov. The genomic-based metabolic reconstructions indicated that members of this new haloarchaeal genus have photoheterotrophic aerobic lifestyle with a typical salt-in signature. 16S rRNA gene sequence reads abundance profiles and genomic recruitment analyses revealed that the Halosegnis genus has a worldwide geographical distribution, reaching high abundance (up to 8%) in habitats with intermediate salinities.

Effects of high-level expression of A1-ATPase on H2 production in Thermococcus kodakarensis

The hyperthermophilic archaeon Thermococcus kodakarensis can grow on pyruvate or maltooligosaccharides through H₂ fermentation. H₂ production levels of members of the Thermococcales are high, and studies to improve their production potential have been reported. Although H₂ production is primary metabolism, here we aimed to partially uncouple cell growth and H₂ production of T. kodakarensis. Additional A₁-type ATPase genes were introduced into T. kodakarensis KU216 under the control of two promoters; the strong constitutive cell surface glycoprotein promoter, P_csg, and the sugar-inducible fructose-1,6-bisphosphate aldolase promoter, P_fba. Whereas cells with the A₁-type ATPase genes under the control of P_csg displayed only trace levels of growth, cells with P_fba (strain KUA-PF) displayed growth sufficient for further analysis. Increased levels of A₁-type ATPase protein were detected in KUA-PF cells grown on pyruvate or maltodextrin, when compared to the levels in the host strain KU216. The growth and H₂ production levels of strain KUA-PF with pyruvate or maltodextrin as a carbon and electron source were analyzed and compared to those of the host strain KU216. Compared to a small decrease in total H₂ production, significantly larger decreases in cell growth were observed, resulting in an increase in cell-specific H₂ production. Quantification of the substrate also revealed that ATPase overexpression led to increased cell-specific pyruvate and maltodextrin consumptions. The results clearly indicate that ATPase production results in partial uncoupling of cell growth and H₂ production in T. kodakarensis.

The ubiquity and coexistence of two FBPases in chloroplasts of photosynthetic eukaryotes and its evolutionary and functional implications

In photosynthetic eukaryotes, there are two well-characterized fructose-1,6-bisphosphatases (FBPases): the redox-insensitive cytosolic FBPase (cyFBPase), which participates in gluconeogenesis, and the redox-sensitive chloroplastic FBPase (cpFBPase1), which is a critical enzyme in the Calvin cycle. Recent studies have identified a new chloroplastic FBPase, cpFBPase2; however, its phylogenetic distribution, evolutionary origin, and physiological function remain unclear. In this study, we identified and characterized these three FBPase isoforms in diverse, representative photosynthetic lineages and analyzed their phylogeny. In contrast to previous hypotheses, we found that cpFBPase2 is ubiquitous in photosynthetic eukaryotes. Additionally, all cpFBPase2s from diverse lineages form a monophyly, suggesting cpFBPase2 is not a recently evolved enzyme restricted to land plants but rather evolved early in the evolution of photosynthetic organisms, and most likely, in the common ancestor of photosynthetic eukaryotes. cyFBPase was probably first duplicated to produce cpFBPase2, and then the latter duplicated to produce cpFBPase1. The ubiquitous coexistence of these two cpFBPases in chloroplasts is most likely the consequence of adaptation to different redox conditions of photosynthesis, especially those caused by recurrent changes in light conditions.

Genomic Evidence for Formate Metabolism by Chloroflexi as the Key to Unlocking Deep Carbon in Lost City Microbial Ecosystems

Primitive forms of life may have originated around hydrothermal vents at the bottom of the ancient ocean. The Lost City hydrothermal vent field, fueled by just rock and water, provides an analog for not only primitive ecosystems but also potential extraterrestrial rock-powered ecosystems. The microscopic life covering the towering chimney structures at the Lost City has been previously documented, yet little is known about the carbon cycling in this ecosystem. These results provide a better understanding of how carbon from the deep subsurface can fuel rich microbial ecosystems on the seafloor.

The Lost City hydrothermal field on the Mid-Atlantic Ridge supports dense microbial life on the lofty calcium carbonate chimney structures. The vent field is fueled by chemical reactions between the ultramafic rock under the chimneys and ambient seawater. These serpentinization reactions provide reducing power (as hydrogen gas) and organic compounds that can serve as microbial food; the most abundant of these are methane and formate. Previous studies have characterized the interior of the chimneys as a single-species biofilm inhabited by the Lost City Methanosarcinales, but they also indicated that this methanogen is unable to metabolize formate. The new metagenomic results presented here indicate that carbon cycling in these Lost City chimney biofilms could depend on the metabolism of formate by Chloroflexi populations. Additionally, we present evidence for metabolically diverse, formate-utilizing Sulfurovum populations and new genomic and phylogenetic insights into the unique Lost City Methanosarcinales.

IMPORTANCE Primitive forms of life may have originated around hydrothermal vents at the bottom of the ancient ocean. The Lost City hydrothermal vent field, fueled by just rock and water, provides an analog for not only primitive ecosystems but also potential extraterrestrial rock-powered ecosystems. The microscopic life covering the towering chimney structures at the Lost City has been previously documented, yet little is known about the carbon cycling in this ecosystem. These results provide a better understanding of how carbon from the deep subsurface can fuel rich microbial ecosystems on the seafloor.

Integrated thermodynamic analysis of electron bifurcating [FeFe]-hydrogenase to inform anaerobic metabolism and H2 production

Electron bifurcating, [FeFe]-hydrogenases are recently described members of the hydrogenase family and catalyze a combination of exergonic and endergonic electron exchanges between three carriers (2 ferredoxin_red^- + NAD(P)H + 3 H⁺ = 2 ferredoxin_ox + NAD(P)⁺ + 2 H₂). A thermodynamic analysis of the bifurcating, [FeFe]-hydrogenase reaction, using electron path-independent variables, quantified potential biological roles of the reaction without requiring enzyme details. The bifurcating [FeFe]-hydrogenase reaction, like all bifurcating reactions, can be written as a sum of two non-bifurcating reactions. Therefore, the thermodynamic properties of the bifurcating reaction can never exceed the properties of the individual, non-bifurcating, reactions. The bifurcating [FeFe]-hydrogenase reaction has three competitive properties: 1) enabling NAD(P)H-driven proton reduction at pH₂ higher than the concurrent operation of the two, non-bifurcating reactions, 2) oxidation of NAD(P)H and ferredoxin simultaneously in a 1:1 ratio, both are produced during typical glucose fermentations, and 3) enhanced energy conservation (~10 kJ mol^-1 H₂) relative to concurrent operation of the two, non-bifurcating reactions. Our analysis demonstrated ferredoxin E°' largely determines the sensitivity of the bifurcating reaction to pH₂, modulation of the reduced/oxidized electron carrier ratios contributed less to equilibria shifts. Hydrogenase thermodynamics data were integrated with typical and non-typical glycolysis pathways to evaluate achieving the 'Thauer limit' (4 H₂ per glucose) as a function of temperature and pH₂. For instance, the bifurcating [FeFe]-hydrogenase reaction permits the Thauer limit at 60 °C if pH ₂ ≤ ~10 mbar. The results also predict Archaea, expressing a non-typical glycolysis pathway, would not benefit from a bifurcating [FeFe]-hydrogenase reaction; interestingly, no Archaea have been observed experimentally with a [FeFe]-hydrogenase enzyme.

Potential regulation by miRNAs on glucose metabolism in liver of common carp (Cyprinus carpio) at different temperatures

Water temperature can affect the metabolism of fish. Common carp (Cyprinus carpio) is a representative eurythermic fish that can survive at a wide range of ambient temperatures, allowing it to live in an extensive geographical range. The goal of this work was to study the glucose metabolism of common carp at different temperatures and determine the miRNAs involved in the regulation of glucose metabolism. We determined the indicators related to glucose metabolism after long-term temperature stress and constructed nine small RNA libraries of livers under different temperature stress (5 °C, 17 °C, and 30 °C, with three biological replicates for each temperature), and subjected these samples to high-throughput sequencing. A positive relationship was observed between weight gain rate (WGR) and temperature increase after 18 days of temperature stress. However, the glucose level in the plasma maintained a gentle decrease. Unexpectedly, liver lactic acid levels were elevated in HTG (high temperature group) and LTG (low temperature group). Six down-regulated miRNAs (miR-122, miR-30b, miR-15b-5p, miR-20a-5p, miR-1, and miR-7b) were identified as involved in the regulation of glycolysis. Twelve genes were predicted as targets of these miRNAs, and these genes are in pathways related to pyruvate metabolism, glycolysis/gluconeogenesis, and the citrate cycle (TCA cycle). The results allowed prediction of a potential regulatory network of miRNAs involved in the regulation of glycolysis. The target genes of six down-regulated miRNAs were up-regulated under temperature stress, including Aldolase C, fructose-bisphosphate, b (ALDOCB), multiple inositol-polyphosphate phosphatase 1 (MINPP1), phosphoenolpyruvate carboxykinase 1 (PCK1), pyruvate dehydrogenase E1 alpha 1 (PDHA1), aldehyde dehydrogenase 9 family member A1a (ALDH9A1A), Acetyl-coenzyme A synthetase (ACSS), lactate dehydrogenase b (LDH-b), and glyoxylate reductase/hydroxypyruvate reductase (GRHPR). Other key genes of glycolysis, glucose transporter 1 (GLUT-1), pyruvate kinase PKM (PKM), and mitochondrial pyruvate carrier (MPC) were significantly up-regulated in LTG and HTG. Overall, the results suggest that miRNAs maintain their energy requirements by regulating glycolysis and play an important role in the molecular response to cold and heat stress of common carp. These data provide the foundation for further studies of the role of miRNAs in environmental adaptation in fish.

Unusual Phosphoenolpyruvate (PEP) Synthetase-Like Protein Crucial to Enhancement of Polyhydroxyalkanoate Accumulation in Haloferax mediterranei Revealed by Dissection of PEP-Pyruvate Interconversion Mechanism

Phosphoenolpyruvate (PEP)/pyruvate interconversion is a major metabolic point in glycolysis and gluconeogenesis and is catalyzed by various sets of enzymes in different Archaea groups. In this study, we report the key enzymes that catalyze the anabolic and catabolic directions of the PEP/pyruvate interconversion in Haloferax mediterranei The in silico analysis showed the presence of a potassium-dependent pyruvate kinase (PYK_Hm [HFX_0773]) and two phosphoenol pyruvate synthetase (PPS) candidates (PPS_Hm [HFX_0782] and a PPS homolog protein named PPS-like [HFX_2676]) in this strain. Expression of the pyk _Hm gene and pps _Hm was induced by glycerol and pyruvate, respectively; whereas the pps-like gene was not induced at all. Similarly, genetic analysis and enzyme activities of purified proteins showed that PYK_Hm catalyzed the conversion from PEP to pyruvate and that PPS_Hm catalyzed the reverse reaction, while PPS-like protein displayed no function in PEP/pyruvate interconversion. Interestingly, knockout of the pps-like gene led to a 70.46% increase in poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) production. The transcriptome sequencing (RNA-Seq) and quantitative reverse transcription-PCR (qRT-PCR) results showed that many genes responsible for PHBV monomer supply and for PHBV synthesis were upregulated in a pps-like gene deletion strain and thereby improved PHBV accumulation. Additionally, our phylogenetic evidence suggested that PPS-like protein diverged from PPS enzyme and evolved as a distinct protein with novel function in haloarchaea. Our findings attempt to fill the gaps in central metabolism of Archaea by providing comprehensive information about key enzymes involved in the haloarchaeal PEP/pyruvate interconversion, and we also report a high-yielding PHBV strain with great future potentials.IMPORTANCE Archaea, the third domain of life, have evolved diversified metabolic pathways to cope with their extreme habitats. Phosphoenol pyruvate (PEP)/pyruvate interconversion during carbohydrate metabolism is one such important metabolic process that is highly differentiated among Archaea However, this process is still uncharacterized in the haloarchaeal group. Haloferax mediterranei is a well-studied haloarchaeon that has the ability to produce polyhydroxyalkanoates (PHAs) under unbalanced nutritional conditions. In this study, we identified the key enzymes involved in this interconversion and discussed their differences with their counterparts from other members of the Archaea and Bacteria domains. Notably, we found a novel protein, phosphoenolpyruvate synthetase-like (PPS-like), which exhibited high homology to PPS enzyme. However, PPS-like protein has evolved some distinct sequence features and functions, and strikingly the corresponding gene deletion helped to enhance poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) synthesis significantly. Overall, we have filled the gap in knowledge about PEP/pyruvate interconversion in haloarchaea and reported an efficient strategy for improving PHBV production in H. mediterranei.

A synthetic glycerol assimilation pathway demonstrates biochemical constraints of cellular metabolism

The engineering of synthetic metabolic routes can provide valuable lessons on the roles of different biochemical constraints in shaping pathway activity. In this study, we designed and engineered a novel glycerol assimilation pathway in Escherichia coli. While the synthetic pathway was based only on well-characterized endogenous reactions, we were not able to establish robust growth using standard concentrations of glycerol. Long-term evolution failed to improve growth via the pathway, indicating that this limitation was not regulatory but rather relates to fundamental aspects of cellular metabolism. We show that the activity of the synthetic pathway is fully controlled by three key physicochemical constraints: thermodynamics, kinetics and metabolite toxicity. Overcoming a thermodynamic barrier at the beginning of the pathway requires high glycerol concentrations. A kinetic barrier leads to a Monod-like growth dependency on substrate concentration, but with a very high substrate saturation constant. Finally, the flat thermodynamic profile of the pathway enforces a pseudoequilibrium between glycerol and the reactive intermediate dihydroxyacetone, which inhibits growth when the feedstock concentration surpasses 1000 mm. Overall, this study serves to demonstrate the use of synthetic biology to elucidate key design principles of cellular metabolism.

Growth Kinetics, Carbon Isotope Fractionation, and Gene Expression in the Hyperthermophile Methanocaldococcus jannaschii during Hydrogen-Limited Growth and Interspecies Hydrogen Transfer

Hyperthermophilic methanogens are often H2 limited in hot subseafloor environments, and their survival may be due in part to physiological adaptations to low H2 conditions and interspecies H2 transfer. The hyperthermophilic methanogen Methanocaldococcus jannaschii was grown in monoculture at high (80 to 83 μM) and low (15 to 27 μM) aqueous H2 concentrations and in coculture with the hyperthermophilic H2 producer Thermococcus paralvinellae The purpose was to measure changes in growth and CH4 production kinetics, CH4 fractionation, and gene expression in M. jannaschii with changes in H2 flux. Growth and cell-specific CH4 production rates of M. jannaschii decreased with decreasing H2 availability and decreased further in coculture. However, cell yield (cells produced per mole of CH4 produced) increased 6-fold when M. jannaschii was grown in coculture rather than monoculture. Relative to high H2 concentrations, isotopic fractionation of CO2 to CH4 (εCO2-CH4) was 16‰ larger for cultures grown at low H2 concentrations and 45‰ and 56‰ larger for M. jannaschii growth in coculture on maltose and formate, respectively. Gene expression analyses showed H2-dependent methylene-tetrahydromethanopterin (H4MPT) dehydrogenase expression decreased and coenzyme F420-dependent methylene-H4MPT dehydrogenase expression increased with decreasing H2 availability and in coculture growth. In coculture, gene expression decreased for membrane-bound ATP synthase and hydrogenase. The results suggest that H2 availability significantly affects the CH4 and biomass production and CH4 fractionation by hyperthermophilic methanogens in their native habitats.IMPORTANCE Hyperthermophilic methanogens and H2-producing heterotrophs are collocated in high-temperature subseafloor environments, such as petroleum reservoirs, mid-ocean ridge flanks, and hydrothermal vents. Abiotic flux of H2 can be very low in these environments, and there is a gap in our knowledge about the origin of CH4 in these habitats. In the hyperthermophile Methanocaldococcus jannaschii, growth yields increased as H2 flux, growth rates, and CH4 production rates decreased. The same trend was observed increasingly with interspecies H2 transfer between M. jannaschii and the hyperthermophilic H2 producer Thermococcus paralvinellae With decreasing H2 availability, isotopic fractionation of carbon during methanogenesis increased, resulting in isotopically more negative CH4 with a concomitant decrease in H2-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression and increase in F420-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression. The significance of our research is in understanding the nature of hyperthermophilic interspecies H2 transfer and identifying biogeochemical and molecular markers for assessing the physiological state of methanogens and possible source of CH4 in natural environments.

The Carbon Switch at the Level of Pyruvate and Phosphoenolpyruvate in Sulfolobus solfataricus P2

Sulfolobus solfataricus P2 grows on different carbohydrates as well as alcohols, peptides and amino acids. Carbohydrates such as D-glucose or D-galactose are degraded via the modified, branched Entner–Doudoroff (ED) pathway whereas growth on peptides requires the Embden–Meyerhof–Parnas (EMP) pathway for gluconeogenesis. As for most hyperthermophilic Archaea an important control point is established at the level of triosephophate conversion, however, the regulation at the level of pyruvate/phosphoenolpyruvate conversion was not tackled so far. Here we describe the cloning, expression, purification and characterization of the pyruvate kinase (PK, SSO0981) and the phosphoenolpyruvate synthetase (PEPS, SSO0883) of Sul. solfataricus. The PK showed only catabolic activity [catalytic efficiency (PEP): 627.95 mM^-1s^-1, 70°C] with phosphoenolpyruvate as substrate and ADP as phosphate acceptor and was allosterically inhibited by ATP and isocitrate (K_i 0.8 mM). The PEPS was reversible, however, exhibited preferred activity in the gluconeogenic direction [catalytic efficiency (pyruvate): 1.04 mM^-1s^-1, 70°C] and showed some inhibition by AMP and α-ketoglutarate. The gene SSO2829 annotated as PEPS/pyruvate:phosphate dikinase (PPDK) revealed neither PEPS nor PPDK activity. Our studies suggest that the energy charge of the cell as well as the availability of building blocks in the citric acid cycle and the carbon/nitrogen balance plays a major role in the Sul. solfataricus carbon switch. The comparison of regulatory features of well-studied hyperthermophilic Archaea reveals a close link and sophisticated coordination between the respective sugar kinases and the kinetic and regulatory properties of the enzymes at the level of PEP-pyruvate conversion.

Pareto Optimality Explanation of the Glycolytic Alternatives in Nature

The Entner-Doudoroff (ED) and Embden-Meyerhof-Parnas (EMP) glycolytic pathways are largely conserved across glycolytic species in nature. Is this a coincidence, convergent evolution or there exists a driving force towards either of the two pathway designs? We addressed this question by first employing a variant of the optStoic algorithm to exhaustively identify over 11,916 possible routes between glucose and pyruvate at different pre-determined stoichiometric yields of ATP. Subsequently, we analyzed the thermodynamic feasibility of all the pathways at physiological metabolite concentrations and quantified the protein cost of the feasible solutions. Pareto optimality analysis between energy efficiency and protein cost reveals that the naturally evolved ED and EMP pathways are indeed among the most protein cost-efficient pathways in their respective ATP yield categories and remain thermodynamically feasible across a wide range of ATP/ADP ratios and pathway intermediate metabolite concentration ranges. In contrast, pathways with higher ATP yield (>2) while feasible, are bound within stringent and often extreme operability ranges of cofactor and intermediate metabolite concentrations. The preponderance of EMP and ED is thus consistent with not only optimally balancing energy yield vs. enzyme cost but also with ensuring operability for wide metabolite concentration ranges and ATP/ADP ratios.

Metatranscriptome analysis of the intestinal microorganisms in Pardosa pseudoannulata in response to cadmium stress

Cadmium (Cd) generates a variety of physiological and ecological toxicity to spiders. However, little is known about the effects of Cd on symbiotic bacteria of spiders. Metatranscriptomics is increasing our knowledge of microorganisms in environment. To better understand the impact of Cd on the symbiotic bacteria of spiders, we generated and compared the metatranscriptomes of the intestinal microorganisms of Pardosa pseudoannulata with and without Cd stress. The community structure of intestinal microorganisms in P. pseudoannulata was composed of 4 kingdoms, namely bacteria, viruses, eukaryotes and archaea, including 46 phyla, 97 classes, 184 orders, 339 families, 470 genera, and 598 species. The abundance of eukaryotes, bacteria and viruses was decreased by 0.14%, 1.22% and 2.52% respectively while the archaea was increased by 99.16% when under Cd stress. We identified 1519 differentially expressed genes (DEGs), including 770 up-regulated and 749 down-regulated genes. The results of KEGG annotation revealed that the expression of genes that are involved in the carbon metabolism, protein and amino acid metabolism and synthesis, glucose metabolism, oxidative phosphorylation, and glutathione metabolism were influenced by Cd. Collectively, these findings showed that Cd significantly impacted the community structure and expression of related functional genes of intestinal microorganisms in P. pseudoannulata.

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

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

Metabolic reconstruction and experimental verification of glucose utilization in Desulfurococcus amylolyticus DSM 16532

Desulfurococcus amylolyticus DSM 16532 is an anaerobic and hyperthermophilic crenarchaeon known to grow on a variety of different carbon sources, including monosaccharides and polysaccharides. Furthermore, D. amylolyticus is one of the few archaea that are known to be able to grow on cellulose. Here, we present the metabolic reconstruction of D. amylolyticus’ central carbon metabolism. Based on the published genome, the metabolic reconstruction was completed by integrating complementary information available from the KEGG, BRENDA, UniProt, NCBI, and PFAM databases, as well as from available literature. The genomic analysis of D. amylolyticus revealed genes for both the classical and the archaeal version of the Embden-Meyerhof pathway. The metabolic reconstruction highlighted gaps in carbon dioxide-fixation pathways. No complete carbon dioxide-fixation pathway such as the reductive citrate cycle or the dicarboxylate-4-hydroxybutyrate cycle could be identified. However, the metabolic reconstruction indicated that D. amylolyticus harbors all genes necessary for glucose metabolization. Closed batch experimental verification of glucose utilization by D. amylolyticus was performed in chemically defined medium. The findings from in silico analyses and from growth experiments are discussed with respect to physiological features of hyperthermophilic organisms.

Exploring the combinatorial space of complete pathways to chemicals

Computational pathway design tools often face the challenges of balancing the stoichiometry of co-metabolites and cofactors, and dealing with reaction rule utilization in a single workflow. To this end, we provide an overview of two complementary stoichiometry-based pathway design tools optStoic and novoStoic developed in our group to tackle these challenges. optStoic is designed to determine the stoichiometry of overall conversion first which optimizes a performance criterion (e.g. high carbon/energy efficiency) and ensures a comprehensive search of co-metabolites and cofactors. The procedure then identifies the minimum number of intervening reactions to connect the source and sink metabolites. We also further the pathway design procedure by expanding the search space to include both known and hypothetical reactions, represented by reaction rules, in a new tool termed novoStoic. Reaction rules are derived based on a mixed-integer linear programming (MILP) compatible reaction operator, which allow us to explore natural promiscuous enzymes, engineer candidate enzymes that are not already promiscuous as well as design de novo enzymes. The identified biochemical reaction rules then guide novoStoic to design routes that expand the currently known biotransformation space using a single MILP modeling procedure. We demonstrate the use of the two computational tools in pathway elucidation by designing novel synthetic routes for isobutanol.

Metabolic versatility of small archaea Micrarchaeota and Parvarchaeota

Small acidophilic archaea belonging to Micrarchaeota and Parvarchaeota phyla are known to physically interact with some Thermoplasmatales members in nature. However, due to a lack of cultivation and limited genomes on hand, their biodiversity, metabolisms, and physiologies remain largely unresolved. Here, we obtained 39 genomes from acid mine drainage (AMD) and hot spring environments around the world. 16S rRNA gene based analyses revealed that Parvarchaeota were only detected in AMD and hot spring habitats, while Micrarchaeota were also detected in others including soil, peat, hypersaline mat, and freshwater, suggesting a considerable higher diversity and broader than expected habitat distribution for this phylum. Despite their small genomes (0.64-1.08 Mb), these archaea may contribute to carbon and nitrogen cycling by degrading multiple saccharides and proteins, and produce ATP via aerobic respiration and fermentation. Additionally, we identified several syntenic genes with homology to those involved in iron oxidation in six Parvarchaeota genomes, suggesting their potential role in iron cycling. However, both phyla lack biosynthetic pathways for amino acids and nucleotides, suggesting that they likely scavenge these biomolecules from the environment and/or other community members. Moreover, low-oxygen enrichments in laboratory confirmed our speculation that both phyla are microaerobic/anaerobic, based on several specific genes identified in them. Furthermore, phylogenetic analyses provide insights into the close evolutionary history of energy related functionalities between both phyla with Thermoplasmatales. These results expand our understanding of these elusive archaea by revealing their involvement in carbon, nitrogen, and iron cycling, and suggest their potential interactions with Thermoplasmatales on genomic scale.

On the Early Evolution of Catabolic Pathways: A Comparative Genomics Approach. I. The Cases of Glucose, Ribose, and the Nucleobases Catabolic Routes

Compared with the large corpus of published work devoted to the study of the origin and early development of anabolism, little attention has been given to the discussion of the early evolution of catabolism in spite of its significance. In the present study, we have used comparative genomics to explore the evolution and phylogenetic distribution of the enzymes that catalyze the extant catabolic pathways of the monosaccharides glucose and ribose, as well as those of the nucleobases adenine, guanine, cytosine, uracil, and thymine. Based on the oxygen dependence of the enzymes, their conservation, and evolution, we speculate on the relative antiquity of the pathways. Our results allow us to suggest which catabolic pathways and enzymes may have already been present in the last common ancestor. We conclude that the enzymatic degradations of ribose, as well as those of purines adenine and guanine, are among the most ancient catabolic pathways which can be traced by protein-based methodologies.

A review on comparative mechanistic studies of antimicrobial peptides against archaea

Archaea was until recently considered as a third domain of life in addition to bacteria and eukarya but recent studies support the existence of only two superphyla (bacteria and archaea). The fundamental differences between archaeal, bacterial, and eukaryal cells are probably the main reasons for the comparatively lower susceptibility of archaeal strains to current antimicrobial agents. The possible emerging pathogenicity of archaea and the role of archaeal methanogens in methane emissions, a potent greenhouse gas, has led many researchers to examine the sensitivity patterns of archaea and make attempts to find agents that have significant anti-archaeal activity. Even though antimicrobial peptides (AMPs) are well known with several published reviews concerning their mode of action against bacteria and eukarya, to our knowledge, to date no reviews are available that focus on the action of these peptides against archaea. Herein, we present a review on all the peptides that have been tested against archaea. In addition, in an attempt to shed more light on possible future work that needs to be performed we have included a brief overview of the chemical characteristics, spectrum of activity, and the known mechanism of action of each of these peptides against bacteria and/or fungi. We also discuss the nature of and key physiological differences between Archaea, Bacteria, and Eukarya that are relevant to the development of anti-archaeal peptides. Despite our relatively limited knowledge about archaea, available data suggest that AMPs have an even broader spectrum of activity than currently recognized.

Two functionally distinct NADP+-dependent ferredoxin oxidoreductases maintain the primary redox balance of Pyrococcus furiosus

Electron bifurcation has recently gained acceptance as the third mechanism of energy conservation in which energy is conserved through the coupling of exergonic and endergonic reactions. A structure-based mechanism of bifurcation has been elucidated recently for the flavin-based enzyme NADH-dependent ferredoxin NADP⁺ oxidoreductase I (NfnI) from the hyperthermophillic archaeon Pyrococcus furiosus. NfnI is thought to be involved in maintaining the cellular redox balance, producing NADPH for biosynthesis by recycling the two other primary redox carriers, NADH and ferredoxin. The P. furiosus genome encodes an NfnI paralog termed NfnII, and the two are differentially expressed, depending on the growth conditions. In this study, we show that deletion of the genes encoding either NfnI or NfnII affects the cellular concentrations of NAD(P)H and particularly NADPH. This results in a moderate to severe growth phenotype in deletion mutants, demonstrating a key role for each enzyme in maintaining redox homeostasis. Despite their similarity in primary sequence and cofactor content, crystallographic, kinetic, and mass spectrometry analyses reveal that there are fundamental structural differences between the two enzymes, and NfnII does not catalyze the NfnI bifurcating reaction. Instead, it exhibits non-bifurcating ferredoxin NADP oxidoreductase-type activity. NfnII is therefore proposed to be a bifunctional enzyme and also to catalyze a bifurcating reaction, although its third substrate, in addition to ferredoxin and NADP(H), is as yet unknown.

Refactoring the Embden-Meyerhof-Parnas Pathway as a Whole of Portable GlucoBricks for Implantation of Glycolytic Modules in Gram-Negative Bacteria

The Embden–Meyerhof–Parnas (EMP) pathway is generally considered to be the biochemical standard for glucose catabolism. Alas, its native genomic organization and the control of gene expression in Escherichia coli are both very intricate, which limits the portability of the EMP pathway to other biotechnologically important bacterial hosts that lack the route. In this work, the genes encoding all the enzymes of the linear EMP route have been individually recruited from the genome of E. coli K-12, edited in silico to remove their endogenous regulatory signals, and synthesized de novo following a standard (GlucoBrick) that enables their grouping in the form of functional modules at the user’s will. After verifying their activity in several glycolytic mutants of E. coli, the versatility of these GlucoBricks was demonstrated in quantitative physiology tests and biochemical assays carried out in Pseudomonas putida KT2440 and P. aeruginosa PAO1 as the heterologous hosts. Specific configurations of GlucoBricks were also adopted to streamline the downward circulation of carbon from hexoses to pyruvate in E. coli recombinants, thereby resulting in a 3-fold increase of poly(3-hydroxybutyrate) synthesis from glucose. Refactoring whole metabolic blocks in the fashion described in this work thus eases the engineering of biochemical processes where the optimization of carbon traffic is facilitated by the operation of the EMP pathway—which yields more ATP than other glycolytic routes such as the Entner–Doudoroff pathway.

Pcal_0111, a highly thermostable bifunctional fructose-1,6-bisphosphate aldolase/phosphatase from Pyrobaculum calidifontis

Pyrobaculum calidifontis genome harbors an open reading frame Pcal_0111 annotated as fructose bisphosphate aldolase. Although the gene is annotated as fructose bisphosphate aldolase, it exhibits a high homology with previously reported fructose-1,6-bisphosphate aldolase/phosphatase from Thermoproteus neutrophilus. To examine the biochemical properties of Pcal_0111, we have cloned and expressed the gene in Escherichia coli. Purified recombinant Pcal_0111 catalyzed both phosphatase and aldolase reactions with specific activity values of 4 U and 1.3 U, respectively. These values are highest among the fructose 1,6-bisphosphatases/aldolases characterized from archaea. The enzyme activity increased linearly with the increase in temperature until 100 °C. Recombinant Pcal_0111 is highly stable with a half-life of 120 min at 100 °C. There was no significant change in the circular dichroism spectra of the protein up to 90 °C. The enzyme activity was not affected by AMP but strongly inhibited by ATP with an IC₅₀ value of 0.75 mM and mildly by ADP. High thermostability and inhibition by ATP make Pcal_0111 a unique fructose 1,6-bisphosphatase/aldolase.

An archaeal ADP-dependent serine kinase involved in cysteine biosynthesis and serine metabolism

Routes for cysteine biosynthesis are still unknown in many archaea. Here we find that the hyperthermophilic archaeon Thermococcus kodakarensis generates cysteine from serine via O-phosphoserine, in addition to the classical route from 3-phosphoglycerate. The protein responsible for serine phosphorylation is encoded by TK0378, annotated as a chromosome partitioning protein ParB. The TK0378 protein utilizes ADP as the phosphate donor, but in contrast to previously reported ADP-dependent kinases, recognizes a non-sugar substrate. Activity is specific towards free serine, and not observed with threonine, homoserine and serine residues within a peptide. Genetic analyses suggest that TK0378 is involved in serine assimilation and clearly responsible for cysteine biosynthesis from serine. TK0378 homologs, present in Thermococcales and Desulfurococcales, are most likely not ParB proteins and constitute a group of kinases involved in serine utilization.

Archaea metabolism has unique adaptations to hostile environments. Here Makino et al. describe an unusual ADP-dependent kinase that phosphorylates free serine to O-phosphoserine and participates in an additional cysteine biosynthetic pathway in the archaeon Thermococcus kodakarensis.

Draft genome sequence of sulfur-reducing archaeon Thermococcus thioreducens DSM 14981T

Thermococcus thioreducens DSM 14981^T, a sulfur-reducing archaeon, was isolated from the rainbow hydrothermal vent site on the Mid-Atlantic Ridge. Herein, we report the draft genome sequence of T. thioreducens DSM 14981^T; we obtained 41 contigs with a genome size of 2,052,483 bp and G + C content of 53.5%. This genome sequence will not only help understand how the archaeon adapts to the deep-sea hydrothermal environment but also aid the development of enzymes that are highly stable under extreme conditions for industrial applications.

Metabolic processes of Methanococcus maripaludis and potential applications

Methanococcus maripaludis is a rapidly growing, fully sequenced, genetically tractable model organism among hydrogenotrophic methanogens. It has the ability to convert CO2 and H2 into a useful cleaner energy fuel (CH4). In fact, this conversion enhances in the presence of free nitrogen as the sole nitrogen source due to prolonged cell growth. Given the global importance of GHG emissions and climate change, diazotrophy can be attractive for carbon capture and utilization applications from appropriately treated flue gases, where surplus hydrogen is available from renewable electricity sources. In addition, M. maripaludis can be engineered to produce other useful products such as terpenoids, hydrogen, methanol, etc. M. maripaludis with its unique abilities has the potential to be a workhorse like Escherichia coli and S. cerevisiae for fundamental and experimental biotechnology studies. More than 100 experimental studies have explored different specific aspects of the biochemistry and genetics of CO2 and N2 fixation by M. maripaludis. Its genome-scale metabolic model (iMM518) also exists to study genetic perturbations and complex biological interactions. However, a comprehensive review describing its cell structure, metabolic processes, and methanogenesis is still lacking in the literature. This review fills this crucial gap. Specifically, it integrates distributed information from the literature to provide a complete and detailed view for metabolic processes such as acetyl-CoA synthesis, pyruvate synthesis, glycolysis/gluconeogenesis, reductive tricarboxylic acid (RTCA) cycle, non-oxidative pentose phosphate pathway (NOPPP), nitrogen metabolism, amino acid metabolism, and nucleotide biosynthesis. It discusses energy production via methanogenesis and its relation to metabolism. Furthermore, it reviews taxonomy, cell structure, culture/storage conditions, molecular biology tools, genome-scale models, and potential industrial and environmental applications. Through the discussion, it develops new insights and hypotheses from experimental and modeling observations, and identifies opportunities for further research and applications.

The nutritional status of Methanosarcina acetivorans regulates glycogen metabolism and gluconeogenesis and glycolysis fluxes

Gluconeogenesis is an essential pathway in methanogens because they are unable to use exogenous hexoses as carbon source for cell growth. With the aim of understanding the regulatory mechanisms of central carbon metabolism in Methanosarcina acetivorans, the present study investigated gene expression, the activities and metabolic regulation of key enzymes, metabolite contents and fluxes of gluconeogenesis, as well as glycolysis and glycogen synthesis/degradation pathways. Cells were grown with methanol as a carbon source. Key enzymes were kinetically characterized at physiological pH/temperature. Active consumption of methanol during exponential cell growth correlated with significant methanogenesis, gluconeogenic flux and steady glycogen synthesis. After methanol exhaustion, cells reached the stationary growth phase, which correlated with the rise in glycogen consumption and glycolytic flux, decreased methanogenesis, negligible acetate production and an absence of gluconeogenesis. Elevated activities of carbon monoxide dehydrogenase/acetyl-CoA synthetase complex and pyruvate: ferredoxin oxidoreductase suggested the generation of acetyl-CoA and pyruvate for glycogen synthesis. In the early stationary growth phase, the transcript contents and activities of pyruvate phosphate dikinase, fructose 1,6-bisphosphatase and glycogen synthase decreased, whereas those of glycogen phosphorylase, ADP-phosphofructokinase and pyruvate kinase increased. Therefore, glycogen and gluconeogenic metabolites were synthesized when an external carbon source was provided. Once such a carbon source became depleted, glycolysis and methanogenesis fed by glycogen degradation provided the ATP supply. Weak inhibition of key enzymes by metabolites suggested that the pathways evaluated were mainly transcriptionally regulated. Because glycogen metabolism and glycolysis/gluconeogenesis are not present in all methanogens, the overall data suggest that glycogen storage might represent an environmental advantage for methanosarcinales when carbon sources are scarce. Also, the understanding of the central carbohydrate metabolism in methanosarcinales may help to optimize methane production.

Hydrogen Production and Enzyme Activities in the Hyperthermophile Thermococcus paralvinellae Grown on Maltose, Tryptone, and Agricultural Waste

Thermococcus may be an important alternative source of H2 in the hot subseafloor in otherwise low H2 environments such as some hydrothermal vents and oil reservoirs. It may also be useful in industry for rapid agricultural waste treatment and concomitant H2 production. Thermococcus paralvinellae grown at 82°C without sulfur produced up to 5 mmol of H2 L(-1) at rates of 5-36 fmol H2 cell(-1) h(-1) on 0.5% (wt vol(-1)) maltose, 0.5% (wt vol(-1)) tryptone, and 0.5% maltose + 0.05% tryptone media. Two potentially inhibiting conditions, the presence of 10 mM acetate and low pH (pH 5) in maltose-only medium, did not significantly affect growth or H2 production. Growth rates, H2 production rates, and cell yields based on H2 production were the same as those for Pyrococcus furiosus grown at 95°C on the same media for comparison. Acetate, butyrate, succinate, isovalerate, and formate were also detected as end products. After 100 h, T. paralvinellae produced up to 5 mmol of H2 L(-1) of medium when grown on up to 70% (vol vol(-1)) waste milk from cows undergoing treatment for mastitis with the bacterial antibiotic Ceftiofur and from untreated cows. The amount of H2 produced by T. paralvinellae increased with increasing waste concentrations, but decreased in P. furiosus cultures supplemented with waste milk above 1% concentration. All mesophilic bacteria from the waste milk that grew on Luria Bertani, Sheep's Blood (selective for Staphylococcus, the typical cause of mastitis), and MacConkey (selective for Gram-negative enteric bacteria) agar plates were killed by heat during incubation at 82°C. Ceftiofur, which is heat labile, was below the detection limit following incubation at 82°C. T. paralvinellae also produced up to 6 mmol of H2 L(-1) of medium when grown on 0.1-10% (wt vol(-1)) spent brewery grain while P. furiosus produced < 1 mmol of H2 L(-1). Twelve of 13 enzyme activities in T. paralvinellae showed significant (p < 0.05) differences across six different growth conditions; however, methyl viologen-dependent membrane hydrogenase activity remained constant across all media types. The results demonstrate the potential of at least some Thermococcus species to produce H2 if protein and α-glucosides are present as substrates.

Preserved hydride transfer mechanism in evolutionarily divergent thymidylate synthases

Thymidylate synthase (TSase) catalyzes a hydride transfer in the last step of the de novo biosynthesis of the DNA nucleotide thymine. We compared two isozymes, namely, TSase from Escherichia coli (ecTSase) and TSase from Bacillus subtilis (bsTSase) that represent a case of divergent evolution. Interestingly, a highly conserved histidine (H147 of ecTSase) was proposed to serve a critical role in catalysis, but in bsTSase it is naturally substituted by valine (Val). Yet, bsTSase is more active than ecTSase, and the intrinsic kinetic isotope effects (KIEs) of both are temperature-independent, suggesting a similarly well-organized transition state (TS) for the catalyzed hydride transfer. To examine the role of that histidine (His) in TSase catalysis, we examined the kinetics of H147V ecTSase, which "bridges" between these two TSases. In contrast to both wild-type TSases, the single mutation results in deficient catalysis. The mutation leads to intrinsic KIEs that are temperature-dependent, indicating a substantial imperfection in its TS. The findings reveal two important features: a direct role of H147 in the hydride transfer step catalyzed by the ecTSase and the evolutionary compensation for its deficiency in bsTSase via extensive polymorphism across the protein. Very different active site residues are observed for these evolutionarily divergent isozymes, which result in a well-organized TS for both. It is suggested that evolutionary pressure compensated for the H to V substitution at the active site of bsTSase by polymorphism leading to a well-organized TS in both enzymes.

Draft genome sequence of Thermococcus sp. EP1, a novel hyperthermophilic archaeon isolated from a deep-sea hydrothermal vent on the East Pacific Rise

Thermococcus sp. strain EP1 is a novel anaerobic hyperthermophilic archaeon isolated from a deep-sea hydrothermal vent on the East Pacific Rise. It grows optimally at 80 °C and can produce industrial enzymes at high temperature. We report here the draft genome of EP1, which contains 1,819,157 bp with a G+C content of 39.3%. The sequence will provide the genetic basis for better understanding of adaptation to hydrothermal environment and the development of novel thermostable enzymes for industrial application.