Deconvolution of Gene Expression Noise into ... - ACS Publications

Mar 29, 2017 - Transcription Factor−Promoter Interplay. Ángel Goñi-Moreno,. † ... Experiments with engineered bacteria in which this distance is...
1 downloads 13 Views 1MB Size
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 CC118pir 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