Control and regulation of acetate overflow in Escherichia coli

  1. Pierre Millard  Is a corresponding author
  2. Brice Enjalbert
  3. Sandrine Uttenweiler-Joseph
  4. Jean-Charles Portais
  5. Fabien Létisse
  1. TBI, Université de Toulouse, CNRS, INRAE, INSA, France
  2. MetaToul-MetaboHUB, National Infrastructure of Metabolomics and Fluxomics, France
  3. RESTORE, Université de Toulouse, INSERM U1031, CNRS 5070, Université Toulouse III - Paul Sabatier, EFS, France
  4. Université Toulouse III - Paul Sabatier, France
6 figures, 3 tables and 2 additional files

Figures

Figure 1 with 3 supplements
Representation of glucose and acetate metabolism in Escherichia coli (A), in Systems Biology Graphical Notation format (http://sbgn.org) (Le Novère et al., 2009).

We performed 13C-labelling experiments to calibrate the model and evaluated the goodness of fit for different topologies (B). The initial model (model 1), which does not include inhibition of the glycolytic pathway and TCA cycle by acetate, did not fit the data satisfactorily. Adding inhibition by acetate of glycolysis (model 2) or of the TCA cycle (model 3) improved the fit, but both pathways had to be inhibited (model 4) for the goodness-of-fit criterion to be satisfied. In (B), the horizontal line represents the 95% confidence threshold for the variance-weighted sum of squared residuals (SSR). The best fits of the experimental data obtained with model 4 are shown in (C), where the shaded areas represent the 95% confidence interval on the fits. The best fits obtained with the alternative models (models 1–3) are shown in Figure 1—figure supplements 13.

Figure 1—source data 1

Experimental data used to calibrate the model.

https://cdn.elifesciences.org/articles/63661/elife-63661-fig1-data1-v2.xlsx
Figure 1—figure supplement 1
Best fit obtained for model 1.
Figure 1—figure supplement 2
Best fit obtained for model 2.
Figure 1—figure supplement 3
Best fit obtained for model 3.
Response of the E. coli transcriptome to changes in acetate concentration (0, 10, 50, or 100 mM) during growth on glucose (15 mM).

The changes in gene expression are shown in (A). Each line represents the expression of a single gene relative to its expression level measured in the absence of acetate. Up- and downregulated genes are shown in red and green, respectively. The Venn diagrams (B) represent the total number of genes upregulated (left) and downregulated (right) by at least a factor of 2 under each condition and during growth on glucose in the absence of acetate but at the same growth rate as in the presence of 100 mM acetate (0.35 hr−1, extrapolated from the data from Esquerré et al., 2014). Biological functions modulated by the presence of acetate (based on Gene Ontology analysis) are shown in (C), with the corresponding p-values. The expression levels of central metabolic genes are shown in (D). The data were obtained from four independent biological replicates for each condition.

Figure 2—source data 1

Response of the E. coli transcriptome to changes in acetate concentration (0, 10, 50, or 100 mM) during growth on glucose (15 mM).

https://cdn.elifesciences.org/articles/63661/elife-63661-fig2-data1-v2.xlsx
Figure 3 with 4 supplements
Comparison of model predictions with experimental data.

We used the model to simulate (i) steady-state glucose and acetate fluxes in glucose-limited chemostat cultures at dilution rates of 0.1–0.5 hr−1 (A), (ii) the growth rates and glucose and acetate fluxes during growth on glucose at various acetate concentrations (B), and (iii) the time courses of the changes in glucose and acetate concentrations during exponential growth on glucose after a pulse of either acetate or water (C). Model predictions are represented by lines and experimental data are shown as dots (the error bars represent one standard deviation), the shaded areas represent the 95% confidence intervals on the predictions. Predictions obtained with the alternative models (models 1–3) are shown in Figure 3—figure supplements 13. The predictive accuracy was compared between models based on the variance-weighted sum of squared residuals between simulated and experimental data (Figure 3—figure supplements 14).

Figure 3—figure supplement 1
Comparison of measured and predicted data (model 1).
Figure 3—figure supplement 2
Comparison of measured and predicted data (model 2).
Figure 3—figure supplement 3
Comparison of measured and predicted data (model 3).
Figure 3—figure supplement 4
Comparison of the predictive accuracy of models 1–4.
Heatmap of flux control coefficients during growth on glucose (15 mM) and acetate (0.1 mM).

Each column represents a controlling reaction (E) and each row, a flux (J). Red and blue cells represent negative and positive flux control coefficients (CEJ), respectively, with darker (lighter) tones indicating stronger (weaker) control.

Figure 5 with 3 supplements
Control of acetate flux over a broad range of acetate concentrations.

The shaded areas represent the 95% confidence intervals. Flux control coefficients calculated with the alternative models (models 1–3) are shown in Figure 5—figure supplements 13.

Figure 5—figure supplement 1
Flux control coefficients obtained with model 1.
Figure 5—figure supplement 2
Flux control coefficients obtained with model 2.
Figure 5—figure supplement 3
Flux control coefficients obtained with model 3.
Regulation of acetate flux in Escherichia coli.

The different routes through which acetate flux can be regulated by the acetate concentration are shown in (A). Dotted lines represent indirect (hierarchical) regulation, and straight lines represent direct (metabolic) regulation. The strengths of the three regulatory routes are respectively shown in (B–D), and their relative contributions are shown in (E). The shaded areas represent the 95% confidence intervals.

Tables

Table 1
Reactions included in the kinetic model of glucose and acetate metabolism of Escherichia coli.
NameReactionRate lawComment
Glucose_feedø → GLCCGlucose inflow and medium outflow to simulate chemostat experiments
Acetate_outflowACEenv → øMA
Biomass_outflowX → øMA
Glucose_outflowGLC → øMA
GlycolysisGLC → 1.4 × ACCOAIMMStoichiometric coefficient taken from Millard et al., 2014
TCA_cycleACCOA → øIMMUtilisation of AcCoA by the TCA cycle
PtaACCOA ↔ ACPRMMRate law from Enjalbert et al., 2017; Millard et al., 2017; Kadir et al., 2010
AckAACP ↔ ACEcellRMMRate law from Enjalbert et al., 2017; Millard et al., 2017; Kadir et al., 2010
Acetate_exchangeACEcell ↔ ACEenvRMMRate law from Millard et al., 2017
GrowthX → 2 × XMARate calculated from the TCA cycle flux, assuming a constant biomass yield (Enjalbert et al., 2017; Pinhal et al., 2019)
  1. Abbreviations: C: constant flux; MA: mass action; RMM: reversible Michaelis-Menten; IMM: irreversible Michaelis-Menten.

Table 3
Values and 95% confidence intervals of the estimated parameters.
ReactionParameterValue95 % CI
 AckAVmax3.4 × 1052.8 × 105 – 5.5 × 105
 PtaVmax9.8 × 1054.9 × 104 – 9.9 × 106
 GlycolysisVmax5.6 × 1035.3 × 103 – 5.9 × 103
Ki_ACE36.730.9 – 46.9
 TCA cycleKm_ACCOA24.88.4 – 615.4
Vmax7.4 × 1052.4 × 105 – 1.7 × 106
Ki_ACE2.31.8 – 3.4
 GrowthY1.0 × 10−49 × 10−5 – 1.1 × 10−4
 Acetate exchangeVmax4.8 × 1058 × 104 – 1.5 × 106
Km_ACE33.21.5 – 99.8

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Pierre Millard
  2. Brice Enjalbert
  3. Sandrine Uttenweiler-Joseph
  4. Jean-Charles Portais
  5. Fabien Létisse
(2021)
Control and regulation of acetate overflow in Escherichia coli
eLife 10:e63661.
https://doi.org/10.7554/eLife.63661