Team:NTU-Singapore/Parts/Parameter Estimation using Characterization Results
From 2008.igem.org
Lalala8585 (Talk | contribs) (→Links) |
Lalala8585 (Talk | contribs) |
||
(10 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: | + | <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 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 | + | [[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=== | ||
Line 87: | Line 87: | ||
<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. | <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: | + | [[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| | + | [[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 101: | 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. | ||
Line 114: | Line 113: | ||
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 129: | 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
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
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
Nmin = | Minimum number of species in the ecosystem |
c = | Adjustment factor |
Logistic Growth Model with Varying Carrying Capacity
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 , and assuming mass action kinetics, the differential equation for C is
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).
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.
By giving our reactor model a signal to "grow", we can try to simulate this logistic growth output as observed in our experiments.
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.
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