Interdependence of Cell Growth and Gene Expression: Origins and Consequences

See allHide authors and affiliations

Science  19 Nov 2010:
Vol. 330, Issue 6007, pp. 1099-1102
DOI: 10.1126/science.1192588


In bacteria, the rate of cell proliferation and the level of gene expression are intimately intertwined. Elucidating these relations is important both for understanding the physiological functions of endogenous genetic circuits and for designing robust synthetic systems. We describe a phenomenological study that reveals intrinsic constraints governing the allocation of resources toward protein synthesis and other aspects of cell growth. A theory incorporating these constraints can accurately predict how cell proliferation and gene expression affect one another, quantitatively accounting for the effect of translation-inhibiting antibiotics on gene expression and the effect of gratuitous protein expression on cell growth. The use of such empirical relations, analogous to phenomenological laws, may facilitate our understanding and manipulation of complex biological systems before underlying regulatory circuits are elucidated.

Systems biology is as an integrative approach to connect molecular-level mechanisms to cell-level behavior (1). Many studies have characterized the impact of molecular circuits and networks on cellular physiology (1, 2), but less is known about the impact of cellular physiology on the functions of molecular networks (35). Endogenous and synthetic genetic circuits can be strongly affected by the physiological states of the organism, resulting in unpredictable outcomes (4, 68). Consequently, both the understanding and implementation of molecular control are predicated on distinguishing global physiological constraints from specific regulatory interactions.

For bacterial cells under steady-state exponential growth, the rate of cell proliferation (the “growth rate”) is an important characteristic of the physiological state. It is well known that the macromolecular composition (e.g., the mass fractions of protein, RNA, and DNA) of bacterial cells under exponential growth depends on the growth medium predominantly through the growth rate allowed by the nutritional content of the medium (9, 10). Such growth rate dependencies inevitably affect the expression of individual genes (4, 11) because protein synthesis is directly dependent on the cell’s ribosome content. The latter is reflected by the RNA/protein ratio. In Escherichia coli, most of the RNA (~85%) is rRNA folded in ribosomes (10, 11). A predictive understanding of the impact of growth physiology on gene expression therefore first requires an understanding of the cell’s allocation of cellular resources to ribosome synthesis (manifested by the RNA/protein ratio) at different growth rates.

For exponentially growing E. coli cells (10, 12), the RNA/protein ratio r is linearly correlated with the specific growth rate λ [ = (ln 2)/doubling time] (Fig. 1A). The correlation is described mathematically as r=r0+λκt(1)where r0 is the vertical intercept and κt is the inverse of the slope (table S1). This linear correlation holds for various E. coli strains growing in medium that supports fast to moderately slow growth [e.g., 20 min to ~2 hours per doubling (11)], and it appears to be quite universal; similar linear correlations have been observed in many other microbes, including slow-growing unicellular eukaryotes (fig. S1). As suggested long ago from mass-balance considerations (11) and elaborated in (13), this linear correlation is expected if the ribosomes are growth-limiting and are engaged in translation at a constant rate, with the phenomenological parameter κt predicted to be proportional to the rate of protein synthesis. Consistent with the prediction, data on RNA/protein ratios from slow-translation mutants of E. coli K-12 (triangles in Fig. 1B) also exhibited linear correlations with the growth rate λ, but with steeper slopes than the parent strain (circles), which have smaller κt. Moreover, the corresponding κt values correlated linearly with the directly measured speed of translational elongation (14) (Fig. 1B, inset). Consequently, we call κt the “translational capacity” of the organism.

Fig. 1

Correlation of the RNA/protein ratio r with growth rate λ for various strains of E. coli. (A) Comparison among E. coli strains grown in minimal medium: Strain B/r [(10), squares], 15τ-bar [(12), diamonds], and EQ2 (this work, solid circles). The growth rate is modulated by changing the quality of nutrients as indicated in the key at lower left. The fraction of total protein devoted to ribosome-affiliated proteins (ϕR) is given by the RNA/protein ratio as ϕR = ρ · r (table S1). (B) The RNA/protein ratio for a family of translational mutants SmR (triangles) and SmP (inverted triangles) and their parent strain Xac (circles) (27), grown with various nutrients (see key at lower left) (table S2). Translational inhibition of the parent Xac strain via exposure to sublethal doses of chloramphenicol (circled numbers; see legend table) gave RNA/protein ratios similar to those of the mutant strains grown in medium with the same nutrient but without chloramphenicol (light blue symbols). Dashed line is a fit to Eq. 2. Inset: Linear correlation of κt values obtained for the Xac, SmR, and SmP strains (table S2) with the measured translation rate of the respective strains (14) (r2 = 0.99).

Translation can be inhibited in a graded manner by exposing cells to sublethal doses of a translation-inhibiting antibiotic. The RNA/protein ratios obtained for wild-type cells grown in medium with a fixed nutrient source and various amounts of chloramphenicol (Fig. 1B, light blue circles) were consistent with data obtained for the isogenic translational mutants grown in medium with the same nutrient but no antibiotic (light blue triangles). Surprisingly, these data revealed another linear correlation between r and λ (Fig. 1B, dashed line), given byr=rmaxλκn(2)where rmax is the vertical intercept and κn is the inverse slope. Such a linear correlation was obtained for cells grown with each of the six nutrient sources studied (Fig. 2A and table S3). The correlation described by Eq. 2 has been observed in cells subjected to numerous other means of imposing translational limitation (fig. S2).

Fig. 2

Effect of translational inhibition. (A) RNA/protein ratio for strain EQ2 grown in different media, each with various levels of chloramphenicol (see key at lower right) (table S3). Solid lines were obtained from fitting data of the same color to Eq. 2. The black line describes the data in the absence of chloramphenicol (as in Fig. 1A). (B) Translational inhibition results in an increased synthesis of R-class proteins (cyan), effectively decreasing the fraction allocable to the P-class (magenta). (C) Mass fraction of constitutively expressed β-galactosidase (strain EQ3) plotted as a function of growth rate. The lines were fit according to Eq. 4. The growth rate dependence of constitutive gene expression due to nutrient limitation found in (4) is also well described by the theory (fig. S5C).

From Fig. 2A and the best-fit values of the parameters rmax and κn (table S4), we observed that the parameter κn exhibited a strong, positive correlation with the growth rate of cells in drug-free medium (fig. S3A). Thus, κn reflects the nutrient quality and is referred to as the “nutritional capacity” of the organism in a medium [see eq. S18 in (13) for a molecular interpretation of κn]. In contrast, the vertical intercept rmax depended only weakly on the composition of the growth medium (fig. S3B). Qualitatively, the increase of the RNA/protein ratio r with increasing degree of translational inhibition can be seen as a compensation for the reduced translational capacity, implemented possibly through the relief of repression of rRNA synthesis by the alarmone ppGpp (15), in response to the buildup of intracellular amino acid pools resulting from slow translation. Because rmax is the (extrapolated) maximal RNA/protein ratio as translation capacity is reduced toward zero, its weak dependence on the quality of the nutrients suggests a common limit in the allocation of cellular resources toward ribosome synthesis.

The simplest model connecting ribosome abundance to gene expression assumes that the total protein content of the cell (called the proteome) is composed of two classes: ribosome-affiliated “class R” proteins (with mass fraction ϕR), and “others” (with mass fraction 1 – ϕR) (5, 16). But the maximum allocation to the R-class proteins as derived from the value of rmax, ϕRmax=ρrmax0.55, is well below 1 [see (13) for the conversion factor ρ]. This suggests that the “other” proteins can be further subdivided minimally into two classes (Fig. 2B): “class Q” of mass fraction ϕQ, which is not affected translational inhibition, and the remainder, “class P” of mass fraction ϕP, with ϕP → 0 as ϕRϕRmax (17). Because ϕP + ϕQ + ϕR = 1, we must have ϕRmax=1ϕQ, withϕP=ϕRmaxϕR=ρ(rmaxr)(3)representing an important constraint between ϕP and ϕR. Together with Eq. 2, the model predictsϕP=ρλ/κn(4)which describes a linear relation between the abundance of the P-class proteins and the growth rate λ for a fixed nutritional capacity κn. The growth rate independence of protein abundance may be maintained through negative autoregulation (4) (fig. S4). Unregulated (or “constitutively expressed”) proteins belong instead to the P-class and can be used to test the prediction of Eq. 4: Expression of β-galactosidase driven by a constitutive promoter (ϕZ, mass of β-galactosidase per total protein mass) in cells grown under different degrees of chloramphenicol inhibition indeed correlated linearly with λ for each nutrient source studied (Fig. 2C), and the slopes of these correlations (colored lines) agree quantitatively with the nutritional capacity κn (fig. S5, A and B) as predicted by Eq. 4.

Although the correlations (Eqs. 2 and 4) were revealed by growth with antibiotics, their forms do not pertain specifically to translational inhibition. Equation 4 may be interpreted as a manifestation of P-class proteins providing the nutrients needed for growth [eqs. S15 to S18 in (13)], just as Eq. 1 is a reflection of R-class proteins providing the protein synthesis needed for growth (Fig. 3A). For different combinations of the nutritional and translational capacities (κn, κt), efficient resource allocation requires that the abundance of P- and R-class proteins be adjusted so that the rate of nutrient influx provided by P (via import or biosynthesis) matches the rate of protein synthesis achievable by R (Fig. 3B), while simultaneously satisfying the constraint of Eq. 3 (Fig. 3C). We can derive the resulting allocation mathematically by postulating that λ, ϕR (or r), and ϕP are analytical functions of the variables κt and κn that respectively capture all molecular details of translation and nutrition (analogous to state variables in thermodynamics). The mathematics is identical to the description of an electric circuit with two resistors (fig. S6), with Eqs. 1 and 4 being analogous to Ohm’s law. Solving these equations simultaneously leads to the Michaelis-Menten relation known empirically for the dependence of cell growth on nutrient level (18) λ(κt,κn)=λc(κt)κnκt+κn(5)The value of the maximal growth rate λct) = κt⋅(rmaxr0) ≈ 2.85 hour–1 (based on the average rmax) corresponds well to the doubling time of ~20 min for typical E. coli strains in rich media. Moreover, Eq. 5 quantitatively accounts for the correlation observed between growth rate λ and nutritional capacity κn (fig. S3A).

Fig. 3

A phenomenological theory of bacterial growth. (A) The growth theory comprises three key ingredients: (i) a three-component partition of the proteome, consisting of a fixed core sector (Q) and two adjustable sectors (R and P) whose fractions (ϕR and ϕP) must add up to a constant (Embedded Image); (ii) a ribosomal fraction ϕR containing all the ribosomal proteins and their affiliates and exerting a positive effect on growth (with growth rate λ ∝ ϕR – ϕ0, where ϕ0 = ρ · r0 corresponds to the vertical intercept in Fig. 1A); (iii) a remaining fraction ϕP exerting a similarly positive effect on growth (with growth rate λ ∝ ϕP) by providing an influx of nutrients (13). (B) During steady-state exponential growth, efficient resource allocation requires that the nutrient influx (κn · ϕP) be flux-matched to the amino acid outflux κt · (ϕR – ϕ0). This can be coordinated by ppGpp, which up-regulates ribosome synthesis and hence amino acid outflux in response to increase in the amino acid pools, and has the opposite effect (down-regulating ribosome synthesis and hence amino acid outflux) in response to decrease in the amino acid pools (15). Changes in ϕR also indirectly regulate nutrient influx through the constraint of Eq. 3, in addition to direct regulatory mechanisms (dashed line). (C) Balancing the demands of protein synthesis and nutrient influx leads to the constraint (ϕR – ϕ0)/ϕP = κnt, sketched as the dashed black line. The other constraint (Eq. 3) is shown as the solid black line. The unique combination (ϕP, ϕR) satisfying both constraints is shown as the solid black circle. Upon increasing κn and/or decreasing κt, ϕR needs to be increased and ϕP decreased to maintain the balance (green line and circle), whereas with decreasing κn and/or increasing κt, ϕP needs to be increased and ϕR decreased (red line and symbols).

This theory can be inverted to predict the effect of protein expression on cell growth. Unnecessary protein expression leads to diminished growth (19). Understanding the origin of this growth inhibition is of value in efforts to increase the yield of heterologous protein in bacteria (20) and to understand the fitness benefit of gene regulation (21, 22). Aside from protein-specific toxicity, several general causes of growth inhibition have been suggested, including diversion of metabolites (23), competition among sigma factors for RNA polymerases (24), and competition among mRNA for ribosomes (19, 25).

We modeled the expression of unnecessary protein (of mass fraction ϕU) as an additional (neutral) component of the proteome that effectively causes a reduction of rmax to rmax – ϕU/ρ (Fig. 4A). Equation 5 then predicts a linear reduction of the growth rate, λ(ϕU)=λ(ϕU=0)[1(ϕU/ϕc)](6)extrapolating toward zero growth at ϕc = ρ · (rmaxr0) ≈ 0.48. The prediction quantitatively described the observed growth defect caused by inducible expression of β-galactosidase (Fig. 4B), as well as previous results obtained for various proteins and expression vectors (Fig. 4C) (19, 26), without any adjustable parameters. These results suggest that growth reduction is a simple consequence of ribosome allocation subject to the constraints of Eqs. 1, 3, and 4.

Fig. 4

Effect of unnecessary protein expression on growth. (A) Expression of an unnecessary protein (U) effectively decreases the fraction allocable to the R and P sectors, leading to a decrease in growth rate. (B) The growth defect caused by unnecessary protein expression in strain EQ30 (13). Data were collected from cells grown in three different media (same symbols as Fig. 2A) with various degrees of LacZ induction. The lines are predictions by Eq. 6, using the growth rate λ(ϕU = 0) and intercept ϕc = 0.48 computed from the average rmax (table S4). (C) Comparison of unnecessary protein expression from various expression systems: β-lactamase driven by the σ70-dependent bla promoter (open green circles) on a medium-copy plasmid (26), β-galactosidase driven by a T7 promoter (open red circles) or truncated EF-Tu driven by the σ70-dependent tac promoter (open black circles), both on high-copy plasmids (19). β-Galactosidase driven by the σ54-dependent Pu promoter (solid blue circles) on a medium-copy plasmid [as in (B)] is shown for comparison.

Robust empirical correlations of the RNA/protein ratio with the growth rate (Figs. 1A and 2A and figs. S1 and S2) revealed underlying constraints of cellular resource allocation and led to the formulation of a simple growth theory that provided quantitative predictions and unifying descriptions of many important but seemingly unrelated aspects of bacterial physiology. Like Ohm’s law, which greatly expedited the design of electrical circuits well before electricity was understood microscopically, the empirical correlations described here may be viewed as microbial “growth laws,” the use of which may facilitate our understanding of the operation and design of complex biological systems well before all the underlying regulatory circuits are elucidated at the molecular level.

Supporting Online Material

Materials and Methods

SOM Text

Figs. S1 to S6

Tables S1 to S7


References and Notes

  1. See supporting material on Science Online.
  2. A particular protein species may belong to multiple classes; sometimes this is a result of expression from multiple promoters that are differently regulated.
  3. We thank H. Bremer, L. Csonka, P. Dennis, M. Ehrenberg, P. Geiduschek, M. Schaechter, A. Tadmor, and members of the Hwa lab for suggestions and discussions; D. Hughes for the Sm mutant strains; and P.-h. Lee and B. Willumsen for the use of unpublished data. Supported by NIH grant RO1GM77298, NSF grant MCB0746581, and the NSF-supported Center for Theoretical Biological Physics (grant PHY0822283). M.S. was supported by a Natural Sciences and Engineering Research Council of Canada fellowship.
View Abstract

Stay Connected to Science

Navigate This Article