A bistable genetic switch based on designable

Full text

Turn on search term navigation

ARTICLE

Received 4 Apr 2014 | Accepted 15 Aug 2014 | Published 29 Sep 2014

Tina Lebar1,2, Urban Bezeljak3, Anja Golob1,3, Miha Jerala3, Lucija Kadunc1,3, Botjan Pir3, Martin Straar3,4, Duan Vuko3, Uro Zupani3, Mojca Benina1,2, Vida Forstneri1, Rok Gaber1,2, Jan Lonzari1,2,Andreja Majerle1,2, Alja Oblak1,2, Ane Smole1 & Roman Jerala1,2

Bistable switches are fundamental regulatory elements of complex systems, ranging from electronics to living cells. Designed genetic toggle switches have been constructed from pairs of natural transcriptional repressors wired to inhibit one another. The complexity of the engineered regulatory circuits can be increased using orthogonal transcriptional regulators based on designed DNA-binding domains. However, a mutual repressor-based toggle switch comprising DNA-binding domains of transcription-activator-like effectors (TALEs) did not support bistability in mammalian cells. Here, the challenge of engineering a bistable switch based on monomeric DNA-binding domains is solved via the introduction of a positive feedback loop composed of activators based on the same TALE domains as their opposing repressors and competition for the same DNA operator site. This design introduces nonlinearity and results in epigenetic bistability. This principle could be used to employ other monomeric DNA-binding domains such as CRISPR for applications ranging from reprogramming cells to building digital biological memory.

DOI: 10.1038/ncomms6007

A bistable genetic switch based on designable DNA-binding domains

1 Department of Biotechnology, National Institute of Chemistry, Hajdrihova 19, Ljubljana 1000, Slovenia. 2 EN-FIST Centre of Excellence, Ljubljana 1000, Slovenia. 3 Slovenian iGEM Team 2012, National Institute of Chemistry and University of Ljubljana, Ljubljana 1000, Slovenia. 4 Faculty of Computer and Information Science, University of Ljubljana, Vena pot 113, Ljubljana 1000, Slovenia. Correspondence and requests for materials should be addressed to R.J. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

Multistability is a widespread natural phenomenon in diverse biological processes, ranging from cell differentiation to the response to environmental signals. The

basic feature of bistable systems is that once an input signal reaches a dened threshold value, the system switches into the corresponding stable state and remains in that state even after this signal is removed. Switching back to the initial state requires the pulse of a second input signal. Flip-op (or SR latch) bistable switches represent the basic elements of nite-state machines composed of sequential logic circuits and the building blocks of memory storage

elements in electronic circuits. In biological systems, genetic toggle switches play important roles in cell cycle progression, embryonal development, cancer, apoptosis and many other processes16. A basic engineered genetic toggle switch consists of two operons encoding a pair of repressors that mutually inhibit one anothers transcription7. This type of design for a genetic toggle switch has been implemented in bacteria8 as well as in mammalian cells9. Although several natural prokaryotic transcription factors have been used to assemble genetic circuits811, the number of orthogonal transcription factors for the toolbox should be

a PCMV

PCMV

TALE-B KRAB

KRAB

TALE-A

PCMV

PCMV PCMV

etr A

PCMV

pir

TagBFP mCitrine

Overlay

Cooperativity 1 Leakage 0.00

Cooperativity 2Leakage 0.05 Leakage 0.00 Leakage 0.05

PI -

900

200

600

100

300

0 0 500

Concentration

0 500

ER -

900

200

600

100

300

0 500

Time (h)

Figure 1 | The mutual repressor switch topology based on TALE repressors does not exhibit bistability. (a) The mutual repressor switch topology, consisting of two TALE repressor genes placed downstream of the constitutive CMV promoter with an operator for the opposing TALE repressor upstream of the promoter, and (b) TALE-based toggle with additional constructs for the external control and two different uorescent protein reporters. The overlay of the microscopy images of the TALE-based mutual repressor switch indicates the simultaneous expression of both uorescent reporter proteins in individual cells in the absence of external regulation. The images are representative of at least three independent experiments. Cells expressing TagBFP or mCitrine reporter genes are displayed in magenta or yellow, respectively. Cells expressing both reports at the same time (53%) are coloured white on the overlay image. Scale bar, 50 mm. The mCitrine and TagBFP reporters were linked to the repressor operons via the t2A peptide linker. For external switch control, the pristinamycin- and erythromycin-inducible transcriptional repressors (PIP:KRAB and E:KRAB, respectively) were added along with the pristinamycin- and erythromycin-inducible transcriptional control of the same pair of TALE repressors. (c) The deterministic modelling of the mutual repressors topology indicates that this setup can achieve state stability only with cooperative elements. The models of the non-cooperative systems did not result in state stability, even with high transcriptional repression as described in Supplementary Note 1. PI and ER stand for pristinamycin and erythromycin induction, respectively. Magenta and green lines indicate the calculated TagBFP and mCitrine reporter gene expression levels, respectively.

2 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

expanded to construct complex designed genetic circuits. Furthermore, the diverse biochemical properties of natural repressors may require adjustments to their strength and orthogonality to make them suitable for designed circuits12. A systematic approach for constructing bistability based on designed regulatory elements is needed.

Modular DNA-binding domains are powerful tools for engineering synthetic cellular regulatory networks because these domains allow the plug-in combinations of DNA-binding and transcriptional effector domains. The DNA-binding domains of zinc nger proteins, transcription-activator-like effectors (TALEs), or CRISPR/Cas nuclease can be genetically fused to transcriptional effector domains to function as designed site-specic transcriptional activators or repressors targeted against almost any DNA sequence1320. The ability to design these synthetic transcription factors to target a large sequential diversity of DNA operators should support the construction of complex designed genetic circuits. Recently, a complete set of static two-input logic gates based on designed TALE repressors was demonstrated in mammalian cells21. It would be desirable to use designable, orthogonal building elements such as TALEs as building blocks in the construction of dynamic, time-dependent logic systems, such as bistable switches. However, one potential problem of designed DNA-binding domains, such as TALEs, zinc nger proteins or CRISPR/Cas nucleases, lies in the monomeric nature of their binding to DNA. Systems theory concludes that a nonlinear response is essential for the functionality of multistable switches9,22. In cases of designed bistable switches implemented up to now, this has been achieved via the application of natural transcription factors that bind as oligomers to the operator DNA. The oligomerization of natural transcription factors evolved separately for each protein, which prevents a systematic modular engineering of orthogonal switches employing the same principles used to construct static logic gates. Therefore, the challenge is to introduce a nonlinear response into a system composed of orthogonal regulators with a linear response introducing a general principle that could be applied across different chassis and for different designable DNA-binding domains.

Here we present a strategy for constructing bistable toggle switches based on designed monomeric DNA-binding domains. We demonstrated that the conventional design of a genetic toggle switch constructed from orthogonal TALE-based repressors did not result in a bistable system. A mathematical model of this design, based on experimentally determined parameters, showed that this was due to the non-cooperativity of TALE domain binding to DNA. Nonlinearity was introduced into the design via the introduction of activators and repressors based on the same DNA-binding domain. Those regulatory elements were wired to introduce competition for the DNA-binding site between an activator and a repressor, and to introduce positive feedback loops. The modelling and experimental validation of this genetic circuit design in human cell lines demonstrate the epigenetic stability of each of the two states of the switch and multiple toggling between the two states via induction by chemical signals, and provides the basis for the scalability for generating multiple switches and more complex logic circuits.

ResultsToggle switch based on two TALE repressors is not bistable. Our initial goal was to prepare a mutual repressor switch comprising a pair of designed orthogonal TALE domain-based repressors and test its bistability in mammalian cells. The designed repressors were constructed via fusion of the selected TALE DNA-binding domains with the KRAB chromatin-silencing domain23. The position of the KRAB repression domains relative to the

DNA-binding domain was investigated to determine the optimal arrangement of the effector domain regarding the DNA-binding domain (Supplementary Fig. 1a,b). Multiple copies of the binding site upstream from the promoter were used21 to decrease the off-target effects and enhance the efciency of repression. This resulted in up to 98% transcriptional repression, while the two selected TALE domains were mutually orthogonal and exhibited virtually no cross-talk (Supplementary Fig. 1a,b). To implement the toggle

c a c a

P TALE-C

TALE-B

d b d b

TALE-D

TALE-A

PI ER

2.0

Concentration

1.6

1.2

0.8

0.4

0.0

0 2 4 6 8 10 12 14 16

Time (days)

TALE-A

TALE-B

TALE-A

2.0

1.6

1.2

0.8

0.4

0.0

Concentration

0.0

0 2 4 6 8 10 12 14 16

PI ER PI

2.0

1.6

0.4

1.2

0.8

0.4

0.0

0 2 4 6 8 10 12 14 16

Time (days)

Figure 2 | Stochastic models demonstrate that competition for the binding sites between an activator and a repressor is an essential feature required for the bistability of the switch based on the monomeric binding of transcription factors. (a) Introduction of a positive feedback loop into the TALE-based mutual repressor switch by the addition of two TALE activator genes. Each of the TALE activators controls expression of its own gene and the corresponding repressor gene, while the latter also inhibits expression of the opposing TALE activator. Each TALE activator or repressor in this setup has a unique DNA-binding domain, eliminating the competition between activators and repressors. (b) Simulations of this genetic circuit demonstrate that the addition of a positive feedback loop is not sufcient for state stability. Details of the model are available in Supplementary Note 2. (c) Introduction of binding competition into the TALE-based toggle switch by the introduction of two TALE domains without regulatory domains. Each of the introduced TALE DNA-binding domains can bind to its own operon and to the operon of the opposing TALE repressor, introducing competition but not exerting direct operon activation or repression. (d) Simulations of this topology reveal that competition for binding sites is sufcient for the switch to assume a stable state in the presence of an inducer molecule. Details on the model are available in Supplementary Note 2.

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

switch, the TALE repressor coding region fused to the carboxy-terminal KRAB domain was positioned under the control of the opposing DNA operator (Fig. 1a). Two uorescent reporters, blue uorescent protein and monomeric citrine protein, were genetically fused to each of the repressor operons via the t2A peptide linker24 to help monitor the expression of each repressor in individual cells. For the external time-dependent selection of each of the two states, we introduced an input layer to regulate the transcription of selected TALE repressors under the control of chemical inducers (that is, pristinamycin and erythromycin)25,26 (Fig. 1b and Supplementary Fig. 1c,e). The pristinamycin- and erythromycin-inducible systems coupled to the TALE repressors exhibited low leakage and an approximately vefold difference in reporter expression between the stimulated and unstimulated states (Supplementary Fig. 1d,f). Experimental results regarding the TALE repressors and the induction layer provided parameters for the mathematical model (Supplementary Notes 1 and 2). The complete setup of the TALE-based mutual repressor switch consisted of the induction layer regulating two mutual repressors, TALA:KRAB and TALB:KRAB (Fig. 1b). However, stimulation of mammalian cells transiently transfected with plasmids for this switch setup with either of the two chemical inducers degenerated into a simultaneous expression of both reporters, suggesting that this system cannot retain the stable state (data not shown). As an additional test, we examined cells with transiently introduced genetic switch circuit plasmids in the

absence of external chemical induction. In a bistable system, each individual cell would be expected to adopt one of the two stable states based on the stochastic variations of each repressor level or based on the previous state of each cell. However, uorescence microscopy analysis demonstrated that the majority of cells expressed both reporters within each cell simultaneously (Fig. 1b). We could not distinguish between the simultaneous transcription of both reporters and dynamic switching between both states in individual cells. Nevertheless, both explanations lead to the conclusion that this design of a mutual TALE-repressor-based genetic circuit is not bistable.

Although most natural transcription factors bind to DNA as oligomers (that is, dimers or tetramers), the crystal structure of the TALE DNA-binding domain clearly shows that it binds to the target DNA as a monomer27,28; thus, the cooperativity of binding should be low. The same argument is also valid for other designable DNA-binding domains such as zinc nger- and CRISPR/Cas-based repressors. Deterministic models of the TALE-based mutual repressor switch demonstrated that this system could not retain the stable state and drifted to the mixed state (Fig. 1c and Supplementary Note 1). Even in the case of the cooperativity coefcient of 2, this type of switch is sensitive to as little as 5% incomplete repression.

Introduction of nonlinearity using monomeric binding domains. The theoretical solution of introducing nonlinearity

0 0 2 4 6 8 10 12 14 16

a c

positive feedback loopconstitutive CMV promoter

PI ER PI

Time (days)

TALE-A KRAB TALE-A VP16

Binding site A Promoter

PCMV

TALE-A

TALE-B

PCMV

TALE-B

TALE-A

PCMV

Activity (RLU) Activity (RLU)

a Pmin

Concentration

TALA:KRAB (ng)

0 2.5 5 10 25 TALA:VP16 (ng)[a]-[Pmin]_fLuc [75 ng]

0 0 2 4 6 8 10 12 14 16

b Pmin

TALB:KRAB (ng)

0 2.5 5 10 25 TALB:VP16 (ng)[b]-[Pmin]_fLuc [75 ng]

Figure 3 | Introduction of nonlinearity into a genetic circuit based on monomeric DNA-binding elements employing competition and positive feedback loops. (a) A TALE repressor and a TALE activator based on the same DNA-binding domain competing for binding to the same DNA target site. (b) The competition of the TALE-based activator and repressor for the same binding site was experimentally tested using the luciferase reporter. Various ratios of plasmids constitutively expressing a TALE activator and a TALE repressor based on the same DNA-binding domain were co-transfected into HEK293T cells, along with 75 ng of the corresponding reporter plasmid, which resulted in different levels of the expression of the reporter. The values represent the mean and s.d. of three replicates and are representative of at least three independent experiments. Fit of the model based on parameters extracted from the experimental data of all experiments overlays the experimental data. (c) Improved bistable switch design incorporating binding competition, along with positive feedback loops. Each TALE activator controls the expression of its own gene and the corresponding repressor gene, while the latter inhibits the expression of the opposing ranch with TALE activator and repressor. (d) A stochastic model of the switch with a positive feedback loop. Modelling determined that the improved design should be able to adopt one of the two stable states and toggle between them on the input signal pulse. The details of the model are available in Supplementary Note 2. PI and ER stand for pristinamycin and erythromycin, respectively. The magenta and green lines indicate the calculated TagBFP and mCitrine reporter gene expression levels, respectively.

4 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

a b

mCitrine (FI)

16%

102

101

etr

PCMV

A A

B B

pir

PCMV

mCitrine

100

15.5%

60%

100 101 102

TagBFP (FI)

TagBFP

Overlay

c d

PI added

PI removed ()

etr

PCMV

A A

B B

pir

PCMV

etr

PCMV

A A

B B

pir

PCMV

PI PI

f g

Both TagBFP mCitrine

3 days

6 days 9 days

100

Cell number (%)

ConcentrationConcentration

TagBFP mCitrine

0 0 2 4 6 8 10 12 14 16

0 3 6 9 6 9days PI PI

Time (days)

PI (3 d) (9 d)

PI (9 d)

mCitrine (FI)

102

101

100

ER added

ER removed ()

100 101 101 TagBFP (FI)

etr

PCMV

A A

B B

pir

PCMV

etr

PCMV

A A

B B

pir

PCMV

ER ER

k l

Both TagBFP mCitrine

3 days

6 days

9 days

6 days 9 days

100

Cell number (%)

TagBFP

mCitrine

0 0 2 4 6 8 10 12 14 16

0 3 6 9 6 9days ERER

Time (days)

ER (3 d) (9 d)

ER (9 d)

mCitrine (FI)

102

101

100

100 101 102 TagBFP (FI)

Figure 4 | The improved design of the TALE-based bistable genetic switch displays bistability. (a) Introduction of binding competition, along with positive feedback loops, into the TALE-based mutual repressor switch. The TALE activators contain the same DNA-binding domains as the TALE repressors of the opposite state. External control was provided in the same fashion as for the mutual repressor switch (Fig. 1b), with the addition of two inducer-dependent TALE activator operons. (b) Confocal microscopy and ow cytometry analyses of the TALE-based competitive feedback switch in the absence of external regulation demonstrate bimodal uorescence distribution. Cells expressing either a TagBFP or mCitrine reporter are marked magenta or yellow, respectively. Scale bar, 50 mm. Schematic presentation of the pristinamycin (PI) (c) and erythromycin-induced (ER) (h) state, and the subsequent stable state on PI (d) and ER removal (i). Confocal microscopy (e,j) and ow cytometry (f,k) analyses of the PI-induced and stable ( ) states (e,f) and of the ER-

induced and stable ( ) states (j,k). Scale bar, 250 mm. The fraction of cells expressing specic reporters obtained from ow cytometry is presented in the

histograms. (g,l) Stochastic modelling determined that the improved design is able to retain a stable state, even after the removal of an inducer: (g) removal of pristinamycin (PI to ) and (l) erythromycin (ER to ) as described in Supplementary Note 2. Magenta and green lines indicate the calculated

TagBFP and mCitrine reporter gene expression levels, respectively. Results are representative of at least three independent experiments.

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

into a system composed of monomeric DNA-binding elements proposed the introduction of a pair of hypothetical transcription factors, each of which should be able to activate its own expression and simultaneously repress the transcription of the opposing transcription factor29,30. Such types of transcription factors would be very difcult to design. However, the modularity of the designed transcription factors allows the construction of a transcriptional activator or repressor, each targeting the same DNA-binding sequence, by coupling the DNA-binding domain to either the VP16 transactivation domain (Supplementary Fig. 2a,b) or the KRAB repression domain, respectively. We reasoned that such a pair of transcriptional regulators based on TALEs could be used to introduce nonlinearity into the designed switch circuit using two strategies: (a) the introduction of positive feedback loops and (b) competition between an activator and a repressor for the same DNA operator binding site. Simulations of a cellular circuit employing TALE activators and repressors wired to bind to separate DNA sites upstream from the promoter, thus without competition effects (Fig. 2a), did not result in bistability (Fig. 2b and Supplementary Note 2). As the designed TALE-based repressor and activator can be based on the same DNA-binding domain, they could compete for binding to the target DNA site (Fig. 3a). The relative occupancy of the binding site and the resulting gene transcription should therefore depend on the ratio of the repressor and the activator, as they comprise the same DNA-binding domain (Fig. 3b). We performed simulations of the designed circuit using competition between TALE-based repressors and TALE DNA-binding domains without an

effector domain (Fig. 2c), which demonstrated the functional bistability of this switch circuit topology (Fig. 2d and Supplementary Note 2). This shows that the competition between the activator and the repressor for the same DNA site is an essential feature for introducing nonlinearity into a genetic circuit based on non-cooperative regulators.

Maintenance of the stability in a TALE-based genetic switch. For the experimental implementation of the bistable switch, we combined the positive feedback loops and competition effects into the same circuit (Fig. 3c). This circuit topology requires only one pair of designed DNA-binding (TALE) domains to construct two repressoractivator pairs, requiring only one DNA operator sequence for each of the two branches of the switch device. Deterministic (Supplementary Fig. 3 and Supplementary Note 1) and stochastic (Fig. 3d) system analysis showed that the switch design that incorporated positive feedback loops with an activator and binding competition performed best. The state of the switch was regulated externally by the same pair of chemical inducers (pristinamycin and erythromycin) used for the mutual repressor switch, with the addition of two TALE activators under the control of chemically inducible transcription factors (Fig. 4a and Supplementary Fig. 4). Both inducible systems, when coupled with TALE activators (Supplementary Fig. 2c,e), exhibited low leakage and high differences in the reporter activity between the stimulated and unstimulated states (Supplementary Fig. 2d,f). Experiments on

a b c

PCMV PCMV

Pmin

PCMV

200

150

100

0 0 50 100 150 200

Repressor

1.00.90.80.70.60.50.40.30.20.10.0

100

101

Stability

Activator

200

150

100

00 50 100 150 200

Repressor

10 98 76 54 32 10

Time until switch

Activator

200

150

100

00 50 100 150 200

0.90.80.70.60.50.40.30.20.10.0

TAL A dissociation

1.00.90.80.70.60.50.40.30.20.10.0

Stability (composite) Stability (composite)

Activator

Stability

102

103

104

10 104

105

TAL B dissociation

103

102

101

100

Repressor

d e f

Pmin

100

Activator

200

150

100

0 0 50 100 150 200

Repressor

1.00.90.80.70.60.50.40.30.20.10.0

Stability

Activator

200

150

100

0 0 50 100 150 200

Repressor Repressor

10 98 76 54 32 1 0

Time until switch

Activator

200

150

100

0 0 50 100 150 200

1.00.90.80.70.60.50.40.30.20.10.0

101

TAL A dissociation

1.00.90.80.70.60.50.40.30.20.10.0

Stability

102

103

104

10 104

105

TAL B dissociation

103

102

101

100

Figure 5 | The minimal promoter competitive feedback switch is more robust to variations in the repressor:activator ratio than the switch based on the constitutive (cytomegalovirus (CMV)) promoter. Parameter scans of stochastic topology models regarding the concentration of activator and repressor. The evaluation criteria were calculated after 17 simulation days. Stability is dened as R2R1, where R2 is the concentration of an induced state reporter and R1 is the concentration of the repressed state reporter. This criterion is normalized for its maximal value across the scanned region. The rst inducer was introduced at t0 0 and removed after t1 3 days. The time until the switch is dened as t2 t1 (in number of days), where t1 represents the

time of the second inducer introduction and t2 the time until the R2 concentration becomes larger than R1. The rst inducer was introduced at t0 0 and the

second inducer at t1 3 days. White areas represent parameter values for which the simulations predict instability (negative stability value) and/or the

inability to switch between the states (R2 never becomes greater than R1). Amounts of activators and repressors were varied as dosage of the constructs PMIN/CMVA0, PMIN/CMVB0, PCMVPIP0, PCMVE0 and PMIN/CMVA1, PMIN/CMVB1, PCMVPIP1, PCMVE1. Details on the constructs are available in

Supplementary Tables 8 and 9. (a,b) Evaluation of the stability and state switching ranges for the CMV promoter competitive feedback switch (Fig. 2c). Simulations indicate that toggling between states is possible only within a narrow range of repressor:activator ratios for topology based on PCMV. (c) Composite image of a and b, displaying the range of robust stability and switching. (d,e) Evaluation of stability and state switching ranges for the minimal promoter competitive feedback switch (Fig. 5a). The range of positive stability and state switching is signicantly wider than for the CMV promoter competitive feedback switch (a,b), displaying higher robustness to variations in the activator:repressor ratios. (f) Composite image of d and e, displaying range robust stability and switching. Details on the model are available in Supplementary Note 2. (g,h) Scan of the effect of the afnities of the TALE-binding domains displaying the similar effect on stability of both types of switches.

6 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

the activation of individual branches of the switch provided parameter values to model the optimal ratio of the activator/ repressor operons (Fig. 5 and Supplementary Note 2) and test the experimental implementation of the competitive feedback switch. In principle, an activator and repressor targeted against different binding sites could also be joined in a single t2A sequence-linked reading frame. However, simulations demonstrated that the separate operons provide an additional tool for tuning the performance of the switch that the optimal ratio of activator versus repressor is different from 1:1.

Cells harbouring the improved bistable switch device assumed a bimodal distribution in the absence of external induction, characterized by the expression of a single uorescent reporter in

the majority (60% TagBFP, 16% mCitrine and only 9% co-expressed both reporters) of the transfected individual cells (Fig. 4b). It has been previously observed that TALE domains can exhibit variations in terms of binding afnity31. The slightly stronger repression of TALB:KRAB compared with TALA:KRAB, as observed in the titration experiment (Fig. 3b), caused a strong bias (1:4) in the state selection in the absence of an inducer. Nevertheless, the overwhelming population of cells expressing a single reporter rather than both reporters simultaneously implied that such a system is bistable. Transient induction via either of the two chemical inducers shifted the system towards the expected state, with a very low level of expression of the opposing state (Fig. 4cf,hk). Despite the imbalance in the repression strength

PCMV

Induction system leakage Pmin

Pmin

1.00.90.80.70.60.50.40.30.20.10.0

100

10 9 8 7 6 5 4 3 2 1 0

Time until switch

101

PIR leakage

Stability

PIR leakage

Time until switch

PIR leakage

1.00.90.80.70.60.50.40.30.20.10.0

101

Stability

PIR leakage

10 9 8 7 6 5 4 3 2 1 0

Time until switch

102

103

104

10 104

104

10 104

104

10 104

104

10 104

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

ETR leakage

Association

Degradation

Effect on transcription

PCMV

100

TAL A association

TAL A degradation

TAL A effect

1.00.90.80.70.60.50.40.30.20.10.0

TAL A association

10 9 8 7 6 5 4 3 2 1 0

TAL A association

10 9 8 7 6 5 4 3 2 1 0

101

Stability

1.00.90.80.70.60.50.40.30.20.10.0

Stability

Time until switch

102

103

104

10 104

104

10 104

104

10 104

104

10 104

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

TAL B association

PCMV

100

1.00.90.80.70.60.50.40.30.20.10.0

10 9 8 7 6 5 4 3 2 1 0

100

TAL A degradation

Time until switch

TAL A degradation

1.00.90.80.70.60.50.40.30.20.10.0

10 9 8 7 6 5 4 3 2 1 0

101

102

Stability

TAL A degradation

Time until switch

102

103

104

10 104

104

10 104

104

10 104

104

10 104

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

105

103

102

101

100

TAL B degradation

100

TAL B degradation

1.00.90.80.70.60.50.40.30.20.10.0

101

Stability

TAL A effect

10 9 8 7 6 5 4 3 2 1 0

Time until switch

TAL A effect

1.00.90.80.70.60.50.40.30.20.10.0

TAL A effect

10 9 8 7 6 5 4 3 2 1 0

Time until switch

100

Stability

101

102

101

100

101

100

101

100

101

100

101

TAL B effect

Figure 6 | Sensitivity of the improved bistable switch to different TALE-regulator-dependent parameters. Parameter scans of stochastic simulation of topology models based on Pmin or PCMV promoters as a function of the TALE regulator-dependent parameters. The simulation setting with evaluation criteria for stability and switching times is dened as in Fig. 4. White areas represent parameter values for which the simulations predict instability (negative stability value) and/or the inability to switch between the states. Parameter labels refer to Supplementary Table 7. (a) Parameter scans of the leakage of the PIP/ETR induction systems (lCMV ). The parameter does not signicantly change any evaluation criterion in either switch model.(b) Parameter scans of the leakage of the TALE association constants (kT). The parameter does not signicantly change any evaluation criterion in either switch model. (c) Parameter scans of the TALE-A/TALE-B degradation rates (dT). The simulations predict a slightly higher stability and switching range of the switch based on PCMV. (d) Parameter scans of TALE-A/TALE-B effect on transcription. To account for variation in the TALE-A/TALE-B effect, the parameters were multiplied by factors, ranging from 10 2 to 101, as shown on the parameter axes. The variation signicantly affects the stability and the switching ability of the PCMV-based switch, while the PMIN variant displays higher robustness in this respect. Details on the models are available in

Supplementary Note 2.

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

of the two TALE repressors, cells retained either of the two selected states equally well for more than 6 days after the chemical inducer was removed (Fig. 4e,f,j,k). This indicates that the improved switch design functions as predicted by the model (Fig. 4g,l). Modelling indicated that the system is insensitive to the several orders of magnitude of difference in DNA-binding strength between the pair of TALE domains due to the positive feedback loop (Fig. 5 and Supplementary Note 2).

Functional switching between stable states. Although the improved switch topology retained the epigenetic state after the chemical inducer was removed, a functional bistable switch should also support controlled reversible toggling between the two stable states. We tested this property via the initial induction of one state, followed by the replacement with the second chemical inducer. Experimental results showed the expected trends in expression of the two reporters 6 days after the inducer

0 2 4 6 8 10 12 14 16

PI added

PI removed ()

etr

PCMV PCMV

PCMV

A A

B B

pir

Pmin Pmin

etr

PCMV PCMV

PCMV

A A

B B

pir

Pmin Pmin

PI PI

3 days

6 days 10 days 13 days

2.0

TagBFP mCitrine

Concentration

1.6

1.2

0.8

0.4

0.0

Time (days)

PI (3 days)

PI (15 days)

100

(15 days)

100

Cell number (%)

mCitrine (FI)

102

101

mCitrine (FI)

101

Both TagBFP mCitrine

100

0 3 6 9 12 15 Time (days)

100 101 102

TagBFP (FI)

100 101 102

TagBFP (FI)

ER added

ER removed ()

etr

PCMV PCMV

PCMV

A A

B B

pir

Pmin Pmin

etr

PCMV PCMV

PCMV

A A

B B

pir

Pmin Pmin

ER ER ER

3 days

6 days 10 days 13 days

2.0

TagBFP mCitrine

Concentration

1.6

1.2

0.8

0.4

0.0

0 2 4 6 8 10 12 14 16

Time (days)

100

(15 days)

100

ER (3 days)

ER (15 days)

mCitrine (FI)

Cell number (%)

mCitrine (FI)

Cell number (%)

102

101

Both TagBFP mCitrine

100

100 101 102 TagBFP (FI)

0 3 6 9 12 15 Time (days)

100 101 102

TagBFP (FI)

0 3 6 9 12 15 Time (days)

Figure 7 | The competitive feedback switch with minimal promoters exhibits state stability in the absence of an inducer. The presence of pristinamycin (PI) or erythromycin (ER) triggers expression of the selected state, as depicted in the schematic representation (a,f), and conrmed with confocal microscopy (c,h) and ow cytometry (d,i). Removal of the inducer molecule 3 days after transfection resulted in stable states as depicted in the schematic representation (b,g), and conrmed even after 13 days with confocal microscopy (c,h) and ow cytometry (d,i). The microscopy and ow cytometry results are representative of at least three independent experiments. Scale bar, 50 mm. The number of cells expressing specic reporters (as obtained from ow cytometry) is presented in the histograms. A stochastic model of the pristinamycin-induced state (e) or the erythromycin-induced state (j) results in state stability, matching the experimental results. Details on the model are available in Supplementary Note 2.

8 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

replacement, but the emergence of the second state was too slow for the reliable conrmation (Supplementary Fig. 5).

To increase the difference between the levels of transcription for the repressed and activated states, the constitutive promoter of the core switch was replaced with a minimal promoter, which exhibits essentially no reporter gene expression in the absence of an activator or in the presence of a repressor. However, the minimal promoter setup could not be used for the basic toggle switch, as it requires repression. Modelling indicated that a system based on minimal promoters is more robust, that is, it is even less sensitive to the incomplete repression than the competitive feedback switch based on the constitutive promoter (Fig. 5 and Supplementary Note 2) and is able to transition between the states more quickly (Fig. 6 and Supplementary Fig. 6).

The experimental results of the switch based on the minimal promoter clearly demonstrated state stability even 12 days after the inducer was removed (Fig. 7). Cells harbouring minimal promoter-based competitive feedback switch plasmids (Fig. 8a) assumed a bimodal distribution in the absence of external induction, characterized by the expression of a single uorescent reporter in the majority (80% TagBFP, 19% mCitrine and only 1% co-expressed both reporters) of the individual cells expressing TagBFP or mCitrine (Fig. 8b). State transition was monitored from both initial conditions via the confocal microscopy and ow cytometry (Fig. 7). However, some of the initial uorescent reporter proteins remained after the transition, presumably due to the slow degradation of uorescent proteins, making them a less suitable choice as reporters for dynamic experiments (Supplementary Fig. 7). To overcome this drawback, we used a

mCitrine (FI)

102

101

B A

etr

PCMV

A A

B B

pir

PCMV

Pmin

mCitrine 100

100 101 102

TagBFP (FI)

TagBFP

Overlay

PI added

ER added

B A

etr

PCMV

A A

B B

pir

PCMV

Pmin

B etr

etr

PCMV

A A

B B

pir

PCMV

Pmin

Switch

PI ER

ER PI

ER PI ER

0.2

2.0

0.2

2.0

Activity (RLU)

Concentration

Activity (RLU)

Concentration

0.2

1.6

1.2

0.1

0.8

0.1

0.8

0.4

0.0

3 6 9 1317 Time (days)

0 2 4 6 8 10 12 14 16

3 6 9 1317Time (days) Time (days)

3 6 9 1317 Time (days)

3 6 9 1317

3 6 9 1317 Time (days)

0 2 4 6 8 10 12 14 16

Time (days)

f g

Switch ER

Switch PI

Switch ER

Switch PI

0.05

Concentration

1.2

0.2

Activity (RLU)

1.2

0.8

Activity (RLU)

Concentration

0.2

0.1

0.8

0.1

0.4

0.0

0.00

0.0

0 2 4 6 8 10 12 14 16

Time (days)

Figure 8 | Regulated transitions between the two stable states of the switch. (a) Introduction of minimal promoters into the TALE-based competitive feedback switch. The mutual TALE repressors and the corresponding TALE activators of the feedback loop are placed under the control of a minimal promoter. (b) Confocal microscopy analysis of the TALE-based competitive feedback switch in the absence of external regulation indicates bimodal distribution of uorescent proteins. Cells expressing either a TagBFP or mCitrine reporter are marked in magenta or yellow, respectively. Cells expressing both reports at the same time (3%) are marked white on the overlay image. Scale bar, 50 mm. (c) Schematic representation of toggling between the pristinamycin (PI)- and erythromycin (ER)-induced states. (dg) Switching between stable states was experimentally tested using luminometry (mean with s.d., n 4). The purple bar graphs represent rey luciferase activity in a switch in which the luciferase gene was used to replace the TagBFP reporter, while

the yellow bar graphs correspond to the version with the luciferase gene replacing the mCitrine reporter. Switching between the states was tested in both directions, (d) pristinamycinerythromycin (PI-ER) and (e) erythromycinpristinamycin (ER-PI), as well as for multiple switch events: (f) pristinamycin erythromycinpristinamycin and (g) erythromycinpristinamycinerythromycin. Experimental results are in agreement with the stochastic models (shown on the right side of the panels) and demonstrate transitions between the on and off states depending on the presence of the selected chemical inducer as described in Supplementary Note 2. Results are representative of at least three independent experiments.

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

PI added

PI removed

pir

PCMV

a Pmin

Pmin

etr B

A C

PCMV

pir PCMV

PCMV

C B

Pmin

tet

PCMV

c Pmin

A B

etr

tet

Pmin

A C

Pmin

A C

etr B

A C

PCMV

c A C

Pmin B

d e f

TE PI ER TE PI ER PI

Concentration

2.4

1.6

2.0

Concentration

1.6

1.2

0.8

0.4

0.0

0 2 4 6 8 10 12 14 16

0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16

0 2 4 6 8 10 12 14 16

Time (days)

Figure 9 | Design of a tristable genetic transcription circuit regulated by three different input signals. (a) Schematic representation of the circuit components comprising three pairs of chemical inducer sensitive regulators, combinations of two repressors that repress transcription of two other constructs and an activator that activates its own transcription; (b) state of the switch in the presence of one chemical inducer (PI). Other two states are regulated in the same manner. (c) State of the switch when the particular inducer has been removed. (d) Stochastic simulation of the switch in the presence of each inducer for 3 days and its removal afterward demonstrates the stability of the circuit. (e) Simulation of switching from the pristinamycinto the erythromycin-induced state, followed by switching to the tetracyclin-induced state (PI-ER-TET). (f) Simulation of switching from the pristinamycin- to the erythromycin-induced state, followed by switching back to the pristinamycin-induced state (PI- ER-PI). The design enables toggling between all three states, regardless of their order and the initial state. Details on the model are available in Supplementary Note 2.

rey luciferase reporter, which has a higher turnover rate instead. This nal setup of the competitive feedback loop switch demonstrated successful multiple toggling of a TALE-based switch between the on and off states in both directions over the course of 17 days (Fig. 8dg).

DiscussionDesignable modular DNA-binding domains provide a powerful enrichment of the toolbox for the rational regulation of biological processes21,32,33. The already demonstrated ability to generate tens of thousands of designed DNA-binding domains targeting different DNA sequences34,35 represents the great potential for the designed regulatory circuits of increased complexity and orthogonality. In this report we show that in addition to the static logic circuits, designable DNA-binding domains can be used to construct dynamic logic devices. Here we present a general strategy for implementing the designed DNA-binding domains to construct bistable switches, which has strong scaling potential due to the modularity and orthogonality. The main feature of our design is the introduction of a nonlinear response, which is required to construct complex dynamic regulatory circuits, through the introduction of positive feedback loops and competition between the activator and the repressor for binding to the same DNA operator site. Competition and feedback loops in transcriptional regulatory circuits have been already analysed previously36,37; however, those designs incorporated cooperative transcriptional regulators or regulator features that would be very difcult to implement experimentally.

The design of a bistable switch based on linear response elements has been experimentally demonstrated in a few special cases, such as using a T7 RNA polymerase that activates its own expression. In that case, the retardation of cell growth due to the circuit activation caused nonlinear dilution of the T7 RNA

polymerase in individual bacteria, which, combined with the feedback, resulted in bistability38. A bistable switch with non-cooperative elements has also been demonstrated based on sequestration and positive feedback using sigma and anti-sigma factors39. In contrast to those designs, which are based on few selected biological phenomena or molecules, our platform is based on an extremely large pool of orthogonal DNA-binding elements that can be applied to the design of multiple orthogonal switches that could operate within the same cell or more complex switches. As a demonstration we designed a tristable switch based on three different inputs, which is an extension of a bistable switch and also involves the positive feedback loop and competition (Fig. 9). In addition, the tools developed in this report allow design of a tristable switch, regulated by different combinations of only two input signals (Fig. 10). This switch device requires four TALE-binding domains to build a topology that can distinguish between the rst, second or both input channel pulses simultaneously and allows reversible switching between any combinations, extending the ability to control cell states. Simulation of those switches based on the parameters from an experimentally veried bistable TALE-based switch indicates that the design of those more complex switches is functional. The designable orthogonal DNA-binding modules allow the construction of parallel orthogonal bistable switches within a single cell, requiring only one pair of orthogonal designed DNA-binding domains per bit of cellular memory. With the ability to massively generate different TALE domains or sgRNA/Cas combinations for the CRISPR system, the design presented here therefore expands the potential regulatory and memory capacity of biological systems. Different input signals from the environment or internal cellular processes (for example, cell type-characteristic transcription factors40) that could be converted into the transcriptional signal could be used to regulate the designed circuits for different applications, for example, sensing or therapeutic.

10 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

pir

PCMV

Pmin

pir

etr

A D

pir etr PCMV B C

A D

pir etr PCMV B C

A D

pir etr PCMV B C

A D

PI added

d pir PCMV B C A

a Pmin B A

2.4

1.8

1.2

0.6

c Pmin B

A C

c Pmin B

A C

c Pmin B

A C

0.0

d etr PCMV B

b Pmin B

A C

b Pmin B

A C

b Pmin B

A C

c ER added

d pir PCMV B C A

2.4

0 2 4 6 8 10 12 14 16

Concentration

1.8

1.2

0.6

0.0

d etr PCMV B

0 2 4 6 8 10 12 14 16

d PI and ER added

PI+ER

a Pmin B A

2.4

1.8

1.2

0.6

0.0

d etr PCMV B

0 2 4 6 8 10 12 14 16

e f

PI ER PI +ER PI+ER ER PI+ER Time (days)

1.6

Concentration

1.2

0.8

0.4

0.0

0 2 4 6 8 10 12 14 16

Time (days) Time (days)

Figure 10 | Design of a tristable genetic transcription circuit regulated by different combinations of two input signals. (a) Schematic representation of the circuit components, comprising four TALE-based activators, repressors and their binding sites. State 3 is induced in the presence of two inducers simultaneously. As the presence of both inducers can potentially induce each of the three states, an additional TALE repressor is included (TALE-D)to repress induction system constructs for pristinamycin- and erythromycin-induced states; (bd) states of the switch in the presence of PI, ER and both PI and ER, respectively. Simulations (right panels) represent states in the presence of the appropriate inducers for three days and their removal afterwards; (e) switching from the pristinamycin- to the erythromycin-induced state, followed by the addition of pristinamycin and erythromycin (PI-ER-PI ER);

(f) switching from the two inducer-induced state to the erythromycin-induced state and back to the two inducer-induced state (PI ER - ER-PI ER).

The design enables toggling between all three possible states, regardless of their order and the initial state. The transition from the two inducer-i nduced state is possible regardless of the additional TALE-D repressor presence. Details on the model are available in Supplementary Note 2.

In conclusion, orthogonally designed DNA-binding domains can be used for the dynamic cell regulation through constructing multistable genetic switches. The availability of a very large number of designable orthogonal DNA-binding domains represents the foundation for the construction of parallel orthogonal switches, each requiring only two DNA-binding domains to assume the selected epigenetic memory states. The introduction of nonlinearity into the system based on linear responseregulatory elements based on the principles used here could be applied to other biological information mediators that exhibit non-cooperative binding, particularly for other modular DNA-binding domains, such as zinc ngers, CRISPR/Cas or helix-coil-helix domains. Finally, although this type of regulation has been demonstrated in mammalian cells, it is transferable to other cell chassis through the use of appropriate activation/repression domain modules.

Methods

Plasmid construction. All plasmids used were constructed with classic molecular cloning methods or the Gibson assembly method41 and are listed in Supplementary Table 1. The nucleotide and amino acid sequences are provided in Supplementary Table 2. The TALE DNA-binding domains (TALA, aa 2804 and TALB, aa 2804) were PCR amplied from TALEN1257 (plasmid 32280, Addgene) and TALEN1297 (plasmid 32279, Addgene)42, respectively. The KRAB domain was obtained from pLVPT-rtTR-KRAB-2SM2 (plasmid 11652, Addgene)43 and the VP16 domain was

obtained from the vector pSGVP, which was provided by Professor Dr Mark Ptashne (Memorial Sloan Kettering Cancer Center, New York, NY, USA)44. TALA and TALB operators21 were designed as multiple repeats of their target sequences, separated by hypervariable linkers, and synthesized by Life Technologies.

Erythromycin- and pristinamycin-dependent DNA-binding proteins E and PIP (plasmids pWW35 and pMF207), as well as the corresponding DNA-binding sites etr and pir (plasmids pWW56 and pMF208), were provided by Professor Dr Martin Fussenegger (Department of Biosystems Science and Engineering, ETH Zrich, Switzerland)25,26. The reporter genes were PCR amplied from commercial plasmids: rey luciferase from pGL4.16 (Promega), TagBFP from pTagBFP-N (Evrogen) and mCitrine45 obtained as a gift from Dr Oliver Griesbeck (the Max Planck Institute of Neurobiology, Munchen, Germany). Renilla luciferase (phRLTK, Promega), mCherry (pmCherry-C1, Clontech) and iRFP670 (ref. 46) (synthesized by Life Technologies) were used as the transfection controls. The SV40 large T-antigen nuclear localization sequence, the hexa-histidine tag and the t2A sequence47 were introduced into the constructs with PCR.

Cell culture. The human embryonic kidney (HEK) 293 T cell line (American Type Culture Collection, Manassas, VA, USA) was cultured in DMEM medium (Invitrogen) supplemented with 10% fetal bovine serum) (BioWhittaker, Walkersville, MD, USA) at 37 C in a 5% CO2 environment.

Inducer molecules. Erythromycin (Sigma-Aldrich) and pristinamycin (Sano-Aventis) were prepared as stock solutions of 50 mg ml 1 (erythromycin in ethanol, pristinamycin in dimethyl sulphoxide) and used at a nal concentration of2 mg ml 1.

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007

Luciferase reporter assay. HEK293T cells (2 104 cells per well) were seeded in

CoStar White 96-well plates (Corning). At 5070% conuence, they were transfected with a mixture of DNA and JetPEI according to the manufacturers instructions (Polyplus-transfection). Every 3 days, the cells were supplemented with fresh medium and (where appropriate) inducer molecules. At various time points, the cells were harvested and lysed with 30 ml of 1 passive lysis buffer

(Promega). Firey luciferase and Renilla luciferase expression were measured using a dual luciferase assay (Promega) and an Orion II microplate reader (Berthold Technologies, Pforzheim, Germany). Relative luciferase activity was calculated by normalizing each samples rey luciferase activity to the constitutive Renilla luciferase activity determined within the same sample. The graphs represent the mean and s.d. from the mean of at least three measurements within an experiment. The data are representative of at least three independent experiments.

Flow cytometry. HEK293T cells were seeded in 24-well plates (1 105 cells per

well) or in 12-well plates (5 105 cells per well) (TPP, Trasadingen, Switzerland).

At 5070% conuence, they were transfected with a mixture of DNA and JetPEI according to the manufacturers instructions (Polyplus-transfection). Every 3 days, the cells were supplemented with fresh medium and (where appropriate) inducer molecules. At various time points, the cells were washed with ice-cold PBS and harvested with 500 ml of FACS buffer (3% fetal bovine serum in PBS). A ow cytometry analysis was performed with a CyFlow space ow cytometer (Partec, Mnster, Germany). A 488-nm diode laser was used to detect mCitrine, a 405-nm diode laser was used to detect TagBFP and a 633-nm diode laser was used to detect the iRFP670 transfection control. In each sample, at least 30,000 cells were analysed and gated to iRFP670-positive cells. The collected data were processed using FlowJo software (TreeStar, Ashland, OR, USA). The data are representative of at least three independent experiments.

Confocal microscopy. HEK293T cells (5 104 cells per well) were seeded in eight-

well tissue culture chambers (m-Slide 8 well, Ibidi Integrated BioDiagnostics, Martinsried Mnchen, Germany). At 3050% conuence, the cells were transfected with a mixture of DNA and JetPEI transfection reagent according to the manufacturers instructions (Polyplus-transfection). Every 3 days, the cells were supplemented with fresh medium and (where appropriate) inducer molecules. Microscopic images were acquired at various time points using the Leica TCS SP5 inverted laser-scanning microscope on a Leica DMI 6000 CS module equipped with a HCX PL Fluotar L 20 , numerical aperture 0.4 (Leica Microsystems,

Wetzlar, Germany). A 514-nm laser line of a 100-mW argon laser with 25% laser power was used for mCitrine excitation, and the emitted light was detected between 520 and 580 nm. A 50-mW 405-nm diode laser was used for TagBFP excitation and the emitted light was detected between 420 and 460 nm. A 1-mW 543-nm HeNe laser was used for mCherry excitation and the emitted light was detected between 560 and 630 nm. Leica LAS AF software was used for acquisition and ImageJ software was used for image processing. The data are representative of at least three independent experiments and ve separate observations within each experiment.

References

1. Yao, G., Lee, T. J., Mori, S., Nevins, J. R. & You, L. A bistable Rb-E2F switch underlies the restriction point. Nat. Cell Biol. 10, 476482 (2008).

2. Tyson, J. J., Chen, K. C. & Novak, B. Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr. Opin. Cell Biol. 15, 221231 (2003).

3. Xiong, W. & Ferrell, J. E. A positive-feedback-based bistable memory module that governs a cell fate decision. Nature 426, 460465 (2003).

4. Ferrell, Jr J. E. The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science 280, 895898 (1998).

5. Kim, D., Rath, O., Kolch, W. & Cho, K.-H. A hidden oncogenic positive feedback loop caused by crosstalk between Wnt and ERK pathways. Oncogene 26, 45714579 (2007).

6. Eissing, T. et al. Bistability analyses of a caspase activation model for receptor-induced apoptosis. J. Biol. Chem. 279, 3689236897 (2004).

7. Monod, J. & Jacob, F. Teleonomic mechanisms in cellular metabolism, growth, and differentiation. Cold Spring Harb. Symp. Quant. Biol. 26, 389401 (1961).

8. Gardner, T. S., Cantor, C. R. & Collins, J. J. Construction of a genetic toggle switch in Escherichia coli. Nature 403, 339342 (2000).

9. Kramer, B. P. et al. An engineered epigenetic transgene switch in mammalian cells. Nat. Biotechnol. 22, 867870 (2004).

10. Elowitz, M. B. & Leibler, S. A synthetic oscillatory network of transcriptional regulators. Nature 403, 335338 (2000).

11. Atkinson, M. R., Savageau, M. a., Myers, J. T. & Ninfa, A. J. Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli. Cell 113, 597607 (2003).

12. Moon, T. S., Lou, C. B., Tamsir, A., Stanton, B. C. & Voigt, C. A. Genetic programs constructed from layered logic gates in single cells. Nature 491, 249253 (2012).

13. Garg, A., Lohmueller, J. J., Silver, P. A. & Armel, T. Z. Engineering synthetic TAL effectors with orthogonal target sites. Nucleic Acids Res. 40, 75847595 (2012).

14. Cong, L., Zhou, R., Kuo, Y. C., Cunniff, M. & Zhang, F. Comprehensive interrogation of natural TALE DNA-binding modules and transcriptional repressor domains. Nat. Commun. 3, 968 (2012).

15. Zhang, F. et al. Efcient construction of sequence-specic TAL effectors for modulating mammalian transcription. Nat. Biotechnol. 29, 149153 (2011).

16. Cermak, T. et al. Efcient design and assembly of custom TALEN and other TAL effector-based constructs for DNA targeting. Nucleic Acids Res. 39, e82 (2011).

17. Gommans, W. M., Haisma, H. J. & Rots, M. G. Engineering zincnger protein transcription factors: the therapeutic relevance of switching endogenous gene expression on or off at command. J. Mol. Biol. 354, 507519 (2005).

18. Qi, L. S. et al. Repurposing CRISPR as an RNA-guided platform for sequence-specic control of gene expression. Cell 152, 11731183 (2013).

19. Mali, P. et al. CAS9 transcriptional activators for target specicity screening and paired nickases for cooperative genome engineering. Nat. Biotechnol. 31, 833838 (2013).

20. Gilbert, L. a. et al. CRISPR-mediated modular RNA-guided regulation of transcription in eukaryotes. Cell 154, 442451 (2013).

21. Gaber, R. et al. Designable DNA-binding domains enable construction of logic circuits in mammalian cells. Nat. Chem. Biol. 10, 203208 (2014).

22. Cherry, J. L. & Adler, F. R. How to make a biological switch. J. Theor. Biol. 203, 117133 (2000).

23. Witzgall, R., OLeary, E., Leaf, A., Onaldi, D. & Bonventre, J. V. The Kruppel-associated box-A (KRAB-A) domain of zinc nger proteins mediates transcriptional repression. Proc. Natl Acad. Sci. USA 91, 45144518 (1994).

24. Kim, J. H. et al. High cleavage efciency of a 2A peptide derived from porcine teschovirus-1 in human cell lines, zebrash and mice. PLoS ONE 6, e18556 (2011).

25. Weber, W. et al. Macrolide-based transgene control in mammalian cells and mice. Nat. Biotechnol. 20, 901907 (2002).

26. Fussenegger, M. et al. Streptogramin-based gene regulation systems for mammalian cells. Nat. Biotechnol. 18, 12031208 (2000).

27. Mak, A. N.-S., Bradley, P., Cernadas, R. A., Bogdanove, A. J. & Stoddard, B. L. The crystal structure of TAL effector PthXo1 bound to its DNA target. Science 335, 716719 (2012).

28. Deng, D. et al. Structural basis for sequence-specic recognition of DNA by TAL effectors. Science 335, 720723 (2012).

29. Andrecut, M., Halley, J. D., Winkler, D. a. & Huang, S. A general model for binary cell fate decision gene circuits with degeneracy: indeterminacy and switch behavior in the absence of cooperativity. PLoS ONE 6, e19358 (2011).

30. Maca, J., Widder, S. & Sol, R. Why are cellular switches Boolean?

General conditions for multistable genetic circuits. J. Theor. Biol. 261, 126135 (2009).

31. Meckler, J. F. et al. Quantitative analysis of TALE-DNA interactions suggests polarity effects. Nucleic Acids Res. 41, 41184128 (2013).

32. Morbitzer, R., Elsaesser, J., Hausner, J. & Lahaye, T. Assembly of custom TALE-type DNA binding domains by modular cloning. Nucleic Acids Res. 39, 57905799 (2011).

33. Cong, L. et al. Multiplex genome engineering using CRISPR/Cas systems. Science 339, 819823 (2013).

34. Weber, E., Gruetzner, R., Werner, S., Engler, C. & Marillonnet, S. Assembly of designer TAL effectors by Golden Gate cloning. PLoS ONE 6, e19722 (2011).

35. Kim, Y. et al. A library of TAL effector nucleases spanning the human genome. Nat. Biotechnol. 31, 251258 (2013).

36. Munteanu, A., Constante, M., Isalan, M. & Sol, R. V. Avoiding transcription factor competition at promoter level increases the chances of obtaining oscillation. BMC Syst. Biol. 4, 66 (2010).

37. Siegal-Gaskins, D., Mejia-Guerra, M. K., Smith, G. D. & Grotewold, E. Emergence of switch-like behavior in a large family of simple biochemical networks. PLoS Comput. Biol. 7, e1002039 (2011).

38. Tan, C., Marguet, P. & You, L. Emergent bistability by a growth-modulating positive feedback circuit. Nat. Chem. Biol. 5, 842848 (2009).

39. Chen, D. & Arkin, A. P. Sequestration-based bistability enables tuning of the switching boundaries and design of a latch. Mol. Syst. Biol. 8, 620 (2012).

40. Nissim, L. & Bar-Ziv, R. H. A tunable dual-promoter integrator for targeting of cancer cells. Mol. Syst. Biol. 6, 444 (2010).

41. Gibson, D. G. et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 6, 343345 (2009).

42. Sander, J. D. et al. Targeted gene disruption in somatic zebrash cells using engineered TALENs. Nat. Biotechnol. 29, 697698 (2011).

12 NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6007 ARTICLE

43. Szulc, J., Wiznerowicz, M., Sauvain, M. O., Trono, D. & Aebischer, P.

A versatile tool for conditional gene expression and knockdown. Nat. Methods 3, 109116 (2006).

44. Sadowski, I., Ma, J., Triezenberg, S. & Ptashne, M. GAL4-VP16 is an unusually potent transcriptional activator. Nature 335, 563564 (1988).

45. Griesbeck, O., Baird, G. S., Campbell, R. E., Zacharias, D. A. & Tsien, R. Y. Reducing the environmental sensitivity of yellow uorescent protein. Mechanism and applications. J. Biol. Chem. 276, 2918829194 (2001).

46. Shcherbakova, D. M. & Verkhusha, V. V. Near-infrared uorescent proteins for multicolor in vivo imaging. Nat. Methods 10, 751754 (2013).

47. Lohmueller, J. J., Armel, T. Z. & Silver, P. A. A tunable zinc nger-based framework for Boolean logic computation in mammalian cells. Nucleic. Acids. Res. 40, 51805187 (2012).

Acknowledgements

This research was supported by grants from the Slovenian Research Agency and in part by EU structural funds to the EN-FIST Centre of Excellence. We thank Zala Lunik, Fedja Pavlovec and Maja Somrak for their contribution to the Slovenian iGEM2012 project, which initiated this work. We thank Prof Dr Martin Fussenegger for providing the constructs that encode E:KRAB, PIP:KRAB and their respective DNA-binding sites, Professor Dr Mark Ptashne for the plasmid containing the VP16 domain and Dr Oliver Griesbeck for the plasmid with the mCitrine uorescent protein.

Author contributions

R.J. designed the study. T.L., U.B., A.G., M.J., L.K., B.P., M.S., U.Z., V.F., R.G., J.L., A.O. and A.S. designed and prepared all the plasmid constructs. T.L., M.B., A.M., A.G., B.P. and U.Z. performed and analysed the confocal microscopy experiments. T.L., M.B., A.O., U.B., M.J. and L.K. performed and analysed the ow cytometry experiments. T.L., A.G., M.J., L.K. and B.P. performed and analysed the luminometry experiments. M.S. and D.V. performed all modelling and simulations of genetic switch topologies, while the other authors advised on the model design and provided experimental results and references for evaluating the parameters. T.L., M.B. and R.J. wrote the manuscript with the help of M.S. and D.V. All authors discussed the results and implications, and commented on the manuscript before submission.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Lebar, T. et al. A bistable genetic switch based on designable DNA-binding domains. Nat. Commun. 5:5007 doi: 10.1038/ncomms6007 (2014).

NATURE COMMUNICATIONS | 5:5007 | DOI: 10.1038/ncomms6007 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 13

Word count: 10602

Show less

Abstract

Translate

Details

Title

A bistable genetic switch based on designable DNA-binding domains

Author

Lebar, Tina; Bezeljak, Urban; Golob, Anja; Jerala, Miha; Kadunc, Lucija; Pirs, Bostjan; Strazar, Martin; Vucko, Dusan; Zupancic, Uros; Bencina, Mojca; Forstneric, Vida; Gaber, Rok; Lonzaric, Jan; Majerle, Andreja; Oblak, Alja; Smole, Anze; Jerala, Roman

Pages

5007

Publication year

2014

Publication date

Sep 2014

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms6007

ProQuest document ID

1566106492

A bistable genetic switch based on designable DNA-binding domains

Jump to:

Full text

Abstract

Details

Suggested sources