Comparison and analysis of objective functions in flux balance analysis

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. The advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

Integration of graph neural networks and genome-scale metabolic models for predicting gene essentiality

Genome-scale metabolic models are powerful tools for understanding cellular physiology. Flux balance analysis (FBA), in particular, is an optimization-based approach widely employed for predicting metabolic phenotypes. In model microbes such as Escherichia coli, FBA has been successful at predicting essential genes, i.e. those genes that impair survival when deleted. A central assumption in this approach is that both wild type and deletion strains optimize the same fitness objective. Although the optimality assumption may hold for the wild type metabolic network, deletion strains are not subject to the same evolutionary pressures and knock-out mutants may steer their metabolism to meet other objectives for survival. Here, we present FlowGAT, a hybrid FBA-machine learning strategy for predicting essentiality directly from wild type metabolic phenotypes. The approach is based on graph-structured representation of metabolic fluxes predicted by FBA, where nodes correspond to enzymatic reactions and edges quantify the propagation of metabolite mass flow between a reaction and its neighbours. We integrate this information into a graph neural network that can be trained on knock-out fitness assay data. Comparisons across different model architectures reveal that FlowGAT predictions for E. coli are close to those of FBA for several growth conditions. This suggests that essentiality of enzymatic genes can be predicted by exploiting the inherent network structure of metabolism. Our approach demonstrates the benefits of combining the mechanistic insights afforded by genome-scale models with the ability of deep learning to infer patterns from complex datasets.

COSMIC-dFBA: A novel multi-scale hybrid framework for bioprocess modeling

Metabolism governs cell performance in biomanufacturing, as it fuels growth and productivity. However, even in well-controlled culture systems, metabolism is dynamic, with shifting objectives and resources, thus limiting the predictive capability of mechanistic models for process design and optimization. Here, we present Cellular Objectives and State Modulation In bioreaCtors (COSMIC)-dFBA, a hybrid multi-scale modeling paradigm that accurately predicts cell density, antibody titer, and bioreactor metabolite concentration profiles. Using machine-learning, COSMIC-dFBA decomposes the instantaneous metabolite uptake and secretion rates in a bioreactor into weighted contributions from each cell state (growth or antibody-producing state) and integrates these with a genome-scale metabolic model. A major strength of COSMIC-dFBA is that it can be parameterized with only metabolite concentrations from spent media, although constraining the metabolic model with other omics data can further improve its capabilities. Using COSMIC-dFBA, we can predict the final cell density and antibody titer to within 10% of the measured data, and compared to a standard dFBA model, we found the framework showed a 90% and 72% improvement in cell density and antibody titer prediction, respectively. Thus, we demonstrate our hybrid modeling framework effectively captures cellular metabolism and expands the applicability of dFBA to model the dynamic conditions in a bioreactor.

Genome-scale metabolic modelling of extremophiles and its applications in astrobiological environments

Metabolic modelling approaches have become the powerful tools in modern biology. These mathematical models are widely used to predict metabolic phenotypes of the organisms or communities of interest, and to identify metabolic targets in metabolic engineering. Apart from a broad range of industrial applications, the possibility of using metabolic modelling in the contexts of astrobiology are poorly explored. In this mini-review, we consolidated the concepts and related applications of applying metabolic modelling in studying organisms in space-related environments, specifically the extremophilic microbes. We recapitulated the current state of the art in metabolic modelling approaches and their advantages in the astrobiological context. Our review encompassed the applications of metabolic modelling in the theoretical investigation of the origin of life within prebiotic environments, as well as the compilation of existing uses of genome-scale metabolic models of extremophiles. Furthermore, we emphasize the current challenges associated with applying this technique in extreme environments, and conclude this review by discussing the potential implementation of metabolic models to explore theoretically optimal metabolic networks under various space conditions. Through this mini-review, our aim is to highlight the potential of metabolic modelling in advancing the study of astrobiology.

Flux sampling in genome-scale metabolic modeling of microbial communities

BACKGROUND

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates.

RESULTS

In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux.

CONCLUSIONS

Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Metabolic regulation boosts bioelectricity generation in Zymomonas mobilis microbial fuel cell, surpassing ethanol production

Zymomonas mobilis (Z. mobilis), a bacterium known for its ethanol production capabilities, can also generate electricity by transitioning from ethanol production to electron generation. The purpose of this study is to investigate the ability of Z. mobilis to produce bioelectricity when utilized as a biocatalyst in a single-chamber microbial fuel cell (MFC). Given the bacterium's strong inclination towards ethanol production, a metabolic engineering strategy was devised to identify key reactions responsible for redirecting electrons from ethanol towards electricity generation. To evaluate the electroactivity of cultured Z. mobilis and its ethanol production in the presence of regulators, the reduction of soluble Fe(III) was utilized. Among the regulators tested, CuCl2 demonstrated superior effectiveness. Consequently, the MFC was employed to analyze the electrochemical properties of Z. mobilis using both a minimal and modified medium. By modifying the bacterial medium, the maximum current and power density of the MFC fed with Z. mobilis increased by more than 5.8- and sixfold, respectively, compared to the minimal medium. These findings highlight the significant impact of metabolic redirection in enhancing the performance of MFCs. Furthermore, they establish Z. mobilis as an active electrogenesis microorganism capable of power generation in MFCs.

Cell-culture process optimization via model-based predictions of metabolism and protein glycosylation

In order to meet the rising demand for biologics and become competitive on the developing biosimilar market, there is a need for process intensification of biomanufacturing processes. Process development of biologics has historically relied on extensive experimentation to develop and optimize biopharmaceutical manufacturing. Experimentation to optimize media formulations, feeding schedules, bioreactor operations and bioreactor scale up is expensive, labor intensive and time consuming. Mathematical modeling frameworks have the potential to enable process intensification while reducing the experimental burden. This review focuses on mathematical modeling of cellular metabolism and N-linked glycosylation as applied to upstream manufacturing of biologics. We review developments in the field of modeling cellular metabolism of mammalian cells using kinetic and stoichiometric modeling frameworks along with their applications to simulate, optimize and improve mechanistic understanding of the process. Interest in modeling N-linked glycosylation has led to the creation of various types of parametric and non-parametric models. Most published studies on mammalian cell metabolism have performed experiments in shake flasks where the pH and dissolved oxygen cannot be controlled. Efforts to understand and model the effect of bioreactor-specific parameters such as pH, dissolved oxygen, temperature, and bioreactor heterogeneity are critically reviewed. Most modeling efforts have focused on the Chinese Hamster Ovary (CHO) cells, which are most commonly used to produce monoclonal antibodies (mAbs). However, these modeling approaches can be generalized and applied to any mammalian cell-based manufacturing platform. Current and potential future applications of these models for Vero cell-based vaccine manufacturing, CAR-T cell therapies, and viral vector manufacturing are also discussed. We offer specific recommendations for improving the applicability of these models to industrially relevant processes.

Genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer

Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment (TME), which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the TME. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

Flux balance analysis-based metabolic modeling of microbial secondary metabolism: Current status and outlook

In microorganisms, different from primary metabolism for cellular growth, secondary metabolism is for ecological interactions and stress responses and an important source of natural products widely used in various areas such as pharmaceutics and food additives. With advancements of sequencing technologies and bioinformatics tools, a large number of biosynthetic gene clusters of secondary metabolites have been discovered from microbial genomes. However, due to challenges from the difficulty of genome-scale pathway reconstruction and the limitation of conventional flux balance analysis (FBA) on secondary metabolism, the quantitative modeling of secondary metabolism is poorly established, in contrast to that of primary metabolism. This review first discusses current efforts on the reconstruction of secondary metabolic pathways in genome-scale metabolic models (GSMMs), as well as related FBA-based modeling techniques. Additionally, potential extensions of FBA are suggested to improve the prediction accuracy of secondary metabolite production. As this review posits, biosynthetic pathway reconstruction for various secondary metabolites will become automated and a modeling framework capturing secondary metabolism onset will enhance the predictive power. Expectedly, an improved FBA-based modeling workflow will facilitate quantitative study of secondary metabolism and in silico design of engineering strategies for natural product production.

Flux Sampling in Genome-scale Metabolic Modeling of Microbial Communities

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model. However, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling may capture additional heterogeneity across cells, especially when cells exhibit sub-maximal growth rates. In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. We find significant differences in the predicted metabolism with sampling, including increased cooperative interactions and pathway-specific changes in predicted flux. Our results suggest the importance of sampling-based and objective function-independent approaches to evaluate metabolic interactions and emphasize their utility in quantitatively studying interactions between cells and organisms.

Ensemble-based genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer

1Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment, which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the tumor microenvironment. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

Context-Specific Genome-Scale Metabolic Modelling and Its Application to the Analysis of COVID-19 Metabolic Signatures

Genome-scale metabolic models (GEMs) have found numerous applications in different domains, ranging from biotechnology to systems medicine. Herein, we overview the most popular algorithms for the automated reconstruction of context-specific GEMs using high-throughput experimental data. Moreover, we describe different datasets applied in the process, and protocols that can be used to further automate the model reconstruction and validation. Finally, we describe recent COVID-19 applications of context-specific GEMs, focusing on the analysis of metabolic implications, identification of biomarkers and potential drug targets.

The choice of the objective function in flux balance analysis is crucial for predicting replicative lifespans in yeast

Flux balance analysis (FBA) is a powerful tool to study genome-scale models of the cellular metabolism, based on finding the optimal flux distributions over the network. While the objective function is crucial for the outcome, its choice, even though motivated by evolutionary arguments, has not been directly connected to related measures. Here, we used an available multi-scale mathematical model of yeast replicative ageing, integrating cellular metabolism, nutrient sensing and damage accumulation, to systematically test the effect of commonly used objective functions on features of replicative ageing in budding yeast, such as the number of cell divisions and the corresponding time between divisions. The simulations confirmed that assuming maximal growth is essential for reaching realistic lifespans. The usage of the parsimonious solution or the additional maximisation of a growth-independent energy cost can improve lifespan predictions, explained by either increased respiratory activity using resources otherwise allocated to cellular growth or by enhancing antioxidative activity, specifically in early life. Our work provides a new perspective on choosing the objective function in FBA by connecting it to replicative ageing.

Role of Dissimilative Pathway of Komagataella phaffii (Pichia pastoris): Formaldehyde Toxicity and Energy Metabolism

Komagataella phaffii (aka Pichia pastoris) is a yeast able to grow in methanol as the sole carbon and energy source. This substrate is converted into formaldehyde, a toxic intermediary that can either be assimilated to biomass or dissimilated to CO₂ through the enzymes formaldehyde dehydrogenase (FLD) and formate dehydrogenase, also producing energy in the form of NADH. The dissimilative pathway has been described as an energy producing and a detoxifying route, but conclusive evidence has not been provided for this. In order to elucidate this theory, we generated mutants lacking the FLD activity (Δfld1) and used flux analysis to evaluate the metabolic impact of this disrupted pathway. Unexpectedly, we found that the specific growth rate of the Δfld1 strain was only slightly lower (92%) than the control. In contrast, the sensitivity to formaldehyde pulses (up to 8mM) was significantly higher in the Δfld1 mutant strain and was associated with a higher maintenance energy. In addition, the intracellular flux estimation revealed a high metabolic flexibility of K. phaffii in response to the disrupted pathway. Our results suggest that the role of the dissimilative pathway is mainly to protect the cells from the harmful effect of formaldehyde, as they were able to compensate for the energy provided from this pathway when disrupted.

Integration of omics data to generate and analyse COVID-19 specific genome-scale metabolic models

COVID-19 presents a complex disease that needs to be addressed using systems medicine approaches that include genome-scale metabolic models (GEMs). Previous studies have used a single model extraction method (MEM) and/or a single transcriptomic dataset to reconstruct context-specific models, which proved to be insufficient for the broader biological contexts. We have applied four MEMs in combination with five COVID-19 datasets. Models produced by GIMME were separated by infection, while tINIT preserved the biological variability in the data and enabled the best prediction of the enrichment of metabolic subsystems. Vitamin D3 metabolism was predicted to be down-regulated in one dataset by GIMME, and in all by tINIT. Models generated by tINIT and GIMME predicted downregulation of retinol metabolism in different datasets, while downregulated cholesterol metabolism was predicted only by tINIT-generated models. Predictions are in line with the observations in COVID-19 patients. Our data indicated that GIMME and tINIT models provided the most biologically relevant results and should have a larger emphasis in further analyses. Particularly tINIT models identified the metabolic pathways that are a part of the host response and are potential antiviral targets. The code and the results of the analyses are available to download from https://github.com/CompBioLj/COVID_GEMs_and_MEMs.

Intelligent host engineering for metabolic flux optimisation in biotechnology

Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.

Uncovering the Role of Metabolism in Oomycete-Host Interactions Using Genome-Scale Metabolic Models

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the "mass-flow" through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host-pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens.

Combining Kinetic and Constraint-Based Modelling to Better Understand Metabolism Dynamics

To understand the phenotypic capabilities of organisms, it is useful to characterise cellular metabolism through the analysis of its pathways. Dynamic mathematical modelling of metabolic networks is of high interest as it provides the time evolution of the metabolic components. However, it also has limitations, such as the necessary mechanistic details and kinetic parameters are not always available. On the other hand, large metabolic networks exhibit a complex topological structure which can be studied rather efficiently in their stationary regime by constraint-based methods. These methods produce useful predictions on pathway operations. In this review, we present both modelling techniques and we show how they bring complementary views of metabolism. In particular, we show on a simple example how both approaches can be used in conjunction to shed some light on the dynamics of metabolic networks.

Systematic evaluation of parameters for genome-scale metabolic models of cultured mammalian cells

Genome-scale metabolic models describe cellular metabolism with mechanistic detail. Given their high complexity, such models need to be parameterized correctly to yield accurate predictions and avoid overfitting. Effective parameterization has been well-studied for microbial models, but it remains unclear for higher eukaryotes, including mammalian cells. To address this, we enumerated model parameters that describe key features of cultured mammalian cells - including cellular composition, bioprocess performance metrics, mammalian-specific pathways, and biological assumptions behind model formulation approaches. We tested these parameters by building thousands of metabolic models and evaluating their ability to predict the growth rates of a panel of phenotypically diverse Chinese Hamster Ovary cell clones. We found the following considerations to be most critical for accurate parameterization: (1) cells limit metabolic activity to maintain homeostasis, (2) cell morphology and viability change dynamically during a growth curve, and (3) cellular biomass has a particular macromolecular composition. Depending on parameterization, models predicted different metabolic phenotypes, including contrasting mechanisms of nutrient utilization and energy generation, leading to varying accuracies of growth rate predictions. Notably, accurate parameter values broadly agreed with experimental measurements. These insights will guide future investigations of mammalian metabolism.

The Metano Modeling Toolbox MMTB: An Intuitive, Web-Based Toolbox Introduced by Two Use Cases

Genome-scale metabolic models are of high interest in a number of different research fields. Flux balance analysis (FBA) and other mathematical methods allow the prediction of the steady-state behavior of metabolic networks under different environmental conditions. However, many existing applications for flux optimizations do not provide a metabolite-centric view on fluxes. Metano is a standalone, open-source toolbox for the analysis and refinement of metabolic models. While flux distributions in metabolic networks are predominantly analyzed from a reaction-centric point of view, the Metano methods of split-ratio analysis and metabolite flux minimization also allow a metabolite-centric view on flux distributions. In addition, we present MMTB (Metano Modeling Toolbox), a web-based toolbox for metabolic modeling including a user-friendly interface to Metano methods. MMTB assists during bottom-up construction of metabolic models by integrating reaction and enzymatic annotation data from different databases. Furthermore, MMTB is especially designed for non-experienced users by providing an intuitive interface to the most commonly used modeling methods and offering novel visualizations. Additionally, MMTB allows users to upload their models, which can in turn be explored and analyzed by the community. We introduce MMTB by two use cases, involving a published model of Corynebacterium glutamicum and a newly created model of Phaeobacter inhibens.

Opportunities at the Interface of Network Science and Metabolic Modeling

Metabolism plays a central role in cell physiology because it provides the molecular machinery for growth. At the genome-scale, metabolism is made up of thousands of reactions interacting with one another. Untangling this complexity is key to understand how cells respond to genetic, environmental, or therapeutic perturbations. Here we discuss the roles of two complementary strategies for the analysis of genome-scale metabolic models: Flux Balance Analysis (FBA) and network science. While FBA estimates metabolic flux on the basis of an optimization principle, network approaches reveal emergent properties of the global metabolic connectivity. We highlight how the integration of both approaches promises to deliver insights on the structure and function of metabolic systems with wide-ranging implications in discovery science, precision medicine and industrial biotechnology.

Personalized Genome-Scale Metabolic Models Identify Targets of Redox Metabolism in Radiation-Resistant Tumors

Redox cofactor production is integral toward antioxidant generation, clearance of reactive oxygen species, and overall tumor response to ionizing radiation treatment. To identify systems-level alterations in redox metabolism that confer resistance to radiation therapy, we developed a bioinformatics pipeline for integrating multi-omics data into personalized genome-scale flux balance analysis models of 716 radiation-sensitive and 199 radiation-resistant tumors. These models collectively predicted that radiation-resistant tumors reroute metabolic flux to increase mitochondrial NADPH stores and reactive oxygen species (ROS) scavenging. Simulated genome-wide knockout screens agreed with experimental siRNA gene knockdowns in matched radiation-sensitive and radiation-resistant cancer cell lines, revealing gene targets involved in mitochondrial NADPH production, central carbon metabolism, and folate metabolism that allow for selective inhibition of glutathione production and H₂O₂ clearance in radiation-resistant cancers. This systems approach represents a significant advancement in developing quantitative genome-scale models of redox metabolism and identifying personalized metabolic targets for improving radiation sensitivity in individual cancer patients.

QTL mapping of modelled metabolic fluxes reveals gene variants impacting yeast central carbon metabolism

The yeast Saccharomyces cerevisiae is an attractive industrial microorganism for the production of foods and beverages as well as for various bulk and fine chemicals, such as biofuels or fragrances. Building blocks for these biosyntheses are intermediates of yeast central carbon metabolism (CCM), whose intracellular availability depends on balanced single reactions that form metabolic fluxes. Therefore, efficient product biosynthesis is influenced by the distribution of these fluxes. We recently demonstrated great variations in CCM fluxes between yeast strains of different origins. However, we have limited understanding of flux modulation and the genetic basis of flux variations. In this study, we investigated the potential of quantitative trait locus (QTL) mapping to elucidate genetic variations responsible for differences in metabolic flux distributions (fQTL). Intracellular metabolic fluxes were estimated by constraint-based modelling and used as quantitative phenotypes, and differences in fluxes were linked to genomic variations. Using this approach, we detected four fQTLs that influence metabolic pathways. The molecular dissection of these QTLs revealed two allelic gene variants, PDB1 and VID30, contributing to flux distribution. The elucidation of genetic determinants influencing metabolic fluxes, as reported here for the first time, creates new opportunities for the development of strains with optimized metabolite profiles for various applications.

Flux sampling is a powerful tool to study metabolism under changing environmental conditions

The development of high-throughput 'omic techniques has sparked a rising interest in genome-scale metabolic models, with applications ranging from disease diagnostics to crop adaptation. Efficient and accurate methods are required to analyze large metabolic networks. Flux sampling can be used to explore the feasible flux solutions in metabolic networks by generating probability distributions of steady-state reaction fluxes. Unlike other methods, flux sampling can be used without assuming a particular cellular objective. We have undertaken a rigorous comparison of several sampling algorithms and concluded that the coordinate hit-and-run with rounding (CHRR) algorithm is the most efficient based on both run-time and multiple convergence diagnostics. We demonstrate the power of CHRR by using it to study the metabolic changes that underlie photosynthetic acclimation to cold of Arabidopsis thaliana plant leaves. In combination with experimental measurements, we show how the regulated interplay between diurnal starch and organic acid accumulation defines the plant acclimation process. We confirm fumarate accumulation as a requirement for cold acclimation and further predict γ-aminobutyric acid to have a key role in metabolic signaling under cold conditions. These results demonstrate how flux sampling can be used to analyze the feasible flux solutions across changing environmental conditions, whereas eliminating the need to make assumptions which introduce observer bias.

Metabolic modeling of bacterial co-culture systems predicts enhanced carbon monoxide-to-butyrate conversion compared to monoculture systems

We used metabolic modeling to computationally investigate the potential of bacterial coculture system designs for CO conversion to the platform chemical butyrate. By taking advantage of the native capabilities of wild-type strains, we developed two anaerobic coculture designs by combining Clostridium autoethanogenum for CO-to-acetate conversion with bacterial strains that offer high acetate-to-butyrate conversion capabilities: the environmental bacterium the human gut bacteriumEubacterium rectale. When grown in continuous stirred tank reactor on a 70/0/30 CO/H₂/N₂ gas mixture, the C. autoethanogenum-C Kluyveri co-culture was predicted to offer no mprovement in butyrate volumetric productivity compared to an engineered C. autoethanogenum monoculture despite utilizing vinyl acetate as a secondary carbon source for C. kluyveri growth enhancement. A coculture consisting of C. autoethanogenum and C. kluyveri engineered in silico to eliminate hexanoate synthesis was predicted to enhance both butyrate productivity and titer. The C. autoethanogenum-E. rectale coculture offered similar improvements in butyrate productivity without the need for metabolic engineering when glucose was provided as a secondary carbon source to enhance E. rectale growth. A bubble column model developed to assess the potential for large-scale butyrate production of the C. autoethanogenum-E. rectale design predicted that a 40/30/30 CO/H₂/N₂ gas mixture and a 5 m column length would be preferred to enhance C. autoethanogenum growth and counteract CO inhibitory effects on E. rectale.

Targeting NAD+ Metabolism to Enhance Radiation Therapy Responses

Nicotinamide adenine dinucleotide (NAD⁺) metabolism is integrally connected with the mechanisms of action of radiation therapy and is altered in many radiation-resistant tumors. This makes NAD⁺ metabolism an ideal target for therapies that increase radiation sensitivity and improve patient outcomes. This review provides an overview of NAD⁺ metabolism in the context of the cellular response to ionizing radiation, as well as current therapies that target NAD⁺ metabolism to enhance radiation therapy responses. Additionally, we summarize state-of-the-art methods for measuring, modeling, and manipulating NAD⁺ metabolism, which are being used to identify novel targets in the NAD⁺ metabolic network for therapeutic interventions in combination with radiation therapy.

How adaptive evolution reshapes metabolism to improve fitness: recent advances and future outlook

Adaptive laboratory evolution (ALE) has emerged as a powerful tool in basic microbial research and strain development. In the context of metabolic science and engineering, it has been applied to study gene knockout responses, expand substrate ranges, improve tolerance to process conditions, and to improve productivity via designed growth coupling. In recent years, advancements in ALE methods and systems biology measurement technologies, particularly genome sequencing and ¹³C metabolic flux analysis (¹³C-MFA), have enabled detailed study of the mechanisms and dynamics of evolving metabolism. In this review, we discuss a range of studies that have applied flux analysis to adaptively evolved strains, as well as modeling frameworks developed to predict and interpret evolved fluxes. These efforts link mutations to fitness-enhanced phenotypes, identify bottlenecks and approaches to resolve them, and address systems concepts such as optimality.

Advances in metabolic modeling of oleaginous microalgae

Production of biofuels and bioenergy precursors by phototrophic microorganisms, such as microalgae and cyanobacteria, is a promising alternative to conventional fuels obtained from non-renewable resources. Several species of microalgae have been investigated as potential candidates for the production of biofuels, for the most part due to their exceptional metabolic capability to accumulate large quantities of lipids. Constraint-based modeling, a systems biology approach that accurately predicts the metabolic phenotype of phototrophs, has been deployed to identify suitable culture conditions as well as to explore genetic enhancement strategies for bioproduction. Core metabolic models were employed to gain insight into the central carbon metabolism in photosynthetic microorganisms. More recently, comprehensive genome-scale models, including organelle-specific information at high resolution, have been developed to gain new insight into the metabolism of phototrophic cell factories. Here, we review the current state of the art of constraint-based modeling and computational method development and discuss how advanced models led to increased prediction accuracy and thus improved lipid production in microalgae.

Genome-wide characterization of Phytophthora infestans metabolism: a systems biology approach

Genome-scale metabolic models (GEMs) provide a functional view of the complex network of biochemical reactions in the living cell. Initially mainly applied to reconstruct the metabolism of model organisms, the availability of increasingly sophisticated reconstruction methods and more extensive biochemical databases now make it possible to reconstruct GEMs for less well-characterized organisms, and have the potential to unravel the metabolism in pathogen-host systems. Here, we present a GEM for the oomycete plant pathogen Phytophthora infestans as a first step towards an integrative model with its host. We predict the biochemical reactions in different cellular compartments and investigate the gene-protein-reaction associations in this model to obtain an impression of the biochemical capabilities of P. infestans. Furthermore, we generate life stage-specific models to place the transcriptomic changes of the genes encoding metabolic enzymes into a functional context. In sporangia and zoospores, there is an overall down-regulation, most strikingly reflected in the fatty acid biosynthesis pathway. To investigate the robustness of the GEM, we simulate gene deletions to predict which enzymes are essential for in vitro growth. This model is an essential first step towards an understanding of P. infestans and its interactions with plants as a system, which will help to formulate new hypotheses on infection mechanisms and disease prevention.

Metabolic adaptation of two in silico mutants of Mycobacterium tuberculosis during infection

BACKGROUND

Up to date, Mycobacterium tuberculosis (Mtb) remains as the worst intracellular killer pathogen. To establish infection, inside the granuloma, Mtb reprograms its metabolism to support both growth and survival, keeping a balance between catabolism, anabolism and energy supply. Mtb knockouts with the faculty of being essential on a wide range of nutritional conditions are deemed as target candidates for tuberculosis (TB) treatment. Constraint-based genome-scale modeling is considered as a promising tool for evaluating genetic and nutritional perturbations on Mtb metabolic reprogramming. Nonetheless, few in silico assessments of the effect of nutritional conditions on Mtb's vulnerability and metabolic adaptation have been carried out.

RESULTS

A genome-scale model (GEM) of Mtb, modified from the H37Rv iOSDD890, was used to explore the metabolic reprogramming of two Mtb knockout mutants (pfkA- and icl-mutants), lacking key enzymes of central carbon metabolism, while exposed to changing nutritional conditions (oxygen, and carbon and nitrogen sources). A combination of shadow pricing, sensitivity analysis, and flux distributions patterns allowed us to identify metabolic behaviors that are in agreement with phenotypes reported in the literature. During hypoxia, at high glucose consumption, the Mtb pfkA-mutant showed a detrimental growth effect derived from the accumulation of toxic sugar phosphate intermediates (glucose-6-phosphate and fructose-6-phosphate) along with an increment of carbon fluxes towards the reductive direction of the tricarboxylic acid cycle (TCA). Furthermore, metabolic reprogramming of the icl-mutant (icl1&icl2) showed the importance of the methylmalonyl pathway for the detoxification of propionyl-CoA, during growth at high fatty acid consumption rates and aerobic conditions. At elevated levels of fatty acid uptake and hypoxia, we found a drop in TCA cycle intermediate accumulation that might create redox imbalance. Finally, findings regarding Mtb-mutant metabolic adaptation associated with asparagine consumption and acetate, succinate and alanine production, were in agreement with literature reports.

CONCLUSIONS

This study demonstrates the potential application of genome-scale modeling, flux balance analysis (FBA), phenotypic phase plane (PhPP) analysis and shadow pricing to generate valuable insights about Mtb metabolic reprogramming in the context of human granulomas.

Fast growth phenotype of E. coli K-12 from adaptive laboratory evolution does not require intracellular flux rewiring

Adaptive laboratory evolution (ALE) is a widely-used method for improving the fitness of microorganisms in selected environmental conditions. It has been applied previously to Escherichia coli K-12 MG1655 during aerobic exponential growth on glucose minimal media, a frequently used model organism and growth condition, to probe the limits of E. coli growth rate and gain insights into fast growth phenotypes. Previous studies have described up to 1.6-fold increases in growth rate following ALE, and have identified key causal genetic mutations and changes in transcriptional patterns. Here, we report for the first time intracellular metabolic fluxes for six such adaptively evolved strains, as determined by high-resolution ¹³C-metabolic flux analysis. Interestingly, we found that intracellular metabolic pathway usage changed very little following adaptive evolution. Instead, at the level of central carbon metabolism the faster growth was facilitated by proportional increases in glucose uptake and all intracellular rates. Of the six evolved strains studied here, only one strain showed a small degree of flux rewiring, and this was also the strain with unique genetic mutations. A comparison of fluxes with two other wild-type (unevolved) E. coli strains, BW25113 and BL21, showed that inter-strain differences are greater than differences between the parental and evolved strains. Principal component analysis highlighted that nearly all flux differences (95%) between the nine strains were captured by only two principal components. The distance between measured and flux balance analysis predicted fluxes was also investigated. It suggested a relatively wide range of similar stoichiometric optima, which opens new questions about the path-dependency of adaptive evolution.

Clostridium butyricum maximizes growth while minimizing enzyme usage and ATP production: metabolic flux distribution of a strain cultured in glycerol

BACKGROUND

The increase in glycerol obtained as a byproduct of biodiesel has encouraged the production of new industrial products, such as 1,3-propanediol (PDO), using biotechnological transformation via bacteria like Clostridium butyricum. However, despite the increasing role of Clostridium butyricum as a bio-production platform, its metabolism remains poorly modeled.

RESULTS

We reconstructed iCbu641, the first genome-scale metabolic (GSM) model of a PDO producer Clostridium strain, which included 641 genes, 365 enzymes, 891 reactions, and 701 metabolites. We found an enzyme expression prediction of nearly 84% after comparison of proteomic data with flux distribution estimation using flux balance analysis (FBA). The remaining 16% corresponded to enzymes directionally coupled to growth, according to flux coupling findings (FCF). The fermentation data validation also revealed different phenotype states that depended on culture media conditions; for example, Clostridium maximizes its biomass yield per enzyme usage under glycerol limitation. By contrast, under glycerol excess conditions, Clostridium grows sub-optimally, maximizing biomass yield while minimizing both enzyme usage and ATP production. We further evaluated perturbations in the GSM model through enzyme deletions and variations in biomass composition. The GSM predictions showed no significant increase in PDO production, suggesting a robustness to perturbations in the GSM model. We used the experimental results to predict that co-fermentation was a better alternative than iCbu641 perturbations for improving PDO yields.

CONCLUSIONS

The agreement between the predicted and experimental values allows the use of the GSM model constructed for the PDO-producing Clostridium strain to propose new scenarios for PDO production, such as dynamic simulations, thereby reducing the time and costs associated with experimentation.

Making life difficult for Clostridium difficile: augmenting the pathogen's metabolic model with transcriptomic and codon usage data for better therapeutic target characterization

BACKGROUND

Clostridium difficile is a bacterium which can infect various animal species, including humans. Infection with this bacterium is a leading healthcare-associated illness. A better understanding of this organism and the relationship between its genotype and phenotype is essential to the search for an effective treatment. Genome-scale metabolic models contain all known biochemical reactions of a microorganism and can be used to investigate this relationship.

RESULTS

We present icdf834, an updated metabolic network of C. difficile that builds on iMLTC806cdf and features 1227 reactions, 834 genes, and 807 metabolites. We used this metabolic network to reconstruct the metabolic landscape of this bacterium. The standard metabolic model cannot account for changes in the bacterial metabolism in response to different environmental conditions. To account for this limitation, we also integrated transcriptomic data, which details the gene expression of the bacterium in a wide array of environments. Importantly, to bridge the gap between gene expression levels and protein abundance, we accounted for the synonymous codon usage bias of the bacterium in the model. To our knowledge, this is the first time codon usage has been quantified and integrated into a metabolic model. The metabolic fluxes were defined as a function of protein abundance. To determine potential therapeutic targets using the model, we conducted gene essentiality and metabolic pathway sensitivity analyses and calculated flux control coefficients. We obtained 92.3% accuracy in predicting gene essentiality when compared to experimental data for C. difficile R20291 (ribotype 027) homologs. We validated our context-specific metabolic models using sensitivity and robustness analyses and compared model predictions with literature on C. difficile. The model predicts interesting facets of the bacterium's metabolism, such as changes in the bacterium's growth in response to different environmental conditions.

CONCLUSIONS

After an extensive validation process, we used icdf834 to obtain state-of-the-art predictions of therapeutic targets for C. difficile. We show how context-specific metabolic models augmented with codon usage information can be a beneficial resource for better understanding C. difficile and for identifying novel therapeutic targets. We remark that our approach can be applied to investigate and treat against other pathogens.

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures-Systems biology-based interpretation using genome-scale metabolic flux balance model and multivariate data analysis

Genome-scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745-753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome-scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163-1173, 2016.

Diversity of flux distribution in central carbon metabolism of S. cerevisiae strains from diverse environments

Background

S. cerevisiae has attracted considerable interest in recent years as a model for ecology and evolutionary biology, revealing a substantial genetic and phenotypic diversity. However, there is a lack of knowledge on the diversity of metabolic networks within this species.

Results

To identify the metabolic and evolutionary constraints that shape metabolic fluxes in S. cerevisiae, we used a dedicated constraint-based model to predict the central carbon metabolism flux distribution of 43 strains from different ecological origins, grown in wine fermentation conditions. In analyzing these distributions, we observed a highly contrasted situation in flux variability, with quasi-constancy of the glycolysis and ethanol synthesis yield yet high flexibility of other fluxes, such as the pentose phosphate pathway and acetaldehyde production. Furthermore, these fluxes with large variability showed multimodal distributions that could be linked to strain origin, indicating a convergence between genetic origin and flux phenotype.

Conclusions

Flux variability is pathway-dependent and, for some flux, a strain origin effect can be found. These data highlight the constraints shaping the yeast operative central carbon network and provide clues for the design of strategies for strain improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0456-0) contains supplementary material, which is available to authorized users.

Cell scale host-pathogen modeling: another branch in the evolution of constraint-based methods

Constraint-based models have become popular methods for systems biology as they enable the integration of complex, disparate datasets in a biologically cohesive framework that also supports the description of biological processes in terms of basic physicochemical constraints and relationships. The scope, scale, and application of genome scale models have grown from single cell bacteria to multi-cellular interaction modeling; host-pathogen modeling represents one of these examples at the current horizon of constraint-based methods. There are now a small number of examples of host-pathogen constraint-based models in the literature, however there has not yet been a definitive description of the methodology required for the functional integration of genome scale models in order to generate simulation capable host-pathogen models. Herein we outline a systematic procedure to produce functional host-pathogen models, highlighting steps which require debugging and iterative revisions in order to successfully build a functional model. The construction of such models will enable the exploration of host-pathogen interactions by leveraging the growing wealth of omic data in order to better understand mechanism of infection and identify novel therapeutic strategies.

Flux balance analysis in the production of clavulanic acid by Streptomyces clavuligerus

In this work, in silico flux balance analysis is used for predicting the metabolic behavior of Streptomyces clavuligerus during clavulanic acid production. To choose the best objective function for use in the analysis, three different optimization problems are evaluated inside the flux balance analysis formulation: (i) maximization of the specific growth rate, (ii) maximization of the ATP yield, and (iii) maximization of clavulanic acid production. Maximization of ATP yield showed the best predictions for the cellular behavior. Therefore, flux balance analysis using ATP as objective function was used for analyzing different scenarios of nutrient limitations toward establishing the effect of limiting the carbon, nitrogen, phosphorous, and oxygen sources on the growth and clavulanic acid production rates. Obtained results showed that ammonia and phosphate limitations are the ones most strongly affecting clavulanic acid biosynthesis. Furthermore, it was possible to identify the ornithine flux from the urea cycle and the α-ketoglutarate flux from the TCA cycle as the most determinant internal fluxes for promoting clavulanic acid production.

Effect of plasmid replication deregulation via inc mutations on E. coli proteome & simple flux model analysis

When the replication of a plasmid based on sucrose selection is deregulated via the inc1 and inc2 mutations, high copy numbers (7,000 or greater) are attained while the growth rate on minimal medium is negligibly affected. Adaptions were assumed to be required in order to sustain the growth rate. Proteomics indicated that indeed a number of adaptations occurred that included increased expression of ribosomal proteins and 2-oxoglutarate dehydrogenase. The operating space prescribed by a basic flux model that maintained phenotypic traits (e.g. growth, byproducts, etc.) within typical bounds of resolution was consistent with the flux implications of the proteomic changes.

MetDFBA: incorporating time-resolved metabolomics measurements into dynamic flux balance analysis

This paper presents MetDFBA, a new approach incorporating experimental metabolomics time-series into constraint-based modeling. The method can be used for hypothesis testing and predicting dynamic flux profiles.