Team:NTU-Singapore/Modeling

From 2008.igem.org

(Difference between revisions)
(E7 prodcution system)
 
(24 intermediate revisions not shown)
Line 3: Line 3:
<div id="header">{{User:Greenbear/sandbox/header}}</div>
<div id="header">{{User:Greenbear/sandbox/header}}</div>
-
<div id="maincontent" style="margin-top:200px;">
+
<div id="maincontent" style="margin-top:100px;">
-
 
+
 +
[[Image:Blackboard.jpg|center|400px]]
=Introduction=
=Introduction=
-
The use of models to describe synthetic biology has its merits. Synthetic biology investigates the use of different biological parts to put together and assemble devices that carry out specific functions. Good mathematical models to describe each part would greatly help not only in the characterization of a part but also facilitate the use of the part by other people when they choose to use the part within their devices or systems.
 
-
Simulations based on modeling can give a first insight on how the system would turn out and provide a rough guide of the system’s behaviour.
+
The primary objective of this section is to assist us in developing an understanding of the dynamic behaviour of the systems that we would create in our iGEM project. A requirement for assessing the dynamic behaviour is a time-dependent mathematical model of the bio-chemical processes that take place.  
-
=System=
 
-
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).
 
-
=ODEs used in modeling=
+
-----------------------------------------------------------------------------------------------------------------
-
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
+
Dictionary Definitions : Model
 +
Model is derived from the Latin word ''modus'', which stands for ''measure''. As a noun it means, " a small representation of a planned or exiting object" (Webster's New World Dictionary)
 +
"A mathematical or physical system, obeying certain specified conditions, whose behaviour is used to understand a physical, biological, or social system to which it is analogous in some way" (McGraw-Hill dictionary of Scientific and Technical Terms.)
 +
-----------------------------------------------------------------------------------------------------------------
-
==Constant synthesis & Linear Synthesis==
+
The use of models to describe synthetic biology has its merits. Synthetic biology investigates the use of different biological parts to put together and assemble devices that carry out specific functions. Good mathematical models to describe each part would greatly help not only in the characterization of a part but also facilitate the use of the part by other people when they choose to use the part within their devices or systems.
-
*Simple ode to describe constant synthesis
+
-
*Gives an explicit analytical solution
+
-
*Unique solution once a IC is posed
+
-
==Linear Degradation==
+
Simulations based on modeling can give a first insight on how the system would turn out and provide a rough guide of the system’s behaviour.
-
*Rate of degradation is proportional to how much of the molecule is present
+
-
*Gives an explicit analytical solution
+
-
*Constant half life
+
-
==Simple Forward Reaction==
+
We hope to use a control engineering approach in this project, and we hope that by understanding the transient behaviours of the systems can we better design our system to suit their purpose.  
-
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
+
The models that are created and found in the subsequent sections can still be improved by comparing the modeling results against actual wet lab results. SIMULINK has proved to be an immensely useful tool for Deterministic modeling and solving ODEs. Cellware had its own limitations when dealing with Stochastic models. As the system got more complicated, the shorter the time span the program was able to simulate and produce useful results.  
-
*Simple, unique solution available with I.C
+
-
==Phosphorylation and Dephosphorylation==
+
We had seen from our models that our AND gate system should have the shorter SupD gene attached to the Fe promoter. This was because in our assumptions, Fe ions would be ion lower concentration in the intestinal environment compared to Ai-2 and to have a similar response rate at the AND gate, the choice of AND gate inputs is critical.
-
Assumptions:
+
-
*Linear kintic rate laws apply only if XT is much less than the Michaelis constants of both kinase and phosphotase
+
 +
We also realised that differing inputs of Lactose or Ai-2 would have an impact on the system via the models and this would be explored in greater detail in our wet lab experiments to see if the models are predicting the outcome correctly. <br>
 +
<br>
 +
<br>
-
*Modeled after simple linear kinetics
+
''"All models are wrong; Some are useful"'' <br>
-
*Gives a hyperbolic signal response curve when X plotted vs S
+
'''George P. Box.''' <br>
 +
<br><br>
 +
<html>
 +
<script language=Javascript1.2>
 +
<!--
-
==Regulated Transcription==
+
var tags_before_clock = "<b>It is now "
 +
var tags_middle_clock = "on"
 +
var tags_after_clock  = "</b>"
-
=ODE system used in model=
+
if(navigator.appName == "Netscape") {
-
==Lactose controlled production of E7 + Imm==
+
document.write('<layer id="clock"></layer><br>');
-
===Variables===
+
}
-
#LacI = A
+
-
#Lactose =B
+
-
#E7 = C
+
-
==Iron and Ai2 controlled production of Lysis==
+
if (navigator.appVersion.indexOf("MSIE") != -1){
-
===Variables===
+
document.write('<span id="clock"></span>');
-
Variables
+
}
-
#Ai-2 : A
+
-
#ai-2-phos : B
+
-
#LsrR : C
+
-
#SupD derivatives : D
+
-
#T7ptag : E
+
-
#Lysis : F
+
-
=Parameter Estimation=
+
DaysofWeek = new Array()
-
Estimation of different parameters
+
  DaysofWeek[0]="Sunday"
-
# Transcription         : 70nt/s
+
  DaysofWeek[1]="Monday"
-
# Translation : 40aa/s
+
  DaysofWeek[2]="Tuesday"
-
# Number of Essential Genes : 297
+
  DaysofWeek[3]="Wednesday"
-
# Number of mRNA per cell : 4000
+
  DaysofWeek[4]="Thursday"
-
# Average mRNA half life : 5min
+
  DaysofWeek[5]="Friday"
-
# Average mRNA length         : 1100
+
  DaysofWeek[6]="Saturday"
-
Assumptions
+
Months = new Array()
-
#Rate of transcription is dependent on length of gene
+
  Months[0]="January"
-
#Number of amino acids is 1/3 of the number of nucleotides in a gene
+
  Months[1]="February"
-
#Rate of Translation is dependent on number of nucleotides
+
  Months[2]="March"
-
#For each gene mRNA = 10 at steady state
+
  Months[3]="April"
-
#Rate of degradation of average mRNA = 1100/ 5 min
+
  Months[4]="May"
-
#Rate of degradation of protein is equivalent to time for cell division i.e. 40 min
+
  Months[5]="June"
 +
  Months[6]="July"
 +
  Months[7]="August"
 +
  Months[8]="September"
 +
  Months[9]="October"
 +
  Months[10]="November"
 +
  Months[11]="December"
-
[[Image:Parameter_Estimation.JPG|450px|Results]]
+
function upclock(){
-
==E7 production system==
+
var dte = new Date();
-
{| align="center" border="1"
+
var hrs = dte.getHours();
-
|Type
+
var min = dte.getMinutes();
-
|Parameter
+
var sec = dte.getSeconds();
-
|Values
+
var day = DaysofWeek[dte.getDay()]
-
|Comments
+
var date = dte.getDate()
-
|-
+
var month = Months[dte.getMonth()]
-
|Transcription Rate of Lac I gene
+
var year = dte.getFullYear()
-
|k1A
+
-
|21
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Transcription Rate of E7 + Imm gene
+
-
|k1C
+
-
|2.470588
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Degradation Rate of Lac I mRNA
+
-
|d1A
+
-
|0.76246
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Transcription Rate of E7 + Imm mRNA
+
-
|d1C
+
-
|0.0897
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Translation Rate of Lac I mRNA
+
-
|k2A
+
-
|36
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Translation Rate of E7 + Imm mRNA
+
-
|k2C
+
-
|4.23539
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Protein Degradation Rate
+
-
|d2A,d2C
+
-
|0.03465
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Hill coefficeint for E7 + Imm
+
-
|nC
+
-
|1
+
-
|This is obtained on the assumption that one Repressor <br>Protein binds to one Lactose molecule complex
+
-
|-
+
-
|Dissociation Constant for E7 + Imm
+
-
|KC
+
-
|0.8
+
-
|[3]
+
-
|-
+
-
|Constitutive Portion for E7 + Imm
+
-
|aC
+
-
|0.5
+
-
|Estimate since a is between 0 and 1<br> Implication that Lactose may not be <br> a very strongly Regulated Promoter
+
-
|-
+
-
|Complex Formation Rate Between <br> Lac I Repressor and Lactose
+
-
|k3AB
+
-
|1
+
-
|Estimate
+
-
|}
+
-
==E7 production system==
+
var col = ":";
-
{| align="center" border="1"
+
var spc = " ";
-
|Type
+
var com = ",";
-
|Parameter
+
var apm;
-
|Values
+
-
|Comments
+
-
|-
+
-
|Phosphorylation Rate of Ai-2
+
-
|kPF
+
-
|1
+
-
|Estimate
+
-
|-
+
-
|Dephosphorylation Rate of Ai-2
+
-
|kPB
+
-
|1
+
-
|Estimate
+
-
|-
+
-
|Transcription Rate of LsrR gene
+
-
|k1C
+
-
|4.402517
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Transcription Rate of SupD gene
+
-
|k1D
+
-
|46.667
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Transcription Rate of t7pTag gene
+
-
|k1E
+
-
|1.5556
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Transcription Rate of Lysis gene
+
-
|k1F
+
-
|28
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Degradation Rate of LsrR mRNA
+
-
|d1C
+
-
|0.159845
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Degradation Rate of SupD mRNA
+
-
|d1D
+
-
|1.694359
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Degradation Rate of t7pTag mRNA
+
-
|d1E
+
-
|0.056478
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Degradation Rate of Lysis mRNA
+
-
|d1F
+
-
|1.0166
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Translation Rate of LsrR Protein
+
-
|k2C
+
-
|7.54716
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Translation Rate of t7 pTag Protein
+
-
|k2E
+
-
|2.6667
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Translation Rate of Lysis Protein
+
-
|k2F
+
-
|48
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Protein Degradation Rate
+
-
|d2C,d2E,d2F
+
-
|0.03465
+
-
|Made using Earlier assumptions
+
-
|-
+
-
|Hill coefficeint for SupD
+
-
|nD
+
-
|50
+
-
|Trial And Error
+
-
|-
+
-
|Hill coefficeint for t7 pTag
+
-
|nE
+
-
|1
+
-
|Trial And Error
+
-
|-
+
-
|Hill coefficeint for Lysis
+
-
|nF
+
-
|1
+
-
|Estimate
+
-
|-
+
-
|Dissociation Constant for SupD
+
-
|KD
+
-
|15
+
-
|[5]
+
-
|-
+
-
|Dissociation Constant for t7 pTag
+
-
|KE
+
-
|1
+
-
|Estimate
+
-
|-
+
-
|Dissociation Constant for Lysis
+
-
|KF
+
-
|0.8
+
-
|Estimate
+
-
|-
+
-
|Constitutive Portion for SupD
+
-
|aD
+
-
|0.01
+
-
|Estimate
+
-
|-
+
-
|Constitutive Portion for t7 pTag
+
-
|aE
+
-
|0.01
+
-
|Estimate
+
-
|-
+
-
|Constitutive Portion for Lysis
+
-
|aF
+
-
|0.0001
+
-
|Estimate
+
-
|-
+
-
|Complex Formation Rate Between <br> LsrR Repressor and Ai-2-Phosphorylated
+
-
|k3BC
+
-
|0.01
+
-
|Estimate
+
-
|-
+
-
|Complex Formation Rate Between <br> SupD tRNA and t7 pTag mRNA
+
-
|k3DE
+
-
|0.000000001
+
-
|Estimate
+
-
|-
+
-
|Amount of Protein Kinase
+
-
|S
+
-
|28.747
+
-
|From Earlier Assumptions
+
-
|-
+
-
|}
+
-
=Results of Modeling=
+
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>"}
-
==E7 production system==
+
else
-
[[Image:Varying_lactose_input.jpg|850px| E7 Production Modeling Results]]
+
  {ender = "<sup>th</sup>"}
-
Modeling of the System shows that Lactose induction is essential to produce E7 and a variation of Lactose input can result in different yields of E7
+
if (12 < hrs) {
 +
apm="<font size='-1'>pm</font>";
 +
hrs-=12;
 +
}
-
==Lysis Production system==
+
else {
-
[[Image:Varying_Lysis.jpg|850px|Lysis Prodcution Modeling Results]]
+
apm="<font size='-1'>am</font>";
 +
}
-
Here we define 1 as a certain threshold e.g. (<250 µM) that when the lysis protein reaches it, lysis in the cell definitely occurs.
+
if (hrs == 0) hrs=12;
-
Although we wish for an absolute AND gate, where 0 will have no lysis production at all, simulation on biological systems shows that such results are impossible. Both addition Fe ions and Ai-2 alone would induce a certain level of lysis production. However when both are available, the lysis protein production would be higher.
+
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();
 +
}
-
In the presence of Ai-2 alone, Lysis will most likely still occur albeit at a slower rate compared to situation when both Fe and Ai-2 are present.
+
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;
 +
}
 +
}
-
With the addition of iron and Ai-2 together the rate of lysis production is still significantly much higher compared to Ai-2 alone.
+
setInterval("upclock()",1000);
 +
//-->
 +
</script>
 +
</html>
-
Therefore, by giving the logical output ‘1’ as a suitable threshold value of Lysis protein production higher than that of Ai-2 induction alone, we would still be able to obtain an AND gate based on our definition.
 
 +
<!---End of main content--->
-
 
-
 
-
 
-
</div>
 
-
<!---End of main content--->
 
-
</div>
 
<!---End of body--->
<!---End of body--->

Latest revision as of 05:42, 27 October 2008

Blackboard.jpg

Introduction

The primary objective of this section is to assist us in developing an understanding of the dynamic behaviour of the systems that we would create in our iGEM project. A requirement for assessing the dynamic behaviour is a time-dependent mathematical model of the bio-chemical processes that take place.



Dictionary Definitions : Model

Model is derived from the Latin word modus, which stands for measure. As a noun it means, " a small representation of a planned or exiting object" (Webster's New World Dictionary)

"A mathematical or physical system, obeying certain specified conditions, whose behaviour is used to understand a physical, biological, or social system to which it is analogous in some way" (McGraw-Hill dictionary of Scientific and Technical Terms.)


The use of models to describe synthetic biology has its merits. Synthetic biology investigates the use of different biological parts to put together and assemble devices that carry out specific functions. Good mathematical models to describe each part would greatly help not only in the characterization of a part but also facilitate the use of the part by other people when they choose to use the part within their devices or systems.

Simulations based on modeling can give a first insight on how the system would turn out and provide a rough guide of the system’s behaviour.

We hope to use a control engineering approach in this project, and we hope that by understanding the transient behaviours of the systems can we better design our system to suit their purpose.

The models that are created and found in the subsequent sections can still be improved by comparing the modeling results against actual wet lab results. SIMULINK has proved to be an immensely useful tool for Deterministic modeling and solving ODEs. Cellware had its own limitations when dealing with Stochastic models. As the system got more complicated, the shorter the time span the program was able to simulate and produce useful results.

We had seen from our models that our AND gate system should have the shorter SupD gene attached to the Fe promoter. This was because in our assumptions, Fe ions would be ion lower concentration in the intestinal environment compared to Ai-2 and to have a similar response rate at the AND gate, the choice of AND gate inputs is critical.

We also realised that differing inputs of Lactose or Ai-2 would have an impact on the system via the models and this would be explored in greater detail in our wet lab experiments to see if the models are predicting the outcome correctly.


"All models are wrong; Some are useful"
George P. Box.