Team:Slovenia/ModelingPositiveFeedbackLoopSwitch

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" />
 +
 +
<!-- back to top -->
 +
<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 24:
-
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 201: Line 212:
#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 217: Line 229:
/* 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*/
ul {
ul {
Line 292: Line 310:
<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 317: Line 336:
<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 345:
<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 353:
<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 361:
  <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 377:
  <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 400:
<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 386: Line 408:
</div>
</div>
<!-- end main menu -->
<!-- end main menu -->
 +
 +
 +
 +
 +
</div> <!-- end menu -->
</div> <!-- end menu -->
<div id="main">
<div id="main">
 +
 +
<!-- 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>
 +
<br/>
<br/>
<h1>Modeling - positive feedback loop switch</h1>
<h1>Modeling - positive feedback loop switch</h1>
Line 414: Line 448:
<td class="summary" style="font-size:110%;">
<td class="summary" style="font-size:110%;">
<p>
<p>
-
Research suggests that bistability is in theory possible without cooperativity (see Cooperativity). Deterministic analysis showed that the positive feedback switch - due to positive feedback loops and competitive binding of activators and repressors – was exhibiting bistability even at very low functional cooperativity values close to 1.
+
Research suggests that bistability is in theory possible without cooperativity (see <a href="https://2012.igem.org/Team:Slovenia/ModelingMethods#cooperativity">Cooperativity</a>). Deterministic analysis showed that the positive feedback switch - due to positive feedback loops and competitive binding of activators and repressors – was exhibiting bistability even at no or very low functional cooperativity (values equal or close to 1).
</p>
</p>
   
   
Line 436: Line 470:
<h2><a name="model">The model</a></h2>
<h2><a name="model">The model</a></h2>
<p>
<p>
-
We can describe the relations for the mutual repressor switch by the following equations. Fractional occupancies of promoters are:
+
A description of the positive feedback loop switch is available <a href="https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch">here</a>.
-
<p><img src="https://static.igem.org/mediawiki/2012/7/71/Svn12_mrs_det_f1.png"/></p>
+
</p>
-
<p><img src="https://static.igem.org/mediawiki/2012/d/d3/Svn12_mrs_det_f2.png"/></p>
+
 
-
<p><img src="https://static.igem.org/mediawiki/2012/a/a3/Svn12_mrs_det_f3.png"/></p>
+
 
-
<p><img src="https://static.igem.org/mediawiki/2012/3/3d/Svn12_mrs_det_f4.png"/></p>
+
<p>
 +
For modeling purposes, the positive feedback loop switch can be represented as a set of constructs:
 +
<p><img src="https://static.igem.org/mediawiki/2012/c/ca/Svn12_modeling_pfs_shema.png"/></p>
 +
</p>
 +
 
 +
<p><i>[A]</i> is TAL-A:KRAB binding site, <i>[B]</i> is TAL-B:KRAB binding site, <i>pir</i> is PIP:KRAB binding site and <i>etr</i> is E:KRAB binding site. The presence of signal 1 causes PIP:KRAB to dissociate from promoter 3 (construct 3) binding site. The presence of signal 2 causes E:KRAB to dissociate from promoter 4 (construct 4) binding site.</p>
 +
 
 +
 
 +
 
 +
<p>
 +
Fractional occupancies of promoters are:
 +
<p><img src="https://static.igem.org/mediawiki/2012/0/0c/Svn12_pfs_det_f_eqs.PNG"/></p>
 +
 
where:
where:
<ul style="margin-left:30px;">
<ul style="margin-left:30px;">
-
<li>f<sub>1</sub>, f<sub>2</sub>, f<sub>3</sub> and f<sub>4</sub> are probabilities of promoters 1 (construct 1), 2 (construct 2), 3 (construct 3) and 4 (construct 4), respectively, being in an active state, resulting in gene expression;</li>
+
<li>f<sub>1</sub>, f<sub>2</sub>, f<sub>3</sub> and f<sub>4</sub> are probabilities of constructs’ promoters 1, 2, 3 and 4, respectively, being in an active state, resulting in gene expression;</li>
-
<li>[TAL-A:KRAB], [TAL-B:KRAB], [PIP:KRAB] and [E:KRAB] are protein concentrations at a given time;</li>
+
<li>[TAL-A:KRAB], [TAL-B:KRAB], [TAL-A:VP16], [TAL-B:VP16], [PIP:KRAB] and [E:KRAB] are protein concentrations at a given time;</li>
-
<li>k<sub>1</sub>, k<sub>2</sub>, k<sub>3</sub> and k<sub>4</sub> are association constants;</li>
+
<li>k<sub>1</sub>, k<sub>2</sub>, k<sub>3</sub>, k<sub>4</sub>, k<sub>5</sub> and k<sub>6</sub> are association constants;</li>
-
<li>n<sub>1</sub>, n<sub>2</sub>, n<sub>3</sub> and n<sub>4</sub> are exponents representing the degree of functional cooperativity;</li>
+
<li>n<sub>1</sub>, n<sub>2</sub>, n<sub>3</sub>, n<sub>4</sub>, n<sub>5</sub> and n<sub>6</sub> are exponents representing the degree of functional cooperativity;</li>
-
<li>K<sub>r</sub> is the amount of repressor required for 50% repression of constitutive promoter (equal to 1 in our simulations);</li>  
+
<li>K<sub>r</sub> is the amount of repressor required for 50% repression of constitutive promoter;</li>
 +
<li>K<sub>a</sub> is the amount of activator required for 50% activation of minimal promoter.</li>
</ul>
</ul>
Line 455: Line 502:
<p>
<p>
-
ODEs representing protein production are described by a set of equations:
+
Protein production is described by a set of equations:
<p><img src="https://static.igem.org/mediawiki/2012/0/09/Svn12_mrs_det_d1_bfp.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/0/09/Svn12_mrs_det_d1_bfp.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/7/74/Svn12_mrs_det_d2_mcit.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/7/74/Svn12_mrs_det_d2_mcit.png"/></p>
-
<p><img src="https://static.igem.org/mediawiki/2012/a/a6/Svn12_mrs_det_d3_talbkrab.png"/></p>
+
<p><img src="https://static.igem.org/mediawiki/2012/e/e8/Svn12_det_pfs_d3_bkrab.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/a/ae/Svn12_mrs_det_d4_talakrab.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/a/ae/Svn12_mrs_det_d4_talakrab.png"/></p>
 +
<p><img src="https://static.igem.org/mediawiki/2012/2/28/Svn12_det_pfs_d5_bvp.png"/></p>
 +
<p><img src="https://static.igem.org/mediawiki/2012/e/e2/Svn12_det_pfs_d6_avp.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/b/bf/Svn12_mrs_det_dpip_mrs_new.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/b/bf/Svn12_mrs_det_dpip_mrs_new.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/6/65/Svn12_mrs_det_dekrab.png"/></p>
<p><img src="https://static.igem.org/mediawiki/2012/6/65/Svn12_mrs_det_dekrab.png"/></p>
Line 467: Line 516:
where:
where:
<ul style="margin-left:30px;">
<ul style="margin-left:30px;">
-
     <li>[BFP], [mCitrine], [TAL-A:KRAB], [TAL-B:KRAB], [PIP:KRAB] and [E:KRAB] are protein concentrations;</li>
+
      
-
    <li>k<sub>BFP</sub> is BFP production rate from construct 1 (i.e. production rate when construct 1 promoter is active);</li>
+
<li>[BFP], [mCitrine], [TAL-A:KRAB], [TAL-B:KRAB], [TAL-A:VP16], [TAL-B:VP16], [PIP:KRAB] and [HetAct] are protein concentrations;
-
    <li>kb<sub>BFP</sub> is basal BFP production rate from construct 1 (i.e. production rate when construct 1 promoter is inactive);</li>
+
</li>
-
     <li>deg<sub>BFP</sub> is BFP degradation rate;</li>
+
 
-
    <li>k<sub>cit</sub> is mCitrine production rate from construct 2 (i.e. production rate when construct 2 promoter is active);</li>
+
 
-
    <li>kb<sub>cit</sub> is basal mCitrine production rate from construct 2 (i.e. production rate when construct 2 promoter in inactive);</li>
+
<li>k<sub>BFP</sub> is BFP production rate from construct 1;
-
    <li>deg<sub>cit</sub> is mCitrine degradation rate;</li>
+
</li>
-
    <li>k<sub>2AKR</sub> is TAL-A:KRAB production rate from construct 2;</li>
+
 
-
    <li>kb<sub>2AKR</sub> is basal TAL-A:KRAB production rate from construct 2;</li>
+
<li>kb<sub>BFP</sub> is leaky BFP production rate from construct 1;
-
    <li>k<sub>4AKR</sub> is TAL-A:KRAB production rate from construct 4;</li>
+
</li>
-
     <li>kb<sub>4AKR</sub> is basal TAL-A:KRAB production rate from construct 4;</li>
+
      
-
     <li>deg<sub>AKR</sub> is TAL-A:KRAB degradation rate;</li>
+
 
-
     <li>k<sub>1BKR</sub> is TAL-B:KRAB production rate from construct 1;</li>
+
<li>deg<sub>BFP</sub> is BFP degradation rate;</li>
-
     <li>kb<sub>1BKR</sub> Is basal TAL-B:KRAB production rate from construct 1;</li>
+
 
-
     <li>k<sub>3BKR</sub> is TAL-B:KRAB production rate from construct 3;</li>
+
 
-
     <li>kb<sub>3BKR</sub> is basal TAL-B:KRAB production rate from construct 3;</li>
+
<li>k<sub>cit</sub> is mCitrine production rate from construct 2;
-
     <li>deg<sub>BKR</sub> is TAL-B:KRAB degradation rate;</li>
+
</li>
-
     <li>k<sub>PIP</sub> is PIP:KRAB production rate;</li>
+
 
-
    <li>deg<sub>PIP</sub> is PIP:KRAB degradation rate;</li>
+
 
-
    <li>k<sub>E</sub> is E:KRAB production rate;</li>
+
<li>kb<sub>cit</sub> is leaky mCitrine production rate from construct 2;
-
    <li>deg<sub>E</sub> is E:KRAB degradation rate.</li>
+
</li>
 +
 
 +
 
 +
<li>deg<sub>cit</sub> is mCitrine degradation rate;
 +
</li>
 +
 
 +
 
 +
<li>k<sub>2AKR</sub> is TAL-A:KRAB production rate from construct 2;
 +
</li>
 +
 
 +
 
 +
<li>kb<sub>2AKR</sub> is leaky TAL-A:KRAB production rate from construct 2;
 +
</li>
 +
 
 +
 
 +
<li>k<sub>4AKR</sub> is TAL-A:KRAB production rate from construct 4;
 +
</li>
 +
      
 +
 
 +
<li>kb<sub>4AKR</sub> is leaky TAL-A:KRAB production rate from construct 4;
 +
</li>
 +
      
 +
 
 +
<li>deg<sub>AKR</sub> is TAL-A:KRAB degradation rate;
 +
</li>
 +
      
 +
 
 +
<li>k<sub>1BKR</sub> is TAL-B:KRAB production rate from construct 1;
 +
</li>
 +
      
 +
<li>kb<sub>1BKR</sub> Is leaky TAL-B:KRAB production rate from construct 1;
 +
</li>
 +
      
 +
 
 +
<li>k<sub>3BKR</sub> is TAL-B:KRAB production rate from construct 3;
 +
</li>
 +
      
 +
 
 +
<li>kb<sub>3BKR</sub> is leaky TAL-B:KRAB production rate from construct 3;
 +
</li>
 +
      
 +
 
 +
<li>deg<sub>BKR</sub> is TAL-B:KRAB degradation rate;
 +
</li>
 +
      
 +
 
 +
<li>k1<sub>AVP</sub> is TAL-A:VP16 production rate from construct 1;
 +
</li>
 +
   
 +
 
 +
<li>kb<sub>1AVP</sub> is leaky TAL-A:VP16 production rate from construct 1;
 +
</li>
 +
   
 +
 
 +
<li>k<sub>3AVP</sub> is TAL-A:VP16 production rate from construct 3;
 +
</li>
 +
   
 +
 
 +
<li>kb<sub>3AVP</sub> is leaky TAL-A:VP16 production rate from construct 3;
 +
</li>
 +
   
 +
 
 +
<li>deg<sub>AVP</sub> is TAL-A:VP16 degradation rate;
 +
</li>
 +
 
 +
 
 +
<li>k<sub>2BVP</sub> is TAL-B:VP16 production rate from construct 2;</li>
 +
 
 +
 
 +
 
 +
<li>kb<sub>2BVP</sub> is leaky TAL-B:VP16 production rate from construct 2;
 +
</li>
 +
 
 +
<li>k<sub>4BVP</sub> is TAL-B:VP16 production rate from construct 4;
 +
</li>
 +
 
 +
 
 +
<li>kb<sub>4BVP</sub> is leaky TAL-B:VP16 production rate from construct 4;
 +
</li>
 +
 
 +
 
 +
<li>deg<sub>BVP</sub> is TAL-B:VP16 degradation rate;
 +
</li>
 +
 
 +
<li>k<sub>PIP</sub> is PIP:KRAB production rate;
 +
</li>
 +
 
 +
<li>deg<sub>PIP</sub> is PIP:KRAB degradation rate;
 +
</li>
 +
 
 +
<li>k<sub>E</sub> is E:KRAB production rate;
 +
</li>
 +
 
 +
<li>deg<sub>E</sub> is E:KRAB degradation rate;
 +
</li>
 +
 
</ul></p>
</ul></p>
<p>
<p>
Line 496: Line 640:
<h2><a name="results">Simulation results</a></h2>
<h2><a name="results">Simulation results</a></h2>
 +
<p>
 +
Similar simulation methods and assumptions as for the mutual repressor switch were used. Again, switching between states was achieved using two signals, each of which was used to induce a different stable state. The presence of signal 1 caused PIP:KRAB to dissociate from promoter 3 binding site. The presence of signal 2 caused E:KRAB to dissociate from promoter 4 binding site.
 +
</p>
 +
 +
<p>
 +
The ratio between the amount of activator required for 50% activation in case of a minimal promoter to the amount of repressor required for 50% repression in case of a constitutive promoter was set to 3.
 +
</p>
 +
 +
<p>
 +
Initial species' concentrations were 0 for all simulations.
 +
</p>
 +
 +
<p>Detailed parameter values for each simulation can be found in MATLAB files <a href="https://2012.igem.org/Team:Slovenia/SourceCode
 +
">here</a>.</p>
 +
 +
<h3>Simulation set 1 - switching between stable states</h3>
 +
<p>
 +
The purpose of this simulation was to show that deterministic conditions exist where the switch would exhibit stable states. In order to compare this switch with the mutual repressor switch, the same parameter values were used: functional cooperativity of each transcription factor (TF) was equal to 2.5; product production-degradation ratio was equal to 100 (maximal possible expression level); zero leaky expression was assumed. Figure 1 shows the result of this simulation, indicating that stable states are possible. The state-switching scenario used was the same as the one in simulation 1 for mutual repressor switch:
 +
<ul style="margin-left:30px;">
 +
    <li>at time = 0, signal 1 was introduced  (signal 1 = 1) to induce the first stable state, depicted by a high level of BFP;</li>
 +
    <li>at time = 100, signal 1 was removed from the system (signal 1 = 0); high-BFP state persisted even after that, meaning a stable state was achieved;</li>
 +
    <li>at time = 400, signal 2 was introduced (signal 2 = 1) to induce the second, high-mCitrine stable state;</li>
 +
    <li>at time = 500, signal 2 was removed from the system (signal 2 = 0); high-mCitrine state persisted, meaning a stable state was achieved;</li>
 +
    <li>at time = 800, signal 1 was re-introduced (signal 1 = 1) to show that a transition to high-BFP state ce be done from a high-mCitrine state;</li>
 +
    <li>at time = 900, signal 1 was removed (signal 1 = 0); high-BFP stable state persisted.</li>
 +
 +
</ul>
 +
</p>
 +
 +
 +
<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/a/ac/Svn12_modeling_PositiveLoopSwitch_test1.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> Positive feedback loop switch transitioning between stable states. Reporter (BFP, mCitrine) concentration is depicted as a function of time. Here, stable states were observed under ideal conditions with no leaky expression and TFs' cooperativity values equal to 2.5.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
 +
<p>
 +
Comparing the result to the mutual repressor switch under similar conditions, it is obvious that here, each steady state is reached sooner. For the mutual repression switch, the transition to a recently-induced state – i.e. rise in its expression level – only began after the level of the previous stable state was very low (close to zero). For the positive feedback loop switch, the rise coincided in time with the introduction of the signal that causes the transition. Hence, stable state was reached more quickly.
 +
</p>
 +
 +
<p>
 +
The mutual repressor switch exhibited no bistability for low cooperativity values, such as 1.15 (when leaky expression was present) or very low stable-state level bistability (when no or extremely little leaky expression was present). The positive feedback loop switch, however, exhibited bistability with high stable-state levels even at this cooperativity (1.15), as shown in Figure 2. This indicated that the positive feedback loop switch was much less dependent on cooperativity and could exhibit high expression-level bistability even at low cooperativity values. Bistability was exhibited with no noticeable change in reached levels even when PIP:KRAB and E:KRAB production rate was only half of TAL production rates.
 +
</p>
 +
<p>
<p>
Line 503: Line 707:
<tr class="invisible">
<tr class="invisible">
<td class="invisible">
<td class="invisible">
-
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/d/dc/Svn12_MutualSwitch_test2b.png"/>
+
<img  class="invisible" src="https://static.igem.org/mediawiki/2012/d/de/Svn12_modeling_PositiveLoopSwitch_test2.png"/>
</td>
</td>
</tr>
</tr>
Line 511: Line 715:
<tbody  class="invisible">
<tbody  class="invisible">
<tr class="normal"><td class="invisible">
<tr class="normal"><td class="invisible">
-
<b>Figure 2.</b> Degraded bistability of the mutual repressor toggle switch for cooperativity of 1.15. In comparison to  Figure 1 - the same state-switching scenario applies - bistability was exhibited, but very low difference between the states was achieved when cooperativity was set to 1.15 (other parameters were identical to parameters of Figure 1 ) and no leaky transcription was present. Increasing the the production rate of E:KRAB and PIP:KRAB slightly increased the levels, but even a 100-fold increase in E:KRAB and PIP:KRAB production rate compared to other protein production rates did not produce significantly higher stable-state levels. Introducing residual TAL expression of as little as 1% resulted in the loss of bistability for the cooperativity of 1.15.
+
<b>Figure 2.</b> Positive feedback loop switch transitioning between stable states for TF cooperativity of 1.15. Bistability was preserved, with low cooperativity values only slightly reducing stable state expression levels. No cooperativity of PIP:KRAB and E:KRAB was required for the same result to be produced.
</td></tr>
</td></tr>
</tbody>
</tbody>
Line 519: Line 723:
 +
<h3>Simulation set 2 - cooperativity and leaky expression</h3>
 +
<p>
 +
As with the mutual repressor switch, higher cooperativity values were found to increase the robustness of the switch. However, the robustness at a given cooperativity was much higher for the positive feedback loop switch. Figure 3 shows that even for a low cooperativity of 1.15, the switch was exhibiting bistability with high-expression levels in the presence of leaky expression of 0.1. Cooperativity of 1.3 was enough to preserve high-level (above 80%) bistability in the presence of leaky expression of 0.5. Increasing cooperativity to a still relatively low value of 1.85 dramatically increased switch’s capability to retain bistable behavior. In such case, as shown in Figure 4, even a leaky expression of 1.2 for each TAL did not prevent bistability from occuring, although the expression levels were reduced slightly.
 +
</p>
 +
<p>This results suggested that positive feedback loop switch was far more robust and leaky expression-tolerant than the mutual repressor switch.
 +
</p>
 +
 +
<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/3/3b/Svn12_modeling_PositiveLoopSwitch_test3.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> Positive feedback loop switch state transitions. Cooperativity of 1.15 was enough to preserve high-level bistability in the presence of leaky expression of each TAL equal to 0.1.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
<br/>
 +
 +
<p>
 +
<!-- 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/c6/Svn12_modeling_PositiveLoopSwitch_test4.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> Positive feedback loop switch state transitions. Cooperativity of 1.85 was enough for the positive feedback loop switch to preserve bistability even for leaky expression as high as 1.2, thus demonstrating superior robustness of the switch compared to the mutual repressor switch.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
 +
 +
<h3>Simulation set 3 - parameter evaluation</h3>
 +
<p>
 +
As for the mutual repressor switch, parameter evaluation was performed to search bistabiliy-exhibiting parameter space more thoroughly. The following state-switching scenario was used as a foundation:
 +
<ul style="margin-left:30px";>
 +
    <li>at time = 0, signal 1 was introduced to induce stable state 1 (high BFP);</li>
 +
    <li>at time = 50, signal 1 was removed;</li>
 +
    <li>at time = 100, signal 2 was introduced to induce stable state 2 (high mCitrine);</li>
 +
    <li>at time = 150, signal 2 was removed.</li>
 +
</ul>
 +
</p>
 +
 +
<p>
 +
Default TF cooperativity was 1.5. Product production-degradation ratio was 100. Default leaky expression of all TALs was 0.25, with zero leaky expression of reporters. Simulation results are shown in Figure 5 and show that bistability was achieved.
 +
</p>
 +
 +
 +
<p>
 +
<!-- figure 5 -->
 +
<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/0/03/Svn12_modeling_2dPositiveLoopSwitch_parameterTest1_2d.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 5.</b> Positive feedback loop switch exhibiting bistability. This simulation scenario was used in parameter evaluation tests.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
<p>
 +
Again, we performed several simulations using this switching scenario, each time varying different parameters, producing 3D graphs that showed how each state was affected.
 +
</p>
 +
 +
<p>
 +
First, we tested the effect of TF cooperativity on each state. Figure 6 shows that, for the first stable state, maximal expression was observed for cooperativity of each TF for values only slightly greater than one. Figure 7 suggested similar conclusion for the second stable state.
 +
</p>
 +
 +
 +
<p>
 +
<!-- figure 6 -->
 +
<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/3/31/Svn12_3dPositiveLoopSwitch_parameterTest1_BFP.png"/><br/>
 +
<img class="invisible" src="https://static.igem.org/mediawiki/2012/8/81/Svn12_3dPositiveLoopSwitch_parameterTest1_BFP_top.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 6.</b> Stable state 1 levels dependence on cooperativity. Switching scenario from Figure 5 was used.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
<br/>
 +
 +
 +
<p>
 +
<!-- figure 7 -->
 +
<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/b/b2/Svn12_3dPositiveLoopSwitch_parameterTest1_Cit.png"/><br/>
 +
<img class="invisible" src="https://static.igem.org/mediawiki/2012/b/b3/Svn12_3dPositiveLoopSwitch_parameterTest1_Cit_top.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 7.</b> Stable state 2 levels dependence on ccoperativity. Switching scenario from figure 5 was used.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
 +
 +
<p>
 +
Bistability was retained for relatively high values of leaky expression (compared to the mutual repressor switch), as shown in Figure 8 (depicting first state) and Figure 9 (depicting second state). Higher cooperativity further improved leaky expression tolerance, increasing expression levels.
 +
</p>
 +
 +
 +
 +
 +
<p>
 +
<!-- figure 8 -->
 +
<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/1/16/Svn12_3d_PositiveLoopSwitch_parameterTest2_BFP_b.png"/><br/>
 +
<img class="invisible" src="https://static.igem.org/mediawiki/2012/f/f1/Svn12_3d_PositiveLoopSwitch_parameterTest2_BFP_top_b.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 8.</b> Stable state 1 levels at different TAL leaky expression rates. Each TAL gene was assigned a leaky expression rate shown. Switching scenario from Figure 5 was used. High tolerance to leaky expression was observed compared to the mutual repressor switch.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
<br/>
 +
 +
 +
 +
 +
 +
<p>
 +
<!-- figure 9 -->
 +
<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/0/09/Svn12_3d_PositiveLoopSwitch_parameterTest2_Cit_b.png"/><br/>
 +
<img class="invisible" src="https://static.igem.org/mediawiki/2012/4/44/Svn12_3d_PositiveLoopSwitch_parameterTest2_Cit_top_b.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 9.</b> Stable state 2 levels at different TAL leaky expression rates. Each TAL gene was assigned a leaky expression rate shown. Switching scenario from figure 5 was used. High tolerance to leaky expression was observed compared to the mutual repressor switch.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
 +
 +
 +
<h3>Simulation set 4</h3>
 +
<p>
 +
A concern in our wet-lab was that the positive feedback loops might be too strong, either rendering the switch incapable of switching from one stable state to another due to a large amount of activator produced, or negatively affecting bistability in some other way. Deterministic modeling revealed that in some cases, it is beneficial for bistability if the positive feedback loops (TAL activators) have lower production rates compared to other transcription factor production rates (Figure 10). One way to achieve this in experiments is to have lower amounts of constructs with positive feedback loops present compared to other constructs. The results also demonstrate that bistability is possible without any transcription factor cooperativity even in deterministic models when protein production to degradation rate ratio is high enough (e.g. 10000) and leaky transcription is low.
 +
</p>
 +
 +
 +
<p>
 +
<!-- figure 10 -->
 +
<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/e/e4/Svn12_pfs_razmerja_gif.gif"/>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<table class="invisible" style="width:90%; text-align:center;">
 +
<tbody  class="invisible">
 +
<tr class="normal"><td class="invisible">
 +
<b>Figure 10.</b> The effect of different positive feedback loops production rates on bistability at constant production rates of other transcription factors, for which the production to degradation rate ratio was equal to 10000. No cooperativity was assumed and leaky transcription was equal to 1%.
 +
</td></tr>
 +
</tbody>
 +
</table>
 +
<!-- end table-->
 +
</p>
Line 526: Line 961:
Next: <a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitchStochastic'>Stochastic model of the positive feedback loop switch >></a>
Next: <a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitchStochastic'>Stochastic model of the positive feedback loop switch >></a>
</b>
</b>
-
 
</div>
</div>

Latest revision as of 21:03, 26 October 2012


Modeling - positive feedback loop switch

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

Deterministic model of the positive feedback loop switch

Research suggests that bistability is in theory possible without cooperativity (see Cooperativity). Deterministic analysis showed that the positive feedback switch - due to positive feedback loops and competitive binding of activators and repressors – was exhibiting bistability even at no or very low functional cooperativity (values equal or close to 1).

The switch proved much more robust than the mutual repressor switch, exhibiting bistable behavior with high stable-state expression levels even when leaky TAL expression was high. As with the mutual repressor switch, increase in cooperativity further improved robustness. However, even at a relatively low cooperativity, the switch was robust enough to suggest that an experimental realization may be possible.

Another advantage of the positive feedback loop switch over the mutual repressor switch, according to our deterministic model, was faster transition from one stable state to another after induction.

Our experimental results showed bistable behavior of the switch, as predicted by the model.

The model

A description of the positive feedback loop switch is available here.

For modeling purposes, the positive feedback loop switch can be represented as a set of constructs:

[A] is TAL-A:KRAB binding site, [B] is TAL-B:KRAB binding site, pir is PIP:KRAB binding site and etr is E:KRAB binding site. The presence of signal 1 causes PIP:KRAB to dissociate from promoter 3 (construct 3) binding site. The presence of signal 2 causes E:KRAB to dissociate from promoter 4 (construct 4) binding site.

Fractional occupancies of promoters are:

where:
  • f1, f2, f3 and f4 are probabilities of constructs’ promoters 1, 2, 3 and 4, respectively, being in an active state, resulting in gene expression;
  • [TAL-A:KRAB], [TAL-B:KRAB], [TAL-A:VP16], [TAL-B:VP16], [PIP:KRAB] and [E:KRAB] are protein concentrations at a given time;
  • k1, k2, k3, k4, k5 and k6 are association constants;
  • n1, n2, n3, n4, n5 and n6 are exponents representing the degree of functional cooperativity;
  • Kr is the amount of repressor required for 50% repression of constitutive promoter;
  • Ka is the amount of activator required for 50% activation of minimal promoter.

Protein production is described by a set of equations:

where:
  • [BFP], [mCitrine], [TAL-A:KRAB], [TAL-B:KRAB], [TAL-A:VP16], [TAL-B:VP16], [PIP:KRAB] and [HetAct] are protein concentrations;
  • kBFP is BFP production rate from construct 1;
  • kbBFP is leaky BFP production rate from construct 1;
  • degBFP is BFP degradation rate;
  • kcit is mCitrine production rate from construct 2;
  • kbcit is leaky mCitrine production rate from construct 2;
  • degcit is mCitrine degradation rate;
  • k2AKR is TAL-A:KRAB production rate from construct 2;
  • kb2AKR is leaky TAL-A:KRAB production rate from construct 2;
  • k4AKR is TAL-A:KRAB production rate from construct 4;
  • kb4AKR is leaky TAL-A:KRAB production rate from construct 4;
  • degAKR is TAL-A:KRAB degradation rate;
  • k1BKR is TAL-B:KRAB production rate from construct 1;
  • kb1BKR Is leaky TAL-B:KRAB production rate from construct 1;
  • k3BKR is TAL-B:KRAB production rate from construct 3;
  • kb3BKR is leaky TAL-B:KRAB production rate from construct 3;
  • degBKR is TAL-B:KRAB degradation rate;
  • k1AVP is TAL-A:VP16 production rate from construct 1;
  • kb1AVP is leaky TAL-A:VP16 production rate from construct 1;
  • k3AVP is TAL-A:VP16 production rate from construct 3;
  • kb3AVP is leaky TAL-A:VP16 production rate from construct 3;
  • degAVP is TAL-A:VP16 degradation rate;
  • k2BVP is TAL-B:VP16 production rate from construct 2;
  • kb2BVP is leaky TAL-B:VP16 production rate from construct 2;
  • k4BVP is TAL-B:VP16 production rate from construct 4;
  • kb4BVP is leaky TAL-B:VP16 production rate from construct 4;
  • degBVP is TAL-B:VP16 degradation rate;
  • kPIP is PIP:KRAB production rate;
  • degPIP is PIP:KRAB degradation rate;
  • kE is E:KRAB production rate;
  • degE is E:KRAB degradation rate;

See model derivation for details.

Simulation results

Similar simulation methods and assumptions as for the mutual repressor switch were used. Again, switching between states was achieved using two signals, each of which was used to induce a different stable state. The presence of signal 1 caused PIP:KRAB to dissociate from promoter 3 binding site. The presence of signal 2 caused E:KRAB to dissociate from promoter 4 binding site.

The ratio between the amount of activator required for 50% activation in case of a minimal promoter to the amount of repressor required for 50% repression in case of a constitutive promoter was set to 3.

Initial species' concentrations were 0 for all simulations.

Detailed parameter values for each simulation can be found in MATLAB files here.

Simulation set 1 - switching between stable states

The purpose of this simulation was to show that deterministic conditions exist where the switch would exhibit stable states. In order to compare this switch with the mutual repressor switch, the same parameter values were used: functional cooperativity of each transcription factor (TF) was equal to 2.5; product production-degradation ratio was equal to 100 (maximal possible expression level); zero leaky expression was assumed. Figure 1 shows the result of this simulation, indicating that stable states are possible. The state-switching scenario used was the same as the one in simulation 1 for mutual repressor switch:

  • at time = 0, signal 1 was introduced (signal 1 = 1) to induce the first stable state, depicted by a high level of BFP;
  • at time = 100, signal 1 was removed from the system (signal 1 = 0); high-BFP state persisted even after that, meaning a stable state was achieved;
  • at time = 400, signal 2 was introduced (signal 2 = 1) to induce the second, high-mCitrine stable state;
  • at time = 500, signal 2 was removed from the system (signal 2 = 0); high-mCitrine state persisted, meaning a stable state was achieved;
  • at time = 800, signal 1 was re-introduced (signal 1 = 1) to show that a transition to high-BFP state ce be done from a high-mCitrine state;
  • at time = 900, signal 1 was removed (signal 1 = 0); high-BFP stable state persisted.

Comparing the result to the mutual repressor switch under similar conditions, it is obvious that here, each steady state is reached sooner. For the mutual repression switch, the transition to a recently-induced state – i.e. rise in its expression level – only began after the level of the previous stable state was very low (close to zero). For the positive feedback loop switch, the rise coincided in time with the introduction of the signal that causes the transition. Hence, stable state was reached more quickly.

The mutual repressor switch exhibited no bistability for low cooperativity values, such as 1.15 (when leaky expression was present) or very low stable-state level bistability (when no or extremely little leaky expression was present). The positive feedback loop switch, however, exhibited bistability with high stable-state levels even at this cooperativity (1.15), as shown in Figure 2. This indicated that the positive feedback loop switch was much less dependent on cooperativity and could exhibit high expression-level bistability even at low cooperativity values. Bistability was exhibited with no noticeable change in reached levels even when PIP:KRAB and E:KRAB production rate was only half of TAL production rates.

Simulation set 2 - cooperativity and leaky expression

As with the mutual repressor switch, higher cooperativity values were found to increase the robustness of the switch. However, the robustness at a given cooperativity was much higher for the positive feedback loop switch. Figure 3 shows that even for a low cooperativity of 1.15, the switch was exhibiting bistability with high-expression levels in the presence of leaky expression of 0.1. Cooperativity of 1.3 was enough to preserve high-level (above 80%) bistability in the presence of leaky expression of 0.5. Increasing cooperativity to a still relatively low value of 1.85 dramatically increased switch’s capability to retain bistable behavior. In such case, as shown in Figure 4, even a leaky expression of 1.2 for each TAL did not prevent bistability from occuring, although the expression levels were reduced slightly.

This results suggested that positive feedback loop switch was far more robust and leaky expression-tolerant than the mutual repressor switch.


Simulation set 3 - parameter evaluation

As for the mutual repressor switch, parameter evaluation was performed to search bistabiliy-exhibiting parameter space more thoroughly. The following state-switching scenario was used as a foundation:

  • at time = 0, signal 1 was introduced to induce stable state 1 (high BFP);
  • at time = 50, signal 1 was removed;
  • at time = 100, signal 2 was introduced to induce stable state 2 (high mCitrine);
  • at time = 150, signal 2 was removed.

Default TF cooperativity was 1.5. Product production-degradation ratio was 100. Default leaky expression of all TALs was 0.25, with zero leaky expression of reporters. Simulation results are shown in Figure 5 and show that bistability was achieved.

Again, we performed several simulations using this switching scenario, each time varying different parameters, producing 3D graphs that showed how each state was affected.

First, we tested the effect of TF cooperativity on each state. Figure 6 shows that, for the first stable state, maximal expression was observed for cooperativity of each TF for values only slightly greater than one. Figure 7 suggested similar conclusion for the second stable state.


Bistability was retained for relatively high values of leaky expression (compared to the mutual repressor switch), as shown in Figure 8 (depicting first state) and Figure 9 (depicting second state). Higher cooperativity further improved leaky expression tolerance, increasing expression levels.


Simulation set 4

A concern in our wet-lab was that the positive feedback loops might be too strong, either rendering the switch incapable of switching from one stable state to another due to a large amount of activator produced, or negatively affecting bistability in some other way. Deterministic modeling revealed that in some cases, it is beneficial for bistability if the positive feedback loops (TAL activators) have lower production rates compared to other transcription factor production rates (Figure 10). One way to achieve this in experiments is to have lower amounts of constructs with positive feedback loops present compared to other constructs. The results also demonstrate that bistability is possible without any transcription factor cooperativity even in deterministic models when protein production to degradation rate ratio is high enough (e.g. 10000) and leaky transcription is low.


Next: Stochastic model of the positive feedback loop switch >>