Investigation of network flexibility and identification of reaction essentiality of human host cells infected by S. flexneri

We have developed a consistently reduced core model for HeLa based on the genome scale human reconstruction. This reduced model includes the central carbon pathways, the electron transport chains, all the necessary transport and exchange reactions and the essential compartments for HeLa. In addition, the pyruvate/glutamate metabolism and the purine/pyrimidine catabolism for HeLa were incorporated based on available experimental observations. The core HeLa model consists of 273 reactions and 202 metabolites. The tFBA analysis of the resulting model, where we integrated experimental information about the metabolites concentrations and fluxes, allowed us to: (i) identify the thermodynamically feasible direction profiles of the system; (ii) classify the reactions according to how close (or far) they operate from their thermodynamic equilibrium. This analysis also shed light on the flexibility of the network in terms of the feasible ranges of the metabolite levels and Gibbs energies of reactions. Furthermore, we used the ORACLE (Optimization and Risk Analysis of Complex Living Entities) framework to develop sets of mechanistic kinetics models. ORACLE allows us, despite incomplete information about kinetic properties of the network, to integrate thermodynamics, and available omics and kinetic data. Using these models, we identified and ranked the metabolites according to their impact on the HeLa physiology when drained from the system. This analysis allowed us to uncover the ‘weak’ and ‘strong’ links in host-pathogen interactions. The identification of such links of the metabolic interface between host and pathogen networks could potentially be validated in vitro and can help us understand the mechanisms of infection and provide insights in the fight against intracellular pathogens.

Kinetic reconstruction and analysis of sphingolipid metabolism

Vassily Hatzimanikatis, Howard Riezman, Georgios Savoglidis

Lipids are major constituents of the cell. They are responsible for major properties of the cellular membranes: hydrophobicity, selective permeability and being the scaffold of signaling proteins. Many diseases are associated with alterations in the lipid distribution in the cell and the composition of membrane domains. Metabolic syndrome, obesity, atherosclerosis, as well as Alzheimer’s, Huntington’s diseases and cancer, have an impact in the levels of lipids, with observed alterations in their concentrations compared to the healthy state. Applying computational techniques along with systematic modeling of lipid metabolism can provide insights that can guide biomedical research and develop potential strategies for prevention and cure. In the present study we developed a comprehensive model of the sphingolipid biosynthesis in the yeast Saccharomyces cerevisiae. Sphingolipids are one of the four major lipid categories, along with (glycero)phospholipids, sterols and fatty acids synthesized in the yeast S. cerevisiae. The importance of sphingolipids present in any higher eukaryote has been demonstrated in many recent studies. For this study, S. cerevisiae has been chosen as a model organism due to its high homology of cellular processes with mammalian cells. The developed model will be an essential part towards the construction of a detailed kinetic model of the whole lipidome of the cell. We first constructed a stoichiometric model that contains all the currently known reactions for the biosynthesis of ceramides and complex sphingolipids in the yeast. Additionally, we have accounted for all five reported hydroxylation states, along with the reactions synthesizing the necessary precursor metabolites from other lipid pathways (i.e. palmitate-CoA from fatty acid synthesis and phosphatidylinositol from the phospholipids metabolism). We next developed a kinetic model and we used a large number of lipidomic measurements of wild type yeast to consistently calibrate our model. Curation of the kinetic information of the model came from the comprehensive mining of references for operation of the enzymes as well as ranges of kinetic parameters from online databases. These datasets created the pool for a sampling technique that accounts for the uncertainty in the parameters of the model. This lead to a robust dynamic model, containing mass balances for all the components of the biochemical network as well as terms that accounted for the dillution of these molecules in the cell due to growth. We performed a thorough kinetic analysis of the system by examining the impact of different assumptions in enzyme operation on the levels of ceramides and complex sphingolipids (e.g. substrate competition, (un)competitive inhibition). By applying the principles of Metabolic Control Analysis (MCA) we were able to quantify the effect of enzyme activities on the lipid profiles and we identified enzymes of the biochemical reaction network as efficient targets for metabolic engineering towards a desired state. We also applied a reverse engineering method and we were able to identify enzyme perturbations responsible for an observed altered state compared to the wild type cellular lipidomic profile. Such an approach could lead to the identification of genetic mutations by imploring the information containd within metabolomic measurements. Although we demonstrate the application of the method in models of lipid metabolism it is possible to be applied to biochemical systems of various parts of metabolism and sizes creating a modular platform for kinetic analysis of cellular operations.

2014

Analysis and Design of Metabolic Engineering and Synthetic Biology Strategies using Large Scale Dynamic Models of Metabolism

Georgios Fengos, Vassily Hatzimanikatis, Ljubisa Miskovic

Systematic analysis for the redirection of carbon flux in metabolite-producing microorganisms requires the comprehensive understanding of their complex metabolic processes. The use of large-scale dynamic models of metabolism plays key role in the understanding of these processes and the study of possible metabolic engineering interventions. However, the generation of such models is hampered by the intrinsic nonlinearities of enzymatic reactions, and the uncertainties at different levels. In particular there is limited knowledge about the exact kinetic mechanisms, and many of the parameters involved in these mechanisms remain largely unknown. In this study we propose a systematic methodology to generate large populations of dynamic non-linear models of metabolism using the ORACLE[1] (Optimization and Risk Analysis of Complex Living Entities) framework. Instead of seeking for an optimal value of the estimated kinetic model parameter values, we integrate thermodynamics, available omics, and kinetic data to construct populations of models that are locally stable, and consistent with the observed physiology. To demonstrate the utility of this methodology we constructed a population of large-scale dynamical models of optimally grown E. coli that involves 283 metabolites and 409 reactions. We used these models to (i) study the response of E. coli metabolism upon large-scale perturbations, such as gene knockouts; (ii) identify and analyze multiple steady states; and (iii) characterize basins of attraction around the identified steady states. The aforementioned analyses provide valuable insight for the design of metabolic engineering strategies towards amplification of desired product-forming pathways. 1. Miskovic, L. & Hatzimanikatis, V. Production of biofuels and biochemicals: in need of an ORACLE. Trends in biotechnology 28, 391–7 (2010).

2016

Metabolic control analysis of the central carbon pathway in optimally grown E. coli

Stefano Andreozzi, Vassily Hatzimanikatis, Ljubisa Miskovic, Keng Cher Soh

The engineering of the metabolism requires the building of reliable kinetic models of the metabolic pathways through the integration of the information from heterogeneous data sources. Unfortunately, the quantitative knowledge about fluxes, mechanisms or kinetic parameters of the constituent reactions of the metabolic network is most of the times uncertain or missing, due to limitations in the measurement techniques of the relevant biological properties of the cell and variability in experimental conditions. In this contribution, we employ ORACLE (Optimization and Risk Analysis of Complex Living Entities), the computational framework for the analysis and optimization of metabolic networks that addresses these issues. ORACLE allows us to generate populations of kinetic models that account explicitly for the physico-chemical and thermodynamic features and constraints of metabolic networks. We show how ORACLE incorporates fluxomics and metabolomics data and constraints into the central carbon core model of E. coli to uncover its possible steady state behaviors, when the cells are cultivated under the physiological condition of optimal growth on different carbon sources, both anaerobically and aerobically. Alternative cases of possible intracellular flux and thermodynamic states, corresponding to the given experimental data, are identified and subsequently the kinetic analysis is performed to evaluate how the control over fluxes and concentrations is distributed over the metabolic network. We demonstrate that ORACLE methodology provides guidance for metabolic engineering in spite of uncertain/missing information in the local/global properties of the metabolic network, integrating a statistical account on the semi-quantitative available data with consistent thermodynamic information and plausible hypotheses.