Pyruvate kinase of the hyperthermophilic crenarchaeote Thermoproteus tenax: physiological role and phylogenetic aspects

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

Thermostable in vitro transcription-translation compatible with microfluidic droplets

BACKGROUND

In vitro expression involves the utilization of the cellular transcription and translation machinery in an acellular context to produce one or more proteins of interest and has found widespread application in synthetic biology and in pharmaceutical biomanufacturing. Most in vitro expression systems available are active at moderate temperatures, but to screen large libraries of natural or artificial genetic diversity for highly thermostable enzymes or enzyme variants, it is instrumental to enable protein synthesis at high temperatures.

OBJECTIVES

Develop an in vitro expression system operating at high temperatures compatible with enzymatic assays and with technologies that enable ultrahigh-throughput protein expression in reduced volumes, such as microfluidic water-in-oil (w/o) droplets.

RESULTS

We produced cell-free extracts from Thermus thermophilus for in vitro translation including thermostable enzymatic cascades for energy regeneration and a moderately thermostable RNA polymerase for transcription, which ultimately limited the temperature of protein synthesis. The yield was comparable or superior to other thermostable in vitro expression systems, while the preparation procedure is much simpler and can be suited to different Thermus thermophilus strains. Furthermore, these extracts have enabled in vitro expression in microfluidic droplets at high temperatures for the first time.

CONCLUSIONS

Cell-free extracts from Thermus thermophilus represent a simpler alternative to heavily optimized or pure component thermostable in vitro expression systems. Moreover, due to their compatibility with droplet microfluidics and enzyme assays at high temperatures, the reported system represents a convenient gateway for enzyme screening at higher temperatures with ultrahigh-throughput.

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family-pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization-including quantitative evaluation of allosteric regulation-of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms. Database URL https://github.com/djparente/PYK-DB.

Genome-Wide Analysis and Functional Characterization of Pyruvate Kinase (PK) Gene Family Modulating Rice Yield and Quality

Pyruvate kinase (PK) is one of the three rate-limiting enzymes of glycolysis, and it plays a pivotal role in energy metabolism. In this study, we have identified 10 PK genes from the rice genome. Initially, these genes were divided into two categories: cytoplasmic pyruvate kinase (PKc) and plastid pyruvate kinase (PKp). Then, an expression analysis revealed that OsPK1, OsPK3, OsPK4, OsPK6, and OsPK9 were highly expressed in grains. Moreover, PKs can form heteropolymers. In addition, it was found that ABA significantly regulates the expression of PK genes (OsPK1, OsPK4, OsPK9, and OsPK10) in rice. Intriguingly, all the genes were found to be substantially involved in the regulation of rice grain quality and yield. For example, the disruption of OsPK3, OsPK5, OsPK7, OsPK8, and OsPK10 and OsPK4, OsPK5, OsPK6, and OsPK10 decreased the 1000-grain weight and the seed setting rate, respectively. Further, the disruption of OsPK4, OsPK6, OsPK8, and OsPK10 through the CRISPR/Cas9 system showed an increase in the content of total starch and a decrease in protein content compared to the WT. Similarly, manipulations of the OsPK4, OsPK8, and OsPK10 genes increased the amylose content. Meanwhile, the grains of all CRISPR mutants and RNAi lines, except ospk6, showed a significant increase in the chalkiness rate compared to the wild type. Overall, this study characterizes the functions of all the genes of the PK gene family and shows their untapped potential to improve rice yield and quality traits.

Odd one out? Functional tuning of Zymomonas mobilis pyruvate kinase is narrower than its allosteric, human counterpart

Various protein properties are often illuminated using sequence comparisons of protein homologs. For example, in analyses of the pyruvate kinase multiple sequence alignment, the set of positions that changed during speciation ("phylogenetic" positions) were enriched for "rheostat" positions in human liver pyruvate kinase (hLPYK). (Rheostat positions are those which, when substituted with various amino acids, yield a range of functional outcomes). However, the correlation was moderate, which could result from multiple biophysical constraints acting on the same position during evolution and/or various sources of noise. To further examine this correlation, we here tested Zymomonas mobilis PYK (ZmPYK), which has <65% sequence identity to any other PYK sequence. Twenty-six ZmPYK positions were selected based on their phylogenetic scores, substituted with multiple amino acids, and assessed for changes in Kapp-PEP . Although we expected to identify multiple, strong rheostat positions, only one moderate rheostat position was detected. Instead, nearly half of the 271 ZmPYK variants were inactive and most others showed near wild-type function. Indeed, for the active ZmPYK variants, the total range of Kapp,PEP values ("tunability") was 40-fold less than that observed for hLPYK variants. The combined functional studies and sequence comparisons suggest that ZmPYK has evolved functional and/or structural attributes that differ from the rest of the family. We hypothesize that including such "orphan" sequences in MSA analyses obscures the correlations used to predict rheostat positions. Finally, results raise the intriguing biophysical question as to how the same protein fold can support rheostat positions in one homolog but not another.

Molecular cloning, expression in Escherichia coli and structural-functional analysis of a pyruvate kinase from Pyrobaculum calidifontis

This manuscript describes recombinant production, characterization and structural analysis of wild-type and mutant Pcal_0029, a pyruvate kinase from Pyrobaculum calidifontis. Recombinant Pcal_0029 was produced in soluble and highly active form in Escherichia coli. Purified protein exhibited divalent metal-dependent activity which increased with the increase in temperature till 85 °C. Recombinant Pcal_0029 was highly thermostable with no significant loss in activity even after an incubation of 120 min at 100 °C. The enzyme exhibited apparent S_0.5 and V_max values of 0.44 ± 0.05 mM and 840 ± 39 units, respectively, towards phosphoenolpyruvate. These values towards adenosine-5'-diphosphate were 0.5 ± 0.07 mM and 870 ± 26 units, respectively. In silico structural analysis and comparison with the characterized enzymes revealed the presence of eight conserved regions. Two substitutions, K130E and S155G, resulted in a 10-fold decrease in activity. Secondary structure analysis indicated similar structures for the wild-type and the mutant enzymes. Bioinformatics analysis revealed disruption of interatomic interactions and hydrogen bond formation, leading to a decreased flexibility and solvent accessibility, which may have led to decrease in activity. To the best of our knowledge, Pcal_0029 is the most thermostable pyruvate kinase reported so far. Moreover, this is the first study on the role of non-catalytic residues in a pyruvate kinase.

Exploring species-specific inhibitors with multiple target sites on S. aureus pyruvate kinase using a computational workflow

Experimental evidence indicated that bacterial pyruvate kinase of glycolysis can be evaluated as an alternative target to eliminate infections, while antibiotic resistance poses a global threat. Here, we use a computational workflow to reveal and investigate the potential allosteric sites of methicillin-resistant S. aureus PK, which can help in designing species-specific drugs to inhibit activity of this organism. Residue interaction networks point to a known allosteric site at the small C-C interface, a potential allosteric site near the small interface (site #1), and a second potential allosteric site at the large interface (site #2). 2 µs-long molecular dynamics (MD) simulations with AMBER16 generate different conformations of one narrow target site. Known and potential allosteric sites on the selected conformers are investigated using ensemble docking with AutoDock Vina and a library of 2447 FDA-approved drugs. We determine 18 hits, comprising ergot-alkaloids, anti-cancer-agents, antivirals, analgesics, cardiac glycosides, all with a high docking z-score for three sites. 5 selected compounds with high, average and low z-scores are subjected to 50 ns-long MD simulations for MM-GBSA calculations. ΔG_bind values up to -49.3 kcal/mol at the C-C interface, up to -32.7 kcal/mol at site #1, and up to -53.3 kcal/mol at site #2 support the docking calculations. We investigate mitapivat and TT-232 as reference compounds under clinical trial, targeting human PK isomers. We suggest 18 FDA-approved hits from the docking calculations and TT-232 as potential inhibitors with multiple target sites on S. aureus PK. This study also proposes pharmacophores models for de novo drug design.Communicated by Ramaswamy H. Sarma.

The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms

Eukarya pyruvate kinases possess glutamate at position 117 (numbering of rabbit muscle enzyme), whereas bacteria have either glutamate or lysine. Those with E117 are K⁺-dependent, whereas those with K117 are K⁺-independent. In a phylogenetic tree, 80% of the sequences with E117 are occupied by T113/K114/T120 and 77% of those with K117 possess L113/Q114/(L,I,V)120. This work aims to understand these residues’ contribution to the K⁺-independent pyruvate kinases using the K⁺-dependent rabbit muscle enzyme. Residues 117 and 120 are crucial in the differences between the K⁺-dependent and -independent mutants. K⁺-independent activity increased with L113 and Q114 to K117, but L120 induced structural differences that inactivated the enzyme. T120 appears to be key in folding the protein and closure of the lid of the active site to acquire its active conformation in the K⁺-dependent enzymes. E117K mutant was K⁺-independent and the enzyme acquired the active conformation by a different mechanism. In the K⁺-independent apoenzyme of Mycobacterium tuberculosis, K72 (K117) flips out of the active site; in the holoenzyme, K72 faces toward the active site bridging the substrates through water molecules. The results provide evidence that two different mechanisms have evolved for the catalysis of this reaction.

The PEP-pyruvate-oxaloacetate node: variation at the heart of metabolism

At the junction between the glycolysis and the tricarboxylic acid cycle-as well as various other metabolic pathways-lies the phosphoenolpyruvate (PEP)-pyruvate-oxaloacetate node (PPO-node). These three metabolites form the core of a network involving at least eleven different types of enzymes, each with numerous subtypes. Obviously, no single organism maintains each of these eleven enzymes; instead, different organisms possess different subsets in their PPO-node, which results in a remarkable degree of variation, despite connecting such deeply conserved metabolic pathways as the glycolysis and the tricarboxylic acid cycle. The PPO-node enzymes play a crucial role in cellular energetics, with most of them involved in (de)phosphorylation of nucleotide phosphates, while those responsible for malate conversion are important redox enzymes. Variations in PPO-node therefore reflect the different energetic niches that organisms can occupy. In this review, we give an overview of the biochemistry of these eleven PPO-node enzymes. We attempt to highlight the variation that exists, both in PPO-node compositions, as well as in the roles that the enzymes can have within those different settings, through various recent discoveries in both bacteria and archaea that reveal deviations from canonical functions.

Biochemical and biophysical characterization of the smallest pyruvate kinase from Entamoeba histolytica

Entamoeba histolytica infection is highly prevalent in developing countries across the globe. The ATP synthesis in this pathogen is solely dependent on the glycolysis pathway where pyruvate kinase (Pyk) catalyzes the final reaction. Here, we have cloned, overexpressed and purified the pyruvate kinase (EhPyk) from E. histolytica. EhPyk is the shortest currently known Pyk till date as it contains only two of the three characterized domains when compared to the other homologues and our phylogenetic analysis places it on a distinct branch from the known type I/II Pyks. Our purification results suggested that it exists as a homodimer in solution. The kinetic characterization showed that EhPyk has maximum activity at pH 7.5 where it exhibited Michaelis-Menten's kinetics for phosphoenolpyruvate with a K_m of 0.23 mM, and it lost its activity at both the acidic pH 4.0 and basic pH 10.0. We also determined the key secondary structural elements of EhPyk at different pH values. MD simulation of EhPyk structure at different pH values suggested that it is most stable at pH 7.0, while least stable at pH 10.0 followed by pH 4.0. Together, our computational simulations correlate well with the experimental studies. In summary, this study expands the current understanding of the EhPyk identified earlier in the amoebic genome and provides the first characterization of this bacterially expressed protein.

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

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

An overview of structure, function, and regulation of pyruvate kinases

In the last step of glycolysis Pyruvate kinase catalyzes the irreversible conversion of ADP and phosphoenolpyruvate to ATP and pyruvic acid, both crucial for cellular metabolism. Thus pyruvate kinase plays a key role in controlling the metabolic flux and ATP production. The hallmark of the activity of different pyruvate kinases is their tight modulation by a variety of mechanisms including the use of a large number of physiological allosteric effectors in addition to their homotropic regulation by phosphoenolpyruvate. Binding of effectors signals precise and orchestrated movements in selected areas of the protein structure that alter the catalytic action of these evolutionarily conserved enzymes with remarkably conserved architecture and sequences. While the diverse nature of the allosteric effectors has been discussed in the literature, the structural basis of their regulatory effects is still not well understood because of the lack of data representing conformations in various activation states. Results of recent studies on pyruvate kinases of different families suggest that members of evolutionarily related families follow somewhat conserved allosteric strategies but evolutionarily distant members adopt different strategies. Here we review the structure and allosteric properties of pyruvate kinases of different families for which structural data are available.

New views on an old enzyme: allosteric regulation and evolution of archaeal pyruvate kinases

Pyruvate kinases (PKs) synthesize ATP as the final step of glycolysis in the three domains of life. PKs from most bacteria and eukarya are allosteric enzymes that are activated by sugar phosphates; for example, the feed-forward regulator fructose-1,6-bisphosphate, or AMP as a sensor of energy charge. Archaea utilize unusual glycolytic pathways, but the allosteric properties of PKs from these species are largely unknown. Here, we present an analysis of 24 PKs from most archaeal clades with respect to allosteric properties, together with phylogenetic analyses constructed using a novel mode of rooting protein trees. We find that PKs from many Thermoproteales, an order of crenarchaeota, are allosterically activated by 3-phosphoglycerate (3PG). We also identify five conserved amino acids that form the binding pocket for 3PG. 3PG is generated via an irreversible reaction in the modified glycolytic pathway of these archaea and therefore functions as a feed-forward regulator. We also show that PKs from hyperthermophilic Methanococcales, an order of euryarchaeota, are activated by AMP. Phylogenetic analyses indicate that 3PG-activated PKs form an evolutionary lineage that is distinct from that of sugar-phosphate activated PKs, and that sugar phosphate-activated PKs originated as AMP-regulated PKs in hyperthermophilic Methanococcales. Since the phospho group of sugar phosphates and 3PG overlap in the allosteric site, our data indicate that the allostery in PKs first started from a progenitor phosphate-binding site that evolved in two spatially distinct directions: one direction generated the canonical site that responds to sugar phosphates and the other gave rise to the 3PG site present in Thermoproteales. Overall, our data suggest an intimate connection between the allosteric properties and evolution of PKs.

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

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

OsPK2 encodes a plastidic pyruvate kinase involved in rice endosperm starch synthesis, compound granule formation and grain filling

Starch is the main form of energy storage in higher plants. Although several enzymes and regulators of starch biosynthesis have been defined, the complete molecular machinery remains largely unknown. Screening for irregularities in endosperm formation in rice represents valuable prospect for studying starch synthesis pathway. Here, we identified a novel rice white-core endosperm and defective grain filling mutant, ospk2, which displays significantly lower grain weight, decreased starch content and alteration of starch physicochemical properties when compared to wild-type grains. The normal starch compound granules were drastically reduced and more single granules filled the endosperm cells of ospk2. Meanwhile, the germination rate of ospk2 seeds after 1-year storage was observably reduced compared with wild-type. Map-based cloning of OsPK2 indicated that it encodes a pyruvate kinase (PK, ATP: pyruvate 2-O-phosphotransferase, EC 2.7.1.40), which catalyses an irreversible step of glycolysis. OsPK2 has a constitutive expression in rice and its protein localizes in chloroplasts. Enzyme assay showed that the protein product from expressed OsPK2 and the crude protein extracted from tissues of wild-type exhibits strong PK activity; however, the mutant presented reduced protein activity. OsPK2 (PKpα1) and three other putative rice plastidic isozymes, PKpα2, PKpβ1 and PKpβ2, can interact to form heteromer. Moreover, the mutation leads to multiple metabolic disorders. Altogether, these results denote new insights into the role of OsPK2 in plant seed development, especially in starch synthesis, compound granules formation and grain filling, which would be useful for genetic improvement of high yield and rice grain quality.

Glycolysis as the Central Core of Fermentation

The increasing concerns of greenhouse gas emissions have increased the interest in dark fermentation as a means of productions for industrial chemicals, especially from renewable cellulosic biomass. However, the metabolism, including glycolysis, of many candidate organisms for cellulosic biomass conversion through consolidated bioprocessing is still poorly understood and the genomes have only recently been sequenced. Because a variety of industrial chemicals are produced directly from sugar metabolism, the careful understanding of glycolysis from a genomic and biochemical point of view is essential in the development of strategies for increasing product yields and therefore increasing industrial potential. The current review discusses the different pathways available for glycolysis along with unexpected variations from traditional models, especially in the utilization of alternate energy intermediates (GTP, pyrophosphate). This reinforces the need for a careful description of interactions between energy metabolites and glycolysis enzymes for understanding carbon and electron flux regulation.

Permanent draft genome sequence of Desulfurococcus mobilis type strain DSM 2161, a thermoacidophilic sulfur-reducing crenarchaeon isolated from acidic hot springs of Hveravellir, Iceland

This report presents the permanent draft genome sequence of Desulfurococcus mobilis type strain DSM 2161, an obligate anaerobic hyperthermophilic crenarchaeon that was isolated from acidic hot springs in Hveravellir, Iceland. D. mobilis utilizes peptides as carbon and energy sources and reduces elemental sulfur to H2S. A metabolic construction derived from the draft genome identified putative pathways for peptide degradation and sulfur respiration in this archaeon. Existence of several hydrogenase genes in the genome supported previous findings that H2 is produced during the growth of D. mobilis in the absence of sulfur. Interestingly, genes encoding glucose transport and utilization systems also exist in the D. mobilis genome though this archaeon does not utilize carbohydrate for growth. The draft genome of D. mobilis provides an additional mean for comparative genomic analysis of desulfurococci. In addition, our analysis on the Average Nucleotide Identity between D. mobilis and Desulfurococcus mucosus suggested that these two desulfurococci are two different strains of the same species.

New insights on the mechanism of the K(+-) independent activity of crenarchaeota pyruvate kinases

Eukarya pyruvate kinases have glutamate at position 117 (numbered according to the rabbit muscle enzyme), whereas in Bacteria have either glutamate or lysine and in Archaea have other residues. Glutamate at this position makes pyruvate kinases K+-dependent, whereas lysine confers K+-independence because the positively charged residue substitutes for the monovalent cation charge. Interestingly, pyruvate kinases from two characterized Crenarchaeota exhibit K+-independent activity, despite having serine at the equivalent position. To better understand pyruvate kinase catalytic activity in the absence of K+ or an internal positive charge, the Thermofilum pendens pyruvate kinase (valine at the equivalent position) was characterized. The enzyme activity was K+-independent. The kinetic mechanism was random order with a rapid equilibrium, which is equal to the mechanism of the rabbit muscle enzyme in the presence of K+ or the mutant E117K in the absence of K+. Thus, the substrate binding order of the T. pendens enzyme was independent despite lacking an internal positive charge. Thermal stability studies of this enzyme showed two calorimetric transitions, one attributable to the A and C domains (Tm of 99.2°C), and the other (Tm of 105.2°C) associated with the B domain. In contrast, the rabbit muscle enzyme exhibits a single calorimetric transition (Tm of 65.2°C). The calorimetric and kinetic data indicate that the B domain of this hyperthermophilic enzyme is more stable than the rest of the protein with a conformation that induces the catalytic readiness of the enzyme. B domain interactions of pyruvate kinases that have been determined in Pyrobaculum aerophilum and modeled in T. pendens were compared with those of the rabbit muscle enzyme. The results show that intra- and interdomain interactions of the Crenarchaeota enzymes may account for their higher B domain stability. Thus the structural arrangement of the T. pendens pyruvate kinase could allow charge-independent catalysis.

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.

An uncharacterized member of the ribokinase family in Thermococcus kodakarensis exhibits myo-inositol kinase activity

Here we performed structural and biochemical analyses on the TK2285 gene product, an uncharacterized protein annotated as a member of the ribokinase family, from the hyperthermophilic archaeon Thermococcus kodakarensis. The three-dimensional structure of the TK2285 protein resembled those of previously characterized members of the ribokinase family including ribokinase, adenosine kinase, and phosphofructokinase. Conserved residues characteristic of this protein family were located in a cleft of the TK2285 protein as in other members whose structures have been determined. We thus examined the kinase activity of the TK2285 protein toward various sugars recognized by well characterized ribokinase family members. Although activity with sugar phosphates and nucleosides was not detected, kinase activity was observed toward d-allose, d-lyxose, d-tagatose, d-talose, d-xylose, and d-xylulose. Kinetic analyses with the six sugar substrates revealed high Km values, suggesting that they were not the true physiological substrates. By examining activity toward amino sugars, sugar alcohols, and disaccharides, we found that the TK2285 protein exhibited prominent kinase activity toward myo-inositol. Kinetic analyses with myo-inositol revealed a greater kcat and much lower Km value than those obtained with the monosaccharides, resulting in over a 2,000-fold increase in kcat/Km values. TK2285 homologs are distributed among members of Thermococcales, and in most species, the gene is positioned close to a myo-inositol monophosphate synthase gene. Our results suggest the presence of a novel subfamily of the ribokinase family whose members are present in Archaea and recognize myo-inositol as a substrate.

The complete genome sequence of Thermoproteus tenax: a physiologically versatile member of the Crenarchaeota

Here, we report on the complete genome sequence of the hyperthermophilic Crenarchaeum Thermoproteus tenax (strain Kra1, DSM 2078^T) a type strain of the crenarchaeotal order Thermoproteales. Its circular 1.84-megabase genome harbors no extrachromosomal elements and 2,051 open reading frames are identified, covering 90.6% of the complete sequence, which represents a high coding density. Derived from the gene content, T. tenax is a representative member of the Crenarchaeota. The organism is strictly anaerobic and sulfur-dependent with optimal growth at 86°C and pH 5.6. One particular feature is the great metabolic versatility, which is not accompanied by a distinct increase of genome size or information density as compared to other Crenarchaeota. T. tenax is able to grow chemolithoautotrophically (CO₂/H₂) as well as chemoorganoheterotrophically in presence of various organic substrates. All pathways for synthesizing the 20 proteinogenic amino acids are present. In addition, two presumably complete gene sets for NADH:quinone oxidoreductase (complex I) were identified in the genome and there is evidence that either NADH or reduced ferredoxin might serve as electron donor. Beside the typical archaeal A₀A₁-ATP synthase, a membrane-bound pyrophosphatase is found, which might contribute to energy conservation. Surprisingly, all genes required for dissimilatory sulfate reduction are present, which is confirmed by growth experiments. Mentionable is furthermore, the presence of two proteins (ParA family ATPase, actin-like protein) that might be involved in cell division in Thermoproteales, where the ESCRT system is absent, and of genes involved in genetic competence (DprA, ComF) that is so far unique within Archaea.

Design of a coupled bioluminescent assay for a recombinant pyruvate kinase from a thermophilic Geobacillus

A simple and rapid method using coupled bioluminescent assay was developed to determine level of ADP. ADP is involved in many biological reactions and ADP assay can be used for assaying some reactions universally by monitoring ADP formation or depletion. ADP analysis involves incubation of ADP or extracts containing ADP with pyruvate kinase (PK) and PEP. The ATP formed by this reaction is determined by measuring the intensity of the initial light flash produced when luciferin-luciferase preparation injected into the reaction mixture. In regard to the main role of the PK in this assay, the gene of PK from a Geobacillus species has been cloned in expression vector pET28a (+), sequenced and overexpressed in Escherichia coli. Recombinant protein was purified using Ni-NTA column and then the purified PK was used in a coupled bioluminescent assay for ADP measurement. Kinetic properties of PK are determined according to a bioluminescent assay using firefly luciferase.

Crystal structure of pyruvate kinase from Geobacillus stearothermophilus

The pyruvate kinase (PK) from a moderate thermophile, Geobacillus stearothermophilus, is an allosteric enzyme activated by AMP and ribose 5-phosphate but not fructose 1, 6-bisphosphate (FBP), which is a common activator of PKs. It has an extra C-terminal sequence (ECTS), which contains a highly conserved phosphoenolpyruvate (PEP) binding motif, but its function and structure remain unclear. To elucidate the structural characteristics of the effector-binding site and the ECTS, the crystal structure of the C9S/C268S mutant of the enzyme was determined at 2.4 A resolution. The crystal belonged to space group P6(2)22, with unit cell parameters a, b = 145.97 A, c = 118.03 A. The enzyme was a homotetramer and its overall domain structure was similar to the previously solved structures except that the ECTS formed a new domain (C' domain). The structure of the C' domain closely resembled that of the PEP binding domain of maize pyruvate phosphate dikinase. A sulphate ion was found in a pocket in the effector-binding C domain. This site corresponds to the 6-phosphate group-binding site in yeast PK bound FBP and seems to be the effector-binding site. Through comparison of the structure of the putative effector-binding site to that of the FBP binding site of the yeast enzyme, the structural basis of the effector specificity of the G. stearothermophilus PK is discussed.

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.

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.

A heteromeric plastidic pyruvate kinase complex involved in seed oil biosynthesis in Arabidopsis

Glycolysis is a ubiquitous pathway thought to be essential for the production of oil in developing seeds of Arabidopsis thaliana and oil crops. Compartmentation of primary metabolism in developing embryos poses a significant challenge for testing this hypothesis and for the engineering of seed biomass production. It also raises the question whether there is a preferred route of carbon from imported photosynthate to seed oil in the embryo. Plastidic pyruvate kinase catalyzes a highly regulated, ATP-producing reaction of glycolysis. The Arabidopsis genome encodes 14 putative isoforms of pyruvate kinases. Three genes encode subunits alpha, beta(1), and beta(2) of plastidic pyruvate kinase. The plastid enzyme prevalent in developing seeds likely has a subunit composition of 4alpha4beta(1), is most active at pH 8.0, and is inhibited by Glu. Disruption of the gene encoding the beta(1) subunit causes a reduction in plastidic pyruvate kinase activity and 60% reduction in seed oil content. The seed oil phenotype is fully restored by expression of the beta(1) subunit-encoding cDNA and partially by the beta(2) subunit-encoding cDNA. Therefore, the identified pyruvate kinase catalyzes a crucial step in the conversion of photosynthate into oil, suggesting a preferred plastid route from its substrate phosphoenolpyruvate to fatty acids.

The glycolytic pathway of Trimastix pyriformis is an evolutionary mosaic

Background

Glycolysis and subsequent fermentation is the main energy source for many anaerobic organisms. The glycolytic pathway consists of ten enzymatic steps which appear to be universal amongst eukaryotes. However, it has been shown that the origins of these enzymes in specific eukaryote lineages can differ, and sometimes involve lateral gene transfer events. We have conducted an expressed sequence tag (EST) survey of the anaerobic flagellate Trimastix pyriformis to investigate the nature of the evolutionary origins of the glycolytic enzymes in this relatively unstudied organism.

Results

We have found genes in the Trimastix EST data that encode enzymes potentially catalyzing nine of the ten steps of the glycolytic conversion of glucose to pyruvate. Furthermore, we have found two different enzymes that in principle could catalyze the conversion of phosphoenol pyruvate (PEP) to pyruvate (or the reverse reaction) as part of the last step in glycolysis. Our phylogenetic analyses of all of these enzymes revealed at least four cases where the relationship of the Trimastix genes to homologs from other species is at odds with accepted organismal relationships. Although lateral gene transfer events likely account for these anomalies, with the data at hand we were not able to establish with confidence the bacterial donor lineage that gave rise to the respective Trimastix enzymes.

Conclusion

A number of the glycolytic enzymes of Trimastix have been transferred laterally from bacteria instead of being inherited from the last common eukaryotic ancestor. Thus, despite widespread conservation of the glycolytic biochemical pathway across eukaryote diversity, in a number of protist lineages the enzymatic components of the pathway have been replaced by lateral gene transfer from disparate evolutionary sources. It remains unclear if these replacements result from selectively advantageous properties of the introduced enzymes or if they are neutral outcomes of a gene transfer 'ratchet' from food or endosymbiotic organisms or a combination of both processes.

Cell-free protein synthesis at high temperatures using the lysate of a hyperthermophile

Systems for cell-free protein synthesis available today are usually based on the lysates of either Escherichia coli, wheat germ or rabbit reticulocyte, and protein synthesis reactions using these extracts are performed at moderate temperatures (20-40 degrees C). We report here the development of a novel system for cell-free protein synthesis that can be operated at high temperatures using a lysate of the hyperthermophilic archaeon, Thermococcus kodakaraensis. With the system, cell-free protein synthesis of ChiADelta4, a derivative of T. kodakaraensis chitinase (ChiA), was observed within a temperature range of 40-80 degrees C, with an optimum at 65 degrees C. Corresponding chitinase activity was also detected in the reaction mixtures after cell-free protein synthesis, indicating that the synthesized ChiADelta4 folded in a proper tertiary structure. The maximum concentration of active ChiADelta4 synthesized was determined to be approximately 1.3 microg/mL. A time course experiment indicated that the amount of synthesized ChiADelta4 saturated within 30 min at 65 degrees C, and energy depletion was suggested to be the main cause of this saturation. We further developed a system for transcription and translation-coupled protein synthesis at high temperatures using a combination of T. kodakaraensis lysate and thermostable T7 RNA polymerase.

Phosphoenolpyruvate synthase plays an essential role for glycolysis in the modified Embden-Meyerhof pathway in Thermococcus kodakarensis

We have carried out a genetic analysis on pyruvate kinase (PykTk) and phosphoenolpyruvate synthase (PpsTk) in the hyperthermophilic archaeon, Thermococcus kodakarensis. In principle, both enzymes can catalyse the final step of the modified Embden-Meyerhof (EM) pathway found in Thermococcales, the conversion of phosphoenolpyruvate (PEP) to pyruvate, with the former utilizing ADP, while the latter is dependent on AMP and phosphate. Enzyme activities and transcript levels of both PykTk and PpsTk increased in T. kodakarensis under glycolytic conditions when compared with cells grown on pyruvate or amino acids. Using KW128, a tryptophan auxotrophic mutant with a trpE gene disruption, as a host strain, we obtained mutant strains with single gene disruptions in either the pykTk (Deltapyk strain) or ppsTk (Deltapps strain) gene. Specific growth rates and cell yields were examined in various media and compared with the host KW128 strain. The results indicated that both enzymes participated in pyruvate metabolism, but were not essential. In the presence of maltooligosaccharides, the Deltapyk strain displayed a 15% decrease in growth rate compared with the host strain, indicating that PykTk does participate in glycolysis. However an even more dramatic effect was observed in the Deltapps strain in that the strain could not grow at all on maltooligosaccharides. The results clearly indicate that, in contrast to the conventional EM pathway dependent on pyruvate kinase, PEP synthase is the essential enzyme for the glycolytic conversion of PEP to pyruvate in T. kodakarensis. The physiological roles of the two enzymes under various growth conditions are discussed.

Dichotomic Phylogenetic Tree of the Pyruvate Kinase Family

K+ dependence was assumed to be a feature of all pyruvate kinases until it was discovered that some enzymes express K+ -independent activity. Almost all the K+ -independent pyruvate kinases have Lys at position 117, instead of the Glu present in the K+ -dependent muscle enzyme. Mutagenesis studies show that the internal positive charge substitutes for the K+ requirement (Laughlin, L. T. & Reed, G. H. (1997) Arch. Biochem. Biophys. 348, 262-267). In this work a phylogenetic analysis of pyruvate kinase was performed to ascertain the abundance of K+ -independent activities and to explore whether the K+ activating effect is related to the evolutionary history of the enzyme. Of the 230 studied sequences, 46% have Lys at position 117, and the rest have Glu. Pyruvate kinases with Lys117 and Glu117 are separated in two clusters. All of the enzymes of the Glu117 cluster that have been characterized are K+ -dependent, whereas those of the Lys117 cluster are K+ -independent. Thus, there is a strict correlation between the dichotomy of the tree and the dependence of activity on K+. 77% of the pyruvate kinases that possess Lys117 have Lys113/Gln114; they also have Ile, Val, or Leu at position 120. These residues are replaced by Glu117 and Thr113/Lys114/Thr120 in 80% of K+ -dependent pyruvate kinases. Structural analysis indicates that these residues are in a hinge region involved in the acquisition of the catalytic conformation of the enzyme. The route of conversion from K+ -independent to K+ -dependent pyruvate kinases is described. A plausible explanation of how enzymes developed K+ dependence is put forth.

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.

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.

The pyruvate kinase of Stigmatella aurantiaca is an indole binding protein and essential for development

Myxospore formation of the myxobacterium Stigmatella aurantiaca can be uncoupled from the cooperative development i.e. fruiting body formation, by low concentrations of indole. Two putative indole receptor proteins were isolated by their capacity to bind indole and identified as pyruvate kinase (PK) and aldehyde dehydrogenase. The PK activity of Stigmatella crude extracts was stimulated by indole. Cloning of the PK gene (pykA) and the construction of a pykA disruption mutant strikingly revealed that PK is essential for multicellular development: Fruiting body formation was abolished in the mutant strain and indole-induced spore formation was delayed. The developmental defects could be complemented by insertion of the pykA gene at the mtaB locus of the Stigmatella genome excluding any polar effects of the pykA disruption.

Pyruvate kinase revisited: the activating effect of K+

For more than 50 years, it has been known that K(+) is an essential activator of pyruvate kinase (Kachmar, J. F., and Boyer, P. D. (1953) J. Biol. Chem. 200, 669-683). However, the role of K(+) in the catalysis by pyruvate kinase has not been totally understood. Previous studies without K(+) showed that the affinity of ADP-Mg(2+) depends on the concentration of phosphoenolpyruvate, although the kinetics of the enzyme at saturating K(+) concentrations show independence in the binding of substrates (Reynard, A. M., Hass, L. F., Jacobsen, D. D. & Boyer, P. D. (1961) J. Biol. Chem. 236, 2277-2283). Here, we explored the kinetics of the enzyme with and without K(+). The results show that without K(+), the kinetic mechanism of pyruvate kinase changes from random to ordered with phosphoenol-pyruvate as first substrate. V(max) with K(+) was about 400 higher than without K(+). In the presence of K(+), the affinities for phosphoenol-pyruvate, ADP-Mg(2+), oxalate, and ADP-Cr(2+) were 2-6-fold higher than in the absence of K(+). This as well as fluorescence data also indicate that K(+) is involved in the acquisition of the active conformation of the enzyme, allowing either phosphoenolpyruvate or ADP to bind independently (random mechanism). In the absence of K(+), ADP cannot bind to the enzyme until phosphoenolpyruvate forms a competent active site (ordered mechanism). We propose that K(+) induces the closure of the active site and the arrangement of the residues involved in the binding of the nucleotide.

Crystallization and preliminary X-ray analysis of pyruvate kinase from Bacillus stearothermophilus

Pyruvate kinase (PK) from a moderate thermophile, Bacillus stearothermophilus (BstPK), is an allosteric enzyme activated by AMP and ribose 5-phosphate but not by fructose 1,6-bisphosphate (FBP). However, almost all other PKs are activated by FBP. The wild-type and W416F/V435W mutant BstPKs were crystallized by the hanging-drop vapour-diffusion method. However, they were unsuitable for structural analysis because their data sets exhibited low completeness. A crystal suitable for structural analysis was obtained using C9S/C268S enzyme. The crystal belonged to space group P6(2)22, with unit-cell parameters a = b = 145.97, c = 118.03 A.

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.

Entamoeba histolytica: kinetic and molecular evidence of a previously unidentified pyruvate kinase

We report the kinetic characterization of a previously unidentified pyruvate kinase (PK) activity in extracts from Entamoeba histolytica trophozoites. This activity was about 74% of the activity of pyruvate phosphate dikinase. EhPK differed from most PKs in that its pH optimum was 5.5-6.5 and was inhibited by high PEP concentrations (1-5mM); these are concentrations at which PK is usually assayed. The optimal temperature was above 40 degrees C with negligible activity below 20 degrees C. EhPK exhibited hyperbolic kinetics with respect to both PEP (K(m) = 0.018 mM) and ADP (K(m) = 1.05 mM). However, it exhibited a sigmoidal behavior with respect to PEP at sub-saturating ADP concentrations. EhPK did not require monovalent cations for activity. Fructose-1,6 bisphosphate was a potent non-essential activator; it increased the affinity for ADP without modification of the V(max) or the affinity for PEP. Phosphate, citrate, malate, and alpha-ketoglutarate significantly inhibited EhPK activity. A putative EhPK gene fragment found in EhDNA was analyzed. The data indicate that E. histolytica trophozoites contain an active PK, which might contribute to the generation of glycolytic ATP for parasite survival.

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.

The unique features of glycolytic pathways in Archaea

An early divergence in evolution has resulted in two prokaryotic domains, the Bacteria and the Archaea. Whereas the central metabolic routes of bacteria and eukaryotes are generally well-conserved, variant pathways have developed in Archaea involving several novel enzymes with a distinct control. A spectacular example of convergent evolution concerns the glucose-degrading pathways of saccharolytic archaea. The identification, characterization and comparison of the glycolytic enzymes of a variety of phylogenetic lineages have revealed a mosaic of canonical and novel enzymes in the archaeal variants of the Embden-Meyerhof and the Entner-Doudoroff pathways. By means of integrating results from biochemical and genetic studies with recently obtained comparative and functional genomics data, the structure and function of the archaeal glycolytic routes, the participating enzymes and their regulation are re-evaluated.

Bifunctional aldehyde/alcohol dehydrogenase (ADHE) in chlorophyte algal mitochondria

Protein profiles of mitochondria isolated from the heterotrophic chlorophyte Polytomella sp. grown on ethanol at pH 6.0 and pH 3.7 were analyzed by Blue Native and denaturing polyacrylamide gel electrophoresis. Steady-state levels of oxidative phosphorylation complexes were influenced by external pH. Levels of an abundant, soluble, mitochondrial protein of 85 kDa and its corresponding mRNA increased at pH 6.0 relative to pH 3.7. N-terminal and internal sequencing of the 85 kDa mitochondrial protein together with the corresponding cDNA identified it as a bifunctional aldehyde/alcohol dehydrogenase (ADHE) with strong similarity to homologues from eubacteria and amitochondriate protists. A mitochondrial targeting sequence of 27 amino acids precedes the N-terminus of the mature mitochondrial protein. A gene encoding an ADHE homologue was also identified in the genome of Chlamydomonas reinhardtii, a photosynthetic relative of Polytomella. ADHE reveals a complex picture of sequence similarity among homologues. The lack of ADHE from archaebacteria indicates a eubacterial origin for the eukaryotic enzyme. Among eukaryotes, ADHE has hitherto been characteristic of anaerobes since it is essential to cytosolic energy metabolism of amitochondriate protists such as Giardia intestinalis and Entamoeba histolytica. Its abundance and expression pattern suggest an important role for ADHE in mitochondrial metabolism of Polytomella under the conditions studied. The current data are compatible with the view that Polytomella ADHE could be involved either in ethanol production or assimilation, or both, depending upon environmental conditions. Presence of ADHE in an oxygen-respiring algal mitochondrion and co-expression at ambient oxygen levels with respiratory chain components is unexpected with respect to the view that eukaryotes acquired ADHE genes specifically as an adaptation to an anaerobic lifestyle.

Comparative analysis of pyruvate kinases from the hyperthermophilic archaea Archaeoglobus fulgidus, Aeropyrum pernix, and Pyrobaculum aerophilum and the hyperthermophilic bacterium Thermotoga maritima: unusual regulatory properties in hyperthermophilic archaea

Pyruvate kinases (PK, EC 2.7.1.40) from three hyperthermophilic archaea (Archaeoglobus fulgidus strain 7324, Aeropyrum pernix, and Pyrobaculum aerophilum) and from the hyperthermophilic bacterium Thermotoga maritima were compared with respect to their thermophilic, kinetic, and regulatory properties. PKs from the archaea are 200-kDa homotetramers composed of 50-kDa subunits. The enzymes required divalent cations, Mg2+ and Mn2+ being most effective, but were independent of K+. Temperature optima for activity were 85 degrees C (A. fulgidus) and above 98 degrees C (A. pernix and P. aerophilum). The PKs were highly thermostable up to 110 degrees C (A. pernix) and showed melting temperatures for thermal unfolding at 93 degrees C (A. fulgidus) or above 98 degrees C (A. pernix and P. aerophilum). All archaeal PKs exhibited sigmoidal saturation kinetics with phosphoenolpyruvate (PEP) and ADP indicating positive homotropic cooperative response with both substrates. Classic heterotropic allosteric regulators of PKs from eukarya and bacteria, e.g. fructose 1,6-bisphosphate or AMP, did not affect PK activity of hyperthermophilic archaea, suggesting the absence of heterotropic allosteric regulation. PK from the bacterium T. maritima is also a homotetramer of 50-kDa subunits. The enzyme was independent of K+ ions, had a temperature optimum of 80 degrees C, was highly thermostable up to 90 degrees C, and had a melting temperature above 98 degrees C. The enzyme showed cooperative response to PEP and ADP. In contrast to its archaeal counterparts, the T. maritima enzyme exhibited the classic allosteric response to the activator AMP and to the inhibitor ATP. Sequences of hyperthermophilic PKs showed significant similarity to characterized PKs from bacteria and eukarya. Phylogenetic analysis of PK sequences of all three domains indicates a distinct archaeal cluster that includes the PK from the hyperthermophilic bacterium T. maritima.

Pyruvate kinase: current status of regulatory and functional properties

Pyruvate kinase (PK) is a key enzyme for the glycolytic pathway and carbon metabolism in general. On the basis of the relevance and enormous diverse properties of this enzyme, this paper describes the results of a current and extensive review that determines the sites of conservation and/or difference in PK sequences, and the differences in the functional and regulatory properties of the enzymes. An alignment and analysis of 50 PK sequences from different sources and a phylogenetic tree are presented. This analysis was performed with reference to crystallographically characterized PK principally from E. coli, cat and rabbit muscle. A number of attributes of the enzyme that make it of particular interest in biomedicine and industry are also discussed.