Team:NTU-Singapore/Parts/Parameter Estimation using Characterization Results

From 2008.igem.org

(Difference between revisions)
(Creating a Logistic Growth Construct)
 
(15 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:180px;">
-
=Parameter Estimation=
+
='''Parameter Estimation'''=
<br>[[Image:GFP_model_sim.JPG|800px|thumb|center| GFP Production Model in SIMULINK]]<br>
<br>[[Image:GFP_model_sim.JPG|800px|thumb|center| GFP Production Model in SIMULINK]]<br>
Line 19: Line 19:
2)The model seems to hit steady state much faster than the actual situation<br>
2)The model seems to hit steady state much faster than the actual situation<br>
-
To test how the parameter estimation works out, we tried to input the RFU values as GFP protein produced, intrinsically assuming that both are directly proportional to one another. The results however are discouraging, prompting us to reconsider what went wrong. The results can be found in [https://2008.igem.org/Team:NTU-Singapore/Parts/Failed_Parameter_Estimation Failed Parameter Estimation Results] link. <br>
+
To test how the parameter estimation works out, we tried to input the RFU values as GFP protein produced, intrinsically assuming that both are directly proportional to one another. The results however are discouraging, prompting us to reconsider what went wrong. The results can be found in <span style="font-size:20px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Failed_Parameter_Estimation Failed Parameter Estimation Results]</span> link. <br>
Line 26: Line 26:
This prompts us to look into a better model such that we can first capture the essentials of the real situation before we decide to use our models with our characterisation results.
This prompts us to look into a better model such that we can first capture the essentials of the real situation before we decide to use our models with our characterisation results.
-
=Creating a Logistic Growth Construct=
+
='''Creating a Logistic Growth Construct'''=
The logistic growth model is typically applied to growth kinetics of Microbes and other natural phenomena. This phenomena is based on the fact that no population can expand exponentially forever and there will be a certain limit to the number of species based on the space that is available.  
The logistic growth model is typically applied to growth kinetics of Microbes and other natural phenomena. This phenomena is based on the fact that no population can expand exponentially forever and there will be a certain limit to the number of species based on the space that is available.  
Line 78: Line 78:
Using [[Image:Complex_reaction.JPG|100px| Complex Formation reaction]] , and assuming mass action kinetics, the differential equation for C is
Using [[Image:Complex_reaction.JPG|100px| Complex Formation reaction]] , and assuming mass action kinetics, the differential equation for C is
-
[[Image:Complex_eq1.JPG|200px|thumb|center| Complex Formation equation]]  
+
[[Image:Complex_eq1.JPG|200px| Complex Formation equation]]  
kf = forward reaction rate
kf = forward reaction rate
-
kr = reverse reaction rate  
+
kr = reverse reaction rate
===Links===
===Links===
-
Do refer to [https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter#GFP_production_system_with_logistic_growth Parameter Estimation of GFP production system] for more details on the above parameters used.
+
Do refer to
 +
<span style="font-size:20px;">[https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter#GFP_production_system_with_logistic_growth Parameter Estimation of GFP production system]</span> for more details on the above parameters used.
-
==Building the Construct==
+
=='''Building the Construct'''==
Revisiting the general idea of our modeling, we have treated each bacteria cell to work as a CSTR (Continuous stirred Reactor).<br>
Revisiting the general idea of our modeling, we have treated each bacteria cell to work as a CSTR (Continuous stirred Reactor).<br>
-
[[Image:Bactoy.JPG|150px| Bacteria]]  [[Image:Equals.jpg|100px| equal]]  [[Image:CSTR.JPG|300px| CSTR]]<br>
+
[[Image:Reactor_ntu_igem.JPG|550px|thumb|center| Bacteria Reactor]]<br>
-
 
+
However, each bacteria is itself a reactor and these "mini reactors" multiply and grow~! Therefore we want our new model to be able to account for both the production of proteins and the growth pattern of the bacteria. So now instead of a simple basal signal into our model, we would send in a signal that codes for a "logistic growth" in the basal transcription of the model.<br>
However, each bacteria is itself a reactor and these "mini reactors" multiply and grow~! Therefore we want our new model to be able to account for both the production of proteins and the growth pattern of the bacteria. So now instead of a simple basal signal into our model, we would send in a signal that codes for a "logistic growth" in the basal transcription of the model.<br>
-
[[Image:Bacteria_reactor_concept.JPG|850px| Bacteria and reactor]]  
+
[[Image:Bacteria_reactor_concept.JPG|850px|thumb|center| Bacterias and reactors Concept]]  
By giving our reactor model a signal to "grow", we can try to simulate this logistic growth output as observed in our experiments.
By giving our reactor model a signal to "grow", we can try to simulate this logistic growth output as observed in our experiments.
Line 100: Line 100:
As seen, the logistic input can still capture varying lactose inputs and this model seems to be a more representative one for the characterization results. Interestingly, as lactose input increases, the model shows that the output changes less. We intend to use these results for our characterization experience and it will be found in the next section.
As seen, the logistic input can still capture varying lactose inputs and this model seems to be a more representative one for the characterization results. Interestingly, as lactose input increases, the model shows that the output changes less. We intend to use these results for our characterization experience and it will be found in the next section.
-
=Correlation section=
+
='''Correlation section'''=
Correlations form an integral part in chemical engineering. Many phenomena observed in chemical processes have been extensively studied and relevant data collected to provide correlations between the variables of interest. This is to allow end users to be able to predict certain outcomes based on the knowledge of a specified variable if the correlation is availble. This practise is common when a process is not fully understood and a quick but efficient method to determine the process output is desired. Herein lies the aim of the next part of our project: to find a correlation between the amount of GFP produced and the fluorescence observed.
Correlations form an integral part in chemical engineering. Many phenomena observed in chemical processes have been extensively studied and relevant data collected to provide correlations between the variables of interest. This is to allow end users to be able to predict certain outcomes based on the knowledge of a specified variable if the correlation is availble. This practise is common when a process is not fully understood and a quick but efficient method to determine the process output is desired. Herein lies the aim of the next part of our project: to find a correlation between the amount of GFP produced and the fluorescence observed.
-
=Model results and RFU data=
+
=='''Model results and RFU data'''==
A simple correlation was done between the model results and the RFU collected.  
A simple correlation was done between the model results and the RFU collected.  
-
[[Image:RFU_protein_37.jpg|800px|RFU versus protein results]]
+
[[Image:RFU_protein_37.jpg|800px|thumb|center|RFU versus protein results]]
For the 10 graphs, correlations were generated for each of them using the [http://www.mathworks.com/products/curvefitting/ Curve Fitting toolbox] found in MATLAB.
For the 10 graphs, correlations were generated for each of them using the [http://www.mathworks.com/products/curvefitting/ Curve Fitting toolbox] found in MATLAB.
The results are found in each of the links below. Each curve is shown together with the prediction bounds to a 95% cofidence level and the equation of the curve.<br>
The results are found in each of the links below. Each curve is shown together with the prediction bounds to a 95% cofidence level and the equation of the curve.<br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve1 Lactose = 1mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve1 Lactose = 1mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve2 Lactose = 2mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve2 Lactose = 2mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve3 Lactose = 3mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve3 Lactose = 3mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve4 Lactose = 4mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve4 Lactose = 4mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve5 Lactose = 5mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve5 Lactose = 5mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve6 Lactose = 6mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve6 Lactose = 6mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve7 Lactose = 7mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve7 Lactose = 7mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve8 Lactose = 8mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve8 Lactose = 8mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve9 Lactose = 9mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve9 Lactose = 9mM]</span><br>
-
[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve10 Lactose = 10mM]<br>
+
<span style="font-size:15px;">[https://2008.igem.org/Team:NTU-Singapore/Parts/Curve10 Lactose = 10mM]</span><br>
The curves were fitted once a R-square value of higer than 0.99 was obtained. Although the curves could be successfully obtained there remains some questions to be asked:<br>
The curves were fitted once a R-square value of higer than 0.99 was obtained. Although the curves could be successfully obtained there remains some questions to be asked:<br>
Line 128: Line 128:
2) Will temperature changes have an impact on the results? However it is still possible that our model can account for temperature effects by tweaking the logistic growth input parameter. Recall that our model accounts for varying carrying capacity, and a different temperature could be represented by different paramter values.<br>
2) Will temperature changes have an impact on the results? However it is still possible that our model can account for temperature effects by tweaking the logistic growth input parameter. Recall that our model accounts for varying carrying capacity, and a different temperature could be represented by different paramter values.<br>
-
=References=
+
='''References'''=
1. A new logistic model for Escherichia coli growth at constant and dynamic temperatures FUJIKAWA Hiroshi ; KAI Akemi ; MOROZUMI Satoshi <br>
1. A new logistic model for Escherichia coli growth at constant and dynamic temperatures FUJIKAWA Hiroshi ; KAI Akemi ; MOROZUMI Satoshi <br>
2. www.oakenstaff.org/munin/biolmodcontrols/pdfs/populationdynamicsr1.pdf <br>
2. www.oakenstaff.org/munin/biolmodcontrols/pdfs/populationdynamicsr1.pdf <br>
 +
<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:51, 27 October 2008


Contents

Parameter Estimation


GFP Production Model in SIMULINK


GFP Production Model Results


RFU Results



In the above two figures, we have a comparison between the model generated by our very own equations versus the RFU results that we have obtained using our characterization procedures. There are certain features that are our outstanding when these two figures are compared.

1)The deterministic model has not captured any of the stochastic features the real situation presented before us.
2)The model seems to hit steady state much faster than the actual situation

To test how the parameter estimation works out, we tried to input the RFU values as GFP protein produced, intrinsically assuming that both are directly proportional to one another. The results however are discouraging, prompting us to reconsider what went wrong. The results can be found in Failed Parameter Estimation Results link.


Careful observation of the results from the characterisation and one would notice that the output takes on the form of a sigmoid curve. Indeed, one of the major flaws in our modeling exercise was that the models don't account for cell growth. Even though we used M9 medium to do our characterisation, it is impossible to assume the cells don't multiply and grow within the 12 hour time period!

This prompts us to look into a better model such that we can first capture the essentials of the real situation before we decide to use our models with our characterisation results.

Creating a Logistic Growth Construct

The logistic growth model is typically applied to growth kinetics of Microbes and other natural phenomena. This phenomena is based on the fact that no population can expand exponentially forever and there will be a certain limit to the number of species based on the space that is available.

Here we explore the results of three Logistic Growth equations and their various outputs to see how well the results turn out to be.

Logistic Growth Model

Logisitic Growth equation
Logisitic Growth Model
N = number of species in the population
t = time
Nmax = Maximum number of species in the ecosystem
r = intrinsic rate of natural increase

Logistic Growth Model with lag phase

Logisitic Growth equation with lag phase
Logisitic Growth Model
Nmin = Minimum number of species in the ecosystem
c = Adjustment factor

Logistic Growth Model with Varying Carrying Capacity

Logisitic sine Growth equation
Logisitic Growth Model

From the model, we can see that population that has an early exponential phase of growth and is not so affected by the variability in the carrying capacity. As the population reaches the limit, the population starts to oscillate at the same frequency as the carrying capacity and some of the population starts to die off.

We would use this as our logisitc equation model to add a bit more variability into our deterministic model.

Complex Formation

In an ODE model, a reaction network is expressed as a set of differential equations with one equation per chemical and with terms that represent the reactions. Our complex formation equation has ignored the fact that the reaction could be reversible. Here we add in a new complex formation equation to make our models more robust.

Using Complex Formation reaction , and assuming mass action kinetics, the differential equation for C is

Complex Formation equation

kf = forward reaction rate kr = reverse reaction rate

Links

Do refer to Parameter Estimation of GFP production system for more details on the above parameters used.

Building the Construct

Revisiting the general idea of our modeling, we have treated each bacteria cell to work as a CSTR (Continuous stirred Reactor).

Bacteria Reactor

However, each bacteria is itself a reactor and these "mini reactors" multiply and grow~! Therefore we want our new model to be able to account for both the production of proteins and the growth pattern of the bacteria. So now instead of a simple basal signal into our model, we would send in a signal that codes for a "logistic growth" in the basal transcription of the model.

Bacterias and reactors Concept

By giving our reactor model a signal to "grow", we can try to simulate this logistic growth output as observed in our experiments.


GFP Production Model in SIMULINK


GFP Production Model Result

As seen, the logistic input can still capture varying lactose inputs and this model seems to be a more representative one for the characterization results. Interestingly, as lactose input increases, the model shows that the output changes less. We intend to use these results for our characterization experience and it will be found in the next section.

Correlation section

Correlations form an integral part in chemical engineering. Many phenomena observed in chemical processes have been extensively studied and relevant data collected to provide correlations between the variables of interest. This is to allow end users to be able to predict certain outcomes based on the knowledge of a specified variable if the correlation is availble. This practise is common when a process is not fully understood and a quick but efficient method to determine the process output is desired. Herein lies the aim of the next part of our project: to find a correlation between the amount of GFP produced and the fluorescence observed.

Model results and RFU data

A simple correlation was done between the model results and the RFU collected.


RFU versus protein results

For the 10 graphs, correlations were generated for each of them using the [http://www.mathworks.com/products/curvefitting/ Curve Fitting toolbox] found in MATLAB.

The results are found in each of the links below. Each curve is shown together with the prediction bounds to a 95% cofidence level and the equation of the curve.
Lactose = 1mM
Lactose = 2mM
Lactose = 3mM
Lactose = 4mM
Lactose = 5mM
Lactose = 6mM
Lactose = 7mM
Lactose = 8mM
Lactose = 9mM
Lactose = 10mM

The curves were fitted once a R-square value of higer than 0.99 was obtained. Although the curves could be successfully obtained there remains some questions to be asked:
1) Is the model results truly representative of the real situation?
2) Will temperature changes have an impact on the results? However it is still possible that our model can account for temperature effects by tweaking the logistic growth input parameter. Recall that our model accounts for varying carrying capacity, and a different temperature could be represented by different paramter values.

References

1. A new logistic model for Escherichia coli growth at constant and dynamic temperatures FUJIKAWA Hiroshi ; KAI Akemi ; MOROZUMI Satoshi
2. www.oakenstaff.org/munin/biolmodcontrols/pdfs/populationdynamicsr1.pdf