Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into GEMs was first enabled by the GECKO method, allowing the study of phenotypes constrained by protein limitations. Here, we upgraded the GECKO toolbox in order to enhance models with enzyme and proteomics constraints for any organism with an available GEM reconstruction. With this, enzyme-constrained models (ecModels) for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus were generated, aiming to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions revealed that upregulation and high saturation of enzymes in amino acid metabolism were found to be common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO was further developed by the implementation of an automated framework for continuous and version-controlled update of ecModels, which was validated by producing additional high-quality ecModels for Escherichia coli and Homo sapiens. These efforts aim to facilitate the utilization of ecModels in basic science, metabolic engineering and synthetic biology purposes.

Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into GEMs was first enabled by the GECKO method, allowing the study of phenotypes constrained by protein limitations. Here, we upgraded the GECKO toolbox in order to enhance models with enzyme and proteomics constraints for any organism with an available GEM reconstruction. With this, enzyme-constrained models (ecModels) for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus were generated, aiming to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions revealed that upregulation and high saturation of enzymes in amino acid metabolism were found to be common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO was further developed by the implementation of an automated framework for continuous and version-controlled update of ecModels, which was validated by producing additional high-quality ecModels for Escherichia coli and Homo sapiens. These efforts aim to facilitate the utilization of ecModels in basic science, metabolic engineering and synthetic biology purposes.

Environmental flexibility does not explain metabolic robustness

Cells show remarkable resilience against genetic and environmental perturbations. However, its evolutionary origin remains obscure. In order to leverage methods of systems biology for examining cellular robustness, a computationally accessible way of quantification is needed. Here, we present an unbiased metric of structural robustness in genome-scale metabolic models based on concepts prevalent in reliability engineering and fault analysis. The probability of failure (PoF) is defined as the (weighted) portion of all possible combinations of loss-of-function mutations that disable network functionality. It can be exactly determined if all essential reactions, synthetic lethal pairs of reactions, synthetic lethal triplets of reactions etc. are known. In theory, these minimal cut sets (MCSs) can be calculated for any network, but for large models the problem remains computationally intractable. Herein, we demonstrate that even at the genome scale only the lowest-cardinality MCSs are required to efficiently approximate the PoF with reasonable accuracy. Building on an improved theoretical understanding, we analysed the robustness of 489 E. coli, Shigella, Salmonella, and fungal genome-scale metabolic models (GSMMs). In contrast to the popular “congruence theory”, which explains the origin of genetic robustness as a byproduct of selection for environmental flexibility, we found no correlation between network robustness and the diversity of growth-supporting environments. On the contrary, our analysis indicates that amino acid synthesis rather than carbon metabolism dominates metabolic robustness.

Environmental flexibility does not explain metabolic robustness

Cells show remarkable resilience against genetic and environmental perturbations. However, its evolutionary origin remains obscure. In order to leverage methods of systems biology for examining cellular robustness, a computationally accessible way of quantification is needed. Here, we present an unbiased metric of structural robustness in genome-scale metabolic models based on concepts prevalent in reliability engineering and fault analysis.The probability of failure (PoF) is defined as the (weighted) portion of all possible combinations of loss-of-function mutations that disable network functionality. It can be exactly determined, if all essential reactions, synthetic lethal pairs of reactions, synthetic lethal triplets of reactions etc., are known. In theory, these minimal cut sets (MCSs) can be calculated for any network, but for large models the problem remains computationally intractable. Herein, we demonstrate that even at the genome scale only the lowest-cardinality MCSs are required to efficiently approximate the PoF with reasonable accuracy.We analysed the robustness of 489 E. coli, Shigella, Salmonella, and fungal genome-scale metabolic models (GSMMs). In contrast to the popular “congruence theory”, which explains the origin of genetic robustness as a byproduct of selection for environmental flexibility, we found no correlation between network robustness and the diversity of growth-supporting environments. On the contrary, our analysis indicates that amino acid synthesis rather than carbon metabolism dominates metabolic robustness.

Redox controls metabolic robustness in the gas-fermenting acetogenClostridium autoethanogenum

Living biological systems display a fascinating ability to self-organize their metabolism. This ability ultimately determines the metabolic robustness that is fundamental to controlling cellular behavior. However, fluctuations in metabolism can affect cellular homeostasis through transient oscillations. For example, yeast cultures exhibit rhythmic oscillatory behavior in high cell-density continuous cultures. Oscillatory behavior provides a unique opportunity for quantitating the robustness of metabolism, as cells respond to changes by inherently compromising metabolic efficiency. Here, we quantify the limits of metabolic robustness in self-oscillating autotrophic continuous cultures of the gas-fermenting acetogenClostridium autoethanogenum. Online gas analysis and high-resolution temporal metabolomics showed oscillations in gas uptake rates and extracellular byproducts synchronized with biomass levels. The data show initial growth on CO, followed by growth on CO and H2. Growth on CO and H2results in an accelerated growth phase, after which a downcycle is observed in synchrony with a loss in H2uptake. Intriguingly, oscillations are not linked to translational control, as no differences were observed in protein expression during oscillations. Intracellular metabolomics analysis revealed decreasing levels of redox ratios in synchrony with the cycles. We then developed a thermodynamic metabolic flux analysis model to investigate whether regulation in acetogens is controlled at the thermodynamic level. We used endo- and exo-metabolomics data to show that the thermodynamic driving force of critical reactions collapsed as H2uptake is lost. The oscillations are coordinated with redox. The data indicate that metabolic oscillations in acetogen gas fermentation are controlled at the thermodynamic level.

Genome-Scale Metabolic Reconstruction and in Silico Perturbation Analysis of the Polar Diatom Fragilariopsis cylindrus Predicts High Metabolic Robustness

Diatoms are major primary producers in polar environments where they can actively grow under extremely variable conditions. Integrative modeling using a genome-scale model (GSM) is a powerful approach to decipher the complex interactions between components of diatom metabolism and can provide insights into metabolic mechanisms underlying their evolutionary success in polar ecosystems. We developed the first GSM for a polar diatom, Fragilariopsis cylindrus, which enabled us to study its metabolic robustness using sensitivity analysis. We find that the predicted growth rate was robust to changes in all model parameters (i.e., cell biochemical composition) except the carbon uptake rate. Constraints on total cellular carbon buffer the effect of changes in the input parameters on reaction fluxes and growth rate. We also show that single reaction deletion of 20% to 32% of active (nonzero flux) reactions and single gene deletion of 44% to 55% of genes associated with active reactions affected the growth rate, as well as the production fluxes of total protein, lipid, carbohydrate, DNA, RNA, and pigments by less than 1%, which was due to the activation of compensatory reactions (e.g., analogous enzymes and alternative pathways) with more highly connected metabolites involved in the reactions that were robust to deletion. Interestingly, including highly divergent alleles unique for F. cylindrus increased its metabolic robustness to cellular perturbations even more. Overall, our results underscore the high robustness of metabolism in F. cylindrus, a feature that likely helps to maintain cell homeostasis under polar conditions.

Lactate and glycerol-3-phosphate metabolism cooperatively regulate growth and redox balance during Drosophila melanogaster larval development

ABSTRACTThe dramatic growth that occurs during Drosophila larval development requires rapid conversion of nutrients into biomass. Many larval tissues respond to these biosynthetic demands by increasing carbohydrate metabolism and lactate dehydrogenase (dLDH) activity. The resulting metabolic program is ideally suited to synthesize macromolecules and mimics the manner by which cancer cells rely on aerobic glycolysis. To explore the potential role of Drosophila dLDH in promoting biosynthesis, we examined how dLdh mutations influence larval development. Our studies unexpectantly found that dLdh mutants grow at a normal rate, indicating that dLDH is dispensable for larval biomass production. However, subsequent metabolomic analyses suggested that dLdh mutants compensate for the inability to produce lactate by generating excess glycerol-3-phosphate (G3P), the production of which also influences larval redox balance. Consistent with this possibility, larvae lacking both dLDH and G3P dehydrogenase (GPDH1) exhibit developmental delays, synthetic lethality, and aberrant carbohydrate metabolism. Considering that human cells also generate G3P upon Lactate Dehydrogenase A (LDHA) inhibition, our findings hint at a conserved mechanism in which the coordinate regulation of lactate and G3P synthesis imparts metabolic robustness upon growing animal tissues.