Team:IBB Pune/Modeling

From 2009.igem.org

(Difference between revisions)
(→Assumptions)
 
(46 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>
-
<br><br><br>
+
<html>
-
<html><p><span style="font-weight:bold; font-size:200%; color:#F07211;">Project 1- A Model Of The Snowdrift Game</span></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]]
-
This is a Model of the Snowdrift game
+
<br>
-
This model assumes use of two different strains as co-operators and defectors.
+
 
 +
 
 +
 
 +
<p>This is a Model of the Snowdrift game. This model assumes use of two different strains as '''co-operators''' and '''defectors'''.</p>
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 expts will be the output of the previous cycle of expts.
+
<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>
 +
 
 +
*<p>There is no intracellular lactose metabolism (only extracellular). </p>
 +
 
 +
*<p>There is no lag in enzyme production and secretion.</p>
 +
 
 +
*<p>Rate of degradation of enzyme is zero</p>
 +
<span style="font-weight:bold; font-size:125%; color:#006600;"><p>Model</p></span>
 +
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:
 +
<br>
 +
L = L - L * Ec * k2                                  .... (5)
 +
 
-
'''Artificial assumptions'''
 
-
<p>Glucose is consumed by all cells. Galactose is also consumed at the same rate  Gc mg/cell/min/ml.
+
The net Glucose present in the medium at time ''''t'''' is given by:
-
The metabolic benefit due to glucose and galactose is same.So effectively each lactose molecule gives rise to 2 glucose molecules</p>
+
<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>
-
<p>There is no intracellular lactose metabolism (only extracellular). </p>
+
[[Image:Pic2.GIF |center|400px|thumbnail|The population distribution at k2=1, Gc = 0.0001 mg/cell/min, R=0.99, c=0.001]]
-
<p>There is no lag in enzyme production and secretion.</p>
+
<br>
-
<p>Rate of degradation of enzyme is zero</p>
+
[[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 ]]

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"