The metabolic and epigenetic effects of the isocitrate dehydrogenase-1 mutation in glioma

Cancer cells have been noted to have an altered metabolic phenotype for over ninety years. In the presence of oxygen, differentiated cells predominately utilise the tricarboxylic acid (TCA) cycle and oxidative phosphorylation to efficiently produce energy and the metabolites necessary for protein and lipid synthesis. However, in hypoxia, this process is altered and cells switch to a higher rate of glycolysis and lactate production to maintain their energy and metabolic needs. In cancer cells, glycolysis is maintained at a high rate, even in the presence of oxygen; a term described as “aerobic glycolysis”. Tumour cells are rapidly dividing and have a much greater need for anabolism compared to normal differentiated cells. Rapid glucose metabolism enables faster ATP production as well as a greater redistribution of carbons to nucleotide, protein, and fatty acid synthesis, thus maximising cell growth. Recently, other metabolic changes, driven by mutations in genes related to the TCA cycle, indicate an alternative role for metabolism in cancer, the “oncometabolite”. This is where a particular metabolite builds up within the cell and contributes to the tumorigenic process. One of these genes is isocitrate dehydrogenase (IDH) IDH is an enzyme that forms part of the tricarboxylic acid (TCA) cycle and converts isocitrate to α-ketoglutarate (α-KG). It exists in three isoforms; IDH1, IDH2 and IDH3 with the former present in the cytoplasm and the latter two in the mitochondria. Point mutations have been identified in the IDH1 and IDH2 genes in glioma which result in a gain of function by converting α-KG to 2-hydroxyglutarate (2HG), an oncometabolite. 2HG acts as a competitive inhibitor of the α-KG dependent dioxygenases, a superfamily of enzymes that are involved in numerous cellular processes such as DNA and histone demethylation. It was hypothesised that the IDH1 mutation would result in other metabolic changes in the cell other than 2HG production, and could potentially identify pathways which could be targeted for therapeutic treatment. In addition, 2HG can act as a potential competitive inhibitor of α-KG dependent dioxygenases, so it was hypothesised that there would be an effect on histone methylation. This may alter gene expression and provide a mechanism for tumourogenesis and potentially identify further therapeutic targets. Metabolic analysis of clinical tumour samples identified changes associated with the IDH1 mutation, which included a reduction in α-KG and an increase in GABA, in addition to the increase in 2HG. This was replicated in several cell models, where 13C labelled metabolomics was also used to identify a possible increase in metabolic flux from glutamate to GABA, as well as from α-KG to 2HG. This may provide a mechanism whereby the cell can bypass the IDH1 mutation as GABA can be metabolised to succinate in the mitochondria by GABA transaminase via the GABA shunt. JMJ histone demethylases are a subset of the α-KG dependent dioxygenases, and are involved in removing methyl groups from histone tails. Changes in histone methylation are associated with changes in gene expression depending on the site and extent of chemical modification. To identify whether the increase in 2HG and fall in α-KG was associated with inhibition of histone demethylases a histone methylation screen was used. The IDH1 mutation was associated with an increase in methylation of H3K4, which is associated with gene activation. ChiP and RNA sequencing identified an increase in H3K4me3 at the transcription start site of the GABRB3 subunit, resulting in an increase in gene expression. The GABRB3 subunit forms part of the GABA-A receptor, a chloride channel, which on activation can reduce cell proliferation. The IDH1 mutation was associated with an increase in GABA and GABRB3 subunit of the GABA-A receptor. This raises the possibility of GABA transaminase as a potential therapeutic target. Inhibition of this enzyme could reduce GABA metabolism, potentially reducing any beneficial effect of the GABA shunt in IDH1 mutant tumours, and increasing activation of the GABA-A receptor by increasing the concentration of GABA in the brain. This in turn may reduce cell proliferation, and could be achieved by using Vigabatrin, a GABA transaminase inhibitor licensed for use in epilepsy.

Unravelling metabolism of Leishmania by metabolomics

The leishmaniases are neglected tropical diseases with an urgent need for effective drugs. Better understanding of the metabolism of the causative parasites will hopefully lead to development of new compounds targeted at critical points of the parasite’s biochemical pathways. In my work I focused on the pentose phosphate pathway of Leishmania, specifically on transketolase, sugar utilisation, and comparison between insect and mammalian infective stages of the parasites. The pentose phosphate pathway (PPP) is the major cellular source of NADPH, an agent critical for oxidative stress defence. The PPP uses glucose, reduces the NADP+ cofactor and produces various sugar phosphates by mutual interconversions. One of the enzymes involved in this latter part is transketolase (TKT). A Leishmania mexicana cell line deleted in transketolase (Δtkt) was assessed regarding viability, sensitivity to a range of drugs, changes in metabolism, and infectivity. The Δtkt cell line had no obvious growth defect in the promastigote stage, but it was more sensitive to an oxidative stress inducing agent and most of the drugs tested. Most importantly, the Δtkt cells were not infective to mice, establishing TKT as a new potential drug target. Metabolomic analyses revealed multiple changes as a consequence of TKT deletion. Levels of the PPP intermediates upstream of TKT increased substantially, and were diverted into additional reactions. The perturbation triggered further changes in metabolism, resembling the ‘stringent metabolic response’ of amastigotes. The Δtkt cells consumed less glucose and glycolytic intermediates were decreased indicating a decrease in flux, and metabolic end products were diminished in production. The decrease in glycolysis was possibly caused by inhibition of fructose-1,6-bisphosphate aldolase by accumulation of the PPP intermediates 6-phosphogluconate and ribose 5-phosphate. The TCA cycle was fuelled by alternative carbon sources, most likely amino acids, instead of glucose. It remains unclear why deletion of TKT is lethal for amastigotes, increased sensitivity to oxidative stress or drop in mannogen levels may contribute, but no definite conclusions can be made. TKT localisation indicated interesting trends too. The WT enzyme is present in the cytosol and glycosomes, whereas a mutant version, truncated by ten amino acids, but retaining a C-terminal targeting sequence, localised solely to glycosomes. Surprisingly, cells expressing purely cytosolic or glycosomal TKT did not have different phenotypes regarding growth, oxidative stress sensitivity or any detected changes in metabolism. Hence, control of the subcellular localisation remains unclear as well as its function. However, these data are in agreement with the presumed semipermeable nature of the glycosome. Further, L. mexicana promastigote cultures were grown in media with different combinations of labelled glucose and ribose and their incorporation into metabolism was followed. Glucose was the preferred carbon source, but when not available, it could be fully replaced with ribose. I also compared metabolic profiles from splenic amastigotes, axenic amastigotes and promastigotes of L. donovani. Metabolomic analysis revealed a substantial drop in amino acids and other indications coherent with a stringent metabolic response in amastigotes. Despite some notable differences, axenic and splenic amastigotes demonstrated fairly similar results both regarding the total metabolic profile and specific metabolites of interest.

Filling in gaps of Drosophila melanogaster urate degradation metabolic pathway using metabolomics approaches: towards the core metabolome of the fruit fly

The primary goal of systems biology is to integrate complex omics data, and data obtained from traditional experimental studies in order to provide a holistic understanding of organismal function. One way of achieving this aim is to generate genome-scale metabolic models (GEMs), which contain information on all metabolites, enzyme-coding genes, and biochemical reactions in a biological system. Drosophila melanogaster GEM has not been reconstructed to date. Constraint-free genome-wide metabolic model of the fruit fly has been reconstructed in our lab, identifying gaps, where no enzyme was identified and metabolites were either only produced or consume. The main focus of the work presented in this thesis was to develop a pipeline for efficient gap filling using metabolomics approaches combined with standard reverse genetics methods, using 5-hydroxyisourate hydrolase (5-HIUH) as an example. 5-HIUH plays a role in urate degradation pathway. Inability to degrade urate can lead to inborn errors of metabolism (IEMs) in humans, including hyperuricemia. Based on sequence analysis Drosophila CG30016 gene was hypothesised to encode 5- HIUH. CG30016 knockout flies were examined to identify Malpighian tubules phenotype, and shortened lifespan might reflect kidney disorders in hyperuricemia in humans. Moreover, LC-MS analysis of mutant tubules revealed that CG30016 is involved in purine metabolism, and specifically urate degradation pathway. However, the exact role of the gene has not been identified, and the complete method for gap filling has not been developed. Nevertheless, thanks to the work presented here, we are a step closer towards the development of a gap-filling pipeline in Drosophila melanogaster GEM. Importantly, the areas that require further optimisation were identified and are the focus of future research. Moreover, LC-MS analysis confirmed that tubules rather than the whole fly were more suitable for metabolomics analysis of purine metabolism. Previously, Dow/Davies lab has generated the most complete tissue-specific transcriptomic atlas for Drosophila – FlyAtlas.org, which provides data on gene expression across multiple tissues of adult fly and larva. FlyAtlas revealed that transcripts of many genes are enriched in specific Drosophila tissues, and that it is possible to deduce the functions of individual tissues within the fly. Based on FlyAtlas data, it has become clear that the fly (like other metazoan species) must be considered as a set of tissues, each 2 with its own distinct transcriptional and functional profile. Moreover, it revealed that for about 30% of the genome, reverse genetic methods (i.e. mutation in an unknown gene followed by observation of phenotype) are only useful if specific tissues are investigated. Based on the FlyAtlas findings, we aimed to build a primary tissue-specific metabolome of the fruit fly, in order to establish whether different Drosophila tissues have different metabolomes and if they correspond to tissue-specific transcriptome of the fruit fly (FlyAtlas.org). Different fly tissues have been dissected and their metabolome elucidated using LC-MS. The results confirmed that tissue metabolomes differ significantly from each other and from the whole fly, and that some of these differences can be correlated to the tissue function. The results illustrate the need to study individual tissues as well as the whole organism. It is clear that some metabolites that play an important role in a given tissue might not be detected in the whole fly sample because their abundance is much lower in comparison to other metabolites present in all tissues, which prevent the detection of the tissue-specific compound.