Molecular basis for the reverse reaction of African human trypanosomes glycerol kinase

Crystal structure of glycerol kinase from Trypanosoma cruzi, a potential molecular target in Chagas disease

Chagas disease is a neglected tropical disease caused by the protozoan parasite Trypanosoma cruzi. It bears a significant global health burden with limited treatment options, thus calling for the development of new and effective drugs. Certain trypanosomal metabolic enzymes have been suggested to be druggable and valid for subsequent inhibition. In this study, the crystal structure of glycerol kinase from T. cruzi, a key enzyme in glycerol metabolism in this parasite, is presented. Structural analysis allowed a detailed description of the glycerol binding pocket, while comparative assessment pinpointed a potential regulatory site which may serve as a target for selective inhibition. These findings advance the understanding of glycerol metabolism in eukaryotes and provide a solid basis for the future treatment of Chagas disease.

Glycerol biosynthetic pathway plays an essential role in proliferation and antioxidative defense in the human enteric protozoan parasite Entamoeba histolytica

Amebiasis is caused by the protozoan parasite Entamoeba histolytica. Treatment options other than metronidazole and its derivatives are few, and their low efficacy against asymptomatic cyst carriers, and experimental evidence of resistance in vitro justify the discovery/repurposing campaign for new drugs against amebiasis. Global metabolic responses to oxidative stress and cysteine deprivation by E. histolytica revealed glycerol metabolism may represent a rational target for drug development. In this study using 14C-labelled glucose, only 11% of the total glucose taken up by E. histolytica trophozoites is incorporated to lipids. To better understand the role of glycerol metabolism in this parasite, we focused on characterizing two important enzymes, glycerol kinase (GK) and glycerol 3-phosphate dehydrogenase (G3PDH). Recombinant GK was biochemically characterized in detail, while G3PDH was not due to failure of protein expression and purification. GK revealed novel characteristics and unprecedented kinetic properties in reverse reaction. Gene silencing revealed that GK is essential for optimum growth, whereas G3PDH is not. Gene silencing of G3PDH caused upregulated GK expression, while that of GK resulted in upregulation of antioxidant enzymes as shown by RNA-seq analysis. Although the precise molecular link between GK and the upregulation of antioxidant enzymes was not demonstrated, the observed increase in antioxidant enzyme expression upon GK gene silencing suggests a potential connection between GK and the cellular response to oxidative stress. Together, these results provide the first direct evidence of the biological importance and coordinated regulation of the glycerol metabolic pathways for proliferation and antioxidative defense in E. histolytica, justifying the exploitation of these enzymes as future drug targets.

Targeted protein degradation might present a novel therapeutic approach in the fight against African trypanosomiasis

African trypanosomiasis (AT) is a hemoparasitic disease caused by infection with African trypanosomes and it is prevalent in many sub-Saharan African countries, affecting both humans and domestic animals. The disease is transmitted mostly by haematophagous insects of the genus Glossina while taking blood meal, in the process spreading the parasites from an infected animal to an uninfected animal. The disease is fatal if untreated, and the available drugs are generally ineffective and resulting in toxicities. Therefore, it is still pertinent to explore novel methods and targets for drug discovery. Proteolysis-targeting chimeras (PROTACs) present a new strategy for development of therapeutic molecules that mimic cellular proteasomal-mediated protein degradation to target proteins involved in different disease types. PROTACs have been used to degrade proteins involved in various cancers, neurodegenerative diseases, and immune disorders with remarkable success. Here, we highlight the problems associated with the current treatments for AT, discuss the concept of PROTACs and associated targeted protein degradation (TPD) approaches, and provide some insights on the future potential for the use of these emerging technologies (PROTACs and TPD) for the development of new generation of anti-Trypanosoma drugs and the first "TrypPROTACs".

Carbohydrate metabolism in trypanosomatids: New insights revealing novel complexity, diversity and species-unique features

The human pathogenic trypanosomatid species collectively called the "TriTryp parasites" - Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp. - have complex life cycles, with each of these parasitic protists residing in a different niche during their successive developmental stages where they encounter diverse nutrients. Consequently, they adapt their metabolic network accordingly. Yet, throughout the life cycles, carbohydrate metabolism - involving the glycolytic, gluconeogenic and pentose-phosphate pathways - always plays a central role in the biology of these parasites, whether the available carbon and free energy sources are saccharides, amino acids or lipids. In this paper, we provide an updated review of the carbohydrate metabolism of the TriTryps, highlighting new data about this metabolic network, the interconnection of its pathways and the compartmentalisation of its enzymes within glycosomes, cytosol and mitochondrion. Differences in the expression of the branches of the metabolic network between the successive life-cycle stages of each of these parasitic trypanosomatids are discussed, as well as differences between them. Recent structural and kinetic studies have revealed unique regulatory mechanisms for some of the network's key enzymes with important species-specific variations. Furthermore, reports of multiple post-translational modifications of trypanosomal glycolytic enzymes suggest that additional mechanisms for stage- and/or environmental cues that regulate activity are operational in the parasites. The detailed comparison of the carbohydrate metabolism of the TriTryps has thus revealed multiple differences and a greater complexity, including for the reduced metabolic network in bloodstream-form T. brucei, than previously appreciated. Although these parasites are related, share many cytological and metabolic features and are grouped within a single taxonomic family, the differences highlighted in this review reflect their separate evolutionary tracks from a common ancestor to the extant organisms. These differences are indicative of their adaptation to the different insect vectors and niches occupied in their mammalian hosts.

Structural Characterization of Glycerol Kinase from the Thermophilic Fungus Chaetomium thermophilum

Glycerol is an organic compound that can be utilized as an alternative source of carbon by various organisms. One of the ways to assimilate glycerol by the cell is the phosphorylative catabolic pathway in which its activation is catalyzed by glycerol kinase (GK) and glycerol-3-phosphate (G3P) is formed. To date, several GK crystal structures from bacteria, archaea, and unicellular eukaryotic parasites have been solved. Herein, we present a series of crystal structures of GK from Chaetomium thermophilum (CtGK) in apo and glycerol-bound forms. In addition, we show the feasibility of an ADP-dependent glucokinase (ADPGK)-coupled enzymatic assay to measure the CtGK activity. New structures described in our work provide structural insights into the GK catalyzed reaction in the filamentous fungus and set the foundation for understanding the glycerol metabolism in eukaryotes.

Discovery of trypanocidal coumarins with dual inhibition of both the glycerol kinase and alternative oxidase of Trypanosoma brucei brucei

African trypanosomiasis, sleeping sickness in humans or nagana in animals, is a potentially fatal neglected tropical disease and a threat to 65 million human lives and 100 million small and large livestock animals in sub-Saharan Africa. Available treatments for this devastating disease are few and have limited efficacy, prompting the search for new drug candidates. Simultaneous inhibition of the trypanosomal glycerol kinase (TGK) and trypanosomal alternative oxidase (TAO) is considered a validated strategy toward the development of new drugs. Our goal is to develop a TGK-specific inhibitor for coadministration with ascofuranone (AF), the most potent TAO inhibitor. Here, we report on the identification of novel compounds with inhibitory potency against TGK. Importantly, one of these compounds (compound 17) and its derivatives (17a and 17b) killed trypanosomes even in the absence of AF. Inhibition kinetics revealed that derivative 17b is a mixed-type and competitive inhibitor for TGK and TAO, respectively. Structural data revealed the molecular basis of this dual inhibitory action, which, in our opinion, will aid in the successful development of a promising drug to treat trypanosomiasis. Although the EC₅₀ of compound 17b against trypanosome cells was 1.77 µM, it had no effect on cultured human cells, even at 50 µM.-Balogun, E. O., Inaoka, D. K., Shiba, T., Tsuge, C., May, B., Sato, T., Kido, Y., Nara, T., Aoki, T., Honma, T., Tanaka, A., Inoue, M., Matsuoka, S., Michels, P. A. M., Watanabe, Y.-I., Moore, A. L., Harada, S., Kita, K. Discovery of trypanocidal coumarins with dual inhibition of both the glycerol kinase and alternative oxidase of Trypanosoma brucei brucei.

Insights into the ubiquinol/dioxygen binding and proton relay pathways of the alternative oxidase

The alternative oxidase (AOX) is a monotopic diiron carboxylate protein which catalyzes the four-electron reduction of dioxygen to water by ubiquinol. Although we have recently determined the crystal structure of Trypanosoma brucei AOX (TAO) in the presence and absence of ascofuranone (AF) derivatives (which are potent mixed type inhibitors) the mechanism by which ubiquinol and dioxygen binds to TAO remain inconclusive. In this article, ferulenol was identified as the first competitive inhibitor of AOX which has been used to probe the binding of ubiquinol. Surface plasmon resonance reveals that AF is a quasi-irreversible inhibitor of TAO whilst ferulenol binding is completely reversible. The structure of the TAO-ferulenol complex, determined at 2.7 Å, provided insights into ubiquinol binding and has also identified a potential dioxygen molecule bound in a side-on conformation to the diiron center for the first time.

The kinetic characteristics of human and trypanosomatid phosphofructokinases for the reverse reaction

Eukaryotic ATP-dependent phosphofructokinases (PFKs) are often considered unidirectional enzymes catalysing the transfer of a phospho moiety from ATP to fructose 6-phosphate to produce ADP and fructose 1,6-bisphosphate. The reverse reaction is not generally considered to occur under normal conditions and has never been demonstrated for any eukaryotic ATP-dependent PFKs, though it does occur in inorganic pyrophosphate-dependent PFKs and has been experimentally shown for bacterial ATP-dependent PFKs. The evidence is provided via two orthogonal assays that all three human PFK isoforms can catalyse the reverse reaction in vitro, allowing determination of kinetic properties. Additionally, the reverse reaction was shown possible for PFKs from three clinically important trypanosomatids; these enzymes are contained within glycosomes in vivo. This compartmentalisation may facilitate reversal, given the potential for trypanosomatids to have an altered ATP/ADP ratio in glycosomes compared with the cytosol. The kinetic properties of each trypanosomatid PFK were determined, including the response to natural and artificial modulators of enzyme activity. The possible physiological relevance of the reverse reaction in trypanosomatid and human PFKs is discussed.

Molecular and biochemical characterization of natural and recombinant phosphoglycerate kinase B from Trypanosoma rangeli

T. rangeli epimastigotes contain only a single detectable phosphoglycerate kinase (PGK) enzyme in their cytosol. Analysis of this parasite's recently sequenced genome showed a gene predicted to code for a PGK with the same molecular mass as the natural enzyme, and with a cytosolic localization as well. In this work, we have partially purified the natural PGK from T. rangeli epimastigotes. Furthermore, we cloned the predicted PGK gene and expressed it as a recombinant active enzyme. Both purified enzymes were kinetically characterized and displayed similar substrate affinities, with Km_ATP values of 0.13 mM and 0.5 mM, and Km_3PGA values of 0.28 mM and 0.71 mM, for the natural and recombinant enzyme, respectively. The optimal pH for activity of both enzymes was in the range of 8-10. Like other PGKs, TrPGK is monomeric with a molecular mass of approximately 44 kDa. The enzyme's kinetic characteristics are comparable with those of cytosolic PGK isoforms from related trypanosomatid species, indicating that, most likely, this enzyme is equivalent with the PGKB that is responsible for generating ATP in the cytosol of other trypanosomatids. This is the first report of a glycolytic enzyme characterization from T. rangeli.

Glycerol kinase of African trypanosomes possesses an intrinsic phosphatase activity

BACKGROUND

In general, glycerol kinases (GKs) are transferases that catalyze phospho group transfer from ATP to glycerol, and the mechanism was suggested to be random bi-bi. The reverse reaction i.e. phospho transfer from glycerol 3-phosphate (G3P) to ADP is only physiologically feasible by the African trypanosome GK. In contrast to other GKs the mechanism of Trypanosoma brucei gambiense glycerol kinase (TbgGK) was shown to be in an ordered fashion, and proceeding via autophosphorylation. From the unique reaction mechanism of TbgGK, we envisaged its potential to possess phosphatase activity in addition to being a kinase.

METHODS

Our hypothesis was tested by spectrophotometric and LC-MS/MS analyses using paranitrophenyl phosphate (pNPP) and TbgGK's natural substrate, G3P respectively. Furthermore, protein X-ray crystallography and site-directed mutagenesis were performed to examine pNPP binding, catalytic residues, and the possible reaction mechanism.

RESULTS

In addition to its widely known and expected phosphotransferase (class II) activity, TbgGK can efficiently facilitate the hydrolytic cleavage of phosphoric anhydride bonds (a class III property). This phosphatase activity followed the classical Michaelis-Menten pattern and was competitively inhibited by ADP and G3P, suggesting a common catalytic site for both activities (phosphatase and kinase). The structure of the TGK-pNPP complex, and structure-guided mutagenesis implicated T276 to be important for the catalysis. Remarkably, we captured a crystallographic molecular snapshot of the phosphorylated T276 reaction intermediate.

CONCLUSION

We conclude that TbgGK has both kinase and phosphatase activities.

GENERAL SIGNIFICANCE

This is the first report on a bifunctional kinase/phosphatase enzyme among members of the sugar kinase family.

Improved Metabolic Models for E. coli and Mycoplasma genitalium from GlobalFit, an Algorithm That Simultaneously Matches Growth and Non-Growth Data Sets

Constraint-based metabolic modeling methods such as Flux Balance Analysis (FBA) are routinely used to predict the effects of genetic changes and to design strains with desired metabolic properties. The major bottleneck in modeling genome-scale metabolic systems is the establishment and manual curation of reliable stoichiometric models. Initial reconstructions are typically refined through comparisons to experimental growth data from gene knockouts or nutrient environments. Existing methods iteratively correct one erroneous model prediction at a time, resulting in accumulating network changes that are often not globally optimal. We present GlobalFit, a bi-level optimization method that finds a globally optimal network, by identifying the minimal set of network changes needed to correctly predict all experimentally observed growth and non-growth cases simultaneously. When applied to the genome-scale metabolic model of Mycoplasma genitalium, GlobalFit decreases unexplained gene knockout phenotypes by 79%, increasing accuracy from 87.3% (according to the current state-of-the-art) to 97.3%. While currently available computers do not allow a global optimization of the much larger metabolic network of E. coli, the main strengths of GlobalFit are already played out when considering only one growth and one non-growth case simultaneously. Application of a corresponding strategy halves the number of unexplained cases for the already highly curated E. coli model, increasing accuracy from 90.8% to 95.4%.

Mathematical models that aim to describe the complete metabolism of a cell help us understand cellular metabolic capabilities and evolution, and aid the biotechnological design of microbial strains with desired properties. Draft models are frequently improved through adjustments that increase the agreement of growth/non-growth predictions with observations from gene knockout experiments. Automated methods for this task typically correct one erroneous prediction after the other. We present GlobalFit, a novel method that can consider all experiments and all possible changes simultaneously to identify model modifications that are globally optimal (i.e., that correct the largest possible number of wrong predictions while introducing sets of changes that are most compatible with existing knowledge). This becomes computationally very hard when considering large metabolic models; however, a reduced application of GlobalFit that only looks at small subsets of experiments simultaneously works very well in practice. Allowing only changes that are conservative (e.g., introducing new reactions only if supported by significant genomic evidence), GlobalFit halves the number of wrong growth/non-growth predictions for the state-of-the-art metabolic models of E. coli and Mycoplasma genitalium, increasing prediction accuracy to 95.4% and 93.0%, respectively. By additionally allowing less conservative changes, we are able to improve accuracy further to 97.3% for the M. genitalium model.