Skip to main content
Advertisement
  • Loading metrics

Plasticity of growth laws tunes resource allocation strategies in bacteria

  • Avik Mukherjee ,

    Contributed equally to this work with: Avik Mukherjee, Yu-Fang Chang

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Yu-Fang Chang ,

    Contributed equally to this work with: Avik Mukherjee, Yu-Fang Chang

    Roles Investigation, Methodology, Resources, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Yanqing Huang,

    Roles Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Nina Catherine Benites,

    Roles Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Leander Ammar,

    Roles Investigation, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Jade Ealy,

    Roles Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Mark Polk,

    Roles Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

  • Markus Basan

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    markus@hms.harvard.edu

    Affiliation Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America

Abstract

Bacteria like E. coli grow at vastly different rates on different substrates, however, the precise reason for this variability is poorly understood. Different growth rates have been attributed to ‘nutrient quality’, a key parameter in bacterial growth laws. However, it remains unclear to what extent nutrient quality is rooted in fundamental biochemical constraints like the energy content of nutrients, the protein cost required for their uptake and catabolism, or the capacity of the plasma membrane for nutrient transporters. Here, we show that while nutrient quality is indeed reflected in protein investment in substrate-specific transporters and enzymes, this is not a fundamental limitation on growth rate, at least for certain ‘poor’ substrates. We show that it is possible to turn mannose, one of the ‘poorest’ substrates of E. coli, into one of the ‘best’ substrates by reengineering chromosomal promoters of the mannose transporter and metabolic enzymes required for mannose degradation. This result falls in line with previous observations of more subtle growth rate improvement for many other carbon sources. However, we show that this faster growth rate comes at the cost of diverse cellular capabilities, reflected in longer lag phases, worse starvation survival and lower motility. We show that addition of cAMP to the medium can rescue these phenotypes but imposes a corresponding growth cost. Based on these data, we propose that nutrient quality is largely a self-determined, plastic property that can be modulated by the fraction of proteomic resources devoted to a specific substrate in the much larger proteome sector of catabolically activated genes. Rather than a fundamental biochemical limitation, nutrient quality reflects resource allocation decisions that are shaped by evolution in specific ecological niches and can be quickly adapted if necessary.

Author summary

Bacteria grow at very different rates on different substrates. Therefore, the substrates themselves are often denoted as ‘rich’ substrates versus ‘poor’ substrates. However, it remains unclear what makes a nutrient a ‘rich’ or a ‘poor’ substrate. A host of different explanations have been suggested, such as the energy content and chemical composition of the nutrient, the amount of protein required for efficient catabolism of the substrate, or limitations of ‘space’ in the plasma membrane for fitting transporters of the substrate. All of these hypotheses are actively debated today. Here, we show instead that, at least for certain substrates, nutrient quality can be a plastic property that can be dialed by evolution. Different nutrients enable bacteria to grow, but they also serve as a major signal that allows microbes to infer information about their environment. We propose that nutrient quality encoded in a combination of regulatory architecture and enzymatic properties serves both as a map of the safety and reliability of the environment and as a regulatory mechanism that implements proteome allocation decisions.

Introduction

Different growth rates of bacteria on different substrates are often attributed to the ‘quality’ of nutrients, which is closely related to carbon hierarchy in catabolite repression [1,2]. ‘Nutrient quality’ is the key parameter in bacterial growth laws [35], originally formulated by Scott et al. [3,4], which elegantly link growth rates and proteome allocation. However, the precise meaning of this growth-rate-determining parameter remains unclear.

As illustrated in Fig 1A, growth laws originate from the balance of two fluxes: First, protein synthesis flux from translating ribosomes generates the increase in total protein proportional to the growth rate λ: (1) where ϕR is the ribosomal fraction of the proteome. Note that ϕR implicitly includes the rest of protein biosynthesis pathways, including amino acid biosynthesis, whose protein cost theoretically scales identical to ribosomes. ϕ0 is a constant, and κt is a parameter denoting the translational capacity, proportional to the speed of ribosomal elongation. Second, biomass building blocks for biomass production rate (proportional to growth rate) are sustained by an equivalent nutrient influx (2) where ϕP is called the P-sector, the proteome fraction that includes ‘catabolic’ enzymes catalyzing this nutrient flux. Here, κn is the key parameter, denoting the ‘nutritional capacity’ or the ‘nutrient quality’ of the substrate. Hence, both nutrient flux and translational flux promote faster growth. However, the proteome sectors ϕ are fractions of the total proteome and are therefore constrained. In other words, the sum of these sectors cannot surpass a maximum fraction, denoted by , which could theoretically by equal to 100% of the proteome, but is found to be smaller than 100% empirically: (3)

thumbnail
Fig 1. Growth laws and protein investment.

a, Overview of bacterial growth law model. Growth rate is determined by the balance of nutrient influx catalyzed by the P-sector and nutrient consumption by protein biosynthesis catalyzed by the R-sector. The Q-sector is a constant growth-rate-independent part of the proteome. Nutrient quality κn determines how much nutrient flux is achieved for a given P-sector fraction. For higher nutrient quality (κn large), a smaller P-sector achieves sufficient nutrient flux and frees up proteome for the R-sector for more biosynthesis, resulting in a higher growth. (panel created with Biorender.) b, Protein copy number of transporters (or first metabolic enzyme) of different substrates plotted against growth rate achieved on these respective substrates. Protein copy numbers were calculated for the same slow carbon-limited growth conditions by combining proteomics [6] and ribosome profiling data sets [12].

https://doi.org/10.1371/journal.pcbi.1011735.g001

Eqs [13] can be combined to give an expression relating the growth rate, the translational capacity, and the ‘quality’ of the substrate: (4) where λmax resembles a maximum growth rate that is related to the growth-rate invariant fraction of the proteome ϕfixed via λmax = (1−ϕfixed)κt.

As neither translational capacity κt, nor maximum growth rate λmax depend on the carbon substrate, according to Eq [4], the key growth rate-determining parameter in the growth laws is the nutrient quality κn. As illustrated in Fig 2B, this parameter determines the size of the P-sector at maximum growth rate. Growth on a ‘good’ substrate (large κn) results in a small P-sector, freeing up more proteomic resources for ribosomes (R-sector) and thereby allowing faster growth (Fig 2B, right side). Conversely, a ‘poor’ substrate (small κn) results in a larger P-sector, constraining the size of the ribosomal R-sector and therefore resulting in slower growth (Fig 2B, left side).

thumbnail
Fig 2. Model summary.

a, Summary of the growth law model and parameters. ϕR is the ribosomal fraction of the proteome. ϕ0 is a constant, and κt is a parameter denoting the translational capacity. ϕP is called the P-sector, the proteome fraction that includes ‘catabolic’ enzymes catalyzing this nutrient flux. κn is the ‘nutritional capacity’ or ‘nutrient quality’. The sum of proteome sectors cannot surpass a maximum fraction, denoted by . λmax resembles a maximum growth rate that is related to the growth-rate invariant fraction of the proteome ϕfixed via λmax = (1−ϕfixedt. b, Illustration how different nutrient quality results in different growth rates. For a low nutrient quality (left), a high expression level of the P-sector is required to achieve sufficient nutrient flux. This leaves fewer resources available for the ribosomal sector ϕR and overall results in slower growth. Conversely, for high nutrient quality κn (right), a higher nutrient flux is catalyzed by a smaller P-sector, freeing up proteomic resources for higher expression of the ribosomal sector ϕR and resulting in faster growth. c, Illustration of nutrient flux and growth rate as a function of cAMP-mediated C-sector expression (left). The C-sector is one of the major components of the P-sector in the growth theory, as illustrated in the pie chart (right). Another large part being the ppGpp-activated protein sector that we denote as the S-sector, where S stands for stress. Most transporters and substrate-specific metabolic genes are part of the C-sector and increasing the C-sector increases nutrient flux (dashed lines, left top panel). Higher nutrient flux leads to an increase in growth rate, but only up to an optimum level, at which flux for biosynthesis balances nutrient flux (left bottom panel). At even higher C-sector expression, growth rate drops because there are insufficient proteomic resources for the R-sector to process nutrient flux generated by the disproportionate C-sector. Nutrient quality κn is determined by how much nutrient flux is achieved per C-sector induction. A “good” substrate results in a steep increase in nutrient flux (blue dashed line, left top panel), whereas a “poor” substrate results in a much flatter induction of nutrient flux (orange dashed line, left bottom panel). Therefore, “poor” substrates result in higher cAMP levels but slower growth rates than “good” substrates (left bottom panel). The steepness of nutrient flux induction, defining nutrient quality, is determined by the catalytic rates of substrate-specific enzymes, but also by the expression level of substrate-specific transporters and enzymes. We denote the core proteome fraction of substrate-specific transporters and enzymes by C*-sector, which is a fraction f in the much larger P-sector, f = ϕC*P. We then denote the core nutrient quality based on fundamental biochemical enzymatic properties by . The effective nutrient quality that emerges in the growth laws can be modulated by changing the expression fraction f, . We hypothesize that because the core enzyme proteome fraction ϕC* is a small fraction of the P-sector ϕP, even for the costliest substrates in terms of protein cost, nutrient quality can be dialed up or down in response to ecological needs by changing the expression fraction f. We denote the P-sector fraction that is not part of substrate-specific metabolism as the “adaptability” sector ϕAD: ϕAD = ϕP−ϕC*. Components of this sector are not important for supporting growth in the current growth conditions, but instead constitute a preparatory response. d, Illustration how the core catabolic fraction determines nutrient quality and growth rate. Within the co-regulated C-sector, low expression and a weak induction of substrate-specific enzymes will result in a lower effective nutrient quality and a higher expression level of the adaptability sector ϕAD (left). Conversely, a large core catabolic fraction results in fast growth and relatively low expression of the adaptability sector ϕAD (right). By dialing the core catabolic fraction, bacteria can convert information about their environment conveyed by the nutrient present, into resource allocation decisions determining their adaptability and preparedness for changing environments or the onset of stress. (illustrations created with Biorender).

https://doi.org/10.1371/journal.pcbi.1011735.g002

Growth laws make a set of elegant predictions that have been quantitatively validated experimentally [3,4] and form the basis of many more detailed models [1,6,7]. However, the basic question how the growth rate determining parameter nutrient quality κn should be interpreted remains. Nutrient quality could be an intrinsic carbon-specific property, such as proteome cost of the catabolic pathways, required for metabolizing the substrate and therefore determined by fundamental enzymatic properties like the maximum catalytic rate Vmax of these metabolic enzymes and transporters. In this case, growth on a ‘poor’ substrate leads to an obligatory increase in the P-sector ϕP that limits growth rate. The other possibility is that nutrient quality κn is a ‘plastic’ parameter, largely determined by the relative expression levels of carbon-specific transporters and enzymes. In this case, bacteria are programmed by evolution to grow slowly on certain substrates. Different species of bacteria and even different strains could set ‘nutrient quality’ at different levels and grow at different rates. It should then be possible to change growth rates by simply changing the expression level of carbon-specific transporters and metabolic enzymes. However, despite the simplicity of this idea, it is challenging to test experimentally because metabolic pathways of individual substrates typically involve many enzymes and transporters that need to be expressed in concert, as well as additional layers of transcriptional and allosteric regulation that affect expression levels and fluxes.

Materials and methods

Preparation of bacterial growth media

For the experiments mentioned in this study, we used single carbon and single nitrogen-source containing minimal media denoted N+C+ minimal medium. In this minimal medium, 20mM NH4Cl was maintained as single nitrogen source, and 20 mM glucose or 20mM mannose were used as carbon sources as required.

To make the minimal media, a 4X N-C- salt stock-solution was first prepared, and the recipe to prepare 1L 4X N-C- solution is shown below.

To prepare 1X N+C+ minimal media, 4X N-C- salt solution stock was diluted to 1X concentration in final media formulation and supplemented with 20 mM NH4Cl (final concentration), and the carbon source (20mM glucose or 20 mM mannose) was added separately. After preparation, each medium was filtered using disposable vacuum filtration system (Corning, PES filter with 0.22 pore size). Carbon sources were added from 1M stock solution, prepared from D (+) Mannose (Sigma Aldrich, M6020-25G, Lot #BCBV4824), and glucose (Sigma, G5146).

YCE119 construction

We replaced the chromosomal promoter driving the D-mannose PTS permease manXYZ with the promoter of the glucose PTS transporter PptsG, and deleted the mlc gene in the same strain.

For strain construction, we used NCM3722 as background.

Strain genotype: NCM3722 FRT-Kan-FRT:rrnBT:PptsG-m5’UTR-manX, Δmlc, FRT-FRT:Ptet-manA.

The mlc gene was knocked out with the mlcΔ p1 lysate [8] generated from Keio collection [9]; the kanamycin cassette was then flipped out with the one-step inactivation method [10]. For integrating the Tet promotor to manA gene chromosomally, the pkD13 plasmid was used to amplify the Ptet promotor with primers that also contain manA target region for replacing manA’s own chromosomal promotor. The PCR product was sequenced before chromosomally integrated into the target strain using pKD46 recombineering method, and the kanamycin cassette was flipped out [10]. The promotor of the glucose PTS transporter PptsG was first cloned into the pKD13 plasmid using Gibson Assembly. The Gibson Assembly product was then used as the PCR template with primers that span from ptsG promotor to manX gene for PCR amplification and included homologous regions for integration. The PCR product was sequenced before being chromosomally integrated into the target strain using pKD46 recombineering method [10] and selected with kanamycin.

Bacterial culture and growth rate measurement

A single isolated colony of WT and YCE119 was picked up from respective LB plates and inoculated in glass tubes containing 5 ml LB and grown for 2–3 Hours. For biological replicates, multiple tubes were inoculated, each from different single colonies. Inoculated LB tubes were incubated in a shaking air incubator (infros HT), set at 37°C with 200 rpm orbital shaking. After 2–3 hours of growth in LB, 1% inoculum was transferred to glass tubes containing fresh minimal media (N+C+glucose or N+C+mannose) and incubated in the shaker incubator overnight. In this way, the same colony can be propagated in two different types of media for downstream experiments. Next morning, from the overnight in minimal media, 1% inoculum was taken out and transferred to glass tubes containing 5 ml fresh minimal media (glucose or mannose respectively). Freshly inoculated tubes were placed in shaker water bath set at 37°C with 200 rpm orbital shaking.

For experiments with cyclic AMP, we added cAMP directly to the glass tubes containing fresh minimal media from a 25 mM cAMP stock solution (cyclic AMP sodium salt, Sigma, A6885), and then the tube was inoculated with desired bacterial strain.

For growth rate measurements, for each biological replicate, samples were taken for at least 4 time points during steady-state growth, and optical density (OD600) was measured using a spectrophotometer (Genesys 30 visible spectrophotometer, Thermo Scientific). Doubling time and growth rates were determined by plotting the measured OD as a function of elapsed time in a semi-log scale.

Range extension assay

Range extension assay was done following the protocol modified from Cremer et. al. [11]. To prepare soft agar plates, we mixed 0.25% agar in 1X N+C- minimal media and autoclaved the agar-containing media. Carbon sources were not added at the time of autoclaving to avoid any denaturation and caramelization of the sugars. 1M glucose and 1M mannose solution were prepared separately and sterilized by passing through a 0.22μm filter. Autoclaved soft agar and carbon stock solutions (glucose and mannose) were taken inside a sterile biosafety cabinet, and we added glucose and mannose separately to the soft agar to obtain a final glucose or mannose concentration of 20 mM. After addition of respective carbon sources, 15 ml soft agar was poured on 100mm sterile plastic petri plates and kept inside a laminar airflow biosafety cabinet to maintain sterility. Note, unlike Cremer et. al. no attractant was added to soft agar plates. For cAMP containing soft agar plates, cAMP was added to still warm soft agar (after addition of respective carbon sources) from a 25mM stock solution (cyclic AMP sodium salt, sigma), as mentioned before, to obtain a final concentration of 3.5 mM. 7 ml warm soft agar was poured on 60 mm plates inside a laminar airflow biosafety cabinet.

After pouring into plates, soft agar was allowed to cool off, and when solidified, plate lids were placed, and plates were taken out of the biosafety cabinet. Due to the soft consistency of agar, these plates cannot be inverted.

For range extension assay, WT (NCM3722) and YCE119 cultures were grown in minimal media with 20mM glucose, 20mM mannose, and 20mM mannose with 3.5 mM cAMP, following the protocol mentioned above. When respective cultures were still in the exponential growth phase and reaching ~0.4–0.5 OD, 2 μl of culture was taken out of each exponentially growing tube and gently placed on the pre-warmed (to 37°C) soft agar plates as a droplet. On each plate, 4 droplets were placed equidistant from one another. Plates were carefully placed in a 37°C air incubator and were not inverted. After overnight incubation, plates were taken out and imaged in a custom-built plate imaging set-up. After imaging, plates were placed back in the incubator and re-imaged on day 3. We kept the incubator humidified by keeping an open water pan to ensure plates did not dry out during long incubation periods. Water droplets inside the lid were removed by soaking with sterile absorbent paper during the incubation period.

For colony size measurement, plate images were analyzed in FIJI. Each colony image was binarized by applying a threshold, and the colony area was calculated from the binary image.

Starvation assay

WT (NCM3722) cells were grown exponentially in minimal media containing 20mM glucose or 20 mM mannose, and YCE119 cells were grown in minimal media containing 20 mM mannose or 20 mM mannose and 3.5 mM cAMP. When cells were still in the exponential phase (0.2–0.4 OD), an aliquot from the exponentially growing culture was taken out for the starvation assay. To keep the cell number roughly similar at the beginning of starvation assay, volume of sample taken was adjusted according to optical density. Cells were centrifuged and pelleted down, and resuspended in carbon free minimal media (N+C-). This step was repeated twice to remove traces of carbon sources. Cells were finally resuspended in 5 ml N+C- starvation media, lacking the carbon source, in glass tubes and placed in a shaking air incubator set at 37°C at 200 rpm. Immediately after resuspension of the cells in carbon free medium at the beginning of the starvation assay, samples were taken out of each tube and diluted 100,000-fold in carbon free medium. 100μl of the diluted samples from each tube were plated on LB-agar plates containing tetrazolium chloride. Metabolically active live cells give red-colored colonies on these plates. After plating, plates were incubated in 37°C air incubator overnight and imaged next morning using a custom-built plate imaging setup. Starving cultures in carbon free media were kept in the shaking incubator for a week, and at the end of experiment, the number of surviving cells were quantified by plating on tetrazolium chloride containing LB plates. On the final day of the experiment, cells were diluted 20,000-fold in carbon free media and plated. Plates were incubated in a 37°C incubator overnight and imaged next morning using our custom-built setup. The number of colonies on each plate was counted using Cell Profiler (open-source image processing software) and normalized to colony number of the same condition in the beginning of the experiment, from which survival percentages were calculated.

Diauxic shift lag time measurement

WT (NCM 3722) and YCE119 strains were grown in LB culture from single isolated colonies for 2–3 hours in a shaking (200 rpm) air incubator set at 37°C. For biological replicates, multiple isolated colonies (n = 3 for NCM 3722 and n = 3 for YCE119) were inoculated in LB. After growing for a few hours in LB, 1% inoculum was transferred to minimal medium containing 1% (approximately 555μM glucose or mannose). WT strain was inoculated in N+C+555μM glucose medium, and YCE119 was inoculated in N+C+555μM mannose medium. Inoculated tubes were placed in a shaking air incubator (Infors HT) maintained at 37°C with 200rpm orbital shaking. The next day, 1% of the overnight cultures were freshly inoculated separately in 5ml fresh N+C+ minimal media with 555 μM glucose (WT) or 555 μM mannose (YCE119). From the overnight culture, the YCE119 strain was also inoculated in minimal media with 555uM mannose and supplemented with 3.5 mM cyclic AMP (added from 25mM aqueous stock solution). To measure WT (NCM3722) in mannose, we inoculated the 1% inoculum from the LB cultures (same 3 biological replicate colonies) to 20mM mannose containing minimal medium early in the morning and grown in a shaking air incubator. When these cultures reached exponential phase of growth and ~0.5 OD, 1% inoculum was transferred to minimal medium with 555 μM mannose. Lag phase measurement for WT in mannose was done on a separate plate.

To measure lag phase and growth in acetate after diauxic shift, 5mM (final concentration) sodium acetate were added to each media condition mentioned above and mixed well by vortex mixing. Freshly inoculated tubes were grown for 2–3 hours in an air incubator, and after that 200 μl samples were transferred to a clear bottom 96 well plate with dark wall (Greiner bio-one, 655090). For each biological replicate, multiple wells of the 96 well plate were used. OD600 from each well was measured over time in a BioTek Synergy H1 microplate reader, maintained at 37°C with continuous double-orbital shaking. Background subtracted averages of the OD600 from replicate wells for each biological replicate were calculated first, and then the average of all biological replicates with standard deviation was used for plotting. To quantify lag time, each condition was normalized to the point when the exponential growth slowed down for each strain, indicating the depletion of glucose and initiation of diauxic shift. Time of normalization for each strain was denoted as time 0, and normalized fluorescence intensity as a function of time was plotted in a semi-log scale.

Calculation of P sector allocation and protein cost

To estimate protein copy numbers of the main carbon transporting enzyme (or the first enzyme in the primary carbon degradation pathway) at slow growth rate (0.45), where the P-sector proteins are highly expressed, we combined protein copy numbers determined via ribosome profiling for E. coli grown in glucose minimal medium by Li et al. [12] (see S1 Table) with fold-changes between different growth rates determined by 15N labeling and mass spectrometry in Hui et al. [6] (see S2 Table) for those proteins that showed up in both datasets. Growth rate in respective substrates (acetate, mannose, lactose, glycerol, and pyruvate) taken from You et. al. [1], except for growth rate of the wildtype strain on mannitol, which needed to be measured in our lab.

Results

First, we wanted to determine to what extent nutrient quality of different substrates was reflected in proteome allocation to substrate-specific transporters and enzymes. If nutrient quality were an intrinsic cost of carbon catabolism, poorer substrates would require a larger proteome fraction allocated to them during growth on these substrates. However, there is currently no dataset available that allows comparison of absolute protein abundances of substrate-specific catabolic proteins on many different substrates. To address this question, we wanted to compare protein copy numbers of different P-sector proteins at the same growth rate. Specifically, we were interested in transporters and metabolic enzymes catabolizing different substrates and their expression level at slow growth rates when the P-sctor is fully activated. Therefore, we combined absolute protein copy numbers determined via ribosome profiling for E. coli grown in glucose minimal medium by Li et al. [12] with relative fold-changes between different growth rates determined by 15N labeling and mass spectrometry in Hui et al. [6]. Strikingly, when we plotted the copy number of carbon transporters (alternatively, the first metabolic enzyme in the degradation pathway) versus the growth rate of E. coli on these carbon sources, we discovered a rank-ordered anticorrelation (Fig 1B). Transporters and metabolic enzymes related to poorer substrates were expressed at higher copy numbers than those of better carbon sources.

Although there is limited coverage in the proteomics dataset, and we only analyzed the transporters or first enzymes in the substrate utilization pathways, this finding suggests a relationship between nutrient quality and protein investment. In fact, if nutrient quality were directly caused by differences in the maximum catalytic rate of enzymes and transporters (Vmax), such an anticorrelation between protein cost and growth rate would follow from the growth laws during growth on these substrates because growth on poor substrates would require more enzymes. The observed anticorrelation is still surprising as these substrates were not present in the experimental growth conditions from which protein abundances were determined. Hence, growth rate on different substrates appears to be reflected in expression levels of transporters and metabolic enzymes of these substrates. We found a similar relation to hold when we converted copy numbers to protein cost by considering the size of the individual proteins (S1 Fig).

Nevertheless, we were skeptical if growth rates on poor substrates were directly limited by protein cost or available capacity on the plasma membrane because proteome fractions dedicated to substrate-specific transporters and enzymes are relatively small, and membrane capacity could be increased by changing cell shape. We therefore wanted to test directly whether growth rates of E. coli on mannose are directly limited by the cost of producing the necessary proteins for substrate transport and metabolism, or if they are limited in some other way by the available capacity for substrate transport. We decided to focus on the substrate mannose, which is one of the slowest carbon sources and exhibits the second highest apparent protein cost in the set of substrates in Fig 1B. Whereas growth on acetate requires complete rerouting of central metabolism and reversal of flux in glycolysis to a gluconeogenic mode [13], mannose is a glycolytic substrate that requires only the expression of only a handful of dedicated transporters and enzymes. Our strategy was to try to simultaneously alleviate all potential bottlenecks in mannose catabolism to achieve a faster growth rate. To change expression levels while avoiding the burden of plasmid expression, which can affect growth rates, we decided to replace the chromosomal promoter driving the D-mannose phosphotransferase system (PTS) permease manXYZ with the promoter of the glucose PTS transporter PptsG because glucose is one of the best substrates of E. coli. To prevent carbon-specific repression of this promoter due to the lack of glucose, we also deleted the mlc gene. Mlc is a transcription factor that represses the expression of the PTS system dependent on the presence of glucose [14,15]. Finally, to prevent a metabolic bottleneck and to ensure that mannose flux enters glycolysis efficiently, we replaced the chromosomal promoter of the mannose-6-phosphate isomerase gene manA gene with the Tet-promoter. Without expression of TetR, this promoter acts like a strong constitutive promoter [16]. The resulting genetic systems in strain denoted YCE119 is illustrated in Fig 3A and will be referred to as the ‘swapped promoter strain’.

thumbnail
Fig 3. Plasticity of nutrient quality.

a, We replaced chromosomal promoters of transporters and metabolic enzymes required for mannose metabolism with the promoter of glucose transporter (P-ptsG). To prevent repression due to lack of glucose, we knocked out the glucose-specific transcriptional regulator Mlc. Finally, to ensure that carbon flux from processed mannose enters the glycolysis pathway, we placed mannose-6-phosphate isomerase under a strong constitutive promoter (P-tet). We refer to this strain as the swapped promoter strain (YCE119). (construct map illustration created with Biorender.) b, Growth rates of wildtype and swapped promoter strain (YCE119) on glucose and mannose. Mannose is one of the slowest substrates for the wildtype (orange circles), and glucose is often considered the best substrate of E. coli with the fastest growth rates in minimal medium (blue circles). Yet the swapped promoter strain grows on mannose as fast as the wildtype on glucose (orange triangles). Therefore, the genetic modifications in the swapped promoter strain have changed the nutrient quality of mannose into the nutrient quality of glucose. This means that nutrient quality is not limited by fundamental biochemical constraints. Unpaired t-test was performed, and following P-values were obtained. WT Glucose:WT mannose P-value<0.0001. WT Mannose:YCE 119 Mannose P-value<0.0001, and WT Glucose:YCE119 Mannose: non-significant (ns) c, Data from Towbin et al. [17], who titrated cAMP levels on different substrates. “Good” substrates have a higher peak growth rate at lower levels of cAMP. d, An inverse correlation of endogenous Crp activity on different substrates with growth rate on the respective substrate from Towbin et al. [17]. Poor substrates require higher cAMP levels for maximum growth. A similar relationship is found when plotting CRP activity vs. growth rate for maximum growth rate (S3 Fig).

https://doi.org/10.1371/journal.pcbi.1011735.g003

Next, we measured growth rates of swapped promoter strain as compared to the wildtype strain in batch culture. If growth rate were truly limited by the protein cost of mannose-specific metabolism or the uptake rate that could be accomplished due to limited capacity of the plasma membrane, overexpression of mannose transporters and metabolic enzymes in the swapped promoter strain should have no positive effect on growth rate. In fact, according to the growth laws, overexpressing too much of the mannose-specific pathway should result in decreased growth rates by resulting in a misbalance of the P-sector and R-sector. Instead, we found that the swapped promoter strain grew on mannose minimal medium at the same rate as the wildtype strain on glucose minimal medium (Fig 3B). While we refer to this strain as the swapped promoter strain for simplicity due to the exchanged promoter of the manXYZ operon, in this strain, the promoter of manA was also exchanged by a strong constitutive promoter (Ptet), and the gene encoding the transcription factor Mlc was knocked out. The effect of subsets of these modifications on growth rate on mannose is shown in S2 Fig. Notably, none of the individual modifications were sufficient to confer the full growth rate advantage. We also found that the strong constitutive promoter (Ptet) driving the manXYZ operon conferred the same growth rate advantage for growth on mannose as the PptsG promoter. These phenotypes vividly demonstrate that mannose is not an inherently ‘poor’ substrate and that the slow wildtype growth rate on mannose is not the result of any fundamental constraint or limitation. Nutrient quality, at least for the case of mannose, appears to be a plastic property, rather than a fundamental biochemical or enzymatic limitation. Indeed, this result is consistent with previous observations about suboptimal growth rates of E. coli on glycerol and fructose [17,18], as well as reports about suboptimal growth rates due to endogenous levels of cAMP [17]. It is also consistent with the observation that substrate preferences can be completely reversed in other bacterial species like Pseudomonas aeruginosa, which grow faster on supposedly ‘poor’ gluconeogenic substrates like acetate than on supposedly ‘rich’ glycolytic substrates [13], despite the high degree of conservation of central metabolism.

How should nutrient quality on different substrates then be interpreted? We propose that different substrates convey information about the safety and reliability of the current growth environment in the context of the ecological niche via their ‘nutrient quality’. Many carbon-specific transporters, metabolic enzymes, and other proteins that are unrelated to the current growth conditions are co-regulated in the cAMP-activated C-sector [1] that is a sub-sector of the P-sector (Fig 2C, red sector). In fact, it is long known that E. coli progressively activates more genes as the ’quality’ of carbon source declines, including uptake systems for alternative carbon sources [19]. Moreover, the P-sector includes other preparatory pathways for the response to stress and adverse conditions that are activated by ppGpp, and that we denote the S(tress)-sector (Fig 2C, green sector). As illustrated in Fig 2C and 2D, although being co-regulated, the core C-sector constituting the substrate-specific transporters and enzymes enabling growth on the current growth conditions (denoted the C*-sector) may constitute only a small fraction of the overall C-sector. According to the growth laws, the entire P-sector is upregulated up to a level, where nutrient flux catalyzed by the core C*-sector balances the flux of biomass precursors consumed by ribosomes for protein translation via the R-sector. Both too low and too high levels of induction of the P-sector result in slower growth rates as proteome allocation is unbalanced (Fig 2C, left panels), which has been vividly demonstrated experimentally by Towbin et al. who titrated cAMP concentrations [17] (Fig 3C). The cAMP level where the maximum growth rate is achieved depends on the nutrient quality, as can be seen by replotting the maximum growth rate as a function of the corresponding cAMP level from the data by Towbin et al. [17] (Figs 2D and S3).

‘Poor’ substrates, as compared to ‘rich’ substrates, require a higher induction of the P-sector to achieve the same nutrient flux and therefore have higher levels of cAMP (Fig 2C, left panel). However, while this flux depends on the catalytic rates of the substrate-specific enzymes and transporters represented by the core nutrient quality , their total expression fraction f within the much larger P-sector is just as important and results in the growth-rate determining nutrient quality given by , with f≤1. Hence, most of the growth-rate-limiting protein cost comes from P-sector components unrelated to current growth conditions (Fig 2D). This is consistent with data from proteomics [6], as well as genome scale metabolic models [20]. In this case, by tuning the expression level of substrate-specific transporters and pathways within the P-sector, nutrient quality of a specific substrate can be freely tuned. This effect is vividly illustrated by the dramatic increase in growth rate on mannose in the swapped promoter strain (YCE 119).

So why would bacteria want to decrease nutrient quality of a given substrate by reducing expression levels of transporters and metabolic enzymes within the P-sector? As illustrated in Fig 4A, lower nutrient quality results in higher levels of the adaptability sector (Fig 2C and 2D). The adaptability sector may not only include many alternative transporters and metabolic enzymes for substrates that are not present in the current medium, but also proteins related to starvation survival, and motility etc. According to this interpretation, ‘nutrient quality’ reflects a resource allocation decision of fast growth versus preparation for changing growth conditions that is consistent with trade-offs between growth rate and adaptability that have recently been characterized [21,22].

thumbnail
Fig 4. Lag and starvation cost of higher nutrient quality.

a, Illustration of resource allocation for high nutrient quality (left) and low nutrient quality (right). Low nutrient quality reflects high expression levels of the adaptability sector ϕAd. (pie chart illustrations created with Biorender.) b, Relative quantification of proteins that have been shown to be related to lag phase in Basan et al. [22]. Proteomics data from Hui et al [6]. (color bar illustration created with Biorender.) c, Relative quantification of proteins that have been shown to be related to improved starvation survival by Schink et al. [23]. Proteomics data from Hui et al [6]. d, Diauxic shifts from glucose or mannose to acetate. The wildtype strain exhibits a lag phase of several hours from glucose to acetate (blue circles), but almost no lag phase from mannose to acetate (orange circles). The swapped promoter strain (YCE 119) grows at a similar growth rate to glucose on mannose (Fig 3B) but has lost the ability to quickly switch to acetate. Upregulation of the C-sector by addition of cAMP (3.5mM) to the growth medium results in a much shorter lag time, however, only at the cost of slower growth (S4 Fig), as predicted by the model (Fig 2). Three biological replicates were used for all conditions. Mean values of 3 biological replicates is plotted, and error bars represent standard deviation. e, Starvation survival after 7 days of starvation relative to initial CFUs. The wildtype survived carbon starvation much better after growing on mannose (bottom bar) than on glucose (top bar). The faster-growing swapped promoter strain (YCE 119) lost this improved survival on mannose (2nd bar from top). With addition of cAMP (3.5mM) to the mannose growth medium but not to the starvation medium, the swapped promoter strain exhibited improved starvation survival (3rd bard from top) at the cost of a slower growth rate (S4 Fig). Three biological replicates were used for YCE119 grown in N+C+Mannose, WT grown in N+C+Glucose, and 4 biological replicated used for YCE 119 grown in N+C+Mannose with 3.5mM cAMP and WT grown in N+C+Mannose. Every biological replicate was plated in triplicate. Mean survival frequency plotted, and error bars represent standard deviation. Unpaired t-test performed, and following P-values obtained: WT Glucose-YCE119 Mannose: non-significant, WT Glucose-YCE119 Mannose+3.5mM cAMP:0.0054, WT Glucose-WT mannose<0.0001.

https://doi.org/10.1371/journal.pcbi.1011735.g004

Indeed, using the dataset by Hui et al. [6], we found that with decreasing nutrient quality abundances of proteins shown to be involved in minimizing lag times [22], as well as proteins shown to help starvation survival [23] strongly increased as shown in Figs 4B and 3C respectively. These changes in protein abundances reflect a larger allocation to the adaptability sector ϕAD (Fig 4A). We wanted to quantify if these changes in protein abundances were also reflected in phenotypes directly, and specifically, if the swapped promoter strain exhibited a loss of capabilities as compared to the wildtype strain growing on mannose, reflecting the cost of faster growth. We first tested the ability to adapt to growth on different substrates as quantified by lag time. One of the longest, but also most frequently encountered lag times comes from the shift from glycolytic substrates to the primary fermentation product acetate [21,22]. The wildtype growing on glucose has a much longer lag time in a diauxic shift to acetate than when growing on mannose (Fig 4C, blue circles vs. orange circles). However, as predicted, the swapped promoter strain growing on mannose lost this ability to quickly shift to acetate and exhibited a long lag phase (Fig 4D, orange triangles), comparable to lag time of the wildtype on glucose.

According to the model outlined in Fig 2, it should be possible to reduce lag times by artificially upregulating the C-sector, that is part of the larger adaptability sector, by increasing the abundances of lag related proteins (Fig 4B). Therefore, we measured lag time in a diauxic shift to acetate of the swapped promoter strain with the addition of cAMP to the growth medium, which activates the C-sector. Indeed, the addition of cAMP to the mannose medium resulted in substantially reduced lag of the swapped promoter strain (Fig 4D, orange inverted triangles). However, this improvement in lag times came at the cost of a slower growth rate in the presence of cAMP (see S4 Fig).

Next, we tested the ability of the swapped promoter strain and the wildtype to survive periods of starvation. Bacteria spend a large fraction of their time coping with periods of starvation and therefore adapting to and surviving starvation is likely a key contributor to bacterial fitness in natural environments. After growing bacteria on different substrates, we washed them and resuspended them in carbon-free medium. We then quantified survival by plating and counting colony-forming units at different timepoints during starvation. As anticipated from changes in the abundance of starvation-related proteins (Fig 4C), after growing more slowly on mannose minimal medium, the wildtype strain survived substantially better than after faster growth on glucose (Fig 4E, 1st bar vs. 4th bar). However, consistent with our expectation, this improved starvation survival ability after growth on mannose was completely lost in the faster-growing swapped promoter strain (Fig 4E, 2nd bar from the top). According to the model outlined in Fig 2, we hoped to rescue the improved starvation survival of the swapped promoter strain by inducing higher expression of the AD-sector with the addition of cAMP to the pre-starvation medium. Indeed, the addition of cAMP to the pre-starvation medium improved survival (Fig 4E, 3rd bar from the top). However, this improved starvation survival came at the cost of a reduced growth rate on mannose minimal medium with cAMP (S4 Fig). According to our model, this slower growth rate is due to the protein expression burden of the increased AD-sector (Fig 2C).

Finally, we investigated cell motility, which is another adaptability phenotype that serves the purpose of escaping the current environment using different motility metrics on soft agar plates (Fig 5A). Motility also requires a substantial investment of proteomic resources in flagellar proteins [6]. The abundances of motility proteins increase with decreasing nutrient quality based on the dataset by Hui et al. [6], as shown in Fig 5B. Indeed, we found that wildtype bacteria growing on mannose were substantially more motile than wildtype bacteria growing on glucose as quantified by the increased fraction of colonies that showed large scale swarming (Fig 5C, blue circles vs. orange circles), as well as by an increase in colony range extension (Fig 4D, blue circles vs. orange circles. In contrast, the swapped promoter strain lost this increased motility when growing on mannose according to both metrics (Figs 5C and 4D, orange triangles). Like the wildtype strain on glucose, the swapped promoter strain on mannose formed small colonies with smooth edges (Fig 5A), and no colonies exhibited swarming (Fig 5C). However, when we artificially induced higher expression of AD-sector with addition of cAMP to the plates, the swapped promoter strain regained high motility (Figs 5C and 4D, inverted triangles). As demonstrated in batch culture, this higher motility comes at the cost of a reduced growth rate of the swapped promoter strain on mannose (S4 Fig). Hence, microbial motility closely matches the nutrient quality dependence of other adaptability phenotypes (Fig 4). We note that while Figs 3 and 4 show that gene expression patterns and phenotypes are consistent with our proposed hypothesis, they should not be considered as causal evidence for the slow growth rate on mannose.

thumbnail
Fig 5. Motility cost of higher nutrient quality.

a, Motility assays on soft agar plates. Colonies on fast growth medium are smaller and show smooth edges (top panel). Colonies on slow growth medium are larger with ruffled edges (middle panel). A fraction of colonies on slow growth medium shows a longer-range swarming phenotype (bottom panel). (pie chart illustrations created with Biorender.) b, Relative quantification of proteins related to motility. Proteomics data from Hui et al. [6]. (color bar illustration created with Biorender.) c, Fraction of colonies exhibiting swarming on soft agar. For wildtype on glucose (blue circle), we observed no swarming, whereas on mannose a substantial number of colonies showed swarming (orange circles). The swapped promoter strain (YCE 119) lost the ability to swarm on mannose (orange triangles), but addition of cAMP (3.5mM) resulted in a large fraction of swarming colonies. 24 colonies were analyzed for WT grown in Glucose, 72 colonies were analyzed for WT grown in Mannose, 24 colonies were analyzed for YCE 119 grown in Glucose, 48 colonies were analyzed for YCE119 grown in mannose, and 24 colonies were analyzed for YSE 119 grown in Mannose with 3.5mM cAMP. d, We tracked growth of colony size over 3 days on minimal medium soft agar plates. Colony sizes of the wildtype increased much more on mannose (orange circles) than on glucose (blue circles) and exhibited rough irregular edges indicative of higher motility. The swapped promoter strain lost this enhanced motility on mannose and only showed an increase in colony size similar to the wildtype on glucose (yellow triangles). However, with the addition of cAMP (3.5mM) to mannose medium, the swapped promoter recovered high motility (yellow inverted triangles).24 colonies were analyzed for WT grown in Glucose, 72 colonies were analyzed for WT grown in Mannose, 24 colonies were analyzed for YCE 119 grown in Glucose, 48 colonies were analyzed for YCE119 grown in Mannose, and 24 colonies were analyzed for YSE 119 grown in Mannose with 3.5mM cAMP. Normalized mean area fold change plotted with standard deviation).

https://doi.org/10.1371/journal.pcbi.1011735.g005

Discussion

Intriguingly, despite ecological complexity, global cellular resource allocation can be quantitatively described by surprisingly simple bacterial growth laws. We hypothesize that the simplicity of the growth laws is in fact a consequence of an elegant regulatory strategy, in which nutrient quality constitutes a convenient regulatory lever controlling resource allocation in specific ecological environments. As illustrated in Fig 6, substrates serve as signals that microbes use to identify the environment in which they find themselves, and substrate quality is the regulatory mediator of resource allocation strategies. Rather than reflecting fundamental biochemical or biophysical limitations like protein cost or limitations in capacity for transporters in the membrane, nutrient quality reflects a probabilistic map of the safety, reliability, and profitability of different ecological environments shaped by evolution. Low nutrient quality substrates thus indicate environments where the carbon source may quickly be depleted, making investment in alternative pathways worthwhile, or environments that tend to co-occur or precede other forms of stress that the bacterium should prepare for. While nutrient quality is correlated with the efficiency of transporters and metabolic enzymes reflected in protein investment, it is not fundamentally limited by these efficiencies in physiological conditions, at least for certain substrates. In this work, we demonstrated this for mannose, the ‘poorest’ substrate of E. coli that we know. We also collected prior evidence from the literature, demonstrating substantial growth rate improvements on the intermediate quality substrates, glycerol, and fructose, with resulting growth rates similar to those on glucose (S5 Fig). While these are isolated examples, and it is impossible to exhaustively demonstrate the plasticity of nutrient quality for all substrates of E. coli, these data suggest that sub-maximum growth on mannose is far from the exception and much more common than commonly appreciated. We propose here that sub-maximum growth rates are actually the norm and that growth rates in most growth conditions can likely be improved substantially.

thumbnail
Fig 6. Nutrient quality as a map of environments and mediator of resource allocation strategies.

Ecological niches in which bacteria like E. coli evolved consist of many strikingly different environments, each coming with their own set of payoffs and hazards. The trajectory in which the bacterium transitions between these environments is partly governed by its natural life cycle, but also has a probabilistic component. To maximize their fitness, bacteria need to identify the environment in which they find themselves and adjust proteome resource allocation strategies accordingly. Nutrients allow bacteria to grow, but they also serve as a major signal that allows microbes to infer information about their environment. Nutrient quality encoded in regulatory architecture and enzymatic properties that were shaped by evolution, serves both as a map of the safety and reliability of the environment and as a regulatory mechanism implementing proteome allocation decisions. The conservation of this simple regulatory architecture across environments is what gives rise to the striking bacterial growth laws and their elegant predictions. (Illustration created with Biorender).

https://doi.org/10.1371/journal.pcbi.1011735.g006

We argue that nutrient quality constitutes a dial that can be quickly adapted by evolution in response to changing ecological conditions. This is consistent with experimental observations that adaptive laboratory evolution can quickly lead to more rapid growth even on ‘good’ carbon sources like glucose [24] and that single point mutations can sometimes lead to substantially enhanced growth rates [25]. Crucially adaptation to changing environments by modifying nutrient quality does not disrupt finely tuned cellular programs coordinating different parts of metabolism. Thereby, this regulatory architecture allows bacteria to be highly adaptable to changing environments yet retain important optimized regulatory programs coordinating central metabolism, adaptability, stress response, and ribosomal flux. Such an integration of substrate-specific and global regulatory systems as a modulator of risk-dependent gene expression has been previously suggested [19]. Here, however, we argue that rather than being implemented as a consequence of ‘energy’-rich or ‘energy’-poor substrates, this strategy may be the central determinant of substrate quality itself.

Supporting information

S1 Fig. Protein cost of substrate-specific transporters and enzymes.

Protein copy number can be converted to protein cost in units of number of amino acids by multiplying copy number with the number of amino acids in each protein. A similar inverse correlation as the one observed for copy number in Fig 1B of the main text, also holds for protein cost.

https://doi.org/10.1371/journal.pcbi.1011735.s001

(DOCX)

S2 Fig. Growth rate effect of individual genetic modifications and combinations.

Growth rates on glucose and mannose minimal medium of different strains. Data points are biological repeats. Unpaired t-test with welch correction was performed. (please see S3 Table for detailed description of strains).

https://doi.org/10.1371/journal.pcbi.1011735.s002

(DOCX)

S3 Fig. Crp activity at maximum growth rates for cAMP titrations.

In Fig 2E of the main text, we plotted endogenous CRP activity measured by Towbin et al. [17] against growth rate on the respective substrate. Endogenous CRP activity and CRP activity where growth rate is maximum are not always the same, as shown by Towbin et al. [17]. Therefore, we plot CRP activity at maximum growth rate against maximum growth rate from Towbin et al. [17] here and find a similar relationship.

https://doi.org/10.1371/journal.pcbi.1011735.s003

(DOCX)

S4 Fig. Growth rate of swapped promoter strain with the addition of cAMP.

As shown in the main text, many adaptability phenotypes that are lost in the swapped promoter strain (YCE 119) can be rescued by the addition of cAMP (3.5mM) to the growth medium, which upregulates the AD-sector. However, these improvements come at the cost of a reduced growth rate due to the substantial protein cost of the AD-sector.

https://doi.org/10.1371/journal.pcbi.1011735.s004

(DOCX)

S5 Fig. Growth rate improvements on other substrates.

Left panel: Point mutations in glycerol kinase glpK result in faster growth rates on glycerol [2527]. Resulting growth rates are effectively identical to growth rates of the wildtype on glucose. Data was replotted from Basan et al. [7]. Note that these experiments were performed in MOPS buffered medium, that results in somewhat faster growth rates in all conditions. Right panel: Knockout of the regulator of gluconeogenesis Cra results in faster growth rates on glucose and mannose. Growth rate of the Cra knockout strain on glucose is effectively identical to the growth rate of the wildtype on glucose. Data was replotted from Basan et al. [22].

https://doi.org/10.1371/journal.pcbi.1011735.s005

(DOCX)

S1 Table. Protein copy number calculation.

Protein copy number calculation of the main carbon transporting enzyme (or the first enzyme in the primary carbon degradation pathway) from Li at. al. [12].

https://doi.org/10.1371/journal.pcbi.1011735.s006

(DOCX)

S2 Table. Protein copy number fold changes across growth rates.

Fold change calculation of the main carbon transporting enzyme (or the first enzyme in the primary carbon degradation pathway) at slow growth rate (0.45) from Hui. et. al. [6].

https://doi.org/10.1371/journal.pcbi.1011735.s007

(DOCX)

S3 Table. Detailed description of the strains used in S2 Fig.

https://doi.org/10.1371/journal.pcbi.1011735.s008

(DOCX)

References

  1. 1. You C, Okano H, Hui S, Zhang Z, Kim M, Gunderson CW, et al. Coordination of bacterial proteome with metabolism by cyclic AMP signalling. Nature. 2013;500: 301–6. pmid:23925119
  2. 2. New AM, Cerulus B, Govers SK, Perez-Samper G, Zhu B, Boogmans S, et al. Different levels of catabolite repression optimize growth in stable and variable environments. PLoS Biol. 2014;12: e1001764. pmid:24453942
  3. 3. Scott M, Gunderson CW, Mateescu EM, Zhang Z, Hwa T. Interdependence of cell growth and gene expression: Origins and consequences. Science (1979). 2010;330: 1099–1102. pmid:21097934
  4. 4. Scott M, Hwa T. Bacterial growth laws and their applications. Curr Opin Biotechnol. 2011;22: 559–65. pmid:21592775
  5. 5. Scott M, Klumpp S, Mateescu EM, Hwa T. Emergence of robust growth laws from optimal regulation of ribosome synthesis. Mol Syst Biol. 2014;10: 747. pmid:25149558
  6. 6. Hui S, Silverman JM, Chen SS, Erickson DW, Basan M, Wang J, et al. Quantitative proteomic analysis reveals a simple strategy of global resource allocation in bacteria. Mol Syst Biol. 2015;11. pmid:25678603
  7. 7. Basan M, Hui S, Okano H, Zhang Z, Shen Y, Williamson JR, et al. Overflow metabolism in Escherichia coli results from efficient proteome allocation. Nature. 2015;528. pmid:26632588
  8. 8. Thomason LC, Costantino N, Court DL. E. coli Genome Manipulation by P1 Transduction. Curr Protoc Mol Biol. 2007;79: 1.17.1–1.17.8. pmid:18265391
  9. 9. Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2: 2006.0008. pmid:16738554
  10. 10. Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A. 2000;97: 6640–5. pmid:10829079
  11. 11. Cremer J, Honda T, Tang Y, Wong-Ng J, Vergassola M, Hwa T. Chemotaxis as a navigation strategy to boost range expansion. Nature 2019 575:7784. 2019;575: 658–663. pmid:31695195
  12. 12. Li G-W, Burkhardt D, Gross C, Weissman JS. Quantifying Absolute Protein Synthesis Rates Reveals Principles Underlying Allocation of Cellular Resources. Cell. 2014;157: 624–635. pmid:24766808
  13. 13. Schink SJ, Christodoulou D, Mukherjee A, Athaide E, Brunner V, Fuhrer T, et al. Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing. Mol Syst Biol. 2022;18: e10704. pmid:34994048
  14. 14. Decker K, Plumbridge J, Boos W. Negative transcriptional regulation of a positive regulator: the expression of malT, encoding the transcriptional activator of the maltose regulon of Escherichia coli, is negatively controlled by Mlc. Mol Microbiol. 1998;27: 381–90. pmid:9484893
  15. 15. Plumbridge J. Regulation of gene expression in the PTS in Escherichia coli: the role and interactions of Mlc. Curr Opin Microbiol. 2002;5: 187–93. pmid:11934616
  16. 16. Klumpp S, Zhang Z, Hwa T. Growth rate-dependent global effects on gene expression in bacteria. Cell. 2009;139: 1366–75. pmid:20064380
  17. 17. Towbin BD, Korem Y, Bren A, Doron S, Sorek R, Alon U. Optimality and sub-optimality in a bacterial growth law. Nat Commun. 2017;8: 14123. pmid:28102224
  18. 18. Basan M, Hui S, Williamson JR. ArcA overexpression induces fermentation and results in enhanced growth rates of E. coli. Sci Rep. 2017;7: 11866. pmid:28928483
  19. 19. Liu M, Durfee T, Cabrera JE, Zhao K, Jin DJ, Blattner FR. Global transcriptional programs reveal a carbon source foraging strategy by Escherichia coli. Journal of Biological Chemistry. 2005;280: 15921–15927. pmid:15705577
  20. 20. O’Brien EJ, Lerman JA, Chang RL, Hyduke DR, Palsson BO. Genome-scale models of metabolism and gene expression extend and refine growth phenotype prediction. Mol Syst Biol. 2014;9: 693–693. pmid:24084808
  21. 21. Basan M. Resource allocation and metabolism: the search for governing principles. Curr Opin Microbiol. 2018;45. pmid:29544124
  22. 22. Basan M, Honda T, Christodoulou D, Hörl M, Chang Y-F, Leoncini E, et al. A universal trade-off between growth and lag in fluctuating environments. Nature. 2020. pmid:32669712
  23. 23. Schink S, Ammar C, Chang Y-F, Zimmer R, Basan M. Analysis of proteome adaptation reveals a key role of the bacterial envelope in starvation survival. Mol Syst Biol. 2022;18: e11160. pmid:36479616
  24. 24. LaCroix RA, Sandberg TE, O’Brien EJ, Utrilla J, Ebrahim A, Guzman GI, et al. Use of Adaptive Laboratory Evolution To Discover Key Mutations Enabling Rapid Growth of Escherichia coli K-12 MG1655 on Glucose Minimal Medium. Appl Environ Microbiol. 2015;81: 17–30. pmid:25304508
  25. 25. Applebee MK, Joyce AR, Conrad TM, Pettigrew DW, Palsson BØ. Functional and Metabolic Effects of Adaptive Glycerol Kinase (GLPK) Mutants in Escherichia coli. Journal of Biological Chemistry. 2011;286: 23150–23159. pmid:21550976
  26. 26. Zwaig N, Lin ECC. Feedback Inhibition of Glycerol Kinase, a Catabolic Enzyme in Escherichia coli. Science (1979). 1966;153: 755–757. pmid:5328677
  27. 27. Pettigrew DW, Liu WZ, Holmes C, Meadow ND, Roseman S. A single amino acid change in Escherichia coli glycerol kinase abolishes glucose control of glycerol utilization in vivo. J Bacteriol. 1996;178: 2846–2852. pmid:8631672