Daniele De Martino,Anna MC Andersson,Tobias Bergmiller,Călin C Guet,Gašper Tkačik
Abstract:Which properties of metabolic networks can be derived solely from stoichiometric information about the network's constituent reactions? Predictive results have been obtained by Flux Balance Analysis (FBA), by postulating that cells set metabolic fluxes within the allowed stoichiometry so as to maximize their growth. Here, we generalize this framework to single cell level using maximum entropy models from statistical physics. We define and compute, for the core metabolism of Escherichia coli, a joint distribution over all fluxes that yields the experimentally observed growth rate. This solution, containing FBA as a limiting case, provides a better match to the measured fluxes in the wild type and several mutants. We find that E. coli metabolism is close to, but not at, the optimality assumed by FBA. Moreover, our model makes a wide range of predictions: (i) on flux variability, its regulation, and flux correlations across individual cells; (ii) on the relative importance of stoichiometric constraints vs. growth rate optimization; (iii) on quantitative scaling relations for singe-cell growth rate distributions. We validate these scaling predictions using data from individual bacterial cells grown in a microfluidic device at different sub-inhibitory antibiotic concentrations. Under mild dynamical assumptions, fluctuation-response relations further predict the autocorrelation timescale in growth data and growth rate adaptation times following an environmental perturbation.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is: **Can the properties of the metabolic network during steady - state growth be deduced solely from the stoichiometric information of its constituent reactions?**
Specifically, the authors are concerned with how to understand the behavior of the metabolic network at the single - cell level and go beyond the traditional Flux Balance Analysis (FBA) method. FBA assumes that the metabolic fluxes within the cell are the result of maximizing the growth rate under the allowed stoichiometric constraints. However, FBA ignores the heterogeneity between cells, that is, there may be differences in the metabolic behavior among different cells.
To answer this question, the authors introduced the Maximum Entropy Model in statistical physics and applied it to the core metabolic network of Escherichia coli. They defined and calculated a joint distribution that can explain the experimentally observed growth rates. This model not only includes FBA as a limiting case but also can better match the measured metabolic fluxes in the wild - type and various mutants.
### Main problems and solutions
1. **Problem description**:
- The traditional FBA method assumes that cells maximize the growth rate by adjusting metabolic fluxes under steady - state conditions, but ignores the heterogeneity between cells.
- How to understand and predict the behavior of the metabolic network at the single - cell level?
2. **Solution**:
- Use the Maximum Entropy Model to construct a joint distribution that can achieve the specified average growth rate under stoichiometric constraints.
- By adjusting the optimization parameter β, the model can interpolate between a uniform distribution (no optimization) and FBA (full optimization).
- Verify the model's prediction of the metabolic fluxes of Escherichia coli wild - type and mutants, and find that the Maximum Entropy Model provides a better match at the intermediate optimization level (rather than full optimization).
### Key points of the model
- **Maximum entropy distribution**: \[ P_\beta(v) = \begin{cases}
\frac{1}{Z} e^{\beta \lambda(v)}, & v \in P \\
0, & v \notin P
\end{cases} \]
where \( Z \) is the normalization constant, \( \lambda(v) \) is the growth rate function, and \( P \) is the feasible metabolic phenotype space.
- **Optimization parameter β**: By adjusting β, the transition of the model from a uniform distribution to the FBA solution can be controlled, thereby capturing the heterogeneity between cells.
- **Prediction and verification**:
- The model predicts the mean, fluctuation, and correlation of metabolic fluxes.
- Experiments verify the model's prediction of growth rate fluctuations, especially the single - cell growth rate distribution under different antibiotic concentrations.
Through this method, the authors not only extended the application range of FBA but also revealed the fact that the Escherichia coli metabolic network operates in a state close to but not completely optimal, and proposed new insights into the cost of metabolic regulation.