Team:HKUSTers/Modelling

From 2008.igem.org

(Difference between revisions)
Line 122: Line 122:
     <td align="center"id="navactive" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Team" >The Team</a> </td>
     <td align="center"id="navactive" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Team" >The Team</a> </td>
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Project" >The Project</a> </td>  
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Project" >The Project</a> </td>  
-
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Parts" >Parts to the Registry</a> </td>
+
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Parts" >Parts</a> </td>
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Modelling">Modelling</a> </td>
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Modelling">Modelling</a> </td>
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Notebook" >Notebook</a> </td>
     <td align="center" ><a class="mainLinks" href="https://2008.igem.org/Team:HKUSTers/Notebook" >Notebook</a> </td>

Revision as of 03:33, 30 October 2008

HK Harbour-1-2.jpg

Simulation Demo

Please download and unzip this executable file, it's fun to play with the parameters!


File:HKUSTers simulation.zip


The parameters are explained bellow. If you click "Simulate All", an "output.txt" file would be generated that piles up key-information of all simulations from a predefined parameter range.


Brief Description

  • our design: based on random diffusion, signal amplified by reciprocal negative feedback loop
  • simulation shows: when requirements for a critical parameter is met, the Randomizer is feasible and robust to other molecular fluctuations.

Assumptions (selected)

  • T7 Polymerase concentration is low and stable over time
  • Protein degradation is insignificant in the time scale of fate-decision (the first few mRNA transcription events has a preponderant effect)
  • Cell devision is inhibited by a DNA-replication inhibitor, so the volume of the Cell remains constant (e.g. "total number" represents concentration).
  • The expression of the Reporters of the two arms of the switch (GFP and RFP, respectively) is ideally coupled to the expression of the repressors at one-to-one ratio.

List of the Main Parameters

a1 - time a2 - time a3 - probability of promoter bias a4 - probability of RepLeft blockage a5 - probability of RepRight blockage

(..under construction..)

flow of events

  1. One T7 RNA Polymerase binds to the bidirectional promotors (one "Pol-Hit"), with a probability of a3 binding to the LeftPromoter, then pauses (abortive transcription and conformational change). The time interval between 2 Pol-Hits is described by the variable "Hit-interval". This parameter fluctuates above a minimum of a1 seconds.
  2. Polymerase proceed into the elongation state but may be blocked by a repressor sitting at the repressor binding site. i.e. at Decision Timepoint, there is a certain probability that the event is terminated at this step.
  3. If the event is not blocked, transcription is followed by translation, and the nascent repressor protein folds, diffuses and may bind to its repressor binding site (on the other arm of the switch). In one Pol-Hit event, we define a variable "bindingDelay" as the time between the Decision Timepoint and timepoint where the first protein transcribed from this mRNA is functionally available. This variable fluctuates above a minimum of a2 seconds.
  4. We assume that the average life-time of this mRNA is around one minute, so in our model, after the first repressor protein is readily diffusing in the cell, a "burst" of proteins (number of molecules fluctuating between 10 and 20) are subsequently generated from this mRNA at an interval of 2 seconds.
  5. The repression at the repressor binding site is concentration dependent, which we simulate with
    Probability of Effective Blockage = 1-(1-Eff)^n(t),


where Eff reflects the inherent efficiency of the repressor (assign parameters Eff = a4 for the RepLeft and Eff = a5 for the RepRight, Eff ∈(0,1)), and n(t) is the total number of available repressors at timepoint t.

  1. This Loop of a single Pol-Hit Event ends here. As mentioned above, the next Pol-Hit occurs after a time-interval described by the variable "Hit-interval"(=a1+fluctuation).
  2. For each cell, we trace 100 successive Pol-Hit Events. 200 cells in total are sampled.


Qualitative Analysis

Evaluation of Randomizer Performance

Two Criteria:

  1. Our "Dice" - each cell - should deliver the binary output by displaying distinct red or green colour, not something yellowish...

for each cell, Define an output as distinct if 1) (# of green proteins + # of red proteins) > 500 (demand sufficient strength of output) AND 2) r<0.05 or r>0.95, with r =(# of green proteins)/(# of green proteins + # of red proteins), corresponding to distinct red and distinct green, respectively for the total pool of cells, Define SuccessRate = (# of distinct cells)/(total # of cells)

  1. The "Dice" should not be biased:

Define Bias = (# of distinct green cells - # of distinct red cells)/ (# of distinct cells)



the cryptic code behind..

File:SimulationDlg.zip