Team:IBB Pune/Modeling

From 2009.igem.org

(Difference between revisions)
(→Model)
 
(42 intermediate revisions not shown)
Line 2: Line 2:
{{Team:IBB_Pune/menu}}
{{Team:IBB_Pune/menu}}
-
<html><p><span style="font-weight:bold; font-size:200%; color:#6600FF;">Modelling</span></p></html>
+
<html>
 +
<br>
 +
<p><span style="font-weight:bold; font-size:200%; color:#006600;">Modelling</span></p></html>
 +
<p>The aim of the modelling aspect of our project is to develop simulations that would enable us to predict the behavior of the export tag with the Snowdrift Game as the backdrop. </p>
 +
 
 +
<html>
 +
<p><span style="font-weight:bold; font-size:200%; color:#006600;">
 +
<br>
 +
Project 1- A Model Of The Snowdrift Game</span>
 +
</p></html>
 +
 
 +
<br><br>
 +
[[Image:Modelling (1).JPG|center|800 px|thumbnail|Model]]
 +
 
 +
<br>
-
<br><br><br>
 
-
<html><p><span style="font-weight:bold; font-size:200%; color:#F07211;">Project 1- A Model Of The Snowdrift Game</span></p></html>
 
Line 11: Line 23:
Co-operators constitutively secrete/produce b-gal. Defectors cannot secrete/produce B-gal.
Co-operators constitutively secrete/produce b-gal. Defectors cannot secrete/produce B-gal.
-
This model is for a Single iteration of the experiment. i.e. the culture is initiated in shaker flask and allowed to grow.
+
This model is for a Single iteration of the experiment i.e. the culture is initiated in shaker flask and allowed to grow. In order to find ESS or stable equilibria, we will have to repeatedly subculture the cells under similar conditions, by replacing the medium or supplying fresh medium. In such a setup, the input for the next cycle of experiments will be the output of the previous cycle of experiments.  
   
   
-
In order to find ESS or stable equilibria, we will have to repeatedly subculture the cells under similar conditions, by replacing the medium or supplying fresh medium. In such a setup, the input for the next cycle of experiments will be the output of the previous cycle of experiments.
+
<span style="font-weight:bold; font-size:125%; color:#006600;"><p>Assumptions</p></span>
-
+
-
===Assumptions===
+
''At time t=0;''
''At time t=0;''
-
There are 'k' co-operators and 'N-k' defectors
+
*There are ''''k'''' co-operators and ''''N-k'''' defectors
-
Medium contains 'L' mg/ml of lactose, glucose conc. (g) = 0
+
*Medium contains ''''L'''' mg/ml of lactose, glucose conc. (g) = 0
-
Culture is well mixed.
+
*Culture is well mixed.
-
Extracellular Enzyme conc (Ec) = 0 units/ml
+
*Extracellular Enzyme conc ('''Ec''') = 0 units/ml
 +
<span style="font-weight:bold; font-size:150%; color:#006600;"><p>Artificial assumptions</p></span>
 +
*<p>Glucose is consumed by all cells. Galactose is also consumed at the same rate  Gc mg/cell/min/ml.
 +
*The metabolic benefit due to glucose and galactose is same.So effectively each lactose molecule gives rise to 2 glucose molecules</p>
-
'''Artificial assumptions'''
+
*<p>There is no intracellular lactose metabolism (only extracellular). </p>
-
<p>Glucose is consumed by all cells. Galactose is also consumed at the same rate  Gc mg/cell/min/ml.
+
*<p>There is no lag in enzyme production and secretion.</p>
-
The metabolic benefit due to glucose and galactose is same.So effectively each lactose molecule gives rise to 2 glucose molecules</p>
+
-
<p>There is no intracellular lactose metabolism (only extracellular). </p>
+
*<p>Rate of degradation of enzyme is zero</p>
-
 
+
<span style="font-weight:bold; font-size:125%; color:#006600;"><p>Model</p></span>
-
<p>There is no lag in enzyme production and secretion.</p>
+
-
 
+
-
<p>Rate of degradation of enzyme is zero</p>
+
-
 
+
-
 
+
-
===Model===
+
In this Model,
In this Model,
-
1. Each co-operator secretes ''''B' units of enzyme/ min and pays cost of 'c' / min'''
+
1. Each co-operator secretes ''''B'''' units of enzyme/ min and pays cost of ''''c''''/min
-
These molecules diffuse freely through the medium
+
2. These molecules diffuse freely through the medium
-
They convert Lactose to Glucose + Galactose given by
+
3. They convert Lactose to Glucose + Galactose given by
-
Rate of conversion of lactose to glucose is given by (k2)*(Ec)*(L)  
+
   
-
g mg/ml/min
+
g = (k2)*(Ec)*(L) mg/ml/min                       .... (1)
-
Absorbed glucose confers growth rate advantage of '''r'''  
+
4. Absorbed glucose confers growth rate advantage of '''r'''  
-
Advantage conferred by Glucose r= (R)*(g)*(Gc)  
+
 +
r = (R)*(g)*(Gc)                                      .... (2)
-
After each min, the population of co-operators and defectors is updated depending upon the growth
+
After each min, the population of co-operators and defectors is updated depending upon their respective growth rates.
-
rates.
+
-
+
D(t) = D(t-1)+ r * D(t-1)                             .... (3)    Defector population (t)
-
D(t) = D(t-1)+r*D(t-1) Defector population
+
k(t) = k(t-1) + (r-c) * k(t-1)                       .... (4)    Co-operator population (t)
-
k(t) = k(t-1) + (r-c)*k(t-1) Co-operator population
+
 
 +
 
 +
 
 +
 
 +
The Lactose that remains at time ''''t'''' is given by:
 +
<br>
 +
L = L - L * Ec * k2                                  .... (5)
 +
 
 +
 
 +
 
 +
The net Glucose present in the medium at time ''''t'''' is given by:
 +
<br>
 +
G = G +(( 2* L* Ec* k2)-( N * Gc)                    .... (6)
 +
<span style="font-weight:bold; font-size:125%; color:#006600;"><p>Results</p></span>
 +
[[Image:Pic1.GIF |center|400px|thumbnail|Plot of lactose and glucose levels at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001  ]]
 +
<br>
 +
 
 +
[[Image:Pic2.GIF |center|400px|thumbnail|The population distribution at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001]]
-
===Equations===
+
<br>
-
C=C0+C*((R*G0*Gc)-cost)
+
[[Image:Pic3.gif |center|400px|thumbnail|The Cooperator to defector ratio at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001 ]]
-
D=D+D*R*Gc*G0
+
-
L=L-L*Ec*k2
+
-
G=G+((2*L*Ec*k2)-(N*Gc)
+

Latest revision as of 00:03, 22 October 2009





Modelling

The aim of the modelling aspect of our project is to develop simulations that would enable us to predict the behavior of the export tag with the Snowdrift Game as the backdrop.


Project 1- A Model Of The Snowdrift Game



Model



This is a Model of the Snowdrift game. This model assumes use of two different strains as co-operators and defectors.

Co-operators constitutively secrete/produce b-gal. Defectors cannot secrete/produce B-gal.

This model is for a Single iteration of the experiment i.e. the culture is initiated in shaker flask and allowed to grow. In order to find ESS or stable equilibria, we will have to repeatedly subculture the cells under similar conditions, by replacing the medium or supplying fresh medium. In such a setup, the input for the next cycle of experiments will be the output of the previous cycle of experiments.

Assumptions

At time t=0;

  • There are 'k' co-operators and 'N-k' defectors
  • Medium contains 'L' mg/ml of lactose, glucose conc. (g) = 0
  • Culture is well mixed.
  • Extracellular Enzyme conc (Ec) = 0 units/ml

Artificial assumptions

  • Glucose is consumed by all cells. Galactose is also consumed at the same rate Gc mg/cell/min/ml.

  • The metabolic benefit due to glucose and galactose is same.So effectively each lactose molecule gives rise to 2 glucose molecules

  • There is no intracellular lactose metabolism (only extracellular).

  • There is no lag in enzyme production and secretion.

  • Rate of degradation of enzyme is zero

Model

In this Model,

1. Each co-operator secretes 'B' units of enzyme/ min and pays cost of 'c'/min

2. These molecules diffuse freely through the medium

3. They convert Lactose to Glucose + Galactose given by 


g = (k2)*(Ec)*(L) mg/ml/min	                       .... (1)

4. Absorbed glucose confers growth rate advantage of r 

r = (R)*(g)*(Gc)                                      .... (2)

After each min, the population of co-operators and defectors is updated depending upon their respective growth rates. 

D(t) = D(t-1)+ r * D(t-1)                             .... (3)     Defector population (t)	
k(t) = k(t-1) + (r-c) * k(t-1)                        .... (4)     Co-operator population (t)



The Lactose that remains at time 't' is given by: 

L = L - L * Ec * k2                                   .... (5)


The net Glucose present in the medium at time 't' is given by: 

G = G +(( 2* L* Ec* k2)-( N * Gc)                     .... (6)

Results
Plot of lactose and glucose levels at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001


The population distribution at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001


The Cooperator to defector ratio at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001
Retrieved from "http://2009.igem.org/Team:IBB_Pune/Modeling"