Team:BGU Israel/Model2.html

From 2013.igem.org

(Difference between revisions)
 
(11 intermediate revisions not shown)
Line 22: Line 22:
Leakage rate was set between 0.1-0.5, and mechanism strengths were given numbers between 1-5, representing the 'a' parameter: protein concentrations functions.</br></br>
Leakage rate was set between 0.1-0.5, and mechanism strengths were given numbers between 1-5, representing the 'a' parameter: protein concentrations functions.</br></br>
-
<b>Model Results</b></br></br>
+
<b>Model Results</b></br>
-
 
+
<ol class="bulletlist">
-
 
+
        <li class="bulletlist">Die-out timeframe: 300-750 minutes.</li>
-
How does a single combination simulation of protein leakage and Mechanism strength look like? </br>Let&rsquo;s see a few representing examples:</br></br>
+
        <li class="bulletlist">Working timeframe: 190-600 minutes.</li>
-
<b>Weakest mechanism, strongest leakage:</b></br>
+
        <li class="bulletlist">Generation range: 13-32 generations.</li>
 +
</ol>
 +
<p>
 +
What does a simulation of one of these combinations of protein leakage and mechanism strength look like? In each example below, 5 experiments out of the 30 are marked in grey, an average curve is red, and standard deviation curves are black.</br></br>
 +
</p>
</p>
</p>
<div style="margin-left:30px;">
<div style="margin-left:30px;">
-
<img src="https://static.igem.org/mediawiki/2013/0/04/BGU_1comb5-5.png" /></br>
+
<img src="https://static.igem.org/mediawiki/2013/d/d9/Bgu_model_1.png" /></br>
-
<h8>Figure 1: 5 curves represent 5 repeats of the experiment with same initial conditions.</h8></br></br>
+
<h8>Strongest mechanism, Weakest leakage rate</h8></br></br>
</div>
</div>
<p>
<p>
-
In this one it is visible that when the mechanism is weak and the protein is leakage rate is high, the population behaves exponentially, variance is visible, though not major. </br></br>
+
Example 1: Strongest mechanism, weakest leakage rate. All experiments end with the entire population dying out. Cells experience a log phase of around 90 minutes, a maximum population of around 230 is reached, and the variance is small.</br></br>
 +
 
-
<b>Strongest mechanism, weakest leakage: </b></br>
 
</p>
</p>
<div style="margin-left:30px;">
<div style="margin-left:30px;">
-
<img src="https://static.igem.org/mediawiki/2013/4/49/BGU_2comb1-1.png" /></br>
+
<img src="https://static.igem.org/mediawiki/2013/b/ba/Bgu_model_2.png" /></br>
-
<h8>Figure 2</h8></br></br>
+
<h8>Intermediate combination. Mechanism strength=4, Leakage rate=0.4</h8></br></br>
</div>
</div>
<p>
<p>
-
We can see that although variance is bigger, all experiments end up with all the population dying out. Cells experience a log phase of around 100 minutes, and at low count some of them survive around 100 minutes more.</br></br>
+
Example 2: Intermediate parameters. Here, the stochasticity of the system takes an important role. Some of the experiments are trending down at the end of the timeframe, while others are trending up. Variance is much bigger. Log phase is about 200 minutes. A maximum population of around 410 is reached.</br></br>
-
<b>Intermediate:</b></br></br>
+
In order to compare all of the experiments, 3D surfaces projecting the data of all combinations were made. In each one of figures 4-6, a specific parameter was chosen to be compared between the combinations. The coefficient of variation (CV) of the 30 repeats was computed and is represented by the surface’s color:</br></br>
</p>
</p>
<div style="margin-left:30px;">
<div style="margin-left:30px;">
-
<img src="https://static.igem.org/mediawiki/2013/0/06/BGU_3comb4-4.png" /></br>
+
<img src="https://static.igem.org/mediawiki/2013/b/b4/Bgu_model_3.png" height="700" width="1100"/></br>
-
<h8>Figure 3</h8></br></br>
+
<h8>Figure 3: Z axis shows the time it took for 95% of the population to die out. X and Y are the values for the tested parameters, the bacteria symbol represents a non- reliable combination where less than 90% of the experiments died out within the model time frame, and the rest lived on. </h8></br></br>
</div>
</div>
<p>
<p>
-
In intermediate parameters the stochasticity of the system takes important role. While 4 experiments are trending down, on is trending up. Variance is much bigger. Log phase is about 200 minutes.</br></br>
+
The graph omits combinations in which the weakest mechanism strength 1 was used, because those combinations never died out. This parameter is robust and has a clear limit of around 2 as defined in the model.</br></br>
-
In order to show a complete comparable picture, 3d surfaces projecting the data of all combinations were made, in each one of figures 4-6 one specific parameter was chosen to be compared for each combination, the variance between the 5 repeats was computed and normalized to the mean and is represented by the surfaces color:</br></br>
+
The leakage is a more modular parameter, giving delicate tuning for high mechanism strengths, thus strengthening our hypothesis that leakage can provide nodularity for the system.</br></br>
 +
A time frame of around 300-750 minutes is seen, as the z-differential between the lowest combination and the highest one.</br></br>
</p>
</p>
<div style="margin-left:30px;">
<div style="margin-left:30px;">
-
<img src="https://static.igem.org/mediawiki/2013/2/25/BGU_4-95_Threshold-Hot.png" /></br>
+
<img src="https://static.igem.org/mediawiki/2013/4/4a/Bgu_model_4.png" height="700" width="1100"/></br>
-
<h8>Figure 4: Z axis shows the time it took for 95% of the population to die out, X and Y are the values for the tested parameters.</h8></br></br></br>
+
<h8>Figure 4: Z axis is the number of generations reached after 1000 minutes of simulation, in combinations that died out - a high number wasn’t reached</h8></br></br></br>
-
<img src="https://static.igem.org/mediawiki/2013/d/d0/BGU_5-generations2-Hot.png" /></br>
+
<p>A generation frame of 13-32 generations within the effective combinations is seen in this model. This effects the chance of the system to transfer genes horizontally, and the chance for a mutation in the mechanism to occur.</br></br>
-
<h8>Figure 5: Z axis is the number of generations reached after 1000 minutes of simulation, for combinations that died out a high number wasn&rsquo;t reached.</h8></br></br></br>
+
</p>
-
<img src="https://static.igem.org/mediawiki/2013/0/00/BGU_6-50_Threshold-Hot.png" /></br>
+
<img src="https://static.igem.org/mediawiki/2013/1/16/Bgu_model_5.png" height="700" width="1100"/></br>
-
<h8>Figure 6: Z axis is the time it took until population reached a 50% mark from initial conditions.</h8></br></br></br>
+
<h8>Figure 5: Z axis is the time it took until population reached a 50% mark from initial conditions.</h8></br></br></br>
 +
 
 +
<p>
 +
This parameter, the measurement of the half-life, can give an indication for the possible working time for the system. A time frame of 190-600 minutes has been calculated.</br></br>
 +
While inducer concentration usually doesn’t provide adequate control of the system because of the promoter’s bimodality[1] (all or nothing behavior when induced), some control can be achieved from different inducer combinations [2]. Therefore it is possible to induce P.A.S.E-containing bacteria to express only a partial percentage of the maximal protein expression capability, allowing for control of the overall population lifetime.</br></br>
 +
</p>
 +
View our <a href="https://static.igem.org/mediawiki/2013/b/b2/Birth_death_iGEM_BGU.m" target="_blank">code.</a>
</div>
</div>
<p>
<p>
-
In the threshold surfaces, combinations reaching the 1000 minutes plateau never died out.</br></br>
+
</br></br><h6>  Continue the journey: read about our <a href="/Team:BGU_Israel/HPOverview">Human Practice</a>.</h6></br></br>
-
It is important to estimate number of generations because the higher it gets, the higher the probability is for a mutation to happen. Another thing to think about is the working time of the GEM, effective combinations reached the 50% threshold after between 200-500 minutes.</br></br>
+
-
While inducer concentration usually don’t provide adequate control of the system because of the promoter’s bimodality<b>[1]</b> (all or nothing behavior when induced), some control can be achieved from different inducer combinations <b>[2]</b>, it is possible to release the P.A.S.E after being induced to express only partial percentage of the maximal protein expression capability.</br></br></br>
+
<hr/>
<hr/>

Latest revision as of 00:52, 29 October 2013

BGU_Israel

Modelling Results

As explained in the theoretical section on the stochastic birth-death model, 25 combinations of protein leakage and mechanism strength were simulated in-silico. For each combination, 30 repeats ("experiments") were performed.

Leakage rate was set between 0.1-0.5, and mechanism strengths were given numbers between 1-5, representing the 'a' parameter: protein concentrations functions.

Model Results

  1. Die-out timeframe: 300-750 minutes.
  2. Working timeframe: 190-600 minutes.
  3. Generation range: 13-32 generations.

What does a simulation of one of these combinations of protein leakage and mechanism strength look like? In each example below, 5 experiments out of the 30 are marked in grey, an average curve is red, and standard deviation curves are black.


Strongest mechanism, Weakest leakage rate

Example 1: Strongest mechanism, weakest leakage rate. All experiments end with the entire population dying out. Cells experience a log phase of around 90 minutes, a maximum population of around 230 is reached, and the variance is small.


Intermediate combination. Mechanism strength=4, Leakage rate=0.4

Example 2: Intermediate parameters. Here, the stochasticity of the system takes an important role. Some of the experiments are trending down at the end of the timeframe, while others are trending up. Variance is much bigger. Log phase is about 200 minutes. A maximum population of around 410 is reached.

In order to compare all of the experiments, 3D surfaces projecting the data of all combinations were made. In each one of figures 4-6, a specific parameter was chosen to be compared between the combinations. The coefficient of variation (CV) of the 30 repeats was computed and is represented by the surface’s color:


Figure 3: Z axis shows the time it took for 95% of the population to die out. X and Y are the values for the tested parameters, the bacteria symbol represents a non- reliable combination where less than 90% of the experiments died out within the model time frame, and the rest lived on.

The graph omits combinations in which the weakest mechanism strength 1 was used, because those combinations never died out. This parameter is robust and has a clear limit of around 2 as defined in the model.

The leakage is a more modular parameter, giving delicate tuning for high mechanism strengths, thus strengthening our hypothesis that leakage can provide nodularity for the system.

A time frame of around 300-750 minutes is seen, as the z-differential between the lowest combination and the highest one.


Figure 4: Z axis is the number of generations reached after 1000 minutes of simulation, in combinations that died out - a high number wasn’t reached


A generation frame of 13-32 generations within the effective combinations is seen in this model. This effects the chance of the system to transfer genes horizontally, and the chance for a mutation in the mechanism to occur.


Figure 5: Z axis is the time it took until population reached a 50% mark from initial conditions.


This parameter, the measurement of the half-life, can give an indication for the possible working time for the system. A time frame of 190-600 minutes has been calculated.

While inducer concentration usually doesn’t provide adequate control of the system because of the promoter’s bimodality[1] (all or nothing behavior when induced), some control can be achieved from different inducer combinations [2]. Therefore it is possible to induce P.A.S.E-containing bacteria to express only a partial percentage of the maximal protein expression capability, allowing for control of the overall population lifetime.

View our code.



Continue the journey: read about our Human Practice.



References

[1] [1] J Biotechnol. 2007 Feb 1;128(2):362-75. Epub 2006 Oct 17. Cell population heterogeneity in expression of a gene-switching network with fluorescent markers of different half-lives. Portle S, Causey TB, Wolf K, Bennett GN, San KY, Mantzaris N. View Source
[2] R. Lutz, H. Bujard. Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements. Nucleic Acids Research 25(6), 1203–1210 (1997). View Source