Skip to main content

MINI REVIEW article

Front. Bioeng. Biotechnol., 25 January 2021
Sec. Synthetic Biology

Opportunities at the Interface of Network Science and Metabolic Modeling

  • 1Department of Life Sciences, Imperial College London, London, United Kingdom
  • 2Department of Mathematics, Imperial College London, London, United Kingdom
  • 3Department of Surgery and Cancer, Imperial College London, London, United Kingdom
  • 4Department of Metabolism, Digestion and Reproduction, Imperial College London, London, United Kingdom
  • 5School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
  • 6School of Informatics, University of Edinburgh, Edinburgh, United Kingdom

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.

1. Introduction

Metabolism comprises the biochemical reactions that convert nutrients into biomolecules and energy to sustain cellular functions. Advances in high-throughput screening technologies have enabled the quantitative characterization of metabolites, proteins and nucleic acids at the genome-scale, revealing previously unknown links between metabolism and many other cellular processes. For example, gene regulation (Chubukov et al., 2014), signal transduction (Tretter et al., 2016), immunity (Loftus and Finlay, 2016), and epigenetic modifications (Reid et al., 2017) have been shown to interact closely with metabolic processes. The increasing availability of data and the fundamental roles of metabolism in various cellular phenotypes (Tonn et al., 2019) have triggered a surge in metabolic research, together with a revived need for computational methods to untangle its complexity.

At the genome scale, metabolism comprises multiple interconnected reactions devoted to the production of energy and synthesis of essential biomolecules (e.g., proteins, lipids, or nucleic acids). The notion of a metabolic pathway is typically employed to organize sets of related reactions into functionally cohesive subsystems. Thus, lipid pathways, for example, are traditionally studied as distinct subsystems from amino acid or aerobic respiration pathways. Although conveniently descriptive, such a priori partitioning can obscure the links between other relevant layers of metabolic organization. Furthermore, metabolic connectivity is not static but actively responds and adapts to extracellular cues. Through various layers of transcriptional, translational, and post-translational regulation, metabolic pathways can be activated or shut down depending on external perturbations. These metabolic shifts drive a number of fundamental biological processes, such as microbial adaptations to growth conditions (Dai et al., 2016; Hartline et al., 2020) or the ability of pathogens to rewire their metabolism and evade the action of antimicrobial drugs (Olive and Sassetti, 2016). Metabolic adaptations are also thought to modulate the onset of complex diseases such as cancer (Hanahan and Weinberg, 2011; Pavlova and Thompson, 2016), diabetes, Alzheimer's, among others (DeBerardinis and Thompson, 2012; Suhre and Gieger, 2012). As a result, there is a growing need for computational tools that go beyond classical pathway definitions and can uncover hidden relations between metabolic components.

The complexity of metabolism has prompted the development of a myriad of methods to analyse its connectivity (Wishart et al., 2018). For specific pathways, kinetic models based on differential equations are widely employed to describe temporal dynamics of metabolites (Steuer et al., 2006; Saa and Nielsen, 2017). At the genome scale, however, the construction of kinetic models faces substantial challenges (Srinivasan et al., 2015). Such models require a large number of parameters, many of which have not been experimentally measured, or their values are subject to large uncertainty. As a result, the majority of genome-scale analyses are based on the metabolic stoichiometry alone. A widely adopted method for genome-scale modeling is Flux Balance Analysis (Palsson, 2015) (FBA), a powerful framework to predict metabolic fluxes on the basis of an optimization principle applied to the network stoichiometry. Alternatively, from the stoichiometry one can build graphs, a computational description of complex systems that has become the cornerstone of network science (Newman, 2010).

In this paper we discuss the relationship between FBA and graph-based analyses of metabolism, and we underline the complementary perspectives they bring to the understanding of metabolic organization. On the one hand, FBA has been shown to predict metabolic activity in various environmental and genetic contexts; on the other, network science can shed light on the emergent properties of global metabolic connectivity. Both approaches share a common root in the genome-scale stoichiometry of cellular metabolism, yet they offer different tools for its analysis. In the following, we discuss their advantages and caveats, highlighting the need and opportunities for integrated methods that combine flux optimization with network science.

2. Flux Balance Analysis

A large number of methods have been developed for the analysis of genome-scale metabolic networks (Lewis et al., 2012); these are generally described as constraint-based methods (Orth et al., 2010), an umbrella term for various techniques focused on the solution of the steady state equation:

Sv=0,    (1)

where S is the n × m stoichiometry matrix for a model with n metabolites and m reactions, and v is a vector containing the m reaction fluxes.

In general, Equation (1) is satisfied by an infinite number of flux vectors. A number of methods aim at probing the geometry of such flux solution space. For example, Elementary Flux Modes (Klamt et al., 2017) and Extreme Pathways (Wagner and Urbanczik, 2005) are two complementary techniques for decomposing the solution space into simpler units (Zanghellini et al., 2013; Muller and Bockmayr, 2014). Other methods for exploring the solution space include random flux sampling with Monte Carlo methods (Wiback et al., 2004), the use of dimensionality reduction techniques (Bhadra et al., 2018), and various structural decompositions of the stoichiometric matrix (Ghaderi et al., 2020).

The most widespread method for genome-scale modeling is Flux Balance Analysis (FBA), which selects a vector of metabolic fluxes v in Equation (1) as a solution to the optimization problem:

                       maxv J(v)subject to: Sv=0                        ViminviVimax, i=1,,m,    (2)

where (Vimin,Vimax) are bounds on each flux. The objective function J(v) is chosen to describe the physiology of a particular organism under study. In microbes, biomass production is the most common choice for the objective function, in which J(v) = cT·v, i.e., the rate of biomass production is assumed to be a linear combination of specific biosynthetic fluxes, defined by the positive vector of weights c. There are many dedicated FBA software packages (Lakshmanan et al., 2014; Lieven et al., 2020) and its popularity has led to a myriad of extensions (Lewis et al., 2012) that account for other complexities of cell physiology such as gene regulation (Covert et al., 2008), dynamic adaptations (Rügen et al., 2015; Waldherr et al., 2015), and many others (Heirendt et al., 2019).

Flux Balance Analysis has found applications in diverse domains, including cell biology (McCloskey et al., 2013), metabolic engineering (Nielsen and Keasling, 2016), microbiome studies (Khandelwal et al., 2013; Manor et al., 2014; Rosario et al., 2018), and personalized medicine (Diener and Resendis-Antonio, 2016; Nielsen, 2017; Raškevičius et al., 2018). A salient feature of FBA is its ability to incorporate various types of omics datasets into its predictions. Various approaches have been developed for this purpose (Becker and Palsson, 2008; Colijn et al., 2009; Jerby et al., 2010; Lee et al., 2012a; Wang et al., 2012; Agren et al., 2014; Nam et al., 2014; Yizhak et al., 2014), most of which incorporate experimental data into the metabolic model through adjustments of the stoichiometric matrix S or the flux bounds Vimin and Vimax in (2).

A popular use case of FBA is the identification of essential genes, i.e., genes that severely impact cellular growth when knocked out. Through simulation of gene deletions, FBA can serve as a systematic tool for in silico screening of lethal mutations, and identification of biomarkers and drug targets in disease (Lehár et al., 2009; Raman and Chandra, 2009; Folger et al., 2011; Gatto et al., 2015; Krueger et al., 2016; Pagliarini et al., 2016; Robinson and Nielsen, 2017). A related application of FBA is the study of metabolic robustness. Since only a fraction of all metabolic reactions are essential in a given environment, knocking out non-essential reactions often has little effect on the phenotype. This is because many reactions have functional backups through other pathways, so as to preserve cellular function in face of perturbations. By providing insights into the reorganization of fluxes under different conditions, FBA can also help improve our understanding of robustness to gene knockouts (Blank et al., 2005; Deutscher et al., 2006; Larhlimi et al., 2011; Palsson, 2015; Ho and Zhang, 2016), gene mutations (Fong and Palsson, 2004), and different growth conditions (Ibarra et al., 2002).

One limitation of FBA is the crucial importance of the objective function to be optimized, which needs to be designed to represent cellular physiology. In microbes, a common choice is maximization of growth rate, but it is questionable whether this is a realistic cellular objective across organisms or in different growth conditions (Schuetz et al., 2007; Feist and Palsson, 2010; Garćıa Sánchez and Torres Sáez, 2014). Although the vast majority of FBA studies rely on the maximization of cellular growth, other objective functions have been proposed, including maximization of ATP production (Nam et al., 2014) and minimization of substrate uptake rate (Raman and Chandra, 2009).

3. Network Science in Metabolic Modeling

Network science represents complex systems as graphs where the nodes are the system components and the edges describe interactions between them. This general description provides a framework for modeling large, interconnected systems across many disciplines, including biology, sociology, economics, and others (Newman, 2010). Numerous works have analyzed metabolism under the lens of network science. Graph-theoretic concepts such as degree distributions and centrality measures (Jeong et al., 2000; Wagner and Fell, 2001; Ma and Zeng, 2003b) can reveal structural features of the connectivity of the overall system, while clustering algorithms can uncover substructures hidden in the network topology. Such tools can be combined with the analysis of perturbations, such as deletions of network nodes or edges (Palumbo et al., 2005; Larhlimi et al., 2011), which can represent changes in the environment, gene knockouts, or therapeutic drugs that target specific metabolic enzymes. Unlike FBA, in which the analysis depends on the choice of a specific objective function, network methods rely on the metabolic stoichiometry alone.

Metabolic modularity is an area where network science has shown promising results. Intuitively, a network module is a subset of the network containing nodes that are more connected among themselves than to the rest of the network. Several studies have focused on the modularity of metabolism, and how network modules can be used to coarse-grain the metabolic connectivity into subunits (Ma and Zeng, 2003b; Tanaka et al., 2005; Zhao et al., 2007; da Silva et al., 2008; Kreimer et al., 2008). The modules identified with network analysis have been found to capture the organization of textbook biochemical pathways while uncovering novel links and relationships between them (Ravasz et al., 2002). A recurring theme in these analyses is the bow-tie topology, whereby a metabolic network can be divided into an input component, an output component and a strongly connected internal component. This architecture aligns well with an intuitive understanding of metabolism, which comprises nutrient uptake, waste production and secretion, and a large number of internal cycles which produce biomass and energy (Ma and Zeng, 2003b; Tanaka et al., 2005; Zhao et al., 2007; da Silva et al., 2008; Kreimer et al., 2008; Cooper and Barahona, 2010).

Despite its promise, however, network science has generally achieved mixed success in metabolic research. For example, from a network perspective it would be natural to expect that essential genes should be associated with high centrality scores (Jeong et al., 2001; Plaimas et al., 2010; Raman et al., 2014; Jalili et al., 2016). This idea draws parallels from other domains, such as the internet and social networks, where highly central nodes are deemed critical for network connectivity. However, correlations between gene essentiality and node centrality have been so far shown to be weak, with various essential metabolites and reactions exhibiting low centrality scores (Mahadevan and Palsson, 2005; Samal et al., 2006). This happens because poorly connected nodes are often the sole route for producing precursors that are essential for growth; in other words, such nodes lack a functional backup that can compensate for their loss. For example, Samal et al. (2006) showed that more than 50% of essential reactions in Escherichia coli, Saccharomyces cerevisiae, and Staphylococcus aureus are involved in such unique pathways, while other works noted that removal of poorly connected metabolites nodes can disrupt subsystems leading to failure of entire networks (Mahadevan and Palsson, 2005; Winterbach et al., 2011). Other studies have attempted to resolve this problem with new network metrics specifically tailored to describe important features of metabolism (Palumbo et al., 2005; Rahman and Schomburg, 2006; Wunderlich and Mirny, 2006; Cooper and Barahona, 2010; Kim et al., 2019; Yeganeh et al., 2020).

A key challenge for the use of network science in metabolic modeling is the lack of consensus on how to build a graph from a metabolic model. For a network with q nodes, the graph is encoded through the q × q adjacency matrix A, which has an entry Aij ≠ 0 if nodes i and j are connected, and Aij = 0 otherwise. As illustrated in Figure 1, depending on how nodes and edges are defined, one can build different graphs for the same metabolic model described by the stoichiometry matrix S in (2). From a metabolite-centric perspective one can build a graph where the nodes are metabolites and the edges corresponds to reactions between them (Ma and Zeng, 2003b; Asgari et al., 2013). In this case the adjacency matrix is

An×n=S^S^T,    (3)

where S^ is the binary version of the stoichiometry matrix S (i.e., Ŝij = 1 when Sij ≠ 0, and Ŝij = 0 otherwise). Conversely, from a reaction-centric perspective we can construct graphs with reactions as nodes and edges describing the sharing of metabolites as reactants or products (Ma et al., 2004; Beguerisse-D́ıaz et al., 2018). Such graph has an adjacency matrix

Am×m=S^TS^.    (4)

One can also build bipartite graphs, where both metabolites and reactions are nodes of different types (Holme, 2009; Beber et al., 2012), or even hypergraphs where an edge connects a set of reactants to a set of products (Cottret et al., 2010; Pearcy et al., 2016). In addition, all of these graphs can be directed/undirected (when the matrix A is asymmetric/symmetric), or weighted/unweighted (where the elements Aij can have weights encoding different properties). Such modeling choices can strongly influence the conclusions drawn from network analyses (Klamt et al., 2009; Bernal and Daza, 2011; Beguerisse-D́ıaz et al., 2018). For example, the existence of power law degree distributions (Jeong et al., 2000) and the small-world property in metabolism (Wagner and Fell, 2001), two cornerstone concepts in network science, have been disputed (Arita, 2004; Lima-Mendez and van Helden, 2009) and attributed to specific ways of constructing the network graph (Montañez et al., 2010; Bernal and Daza, 2011).

FIGURE 1
www.frontiersin.org

Figure 1. Graph constructions for metabolic networks. Depending on how nodes and edges are defined, several graphs can be built from a single metabolic model (Palsson, 2015). The conclusions drawn from graph analyses depend strongly on the choice of graph. The lack of consensus on the construction of such graphs is a key challenge for the use of network science in metabolic modeling.

A further limitation of graph-based analyses is their ad hoc treatment of pool metabolites, e.g., H2O, ATP, NADH, and other enzymatic co-factors. Because pool metabolites participate in a large number of reactions, they can distort and dominate the topological properties of reaction-centric graphs (Ma and Zeng, 2003b). A common approach to mitigate this problem is pruning the pool metabolites from the graph; yet there are no established best practices on how to choose which pool metabolites to prune, or how to mitigate the potential loss of information in so doing (Ma and Zeng, 2003a; Gerlee et al., 2009).

Another challenge arises from the reversibility of metabolic reactions in the graph. Although all biochemical reactions are reversible, they take one direction depending on the physiological conditions. The analysis of reaction-centric graphs typically prescribe a direction for reaction flux, or they split them into forward and backward components (Wagner and Fell, 2001; Helden et al., 2002). Neither of these approaches is ideal: assigning the direction of a reaction based on one condition may not generalize across other conditions, whereas incorporating bi-directional edges increases the complexity of the analysis.

4. Flux-Weighted Graphs: Integration of FBA and Network Science

As discussed in previous sections, both FBA and network science require modeling choices that can shape the conclusions drawn from their analyses. Tools from network science have already been employed to improve FBA pipelines in various ways (Lewis et al., 2012). Here we argue that the converse, i.e., using FBA to enrich the metabolic graphs, offers promising avenues to overcome some of their individual shortcomings. Flux information obtained from FBA solutions can be employed to assign direction and strength to the interactions between nodes in a graph. Such flux-weighted graphs allow to constrain their connectivity to various growth conditions, resulting in graphs that do not represent one universal network but are rather tailored to specific environmental or physiological contexts. As illustrated in Figure 2A, the integration of FBA and graph construction can thus result in a highly flexible pipeline to study metabolic connectivity in different functional states of an organism.

FIGURE 2
www.frontiersin.org

Figure 2. Integration of Flux Balance Analysis and graph theory. (A) Construction of flux-weighted graphs. Starting from a network graph, FBA solutions can be used to weight the graph edges, leading to different instances of the same graph for various perturbations such as gene knockouts, heterologous gene expression, or changes in growth conditions. (B) Clustering of flux-weighted graphs can reveal hidden groupings in the structure of metabolic networks (Yoon et al., 2007; Beguerisse-D́ıaz et al., 2018). (C) Essentiality and node centrality are two areas where flux-weighted graphs offer promising potential. (Left) New essentiality scores can be defined to rank metabolites according to their impact on the phenotype (Koschützki et al., 2010; Riemer et al., 2013; Laniau et al., 2017). (Right) Changes in node centrality between wild type (WT) and deletion (WT-) networks can reveal the molecular players associated to disease phenotypes; orange nodes denote reactions that undergo substantial changes in centrality upon gene deletions (Li et al., 2013; Beguerisse-D́ıaz et al., 2018).

Although the literature on this subject is still scarce, a number of studies have demonstrated the potential of the integration of FBA into graph analyses. These studies cover a range of methodologies and applications, including e.g., the identification of biomarkers (Li et al., 2013), detection of metabolic drug targets (Li et al., 2010), and quantification of metabolite essentiality (Riemer et al., 2013; Laniau et al., 2017). In one of the early works in the subject, Smart et al. (2008) proposed an adaptation of FBA that takes into account the connectivity of individual nodes. This idea revealed new insights on how the connectivity of specific metabolites provides robustness to metabolic networks.

Other studies have explored the use of FBA to construct flux-weighted graphs with either metabolites as nodes (Yoon et al., 2007; Koschützki et al., 2010; Riemer et al., 2013) or reactions as nodes (Kelk et al., 2012; Li et al., 2013; Beguerisse-D́ıaz et al., 2018). An alternative approach defined the concept of flux similarity (Li et al., 2010) to build reaction-drug graphs for detection of drug targets in cancer. Most recently, Hari and Lobo (2020) developed Fluxer, a web tool for visualization and analysis of flux-weighted metabolite graphs. The software allows the inclusion of customizable edge weights based on reaction fluxes and can perform multi-reaction knockout simulations.

In terms of applications, most studies have focused on flux-weighted graphs for the analysis of metabolic modularity and essentiality. Next we briefly discuss some of the approaches so far in these two application domains.

4.1. Network Clustering

A promising application of flux-weighted graphs is the detection of modular subunits within genome-scale metabolic models (Figure 2B). The idea is that flux-weighted graphs can encode information on the strengths on interactions between graph nodes that are specific to a particular physiological state, as modeled by the FBA solution. This can potentially reveal hidden groupings within metabolism, or how known groupings change across different contexts. For example, Yoon et al. (2007) employed experimentally determined fluxes to build flux-weighted graphs with metabolites as nodes. Using clustering algorithms on the graphs for energy metabolism of rat liver and adipose tissue formation, the approach revealed changes in cluster membership under different physiological flux distributions.

Another promising approach is the “mass flow graph” proposed by Beguerisse-D́ıaz et al. (2018), which uses FBA solutions to weight the edges of graph with reactions as nodes. In this approach, if reaction Ri produces a metabolite xk that is consumed by Rj, then the weight of the edge between both reactions is

wij=k(mass flow ofxk from Rito Rj),    (5)

where the sum acts on all the metabolites that are produced by Ri and consumed by Rj. The mass flows in (5) are directly computed from the stoichiometric matrix S and a flux vector obtained with FBA. Different mass flow graphs can be then computed for FBA solutions corresponding to specific environmental conditions. Thanks to the flux weighting, mass flow graphs avoid the need to prune pool metabolites, a common limitation of reaction graphs (Gerlee et al., 2009). Although pool metabolites do create many connections between functionally unrelated reactions, in mass flow graphs such connections are weak as a result of the flux weighting. This feature allowed the use of multiscale community detection algorithms to study changes in the modular structure of E. coli metabolism in various growth media (Beguerisse-D́ıaz et al., 2018).

4.2. Centrality and Essentiality

Flux-graph integration has also provided opportunities to explore centrality scores for quantifying essentiality of reactions and metabolites (Figure 2C). One example of this approach (Li et al., 2013) demonstrated that the combination of PageRank centrality (Newman, 2010) with flux information can help identifying candidate biomarker genes in disease. The use of flux-weighted graphs also allows to compare their connectivity between models that lack specific metabolic genes, e.g., in the case of mutants or genetic deficiencies found in metabolic disorders. For example, PageRank centrality was employed in conjunction with mass flow graphs (Beguerisse-D́ıaz et al., 2018) to study structural changes in hepatocyte metabolism in primary hyperoxaluria type 1, a rare metabolic disease characterized by the lack of the agt gene involved in glyoxylate breakdown (Pagliarini et al., 2016). This approach showed that reactions which underwent the highest PageRank changes between healthy and diseased states were directly related to the PH1 phenotype (Figure 2C). Importantly, some of the changes in PageRank centrality did not correlate with changes in flux, providing strong evidence that metabolic graphs can encode information that cannot inferred from FBA alone.

A number of other works have sought to define new, metabolism-specific, centrality scores that can reveal new information on the topology of metabolic networks. For example, Koschützki et al. (2010) built a novel “flux centrality” score for metabolites in networks where only the carbon exchanges are modeled as edges. This metric emphasizes the role that a metabolite plays in biomass formation based on both topology and flux, penalizing the impact of highly connected pool metabolites. Riemer et al. (2013) combined the classic notion of metabolic branch points, i.e., metabolites that are substrates to multiple downstream pathways, with reaction fluxes so as to rank metabolites according to various metrics of essentiality. A similar approach to establish metabolite essentiality was presented by Laniau et al. (2017), where they classify metabolites on the basis of their capacity to influence the activation of a target objective function.

5. Discussion

Recent discoveries have led to a renewed interest in the interplay of metabolism with other layers of the cellular machinery (Chubukov et al., 2014; Loftus and Finlay, 2016; Tretter et al., 2016; Reid et al., 2017; Tonn et al., 2020). Due to the complexity and scale of metabolic reaction networks, computational methods are essential to tease apart the influence of metabolic architectures on cellular function. Here we have discussed the complementary roles of Flux Balance analysis and network science in the analysis of metabolism at the genome scale. Although both approaches start from the metabolic stoichiometry, they differ in their mathematical foundations and the type of predictions they produce. FBA predictions can be accurate but their effectiveness requires high quality omics datasets. Network science, in contrast, requires nothing more than the metabolic stoichiometry, yet can lead to misleading predictions depending on how the network graph is built. As a result, so far FBA has led to more successful connections with experimental results than network science.

When used in isolation, both FBA and network science can be insufficient to understand changes in metabolic connectivity triggered by physiological or environmental perturbations. Here we argue that the use of flux-weighted graphs (Figure 2A) allows for a natural integration of FBA and network science, applicable in many subject domains. For example, with the rise of big data in the life sciences, there is a growing interest in using patient metabolic signatures to tailor treatments (O'Day et al., 2018). Computational methods can play a key role in detecting drug targets involved in metabolic activity, and how their targeting can disrupt metabolic connectivity. A particularly promising area is cancer treatment, where there is considerable interest on drugs that target specific metabolic enzymes (Neradil et al., 2012; Nishi et al., 2016). Moreover, novel data-driven approaches based on machine learning can also be integrated with FBA (Zampieri et al., 2019; Kavvas et al., 2020) and network science to extend their capabilities into novel applications.

Another promising application domain is industrial biotechnology (de Lorenzo et al., 2018), where microbial cell factories are engineered for production of commodity chemicals and fine products (Lee et al., 2012b). In this field, FBA is widely employed for strain design, with the goal of finding combinations of genetic interventions that maximize production of a desired metabolite. A recent trend is to increase production with synthetic biology tools and dynamic control of gene expression (Brockman and Prather, 2015; Liu et al., 2018). This approach needs computational methods that capture the dynamic reallocation of metabolic flux (Oyarzún and Stan, 2013). Integrating FBA solutions with network models can provide a versatile tool to identify suitable genetic modifications for microbial strains with increased production.

Further developments at the interface of FBA and network science offer a novel way to explore the impact of perturbations on metabolic connectivity. The flexibility of FBA allows for the modeling of metabolic perturbations of various kinds, including changes in growth conditions, deletion of metabolic genes or the action of enzyme inhibitors, whereas the application of network-theoretical tools can bring a broadened understanding of emergent properties of the overall system. This flexibility offers promising potential to deploy network science tools across a range of questions in basic science, biomedicine, and industrial biotechnology.

Author Contributions

VD and DO conceived the study. All authors contributed to writing the manuscript.

Funding

This work was supported by Cancer Research UK (C24523/A27435), the Cancer Research UK Imperial Centre, and the EPSRC Centre for Mathematics of Precision Healthcare (EP/N014529/1).

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.

The handling editor declared a past co-authorship with one of the author, DO.

References

Agren, R., Mardinoglu, A., Asplund, A., Kampf, C., Uhlen, M., and Nielsen, J. (2014). Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling. Mol. Syst. Biol. 10:721. doi: 10.1002/msb.145122

CrossRef Full Text | Google Scholar

Arita, M. (2004). The metabolic world of Escherichia coli is not small. Proc. Natl. Acad. Sci. U. S. A. 101, 1543–1547. doi: 10.1073/pnas.0306458101

PubMed Abstract | CrossRef Full Text | Google Scholar

Asgari, Y., Salehzadeh-Yazdi, A., Schreiber, F., and Masoudi-Nejad, A. (2013). Controllability in cancer metabolic networks according to drug targets as driver nodes. PLoS ONE 8:e79397. doi: 10.1371/journal.pone.0079397

PubMed Abstract | CrossRef Full Text | Google Scholar

Beber, M. E., Fretter, C., Jain, S., Sonnenschein, N., Muller-Hannemann, M., and Hutt, M. T. (2012). Artefacts in statistical analyses of network motifs: general framework and application to metabolic networks. J. R. Soc. Interface 9, 3426–3435. doi: 10.1098/rsif.2012.0490

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, S. A., and Palsson, B. O. (2008). Context-specific metabolic networks are consistent with experiments. PLoS Comput. Biol. 4:e1000082. doi: 10.1371/journal.pcbi.1000082

PubMed Abstract | CrossRef Full Text | Google Scholar

Beguerisse-Díaz, M., Bosque, G., Oyarzún, D., Picó, J., and Barahona, M. (2018). Flux-dependent graphs for metabolic networks. NPJ Syst. Biol. Appl. 4:32. doi: 10.1038/s41540-018-0067-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernal, A., and Daza, E. (2011). Metabolic networks: beyond the graph. Curr. Comput. Aided Drug Des. 7, 122–132. doi: 10.2174/157340911795677611

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhadra, S., Blomberg, P., Castillo, S., and Rousu, J. (2018). Principal metabolic flux mode analysis. Bioinformatics 34, 2409–2417. doi: 10.1093/bioinformatics/bty049

PubMed Abstract | CrossRef Full Text | Google Scholar

Blank, L. M., Kuepfer, L., and Sauer, U. (2005). Large-scale 13c-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast. Genome Biol. 6:R49. doi: 10.1186/gb-2005-6-6-r49

PubMed Abstract | CrossRef Full Text | Google Scholar

Brockman, I. M., and Prather, K. L. (2015). Dynamic metabolic engineering: new strategies for developing responsive cell factories. Biotechnol. J. 10, 1360–1369. doi: 10.1002/biot.201400422

PubMed Abstract | CrossRef Full Text | Google Scholar

Chubukov, V., Gerosa, L., Kochanowski, K., and Sauer, U. (2014). Coordination of microbial metabolism. Nat. Rev. Microbiology 12, 327–340. doi: 10.1038/nrmicro3238

PubMed Abstract | CrossRef Full Text | Google Scholar

Colijn, C., Brandes, A., Zucker, J., Lun, D. S., Weiner, B., Farhat, M. R., et al. (2009). Interpreting expression data with metabolic flux models: predicting mycobacterium tuberculosis mycolic acid production. PLoS Comput. Biol. 5:e1000489. doi: 10.1371/journal.pcbi.1000489

PubMed Abstract | CrossRef Full Text | Google Scholar

Cooper, K., and Barahona, M. (2010). Role-based similarity in directed networks. arXiv. arXiv:1012.2726.

Google Scholar

Cottret, L., Milreu, P. V., Acuna, V., Marchetti-Spaccamela, A., Stougie, L., Charles, H., et al. (2010). Graph-based analysis of the metabolic exchanges between two co-resident intracellular symbionts, Baumannia cicadellinicola and Sulcia muelleri, with their insect host, Homalodisca coagulata. PLoS Comput. Biol. 6:e1000904. doi: 10.1371/journal.pcbi.1000904

PubMed Abstract | CrossRef Full Text | Google Scholar

Covert, M. W., Xiao, N., Chen, T. J., and Karr, J. R. (2008). Integrating metabolic, transcriptional regulatory and signal transduction models in Escherichia coli. Bioinformatics 24, 2044–2050. doi: 10.1093/bioinformatics/btn352

CrossRef Full Text | Google Scholar

da Silva, M. R., and Hongwu Ma, An-Ping Zeng (2008). Centrality, network capacity, and modularity as parameters to analyze the Core-Periphery structure in metabolic networks. Proc. IEEE 96, 1411–1420. doi: 10.1109/JPROC.2008.925418

CrossRef Full Text | Google Scholar

Dai, X., Yan, H., Li, N., He, J., Ding, Y., Dai, L., et al. (2016). Metabolic adaptation of microbial communities to ammonium stress in a high solid anaerobic digester with dewatered sludge. Sci. Rep. 6:28193. doi: 10.1038/srep28193

PubMed Abstract | CrossRef Full Text | Google Scholar

de Lorenzo, V., Prather, K. L., Chen, G., O'Day, E., von Kameke, C., Oyarzún, D. A., et al. (2018). The power of synthetic biology for bioproduction, remediation and pollution control. EMBO Rep. 19:e45658. doi: 10.15252/embr.201745658

PubMed Abstract | CrossRef Full Text | Google Scholar

DeBerardinis, R. J., and Thompson, C. B. (2012). Cellular metabolism and disease: what do metabolic outliers teach us? Cell 148, 1132–1144. doi: 10.1016/j.cell.2012.02.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Deutscher, D., Meilijson, I., Kupiec, M., and Ruppin, E. (2006). Multiple knockout analysis of genetic robustness in the yeast metabolic network. Nat. Genet. 38, 993–998. doi: 10.1038/ng1856

PubMed Abstract | CrossRef Full Text | Google Scholar

Diener, C., and Resendis-Antonio, O. (2016). Personalized prediction of proliferation rates and metabolic liabilities in cancer biopsies. Front. Physiol. 7:644. doi: 10.3389/fphys.2016.00644

PubMed Abstract | CrossRef Full Text | Google Scholar

Feist, A. M., and Palsson, B. O. (2010). The biomass objective function. Curr. Opin. Microbiol. 13, 344–349. doi: 10.1016/j.mib.2010.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Folger, O., Jerby, L., Frezza, C., Gottlieb, E., Ruppin, E., and Shlomi, T. (2011). Predicting selective drug targets in cancer through metabolic networks. Mol. Syst. Biol. 7:501. doi: 10.1038/msb.2011.35

PubMed Abstract | CrossRef Full Text | Google Scholar

Fong, S. S., and Palsson, B. Ø. (2004). Metabolic gene-deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes. Nat. Genet. 36, 1056–1058. doi: 10.1038/ng1432

PubMed Abstract | CrossRef Full Text | Google Scholar

García Sánchez, C. E., and Torres Sáez, R. G. (2014). Comparison and analysis of objective functions in flux balance analysis. Biotechnol. Prog. 30, 985–991. doi: 10.1002/btpr.1949

PubMed Abstract | CrossRef Full Text | Google Scholar

Gatto, F., Miess, H., Schulze, A., and Nielsen, J. (2015). Flux balance analysis predicts essential genes in clear cell renal cell carcinoma metabolism. Sci. Rep. 5:10738. doi: 10.1038/srep10738

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerlee, P., Lizana, L., and Sneppen, K. (2009). Pathway identification by network pruning in the metabolic network of Escherichia coli. Bioinformatics 25, 3282–3288. doi: 10.1093/bioinformatics/btp575

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghaderi, S., Haraldsdóttir, H. S., Ahookhosh, M., Arreckx, S., and Fleming, R. M. (2020). Structural conserved moiety splitting of a stoichiometric matrix. J. Theor. Biol. 499:110276. doi: 10.1016/j.jtbi.2020.110276

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. Cell 144, 646–674. doi: 10.1016/j.cell.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hari, A., and Lobo, D. (2020). Fluxer: a web application to compute, analyze and visualize genome-scale metabolic flux networks. Nucleic Acids Res. 48, W427–W435. doi: 10.1093/nar/gkaa409

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartline, C. J., Mannan, A. A., Liu, D., Zhang, F., and Oyarzún, D. A. (2020). Metabolite sequestration enables rapid recovery from fatty acid depletion in Escherichia coli. mBio 11:e03112-19. doi: 10.1128/mBio.03112-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Heirendt, L., Arreckx, S., Pfau, T., Mendoza, S. N., Richelle, A., Heinken, A., et al. (2019). Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nat. Protoc. 14, 639–702. doi: 10.1038/s41596-018-0098-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Helden, J., Wernisch, L., Gilbert, D., and Wodak, S. J. (2002). “Graph-Based analysis of metabolic networks,” in Bioinformatics and Genome Analysis, eds H.-W. Mewes, H. Seidel, and B. Weiss (Berlin; Heidelberg: Springer Berlin Heidelberg), 245–274. doi: 10.1007/978-3-662-04747-7_12

CrossRef Full Text | Google Scholar

Ho, W.-C., and Zhang, J. (2016). Adaptive genetic robustness of Escherichia coli metabolic fluxes. Mol. Biol. Evol. 33, 1164–1176. doi: 10.1093/molbev/msw002

PubMed Abstract | CrossRef Full Text | Google Scholar

Holme, P. (2009). Model validation of simple-graph representations of metabolism. J. R. Soc. Interface 6, 1027–1034. doi: 10.1098/rsif.2008.0489

PubMed Abstract | CrossRef Full Text | Google Scholar

Ibarra, R. U., Edwards, J. S., and Palsson, B. O. (2002). Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth. Nature 420, 186–189. doi: 10.1038/nature01149

PubMed Abstract | CrossRef Full Text | Google Scholar

Jalili, M., Salehzadeh-Yazdi, A., Gupta, S., Wolkenhauer, O., Yaghmaie, M., Resendis-Antonio, O., et al. (2016). Evolution of centrality measurements for the detection of essential proteins in biological networks. Front. Physiol. 7:375. doi: 10.3389/fphys.2016.00375

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeong, H., Mason, S. P., Barabási, A. L., and Oltvai, Z. N. (2001). Lethality and centrality in protein networks. Nature 411, 41–42. doi: 10.1038/35075138

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeong, H., Tombor, B., Albert, R., Oltvai, Z. N., and Barabási, A. L. (2000). The large-scale organization of metabolic networks. Nature 407, 651–654. doi: 10.1038/35036627

PubMed Abstract | CrossRef Full Text | Google Scholar

Jerby, L., Shlomi, T., and Ruppin, E. (2010). Computational reconstruction of tissue-specific metabolic models: application to human liver metabolism. Mol. Syst. Biol. 6:401. doi: 10.1038/msb.2010.56

PubMed Abstract | CrossRef Full Text | Google Scholar

Kavvas, E. S., Yang, L., Monk, J. M., Heckmann, D., and Palsson, B. O. (2020). A biochemically-interpretable machine learning classifier for microbial GWAS. Nat. Commun. 11:2580. doi: 10.1038/s41467-020-16310-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelk, S. M., Olivier, B. G., Stougie, L., and Bruggeman, F. J. (2012). Optimal flux spaces of genome-scale stoichiometric models are determined by a few subnetworks. Sci. Rep. 2:580. doi: 10.1038/srep00580

PubMed Abstract | CrossRef Full Text | Google Scholar

Khandelwal, R. A., Olivier, B. G., Roling, W. F., Teusink, B., and Bruggeman, F. J. (2013). Community flux balance analysis for microbial consortia at balanced growth. PLoS ONE 8:e64567. doi: 10.1371/journal.pone.0064567

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, E.-Y., Ashlock, D., and Yoon, S. H. (2019). Identification of critical connectors in the directed reaction-centric graphs of microbial metabolic networks. BMC Bioinformatics 20:328. doi: 10.1186/s12859-019-2897-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Klamt, S., Haus, U.-U., and Theis, F. (2009). Hypergraphs and cellular networks. PLoS Comput. Biol. 5:e1000385. doi: 10.1371/journal.pcbi.1000385

PubMed Abstract | CrossRef Full Text | Google Scholar

Klamt, S., Regensburger, G., Gerstl, M. P., Jungreuthmayer, C., Schuster, S., Mahadevan, R., et al. (2017). From elementary flux modes to elementary flux vectors: Metabolic pathway analysis with arbitrary linear flux constraints. PLoS Comput. Biol. 13:e1005409. doi: 10.1371/journal.pcbi.1005409

PubMed Abstract | CrossRef Full Text | Google Scholar

Koschützki, D., Junker, B. H., Schwender, J., and Schreiber, F. (2010). Structural analysis of metabolic networks based on flux centrality. J. Theor. Biol. 265, 261–269. doi: 10.1016/j.jtbi.2010.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreimer, A., Borenstein, E., Gophna, U., and Ruppin, E. (2008). The evolution of modularity in bacterial metabolic networks. Proc. Natl. Acad. Sci. U. S. A. 105, 6976–6981. doi: 10.1073/pnas.0712149105

PubMed Abstract | CrossRef Full Text | Google Scholar

Krueger, A. S., Munck, C., Dantas, G., Church, G. M., Galagan, J., Lehár, J., et al. (2016). Simulating Serial-Target antibacterial drug synergies using flux balance analysis. PLoS ONE 11:e0147651. doi: 10.1371/journal.pone.0147651

PubMed Abstract | CrossRef Full Text | Google Scholar

Lakshmanan, M., Koh, G., Chung, B. K. S., and Lee, D.-Y. (2014). Software applications for flux balance analysis. Brief. Bioinformatics 15, 108–122. doi: 10.1093/bib/bbs069

PubMed Abstract | CrossRef Full Text | Google Scholar

Laniau, J., Frioux, C., Nicolas, J., Baroukh, C., Cortes, M.-P., Got, J., et al. (2017). Combining graph and flux-based structures to decipher phenotypic essential metabolites within metabolic networks. PeerJ 5:e3860. doi: 10.7717/peerj.3860

PubMed Abstract | CrossRef Full Text | Google Scholar

Larhlimi, A., Blachon, S., Selbig, J., and Nikoloski, Z. (2011). Robustness of metabolic networks: a review of existing definitions. Biosystems 106, 1–8. doi: 10.1016/j.biosystems.2011.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, D., Smallbone, K., Dunn, W. B., Murabito, E., Winder, C. L., Kell, D. B., et al. (2012a). Improving metabolic flux predictions using absolute gene expression data. BMC Syst. Biol. 6:73. doi: 10.1186/1752-0509-6-73

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. W., Na, D., Park, J. M., Lee, J., Choi, S., and Lee, S. Y. (2012b). Systems metabolic engineering of microorganisms for natural and non-natural chemicals. Nat. Chem. Biol. 8, 536–546. doi: 10.1038/nchembio.970

CrossRef Full Text | Google Scholar

Lehár, J., Krueger, A. S., Avery, W., Heilbut, A. M., Johansen, L. M., Price, E. R., et al. (2009). Synergistic drug combinations tend to improve therapeutically relevant selectivity. Nat. Biotechnol. 27, 659–666. doi: 10.1038/nbt.1549

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, N., Nagarajan, H., and Palsson, B. Ø. (2012). Constraining the metabolic genotype-phenotype relationship using a phylogeny of in silico methods. Nat. Rev. Microbiol. 10, 291–305. doi: 10.1038/nrmicro2737

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Jiang, H., Qiu, Y., Ching, W. K., and Vassiliadis, V. S. (2013). Discovery of metabolite biomarkers: flux analysis and reaction-reaction network approach. BMC Syst. Biol. 7(Suppl. 2):S13. doi: 10.1186/1752-0509-7-S2-S13

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Zhou, X., Ching, W.-K., and Wang, P. (2010). Predicting enzyme targets for cancer drugs by profiling human metabolic reactions in NCI-60 cell lines. BMC Bioinformatics 11:501. doi: 10.1186/1471-2105-11-501

PubMed Abstract | CrossRef Full Text | Google Scholar

Lieven, C., Beber, M. E., Olivier, B. G., Bergmann, F. T., Ataman, M., Babaei, P., et al. (2020). MEMOTE for standardized genome-scale metabolic model testing. Nat. Biotechnol. 38, 272–276. doi: 10.1038/s41587-020-0446-y

CrossRef Full Text | Google Scholar

Lima-Mendez, G., and van Helden, J. (2009). The powerful law of the power law and other myths in network biology. Mol. Biosyst. 5, 1482–1493. doi: 10.1039/b908681a

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D., Mannan, A. A., Han, Y., Oyarzún, D. A., and Zhang, F. (2018). Dynamic metabolic control: towards precision engineering of metabolism. J. Indus. Microbiol. Biotechnol. 45, 535–543. doi: 10.1007/s10295-018-2013-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Loftus, R. M., and Finlay, D. K. (2016). Immunometabolism: cellular metabolism turns immune regulator. J. Biol. Chem. 291, 1–10. doi: 10.1074/jbc.R115.693903

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, H., and Zeng, A.-P. (2003a). Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms. Bioinformatics 19, 270–277. doi: 10.1093/bioinformatics/19.2.270

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, H.-W., and Zeng, A.-P. (2003b). The connectivity structure, giant strong component and centrality of metabolic networks. Bioinformatics 19, 1423–1430. doi: 10.1093/bioinformatics/btg177

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, H.-W., Zhao, X.-M., Yuan, Y.-J., and Zeng, A.-P. (2004). Decomposition of metabolic network into functional modules based on the global connectivity structure of reaction graph. Bioinformatics 20, 1870–1876. doi: 10.1093/bioinformatics/bth167

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahadevan, R., and Palsson, B. O. (2005). Properties of metabolic networks: structure versus function. Biophys. J. 88, L07–L09. doi: 10.1529/biophysj.104.055723

PubMed Abstract | CrossRef Full Text | Google Scholar

Manor, O., Levy, R., and Borenstein, E. (2014). Mapping the inner workings of the microbiome: genomic- and metagenomic-based study of metabolism and metabolic interactions in the human microbiome. Cell Metab. 20, 742–752. doi: 10.1016/j.cmet.2014.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

McCloskey, D., Palsson, B., and Feist, A. M. (2013). Basic and applied uses of genome-scale metabolic network reconstructions of Escherichia coli. Mol. Syst. Biol. 9:661. doi: 10.1038/msb.2013.18

PubMed Abstract | CrossRef Full Text | Google Scholar

Monta nez, R., Medina, M. A., Solé, R. V., and Rodríguez-Caso, C. (2010). When metabolism meets topology: reconciling metabolite and reaction networks. Bioessays 32, 246–256. doi: 10.1002/bies.200900145

CrossRef Full Text | Google Scholar

Muller, A. C., and Bockmayr, A. (2014). Flux modules in metabolic networks. J. Math. Biol. 69, 1151–1179. doi: 10.1007/s00285-013-0731-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Nam, H., Campodonico, M., Bordbar, A., Hyduke, D. R., Kim, S., Zielinski, D. C., et al. (2014). A systems approach to predict oncometabolites via context-specific genome-scale metabolic networks. PLoS Comput. Biol. 10:e1003837. doi: 10.1371/journal.pcbi.1003837

PubMed Abstract | CrossRef Full Text | Google Scholar

Neradil, J., Pavlasova, G., and Veselska, R. (2012). New mechanisms for an old drug: DHFR- and non-DHFR-mediated effects of methotrexate in cancer cells. Klin Onkol. 25(Suppl 2), 2S87–2S92.

PubMed Abstract | Google Scholar

Newman, M. (2010). Networks: An Introduction. Oxford University Press. doi: 10.1093/acprof:oso/9780199206650.003.0001

CrossRef Full Text | Google Scholar

Nielsen, J. (2017). Systems biology of metabolism: a driver for developing personalized and precision medicine. Cell Metab. 25, 572–579. doi: 10.1016/j.cmet.2017.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, J., and Keasling, J. D. (2016). Engineering cellular metabolism. Cell 164, 1185–1197. doi: 10.1016/j.cell.2016.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishi, K., Suzuki, K., Sawamoto, J., Tokizawa, Y., Iwase, Y., Yumita, N., et al. (2016). Inhibition of fatty acid synthesis induces apoptosis of human pancreatic cancer cells. Anticancer Res. 36, 4655–4660. doi: 10.21873/anticanres.11016

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Day, E., Hosta-Rigau, L., Oyarzún, D. A., Okano, H., de Lorenzo, V., von Kameke, C., et al. (2018). Are we there yet? How and when specific biotechnologies will improve human health. Biotechnol. J. 14:1800195. doi: 10.1002/biot.201800195

PubMed Abstract | CrossRef Full Text | Google Scholar

Olive, A. J., and Sassetti, C. M. (2016). Metabolic crosstalk between host and pathogen: sensing, adapting and competing. Nat. Rev. Microbiol. 14, 221–234. doi: 10.1038/nrmicro.2016.12

PubMed Abstract | CrossRef Full Text | Google Scholar

Orth, J. D., Thiele, I., and Palsson, B. Ø. (2010). What is flux balance analysis? Nat. Biotechnol. 28, 245–248. doi: 10.1038/nbt.1614

PubMed Abstract | CrossRef Full Text | Google Scholar

Oyarzún, D. A., and Stan, G.-B. V. (2013). Synthetic gene circuits for metabolic control: design trade-offs and constraints. J. R. Society Interface 10:671. doi: 10.1098/rsif.2012.0671

PubMed Abstract | CrossRef Full Text | Google Scholar

Pagliarini, R., Castello, R., Napolitano, F., Borzone, R., Annunziata, P., Mandrile, G., et al. (2016). In Silico modeling of liver metabolism in a human disease reveals a key enzyme for histidine and histamine homeostasis. Cell Rep. 15, 2292–2300. doi: 10.1016/j.celrep.2016.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Palsson, B. Ø. (2015). Systems Biology: Constraint-based Reconstruction and Analysis. Cambridge: Cambridge University Press. doi: 10.1017/CBO9781139854610

CrossRef Full Text | Google Scholar

Palumbo, M. C., Colosimo, A., Giuliani, A., and Farina, L. (2005). Functional essentiality from topology features in metabolic networks: a case study in yeast. FEBS Lett. 579, 4642–4646. doi: 10.1016/j.febslet.2005.07.033

CrossRef Full Text | Google Scholar

Pavlova, N. N., and Thompson, C. B. (2016). The emerging hallmarks of cancer metabolism. Cell Metab. 23, 27–47. doi: 10.1016/j.cmet.2015.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearcy, N., Chuzhanova, N., and Crofts, J. J. (2016). Complexity and robustness in hypernetwork models of metabolism. J. Theor. Biol. 406, 99–104. doi: 10.1016/j.jtbi.2016.06.032

CrossRef Full Text | Google Scholar

Plaimas, K., Eils, R., and Konig, R. (2010). Identifying essential genes in bacterial metabolic networks with machine learning methods. BMC Syst. Biol. 4:56. doi: 10.1186/1752-0509-4-56

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, S. A., and Schomburg, D. (2006). Observing local and global properties of metabolic pathways: 'load points' and 'choke points' in the metabolic networks. Bioinformatics 22, 1767–1774. doi: 10.1093/bioinformatics/btl181

PubMed Abstract | CrossRef Full Text | Google Scholar

Raman, K., and Chandra, N. (2009). Flux balance analysis of biological systems: applications and challenges. Brief. Bioinform. 10, 435–449. doi: 10.1093/bib/bbp011

PubMed Abstract | CrossRef Full Text | Google Scholar

Raman, K., Damaraju, N., and Joshi, G. K. (2014). The organisational structure of protein networks: revisiting the centrality-lethality hypothesis. Syst. Synth. Biol. 8, 73–81. doi: 10.1007/s11693-013-9123-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Raškevičius, V., Mikalayeva, V., Antanavičiūtė, I., Ceslevičienė, I., Skeberdis, V. A., Kairys, V., et al. (2018). Genome scale metabolic models as tools for drug design and personalized medicine. PLoS ONE 13:e0190636. doi: 10.1371/journal.pone.0190636

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N., and Barabási, A. L. (2002). Hierarchical organization of modularity in metabolic networks. Science 297, 1551–1555. doi: 10.1126/science.1073374

PubMed Abstract | CrossRef Full Text | Google Scholar

Reid, M. A., Dai, Z., and Locasale, J. W. (2017). The impact of cellular metabolism on chromatin dynamics and epigenetics. Nat. Cell Biol. 19, 1298–1306. doi: 10.1038/ncb3629

PubMed Abstract | CrossRef Full Text | Google Scholar

Riemer, S. A., Rex, R., and Schomburg, D. (2013). A metabolite-centric view on flux distributions in genome-scale metabolic models. BMC Syst. Biol. 7:33. doi: 10.1186/1752-0509-7-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, J. L., and Nielsen, J. (2017). Anticancer drug discovery through genome-scale metabolic modeling. Curr. Opin. Syst. Biol. 4, 1–8. doi: 10.1016/j.coisb.2017.05.007

CrossRef Full Text | Google Scholar

Rosario, D., Benfeitas, R., Bidkhori, G., Zhang, C., Uhlen, M., Shoaie, S., et al. (2018). Understanding the representative gut microbiota dysbiosis in metformin-treated type 2 diabetes patients using genome-scale metabolic modeling. Front. Physiol. 9:775. doi: 10.3389/fphys.2018.00775

PubMed Abstract | CrossRef Full Text | Google Scholar

Rügen, M., Bockmayr, A., and Steuer, R. (2015). Elucidating temporal resource allocation and diurnal dynamics in phototrophic metabolism using conditional FBA. Sci. Rep. 5:15247. doi: 10.1038/srep15247

PubMed Abstract | CrossRef Full Text | Google Scholar

Saa, P. A., and Nielsen, L. K. (2017). Formulation, construction and analysis of kinetic models of metabolism: a review of modelling frameworks. Biotechnol. Adv. 35, 981–1003. doi: 10.1016/j.biotechadv.2017.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Samal, A., Singh, S., Giri, V., Krishna, S., Raghuram, N., and Jain, S. (2006). Low degree metabolites explain essential reactions and enhance modularity in biological networks. BMC Bioinformatics 7:118. doi: 10.1186/1471-2105-7-118

PubMed Abstract | CrossRef Full Text | Google Scholar

Schuetz, R., Kuepfer, L., and Sauer, U. (2007). Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol. Syst. Biol. 3:119. doi: 10.1038/msb4100162

PubMed Abstract | CrossRef Full Text | Google Scholar

Smart, A. G., Amaral, L. A. N., and Ottino, J. M. (2008). Cascading failure and robustness in metabolic networks. Proc. Natl. Acad. Sci. U. S. A. 105, 13223–13228. doi: 10.1073/pnas.0803571105

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivasan, S., Cluett, W. R., and Mahadevan, R. (2015). Constructing kinetic models of metabolism at genome-scales: a review. Biotechnol. J. 10, 1345–1359. doi: 10.1002/biot.201400522

PubMed Abstract | CrossRef Full Text | Google Scholar

Steuer, R., Gross, T., Selbig, J., and Blasius, B. (2006). Structural kinetic modeling of metabolic networks. Proc. Natl. Acad. Sci. U. S. A. 103, 11868–11873. doi: 10.1073/pnas.0600013103

PubMed Abstract | CrossRef Full Text | Google Scholar

Suhre, K., and Gieger, C. (2012). Genetic variation in metabolic phenotypes: study designs and applications. Nat. Rev. Genet. 13, 759–769. doi: 10.1038/nrg3314

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, R., Csete, M., and Doyle, J. (2005). Highly optimised global organisation of metabolic networks. Syst. Biol. 152, 179–184. doi: 10.1049/ip-syb:20050042

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonn, M. K., Thomas, P., Barahona, M., and Oyarzún, D. A. (2019). Stochastic modelling reveals mechanisms of metabolic heterogeneity. Commun. Biol. 2:108. doi: 10.1038/s42003-019-0347-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonn, M. K., Thomas, P., Barahona, M., and Oyarzún, D. A. (2020). Computation of single-cell metabolite distributions using mixture models. Cell Dev. Biol. 8, 1–11. doi: 10.3389/fcell.2020.614832

PubMed Abstract | CrossRef Full Text | Google Scholar

Tretter, L., Patocs, A., and Chinopoulos, C. (2016). Succinate, an intermediate in metabolism, signal transduction, ROS, hypoxia, and tumorigenesis. Biochim. Biophys. Acta 1857, 1086–1101. doi: 10.1016/j.bbabio.2016.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, A., and Fell, D. A. (2001). The small world inside large metabolic networks. Proc. Biol. Sci. 268, 1803–1810. doi: 10.1098/rspb.2001.1711

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, C., and Urbanczik, R. (2005). The geometry of the flux cone of a metabolic network. Biophys. J. 89, 3837–3845. doi: 10.1529/biophysj.104.055129

PubMed Abstract | CrossRef Full Text | Google Scholar

Waldherr, S., Oyarzún, D. A., and Bockmayr, A. (2015). Dynamic optimization of metabolic networks coupled with gene expression. J. Theor. Biol. 365, 469–485. doi: 10.1016/j.jtbi.2014.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Eddy, J. A., and Price, N. D. (2012). Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE. BMC Syst. Biol. 6:153. doi: 10.1186/1752-0509-6-153

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiback, S. J., Famili, I., Greenberg, H. J., and Palsson, B. O. (2004). Monte carlo sampling can be used to determine the size and shape of the steady-state flux space. J. Theor. Biol. 228, 437–447. doi: 10.1016/j.jtbi.2004.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Winterbach, W., Wang, H., Reinders, M., Van Mieghem, P., and de Ridder, D. (2011). Metabolic network destruction: relating topology to robustness. Nano Commun. Netw. 2, 88–98. doi: 10.1016/j.nancom.2011.05.001

CrossRef Full Text | Google Scholar

Wishart, D. S., Feunang, Y. D., Marcu, A., Guo, A. C., Liang, K., Vazquez-Fresno, R., et al. (2018). HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. 46, D608–D617. doi: 10.1093/nar/gkx1089

CrossRef Full Text | Google Scholar

Wunderlich, Z., and Mirny, L. A. (2006). Using the topology of metabolic networks to predict viability of mutant strains. Biophys. J. 91, 2304–2311. doi: 10.1529/biophysj.105.080572

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeganeh, P. N., Richardson, C., Saule, E., Loraine, A., and Mostafavi, M. T. (2020). Revisiting the use of graph centrality models in biological pathway analysis. BioData Mining 13. doi: 10.1186/s13040-020-00214-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yizhak, K., Gaude, E., Le Dévédec, S., Waldman, Y. Y., Stein, G. Y., van de Water, B., et al. (2014). Phenotype-based cell-specific metabolic modeling reveals metabolic liabilities of cancer. Elife 3:e03641. doi: 10.7554/eLife.03641.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, J., Si, Y., Nolan, R., and Lee, K. (2007). Modular decomposition of metabolic reaction networks based on flux analysis and pathway projection. Bioinformatics 23, 2433–2440. doi: 10.1093/bioinformatics/btm374

PubMed Abstract | CrossRef Full Text | Google Scholar

Zampieri, G., Vijayakumar, S., Yaneske, E., and Angione, C. (2019). Machine and deep learning meet genome-scale metabolic modeling. PLoS Comput. Biol. 15:e1007084. doi: 10.1371/journal.pcbi.1007084

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanghellini, J., Ruckerbauer, D. E., Hanscho, M., and Jungreuthmayer, C. (2013). Elementary flux modes in a nutshell: properties, calculation and applications. Biotechnol. J. 8, 1009–1016. doi: 10.1002/biot.201200269

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Ding, G.-H., Tao, L., Yu, H., Yu, Z.-H., Luo, J.-H., et al. (2007). Modular co-evolution of metabolic networks. BMC Bioinformatics 8:311. doi: 10.1186/1471-2105-8-311

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: genome scale metabolic modeling, network science, systems biology, flux balance analysis, machine learning, synthetic biology

Citation: Dusad V, Thiel D, Barahona M, Keun HC and Oyarzún DA (2021) Opportunities at the Interface of Network Science and Metabolic Modeling. Front. Bioeng. Biotechnol. 8:591049. doi: 10.3389/fbioe.2020.591049

Received: 03 August 2020; Accepted: 22 December 2020;
Published: 25 January 2021.

Edited by:

Pablo Carbonell, Universitat Politècnica de València, Spain

Reviewed by:

Meiyappan Lakshmanan, Bioprocessing Technology Institute (A*STAR), Singapore
Hidde De Jong, INRIA, France

Copyright © 2021 Dusad, Thiel, Barahona, Keun and Oyarzún. 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) and the copyright owner(s) 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: Diego A. Oyarzún, d.oyarzun@ed.ac.uk

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.