"
Team:Freiburg/Project/modeling
From 2013.igem.org
| Line 207: | Line 207: | ||
</tr> | </tr> | ||
</tbody></table> | </tbody></table> | ||
| - | |||
| - | |||
| - | |||
| Line 267: | Line 264: | ||
For more detailed information refer our <a id="link" href="https://2013.igem.org/Team:Freiburg/Notebook/modeling">modeling notebook</a>. </p> | For more detailed information refer our <a id="link" href="https://2013.igem.org/Team:Freiburg/Notebook/modeling">modeling notebook</a>. </p> | ||
| - | + | ||
<p id="h3"> Fitting Procedure and Results </p> | <p id="h3"> Fitting Procedure and Results </p> | ||
| + | <p id="h4"> dCas-VP16 </p> | ||
<p>Assuming the given ODE and using the fminsearch-function implemented in matlab with various initial parameter vectors the fitting process results in one optimal parameter composition (Table.1). During the fitting process not only the kinetic parameters were assumed as variable, but also the initial concentrations, because our first measured point is at time point 6 h. However the model is written in a way that the change between flexible and variable parameters is easily to perform, because of an additional vector (qfit). This vector contains boolean values depending on wether the parameter is fixed or flexible during the fitting process. Moreover another parameter is required for adjusting the absolute values of the duplicates.</p> | <p>Assuming the given ODE and using the fminsearch-function implemented in matlab with various initial parameter vectors the fitting process results in one optimal parameter composition (Table.1). During the fitting process not only the kinetic parameters were assumed as variable, but also the initial concentrations, because our first measured point is at time point 6 h. However the model is written in a way that the change between flexible and variable parameters is easily to perform, because of an additional vector (qfit). This vector contains boolean values depending on wether the parameter is fixed or flexible during the fitting process. Moreover another parameter is required for adjusting the absolute values of the duplicates.</p> | ||
| Line 301: | Line 299: | ||
</tbody></table> | </tbody></table> | ||
</div> | </div> | ||
| - | |||
| + | |||
| + | <p id="h4"> dCas-KRAB </p> | ||
| + | |||
| + | <p></p> | ||
| + | <p></p> | ||
| + | |||
| + | |||
| + | <div> | ||
| + | <table class="gelpic"> | ||
| + | <tbody><tr> | ||
| + | <td> <img id="bild" src=https://static.igem.org/mediawiki/2013/4/47/Freiburg2013_Rep100_1.png"> </td> | ||
| + | </tr> | ||
| + | <tr> | ||
| + | <td> <b id="Fig6" >Fig. 6: Modeling Result:</b><br> Shown are the experimental results (purple square) in comparison to the model prediction values (cyan cross) for SEAP and Cas, as well as the model prediction for not measured components </td> | ||
| + | </tr> | ||
| + | </tbody></table> | ||
| + | </div> | ||
| + | |||
| + | <p> This time we also fitted the different error values. <span id="refer"> <a href="#Fig7"> Fig.7 </a></span>.</p> | ||
| + | |||
| + | <div> | ||
| + | <table class="gelpic"> | ||
| + | <tbody><tr> | ||
| + | <td> <img id="bild" src="https://static.igem.org/mediawiki/2013/d/d9/Rep100_Fehler1.png"> </td> | ||
| + | </tr> | ||
| + | <tr> | ||
| + | <td> <b id="Fig7" >Fig. 7: Different error values plotted in increasing order.</b><br> </td> | ||
| + | </tr> | ||
| + | </tbody></table> | ||
| + | </div> | ||
| + | |||
| + | --> | ||
<p><small><b>Sources</b> <br> | <p><small><b>Sources</b> <br> | ||
</small></p> | </small></p> | ||
Revision as of 15:48, 3 October 2013
Modeling our dCAS
Introduction
We used a thermodynamic approach to model and characterize our system. It is based on various ordinary differential equations (ODE) that describe the behaviour of our network. Due to the limited measurment possibilities and the unwritten law, that you should at least measure half of the number of components of your network we started by using a small network with a limited amount of different components.
The Networks
1. dCAS-VP16
Our network includes four different components dCas-VP16, a RNA complex (tracr/cr RNA), a RNA-dCas-VP16 complex and the Secreted alkaline phosphatase (SEAP). DCas-VP16 binds the RNA-complex and the whole complex binds the DNA, which leads to the production of SEAP.
 |
| Fig. 1: Transcriptional Activation via dCAS-VP16: The dCAS-VP16 fusion protein is guided to the desired DNA sequence by a co-expressed crRNA and tracrRNA. The binding of the gene recognition complex leads to an expression of SEAP. |
Setting up the ODE
According to the graphical reaction network the ODE can be set up.
Cas9 is constitutively expressed by the CBh promoter and degraded proportional to the current concentration. It is used to build the DNA recognition complex and produced during complex decay.
The RNA-complex is build linearly. The production constant can be seen as production constant of the lower expressed RNA, because this expression limits the complex building. It is assumed that the RNA is degraded after DNA recognition complex decay and therefore the complex decay does not lead to more RNA.
The DNA recognition complex is built, when Cas9 and RNA meets and degraded proportional to the current DNA recognition complex concentration.
There is a leaky SEAP production and one that depends on the current concentration of the Cas9/RNA Complex. This dependency is assumed to follow the Monod-kinetic. Because of the long half time (T2 > 500 h) of SEAP we can neglect the SEAP decay [2, 3].
| The parameters are: | |
|---|---|
| k1: linear production rate of Cas9 k2: Cas9 degradation rate k3: tracr/crRNA production rate k4: tracr/crRNA degradation rate k5: gene recognition complex building rate |
k6: cr/trRNA /Cas9 degradation rate k7: SEAPs leaky production rate k8: Complex dependent SEAP production rate k9: |
Sources
The Code Files
