NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase from Thermoproteus tenax. The first identified archaeal member of the aldehyde dehydrogenase superfamily is a glycolytic enzyme with unusual regulatory properties

Adaptive laboratory evolution recruits the promiscuity of succinate semialdehyde dehydrogenase to repair different metabolic deficiencies

Promiscuous enzymes often serve as the starting point for the evolution of novel functions. Yet, the extent to which the promiscuity of an individual enzyme can be harnessed several times independently for different purposes during evolution is poorly reported. Here, we present a case study illustrating how NAD(P)+-dependent succinate semialdehyde dehydrogenase of Escherichia coli (Sad) is independently recruited through various evolutionary mechanisms for distinct metabolic demands, in particular vitamin biosynthesis and central carbon metabolism. Using adaptive laboratory evolution (ALE), we show that Sad can substitute for the roles of erythrose 4-phosphate dehydrogenase in pyridoxal 5'-phosphate (PLP) biosynthesis and glyceraldehyde 3-phosphate dehydrogenase in glycolysis. To recruit Sad for PLP biosynthesis and glycolysis, ALE employs various mechanisms, including active site mutation, copy number amplification, and (de)regulation of gene expression. Our study traces down these different evolutionary trajectories, reports on the surprising active site plasticity of Sad, identifies regulatory links in amino acid metabolism, and highlights the potential of an ordinary enzyme as innovation reservoir for evolution.

A vesicular Warburg effect: Aerobic glycolysis occurs on axonal vesicles for local NAD+ recycling and transport

In neurons, fast axonal transport (FAT) of vesicles occurs over long distances and requires constant and local energy supply for molecular motors in the form of adenosine triphosphate (ATP). FAT is independent of mitochondrial metabolism. Indeed, the glycolytic machinery is present on vesicles and locally produces ATP, as well as nicotinamide adenine dinucleotide bonded with hydrogen (NADH) and pyruvate, using glucose as a substrate. It remains unclear whether pyruvate is transferred to mitochondria from the vesicles as well as how NADH is recycled into NAD+ on vesicles for continuous glycolysis activity. The optimization of a glycolytic activity test for subcellular compartments allowed the evaluation of the kinetics of vesicular glycolysis in the brain. This revealed that glycolysis is more efficient on vesicles than in the cytosol. We also found that lactate dehydrogenase (LDH) enzymatic activity is required for effective vesicular ATP production. Indeed, inhibition of LDH or the forced degradation of pyruvate inhibited ATP production from axonal vesicles. We found LDHA rather than the B isoform to be enriched on axonal vesicles suggesting a preferential transformation of pyruvate to lactate and a concomitant recycling of NADH into NAD+ on vesicles. Finally, we found that LDHA inhibition dramatically reduces the FAT of both dense-core vesicles and synaptic vesicle precursors in a reconstituted cortico-striatal circuit on-a-chip. Together, this shows that aerobic glycolysis is required to supply energy for vesicular transport in neurons, similar to the Warburg effect.

Encapsulins

Subcellular compartmentalization is a defining feature of all cells. In prokaryotes, compartmentalization is generally achieved via protein-based strategies. The two main classes of microbial protein compartments are bacterial microcompartments and encapsulin nanocompartments. Encapsulins self-assemble into proteinaceous shells with diameters between 24 and 42 nm and are defined by the viral HK97-fold of their shell protein. Encapsulins have the ability to encapsulate dedicated cargo proteins, including ferritin-like proteins, peroxidases, and desulfurases. Encapsulation is mediated by targeting sequences present in all cargo proteins. Encapsulins are found in many bacterial and archaeal phyla and have been suggested to play roles in iron storage, stress resistance, sulfur metabolism, and natural product biosynthesis. Phylogenetic analyses indicate that they share a common ancestor with viral capsid proteins. Many pathogens encode encapsulins, and recent evidence suggests that they may contribute toward pathogenicity. The existing information on encapsulin structure, biochemistry, biological function, and biomedical relevance is reviewed here.

Structural and biochemical evidence that ATP inhibits the cancer biomarker human aldehyde dehydrogenase 1A3

Human aldehyde dehydrogenase (ALDH) participates in the oxidative stress response and retinoid metabolism, being involved in several diseases, including cancer, diabetes and obesity. The ALDH1A3 isoform has recently elicited wide interest because of its potential use as a cancer stem cell biomarker and drug target. We report high-resolution three-dimensional ALDH1A3 structures for the apo-enzyme, the NAD+ complex and a binary complex with ATP. Each subunit of the ALDH1A3-ATP complex contains one ATP molecule bound to the adenosine-binding pocket of the cofactor-binding site. The ATP complex also shows a molecule, putatively identified as a polyethylene glycol aldehyde, covalently bound to the active-site cysteine. This mimics the thioacyl-enzyme catalytic intermediate, which is trapped in a dead enzyme lacking an active cofactor. At physiological concentrations, ATP inhibits the dehydrogenase activity of ALDH1A3 and other isoforms, with a Ki value of 0.48 mM for ALDH1A3, showing a mixed inhibition type against NAD+. ATP also inhibits esterase activity in a concentration-dependent manner. The current ALDH1A3 structures at higher resolution will facilitate the rational design of potent and selective inhibitors. ATP binding to ALDH1A3 enables activity modulation by the energy status of the cell and metabolic reprogramming, which may be relevant in several disease conditions.

The glycerol-3-phosphate dehydrogenases GpsA and GlpD constitute the oxidoreductive metabolic linchpin for Lyme disease spirochete host infectivity and persistence in the tick

We have identified GpsA, a predicted glycerol-3-phosphate dehydrogenase, as a virulence factor in the Lyme disease spirochete Borrelia (Borreliella) burgdorferi: GpsA is essential for murine infection and crucial for persistence of the spirochete in the tick. B. burgdorferi has a limited biosynthetic and metabolic capacity; the linchpin connecting central carbohydrate and lipid metabolism is at the interconversion of glycerol-3-phosphate and dihydroxyacetone phosphate, catalyzed by GpsA and another glycerol-3-phosphate dehydrogenase, GlpD. Using a broad metabolomics approach, we found that GpsA serves as a dominant regulator of NADH and glycerol-3-phosphate levels in vitro, metabolic intermediates that reflect the cellular redox potential and serve as a precursor for lipid and lipoprotein biosynthesis, respectively. Additionally, GpsA was required for survival under nutrient stress, regulated overall reductase activity and controlled B. burgdorferi morphology in vitro. Furthermore, during in vitro nutrient stress, both glycerol and N-acetylglucosamine were bactericidal to B. burgdorferi in a GlpD-dependent manner. This study is also the first to identify a suppressor mutation in B. burgdorferi: a glpD deletion restored the wild-type phenotype to the pleiotropic gpsA mutant, including murine infectivity by needle inoculation at high doses, survival under nutrient stress, morphological changes and the metabolic imbalance of NADH and glycerol-3-phosphate. These results illustrate how basic metabolic functions that are dispensable for in vitro growth can be essential for in vivo infectivity of B. burgdorferi and may serve as attractive therapeutic targets.

Lyme disease (borreliosis) is the most common tick-borne disease in the Northern hemisphere and its prevalence is increasing. Borrelia burgdorferi, the etiological agent of Lyme disease, is an enzootic pathogen that alternates between a tick vector and vertebrate host. Humans are considered an incidental host after transmission of B. burgdorferi following the bite of an infected tick. The mechanisms by which B. burgdorferi persists in the Ixodid tick, transmits to a vertebrate host and establishes infection are not well understood. Therefore, identifying virulence factors and uncovering the pathogenic strategies in the spirochete remain important to address the public health concerns of Lyme disease. In this study, we identify an enzyme involved in three-carbon metabolism, GpsA, as a new virulence factor with an effect on persistence in ticks. GpsA and GlpD, another enzyme, constitute a bidirectional metabolic node connecting lipid biosynthesis and glycolysis, which serves as the linchpin for regulating carbon utilization for B. burgdorferi throughout its enzootic cycle. Disruption of this node causes a lethal metabolic imbalance revealing a potential therapeutic target for the treatment of Lyme disease.

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.

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.

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.

Insights into stress responses in mandarins triggered by Bacillus subtilis cyclic lipopeptides and exogenous plant hormones upon Penicillium digitatum infection

Bacillus subtilis CLP extract activates defense gene expression and increases the unique protein production involving in pathways of ISR, SAR, ubiquitin-proteasome system, and glycolysis for stress responses in flavedo tissues. Cyclic lipopeptides (CLPs) of Bacillus subtilis ABS-S14 had ability to activate plant defensive pathways, increase resistance and control green mold rot caused by Penicillium digitatum in mandarin fruit. The current study investigated transcriptional and proteomic data to highlight the unique induction effect of CLPs produced by B. subtilis ABS-S14 on the defense mechanism of mandarins in response to P. digitatum attack, and their differences from those following the exogenous plant hormone application. The proteomic patterns of the flavedo tissues as affected by Bacillus CLP extract, salicylic acid (SA), methyl jasmonate (MeJA), and ethephon (Et) were explored. qPCR analysis revealed the great effects of CLP extract in enhancing the transcription of PAL, ACS1, GLU, POD, and PR1. Tryptic peptides by LC-MS analysis between treatments with and without fungal infection were compared. B. subtilis CLP extract empowered the plant's immune response to wound stress by the significant production of calmodulin-binding receptor-like cytoplasmic kinase 2, molybdenum cofactor sulfurase, and NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase. Ubiquitin carrier protein abundance was developed only in the treated flavedo with CLP extract coupled with P. digitatum infection. The gene expression and overall proteome findings involving pathways of ubiquitin proteasome system, ISR, SAR, and energy production provide a new insight into the molecular mechanisms of the antagonist B. subtilis ABS-S14 inducing resistance against green mold in mandarins.

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.

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.

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.

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.

Bacterial Anaerobic Synthesis Gas (Syngas) and CO2+H2 Fermentation

Anaerobic bacterial gas fermentation gains broad interest in various scientific, social, and industrial fields. This microbial process is carried out by a specific group of bacterial strains called acetogens. All these strains employ the Wood-Ljungdahl pathway but they belong to different taxonomic groups. Here we provide an overview of the metabolism of acetogens and naturally occurring products. Characteristics of 61 strains were summarized and selected acetogens described in detail. Acetobacterium woodii, Clostridium ljungdahlii, and Moorella thermoacetica serve as model organisms. Results of approaches such as genome-scale modeling, proteomics, and transcriptomics are discussed. Metabolic engineering of acetogens can be used to expand the product portfolio to platform chemicals and to study different aspects of cell physiology. Moreover, the fermentation of gases requires specific reactor configurations and the development of the respective technology, which can be used for an industrial application. Even though the overall process will have a positive effect on climate, since waste and greenhouse gases could be converted into commodity chemicals, some legislative barriers exist, which hamper successful exploitation of this technology.

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.

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.

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.

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.

Exploiting protein flexibility to predict the location of allosteric sites

BACKGROUND

Allostery is one of the most powerful and common ways of regulation of protein activity. However, for most allosteric proteins identified to date the mechanistic details of allosteric modulation are not yet well understood. Uncovering common mechanistic patterns underlying allostery would allow not only a better academic understanding of the phenomena, but it would also streamline the design of novel therapeutic solutions. This relatively unexplored therapeutic potential and the putative advantages of allosteric drugs over classical active-site inhibitors fuel the attention allosteric-drug research is receiving at present. A first step to harness the regulatory potential and versatility of allosteric sites, in the context of drug-discovery and design, would be to detect or predict their presence and location. In this article, we describe a simple computational approach, based on the effect allosteric ligands exert on protein flexibility upon binding, to predict the existence and position of allosteric sites on a given protein structure.

RESULTS

By querying the literature and a recently available database of allosteric sites, we gathered 213 allosteric proteins with structural information that we further filtered into a non-redundant set of 91 proteins. We performed normal-mode analysis and observed significant changes in protein flexibility upon allosteric-ligand binding in 70% of the cases. These results agree with the current view that allosteric mechanisms are in many cases governed by changes in protein dynamics caused by ligand binding. Furthermore, we implemented an approach that achieves 65% positive predictive value in identifying allosteric sites within the set of predicted cavities of a protein (stricter parameters set, 0.22 sensitivity), by combining the current analysis on dynamics with previous results on structural conservation of allosteric sites. We also analyzed four biological examples in detail, revealing that this simple coarse-grained methodology is able to capture the effects triggered by allosteric ligands already described in the literature.

CONCLUSIONS

We introduce a simple computational approach to predict the presence and position of allosteric sites in a protein based on the analysis of changes in protein normal modes upon the binding of a coarse-grained ligand at predicted cavities. Its performance has been demonstrated using a newly curated non-redundant set of 91 proteins with reported allosteric properties. The software developed in this work is available upon request from the authors.

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.

Biochemical characterization of glyceraldehyde-3-phosphate dehydrogenase from Thermococcus kodakarensis KOD1

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) plays an essential role in glycolysis by catalyzing the conversion of D-glyceraldehyde 3-phosphate (D-G3P) to 1,3-diphosphoglycerate using NAD(+) as a cofactor. In this report, the GAPDH gene from the hyperthermophilic archaeon Thermococcus kodakarensis KOD1 (GAPDH-tk) was cloned and the protein was purified to homogeneity. GAPDH-tk exists as a homotetramer with a native molecular mass of 145 kDa; the subunit molecular mass was 37 kDa. GAPDH-tk is a thermostable protein with a half-life of 5 h at 80-90°C. The apparent K (m) values for NAD(+) and D-G3P were 77.8 ± 7.5 μM and 49.3 ± 3.0 μM, respectively, with V (max) values of 45.1 ± 0.8 U/mg and 59.6 ± 1.3 U/mg, respectively. Transmission electron microscopy (TEM) and image processing confirmed that GAPDH-tk has a tetrameric structure. Interestingly, GAPDH-tk migrates as high molecular mass forms (~232 kDa and ~669 kDa) in response to oxidative stress.

Regulation of autotrophic CO2 fixation in the archaeon Thermoproteus neutrophilus

Thermoproteus neutrophilus, a hyperthermophilic, chemolithoautotrophic, anaerobic crenarchaeon, uses a novel autotrophic CO(2) fixation pathway, the dicarboxylate/hydroxybutyrate cycle. The regulation of the central carbon metabolism was studied on the level of whole cells, enzyme activity, the proteome, transcription, and gene organization. The organism proved to be a facultative autotroph, which prefers organic acids as carbon sources that can easily feed into the metabolite pools of this cycle. Addition of the preferred carbon sources acetate, pyruvate, succinate, and 4-hydroxybutyrate to cultures resulted in stimulation of the growth rate and a diauxic growth response. The characteristic enzyme activities of the carbon fixation cycle, fumarate hydratase, fumarate reductase, succinyl coenzyme A (CoA) synthetase, and enzymes catalyzing the conversion of succinyl-CoA to crotonyl-CoA, were differentially downregulated in the presence of acetate and, to a lesser extent, in the presence of other organic substrates. This regulation pattern correlated well with the differential expression profile of the proteome as well as with the transcription of the encoding genes. The genes encoding phosphoenolpyruvate (PEP) carboxylase, fumarate reductase, and four enzymes catalyzing the conversion of succinyl-CoA to crotonyl-CoA are clustered. Two putative operons, one comprising succinyl-CoA reductase plus 4-hydroxybutyrate-CoA ligase genes and the other comprising 4-hydroxybutyryl-CoA dehydratase plus fumarate reductase genes, were divergently transcribed into leaderless mRNAs. The promoter regions were characterized and used for isolating DNA binding proteins. Besides an Alba protein, a 18-kDa protein characteristic for autotrophic Thermoproteales that bound specifically to the promoter region was identified. This system may be suitable for molecular analysis of the transcriptional regulation of autotrophy-related genes.

Complete genome sequence of Thermosphaera aggregans type strain (M11TL)

Thermosphaera aggregans Huber et al. 1998 is the type species of the genus Thermosphaera, which comprises at the time of writing only one species. This species represents archaea with a hyperthermophilic, heterotrophic, strictly anaerobic and fermentative phenotype. The type strain M11TL(T) was isolated from a water-sediment sample of a hot terrestrial spring (Obsidian Pool, Yellowstone National Park, Wyoming). Here we describe the features of this organism, together with the complete genome sequence and annotation. The 1,316,595 bp long single replicon genome with its 1,410 protein-coding and 47 RNA genes is a part of the Genomic Encyclopedia of Bacteria and Archaea project.

Inhibition kinetics of catabolic dehydrogenases by elevated moieties of ATP and ADP--implication for a new regulation mechanism in Lactococcus lactis

ATP and ADP inhibit, in varying degrees, several dehydrogenases of the central carbon metabolism of Lactococcus lactis ATCC 19435 in vitro, i.e. glyceraldehyde-3-phosphate dehydrogenase (GAPDH), lactate dehydrogenase (LDH) and alcohol dehydrogenase (ADH). Here we demonstrate mixed inhibition for GAPDH and competitive inhibition for LDH and ADH by adenine nucleotides in single inhibition studies. The nonlinear negative co-operativity was best modelled with Hill-type kinetics, showing greater flexibility than the usual parabolic inhibition equation. Because these natural inhibitors are present simultaneously in the cytoplasm, multiple inhibition kinetics was determined for each dehydrogenase. For ADH and LDH, the inhibitor combinations ATP plus NAD and ADP plus NAD are indifferent to each other. Model discrimination suggested that the weak allosteric inhibition of GAPDH had no relevance when multiple inhibitors are present. Interestingly, with ADH and GAPDH the combination of ATP and ADP exhibits lower dissociation constants than with either inhibitor alone. Moreover, the concerted inhibition of ADH and GAPDH, but not of LDH, shows synergy between the two nucleotides. Similar kinetics, but without synergies, were found for horse liver and yeast ADHs, indicating that dehydrogenases can be modulated by these nucleotides in a nonlinear manner in many organisms. The action of an elevated pool of ATP and ADP may effectively inactivate lactococcal ADH, but not GAPDH and LDH, providing leverage for the observed metabolic shift to homolactic acid formation in lactococcal resting cells on maltose. Therefore, we interpret these results as a regulation mechanism contributing to readjusting the flux of ATP production in L. lactis.

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.

A global transcriptional regulator in Thermococcus kodakaraensis controls the expression levels of both glycolytic and gluconeogenic enzyme-encoding genes

We identified a novel regulator, Thermococcales glycolytic regulator (Tgr), functioning as both an activator and a repressor of transcription in the hyperthermophilic archaeon Thermococcus kodakaraensis KOD1. Tgr (TK1769) displays similarity (28% identical) to Pyrococcus furiosus TrmB (PF1743), a transcriptional repressor regulating the trehalose/maltose ATP-binding cassette transporter genes, but is more closely related (67%) to a TrmB paralog in P. furiosus (PF0124). Growth of a tgr disruption strain (Deltatgr) displayed a significant decrease in growth rate under gluconeogenic conditions compared with the wild-type strain, whereas comparable growth rates were observed under glycolytic conditions. A whole genome microarray analysis revealed that transcript levels of almost all genes related to glycolysis and maltodextrin metabolism were at relatively high levels in the Deltatgr mutant even under gluconeogenic conditions. The Deltatgr mutant also displayed defects in the transcriptional activation of gluconeogenic genes under these conditions, indicating that Tgr functions as both an activator and a repressor. Genes regulated by Tgr contain a previously identified sequence motif, the Thermococcales glycolytic motif (TGM). The TGM was positioned upstream of the Transcription factor B-responsive element (BRE)/TATA sequence in gluconeogenic promoters and downstream of it in glycolytic promoters. Electrophoretic mobility shift assay indicated that recombinant Tgr protein specifically binds to promoter regions containing a TGM. Tgr was released from the DNA when maltotriose was added, suggesting that this sugar is most likely the physiological effector. Our results strongly suggest that Tgr is a global transcriptional regulator that simultaneously controls, in response to sugar availability, both glycolytic and gluconeogenic metabolism in T. kodakaraensis via its direct binding to the TGM.

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.

Cloning, gene expression and characterization of a novel bacterial NAD-dependent non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase from Neisseria meningitidis strain Z2491

Alignment of the amino acid sequence of some archaeal, bacterial and eukaryotic non-phosphorylating glyceraldehydes-3-phosphate dehydrogenases (GAPNs) and aldehyde dehydrogenases (ALDHs) with the sequence of a putative GAPN present in the genome of the Gram-negative bacterium Neisseria meningitidis strain Z2491 demonstrated the conservation of residues involved in the catalytic activity. The predicted coding sequence of the N. meningitidis gapN gene was cloned in Escherichia coli XL1-blue under the expression of an inducible promoter. The IPTG-induced GAPN was purified ca. 48-fold from E. coli cells using a procedure that sequentially employed conventional ammonium sulfate fractionation as well as anion-exchange and affinity chromatography. The purified recombinant enzyme was thoroughly characterized. The protein is a homotetramer with a 50-kDa subunit, exhibiting absolute specificity for NAD and a broad spectrum of aldehyde substrates. Isoelectric focusing analysis with the purified fraction showed the presence of an acidic polypeptide with an isoelectric point of 6.3. The optimum pH of the purified enzyme was between 9 and 10. Studies on the effect of increasing temperatures on the enzyme activity revealed an optimal value ca. 64 degrees C. Molecular phylogenetic data suggest that N. meningitidis GAPN has a closer relationship with archaeal GAPNs and glyceraldehyde dehydrogenases than with the typical NADP-specific GAPNs from Gram-positive bacteria and photosynthetic eukaryotes.

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

Archaea utilize a branched modification of the classical Entner-Doudoroff (ED) pathway for sugar degradation. The semi-phosphorylative branch merges at the level of glyceraldehyde 3-phosphate (GAP) with the lower common shunt of the Emden-Meyerhof-Parnas pathway. In Sulfolobus solfataricus two different GAP converting enzymes-classical phosphorylating GAP dehydrogenase (GAPDH) and the non-phosphorylating GAPDH (GAPN)-were identified. In Sulfolobales the GAPN encoding gene is found adjacent to the ED gene cluster suggesting a function in the regulation of the semi-phosphorylative ED branch. The biochemical characterization of the recombinant GAPN of S. solfataricus revealed that-like the well-characterized GAPN from Thermoproteus tenax-the enzyme of S. solfataricus exhibits allosteric properties. However, both enzymes show some unexpected differences in co-substrate specificity as well as regulatory fine-tuning, which seem to reflect an adaptation to the different lifestyles of both organisms. Phylogenetic analyses and database searches in Archaea indicated a preferred distribution of GAPN (and/or GAP oxidoreductase) in hyperthermophilic Archaea supporting the previously suggested role of GAPN in metabolic thermoadaptation. This work suggests an important role of GAPN in the regulation of carbon degradation via modifications of the EMP and the branched ED pathway in hyperthermophilic Archaea.

The phosphofructokinase-B (MJ0406) from Methanocaldococcus jannaschii represents a nucleoside kinase with a broad substrate specificity

Recently, unusual non-regulated ATP-dependent 6-phosphofructokinases (PFK) that belong to the PFK-B family have been described for the hyperthermophilic archaea Desulfurococcus amylolyticus and Aeropyrum pernix. Putative homologues were found in genomes of several archaea including the hyperthermophilic archaeon Methanocaldococcus jannaschii. In this organism, open reading frame MJ0406 had been annotated as a PFK-B sugar kinase. The gene encoding MJ0406 was cloned and functionally expressed in Escherichia coli. The purified recombinant enzyme is a homodimer with an apparent molecular mass of 68 kDa composed of 34 kDa subunits. With a temperature optimum of 85 degrees C and a melting temperature of 90 degrees C, the M. jannaschii nucleotide kinase represents one of the most thermoactive and thermostable members of the PFK-B family described so far. The recombinant enzyme was characterized as a functional nucleoside kinase rather than a 6-PFK. Inosine, guanosine, and cytidine were the most effective phosphoryl acceptors. Besides, adenosine, thymidine, uridin and xanthosine were less efficient. Extremely low activity was found with fructose-6-phosphate. Further, the substrate specificity of closely related PFK-Bs from D. amylolyticus and A. pernix were reanalysed.

Identification of lactaldehyde dehydrogenase in Methanocaldococcus jannaschii and its involvement in production of lactate for F420 biosynthesis

One of the early steps in the biosynthesis of coenzyme F(420) in Methanocaldococcus jannaschii requires generation of 2-phospho-L-lactate, which is formed by the phosphorylation of L-lactate. Preliminary studies had shown that L-lactate in M. jannaschii is not derived from pyruvate, and thus an alternate pathway(s) for its formation was examined. Here we report that L-lactate is formed by the NAD(+)-dependent oxidation of l-lactaldehyde by the MJ1411 gene product. The lactaldehyde, in turn, was found to be generated either by the NAD(P)H reduction of methylglyoxal or by the aldol cleavage of fuculose-1-phosphate by fuculose-1-phosphate aldolase, the MJ1418 gene product.

Phosphoenolpyruvate synthetase and pyruvate, phosphate dikinase of Thermoproteus tenax: key pieces in the puzzle of archaeal carbohydrate metabolism

The interconversion of phosphoenolpyruvate and pyruvate represents an important control point of the Embden-Meyerhof-Parnas (EMP) pathway in Bacteria and Eucarya, but little is known about this site of regulation in Archaea. Here we report on the coexistence of phosphoenolpyruvate synthetase (PEPS) and the first described archaeal pyruvate, phosphate dikinase (PPDK), which, besides pyruvate kinase (PK), are involved in the catalysis of this reaction in the hyperthermophilic crenarchaeote Thermoproteus tenax. The genes encoding T. tenax PEPS and PPDK were cloned and expressed in Escherichia coli, and the enzymic and regulatory properties of the recombinant gene products were analysed. Whereas PEPS catalyses the unidirectional conversion of pyruvate to phosphoenolpyruvate, PPDK shows a bidirectional activity with a preference for the catabolic reaction. In contrast to PK of T. tenax, which is regulated on transcript level but exhibits only limited regulatory potential on protein level, PEPS and PPDK activities are modulated by adenosine phosphates and intermediates of the carbohydrate metabolism. Additionally, expression of PEPS is regulated on transcript level in response to the offered carbon source as revealed by Northern blot analyses. The combined action of the differently regulated enzymes PEPS, PPDK and PK represents a novel way of controlling the interconversion of phosphoenolpyruvate and pyruvate in the reversible EMP pathway, allowing short-term and long-term adaptation to different trophic conditions. Comparative genomic analyses indicate the coexistence of PEPS, PPDK and PK in other Archaea as well, suggesting a similar regulation of the carbohydrate metabolism in these organisms.

Reconstruction of central carbon metabolism inSulfolobus solfataricus using a two-dimensional gel electrophoresis map, stable isotope labelling and DNA microarray analysis

In the last decade, an increasing number of sequenced archaeal genomes have become available, opening up the possibility for functional genomic analyses. Here, we reconstructed the central carbon metabolism in the hyperthermophilic crenarchaeon Sulfolobus solfataricus (glycolysis, gluconeogenesis and tricarboxylic acid cycle) on the basis of genomic, proteomic, transcriptomic and biochemical data. A 2-DE reference map of S. solfataricus grown on glucose, consisting of 325 unique ORFs in 255 protein spots, was created to facilitate this study. The map was then used for a differential expression study based on (15)N metabolic labelling (yeast extract + tryptone-grown cells (YT) vs. glucose-grown cells (G)). In addition, the expression ratio of the genes involved in carbon metabolism was studied using DNA microarrays. Surprisingly, only 3 and 14% of the genes and proteins, respectively, involved in central carbon metabolism showed a greater than two-fold change in expression level. All results are discussed in the light of the current understanding of central carbon metabolism in S. solfataricus and will help to obtain a system-wide understanding of this organism.

Effects of lactate dehydrogenase suppression and glycerol-3-phosphate dehydrogenase overexpression on cellular metabolism

In order to conduct a physiological functional study of lactate dehydrogenase (LDH) and glycerol-3-phosphate dehydrogenase (GPDH), we engineered a CHO dhfr(-) cell, by overexpressing either the anti-sense LDH-A RNA (anti-LDH cells) or GPDH (GP3 cells), or both (GP3/anti-LDH cells). LDH activity in the cell cytosol, and lactate content and pHe change in the growth media were found to decrease according to the order: cell lines GP3/anti-LDH > anti-LDH > GP3 > CHO. Intracellular ATP contents, representing the extent of respiration rate, also decreased, according to a rank order as follows: GP3 > CHO > GP3/anti-LDH > anti-LDH. We also attempted to identify and characterize any physiological changes occurring in the cells which harbored diverse metabolic pathways. First, anti-LDH cells with heightened respiration rates were found to display a higher degree of sensitivity to the prooxidant tert-butyl hydroperoxide (tBOOH), and the mitochondrial complex III inhibitor, antimycin A, than the GPDH-expressing cells (GP3 and GP3/anti-LDH), which have a lower respiration rate. Second, the anti-sense LDH-A RNA-expressing cells (anti-LDH and GP3/anti-LDH) evidenced a higher degree of resistance to apoptosis by cell-cell contact inhibition, and a faster doubling time ( approximately 19 h compared with approximately 26 h) than the CHO and GP3 cells. Additionally, cell growth in an extended culture under HCO(3) (-)-free conditions to induce a steep acidification could be maintained with the anti-sense LDH-A RNA-expressing cells, but could not be maintained with the CHO and GP3 cells. Third, we observed that the most appropriate cell line for the optical production of a certain therapeutic protein (Tissue-Plasminogen Activator) was the GP3/anti-LDH cells. Collectively, our data indicate a variety of physiological roles for LDH and GPDH, including cellular acidosis, oxidoresistance, apoptosis by both acidosis and cell-cell contact inhibition, cell growth, and the generation of recombinant proteins.

The semi-phosphorylative Entner-Doudoroff pathway in hyperthermophilic archaea: a re-evaluation

Biochemical studies have suggested that, in hyperthermophilic archaea, the metabolic conversion of glucose via the ED (Entner-Doudoroff) pathway generally proceeds via a non-phosphorylative variant. A key enzyme of the non-phosphorylating ED pathway of Sulfolobus solfataricus, KDG (2-keto-3-deoxygluconate) aldolase, has been cloned and characterized previously. In the present study, a comparative genomics analysis is described that reveals conserved ED gene clusters in both Thermoproteus tenax and S. solfataricus. The corresponding ED proteins from both archaea have been expressed in Escherichia coli and their specificity has been identified, revealing: (i) a novel type of gluconate dehydratase (gad gene), (ii) a bifunctional 2-keto-3-deoxy-(6-phospho)-gluconate aldolase (kdgA gene), (iii) a 2-keto-3-deoxygluconate kinase (kdgK gene) and, in S. solfataricus, (iv) a GAPN (non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase; gapN gene). Extensive in vivo and in vitro enzymatic analyses indicate the operation of both the semi-phosphorylative and the non-phosphorylative ED pathway in T. tenax and S. solfataricus. The existence of this branched ED pathway is yet another example of the versatility and flexibility of the central carbohydrate metabolic pathways in the archaeal domain.

Unusual pathways and enzymes of central carbohydrate metabolism in Archaea

Sugar-utilizing hyperthermophilic and halophilic Archaea degrade glucose and glucose polymers to acetate or to CO2 using O2, nitrate, sulfur or sulfate as electron acceptors. Comparative analyses of glycolytic pathways in these organisms indicate a variety of differences from the classical Emden-Meyerhof and Entner-Doudoroff pathways that are operative in Bacteria and Eukarya, respectively. The archaeal pathways are characterized by the presence of numerous novel enzymes and enzyme families that catalyze, for example, the phosphorylation of glucose and of fructose 6-phosphate, the isomerization of glucose 6-phosphate, the cleavage of fructose 1,6-bisphosphate, the oxidation of glyceraldehyde 3-phosphate and the conversion of acetyl-CoA to acetate. Recent major advances in deciphering the complexity of archaeal central carbohydrate metabolism were gained by combination of classical biochemical and genomic-based approaches.

Corynebacterium glutamicum glyceraldehyde-3-phosphate dehydrogenase isoforms with opposite, ATP-dependent regulation

Corynebacterium glutamicum gapA and gapB encode glyceraldehyde-3-phosphate dehydrogenases (GAPDHs) that differ in molecular weight and activity in the presence of ATP. Comparative genome analysis revealed that GapA, the product of gapA, represented the canonical GAPDH that is highly conserved across the three major life forms. GapB, with an additional 110-residue-long sequence upstream of its GAPDH-specific domain, was homologous only to select microbial putative GAPDHs. Upon gene disruption, the initial growth rates of the wild-type, DeltagapA and DeltagapB strains on glucose (0.77, 0.00 and 0.76 h(-1), respectively), lactate (0.20, 0.18 and 0.15 h(-1), respectively), pyruvate (0.39, 0.29 and 0.20 h(-1), respectively), and acetate (0.06, 0.06 and 0.04 h(-1), respectively), implied that GapA was indispensable for growth on glucose, that GapB, but not GapA, affected early growth on acetate, and that GapB had a greater influence on growth under gluconeogenic conditions than GapA. The disruption of either gapA or gapB showed no significant effect on the transcription of any of the other gap cluster genes although it led to reduced triosephosphate isomerase (TPI) activities. Glycolytic GAPDH activity at low in vitro ATP concentrations was solely attributed to the 35.9-kDa GapA. At higher ATP concentrations, the same activity was attributed to the 51.2-kDa GapB. Both enzymes, however, exhibited similar NADP-dependent GAPDH activities at the higher ATP concentrations. In effect therefore, the GAPDH-catalyzed reaction at low ATP concentrations was irreversible, with all the glycolytic activity strictly NAD-dependent and attributed to GapA. At higher ATP concentrations, the reaction was reversible, with glycolytic activity NAD- or NADP-dependent and attributed to GapB, while gluconeogenic activity was attributable to both GapA and GapB.

Transaldolase of Methanocaldococcus jannaschii

The Methanocaldococcus jannaschii genome contains putative genes for all four nonoxidative pentose phosphate pathway enzymes. Open reading frame (ORF) MJ0960 is a member of the mipB/talC family of 'transaldolase-like' genes, so named because of their similarity to the well-characterized transaldolase B gene family. However, recently, it has been reported that both the mipB and the talC genes from Escherichia coli encode novel enzymes with fructose-6-phosphate aldolase activity, not transaldolase activity (Schürmann and Sprenger 2001). The same study reports that other members of the mipB/talC family appear to encode transaldolases. To confirm the function of MJ0960 and to clarify the presence of a nonoxidative pentose phosphate pathway in M. jannaschii, we have cloned ORF MJ0960 from M. jannaschii genomic DNA and purified the recombinant protein. MJ0960 encodes a transaldolase and displays no fructose-6-phosphate aldolase activity. It etained full activity for 4 h at 80 degrees C, and for 3 weeks at 25 degrees C. Methanocaldococcus jannaschii transaldolase has a maximal velocity (Vmax) of 1.0 +/- 0.2 micromol min(-1) mg(-1) at 25 degrees C, whereas Vmax = 12.0 +/- 0.5 micromol min(-1) mg(-1) at 50 degrees C. Apparent Michaelis constants at 50 degrees C were Km = 0.65 +/- 0.09 mM for fructose-6-phosphate and Km = 27.8 +/- 4.3 microM for erythrose-4-phosphate. When ribose-5-phosphate replaced erythrose-4-phosphate as an aldose acceptor, Vmax decreased twofold, whereas the Km was 150-fold higher. The molecular mass of the active enzyme is 271 +/- 27 kDa as estimated by gel filtration, whereas the predicted monomer size is 23.96 kDa, suggesting that the native form of the protein is probably a decamer. A readily available source of thermophilic pentose phosphate pathway enzymes including transaldolase may have direct application in enzymatic biohydrogen production.

Distribution and phylogenies of enzymes of the Embden-Meyerhof-Parnas pathway from archaea and hyperthermophilic bacteria support a gluconeogenic origin of metabolism

Enzymes of the gluconeogenic/glycolytic pathway (the Embden-Meyerhof-Parnas (EMP) pathway), the reductive tricarboxylic acid cycle, the reductive pentose phosphate cycle and the Entner-Doudoroff pathway are widely distributed and are often considered to be central to the origins of metabolism. In particular, several enzymes of the lower portion of the EMP pathway (the so-called trunk pathway), including triosephosphate isomerase (TPI; EC 5.3.1.1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH; EC 1.2.1.12/13), phosphoglycerate kinase (PGK; EC 2.7.2.3) and enolase (EC 4.2.1.11), are extremely well conserved and universally distributed among the three domains of life. In this paper, the distribution of enzymes of gluconeogenesis/glycolysis in hyperthermophiles--microorganisms that many believe represent the least evolved organisms on the planet--is reviewed. In addition, the phylogenies of the trunk pathway enzymes (TPIs, GAPDHs, PGKs and enolases) are examined. The enzymes catalyzing each of the six-carbon transformations in the upper portion of the EMP pathway, with the possible exception of aldolase, are all derived from multiple gene sequence families. In contrast, single sequence families can account for the archaeal and hyperthermophilic bacterial enzyme activities of the lower portion of the EMP pathway. The universal distribution of the trunk pathway enzymes, in combination with their phylogenies, supports the notion that the EMP pathway evolved in the direction of gluconeogenesis, i.e., from the bottom up.

Structural Basis of allosteric regulation and substrate specificity of the non-phosphorylating glyceraldehyde 3-Phosphate dehydrogenase from Thermoproteus tenax

The non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN) of the hyperthermophilic Archaeum Thermoproteus tenax is a member of the superfamily of aldehyde dehydrogenases (ALDH). GAPN catalyses the irreversible oxidation of glyceraldehyde 3-phosphate (GAP) to 3-phosphoglycerate in the modified glycolytic pathway of this organism. In contrast to other members of the ALDH superfamily, GAPN from T.tenax (Tt-GAPN) is regulated by a number of intermediates and metabolites. In the NAD-dependent oxidation of GAP, glucose 1-phosphate, fructose 6-phosphate, AMP and ADP increase the affinity for the cosubstrate, whereas ATP, NADP, NADPH and NADH decrease it leaving, however, the catalytic rate virtually unaltered. As we show here, the enzyme also uses NADP as a cosubstrate, displaying, however, unusual discontinuous saturation kinetics indicating different cosubstrate affinities and/or reactivities of the four active sites of the protein tetramer caused by cooperative effects. Furthermore, in the NADP-dependent reaction the presence of activators decreases the overall S0.5 and increases Vmax by a factor of 3. To explore the structural basis for the different effects of both pyridine nucleotides we solved the crystal structure of Tt-GAPN in complex with NAD at 2.2 A resolution and compared it to the binary Tt-GAPN-NADPH structure. Although both pyridine nucleotides show a similar binding mode, NADPH appears to be more tightly bound to the protein via the 2' phosphate moiety. Moreover, we present four co-crystal structures with the activating molecules glucose 1-phosphate, fructose 6-phosphate, AMP and ADP determined at resolutions ranging from 2.3 A to 2.6 A. These crystal structures reveal a common regulatory site able to accommodate the different activators. A phosphate-binding pocket serves as an anchor point ensuring similar binding geometry. The observed conformational changes upon activator binding are discussed in terms of allosteric regulation. Furthermore, we present a crystal structure of Tt-GAPN in complex with the substrate D-GAP at 2.3 A resolution, which allows us to analyse the structural basis for substrate binding, the mechanism of catalysis as well as the stereoselectivity of the enzymatic reaction.

Structure and Function of a Regulated Archaeal Triosephosphate Isomerase Adapted to High Temperature

Triosephophate isomerase (TIM) is a dimeric enzyme in eucarya, bacteria and mesophilic archaea. In hyperthermophilic archaea, however, TIM exists as a tetramer composed of monomers that are about 10% shorter than other eucaryal and bacterial TIM monomers. We report here the crystal structure of TIM from Thermoproteus tenax, a hyperthermophilic archaeon that has an optimum growth temperature of 86 degrees C. The structure was determined from both a hexagonal and an orthorhombic crystal form to resolutions of 2.5A and 2.3A, and refined to R-factors of 19.7% and 21.5%, respectively. In both crystal forms, T.tenax TIM exists as a tetramer of the familiar (betaalpha)(8)-barrel. In solution, however, and unlike other hyperthermophilic TIMs, the T.tenax enzyme exhibits an equilibrium between inactive dimers and active tetramers, which is shifted to the tetramer state through a specific interaction with glycerol-1-phosphate dehydrogenase of T.tenax. This observation is interpreted in physiological terms as a need to reduce the build-up of thermolabile metabolic intermediates that would be susceptible to destruction by heat. A detailed structural comparison with TIMs from organisms with growth optima ranging from 15 degrees C to 100 degrees C emphasizes the importance in hyperthermophilic proteins of the specific location of ionic interactions for thermal stability rather than their numbers, and shows a clear correlation between the reduction of heat-labile, surface-exposed Asn and Gln residues with thermoadaptation. The comparison confirms the increase in charged surface-exposed residues at the expense of polar residues.

Reconstruction of the central carbohydrate metabolism of Thermoproteus tenax by use of genomic and biochemical data

The hyperthermophilic, facultatively heterotrophic crenarchaeum Thermoproteus tenax was analyzed using a low-coverage shotgun-sequencing approach. A total of 1.81 Mbp (representing 98.5% of the total genome), with an average gap size of 100 bp and 5.3-fold coverage, are reported, giving insights into the genome of T. tenax. Genome analysis and biochemical studies enabled us to reconstruct its central carbohydrate metabolism. T. tenax uses a variant of the reversible Embden-Meyerhof-Parnas (EMP) pathway and two different variants of the Entner-Doudoroff (ED) pathway (a nonphosphorylative variant and a semiphosphorylative variant) for carbohydrate catabolism. For the EMP pathway some new, unexpected enzymes were identified. The semiphosphorylative ED pathway, hitherto supposed to be active only in halophiles, is found in T. tenax. No evidence for a functional pentose phosphate pathway, which is essential for the generation of pentoses and NADPH for anabolic purposes in bacteria and eucarya, is found in T. tenax. Most genes involved in the reversible citric acid cycle were identified, suggesting the presence of a functional oxidative cycle under heterotrophic growth conditions and a reductive cycle for CO2 fixation under autotrophic growth conditions. Almost all genes necessary for glycogen and trehalose metabolism were identified in the T. tenax genome.