Noisy metabolism can promote microbial cross-feeding

  1. Jaime G Lopez
  2. Ned S Wingreen  Is a corresponding author
  1. Lewis-Sigler Institute for Integrative Genomics, Princeton University, United States
  2. Department of Molecular Biology, Princeton University, United States
13 figures and 1 additional file

Figures

Isolated bacterial cells suffer negative growth effects from noisy enzyme regulation and metabolite leakage.

(A) Schematic of modeled intracellular dynamics. Cells import an external nutrient (yellow hexagon) that can be converted by enzymes (magenta and blue) into two essential metabolites. Metabolites passively exchange with the extracellular medium (‘leakage’), and degrade at a fixed rate. The two metabolites are used for growth in accord with Liebig’s law of the minimum. (B) Schematic of dynamic enzyme regulation scheme: the type of enzyme produced is always the one associated with the lower metabolite pool. (C) Metabolite timecourse of a cell that produces enzymes in small bursts. See Equations 1; 2 for details; parameters β=2, γ=50, κ=100, rV=10, P=1, δ=1, g*=1×10-5. Time is normalized by the inverse of the maximum steady-state growth rate (t=tκγg*2(1+δ)). (D) Metabolite timecourse of a cell that produces enzymes in large bursts. Same parameters as in C but with β=20. (E) Average growth rate of an isolated cell for differing values of permeability; parameters as in C and D except as specified. Growth rate is normalized to the maximum possible growth rate, i.e. with perfect regulation and zero permeability. (F) Average growth rate of an isolated cell for differing values of burst size. Parameters as in C and D except as specified.

Bacterial cells can compensate for noisy enzyme regulation and increase growth rate by exchanging metabolites within a clonal community.

(A) Schematic of multi-cell metabolism model. Individual cells regulate their own enzyme levels, but metabolites leak into the local medium and can be used by other cells in the community. (B) Intracellular metabolite coefficient of variation (CV) for a community of 10 cells as a function of cell permeability. (C) Average growth rate for community of 10 cells as a function of cell permeability. Parameters in B and C same as in Figure 1E.

Sharing multiple metabolites can generically reduce noise and improve overall colony growth rate.

(A) Implication of Liebig’s law of the minimum for fluctuating metabolites: the growth rate at any time is set by the lowest metabolite level, hence the average growth rate is lower than at the average metabolite level. The magnitude of this decrease grows with increasing metabolite variance. (B) Distributions of growth rates, set by minimum metabolite level Min(mi), for varying numbers of essential metabolites. (C) Mean growth rate as a function of number of metabolites and of metabolite CV. The dashed curve corresponds to CV =0.4, as measured for essential proteins in E. coli (Taniguchi et al., 2010). (D) Mean growth rate as a function of number of metabolites and of the correlation coefficient between metabolites, for CV=0.4. The dashed curve indicates a metabolite correlation of 0.7, approximately that observed for proteins in E. coli (Taniguchi et al., 2010).

Appendix 9—figure 1
Enzyme timecourses corresponding to Figure 1CD.

(A) Enzyme timecourse of a cell that produces enzymes in small bursts, same parameters as in Figure 1C. (B) Enzyme timecourse of a cell that produces enzymes in large bursts, same parameters as in Figure 1D. Note that the enzyme levels are substantially correlated with each other (r=0.89 for β=2 and r=0.53 for β=20). This correlation between enzyme levels also results in correlation between metabolite levels (r=0.8 for β=2 and r=0.25 for β=20).

Appendix 9—figure 2
Metabolite CV corresponding to simulations in Figure 1E.

All simulation parameters identical to those in Figure 1E. Data corresponding to β=100 not shown as metabolite CV cannot be meaningfully estimated for cells with arrested metabolism.

Appendix 9—figure 3
Impact of extracellular volume on growth rate of isolated cells.

As extracellular volume is increased, cells are increasingly less likely to recover their leaked metabolites, leading to decreased growth. This effect plateaus once the extracellular volume becomes effectively infinite such that cells never recover their leaked metabolites. The dashed black line represents the maximum growth rate of cells at infinite volume, calculated by treating the permeability as an additional intracellular metabolite degradation term (i.e. g=κγg*2(1+δ+P)). Parameters used are the same as in 1E with a permeability of P=1. Values presented are averaged from 100 replicate simulations.

Appendix 9—figure 4
Comparison of Langevin equation predictions to full metabolism model.

The full model is simulated with the same parameters as in Figure 1EF, except with V=0 and P=100. We apply the following substitutions to the Langevin model in order to directly compare it to the full model: Γ=γgmax, μE=gmax, and μm=δ+1, where gmax=κγg*2(1+δ). (A) Comparison of metabolite CV in the full model (solid curves) and the Langevin model (dashed curves) as a function of community size at different burst sizes. (B) Comparison of metabolite CV in the full model (solid curve) and the Langevin model (dashed curve) as a function of burst size. As expected, the agreement between the full model and the Langevin model is best at high N and low β - conditions that lead to low levels of noise.

Appendix 9—figure 5
Estimated fractional production losses of amino acids in E. coli due to leakage.

E. coli due to leakage. In order to determine the significance of amino-acid leakage in E. coli, we estimated the fraction of E. coli’s amino-acid production that is lost to leakage. Formally, we define the fraction of production lost as the ratio of the number of amino acids lost via leakage over the period of one division to the number of amino acids required to produce a daughter cell. For each amino acid, we require three experimentally measured quantities for this calculation: (1) the rate of leakage through the cell membrane, (2) the intracellular concentration of the amino acid, and (3) the number of amino acids required to produce a daughter cell. For leakage rates, we use data from artificial liposomes (Chakrabarti and Deamer, 1992). In cases where multiple pH conditions were tested, we used data measured at pH 7 (though leakage rates did not vary substantially with pH). This study measured data for only a limited set of amino acids. For other amino acids, we estimated their leakage rates using a linear regression of leakage rate versus log octanol/water partition coefficient (r2=0.96). Leakage rates were also adjusted for the differing size of the liposomes and typical E. coli cells, assuming a liposome radius of 100 nm and an E. coli radius of 400 nm (Grossman et al., 1982). Intracellular concentrations were taken from (Bennett et al., 2009) and per-cell pool sizes were calculated assuming a cell volume of 1.8×10-15L (Outten and O’Halloran, 2001). A cell’s amino-acid production was assumed to be the number of amino acids required to produce a daughter cell, taken from (Mee et al., 2014). With all of these experimental values, the fraction of production lost is f=kNIτNtot where k is the leakage rate, NI is the intracellular molecule count, τ is the doubling time (assumed to be 24 minutes), and Ntot is the number of amino acids required to produce a daughter cell.

Appendix 9—figure 6
Version of Figure 3C with an alternative growth function.

To demonstrate that our findings are not limited to Liebig’s law of the minimum, we repeat the analysis in Figure 3C with an alternative growth function from (Goyal and Wingreen, 2007). The function is g=gmax1Ni=1Nmi+mi*mi, where gmax is the maximum growth rate, N is the total number of unique metabolites, and mi* is the half-substrate constant of each metabolite. In this analysis we assume gmax=1 and mi=E[mi].

Appendix 9—figure 7
Correlation between initial conditions and later times in the full metabolism model.

To compute correlations to the initial condition, we perform 2,500 replicate simulations of single-cell metabolism with two different values of burst size β, each with the same initial conditions except for the value of E1, which is uniformly distributed in the range [0,γ]. Simulations are run for 12.5 time units, where time is normalized to the maximum steady-state growth rate, as in Figure 1. We bin the resulting timeseries with a window size of 0.025 time units, computing Pearson correlations between the initial conditions and all data points within the window. From each bin, we also compute the mean time of datapoints within the bin, which we use as the x-axis data in this figure. As can be seen, the system rapidly forgets its initial conditions, nearing zero correlation after a short time lag.

Appendix 9—figure 8
Examples of non-globally concave growth functions that still allow NAC to be potentially beneficial.

The black curve is g(x)=x/(0.3+x)+0.02sin(100x) and the red curve is g(x)=1/(1+e-5x+2.1)-1/(1+e2.1). For metabolites uniformly distributed in the domain [0,1], both functions have negative Jensen gaps and are thus negatively impacted by noise.

Appendix 9—figure 9
Spatial profiles of NAC metabolites within a biofilm.

Here, we plot solutions of the biofilm model analyzed and parameterized in Appendix 7. The dashed lines represent the boundary between the NAC and non-NAC sectors of the biofilm, with the NAC sector beginning at r^=0. Left: Metabolite profiles with a NAC sector of radius a^=5μm. Right: Metabolite profiles with a NAC sector of radius a^=20μm.

Appendix 9—figure 10
Analysis of metabolism model without metabolite degradation.

(A–B) Versions of Figure 1EF with δ=0. Panel A shows an extended range of permeability values to show the non-monotonic growth behavior. Data in A is computed from 100 replicates and the data from B is computed from 200 replicates. All parameters identical to the original figure with the exception of the degradation rate, δ. Simulations are initialized with zero extracellular metabolites. (C–D) Timecourses of metabolite and enzyme levels in an isolated cell exhibiting build-up and crash behavior. Parameters the same as in Figure 1D with the exception of P=0 and δ=0. Time normalized by inverse of maximum steady-state growth rate. (E–F) Versions of Figure 2BC with δ=0. Data is computed from 40 replicates. All parameters identical to the original figure with the exception of the degradation rate, δ. Simulations are initialized with zero extracellular metabolites.

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. Jaime G Lopez
  2. Ned S Wingreen
(2022)
Noisy metabolism can promote microbial cross-feeding
eLife 11:e70694.
https://doi.org/10.7554/eLife.70694