Team:DTU-Denmark/Modelling

From 2010.igem.org

(Difference between revisions)
 
(21 intermediate revisions not shown)
Line 1: Line 1:
-
__NOTITLE__
 
-
__NOTOC__
 
-
[[image:wiki banner_967px.png]]
 
-
 
<html>
<html>
 +
<img src="https://static.igem.org/mediawiki/2010/6/6b/Wiki_banner_967px.png"  align="center"> </img>
<head>
<head>
<title>Welcome to the DTU iGEM wiki!</title>
<title>Welcome to the DTU iGEM wiki!</title>
Line 44: Line 41:
.links tr {
.links tr {
-
  background-color:#660066;
+
  background-color:#990000;
  height:30px;
  height:30px;
}
}
Line 68: Line 65:
.links a:hover {
.links a:hover {
-
background-color: #660066;
+
background-color: #999999;
}
}
Line 105: Line 102:
  padding: .8ex;
  padding: .8ex;
}
}
 +
 +
#denav dd{
#denav dd{
Line 123: Line 122:
  letter-spacing: 1.2;
  letter-spacing: 1.2;
}
}
 +
caption {font-size:85%}
.greybox {
.greybox {
Line 133: Line 133:
</head>
</head>
<body>
<body>
-
 
<!-- Table displaying top row of links -->
<!-- Table displaying top row of links -->
Line 142: Line 141:
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Project" >The Project</a> </font></td>
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Project" >The Project</a> </font></td>
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Parts" >Parts submitted</a> </font></td>
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Parts" >Parts submitted</a> </font></td>
-
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Modelling">Modelling</a></font> </td>
+
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Results">Results</a></font> </td>
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Notebook" title="Day to day lab activity">Notebook</a>  
     <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Notebook" title="Day to day lab activity">Notebook</a>  
   <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Blog">Blog</a></font> </td>
   <td align="center" ><font face="arial" size="3"><a class="mainLinks" href="https://2010.igem.org/Team:DTU-Denmark/Blog">Blog</a></font> </td>
Line 153: Line 152:
<tr>
<tr>
<td width="163px" height="100%" valign="top">
<td width="163px" height="100%" valign="top">
 +
  <font color="#990000" face="arial" size="3">
 +
<ul type="circle">
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Background">Synthetic Biology</a></li><br>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Regulatory_sytems">Regulatory Systems</a></li>
 +
<ul><font size="2">
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Regulatory_sytems#lambda">Lambda Phage</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Regulatory_sytems#gifsy">Gifsy Phages</a></li>
 +
</font></ul>
 +
<br>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch">The Switch</a></li>
 +
<ul><font size="2">
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch#Biological_Switch">Biological Switches</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch#Bistable_Switches">Bistable Switches</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch#Design">Design of our Bi[o]stable Switch</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch#Engineering">Step-wise Engineering of the Switch</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Switch#Applications">Applications</a></li>
 +
</font></ul>
 +
<br><li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL">Synthetic Promoter Library</a>
 +
<ul><font size="2">
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL#standard">The DTU SPL Standard</a>
 +
<ul>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL#strategy">Strategy</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL#design">Primer Design</a></li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL#protocol">Protocol</a></li>
 +
</ul>
 +
</li>
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/SPL#advantages">Advantages</a></li>
 +
</ul></font>
 +
</li><br>
 +
<li ><a href="https://2010.igem.org/Team:DTU-Denmark/Modelling">Modeling</a>
 +
<ul><font size="2">
 +
<li><a href="https://2010.igem.org/Team:DTU-Denmark/Modelling#Approach">Modeling Approach</a></li>
 +
</ul></font>
 +
</li><br>
 +
</ul>
 +
</font>
</td>
</td>
<td>
<td>
-
<p align="center"><b>UNDER CONSTRUCTION</b></p>
+
 
-
<h1>Introduction</h1>
+
  <font color="#990000" face="arial" size="5">
 +
<br> 
 +
<b>Introduction</b><br><br>
 +
  </font>
 +
 
 +
 
 +
 
<p align="justify">One of the ideas behind synthetic biology is to engineer artificial systems in biological cells. This can be achieved by using well-characterized and simplified parts that can be found in nature. We will be using parts found in different systems in nature to design a bistable switch, an idea stolen from the world of electronics. A simplified version of our bistable switch is illustrated in Figure 1. </p>
<p align="justify">One of the ideas behind synthetic biology is to engineer artificial systems in biological cells. This can be achieved by using well-characterized and simplified parts that can be found in nature. We will be using parts found in different systems in nature to design a bistable switch, an idea stolen from the world of electronics. A simplified version of our bistable switch is illustrated in Figure 1. </p>
-
<p align="center"><img src="https://static.igem.org/mediawiki/2010/9/9e/DTU_modeling1.png" width="400px"  align="center"> </img></p>
+
<table class="https://static.igem.org/mediawiki/2010/9/9e/DTU_modeling1.png" align="center">
-
<p align="center">Figure 1: Simplified version of our bistable swicth </img></p><br>
+
  <caption align="bottom"><p align="justify"><b>Figure 1</b>: Simplified version of our bistable swicth</p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/9/9e/DTU_modeling1.png"  width="400px"></td></tr>
 +
</table><br>
<p align="justify">This simplified version of our switch, the repressor-repressor switch is analogous SR flip-flop of digital electronics. The SR (set-reset) flip-flop is constructed from a pair of cross-linked NAND or NOR logic gates. An SR flip-flop circuit built with NOR logic gates is illustrated in Figure 2.</p>
<p align="justify">This simplified version of our switch, the repressor-repressor switch is analogous SR flip-flop of digital electronics. The SR (set-reset) flip-flop is constructed from a pair of cross-linked NAND or NOR logic gates. An SR flip-flop circuit built with NOR logic gates is illustrated in Figure 2.</p>
-
<p align="center"><img src="https://static.igem.org/mediawiki/2010/5/53/SRflipflop.png" width="200px"  align="center"> </img></p>
+
<table class="https://static.igem.org/mediawiki/2010/5/53/SRflipflop.png" align="center">
-
<p align="center">Figure 2: SR flip-flop </img></p><br>
+
  <caption align="bottom"><p align="justify"><b>Figure 2</b>: SR flip-flop</p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/5/53/SRflipflop.png"  width="200px"></td></tr>
 +
</table><br>
 +
 
<p align="justify">The circuit, built with the pair of cross-linked NOR logic gates, functions such that output from a gate will only be true if no input is given as demonstrated by the truth table in Table 1. The biological interpretation of this logical operation shows is that the repression by a protein can be prevented either by an inducer (or anti-repressor protein as in our case) (Input 1/2) or by repression of its synthesis by another repressor-protein (Output 2/1). The elements in the circuit do not possess memory in them, but by cross-linking the elements, the circuit is capable of “remembering” and holding its current state. This means that when neither Input 1 nor Input 2 (no inducers) is present, the circuit is bistable with either Output 2 or Output 1 being active (only one of the repressor-proteins are expressed), respectively.</p>
<p align="justify">The circuit, built with the pair of cross-linked NOR logic gates, functions such that output from a gate will only be true if no input is given as demonstrated by the truth table in Table 1. The biological interpretation of this logical operation shows is that the repression by a protein can be prevented either by an inducer (or anti-repressor protein as in our case) (Input 1/2) or by repression of its synthesis by another repressor-protein (Output 2/1). The elements in the circuit do not possess memory in them, but by cross-linking the elements, the circuit is capable of “remembering” and holding its current state. This means that when neither Input 1 nor Input 2 (no inducers) is present, the circuit is bistable with either Output 2 or Output 1 being active (only one of the repressor-proteins are expressed), respectively.</p>
 +
<table cellpadding="2" cellspacing="0" border="1px" align="center" >
 +
<font size="2">
 +
<caption><p align="justify"><b>Table 1</b></p></caption>
 +
<thead>
 +
<td align="center"><b>Input 1<br>(Anti-Repressor 1)</td><td align="center"><b>Input 2<br>(Anti-Repressor 2)</td><td align="center"><b>Output 1<br>(Reporter 1)</td><td align="center"><b>Output 2<br>(Reporter 2)</td>
 +
</thead>
 +
<tr>
 +
<td align="center">0</td><td align="center">0</td><td align="center">No change</td><td align="center">No change</td>
 +
</tr>
 +
<tr>
 +
<td align="center">0</td><td align="center">1</td><td align="center">1</td><td align="center">0</td>
 +
</tr>
 +
<tr>
 +
<td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">1</td>
 +
</tr>
 +
<tr>
 +
<td align="center">1</td><td align="center">1</td><td align="center">Invalid</td><td align="center">Invalid</td>
 +
</tr>
 +
</font>
 +
</table>
 +
 +
 +
  <font color="#990000" face="arial" size="5">
 +
<br> 
 +
<b>Modeling approach</b><br><br>
 +
  </font>
 +
 +
 +
<p align="justify">The simplified version of our <a href="https://2010.igem.org/Team:DTU-Denmark/Switch#final" target="_blank">Bi[o]stable switch</a> illustrated in Figure 1 was used in the modeling. A system must possess two or more stable states for it to function as a switch. The essential components of the switch are the two repressors and two repressible promoters. Each promoter is inhibited by the repressor expressed from the opposing promoter. The aim of our project is to design a robust bistable switch that exhibits bi-stability over a wide-range of parameters and that the two states are tolerant of the fluctuations inherent in gene expression. This will allow the latest induced state to remain constant over cell-generations. The bi-stability will arise as a result of the mutually inhibitory arrangement. As denoted in Table 1, the absence of inducers will leave the cell in its current state, where either promoter 1 expresses repressor 1 or promoter 2 expresses repressor 2. The expression of the repressor will prevent expression from the opposing promoter, preventing the spontaneous change of state.</p>
 +
<p align="justify">The change of state is accomplished by transiently introducing an inducer (or an anti-repressor protein as in our case) of the currently active repressor. Maximal expression from the opposing promoter and establishment of the alternative stable state will ensue (Gardner, T.S. et al, 2000)(Cherry, J., 2000).</p>
 +
 +
<h3>Michaelis-Menten Function vs Hill Function</h3>
 +
<p align="justify">In our repressor-repressor switch, the simplest switch is one where each repressor binds to its operator by mass action and the rate of transcription is proportional to the amount of unbound operator. A Michaelis-Menten equation can approximate the fraction of bound operator when the molar quantity of the repressor is greater than the operator. The rate of production of repressor protein <i>x</i> is given by:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/1/18/DTU_modeling_form1.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/1/18/DTU_modeling_form1.png" width="110px"></td></tr>
 +
</table><br>
 +
<p align="justify">Where <i>k<sub>1</sub></i> is the maximal rate of protein production, <i>K<sub>D</sub></i> is the equilibrium constant for dissociation of repressor from its operator site, <i>y</i> is the concentration of repressor protein y</p>
 +
<p align="justify">By expressing K<sub><i>d</i></sub> in terms of repressor y:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/c/c7/DTU_modeling_form2.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/c/c7/DTU_modeling_form2.png" width="100px"></td></tr>
 +
</table><br>
 +
<p align="justify">And conversely, the rate of production for repressor protein 2 is given by:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/1/12/DTU_modeling_form3.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/1/12/DTU_modeling_form3.png" width="100px"></td></tr>
 +
</table><br>
 +
<p align="justify">Where similarly, <i>k<sub>2</sub></i> is the maximal rate of protein production of repressor y, <i>K<sub>D</sub></i> is the equilibrium constant for dissociation of repressor x from its operator site, <i>x</i> is the concentration of repressor protein x</p>
 +
<p align="justify">This results in the following repressor-repressor system:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" width="120px"></td></tr>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/7/77/DTU_modeling_form5.png" width="120px"></td></tr>
 +
</table><br>
 +
 +
 +
<p align="justify">The resulting nullclines are defined by:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/0/00/DTU_modeling_form6.png" width="120px"></td></tr>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/3/3e/DTU_modeling_form7.png" width="120px"></td></tr>
 +
</table><br>
 +
<p align="justify">At equilibrium:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/a/a8/DTU_modeling_form8.png" width="100px"></td></tr>
 +
</table><br>
 +
 +
<p align="justify">This results in a linear fractional transformation with a negative second derivative everywhere and is unable to have multiple intersections with the diagonal. As the multiple intersections of the second derivative with the diagonal is a requirement for a switch, the repression described can therefore not support a switch. However, the inclusion of the depletion term and more than one operator site per cell makes a switch possible. The Hill function that takes cooperativity of binding when describing the level of repression is:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/4/4f/DTU_modeling_form9.png" width="100px" ></td></tr>
 +
</table><br>
 +
<p align="justify">The behavior and conditions for bistability can be elucidated with the following dimensionless model by replacing the Hill equation into the repressor-repressor system, where the depletion term is also taken into account:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/c/cd/DTU_modeling_form10.png" width="120px"></td></tr>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/8/8d/DTU_11.png"  width="120px"></td></tr>
 +
</table><br>
 +
<p align="justify">The parameters k<sub>1</sub> and k<sub>2</sub> describes the net effect of RNA polymerase binding, open-complex formation, transcription elongation, transcription termination, repressor binding, ribosome binding and translation. The cooperativity, described by n, can be due to the multimerization of the repressor proteins as well as the cooperative binding of repressor multimers to multiple binding sites within the operator. The nullclines of the system are:</p>
 +
<table class="https://static.igem.org/mediawiki/2010/3/3f/DTU_modeling_form4.png" align="center">
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/3/39/DTU_modeling_form12.png" width="120px"></td></tr>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/d/dd/DTU_13.png" width="120px"></td></tr>
 +
</table><br>
 +
 +
<table class="https://static.igem.org/mediawiki/2010/6/62/DTU_modeling_fi3.png" align="center">
 +
<caption align="bottom"><p align="justify"><b>Figure 3</b>: Shows the two stable states and the unstable state of the switch.</p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/6/62/DTU_modeling_fi3.png"  width="200px"></td></tr>
 +
</table><br>
 +
<p align="justify">The nullclines intersect at three points, resulting in the two steady states and one unstable state shown in Figure 3. From this it is clear that three intersection points occur due to the sigmoid shape of the graph, which arises when n > 1 showing that cooperativity of binding is necessary. Another key that can be extracted is the importance of similar rates of synthesis of both repressors. If the rates are unbalanced, the nullclines will only intersect once producing only a single stable state. An increase in the rate of repressor synthesis will result in a more robust switch (Gardner, T.S. et al, 2000)(Cherry, J., 2000). Figure 4 illustrates the functions of a good switch.</p>
 +
<table class="https://static.igem.org/mediawiki/2010/1/14/DTU_modeling_goodswitch.png" align="center">
 +
<caption align="bottom"><p align="center"><b>Figure 4</b></p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/1/14/DTU_modeling_goodswitch.png"  width="200px"></td></tr>
 +
</table><br>
 +
 +
<h3>Robustness of the switch</h3>
 +
<p align="justify">Real switches are prone to stochastic fluctuations or noise, this is a consequence of the lack of separation of the states shown in Figure 5. Poor transversality of the nullclines, illustrated in Figure 6, due to real world fluctuations in parameters can lead to spontaneous change of state and thereby also loss of the information stored by the switch. Systems with highly transverse nullclines are able to withstand movements of the nullclines without spontaneous change of state. Stochastic fluctuations in protein levels occur due to the fact that the number of protein and mRNA molecules are finite. This means that it is more beneficial for a switch to be based on high concentrations of low-affinity repressors compared to a low concentration of high-affinity repressors. The cost of additional protein production is balanced out by the reduction of noise in the system (Gardner, T.S. et al, 2000)(Cherry, J., 2000).</p>
 +
<table class="https://static.igem.org/mediawiki/2010/2/2b/DTU_modeling_poorsep.png" align="center">
 +
<caption align="bottom"><p align="center"><b>Figure 5</b></p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/2/2b/DTU_modeling_poorsep.png"  width="200px"></td></tr>
 +
</table><br>
 +
<table class="https://static.igem.org/mediawiki/2010/7/70/DTU_modeling_poortrans.png" align="center">
 +
<caption align="bottom"><p align="center"><b>Figure 6</b></p></caption>
 +
<tr><td><img src="https://static.igem.org/mediawiki/2010/7/70/DTU_modeling_poortrans.png"  width="200px"></td></tr>
 +
</table><br>
 +
 +
<h3>References</h3>
 +
<ul>
 +
<li>Gardner, T.S., Cantor, C.R. & Collins, J.J. Construction of a genetic toggle switch in Escherichia coli. Nature 403, 339-342 (2000).
 +
<li>Cherry, J. How to make a Biological Switch. Journal of Theoretical Biology 203, 117-133 (2000).
 +
</ul>
</td>
</td>
<td width="163px" height="100%" valign="top">
<td width="163px" height="100%" valign="top">

Latest revision as of 03:30, 28 October 2010

Welcome to the DTU iGEM wiki!



Introduction

One of the ideas behind synthetic biology is to engineer artificial systems in biological cells. This can be achieved by using well-characterized and simplified parts that can be found in nature. We will be using parts found in different systems in nature to design a bistable switch, an idea stolen from the world of electronics. A simplified version of our bistable switch is illustrated in Figure 1.

Figure 1: Simplified version of our bistable swicth


This simplified version of our switch, the repressor-repressor switch is analogous SR flip-flop of digital electronics. The SR (set-reset) flip-flop is constructed from a pair of cross-linked NAND or NOR logic gates. An SR flip-flop circuit built with NOR logic gates is illustrated in Figure 2.

Figure 2: SR flip-flop


The circuit, built with the pair of cross-linked NOR logic gates, functions such that output from a gate will only be true if no input is given as demonstrated by the truth table in Table 1. The biological interpretation of this logical operation shows is that the repression by a protein can be prevented either by an inducer (or anti-repressor protein as in our case) (Input 1/2) or by repression of its synthesis by another repressor-protein (Output 2/1). The elements in the circuit do not possess memory in them, but by cross-linking the elements, the circuit is capable of “remembering” and holding its current state. This means that when neither Input 1 nor Input 2 (no inducers) is present, the circuit is bistable with either Output 2 or Output 1 being active (only one of the repressor-proteins are expressed), respectively.

Table 1

Input 1
(Anti-Repressor 1)
Input 2
(Anti-Repressor 2)
Output 1
(Reporter 1)
Output 2
(Reporter 2)
00No changeNo change
0110
1001
11InvalidInvalid

Modeling approach

The simplified version of our Bi[o]stable switch illustrated in Figure 1 was used in the modeling. A system must possess two or more stable states for it to function as a switch. The essential components of the switch are the two repressors and two repressible promoters. Each promoter is inhibited by the repressor expressed from the opposing promoter. The aim of our project is to design a robust bistable switch that exhibits bi-stability over a wide-range of parameters and that the two states are tolerant of the fluctuations inherent in gene expression. This will allow the latest induced state to remain constant over cell-generations. The bi-stability will arise as a result of the mutually inhibitory arrangement. As denoted in Table 1, the absence of inducers will leave the cell in its current state, where either promoter 1 expresses repressor 1 or promoter 2 expresses repressor 2. The expression of the repressor will prevent expression from the opposing promoter, preventing the spontaneous change of state.

The change of state is accomplished by transiently introducing an inducer (or an anti-repressor protein as in our case) of the currently active repressor. Maximal expression from the opposing promoter and establishment of the alternative stable state will ensue (Gardner, T.S. et al, 2000)(Cherry, J., 2000).

Michaelis-Menten Function vs Hill Function

In our repressor-repressor switch, the simplest switch is one where each repressor binds to its operator by mass action and the rate of transcription is proportional to the amount of unbound operator. A Michaelis-Menten equation can approximate the fraction of bound operator when the molar quantity of the repressor is greater than the operator. The rate of production of repressor protein x is given by:


Where k1 is the maximal rate of protein production, KD is the equilibrium constant for dissociation of repressor from its operator site, y is the concentration of repressor protein y

By expressing Kd in terms of repressor y:


And conversely, the rate of production for repressor protein 2 is given by:


Where similarly, k2 is the maximal rate of protein production of repressor y, KD is the equilibrium constant for dissociation of repressor x from its operator site, x is the concentration of repressor protein x

This results in the following repressor-repressor system:


The resulting nullclines are defined by:


At equilibrium:


This results in a linear fractional transformation with a negative second derivative everywhere and is unable to have multiple intersections with the diagonal. As the multiple intersections of the second derivative with the diagonal is a requirement for a switch, the repression described can therefore not support a switch. However, the inclusion of the depletion term and more than one operator site per cell makes a switch possible. The Hill function that takes cooperativity of binding when describing the level of repression is:


The behavior and conditions for bistability can be elucidated with the following dimensionless model by replacing the Hill equation into the repressor-repressor system, where the depletion term is also taken into account:


The parameters k1 and k2 describes the net effect of RNA polymerase binding, open-complex formation, transcription elongation, transcription termination, repressor binding, ribosome binding and translation. The cooperativity, described by n, can be due to the multimerization of the repressor proteins as well as the cooperative binding of repressor multimers to multiple binding sites within the operator. The nullclines of the system are:


Figure 3: Shows the two stable states and the unstable state of the switch.


The nullclines intersect at three points, resulting in the two steady states and one unstable state shown in Figure 3. From this it is clear that three intersection points occur due to the sigmoid shape of the graph, which arises when n > 1 showing that cooperativity of binding is necessary. Another key that can be extracted is the importance of similar rates of synthesis of both repressors. If the rates are unbalanced, the nullclines will only intersect once producing only a single stable state. An increase in the rate of repressor synthesis will result in a more robust switch (Gardner, T.S. et al, 2000)(Cherry, J., 2000). Figure 4 illustrates the functions of a good switch.

Figure 4


Robustness of the switch

Real switches are prone to stochastic fluctuations or noise, this is a consequence of the lack of separation of the states shown in Figure 5. Poor transversality of the nullclines, illustrated in Figure 6, due to real world fluctuations in parameters can lead to spontaneous change of state and thereby also loss of the information stored by the switch. Systems with highly transverse nullclines are able to withstand movements of the nullclines without spontaneous change of state. Stochastic fluctuations in protein levels occur due to the fact that the number of protein and mRNA molecules are finite. This means that it is more beneficial for a switch to be based on high concentrations of low-affinity repressors compared to a low concentration of high-affinity repressors. The cost of additional protein production is balanced out by the reduction of noise in the system (Gardner, T.S. et al, 2000)(Cherry, J., 2000).

Figure 5


Figure 6


References

  • Gardner, T.S., Cantor, C.R. & Collins, J.J. Construction of a genetic toggle switch in Escherichia coli. Nature 403, 339-342 (2000).
  • Cherry, J. How to make a Biological Switch. Journal of Theoretical Biology 203, 117-133 (2000).