Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Res Notes**|**v.9; 2016**|**PMC4915090

Formats

Article sections

Authors

Related links

BMC Res Notes. 2016; 9: 317.

Published online 2016 June 21. doi: 10.1186/s13104-016-2117-0

PMCID: PMC4915090

Systems Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant Physiology, Am Mühlenweg 1, 14476 Potsdam, Germany

Max Sajitz-Hermstein, Email: ed.gpm.mlog-pmipm@ztijas.

Received 2015 September 29; Accepted 2016 June 2.

Copyright © The Author(s) 2016

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.

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 control-effective fluxes, another measure aiming at capturing metabolic control using structural information.

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.

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 constraint-based 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 large-scale 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 control-effective 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 control-effective fluxes in predicting changes in metabolic flux in the central carbon metabolism of *Escherichia coli*.

We examine the capability to predict flux changes using functional centralities as well as control-effective fluxes, which were found to be the two best-performing 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 control-effective 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 control-effective 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 medium-size metabolic networks [1].

The ability of control-effective 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 non-linear relationship in the case of both measures. We utilize logarithmic transformation of the predictor’s and of the regressor’s values, respectively, to discount non-linearity and to render the distribution of the predicted quantity closer to normal. A linear regression model is fitted to the double-log-transformed 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

$$\begin{array}{c}\hfill {R}^{2}=1-\frac{{\sum}_{i}{R}_{i}^{2}}{{\sum}_{i}{({y}_{i}-\overline{y})}^{2}},\end{array}$$

1

whereby *y*_{i} and ${y}_{i}^{\ast}$ denote the *i*th (log-transformed) data point and the (log-transformed) predicted value, respectively; ${R}_{i}={y}_{i}-{y}_{i}^{\ast}$ denotes the *i*th residual and $\overline{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 control-effective fluxes (see Figs. Figs.1,1, ,2).2). While control-effective 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.

Prediction of flux changes by control-effective fluxes. Calculated ratios between flux during growth on glucose (GLC) under conditions of aerobic respiration and fermentative conditions based on control-effective fluxes versus experimentally determined **...**

Prediction of flux changes by functional centralities. Calculated ratios between flux during growth on glucose (GLC) under conditions of aerobic respiration and fermentative conditions based on functional centralities versus experimentally determined **...**

Functional centrality, like control-effective 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 control-effective 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.

Flux balance analysis is a structural modeling framework enabling prediction of steady-state 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* ∈ 𝒩, with 𝒩 denoting the set of reactions forming the metabolic network). Optimization can then be formulated in terms of a linear programming problem:

$$\begin{array}{c}\hfill \text{min(max)}z=\sum _{i}{c}_{i}{f}_{i},\phantom{\rule{0.333333em}{0ex}}\text{s.t.}\\ \hfill \mathbf{S}\xb7\mathbf{f}=\mathbf{0},\\ \hfill {f}_{min}(i)\le {f}_{i}\le {f}_{max}(i),\end{array}$$

2

with *z* representing the objective to be optimized, **c** being a vector of coefficients quantifying the contribution of each flux to this objective, and matrix **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 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 Sajitz-Hermstein 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*_{𝒮} = (𝒮, 𝒜) be a directed graph. The set of nodes 𝒮 encompasses all subsets *S* ⊆ 𝒩 which correspond to functional subnetworks, i.e., exhibiting positive outcome of FBA, and the empty set ∅. The set of arcs 𝒜 consists of all (*S*, *S*^{′}) with *S* ⊊ *S*^{′} for which it holds that there exists no *S*^{′′} ∈ 𝒮 with *S* ⊊ *S*^{′′} ⊊ *S*^{′}. In *G*_{𝒮}, every path from the empty set to the set 𝒩 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}, …, *S*_{l(W)}) with *S*_{0} = ∅ and *S*_{l(W)} = 𝒩 be a maximal chain for inclusion in 𝒮, implying *S*_{0} ⊊ *S*_{1} ⊊ ⋯ ⊊ *S*_{l(W)}. Then, the contribution of reaction *i* in maximal chain *W* is given by

$$\begin{array}{c}\hfill {\mathit{\psi}}_{W,i}(v)=\sum _{j=1}^{l(W)}\frac{v({S}_{j})-v({S}_{j-1})}{|{S}_{j}|-|{S}_{j-1}|}\xb7\mathit{\chi}({S}_{j}\backslash {S}_{j-1},i)\phantom{\rule{1em}{0ex}}\phantom{\rule{0.333333em}{0ex}}\text{for}\phantom{\rule{0.333333em}{0ex}}i\in \mathcal{N},\end{array}$$

3

$$\begin{array}{c}\hfill \mathit{\chi}(S,i)=\left\{\begin{array}{cc}1\hfill & \phantom{\rule{1em}{0ex}}\text{if}\phantom{\rule{0.333333em}{0ex}}\phantom{\rule{0.166667em}{0ex}}i\in S\hfill \\ 0\hfill & \phantom{\rule{1em}{0ex}}\text{otherwise},\hfill \end{array}\right.\end{array}$$

4

with *l*(*W*) being the length of the maximal chain *W* and | · | denoting the cardinality of a set. Let the set of all maximal chains be denoted by 𝒲. The functional centrality Φ_{i}(*v*) of reaction *i* with respect to the objective function *v* is

$$\begin{array}{c}\hfill {\mathit{\varphi}}_{i}(v)=\sum _{W\in \mathcal{W}}\frac{1}{|\mathcal{W}|}{\mathit{\psi}}_{W,i}(v).\end{array}$$

5

Control-effective fluxes are defined by *efficiencies**ε*_{i} of the individual elementary flux modes (EFM) *e*_{i} with respect to a substrate *S*_{k} and the production of biomass (*μ*) and ATP (EFMs are normalized by substrate uptake in advance) [10]:

$$\begin{array}{c}\hfill {\mathit{\epsilon}}_{i}({S}_{k},\mathit{\mu})=\frac{{e}_{i}^{\mathit{\mu}}}{{\sum}_{l}|{e}_{i}^{l}|},\end{array}$$

6

$$\begin{array}{c}\hfill {\mathit{\epsilon}}_{i}({S}_{k},\mathrm{ATP})=\frac{{e}_{i}^{\mathrm{ATP}}}{{\sum}_{l}|{e}_{i}^{l}|},\end{array}$$

7

whereby ${e}_{i}^{l}$ denotes the flux through reaction *l* in the EFM *e*_{i}. The control-effective fluxes are then given by

$$\begin{array}{cc}\hfill {C}_{l}({S}_{k})& =\frac{1}{{Y}_{\mathit{\mu}/{S}_{k}}^{\mathrm{max}}}\xb7\frac{{\sum}_{i}{\mathit{\epsilon}}_{i}({S}_{k},\mathit{\mu})\xb7|{e}_{i}^{l}|}{{\sum}_{i}{\mathit{\epsilon}}_{i}({S}_{k},\mathit{\mu})}+\frac{1}{{Y}_{\mathrm{ATP}/{S}_{k}}^{\mathrm{max}}}\hfill \\ \hfill & \phantom{\rule{1em}{0ex}}\xb7\frac{{\sum}_{i}{\mathit{\epsilon}}_{i}({S}_{k},\mathrm{ATP})\xb7|{e}_{i}^{l}|}{{\sum}_{i}{\mathit{\epsilon}}_{i}({S}_{k},\mathrm{ATP})},\hfill \end{array}$$

8

with ${Y}_{(\mathit{\mu}/ATP)/{S}_{k}}^{\mathrm{max}}$ being the maximum yield of biomass or ATP production, respectively, for substrate *S*_{k}.

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).

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.

The authors declare that they have no competing interests.

The utilized metabolic network model as well as the flux data is available as supplementary material of [13].

Not applicable.

The study did not involve human participants, human data or human tissue.

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.

- FBA
- flux balance analysis
- EFM
- elementary flux modes

Max Sajitz-Hermstein, Email: ed.gpm.mlog-pmipm@ztijas.

Zoran Nikoloski, Email: ed.gpm.mlog-pmipm@iksolokin.

1. Sajitz-Hermstein M, Nikoloski Z. Structural control of metabolic flux. PLoS Comput Biol. 2013;9(12):e1003368. doi: 10.1371/journal.pcbi.1003368. [PMC free article] [PubMed] [Cross Ref]

2. Nielsen J. It is all about metabolic fluxes. J Bacteriol. 2003;185(24):7031–7035. doi: 10.1128/JB.185.24.7031-7035.2003. [PMC free article] [PubMed] [Cross Ref]

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. doi: 10.1007/BF00871629. [PubMed] [Cross Ref]

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–15117. doi: 10.1073/pnas.232349399. [PubMed] [Cross Ref]

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–7700. doi: 10.1073/pnas.0406346102. [PubMed] [Cross Ref]

6. Liu YY, Slotine JJ, Barabási AL. Controllability of complex networks. Nature. 2011;473:167–173. doi: 10.1038/nature10011. [PubMed] [Cross Ref]

7. Aguilera NE, Di Marco SC, Escalante MS. The Shapley value for arbitrary families of coalitions. Eur J Oper Res. 2010;204:125–138. doi: 10.1016/j.ejor.2009.09.031. [Cross Ref]

8. Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010;28(3):245–248. doi: 10.1038/nbt.1614. [PMC free article] [PubMed] [Cross Ref]

9. Sajitz-Hermstein M, Nikoloski Z. Restricted cooperative games on metabolic networks reveal functionally important reactions. J Theor Biol. 2012;314:192–203. doi: 10.1016/j.jtbi.2012.08.018. [PubMed] [Cross Ref]

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–193. doi: 10.1038/nature01166. [PubMed] [Cross Ref]

11. Kruger N, Ratcliffe RG. Fluxes through plant metabolic networks: measurements, predictions, insights and challenges. Biochem J. 2015;465:27–38. doi: 10.1042/BJ20140984. [PubMed] [Cross Ref]

12. Fischer E, Zamboni N, Sauer U. High-throughput metabolic flux analysis based on gas chromatography–mass spectrometry derived ^{13}C constraints. Anal Biochem. 2004;325:308–316. doi: 10.1016/j.ab.2003.10.036. [PubMed] [Cross Ref]

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. doi: 10.1038/msb4100162. [PMC free article] [PubMed] [Cross Ref]

Articles from BMC Research Notes are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |