Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)
Article
Deconvolution of gene expression noise into spatial dynamics of transcription factor-promoter interplay Angel Goñi-Moreno, Ilaria Benedetti, Juhyun Kim, and Victor de Lorenzo ACS Synth. Biol., Just Accepted Manuscript • DOI: 10.1021/acssynbio.6b00397 • Publication Date (Web): 29 Mar 2017 Downloaded from http://pubs.acs.org on March 30, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
ACS Synthetic Biology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
39x19mm (300 x 300 DPI)
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2 3
Deconvolution of gene expression noise into spatial dynamics of transcription factorpromoter interplay
4
by
5 Ángel Goñi-Moreno1,$, Ilaria Benedetti1, Juhyun Kim1, Víctor de Lorenzo1,*
6 7 8
1
Systems Biology Program, Centro Nacional de Biotecnología, Cantoblanco-Madrid, Spain.
9 10 11
Keywords:
12 13
Transcriptional noise, bet hedging, Pseudomonas putida, intracellular heterogeneity, TOL plasmid
Running title:
Intracellular microgranularity and intrinsic noise
14 15 16 17 18 19
*Correspondence to:
Víctor de Lorenzo (
[email protected])
20
Centro Nacional de Biotecnología (CNB-CSIC)
21
Campus de Cantoblanco, 28049 Madrid (Spain)
22
Tel: (+34 91) 585 45 73; Fax: (+34 91) 585 45 06
23
______________________________________________________________________________
24
$
25 26
United Kingdom
Present address: School of Computing Science, Newcastle University, Newcastle upon Tyne,
1
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
27
Abstract
28
Gene expression noise is not only the mere consequence of stochasticity, but also a signal that
29
reflects the upstream physical dynamics of the cognate molecular machinery. Soil bacteria facing
30
recalcitrant pollutants exploit noise of catabolic promoters to deploy beneficial phenotypes such
31
as metabolic bet-hedging and/or division of biochemical labour. Although the role of upstream
32
promoter-regulator interplay in the origin of this noise is little understood, its specifications are
33
probably ciphered in flow cytometry data patterns. We studied Pm promoter activity of the
34
environmental bacterium Pseudomonas putida and its cognate regulator XylS by following
35
expression of Pm-gfp fusions in single cells. Using mathematical modelling and computational
36
simulations, we determined the kinetic properties of the system and used them as a baseline code
37
to interpret promoter activity in terms of upstream regulator dynamics. Transcriptional noise was
38
predicted to depend on the intracellular physical distance between regulator source (where XylS
39
is produced) and the target promoter. Experiments with engineered bacteria in which this distance
40
is minimised or enlarged confirmed the predicted effects of source/target proximity on noise
41
patterns. This approach allowed deconvolution of cytometry data into mechanistic information on
42
gene expression flow. It also provided a basis for selecting programmable noise levels in
43
synthetic regulatory circuits.
44
Introduction
45
Information processing inside bacterial cells in response to physicochemical stimuli requires
46
regulatory cascades to propagate input/output signals effectively. This process typically involves
47
several steps in which transcription factors (TF) interact with promoters to trigger gene
48
expression responses. These interactions occur stochastically rather than deterministically1-5,
49
leading to specific and variable noisy signals. The customary view considers this effect the
50
necessary result of random fluctuations of regulatory elements present in short supply in
51
individual cells6. The dynamic properties of promoter activation have a determining influence on
52
expression noise range and intensity7-9.
53
Single-cell technologies10-12 has helped to clarify various mechanisms behind noise generation. 2
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
54
A major noise source in virtually every prokaryotic promoter is the so-called bursting effect13, 14,
55
a pulse-like activity that results largely from discontinuous topological changes in DNA caused
56
by RNA polymerase progression through DNA15, 16. Different in vivo noise generators can be
57
measured by fluorescence distribution in single cells, which can be followed by cytometry17.
58
Cell cytometry profiles bear embedded information on the mechanistic origin of the gene
59
expression noise, but it remains unclear how these data can be retraced to the physical TF-
60
promoter dynamics that produces this fluorescence distribution.
61
The Gram-negative soil bacterium Pseudomonas putida mt-2 provides an excellent model for
62
responding to these questions. This microorganism can thrive in sites polluted with aromatic
63
chemicals18 such as m-xylene (m-xyl), because of a complex metabolic and regulatory network
64
encoded in its single-copy TOL plasmid pWW019 (Figure 1A). The noise from each of the four
65
promoters seems to be precisely controlled, giving rise to metabolic diversification20. This allows
66
a fraction of the cells in a population to explore new nutritional landscapes without risking
67
communal collapse20, 21.
68 69
The noise of the TOL network Pu and Ps promoters is explained by the small number of
70
molecules of their regulatory protein XylR22; however, that of Pm, which controls the lower
71
operon of the TOL pathway in response to 3-methylbenzoate23 (3MBz) is puzzling. This
72
promoter can be activated by two separate mechanisms, [i] low intracellular concentration of its
73
regulator, XylS, bound to its effector, 3MBz, or [ii] m-xyl-induced XylS overproduction with no
74
3MBz involvement (Figure 1A). When these conditions coincide (i.e., high XylS levels and
75
presence of 3MBz), Pm activity is very high24, 25. Therefore, Pm function is that of an OR logic
76
gate (Figure 1B) where either input, 3MBz or m-xyl, can trigger its activity. The puzzling feature
77
of the Pm regulatory node is that the Pm noise pattern varies greatly depending on the induction
78
mode. We analysed whether cell cytometry data for transcriptional Pm-gfp fusions could be
79
decoded into information on the physical dynamics of promoter activation, including clues to
80
XylS/Pm regulatory node spatial arrangement.
81 82
A combination of modelling and wet experiments shows these noise regimes can be changed by 3
ACS Paragon Plus Environment
Page 4 of 32
Page 5 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
83
regulator numbers and the spatial arrangement of regulatory components. These predictions were
84
validated in cells engineered to minimise the distance between the XylS source site and Pm
85
location.
86 87
Results
88 89
Two distinct noise regimes define Pm output. The activity of the TOL plasmid promoter Pm
90
(Figure 1) can be induced by exposing P. putida mt-2 to one of two inputs. The first input, 3MBz,
91
activates the XylS molecules leaked from a non-active Ps promoter. The second, m-xyl, is
92
converted metabolically inside cells by the upper TOL pathway enzymes to 3MBz. In addition,
93
the latter scenario leads to XylS overproduction since m-xyl triggers the activity of the Ps
94
promoter via the XylR regulator. Therefore, the input m-xyl leads to a higher concentration of
95
activated XylS molecules (XylSa) than 3MBz (Figure 1C). Although the Pm/XylS node of the
96
TOL network is often abstracted as a binary switch with only ON/OFF states, the unique noisy
97
nature of Pm output invalidates this view and highlights the role of signal variability26.
98 99
The P. putida mt-2-Pm, strain derives from the natural P. putida mt-2 isolate bearing the TOL
100
plasmid pWW0, but was engineered to bear a transcriptional Pm-gfp fusion in its chromosome
101
(Table 1 and Fig. 2A). This strain has its sole XylS source in the pWW0-encoded gene, which is
102
expressed through the TOL plasmid Ps promoter (see Figure 1 and Methods), and the XylS target,
103
the encoded Pm-gfp, in the chromosome. Because of this arrangement, the DNA region that
104
encodes and supplies XylS transcriptional regulators is not adjacent to the target promoter, from
105
which it is physically separated in trans. Differences in Pm expression regimes depend on
106
whether the XylS/Pm regulatory node is induced with m-xyl or 3MBz (Figure 2A), as shown by
107
flow cytometry results of promoter activity (Figure 2B, C). Induction with m-xyl (XylS
108
overproduction and intracellular 3MBz production) leads to a situation in which the output signal
109
noise range allows a null overlap between ON (uninduced) and OFF (induced) states. In contrast,
110
induction with exogenous 3MBz (low XylS) left the OFF state unaffected while producing an ON
111
state in which the noise regime consisted of a broad, plateau-like distribution from the lowest to
112
the highest intensity value. The output ON signals are patently different in both cases, suggesting 4
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
113
that they originate in a different type of TF-promoter interplay beyond randomness. We measured
114
the noise produced by leaky XylS expression (no Ps activity, Figure 2C) and that of the full XylS
115
production (Ps activity, Figure 2B). The former scenario, while very noisy, does not depend of
116
upstream nodes. The latter does depend on upstream nodes, but does not generate noisy patterns.
117
Altogether, we can assume that upstream Pm nodes are not a source of noise in our setup.
118
Therefore, we focus our study in Pm-gfp kinetics and TF-promoter dynamics.
119 120
Noise deconvolution and rate optimisation. The regulatory node formed by the pair XylS-Pm
121
can be modelled according to the kinetic rates in Figure 3A. The regulator in its active form,
122
XylSa, binds Pm (k1, molecules-1hour-1) to fire its activity and unbinds it (k-1, hour-1) back to the
123
default silent state. When bound, mRNA molecules are transcribed (k2, hour-1) from the
124
downstream gene, gfp, which produce proteins through translation (k3, hour-1). Even when the
125
regulator is not bound, there is some leakage of basal Pm transcription (k6, hour-1). To complete
126
the model, we included degradation rates for mRNA (k4, hour-1) and GFP (k5, hour-1).
127 128
Our goal was to find those values for the rates that allow Pm to produce different noise regimes
129
depending on the inducer used. To this end, we considered a training vector θ with the rates
130
mainly responsible for Pm dynamics, defined as follows: θ = (k1, k-1, k2, k3, [XylSa]), where [x]
131
(molecules) denotes number of molecule x. Basal activity (k6 = 15) and degradations (k4 = 10, k5
132
= 2) were specified within standard ranges for mathematical analysis27-31 (see Methods). To
133
identify the set of values that best simulated the experimental output, we defined two fitness
134
parameters (f) based on the ON state produced by 3MBz (Figure 2C), wide-range signal (f1) and
135
plateau-like surface (f2) (Figure 3B). An optimisation process (see Methods) yielded vector θf =
136
(0.004, 1.5, 900, 80, 200). Strong expression kinetics (transcription + translation) are needed to
137
produce high signal intensity (requirement for f1), while low binding/unbinding rates guarantees
138
affinity instability and thus helps generate the final plateau-shaped distribution (f2). These values
139
produce the broadest range signal possible while assuring flatness which resembles the ON state
140
induced by 3MBz (Figure 2B). To test whether the noise regime observed in the ON state during
141
m-xyl induction could be reproduced, we carried out stochastic simulations (see Methods) in
142
which the number of XylSa molecules was increased while vector θf rates were untouched; this 5
ACS Paragon Plus Environment
Page 6 of 32
Page 7 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
143
concentration was fixed at [XylSa] = 3000 molecules, which reproduced the m-xyl-induced ON
144
state (simulations in Figure 3C). The balanced relationship of the two quantities, 200 and 3000,is
145
based on qualitative observations32. Although some rates may have unusual values e.g. low
146
binding, it is important to take into account that these are the output of the optimisation process.
147
Therefore, they are information-rich from a mathematical standpoint since they indicate where
148
the model needs further analysis. During the present study, more realistic rates33,
149
investigated (while maintaining system output) under more complex scenarios than mere time-
150
based approaches.
34
will be
151 152
We observed that Pm activity or its in-silico counterpart θf, is very specific; in other words,
153
changes in certain rates can cause incorrect promoter function. As an example, the stability of the
154
system at variable binding and unbinding rates is shown in Figure 3D. The vertically aligned
155
graphs in the figure shared the value for k-1, while that for k1 was increased to 250% its original
156
value (top graph) or decreased to 40% (bottom). The same ratios were applied to the changes in
157
k-1 in the horizontal simulations, for which k1 was constant (left, right). When these key rates
158
were altered, the ideal behaviour (centre, maximum differential variability of DV =11.7) was no
159
longer maintained.
160 161
A stress analysis on the model tested its robustness to rate variation. GFP degradation rate (k5)
162
was gradually decreased from 2.0 to 0.7 to render the protein more stable; affinity rates (k1 and k-1)
163
were increased (up to 0.64 and 24.0, respectively) to analyse different TF-to-promoter behaviour
164
(Supplementary Text S1). In all cases an adjustment of XylSa numbers for high and low induction
165
(maintaining the 3000-to-200 ratio previously optimised) was sufficient to restore system
166
function. Additionally, we performed a multi-agent simulation (Supplementary Figure S1) to
167
include [i] XylSa entry and degradation rates, [ii] molecule dilution upon cell division and [iii]
168
extrinsic noise (see Methods). A decrease in the unbinding rate was enough to reproduce the
169
noise patterns under this new scenario. That indicates that dilution is not a decisive source of
170
noise. Moreover, both Ps and Pm double their copy number (from one to two) just before
171
division, so potential effects due to dilution are mitigated by a synchronous increase/decrease in
172
both regulators and promoters. Taken together, the above information, plus a more general 6
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
173
simulation of the full TOL network (Supplementary Text S2) suggested that for a given possible
174
set of rates, as θf, the model for Pm expression depends crucially on the number of TF molecules
175
to produce the observed behaviour. We now analyse these considerations in detail.
176 177
System sensitivity to alterations in the number of transcription factor. Simulations of the Pm
178
response to gradual changes in regulator numbers showed our estimated XylSa figures for both
179
ON states (200 and 3000 for 3MBz and m-xyl, respectively) were optimal for maximising
180
differences between the two expression noise regimes (under θf). The simulated Pm transfer
181
function indicated the signal range and mean value at a given XylSa concentration (Figure 4A).
182
Unlike other reported promoter transfer functions35-37, by which transcriptional activity produces
183
similar noise (error bars in graphs) regardless of regulator concentration, here the central section
184
of the curve shows a wider noise profile than the remainder. Indeed, the noise ranges reach
185
maximum and minimum levels at ~200 and 3000 XylSa molecules, respectively. Such wider
186
noise at intermediate induction has been observed38, 39 for gradients of one single input. However,
187
the noise patterns of the Pm promoter are controlled by two distinct inputs, yet one regulator,
188
which makes it dynamically unique. The simulation in Figure 4B shows the system tested in
189
continuous function in which the inducer is changed sequentially. Pm activity suggests a trinary
190
(rather than binary) signal, with three states: one OFF and two ON, each state with a distinct
191
shape that unequivocally recalls its input. There is thus a direct correlation between the time
192
intervals of the bursting effect and the max-min distance (amplitude) of the resulting gene
193
expression levels. Furthermore, such a signal could be of potential use for multi-valued genetic
194
logic circuits beyond the mere binary (0/1) abstraction40, 41.
195 196
Two further analyses that link regulator dynamics with output noise are shown in Figure 4C.
197
They must be interpreted in terms of the pulsing transcriptional bursts that frame the prokaryotic
198
promoter activity. Two measurements, cumulative pulse duration and signal amplitude, were
199
monitored in 24-h simulations of the system at several XylSa numbers (Figure 4C). Cumulative
200
pulse refers to the core of the bursting effect, the total time that Pm is in its active state when
201
XylSa is bound to it. Total XylSa residence time on Pm increased with the number of regulator
202
molecules, which indicated that the ON state produced during m-xyl induction corresponds to 7
ACS Paragon Plus Environment
Page 8 of 32
Page 9 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
203
large numbers of XylSa molecules. Signal amplitude, the difference between the highest and
204
lowest output values of a single simulation run, decreased except in the interval in which the
205
number of XylSa molecules was in the range ∈ [0, ≃150], when the distance between the
206
uppermost/lowermost signal increases. The cognate inflexion point can thus be explained as the
207
number of XylSa molecules that brings about the ON state during 3MBz induction. The fact that
208
such inflection point occurs at low regulator levels (direct consequence of the rate values used)
209
matches our intuition since the inducer 3MBz will activate the only XylS molecules leaked from
210
a silent Ps promoter.
211 212
Influence of intracellular regulator-promoter proximity on transcriptional output. Based
213
on these findings, we asked whether the low-noise regime produced by m-xyl induction would be
214
generated if 3MBz were the only input. When we interrogated the model with this question, the
215
answer was positive when and only when 3MBz co-occurs with large numbers (3000) of the
216
regulator XylSa. Although this condition would appear impossible to achieve since m-xyl, and
217
not 3MBz, is needed to stimulate XylSa production (see Figure 1), we must consider the possible
218
spatial effects of XylSa molecules. In the initial zero-dimensional model, it is assumed that each
219
regulator is able to bind its target promoter at a given fixed rate (thus only time-based), as if it
220
were a pure chemical reaction. When measuring living cells (as in Figure 2), one must consider
221
that due to imperfect diffusion caused by molecular crowding and non-homogenous
222
microviscosity42, 43, not all regulators are equally effective in reaching and binding cognate target
223
DNA sequences. A specific regulator will not interact with its promoter if they cannot meet, since
224
its access to the target promoter will be limited by the ease of diffusion towards the physical Pm
225
location.
226 227
To examine this possibility, we simulated protein trajectories33,
228
Brownian motion43, 46, 47. Figures 5A and B record the trajectories of regulator molecules (XylSa)
229
being expressed from what we term the source region: the physical Ps promoter location from
230
which xylS is expressed. Given the coupling of transcription/translation processes in prokaryotic
231
gene expression48-50, it is safe to assume that the TF protein is produced in close proximity to the
232
Ps-xylS promoter gene pair (Figure 1). To trigger transcription, XylSa must migrate to a 8
ACS Paragon Plus Environment
44, 45
inside a cell, following
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 32
233
physically separate target site where Pm is located; Figure 5A and B illustrate two possible
234
scenarios, which diverge only in the number of proteins stemming from the source region. If
235
regulator numbers are low, there are necessarily empty locations within the cell that XylSa may
236
not encounter easily. Should Pm be located in one such regulator-empty sector, a productive
237
contact is physically impossible, which leads to an OFF promoter state.
238 239
When this situation is scaled up to several thousand bacteria, each individual target region could
240
accommodate a different number of regulators ranging from all to none, leading to pronounced
241
cell-to-cell variability. In contrast, when many regulatory proteins originate in the source region,
242
there are few empty areas in the intracellular space and the variability range narrows. Protein
243
distribution at any given time does not necessarily match trajectory distribution (Figure 5A),
244
meaning that the trails of the regulators are more space-dependent than their spread. As a result,
245
the apparent binding rate, k1, in reality combines promoter-TF affinity proper (the ability of the
246
two molecular partners to interact physically) with the probability that the regulator is located
247
near the promoter (availability). In our initial time-based kinetic model, it was not possible to
248
obtain both expression noises with usual parameter values, for instance, k1=0.5 and k-1=50. Using
249
such parameters in a simulation with 200 TFs (corresponding to 3MBz induction) returned a
250
unimodal distribution (Figure 5B, right) that was far from the noisy pattern observed
251
experimentally (Figure 2C). The initial optimisation process returned a binding value of k1 =
252
0.004, unrealistic yet necessary to reproduce the observations. That was because all TFs were
253
available to bind at any time i.e. the number of XylSa was fixed for a single simulation run.
254
However, by taking a spatial-based approach, where the number of TFs was variable during the
255
stochastic simulation, we restored system function even with k1=0.5 and k-1=50 (Figure 5B, right).
256
To that end, the spatial Brownian-motion simulation updated the Gillespie algorithm with the
257
number of TFs that crossed a specific target region at given time-points (see Methods). This
258
generated a crucial fluctuation in the availability of TFs that could bind Pm. The source-target
259
distance becomes then a decisive parameter.
260 261
Upon the inclusion of spatial dynamics in time-based kinetics, the distinction between the
262
absolute number of TFs and the trajectory points becomes unavoidable. Since a given TF would 9
ACS Paragon Plus Environment
Page 11 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
263
bind/unbind more than once51, it seems coherent to talk about TF crossings in a particular region
264
when referring to the availability parameter described above. That value would be obviously
265
higher than the actual number of TF molecules present in the cell. In our study, we match time-
266
based absolute TF numbers (Figure 3C) with spatial-based trajectory points (Figure 5B).
267 268
Physical proximity between Pm and XylS decreases transcriptional noise.
In the
269
experimental setting that showed the differences in Pm noise (Figure 2), the promoter and its
270
regulator were placed at distant locations within the cell. This was done by placing the Ps
271
promoter for xylS expression and the reporter Pm-gfp fusion in different replicons, the P. putida
272
mt-2-Pm TOL plasmid and the chromosome, respectively (see Methods). A key interpretation of
273
the simulations is that physical proximity between genomic sites bearing the Ps-xylS and Pm-gfp
274
DNA segments would result in better Pm occupation at lower XylS concentrations, and thus in
275
reduced GFP expression noise.
276 277
To test this prediction, we positioned the Ps-xylS and Pm-gfp sequences within the frame of a
278
mini-Tn7 transposon vector (see Methods), which was delivered to the single attTn7 site of P.
279
putida KT2440 (identical to P. putida mt-2 without the TOL plasmid) to generate a strain termed
280
P. putida KT-BGS (Table 1). In these engineered bacteria, the two components of the regulatory
281
device (Ps-xylS/Pm-gfp) were designed to be adjacent, in monocopy and at a fixed chromosomal
282
site, with an artificially minimised distance between TF source and promoter target regions
283
(Figure 6A).
284 285
We then carried out flow cytometry measurements of 3MBz induction in the P. putida KT-BGS
286
strain as for the reference P. putida mt-2-Pm strain (Figure 6A, top). For the sake of comparison,
287
Figure 6A (bottom) reproduces the information for 3MBz-induced P. putida mt-2-Pm (from
288
Figure 2C). The proximity of Ps-xylS to Pm-gfp in P. putida KT-BGS results in a 3MBz
289
response that delivers a much narrower noise regime at high GFP intensity values.
290
fluorescent signals of 3MBz-induced P. putida KT-BGS (in which xylS expression is low but
291
spatially proximal to Pm, Figure 6A) were indistinguishable from those of m-xyl-induced P.
292
putida mt-2-Pm (in which XylS expression is high but distant from the Pm target promoter). 10
ACS Paragon Plus Environment
The
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 32
293
According to simulations, for the P. putida KT-BGS strain, Ps-xylS proximity to Pm-gfp yields a
294
larger number of regulators in the local molecular environment of Pm, which are thus available
295
for binding. The uninduced performance is null in both strains (Supplementary Figure S2)
296
pointing out that XylS molecules alone do not trigger Pm activity. It has been reported that a high
297
concentration of non-active XylS molecules would, in principle, induce Pm52. However, that
298
scenario is not applicable to our study. Not even in KT-BGS where XylS molecules are more
299
abundant in the proximity of Pm.
300 301
To ensure that the modifications in P. putida KT-BGS did not distort the physical structure of the
302
bacteria53, we compared the size and complexity of individual cells to those of the P. putida mt-
303
2-Pm counterpart, in which Ps-xylS and Pm-gfp are separated. The two strains were virtually
304
indistinguishable, with no important differences in physical quality (Supplementary Figure S3).
305
We compared Ps promoter activity in mt-2-Pm and KT-BGS strains by quantitative PCR, which
306
showed similar levels of XylS molecules in response to 3MBz, (Figure 6B); growth curves
307
validated strain comparability (Supplementary Figure S4). Visualisation of the TOL plasmid
308
(Figure 6C, see Methods) allowed us to confirm its single-copy nature, thus comparable to the
309
single copy Ps insertion in the chromosome.
310 311
It has been recently suggested that molecular crowding effects inside bacteria lead to a non-
312
homogenous intracellular space54. This spatial heterogeneity adds complexity to the abstract cell
313
compartment of Figure 5 and would cause an uneven diffusion of molecules. To analyse the
314
simulation effects resulting from local environments with different diffusion specifications (see
315
Methods), we included regions where the Brownian motion was slowed down (Supplementary
316
Figure S5). The presence of low mobility regions facilitated TF accumulation and help regulator
317
numbers increase in highly-condensed areas. This scenario matched the physical architecture of P.
318
putida KT-BGS, in which the target region was inserted into the chromosome, presumably a
319
crowded environment with restricted mobility.
320 321
The noise-dependence of promoter-to-regulator distance is likely to increase if the TF is very
322
unstable, as appears to be the case for XylS55. These data identify a function within the intricate 11
ACS Paragon Plus Environment
Page 13 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
323
architecture of the regulatory network that governs biodegradation of m-xylene in Pseudomonas.
324
Future work on this line would focus on re-arranging network components across different
325
genome locations and/or vary their copy number to establish universal spatial-dependent design
326
rules. Moreover, it is worth studying in detail the complex motility of molecules inside crowded
327
regions of the cell56 in order to develop more accurate spatial predictions.
328 329
Discussion
330 331
Intracellular signals are transmitted according to specific dynamics of the components involved
332
in their transfer. These communications are therefore endowed with precise information, whose
333
decodification promises valuable insights into cellular kinetic and structural properties9. Signal
334
variability, commonly referred to as gene expression noise4, 5, constitutes the fingerprint of such
335
transmission, and thus the target data to be interpreted. In the case documented here, the
336
expression signals produced initially by the Pm promoter activity in P. putida19, 22, 23 lead to
337
highly specific, stable noise patterns depending on the stimulus to which the cells were exposed.
338 339
Using mathematical modelling and computational analysis, we deconvoluted the flow cytometry
340
data for each setting to describe the kinetics that could reproduce that behaviour. As a result, the
341
kinetic values that fit the experimental observations highlight the importance of the bursting-
342
specific rates, binding and unbinding13-16, where each of these values influences the final
343
promoter activity distinctly. We pinpointed how the dynamics of Pm-regulator interplay
344
determines gene expression by including spatial effects, in particular protein distribution, within a
345
cell. Our model, validated by the experiments shown above, indicates that the physical distance
346
between the regulator source and the target promoter is translated into specific noise patterns that
347
change radically depending on promoter-TF proximity. This is due to the fact that regulators, or
348
rather their trajectories (Figure 5), are not distributed homogeneously57 and TF are thus more
349
likely to meet the promoters they regulate if located near the source58. This concept was
350
hypothesized by Ten Wolde to explain the frequent genomic association of TF and target
351
promoters as an evolutionary remedy to an excess of noise59, 60. In contrast, our analyses raise
352
questions as to whether gene expression noise caused by a non-homogeneous intracellular matrix 12
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 32
353
is an adaptive trait that endows regulatory networks with specific properties. We show that
354
changing the spatial positioning of components, Pm noise patterns can be altered, which opens
355
the opportunity to use one expression profile or another depending on needs. As one Pm activity
356
regime is much more variable than the other, it might well have been co-opted evolutionarily to
357
create phenotypic heterogeneity within a population to increase its metabolic or else fitness61, 62.
358
The other regime, full-expression and low variability could be used, for example for decreasing
359
phenotypic diversity of a clonal population of productive cells63.
360 361
Our data also offer a new challenge for engineering non-native regulatory circuits or, more
362
generally, synthetic genetic implants in the genomic and biochemical chassis of a bacterial cell36,
363
41, 64
364
in the spatial frame of a cell for optimal performance, a question that is rarely considered in
365
contemporary synthetic biology and which deserves more attention. Finally, our results add a new
366
perspective to the much-debated generality of transcription/translation coupling in prokaryotes48-
367
50
368
generated near the genes that encode them or whether they must migrate to other cell sites.
. Each gene sequence and each protein (including TF) might need a specific physical address
, as the noise regime of promoters is certainly influenced by whether their cognate TF are
369 370
Materials and Methods
371 372
Bacterial strains, growth conditions and genetic constructs. Bacterial strains and plasmids
373
used are listed in Table 1. Escherichia coli cells were grown at 37ºC in LB medium and used as
374
hosts for cloning procedures. Pseudomonas putida cells were incubated at 30ºC in M9 minimal
375
medium supplemented with 2 mM MgSO4 and 20 mM citrate as sole carbon source65. When
376
needed, gentamycin (Gm; 10 µg mL-1), kanamycin (Km; 50 µg mL-1), ampicillin (Ap; 150 µg
377
mL-1) and chloramphenicol (Cm; 30 µg mL-1) were added to growth media. Reporter strain P.
378
putida mt-2-Pm is the original TOL plasmid pWW0-containing P. putida mt-2, which was
379
inserted in the single attTn7 site of its genome with a Pm-gfp transcriptional fusion. The Pm
380
promoter sequence was amplified from plasmid pSEVA22866 as a 122 base pair (bp) PacI/AvrII
381
fragment with primers 5’TTAATTAAGGTTTGATAGGGATAAGTCC3’ and 5’CCTAGGT
382
CTGTTGCATAAAGCCTAA3’, and cloned into the mini-Tn7 promoter-calibrating vector pBG 13
ACS Paragon Plus Environment
Page 15 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
383
(Zoebel et al., 2015). The organization of this vector (Supplementary Figure S6A) is such that
384
inserting promoter-bearing PacI/AvrII originates a standardised translation/transcription fusion
385
that minimises any effect of the non-translated 5’ end of the reporter transcript in the final GFP
386
readout. Cloning the Pm promoter in pBG generated mini-Tn7 delivery vector pBG-Pm
387
(Supplementary Figure S6B). This construct was then mobilised to pWW0-containing P. putida
388
mt-2 strain by tetra-parental mating67. Finally, GmR exconjugants were verified for insertion of
389
the hybrid mini-Tn7 transposon (bearing the Pm-gfp fusion) in a specific orientation at the attTn7
390
site by amplifying the genomic region of interest with diagnostic PCR using primer pairs 5-Pput-
391
glmS
392
ACACCC3’
393
5’CACAGCATAACTGGACTGATTTC3’. One of these clones yielding DNA products of 400
394
and 200 bp68, 69, was designated as P. putida mt-2-Pm and used for the experiments discussed
395
above. To obtain an entirely equivalent P. putida strain with a physically rearranged XylS/Pm
396
regulatory node, a 1088 bp DNA segment containing the array of regulatory parts xylS ← Ps -
397
Pm → was excised from plasmid pSEVA22866 as a PacI/AvrII fragment and cloned in the
398
corresponding sites of the pBG vector (Supplementary Figure S6C). The resulting construct
399
(pBGS) was mobilised to the genome of the pWW0-less strain P. putida KT2440, and GmR
400
exconjugants were tested for insertion of the mini-Tn7 transposon (with the Pm-gfp fusion
401
adjacent to the xylS gene) in the same genomic site and orientation as before. One of these clones,
402
termed P. putida KT-BGS, was chosen to test the effects of XylS/Pm proximity. This genetic
403
strategy allowed a faithful comparison between the expression noise produced by the Pm-gfp
404
fusion borne by either P. putida mt-2-Pm (Ps → xylS and Pm in non-adjacent, separate replicons)
405
or P. putida KT-BGS (Ps → xylS and Pm in close genomic proximity).
UP
5’AGTCAGAGTTACGGAATTGTAGG3’/3-Tn7L and
5-PpuglmS
DOWN
(5’ATTAGCTTACGACGCT
5’TTACGTGGCCGTGCTAAAGGG3’/3-Tn7R
406 407
Single cell analysis by flow cytometry. Single-cell experiments were performed with a Gallios
408
(Beckam Coulter) flow cytometer. GFP was excited at 488 nm, and the fluorescence signal
409
recovered with a 525(40) BP filter. Strains grown overnight were diluted 1/100 and allowed to
410
grow at 30ºC in pre-filtered M9 citrate medium and incubated (3-4 h). After pre-incubation, cells
411
in the late exponential phase (OD600nm = 0.4) were treated with the inducer 3MBz (1 mM);
412
cultures were incubated with aeration (30ºC). At 3 h post-induction, an aliquot of each sample 14
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
413
Page 16 of 32
was analysed by flow cytometry; 20,000 events were analysed for each sample.
414 415
RNA purification and real-time q-PCR. P. putida strains KT-BGS and mt-2-Pm were grown
416
overnight (30ºC) in citrate-supplemented M9 with in aerated flasks. Cultures were diluted 1/100
417
in the same medium, alone or with 1 mM 3MBz and grown to exponential phase (OD600 = 0.3-
418
0.5). A 3 ml aliquot of each sample was treated with 300 µl ice-top solution (5% water-saturated
419
phenol in ethanol) and centrifuged; pellets were frozen in liquid nitrogen and stored at -80ºC.
420
RNA was extracted using the miRNeasy kit (Qiagen) with some modifications to optimise
421
isolation of high quality RNA from P. putida. The quantity of total RNA was determined in a
422
Nanovue Plus spectrophotometer (GE Healthcare Life Sciences), and RNA integrity analysed by
423
agarose gel electrophoresis. The absence of DNA was confirmed using primers for rpoN
424
(5’TCGACCCGGAGCTGGATA
425
(5’CGAGTTGCTGGAGATTGTGT and 3’TCGTTAAATTGCCCTCAGTG). Real-time RT-
426
PCR was performed using total RNA preparations from three independent cultures (three
427
biological replicates). The process was monitored by q-PCR in an ABI pRISM 7900HT Fast
428
RealTime PCR system (Applied Biosystems). To calculate the relative amount of xylS transcript,
429
we used the ∆∆Ct method70, which is designed to compare levels of a given RNA in two
430
conditions (induced and uninduced cells). The primers used for xylS were (xylS 141: 5’TAAT
431
CCAGGCGAGATTACCC and 3’AACCAGTATGTCGGTACGCA; xylS 108: 5’CGAGTT
432
GCTGGAGATTGTGT and 3’TCGTTAAATTGCCCTCAGTG; rpoN: 5’TCGACCCGGAG
433
CTGGATA and 3’CGGCTCGAACTGCTGGAT). Results were normalised relative to for the
434
rpoN gene, the expression of which remains constant throughout the growth curve.
and
3’CGGCTCGAACTGCTGGAT)
and
xylS
435 436
Kinetic reactions, simulation and optimisation. The kinetic reactions that describe the model
437
depicted in Figure 3 are:
15
ACS Paragon Plus Environment
Page 17 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
438 439
where Pma and Pm are the promoter with and without bound XylSa, respectively, XylSa denotes
440
the regulator in its active form, mRNA is the output of the transcription process and GFP is the
441
final green fluorescent protein. The description of the rates is as follows: k1 is the binding rate of
442
XylSa to Pm (molecules-1 hour-1), k-1, the unbinding rate of XylSa from the promoter (hour-1), k2
443
and k3, the transcription and translation rates (hour-1), k4 and k5 are the degradation rates of mRNA
444
and GFP (hour-1) and k6, the basal transcription of the promoter, which is Pm activity with no
445
regulator bound (hour-1). Stochastic simulations were performed using the Gillespie algorithm71.
446
To obtain cytometry-like graphs (Figure 3C), the values Gillespie’s algorithm returns must be
447
converted into a time-course array in which time intervals are fixed, τa, and are small enough to
448
have cells (each time point) that correctly represent (in terms of frequency) all possible molecular
449
levels72. Here we used τa = 0.01 h (Supplementary Figure S7). The values in Figure 3C where
450
optimised within the ranges k1 ∈[0.001 - 1.2], k-1 ∈[0.2 - 80], k2 ∈[100 - 1000], k3 ∈[10 - 120]
451
and XylSa ∈ [20 – 3000], were each combination was weighted following the fitness parameters
452
width and flatness (Figure 3B). Each combination of rates was assigned a fitness value. The
453
width of the output distribution was measured as the difference between maximum and minim
454
protein values. The flatness parameter corresponds to the minimum probability that represents
455
any given protein value of a given run. Such value would be 0 in a one-peak Gaussian-like
456
distribution while higher for a flat surface. Therefore, the goal (for 3MBz induction) was to
457
maximize width while at the same time maximizing flatness. Five runs per combination of rates
458
assisted the process to discard outlier values (not representative of the final distribution).
459
Different degradation rates were tested, as well as affinity values, in the stress analysis of
460
Supplementary Text S1 to test system robustness. The full TOL network was formalized,
461
including reactions and rates, and simulated in Supplementary Text S2.
462 463
The multi-agent simulation of Supplementary Figure S2 was performed using our software
464
DiSCUS72, 73. XylSa entry rate (k7 = 400 molecules hour-1) and degradation (k8 = 2 hour-1) were
465
added to reactions 1-6. Extrinsic noise was simulated by changing the rates vector after division
466
to reflect fluctuations in environmental conditions. Every new rate was the result of a Gaussian
467
distribution where the mean was the previous rate value and the standard deviation a 20% of the 16
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 32
468
original rate. Dilution was simulated by dividing by two the number of molecules in each
469
daughter cell compared with the mother.
470 471
Differential variability (DV), the relation between the variance of the noise under two different
472
conditions74, was used to measure the simulations of Figure 3D. This value was defined as f =
473
σ12/ σ22 where σ12 and σ22 are the variances of the signal at low induction (3MBz) and high
474
induction (m-xyl), respectively.
475 476
Spatial protein trajectories. For the spatial simulation of regulators shown in Figures 5 we
477
implemented a two-dimensional Brownian motion instance, written as an iteration scheme as
478
follows:
479
480 481
where t identifies the last time event, dt is the time step (dt = T⁄N with T the total time per
482
iteration and N the number of steps, 15.0 and 1.0, respectively, in this case) and δ the so-called
483
Wiener process parameter (here, 0.25). Each protein ran for 400 iterations. Time parameters were
484
dimensionless and the simulated cell area was based on a 60 x 20 2D lattice. The target region
485
monitored to obtain the graphs of Figure 5B (right) was a 10 x 5 lattice situated on a pole (source
486
region at the middle). The numbers of trajectory points (not TF numbers, but TF crossings) in the
487
target region for the ‘variable TF’ simulation were obtained during 16 spatial simulations: 147,
488
87, 47, 128, 103, 71, 12, 26, 91, 35, 47, 0, 214, 21, 29, 152. These numbers updated the TF
489
variable within a single Gillespie simulation to generate the wide-rage plateau-like signal. To
490
simulate the low mobility areas of the heterogeneous spatial distribution (Figure S5) the time step
491
was set to 1.0 in such areas in contrast with the 15.0 of the normal mobility regions.
492 493
All computational simulations were written in Python. The source code is available in
494
Supplementary File S1.
495
17
ACS Paragon Plus Environment
Page 19 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Synthetic Biology
496
Visualization of the TOL plasmid (Figure 6C). The plasmid was labelled using the fluorescent
497
operator repressor system tetO-TetR-EYFP75 to investigate its localization within the cell
498
P.putida mt-2. Tandemly repeated tetO sequences were introduced into specific locus of the TOL
499
plasmid. Given TetR fused EYFP chimera, the labelled DNA was visualized using epifluorescent
500
microscopy. Cells were grown on succinate amended agarose pad.
501 502
ASSOCIATED CONTENT
503 504
This material is available free of charge via the Internet at http://pubs.acs.org:
505
Supplementary Text S1. Stress analysis document. Effects of rate variation on circuit
506
performance.
507
Supplementary Figure S1. System variances due to the inclusion of dilution dynamics and
508
extrinsic noise in the model.
509
Supplementary Figure S2. Side-by-side comparison of Pm activity in strains mt-2-Pm and KT-
510
BGS induced with 3MBz.
511
Supplementary Figure S3. Comparison of cell complexity between mt-2-Pm and KT-BGS
512
strains.
513
Supplementary Figure S4. Growth curves for P. putida KT-BGS and mt-2-Pm strains.
514
Supplementary Figure S5. Local TF density affected by heterogeneous diffusion areas within
515
the cell.
516
Supplementary Figure S6. Organisation of genetic constructs for inspection of transcriptional
517
noise of the Pm promoter using GFP cytometry.
518
Supplementary Figure S7. Time-course conversion from the Gillespie output to a frequency-
519
constant line.
520
Supplementary Text S2. Full TOL simulation
521
Supplementary File S1. Compressed file containing all Python source code used in this study.
522 523
COMPETING INTERESTS 18
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
524
Page 20 of 32
The authors declare that there are no competing interests.
525 526
AUTHORS' CONTRIBUTIONS
527
AGM and VDL conceived the whole study and wrote the article. IB and JK carried out the
528
experimental parts of the work. All the authors contributed to the discussion of the research and
529
interpretation of the data.
530 531
ACKNOWLEDGEMENTS
532
This work was funded by the CAMBIOS Project of the Spanish Ministry of Economy and
533
Competitiveness RTC-2014-1777-3 (MINECO), HELIOS Project of the Spanish Ministry of
534
Economy and Competitiveness BIO 2015-66960-C3-2-R (MINECO/FEDER). ARISYS (ERC-
535
2012-ADG-322797), EmPowerPutida (EU-H2020-BIOTEC-2014-2015-6335536) and FUTURE
536
(704410-H2020-MSCA-IF-15) Contracts of the European Union.
537 538
References
539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558
1. Kærn, M., Elston, T. C., Blake, W. J., and Collins, J. J. (2005) Stochasticity in gene expression: from theories to phenotypes, Nature Reviews Genetics 6, 451-464. 2. Golding, I., Paulsson, J., Zawilski, S. M., and Cox, E. C. (2005) Real-time kinetics of gene activity in individual bacteria, Cell 123, 1025-1036. 3. Raj, A., and van Oudenaarden, A. (2008) Nature, nurture, or chance: stochastic gene expression and its consequences, Cell 135, 216-226. 4. Eldar, A., and Elowitz, M. B. (2010) Functional roles for noise in genetic circuits, Nature 467, 167-173. 5. Rinott, R., Jaimovich, A., and Friedman, N. (2011) Exploring transcription regulation through cell-to-cell variability, Proceedings of the National Academy of Sciences 108, 6329-6334. 6. McAdams, H. H., and Arkin, A. (1999) It’sa noisy business! Genetic regulation at the nanomolar scale, Trends in genetics 15, 65-69. 7. Hansen, A. S., and O'Shea, E. K. (2013) Promoter decoding of transcription factor dynamics involves a trade‐off between noise and control of gene expression, Molecular systems biology 9, 704. 8. Munsky, B., Neuert, G., and van Oudenaarden, A. (2012) Using gene expression noise to understand gene regulation, Science 336, 183-187. 9. Purvis, J. E., and Lahav, G. (2013) Encoding and decoding cellular information through signaling dynamics, Cell 152, 945-956. 19
ACS Paragon Plus Environment
Page 21 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603
ACS Synthetic Biology
10. Brehm-Stecher, B. F., and Johnson, E. A. (2004) Single-cell microbiology: tools, technologies, and applications, Microbiology and molecular biology reviews 68, 538-559. 11. Czechowska, K., Johnson, D. R., and van der Meer, J. R. (2008) Use of flow cytometric methods for single-cell analysis in environmental microbiology, Current opinion in microbiology 11, 205-212. 12. Kortmann, H., Blank, L. M., and Schmid, A. (2010) Single cell analytics: An overview, In High Resolution Microbial Single Cell Analytics, pp 99-122, Springer. 13. Taniguchi, Y., Choi, P. J., Li, G.-W., Chen, H., Babu, M., Hearn, J., Emili, A., and Xie, X. S. (2010) Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells, Science 329, 533-538. 14. Zong, C., So, L. h., Sepúlveda, L. A., Skinner, S. O., and Golding, I. (2010) Lysogen stability is determined by the frequency of activity bursts from the fate‐determining gene, Molecular systems biology 6, 440. 15. Chong, S., Chen, C., Ge, H., and Xie, X. S. (2014) Mechanism of transcriptional bursting in bacteria, Cell 158, 314-326. 16. So, L.-h., Ghosh, A., Zong, C., Sepúlveda, L. A., Segev, R., and Golding, I. (2011) General properties of transcriptional time series in Escherichia coli, Nature genetics 43, 554-560. 17. Nikel, P. I., Silva‐Rocha, R., Benedetti, I., and Lorenzo, V. (2014) The private life of environmental bacteria: pollutant biodegradation at the single cell level, Environmental microbiology 16, 628-642. 18. Nikel, P. I., Martínez-García, E., and de Lorenzo, V. (2014) Biotechnological domestication of pseudomonads using synthetic biology, Nature Reviews Microbiology 12, 368-379. 19. De Las Heras, A., Fraile, S., and de Lorenzo, V. (2012) Increasing signal specificity of the TOL network of Pseudomonas putida mt-2 by rewiring the connectivity of the master regulator XylR, PLoS Genet 8, e1002963. 20. Nikel, P. I., Romero-Campero, F. J., Zeidman, J. A., Goñi-Moreno, Á., and de Lorenzo, V. (2015) The glycerol-dependent metabolic persistence of Pseudomonas putida KT2440 reflects the regulatory logic of the GlpR repressor, mBio 6, e00340-00315. 21. Silva-Rocha, R., Pérez-Pantoja, D., and de Lorenzo, V. (2013) Decoding the genetic networks of environmental bacteria: regulatory moonlighting of the TOL system of Pseudomonas putida mt-2, The ISME journal 7, 229. 22. Ramos, J. L., Marqués, S., and Timmis, K. N. (1997) Transcriptional control of the Pseudomonas TOL plasmid catabolic operons is achieved through an interplay of host factors and plasmid-encoded regulators, Annual Reviews in Microbiology 51, 341-373. 23. González-Pérez, M., Marqués, S., Domı́nguez-Cuevas, P., and Ramos, J. L. (2002) XylS activator and RNA polymerase binding sites at the Pm promoter overlap, FEBS letters 519, 117-122. 24. Pérez‐Pantoja, D., Kim, J., Silva‐Rocha, R., and Lorenzo, V. (2015) The differential response of the Pben promoter of Pseudomonas putida mt‐2 to BenR and XylS prevents metabolic conflicts in m‐xylene biodegradation, Environmental microbiology 17, 64-75. 25. González-Pérez, M. M., Ramos, J. L., and Marqués, S. (2004) Cellular XylS levels are a function of transcription of xylS from two independent promoters and the differential efficiency of translation of the two mRNAs, Journal of bacteriology 186, 1898-1901. 26. Goñi-Moreno, A. (2014) On genetic logic circuits: forcing digital electronics standards?, Memetic Computing 6, 149-155. 20
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648
Page 22 of 32
27. Dublanche, Y., Michalodimitrakis, K., Kümmerer, N., Foglierini, M., and Serrano, L. (2006) Noise in transcription negative feedback loops: simulation and experimental analysis, Molecular systems biology 2, 41. 28. Goñi-Moreno, A., and Amos, M. (2012) Continuous computation in engineered gene circuits, Biosystems 109, 52-56. 29. Balagaddé, F. K., Song, H., Ozaki, J., Collins, C. H., Barnet, M., Arnold, F. H., Quake, S. R., and You, L. (2008) A synthetic Escherichia coli predator–prey ecosystem, Molecular systems biology 4, 187. 30. Andersen, J. B., Sternberg, C., Poulsen, L. K., Bjørn, S. P., Givskov, M., and Molin, S. (1998) New unstable variants of green fluorescent protein for studies of transient gene expression in bacteria, Applied and environmental microbiology 64, 2240-2246. 31. de-Leon, S. B.-T., and Davidson, E. H. (2009) Modeling the dynamics of transcriptional gene regulatory networks for animal development, Developmental biology 325, 317-328. 32. Velázquez, F., Parro, V., and de Lorenzo, V. (2005) Inferring the genetic network of m‐xylene metabolism through expression profiling of the xyl genes of Pseudomonas putida mt‐2, Molecular microbiology 57, 1557-1569. 33. Elf, J., Li, G.-W., and Xie, X. S. (2007) Probing transcription factor dynamics at the singlemolecule level in a living cell, Science 316, 1191-1194. 34. Miró-Bueno, J. M., and Rodríguez-Patón, A. (2011) A simple negative interaction in the positive transcriptional feedback of a single gene is sufficient to produce reliable oscillations, PloS one 6, e27414. 35. Moon, T. S., Lou, C., Tamsir, A., Stanton, B. C., and Voigt, C. A. (2012) Genetic programs constructed from layered logic gates in single cells, Nature 491, 249-253. 36. Wang, B., Kitney, R. I., Joly, N., and Buck, M. (2011) Engineering modular and orthogonal genetic logic gates for robust digital-like synthetic biology, Nature communications 2, 508. 37. Bonnet, J., Yin, P., Ortiz, M. E., Subsoontorn, P., and Endy, D. (2013) Amplifying genetic logic gates, Science 340, 599-603. 38. Murphy, K. F., Balázsi, G., and Collins, J. J. (2007) Combinatorial promoter design for engineering noisy gene expression, Proceedings of the National Academy of Sciences 104, 12726-12731. 39. Murphy, K. F., Adams, R. M., Wang, X., Balazsi, G., and Collins, J. J. (2010) Tuning and controlling gene expression noise in synthetic gene networks, Nucleic acids research, gkq091. 40. Goñi-Moreno, A., and Amos, M. (2012) A reconfigurable NAND/NOR genetic logic gate, BMC systems biology 6, 126. 41. Daniel, R., Rubens, J. R., Sarpeshkar, R., and Lu, T. K. (2013) Synthetic analog computation in living cells, Nature 497, 619-623. 42. Miklos, A. C., Sarkar, M., Wang, Y., and Pielak, G. J. (2011) Protein crowding tunes protein stability, Journal of the American Chemical Society 133, 7116-7120. 43. Parry, B. R., Surovtsev, I. V., Cabeen, M. T., O’Hern, C. S., Dufresne, E. R., and JacobsWagner, C. (2014) The bacterial cytoplasm has glass-like properties and is fluidized by metabolic activity, Cell 156, 183-194. 44. Gahlmann, A., and Moerner, W. E. (2014) Exploring bacterial cell biology with singlemolecule tracking and super-resolution imaging, Nature Reviews Microbiology 12, 9-22. 21
ACS Paragon Plus Environment
Page 23 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693
ACS Synthetic Biology
45. Uphoff, S., and Kapanidis, A. N. (2014) Studying the organization of DNA repair by singlecell and single-molecule imaging, DNA repair 20, 32-40. 46. Uhlenbeck, G. E., and Ornstein, L. S. (1930) On the theory of the Brownian motion, Physical review 36, 823. 47. Saffman, P. G., and Delbrück, M. (1975) Brownian motion in biological membranes, Proceedings of the National Academy of Sciences 72, 3111-3113. 48. Gowrishankar, J., and Harinarayanan, R. (2004) Why is transcription coupled to translation in bacteria?, Molecular microbiology 54, 598-603. 49. Burmann, B. M., Schweimer, K., Luo, X., Wahl, M. C., Stitt, B. L., Gottesman, M. E., and Rösch, P. (2010) A NusE: NusG complex links transcription and translation, Science 328, 501-504. 50. Miller, O. L., Hamkalo, B. A., and Thomas, C. A. (1970) Visualization of bacterial genes in action, Science 169, 392-395. 51. Balleza, E., Lopez-Bojorquez, L. N., Martínez-Antonio, A., Resendis-Antonio, O., LozadaChávez, I., Balderas-Martínez, Y. I., Encarnación, S., and Collado-Vides, J. (2009) Regulation by transcription factors in bacteria: beyond description, FEMS microbiology reviews 33, 133-151. 52. Silva-Rocha, R., De Jong, H., Tamames, J., and De Lorenzo, V. (2011) The logic layout of the TOL network of Pseudomonas putida pWW0 plasmid stems from a metabolic amplifier motif (MAM) that optimizes biodegradation of m-xylene, BMC systems biology 5, 191. 53. Weng, X., and Xiao, J. (2014) Spatial organization of transcription in bacterial cells, Trends in Genetics 30, 287-297. 54. Yu, I., Mori, T., Ando, T., Harada, R., Jung, J., Sugita, Y., and Feig, M. (2016) Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm, eLife 5, e19274. 55. González-Pérez, M. M., Ramos, J. L., Gallegos, M. a.-T., and Marqués, S. (1999) Critical nucleotides in the upstream region of the XylS-dependent TOL meta-cleavage pathway operon promoter as deduced from analysis of mutants, Journal of Biological Chemistry 274, 2286-2290. 56. Halford, S. E. (2009) An end to 40 years of mistakes in DNA–protein association kinetics?, Portland Press Limited. 57. Ishihama, A., Kori, A., Koshio, E., Yamada, K., Maeda, H., Shimada, T., Makinoshima, H., Iwata, A., and Fujita, N. (2014) Intracellular concentrations of transcription factors in Escherichia coli: 65 species with known regulatory functions, Journal of Bacteriology, JB. 01579-01514. 58. Dröge, P., and Müller‐Hill, B. (2001) High local protein concentrations at promoters: strategies in prokaryotic and eukaryotic cells, Bioessays 23, 179-183. 59. Warren, P. B., and Ten Wolde, P. R. (2004) Statistical analysis of the spatial distribution of operons in the transcriptional regulation network of Escherichia coli, Journal of molecular biology 342, 1379-1390. 60. van Zon, J. S., Morelli, M. J., Tǎnase-Nicola, S., and ten Wolde, P. R. (2006) Diffusion of transcription factors can drastically enhance the noise in gene expression, Biophysical journal 91, 4350-4367. 61. de Jong, I. G., Haccou, P., and Kuipers, O. P. (2011) Bet hedging or not? A guide to proper 22
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738
Page 24 of 32
classification of microbial survival strategies, Bioessays 33, 215-223. 62. Veening, J.-W., Smits, W. K., and Kuipers, O. P. (2008) Bistability, epigenetics, and bethedging in bacteria, Annu. Rev. Microbiol. 62, 193-210. 63. Delvigne, F., and Goffin, P. (2014) Microbial heterogeneity affects bioprocess robustness: Dynamic single‐cell analysis contributes to understanding of microbial populations, Biotechnology journal 9, 61-72. 64. Lou, C., Liu, X., Ni, M., Huang, Y., Huang, Q., Huang, L., Jiang, L., Lu, D., Wang, M., and Liu, C. (2010) Synthesizing a novel genetic sequential logic circuit: a push‐on push‐off switch, Molecular systems biology 6, 350. 65. Abril, M. A., Michan, C., Timmis, K. N., and Ramos, J. L. (1989) Regulator and enzyme specificities of the TOL plasmid-encoded upper pathway for degradation of aromatic hydrocarbons and expansion of the substrate range of the pathway, Journal of Bacteriology 171, 6782-6790. 66. Martínez-García, E., Aparicio, T., Goñi-Moreno, A., Fraile, S., and de Lorenzo, V. (2014) SEVA 2.0: an update of the Standard European Vector Architecture for de-/reconstruction of bacterial functionalities, Nucleic acids research, gku1114. 67. Keen, N. T., Tamaki, S., Kobayashi, D., and Trollinger, D. (1988) Improved broad-host-range plasmids for DNA cloning in gram-negative bacteria, Gene 70, 191-197. 68. Schweizer, H. P. (2001) Vectors to express foreign genes and techniques to monitor gene expression in Pseudomonads, current opinion in biotechnology 12, 439-445. 69. Bao, Y., Lies, D. P., Fu, H., and Roberts, G. P. (1991) An improved Tn7-based system for the single-copy insertion of cloned genes into chromosomes of gram-negative bacteria, Gene 109, 167-168. 70. Schmittgen, T. D., and Livak, K. J. (2008) Analyzing real-time PCR data by the comparative CT method, Nature protocols 3, 1101-1108. 71. Gillespie, D. T. (1977) Exact stochastic simulation of coupled chemical reactions, The journal of physical chemistry 81, 2340-2361. 72. Goñi-Moreno, A., Carcajona, M., Kim, J., Martinez-García, E., Amos, M., and de Lorenzo, V. (2016) An implementation-focused bio/algorithmic workflow for synthetic biology, ACS Synthetic Biology. 73. Goni-Moreno, A., and Amos, M. (2015) DiSCUS: A simulation platform for conjugation computing, In International Conference on Unconventional Computation and Natural Computation, pp 181-191, Springer. 74. Ho, J. W. K., Stefani, M., dos Remedios, C. G., and Charleston, M. A. (2008) Differential variability analysis of gene expression and its application to human diseases, Bioinformatics 24, i390-i398. 75. Vallet-Gely, I., and Boccard, F. (2013) Chromosomal organization and segregation in Pseudomonas aeruginosa, PLoS Genet 9, e1003492. 76. Herrero, M., de Lorenzo, V., and Timmis, K. N. (1990) Transposon vectors containing nonantibiotic resistance selection markers for cloning and stable chromosomal insertion of foreign genes in gram-negative bacteria, Journal of bacteriology 172, 6557-6567. 77. Boyer, H. W., and Roulland-dussoix, D. (1969) A complementation analysis of the restriction and modification of DNA in Escherichia coli, Journal of molecular biology 41, 459-472. 78. Nelson, K. E., Weinel, C., Paulsen, I. T., Dodson, R. J., Hilbert, H., Martins dos Santos, V. A. P., Fouts, D. E., Gill, S. R., Pop, M., and Holmes, M. (2002) Complete genome sequence 23
ACS Paragon Plus Environment
Page 25 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757
ACS Synthetic Biology
and comparative analysis of the metabolically versatile Pseudomonas putida KT2440, Environmental microbiology 4, 799-808. 79. Bagdasarian, M., Lurz, R., Rückert, B., Franklin, F. C. H., Bagdasarian, M. M., Frey, J., and Timmis, K. N. (1981) Specific-purpose plasmid cloning vectors II. Broad host range, high copy number, RSF 1010-derived vectors, and a host-vector system for gene cloning in Pseudomonas, Gene 16, 237-247. 80. Worsey, M. J., and Williams, P. A. (1975) Metabolism of toluene and xylenes by Pseudomonas (putida (arvilla) mt-2: evidence for a new function of the TOL plasmid, Journal of Bacteriology 124, 7-13. 81. Kessler, B., Herrero, M., Timmis, K. N., and De Lorenzo, V. (1994) Genetic evidence that the XylS regulator of the Pseudomonas TOL meta operon controls the Pm promoter through weak DNA-protein interactions, Journal of bacteriology 176, 3171-3176. 82. Choi, K.-H., Gaynor, J. B., White, K. G., Lopez, C., Bosio, C. M., Karkhoff-Schweizer, R. R., and Schweizer, H. P. (2005) A Tn7-based broad-range bacterial cloning and expression system, Nature methods 2, 443-448. 83. Zobel, S., Benedetti, I., Eisenbach, L., de Lorenzo, V., Wierckx, N., and Blank, L. M. (2015) Tn7-Based device for calibrated heterologous gene expression in Pseudomonas putida, ACS synthetic biology 4, 1341-1351.
758
24
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
759
Page 26 of 32
Figures
760 761 762 763 764
765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784
Figure 1. Structure of the TOL network borne by plasmid pWW0of Pseudomonas putida mt-2 with Pm promoter as output.
A. As shown, m-xylene is first converted to 3-methylbenzoate (3MBz) through the action of the enzymes encoded by the upper TOL pathway; this compound is further processed by the activity of the lower, and the resulting intermediates are metabolised in the tricarboxylic acid (TCA) cycle. XylR and XylS are transcriptional regulators, whereas Pu, Pm, Ps and Pr are promoters. The master regulatory gene xylR controls expression of both the upper pathway and the second transcription factor, XylS, which is encoded in a location adjacent to the end of the lower operon.. This regulatory architecture has a decisive role in Pm activation dynamics, as the levels of its cognate activator (XylS) vary depending on the inducer. In one case, 3MBz activates those XylS molecules found in low numbers in the cell due to leaky Ps promoter expression. This results in the active form of the protein, termed XylSa, which is able to bind and activate Pm. In the second case, m-xylene (m-xyl) causes both XylS overexpression (due to Ps activation by XylR) and intracellular production of metabolic 3MBz (due to Pu-driven activity of the upper pathway operon)... B. Abstracted Pm activation as an OR logic gate. Pm activity will be triggered by using either 3MBz or m-xyl as the inducer. C. m-xyl leads to a higher concentration of XylSa (XylS overexpression and intracellular generation of 3MBz) than externally added 3MBz (use of only leaked XylS moecules). This difference is the key feature for decoding Pm output
25
ACS Paragon Plus Environment
Page 27 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
785 786 787
788 789 790 791 792 793 794 795 796 797 798 799 800 801 802
ACS Synthetic Biology
Figure 2. Variable noise patterns depending on input signal in P. putida mt-2-Pm strain
A. In these experiments, XylS molecules are produced by the TOL plasmid borne by the P. putida mt-2-Pm, whereas the target Pm-gfp reporter fusion is inserted in the chromosome (see Methods); that is, the source of the TF and its target promoter are non-adjacent and encoded in separate mono-copy replicons (i.e., TOL plasmid and chromosome). B. In cells alone or in the presence of m-xyl (m-xylene), the Pm promoter activity recorded (based on green fluorescent protein intensity) can be abstracted as a binary switch with a 0 or OFF state and a 1 or ON state. Flow cytometry results show this behaviour, where the noise range allows a null overlap between 0 and 1 (termed 1a to differentiate it from the following). C. Using 3MBz as the inducer again provokes switch-like behaviour in Pm, with 0/OFF and 1/ON states. As the cytometry results show, the noise range is much wider here, from maximum expression to the minimum (ON state thus called 1b).
26
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
803 804
805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823
Page 28 of 32
Figure 3. Rate optimisation according to output-state fitness and binding sensitivity
A. The Pm promoter studied here and the rates involved in the model. XylSa in its active form is the activator of the inducible promoter. The Pma complex refers to the promoter with the bound regulator. Rates k1 and k-1 correspond to binding and unbinding events, respectively. Transcription, k2, translation, k3, degradation rates, k4 and k5; the basal transcription rate is represented by k6. B. The two fitness parameters (conditions) used in the optimisation process to select values for the rates: wide-rage signal going from basal to full expression, and plateaushaped distribution with roughly the same number of cells representing each value in between. C. Simulated signals under the same rate values but different TF numbers: 200 (3MBz case, up) and 3000 (m-xyl case, bottom) with flow cytometry signals shadowed for comparison. D. Several time-course simulations are shown in which the Pm promoter is exposed to three concentrations of its regulator, XylSa: 10 (null induction, thus basal, yellow line), 200 (low induction, green line) and 3000 molecules (high induction, purple line). Centre; the graph corresponds to the rates established in Figure 3C, with k1 = 0.004 and k-1 = 1.5. Top, k1 increased to 250%. Bottom, k1 reduced to 40% its original value. Left, k-1 at 250%. Right, k-1 at 40%.
27
ACS Paragon Plus Environment
Page 29 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
824 825 826
827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843
ACS Synthetic Biology
Figure 4. Analysis of Pm activity noise relative to regulator dynamics.
Study of simulation sets in which regulator number is the only parameter changed. A. Pm transfer function, which measures the output level resulting from different input values. Whereas the average value (red line) shows no additional relevant information, the noise produced by the signal (blue error bars, which denote max-min signal values) displays distinct behaviour according to input numbers, with wider range at middle concentrations. B. Number of GFP molecules over time while input changes (3MBz, m-xyl or none). The three logic values of the signal (‘0’, ‘1a’ and ‘1b’ in Figure 2B, C) are wide-range noise (0-20 h and 40-60 h), small-range high-level noise (20-40 h, 60-80 h) and small-range low-level noise (80 h onwards). C. 24 h simulations at different XylSa levels (from 0 to 3000 molecules) were used to measure 1) the cumulative pulse duration, which corresponds to the length of time that the Pm promoter is in the ON state (regulator bound to the DNA) and 2) signal amplitude, defined here as the the maximum difference (in molecules) between the highest and lowest values achieved during the simulation (measured in steady-state).
28
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
844 845
846 847 848 849 850 851 852 853 854 855 856 857 858 859 860
Page 30 of 32
Figure 5. Effects on promoter activity of its distance from regulator source
A. Left: spatial distribution of simulated proteins following Brownian movement in a cell-like compartment with high protein occupancy. Each coloured line inside the cell represents a protein trajectory from its source (in the middle of the space; labelled S) to its final position at a given time. Density of trajectory positions in each section (longitudinal and transversal) is shown in side graphs. Right: final protein position and their distribution in side graphs. B. Left: simulated spatial trajectories (colour lines) under low protein occupancy. Desnsity of trajectory positions in each section (longitudinal and transversal) is shown in side graphs. Two zoom-in regions (source S and target T) with distinct trajectory points occupation are displayed in detailed. Right: stochastic simulations of Pm-gfp activity with k1 = 0.5 and k-1 = 50 under basal, 200 TFs (blue) and variable TFs (green) conditions. The latter reproduce the noise pattern of 3MBz induction with the new rates and spatial-based constraints.
29
ACS Paragon Plus Environment
Page 31 of 32
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
861 862 863
864 865 866 867 868 869 870 871 872 873 874 875 876 877 878
ACS Synthetic Biology
Figure 6. Pm promoter activity influenced by its distance from regulator source
A. Left: physically rearranged XylS/Pm regulatory node engineered in P. putida KT-BGS (Table 1) to maximise proximity between source (XylS production via Ps promoter) and target (Pm). Both promoters were inserted adjacent to each other into the chromosome of strain KT2440, from which the TOL plasmid was removed (see Methods). Right: flow cytometry results with P. putida KT-BGS cells. As predicted by the model, use of 3MBz as inducer with minimal distance between source and target (top, in green) gives results similar to use of m-xyl as inducer of the reference P. putida mt-2-Pm strain (Figure 2B). Bottom, 3MBz induction in P. putida mt-2-Pm cells, where the TF source and the target promoter are not adjacent. B. Quantitative PCR to measure mRNA molecules transcribed from xylS shows similar Ps promoter activity in both strains, uninduced or 3MBz-induced. C. Visualisation of the single-copy TOL plasmid (J. Kim) At division, the TOL plasmid replicates itself using its own machinery.
30
ACS Paragon Plus Environment
ACS Synthetic Biology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 32 of 32
879 880
Table I. Bacterial strains and plasmids used in this study Bacterial strain or plasmid Strains Escherichia coli CC118pir HB101
Pseudomonas putida KT2440
Relevant characteristics
Reference or source
Cloning host; ∆(ara-leu) araD ∆lacX174 galE galK phoA thiE1 rpsE rpoB(Rif) agrE(Am) recA1 λpir lysogen Helper strain; F- λ- hsdS20(rB- mB-) recA13 leuB6(Am) araC14 ∆(gpt-proA)62 lacY1 galK2(Oc) xyl-5 mtl-1 thiE1 rpsL20(Smr) glnX44(AS)
(Herrero et al., 1990)76
Wild-type strain; mt2 derivative cured of the TOL plasmid pWW0
(Boyer and RoullandDussoix, 1969)77
(Nelson et al., 2002; Bagdasarian et al., 1981)78, 79
mt-2
Wild-type strain bearing pWW0 plasmid
mt-2-Pm
Gmr. P. putida KT2440 inserted in its genomic attTn7 with the hybrid mini-Tn7 delivered by plasmid pBGPm (Supplementary Figure S6) Gmr. P. putida KT2440 inserted in its genomic attTn7 with the hybrid mini-Tn7 delivered by plasmid pBGS (Supplementary Figure S6)
KT-BGS Plasmids pRK600 pTnS-1 pBG
pBG-Pm pSEVA228 pBGS
881 882
Cmr. Helper plasmid used for conjugation; oriV ColE1, RK2(mob+ tra+) Apr, oriR6K, TnSABC+D (Tn7 transposase) operon Kmr, Gmr, oriR6K, mini-Tn7 delivery vector; Tn7L and TnR bracketting a mobile DNA segment for engineering standardised BCD2-msf GFP reporter fusions (Supplementary Figure S6) Kmr, Gmr, oriR6K. pBG inserted with Pm promoter and thus bearing a standardised Pm-gfp fusion Kmr, oriRK2, xylS/Pm expression system Kmr, Gmr, oriR6K, pBG inserted with the xylS/Pm module of pSEVA228 and thus bearing the TF adjacent to the same standardised Pm-gfp fusion as in pBG-Pm
883 884 885
31
ACS Paragon Plus Environment
(Worsey and Williams, 1975)80 This study This study
(Kessler et al., 1994; Keen et al., 1988)67, 81 (Choi et al., 2005)82 (Zoebel et al., 2015)83
This study (Martínez-García et al., 2014)66 This study