Team:Groningen/Modelling/Characterization
From 2009.igem.org
(Difference between revisions)
(Ignoring most experiments and focussing on just a few variables.) |
(Enabled time<Infinity measurements and adjusted the starting values a bit.) |
||
Line 99: | Line 99: | ||
|id="K7sol"| | |id="K7sol"| | ||
|- | |- | ||
- | | | + | |tauBbetaB |
|id="tauBbeta4"| | |id="tauBbeta4"| | ||
|id="tauBbeta4cur"| | |id="tauBbeta4cur"| | ||
Line 111: | Line 111: | ||
|id="tauBsol"| | |id="tauBsol"| | ||
|- | |- | ||
- | | | + | |betaB |
|id="beta4"| | |id="beta4"| | ||
|id="beta4cur"| | |id="beta4cur"| | ||
Line 123: | Line 123: | ||
|id="tauRsol"| | |id="tauRsol"| | ||
|- | |- | ||
- | | | + | |betaRN |
|id="beta1"| | |id="beta1"| | ||
|id="beta1cur"| | |id="beta1cur"| | ||
Line 218: | Line 218: | ||
{constants:{Vc:0.006666667,Vs:(1-0.006666667),pro:0,ars2T:0,proK:1.6605e-9},AsT:20e-6, | {constants:{Vc:0.006666667,Vs:(1-0.006666667),pro:0,ars2T:0,proK:1.6605e-9},AsT:20e-6, | ||
data:{AsinT:[2.25e-4,3.47e-4,4.19e-4,3.93e-4,4.19e-4,4.82e-4,4.82e-4,4.95e-4], | data:{AsinT:[2.25e-4,3.47e-4,4.19e-4,3.93e-4,4.19e-4,4.82e-4,4.82e-4,4.95e-4], | ||
- | time:[582,1212,1890,2514,3144,3828,4260,6036]}}, | + | time:[582,1212,1890,2514,3144,3828,4260,6036]}},*/ |
pSB1A2con: // concentration mode this is our first icps measerment wild type | pSB1A2con: // concentration mode this is our first icps measerment wild type | ||
- | {constants:{Vc:0.000808081,Vs:(1-0.000808081),pro:0,ars2T:0},time: | + | {constants:{Vc:0.000808081,Vs:(1-0.000808081),pro:0,ars2T:0},time:3600, |
data:{AsinT:[207.0208222e-6,229.0443139e-6,493.3262146e-6,585.8248799e-6], | data:{AsinT:[207.0208222e-6,229.0443139e-6,493.3262146e-6,585.8248799e-6], | ||
- | AsT:[10e-6,20e-6,50e-6,100e-6]}}, | + | AsT:[10e-6,20e-6,50e-6,100e-6]}}, |
pSB1A2time: // concentration mode this is our first icps measerment wild type | pSB1A2time: // concentration mode this is our first icps measerment wild type | ||
{constants:{Vc:0.002320346,Vs:(1-0.002320346),pro:0,ars2T:0},AsT:10e-6, | {constants:{Vc:0.002320346,Vs:(1-0.002320346),pro:0,ars2T:0},AsT:10e-6, | ||
Line 255: | Line 255: | ||
beta5: function(v){return v.beta5;}};*/ | beta5: function(v){return v.beta5;}};*/ | ||
- | var varsToMutate = [/*'v5_K5','v5',*/'k8_K7','k8','tauBbeta4','beta4',' | + | var varsToMutate = [/*'v5_K5','v5',*/'k8_K7','k8','tauBbeta4','beta4', |
+ | 'tauRbeta1_tauBbeta4','beta1_beta4'/*,'tauFbetaF','betaF', | ||
'tauKbetaK','betaK','tauGbeta5','beta5','tauF','betaF','tauK','betaK','tauG','beta5','ars2T'*/]; | 'tauKbetaK','betaK','tauGbeta5','beta5','tauF','betaF','tauK','betaK','tauG','beta5','ars2T'*/]; | ||
var mutateFuncs = {//v5: function(v){return v.v5;}, | var mutateFuncs = {//v5: function(v){return v.v5;}, | ||
Line 262: | Line 263: | ||
K7: function(v){return v.k8/v.k8_K7;}, | K7: function(v){return v.k8/v.k8_K7;}, | ||
tauB: function(v){return v.tauBbeta4/v.beta4;}, | tauB: function(v){return v.tauBbeta4/v.beta4;}, | ||
- | |||
beta4: function(v){return v.beta4;}, | beta4: function(v){return v.beta4;}, | ||
+ | tauR: function(v){return v.tauRbeta1_tauBbeta4*v.tauBbeta4/(v.beta4*v.beta1_beta4);}, | ||
+ | beta1: function(v){return v.beta4*v.beta1_beta4;}/*, | ||
//tauF: function(v){return v.tauF;}, | //tauF: function(v){return v.tauF;}, | ||
//betaF: function(v){return v.betaF;}, | //betaF: function(v){return v.betaF;}, | ||
Line 271: | Line 273: | ||
//ars2T: function(v){return v.ars2T;}, | //ars2T: function(v){return v.ars2T;}, | ||
//beta5: function(v){return v.beta5;}, | //beta5: function(v){return v.beta5;}, | ||
- | |||
tauF: function(v){return v.tauFbetaF/v.betaF;}, | tauF: function(v){return v.tauFbetaF/v.betaF;}, | ||
betaF: function(v){return v.betaF;}, | betaF: function(v){return v.betaF;}, | ||
Line 297: | Line 298: | ||
xt = simulate(x0,e[i].data.time,function(t,d){return arsenicModelGradient(nc,d);}); | xt = simulate(x0,e[i].data.time,function(t,d){return arsenicModelGradient(nc,d);}); | ||
- | // Sum (squares of) errors, divided by the average value | + | // Sum of (squares of) errors, divided by the average value |
var curcost = 0, n = 0; | var curcost = 0, n = 0; | ||
for(var xn in e[i].data) { | for(var xn in e[i].data) { | ||
Line 334: | Line 335: | ||
} | } | ||
- | // Sum (squares of) errors, divided by the average value | + | // Sum of (squares of) errors, divided by the average value |
for(var xn in e[i].data) { | for(var xn in e[i].data) { | ||
if (xn=='AsT') continue; | if (xn=='AsT') continue; | ||
curcost += Math.pow((e[i].data[xn][j]-xt[xn])/avgv[xn],2); | curcost += Math.pow((e[i].data[xn][j]-xt[xn])/avgv[xn],2); | ||
+ | n++; | ||
+ | } | ||
+ | } | ||
+ | cost += Math.sqrt(curcost/n); // Compute the square root of the average of the squares (RMS) | ||
+ | weight++; | ||
+ | e[i].solution.cost = Math.sqrt(curcost/n); | ||
+ | } else if (!isNaN(e[i].time)) { // Vary AsT, with t = e[i].time | ||
+ | var avgv = {}; | ||
+ | for(var xn in e[i].data) { | ||
+ | avgv[xn] = 0; | ||
+ | for(var j in e[i].data[xn]) avgv[xn] += e[i].data[xn][j]; | ||
+ | avgv[xn] /= e[i].data[xn].length; | ||
+ | } | ||
+ | e[i].solution = {'xt':{'AsT':[]}}; | ||
+ | var curcost = 0, n = 0; | ||
+ | for(var j in e[i].data.AsT) { | ||
+ | // Simulate | ||
+ | x0 = arsenicModelInitialization(nc,e[i].data.AsT[j]); | ||
+ | xt = simulate(x0,e[i].time,function(t,d){return arsenicModelGradient(nc,d);}); | ||
+ | e[i].solution.xt.AsT[j] = e[i].data.AsT[j]; | ||
+ | |||
+ | // Fill solution | ||
+ | for(var xn in xt) { | ||
+ | if (e[i].solution.xt[xn]==undefined) e[i].solution.xt[xn] = []; | ||
+ | e[i].solution.xt[xn][j] = xt[xn][xt[xn].length-1]; | ||
+ | } | ||
+ | |||
+ | // Sum of (squares of) errors, divided by the average value | ||
+ | for(var xn in e[i].data) { | ||
+ | if (xn=='AsT') continue; | ||
+ | curcost += Math.pow((e[i].data[xn][j]-xt[xn][xt[xn].length-1])/avgv[xn],2); | ||
n++; | n++; | ||
} | } | ||
Line 382: | Line 414: | ||
//cSol.v5_K5 = orgC.v5/orgC.K5; | //cSol.v5_K5 = orgC.v5/orgC.K5; | ||
//cSol.v5 = orgC.v5; | //cSol.v5 = orgC.v5; | ||
- | cSol.k8 = | + | cSol.k8 = 10; |
- | cSol.k8_K7 = | + | cSol.k8_K7 = 2e5; |
+ | cSol.tauBbeta4 = 55; | ||
+ | cSol.beta4 = 18; | ||
+ | cSol.tauRbeta1_tauBbeta4 = 400; | ||
+ | cSol.beta1_beta4 = 2; | ||
// cSol.tauBbeta4 = 180000; | // cSol.tauBbeta4 = 180000; | ||
// cSol.tauB = 180; | // cSol.tauB = 180; | ||
Line 405: | Line 441: | ||
//c.v5_K5 = orgC.v5/orgC.K5; | //c.v5_K5 = orgC.v5/orgC.K5; | ||
//c.v5 = orgC.v5; | //c.v5 = orgC.v5; | ||
- | c.k8 = | + | c.k8 = 10; |
- | c.k8_K7 = | + | c.k8_K7 = 2e5; |
+ | c.tauBbeta4 = 55; | ||
+ | c.beta4 = 18; | ||
+ | c.tauRbeta1_tauBbeta4 = 400; | ||
+ | c.beta1_beta4 = 2; | ||
// cSol.tauBbeta4 = 180000; | // cSol.tauBbeta4 = 180000; | ||
// cSol.tauB = 180; | // cSol.tauB = 180; | ||
Line 477: | Line 517: | ||
//document.getElementById('pArsRRFPtimeGraph').refresh(); | //document.getElementById('pArsRRFPtimeGraph').refresh(); | ||
//document.getElementById('Kostal2004fig3AGraph').refresh(); | //document.getElementById('Kostal2004fig3AGraph').refresh(); | ||
- | + | document.getElementById('pSB1A2conGraph').refresh(); | |
//document.getElementById('pArsRRFPconGraph').refresh(); | //document.getElementById('pArsRRFPconGraph').refresh(); | ||
} | } | ||
Line 596: | Line 636: | ||
|{{graph|Team:Groningen/Graphs/Characterization/pArsRRFPtime|id=pArsRRFPtimeGraph}} | |{{graph|Team:Groningen/Graphs/Characterization/pArsRRFPtime|id=pArsRRFPtimeGraph}} | ||
|{{graph|Team:Groningen/Graphs/Characterization/Kostal2004fig3A|id=Kostal2004fig3AGraph}} | |{{graph|Team:Groningen/Graphs/Characterization/Kostal2004fig3A|id=Kostal2004fig3AGraph}} | ||
- | |- | + | |- --> |
|{{graph|Team:Groningen/Graphs/Characterization/pSB1A2con|id=pSB1A2conGraph}} | |{{graph|Team:Groningen/Graphs/Characterization/pSB1A2con|id=pSB1A2conGraph}} | ||
- | |{{graph|Team:Groningen/Graphs/Characterization/pArsRRFPcon|id=pArsRRFPconGraph}}--> | + | <!--|{{graph|Team:Groningen/Graphs/Characterization/pArsRRFPcon|id=pArsRRFPconGraph}}--> |
|} | |} | ||
<!-- Don't forget to update the refreshGraphs function above! --> | <!-- Don't forget to update the refreshGraphs function above! --> |
Revision as of 14:56, 18 October 2009
[http://2009.igem.org/Team:Groningen http://2009.igem.org/wiki/images/f/f1/Igemhomelogo.png]
|
---|
- Modelling
- DetailedModel
- Characterization
- Downloads
TODO: Talk about the devices we have and in what way we want to characterize them.
Uptake measurements
Time (min) | ||||||
---|---|---|---|---|---|---|
0 | 10 | 20 | 40 | 60 | ||
As(III)exT(0) (µM) | 0 | x | ||||
10 | x | x | x | x | x | |
20 | x | |||||
50 | x | |||||
100 | x |
To efficiently look at both time and concentration dependent processes we took samples as in the table on the right. Below we list all results, which have been used for fitting all necessary parameters.
TODO: List results. Take conversion from nmol/mg and mg/ml to µM and Vc/Vs into account.
best | cur | gradient | solved | |
---|---|---|---|---|
v5/K5 | ||||
v5 | ||||
K5 | ||||
k8/K7 | ||||
k8 | ||||
K7 | ||||
tauBbetaB | ||||
tauB | ||||
betaB | ||||
tauR | ||||
betaRN | ||||
tauFbetaF | ||||
tauF | ||||
betaF | ||||
tauKbetaK | ||||
tauK | ||||
betaK | ||||
tauGbeta5 | ||||
tauG | ||||
beta5 | ||||
ars2T | ||||
E |
|
|
Loading graph...
| Loading graph...
|