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

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.

Control and regulation of the pyrophosphate-dependent glucose metabolism in Entamoeba histolytica

Entamoeba histolytica has neither Krebs cycle nor oxidative phosphorylation activities; therefore, glycolysis is the main pathway for ATP supply and provision of carbon skeleton precursors for the synthesis of macromolecules. Glucose is metabolized through fermentative glycolysis, producing ethanol as its main end-product as well as some acetate. Amoebal glycolysis markedly differs from the typical Embden-Meyerhof-Parnas pathway present in human cells: (i) by the use of inorganic pyrophosphate, instead of ATP, as the high-energy phospho group donor; (ii) with one exception, the pathway enzymes can catalyze reversible reactions under physiological conditions; (iii) there is no allosteric regulation and sigmoidal kinetic behavior of key enzymes; and (iv) the presence of some glycolytic and fermentation enzymes similar to those of anaerobic bacteria. These peculiarities bring about alternative mechanisms of control and regulation of the PPi-dependent fermentative glycolysis in the parasite in comparison to the ATP-dependent and allosterically regulated glycolysis in many other eukaryotic cells. In this review, the current knowledge of the carbohydrate metabolism enzymes in E. histolytica is analyzed. Thermodynamics and stoichiometric analyses indicate 2 to 3.5 ATP yield per glucose metabolized, instead of the often presumed 5 ATP/glucose ratio. PPi derived from anabolism seems insufficient for PPi-glycolysis; hence, alternative ways of PPi supply are also discussed. Furthermore, the underlying mechanisms of control and regulation of the E. histolytica carbohydrate metabolism, analyzed by applying integral and systemic approaches such as Metabolic Control Analysis and kinetic modeling, contribute to unveiling alternative and promising drug targets.

Sequence analysis and characterization of pyruvate kinase from Clonorchis sinensis, a 53.1-kDa homopentamer, implicated immune protective efficacy against clonorchiasis

BACKGROUND

Clonorchis sinensis, the causative agent of clonorchiasis, is classified as one of the most neglected tropical diseases and affects more than 15 million people globally. This hepatobiliary disease is highly associated with cholangiocarcinoma. As key molecules in the infectivity and subsistence of trematodes, glycolytic enzymes have been targets for drug and vaccine development. Clonorchis sinensis pyruvate kinase (CsPK), a crucial glycolytic enzyme, was characterized in this research.

RESULTS

Differences were observed in the sequences and spatial structures of CsPK and PKs from humans, rats, mice and rabbits. CsPK possessed a characteristic active site signature (IKLIAKIENHEGV) and some unique sites but lacked the N-terminal domain. The predicted subunit molecular mass (Mr) of CsPK was 53.1 kDa. Recombinant CsPK (rCsPK) was a homopentamer with a Mr. of approximately 290 kDa by both native PAGE and gel filtration chromatography. Significant differences in the protein and mRNA levels of CsPK were observed among four life stages of C. sinensis (egg, adult worm, excysted metacercaria and metacercaria), suggesting that these developmental stages may be associated with diverse energy demands. CsPK was widely distributed in adult worms. Moreover, an intense Th1-biased immune response was persistently elicited in rats immunized with rCsPK. Also, rat anti-rCsPK sera suppressed C. sinensis adult subsistence both in vivo and in vitro.

CONCLUSIONS

The sequences and spatial structures, molecular mass, and expression profile of CsPK have been characterized. rCsPK was indicated to be a homopentamer. Rat anti-rCsPK sera suppressed C. sinensis adult subsistence both in vivo and in vitro. CsPK is worthy of further study as a promising target for drug and vaccine development.

In Vitro Testing of Potential Entamoeba histolytica Pyruvate Phosphate Dikinase Inhibitors

Adverse effects and resistance to metronidazole have motivated the search for new antiamoebic agents against Entamoeba histolytica. Control of amoeba growth may be achieved by inhibiting the function of the glycolytic enzyme and pyruvate phosphate dikinase (PPDK). In this study, we screened 10 compounds using an in vitro PPDK enzyme assay. These compounds were selected from a virtual screening of compounds in the National Cancer Institute database. The antiamoebic activity of the selected compounds was also evaluated by determining minimal inhibitory concentrations (MICs) and IC₅₀ values using the nitro-blue tetrazolium reduction assay. Seven of the 10 compounds showed inhibitory activities against the adenosine triphosphate (ATP)/inorganic phosphate binding site of the ATP-grasp domain. Two compounds, NSC349156 (pancratistatin) and NSC228137 (7-ethoxy-4-[4-methylphenyl] sulfonyl-3-oxido-2, 1, 3-benzoxadiazol-3-ium), exhibited inhibitory effects on the growth of E. histolytica trophozoites with MIC values of 25 and 50 μM, and IC₅₀ values of 14 and 20.7 μM, respectively.

Proteomic analysis of Clostridium thermocellum ATCC 27405 reveals the upregulation of an alternative transhydrogenase-malate pathway and nitrogen assimilation in cells grown on cellulose

Clostridium thermocellum is a Gram-positive thermophilic anaerobic bacterium with the ability to directly convert cellulosic biomass into useful products such as ethanol and hydrogen. In this study, a quantitative comparative proteomic analysis of the organism was performed to identify proteins and biochemical pathways that are differentially utilized by the organism after growth on cellobiose or cellulose. The cytoplasmic and membrane proteomes of C. thermocellum grown on cellulose or cellobiose were quantitatively compared using a metabolic (15)N isotope labelling method in conjunction with nanoLC-ESI-MS/MS (liquid chromatography - electrospray ionization - tandem mass spectrometry). In total, 1255 proteins were identified in the study, and 129 of those were able to have their relative abundance per cell compared in at least one cellular compartment in response to the substrate provided. This study reveals that cells grown on cellulose increase their abundance of phosphoenolpyruvate carboxykinase while decreasing the abundance of pyruvate dikinase and oxaloacetate decarboxylase, suggesting that the organism diverts carbon flow into a transhydrogenase-malate pathway that can increase the production of the biosynthetic intermediates NADPH and GTP. Glutamate dehydrogenase was also found to have increased abundance in cellulose-grown cells, suggesting that the assimilation of ammonia is upregulated in cells grown on the cellulosic substrates. The results illustrate a mechanism by which C. thermocellum can divert carbon into alternative pathways for the purpose of producing biosynthetic intermediates necessary to respond to growth on cellulose, including transhydrogenation of NADH to NADPH and increased nitrogen assimilation.

Biochemistry and evolution of anaerobic energy metabolism in eukaryotes

Major insights into the phylogenetic distribution, biochemistry, and evolutionary significance of organelles involved in ATP synthesis (energy metabolism) in eukaryotes that thrive in anaerobic environments for all or part of their life cycles have accrued in recent years. All known eukaryotic groups possess an organelle of mitochondrial origin, mapping the origin of mitochondria to the eukaryotic common ancestor, and genome sequence data are rapidly accumulating for eukaryotes that possess anaerobic mitochondria, hydrogenosomes, or mitosomes. Here we review the available biochemical data on the enzymes and pathways that eukaryotes use in anaerobic energy metabolism and summarize the metabolic end products that they generate in their anaerobic habitats, focusing on the biochemical roles that their mitochondria play in anaerobic ATP synthesis. We present metabolic maps of compartmentalized energy metabolism for 16 well-studied species. There are currently no enzymes of core anaerobic energy metabolism that are specific to any of the six eukaryotic supergroup lineages; genes present in one supergroup are also found in at least one other supergroup. The gene distribution across lineages thus reflects the presence of anaerobic energy metabolism in the eukaryote common ancestor and differential loss during the specialization of some lineages to oxic niches, just as oxphos capabilities have been differentially lost in specialization to anoxic niches and the parasitic life-style. Some facultative anaerobes have retained both aerobic and anaerobic pathways. Diversified eukaryotic lineages have retained the same enzymes of anaerobic ATP synthesis, in line with geochemical data indicating low environmental oxygen levels while eukaryotes arose and diversified.

Experimental validation of metabolic pathway modeling

In the search for new drug targets in the human parasite Entamoeba histolytica, metabolic control analysis was applied to determine, experimentally, flux control distribution of amebal glycolysis. The first (hexokinase, hexose-6-phosphate isomerase, pyrophosphate-dependent phosphofructokinase (PP(i)-PFK), aldolase and triose-phosphate isomerase) and final (3-phosphoglycerate mutase, enolase and pyruvate phosphate dikinase) glycolytic segments were reconstituted in vitro with recombinant enzymes under near-physiological conditions of pH, temperature and enzyme proportion. Flux control was determined by titrating flux with each enzyme component. In parallel, both glycolytic segments were also modeled by using the rate equations and kinetic parameters previously determined. Because the flux control distribution predicted by modeling and that determined by reconstitution were not similar, kinetic interactions among all the reconstituted components were experimentally revised to unravel the causes of the discrepancy. For the final segment, it was found that 3-phosphoglycerate was a weakly competitive inhibitor of enolase, whereas PP(i) was a moderate inhibitor of 3-phosphoglycerate mutase and enolase. For the first segment, PP(i) was both a strong inhibitor of aldolase and a nonessential mixed-type activator of amebal hexokinase; in addition, lower V(max) values for hexose-6-phosphate isomerase, PP(i)-PFK and aldolase were induced by PP(i) or ATP inhibition. It should be noted that PP(i) and other metabolites were absent from the 3-phosphoglycerate mutase and enolase or aldolase and hexokinase kinetics experiments, but present in reconstitution experiments. Only by incorporating these modifications in the rate equations, modeling predicted values of flux control distribution, flux rate and metabolite concentrations similar to those experimentally determined. The experimentally validated segment models allowed 'in silico experimentation' to be carried out, which is not easy to achieve in in vivo or in vitro systems. The results predicted a nonsignificant effect on flux rate and flux control distribution by adding parallel routes (pyruvate kinase for the final segment and ATP-dependent PFK for the first segment), because of the much lower activity of these enzymes in the ameba. Furthermore, modeling predicted full flux-control by 3-phosphoglycerate mutase and hexokinase, in the presence of low physiological substrate and product concentrations. It is concluded that the combination of in vitro pathway reconstitution with modeling and enzyme kinetics experimentation permits a more comprehensive understanding of the pathway behavior and control properties.

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.

Structure and content of the Entamoeba histolytica genome

The intestinal parasite Entamoeba histolytica is one of the first protists for which a draft genome sequence has been published. Although the genome is still incomplete, it is unlikely that many genes are missing from the list of those already identified. In this chapter we summarise the features of the genome as they are currently understood and provide previously unpublished analyses of many of the genes.

The catalyzing role of PPDK in Giardia lamblia

Giardia lamblia is an early branching eukaryotic microorganism that derives its metabolic energy primarily from anaerobic glycolysis. In most organisms, glycolysis is catalyzed by pyruvate kinase (PK), allowing the generation of two ATP molecules from one molecule of pyruvate. Giardia has both PK and pyrophosphate-dependent pyruvate phosphate dikinase (PPDK), which catalyzes the generation of five ATP molecules from pyruvate by pyrophosphate-dependent glycolysis and offers a potential selective advantage. In order to evaluate the importance of pyrophosphate-dependent glycolysis, we used ribozyme-mediated cleavage of the PPDK transcript to decrease PPDK transcript levels to 20% of normal. The accompanying decrease in PPDK enzyme activity decreased ATP levels to 3% of normal and increased glycogen deposition, confirming the importance pyrophosphate-mediated glycolysis that was previously suggested by cell lysate studies. PPDK is not found in vertebrates, so specific inhibitors may be useful for treatment of infections caused by anaerobic protists that depend on pyrophosphate-dependent glycolysis.

Kinetic modeling can describe in vivo glycolysis in Entamoeba histolytica

Glycolysis in the human parasite Entamoeba histolytica is characterized by the absence of cooperative modulation and the prevalence of pyrophosphate-dependent (over ATP-dependent) enzymes. To determine the flux-control distribution of glycolysis and understand its underlying control mechanisms, a kinetic model of the pathway was constructed by using the software gepasi. The model was based on the kinetic parameters determined in the purified recombinant enzymes, and the enzyme activities, and steady-state fluxes and metabolite concentrations determined in amoebal trophozoites. The model predicted, with a high degree of accuracy, the flux and metabolite concentrations found in trophozoites, but only when the pyrophosphate concentration was held constant; at variable pyrophosphate, the model was not able to completely account for the ATP production/consumption balance, indicating the importance of the pyrophosphate homeostasis for amoebal glycolysis. Control analysis by the model revealed that hexokinase exerted the highest flux control (73%), as a result of its low cellular activity and strong AMP inhibition. 3-Phosphoglycerate mutase also exhibited significant flux control (65%) whereas the other pathway enzymes showed little or no control. The control of the ATP concentration was also mainly exerted by ATP consuming processes and 3-phosphoglycerate mutase and hexokinase (in the producing block). The model also indicated that, in order to diminish the amoebal glycolytic flux by 50%, it was required to decrease hexokinase or 3-phosphoglycerate mutase by 24% and 55%, respectively, or by 18% for both enzymes. By contrast, to attain the same reduction in flux by inhibiting the pyrophosphate-dependent enzymes pyrophosphate-phosphofructokinase and pyruvate phosphate dikinase, they should be decreased > 70%. On the basis of metabolic control analysis, steps whose inhibition would have stronger negative effects on the energy metabolism of this parasite were identified, thus becoming alternative targets for drug design.

Current therapeutics, their problems, and sulfur-containing-amino-acid metabolism as a novel target against infections by "amitochondriate" protozoan parasites

The "amitochondriate" protozoan parasites of humans Entamoeba histolytica, Giardia intestinalis, and Trichomonas vaginalis share many biochemical features, e.g., energy and amino acid metabolism, a spectrum of drugs for their treatment, and the occurrence of drug resistance. These parasites possess metabolic pathways that are divergent from those of their mammalian hosts and are often considered to be good targets for drug development. Sulfur-containing-amino-acid metabolism represents one such divergent metabolic pathway, namely, the cysteine biosynthetic pathway and methionine gamma-lyase-mediated catabolism of sulfur-containing amino acids, which are present in T. vaginalis and E. histolytica but absent in G. intestinalis. These pathways are potentially exploitable for development of drugs against amoebiasis and trichomoniasis. For instance, L-trifluoromethionine, which is catalyzed by methionine gamma-lyase and produces a toxic product, is effective against T. vaginalis and E. histolytica parasites in vitro and in vivo and may represent a good lead compound. In this review, we summarize the biology of these microaerophilic parasites, their clinical manifestation and epidemiology of disease, chemotherapeutics, the modes of action of representative drugs, and problems related to these drugs, including drug resistance. We further discuss our approach to exploit unique sulfur-containing-amino-acid metabolism, focusing on development of drugs against E. histolytica.

Glycolysis in Entamoeba histolytica. Biochemical characterization of recombinant glycolytic enzymes and flux control analysis

The synthesis of ATP in the human parasite Entamoeba histolytica is carried out solely by the glycolytic pathway. Little kinetic and structural information is available for most of the pathway enzymes. We report here the gene cloning, overexpression and purification of hexokinase, hexose-6-phosphate isomerase, inorganic pyrophosphate-dependent phosphofructokinase, fructose-1,6 bisphosphate aldolase (ALDO), triosephosphate isomerase, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase, phosphoglycerate mutase (PGAM), enolase, and pyruvate phosphate dikinase (PPDK) enzymes from E. histolytica. Kinetic characterization of these 10 recombinant enzymes was made, establishing the kinetic constants at optimal and physiological pH values, analyzing the effect of activators and inhibitors, and investigating the storage stability and oligomeric state. Determination of the catalytic efficiencies at the pH optimum and at pH values that resemble those of the amoebal trophozoites was performed for each enzyme to identify possible controlling steps. This analysis suggested that PGAM, ALDO, GAPDH, and PPDK might be flux control steps, as they showed the lowest catalytic efficiencies. An in vitro reconstruction of the final stages of glycolysis was made to determine their flux control coefficients. Our results indicate that PGAM and PPDK exhibit high control coefficient values at physiological pH.

Functional analysis, overexpression, and kinetic characterization of pyruvate kinase from Plasmodium falciparum

The important role of pyruvate kinase during malarial infection has prompted the cloning of a cDNA encoding Plasmodium falciparum pyruvate kinase (pfPyrK), using mRNA from intraerythrocytic-stage malaria parasites. The full-length cDNA encodes a protein with a computed molecular weight of 55.6 kDa and an isoelectric point of 7.5. The purified recombinant pfPyrK is enzymatically active and exists as a homotetramer in its active form. The enzyme exhibits hyperbolic kinetics with respect to phosphoenolpyruvate and ADP, with K(m) of 0.19 and 0.12 mM, respectively. pfPyrK is not affected by fructose-1,6-bisphosphate, a general activating factor of pyruvate kinase for most species. Glucose-6-phosphate, an activator of the Toxoplasma gondii enzyme, does not affect pfPyrK activity. Similar to rabbit pyruvate kinase, pfPyrK is susceptible to inactivation by 1mM pyridoxal-5'-phosphate, but to a lesser extent. A screen for inhibitors to pfPyrK revealed that it is markedly inhibited by ATP and citrate. Detailed kinetic analysis revealed a transition from hyperbolic to sigmoidal kinetics for PEP in the presence of citrate, as well as competitive inhibitory behavior for ATP with respect to PEP. Citrate exhibits non-competitive inhibition with respect to ADP with a K(i) of 0.8mM. In conclusion, P. falciparum expresses an active pyruvate kinase during the intraerythrocytic-stage of its developmental cycle that may play important metabolic roles during infection.

Kinetic mechanism and metabolic role of pyruvate phosphate dikinase from Entamoeba histolytica

The kinetic mechanism and the metabolic role of pyruvate phosphate dikinase from Entamoeba histolytica were investigated. The initial velocity patterns in double reciprocal plots were parallel for the phosphoenolpyruvate/AMP and phosphoenolpyruvate/pyrophosphate substrate pairs and intersecting for the AMP/pyrophosphate pair. This suggests a kinetic mechanism with two independent reactions. The rate of ATP synthesis at saturating and equimolar concentrations of phosphoenolpyruvate, AMP, and pyrophosphate was inhibited by phosphate, which is consistent with an ordered steady-state mechanism. Enzyme phosphorylation by [(32)P(i)]pyrophosphate depends on the formation of a ternary complex between AMP, pyrophosphate, and pyruvate phosphate dikinase. In consequence, the reaction that involves the AMP/pyrophosphate pair follows a sequential steady-state mechanism. The product inhibition patterns of ATP and phosphate versus phosphoenolpyruvate were noncompetitive and uncompetitive, respectively, suggesting that these products were released in an ordered process (phosphate before ATP). The ordered release of phosphate and ATP and the noncompetitive inhibition patterns of pyruvate versus AMP and versus pyrophosphate also supported the sequential kinetic mechanism between AMP and pyrophosphate. Taken together, our data provide evidence for a uni uni bi bi pingpong mechanism for recombinant pyruvate phosphate dikinase from E. histolytica. The Delta G value for the reaction catalyzed by pyruvate phosphate dikinase (+2.7 kcal/mol) determined under near physiological conditions indicates that the synthesis of ATP is not thermodynamically favorable in trophozoites of E. histolytica.