Current status and applications of genome-scale metabolic models

Genome-scale metabolic models (GEMs) computationally describe gene-protein-reaction associations for entire metabolic genes in an organism, and can be simulated to predict metabolic fluxes for various systems-level metabolic studies. Since the first GEM for Haemophilus influenzae was reported in 1999, advances have been made to develop and simulate GEMs for an increasing number of organisms across bacteria, archaea, and eukarya. Here, we review current reconstructed GEMs and discuss their applications, including strain development for chemicals and materials production, drug targeting in pathogens, prediction of enzyme functions, pan-reactome analysis, modeling interactions among multiple cells or organisms, and understanding human diseases.

In silico identification of metabolic engineering strategies for improved lipid production in Yarrowia lipolytica by genome-scale metabolic modeling

BACKGROUND

Yarrowia lipolytica, an oleaginous yeast, is a promising platform strain for production of biofuels and oleochemicals as it can accumulate a high level of lipids in response to nitrogen limitation. Accordingly, many metabolic engineering efforts have been made to develop engineered strains of Y. lipolytica with higher lipid yields. Genome-scale model of metabolism (GEM) is a powerful tool for identifying novel genetic designs for metabolic engineering. Several GEMs for Y. lipolytica have recently been developed; however, not many applications of the GEMs have been reported for actual metabolic engineering of Y. lipolytica. The major obstacle impeding the application of Y. lipolytica GEMs is the lack of proper methods for predicting phenotypes of the cells in the nitrogen-limited condition, or more specifically in the stationary phase of a batch culture.

RESULTS

In this study, we showed that environmental version of minimization of metabolic adjustment (eMOMA) can be used for predicting metabolic flux distribution of Y. lipolytica under the nitrogen-limited condition and identifying metabolic engineering strategies to improve lipid production in Y. lipolytica. Several well-characterized overexpression targets, such as diglyceride acyltransferase, acetyl-CoA carboxylase, and stearoyl-CoA desaturase, were successfully rediscovered by our eMOMA-based design method, showing the relevance of prediction results. Interestingly, the eMOMA-based design method also suggested non-intuitive knockout targets, and we experimentally validated the prediction with a mutant lacking YALI0F30745g, one of the predicted targets involved in one-carbon/methionine metabolism. The mutant accumulated 45% more lipids compared to the wild-type.

CONCLUSION

This study demonstrated that eMOMA is a powerful computational method for understanding and engineering the metabolism of Y. lipolytica and potentially other oleaginous microorganisms.

In silico model-guided identification of transcriptional regulator targets for efficient strain design

BACKGROUND

Cellular metabolism is tightly regulated by hard-wired multiple layers of biological processes to achieve robust and homeostatic states given the limited resources. As a result, even the most intuitive enzyme-centric metabolic engineering endeavours through the up-/down-regulation of multiple genes in biochemical pathways often deliver insignificant improvements in the product yield. In this regard, targeted engineering of transcriptional regulators (TRs) that control several metabolic functions in modular patterns is an interesting strategy. However, only a handful of in silico model-added techniques are available for identifying the TR manipulation candidates, thus limiting its strain design application.

RESULTS

We developed hierarchical-Beneficial Regulatory Targeting (h-BeReTa) which employs a genome-scale metabolic model and transcriptional regulatory network (TRN) to identify the relevant TR targets suitable for strain improvement. We then applied this method to industrially relevant metabolites and cell factory hosts, Escherichia coli and Corynebacterium glutamicum. h-BeReTa suggested several promising TR targets, many of which have been validated through literature evidences. h-BeReTa considers the hierarchy of TRs in the TRN and also accounts for alternative metabolic pathways which may divert flux away from the product while identifying suitable metabolic fluxes, thereby performing superior in terms of global TR target identification.

CONCLUSIONS

In silico model-guided strain design framework, h-BeReTa, was presented for identifying transcriptional regulator targets. Its efficacy and applicability to microbial cell factories were successfully demonstrated via case studies involving two cell factory hosts, as such suggesting several intuitive targets for overproducing various value-added compounds.

Genome-scale model-driven strain design for dicarboxylic acid production in Yarrowia lipolytica

Background

Recently, there have been several attempts to produce long-chain dicarboxylic acids (DCAs) in various microbial hosts. Of these, Yarrowia lipolytica has great potential due to its oleaginous characteristics and unique ability to utilize hydrophobic substrates. However, Y. lipolytica should be further engineered to make it more competitive: the current approaches are mostly intuitive and cumbersome, thus limiting its industrial application.

Results

In this study, we proposed model-guided metabolic engineering strategies for enhanced production of DCAs in Y. lipolytica. At the outset, we reconstructed genome-scale metabolic model (GSMM) of Y. lipolytica (iYLI647) by substantially expanding the previous models. Subsequently, the model was validated using three sets of published culture experiment data. It was finally exploited to identify genetic engineering targets for overexpression, knockout, and cofactor modification by applying several in silico strain design methods, which potentially give rise to high yield production of the industrially relevant long-chain DCAs, e.g., dodecanedioic acid (DDDA). The resultant targets include (1) malate dehydrogenase and malic enzyme genes and (2) glutamate dehydrogenase gene, in silico overexpression of which generated additional NADPH required for fatty acid synthesis, leading to the increased DDDA fluxes by 48% and 22% higher, respectively, compared to wild-type. We further investigated the effect of supplying branched-chain amino acids on the acetyl-CoA turn-over rate which is key metabolite for fatty acid synthesis, suggesting their significance for production of DDDA in Y. lipolytica.

Conclusion

In silico model-based strain design strategies allowed us to identify several metabolic engineering targets for overproducing DCAs in lipid accumulating yeast, Y. lipolytica. Thus, the current study can provide a methodological framework that is applicable to other oleaginous yeasts for value-added biochemical production.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0542-5) contains supplementary material, which is available to authorized users.

In silico and in vivo stability analysis of a heterologous biosynthetic pathway for 1,4-butanediol production in metabolically engineered E. coli

Recently, several approaches have been published in order to develop a functional biosynthesis route for the non-natural compound 1,4-butanediol (BDO) in E. coli using glucose as a sole carbon source or starting from xylose. Among these studies, there was reported as high as 18 g/L product concentration achieved by industrial strains, however BDO production varies greatly in case of the reviewed studies. Our motivation was to build a simple heterologous pathway for this compound in E. coli and to design an appropriate cellular chassis based on a systemic biology approach, using constraint-based flux balance analysis and bi-level optimization for gene knock-out prediction. Thus, the present study reports, at the "proof-of concept" level, our findings related to model-driven development of a metabolically engineered E. coli strain lacking key genes for ethanol, lactate and formate production (ΔpflB, ΔldhA and ΔadhE), with a three-step biosynthetic pathway. We found this strain to produce a limited quantity of 1,4-BDO (.89 mg/L BDO under microaerobic conditions and .82 mg/L under anaerobic conditions). Using glycerol as carbon source, an approach, which to our knowledge has not been tackled before, our results suggest that further metabolic optimization is needed (gene-introductions or knock-outs, promoter fine-tuning) to address the redox potential imbalance problem and to achieve development of an industrially sustainable strain. Our experimental data on culture conditions, growth dynamics and fermentation parameters can consist a base for ongoing research on gene expression profiles and genetic stability of such metabolically engineered E. coli strains.

Flux-sum analysis identifies metabolite targets for strain improvement

Background

Rational design of microbial strains for enhanced cellular physiology through in silico analysis has been reported in many metabolic engineering studies. Such in silico techniques typically involve the analysis of a metabolic model describing the metabolic and physiological states under various perturbed conditions, thereby identifying genetic targets to be manipulated for strain improvement. More often than not, the activation/inhibition of multiple reactions is necessary to produce a predicted change for improvement of cellular properties or states. However, as it is more computationally cumbersome to simulate all possible combinations of reaction perturbations, it is desirable to consider alternative techniques for identifying such metabolic engineering targets.

Results

In this study, we present the modified version of previously developed metabolite-centric approach, also known as flux-sum analysis (FSA), for identifying metabolic engineering targets. Utility of FSA was demonstrated by applying it to Escherichia coli, as case studies, for enhancing ethanol and succinate production, and reducing acetate formation. Interestingly, most of the identified metabolites correspond to gene targets that have been experimentally validated in previous works on E. coli strain improvement. A notable example is that pyruvate, the metabolite target for enhancing succinate production, was found to be associated with multiple reaction targets that were only identifiable through more computationally expensive means. In addition, detailed analysis of the flux-sum perturbed conditions also provided valuable insights into how previous metabolic engineering strategies have been successful in enhancing cellular physiology.

Conclusions

The application of FSA under the flux balance framework can identify novel metabolic engineering targets from the metabolite-centric perspective. Therefore, the current approach opens up a new research avenue for rational design and engineering of industrial microbes in the field of systems metabolic engineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0198-3) contains supplementary material, which is available to authorized users.

In silico model-driven cofactor engineering strategies for improving the overall NADP(H) turnover in microbial cell factories

Optimizing the overall NADPH turnover is one of the key challenges in various value-added biochemical syntheses. In this work, we first analyzed the NADPH regeneration potentials of common cell factories, including Escherichia coli, Saccharomyces cerevisiae, Bacillus subtilis, and Pichia pastoris across multiple environmental conditions and determined E. coli and glycerol as the best microbial chassis and most suitable carbon source, respectively. In addition, we identified optimal cofactor specificity engineering (CSE) enzyme targets, whose cofactors when switched from NAD(H) to NADP(H) improve the overall NADP(H) turnover. Among several enzyme targets, glyceraldehyde-3-phosphate dehydrogenase was recognized as a global candidate since its CSE improved the NADP(H) regeneration under most of the conditions examined. Finally, by analyzing the protein structures of all CSE enzyme targets via homology modeling, we established that the replacement of conserved glutamate or aspartate with serine in the loop region could change the cofactor dependence from NAD(H) to NADP(H).

Applications of genome-scale metabolic network model in metabolic engineering

Genome-scale metabolic network model (GEM) is a fundamental framework in systems metabolic engineering. GEM is built upon extensive experimental data and literature information on gene annotation and function, metabolites and enzymes so that it contains all known metabolic reactions within an organism. Constraint-based analysis of GEM enables the identification of phenotypic properties of an organism and hypothesis-driven engineering of cellular functions to achieve objectives. Along with the advances in omics, high-throughput technology and computational algorithms, the scope and applications of GEM have substantially expanded. In particular, various computational algorithms have been developed to predict beneficial gene deletion and amplification targets and used to guide the strain development process for the efficient production of industrially important chemicals. Furthermore, an Escherichia coli GEM was integrated with a pathway prediction algorithm and used to evaluate all possible routes for the production of a list of commodity chemicals in E. coli. Combined with the wealth of experimental data produced by high-throughput techniques, much effort has been exerted to add more biological contexts into GEM through the integration of omics data and regulatory network information for the mechanistic understanding and improved prediction capabilities. In this paper, we review the recent developments and applications of GEM focusing on the GEM-based computational algorithms available for microbial metabolic engineering.

Transhydrogenase promotes the robustness and evolvability of E. coli deficient in NADPH production

Metabolic networks revolve around few metabolites recognized by diverse enzymes and involved in myriad reactions. Though hub metabolites are considered as stepping stones to facilitate the evolutionary expansion of biochemical pathways, changes in their production or consumption often impair cellular physiology through their system-wide connections. How does metabolism endure perturbations brought immediately by pathway modification and restore hub homeostasis in the long run? To address this question we studied laboratory evolution of pathway-engineered Escherichia coli that underproduces the redox cofactor NADPH on glucose. Literature suggests multiple possibilities to restore NADPH homeostasis. Surprisingly, genetic dissection of isolates from our twelve evolved populations revealed merely two solutions: (1) modulating the expression of membrane-bound transhydrogenase (mTH) in every population; (2) simultaneously consuming glucose with acetate, an unfavored byproduct normally excreted during glucose catabolism, in two subpopulations. Notably, mTH displays broad phylogenetic distribution and has also played a predominant role in laboratory evolution of Methylobacterium extorquens deficient in NADPH production. Convergent evolution of two phylogenetically and metabolically distinct species suggests mTH as a conserved buffering mechanism that promotes the robustness and evolvability of metabolism. Moreover, adaptive diversification via evolving dual substrate consumption highlights the flexibility of physiological systems to exploit ecological opportunities.

Genome engineering for improved recombinant protein expression in Escherichia coli

A metabolic engineering perspective which views recombinant protein expression as a multistep pathway allows us to move beyond vector design and identify the downstream rate limiting steps in expression. In E.coli these are typically at the translational level and the supply of precursors in the form of energy, amino acids and nucleotides. Further recombinant protein production triggers a global cellular stress response which feedback inhibits both growth and product formation. Countering this requires a system level analysis followed by a rational host cell engineering to sustain expression for longer time periods. Another strategy to increase protein yields could be to divert the metabolic flux away from biomass formation and towards recombinant protein production. This would require a growth stoppage mechanism which does not affect the metabolic activity of the cell or the transcriptional or translational efficiencies. Finally cells have to be designed for efficient export to prevent buildup of proteins inside the cytoplasm and also simplify downstream processing. The rational and the high throughput strategies that can be used for the construction of such improved host cell platforms for recombinant protein expression is the focus of this review.