Journal Pcbi 1010087
Journal Pcbi 1010087
RESEARCH ARTICLE
circuit, we examined gene expression in single cells by using a double reporter system to moni-
tor both the regulator (MS2CP) and the regulated gene, and we also developed a mathematical
model to provide a predictive quantitative foundation of the system.
Results
Regulation of translation with an RNA-binding protein in single cells
We engineered a synthetic genetic system in E. coli in which the RNA-binding protein MS2CP
acts as a protein translation factor (Fig 1A). MS2CP was expressed from a synthetic PL-based
promoter repressed by LacI (named as PLlac) [21] in a medium copy number plasmid (about
80 copies/cell). This allowed controlling the expression of the regulator (at the transcriptional
level) with isopropyl β-D-1-thiogalactopyranoside (IPTG). In addition, we fused the enhanced
blue fluorescent protein 2 (eBFP2) [22] to the N terminus of MS2CP (leading to
eBFP2-MS2CP) in order to monitor its expression. As a regulated element, here we used the
superfolder green fluorescent protein (sfGFP) [23], which was expressed from a constitutive
promoter in a low copy number plasmid (about 15 copies/cell). The wild-type RNA motif rec-
ognized by MS2CP (with a dissociation constant of about 3 nM) [24] was placed in frame just
after the start codon of sfGFP. In this way, MS2CP can block the progression of the ribosome
on the regulated gene in the initial phase [15]. This mode of action differs from the natural
one, in which MS2CP prevents translation initiation rather than elongation [16]. The resulting
circuit behaves like an inverter considering IPTG as input and sfGFP as output, MS2CP being
an internal regulator that operates at the level of translation.
We performed single-cell measurements of blue and green fluorescence by flow cytometry
for a concentration gradient of IPTG (9 conditions) in order to quantitatively study the sto-
chastic regulatory dynamics of this engineered system (Fig 1B and 1C). We found a substantial
down-regulation of sfGFP (about 50-fold in expression) as a consequence of the action of
MS2CP on the cognate mRNA. From these data, we calculated the mean and the noise of
expression for both eBFP2-MS2CP and sfGFP (the noise as the square of the coefficient of var-
iation) [3], which were represented as a function of IPTG (Fig 1D and 1E). The mean gives the
average position of the population, and the noise is a measure of the cell-to-cell variability.
These measurements were repeated for different populations, finding consistency in the results
(S1 Fig). We then constructed a mathematical model relying on a series of algebraic equations
from basics on the biochemistry of gene expression and molecular noise propagation [7]. To
derive these mathematical expressions for the mean and the noise, we constructed a system of
stochastic differential equations for mRNA and protein expression following the Langevin for-
malism. The rates of concentration changes were subject to stochastic fluctuations of intrinsic
and extrinsic nature. This system was analytically solved in steady state with the mean-field
approximation for the fluctuations. With a suitable parameterization, our model was able to
recapitulate with reasonable agreement the values of mean expression and noise for both
eBFP2-MS2CP and sfGFP, highlighting the functional form of the different dose-response
curves. In particular, the mean expressions follow Hill-Langmuir equations and the noises
non-monotonous curves presenting a maximum at an intermediate IPTG concentration.
Indeed, the peak-like noise curve is a consequence of a sigmoidal dynamics at the population
level. We also observed that the noise levels in sfGFP are lower than in eBFP2 for all IPTG con-
centrations. We also performed numerical simulations of the stochastic differential equations
(S2 Fig), finding good agreement with the analytical results, as well as sensitivity analyses to
reveal the effect of perturbations in the adjusted parameters (S3 Fig), highlighting how the
curves of mean expression and noise shift in one direction and even change in form. In addi-
tion, we represented the noise versus the mean to show the stochastic expression scaling laws
Fig 1. Regulation with a protein translation factor. a) Schematics of the gene regulatory system implemented in a bacterial cell. IPTG is the external molecule
that controls the expression of the protein translation factor (eBFP2-MS2CP). sfGFP is the final output of the system. b) Histograms of single-cell fluorescence
for eBFP2 (fused to the regulatory protein) for different induction conditions with IPTG. On top, sequence details of the cis-regulatory region (DNA level) for
transcriptional regulation (PLlac promoter). c) Histograms of single-cell fluorescence for sfGFP (the regulated protein) for different induction conditions with
IPTG. On top, sequence details of the cis-regulatory region (RNA level) for post-transcriptional regulation (MS2CP RNA motif). d) Mean and noise of
expression for eBFP2 as a function of IPTG. e) Mean and noise of expression for sfGFP as a function of IPTG. f) Noise for eBFP2 as a function of mean
expression. g) Noise for sfGFP as a function of mean expression. h) Transfer function of the post-transcriptional regulation in terms of mean expression. i)
Transfer function of the post-transcriptional regulation in terms of noise. In plots d-i), points correspond to calculations from the experimental data, while
solid lines to predictions with the mathematical model.
https://doi.org/10.1371/journal.pcbi.1010087.g001
of the system (Fig 1F and 1G). The model was also explicative about the nonlinear transfer
functions in terms of mean expression regulation (Fig 1H) and noise propagation (i.e., how
the noise of eBFP2-MS2CP impacts on the noise of sfGFP; Fig 1I). Together, these results indi-
cated that the protein translation factor is a suitable element to control expression and that the
cell-to-cell variability emerged at this level can be predicted with certain accuracy.
Fig 2. Detailed analysis of stochastic gene expression. a) Model-based calculation of the total noise in eBFP2 expression as a function of IPTG. On top,
schematics of the amplification effect by the transcriptional regulation. b-d) Decomposition of the total eBFP2 noise into extrinsic, intrinsic, and regulation
noise components. e) Model-based calculation of the total noise in sfGFP expression as a function of IPTG. On top, schematics of the buffering effect by the
post-transcriptional regulation. f-h) Decomposition of the total sfGFP noise into extrinsic, intrinsic, and regulation noise components. Insets in c,g) show the
scaling of the intrinsic noise with the mean expression. i) Predicted eBFP2 fluorescence distributions for different induction conditions with IPTG (Gamma
distributions). j) Gamma shape and scale parameters for eBFP2. k) Predicted sfGFP fluorescence distributions for different induction conditions with IPTG
(Gamma distributions). l) Gamma shape and scale parameters for sfGFP. In plots j,l), points correspond to calculations from the experimental values of mean
and noise, while solid lines to predictions with the mathematical model.
https://doi.org/10.1371/journal.pcbi.1010087.g002
IPTG concentrations. Overall, this highlighted the generality of the Gamma distribution to
describe genetic systems regulated at both transcriptional and translational levels.
global response mechanism in bacteria against this type of antibiotics) [31,32]. That is, the cell
is able to sense that a substantial amount of ribosomes is being inhibited upon binding to TC
and produces more. In particular, TC binds to the 30S subunit and interferes with the transfer
RNAs (tRNAs). In turn, the cell growth rate is compromised due to the action of TC. Impor-
tantly, this parameter has been shown to modulate the mean and noise of gene expression
[33,34], so we decided to exploit it as a predictor variable. Over a two-dimensional concentra-
tion gradient of IPTG and TC (81 conditions), we first generated growth curves (S5 Fig). Basi-
cally, only TC showed a significant impact on growth rate (Fig 3B), with a maximal reduction
of almost 3-fold, which was well explained by a Michaelis-Menten function (Fig 3C).
In parallel, we performed single-cell measurements of blue and green fluorescence for each
condition (Fig 3D). We observed that the mean expression levels of both eBFP2-MS2CP and
sfGFP remained almost constant at low TC concentrations, but they increased significantly
from 500 ng/mL TC, irrespective of the induction with IPTG (Fig 3E and 3F). Because protein
expression comes from the ratio between the protein synthesis rate (accounting for both tran-
scription and translation) and the growth rate (in the case of stable proteins, as it is the case
here), this indicated that the protein synthesis rate of both eBFP2-MS2CP and sfGFP scales
with the growth rate (Fig 3G and 3H). It was interesting to note here the logical NOR behavior
of the sfGFP synthesis rate and the difference between protein expression and synthesis rate.
In addition, we calculated the noise levels for each condition (Fig 3I and 3J). We observed that
the regulation noise decreases with TC for both eBFP2-MS2CP and sfGFP, as well as that TC
leads to a substantial increase in the sfGFP noise when this gene is fully repressed by MS2CP.
Fig 3. Growth-dependent regulation with a protein translation factor. a) Extended schematics of the gene regulatory
system in which TC further modulates it through its negative impact on growth rate and positive impact on translation
rate (global effects). b) Heatmap of the mean growth rate as a function of IPTG and TC. c) Dose-response curve between
growth rate and TC. Points correspond to experimental data, while solid line comes from the mathematical model. d)
Projected two-dimensional histograms of single-cell fluorescence for eBFP2 and sfGFP for different induction
conditions with IPTG and TC. e) Heatmap of the mean eBFP2 fluorescence as a function of IPTG and TC. f) Heatmap of
the mean sfGFP fluorescence as a function of IPTG and TC. g) Heatmap of the mean eBFP2 synthesis rate as a function
of IPTG and TC. h) Heatmap of the mean sfGFP synthesis rate as a function of IPTG and TC. i) Heatmap of the eBFP2
noise as a function of IPTG and TC. j) Heatmap of the sfGFP noise as a function of IPTG and TC.
https://doi.org/10.1371/journal.pcbi.1010087.g003
the number of MS2CP molecules per cell changes with IPTG, so there is a coupling between the
translational regulation and the effect of growth rate on expression. While at low IPTG concentra-
tions the relative sfGFP synthesis rate follows the aforementioned trend for eBFP2-MS2CP, at
high IPTG concentrations there is a maximum at a growth rate of about 0.65 h-1 (it was particu-
larly pronounced at the intermediate value of 100 μM). Fig 4F illustrates how the transfer function
in terms of mean expression varies with the TC concentration (i.e., by increasing the maximal
expression level and shifting the inflexion point towards the right).
Finally, we applied the model to project the noise in protein expression. In this case, we
needed to introduce a phenomenological dependence with the growth rate on three noise-
related parameters to explain the data. In particular, we set that the noise in LacI expression
scales with the square of the growth rate (i.e., LacI expression varies from cell to cell in greater
extent when cells grow faster) and that the extrinsic noise of both eBFP2-MS2CP and sfGFP
scales with the inverse of the growth rate by following the translation rate (i.e., the extrinsic
noise is higher at lower growth rates, which seems in tune with recent experiments characteriz-
ing genome-wide noisy expression levels [36]). While for eBFP2-MS2CP the noise decreases
or remains constant with the growth rate (Fig 4G), for sfGFP the noise presents a more com-
plex trend (almost constant at low IPTG concentrations and with a maximum at high IPTG
concentrations; Fig 4H). In turn, Fig 4I illustrates how the transfer function in terms of noise
varies with the TC concentration, showing how the belly shape is reduced with TC, which indi-
cates that noise propagation through the translation factor is less significant (i.e., it is masked)
when the growth rate is low. That is, the regulation noise term, which quantifies how much
fluctuation sfGFP perceives from MS2CP, becomes smaller than the other noise terms (intrin-
sic and extrinsic) with TC. Arguably, at high growth rates, when global perturbations are
small, a translation factor is superior to a transcription factor because it is able to regulate gene
expression without transmitting much noise. However, this is not the case at low growth rates,
when global perturbations become substantial, due to a poor signal-to-noise ratio. Together,
these results highlight the complex impact of growth rate in the system and serve to appreciate
how global and local regulatory mechanisms interplay in the cell.
Discussion
Our development follows previous work on exploiting RNA-binding proteins as translation fac-
tors to engineer gene circuits [15]. In this work, we focused on quantitatively studying the sto-
chastic behavior of this type of circuits (i.e., noise generation and propagation in gene expression
when it is regulated at the level of translation). Our results show that the protein-RNA interaction
in this case leads to a significant down-regulation in expression of about 50-fold by blocking the
progression of the ribosome on the target mRNA, which is comparable to transcriptional fold-
changes. In addition, a general mathematical framework was shown suitable to describe the sto-
chastic behavior in regulations exerted by both transcription and translation factors. Noise propa-
gation from gene to gene is buffered when the regulator acts at the level of translation, as the
amplification process by transcription is avoided, and a Gamma distribution properly
Fig 4. Detailed analysis of stochastic gene expression modulated by growth rate. a) As TC increases, cells grow slower and bigger. b) Mean eBFP2 expression as a
function of growth rate for each IPTG condition. c) Mean sfGFP expression as a function of growth rate for each IPTG condition. d) Relative mean eBFP2 synthesis
rate as a function of growth rate for each IPTG condition. e) Relative mean sfGFP synthesis rate as a function of growth rate for each IPTG condition. In d,e), the
values are relative to the case TC = 0, and the dashed line corresponds to a linear dependence (which comes from a null model in which the translation rate is not
affected by the growth rate). f) Transfer function of the post-transcriptional regulation in terms of mean expression for each TC condition. The inset shows the effect of
TC on the sfGFP mRNA synthesis rate (included in αy, as the growth rate decreases) and the effective dissociation constant between eBFP2-MS2CP and sfGFP mRNA
(θx, as the volume increases) according to the mathematical model. g) eBFP2 noise as a function of growth rate for each IPTG condition. h) sfGFP noise as a function
of growth rate for each IPTG condition. i) Transfer function of the post-transcriptional regulation in terms of noise for each TC condition. In plots b-i), points
correspond to calculations from the experimental data, while solid lines to predictions with the mathematical model.
https://doi.org/10.1371/journal.pcbi.1010087.g004
parameterized can provide deep analytical explanations about the resulting cell-to-cell variability
[28]. By modulating the cellular growth rate, we also reported an interplay between global and
local regulatory mechanisms in the cell that affect both the mean expression and noise levels.
It is important to notice that the growth rate that we measured here corresponds to an average
of the population (as we calculate it having monitored absorbance with time for a culture). Never-
theless, each cell grows differently, especially when the culture is under the effect of TC. In this
regard, a mathematical model incorporating such heterogeneous growth might explain better the
observed noise patterns [37]. Another limitation of our work is the assumption that the noise
sources are independent, which allowed us to derive a compact mathematical expression for the
noise of sfGFP. This is not strictly true since eBFP2-MS2CP and sfGFP share extrinsic noise
sources [7]. Yet, we expect reliability in the conclusions derived from this study.
Here, we exploited the viral protein MS2CP to implement the regulatory system, but in
principle other RNA-binding proteins might be used. For example, the bacteriophage PP7
coat protein or the Mycobacterium enzyme PyrR [16] are suitable elements from which to
engineer orthogonal systems. In fact, given the plethora of RNA-binding proteins in nature,
especially in eukaryotes [38], and noting that the regulatory mechanism only requires a tight
protein-RNA interaction to interfere with the ribosome, multiple implementations might be
achieved. Our mathematical model is general enough to describe these eventual implementa-
tions. We only expect to change the kinetic parameters for each particular protein, preserving
the functional form. In principle, each RNA motif will lead to a different translation rate. In
this regard, the predictability of the system might be strengthened by using the RBS calculator
[39]. We also anticipate that the use of tandem repeats of the RNA motif might enhance the
regulatory fold-change of the system. Furthermore, since RNA is a very versatile molecule,
RNA-binding proteins can regulate gene expression through a variety of mechanisms acting
post-transcriptionally, including the regulation of translation initiation, translation elongation,
transcription termination, and RNA stability [40]. In prokaryotes, the regulation of translation
initiation by controlling RBS accessibility is a widespread mechanism, but in eukaryotes the
blockage of translation elongation has been observed in the case of Argonaute proteins [41].
Arguably, a blockage in the initial phase by MS2CP is key in our synthetic system. Further
work should analyze how those other mechanisms generate and propagate noise.
In sum, our work provides new quantitative insights about the stochastic behavior in
genes regulated translationally by RNA-binding proteins. Protein translation factors can
integrate some advantages distinctively attributed to proteins (as transcription factors), such
as the ability to transduce small signals and achieve large dynamic ranges, or to small RNAs,
such as the ability to produce rapid responses and buffer transcriptional noise [27]. Further-
more, our work paves the way for engineering gene regulatory circuits with greater integra-
bility and then sophistication. Certainly, the combination of different layers within the
genetic information flow (i.e., transcription and translation) leads to an easier integration of
signals to achieve a given function [42]. Therefore, we envision that RNA-binding proteins
will be of great utility in synthetic biology in the close future to face biotechnological and bio-
medical challenges.
the translational control of the MS2CP-recognizing RNA motif. The genetic cassettes were
synthesized by IDT. LB medium was used for both overnight and characterization cultures.
Kanamycin and chloramphenicol were used at the concentration of 50 μg/mL and 34 μg/mL,
respectively. IPTG and TC were used as inducers of the system. The concentration gradient of
IPTG that we tested was 0, 5, 10, 20, 50, 100, 200, 500, and 1000 μM, and the concentration
gradient of TC was 0, 10, 100, 200, 300, 400, 500, 700, and 1000 ng/mL. Compounds provided
by Sigma.
Note that LacI is overexpressed in MG1655-Z1 to efficiently regulate the PLlac promoter in
the pRKFR2 plasmid. This overexpression adds to the wild-type expression of LacI. In addi-
tion, TetR is expressed in E. coli MG1655-Z1 (although it does not play any regulatory role), so
it will bind to TC when this inducer is used. The titration effect will be more relevant at low
concentrations of TC, although overall this will only mean an effective TC concentration
slightly lower. The maximal TC concentration used here reduced significantly the growth rate
of the cells, suggesting a marginal effect of TetR in this case.
Growth curves
Cultures (2 mL) inoculated from single colonies (three replicates) were grown overnight in
LB medium at 37˚C and 200 rpm. Cultures were then diluted 1:100 in fresh LB medium (2
mL) and were grown for 3 h at the same conditions to reach exponential phase (OD600
around 0.5). Cultures were then diluted 1:50 in fresh LB medium (200 μL) to load a micro-
plate (96 wells, black, clear bottom; Corning) with appropriate concentrations of IPTG and
TC. The microplate was then incubated for 10 h at 37˚C and 1,000 rpm in a PST-60HL
plate shaker (Biosan). Absorbance (600 nm) was measured every hour in a Varioskan Lux
fluorometer (Thermo). The growth rate was calculated as the slope between absorbance (in
log scale) and time during the exponential phase. Data analysis performed with MATLAB
(MathWorks) and Python.
Flow cytometry
Cultures (2 mL) inoculated from single colonies (four replicates) were grown overnight in LB
medium at 37˚C and 200 rpm. Cultures were then diluted 1:100 in fresh LB medium (2 mL)
and were grown for 3 h at the same conditions to reach exponential phase (OD600 around 0.5).
Cultures were then diluted 1:50 in fresh LB medium (200 μL) to load a microplate (96 wells,
black, clear bottom; Corning) with appropriate concentrations of IPTG and TC. The micro-
plate was then incubated at 37˚C and 1,000 rpm in a PST-60HL plate shaker (Biosan) until cul-
tures reached a sufficient OD600 (a different incubation time for each TC concentration).
Cultures (6 μL) were then diluted in PBS (1 mL). Fluorescence was measured in an LSRFor-
tessa flow cytometer (BD); a 405 nm laser and a 450 nm filter for blue fluorescence, and a 488
nm laser and a 530 nm filter for green fluorescence. Events were gated by using the forward
and side scatter signals and compensated (~104 events after this process). The mean value of
the autofluorescence of the cells was subtracted in each channel to obtain a final estimate of
expression (S1 Data). Data analysis performed with MATLAB and Python. The mean and the
variance were calculated for each distribution after removing outliers, which served to com-
pute the noise in gene expression.
molecules of control, the Hill-Langmuir equations that dictate average protein expression
(population measure) are
� �n
IPTG i
rx þ
y
heBFP2i ¼ ax � i �ni
IPTG
1þ
yi
� �n ð1Þ
heBFP2i x
1 þ ry
yx
hsfGFPi ¼ ay � �n ;
heBFP2i x
1þ
yx
where αx is the maximal protein level from the PLlac promoter (in presence of IPTG), ρx the
transcriptional repression fold by LacI, αy the maximal protein level from the constitutive
J23119 promoter, and ρy the translational repression fold by eBFP2-MS2CP. In addition, θi is
the effective dissociation constant between LacI and IPTG, ni the effective degree of coopera-
tivity of LacI, θx the effective dissociation constant between eBFP2-MS2CP and the cognate
RNA motif embedded within the sfGFP mRNA, and nx the effective degree of cooperativity of
eBFP2-MS2CP. By using our data, the adjusted parameter values are αx = 1,890 AU, ρx =
0.225, θi = 116 μM, ni = 2.38, αy = 29,000 AU, ρy = 0.016, θx = 610 AU, and nx = 5 (upon vary-
ing IPTG, with no TC).
Moreover, if μ denotes the actual cell growth rate, which is modulated by TC and
dictates the dilution rate of the proteins, it turns out that the following Michaelis-Menten
equation
m0
m¼ ; ð2Þ
1 þ TC
yc
where μ0 is the maximal cell growth rate (in absence of TC) and θc the half maximal inhibi-
tory concentration of TC. From the data, we obtained μ0 = 1.2 h-1 and θc = 526 ng/mL
(upon varying TC, with no IPTG).
To consider the impact of growth rate on protein expression, we noticed that protein
expression is the ratio between the protein synthesis rate (SR) and the growth rate, and that
the protein synthesis rate is the product between the mRNA amount and the translation rate
(λi, for gene i). That is,
hSReBFP2 i lx hmRNAeBFP2 i
heBFP2i ¼ ¼
m m
ð3Þ
hSRsfGFP i ly hmRNAeBFP2 i
hsfGFPi ¼ ¼ :
m m
Because it is known that the mRNA amount per cell is proportional to the growth rate (i.e., the
transcription rate increases as long as the cell grows faster, hmRNAi/μ) [33], the dependence
of the protein synthesis rate on the growth rate is just given by the effect of the growth rate on
the translation rate. Interestingly, the translation rate can be described by a Michaelis function
of the number of active ribosomes [32], which is known to increase linearly when the growth
rate decreases as a consequence of TC. We can then define two effective parameters (ε1i and
Here, we adjusted ε1x = 0.196 h-1, ε1x = 0.299 h-1, ε1y = 0.246 h-1, and ε1y = 0.349 h-1.
In addition, we noticed that the cellular volume changes with the growth rate as a conse-
quence of TC. By using the cube of the forward scatter signal (median of the population) as a
proxy of the volume, a negative exponential trend was identified (i.e., the volume increases
with TC). Because the model considers as variables the amount of fluorescent proteins per cell
and not the concentrations, the parameter θx needs to be corrected as
yx ðmÞ ¼ edðm0 mÞ
yx ðm0 Þ; ð5Þ
where here we took δ = 0.85 h. In essence, when the volume increases, the number of proteins
required to regulate the target gene are higher. Besides, the impact of growth rate on the tran-
scriptional regulation was assumed negligible, as LacI is a highly expressed protein from the
chromosome whose activity is modulated by a chemical inducer, which enters into the model
in terms of concentration (i.e., θi independent of μ).
0 � �nx 1 12
heBFP2i
by 1B y a ð1 r y Þn x yx C gy heBFP2i
2
where Z2x and Z2y are two empirical constants that quantify the levels of noise of extrinsic nature
on eBFP2-MS2CP and sfGFP, respectively. Also, Z2lac is a constant that measures the noise in
LacI expression, βx and βy the Fano factors of noise of intrinsic nature for eBFP2-MS2CP and
sfGFP, respectively, and γy a constant that accounts for the difference between fluorescence
and number of molecules. Note that heBFP2i = hMS2CPi and CV2eBFP2 ¼ CV2MS2CP . Note also
that when a strong repression occurs at the level of translation the transcriptional noise can be
neglected and then βy can be considered constant (i.e., the Fano factor, in number of molecules
per cell, can be approached by 1). By using our data, the adjusted parameter values are Z2x ¼
0:246; βx = 45.6 AU, Z2lac ¼ 6; 470 μM2, Z2y ¼ 0:127; βy = 61.9 AU, and γy = 0.0233 (upon vary-
ing IPTG, with no TC).
It is important to recall that the parameters αx, αy, and θx depend on the growth rate. More-
over, the intrinsic noise Fano factor for eBFP2-MS2CP is proportional to the translation rate
(βx/λx),
lx ðmÞ
b x ðm Þ ¼ b ðm Þ; ð7Þ
lx ðm0 Þ x 0
but not the factor for sfGFP (βy nearly independent of λy), as sfGFP is regulated at the level of
translation and then its mRNA is constitutively expressed. To explain our data, we introduced
the following phenomenological expressions
� �2
2
m
Zlac ðmÞ ¼ Z2lac ðm0 Þ
m0
lx ðmÞ 2
Z2x ðmÞ ¼ Z ðm Þ ð8Þ
lx ðm0 Þ x 0
ly ðmÞ 2
Z2y ðmÞ ¼ Z ðm Þ:
ly ðm0 Þ y 0
In essence, this indicates that LacI expression varies from cell to cell in greater extent when
cells grow faster, and that the extrinsic noise increases when the growth rate is very low.
Finally, we assumed that the stochastic gene expression follows a Gamma distribution (see S2
Appendix for a basic derivation) [28]. Then, the probability for a given expression level reads
eBFP2ax 1 e eBFP2=bx
PðeBFP2Þ ¼
Gðax Þbx ax
ð9Þ
sfGFPay 1 e sfGFP=by
PðsfGFPÞ ¼ ;
Gðay Þby ay
where ax and ay are the Gamma shape parameters for eBFP2-MS2CP and sfGFP, respectively,
and bx and by the Gamma scale parameters. Importantly, by knowing that for a Gamma distribu-
tion ab is the mean and ab2 the variance, these two parameters can be defined as
1
ax ¼
CV2eBFP2
1
ay ¼ ð10Þ
CV2sfGFP
bx ¼ heBFP2iCV2eBFP2
by ¼ hsfGFPiCV2sfGFP :
This means that the Gamma shape parameter is directly the inverse of the noise, and that the
Gamma scale parameter depends on the translation rate in the case of eBFP2-MS2CP (transcrip-
tion regulation) and is nearly independent of it in the case of sfGFP (translation regulation).
Numerical simulations
The system of stochastic differential equations (see S1 Appendix) was solved numerically to
obtain stochastic trajectories of mRNA and protein concentrations. For that, we followed an inte-
gration scheme previously described [43]. The colored stochastic processes (for extrinsic and reg-
ulation noise) were obtained from independent white stochastic processes. The system was solved
in one time interval with the routine ode45s from MATLAB, considering constant the stochastic
fluctuations in that interval. The values of the fluctuations were updated in each interval with the
previous mRNA and protein concentrations. Negative concentration values were avoided.
Supporting information
S1 Fig. Reliability of the dose-response curve. a) Mean of eBFP2 expression as a function of
IPTG. b) Mean of sfGFP expression as a function of IPTG. c) Noise of eBFP2 expression as a
function of IPTG. d) Noise of sfGFP expression as a function of IPTG. Points correspond to
the values of the population shown in the main figures. Error bars correspond to standard
errors calculated from four different populations. Solid lines correspond to predictions with
the mathematical model.
(TIF)
S2 Fig. Numerical simulations of stochastic dynamics. a-d) Stochastic trajectories with time
of eBFP2 and sfGFP for two different IPTG concentrations. In red, deterministic trajectories.
The initial condition corresponds to the uninduced state in all cases. e-h) Histograms of pro-
tein expression computed from long trajectories. The Gamma distributions fitted against the
experimental data (blue lines) were also represented.
(TIF)
S3 Fig. Sensitivity analysis of the model parameters. Plots of mean and noise of expression
as a function of IPTG, where solid lines correspond to the dynamics predicted with the
adjusted parameter, dotted lines to the dynamics if the parameter increases 2-fold, and dashed
lines to the dynamics if the parameter decreases 2-fold.
(TIF)
S4 Fig. Stochastic gene expression described by a Gamma distribution. Histograms of
experimental single-cell fluorescence for both a) eBFP2 and b) sfGFP for different induction
conditions with IPTG, together with fitted Gamma distributions against the data (blue lines)
and predicted Gamma distributions obtained by using the model values of mean and noise
(red lines).
(TIF)
S5 Fig. Growth curves. Three different populations (blue, red, and green) were monitored
with time. Points correspond to absorbance values, while solid lines come from fitted exponen-
tial trends.
(TIF)
S6 Fig. Relationship between cellular growth rate and volume. a) Schematics to show that as
TC increases, cells grow slower and are bigger. b) Scatter plot between the cube of the forward
scattering signal (proxy of cellular volume) and the growth rate for the 81 IPTG and TC condi-
tions (colored by TC condition). An exponential trend was adjusted (solid line).
(TIF)
Author Contributions
Conceptualization: Guillermo Rodrigo.
Investigation: Roswitha Dolcemascolo, Lucas Goiriz, Roser Montagud-Martı́nez, Guillermo
Rodrigo.
Methodology: Roswitha Dolcemascolo, Lucas Goiriz, Roser Montagud-Martı́nez, Guillermo
Rodrigo.
Supervision: Roser Montagud-Martı́nez, Guillermo Rodrigo.
Visualization: Lucas Goiriz, Guillermo Rodrigo.
Writing – original draft: Roswitha Dolcemascolo, Lucas Goiriz, Guillermo Rodrigo.
References
1. Wagner GP, Zhang J (2011) The pleiotropic structure of the genotype-phenotype map: the evolvability
of complex organisms. Nat Rev Genet 12: 204–213. https://doi.org/10.1038/nrg2949 PMID: 21331091
2. Balázsi G, van Oudenaarden A, Collins JJ (2011) Cellular decision-making and biological noise: from
microbes to mammals. Cell 144: 910–925. https://doi.org/10.1016/j.cell.2011.01.030 PMID: 21414483
3. Swain PS, Elowitz MB, Siggia ED (2002) Intrinsic and extrinsic contributions to stochasticity in gene
expression. Proc Natl Acad Sci USA 99: 12795–12800. https://doi.org/10.1073/pnas.162041399
PMID: 12237400
4. Fraser HB, Hirsh AE, Giaever G, Kumm J, Eisen MB (2004) Noise minimization in eukaryotic gene
expression. PLoS Biol 2: e137. https://doi.org/10.1371/journal.pbio.0020137 PMID: 15124029
5. Eldar A, Elowitz MB (2010) Functional roles for noise in genetic circuits. Nature 467: 167–173. https://
doi.org/10.1038/nature09326 PMID: 20829787
6. Beckett D (2001) Regulated assembly of transcription factors and control of transcription initiation. J
Mol Biol 314: 335–352. https://doi.org/10.1006/jmbi.2001.5134 PMID: 11846548
7. Pedraza JM, van Oudenaarden A (2005) Noise propagation in gene networks. Science 307: 1965–
1969. https://doi.org/10.1126/science.1109090 PMID: 15790857
8. Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S (2006) Stochastic mRNA synthesis in mammalian
cells. PLoS Biol 4: e309. https://doi.org/10.1371/journal.pbio.0040309 PMID: 17048983
9. Libby E, Perkins TJ, Swain PS (2007) Noisy information processing through transcriptional regulation.
Proc Natl Acad Sci USA 104: 7151–7156. https://doi.org/10.1073/pnas.0608963104 PMID: 17420464
10. Dunlop MJ, Cox RS, Levine JH, Murray RM, Elowitz MB (2008) Regulatory activity revealed by dynamic
correlations in gene expression noise. Nat Genet 40: 1493–1498. https://doi.org/10.1038/ng.281
PMID: 19029898
11. Siciliano V, Garzilli I, Fracassi C, Criscuolo S, Ventre S, di Bernardo D (2013) MiRNAs confer pheno-
typic robustness to gene networks by suppressing biological noise. Nat Commun 4: 2364. https://doi.
org/10.1038/ncomms3364 PMID: 24077216
12. Schmiedel JM, Klemm SL, Zheng Y, Sahay A, Blüthgen N, Marks DS, van Oudenaarden A (2015)
MicroRNA control of protein expression noise. Science 348: 128–132. https://doi.org/10.1126/science.
aaa1738 PMID: 25838385
13. Dacheux E, Malys N, Meng X, Ramachandran V, Mendes P, McCarthy JE (2017) Translation initiation
events on structured eukaryotic mRNAs generate gene expression noise. Nucleic Acids Res 45: 6981–
6992. https://doi.org/10.1093/nar/gkx430 PMID: 28521011
14. Rodrigo G (2018) Post-transcriptional bursting in genes regulated by small RNA molecules. Phys Rev E
97: 032401. https://doi.org/10.1103/PhysRevE.97.032401 PMID: 29776125
15. Katz N, Cohen R, Solomon O, Kaufmann B, Atar O, Yakhini Z, Goldberg S, Amit R (2019) Synthetic 5’
UTRs can either up- or downregulate expression upon RNA-binding protein binding. Cell Syst 9: 93–
106. https://doi.org/10.1016/j.cels.2019.04.007 PMID: 31129060
16. Babitzke P, Baker CS, Romeo T (2009) Regulation of translation initiation by RNA binding proteins.
Annu Rev Microbiol 63: 27–44. https://doi.org/10.1146/annurev.micro.091208.073514 PMID:
19385727
17. Bertrand E, Chartrand P, Schaefer M, Shenoy SM, Singer RH, Long RM (1998) Localization of ASH1
mRNA particles in living yeast. Mol Cell 2: 437–445. https://doi.org/10.1016/s1097-2765(00)80143-4
PMID: 9809065
18. Graindorge A, Pinheiro I, Nawrocka A, Mallory AC, Tsvetkov P, Gil N, Carolis C, Buchholz F, Ulitsky I,
Heard E, Taipale M, Shkumatava A (2019) In-cell identification and measurement of RNA-protein inter-
actions. Nat Commun 10: 5317. https://doi.org/10.1038/s41467-019-13235-w PMID: 31757954
19. Zalatan JG, Lee ME, Almeida R, Gilbert LA, Whitehead EH, La Russa M, Tsai JC, Weissman JS, Due-
ber JE, Qi LS, Lim WA (2015) Engineering complex synthetic transcriptional programs with CRISPR
RNA scaffolds. Cell 160: 339–350. https://doi.org/10.1016/j.cell.2014.11.052 PMID: 25533786
20. Delebecque CJ, Lindner AB, Silver PA, Aldaye FA (2011) Organization of intracellular reactions with
rationally designed RNA assemblies. Science 333: 470–474. https://doi.org/10.1126/science.1206938
PMID: 21700839
21. Lutz R, Bujard H (1997) Independent and tight regulation of transcriptional units in Escherichia coli via
the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements. Nucleic Acids Res 25: 1203–1210. https://
doi.org/10.1093/nar/25.6.1203 PMID: 9092630
22. Ai HW, Shaner NC, Cheng Z, Tsien RY, Campbell RE (2007) Exploration of new chromophore struc-
tures leads to the identification of improved blue fluorescent proteins. Biochemistry 46: 5904–5910.
https://doi.org/10.1021/bi700199g PMID: 17444659
23. Pédelacq JD, Cabantous S, Tran T, Terwilliger TC, Waldo GS (2006) Engineering and characterization
of a superfolder green fluorescent protein. Nat Biotechnol 24: 79–88. https://doi.org/10.1038/nbt1172
PMID: 16369541
24. Lim F, Peabody DS (1994) Mutations that increase the affinity of a translational repressor for RNA.
Nucleic Acids Res 22: 3748–3752. https://doi.org/10.1093/nar/22.18.3748 PMID: 7937087
25. Volfson D, Marciniak J, Blake WJ, Ostroff N, Tsimring LS, Hasty J (2006) Origins of extrinsic variability
in eukaryotic gene expression. Nature 439: 861–864. https://doi.org/10.1038/nature04281 PMID:
16372021
26. Kaern M, Elston TC, Blake WJ, Collins JJ (2005) Stochasticity in gene expression: from theories to phe-
notypes. Nat Rev Genet 6: 451–464. https://doi.org/10.1038/nrg1615 PMID: 15883588
27. Mehta P, Goyal S, Wingreen N (2008) A quantitative comparison of sRNA-based and protein-based
gene regulation. Mol Syst Biol 4: 221. https://doi.org/10.1038/msb.2008.58 PMID: 18854820
28. Taniguchi Y, Choi PJ, Li GW, Chen H, Babu M, Hearn J, Emili A, Xie XS (2010) Quantifying E. coli pro-
teome and transcriptome with single-molecule sensitivity in single cells. Science 329: 533–538. https://
doi.org/10.1126/science.1188308 PMID: 20671182
29. Friedman N, Cai L, Xie XS (2006) Linking stochastic dynamics to population distribution: an analytical
framework of gene expression. Phys Rev Lett 97: 168302. https://doi.org/10.1103/PhysRevLett.97.
168302 PMID: 17155441
30. Wilson DN (2014) Ribosome-targeting antibiotics and mechanisms of bacterial resistance. Nat Rev
Microbiol 12: 35–48. https://doi.org/10.1038/nrmicro3155 PMID: 24336183
31. Scott M, Gunderson CW, Mateescu EM, Zhang Z, Hwa T (2010) Interdependence of cell growth and
gene expression: origins and consequences. Science 330: 1099–1102. https://doi.org/10.1126/
science.1192588 PMID: 21097934
32. Borkowski O, Goelzer A, Schaffer M, Calabre M, Mäder U, Aymerich S, Jules M, Fromion V (2016)
Translation elicits a growth rate-dependent, genome-wide, differential protein production in Bacillus
subtilis. Mol Syst Biol 12: 870. https://doi.org/10.15252/msb.20156608 PMID: 27193784
33. Klumpp S, Zhang Z, Hwa T (2009) Growth rate-dependent global effects on gene expression in bacte-
ria. Cell 139: 1366–1375. https://doi.org/10.1016/j.cell.2009.12.001 PMID: 20064380
34. Keren L, van Dijk D, Weingarten-Gabbay S, Davidi D, Jona G, Weinberger A, Milo R, Segal E (2015)
Noise in gene expression is coupled to growth rate. Genome Res 25: 1893–1902. https://doi.org/10.
1101/gr.191635.115 PMID: 26355006
35. Harris LK, Theriot JA (2016) Relative rates of surface and volume synthesis set bacterial cell size. Cell
165: 1479–1492. https://doi.org/10.1016/j.cell.2016.05.045 PMID: 27259152
36. Urchueguı́a A, Galbusera L, Chauvin D, Bellement G, Julou T, van Nimwegen E (2021) Genome-wide
gene expression noise in Escherichia coli is condition-dependent and determined by propagation of
noise through the regulatory network. PLoS Biol 19: e3001491. https://doi.org/10.1371/journal.pbio.
3001491 PMID: 34919538
37. Ray JCJ, Wickersheim ML, Jalihal AP, Adeshina YO, Cooper TF, Balázsi G (2016) Cellular growth
arrest and persistence from enzyme saturation. PLoS Comput Biol 12: e1004825. https://doi.org/10.
1371/journal.pcbi.1004825 PMID: 27010473
38. Beckmann BM, Horos R, Fischer B, Castello A, Eichelbaum K, Alleaume AM, Schwarzl T, Curk T,
Foehr S, Huber W, Krijgsveld J, Hentze MW (2015) The RNA-binding proteomes from yeast to man har-
bour conserved enigmRBPs. Nat Commun 6: 10127. https://doi.org/10.1038/ncomms10127 PMID:
26632259
39. Salis HM, Mirsky EA, Voigt CA (2009) Automated design of synthetic ribosome binding sites to control
protein expression. Nat Biotechnol 27: 946–50. https://doi.org/10.1038/nbt.1568 PMID: 19801975
40. van Assche E, van Puyvelde S, Vanderleyden J, Steenackers PH (2015) RNA-binding proteins involved
in post-transcriptional regulation in bacteria. Front Microbiol 6: 141. https://doi.org/10.3389/fmicb.2015.
00141 PMID: 25784899
41. Iwakawa HO, Tomari Y (2013) Molecular insights into microRNA-mediated translational repression in
plants. Mol Cell 52: 591–601. https://doi.org/10.1016/j.molcel.2013.10.033 PMID: 24267452
42. Rosado A, Cordero T, Rodrigo G (2018) Binary addition in a living cell based on riboregulation. PLoS
Genet 14: e1007548. https://doi.org/10.1371/journal.pgen.1007548 PMID: 30024870
43. Rodrigo G, Carrera J, Jaramillo A (2011) Computational design of synthetic regulatory networks from a
genetic library to characterize the designability of dynamical behaviors. Nucleic Acids Res 39: e138.
https://doi.org/10.1093/nar/gkr616 PMID: 21865275