Abstract
Genome-scale metabolic models (GEMs) have become a popular tool for systems biology, and they have been used in many fields such as industrial biotechnology and systems medicine. Since more and more studies are being conducted using GEMs, they have recently received considerable attention. In this review, we introduce the basic concept of GEMs and provide an overview of their applications in biotechnology, systems medicine, and some other fields. In addition, we describe the general principle of the applications and analyses built on GEMs. The purpose of this review is to introduce the application of GEMs in biological analysis and to promote its wider use by biologists.
Introduction
Genome-scale metabolic models (GEMs) are reconstructions of the metabolic networks of many kinds of cells, including those of microorganisms, plants, and mammals. In some cases, GEMs could represent the whole tissue or body of a multicellular organism. In these metabolic networks, the gene-protein-reaction (GPR) relationships are annotated. In addition, all the reactions in GEMs are mass- and energy-balanced, ensuring stoichiometric balance. Thus, GEMs enable researchers to conduct system-level metabolic response analysis and flux simulation, which is not possible using general metabolic pathway databases such as KEGG. Furthermore, since GPR relationships are included in GEMs, other omics data such as transcriptomic and proteomic data could be systematically integrated into GEMs. Thus, GEM-based multi-omic analyses are more informative with stoichiometric balance and could possibly provide deeper biological insights.
In the past 15 years, GEMs have garnered considerable research attention. In 2000, the first GEM, a model of Escherichia coli MG1655, was reported (Edwards and Palsson, ). A few years later, a yeast GEM was published (Doerks et al., ), thus initiating a new era for systems biology. In the beginning, researchers tried to use GEM-based in silico simulations to guide the rational design of industrial microorganisms (hereafter referred to as in silico metabolic engineering). In 2003, a method called OptKnock (Burgard et al., ) was published and it employed a bi-level optimization program to search for reaction knockout targets that would yield overproduction of a desired biochemical while maintaining optimal growth. Following that, a series of in silico metabolic engineering methods were developed for various gene manipulations other than knock-out (Pharkya et al., ; Pharkya and Maranas, ; Choi et al., ; Ranganathan et al., ; Park et al., ; Chowdhury et al., ; Mahalik et al., ), leading to a marked expansion in the usage of GEMs. Furthermore, many of the in silico metabolic engineering methods were experimentally validated (Fong et al., ; Izallalen et al., ; Asadollahi et al., ; Brochado et al., ; Choi et al., ; Yim et al., 2011; Xu et al., 2011; Park et al., ; Ranganathan et al., ; Otero et al., ; Kim et al., ), which showed the power of GEM-based applications. With the development of systems biology, GEMs were also used as scaffolds for systematic integration of omics data because GEMs could be used to reconstruct the relationship among genes, enzymes, and metabolism. Numerous algorithms have been developed to integrate various types of omics data such as thermodynamics (Henry et al., ), transcriptomics/proteomics (Becker and Palsson, ; Colijn et al., ; Zur et al., 2010), fluxomics (Wiback et al., 2004), and metabolomics (Cakir et al., ). In return, the integration of omics data could improve the prediction of GEMs. More recently, GEM has been applied to systems medicine. Since the reconstruction of the first global GEM for humans, Recon 1, which was established in 2007 (Duarte et al., ), researchers have started to explore the possibility of clinical applications of GEMs and have reported several successful cases (Agren et al., ; Gatto et al., ; Jerby-Arnon et al., ). In fact, GEMs and their applications have received considerable attention recently.
Although GEMs are becoming increasingly popular, they are not easy to understand or use by non-experts. The complex code and script usually used for GEM-based computational applications and analyses are not readily available to the community of biologists, greatly hampering the wide usage of GEMs. In this review, we describe the key concepts and assumptions of GEMs. In addition, we describe the general principle of the applications and analyses built on GEMs. The information presented here is expected to promote the spread of GEM usage by biologists.
Basic concept of GEMs
As mentioned above, GEMs are metabolic networks. Figure 1A shows a partly visualized glycolysis pathway in a GEM of E. coli, and within this part, we can see that metabolites are linked with each other by reactions, which are associated with enzymes, which are encoded by genes. It should be noted that the stoichiometric coefficient in metabolic reactions in Figure 1A (as shown in Figure 1B) could not be visualized in a graph. Therefore, GEMs employ a stoichiometric matrix (S matrix) to represent all the coefficients in metabolic reactions (Figure 1C). In the S matrix, the ijth element represents the stoichiometric coefficient of the ith metabolite in the jth reaction in the GEM. If the coefficient is positive, the metabolite is produced; otherwise, it's consumed. In addition, the GPR relationships in GEMs are simplified into a two-dimensional binary matrix showing the association between genes and reactions (Figure 1D), in which the ijth element is one if the ith reaction is associated with the jth gene, and it's zero if they aren't associated.
Figure 1
GEMs have several notable features: (1) They are collections of existing knowledge of the metabolism of a specific organism, and in most GEM-based applications, it's assumed that the metabolic network is complete, with very few exceptions, such as for gap finding and gap filling (Latendresse, ). (2) They are stoichiometric-balanced networks, which means mass as well as energy balance, reduction, and proton balance are well considered. (3) GPR relationships are annotated in GEMs, but the interactions are not quantitatively described. (4) Even though GEMs describe the metabolism, concentrations of metabolites are not directly included and flux balance analysis (FBA; Orth et al., ) is employed for flux simulations, which assumes that there is no (unexpected) accumulation of metabolites within GEMs.
Using GEMs for essentiality and synthetic lethality analysis
As mentioned above, since GEMs are complete metabolic networks, they can be used for gene/reaction essentiality analysis (EA; Edwards and Palsson, ). In general, EA identifies all essential genes or reactions whose knockout will disable a specific biological function through FBA. EA could be easily implemented in silico using GEMs by enumerating all single gene/reaction knockouts and testing whether their biological objectives are still functioning. In addition, synthetic lethality analysis (SLA), which scans for combinatory knockouts of multiple reactions/genes that lead to blocking of the target biological function, could also be implemented in a similar way. And recently, several methods have been developed to perform advanced SLA efficiently (Suthers et al., ; von Kamp and Klamt, 2014; Pratapa et al., ; Zhang et al., 2015).
It's generally believed that gene/reaction EA could be performed by topologic analysis of the metabolic network. However, since the stoichiometric coefficients are absent in topologic metabolic networks, they're less accurate. For example, Figure 2 shows the topologic network of the toy model from Figure 1. Based on its topologic properties, this metabolic work can use D-glucose-6-phosphate, NAD, and phosphate as substrates and produce 3-phospho-D-glycerate, NADH, and a proton. However, this pathway always consumes more ADP than it produces, and produces more ATP than it consumes. Therefore, this pathway will be blocked without ADP supplementation and this finding was not possible by topologic analysis.
Figure 2
Essentially, if a GEM is well established, its EA and SLA results could be very accurate. For example, in the most used E. coli and S. cerevisiae GEMs, around 90% of the predicted essential genes have been validated in vivo (true-negative; Feist et al., ; Heavner et al., ). This is within expectation, because if a function is blocked in silico, it's very unlikely that there could be a complimentary solution in vivo to recover it. The explanation for the very few false-negative predictions (negative growth in silico and positive growth in vivo) is that there's a knowledge gap, such as unknown enzyme or unknown function of an existing enzyme, which leads to the underestimation of the capability of the GEM. On the other hand, even if the GEMs are 100% complete, there may still be false-positive predictions since the missing information of regulation and protein (enzyme) efficiency could lead to extra constraints to GEMs, thereby rendering a nonessential reaction/gene in silico essential in vivo. It's worth mentioning that, after a certain period of adaptive evolution, a false-positive knockout could become nonessential in vivo again (Patil et al., ). EA and SLA have mainly been used to validate newly constructed GEMs and in recent years, EA and SLA were applied to study of systems medicine (see Section Using GEMs in Studies of Systems Medicine).
Using GEMs as scaffolds for multi omics data integration and interpretation
Recently, increasing volumes of transcriptomic, proteomic, and metabolomics data are becoming publically available, and it's believed that GEMs are good scaffolds to make use of these multi omics data. In GEMs, omics data could be quantitatively integrated as constraints for metabolic fluxes, thereby allowing systematic and quantitative evaluation of these data, which was not possible using traditional metabolic networks. This is the most significant advantage of using GEMs as scaffolds.
Although, GEMs are metabolic networks, the most used omic data for GEMs are transcriptomic and proteomic. This is because the technic is really advancing in the field and makes large number of high quality transcriptomic and proteomic data available. However, since the GPR relationships are qualitative in GEMs (Figure 1C), one needs to make assumptions to define the quantitative relationship between gene/protein expression and metabolic fluxes when integrating transcriptomic or proteomic data into GEMs. This is problematic because the complicated relation between fluxes and expression level of genes and enzymes in vivo are unlikely to be captured by a general assumption (MacHado and HerrgÃ¥rd, ). On the other hand, there're many well-defined approaches to integrate fluxomics and metabolomics, data (Khodayari et al., ; MartÃn et al., ; Miskovic et al., ). However, it's very difficult (if not impossible) to get genome scale data of them. Hence, we suggest that even though omics data are integrated, one should be skeptical about the quantitative results of simulations or predictions from GEMs.
Nonetheless, we believe that it is better to qualitatively interpret the omics data using GEMs. For instance, it would be much more reliable to use omics data to determine the presence or absence of reactions and to construct high-quality and specific GEMs (Zur et al., 2010; Agren et al., , ; Mardinoglu et al., ; Yizhak et al., 2014). In addition, many researchers started to integrate significance of differential expression of genes with GEMs rather than their quantitative expression to interpret the biological information behind omic data (Patil and Nielsen, ; Cakir et al., ; Jensen and Papin, ; Fang et al., ; Navid and Almaas, ). Moreover, qualitative interpretation of omics data with GEMs have recently been applied to systems medicine (see Section Using GEMs in Studies of Systems Medicine). These studies demonstrated the usefulness of GEMs as scaffolds. In short, we suggest that GEMs are powerful platforms for integration of omics data for gaining biological insights rather than quantitative results.
Using GEMs for in silico metabolic engineering
Using GEMs for in silico metabolic engineering has been a widely discussed topic for years. It's generally believed that GEM-based methods could predict gene modification strategies for overproduction of desired biochemicals and thus, accelerate the overall metabolic engineering process. In the last decade, various kinds of in silico metabolic engineering methods had been developed and many of them were applied experimentally (Kim et al., ; Long et al., ; MacHado and Herrgård, ).
Although in silico metabolic engineering methods seemed quite different from each other, they follow a similar procedure: (1) they define what a desired strain is and (2) identify approaches that push the wild-type strain to become the desired one. So far, a variety of approaches were used in in silico metabolic engineering, such as reaction/gene knock-out (Burgard et al., ; Patil et al., ; Kim et al., ; Ren et al., ; Ruckerbauer et al., ; Zhang et al., 2015), overexpression/suppression (Pharkya and Maranas, ; Choi et al., ; Ranganathan et al., ; Park et al., ; Chowdhury et al., ), foreign pathway knock-in (Pharkya et al., ), and swapping the co-factor for a target enzyme (NADH to NADPH or vice versa; King and Feist, ). However, the methods for knock-out identification are the majority since a knockout is much easier to define in silico than up-/down-regulation of genes as mentioned before. On the other hand, different methods could have independent definition of desired strains. For instance, some of the methods define the desired strain by simply setting thresholds for growth and production, respectively, and others could define the desired strain following some biological assumptions (Edwards et al., ; Segrè et al., ).
Interestingly, methods pursuing different type of desired strains could all lead to experimentally valid strategies for metabolic engineering (Fong et al., ; Trinh et al., ; Fowler et al., ; Choi et al., ; Yim et al., 2011; Ng et al., ; Nocon et al., ), but the production of target products predicted in silico seldom achieved in vivo. The explanation to this is complicated, and could come from both the computational and experimental side. However, one of the key reasons should be that GEM with only metabolic network is not enough to quantitatively predict the behavior of strains in vivo. In conclusion, we suggested that all kinds of in silico metabolic engineering methods are instructive, but it's better to use them for gaining information rather than to develop exact strategies.
Using GEMs in studies of systems medicine
Using GEMs for systems medicine studies have recently been highlighted (Mardinoglu and Nielsen, ; Yizhak et al., 2015). GEMs simulate the human metabolism in a holistic way, and this greatly advances systems medicine studies by enabling systematic evaluation of metabolic feature of human disease. Great efforts had been made in reconstructing GEMs of human, and there're now several publically available generic human metabolic networks such as Recon 1, Recon 2, EHMN, and HMR (Duarte et al., ; Ma et al., ; Agren et al., ; Thiele et al., ). In addition, since the technology is advancing, tissue specific or cell specific genomic, proteomic and transcriptomic data are becoming available (Cancer Genome Atlas Research Network, ; Uhlén et al., 2015). These led to rapid development in reconstruction of high quality tissue or cancer specific GEMs (Zur et al., 2010; Agren et al., , ; Mardinoglu et al., ) and, therefore, enabled more confident interpretation of metabolism of diseases.
For instance, cancer specific GEMs together with EA and SLA analysis were recently used for identification of oncogenes/metabolites and biomarkers for diagnosing specific cancer (Agren et al., ; Jerby-Arnon et al., ; Gatto et al., ; Gatto and Nielsen, ). Since this procedure mainly uses the true-negative part of EA and SLA, the analysis could be highly reliable. For example, (Agren et al., ) identified 101 drug targets for liver cancer treatment; and 83 of them are currently in use or have shown strong correlation with cancer progression. In addition, together with multi-omic data, GEMs were used to find the mechanistic explanation of various diseases. By interpreting clinical omic data with GEMs, the mechanistic understanding of non-alcoholic fat liver disease and type two diabetes were reported (Mardinoglu et al., ; Väremo et al., 2015). Moreover, GEMs were also used to explore the effect of microbiota (Ji and Nielsen, ). By simulate and predict the interaction of gut microbiota and their effect on hosts, several recent studies revealed that microbiota modulate the amino acid and glutathione metabolism of their host (Shoaie et al., , ; Mardinoglu et al., ). These exciting studies exhibited the great potential of GEMs in the field of systems medicine, and hopefully there would be much more excellent works coming out.
Discussion
GEMs are very useful platforms and tools for systems biology, but they're still very young compared to traditional ones. Fluxes of reactions could be quantitatively simulated using GEMs, although caution should be exercised before drawing conclusions based on simulated fluxes owing to the huge solution space of GEMs (Reed, ). Although solution space could be reduced by adding constraints through integration of omics data, it would be better to gain biological insights by qualitative interpretation of omics data rather than quantitative fluxes.
In order to achieve accurate quantitative prediction, the scope of GEMs should be expanded. The establishment of ME-models set a good example for this (Thiele et al., , ). In ME-models, the interaction of genes (mRNA), enzymes, and metabolic fluxes are quantitatively expressed, enabling proper integration of transcriptomic and proteomic data. However, it is still difficult to integrate metabolomics data into ME-models. A potential option to integrate metabolite concentration into GEMs is cybernetic modeling. However, to date, there has been no study on genome-scale cybernetic modeling because there are too many parameters to simulate, making it computationally infeasible.
In general, no model is perfect. Genome-scale modeling methods are still under development and have several drawbacks. In addition, it has been recently reported that many published GEMs are of low qualities (Chindelevitch et al., ; Ravikrishnan and Raman, ). Therefore, they should be used with caution. As concluded in this review, GEMs are more suitable for qualitative applications at this stage, such as EA and SLA analysis. When using GEMs for quantitative applications such as in silico metabolic engineering, one should be aware of the key assumption behind the method and take the results as instructions. However, it should also be noted that, GEMs are open platforms and have great potential in a wide array of applications. Currently, GEMs are used for simulating the interactions between multiple organisms, multiple tissues (Bordbar et al., ), and even between microbiota and human tissues. On the other hand, EA and SLA were developed years ago, but they were not used in the discovery of anti-cancer drugs until recent years. These are good examples of how to explore novel applications based on classical methods. Thus, in future, GEMs can be expected to be more widely used in biotechnology, bioengineering, and many other fields.
Statements
Author contributions
CZ conducted the writing of this paper, and QH modified and edited it.
Acknowledgments
This work was funded by National Basic Research Program of China (973 Program) (2012CB721101), and China Scholarship Council.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AgrenR.BordelS.MardinogluA.PornputtapongN.NookaewN.NielsenJ. (2012). Reconstruction of genome-scale active metabolic networks for 69 human cell types and 16 cancer types using INIT. PLoS Comput. Biol.8:e1002518. 10.1371/journal.pcbi.1002518
2
AgrenR.MardinogluA.AsplundA.KampfC.UhlenM.NielsenJ. (2014). Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling. Mol. Syst. Biol. 10, 721. 10.1002/msb.145122
3
AsadollahiM. A.MauryJ.PatilK. R.SchalkM.ClarkA.NielsenJ. (2009). Enhancing sesquiterpene production in Saccharomyces cerevisiae through in silico driven metabolic engineering. Metab. Eng.11, 328–334. 10.1016/j.ymben.2009.07.001
4
BeckerS. A.PalssonB. O. (2008). Context-specific metabolic networks are consistent with experiments. PLoS Comput. Biol. 4:e1000082. 10.1371/journal.pcbi.1000082
5
BordbarA.FeistA.Usaite-BlackR.WoodcockJ.PalssonB.FamiliI. (2011). A multi-tissue type genome-scale metabolic network for analysis of whole-body systems physiology. BMC Syst. Biol.5:180. 10.1186/1752-0509-5-180
6
BrochadoA.MatosC.MøllerB.HansenJ.MortensenU.PatilK. (2010). Improved vanillin production in baker's yeast through in silico design. Microb. Cell Factories9:84. 10.1186/1475-2859-9-84
7
BurgardA. P.PharkyaP.MaranasC. D. (2003). Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization. Biotechnol. Bioeng.84, 647–657. 10.1002/bit.10803
8
CakirT.PatilK. R.OnsanZ. I.UlgenK. O.KirdarB.NielsenJ. (2006). Integration of metabolome data with metabolic networks reveals reporter reactions. Mol. Syst. Biol. 2:50. 10.1038/msb4100085
9
Cancer Genome Atlas Research Network, (2008). Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature455, 1061–1068. 10.1038/nature07385
10
ChindelevitchL.TriggJ.RegevA.BergerB. (2014). An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models. Nat. Commun.5:4893. 10.1038/ncomms5893
11
ChoiH. S.LeeS. Y.KimT. Y.WooH. M. (2010). In silico identification of gene amplification targets for improvement of lycopene production. Appl. Environ. Microbiol.76, 3097–3105. 10.1128/AEM.00115-10
12
ChowdhuryA.ZomorrodiA. R.MaranasC. D. (2014). k-OptForce: integrating kinetics with flux balance analysis for strain design. PLoS Comput. Biol.10:e1003487. 10.1371/journal.pcbi.1003487
13
ColijnC.BrandesA.ZuckerJ.LunD. S.WeinerB.FarhatM. R.et al. (2009). Interpreting expression data with metabolic flux models: predicting mycobacterium tuberculosis mycolic acid production. PLoS Comput. Biol.5:e1000489. 10.1371/journal.pcbi.1000489
14
DoerksT.CopleyR. R.SchultzJ.PontingC. P.BorkP. (2002). Systematic identification of novel protein domain families associated with nuclear functions. Genome Res.12, 47–56. 10.1101/gr.203201
15
DuarteN. C.BeckerS. A.JamshidiN.ThieleI.MoM. L.VoT. D.et al. (2007). Global reconstruction of the human metabolic network based on genomic and bibliomic data. Proc. Natl. Acad. Sci. U.S.A.104, 1777–1782. 10.1073/pnas.0610772104
16
EdwardsJ. S.IbarraR. U.PalssonB. O. (2001). In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat. Biotech. 19, 125–130. 10.1038/84379
17
EdwardsJ. S.PalssonB. O. (2000). The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. Proc.Natl Acad. Sci. U.S.A.97, 5528–5533. 10.1073/pnas.97.10.5528
18
FongS. S.BurgardA. P.HerringC. D.KnightE. M.BlattnerF. R.MaranasC. D.et al. (2005). In silico design and adaptive evolution of Escherichia coli for production of lactic acid. Biotechnol. Bioeng.91, 643–648. 10.1002/bit.20542
19
FangX.WallqvistA.ReifmanJ. (2012). Modeling phenotypic metabolic adaptations of Mycobacterium tuberculosis H37Rv under Hypoxia. PLoS Comput. Biol.8:e1002688. 10.1371/journal.pcbi.1002688
20
FeistA. M.HenryC. S.ReedJ. L.KrummenackerM.JoyceA. R.KarpP. D.et al. (2007). A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol. Syst. Biol.3, 121. 10.1038/msb4100155
21
FowlerZ. L.GikandiW. W.KoffasM. A. (2009). Increased malonyl coenzyme a biosynthesis by tuning the Escherichia coli metabolic network and its application to flavanone production. Appl. Environ. Microbiol.75, 5831–5839. 10.1128/AEM.00270-09
22
GattoF.MiessH.SchulzeA.NielsenJ. (2015). Flux balance analysis predicts essential genes in clear cell renal cell carcinoma metabolism.Sci. Rep.5:10738. 10.1038/srep10738
23
GattoF.NielsenJ. (2015). In search for symmetries in the metabolism of cancer. Wiley Interdiscip. Rev. Syst. Biol. Med. [Epub ahead of print]. 10.1002/wsbm.1321
24
GattoF.NookaewI.NielsenJ. (2014). Chromosome 3p loss of heterozygosity is associated with a unique metabolic network in clear cell renal carcinoma. Proc. Natl. Acad. Sci. U.S.A.111, E866–E875. 10.1073/pnas.1319196111
25
HeavnerB. D.SmallboneK.PriceN. D.WalkerL. P. (2013). Version 6 of the consensus yeast metabolic network refines biochemical coverage and improves model performance. Database (Oxford)2013:bat059. 10.1093/database/bat059
26
HenryC. S.BroadbeltL. J.HatzimanikatisV. (2007). Thermodynamics-based metabolic flux analysis. Biophys. J.92, 1792–1805. 10.1529/biophysj.106.093138
27
IzallalenM.MahadevanR.BurgardA.PostierB.DidonatoR.Jr.SunJ.et al. (2008). Geobacter sulfurreducens strain engineered for increased rates of respiration. Metab. Eng. 10, 267–275. 10.1016/j.ymben.2008.06.005
28
JensenP. A.PapinJ. A. (2011). Functional integration of a metabolic network model and expression data without arbitrary thresholding. Bioinformatics27, 541–547. 10.1093/bioinformatics/btq702
29
Jerby-ArnonL.PfetzerN.WaldmanY. Y.McGarryL.JamesD.ShanksE.et al. (2014). Predicting cancer-specific vulnerability via data-driven detection of synthetic lethality. Cell158, 1199–1209. 10.1016/j.cell.2014.07.027
30
JiB.NielsenJ. (2015). From next-generation sequencing to systematic modeling of the gut microbiome. Front. Genet. 6:219. 10.3389/fgene.2015.00219
31
KhodayariA.ZomorrodiA. R.LiaoJ. C.MaranasC. D. (2014). A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data. Metab. Eng.25, 50–62. 10.1016/j.ymben.2014.05.014
32
KimB.KimW. J.KimD. I.LeeS. Y. (2015). Applications of genome-scale metabolie network model in metabolic engineering. J. Ind. Microbiol. Biotechnol. 42, 339–348. 10.1007/s10295-014-1554-9
33
KimJ.ReedJ. L.MaraveliasC. T. (2012). Large-scale bi-level strain design approaches and mixed-integer programming solution techniques. PLoS ONE6:e24162. 10.1371/journal.pone.0024162
34
KimM.Sang YiJ.KimJ.KimJ. N.KimM. W.KimB. G. (2014). Reconstruction of a high-quality metabolic model enables the identification of gene overexpression targets for enhanced antibiotic production in Streptomyces coelicolor A3(2). Biotechnol. J.9, 1185–1194. 10.1002/biot.201300539
35
KingZ. A.FeistA. M. (2013). Optimizing cofactor specificity of oxidoreductase enzymes for the generation of microbial production strains—OptSwap. Ind. Biotechnol. 9, 236–246. 10.1089/ind.2013.0005
36
LatendresseM. (2014). Efficiently gap-filling reaction networks. BMC Bioinformatics15:225. 10.1186/1471-2105-15-225
37
LongM. R.OngW. K.ReedJ. L. (2015). Computational methods in metabolic engineering for strain design. Curr. Opin. Biotech. 34, 135–141. 10.1016/j.copbio.2014.12.019
38
MaH.SorokinA.MazeinA.SelkovA.SelkovE.DeminO.et al. (2007). The Edinburgh human metabolic network reconstruction and its functional analysis. Mol. Syst. Biol. 3:135. 10.1038/msb4100177
39
MacHadoD.HerrgårdM. (2014). Systematic evaluation of methods for integration of transcriptomic data into constraint-based models of metabolism. PLoS Comput. Biol.10:e1003580. 10.1371/journal.pcbi.1003580
40
MacHadoD.HerrgårdM. (2015). Co-evolution of strain design methods based on flux balance and elementary mode analysis. Metab. Eng. Commun. 2, 85–92. 10.1016/j.meteno.2015.04.001
41
MahalikS.SharmaA. K.MukherjeeK. J. (2014). Genome engineering for improved recombinant protein expression in Escherichia coli. Microb. Cell Factories13, 177. 10.1186/s12934-014-0177-1
42
MardinogluA.AgrenR.KampfC.AsplundA.NookaewI.JacobsonP.et al. (2013). Integration of clinical data with a genome-scale metabolic model of the human adipocyte.Mol. Syst. Biol.9. 10.1038/msb.2013.5
43
MardinogluA.AgrenR.KampfC.AsplundA.UhlenM.NielsenJ. (2014). Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease.Nat. Commun.5. 10.1038/ncomms4083
44
MardinogluA.NielsenJ. (2015). New paradigms for metabolic modeling of human cells.Curr. Opin. Biotechnol.34, 91–97. 10.1016/j.copbio.2014.12.013
45
MardinogluA.ShoaieS.BergentallM.GhaffariP.ZhangC.LarssonE.et al. (2015). The gut microbiota modulates host amino acid and glutathione metabolism in mice. Mol. Syst. Biol.11, 834. 10.15252/msb.20156487
46
MartÃnH. G.KumarV. S.WeaverD.GhoshA.ChubukovV.MukhopadhyayA.et al. (2015). A method to constrain genome-scale models with 13C labeling data. PLoS Comput. Biol.11:e1004363. 10.1371/journal.pcbi.1004363
47
MiskovicL.TokicM.FengosG.HatzimanikatisV. (2015). Rites of passage: requirements and standards for building kinetic models of metabolic phenotypes. Curr. Opin. Biotech. 36, 146–153. 10.1016/j.copbio.2015.08.019
48
NavidA.AlmaasE. (2012). Genomic-level transcription data of Yersinia pestis analyzed with a new metabolic constraint-based approarch. BMC Syst. Biol.6:150. 10.1186/1752-0509-6-150
49
NgC. Y.JungM.LeeJ.OhM. K. (2012). Production of 2, 3-butanediol in Saccharomyces cerevisiae by in silico aided metabolic engineering. Microb. Cell Factories11:68. 10.1186/1475-2859-11-68
50
NoconJ.SteigerM. G.PfefferM.SohnS. B.KimT. Y.MaurerM.et al. (2014). Model based engineering of Pichia pastoris central metabolism enhances recombinant protein production. Metab. Eng. 24, 129–138. 10.1016/j.ymben.2014.05.011
51
OrthJ. D.ThieleI.PalssonB. O. (2010). What is flux balance analysis?Nat. Biotechnol.28, 245–248. 10.1038/nbt.1614
52
OteroJ. M.CiminiD.PatilK. R.PoulsenS. G.OlssonL.NielsenJ. (2013). Industrial systems biology of Saccharomyces cerevisiae enables novel succinic acid cell factory. PLoS ONE8:e54144. 10.1371/journal.pone.0054144
53
ParkJ. M.ParkH. M.KimW. J.KimH. U.KimT. Y.LeeS. Y. (2012). Flux variability scanning based on enforced objective flux for identifying gene amplification targets. BMC Syst. Biol.6, 106–106. 10.1186/1752-0509-6-106
54
PatilK. R.NielsenJ. (2005). Uncovering transcriptional regulation of metabolism by using metabolic network topology. Proc. Natl. Acad. Sci. U.S.A.102, 2685–2689. 10.1073/pnas.0406811102
55
PatilK. R.RochaI.FörsterJ.NielsenJ. (2005). Evolutionary programming as a platform for in silico metabolic engineering. BMC Bioinformatics6:308. 10.1186/1471-2105-6-308
56
PharkyaP.BurgardA. P.MaranasC. D. (2004). OptStrain: a computational framework for redesign of microbial production systems. Genome Res.14, 2367–2376. 10.1101/gr.2872004
57
PharkyaP.MaranasC. D. (2006). An optimization framework for identifying reaction activation/inhibition or elimination candidates for overproduction in microbial systems. Metab. Eng.8, 1–13. 10.1016/j.ymben.2005.08.003
58
PratapaA.BalachandranS.RamanK. (2015). Fast-SL: an efficient algorithm to identify synthetic lethal sets in metabolic networks. Bioinformatics31, 3299–3305. 10.1093/bioinformatics/btv352
59
RanganathanS.SuthersP. F.MaranasC. D. (2010). OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput. Biol.6:e1000744. 10.1371/journal.pcbi.1000744
60
RanganathanS.TeeT. W.ChowdhuryA.ZomorrodiA. R.YoonJ. M.FuY.et al. (2012). An integrated computational and experimental study for overproducing fatty acids in Escherichia coli. Metab. Eng.14, 687–704. 10.1016/j.ymben.2012.08.008
61
RavikrishnanA.RamanK. (2015). Critical assessment of genome-scale metabolic networks: the need for a unified standard. Brief. Bioinform. 16, 1057–1068. 10.1093/bib/bbv003
62
ReedJ. L. (2012). Shrinking the metabolic solution space using experimental datasets. PLoS Comput. Biol.8:e1002662. 10.1371/journal.pcbi.1002662
63
RenS.ZengB.QianX. (2013). Adaptive bi-level programming for optimal gene knockouts for targeted overproduction under phenotypic constraints.BMC Bioinformatics14(Suppl. 2):S17. 10.1186/1471-2105-14-s2-s17
64
RuckerbauerD. E.JungreuthmayerC.ZanghelliniJ. (2014). Design of optimally constructed metabolic networks of minimal functionality. PLoS ONE9:e92583. 10.1371/journal.pone.0092583
65
SegrèD.VitkupD.ChurchG. M. (2002). Analysis of optimality in natural and perturbed metabolic networks. Proc. Natl. Acad. Sci. U.S.A.99, 15112–15117. 10.1073/pnas.232349399
66
ShoaieS.GhaffariP.Kovatcheva-DatcharyP.MardinogluA.SenP.Pujos-GuillotE.et al. (2015). Quantifying diet-induced metabolic changes of the human gut microbiome. Cell Metab.22, 320–331. 10.1016/j.cmet.2015.07.001
67
ShoaieS.KarlssonF.MardinogluA.NookaewI.BordelS.NielsenJ. (2013). Understanding the interactions between bacteria in the human gut through metabolic modeling. Sci. Rep. 3, 2532. 10.1038/srep02532
68
SuthersP. F.ZomorrodiA.MaranasC. D. (2009). Genome-scale gene/reaction essentiality and synthetic lethality analysis. Mol. Syst. Biol.5, 301. 10.1038/msb.2009.56
69
ThieleI.FlemingR. M. T.QueR.BordbarA.DiepD.PalssonB. O. (2012). Multiscale modeling of metabolism and macromolecular synthesis in E. coli and its application to the evolution of codon usage. PLoS ONE7:e45635. 10.1371/journal.pone.0045635
70
ThieleI.JamshidiN.FlemingR. M.PalssonB. O. (2009). Genome-scale reconstruction of Escherichia coli's transcriptional and translational machinery: a knowledge base, its mathematical formulation, and its functional characterization. PLoS Comput. Biol. 5:e1000312. 10.1371/journal.pcbi.1000312
71
ThieleI.SwainstonN.FlemingR. M.HoppeA.SahooS.AurichM. K.et al. (2013). A community-driven global reconstruction of human metabolism.Nat. Biotech.31, 419–425. 10.1038/nbt.2488
72
TrinhC. T.UnreanP.SriencF. (2008). Minimal Escherichia coli cell for the most efficient production of ethanol from hexoses and pentoses. Appl. Environ. Microbiol.74, 3634–3643. 10.1128/AEM.02708-07
73
UhlénM.FagerbergL.HallströmB. M.LindskogC.OksvoldP.MardinogluA.et al. (2015). Tissue-based map of the human proteome. Science347:1260419. 10.1126/science.1260419
74
VäremoL.ScheeleC.BroholmC.MardinogluA.KampfC.AsplundA.et al. (2015). Proteome- and transcriptome-driven reconstruction of the human myocyte metabolic network and its use for identification of markers for diabetes. Cell Rep.11, 921–933. 10.1016/j.celrep.2015.04.010
75
von KampA.KlamtS. (2014). Enumeration of smallest intervention strategies in genome-scale metabolic networks. PLoS Comput. Biol.10:e1003378. 10.1371/journal.pcbi.1003378
76
WibackS. J.MahadevanR.PalssonB. O. (2004). Using metabolic flux data to further constrain the metabolic solution space and predict internal flux patterns: the Escherichia coli spectrum. Biotechnol. Bioeng.86, 317–331. 10.1002/bit.20011
77
XuP.RanganathanS.FowlerZ. L.MaranasC. D.KoffasM. A. G. (2011). Genome-scale metabolic network modeling results in minimal interventions that cooperatively force carbon flux towards malonyl-CoA. Metab. Eng.13, 578–587. 10.1016/j.ymben.2011.06.008
78
YimH.HaselbeckR.NiuW.Pujol-Baxley CC.BurgardA.BoldtJ.et al. (2011). Metabolic engineering of Escherichia coli for direct production of 1, 4-butanediol. Nat. Chem. Biol. 7, 445–452. 10.1038/nchembio.580
79
YizhakK.ChanetonB.GottliebE.RuppinE. (2015). Modeling cancer metabolism on a genome-scale. Mol. Syst. Biol.11:817. 10.15252/msb.20145307
80
YizhakK.DévédecS. E. L.RogkotiV. M.BaenkeF.de BoerV. C.FrezzaC.et al. (2014). A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration. Mol. Syst. Biol.10:744. 10.15252/msb.20145746
81
ZhangC.JiB.MardinogluA.NielsenJ.HuaQ. (2015). Logical transformation of genome-scale metabolic models for gene level applications and analysis. Bioinformatics31, 2324–2331. 10.1093/bioinformatics/btv134
82
ZurH.RuppinE.ShlomiT. (2010). iMAT: an integrative metabolic analysis tool. Bioinformatics. 26, 3140–3142. 10.1093/bioinformatics/btq602
Summary
Keywords
genome-scale metabolic models, systems biology, metabolic capability analysis, in silico metabolic engineering, systems medicine
Citation
Zhang C and Hua Q (2016) Applications of Genome-Scale Metabolic Models in Biotechnology and Systems Medicine. Front. Physiol. 6:413. doi: 10.3389/fphys.2015.00413
Received
29 October 2015
Accepted
15 December 2015
Published
07 January 2016
Volume
6 - 2015
Edited by
Xiaogang Wu, Institute for Systems Biology, USA
Reviewed by
Stephen Fong, Virginia Commonwealth University, USA; Kazuyuki Shimizu, Keio University, Japan; Hiroshi Shimizu, Osaka University, Japan
Updates
Copyright
© 2016 Zhang and Hua.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Hua qhua@ecust.edu.cn
This article was submitted to Systems Biology, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.