Team:Slovenia/ModelingMutualRepressorSwitchStochastic

From 2012.igem.org

(Difference between revisions)
 
(14 intermediate revisions not shown)
Line 2: Line 2:
<head>
<head>
<meta http-equiv="X-UA-Compatible" content="IE=edge" />
<meta http-equiv="X-UA-Compatible" content="IE=edge" />
 +
 +
<div style="position:fixed; bottom:45px; right:30px; width:100px; height:66px; background-color:transparent;">
 +
<a href="#topofthepage">
 +
<table style="background-color:transparent;" onclick="window.location = '#topofthepage'" class="invisible" style="height:100%;">
 +
<tr class="invisible" style="background-color:transparent;">
 +
<td class="invisible" style="background-color:transparent;" valign="center">
 +
<img width="100px" src ="https://static.igem.org/mediawiki/2012/1/14/Svn12_hp_btt.png"></img>
 +
</td></tr></table>
 +
</a>
 +
</div>
 +
 +
<style type="text/css">
<style type="text/css">
Line 13: Line 25:
-
p {line-height:1.5em; margin:0 0 15px; text-align:left;}
+
p {line-height:1.5em; margin:0 0 15px; text-align:justify;}
h2 {font-size:1.8em; font-weight:400; margin:0 0 12px;}
h2 {font-size:1.8em; font-weight:400; margin:0 0 12px;}
Line 182: Line 194:
 +
p.inliner{text-align:justify; }
 +
img.inliner{width:100%; border-width:0px;}
 +
table.inliner{border-width:0px; float:right; width:60%; margin-bottom:15px; margin-left:15 px; margin-right:15 px; padding-left:15px; }
 +
td.inliner, th.normal { border-width:0px;}
 +
thead.inliner{}
 +
tbody .inliner{}
 +
 +
/* CSS navigation menu (blue) */
/* CSS navigation menu (blue) */
/* CSS navigation menu (blue) */
#cssmenu{ height:37px; display:block; padding:0; margin:0;  border:1px solid; border-radius:0px; }  
#cssmenu{ height:37px; display:block; padding:0; margin:0;  border:1px solid; border-radius:0px; }  
Line 201: Line 221:
#cssmenu ul li > ul li{display:block; list-style:inside none; padding:0; margin:0; position:relative;}  
#cssmenu ul li > ul li{display:block; list-style:inside none; padding:0; margin:0; position:relative;}  
#cssmenu ul li > ul li a{ outline:none; display:block; position:relative; margin:0; padding:8px 20px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }  
#cssmenu ul li > ul li a{ outline:none; display:block; position:relative; margin:0; padding:8px 20px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }  
 +
#cssmenu ul li > ul li a table tr td span{ outline:none; display:block; position:relative; margin:0; padding:0px 0px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }
#cssmenu, #cssmenu > ul > li > ul > li a:hover
#cssmenu, #cssmenu > ul > li > ul > li a:hover
{ background:#043A6B;  
{ background:#043A6B;  
Line 216: Line 237:
#cssmenu > ul > li > a:hover{background:#0C5DA5;}  
#cssmenu > ul > li > a:hover{background:#0C5DA5;}  
/* end CSS navigation menu (blue) */
/* end CSS navigation menu (blue) */
 +
 +
/*new table start*/
 +
table.newtable {background-color:transparent;}
 +
td.newtable, th.newtable {background-color:transparent;}
 +
thead.newtable{ }
 +
tbody .newtable{}
 +
/*new table start*/
Line 272: Line 300:
</style>
</style>
-
 
<!-- google analytics -->
<!-- google analytics -->
Line 292: Line 319:
<body>
<body>
 +
<a name="topofthepage" style="background-color:transparent;"></a>
 +
<div id="banner">
<div id="banner">
-
<a style="position:absolute; top:0px; left:490px;" href="https://2012.igem.org/Main_Page"><b>iGEM 2012</b></a>
 
</div>
</div>
Line 300: Line 328:
<!--<div id="container">-->
<!--<div id="container">-->
<div id="mainmenu">
<div id="mainmenu">
-
 
<!-- start main menu -->
<!-- start main menu -->
<div id='cssmenu'>
<div id='cssmenu'>
Line 317: Line 344:
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchDesignedTALregulators'><span>Designed TAL regulators</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchDesignedTALregulators'><span>Designed TAL regulators</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>  
+
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch';" class="newtable"><tr class="newtable"><td class="newtable"><span>Positive feedback loop switch</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
 +
    <li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchControls'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TheSwitchControls';" class="newtable"><tr class="newtable"><td class="newtable"><span>Controls</span></td><td class="newtable"><img style="margin-right:-81px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
  </ul>
  </ul>
</li>
</li>
Line 325: Line 353:
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsEscapeTag'><span>Escape tag</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsEscapeTag'><span>Escape tag</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsTermination'><span>Termination</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsTermination'><span>Termination</span></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation'><span>Microcapsule degradation</span></a></li>  
+
    <li><a href="https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation"><table  onclick="window.location = 'https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation';" class="newtable"><tr class="newtable"><td class="newtable"><span>Microcapsule degradation</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
  </ul>
  </ul>
</li>
</li>
Line 333: Line 361:
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationHepatitisC'><span>Hepatitis C</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationHepatitisC'><span>Hepatitis C</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationIschaemicHeartDisease'><span>Ischaemic heart disease</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationIschaemicHeartDisease'><span>Ischaemic heart disease</span></a></li>  
 +
    <li><a href='https://2012.igem.org/Team:Slovenia/ImplementationImpact'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ImplementationImpact';" class="newtable"><tr class="newtable"><td class="newtable"><span>Impact</span></td><td class="newtable"><img style="margin-right:-86px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
 
 
  </ul>
  </ul>
Line 340: Line 369:
  <ul>
  <ul>
<li><a href='https://2012.igem.org/Team:Slovenia/Modeling'><span>Overview</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Modeling'><span>Overview</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPK'><span>Pharmacokinetics</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/ModelingPK'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingPK';" class="newtable"><tr class="newtable"><td class="newtable"><span>Pharmacokinetics</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMethods'><span>Modeling methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMethods'><span>Modeling methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel'><span>Quantitative and stability model</span></a></li>  
+
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel';" class="newtable"><tr class="newtable"><td class="newtable"><span>Experimental model</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations'><span>Interactive simulations</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations';" class="newtable"><tr class="newtable"><td class="newtable"><span>Interactive simulations</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
  </ul>
  </ul>
</li>
</li>
Line 356: Line 385:
  <ul>
  <ul>
<li><a href='https://2012.igem.org/Team:Slovenia/Notebook'><span>Experimental methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Notebook'><span>Experimental methods</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLablog'><span>Lablog</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/NotebookLablog'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/NotebookLablog';" class="newtable"><tr class="newtable"><td class="newtable"><span>Lablog</span></td><td class="newtable"><img style="margin-right:-90px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLabSafety'><span>Lab safety</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLabSafety'><span>Lab safety</span></a></li>  
  </ul>
  </ul>
Line 379: Line 408:
<li><a href='https://2012.igem.org/Team:Slovenia/Team'><span>Team members</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Team'><span>Team members</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamAttributions'><span>Attributions</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamAttributions'><span>Attributions</span></a></li>
 +
<li><a href='https://2012.igem.org/Team:Slovenia/TeamCollaborations'><table  onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TeamCollaborations';" class="newtable"><tr class="newtable"><td class="newtable"><span>Collaborations</span></td><td class="newtable"><img style="margin-right:-20px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamGallery'><span>Gallery</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamGallery'><span>Gallery</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamSponsors'><span>Sponsors</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamSponsors'><span>Sponsors</span></a></li>  
Line 391: Line 421:
<div id="main">
<div id="main">
<br/>
<br/>
 +
 +
 +
<!-- dummy link na bannerju -->
 +
<a href="https://2012.igem.org/Main_Page">
 +
<div id="dummy" style="background-color:transparent; position:absolute; left:870px; top:25px; width:115px; height:80px; z-index:100; opacity:0.0;">
 +
</div>
 +
</a>
 +
<h1>Modeling - mutual repressor switch</h1>
<h1>Modeling - mutual repressor switch</h1>
<p>
<p>
Line 416: Line 454:
<p>
<p>
Stochastic simulation revealed that the mutual repressor switch exhibited bistability only as long as there was no leaky expression present. Otherwise, cooperativity higher than 1 was required. Higher cooperativity improved leaky expression tolerance, but only to a certain threshold. Higher leaky expression resulted in lower stable-state levels.
Stochastic simulation revealed that the mutual repressor switch exhibited bistability only as long as there was no leaky expression present. Otherwise, cooperativity higher than 1 was required. Higher cooperativity improved leaky expression tolerance, but only to a certain threshold. Higher leaky expression resulted in lower stable-state levels.
-
</p>
 
-
 
-
<p>
 
-
Bistable behavior of the switch became inconsistent when a delay between transcription and translation was taken into account.
 
</p>
</p>
</td>
</td>
Line 435: Line 469:
<p><img src="https://static.igem.org/mediawiki/2012/d/da/Svn12_mrs_stoh4.PNG" /></p>
<p><img src="https://static.igem.org/mediawiki/2012/d/da/Svn12_mrs_stoh4.PNG" /></p>
<p><img src="https://static.igem.org/mediawiki/2012/1/19/Svn12_mrs_stoh5.PNG" /></p>
<p><img src="https://static.igem.org/mediawiki/2012/1/19/Svn12_mrs_stoh5.PNG" /></p>
-
<p><img src=https://static.igem.org/mediawiki/2012/7/73/Svn12_mrs_stoh6.PNG"" /></p>
+
<p><img src="https://static.igem.org/mediawiki/2012/7/73/Svn12_mrs_stoh6.PNG" /></p>
<p><img src="https://static.igem.org/mediawiki/2012/1/10/Svn12_mrs_stoh7.PNG" /></p>
<p><img src="https://static.igem.org/mediawiki/2012/1/10/Svn12_mrs_stoh7.PNG" /></p>
 +
 +
<p>
 +
Here:
 +
<ul style="margin-left:30px;">
 +
    <li>Pro1 is construct 1 promoter (i.e. promoter 1 - constitutive);</li>
 +
    <li>Pro2 is construct 2 promoter (i.e. promoter 2 - constitutive);</li>
 +
    <li>Pro3 s construct 3 promoter (i.e. promoter 3 - constitutive);</li> 
 +
    <li>Pro4 is construct 4 promoter (i.e. promoter 4 - constitutive);</li>       
 +
    <li>Pro5 is construct 5 promoter (i.e. promoter 5 - constitutive);</li>
 +
    <li>Ind1 Is inducer 1 (i.e. signal 1), used to induce stable state 1;</li>
 +
    <li>Ind2 is inducer 2 (i.e. signal 2), used to induce stable state 2.</li>
 +
</ul>
 +
</p>
 +
 +
<p>
 +
[multi] is a variable, equal to a degree of oligomerization and is used to model cooperativity. Normal gene expression means constitutive gene expression, occuring when no repressor is bound.
 +
</p>
 +
 +
 +
<p>
 +
When inducer 1 is present, it binds PIP:KRAB and forms a complex, denoted as Ind1.PIP:KRAB. This complex is referred to as inactive PIP:KRAB, meaning PIP:KRAB that cannot bind to the promoter 3. When inducer 2 is present, it binds E:KRAB and forms a complex, denoted as Ind2.E:KRAB. This complex is referred to as inactive E:KRAB, meaning E:KRAB that cannot bind to the promoter 4.
 +
</p>
Line 442: Line 498:
<h2><a name="results">Simulation results</a></h2>
<h2><a name="results">Simulation results</a></h2>
 +
<p>
 +
Simulation results are shown as concentrations of reporter for different states (BFP, mCitrine) indicating one of the two states as a function of time. No specific units were used, hence no absolute interpretation of the results' values in terms of units is in place. Switching between states was achieved using two signals (inducers) that were introduced into the system at appropriate time. The signals were modeled as a high number of molecules that triggered state induction and were removed at an appropriate time.
 +
</p>
 +
<p>
 +
Cooperativity was modeled as a degree of oligomerization, with repressor oligomers binding to a promoter.
 +
</p>
 +
<p>
 +
Protein production:degradation rate ratio was set to 10. Promoter unbinding (e.g. bound transcription factor unbinding from the target binding site) reactions were assumed to be 100 times slower than promoter binding reactions, which is in the range of typical transcription factors. Initial concentrations of proteins were 0 for all simulations.
 +
</p>
 +
<p>Detailed parameter values for each simulation can be found in the corresponding simulation files that can be found <a href="https://2012.igem.org/Team:Slovenia/SourceCode">here</a>.</p>
 +
 +
<p>
 +
The purpose of the first simulation was to show that the mutual repressor switch can exhibit bistability in a stochastic environment. No leaky expression was assumed, and no cooperativity (i.e. no TAL reprssor oligomerization). The following state-switching scenario was used in all stochastic simulations:
 +
<ul style="margin-left:30px;">
 +
    <li>at time = 500, signal 2 was introduced, inducing stable state 2 (high mCitrine state); the signal gradually degraded, but stable state was preserved even after degradation;</li>
 +
    <li>at time = 2500, signal 1 was introduced to induce stable state 1 (high BFP state); the switch remained in a stable state after signal degradation;</li>
 +
    <li>at time = 4500, signal 2 was introduced again, switching the system back to stable state 2;at time = 6500, signal 1 was introduced again, switching the system to stable state 1.</li>
 +
    <li>at time = 6500, signal 1 was introduced again, switching the system to stable state 1.</li>
 +
</ul>
 +
</p>
 +
 +
<p>
 +
As shown in Figure 1, bistability was exhibited for this scenario, despite no cooperativity. However, as soon as sufficient leaky expression (e.g. 0.03 compared to protein production rate of 1) was introduced for each gene, bistability was lost (Figure 2).</p>
 +
 +
<!-- figure 1 -->
 +
<table class="invisible" style="width:90%;">
 +
<tbody  class="invisible">
 +
<tr class="invisible">
 +
<td class="invisible">
 +
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/d/d6/Svn12_MutualRepressorSwitch_test1_mtl.png"/>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<table class="invisible" style="width:90%; text-align:center;">
 +
<tbody  class="invisible">
 +
<tr class="normal"><td class="invisible">
 +
<b>Figure 1.</b> Mutual repressor switch exhibiting transitions between stable states. No cooperativity and zero leaky expression were assumed here.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
 +
 +
<!-- figure 2 -->
 +
<table class="invisible" style="width:90%;">
 +
<tbody  class="invisible">
 +
<tr class="invisible">
 +
<td class="invisible">
 +
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/8/80/Svn12_MutualRepressorSwitch_test2_mtl.png"/>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<table class="invisible" style="width:90%; text-align:center;">
 +
<tbody  class="invisible">
 +
<tr class="normal"><td class="invisible">
 +
<b>Figure 2.</b> Mutual repressor switch exhibited no bistability when no cooperativity was used together with a relatively low leaky expression of 0.03.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
 +
 +
<p>
 +
For the switch to exhibit bistability in the presence of leaky expression, higher cooperativity values were required. Figure 3 shows bistable behavior of the switch for the cooperativity of 3 in the presence of leaky expression. Higher cooperativity led to higher leaky expression tolerance, but only to a certain leaky expression threshold, depending on other parameter values, such as production and degradation rates. Higher leaky expression also lowered the reached stable-state levels. The threshold for maximal cooperativity allowing bistability depended, just like leaky expression threshold, on other parameter values, such as production and degradation rates. Figure 4 shows that for very high cooperativity of 10, no bistability was observed at leaky expression rate of 0.05.
 +
</p>
 +
 +
 +
 +
 +
<!-- figure 3 -->
 +
<table class="invisible" style="width:90%;">
 +
<tbody  class="invisible">
 +
<tr class="invisible">
 +
<td class="invisible">
 +
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/7/77/Svn12_MutualRepressorSwitch_test3_mtl.png"/>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<table class="invisible" style="width:90%; text-align:center;">
 +
<tbody  class="invisible">
 +
<tr class="normal"><td class="invisible">
 +
<b>Figure 3.</b> Bistability was obtained for cooperativity equal to 3 and leaky expression equal to 0.05.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
 +
 +
<!-- figure 4 -->
 +
<table class="invisible" style="width:90%;">
 +
<tbody  class="invisible">
 +
<tr class="invisible">
 +
<td class="invisible">
 +
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/c/ce/Svn12_MutualRepressorSwitch_test4_mtl.png"/>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<table class="invisible" style="width:90%; text-align:center;">
 +
<tbody  class="invisible">
 +
<tr class="normal"><td class="invisible">
 +
<b>Figure 4.</b> Bistability was lost for cooperativity equal to 10 and leaky expression equal to 0.05.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->

Latest revision as of 21:16, 26 October 2012


Modeling - mutual repressor switch

  1. Deterministic model
  2. Stochastic model
  3. C#Sim model

Stochastic model of the mutual repressor switch

Stochastic simulation revealed that the mutual repressor switch exhibited bistability only as long as there was no leaky expression present. Otherwise, cooperativity higher than 1 was required. Higher cooperativity improved leaky expression tolerance, but only to a certain threshold. Higher leaky expression resulted in lower stable-state levels.

The model

The basis for the stochastic simulation of the mutual repressor switch was the following set of reactions that describe the dynamics of the switch:

Here:

  • Pro1 is construct 1 promoter (i.e. promoter 1 - constitutive);
  • Pro2 is construct 2 promoter (i.e. promoter 2 - constitutive);
  • Pro3 s construct 3 promoter (i.e. promoter 3 - constitutive);
  • Pro4 is construct 4 promoter (i.e. promoter 4 - constitutive);
  • Pro5 is construct 5 promoter (i.e. promoter 5 - constitutive);
  • Ind1 Is inducer 1 (i.e. signal 1), used to induce stable state 1;
  • Ind2 is inducer 2 (i.e. signal 2), used to induce stable state 2.

[multi] is a variable, equal to a degree of oligomerization and is used to model cooperativity. Normal gene expression means constitutive gene expression, occuring when no repressor is bound.

When inducer 1 is present, it binds PIP:KRAB and forms a complex, denoted as Ind1.PIP:KRAB. This complex is referred to as inactive PIP:KRAB, meaning PIP:KRAB that cannot bind to the promoter 3. When inducer 2 is present, it binds E:KRAB and forms a complex, denoted as Ind2.E:KRAB. This complex is referred to as inactive E:KRAB, meaning E:KRAB that cannot bind to the promoter 4.

Simulation results

Simulation results are shown as concentrations of reporter for different states (BFP, mCitrine) indicating one of the two states as a function of time. No specific units were used, hence no absolute interpretation of the results' values in terms of units is in place. Switching between states was achieved using two signals (inducers) that were introduced into the system at appropriate time. The signals were modeled as a high number of molecules that triggered state induction and were removed at an appropriate time.

Cooperativity was modeled as a degree of oligomerization, with repressor oligomers binding to a promoter.

Protein production:degradation rate ratio was set to 10. Promoter unbinding (e.g. bound transcription factor unbinding from the target binding site) reactions were assumed to be 100 times slower than promoter binding reactions, which is in the range of typical transcription factors. Initial concentrations of proteins were 0 for all simulations.

Detailed parameter values for each simulation can be found in the corresponding simulation files that can be found here.

The purpose of the first simulation was to show that the mutual repressor switch can exhibit bistability in a stochastic environment. No leaky expression was assumed, and no cooperativity (i.e. no TAL reprssor oligomerization). The following state-switching scenario was used in all stochastic simulations:

  • at time = 500, signal 2 was introduced, inducing stable state 2 (high mCitrine state); the signal gradually degraded, but stable state was preserved even after degradation;
  • at time = 2500, signal 1 was introduced to induce stable state 1 (high BFP state); the switch remained in a stable state after signal degradation;
  • at time = 4500, signal 2 was introduced again, switching the system back to stable state 2;at time = 6500, signal 1 was introduced again, switching the system to stable state 1.
  • at time = 6500, signal 1 was introduced again, switching the system to stable state 1.

As shown in Figure 1, bistability was exhibited for this scenario, despite no cooperativity. However, as soon as sufficient leaky expression (e.g. 0.03 compared to protein production rate of 1) was introduced for each gene, bistability was lost (Figure 2).

For the switch to exhibit bistability in the presence of leaky expression, higher cooperativity values were required. Figure 3 shows bistable behavior of the switch for the cooperativity of 3 in the presence of leaky expression. Higher cooperativity led to higher leaky expression tolerance, but only to a certain leaky expression threshold, depending on other parameter values, such as production and degradation rates. Higher leaky expression also lowered the reached stable-state levels. The threshold for maximal cooperativity allowing bistability depended, just like leaky expression threshold, on other parameter values, such as production and degradation rates. Figure 4 shows that for very high cooperativity of 10, no bistability was observed at leaky expression rate of 0.05.


Next: C#Sim model of the mutual repressor switch >>