- Open Access
Dilution and the theoretical description of growth-rate dependent gene expression
Journal of Biological Engineeringvolume 7, Article number: 22 (2013)
Expression of a gene is not only tuned by direct regulation, but also affected by the global physiological state of the (host) cell. This global dependence complicates the quantitative understanding of gene regulation and the design of synthetic gene circuits. In bacteria these global effects can often be described as a dependence on the growth rate. Here we discuss how growth-rate dependence can be incorporated in mathematical models of gene expression by comparing data for unregulated genes with the predictions of different theoretical descriptions of growth-rate dependence. We argue that a realistic description of growth effects requires a growth-rate dependent protein synthesis rate in addition to dilution by growth.
Genetic circuits are unavoidably coupled to the global state of their host cell, which provides the machinery for gene expression . Rather than being a rigid “chassis” for these gene circuits, the host cell itself is dynamic and adapts to external conditions, complicating the predictive design of synthetic gene circuits [2, 3]. In bacterial cells, the most important characteristic of the global state of the cell is the growth rate, and pmeters of the host cell (such as the ribosome and RNA polymerase content) can often be described as growth-rate dependent . Aiming at a quantitative understanding of genetic circuits, a number of recent studies have therefore addressed growth-rate dependent effects [1, 5–9]. In experiments, such global effects can be studied by comparing the expression of regulated genes with an unregulated reference gene that only reflects global changes in gene expression. To include the global effects in mathematical models of gene circuits, a mathematical description of the growth-rate dependence of an unregulated reference gene is required. In this letter we discuss and compare several version of such mathematical description and compare their predictions with a compilation of experimental data .
In many modeling approaches, the dynamics of the concentration (p) of a protein is described by a balance of synthesis and (effective) degradation with rates α and β, respectively,
The synthesis rate α includes the rates of transcription and translation as well as the degradation rate of the corresponding mRNA and the concentration of the corresponding gene. For regulated genes, α is furthermore a function of the concentrations of the corresponding transcription factors or other regulators. Here, to septe global growth effects and specific regulation, we consider unregulated (constitutive) gene expression, without such dependencies. The degradation rate in this equation is an effective degradation rate, β=β0+λ, given by the sum of the rates of actual degradation by proteolysis (β0) and of dilution due to cell growth with the dilution rate or growth rate λ (which is related to the doubling rate μ via λ=μ ln2). Most proteins in E. coli are not actively degraded during balanced growth , so typically dilution by cell growth is dominant, β≈λ.
When effects of growth are considered in models of gene circuits, they are often identified with this dilution effect and the synthesis rate α is assumed to be constant, independent of the growth rate. As a result, the steady state protein concentration is proportional to the inverse of the growth rate, p(λ)=α/λ. Thus, if these assumptions are correct, the product of growth rate and the concentration of a constitutively expressed protein, p(λ)×λ, should be constant (black lines in Figure 1). In Figure 1 we plot this product (which corresponds to an effective synthesis rate) using a compilation of data from ref.  for the growth rate dependence of protein concentrations for unregulated genes in E. coli. Despite the scatter in the data at the highest growth rates, the product can be seen to be approximately constant for moderate to rapid growth, but it clearly decreases when the growth rate is reduced below 1 doubling/hour.
This additional growth dependence could be explained by a growth dependent synthesis rate or be due to a different growth dependence of the degradation rate. The next simplest possibility is to include active degradation with a constant degradation rate β0. Then the steady state of a constitutive gene reads p(λ)=α/(β0+λ) and the product p×λ is given by αλ/(β0+λ). A least square fit of this hyperbolic function to the data leads to the blue curve in Figure 1 with a degradation rate β0 = 0.29 hr−1, corresponding to a protein half-life of ~2.4 hours. This value is much smaller than the half-lifes measured for the bulk of proteins in exponentially growing E. coli, which has been determined to be ≥10 hours [10, 11]. Likewise, if we specifically fit the data with a beta-galactosidase reporter (Ptet:lacZ data from ref. ), we obtain a half-life of ~1.4 hours, also much smaller than the value of ≥10 hours reported for this protein . We therefore conclude that this approach to modeling growth-rate dependence, while providing a good fit to the data, overestimates protein degradation.
Thus, the growth-rate dependence of p×λ has to be attributed either to a growth-rate dependent protein synthesis rate α or to a growth-rate dependent degradation rate β0. While the latter case cannot be excluded, the extent of growth-rate dependent effects in protein degradation is not very clear. Increased degradation rates have been reported for non-growing or slowly growing cells, but either with small degradation rates (≈0.05 hr-1)  or with only a small fraction of protein (2–5%) that is rapidly degraded a. On the other hand, a decrease of the protein synthesis rate at slow growth has recently been reported , likely reflecting the reduced availability of RNA polymerases under these conditions b. Assuming again stable proteins, i.e. β=λ, an excellent fit to the data in Figure 1 is obtained with the exponential dependence
The resulting curve with λ* = 0.33 hr-1 is shown in red in Figure 1. We note that protein degradation cannot be neglected for arbitrarily slow growth, but including a degradation rate β0 by fitting with pλ=α(λ)×λ/(λ+β0) has only a minor effect on the fit (dashed red line in Figure 1, with λ* = 0.24 hr-1,β0 = 0.06 hr-1).
In summary, we propose to describe an unregulated reference gene by dilution and a growth rate dependent protein synthesis rate (Equation 2). The two simplest and often-used descriptions of growth effects are not sufficient: A description with constant synthesis and degradation rates and dilution is only consistent with the data if an unrealistically high degradation rate is assumed. The data is clearly inconsistent with models where growth only affects the dilution of stable proteins.
aFor comparison, attributing all growth rate dependence of p(λ)×λ to the degradation rate, β0 would be given by λ* from Equation 2 at slow growth.
bThe reported synthesis rates are for a constant gene concentration. This condition is expected to be approximately fulfilled for chromosomal genes at slow growth, where the gene copy number is given by the copy number on the chromosome and the cell volume varies only weakly .
Klumpp S, Zhang Z, Hwa T: Growth rate-dependent global effects on gene expression in bacteria. Cell 2009, 139: 1366-1375. 10.1016/j.cell.2009.12.001
de Lorenzo V: Beware of metaphors: chasses and orthogonality in synthetic biology. Bioengineered bugs 2011, 2: 3-7. 10.4161/bbug.2.1.13388
Kwok R: Five hard truths for synthetic biology. Nature 2010, 463: 288-290. 10.1038/463288a
Bremer H, Dennis PP: Modulation of chemical composition and other pmeters of the cell by growth rate. In Escherichia coli and Salmonella. Washington D.C: Neidhardt FC; 1996:1553-1569.
Tan C, Marguet P, You L: Emergent bistability by a growth-modulating positive feedback circuit. Nat Chem Biol 2009, 5: 842-848. 10.1038/nchembio.218
Scott M, Gunderson CW, Mateescu EM, Zhang Z, Hwa T: Interdependence of cell growth and gene expression: origins and consequences. Science 2010, 330: 1099-1102. 10.1126/science.1192588
Klumpp S: Growth-rate dependence reveals design principles of plasmid copy number control. PLoS One 2011, 6: e20403. 10.1371/journal.pone.0020403
Berthoumieux S, de Jong H, Baptist G, Pinel C, Ranquet C, Ropers D, Geiselmann J: Shared control of gene expression in bacteria by transcription factors and global physiology of the cell. Mol Syst Biol 2013, 9: 634.
Gerosa L, Kochanowski K, Heinemann M, Sauer U: Dissecting specific and global transcriptional regulation of bacterial gene expression. Mol Syst Biol 2013, 9: 658.
Nath K, Koch AL: Protein degradation in Escherichia coli. I. Measurement of rapidly and slowly decaying components. J Biol Chem 1970, 245: 2889-2900.
Larrabee KL, Phillips JO, Williams GJ, Larrabee AR: The relative rates of protein synthesis and degradation in a growing culture of Escherichia coli. J Biol Chem 1980, 255: 4125-4130.
Varshavsky A: The N-end rule: functions, mysteries, uses. Proc Natl Acad Sci USA 1996, 93: 12142-12149. 10.1073/pnas.93.22.12142
Mandelstam J: The intracellular turnover of protein and nucleic acids and its role in biochemical differentiation. Bacteriological reviews 1960, 24: 289-308.
Klumpp S, Hwa T: Growth-rate-dependent partitioning of RNA polymerases in bacteria. Proc Natl Acad Sci USA 2008, 105: 20245-20250. 10.1073/pnas.0804953105
The authors declare that they have no competing interest.
MH and SK conceived and performed the analysis. SK wrote the paper. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.