 Short Report
 Open Access
 Published:
Functional centrality as a predictor of shifts in metabolic flux states
BMC Research Notesvolume 9, Article number: 317 (2016)
Abstract
Background
The flux phenotype describes the entirety of biochemical conversions in a cell, which renders it a key characteristic of metabolic function. To quantify the functional relevance of individual biochemical reactions, functional centrality has been introduced based on cooperative game theory and structural modeling. It was shown to be capable to determine metabolic control properties utilizing only structural information. Here, we demonstrate the capability of functional centrality to predict changes in the flux phenotype.
Results
We use functional centrality to successfully predict changes of metabolic flux triggered by switches in the environment. The predictions via functional centrality improve upon predictions using controleffective fluxes, another measure aiming at capturing metabolic control using structural information.
Conclusions
The predictions of flux changes via functional centrality corroborate the capability of the measure to gain a mechanistic understanding of metabolic control from the structure of metabolic networks.
Background
The metabolic state of a cellular system is characterized by the fluxes of the underlying biochemical reactions and the concentrations of the involved metabolites [1]. In the analysis of metabolism, the flux phenotype is one of the most important cellular observables as it directly relates to metabolic functionality, i.e., the conversion of individual biochemical compounds [2], which also naturally accounts for changes in metabolic state.
Alterations in the flux phenotype are triggered by internal or external shifts and imply rerouting of metabolic flux, for example from respiration towards fermentation in facultative aerobic bacteria [3], which is the result of metabolic control [1]. Several constraintbased approaches were proposed to predict flux phenotypes in unicellular organisms resulting from genetic interventions by minimizing their effect on flux distributions [4, 5]. Nevertheless, elucidating the principles of how largescale metabolic systems achieve efficient metabolic rerouting remains one of the key challenges in developing a mechanistic understanding of metabolic control, in accordance with the general case of controlling complex networks [6]. Here, we demonstrate the ability of our recently introduced measure of functional centrality to predict changes in experimentally determined metabolic fluxes based solely on the structure of the metabolic network captured by its stoichiometric matrix.
Functional centrality combines cooperative game theory, in particular a modified version of the Shapley value for arbitrary restricted games [7], with flux balance analysis (FBA) [8] to determine the functional relevance of biochemical reactions [1, 9]. It quantifies the contribution of individual biochemical reactions to metabolic functionality, e.g., to biomass production, based on the structure of a metabolic network together with constraints on internal and exchange fluxes, which mimic internal and external conditions.
In our previous study we examined the relationship of functional relevance and transcriptional readout [1], whereby we examined four measures based on structural modeling and associated with metabolic control. Our comparative analysis showed that controleffective fluxes [10] were superior in predicting changes of gene expression. However, gene expression may not directly be manifested in changes of flux; a more suitable comparison necessitates the usage of changes in fluxes inferred based on labeling data [11]. In this brief research note, we demonstrate that functional centrality performs better than controleffective fluxes in predicting changes in metabolic flux in the central carbon metabolism of Escherichia coli.
Results
We examine the capability to predict flux changes using functional centralities as well as controleffective fluxes, which were found to be the two bestperforming contending predictors of transcriptional change in a recent study of structural control of metabolism [1]. To this end, we analyze the relationship of flux changes in the central carbon metabolism of E. coli resulting from a shift between aerobic and fermentative conditions and the corresponding changes of the two investigated measures. The flux data were obtained via stationary ^{13}C labeling experiments and metabolic flux analysis by Fischer et al. [12].
In accordance with the definition of controleffective fluxes [10], we utilize the mean of normalized functional centralities obtained for the metabolic functions of ATP and biomass production, which are the two dominant metabolic functions shaping the flux phenotype [13]. The controleffective fluxes and the mean functional centralities are calculated for aerobic as well as for fermentative conditions using the metabolic network of E. coli’s central carbon metabolism provided by Schuetz et al. [13]. For both measures, we calculate the ratios of values obtained for the two conditions and use them as predictors of ratios of the corresponding fluxes. Estimates of functional centralities are used, because exact computation is computationally intractable for large metabolic network. We demonstrated that robust estimates of functional centralities can be obtained via Monte Carlo sampling utilizing the set of elementary flux modes in mediumsize metabolic networks [1].
The ability of controleffective fluxes and functional centrality to predict the observed flux changes (23 data points) is assessed by employing two statistical procedures. In the first, we fit a regression model. Inspection of the dataset indicates a nonlinear relationship in the case of both measures. We utilize logarithmic transformation of the predictor’s and of the regressor’s values, respectively, to discount nonlinearity and to render the distribution of the predicted quantity closer to normal. A linear regression model is fitted to the doublelogtransformed data (residuals do not show significant deviations from a normal distribution). The slope of the linear regression line then indicates the relative (percental) change of flux resulting from a relative (percental) change of the predictor. The goodness of fit is determined by the coefficient of determination
whereby \(y_{i}\) and \(y_{i}^{*}\) denote the ith (logtransformed) data point and the (logtransformed) predicted value, respectively; \(R_{i}=y_{i}y_{i}^{*}\) denotes the ith residual and \(\bar{y}\) the mean of data points. In the second, we use Kendall’s rank correlation coefficient. We find that functional centrality improves prediction of metabolic fluxes in comparison to controleffective fluxes (see Figs. 1, 2). While controleffective fluxes can only explain \(35.1\,\%\) of the variance of measured fluxes, functional centrality is capable to explain \(44.21\,\%\) of the variance. The association found by utilizing Kendall’s rank correlation coefficient is even more pronounced as shown in Table 1, which also provides values of Pearson’s and Spearman’s correlation coefficient.
Conclusions
Functional centrality, like controleffective fluxes, does not require any further information besides reaction stoichiometry, characterization of exchange fluxes and the choice of a metabolic function. Here we demonstrated that functional centrality enables prediction of changes in metabolic flux with improved accuracy, exceeding the results of controleffective fluxes. The functional centrality measure is based on the rich axiomatic framework of power indices in cooperative games, and is suitable for the analysis of structural control of metabolic networks, as supported by the analyzed data. Our results provide further indication that concepts from cooperative game theory can help in obtaining a better understanding of control in metabolic networks.
Methods
Flux balance analysis
Flux balance analysis is a structural modeling framework enabling prediction of steadystate metabolic flux [8]. It relies on the assumption that metabolism is governed by an optimization principle with respect to some cellular objective. Usually, the objective is expressed as a linear combination of individual fluxes \(f_{i}\) (\(i\in \mathcal {N},\) with \(\mathcal {N}\) denoting the set of reactions forming the metabolic network). Optimization can then be formulated in terms of a linear programming problem:
with z representing the objective to be optimized, \(\mathbf {c}\) being a vector of coefficients quantifying the contribution of each flux to this objective, and matrix \(\mathbf {S}\) capturing reaction stoichiometries. The bounds \(f_{min}(i)\) and \(f_{max}(i)\), denote minimum and maximum values of the fluxes and, therefore, determine reaction reversibility.
Functional centrality
Functional centrality quantifies the contribution of individual reactions to a metabolic function, e.g., biomass production, by utilizing a formulation of the Shapley value for restricted cooperative games [9]. The metabolic function of interest is formulated in terms of a linear objective function accessible via FBA. We provide briefly the definition of functional centrality, an extensive derivation and an algorithm for estimation of functional centrality in large networks is provided by SajitzHermstein and Nikoloski [1, 9].
The optimal value of the objective function v in a subnetwork formed by the reaction set S is determined by FBA. Let \(G_{\mathcal {S}}=(\mathcal {S},\mathcal {A})\) be a directed graph. The set of nodes \(\mathcal {S}\) encompasses all subsets \(S\subseteq \mathcal {N}\) which correspond to functional subnetworks, i.e., exhibiting positive outcome of FBA, and the empty set \(\emptyset\). The set of arcs \(\mathcal {A}\) consists of all \((S,S')\) with \(S\subsetneq S'\) for which it holds that there exists no \(S''\in \mathcal {S}\) with \(S\subsetneq S''\subsetneq S'\). In \(G_{\mathcal {S}}\), every path from the empty set to the set \(\mathcal {N}\) represents one possibility to add reactions successively in such a way that the corresponding subnetworks belong to the family of functional subnetworks. These paths are called maximal chains. The calculation of functional centrality is as follows:
Let \(W=(S_{0},S_{1},\ldots ,S_{l(W)})\) with \(S_{0}=\emptyset\) and \(S_{l(W)}=\mathcal {N}\) be a maximal chain for inclusion in \(\mathcal {S}\), implying \(S_{0}\subsetneq S_{1}\subsetneq \cdots \subsetneq S_{l(W)}.\) Then, the contribution of reaction i in maximal chain W is given by
with l(W) being the length of the maximal chain W and \(\cdot \) denoting the cardinality of a set. Let the set of all maximal chains be denoted by \(\mathcal {W}.\) The functional centrality \(\Phi _{i}(v)\) of reaction i with respect to the objective function v is
Controleffective fluxes
Controleffective fluxes are defined by efficiencies \(\varepsilon _{i}\) of the individual elementary flux modes (EFM) \(e_{i}\) with respect to a substrate \(S_{k}\) and the production of biomass (\(\mu\)) and ATP (EFMs are normalized by substrate uptake in advance) [10]:
whereby \(e_{i}^{l}\) denotes the flux through reaction l in the EFM \(e_{i}.\) The controleffective fluxes are then given by
with \(Y_{(\mu /ATP)/S_{k}}^{\mathrm {max}}\) being the maximum yield of biomass or ATP production, respectively, for substrate \(S_{k}\).
Model
We utilize a model of central carbon metabolism of E. coli originally published by Schuetz et al. [13]. Import of acetate and ethanol is disabled in our study. The sum of glucose import by the reactions mglABC and ptsGHI is constrained from above arbitrarily by one. Isozymes in the model were deleted. We examine two environmental conditions: (i) aerobic respiration and (ii), fermentation (no oxygen and no nitrate import). The utilized objective functions are the fluxes through the reactions: (i) maint (ATP production) and (ii), biomass (biomass production).
Abbreviations
 FBA:

flux balance analysis
 EFM:

elementary flux modes
References
 1.
SajitzHermstein M, Nikoloski Z. Structural control of metabolic flux. PLoS Comput Biol. 2013;9(12):e1003368.
 2.
Nielsen J. It is all about metabolic fluxes. J Bacteriol. 2003;185(24):7031–5.
 3.
Unden G, Becker S, Bongaerts J, Schirawski J, Six S. Oxygen regulated gene expression in facultatively anaerobic bacteria. Antonie van Leeuwenhoek. 1994;66:3–22.
 4.
Segrè D, Vitkup D, Church GM. Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002;99(23):15112–7.
 5.
Shlomi T, Berkman O, Ruppin E. Regulatory on/off minimization of metabolic flux changes after genetic perturbations. Proc Natl Acad Sci USA. 2005;102(21):7695–700.
 6.
Liu YY, Slotine JJ, Barabási AL. Controllability of complex networks. Nature. 2011;473:167–73.
 7.
Aguilera NE, Di Marco SC, Escalante MS. The Shapley value for arbitrary families of coalitions. Eur J Oper Res. 2010;204:125–38.
 8.
Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010;28(3):245–8.
 9.
SajitzHermstein M, Nikoloski Z. Restricted cooperative games on metabolic networks reveal functionally important reactions. J Theor Biol. 2012;314:192–203.
 10.
Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED. Metabolic network structure determines key aspects of functionality and regulation. Nature. 2002;420:190–3.
 11.
Kruger N, Ratcliffe RG. Fluxes through plant metabolic networks: measurements, predictions, insights and challenges. Biochem J. 2015;465:27–38.
 12.
Fischer E, Zamboni N, Sauer U. Highthroughput metabolic flux analysis based on gas chromatography–mass spectrometry derived ^{13}C constraints. Anal Biochem. 2004;325:308–16.
 13.
Schuetz R, Kuepfer L, Sauer U. Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol Syst Biol. 2007;3:119.
Authors’ contributions
MSH and ZN designed and conceived the study. MSH carried out the computational analysis. MSH and ZN wrote the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The utilized metabolic network model as well as the flux data is available as supplementary material of [13].
Consent for publication
Not applicable.
Ethics approval and consent to participate
The study did not involve human participants, human data or human tissue.
Funding
MSH and ZN would like to acknowledge the financial support of PROMICS, a Research Unit 1186 of the German Research Foundation, as well as financial support by the Max Planck Society.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Functional centrality
 Metabolic network
 Flux balance analysis
 Elementary flux mode
 Reaction stoichiometry
 Structural modeling
 Metabolic control
 Metabolic flux
 Cooperative game theory
 Shapley value