Computational approaches in metabolic engineering

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

FastKnock: An efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using three Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more useful and important practical solutions. The availability of all the solutions provides the opportunity to further characterize and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock.

Improved 13C metabolic flux analysis in Escherichia coli metabolism: application of a high-resolution MS (GC-EI-QTOF) for comprehensive assessment of MS/MS fragments

Gas chromatography-tandem mass spectrometry with electron ionization (GC-EI-MS/MS) provides rich information on stable-isotope labeling for 13C-metabolic flux analysis (13C-MFA). To pave the way for the routine application of tandem MS data for metabolic flux quantification, we aimed to compile a comprehensive library of GC-EI-MS/MS fragments of tert-butyldimethylsilyl (TBDMS) derivatized proteinogenic amino acids. First, we established an analytical workflow that combines high-resolution gas chromatography-quadrupole time-of-flight mass spectrometry and fully 13C-labeled biomass to identify and structurally elucidate tandem MS amino acid fragments. Application of the high-mass accuracy MS procedure resulted into the identification of 129 validated precursor-product ion pairs of 13 amino acids with 30 fragments being accepted for 13C-MFA. The practical benefit of the novel tandem MS data was demonstrated by a proof-of-concept study, which confirmed the importance of the compiled library for high-resolution 13C-MFA.

ONE SENTENCE SUMMARY

An analytical workflow that combines high-resolution mass spectrometry (MS) and fully 13C-labeled biomass to identify and structurally elucidate tandem MS amino acid fragments, which provide positional information and therefore offering significant advantages over traditional MS to improve 13C-metabolic flux analysis.

The metabolic origins of non-photorespiratory CO2 release during photosynthesis: a metabolic flux analysis

Respiration in the light (RL) releases CO2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since RL lowers net carbon fixation, understanding RL could help improve plant carbon-use efficiency and models of crop photosynthesis. Although RL was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to RL, we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13CO2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N2; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that RL releases 5.2 μmol CO2 g-1 FW h-1 of CO2, which is relatively consistent with a value of 9.3 μmol CO2 g-1 FW h-1 measured by CO2 gas exchange. The results indicated that ≤10% of RL results from TCA cycle reactions, which are widely considered to dominate RL. Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >93% of CO2 released by RL.

Metabolic Reprogramming of Fibroblasts as Therapeutic Target in Rheumatoid Arthritis and Cancer: Deciphering Key Mechanisms Using Computational Systems Biology Approaches

Fibroblasts, the most abundant cells in the connective tissue, are key modulators of the extracellular matrix (ECM) composition. These spindle-shaped cells are capable of synthesizing various extracellular matrix proteins and collagen. They also provide the structural framework (stroma) for tissues and play a pivotal role in the wound healing process. While they are maintainers of the ECM turnover and regulate several physiological processes, they can also undergo transformations responding to certain stimuli and display aggressive phenotypes that contribute to disease pathophysiology. In this review, we focus on the metabolic pathways of glucose and highlight metabolic reprogramming as a critical event that contributes to the transition of fibroblasts from quiescent to activated and aggressive cells. We also cover the emerging evidence that allows us to draw parallels between fibroblasts in autoimmune disorders and more specifically in rheumatoid arthritis and cancer. We link the metabolic changes of fibroblasts to the toxic environment created by the disease condition and discuss how targeting of metabolic reprogramming could be employed in the treatment of such diseases. Lastly, we discuss Systems Biology approaches, and more specifically, computational modeling, as a means to elucidate pathogenetic mechanisms and accelerate the identification of novel therapeutic targets.

A Method to Constrain Genome-Scale Models with 13C Labeling Data

Current limitations in quantitatively predicting biological behavior hinder our efforts to engineer biological systems to produce biofuels and other desired chemicals. Here, we present a new method for calculating metabolic fluxes, key targets in metabolic engineering, that incorporates data from ¹³C labeling experiments and genome-scale models. The data from ¹³C labeling experiments provide strong flux constraints that eliminate the need to assume an evolutionary optimization principle such as the growth rate optimization assumption used in Flux Balance Analysis (FBA). This effective constraining is achieved by making the simple but biologically relevant assumption that flux flows from core to peripheral metabolism and does not flow back. The new method is significantly more robust than FBA with respect to errors in genome-scale model reconstruction. Furthermore, it can provide a comprehensive picture of metabolite balancing and predictions for unmeasured extracellular fluxes as constrained by ¹³C labeling data. A comparison shows that the results of this new method are similar to those found through ¹³C Metabolic Flux Analysis (¹³C MFA) for central carbon metabolism but, additionally, it provides flux estimates for peripheral metabolism. The extra validation gained by matching 48 relative labeling measurements is used to identify where and why several existing COnstraint Based Reconstruction and Analysis (COBRA) flux prediction algorithms fail. We demonstrate how to use this knowledge to refine these methods and improve their predictive capabilities. This method provides a reliable base upon which to improve the design of biological systems.

While metabolic fluxes constitute the most direct window into a cell’s metabolism, their accurate measurement is non trivial. The gold standard for flux measurement involves providing a labeled feed where some of the carbon atoms have been substituted by isotopes with higher atomic mass (¹³C instead of ¹²C). The ensuing labeling found in intracellular metabolites is then used to computationally infer the metabolic fluxes that produced the observed pattern. However, this procedure is typically performed with small metabolic models encompassing only central carbon metabolism. The genomic revolution has afforded us easily available genomes and, with them, comprehensive genome-scale models of cellular metabolism. It would be desirable to use the ¹³C labeling experimental data to constrain genome-scale models: these data constrain fluxes very effectively and provide in the labeling data fit an obvious proof that the underlying model correctly explains measured quantities. Here, we introduce a rigorous, self-consistent method that uses the full amount of information contained in ¹³C labeling data to constrain fluxes for a genome-scale model where underlying assumptions are explicitly stated.

Publishing 13C metabolic flux analysis studies: a review and future perspectives

(13)C-Metabolic flux analysis ((13)C-MFA) is a powerful model-based analysis technique for determining intracellular metabolic fluxes in living cells. It has become a standard tool in many labs for quantifying cell physiology, e.g., in metabolic engineering, systems biology, biotechnology, and biomedical research. With the increasing number of (13)C-MFA studies published each year, it is now ever more important to provide practical guidelines for performing and publishing (13)C-MFA studies so that quality is not sacrificed as the number of publications increases. The main purpose of this paper is to provide an overview of good practices in (13)C-MFA, which can eventually be used as minimum data standards for publishing (13)C-MFA studies. The motivation for this work is two-fold: (1) currently, there is no general consensus among researchers and journal editors as to what minimum data standards should be required for publishing (13)C-MFA studies; as a result, there are great discrepancies in terms of quality and consistency; and (2) there is a growing number of studies that cannot be reproduced or verified independently due to incomplete information provided in these publications. This creates confusion, e.g. when trying to reconcile conflicting results, and hinders progress in the field. Here, we review current status in the (13)C-MFA field and highlight some of the shortcomings with regards to (13)C-MFA publications. We then propose a checklist that encompasses good practices in (13)C-MFA. We hope that these guidelines will be a valuable resource for the community and allow (13)C-flux studies to be more easily reproduced and accessed by others in the future.

COMPLETE-MFA: complementary parallel labeling experiments technique for metabolic flux analysis

We have developed a novel approach for measuring highly accurate and precise metabolic fluxes in living cells, termed COMPLETE-MFA, short for complementary parallel labeling experiments technique for metabolic flux analysis. The COMPLETE-MFA method is based on combined analysis of multiple isotopic labeling experiments, where the synergy of using complementary tracers greatly improves the precision of estimated fluxes. In this work, we demonstrate the COMPLETE-MFA approach using all singly labeled glucose tracers, [1-(13)C], [2-(13)C], [3-(13)C], [4-(13)C], [5-(13)C], and [6-(13)C]glucose to determine precise metabolic fluxes for wild-type Escherichia coli. Cells were grown in six parallel cultures on defined medium with glucose as the only carbon source. Mass isotopomers of biomass amino acids were measured by gas chromatography-mass spectrometry (GC-MS). The data from all six experiments were then fitted simultaneously to a single flux model to determine accurate intracellular fluxes. We obtained a statistically acceptable fit with more than 300 redundant measurements. The estimated flux map is the most precise flux result obtained thus far for E. coli cells. To our knowledge, this is the first time that six isotopic labeling experiments have been successfully integrated for high-resolution (13)C-flux analysis.

Rational design of ¹³C-labeling experiments for metabolic flux analysis in mammalian cells

Background

¹³C-Metabolic flux analysis (¹³C-MFA) is a standard technique to probe cellular metabolism and elucidate in vivo metabolic fluxes. ¹³C-Tracer selection is an important step in conducting ¹³C-MFA, however, current methods are restricted to trial-and-error approaches, which commonly focus on an arbitrary subset of the tracer design space. To systematically probe the complete tracer design space, especially for complex systems such as mammalian cells, there is a pressing need for new rational approaches to identify optimal tracers.

Results

Recently, we introduced a new framework for optimal ¹³C-tracer design based on elementary metabolite units (EMU) decomposition, in which a measured metabolite is decomposed into a linear combination of so-called EMU basis vectors. In this contribution, we applied the EMU method to a realistic network model of mammalian metabolism with lactate as the measured metabolite. The method was used to select optimal tracers for two free fluxes in the system, the oxidative pentose phosphate pathway (oxPPP) flux and anaplerosis by pyruvate carboxylase (PC). Our approach was based on sensitivity analysis of EMU basis vector coefficients with respect to free fluxes. Through efficient grouping of coefficient sensitivities, simple tracer selection rules were derived for high-resolution quantification of the fluxes in the mammalian network model. The approach resulted in a significant reduction of the number of possible tracers and the feasible tracers were evaluated using numerical simulations. Two optimal, novel tracers were identified that have not been previously considered for ¹³C-MFA of mammalian cells, specifically [2,3,4,5,6-¹³C]glucose for elucidating oxPPP flux and [3,4-¹³C]glucose for elucidating PC flux. We demonstrate that ¹³C-glutamine tracers perform poorly in this system in comparison to the optimal glucose tracers.

Conclusions

In this work, we have demonstrated that optimal tracer design does not need to be a pure simulation-based trial-and-error process; rather, rational insights into tracer design can be gained through the application of the EMU basis vector methodology. Using this approach, rational labeling rules can be established a priori to guide the selection of optimal ¹³C-tracers for high-resolution flux elucidation in complex metabolic network models.

Selection of tracers for 13C-metabolic flux analysis using elementary metabolite units (EMU) basis vector methodology

Metabolic flux analysis (MFA) is a powerful technique for elucidating in vivo fluxes in microbial and mammalian systems. A key step in (13)C-MFA is the selection of an appropriate isotopic tracer to observe fluxes in a proposed network model. Despite the importance of MFA in metabolic engineering and beyond, current approaches for tracer experiment design are still largely based on trial-and-error. The lack of a rational methodology for selecting isotopic tracers prevents MFA from achieving its full potential. Here, we introduce a new technique for tracer experiment design based on the concept of elementary metabolite unit (EMU) basis vectors. We demonstrate that any metabolite in a network model can be expressed as a linear combination of so-called EMU basis vectors, where the corresponding coefficients indicate the fractional contribution of the EMU basis vector to the product metabolite. The strength of this approach is the decoupling of substrate labeling, i.e. the EMU basis vectors, from the dependence on free fluxes, i.e. the coefficients. In this work, we demonstrate that flux observability inherently depends on the number of independent EMU basis vectors and the sensitivities of coefficients with respect to free fluxes. Specifically, the number of independent EMU basis vectors places hard limits on how many free fluxes can be determined in a model. This constraint is used as a guide for selecting feasible substrate labeling. In three example models, we demonstrate that by maximizing the number of independent EMU basis vectors the observability of a system is improved. Inspection of sensitivities of coefficients with respect to free fluxes provides additional constraints for proper selection of tracers. The present contribution provides a fresh perspective on an important topic in metabolic engineering, and gives practical guidelines and design principles for a priori selection of isotopic tracers for (13)C-MFA studies.

Dynamic metabolic flux analysis (DMFA): a framework for determining fluxes at metabolic non-steady state

Metabolic flux analysis (MFA) is a key tool for measuring in vivo metabolic fluxes in systems at metabolic steady state. Here, we present a new method for dynamic metabolic flux analysis (DMFA) of systems that are not at metabolic steady state. The advantages of our DMFA method are: (1) time-series of metabolite concentration data can be applied directly for estimating dynamic fluxes, making data smoothing and estimation of average extracellular rates unnecessary; (2) flux estimation is achieved without integration of ODEs, or iterations; (3) characteristic metabolic phases in the fermentation data are identified automatically by the algorithm, rather than selected manually/arbitrarily. We demonstrate the application of the new DMFA framework in three example systems. First, we evaluated the performance of DMFA in a simple three-reaction model in terms of accuracy, precision and flux observability. Next, we analyzed a commercial glucose-limited fed-batch process for 1,3-propanediol production. The DMFA method accurately captured the dynamic behavior of the fed-batch fermentation and identified characteristic metabolic phases. Lastly, we demonstrate that DMFA can be used without any assumed metabolic network model for data reconciliation and detection of gross measurement errors using carbon and electron balances as constraints.

Large-scale bi-level strain design approaches and mixed-integer programming solution techniques

The use of computational models in metabolic engineering has been increasing as more genome-scale metabolic models and computational approaches become available. Various computational approaches have been developed to predict how genetic perturbations affect metabolic behavior at a systems level, and have been successfully used to engineer microbial strains with improved primary or secondary metabolite production. However, identification of metabolic engineering strategies involving a large number of perturbations is currently limited by computational resources due to the size of genome-scale models and the combinatorial nature of the problem. In this study, we present (i) two new bi-level strain design approaches using mixed-integer programming (MIP), and (ii) general solution techniques that improve the performance of MIP-based bi-level approaches. The first approach (SimOptStrain) simultaneously considers gene deletion and non-native reaction addition, while the second approach (BiMOMA) uses minimization of metabolic adjustment to predict knockout behavior in a MIP-based bi-level problem for the first time. Our general MIP solution techniques significantly reduced the CPU times needed to find optimal strategies when applied to an existing strain design approach (OptORF) (e.g., from ∼10 days to ∼5 minutes for metabolic engineering strategies with 4 gene deletions), and identified strategies for producing compounds where previous studies could not (e.g., malate and serine). Additionally, we found novel strategies using SimOptStrain with higher predicted production levels (for succinate and glycerol) than could have been found using an existing approach that considers network additions and deletions in sequential steps rather than simultaneously. Finally, using BiMOMA we found novel strategies involving large numbers of modifications (for pyruvate and glutamate), which sequential search and genetic algorithms were unable to find. The approaches and solution techniques developed here will facilitate the strain design process and extend the scope of its application to metabolic engineering.