Team:Slovenia/ModelingMutualRepressorSwitchStochastic

From 2012.igem.org

(Difference between revisions)
 
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 307: Line 319:
<body>
<body>
 +
<a name="topofthepage" style="background-color:transparent;"></a>
 +
<div id="banner">
<div id="banner">
Line 314: 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 331: 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'><table 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/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 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>  
+
     <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 340: 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'><table 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>  
+
     <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 348: 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 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>  
+
     <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 356: 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'><table 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/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'><table 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/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'><table 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>  
+
     <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 372: 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'><table 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/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 395: 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 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/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 412: Line 425:
<!-- dummy link na bannerju -->
<!-- dummy link na bannerju -->
<a href="https://2012.igem.org/Main_Page">
<a href="https://2012.igem.org/Main_Page">
-
<div id="dummy" style="background-color:white; position:absolute; left:870px; top:25px; width:115px; height:80px; z-index:100; opacity:0.0;">
+
<div id="dummy" style="background-color:transparent; position:absolute; left:870px; top:25px; width:115px; height:80px; z-index:100; opacity:0.0;">
</div>
</div>
</a>
</a>

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 >>