Team:BGU Israel/Model1.html

From 2013.igem.org

(Difference between revisions)
(Created page with "{{Team:BGU_Israel/header}} <html> <!-- MAIN --> <div id="main"> <!-- pcontent --> <div id="pcontent"> <!-- title --> <div id="page-title"> <span class="...")
Line 19: Line 19:
         <div id="page-content">
         <div id="page-content">
   
   
-
                   <p>Pellentesque habitant morbi tristique senectus et netus et malesuada fames ac turpis egestas. Vestibulum tortor quam, feugiat vitae, ultricies eget, tempor sit amet, ante. Donec eu libero sit amet quam egestas semper. Aenean ultricies mi vitae est. Mauris placerat eleifend leo.<span class="pullquote-right">THIS IS A PULL QUOTE RIGHT, BGU_Israel AMET</span> Quisque sit amet est et sapien ullamcorper pharetra.  Ut felis. Praesent dapibus, neque id cursus faucibus,  Vestibulum erat wisi, condimentum sed, commodo vitae, ornare sit amet, wisi. Aenean fermentum, elit eget tincidunt condimentum, eros ipsum rutrum orci, sagittis tempus lacus enim ac dui. Donec non enim in turpis pulvinar facilisis. Ut felis. Praesent dapibus, neque id cursus faucibus, tortor neque egestas augue, eu vulputate magna eros eu erat. Aliquam erat volutpat. Nam dui mi, tincidunt quis, accumsan porttitor, facilisis luctus, metus</p>
+
                   <p>It is possible to describe the behavior of a bacterial population by looking at a specific generation k of bacteria and examining its inter-division rate and death rate:</br></br><img src="https://static.igem.org/mediawiki/2013/f/fb/BGU_eq1.png"/> </br>
 +
Where &beta; is inter-division rate and µ is death rate. </br>
 +
A deterministic equation to describe such a population will look like this: </br>
 +
</br>
-
 
-
<img src="https://static.igem.org/mediawiki/2013/f/fb/BGU_eq1.png"/> </br>
 
<img src="https://static.igem.org/mediawiki/2013/f/f7/BGU_eq2.png"/> </br>
<img src="https://static.igem.org/mediawiki/2013/f/f7/BGU_eq2.png"/> </br>
 +
This deterministic approach assumes that the doubling time equals the inter-division time, meaning
<img src="https://static.igem.org/mediawiki/2013/c/c3/BGU_eq3.png"/> </br>
<img src="https://static.igem.org/mediawiki/2013/c/c3/BGU_eq3.png"/> </br>
 +
The stochastic approach for the same issue suggests that there are different inter-division and death time distributions that will fit a specific population behavior, and not just an average time. </br> </br>
 +
 +
One can replace the use of the rates &beta;  and &mu;  with the use of inter-division and death time distributions instead <b>[1]</b>.</br></br>
 +
 +
There is evidence <b>[2][1]</b> and it is fair to assume that the distributions for cell division and death times are similar to the Gamma distribution. Right after an individual cell has divided, in order for a second proliferation to happen, a set of events in the cell cycle need to occur, therefore proliferation probability is low right after the division, is rising as time advances, and is falling again asymmetrically as time goes by after avg. cell cycle time. Same can be assumed for death time distribution because the cell cycle can be considered as a sensitive time with high death possibility. </br></br>
 +
 +
As described in the overview, we defined our problematic parameters that we are modeling as the leakage of the system, and the strength of the mechanism influencing the cell. When a bacterium is released and no more induction occurs, concentration will behave like this: </br>
 +
</br>
 +
<img src="https://static.igem.org/mediawiki/2013/f/f5/BGU_eq4.png"/> </br>
<img src="https://static.igem.org/mediawiki/2013/f/f5/BGU_eq4.png"/> </br>
 +
Where k is generation (k=1, [protein]=100%).</br>
 +
 +
In this model we chose distributions for inter-division and death times such that they will describe the behavior. Those distributions were initially configured to represent a regular genetically modified E.coli, the initial setting corresponded to a P.A.S.E 1 with 100% CI concentration, or a P.A.S.E 2 with 100% tyrosine synthetase concentration (the state of the organism just after its release when it is fully induced). The algorithm, which is based on a Gillespie algorithm, advanced each individual in the population with time according to those distributions (a decision is being made each iteration for each cell if it will die or proliferate). With each update, individual protein concentration was updated as well. </br>
 +
</br>
 +
 +
 +
So how did we introduce protein concentration differentiation and mechanism strength to the system? </br>
 +
</br>
 +
 +
The Gamma(a,b) distribution is defined by its parameters- ‘a’ is shape (dimensionless) , ‘b’ is scale (minutes in our case). The Mean of the distribution equals a*b and represents the average inter division time, it was initially set for a*b=35 min, similar to a genetically modified E.coli. Control of the distribution means control over the shape parameter ‘a’ , therefore we chose the death time shape parameter to be linearly depended on the concentration of protein as follows: </br>
 +
</br>
 +
 +
a=slope*[protein]+n. </br>
 +
</br>
 +
 +
Mechanism strength was simply introduced by changing the ‘slope’ and ‘n’ parameters to represent a stronger influence of [protein] on ‘a’. </br>
 +
</br>
 +
 +
Protein leakage was changed by changing the function that updates the concentration each iteration. </br>
 +
</br>
 +
 +
25 combinations of leakages and mechanism strengths were simulated with repeats, and means and deviation of different parameters were analyzed. </br>
 +
</br>
 +
 +
<img src="https://static.igem.org/mediawiki/2013/e/ee/BGU_eq5.png"/> </br>
<img src="https://static.igem.org/mediawiki/2013/e/ee/BGU_eq5.png"/> </br>
 
 
 +
</p>
 +
 +
 +
 +
   
   

Revision as of 16:06, 28 September 2013

BGU_Israel

A Stochastic Birth-Death Model

It is possible to describe the behavior of a bacterial population by looking at a specific generation k of bacteria and examining its inter-division rate and death rate:


Where β is inter-division rate and µ is death rate.
A deterministic equation to describe such a population will look like this:


This deterministic approach assumes that the doubling time equals the inter-division time, meaning
The stochastic approach for the same issue suggests that there are different inter-division and death time distributions that will fit a specific population behavior, and not just an average time.

One can replace the use of the rates β and μ with the use of inter-division and death time distributions instead [1].

There is evidence [2][1] and it is fair to assume that the distributions for cell division and death times are similar to the Gamma distribution. Right after an individual cell has divided, in order for a second proliferation to happen, a set of events in the cell cycle need to occur, therefore proliferation probability is low right after the division, is rising as time advances, and is falling again asymmetrically as time goes by after avg. cell cycle time. Same can be assumed for death time distribution because the cell cycle can be considered as a sensitive time with high death possibility.

As described in the overview, we defined our problematic parameters that we are modeling as the leakage of the system, and the strength of the mechanism influencing the cell. When a bacterium is released and no more induction occurs, concentration will behave like this:


Where k is generation (k=1, [protein]=100%).
In this model we chose distributions for inter-division and death times such that they will describe the behavior. Those distributions were initially configured to represent a regular genetically modified E.coli, the initial setting corresponded to a P.A.S.E 1 with 100% CI concentration, or a P.A.S.E 2 with 100% tyrosine synthetase concentration (the state of the organism just after its release when it is fully induced). The algorithm, which is based on a Gillespie algorithm, advanced each individual in the population with time according to those distributions (a decision is being made each iteration for each cell if it will die or proliferate). With each update, individual protein concentration was updated as well.

So how did we introduce protein concentration differentiation and mechanism strength to the system?

The Gamma(a,b) distribution is defined by its parameters- ‘a’ is shape (dimensionless) , ‘b’ is scale (minutes in our case). The Mean of the distribution equals a*b and represents the average inter division time, it was initially set for a*b=35 min, similar to a genetically modified E.coli. Control of the distribution means control over the shape parameter ‘a’ , therefore we chose the death time shape parameter to be linearly depended on the concentration of protein as follows:

a=slope*[protein]+n.

Mechanism strength was simply introduced by changing the ‘slope’ and ‘n’ parameters to represent a stronger influence of [protein] on ‘a’.

Protein leakage was changed by changing the function that updates the concentration each iteration.

25 combinations of leakages and mechanism strengths were simulated with repeats, and means and deviation of different parameters were analyzed.