The non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN) of Sulfolobus solfataricus: a key-enzyme of the semi-phosphorylative branch of the Entner-Doudoroff pathway

Graphene oxide stimulated low-temperature denitrification activity of microbial communities in lake sediments by enhancing anabolism and inhibiting cellular respiration

Nitrate pollution in fresh water is becoming increasingly serious. In this study, the effects of temperature and graphene oxide materials on the potential functions of denitrification communities in lake sediments were investigated by metagenome. The addition of graphene oxide significantly affected the abundance of denitrification genes such as Nap, Nos, and enhanced the contribution of Pseudomonas, making low temperature and material addition conducive to the denitrification process. Module network implied that low temperature increased the centrality of denitrification in community functions. At low temperatures, graphene oxide enhanced community anabolism by stimulation organic carbon consumption and regulating the gene abundance in the citric acid cycle and the semi-phosphorylation Entner-Doudoroff, thus possibly stimulating extracellular polymeric substances (EPS) synthesis and secretion. In addition, graphene oxide may also regulate the transfer of reducing electrons from NADH to denitrifying enzymes by affecting the gene abundances of complex I and complex IV.

The Non-phosphorylating Glyceraldehyde-3-Phosphate Dehydrogenase GapN Is a Potential New Drug Target in Streptococcus pyogenes

The strict human pathogen Streptococcus pyogenes causes infections of varying severity, ranging from self-limiting suppurative infections to life-threatening diseases like necrotizing fasciitis or streptococcal toxic shock syndrome. Here, we show that the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase GapN is an essential enzyme for S. pyogenes. GapN converts glyceraldehyde 3-phosphate into 3-phosphoglycerate coupled to the reduction of NADP to NADPH. The knock-down of gapN by antisense peptide nucleic acids (asPNA) significantly reduces viable bacterial counts of S. pyogenes laboratory and macrolide-resistant clinical strains in vitro. As S. pyogenes lacks the oxidative part of the pentose phosphate pathway, GapN appears to be the major NADPH source for the bacterium. Accordingly, other streptococci that carry a complete pentose phosphate pathway are not prone to asPNA-based gapN knock-down. Determination of the crystal structure of the S. pyogenes GapN apo-enzyme revealed an unusual cis-peptide in proximity to the catalytic binding site. Furthermore, using a structural modeling approach, we correctly predicted competitive inhibition of S. pyogenes GapN by erythrose 4-phosphate, indicating that our structural model can be used for in silico screening of specific GapN inhibitors. In conclusion, the data provided here reveal that GapN is a potential target for antimicrobial substances that selectively kill S. pyogenes and other streptococci that lack the oxidative part of the pentose phosphate pathway.

Study of ALDH from Thermus thermophilus-Expression, Purification and Characterisation of the Non-Substrate Specific, Thermophilic Enzyme Displaying Both Dehydrogenase and Esterase Activity

Aldehyde dehydrogenases (ALDH), found in all kingdoms of life, form a superfamily of enzymes that primarily catalyse the oxidation of aldehydes to form carboxylic acid products, while utilising the cofactor NAD(P)⁺. Some superfamily members can also act as esterases using p-nitrophenyl esters as substrates. The ALDH_Tt from Thermus thermophilus was recombinantly expressed in E. coli and purified to obtain high yields (approximately 15–20 mg/L) and purity utilising an efficient heat treatment step coupled with IMAC and gel filtration chromatography. The use of the heat treatment step proved critical, in its absence decreased yield of 40% was observed. Characterisation of the thermophilic ALDH_Tt led to optimum enzymatic working conditions of 50 °C, and a pH of 8. ALDH_Tt possesses dual enzymatic activity, with the ability to act as a dehydrogenase and an esterase. ALDH_Tt possesses broad substrate specificity, displaying activity for a range of aldehydes, most notably hexanal and the synthetic dialdehyde, terephthalaldehyde. Interestingly, para-substituted benzaldehydes could be processed efficiently, but ortho-substitution resulted in no catalytic activity. Similarly, ALDH_Tt displayed activity for two different esterase substrates, p-nitrophenyl acetate and p-nitrophenyl butyrate, but with activities of 22.9% and 8.9%, respectively, compared to the activity towards hexanal.

The biology of thermoacidophilic archaea from the order Sulfolobales

Thermoacidophilic archaea belonging to the order Sulfolobales thrive in extreme biotopes, such as sulfuric hot springs and ore deposits. These microorganisms have been model systems for understanding life in extreme environments, as well as for probing the evolution of both molecular genetic processes and central metabolic pathways. Thermoacidophiles, such as the Sulfolobales, use typical microbial responses to persist in hot acid (e.g. motility, stress response, biofilm formation), albeit with some unusual twists. They also exhibit unique physiological features, including iron and sulfur chemolithoautotrophy, that differentiate them from much of the microbial world. Although first discovered >50 years ago, it was not until recently that genome sequence data and facile genetic tools have been developed for species in the Sulfolobales. These advances have not only opened up ways to further probe novel features of these microbes but also paved the way for their potential biotechnological applications. Discussed here are the nuances of the thermoacidophilic lifestyle of the Sulfolobales, including their evolutionary placement, cell biology, survival strategies, genetic tools, metabolic processes and physiological attributes together with how these characteristics make thermoacidophiles ideal platforms for specialized industrial processes.

Large-scale computational discovery and analysis of virus-derived microbial nanocompartments

Encapsulins are a class of microbial protein compartments defined by the viral HK97-fold of their capsid protein, self-assembly into icosahedral shells, and dedicated cargo loading mechanism for sequestering specific enzymes. Encapsulins are often misannotated and traditional sequence-based searches yield many false positive hits in the form of phage capsids. Here, we develop an integrated search strategy to carry out a large-scale computational analysis of prokaryotic genomes with the goal of discovering an exhaustive and curated set of all HK97-fold encapsulin-like systems. We find over 6,000 encapsulin-like systems in 31 bacterial and four archaeal phyla, including two novel encapsulin families. We formulate hypotheses about their potential biological functions and biomedical relevance, which range from natural product biosynthesis and stress resistance to carbon metabolism and anaerobic hydrogen production. An evolutionary analysis of encapsulins and related HK97-type virus families shows that they share a common ancestor, and we conclude that encapsulins likely evolved from HK97-type bacteriophages.

Utilization of Phenol as Carbon Source by the Thermoacidophilic Archaeon Saccharolobus solfataricus P2 Is Limited by Oxygen Supply and the Cellular Stress Response

Present in many industrial effluents and as common degradation product of organic matter, phenol is a widespread compound which may cause serious environmental problems, due to its toxicity to animals and humans. Degradation of phenol from the environment by mesophilic bacteria has been studied extensively over the past decades, but only little is known about phenol biodegradation at high temperatures or low pH. In this work we studied phenol degradation in the thermoacidophilic archaeon Saccharolobus solfataricus P2 (basonym: Sulfolobus solfataricus) under extreme conditions (80°C, pH 3.5). We combined metabolomics and transcriptomics together with metabolic modeling to elucidate the organism’s response to growth with phenol as sole carbon source. Although S. solfataricus is able to utilize phenol for biomass production, the carbon source induces profound stress reactions, including genome rearrangement as well as a strong intracellular accumulation of polyamines. Furthermore, computational modeling revealed a 40% higher oxygen demand for substrate oxidation, compared to growth on glucose. However, only 16.5% of oxygen is used for oxidation of phenol to catechol, resulting in a less efficient integration of carbon into the biomass. Finally, our data underlines the importance of the phenol meta-degradation pathway in S. solfataricus and enables us to predict enzyme candidates involved in the degradation processes downstream of 2-hydroxymucconic acid.

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.

The genomes of polyextremophilic cyanidiales contain 1% horizontally transferred genes with diverse adaptive functions

The role and extent of horizontal gene transfer (HGT) in eukaryotes are hotly disputed topics that impact our understanding of the origin of metabolic processes and the role of organelles in cellular evolution. We addressed this issue by analyzing 10 novel Cyanidiales genomes and determined that 1% of their gene inventory is HGT-derived. Numerous HGT candidates share a close phylogenetic relationship with prokaryotes that live in similar habitats as the Cyanidiales and encode functions related to polyextremophily. HGT candidates differ from native genes in GC-content, number of splice sites, and gene expression. HGT candidates are more prone to loss, which may explain the absence of a eukaryotic pan-genome. Therefore, the lack of a pan-genome and cumulative effects fail to provide substantive arguments against our hypothesis of recurring HGT followed by differential loss in eukaryotes. The maintenance of 1% HGTs, even under selection for genome reduction, underlines the importance of non-endosymbiosis related foreign gene acquisition.

Correlating enzyme annotations with a large set of microbial growth temperatures reveals metabolic adaptations to growth at diverse temperatures

Background

The ambient temperature of all habitats is a key physical property that shapes the biology of microbes inhabiting them. The optimal growth temperature (OGT) of a microbe, is therefore a key piece of data needed to understand evolutionary adaptations manifested in their genome sequence. Unfortunately there is no growth temperature database or easily downloadable dataset encompassing the majority of cultured microorganisms. We are thus limited in interpreting genomic data to identify temperature adaptations in microbes.

Results

In this work I significantly contribute to closing this gap by mining data from major culture collection centres to obtain growth temperature data for a nonredundant set of 21,498 microbes. The dataset (10.5281/zenodo.1175608) contains mainly bacteria and archaea and spans psychrophiles, mesophiles, thermophiles and hyperthermophiles. Using this data a full 43% of all protein entries in the UniProt database can be annotated with the growth temperature of the species from which they originate. I validate the dataset by showing a Pearson correlation of up to 0.89 between growth temperature and mean enzyme optima, a physiological property directly influenced by the growth temperature. Using the temperature dataset I correlate the genomic occurance of enzyme functional annotations with growth temperature. I identify 319 enzyme functions that either increase or decrease in occurrence with temperature. Eight metabolic pathways were statistically enriched for these enzyme functions. Furthermore, I establish a correlation between 33 domains of unknown function (DUFs) with growth temperature in microbes, four of which (DUF438, DUF1524, DUF1957 and DUF3458_C) were significant in both archaea and bacteria.

Conclusions

The growth temperature dataset enables large-scale correlation analysis with enzyme function- and domain-level annotations. Growth-temperature dependent changes in their occurrence highlight potential evolutionary adaptations. A few of the identified changes are previously known, such as the preference for menaquinone biosynthesis through the futalosine pathway in bacteria growing at high temperatures. Others represent important starting points for future studies, such as DUFs where their occurrence change with temperature. The growth temperature dataset should become a valuable community resource and will find additional, important, uses in correlating genomic, transcriptomic, proteomic, metabolomic, phenotypic or taxonomic properties with temperature in future studies.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1320-7) contains supplementary material, which is available to authorized users.

Profiling of glucose-induced transcription in Sulfolobus acidocaldarius DSM 639

Sulfolobus species can grow on a variety of organic compounds as carbon and energy sources. These species degrade glucose to pyruvate by the modified branched Entner-Doudoroff pathway. We attempted to determine the differentially expressed genes (DEGs) under sugar-limited and sugar-rich conditions. RNA sequencing (RNA-seq) was used to quantify the expression of the genes and identify those DEGs between the S. acidocaldarius cells grown under sugar-rich (YT with glucose) and sugar-limited (YT only) conditions. The functions and pathways of the DEGs were examined using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Quantitative real-time PCR (qRT-PCR) was performed to validate the DEGs. Transcriptome analysis of the DSM 639 strain grown on sugar-limited and sugar-rich media revealed that 853 genes were differentially expressed, among which 481 were upregulated and 372 were downregulated under the glucose-supplemented condition. In particular, 70 genes showed significant changes in expression levels of ≥ twofold. GO and KEGG enrichment analyses revealed that the genes encoding components of central carbon metabolism, the respiratory chain, and protein and amino acid biosynthetic machinery were upregulated under the glucose condition. RNA-seq and qRT-PCR analyses indicated that the sulfur assimilation genes (Saci_2197-2204) including phosphoadenosine phosphosulfate reductase and sulfite reductase were significantly upregulated in the presence of glucose. The present study revealed metabolic networks in S. acidocaldarius that are induced in a glucose-dependent manner, improving our understanding of biomass production under sugar-rich conditions.

The quaternary structure of Thermus thermophilus aldehyde dehydrogenase is stabilized by an evolutionary distinct C-terminal arm extension

Aldehyde dehydrogenases (ALDH) form a superfamily of dimeric or tetrameric enzymes that catalyze the oxidation of a broad range of aldehydes into their corresponding carboxylic acids with the concomitant reduction of the cofactor NAD(P) into NAD(P)H. Despite their varied polypeptide chain length and oligomerisation states, ALDHs possess a conserved architecture of three domains: the catalytic domain, NAD(P)⁺ binding domain, and the oligomerization domain. Here, we describe the structure and function of the ALDH from Thermus thermophilus (ALDH_Tt) which exhibits non-canonical features of both dimeric and tetrameric ALDH and a previously uncharacterized C-terminal arm extension forming novel interactions with the N-terminus in the quaternary structure. This unusual tail also interacts closely with the substrate entry tunnel in each monomer providing further mechanistic detail for the recent discovery of tail-mediated activity regulation in ALDH. However, due to the novel distal extension of the tail of ALDH_Tt and stabilizing termini-interactions, the current model of tail-mediated substrate access is not apparent in ALDH_Tt. The discovery of such a long tail in a deeply and early branching phylum such as Deinococcus-Thermus indicates that ALDH_Tt may be an ancestral or primordial metabolic model of study. This structure provides invaluable evidence of how metabolic regulation has evolved and provides a link to early enzyme regulatory adaptations.

A Phosphofructokinase Homolog from Pyrobaculum calidifontis Displays Kinase Activity towards Pyrimidine Nucleosides and Ribose 1-Phosphate

The genome of the hyperthermophilic archaeon Pyrobaculum calidifontis contains an open reading frame, Pcal_0041, annotated as encoding a PfkB family ribokinase, consisting of phosphofructokinase and pyrimidine kinase domains. Among the biochemically characterized enzymes, the Pcal_0041 protein was 37% identical to the phosphofructokinase (Ape_0012) from Aeropyrum pernix, which displayed kinase activity toward a broad spectrum of substrates, including sugars, sugar phosphates, and nucleosides, and 36% identical to a phosphofructokinase from Desulfurococcus amylolyticus To examine the biochemical function of the Pcal_0041 protein, we cloned and expressed the gene and purified the recombinant protein. Although the Pcal_0041 protein contained a putative phosphofructokinase domain, it exhibited only low levels of phosphofructokinase activity. The recombinant enzyme catalyzed the phosphorylation of nucleosides and, to a lower extent, sugars and sugar phosphates. Surprisingly, among the substrates tested, the highest activity was detected with ribose 1-phosphate (R1P), followed by cytidine and uridine. The catalytic efficiency (k _cat/K_m ) toward R1P was 11.5 mM^-1 · s^-1 ATP was the most preferred phosphate donor, followed by GTP. Activity measurements with cell extracts of P. calidifontis indicated the presence of nucleoside phosphorylase activity, which would provide the means to generate R1P from nucleosides. The study suggests that, in addition to the recently identified ADP-dependent ribose 1-phosphate kinase (R1P kinase) in Thermococcus kodakarensis that functions in the pentose bisphosphate pathway, R1P kinase is also present in members of the Crenarchaeota.IMPORTANCE The discovery of the pentose bisphosphate pathway in Thermococcus kodakarensis has clarified how this archaeon can degrade nucleosides. Homologs of the enzymes of this pathway are present in many members of the Thermococcales, suggesting that this metabolism occurs in these organisms. However, this is not the case in other archaea, and degradation mechanisms for nucleosides or ribose 1-phosphate are still unknown. This study reveals an important first step in understanding nucleoside metabolism in Crenarchaeota and identifies an ATP-dependent ribose 1-phosphate kinase in Pyrobaculum calidifontis The enzyme is structurally distinct from previously characterized archaeal members of the ribokinase family and represents a group of proteins found in many crenarchaea.

Dihydrodipicolinate synthase is absent in fungi

The class I aldolase dihydrodipicolinate synthase (DHDPS) catalyzes the first committed step of the diaminopimelate (DAP) lysine biosynthesis pathway in bacteria, archaea and plants. Despite the existence, in databases, of numerous fungal sequences annotated as DHDPS, its presence in fungi has been the subject of contradictory claims. We report the characterization of DHDPS candidates from fungi. Firstly, the putative DHDPS from Coccidioides immitis (PDB ID: 3QFE) was shown to have negligible enzyme activity. Sequence analysis of 3QFE showed that three out of the seven amino acid residues critical for DHDPS activity are absent; however, exact matches to catalytic residues from two other class I aldolases, 2-keto-3-deoxygluconate aldolase (KDGA), and 4-hydroxy-2-oxoglutarate aldolase (HOGA), were identified. The presence of both KDGA and HOGA activity in 3QFE was confirmed in vitro using enzyme assays, the first report of such dual activity. Subsequent analyses of all publically available fungal sequences revealed that no entry contains all seven residues important for DHDPS function. The candidate with the highest number of identities (6 of 7), KIW77228 from Fonsecaea pedrosoi, was shown to have trace DHDPS activity in vitro, partially restored by substitution of the seventh critical residue, and to be incapable of complementing DHDPS-deficient E. coli cells. Combined with the presence of all seven sequences for the alternative α-aminoadipate (AAA) lysine biosynthesis pathway in C. immitis and F. pedrosoi, we believe that DHDPS and the DAP pathway are absent in fungi, and further, that robust informed methods for annotating genes need to be implemented.

Molecular basis of autotrophic vs mixotrophic growth in Chlorella sorokiniana

In this work, we investigated the molecular basis of autotrophic vs. mixotrophic growth of Chlorella sorokiniana, one of the most productive microalgae species with high potential to produce biofuels, food and high value compounds. To increase biomass accumulation, photosynthetic microalgae are commonly cultivated in mixotrophic conditions, adding reduced carbon sources to the growth media. In the case of C. sorokiniana, the presence of acetate enhanced biomass, proteins, lipids and starch productivity when compared to autotrophic conditions. Despite decreased chlorophyll content, photosynthetic properties were essentially unaffected while differential gene expression profile revealed transcriptional regulation of several genes mainly involved in control of carbon flux. Interestingly, acetate assimilation caused upregulation of phosphoenolpyruvate carboxylase enzyme, enabling potential recovery of carbon atoms lost by acetate oxidation. The obtained results allowed to associate the increased productivity observed in mixotrophy in C. sorokiniana with a different gene regulation leading to a fine regulation of cell metabolism.

Two distinct glyceraldehyde-3-phosphate dehydrogenases in glycolysis and gluconeogenesis in the archaeon Haloferax volcanii

The halophilic archaeon Haloferax volcanii degrades glucose via the semiphosphorylative Entner-Doudoroff pathway and can also grow on gluconeogenic substrates. Here, the enzymes catalysing the conversion of glyceraldehyde-3-phosphate (GAP) to 3-phosphoglycerate were analysed. The genome contains the genes gapI and gapII encoding two putative GAP dehydrogenases, and pgk encoding phosphoglycerate kinase (PGK). We show that gapI is functionally involved in sugar catabolism, whereas gapII is involved in gluconeogenesis. For pgk, an amphibolic function is indicated. This is the first report of the functional involvement of a phosphorylating glyceraldehyde-3-phosphate dehydrogenase and PGK in sugar catabolism in archaea. Phylogenetic analyses indicate that the catabolic gapI from H. volcanii is acquired from bacteria via lateral genetransfer, whereas the anabolic gapII as well as pgk are of archaeal origin.

Phosphoglycerate kinase acts as a futile cycle at high temperature

In (hyper)thermophilic organisms metabolic processes have to be adapted to function optimally at high temperature. We compared the gluconeogenic conversion of 3-phosphoglycerate via 1,3-bisphosphoglycerate to glyceraldehyde-3-phosphate at 30 °C and at 70 °C. At 30 °C it was possible to produce 1,3-bisphosphoglycerate from 3-phosphoglycerate with phosphoglycerate kinase, but at 70 °C, 1,3-bisphosphoglycerate was dephosphorylated rapidly to 3-phosphoglycerate, effectively turning the phosphoglycerate kinase into a futile cycle. When phosphoglycerate kinase was incubated together with glyceraldehyde 3-phosphate dehydrogenase it was possible to convert 3-phosphoglycerate to glyceraldehyde 3-phosphate, both at 30 °C and at 70 °C, however, at 70 °C only low concentrations of product were observed due to thermal instability of glyceraldehyde 3-phosphate. Thus, thermolabile intermediates challenge central metabolic reactions and require special adaptation strategies for life at high temperature.

Systems biology of the modified branched Entner-Doudoroff pathway in Sulfolobus solfataricus

Sulfolobus solfataricus is a thermoacidophilic Archaeon that thrives in terrestrial hot springs (solfatares) with optimal growth at 80°C and pH 2–4. It catabolizes specific carbon sources, such as D-glucose, to pyruvate via the modified Entner-Doudoroff (ED) pathway. This pathway has two parallel branches, the semi-phosphorylative and the non-phosphorylative. However, the strategy of S.solfataricus to endure in such an extreme environment in terms of robustness and adaptation is not yet completely understood. Here, we present the first dynamic mathematical model of the ED pathway parameterized with quantitative experimental data. These data consist of enzyme activities of the branched pathway at 70°C and 80°C and of metabolomics data at the same temperatures for the wild type and for a metabolic engineered knockout of the semi-phosphorylative branch. We use the validated model to address two questions: 1. Is this system more robust to perturbations at its optimal growth temperature? 2. Is the ED robust to deletion and perturbations? We employed a systems biology approach to answer these questions and to gain further knowledge on the emergent properties of this biological system. Specifically, we applied deterministic and stochastic approaches to study the sensitivity and robustness of the system, respectively. The mathematical model we present here, shows that: 1. Steady state metabolite concentrations of the ED pathway are consistently more robust to stochastic internal perturbations at 80°C than at 70°C; 2. These metabolite concentrations are highly robust when faced with the knockout of either branch. Connected with this observation, these two branches show different properties at the level of metabolite production and flux control. These new results reveal how enzyme kinetics and metabolomics synergizes with mathematical modelling to unveil new systemic properties of the ED pathway in S.solfataricus in terms of its adaptation and robustness.

The Peculiar Glycolytic Pathway in Hyperthermophylic Archaea: Understanding Its Whims by Experimentation In Silico

Mathematical models are key to systems biology where they typically describe the topology and dynamics of biological networks, listing biochemical entities and their relationships with one another. Some (hyper)thermophilic Archaea contain an enzyme, called non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN), which catalyzes the direct oxidation of glyceraldehyde-3-phosphate to 3-phosphoglycerate omitting adenosine 5′-triphosphate (ATP) formation by substrate-level-phosphorylation via phosphoglycerate kinase. In this study we formulate three hypotheses that could explain functionally why GAPN exists in these Archaea, and then construct and use mathematical models to test these three hypotheses. We used kinetic parameters of enzymes of Sulfolobus solfataricus (S. solfataricus) which is a thermo-acidophilic archaeon that grows optimally between 60 and 90 °C and between pH 2 and 4. For comparison, we used a model of Saccharomyces cerevisiae (S. cerevisiae), an organism that can live at moderate temperatures. We find that both the first hypothesis, i.e., that the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) plus phosphoglycerate kinase (PGK) route (the alternative to GAPN) is thermodynamically too much uphill and the third hypothesis, i.e., that GAPDH plus PGK are required to carry the flux in the gluconeogenic direction, are correct. The second hypothesis, i.e., that the GAPDH plus PGK route delivers less than the 1 ATP per pyruvate that is delivered by the GAPN route, is only correct when GAPDH reaction has a high rate and 1,3-bis-phosphoglycerate (BPG) spontaneously degrades to 3PG at a high rate.

Genome Sequencing of Sulfolobus sp. A20 from Costa Rica and Comparative Analyses of the Putative Pathways of Carbon, Nitrogen, and Sulfur Metabolism in Various Sulfolobus Strains

The genome of Sulfolobus sp. A20 isolated from a hot spring in Costa Rica was sequenced. This circular genome of the strain is 2,688,317 bp in size and 34.8% in G+C content, and contains 2591 open reading frames (ORFs). Strain A20 shares ~95.6% identity at the 16S rRNA gene sequence level and <30% DNA-DNA hybridization (DDH) values with the most closely related known Sulfolobus species (i.e., Sulfolobus islandicus and Sulfolobus solfataricus), suggesting that it represents a novel Sulfolobus species. Comparison of the genome of strain A20 with those of the type strains of S. solfataricus, Sulfolobus acidocaldarius, S. islandicus, and Sulfolobus tokodaii, which were isolated from geographically separated areas, identified 1801 genes conserved among all Sulfolobus species analyzed (core genes). Comparative genome analyses show that central carbon metabolism in Sulfolobus is highly conserved, and enzymes involved in the Entner-Doudoroff pathway, the tricarboxylic acid cycle and the CO₂ fixation pathways are predominantly encoded by the core genes. All Sulfolobus species encode genes required for the conversion of ammonium into glutamate/glutamine. Some Sulfolobus strains have gained the ability to utilize additional nitrogen source such as nitrate (i.e., S. islandicus strain REY15A, LAL14/1, M14.25, and M16.27) or urea (i.e., S. islandicus HEV10/4, S. tokodaii strain7, and S. metallicus DSM 6482). The strategies for sulfur metabolism are most diverse and least understood. S. tokodaii encodes sulfur oxygenase/reductase (SOR), whereas both S. islandicus and S. solfataricus contain genes for sulfur reductase (SRE). However, neither SOR nor SRE genes exist in the genome of strain A20, raising the possibility that an unknown pathway for the utilization of elemental sulfur may be present in the strain. The ability of Sulfolobus to utilize nitrate or sulfur is encoded by a gene cluster flanked by IS elements or their remnants. These clusters appear to have become fixed at a specific genomic site in some strains and lost in other strains during the course of evolution. The versatility in nitrogen and sulfur metabolism may represent adaptation of Sulfolobus to thriving in different habitats.

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.

Stability of Ensemble Models Predicts Productivity of Enzymatic Systems

Stability in a metabolic system may not be obtained if incorrect amounts of enzymes are used. Without stability, some metabolites may accumulate or deplete leading to the irreversible loss of the desired operating point. Even if initial enzyme amounts achieve a stable steady state, changes in enzyme amount due to stochastic variations or environmental changes may move the system to the unstable region and lose the steady-state or quasi-steady-state flux. This situation is distinct from the phenomenon characterized by typical sensitivity analysis, which focuses on the smooth change before loss of stability. Here we show that metabolic networks differ significantly in their intrinsic ability to attain stability due to the network structure and kinetic forms, and that after achieving stability, some enzymes are prone to cause instability upon changes in enzyme amounts. We use Ensemble Modelling for Robustness Analysis (EMRA) to analyze stability in four cell-free enzymatic systems when enzyme amounts are changed. Loss of stability in continuous systems can lead to lower production even when the system is tested experimentally in batch experiments. The predictions of instability by EMRA are supported by the lower productivity in batch experimental tests. The EMRA method incorporates properties of network structure, including stoichiometry and kinetic form, but does not require specific parameter values of the enzymes.

NADPH-generating systems in bacteria and archaea

Reduced nicotinamide adenine dinucleotide phosphate (NADPH) is an essential electron donor in all organisms. It provides the reducing power that drives numerous anabolic reactions, including those responsible for the biosynthesis of all major cell components and many products in biotechnology. The efficient synthesis of many of these products, however, is limited by the rate of NADPH regeneration. Hence, a thorough understanding of the reactions involved in the generation of NADPH is required to increase its turnover through rational strain improvement. Traditionally, the main engineering targets for increasing NADPH availability have included the dehydrogenase reactions of the oxidative pentose phosphate pathway and the isocitrate dehydrogenase step of the tricarboxylic acid (TCA) cycle. However, the importance of alternative NADPH-generating reactions has recently become evident. In the current review, the major canonical and non-canonical reactions involved in the production and regeneration of NADPH in prokaryotes are described, and their key enzymes are discussed. In addition, an overview of how different enzymes have been applied to increase NADPH availability and thereby enhance productivity is provided.

Residues that influence coenzyme preference in the aldehyde dehydrogenases

To find out the residues that influence the coenzyme preference of aldehyde dehydrogenases (ALDHs), we reviewed, analyzed and correlated data from their known crystal structures and amino-acid sequences with their published kinetic parameters for NAD(P)(+). We found that the conformation of the Rossmann-fold loops participating in binding the adenosine ribose is very conserved among ALDHs, so that coenzyme specificity is mainly determined by the nature of the residue at position 195 (human ALDH2 numbering). Enzymes with glutamate or proline at 195 prefer NAD(+) because the side-chains of these residues electrostatically and/or sterically repel the 2'-phosphate group of NADP(+). But contrary to the conformational rigidity of proline, the conformational flexibility of glutamate may allow NADP(+)-binding in some enzymes by moving the carboxyl group away from the 2'-phosphate group, which is possible if a small neutral residue is located at position 224, and favored if the residue at position 53 interacts with Glu195 in a NADP(+)-compatible conformation. Of the residues found at position 195, only glutamate interacts with the NAD(+)-adenosine ribose; glutamine and histidine cannot since their side-chain points are opposite to the ribose, probably because the absence of the electrostatic attraction by the conserved nearby Lys192, or its electrostatic repulsion, respectively. The shorter side-chains of other residues-aspartate, serine, threonine, alanine, valine, leucine, or isoleucine-are distant from the ribose but leave room for binding the 2'-phosphate group. Generally, enzymes having a residue different from Glu bind NAD(+) with less affinity, but they can also bind NADP(+) even sometimes with higher affinity than NAD(+), as do enzymes containing Thr/Ser/Gln195. Coenzyme preference is a variable feature within many ALDH families, consistent with being mainly dependent on a single residue that apparently has no other structural or functional roles, and therefore can easily be changed through evolution and selected in response to physiological needs.

Tracking difference in gene expression in a time-course experiment using gene set enrichment analysis

Fistulifera sp. strain JPCC DA0580 is a newly sequenced pennate diatom that is capable of simultaneously growing and accumulating lipids. This is a unique trait, not found in other related microalgae so far. It is able to accumulate between 40 to 60% of its cell weight in lipids, making it a strong candidate for the production of biofuel. To investigate this characteristic, we used RNA-Seq data gathered at four different times while Fistulifera sp. strain JPCC DA0580 was grown in oil accumulating and non-oil accumulating conditions. We then adapted gene set enrichment analysis (GSEA) to investigate the relationship between the difference in gene expression of 7,822 genes and metabolic functions in our data. We utilized information in the KEGG pathway database to create the gene sets and changed GSEA to use re-sampling so that data from the different time points could be included in the analysis. Our GSEA method identified photosynthesis, lipid synthesis and amino acid synthesis related pathways as processes that play a significant role in oil production and growth in Fistulifera sp. strain JPCC DA0580. In addition to GSEA, we visualized the results by creating a network of compounds and reactions, and plotted the expression data on top of the network. This made existing graph algorithms available to us which we then used to calculate a path that metabolizes glucose into triacylglycerol (TAG) in the smallest number of steps. By visualizing the data this way, we observed a separate up-regulation of genes at different times instead of a concerted response. We also identified two metabolic paths that used less reactions than the one shown in KEGG and showed that the reactions were up-regulated during the experiment. The combination of analysis and visualization methods successfully analyzed time-course data, identified important metabolic pathways and provided new hypotheses for further research.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Improvement of NADPH bioavailability in Escherichia coli by replacing NAD(+)-dependent glyceraldehyde-3-phosphate dehydrogenase GapA with NADP (+)-dependent GapB from Bacillus subtilis and addition of NAD kinase

Enzymatic synthesis of some industrially important compounds depends heavily on cofactor NADPH as the reducing agent. This is especially true in the synthesis of chiral compounds that are often used as pharmaceutical intermediates to generate the correct stereochemistry in bioactive products. The high cost and technical difficulty of cofactor regeneration often pose a challenge for such biocatalytic reactions. In this study, to increase NADPH bioavailability, the native NAD(+)-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gapA gene in Escherichia coli was replaced with a NADP(+)-dependent gapB from Bacillus subtilis. To overcome the limitation of NADP(+) availability, E. coli NAD kinase, nadK was also coexpressed with gapB. The recombinant strains were then tested in three reporting systems: biosynthesis of lycopene, oxidation of cyclohexanone with cyclohexanone monooxygenase (CHMO), and an anaerobic system utilizing 2-haloacrylate reductase (CAA43). In all the reporting systems, replacing NAD(+)-dependent GapA activity with NADP(+)-dependent GapB activity increased the synthesis of NADPH-dependent compounds. The increase was more pronounced when NAD kinase was also overexpressed in the case of the one-step reaction catalyzed by CAA43 which approximately doubled the product yield. These results validate this novel approach to improve NADPH bioavailability in E. coli and suggest that the strategy can be applied in E. coli or other bacterium-based production of NADPH-dependent compounds.

Intermediate instability at high temperature leads to low pathway efficiency for an in vitro reconstituted system of gluconeogenesis in Sulfolobus solfataricus

Four enzymes of the gluconeogenic pathway in Sulfolobus solfataricus were purified and kinetically characterized. The enzymes were reconstituted in vitro to quantify the contribution of temperature instability of the pathway intermediates to carbon loss from the system. The reconstituted system, consisting of phosphoglycerate kinase, glyceraldehyde 3-phosphate dehydrogenase, triose phosphate isomerase and the fructose 1,6-bisphosphate aldolase/phosphatase, maintained a constant consumption rate of 3-phosphoglycerate and production of fructose 6-phosphate over a 1-h period. Cofactors ATP and NADPH were regenerated via pyruvate kinase and glucose dehydrogenase. A mathematical model was constructed on the basis of the kinetics of the purified enzymes and the measured half-life times of the pathway intermediates. The model quantitatively predicted the system fluxes and metabolite concentrations. Relative enzyme concentrations were chosen such that half the carbon in the system was lost due to degradation of the thermolabile intermediates dihydroxyacetone phosphate, glyceraldehyde 3-phosphate and 1,3-bisphosphoglycerate, indicating that intermediate instability at high temperature can significantly affect pathway efficiency.

Archaeal aldehyde dehydrogenase ST0064 from Sulfolobus tokodaii, a paralog of non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase, is a succinate semialdehyde dehydrogenase

Aldehyde dehydrogenase ST0064, the closest paralog of previously characterized allosteric non-phosphorylating glyceraldehyde-3-phosphate (GAP) dehydrogenase (GAPN, ST2477) from a thermoacidophilic archaeon, Sulfolobus tokodaii, was expressed heterologously and characterized in detail. ST0064 showed remarkable activity toward succinate semialdehyde (SSA) (Km of 0.0029 mM and kcat of 30.0 s(-1)) with no allosteric regulation. Activity toward GAP was lower (Km of 4.6 mM and kcat of 4.77 s(-1)), and previously predicted succinyl-CoA reductase activity was not detected, suggesting that the enzyme functions practically as succinate semialdehyde dehydrogenase (SSADH). Phylogenetic analysis indicated that archaeal SSADHs and GAPNs are closely related within the aldehyde dehydrogenase superfamily, suggesting that they are of the same origin.

An unusual effect of NADP+ on the thermostability of the nonphosphorylating glyceraldehyde-3-phosphate dehydrogenase from Streptococcus mutans

Adiabatic differential scanning calorimetry was used to investigate the effect of NADP+ on the irreversible thermal denaturation of the nonphosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN) from Streptococcus mutans. The GAPN-NADP+ binary complex showed a strongly decreased thermal stability, with a difference of about 20 °C between the temperatures of the thermal transition peak maxima of the complex and the free protein. This finding was similar to the previously described thermal destabilization of GAPN upon binding of inorganic phosphate to the substrate binding site and can be interpreted as the shift of the equilibrium between 2 conformers of tetrameric GAPN upon addition of the coenzyme. Single amino acid substitution, known to abolish the NADP+ binding, cancelled the calorimetric effect of the coenzyme. GAPN thermal inactivation was considerably decelerated in the presence of NADP+ showing that the apparent change in stability of the active centre can be the opposite to that of the whole protein molecule. NADP+ could also reactivate the inactive GAPN* species, obtained by the heating of the apoenzyme below the thermal denaturation transition temperature. These effects may reflect a mechanism that provides GAPN the sufficient flexibility for the earlier observed profound active site reorganizations required during the catalytic cycle. The elevated thermal stability of the apoenzyme may, in turn, be important for maintaining a constant level of active GAPN--an enzyme that is known to be crucial for the effective supply of the reducing equivalents in S. mutans and its ability to grow under aerobic conditions.

Unraveling the function of paralogs of the aldehyde dehydrogenase super family from Sulfolobus solfataricus

Aldehyde dehydrogenases (ALDHs) have been well established in all three domains of life and were shown to play essential roles, e.g., in intermediary metabolism and detoxification. In the genome of Sulfolobus solfataricus, five paralogs of the aldehyde dehydrogenases superfamily were identified, however, so far only the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN) and α-ketoglutaric semialdehyde dehydrogenase (α-KGSADH) have been characterized. Detailed biochemical analyses of the remaining three ALDHs revealed the presence of two succinic semialdehyde dehydrogenase (SSADH) isoenzymes catalyzing the NAD(P)(+)-dependent oxidation of succinic semialdehyde. Whereas SSO1629 (SSADH-I) is specific for NAD(+), SSO1842 (SSADH-II) exhibits dual cosubstrate specificity (NAD(P)(+)). Physiological significant activity for both SSO-SSADHs was only detected with succinic semialdehyde and α-ketoglutarate semialdehyde. Bioinformatic reconstructions suggest a major function of both enzymes in γ-aminobutyrate, polyamine as well as nitrogen metabolism and they might additionally also function in pentose metabolism. Phylogenetic studies indicated a close relationship of SSO-SSALDHs to GAPNs and also a convergent evolution with the SSADHs from E. coli. Furthermore, for SSO1218, methylmalonate semialdehyde dehydrogenase (MSDH) activity was demonstrated. The enzyme catalyzes the NAD(+)- and CoA-dependent oxidation of methylmalonate semialdehyde, malonate semialdehyde as well as propionaldehyde (PA). For MSDH, a major function in the degradation of branched chain amino acids is proposed which is supported by the high sequence homology with characterized MSDHs from bacteria. This is the first report of MSDH as well as SSADH isoenzymes in Archaea.

Coping with complexity in metabolic engineering

In the past decade, systems biology has revealed great metabolic and regulatory complexity even in seemingly simple microbial systems. Metabolic engineering aims to control this complexity in order to establish sustainable and economically viable production routes for valuable chemicals. Recent advances in systems-level data generation and modeling of cellular metabolism and regulation together with tremendous progress in synthetic biology will provide the tools to put biotechnologists on the fast track for implementing novel production processes. Great potential lies in the reduction of cellular complexity by orthogonalization of metabolic modules. Here, we review recent advances that will eventually enable metabolic engineers to predict, design, and build streamlined microbial cell factories with reduced time and effort.

Mechanism of oxygen detoxification by the surprisingly oxygen-tolerant hyperthermophilic archaeon, Pyrococcus furiosus

The anaerobic archaeon Pyrococcus furiosus grows by fermenting carbohydrates producing H(2), CO(2), and acetate. We show here that it is surprisingly tolerant to oxygen, growing well in the presence of 8% (vol/vol) O(2). Although cell growth and acetate production were not significantly affected by O(2), H(2) production was reduced by 50% (using 8% O(2)). The amount of H(2) produced decreased in a linear manner with increasing concentrations of O(2) over the range 2-12% (vol/vol), and for each mole of O(2) consumed, the amount of H(2) produced decreased by approximately 2 mol. The recycling of H(2) by the two cytoplasmic hydrogenases appeared not to play a role in O(2) resistance because a mutant strain lacking both enzymes was not more sensitive to O(2) than the parent strain. Decreased H(2) production was also not due to inactivation of the H(2)-producing, ferredoxin-dependent membrane-bound hydrogenase because its activity was unaffected by O(2) exposure. Electrons from carbohydrate oxidation must therefore be diverted to relieve O(2) stress at the level of reduced ferredoxin before H(2) production. Deletion strains lacking superoxide reductase (SOR) and putative flavodiiron protein A showed increased sensitivity to O(2), indicating that these enzymes play primary roles in resisting O(2). However, a mutant strain lacking the proposed electron donor to SOR, rubredoxin, was unaffected in response to O(2). Hence, electrons from sugar oxidation normally used to produce H(2) are diverted to O(2) detoxification by SOR and putative flavodiiron protein A, but the electron flow pathway from ferredoxin does not necessarily involve rubredoxin.|

Change of carbon source causes dramatic effects in the phospho-proteome of the archaeon Sulfolobus solfataricus

Protein phosphorylation is known to occur in Archaea. However, knowledge of phosphorylation in the third domain of life is rather scarce. Homology-based searches of archaeal genome sequences reveals the absence of two-component systems in crenarchaeal genomes but the presence of eukaryotic-like protein kinases and protein phosphatases. Here, the influence of the offered carbon source (glucose versus tryptone) on the phospho-proteome of Sulfolobus solfataricus P2 was studied by precursor acquisition independent from ion count (PAcIFIC). In comparison to previous phospho-proteome studies, a high number of phosphorylation sites (1318) located on 690 phospho-peptides from 540 unique phospho-proteins were detected, thus increasing the number of currently known archaeal phospho-proteins from 80 to 621. Furthermore, a 25.8/20.6/53.6 Ser/Thr/Tyr percentage ratio with an unexpectedly high predominance of tyrosine phosphorylation was detected. Phospho-proteins in most functional classes (21 out of 26 arCOGs) were identified, suggesting an important regulatory role in S. solfataricus. Focusing on the central carbohydrate metabolism in response to the offered carbon source, significant changes were observed. The observed complex phosphorylation pattern hints at an important physiological function of protein phosphorylation in control of the central carbohydrate metabolism, which might particularly operate in channeling carbon flux into the respective metabolic pathways.

Synthetic metabolic engineering-a novel, simple technology for designing a chimeric metabolic pathway

BACKGROUND

The integration of biotechnology into chemical manufacturing has been recognized as a key technology to build a sustainable society. However, the practical applications of biocatalytic chemical conversions are often restricted due to their complexities involving the unpredictability of product yield and the troublesome controls in fermentation processes. One of the possible strategies to overcome these limitations is to eliminate the use of living microorganisms and to use only enzymes involved in the metabolic pathway. Use of recombinant mesophiles producing thermophilic enzymes at high temperature results in denaturation of indigenous proteins and elimination of undesired side reactions; consequently, highly selective and stable biocatalytic modules can be readily prepared. By rationally combining those modules together, artificial synthetic pathways specialized for chemical manufacturing could be designed and constructed.

RESULTS

A chimeric Embden-Meyerhof (EM) pathway with balanced consumption and regeneration of ATP and ADP was constructed by using nine recombinant E. coli strains overproducing either one of the seven glycolytic enzymes of Thermus thermophilus, the cofactor-independent phosphoglycerate mutase of Pyrococcus horikoshii, or the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase of Thermococcus kodakarensis. By coupling this pathway with the Thermus malate/lactate dehydrogenase, a stoichiometric amount of lactate was produced from glucose with an overall ATP turnover number of 31.

CONCLUSIONS

In this study, a novel and simple technology for flexible design of a bespoke metabolic pathway was developed. The concept has been testified via a non-ATP-forming chimeric EM pathway. We designated this technology as "synthetic metabolic engineering". Our technology is, in principle, applicable to all thermophilic enzymes as long as they can be functionally expressed in the host, and thus would be potentially applicable to the biocatalytic manufacture of any chemicals or materials on demand.

Comparative analysis of two glyceraldehyde-3-phosphate dehydrogenases from a thermoacidophilic archaeon, Sulfolobus tokodaii

Sulfolobus tokodaii, a thermoacidophilic archaeon, possesses two structurally and functionally different enzymes that catalyze the oxidation of glyceraldehyde-3-phosphate (GAP): non-phosphorylating GAP dehydrogenase (St-GAPN) and phosphorylating GAP dehydrogenase (St-GAPDH). In contrast to previously characterized GAPN from Sulfolobus solfataricus, which exhibits V-type allosterism, St-GAPN showed K-type allosterism in which the positive cooperativity was abolished with concomitant activation by glucose 1-phosphate (G1P). St-GAPDH catalyzed the reversible oxidation of GAP to 1,3-bisphosphoglycerate (1,3-BPG) with high gluconeogenic activity, which was specific for NADPH, while both NAD(+) and NADP(+) were utilized in the glycolytic direction.

Genome sequencing of a genetically tractable Pyrococcus furiosus strain reveals a highly dynamic genome

The model archaeon Pyrococcus furiosus grows optimally near 100°C on carbohydrates and peptides. Its genome sequence (NCBI) was determined 12 years ago. A genetically tractable strain, COM1, was very recently reported, and here we describe its genome sequence. Of 1,909,827 bp in size, it is 1,571 bp longer (0.1%) than the reference NCBI sequence. The COM1 genome contains numerous chromosomal rearrangements, deletions, and single base changes. COM1 also has 45 full or partial insertion sequences (ISs) compared to 35 in the reference NCBI strain, and these have resulted in the direct deletion or insertional inactivation of 13 genes. Another seven genes were affected by chromosomal deletions and are predicted to be nonfunctional. In addition, the amino acid sequences of another 102 of the 2,134 predicted gene products are different in COM1. These changes potentially impact various cellular functions, including carbohydrate, peptide, and nucleotide metabolism; DNA repair; CRISPR-associated defense; transcriptional regulation; membrane transport; and growth at 72°C. For example, the IS-mediated inactivation of riboflavin synthase in COM1 resulted in a riboflavin requirement for growth. Nevertheless, COM1 grew on cellobiose, malto-oligosaccharides, and peptides in complex and minimal media at 98 and 72°C to the same extent as did both its parent strain and a new culture collection strain (DSMZ 3638). This was in spite of COM1 lacking several metabolic enzymes, including nonphosphorylating glyceraldehyde-3-phosphate dehydrogenase and beta-glucosidase. The P. furiosus genome is therefore of high plasticity, and the availability of the COM1 sequence will be critical for the future studies of this model hyperthermophile.

Biochemical and genetic characterization of the three metabolic routes in Thermococcus kodakarensis linking glyceraldehyde 3-phosphate and 3-phosphoglycerate

In the classical Embden-Meyerhof (EM) pathway for glycolysis, the conversion between glyceraldehyde 3-phosphate (GAP) and 3-phosphoglycerate (3-PGA) is reversibly catalysed by phosphorylating GAP dehydrogenase (GAPDH) and phosphoglycerate kinase (PGK). In the Euryarchaeota Thermococcus kodakarensis and Pyrococcus furiosus, an additional gene encoding GAP:ferredoxin oxidoreductase (GAPOR) and a gene similar to non-phosphorylating GAP dehydrogenase (GAPN) are present. In order to determine the physiological roles of the three routes that link GAP and 3-PGA, we individually disrupted the GAPOR, GAPN, GAPDH and PGK genes (gor, gapN, gapDH and pgk respectively) of T. kodakarensis. The Δgor strain displayed no growth under glycolytic conditions, confirming its proposed function to generate reduced ferredoxin for energy generation in glycolysis. Surprisingly, ΔgapN cells also did not grow under glycolytic conditions, suggesting that GAPN plays a key role in providing NADPH under these conditions. Disruption of gor and gapN had no effect on gluconeogenic growth. Growth experiments with the ΔgapDH and Δpgk strains indicated that, unlike their counterparts in the classical EM pathway, GAPDH/PGK play a major role only in gluconeogenesis. Biochemical analyses indicated that T. kodakarensis GAPN did not recognize aldehyde substrates other than d-GAP, preferred NADP(+) as cofactor and was dramatically activated with glucose 1-phosphate.

Novel metabolic pathways in Archaea

The Archaea harbor many metabolic pathways that differ to previously recognized classical pathways. Glycolysis is carried out by modified versions of the Embden-Meyerhof and Entner-Doudoroff pathways. Thermophilic archaea have recently been found to harbor a bi-functional fructose-1,6-bisphosphate aldolase/phosphatase for gluconeogenesis. A number of novel pentose-degrading pathways have also been recently identified. In terms of anabolic metabolism, a pathway for acetate assimilation, the methylaspartate cycle, and two CO2-fixing pathways, the 3-hydroxypropionate/4-hydroxybutyrate cycle and the dicarboxylate/4-hydroxybutyrate cycle, have been elucidated. As for biosynthetic pathways, recent studies have clarified the enzymes responsible for several steps involved in the biosynthesis of inositol phospholipids, polyamine, coenzyme A, flavin adeninedinucleotide and heme. By examining the presence/absence of homologs of these enzymes on genome sequences, we have found that the majority of these enzymes and pathways are specific to the Archaea.

"Hot standards" for the thermoacidophilic archaeon Sulfolobus solfataricus

Within the archaea, the thermoacidophilic crenarchaeote Sulfolobus solfataricus has become an important model organism for physiology and biochemistry, comparative and functional genomics, as well as, more recently also for systems biology approaches. Within the Sulfolobus Systems Biology (“SulfoSYS”)-project the effect of changing growth temperatures on a metabolic network is investigated at the systems level by integrating genomic, transcriptomic, proteomic, metabolomic and enzymatic information for production of a silicon cell-model. The network under investigation is the central carbohydrate metabolism. The generation of high-quality quantitative data, which is critical for the investigation of biological systems and the successful integration of the different datasets, derived for example from high-throughput approaches (e.g., transcriptome or proteome analyses), requires the application and compliance of uniform standard protocols, e.g., for growth and handling of the organism as well as the “–omics” approaches. Here, we report on the establishment and implementation of standard operating procedures for the different wet-lab and in silico techniques that are applied within the SulfoSYS-project and that we believe can be useful for future projects on Sulfolobus or (hyper)thermophiles in general. Beside established techniques, it includes new methodologies like strain surveillance, the improved identification of membrane proteins and the application of crenarchaeal metabolomics.

SulfoSYS (Sulfolobus Systems Biology): towards a silicon cell model for the central carbohydrate metabolism of the archaeon Sulfolobus solfataricus under temperature variation

SulfoSYS (Sulfolobus Systems Biology) focuses on the study of the CCM (central carbohydrate metabolism) of Sulfolobus solfataricus and its regulation under temperature variation at the systems level. In Archaea, carbohydrates are metabolized by modifications of the classical pathways known from Bacteria or Eukarya, e.g. the unusual branched ED (Entner–Doudoroff) pathway, which is utilized for glucose degradation in S. solfataricus. This archaeal model organism of choice is a thermoacidophilic crenarchaeon that optimally grows at 80°C (60–92°C) and pH 2–4. In general, life at high temperature requires very efficient adaptation to temperature changes, which is most difficult to deal with for organisms, and it is unclear how biological networks can withstand and respond to such changes. This integrative project combines genomic, transcriptomic, proteomic and metabolomic, as well as kinetic and biochemical information. The final goal of SulfoSYS is the construction of a silicon cell model for this part of the living cell that will enable computation of the CCM network. In the present paper, we report on one of the first archaeal systems biology projects.

DNA microarray analysis of central carbohydrate metabolism: glycolytic/gluconeogenic carbon switch in the hyperthermophilic crenarchaeum Thermoproteus tenax

In order to unravel the role of regulation on transcript level in central carbohydrate metabolism (CCM) of Thermoproteus tenax, a focused DNA microarray was constructed by using 85 open reading frames involved in CCM. A transcriptional analysis comparing heterotrophic growth on glucose versus autotrophic growth on CO2-H2 was performed.

Glycerate kinase of the hyperthermophilic archaeon Thermoproteus tenax: new insights into the phylogenetic distribution and physiological role of members of the three different glycerate kinase classes

Background

The presence of the branched Entner-Doudoroff (ED) pathway in two hyperthermophilic Crenarchaea, the anaerobe Thermoproteus tenax and the aerobe Sulfolobus solfataricus, was suggested. However, so far no enzymatic information of the non-phosphorylative ED branch and especially its key enzyme – glycerate kinase – was available. In the T. tenax genome, a gene homolog with similarity to putative hydroxypyruvate reductase/glycerate dehydrogenase and glycerate kinase was identified.

Results

The encoding gene was expressed in E. coli in a recombinant form, the gene product purified and the glycerate kinase activity was confirmed by enzymatic studies. The enzyme was active as a monomer and catalyzed the ATP-dependent phosphorylation of D-glycerate forming exclusively 2-phosphoglycerate. The enzyme was specific for glycerate and highest activity was observed with ATP as phosphoryl donor and Mg²⁺ as divalent cation. ATP could be partially replaced by GTP, CTP, TTP and UTP. The enzyme showed high affinity for D-glycerate (K_m 0.02 ± 0.01 mM, V_max of 5.05 ± 0.52 U/mg protein) as well as ATP (K_m of 0.03 ± 0.01 mM, V_max of 4.41 ± 0.04 U/mg protein), although at higher glycerate concentrations, substrate inhibition was observed. Furthermore, the enzyme was inhibited by its product ADP via competitive inhibition. Data bank searches revealed that archaeal glycerate kinases are members of the MOFRL (multi-organism fragment with rich leucine) family, and homologs are found in all three domains of life.

Conclusion

A re-evaluation of available genome sequence information as well as biochemical and phylogenetic studies revealed the presence of the branched ED pathway as common route for sugar degradation in Archaea that utilize the ED pathway. Detailed analyses including phylogenetic studies demonstrate the presence of three distinct glycerate kinase classes in extant organisms that share no common origin. The affiliation of characterized glycerate kinases with the different enzyme classes as well as their physiological/cellular function reveals no association with particular pathways but a separate phylogenetic distribution. This work highlights the diversity and complexity of the central carbohydrate metabolism. The data also support a key function of the conversion of glycerate to 2- or 3-phosphoglycerate via glycerate kinase in funneling various substrates into the common EMP pathway for catabolic and anabolic purposes.