Team:NTU-Singapore/Modelling/ODE

From 2008.igem.org

(Difference between revisions)
(→Regulated Transcription)
 
(7 intermediate revisions not shown)
Line 5: Line 5:
<div id="maincontent" style="margin-top:200px;">
<div id="maincontent" style="margin-top:200px;">
-
=ODEs used in modeling=
+
='''ODEs used in modeling'''=
In this section, we will have a short discussion on the ODEs used for Deterministic Modeling in our iGem project.  
In this section, we will have a short discussion on the ODEs used for Deterministic Modeling in our iGem project.  
Line 14: Line 14:
-
The following equations shows the break down of the different equations that will be used in this modeling exercise. By understanding this section, it would make the understanding of the system of ODEs used in our iGEm modeling exercise.
+
The following equations show the break down of the different equations that will be used in this modeling exercise. By understanding this section, it would make the understanding of the system of ODEs used in our iGEm modeling exercise.
 +
[[Image:Modeling cartoon.gif|center]]
---------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------
Line 101: Line 102:
The system can be viewed as two parts. The first part comprises of lactose induced production of colicin E7 and the immunity protein. The second part comprises of a detection mechanism that produces the lysis protein upon the detection of both Iron ions and Ai-2 ( Autoinducer 2).
The system can be viewed as two parts. The first part comprises of lactose induced production of colicin E7 and the immunity protein. The second part comprises of a detection mechanism that produces the lysis protein upon the detection of both Iron ions and Ai-2 ( Autoinducer 2).
==Lactose controlled production of E7 + Imm==
==Lactose controlled production of E7 + Imm==
-
===Variables===
+
'''Variables'''
-
*LacI = A
+
{|border="0"|
-
*Lactose =B
+
|LacI =
-
*E7 = C
+
|A
 +
|-
 +
|Lactose =
 +
|B
 +
|-
 +
|E7 =
 +
|C
 +
|}
-
# LacI production
+
===LacI production===
-
#* Transcription of LacI gene mRNA
+
* '''Transcription of LacI gene mRNA'''
[[Image:Lac1.JPG‎ |200px|Lac1]]<br>
[[Image:Lac1.JPG‎ |200px|Lac1]]<br>
-
[MA] : LacI mRNA concentration<br>
+
{|border="0"|
-
k1A : kinetic constant of transcription<br>
+
|[MA]:
-
U : system input<br>
+
|LacI mRNA concentration
-
d1A : degradation constant for mRNA<br>
+
|-
 +
|k1A:
 +
|kinetic constant of transcription
 +
|-
 +
|U:
 +
|system input
 +
|-
 +
|d1A:
 +
|degradation constant for mRNA
 +
|}
-
#* Translation of LacI Protein
+
* '''Translation of LacI Protein'''
[[Image:Lac2.JPG‎ |200px|Lac2]]<br>
[[Image:Lac2.JPG‎ |200px|Lac2]]<br>
-
[PA] : Protein concentration<br>
+
{|border="0"|
-
k2A : kinetic constant of translation<br>
+
|[PA]:
-
d2A : degradation constant for Protein<br>
+
|Protein concentration
-
 
+
|-
-
#Complex formation between LacI and Lactose
+
|k2A:
 +
|kinetic constant of translation
 +
|-
 +
|d2A:
 +
|degradation constant for Protein
 +
|}
 +
===Complex formation between LacI and Lactose===
[[Image:Complex1.JPG‎ |300px|Complex]]<br>
[[Image:Complex1.JPG‎ |300px|Complex]]<br>
-
[ComAB] : Complex of LacI and Lactose<br>
+
{|border="0"|
-
K3AB : Rate constant of complex formation<br>
+
|[ComAB] :
-
[PA]o : Protein concentration of LacI at the time of Lactose addition<br>
+
|Complex of LacI and Lactose
-
[B]o : Initial concentration of Lactose<br>
+
|-
 +
|K3AB :
 +
|Rate constant of complex formation
 +
|-
 +
|[PA]o :
 +
|Protein concentration of LacI at the time of Lactose addition
 +
|-
 +
|[B]o :
 +
|Initial concentration of Lactose
 +
|}
-
#E7 production
+
===E7 production===
-
#* Transcription of E7 gene mRNA
+
* '''Transcription of E7 gene mRNA'''
[[Image:E71.JPG‎ |300px|E71]]<br>
[[Image:E71.JPG‎ |300px|E71]]<br>
-
[MC]: E7 mRNA Concentration<br>
+
{|border="0"|
-
k1C: kinetic constant of transcription<br>
+
|[MC]:
-
a : constitutive portion , 0<a<1<br>
+
|E7 mRNA Concentration
-
KC: Hill constant<br>
+
|-
-
nC: Hill coefficient<br>
+
|k1C:
-
d1C: degradation constant for mRNA<br>
+
|kinetic constant of transcription
 +
|-
 +
|a :
 +
|constitutive portion , 0<a<1
 +
|-
 +
|KC:
 +
|Hill constant
 +
|-
 +
|nC:
 +
|Hill coefficient
 +
|-
 +
|d1C:
 +
|degradation constant for mRNA
 +
|}
-
#* Translation of E7 Protein
+
* '''Translation of E7 Protein'''
[[Image:E72.JPG‎ |200px|E72]]<br>
[[Image:E72.JPG‎ |200px|E72]]<br>
-
[PC] : Protein concentration<br>
+
{|border="0"|
-
k2C : kinetic constant of translation<br>
+
|[PC] :
-
d2C : degradation constant for Protein<br>
+
|Protein concentration
 +
|-
 +
|k2C :
 +
|kinetic constant of translation
 +
|-
 +
|d2C :
 +
|degradation constant for Protein
 +
|}
-
==Iron and Ai2 controlled production of Lysis==
+
==Variables of Iron and Ai2 controlled production of Lysis==
-
===Variables===
+
 
-
Variables
+
'''Variables'''
-
#Ai-2 : A
+
{|border="0"|
-
#ai-2-phos : B
+
|#Ai-2:
-
#LsrR : C
+
|A
-
#SupD derivatives : D
+
|-
-
#T7ptag : E
+
|#ai-2-phos:
-
#Lysis : F
+
|B
 +
|-
 +
|#LsrR:
 +
|C
 +
|-
 +
|#SupD derivatives:
 +
|D
 +
|-
 +
|#T7ptag:
 +
|E
 +
|-
 +
|#Lysis:
 +
|F
 +
|}
-
#Phosphorylation of ai-2
+
===Phosphorylation of Ai-2===
[[Image:Phos.JPG‎ |200px|Phos]]<br>
[[Image:Phos.JPG‎ |200px|Phos]]<br>
-
A: Ai-2<br>
+
{|border="0"|
-
B: Ai-2-phos<br>
+
|A:
-
kPF: Forward Phosphorylation<br>
+
|Ai-2
-
S: Concentration of Protein Kinase<br>
+
|-
-
kPB: Backward dephosphorylation<br>
+
|B:
 +
|Ai-2-phos
 +
|-
 +
|kPF:
 +
|Forward Phosphorylation
 +
|-
 +
|S:
 +
|Concentration of Protein Kinase
 +
|-
 +
|kPB:
 +
|Backward dephosphorylation
 +
|}
-
# LsrR production
+
===LsrR production===
-
#* Transcription of LsrR gene mRNA
+
*'''Transcription of LsrR gene mRNA'''
[[Image:LsrR1.JPG‎ |200px|LsrR1]]<br>
[[Image:LsrR1.JPG‎ |200px|LsrR1]]<br>
-
[MC]: LacI mRNA concentration<br>
+
{|border="0"|
-
k1C : kinetic constant of transcription<br>
+
|[MC]:
-
U : system input<br>
+
|LacI mRNA concentration
-
d1C: degradation constant for mRNA,br>
+
|-
 +
|k1C:
 +
|kinetic constant of transcription
 +
|-
 +
|U:
 +
|system input
 +
|-
 +
|d1C:
 +
|degradation constant for mRNA
 +
|}
-
#* Translation of LsrR protein
+
* '''Translation of LsrR protein'''
[[Image:LsrR2.JPG‎ |200px|LsrR2]]<br>
[[Image:LsrR2.JPG‎ |200px|LsrR2]]<br>
-
[PC] : Protein concentration<br>
+
{|border="0"|
-
k2C : kinetic constant of translation<br>
+
|[PC] :
-
d2C : degradation constant for Protein<br>
+
|Protein concentration
 +
|-
 +
|k2C:
 +
|kinetic constant of translation
 +
|-
 +
|d2C:
 +
|degradation constant for Protein
 +
|}
-
#Complex formation between Ai-2-phos and LsrR
+
===Complex formation between Ai-2-phos and LsrR===
-
[[Image:Complex2.JPG‎ |200px|Complex2]]<br>
+
[[Image:Complex2.JPG‎ |350px|Complex2]]<br>
-
[ComAB] : Complex of LsrR and Ai-2-phos<br>
+
{|border="0"|
-
k3BC : Rate constant of complex formation<br>
+
|[ComAB] :
-
[PC]o: Protein concentration of LsrR at the time of start of complex formation <br>
+
|Complex of LsrR and Ai-2-phos
-
[B]o: Initial concentration of Ai-2-phos<br>
+
|-
 +
|k3BC :
 +
|Rate constant of complex formation
 +
|-
 +
|[PC]o:
 +
|Protein concentration of LsrR at the time of start of complex formation
 +
|-
 +
|[B]o:
 +
|Initial concentration of Ai-2-phos
 +
|}
-
The hard part here is to determine what is the initial concentration of the protein LsrR and Ai-2-phos. A way out of this mess is to simplify the equation even more.  
+
The hard part here is to determine what is the initial concentration of the protein LsrR and Ai-2-phos. A way out of this mess is to simplify the equation even more.
-
#AND GATE odes
+
===AND GATE odes===
-
#*Transcription of SupD gene mRNA
+
*'''Transcription of SupD gene mRNA'''
[[Image:SupD1.JPG‎ |300px|SupD]]<br>
[[Image:SupD1.JPG‎ |300px|SupD]]<br>
-
[MD]: t7 mRNA Concentration <br>
+
{|border="0"|
-
k1D: kinetic constant of transcription<br>
+
|[MD]:
-
a : constitutive portion , 0<a<1<br>
+
|t7 mRNA Concentration  
-
KD: Hill constant<br>
+
|-
-
nD: Hill coefficient<br>
+
|k1D:
-
d1D: degradation constant for mRNA<br>
+
|kinetic constant of transcription
 +
|a :
 +
|constitutive portion , 0<a<1
 +
|-
 +
|KD:
 +
|Hill constant
 +
|-
 +
|nD:
 +
|Hill coefficient
 +
|-
 +
|d1D:
 +
|degradation constant for mRNA
 +
|}
-
#*Transcription of T7ptag gene mRNA
+
*'''Transcription of T7ptag gene mRNA'''
[[Image:T71.JPG‎ |300px|T71]]<br>
[[Image:T71.JPG‎ |300px|T71]]<br>
-
[ME]: supD mRNA Concentration <br>
+
{|border="0"|
-
k1E: kinetic constant of transcription<br>
+
|[ME]:
-
a : constitutive portion , 0<a<1<br>
+
|supD mRNA Concentration
-
KE: Hill constant<br>
+
|-
-
nE: Hill coefficient<br>
+
|k1E:
-
d1E: degradation constant for mRNA<br>
+
|kinetic constant of transcription
-
 
+
|-
-
#*Complex formation between tRNA and mRNA
+
|a :
-
[[Image:Complex3.JPG‎ |300px|Complex3]]<br>
+
|constitutive portion , 0<a<1
-
[ComDE] : Complex of tRNA and mRNA t7<br>
+
|-
-
k3DE : Rate constant of complex formation<br>
+
|KE:
-
[MD]o : concentration of tRNA at the time of start of complex formation<br>
+
|Hill constant
-
[ME]o : Initial concentration of mRNA t7<br>
+
|-
 +
|nE:
 +
|Hill coefficient
 +
|-
 +
|d1E:
 +
|degradation constant for mRNA
 +
|}
 +
*'''Complex formation between tRNA and mRNA'''
 +
[[Image:Complex3.JPG‎ |350px|Complex3]]<br>
 +
{|border="0"|
 +
|[ComDE] :
 +
|Complex of tRNA and mRNA t7
 +
|-
 +
|k3DE :
 +
|Rate constant of complex formation
 +
|-
 +
|[MD]o :
 +
|concentration of tRNA at the time of start of complex formation
 +
|-
 +
|[ME]o :
 +
|Initial concentration of mRNA t7
 +
|}
-
#*Translation of t7 protein
+
*'''Translation of t7 protein'''
[[Image:T72.JPG‎ |200px|Complex3]]<br>
[[Image:T72.JPG‎ |200px|Complex3]]<br>
-
[PF] : Protein concentration<br>
+
{|border="0"|
-
k2F : kinetic constant of translation<br>
+
|[PF] :
-
d2F : degradation constant for Protein<br>
+
|Protein concentration
 +
|-
 +
|k2F :
 +
|kinetic constant of translation
 +
|-
 +
|d2F :
 +
|degradation constant for Protein
 +
|}
-
#Lysis
+
===Lysis production===
-
#*Transcription of Lysis gene mRNA
+
*'''Transcription of Lysis gene mRNA'''
[[Image:Lysis1.JPG‎ |300px|Lysis]]<br>
[[Image:Lysis1.JPG‎ |300px|Lysis]]<br>
-
[MF]: Lysis mRNA Concentration <br>
+
{|border="0"|
-
k1F: kinetic constant of transcription<br>
+
|[MF]:
-
a : constitutive portion , 0<a<1<br>
+
|Lysis mRNA Concentration
-
KF: Hill constant<br>
+
|-
-
nF: Hill coefficient<br>
+
|k1F:
-
d1F: degradation constant for mRNA<br>
+
|kinetic constant of transcription
 +
|a :
 +
|constitutive portion , 0<a<1
 +
|-
 +
|KF:
 +
|Hill constant
 +
|-
 +
|nF:
 +
|Hill coefficient
 +
|-
 +
|d1F:
 +
|degradation constant for mRNA
 +
|}
-
 
+
*'''Translation of Lysis protein'''
-
#*Translation of Lysis protein
+
[[Image:Lysis2.JPG‎ |200px|Lysis2]]<br>
[[Image:Lysis2.JPG‎ |200px|Lysis2]]<br>
-
[PF] : Protein concentration<br>
+
{|border="0"|
-
k2F : kinetic constant of translation<br>
+
|[PF] :
-
d2F : degradation constant for Protein<br>
+
|Protein concentration
 +
|-
 +
|k2F :
 +
|kinetic constant of translation
 +
|-
 +
|d2F :
 +
|degradation constant for Protein
 +
|}
 +
-----------------------------------------------------------------------------------------------------------------
 +
<br><br>
 +
<html>
 +
<script language=Javascript1.2>
 +
<!--
 +
var tags_before_clock = "<b>It is now "
 +
var tags_middle_clock = "on"
 +
var tags_after_clock  = "</b>"
-
-----------------------------------------------------------------------------------------------------------------
+
if(navigator.appName == "Netscape") {
 +
document.write('<layer id="clock"></layer><br>');
 +
}
 +
 
 +
if (navigator.appVersion.indexOf("MSIE") != -1){
 +
document.write('<span id="clock"></span>');
 +
}
 +
 
 +
DaysofWeek = new Array()
 +
  DaysofWeek[0]="Sunday"
 +
  DaysofWeek[1]="Monday"
 +
  DaysofWeek[2]="Tuesday"
 +
  DaysofWeek[3]="Wednesday"
 +
  DaysofWeek[4]="Thursday"
 +
  DaysofWeek[5]="Friday"
 +
  DaysofWeek[6]="Saturday"
 +
 
 +
Months = new Array()
 +
  Months[0]="January"
 +
  Months[1]="February"
 +
  Months[2]="March"
 +
  Months[3]="April"
 +
  Months[4]="May"
 +
  Months[5]="June"
 +
  Months[6]="July"
 +
  Months[7]="August"
 +
  Months[8]="September"
 +
  Months[9]="October"
 +
  Months[10]="November"
 +
  Months[11]="December"
 +
 
 +
function upclock(){
 +
var dte = new Date();
 +
var hrs = dte.getHours();
 +
var min = dte.getMinutes();
 +
var sec = dte.getSeconds();
 +
var day = DaysofWeek[dte.getDay()]
 +
var date = dte.getDate()
 +
var month = Months[dte.getMonth()]
 +
var year = dte.getFullYear()
 +
 
 +
var col = ":";
 +
var spc = " ";
 +
var com = ",";
 +
var apm;
 +
 
 +
if (date == 1 || date == 21 || date == 31)
 +
  {ender = "<sup>st</sup>"}
 +
else
 +
if (date == 2 || date == 22)
 +
  {ender = "<sup>nd</sup>"}
 +
else
 +
if (date == 3 || date == 23)
 +
  {ender = "<sup>rd</sup>"}
 +
 
 +
else
 +
  {ender = "<sup>th</sup>"}
 +
 
 +
if (12 < hrs) {
 +
apm="<font size='-1'>pm</font>";
 +
hrs-=12;
 +
}
 +
 
 +
else {
 +
apm="<font size='-1'>am</font>";
 +
}
 +
 
 +
if (hrs == 0) hrs=12;
 +
if (hrs<=9) hrs="0"+hrs;
 +
if (min<=9) min="0"+min;
 +
if (sec<=9) sec="0"+sec;
 +
 
 +
if(navigator.appName == "Netscape") {
 +
document.clock.document.write(tags_before_clock+hrs+col+min+col+sec+apm+spc+tags_middle_clock+spc+day+com+spc+date+ender+spc+month+com+spc+year+tags_after_clock);
 +
document.clock.document.close();
 +
}
 +
 
 +
if (navigator.appVersion.indexOf("MSIE") != -1){
 +
clock.innerHTML = tags_before_clock+hrs+col+min+col+sec+apm+spc+tags_middle_clock+spc+day+com+spc+date+ender+spc+month+com+spc+year+tags_after_clock;
 +
}
 +
}
 +
 
 +
setInterval("upclock()",1000);
 +
//-->
 +
</script>
 +
</html>

Latest revision as of 05:42, 27 October 2008

Contents

ODEs used in modeling

In this section, we will have a short discussion on the ODEs used for Deterministic Modeling in our iGem project.

For the modeling exercise, we use mainly simple linear ODEs to describe the system we have at hand. These ODEs are used to describe processes such as constitutive transcription/ translation etc. Nonlinear equations are used for more complex processes such as regulated transcription. The choice of these equations are in hope that they can describe the system to a relative degree of accuracy and yet be simple enough to be solved. This is very much an engineering compromise between a rigorous description and getting an answer that is good enough which is required.

Our systems are based on the assumptions that Lumped parameter models are sufficient to describe them.
Lumped parameter models are those in which spatial variations are ignored and the dependent variables are uniform throughout the whole system. The only variations are those with respect to time leading to mathematical description in terms of the ODEs used.


The following equations show the break down of the different equations that will be used in this modeling exercise. By understanding this section, it would make the understanding of the system of ODEs used in our iGEm modeling exercise.

Modeling cartoon.gif

Constant synthesis & Linear Synthesis

Linear Synthesis

  • Simple ode to describe constant synthesis
  • Gives an explicit analytical solution
  • Unique solution once a IC is posed

Linear Degradation

Linear Degradation

  • Rate of degradation is proportional to how much of the molecule is present
  • Gives an explicit analytical solution
  • Constant half life

Simple Forward Reaction

Complex Formation

[C] : Complex
kc : Rate constant of complex formation

This equation ignores the fact that dissociation of the complex occurs. We can do so if the dissociation is much slower than the formation.

  • Single solvable equation for the unknown C
  • Simple, unique solution available with I.C

Phosphorylation and Dephosphorylation

Assumptions:

  • Linear kinetic rate laws apply only if XT is much less than the Michaelis constants of both kinase and phosphotase

Phosphorylation and Dephosphorylation

XT : total cost of X protein in phosphorylated and unphosphorylated form
S : protein kinase concentration
k2 : accounts for protein phosphotase
  • Modeled after simple linear kinetics
  • Gives a hyperbolic signal response curve when X plotted vs S

Regulated Transcription

Regulated1

[P]: Protein Formed
µ: Repression, µ=0;
Activation, µ=1
K: Hill Constant Value of input that gives 50% response
n: Hill coefficient Slope of signal-response curve at this input signal
d: degradation of protein
k1: basal gene expression
k: signal-dependent gene expression
a: correlation between k1 and k, 0<a<1

Regulated2
This ODE attempts to capture characteristics of the mRNA dynamics
For our modeling, all our detection systems activates some form of transcription. Therefore µ=1 in all cases for our modeling exercise. Regulated3

ODE system used in model

The system can be viewed as two parts. The first part comprises of lactose induced production of colicin E7 and the immunity protein. The second part comprises of a detection mechanism that produces the lysis protein upon the detection of both Iron ions and Ai-2 ( Autoinducer 2).

Lactose controlled production of E7 + Imm

Variables

LacI = A
Lactose = B
E7 = C

LacI production

  • Transcription of LacI gene mRNA

Lac1

[MA]: LacI mRNA concentration
k1A: kinetic constant of transcription
U: system input
d1A: degradation constant for mRNA
  • Translation of LacI Protein

Lac2

[PA]: Protein concentration
k2A: kinetic constant of translation
d2A: degradation constant for Protein

Complex formation between LacI and Lactose

Complex

[ComAB] : Complex of LacI and Lactose
K3AB : Rate constant of complex formation
[PA]o : Protein concentration of LacI at the time of Lactose addition
[B]o : Initial concentration of Lactose

E7 production

  • Transcription of E7 gene mRNA

E71

[MC]: E7 mRNA Concentration
k1C: kinetic constant of transcription
a : constitutive portion , 0<a<1
KC: Hill constant
nC: Hill coefficient
d1C: degradation constant for mRNA
  • Translation of E7 Protein

E72

[PC] : Protein concentration
k2C : kinetic constant of translation
d2C : degradation constant for Protein

Variables of Iron and Ai2 controlled production of Lysis

Variables

#Ai-2: A
#ai-2-phos: B
#LsrR: C
#SupD derivatives: D
#T7ptag: E
#Lysis: F

Phosphorylation of Ai-2

Phos

A: Ai-2
B: Ai-2-phos
kPF: Forward Phosphorylation
S: Concentration of Protein Kinase
kPB: Backward dephosphorylation

LsrR production

  • Transcription of LsrR gene mRNA

LsrR1

[MC]: LacI mRNA concentration
k1C: kinetic constant of transcription
U: system input
d1C: degradation constant for mRNA
  • Translation of LsrR protein

LsrR2

[PC] : Protein concentration
k2C: kinetic constant of translation
d2C: degradation constant for Protein

Complex formation between Ai-2-phos and LsrR

Complex2

[ComAB] : Complex of LsrR and Ai-2-phos
k3BC : Rate constant of complex formation
[PC]o: Protein concentration of LsrR at the time of start of complex formation
[B]o: Initial concentration of Ai-2-phos

The hard part here is to determine what is the initial concentration of the protein LsrR and Ai-2-phos. A way out of this mess is to simplify the equation even more.

AND GATE odes

  • Transcription of SupD gene mRNA

SupD

[MD]: t7 mRNA Concentration
k1D: kinetic constant of transcription a : constitutive portion , 0<a<1
KD: Hill constant
nD: Hill coefficient
d1D: degradation constant for mRNA
  • Transcription of T7ptag gene mRNA

T71

[ME]: supD mRNA Concentration
k1E: kinetic constant of transcription
a : constitutive portion , 0<a<1
KE: Hill constant
nE: Hill coefficient
d1E: degradation constant for mRNA
  • Complex formation between tRNA and mRNA

Complex3

[ComDE] : Complex of tRNA and mRNA t7
k3DE : Rate constant of complex formation
[MD]o : concentration of tRNA at the time of start of complex formation
[ME]o : Initial concentration of mRNA t7
  • Translation of t7 protein

Complex3

[PF] : Protein concentration
k2F : kinetic constant of translation
d2F : degradation constant for Protein

Lysis production

  • Transcription of Lysis gene mRNA

Lysis

[MF]: Lysis mRNA Concentration
k1F: kinetic constant of transcription a : constitutive portion , 0<a<1
KF: Hill constant
nF: Hill coefficient
d1F: degradation constant for mRNA
  • Translation of Lysis protein

Lysis2

[PF] : Protein concentration
k2F : kinetic constant of translation
d2F : degradation constant for Protein



Retrieved from "http://2008.igem.org/Team:NTU-Singapore/Modelling/ODE"