Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

Metabolic flux analysis for directing metabolism of Streptomyces olindensis from growth to anti-tumor drug (cosmomycin) production

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, CAPES, Brazil

Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

Magnetic flux analysis in superconducting Fault Current Limiters under short-circuit conditions

In order to obtain a more compact Superconducting Fault Current limiter (SFCL), a special geometry of core and AC coil is required. This results in a unique magnetic flux pattern which differs from those associated with conventional round core arrangements. In this paper the magnetic flux density within a Fault Current Limiter (FCL) is described. Both experimental and analytical approaches are considered. A small scale prototype of an FCL was constructed in order to conduct the experiments. This prototype comprises a single phase. The analysis covers both the steady state and the short-circuit condition. Simulation results were obtained using commercial software based on the Finite Element Method (FEM). The magnetic flux saturating the cores, leakage magnetic flux giving rise to electromagnetic forces and leakage magnetic flux flowing in the enclosing tank are computed.

Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

Transcriptome and proteome analysis of xylose-metabolising Saccharomyces cerevisiae

Increasing concern about global climate warming has accelerated research into renewable energy sources that could replace fossil petroleum-based fuels and materials. Bioethanol production from cellulosic biomass by fermentation with baker s yeast Saccharomyces cerevisiae is one of the most studied areas in this field. The focus has been on metabolic engineering of S. cerevisiae for utilisation of the pentose sugars, in particular D-xylose that is abundant in the hemicellulose fraction of biomass. Introduction of a heterologous xylose-utilisation pathway into S. cerevisiae enables xylose fermentation, but ethanol yield and productivity do not reach the theoretical level. In the present study, transcription, proteome and metabolic flux analyses of recombinant xylose-utilising S. cerevisiae expressing the genes encoding xylose reductase (XR) and xylitol dehydrogenase (XDH) from Pichia stipitis and the endogenous xylulokinase were carried out to characterise the global cellular responses to metabolism of xylose. The aim of these studies was to find novel ways to engineer cells for improved xylose fermentation. The analyses were carried out from cells grown on xylose and glucose both in batch and chemostat cultures. A particularly interesting observation was that several proteins had post-translationally modified forms with different abundance in cells grown on xylose and glucose. Hexokinase 2, glucokinase and both enolase isoenzymes 1 and 2 were phosphorylated differently on the two different carbon sources studied. This suggests that phosphorylation of glycolytic enzymes may be a yet poorly understood means to modulate their activity or function. The results also showed that metabolism of xylose affected the gene expression and abundance of proteins in pathways leading to acetyl-CoA synthesis and altered the metabolic fluxes in these pathways. Additionally, the analyses showed increased expression and abundance of several other genes and proteins involved in cellular redox reactions (e.g. aldo-ketoreductase Gcy1p and 6-phosphogluconate dehydrogenase) in cells grown on xylose. Metabolic flux analysis indicated increased NADPH-generating flux through the oxidative part of the pentose phosphate pathway in cells grown on xylose. The most importantly, results indicated that xylose was not able to repress to the same extent as glucose the genes of the tricarboxylic acid and glyoxylate cycles, gluconeogenesis and some other genes involved in the metabolism of respiratory carbon sources. This suggests that xylose is not recognised as a fully fermentative carbon source by the recombinant S. cerevisiae that may be one of the major reasons for the suboptimal fermentation of xylose. The regulatory network for carbon source recognition and catabolite repression is complex and its functions are only partly known. Consequently, multiple genetic modifications and also random approaches would probably be required if these pathways were to be modified for further improvement of xylose fermentation by recombinant S. cerevisiae strains.

Energy metabolism, enzymatic flux capacities, and metabolic flux rates in flying honeybees.

Honeybees rely primarily on the oxidation of hexose sugars to provide the energy required for flight. Measurement of VCO2 (equal to VO2, because VCO2/VO2 = 1.0 during carbohydrate oxidation) during flight allowed estimation of steady-state flux rates through pathways of flight muscle energy metabolism. Comparison of Vmax values for flight muscle hexokinase, phosphofructokinase, citrate synthase, and cytochrome c oxidase with rates of carbon and O2 flux during flight reveal that these enzymes operate closer to Vmax in the flight muscles of flying honeybees than in other muscles previously studied. Possible mechanistic and evolutionary implications of these findings are discussed.

Constraint-based modelling of mixed microbial populations: Application to polyhydroxyalkanoates production

Dissertação para obtenção do Grau de Doutor em Engenharia Química e Bioquímica

Análise de fluxos metabólicos para otimização da síntese do antibiótico cosmomicina por Streptomyces olindensis ICB20

Metabolic flux analysis (MFA) is a powerful tool for analyzing cellular metabolism. In order to control the growth conditions of a specific organism, it is important to have a complete understanding of its MFA. This would allowed us to improve the processes for obtaining products of interest to human and also to understand how to manipulate the genome of a cell, allowing optimization process for genetic engineering. Streptomyces olindensis ICB20 is a promising producer of the antibiotic cosmomycin, a powerful antitumor drug. Several Brazilian researchers groups have been developing studies in order to optimize cosmomycin production in bioreactors. However, to the best of our knowledge, nothing has been done on metabolic fluxes analysis field. Therefore, the aim of this work is to identify several factors that can affect the metabolism of Streptomyces olindensis ICB20, through the metabolic flux analysis. As a result, the production of the secondary metabolite, cosmomycin, can be increased. To achieve this goal, a metabolic model was developed which simulates a distribution of internal cellular fluxes based on the knowledge of metabolic pathways, its interconnections, as well as the constraints of microorganism under study. The validity of the proposed model was verified by comparing the computational data obtained by the model with the experimental data obtained from the literature. Based on the analysis of intracellular fluxes, obtained by the model, an optimal culture medium was proposed. In addition, some key points of the metabolism of Streptomyces olindensis were identified, aiming to direct its metabolism to a greater cosmomycin production. In this sense it was found that by increasing the concentration of yeast extract, the culture medium could be optimized. Furthermore, the inhibition of the biosynthesis of fatty acids was found to be a interesting strategy for genetic manipulation. Based on the metabolic model, one of the optimized medium conditions was experimentally tested in order to demonstrate in vitro what was obtained in silico. It was found that by increasing the concentration of yeast extract in the culture medium would induce to an increase of the cosmomycin production

Early nutritional programming in fish: tailoring the metabolic use of dietary carbohydrates

Tese de doutoramento, Ciências da Vida, do Mar, da Terra e do Ambiente (Nutrição), Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

Enzymatic phosphorylation of muscle glycogen synthase: a mechanism for maintenance of metabolic homeostasis.

We recently analyzed experimental studies of mammalian muscle glycogen synthesis using metabolic control analysis and concluded that glycogen synthase (GSase) does not control the glycogenic flux but rather adapts to the flux which is controlled bv the activity of the proximal glucose transport and hexokinase steps. This model did not provide a role for the well established relationship between GSase fractional activity, determined by covalent phosphorylation, and the rate of glycogen synthesis. Here we propose that the phosphorylation of GSase, which alters the sensitivity to allosteric activation by glucose 6-phosphate (G6P), is a mechanism for controlling the concentration of G6P instead of controlling the flux. When the muscle cell is exposed to conditions which favor glycogen synthesis such as high plasma insulin and glucose concentrations the fractional activity of GSase is increased in coordination with increases in the activity of glucose transport and hexokinase. This increase in GSase fractional activity helps to maintain G6P homeostasis by reducing the G6P concentration required to activate GSase allosterically to match the flux determined by the proximal reactions. This role for covalent phosphorylation also provides a novel solution to the Kacser and Acarenza paradigm which requires coordinated activity changes of the enzymes proximal and distal to a shared intermediate, to avoid unwanted flux changes.

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.