Cellular perception of growth rate and the mechanistic origin of bacterial growth law.

Many cellular activities in bacteria are organized according to their growth rate. The notion that ppGpp measures the cell’s growth rate is well accepted in the field of bacterial physiology. However, despite decades of interrogation and the identification of multiple molecular interactions that connects ppGpp to some aspects of cell growth, we lack a system-level, quantitative picture of how this alleged “measurement” is performed. Through quantitative experiments, we show that the ppGpp pool responds inversely to the rate of translational elongation in Escherichia coli. Together with its roles in inhibiting ribosome biogenesis and activity, ppGpp closes a key regulatory circuit that enables the cell to perceive and control the rate of its growth across conditions. The celebrated linear growth law relating the ribosome content and growth rate emerges as a consequence of keeping a supply of ribosome reserves while maintaining elongation rate in slow growth conditions. Further analysis suggests the elongation rate itself is detected by sensing the ratio of dwelling and translocating ribosomes, a strategy employed to collapse the complex, high-dimensional dynamics of the molecular processes underlying cell growth to perceive the physiological state of the whole.

In the past decade, many efforts have been devoted toward characterizing and understanding the allocation of the bacterial proteome across different growth conditions (1–6). Central to the bacterial proteome allocation strategy is the approximate linear increase of the ribosome content with growth rate (7–9), when growth is varied by using different nutrients. This classic bacterial growth law is rationalized by Maaløe in terms of the need of more ribosomes to synthesize proteins to achieve faster growth rate, when the rate of translational elongation by ribosomes is saturated (9). This strategy of producing ribosomes as needed in different growth conditions forms the basis of the optimal resource allocation strategy, which posits that cells allocate their resources (the proteome in this case) in such a way to maximize its growth (5, 10). However, it is actually long known that the translational elongation rate (ER) itself varies across growth conditions (11), which poses a challenge to the rationalization by Maaløe. Moreover, it is known that in conditions where cells are hardly growing, a significant pool of ribosomes (the ribosome reserve) is kept idle, presumably for rapid transition to fast growth when favorable growth conditions return (12–14). Intriguingly, the ribosome reserve kept by cells is not limited to slow growth but maintained at a constant amount above the minimum needed across growth rates (2, 13, 15). This behavior again challenges the notion of optimal resource allocation for the current growth condition.

One approach toward understanding the bacterial proteome allocation strategy is to follow its regulatory mechanisms, to see how the linear growth law is implemented mechanistically. This involves the sensing and control of the cell’s growth rate since the proteome allocation strategy is strongly dependent on the growth rate. Toward this end, we note that Guanosine tetraphosphate (ppGpp) (16) is a key signaling molecule involved in bacterial response to environmental changes and in coordinating growth rate–dependent responses (17–20). ppGpp signaling has been extensively studied (18–24), both for mechanisms contributing to its synthesis and degradation and for its downstream effects on hundreds of genes (25–28), including the synthesis (29, 30) and activity (19, 23) of the translation machinery. Without a doubt, ppGpp signaling plays a crucial role in responding to the cellular growth rate. Yet, despite the wealth of information at the molecular level, quantitative understanding of how bacteria perceive the state of cell growth is lacking. Here we reveal the underlying signaling strategies employed by Escherichia coli to perceive and respond to growth, established through a series of experiments in which ppGpp and other key physiological variables are quantitatively measured. These strategies provide important insight on the initial question on bacterial proteome allocation strategy.

Results

A Simple, Robust Relation between ppGpp and Translational ER.

During environmental changes such as diauxic shifts, E. coli responds by producing ppGpp (31). Fig. 1 shows a typical diauxic growth curve in minimal medium containing glycerol and a small amount of glucose as the only carbon sources: Cells grow exponentially on glucose without utilizing glycerol until glucose is depleted (32), followed by a period of growth arrest (∼40 to 50 min in this case), before fully resuming growth on glycerol. We followed the kinetics of ppGpp accumulation by performing such growth transition experiments in the presence of 32P-orthophosphate. Throughout the transition, labeled nucleotides were extracted and resolved by thin-layer chromatography (Fig. 1). The ppGpp level relative to that of steady-state growth in glucose, denoted as , increased by over eightfold within the first 10 min of glucose depletion before relaxing to a new steady-state level (Fig. 1).

Fig. 1.

Relation between ppGpp and the translational ER during growth transition. (A) Growth kinetics of E. coli K-12 NCM3722 monitored by measuring the optical density at 600 nm (OD600) during the diauxic transition from glucose to glycerol. The same color scheme is used across the panels to match different measured quantities to samples taken at different time during the growth transition. (B) Resolution of ppGpp in cells sampled at different time during the growth transition by TLC. The spots at the bottom correspond to sample loading, and the migrated ppGpp are indicated alongside the chromatogram. Signal intensity from the chromatograph was used as the measure of the ppGpp level. (C) Relative ppGpp concentration, denoted by , was obtained as signal intensity from B, normalized to that in steady-state growth on glucose before the shift at . (D) The translational ER was obtained as described in , and plotted against the time at which sample was taken. (E) ppGpp levels are plotted against the reciprocal of the corresponding ERs normalized by the maximum ER (), defined as the value of the ER extrapolated to (). The line is the best fit of the data to Eq. , with and aa/s.

We also characterized changes in the translational ER [denoted by ] during the growth recovery period by assaying for the delays in LacZ induction as previous studies have established that ER determined from LacZ is representative of that of typical proteins (33, 34), and single-molecule study of translation kinetics in vivo suggested little heterogeneity in ER across codons and mRNAs in the absence of antibiotics (35, 36). Using the induction time obtained at various time () and taking into account the initiation time which showed little variation (), the instantaneous ER, , was deduced throughout the transition period (Fig. 1 and ). The data show an abrupt drop in ER immediately following glucose depletion and a gradual recovery before growth resumed. During the period of growth arrest, the time course of strikingly mirrored that of the relative ppGpp level (compare Fig. 1 ). The scatterplot of the ppGpp level with the reciprocal of ER exhibits a striking linear relation (). Defining the value of the extrapolated ER at to be , the maximum ER (), the empirical relation between the relative ppGpp level and ER can be expressed aswhere and (Fig. 1).

Since the ppGpp level and ER are also known to change with the cellular growth rate during exponential growth (33, 37), we further examined their mutual relationship during steady exponential growth. We grew E. coli with different nutrient composition at growth rates ranging from 0.13 to 0.96 h−1 () and measured the steady-state ppGpp levels relative to that in glucose, as well as the corresponding translational ERs (). As growth rate was reduced, ppGpp levels increased while ER decreased (Fig. 2 , red squares), consistent with earlier reports (33, 37). Additionally, ER has recently been shown to increase in the presence of sublethal amounts of chloramphenicol (Cm) (33). Accordingly, we observed ppGpp levels to decrease and ER to increase during steady-state growth in the presence of increasing doses of Cm (Fig. 2 , green triangles). Owing to the difficulty of detecting low ppGpp levels, we used a ΔptsG strain (NQ1261) (33) which has reduced glucose intake and thereby shows elevated ppGpp levels in the absence of Cm. This strain allowed us to quantify changes in ppGpp level and compare it to the changes in ER under Cm treatment. The scatterplot of the steady-state ppGpp level with the reciprocal of ER under carbon limitation again exhibited a linear relation (red squares, Fig. 2). Moreover, those from Cm-inhibited cells fell on the same linear relationship (green triangles). Strikingly, this is the same relationship as the one observed during the diauxic shift (compare with blue circles in Fig. 2), i.e., Eq. with the same intercept and proportionality constant.

Fig. 2.

Relation between ppGpp and translational ER during steady-state growth. (A) The average ppGpp concentrations relative to that in the glucose minimal medium, , are plotted against the corresponding growth rates for cells under steady-state carbon-limited growth (red) and translation-limited growth (green) for cells treated with sublethal doses of chloramphenicol (Cm) (). For each growth medium, ppGpp level was obtained by measuring four samples taken from exponentially growing cells at different ODs and using linear regression ( and ). Error bars represent the uncertainty in the linear fit. (B) Translation ERs are plotted against the steady-state growth rates for carbon-limited (red) and Cm-treated (green) cells. (C) Scatterplot of reciprocal ERs (or the step time for ribosome advancement) in milliseconds and the relative ppGpp concentrations () measured during steady-state growth for wild-type E. coli under carbon limitation (red) and translation inhibition (green). (D) The same measurements from steady-state growth (C) are replotted by normalizing the ER to together with the data collected under growth transition from Fig. 1 (blue symbols) for comparison. The line is the same as that in Fig. 1.

Regulatory Circuit Mediated by Translation Rate Links ppGpp Quantitatively to Growth Rate.

A steady-state relationship between ER and ppGpp level allows the cell to link the ppGpp level uniquely to the steady-state growth rate via a simple regulatory circuit (see Fig. 4): Due to negligible rate of protein turnover (41, 42), the rate of protein synthesis is given by the product of ER and the total number of active (translating) ribosomes per cell, . During exponential growth at rate , the total peptide synthesis rate is , where is the total protein mass per cell (in unit of the mass of an amino acid). The number of active ribosomes is the difference between the total number of ribosomes per cell () and the number of inactive ribosomes per cell (). Thus,

The ratio , which is proportional to the cellular ribosome concentration [since is proportional to the cell volume (14)], is set by the ppGpp level through regulation of rRNA expression (30). We take this regulatory function to be

(with an unknown constant ) since the RNA/protein mass ratio, which is proportional to , scales linearly as (Fig. 3). The inactive ribosome concentration, which is proportional to , is more difficult to quantify directly due to a multitude of ppGpp-dependent effects, e.g., the binding of ribosomes to various ribosome hibernation factors including Rmf, Hpf, and RaiA (43–45), which all increase linearly with the ppGpp level as growth rate is reduced upon limiting carbon uptake (). We describe this effect by the form

Fig. 3.

Model data comparison. (A) RNA–protein ratio (circles, left vertical axis) is proportional to the reciprocal of ppGpp level. This ratio is taken to be proportional to the ribosome content, (right vertical axis), and is used as an input to the model (Eq. 3). The proportionality constant between the RNA–protein ratio and , denoted by , is one of the three fitting parameters of the model. The best fit is shown by the line. (B) The approximate linear relation between growth rate and the reciprocal of ppGpp level is well accounted for by the model (line), as is (C) the approximate linear relation between the growth rate and RNA/protein ratio, (D) the Michaelis relation between RNA-protein ratio, and (E) the weak relation between the ER and the growth rate. RNA/protein data in A, C, and D represent mean of at least three independent biological replicates, and the error bars indicate 95% confidence interval determined from SDs. The model is described by Eqs. and . The values of the best-fit parameters were , , and .

for simplicity, with another unknown constant . In principle, the inhibitory effect of ppGpp on translation initiation (46–48) would also be included in the above. However, the magnitude of this effect is not expected to be large, as no strong delay in translational initiation was detected following transient growth arrest during the diauxic shift ().

Putting together the form of the regulatory factors in Eqs. and into Eq. leads us to a relationship between the growth rate and the ppGpp level for exponentially growing cells:where is obtained from the steady-state relation between and (Fig. 2), which is mathematically the same as inverting Eq. . The two constants and in Eq. specify the magnitudes of the two regulatory interactions. With appropriate choices of these two constants, the simple model defined by Eq. is able to quantitatively capture all the observed correlations among the growth rate , the ppGpp level (), the ribosome content (), and the ER () under nutrient limitation (Fig. 3 ), with model predictions based on best-fitted values of and shown as solid lines. In particular, the model recapitulated the well-known inverse relation between the growth rate and ppGpp level (Fig. 3) (37, 38). This illustrates the general principle that the cell can perceive its own growth rate by incorporating the sensing of ER (Eq. into a simple regulatory circuit (Fig. 4) that controls the active ribosome content by the sensor. Equally importantly, the simple correspondence between the growth rate and ppGpp level enables the cell to implement growth rate–dependent control of many cellular functions, ranging from metabolism to cell division control, by simply using ppGpp to control the expression of the relevant genes (17, 25, 28, 39, 49). In a previous study on bacterial growth control by Erickson et al. (50), an ansatz was introduced in which the translational activity was used to successfully predict gene expression dynamics during diauxic shifts. The results here establish a one-to-one relation between and (), thereby justifying the ansatz used in ref. 50.

Fig. 4.

Regulatory circuit connecting ppGpp to growth rate. The steady-state protein synthesis flux is given by the product of the ER and the number of active ribosomes, . Because is simply connected to the ppGpp level (Fig. 2 as summarized by Eq. , and the active ribosome content is given by the difference between the total ribosome content (Fig. 3 and Eq. and the content of the ribosome-sequestering elements (, and Eq. , each of which is a function of due to ppGpp-mediated regulation, it follows that is a function of . This gives rise to the correlation between the ppGpp level and growth rate (37, 38) (Fig. 3). By using ppGpp to regulate a spectrum of cellular processes (19, 22, 24), the cell thus manages to link the regulation of these processes to the growth rate, leading to the appearance of growth rate–dependent control (17, 39, 40).

Linear Bacterial Growth Law Obtained with a Special Condition on ER.

At a quantitative level, our model captured the approximate linear relation between the ribosome content and the growth rate (Fig. 3), the celebrated growth law discovered long ago (7–9). Additionally, the model captured the Michaelis-like relation between the ER and the ribosome content (Fig. 3), substantiated with extensive data collected from many conditions as reported in Dai et al. (33). Notably, a fit of the data to the Michaelis–Menten relation () recovers a maximum ER ( aa/s) that is indistinguishable from aa/s defined by taking in Eq. (). Finally, the model captures the weak dependence of the ER on growth rate (Fig. 3), with the minimal ER in slow growth condition, denoted by , whose value is close to one-half of .

The emergence of a simple linear relation between the ribosome content R and the growth rate is surprising, given the nonlinear regulatory effects exerted by ppGpp (Eq. . In fact, a variety of relationships among these quantities is possible for generic values of and (). However, if the regulatory parameters and are such that the ratio , is exactly one-half, then mathematically, the model yields an exact linear relation between and , with the slope given by , and an exact Michaelis–Menten relation between the ER and the ribosome content, with the maximal ER being (). Thus, to the extent that Eqs. and capture the forms of the regulatory functions, prescribing the appropriate regulatory parameters and to enforce being approximately one-half of is required for the emergence of the approximate linear growth relation between the ribosome content and the growth rate. [We have separately shown that adding offsets to the simplest forms of the regulatory functions used in Eqs. and does not significantly affect the quality of the fit ().]

Maintenance of ER above a minimal level is clearly of physiological importance as too low an ER would lead to problems in the processivity of protein synthesis (51, 52). Another physiological requirement is the maintenance of a sufficient ribosome reserve at slow growth, denoted by , needed for rapid growth recovery when favorable nutrient conditions return (14). Both physiological requirements are satisfied by employing hibernation factors to inactivate ribosomes. By employing both positive and negative regulation through distinct promoters (Fig. 4), the cell can readily attain the required values of and by simply prescribing the regulatory parameters and (). To keep high while also maintaining a ribosome reserve is possible in principle (). With high ER (e.g., above 90% of ; dashed line in ), the ribosome content would even be moderately reduced at fast growth rate (), a fitness benefit from the proteome allocation perspective (2, 53). However, this strategy would also require exquisite mechanism for detecting very small changes in ER (). Thus, the choice of using may reflect a compromise between the physiological demand for keeping high and the molecular constraint for detecting small changes in ER in order to sense slow growth and enforce growth rate–dependent regulation.

The Mechanism of Sensing Translational ER.

It is important to point out that the steady-state results presented here (Eq. and Fig. 3) are predicated on the existence of the empirical relation given by Eq. . We now return to discuss the causal link and the mechanistic origin of this empirical relation. Toward this end, the first question to address is whether ppGpp or ER is the primary driver of this response. One scenario is that ppGpp rises in response to some unknown starvation signal as glucose runs out, and the resulting increase in ppGpp then reduces the ER. However, we will show shortly below that a mutant in which ppGpp does not rise instantaneously still exhibits a strong immediate drop in ER as glucose runs out. An alternative scenario is that the drop of ER occurs first, and this drop is itself the signal that drives up ppGpp. The latter scenario is supported by metabolomic study which found the amino acid pools (and particularly the glutamate pool) drop sharply and immediately following glucose runout (54), thus imposing an obligatory reduction in ER. Sensing the drop in ER could therefore be an effective strategy to sense the nutritional status of the cell.

We next examine the form of the response [] in terms of the known mechanisms of ppGpp synthesis and degradation. It will be convenient to reexpress ER and in Eq. in terms of the elemental steps of the translation cycle (Fig. 5): a time where the ribosome dwells on the A site waiting for the cognate-charged tRNA and a time for peptidyl transfer and translocation to the next codon. This changes Eq. to , with the maximal ER identified with the case where . (Incidentally, the special limit of ER at slow growth being one-half of corresponds to at slow growth.)

Fig. 5.

Sensing of the translational ER by ppGpp. (A) A cycle of translation elongation constitutes the loading of a cognate charged tRNA to the A site (taking time ), followed by peptidyl transfer and the translocation of mRNA/tRNA (taking time ). The total time for one cycle, given by the reciprocal of the translational ER , is thus given by . depends on the molecular properties of the translation machinery, and depends on the concentration of uncharged tRNAs. Hence, long dwell times would lead to slow elongation speeds after a nutrient downshift owing to the increased uncharged tRNA levels, while the fastest elongation speed is obtained when , in which case . The empirical relation observed between ppGpp and the ER , Eq. , can thus be alternatively written as . According to the analysis of , the ratio of dwelling and translocation time is given by the ratio of the dwelling to translocating ribosomes, whose concentrations are and as indicated. As the transition from the dwelling state to the translocation state occurs with rate while the transition from the translocation to the dwelling state occurs with the rate , the condition of detailed balance imposes that regardless of the values of hundreds of molecular parameters underlying the translation process (). It then follows that the empirical relation [] can be obtained if, e.g., the synthesis of ppGpp is proportional to and the hydrolysis of ppGpp is proportional to . (B) According to biochemical and structural studies (24, 55, 56, 57, 58), ppGpp synthesis is activated when the A site of the ribosome is loaded by a RelA-bound uncharged tRNA. This provides a mechanistic model for the control of ppGpp synthesis rate by , the concentration of dwelling ribosomes. Although this RelA-mediated synthesis activity would provide elevated levels of ppGpp in poor nutrient conditions, it is insufficient to generate the empirical relation described in Eq. ; see , with details in . Instead, the involvement of state-dependent ppGpp hydrolysis is predicted for a full account of the empirical relation. (C) As charged tRNAs must be delivered to the ribosomal A site to complete each step of translation, the distribution of dwelling and translocating ribosomes is dependent on the metabolic fluxes directed toward all the biosynthetic precursors needed for protein synthesis, represented here by the orange cloud: triangles, amino acids; clover leaves, tRNA; colored arrows, fluxes. The translational ER is therefore a single quantity that reflects the combined status of the hundreds of diverse metabolic reactions underlying protein synthesis and cell growth (). In this sense, the mechanism of ER sensing described in B is a dimensional reduction technique employed by E. coli to convey the nutrient status by a single quantity, the level of ppGpp. The latter is further connected to the growth rate via the regulatory circuit of Fig. 4.

Next, detailed analysis based on flux balance () establishes a simple relation between two pools of actively translating ribosomes, those in the dwelling state (of concentration ) and those in the process translocation (of concentration ):

Eq. is simply a condition of detailed balance between the flux of ribosomes transitioning from the dwelling state to the translocation state and the flux transition from the translocation state back to the dwelling state, with being the total concentration of actively translating ribosomes. In terms of these ribosome pools, Eq. then becomes

i.e., the ratio of the two pools of ribosomes.

In a simple model of ppGpp composed of rapid equilibration between synthesis and degradation such that , the ppGpp level is given by the ratio of the synthesis rate () and the specific degradation rates (). A simple scenario giving rise to the empirical relation [] or [] is therefore to have the synthesis rate and the degradation rate . The effect of dwelling ribosomes on ppGpp synthesis is well supported molecularly based on the known structure of the ppGpp synthetase RelA in complex with the ribosome (55, 56) as well as earlier biochemical studies (57, 59), which find that RelA is activated (i.e., synthesizing ppGpp) only when it is complexed with an uncharged tRNA together at the A site, i.e., when the ribosome is in the dwelling state (Fig. 5 , Left). Less is known about ppGpp degradation, which is solely catalyzed by SpoT in E. coli (20, 60). The empirical relation [] excludes a model with constitutive hydrolysis of ppGpp by SpoT () and instead predicts regulation of SpoT hydrolysis activity, e.g., its stimulation by the translocating ribosomes (Fig. 5 , Right) or more complex regulation of SpoT activity depending on the ppGpp level itself as discussed in . Consistent with these ideas, we find deletion of relA to disrupt the linear relationship between ppGpp and ER during the diauxic transition (), with the remaining nontrivial ppGpp dynamics attributed to the response of SpoT, which is also the only other enzyme capable of ppGpp synthesis in E. coli (20, 24). [In this strain, a much slower accumulation of ppGpp occurred following glucose runout compared to the wild type (), yet ER dropped more and for a longer period, thus negating the aforementioned scenario that the drop of ER resulted from ppGpp accumulation.]

Intriguingly, relA deletion has no effect on either the ppGpp level or the ER in steady-state growth (), consistent with the knowledge that RelA is not essential during steady-state growth (37). The maintenance of the relationship [] by the ΔrelA strain can only be attributed to SpoT. Given Eq. , our data thus suggest that in steady state, SpoT has acquired the ability to sense both the dwelling and translocating ribosomes, in ways that it is not capable of during transient shift. How SpoT may sense different states of the ribosome is, however, not known. While there is some evidence of SpoT associated with the ribosome (61), the functional significance of this association is unclear. It is likely that this process is aided by some unknown mediators which interact with the ribosomes in ways analogous to RelA (Fig. 5) and convey that information to SpoT.

Discussion

Coordination of bacterial behaviors with growth rate is widely observed (3, 40, 49, 62). While the ppGpp signaling system is known to be central to the growth rate–dependent responses (17, 18), how growth rate is perceived and used for regulation is not known at the quantitative level. Our work establishes a missing central element in E. coli’s strategy to perceive its state of growth and respond to it: Changes in nutrient environment are immediately reflected in the rate of translational elongation, and this in turn sets the ppGpp level within a time scale of 5 to 10 min (Fig. 1). ppGpp’s well-established regulation of the ribosome content and activity, together with its robust relation with the ER (Fig. 5), then forges a unique relation between ppGpp and growth rate since the latter is given quantitatively by the product of active ribosomes and their ER.

As translational elongation is affected by many metabolic processes in the cell, monitoring the ER is an effective strategy to diagnose the state of cell growth independently of specific nutrient bottlenecks. This is the origin of the generality of the phenomenological growth laws, i.e., why the same quantitative relation between the ribosome content and the growth rate is sustained regardless of whether cell growth is limited by carbon supply (2, 7, 63), nitrogen supply (64), partial auxotrophy (3, 63), or drugs which inhibits tRNA charging (65). This mechanism also predicts generally that for perturbations not captured by a repartition of the two ribosomal states, including antibiotics (33, 66), phosphate limitation (64), and lipid stress (67), growth rate perception is distorted and the linear relation between ribosome content and growth rate is altered.

One surprising finding of our analysis is the important role played by factors that inactivate the ribosome. These factors control the amount of ribosome reserve kept by cells at slow growth (68), and the ribosome reserve is important because it affects the rapidity of growth recovery when good growth condition returns (12–14, 68). However, the necessity of employing ribosome inactivating factors can only be appreciated in light of cell physiology at slow growth: If inactivating factors are not used, then keeping a finite ribosome reserve would require the translational ER to drop to low levels in slow growth conditions, leading to problems in protein synthesis (69). Thus, these inactivating factors serve as an effective tool to maintain ER while setting aside ribosome reserve at slow growth.

Nevertheless, the deployment of these factors complicates the overall regulation of active ribosomes, making it difficult to understand how the well-known linear relation between ribosome content and growth rate arises. Our analysis shows that the linear relation emerges for a special choice of regulatory parameters such that the ER at very slow growth is one-half of the maximal ER, or alternatively, the time a ribosome spends in the dwelling state of the translation cycle is not longer than the time it spends in the translocating state, i.e., . It should be noted that this implementation of the linear ribosome–growth rate relation is very different from the existing ideas, based, e.g., on a Michaelis relation between translational factors and the ER (15) or on optimizing the steady-state growth rate (5, 10). Regarding the former, this work shows that the appearance of a Michaelis relation between the translational factors and the ER (33) is actually a consequence of regulation involving the ribosome inactivating factors. Regarding the latter, the existence of ribosome reserve which is detrimental to steady-state growth but needed for rapid transition from slow to fast growth is typically glossed over in optimization theories. However, the specific mechanism used to set the reserve, i.e., the use of ribosome inactivating factors which are needed for maintaining a finite ER in the presence of inactive ribosomes and which produces a constant offset () in the linear ribosome–growth rate relation across the entire growth rate range, calls into question the popular notion that the linear ribosome–growth rate relation is predominantly driven by the optimization of steady-state growth. Instead, it suggests that setting aside a pool of ribosomes as a reserve is something specific that the cell intends to accomplish in its proteome allocation strategy, despite having a cost to the steady-state growth rate.

Turning to the mechanisms for sensing the translational ER, our data and analysis show that it is based on the ratio of translating ribosomes in their two alternating states (Fig. 5). Sensing of the dwelling ribosomes fits well with the elaborate molecular design known for RelA (24, 57). However, our analysis suggests that this RelA-based mechanism is insufficient by itself: As RelA is not required in steady state, ppGpp synthesis activity by SpoT must somehow also be able to respond to dwelling ribosomes. The synthesis activity of SpoT was recently shown to be correlated with the levels of acetyl phosphate, a glycolytic intermediate (70). However, glycolytic flux is not necessarily a proxy for dwelling ribosomes. While SpoT has been shown to interact with ribosome-associated proteins (71), currently, it remains unclear what interactions can enable SpoT to sense dwelling ribosomes. Moreover, even in the presence of RelA, the empirical relation observed between ppGpp and ER during growth transition requires additional regulation by SpoT (). Finer details of the ppGpp signaling system may be revealed by quantifying how the ppGpp–ER relation is modified for various combinations of RelA and SpoT mutants in future studies.

At a broader level, this study provides a rare, trackable example of how cells perform dimensional reduction at the molecular level to attain crucial physiological information at the cellular level (72). The key trick E. coli uses to collapse the high-dimensional complexity of the metabolic state of a cell, e.g., involving 20 amino acid synthesis pathways and the charging of over 60 tRNAs (Fig. 5), is to take advantage of detailed balance between the two alternating states of the elongation ribosome, so that the ratio of the ribosome dwelling and translocation time, which reflects a weighted average of the tRNA charging ratios, can be deduced from the ratio of the dwelling and translocation pools of ribosomes regardless of molecular details (Fig. 5 and ). Identifying and elucidating further instances of such strategies of dimensional reduction employed by cells will be instrumental in fundamentally advancing our understanding of the connection between molecular interaction and cellular physiology for prokaryotes as well as for eukaryotes.

Methods and Materials

Growth Media Composition and Culture Conditions.

Steady-state and growth transitioning cultures were grown in MOPS (3-(N-morpholino)propanesulfonic acid)-based minimal media (73) supplemented with various carbon sources and chloramphenicol as indicated in . All cultures were grown at 37 °C with shaking at 250 rpm. For every experiment, culturing was carried out through sequential propagation of seed cultures in LB, precultures in the experimental medium, and the experimental cultures. Single colonies from fresh LB agar plates were first grown in LB broth for 6 h as the seeding culture. In the preculturing step, depending on the experiment, cells from seeding cultures were diluted into appropriate media such that the precultures would still be in exponential growth phase after overnight growth. Media used for preculturing in steady-state experiments were the same as the experimental media (). For the glucose to glycerol transitions, optical density was monitored at 600 nm (OD600) to follow the growth transition kinetics. Precultures were grown in medium supplemented with 20 mM glucose and 20 mM glycerol, to avoid glucose run-out during the preculturing step. Exponentially growing precultures were then diluted in the appropriate experimental medium (prewarmed) at an initial OD600 of ∼0.005, and various measurements were carried in the OD600 range of 0.1 to 0.4.

Strain Construction.

Wild-type E. coli K-12 NCM3722 (74, 75) and its derivatives were used in this work. HE838 (ΔrelA) was constructed using the λ-Red recombinase method (76) as follows. The km resistance gene was amplified from pKD13 using chimeric oligos relA1-P1 and relA2-P2 (). The PCR products were electroporated into NCM3722 cells expressing Lambada-Red proteins encoded by pKD46. The Km-resistant colonies were confirmed by PCR and sequencing for the replacement of the region harboring relA by the Km gene.

Translation ER Measurement.

ER was measured using LacZ as a reporter as described in Dai et al. (33) with modifications. Depending on the experiment, 10 mL cultures were grown either in different steady-state conditions or as undergoing glucose to glycerol growth transition. When cultures reached OD600 = 0.4 (for steady-state growth) or at specific time points during the growth transition, 5 mM isopropyl-β-D-thiogalactoside was added to induce the lac operon. Immediately after induction, 500 mL samples were taken at 10- or 15-s intervals to precooled ) tubes with 20 mL of 0.1 M chloramphenicol and then rapidly frozen by dry ice. Samples were stored at −80 before beta-galactosidase assay. A sensitive fluorescence substrate, 4-methylumbelliferyl-D-galactopyranoside (MUG), was used to measure LacZ activities in this work. Briefly, each sample was diluted by Z-buffer by fivefold and added to a 96-well plate to a volume of 200 mL. The plate was warmed at 37 for 10 min before adding MUG. A Tecan (SPARK) plate reader was used for MUG injection and fluorescence readings. Twenty milliliters of 2 mg/mL MUG was injected to each well, and fluorescence intensity (365-nm excitation filter, 450-nm emission filter) was measured every 4 min for 2 h. In the linear range of a fluorescence intensity vs. time plot, a linear fit was applied to obtain the slope as the relative LacZ activity for each sample. By plotting the square root of the relative LacZ activity above basal level against time (77) (Schleif plot), the lag time for the synthesis of the first LacZ molecule () was obtained for each sample (). Similar measurements using the -complement of LacZ (33) (strain NQ1468) allows us to estimate the translational initiation time (). The ER measured in this work was found to be slightly higher than that reported by Dai et al. (33), likely due to the higher sensitivity of the substrate MUG compared to ONPG (ortho-Nitrophenyl-β-galactoside) as used previously (33).

ppGpp Measurement.

ppGpp measurements were carried out as described by Cashel (78) with minor modifications. Typically, experimental cultures were grown in 3-mL volumes. Labeling was carried out when the experimental cultures grew to OD600 = 0.02 by adding 0.1 mCi 32P-orthophosphate (Perkin–Elmer) per milliliter culture. For steady-state growth, 20-μL aliquots were drawn at various OD600 values between the range 0.1 to 0.4 () and added to an equal volume of prechilled 10% formic acid. For cultures undergoing diauxic shift, 20-μL aliquots were drawn at various time points during the transition and added to an equal volume of prechilled 10% formic acid. Formic acid extracts were spun down at 13,000 rpm for 10 min, and a total of 2 μL supernatant was spotted 0.5 μL at a time near the base of a Polyethylenimine–Cellulose F thin layer chromatography (TLC) plate (Millipore). The spots were dried and nucleotides were resolved using freshly prepared 1.5 M KH2PO4 (pH 3.4). The TLC plates were dried and exposed to a phosphorimaging screen for 24 to 36 h. Chromatograms were imaged using a Typhoon FLA 9500 scanner (GE) and analyzed using Fiji software. For steady-state growth conditions, the slope of ppGpp signal intensities versus OD600 () was obtained as “ppGpp level”; this quantity is proportional to the average cellular concentration of ppGpp since cell volume per OD600 was found to depend weakly on growth conditions (79, 80). To control for changes in radioactivity between batches, the ppGpp level from a sample of NCM3722 growing steadily in MOPS glucose was always measured as a reference. All ppGpp levels in that batch were normalized to that of the glucose-grown reference in the same batch, and the result is reported as relative ppGpp levels throughout the manuscript.

Total RNA and Protein Measurement.

Total RNA was measured using the method of Benthin et al. (81), and protein was measured using the Biuret method (82), with minor modifications as described by You et al. (63).