The Control of Pyruvate Kinases of Escherichia coli

Structural Basis of Nucleotide Selectivity in Pyruvate Kinase

Nucleoside triphosphates are indispensable in numerous biological processes, with enzymes involved in their biogenesis playing pivotal roles in cell proliferation. Pyruvate kinase (PYK), commonly regarded as the terminal glycolytic enzyme that generates ATP in tandem with pyruvate, is also capable of synthesizing a wide range of nucleoside triphosphates from their diphosphate precursors. Despite their substrate promiscuity, some PYKs show preference towards specific nucleotides, suggesting an underlying mechanism for differentiating nucleotide bases. However, the thorough characterization of this mechanism has been hindered by the paucity of nucleotide-bound PYK structures. Here, we present crystal structures of Streptococcus pneumoniae PYK in complex with four different nucleotides. These structures facilitate direct comparison of the protein-nucleotide interactions and offer structural insights into its pronounced selectivity for GTP synthesis. Notably, this selectivity is dependent on a sequence motif in the nucleotide recognition site that is widely present among prokaryotic PYKs, particularly in Firmicutes species. We show that pneumococcal cell growth is significantly impaired when expressing a PYK variant with compromised GTP and UTP synthesis activity, underscoring the importance of PYK in maintaining nucleotide homeostasis. Our findings collectively advance our understanding of PYK biochemistry and prokaryotic metabolism.

Pyruvate kinase 2 from Synechocystis sp. PCC 6803 increased substrate affinity via glucose-6-phosphate and ribose-5-phosphate for phosphoenolpyruvate consumption

Pyruvate kinase (Pyk, EC 2.7.1.40) is a glycolytic enzyme that generates pyruvate and adenosine triphosphate (ATP) from phosphoenolpyruvate (PEP) and adenosine diphosphate (ADP), respectively. Pyk couples pyruvate and tricarboxylic acid metabolisms. Synechocystis sp. PCC 6803 possesses two pyk genes (encoded pyk1, sll0587 and pyk2, sll1275). A previous study suggested that pyk2 and not pyk1 is essential for cell viability; however, its biochemical analysis is yet to be performed. Herein, we biochemically analyzed Synechocystis Pyk2 (hereafter, SyPyk2). The optimum pH and temperature of SyPyk2 were 7.0 and 55 °C, respectively, and the Km values for PEP and ADP under optimal conditions were 1.5 and 0.053 mM, respectively. SyPyk2 is activated in the presence of glucose-6-phosphate (G6P) and ribose-5-phosphate (R5P); however, it remains unaltered in the presence of adenosine monophosphate (AMP) or fructose-1,6-bisphosphate. These results indicate that SyPyk2 is classified as PykA type rather than PykF, stimulated by sugar monophosphates, such as G6P and R5P, but not by AMP. SyPyk2, considering substrate affinity and effectors, can play pivotal roles in sugar catabolism under nonphotosynthetic conditions.

The 2.4 Å structure of Zymomonas mobilis pyruvate kinase: Implications for stability and regulation

Human liver pyruvate kinase (hlPYK) catalyzes the final step in glycolysis, the formation of pyruvate (PYR) and ATP from phosphoenolpyruvate (PEP) and ADP. Fructose 1,6-bisphosphate (FBP), a pathway intermediate of glycolysis, serves as an allosteric activator of hlPYK. Zymomonas mobilis pyruvate kinase (ZmPYK) performs the final step of the Entner-Doudoroff pathway, which is similar to glycolysis in that energy is harvested from glucose and pyruvate is generated. The Entner-Doudoroff pathway does not have FBP as a pathway intermediate, and ZmPYK is not allosterically activated. In this work, we solved the 2.4 Å X-ray crystallographic structure of ZmPYK. The protein is dimeric in solution as determined by gel filtration chromatography, but crystallizes as a tetramer. The buried surface area of the ZmPYK tetramerization interface is significantly smaller than that of hlPYK, and yet tetramerization using the standard interfaces from higher organisms provides an accessible low energy crystallization pathway. Interestingly, the ZmPYK structure showed a phosphate ion in the analogous location to the 6-phosphate binding site of FBP in hlPYK. Circular Dichroism (CD) was used to measure melting temperatures of hlPYK and ZmPYK in the absence and presence of substrates and effectors. The only significant difference was an additional phase of small amplitude for the ZmPYK melting curves. We conclude that the phosphate ion plays neither a structural or allosteric role in ZmPYK under the conditions tested. We hypothesize that ZmPYK does not have sufficient protein stability for activity to be tuned by allosteric effectors as described for rheostat positions in the allosteric homologues.

Functional and structural characterization of Streptococcus pneumoniae pyruvate kinase involved in fosfomycin resistance

Glycolysis is the primary metabolic pathway in the strictly fermentative Streptococcus pneumoniae, which is a major human pathogen associated with antibiotic resistance. Pyruvate kinase (PYK) is the last enzyme in this pathway that catalyzes the production of pyruvate from phosphoenolpyruvate (PEP) and plays a crucial role in controlling carbon flux; however, while S. pneumoniae PYK (SpPYK) is indispensable for growth, surprisingly little is known about its functional properties. Here, we report that compromising mutations in SpPYK confer resistance to the antibiotic fosfomycin, which inhibits the peptidoglycan synthesis enzyme MurA, implying a direct link between PYK and cell wall biogenesis. The crystal structures of SpPYK in the apo and ligand-bound states reveal key interactions that contribute to its conformational change as well as residues responsible for the recognition of PEP and the allosteric activator fructose 1,6-bisphosphate (FBP). Strikingly, FBP binding was observed at a location distinct from previously reported PYK effector binding sites. Furthermore, we show that SpPYK could be engineered to become more responsive to glucose 6-phosphate instead of FBP by sequence and structure-guided mutagenesis of the effector binding site. Together, our work sheds light on the regulatory mechanism of SpPYK and lays the groundwork for antibiotic development that targets this essential enzyme.

Rewiring cell-free metabolic flux in E. coli lysates using a block-push-pull approach

Cell-free systems can expedite the design and implementation of biomanufacturing processes by bypassing troublesome requirements associated with the use of live cells. In particular, the lack of survival objectives and the open nature of cell-free reactions afford engineering approaches that allow purposeful direction of metabolic flux. The use of lysate-based systems to produce desired small molecules can result in competitive titers and productivities when compared to their cell-based counterparts. However, pathway crosstalk within endogenous lysate metabolism can compromise conversion yields by diverting carbon flow away from desired products. Here, the 'block-push-pull' concept of conventional cell-based metabolic engineering was adapted to develop a cell-free approach that efficiently directs carbon flow in lysates from glucose and toward endogenous ethanol synthesis. The approach is readily adaptable, is relatively rapid and allows for the manipulation of central metabolism in cell extracts. In implementing this approach, a block strategy is first optimized, enabling selective enzyme removal from the lysate to the point of eliminating by-product-forming activity while channeling flux through the target pathway. This is complemented with cell-free metabolic engineering methods that manipulate the lysate proteome and reaction environment to push through bottlenecks and pull flux toward ethanol. The approach incorporating these block, push and pull strategies maximized the glucose-to-ethanol conversion in an Escherichia coli lysate that initially had low ethanologenic potential. A 10-fold improvement in the percent yield is demonstrated. To our knowledge, this is the first report of successfully rewiring lysate carbon flux without source strain optimization and completely transforming the consumed input substrate to a desired output product in a lysate-based, cell-free system.

Glucose Phosphotransferase System Modulates Pyruvate Metabolism, Bacterial Fitness, and Microbial Ecology in Oral Streptococci

Spontaneous mutants with defects in the primary glucose phosphotransferase permease (manLMNO) of Streptococcus sanguinis SK36 showed enhanced fitness at low pH. Transcriptomics and metabolomics with a manL deletion mutant (SK36/manL) revealed redirection of pyruvate to production of acetate and formate, rather than lactate. These observations were consistent with measurements of decreased lactic acid accumulation and increased excretion of acetate, formate, pyruvate, and H2O2. Genes showing increased expression in SK36/manL included those encoding carbohydrate transporters, extracellular glycosidases, intracellular polysaccharide metabolism, and arginine deiminase and pathways for metabolism of acetoin, ethanolamine, ascorbate, and formate, along with genes required for membrane biosynthesis and adhesion. Streptococcus mutans UA159 persisted much better in biofilm cocultures with SK36/manL than with SK36, an effect that was further enhanced by culturing the biofilms anaerobically but dampened by adding arginine to the medium. We posited that the enhanced persistence of S. mutans with SK36/manL was in part due to excess excretion of pyruvate by the latter, as addition of pyruvate to S. mutans-S. sanguinis cocultures increased the proportions of UA159 in the biofilms. Reducing the buffer capacity or increasing the concentration of glucose benefited UA159 when cocultured with SK36, but not with SK36/manL, likely due to the altered metabolism and enhanced acid tolerance of the mutant. When manL was deleted in S. mutans or Streptococcus gordonii, the mutants presented altered fitness characteristics. Our study demonstrated that phosphotransferase system (PTS)-dependent modulation of central metabolism can profoundly affect streptococcal fitness and metabolic interactions, revealing another dimension in commensal-pathogen relationships influencing dental caries development. IMPORTANCE Dental caries is underpinned by a dysbiotic microbiome and increased acid production. As beneficial bacteria that can antagonize oral pathobionts, oral streptococci such as S. sanguinis and S. gordonii can ferment many carbohydrates, despite their relative sensitivity to low pH. We characterized the molecular basis for why mutants of glucose transporter ManLMNO of S. sanguinis showed enhanced production of hydrogen peroxide and ammonia and improved persistence under acidic conditions. A metabolic shift involving more than 300 genes required for carbohydrate transport, energy production, and envelope biogenesis was observed. Significantly, manL mutants engineered in three different oral streptococci displayed altered capacities for acid production and interspecies antagonism, highlighting the potential for targeting the glucose-PTS to modulate the pathogenicity of oral biofilms.

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.

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.

Structure, Function and Regulation of a Second Pyruvate Kinase Isozyme in Pseudomonas aeruginosa

Pseudomonas aeruginosa (PA) depends on the Entner-Doudoroff pathway (EDP) for glycolysis. The main enzymatic regulator in the lower half of the EDP is pyruvate kinase. PA contains genes that encode two isoforms of pyruvate kinase, denoted PykAPA and PykFPA. In other well-characterized organisms containing two pyruvate kinase isoforms (such as Escherichia coli) each isozyme is differentially regulated. The structure, function and regulation of PykAPA has been previously characterized in detail, so in this work, we set out to assess the biochemical and structural properties of the PykFPA isozyme. We show that pykF PA expression is induced in the presence of the diureide, allantoin. In spite of their relatively low amino acid sequence identity, PykAPA and PykFPA display broadly comparable kinetic parameters, and are allosterically regulated by a very similar set of metabolites. However, the x-ray crystal structure of PykFPA revealed significant differences compared with PykAPA. Notably, although the main allosteric regulator binding-site of PykFPA was empty, the "ring loop" covering the site adopted a partially closed conformation. Site-directed mutation of the proline residues flanking the ring loop yielded apparent "locked on" and "locked off" allosteric activation phenotypes, depending on the residue mutated. Analysis of PykFPA inter-protomer interactions supports a model in which the conformational transition(s) accompanying allosteric activation involve re-orientation of the A and B domains of the enzyme and subsequent closure of the active site.

Structural and kinetic characterization of Trypanosoma congolense pyruvate kinase

Trypanosoma are blood-borne parasites and are the causative agents of neglected tropical diseases (NTDs) affecting both humans and animals. These parasites mainly rely on glycolysis for their energy production within the mammalian host, which is why trypanosomal glycolytic enzymes have been pursued as interesting targets for the development of trypanocidal drugs. The structure-function relationships of pyruvate kinases (PYKs) from trypanosomatids (Trypanosoma and Leishmania) have been well-studied within this context. In this paper, we describe the structural and enzymatic characterization of PYK from T. congolense (TcoPYK), the main causative agent of Animal African Trypanosomosis (AAT), by employing a combination of enzymatic assays, thermal unfolding studies and X-ray crystallography.

Acidogenesis, solventogenesis, metabolic stress response and life cycle changes in Clostridium beijerinckii NRRL B-598 at the transcriptomic level

Clostridium beijerinckii NRRL B-598 is a sporulating, butanol and hydrogen producing strain that utilizes carbohydrates by the acetone-butanol-ethanol (ABE) fermentative pathway. The pathway consists of two metabolic phases, acidogenesis and solventogenesis, from which the latter one can be coupled with sporulation. Thorough transcriptomic profiling during a complete life cycle and both metabolic phases completed with flow cytometry, microscopy and a metabolites analysis helped to find out key genes involved in particular cellular events. The description of genes/operons that are closely involved in metabolism or the cell cycle is a necessary condition for metabolic engineering of the strain and will be valuable for all C. beijerinckii strains and other Clostridial species. The study focused on glucose transport and catabolism, hydrogen formation, metabolic stress response, binary fission, motility/chemotaxis and sporulation, which resulted in the composition of the unique image reflecting clostridial population changes. Surprisingly, the main change in expression of individual genes was coupled with the sporulation start and not with the transition from acidogenic to solventogenic metabolism. As expected, solvents formation started at pH decrease and the accumulation of butyric and acetic acids in the cultivation medium.

Allosteric pyruvate kinase-based "logic gate" synergistically senses energy and sugar levels in Mycobacterium tuberculosis

Pyruvate kinase (PYK) is an essential glycolytic enzyme that controls glycolytic flux and is critical for ATP production in all organisms, with tight regulation by multiple metabolites. Yet the allosteric mechanisms governing PYK activity in bacterial pathogens are poorly understood. Here we report biochemical, structural and metabolomic evidence that Mycobacterium tuberculosis (Mtb) PYK uses AMP and glucose-6-phosphate (G6P) as synergistic allosteric activators that function as a molecular "OR logic gate" to tightly regulate energy and glucose metabolism. G6P was found to bind to a previously unknown site adjacent to the canonical site for AMP. Kinetic data and structural network analysis further show that AMP and G6P work synergistically as allosteric activators. Importantly, metabolome profiling in the Mtb surrogate, Mycobacterium bovis BCG, reveals significant changes in AMP and G6P levels during nutrient deprivation, which provides insights into how a PYK OR gate would function during the stress of Mtb infection.

Chemoenzymatic synthesis of 3'-phosphoadenosine-5'-phosphosulfate coupling with an ATP regeneration system

3'-Phosphoadenosine-5'-phosphosulfate (PAPS) is the obligate cosubstrate and source of the sulfonate group in the chemoenzymatic synthesis of heparin, a commonly used anticoagulant drug. Previously, using ATP as the substrate, we had developed a one-pot synthesis to prepare PAPS with 47% ATP conversion efficiency. During the reaction, 47% of ATP was converted into the by-product, ADP. Here, to increase the conversion ratio of ATP to PAPS, an ATP regeneration system was developed to couple with PAPS synthesis. In the ATP regeneration system, the chemical compound, monopotassium phosphoenolpyruvate (PEP^-K⁺), was synthesized and used as the phospho-donor. By using 3-bromopyruvic acid as the starting material, the total yield of PEP^-K⁺ synthesis was over 50% at low cost. Then, the enzyme PykA from Escherichia coli was overexpressed, purified, and used to convert the by-product ADP into ATP. When coupled the ATP regeneration system with PAPS synthesis, the higher ratio of PEP^-K⁺ to ADP was associated with higher ATP conversion efficiency. By using the ATP regeneration system, the conversion ratio of ATP to PAPS was increased to 98% as determined by PAMN-HPLC analysis, and 5 g of PAPS was produced in 1 L of the reaction mixture. Furthermore, the chemoenzymatic synthesized PAPS was purified and freeze-dried without observed decomposition. However, the powdery PAPS was more unstable than the PAPS sodium salt in aqueous solution at ambient temperature. This developed chemoenzymatic approach of PAPS production will contribute to the synthesis of heparin, in which PAPS is necessary as the individual sulfo-donor.

The contribution of two isozymes to the pyruvate kinase activity of Vibrio cholerae: One K+-dependent constitutively active and another K+-independent with essential allosteric activation

In a previous phylogenetic study of the family of pyruvate kinase EC (2.7.1.40), a cluster with Glu117 and another with Lys117 were found (numbered according to the rabbit muscle enzyme). The sequences with Glu117 have been found to be K+-dependent, whereas those with Lys117 were K+-independent. Interestingly, only γ-proteobacteria exhibit sequences in both branches of the tree. In this context, it was explored whether these phylogenetically distinct pyruvate kinases were both expressed and contribute to the pyruvate kinase activity in Vibrio cholerae. The main findings of this work showed that the isozyme with Glu117 is an active K+-dependent enzyme. At the same substrate concentration, its Vmax in the absence of fructose 1,6 bisphosphate was 80% of that with its effector. This result is in accordance with the non-essential activation described by allosteric ligands for most pyruvate kinases. In contrast, the pyruvate kinase with Lys117 was a K+-independent enzyme displaying an allosteric activation by ribose 5-phosphate. At the same substrate concentration, its activity without the effector was 0.5% of the one obtained in the presence of ribose 5-phosphate, indicating that this sugar monophosphate is a strong activator of this enzyme. This absolute allosteric dependence is a novel feature of pyruvate kinase activity. Interestingly, in the K+-independent enzyme, Mn2+ may "mimic" the allosteric effect of Rib 5-P. Despite their different allosteric behavior, both isozymes display a rapid equilibrium random order kinetic mechanism. The intracellular concentrations of fructose 1,6-bisphosphate and ribose 5-phosphate in Vibrio cholerae have been experimentally verified to be sufficient to induce maximal activation of both enzymes. In addition, Western blot analysis indicated that both enzymes were co-expressed. Therefore, it is concluded that VcIPK and VcIIPK contribute to the activity of pyruvate kinase in this γ-proteobacterium.

Reexamination of the Physiological Role of PykA in Escherichia coli Revealed that It Negatively Regulates the Intracellular ATP Levels under Anaerobic Conditions

Pyruvate kinase is one of the three rate-limiting glycolytic enzymes that catalyze the last step of glycolysis, conversion of phosphoenolpyruvate (PEP) into pyruvate, which is associated with ATP generation. Two isozymes of pyruvate kinase, PykF and PykA, are identified in Escherichia coli PykF is considered important, whereas PykA has a less-defined role. Prior studies inactivated the pykA gene to increase the level of its substrate, PEP, and thereby increased the yield of end products derived from PEP. We were surprised when we found a pykA::Tn5 mutant in a screen for increased yield of an end product derived from pyruvate (n-butanol), suggesting that the role of PykA needs to be reexamined. We show that the pykA mutant exhibited elevated intracellular ATP levels, biomass concentrations, glucose consumption, and n-butanol production. We also discovered that the pykA mutant expresses higher levels of a presumed pyruvate transporter, YhjX, permitting the mutant to recapture and metabolize excreted pyruvate. Furthermore, we demonstrated that the nucleotide diphosphate kinase activity of PykA leads to negative regulation of the intracellular ATP levels. Taking the data together, we propose that inactivation of pykA can be considered a general strategy to enhance the production of pyruvate-derived metabolites under anaerobic conditions.IMPORTANCE This study showed that knocking out pykA significantly increased the intracellular ATP level and thus significantly increased the levels of glucose consumption, biomass formation, and pyruvate-derived product formation under anaerobic conditions. pykA was considered to be encoding a dispensable pyruvate kinase; here we show that pykA negatively regulates the anaerobic glycolysis rate through regulating the energy distribution. Thus, knocking out pykA can be used as a general strategy to increase the level of pyruvate-derived fermentative products.

In Vitro Reconstitution and Optimization of the Entire Pathway to Convert Glucose into Fatty Acid

Glucose and fatty acids play essential physiological roles in nearly all living organisms, and the pathway that converts glucose into fatty acid is pivotal to the central metabolic network. We have successfully reconstituted a pathway that converts glucose to fatty acid in vitro using 30 purified proteins. Through systematic titration and optimization of the glycolytic pathway and pyruvate dehydrogenase, we increased the yield of free fatty acid from nondetectable to a level that exceeded 9% of the theoretical yield. We also reconstituted the entire pentose-phosphate pathway of Escherichia coli and established a pentose phosphate-glycolysis hybrid pathway, replacing GAPDH to enhance NADPH availability. Our efforts provide a useful platform for research involving these core biochemical transformations.

A probabilistic framework for the exploration of enzymatic capabilities based on feasible kinetics and control analysis

BACKGROUND

Analysis of limiting steps within enzyme-catalyzed reactions is fundamental to understand their behavior and regulation. Methods capable of unravelling control properties and exploring kinetic capabilities of enzymatic reactions would be particularly useful for protein and metabolic engineering. While single-enzyme control analysis formalism has previously been applied to well-studied enzymatic mechanisms, broader application of this formalism is limited in practice by the limited amount of kinetic data and the difficulty of describing complex allosteric mechanisms.

METHODS

To overcome these limitations, we present here a probabilistic framework enabling control analysis of previously unexplored mechanisms under uncertainty. By combining a thermodynamically consistent parameterization with an efficient Sequential Monte Carlo sampler embedded in a Bayesian setting, this framework yields insights into the capabilities of enzyme-catalyzed reactions with modest kinetic information, provided that the catalytic mechanism and a thermodynamic reference point are defined.

RESULTS

The framework was used to unravel the impact of thermodynamic affinity, substrate saturation levels and effector concentrations on the flux control and response coefficients of a diverse set of enzymatic reactions.

CONCLUSIONS

Our results highlight the importance of the metabolic context in the control analysis of isolated enzymes as well as the use of statistically sound methods for their interpretation.

GENERAL SIGNIFICANCE

This framework significantly expands our current capabilities for unravelling the control properties of general reaction kinetics with limited amount of information. This framework will be useful for both theoreticians and experimentalists in the field.

Mitochondrial substrates in cancer: drivers or passengers?

The majority of cancers demonstrate various tumor-specific metabolic aberrations, such as increased glycolysis even under aerobic conditions (Warburg effect), whereas mitochondrial metabolic activity and their contribution to cellular energy production are restrained. One of the most important mechanisms for this metabolic switch is the alteration in the abundance, utilization, and localization of various mitochondrial substrates. Numerous lines of evidence connect disturbances in mitochondrial metabolic pathways with tumorigenesis and provide an intriguing rationale for utilizing mitochondria as targets for anti-cancer therapy.

Organism-adapted specificity of the allosteric regulation of pyruvate kinase in lactic acid bacteria

Pyruvate kinase (PYK) is a critical allosterically regulated enzyme that links glycolysis, the primary energy metabolism, to cellular metabolism. Lactic acid bacteria rely almost exclusively on glycolysis for their energy production under anaerobic conditions, which reinforces the key role of PYK in their metabolism. These organisms are closely related, but have adapted to a huge variety of native environments. They include food-fermenting organisms, important symbionts in the human gut, and antibiotic-resistant pathogens. In contrast to the rather conserved inhibition of PYK by inorganic phosphate, the activation of PYK shows high variability in the type of activating compound between different lactic acid bacteria. System-wide comparative studies of the metabolism of lactic acid bacteria are required to understand the reasons for the diversity of these closely related microorganisms. These require knowledge of the identities of the enzyme modifiers. Here, we predict potential allosteric activators of PYKs from three lactic acid bacteria which are adapted to different native environments. We used protein structure-based molecular modeling and enzyme kinetic modeling to predict and validate potential activators of PYK. Specifically, we compared the electrostatic potential and the binding of phosphate moieties at the allosteric binding sites, and predicted potential allosteric activators by docking. We then made a kinetic model of Lactococcus lactis PYK to relate the activator predictions to the intracellular sugar-phosphate conditions in lactic acid bacteria. This strategy enabled us to predict fructose 1,6-bisphosphate as the sole activator of the Enterococcus faecalis PYK, and to predict that the PYKs from Streptococcus pyogenes and Lactobacillus plantarum show weaker specificity for their allosteric activators, while still having fructose 1,6-bisphosphate play the main activator role in vivo. These differences in the specificity of allosteric activation may reflect adaptation to different environments with different concentrations of activating compounds. The combined computational approach employed can readily be applied to other enzymes.

Recombinant production of Yersinia enterocolitica pyruvate kinase isoenzymes PykA and PykF

The glycolytic enzyme pyruvate kinase (PK) generates ATP from ADP through substrate-level phosphorylation powered by the conversion of phosphoenolpyruvate to pyruvate. In contrast to other bacteria, Enterobacteriaceae, such as pathogenic yersiniae, harbour two pyruvate kinases encoded by pykA and pykF. The individual roles of these isoenzymes are poorly understood. In an attempt to make the Yersinia enterocolitica pyruvate kinases PykA and PykF amenable to structural and functional characterisation, we produced them untagged in Escherichia coli and purified them to near homogeneity through a combination of ion exchange and size exclusion chromatography, yielding more than 180 mg per litre of batch culture. The solution structure of PykA and PykF was analysed through small angle X-ray scattering which revealed the formation of PykA and PykF tetramers and confirmed the binding of the allosteric effector fructose-1,6-bisphosphate (FBP) to PykF but not to PykA.

Analysis and Design of Stimulus Response Curves of E. coli

Metabolism and signalling are tightly coupled in bacteria. Combining several theoretical approaches, a core model is presented that describes transcriptional and allosteric control of glycolysis in Escherichia coli. Experimental data based on microarrays, signaling components and extracellular metabolites are used to estimate kinetic parameters. A newly designed strain was used that adjusts the incoming glucose flux into the system and allows a kinetic analysis. Based on the results, prediction for intracelluar metabolite concentrations over a broad range of the growth rate could be performed and compared with data from literature.

KINETIC MODEL OF PHOSPHOFRUCTOKINASE-1 FROMESCHERICHIA COLI

This paper presents a kinetic model of phosphofructokinase-1 from Escherichia coli. A complete catalytic cycle has been reconstructed based on available information on the oligomeric structure of the enzyme and kinetic mechanism of its monomer. Applying the generalization of the Monod–Wyman–Changeux approach proposed by Popova and Sel'kov35–37to the reconstructed catalytic cycle rate equation has been derived. Dependence of the reaction rate on pH , magnesium, and effectors has been taken into account. Kinetic parameters have been estimated via fitting the rate equation against experimentally measured dependencies of initial rate on substrates, products, effectors, and pH available from the literature. The model of phosphofructokinase-1 predicts (1) cooperativity of binding both fructose-6-phosphate and ATPMg2-, (2) significant inhibition of the enzyme resulting from an increase in total concentration of ATP under the condition of fixed concentration of Mg2+ions, and (3) dual effect of ADP consisting of allosteric activation and product inhibition of the enzyme. Moreover, the model developed can be used in the kinetic modeling of biochemical pathways containing phosphofructokinase-1.

Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0

Over the past decade, a growing community of researchers has emerged around the use of constraint-based reconstruction and analysis (COBRA) methods to simulate, analyze and predict a variety of metabolic phenotypes using genome-scale models. The COBRA Toolbox, a MATLAB package for implementing COBRA methods, was presented earlier. Here we present a substantial update of this in silico toolbox. Version 2.0 of the COBRA Toolbox expands the scope of computations by including in silico analysis methods developed since its original release. New functions include (i) network gap filling, (ii) (13)C analysis, (iii) metabolic engineering, (iv) omics-guided analysis and (v) visualization. As with the first version, the COBRA Toolbox reads and writes systems biology markup language-formatted models. In version 2.0, we improved performance, usability and the level of documentation. A suite of test scripts can now be used to learn the core functionality of the toolbox and validate results. This toolbox lowers the barrier of entry to use powerful COBRA methods.

Functional analysis, overexpression, and kinetic characterization of pyruvate kinase from methicillin-resistant Staphylococcus aureus

Novel antimicrobial targets are urgently needed to overcome rising antibiotic resistance of important human pathogens including methicillin-resistant Staphylococcus aureus (MRSA). Here we report the essentiality and kinetic properties of MRSA pyruvate kinase (PK). Targetron-mediated gene disruption demonstrated PK is essential for S. aureus growth and survival, suggesting that this protein may be a potential drug target. The presence of the pfk (6-phosphofructokinase)-pyk operon in MRSA252, and the nonessential nature of PFK shown by targetron, further emphasized the essential role of PK in cell viability. The importance of PK in bacterial growth was confirmed by showing that its enzymatic activity peaked during the logarithmic phase of S. aureus growth. PK from Staphylococcus and several other species of bacteria have an extra C-terminal domain (CT) containing a phosphoenolpyruvate (PEP) binding motif. To elucidate the possible structure and function of this sequence, the quaternary structures and kinetic properties of the full-length MRSA PK and truncated MRSA PK lacking the CT domain were characterized. Our results showed that (1) MRSA PK is an allosteric enzyme with homotetramer architecture activated by AMP or ribose 5-phosphate (R5P), but not by fructose 1,6-bisphosphate (FBP), which suggests a different mode of allosteric regulation when compared with human isozymes, (2) the CT domain is not required for the tetramerization of the enzyme; homotetramerization occurred in a truncated PK lacking the domain, (3) truncated enzyme exhibited high affinity toward both PEP and ADP and exhibited hyperbolic kinetics toward PEP in the presence of activators (AMP and R5P) consistent with kinetic properties of full-length enzyme, indicating that the CT domain is not required for substrate binding or allosteric regulation observed in the holoenzyme, (4) the kinetic efficiency (k(cat)/S(0.5)) of truncated enzyme was decreased by 24- and 16-fold, in ligand-free state, toward PEP and ADP, respectively, but was restored by 3-fold in AMP-bound state, suggesting that the sequence containing the CT domain (Gly(473)-Leu(585)) plays a substantial role in enzyme activity and comformational stability, and (5) full-length MRSA PK activity was stimulated at low concentrations of ATP (e.g., 1 mM) and inhibited by inorganic phosphate and high concentrations of FBP (10 mM) and ATP (e.g., >2.5 mM), whereas for truncated enzyme, stimulation at low concentrations of ATP was lost. These findings suggest that the CT domain is involved in maintaining the specificity of allosteric regulation of MRSA PK by AMP, R5P, and ATP. The CT extension also encodes a protein domain with homology to enzyme I of the Escherichia coli sugar-PTS system, suggesting that MRSA PK may also exert an important regulatory role in sugar transport metabolism. These findings yield new insights into MRSA PK function and mode of allosteric regulation which may aid in the development of clinically important drugs targeting this enzyme and further define the role of the extra C-terminal domain in modulating the enzyme's activity.

In silico strategy to rationally engineer metabolite production: A case study for threonine in Escherichia coli

Genetic engineering of metabolic pathways is a standard strategy to increase the production of metabolites of economic interest. However, such flux increases could very likely lead to undesirable changes in metabolite concentrations, producing deleterious perturbations on other cellular processes. These negative effects could be avoided by implementing a balanced increase of enzyme concentrations according to the Universal Method [Kacser and Acerenza (1993) Eur J Biochem 216:361-367]. Exact application of the method usually requires modification of many reactions, which is difficult to achieve in practice. Here, improvement of threonine production via pyruvate kinase deletion in Escherichia coli is used as a case study to demonstrate a partial application of the Universal Method, which includes performing sensitivity analysis. Our analysis predicts that manipulating a few reactions is sufficient to obtain an important increase in threonine production without major perturbations of metabolite concentrations.

Specific ethanol production rate in ethanologenic Escherichia coli strain KO11 Is limited by pyruvate decarboxylase

Modification of ethanol productivity and yield, using mineral medium supplemented with glucose or xylose as carbon sources, was studied in ethanologenic Escherichia coli KO11 by increasing the activity of five key carbon metabolism enzymes. KO11 efficiently converted glucose or xylose to ethanol with a yield close to 100% of the theoretical maximum when growing in rich medium. However, when KO11 ferments glucose or xylose in mineral medium, the ethanol yields decreased to only 70 and 60%, respectively. An increase in GALP(Ec) (permease of galactose-glucose-xylose) or PGK(Ec) (phosphoglycerate kinase) activities did not change xylose or glucose and ethanol flux. However, when PDC(Zm) (pyruvate decarboxylase from Zymomonas mobilis) activity was increased 7-fold, the yields of ethanol from glucose or xylose were increased to 85 and 75%, respectively, and organic acid formation rates were reduced. Furthermore, as a response to a reduction in acetate and ATP yield, and a limited PDC(Zm) activity, an increase in PFK(Ec) (phosphofructokinase) or PYK(Bs) (pyruvate kinase from Bacillus stearothermophilus) activity drastically reduced glucose or xylose consumption and ethanol formation flux. This experimental metabolic control analysis showed that ethanol flux in KO11 is negatively controlled by phosphofructokinase and pyruvate kinase, and positively influenced by the PDC(Zm) activity level.

A novel GDP-dependent pyruvate kinase isozyme from Toxoplasma gondii localizes to both the apicoplast and the mitochondrion

We previously reported a cytosolic pyruvate kinase (EC 2.7.1.40) from Toxoplasma gondii (TgPyKI) that differs from most eukaryotic pyruvate kinases in being regulated by glucose 6-phosphate rather than fructose 1,6-diphosphate. Another putative pyruvate kinase (TgPyKII) was identified from parasite genome, which exhibits 32% amino acid sequence identity to TgPyKI and retains pyruvate kinase signature motifs and amino acids essential for substrate binding and catalysis. Whereas TgPyKI is most closely related to plant/algal enzymes, phylogenetic analysis suggests a proteobacterial origin for TgPyKII. Enzymatic characterization of recombinant TgPyKII shows a high pH optimum at 8.5, and a preference for GDP as a phosphate recipient. Catalytic activity is independent of K+, and no allosteric or regulatory effects were observed in the presence of fructose 1,6-diphosphate, fructose 2,6-diphosphate, glucose 6-phosphate, ribose 5-phosphate, AMP, or ATP. Unlike TgPyKI, native TgPyKII activity was exclusively associated with the membranous fraction of a T. gondii tachyzoite lysate. TgPyKII possesses a long N-terminal extension containing five putative start codons before the conserved region and localizes to both apicoplast and mitochondrion by immunofluorescence assay using native antibody and fluorescent protein fusion to the N-terminal extension. Further deletional and site-directed mutagenesis suggests that a translation product from 1st Met is responsible for the localization to the apicoplast, whereas one from 3rd Met is for the mitochondrion. This is the first study of a potential mitochondrial pyruvate kinase in any system.

A feed-forward loop guarantees robust behavior in Escherichia coli carbohydrate uptake

MOTIVATION

In Escherichia coli, the phosphoenolpyruvate: carbohydrate phosphotransferase system acts like a sensory element which is able to measure the flux through glycolysis. Since the output of the sensor, the phosphorylated form of protein EIIA, is connected to the activity of the global transcription factor Crp, the kinetic and structural properties of the system are important for the understanding of the overall cellular behavior.

RESULTS

A family of mathematical models is presented, varying with respect to their degree of complexity (number of reactions that are taken into account, number of parameters) that show a structurally and quantitatively robust behavior. The models describe a set of experimental data that relates the output of the sensor to the specific growth rate. A central element that is responsible for the structural robustness is a feed-forward loop in the glycolysis, namely the activation of the pyruvate kinase reaction by a metabolite of the upper part of the glycolysis. The robustness is shown for variations of the measured data as well as for variations of the parameters.

AVAILABILITY

MATLAB files for model simulations are available on http://www.mpi-magdeburg.mpg.de/people/kre/robust/ A short description of the files provided on this site can be found in the Supporting information.

Glycolysis and Flux Control

Central metabolism of carbohydrates uses the Embden-Meyerhof-Parnas (EMP), pentose phosphate (PP), and Entner-Doudoroff (ED) pathways. This review reviews the biological roles of the enzymes and genes of these three pathways of E. coli. Glucose, pentoses, and gluconate are primarily discussed as the initial substrates of the three pathways, respectively. The genetic and allosteric regulatory mechanisms of glycolysis and the factors that affect metabolic flux through the pathways are considered here. Despite the fact that a lot of information on each of the reaction steps has been accumulated over the years for E. coli, surprisingly little quantitative information has been integrated to analyze glycolysis as a system. Therefore, the review presents a detailed description of each of the catalytic steps by a systemic approach. It considers both structural and kinetic aspects. Models that include kinetic information of the reaction steps will always contain the reaction stoichiometry and therefore follow the structural constraints, but in addition to these also kinetic rate laws must be fulfilled. The kinetic information obtained on isolated enzymes can be integrated using computer models to simulate behavior of the reaction network formed by these enzymes. Successful examples of such approaches are the modeling of glycolysis in S. cerevisiae, the parasite Trypanosoma brucei, and the red blood cell. With the rapid developments in the field of Systems Biology many new methods have been and will be developed, for experimental and theoretical approaches, and the authors expect that these will be applied to E. coli glycolysis in the near future.

Identification of factors involved in recovery of heat-injured Salmonella Enteritidis

Proteins and genes involved in the recovery of heat-injured Salmonella Enteritidis were investigated. Salmonella Enteritidis cells cultured overnight in tryptic soy broth (TSB; nonselective medium) were suspended in citric acid-disodium hydrogen phosphate buffer (pH 6). After heat treatment at 55 degrees C for 15 min, the culturable counts measured by tryptic soy agar (TSA; nonselective medium) decreased from 10(8) to 10(7) CFU/ml. On the other hand, culturable counts measured by desoxycholate-hydrogen sulfite-lactose (DHL) agar (selective medium) were decreased from 10(8) to 10(4) CFU/ml by the same treatment. The results suggest that 99.9% of Salmonella Enteritidis detected on TSA were injured but recoverable. When injured Salmonella Enteritidis was incubated in TSB, the culturable count measured by TSA did not increase for 2 h, whereas that by DHL agar increased after incubation for 30 min. After incubation for 2 h, the culturable count measured by DHL agar reached a similar level with that by TSA, indicating that Salmonella Enteritidis had recovered. The two-dimensional polyacrylamide gel electrophoresis analysis revealed that elongation factor G (FusA) and pyruvate kinase (PykF) specifically increased in the cells just after heat treatment and in the recovery cells. The levels of transcription of 86 stress-inducible genes were also investigated by reverse transcription PCR. Nineteen heat-inducible (clpB, clpX, degP, dnaJ, fkpA, ftsJ, gapA, hflB, hslJ, hslU, hslV, htpG, htrA, lon, mopA, mopB, mreB, rpoE, and ppiD), and 12 oxidative-stress and DNA damage-inducible (ahpC, ahpF, fldB, fur, grxA, dinF, katG, mutM, recA, soxR, trxC, and zwf) genes were transcribed extensively during recovery in TSB. The results obtained in this study will be used to develop the media or culture conditions that will promote recovery for the detection of food poisoning bacteria, including injured cells from food products.

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.

Pyruvate kinase from Chlamydia trachomatis is activated by fructose-2,6-bisphosphate

Pyruvate kinase is the final regulatory point in the catabolic Embden-Meyerhoff-Parnas pathway, which controls the carbon flux of glycolytic intermediates and regulates the level of ATP in the cell. In a previous study, we identified, cloned and sequenced pyruvate kinase from the obligate intracellular bacterium Chlamydia trachomatis and demonstrated that the enzyme was active in crude extract. Here, we report the kinetic properties of highly purified C. trachomatis pyruvate kinase. The results indicate that C. trachomatis pyruvate kinase is 53.5 kDa with a pH optima of 7.3. Kinetic studies show that C. trachomatis pyruvate kinase requires both K+ and Mg2+ ions for activity, exhibits sigmoidal kinetics with respect to phosphoenolpyruvate and Michaelis-Menten kinetics with respect to ADP. In addition, C. trachomatis pyruvate kinase is able to use alternative nucleoside diphosphates as phosphate acceptors, although it shows the greatest activity with ADP. In contrast to other bacterial pyruvate kinases that are activated by AMP, our data show that AMP, in addition to ATP and GTP, inhibits C. trachomatis pyruvate kinase. Surprisingly, unlike any other known bacterial pyruvate kinase, C. trachomatis pyruvate kinase was allosterically activated by fructose-2,6-bisphosphate, an important regulatory metabolite that has only been reported in eukaryotes.

Micellar electrokinetic chromatography as a complementary method to sodium dodecyl sulfate-polyacrylamide gel electrophoresis for studying limited proteolysis of proteins

Micellar electrokinetic chromatography (MEKC) has been utilized as an analytical method to perform investigations on limited proteolysis of proteins. To this purpose partial proteolysis experiments with a series of proteinases were performed, utilizing as model protein pyruvate kinase (PK) from Escherichia coli, an enzyme that is regulated allosterically by fructose 1,6-bisphosphate (FBP). Data obtained with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and MEKC were compared; the profiles generated by submitting digests of PK treated with different proteinases in the presence and absence of FBP to electrophoretic analysis provided a useful adjunct for a better understanding of the effects of the allosteric activator on the conformation of the model enzyme. MEKC was also found to be a convenient technique for determining the kinetics of substrate proteolysis.

Incorporating qualitative knowledge in enzyme kinetic models using fuzzy logic

Modeling of metabolic pathway dynamics requires detailed kinetic equations at the enzyme level. In particular, the kinetic equations must account for metabolite effectors that contribute significantly to the pathway regulation in vivo. Unfortunately, most kinetic rate laws available in the literature do not consider all the effectors simultaneously, and much kinetic information exists in a qualitative or semiquantitative form. In this article, we present a strategy to incorporate such information into the kinetic equation. This strategy uses fuzzy logic-based factors to modify algebraic rate laws that account for partial kinetic characteristics. The parameters introduced by the fuzzy factors are then optimized by use of a hybrid of simplex and genetic algorithms. The resulting model provides a flexible form that can simulate various kinetic behaviors. Such kinetic models are suitable for pathway modeling without complete enzyme mechanisms. Three enzymes in Escherichia coli central metabolism are used as examples: phosphoenolpyruvate carboxylase; phosphoenolpyruvate carboxykinase; and pyruvate kinase I. Results show that, with fuzzy logic-augmented models, the kinetic data can be much better described. In particular, complex behavior, such as allosteric inhibition, can be captured using fuzzy rules. The resulting models, even though they do not provide additional physical meaning in enzyme mechanisms, allow the model to incorporate semiquantitative information in metabolic pathway models.

Pyruvate kinase of Trypanosoma brucei: overexpression, purification, and functional characterization of wild-type and mutated enzyme

A procedure was developed for overexpression of Trypanosoma brucei pyruvate kinase in Escherichia coli. The enzyme was purified to near-homogeneity from the bacterial lysate by first removing nucleic acids and contaminating proteins by protamine sulfate precipitation and subsequent passage over a phosphocellulose column. The purified protein is essentially indistinguishable in its physicochemical and kinetic properties from the enzyme purified from trypanosomes. Furthermore, experiments were undertaken to locate the binding site of the allosteric effector fructose 2,6-bisphosphate. Regulation of pyruvate kinase by this effector is unique to trypanosomes and related protozoan organisms. Therefore, a three-dimensional structure model of the enzyme was made, and a putative effector-binding site could be identified in an interdomain cleft. Four residues in this cleft were mutated, and the mutant proteins were produced and purified, using the same methodology as for the wild-type pyruvate kinase. Some mutants showed only minor changes in the activation by the effector. However, substitution of Arg22 by Gly resulted in a 9.2-fold higher S(0.5) for phosphoenolpyruvate and a significantly smaller kcat than the wild-type enzyme. Furthermore, the apparent affinity of this mutant for the allosteric effectors fructose 1,6-bisphosphate and fructose 2,6-bisphosphate was 8.2- and 5.2-fold lower than that of its wild-type counterpart. Effector binding was also affected, although to a lesser extent, in a mutant Phe463Val. These data indicate that particularly residue Arg22, but also Phe463, are somehow involved in the binding of the allosteric effectors.

The allosteric regulation of pyruvate kinase by fructose-1,6-bisphosphate

BACKGROUND

Yeast pyruvate kinase (PK) catalyzes the final step in glycolysis. The enzyme therefore represents an important control point and is allosterically activated by fructose-1,6-bisphosphate (FBP). In mammals the enzyme is found as four different isozymes with different regulatory properties: two of these isozymes are produced by alternate splicing. The allosteric regulation of PK is directly related to proliferation of certain cell types, as demonstrated by the expression of an allosterically regulated isozyme in tumor cells. A model for the allosteric transition from the inactive (T) state to the active (R) state has been proposed previously, but until now the FBP-binding site had not been identified.

RESULTS

We report here the structures of PK from yeast complexed with a substrate analog and catalytic metal ions in the presence and absence of bound FBP. The allosteric site is located 40 A from the active site and is entirely located in the enzyme regulatory (C) domain. A phosphate-binding site for the allosteric activator is created by residues encoded by a region of the gene corresponding to the alternately spliced exon of mammalian isozymes. FBP activation appears to induce several conformational changes among active-site sidechains through a mechanism that is most likely to involve significant domain motions, as previously hypothesized.

CONCLUSIONS

The structure and location of the allosteric activator site agrees with the pattern of alternate genetic splicing of the PK gene in multicellular eukaryotes that distinguishes between a non-regulated isozyme and the regulated fetal isozymes. The conformational differences observed between the active sites of inactive and fully active PK enzymes is in agreement with the recently determined thermodynamic mechanism of allosteric activation through a 'metal relay' that increases the affinity of the enzyme for its natural phosphoenolpyruvate substrate.

Metal-ion-mediated allosteric triggering of yeast pyruvate kinase. 1. A multidimensional kinetic linked-function analysis

Regulation of the glycolytic pathway is considered to be primarily achieved by the carbon metabolites resulting from glucose metabolism [e.g., fructose 1,6-diphosphate (FDP), phosphoenolpyruvate (PEP), and citrate] and by the ATP charge of the cell. The divalent cations (e.g., Mg2+ and Mn2+) have not been considered as having regulatory roles in glycolysis, although they are involved in almost every enzyme-catalyzed reaction in the pathway. Using a kinetic linked-function analysis of steady-state kinetic data for the interactions of PEP, FDP, and Mn2+ with yeast pyruvate kinase (YPK), we have found that the divalent metal is the principal trigger of the allosteric responses observed with this enzyme. The interaction of Mn2+ to YPK enhances the interaction of FDP by -1.6 kcal/mol and the interaction of PEP by -2.8 kcal/mol. The simultaneous interaction of all three of these ligands to YPK is favored by -4.3 kcal/mol over the sum of their independent binding free energies. Surprisingly, the binding of the allosteric activator FDP does not directly influence the binding of the substrate PEP since a coupling free energy near zero was calculated for these two ligands. Thus, communication between the PEP and FDP sites occurs structurally through the metal by an allosteric relay mechanism. These conclusions are supported by results of a thermodynamic linked-function analysis of direct binding data for the interactions of PEP, FDP, and Mn2+ with YPK [Mesecar, A. D., & Nowak, T. (1997) Biochemistry (following paper in this series)]. Our findings raise important questions as to the possible roles of divalent metals in modulating multiligand interactions with YPK and in the regulation of the glycolytic pathway.

Affinity labelling of the catalytic and allosteric ATP binding sites on pyruvate kinase type I from Escherichia coli

The allosterically regulated pyruvate kinase type I (PKI) from E. coli was inactivated by the ATP analog 2',3'-dialdehyde ATP (o-ATP) with a Ki of 3.6 mM. ATP and phosphoenolpyruvate protected the enzyme activity while the allosteric activator fructose 1,6-bisphosphate enhanced the rate of inactivation. Incubation with o-ATP, followed by reduction of the formed Schiff bases with radioactive sodium borohydride, was employed to determine the ATP binding sites of PKI. After tryptic digestion, the purification of the labelled peptides and the sequence analysis allowed to identify four modified lysyl residues, namely Lys173, Lys175, Lys272, and Lys317 of the known DNA-deduced sequence of PKI. The close lysines 173 and 175 reacted with o-ATP in a mutually exclusive way and accounted together for 53% of the recovered radioactivity, the rest being distributed on Lys272 (31%) and Lys317 (16%). When fitted on the available three-dimensional structure of muscle pyruvate kinase, the position of the modified lysines defines both the catalytic and the allosteric ATP binding sites on PKI.

Dynamic behavior in glycolytic oscillations with phase shifts

Practically all of the studies of glycolytic oscillations in homogeneous spatial mediums have been performed through the construction of systems of ordinary differential equations and the search for their solutions. In this kind of modelling, the system dynamic behavior is considered to depend only on the values adopted by the parameters related to the dependent variables. In the present work, the modeling of a biochemical system through a system of functional differential equations with delay allows us to analyse the consequences that the variations in the parametric values linked to the independent variable (time) have upon the integral solutions of the system. In our model, the delays correspond with phase shifts in the initial functions for two dependent variables. The results of our researches show that when a instability-generating multienzymatic mechanism suffers variations of the delay time in any of its variables, a wide range of different dynamic responses can be produced. Our work is presented as an enlargement on the dynamic study of biochemical oscillations in general and, particularly, the glycolytic oscillations, under the consideration of the existence of variations in the phase shifts during the oscillations of metabolites involved in the studied reactive processes.

Alteration of the biochemical valves in the central metabolism of Escherichia coli

Although E. coli central metabolism has been studied for several decades, many regulatory features are still unknown. To achieve the goal of rational manipulation of cellular metabolism, it is important to understand how E. coli responds to overexpressed enzymes. By studying the biochemical control of fluxes between PEP, pyruvate, and OAA, we have addressed some fundamental questions that may prove to be essential for applications in metabolic engineering. First, we found that simultaneous overexpression of Pck and Ppc, or Pps alone in the presence of glucose leads to phenotypes consistent with futile cycline. In contrast to our expectation, futile cycling per se does not affect the growth rate significantly. However, excessive futile cycling may cause competitive disadvantage in the natural environment. Overexpression of Pck caused growth inhibition but no futile cycling. Therefore, E. coli controls the expression of gluconeogenic enzymes not only to avoid excessive futile cycling, but also to prevent toxicity effects. In metabolic engineering, futile cycling may be used as a strategy to stimulate metabolism for either production of metabolites or digestion of toxic wastes. Second, we found that the expression levels of Pps and Pck in E. coli are not optimal for growth on pyruvate and succinate, respectively. Overexpression of these enzymes increases the growth rate on pyruvate and on succinate, respectively, indicating that the slow growth rates on these substrates are at least partially caused by the insufficient supply of PEP and its derivatives. Moreover, E. coli also has not optimized the Ppc level for optimal growth yield on glucose in uncontrolled batch cultures. These results demonstrate that the central metabolism is not optimized for growth under defined laboratory conditions. Thus, the possibility exists that adjustment of native enzyme levels in the central metabolism can improve bioreactor performance. Third, we found that overexpression of Pck affects the transcriptional levels of unrelated genes. This example indicates that physiological responses to enzyme (over)expression should be interpreted cautiously, as changing the expression level of a specific enzyme may affect many unlinked genes. Similar results have also been obtained by use of two-dimensional electrophoresis of proteins from E. coli. Although more questions remain to be answered, fast progress in the area of metabolic engineering can be expected in the near future.

Pyruvate kinase from Trichomonas vaginalis, an allosteric enzyme stimulated by ribose 5-phosphate and glycerate 3-phosphate

Trichomonas vaginalis pyruvate kinase was purified over 1750 fold to a specific activity greater than 100 mumol min-1 (mg protein)-1. The enzyme is a tetramer of M(r) 266,000, consisting of subunits of M(r) 53,000 and 56,000 in equivalent amounts. Its activity was dependent on the presence of magnesium but was not stimulated by potassium or ammonium. The enzyme exhibited positive cooperativity towards phosphoenolpyruvate and was inhibited by inorganic phosphate, which increased the sigmoidicity of the saturation curve for phosphoenolpyruvate without affecting maximal activity. It was heterotropically stimulated by ribose 5-phosphate and glycerate 3-phosphate, not previously known to act on eukaryotic pyruvate kinases, but was unaffected by known effectors of most pyruvate kinases, including fructose 1,6-bisphosphate and fructose 2,6-bisphosphate.

Some kinetic properties of pyruvate kinase from Trypanosoma brucei

We have studied the kinetics of the allosteric interactions of pyruvate kinase from Trypanosoma brucei. The kinetics for phosphoenolpyruvate depended strongly on the nature of the bivalent metal ions. Pyruvate kinase activated by Mg2+ had the highest catalytic activity, but also the highest S0.5 for phosphoenolpyruvate, while the opposite was true for pyruvate kinase activated by Mn2+. The reaction rates of Mg(2+)-pyruvate kinase and Mn(2+)-pyruvate kinase were clearly allosteric with respect to phosphoenolpyruvate, while the kinetics with Co(2+)-pyruvate kinase were hyperbolic. However, Co(2+)-pyruvate kinase was still sensitive to heterotropic activation. Trypanosomal pyruvate kinase is unique in that the best activator was fructose 2,6-bisphosphate. Ribulose 1,5-bisphosphate and 5-phosphorylribose 1-pyrophosphate were also strong heterotropic activators, which were much more effective than fructose 1,6-bisphosphate and glucose 1,6-bisphosphate. In the presence of the heterotropic activators, the sigmoidal kinetics with respect to phosphoenolpyruvate and the bivalent metal ions were modified as were the concentrations of phosphoenolpyruvate and the bivalent metal ions needed to attain the maximal activity. Maximal activities were not significantly changed with Mg2+ and Mn2+ as the activating metal ions. Moreover, with Co2+ and fructose 2,6-bisphosphate or ribulose 1,5-bisphosphate or 5-phosphorylribose 1-pyrophosphate, the maximal activity was significantly reduced. Ribulose 1,5-bisphosphate and 5-phosphorylribose 1-pyrophosphate resembled fructose 2,6-bisphosphate rather than fructose 1,6-bisphosphate and glucose 1,6-bisphosphate in their action in that the K0.5 values for the former 3 compounds increased when Mg2+ was replaced by Co2+, while the K0.5 for fructose 1,6-bisphosphate and glucose 1,6-bisphosphate increased.(ABSTRACT TRUNCATED AT 250 WORDS)

Characterization of pyruvate kinase of Trypanosoma brucei and its role in the regulation of carbohydrate metabolism

Pyruvate kinase from Trypanosoma brucei is a labile enzyme, losing its activity within several hours. In mixtures containing 50 mM triethanolamine buffer, pH 7.2, 25% glycerol and 0.5 mM inorganic phosphate the enzyme remained active and could be purified to homogeneity with a specific activity of 417 units mg-1 and a yield of 65%. The enzyme has an activation energy of 31.9 kJ mol-1. Magnesium and potassium ions are essential for activity. Cobalt or manganese ions replace Mg2+ but this leads to a decrease in maximal velocity. Potassium ions can be substituted by ammonium ions, while sodium ions behave as a competitive inhibitor with respect to both K+ and NH4+. All metal ions studied displayed sigmoidal kinetics. The enzyme is activated, with decreasing efficiency by fructose 2-phosphorothioate 6-phosphate, fructose 2,6-bisphosphate, fructose 1,6-bisphosphate and glucose 1,6-bisphosphate. They all display hyperbolic kinetics. Glycerate 2,3-bisphosphate, glyceraldehyde 3-phosphate, CoASAc, oxalate, AMP, ADP, and ATP inhibit the enzyme. At substrate saturation PK was activated by Pi up to a concentration of 0.8 mM. At higher Pi concentrations the enzyme is inhibited. The enzyme is unaffected by most amino acids, only phenylalanine stimulates and tyrosine inhibits.

Fructose-1,6-bisphosphate-activated pyruvate kinase from Escherichia coli. Nature of bonds involved in the allosteric mechanism

The allosteric properties of the fructose-1,6-bis-phosphate-activated pyruvate kinase from Escherichia coli were examined in the presence of a number of fructose bisphosphate analogues, as well as of increased ionic strength (NaCl) and of the hydrogen-bond-breaking agent, formamide. Fructose 2,6-bisphosphate, ribulose 1,5-bisphosphate and 5-phosphorylribose 1-pyrophosphate gave allosteric activation (additive to that of fructose 1,6-bisphosphate). Formamide always decreased Vmax, but left unchanged the Km for phosphoenolpyruvate, while it decreased the concentration of fructose bisphosphate required to give half-maximal activity (K0.5). NaCl increased the K0.5 for both phosphoenolpyruvate and fructose bisphosphate, leaving Vmax unchanged. These results are consistent with ionic binding of fructose bisphosphate through phosphates and with a critical role of hydrogen bonds in stabilizing both the inactive and the active enzyme conformers.