(11 intermediate revisions by 2 users not shown) | |||
Line 84: | Line 84: | ||
synthetic biology to use modeling tools to predict behaviors and thus make more efficient the construction of | synthetic biology to use modeling tools to predict behaviors and thus make more efficient the construction of | ||
genetic circuits.</p> | genetic circuits.</p> | ||
− | <p>We have used these tools to model two of our constructs: pFadBA_34_RFP and | + | <p>We have used these tools to model two of our constructs: pFadBA_34_RFP and pFadBA_34_luxR_plux_32_RFP.</p> |
− | <p>Regarding the model for pFadBA_34_RFP:</p> | + | <p class="subapart3">Regarding the model for pFadBA_34_RFP:</p> |
<p>To analyze the response of pFadBA promoter to palmitic acid (PA) presence, we have used the RFP reporter protein to | <p>To analyze the response of pFadBA promoter to palmitic acid (PA) presence, we have used the RFP reporter protein to | ||
obtain fluorescence measurements. The issue of those experiments was that the PA has effects in the fluorescence | obtain fluorescence measurements. The issue of those experiments was that the PA has effects in the fluorescence | ||
Line 91: | Line 91: | ||
the complete transfer function of this construct. </p> | the complete transfer function of this construct. </p> | ||
<p>Nevertheless, we have made parameter optimization with weighted non-linear regression functions in order to obtain | <p>Nevertheless, we have made parameter optimization with weighted non-linear regression functions in order to obtain | ||
− | the curve that fitted the best our | + | the curve that fitted the best our less noisy data. With this, we were able to estimate the whole |
transfer function of our construct. </p> | transfer function of our construct. </p> | ||
− | <p>Regarding the model for pFadba_34_luxR_plux_32_RFP:</p> | + | <p class="subapart3">Regarding the model for pFadba_34_luxR_plux_32_RFP:</p> |
<p>In this case, even if we had the same problem that with pFadBA, since luxR/plux amplifies considerably the signal, | <p>In this case, even if we had the same problem that with pFadBA, since luxR/plux amplifies considerably the signal, | ||
we could obtain the whole curve with acceptable noise levels. Here, we also adapted our parameters to obtain a | we could obtain the whole curve with acceptable noise levels. Here, we also adapted our parameters to obtain a | ||
Line 99: | Line 99: | ||
<p>pFadBA model was created to complement and refine the experimental data, in contrast to pFadba_34_luxR_plux_32_RFP | <p>pFadBA model was created to complement and refine the experimental data, in contrast to pFadba_34_luxR_plux_32_RFP | ||
whose model was created for a different aim.The pFadba_34_luxR_plux_32_RFP system response is dependent on an | whose model was created for a different aim.The pFadba_34_luxR_plux_32_RFP system response is dependent on an | ||
− | eternal inducer: lactone. | + | eternal inducer: lactone. Having a model that predicts how much lactone is necessary for the |
modulation we want will facilitate its use in the future. </p> | modulation we want will facilitate its use in the future. </p> | ||
+ | <a href="https://2018.igem.org/Team:UPF_CRG_Barcelona/Measurement"><button class="button">Read more...</button></a> | ||
<p>The ODE for the pFadBA_34_RFP and pFadba_34_luxR_plux_32_RFP systems are given below:</p> | <p>The ODE for the pFadBA_34_RFP and pFadba_34_luxR_plux_32_RFP systems are given below:</p> | ||
+ | $$ | ||
+ | X_{FadR}= [FadR]\left(\frac{Kd_2^{n_1}}{[FA]^{n_1}+Kd_1^{n_1}}\right) | ||
+ | $$ | ||
+ | $$ | ||
+ | \frac{d[RFP]}{dt}=\alpha_{RFP} [RFP] \left(\frac{Kd_5^{n_5}}{Kd_5^{n_5}+X_{FadR}^{n_5}}\right) - \gamma_{RFP} [RFP] | ||
+ | $$ | ||
+ | $$ | ||
+ | \frac{d[luxR]}{dt}=\alpha_{luxR} [luxR] \left(\frac{Kd_2^{n_2}}{Kd_2^{n_2}+X_{FadR}^{n_2}}\right) - \gamma_{luxR} [luxR] | ||
+ | $$ | ||
+ | $$ | ||
+ | luxR_{AHL} = [luxR] \left(\frac{[AHL]^{n_3}}{[AHL]^{n_3}+Kd_3^{n_3}}\right) | ||
+ | $$ | ||
+ | $$ | ||
+ | \frac{d[RFP]}{dt} = \alpha_{RFP} \left(\frac{[luxR_{AHL}]^{n_4}}{Kd_4^{n_4}+[luxR_{AHL}]^{n_4}}\right) - \gamma_{RFP} [RFP] | ||
+ | $$ | ||
<p>Both sets of experimental data were convoluted in order to obtain a general deterministic model for pFadBA promoter. | <p>Both sets of experimental data were convoluted in order to obtain a general deterministic model for pFadBA promoter. | ||
</p> | </p> | ||
− | <p> | + | <p style="padding-left:5em">pfadBA_34_RFP experimental fitting to the model.</p> |
− | <p> | + | <img style="width: 30%; display: block; margin: auto" src="https://static.igem.org/mediawiki/2018/f/fd/T--UPF_CRG_Barcelona--transfer_function.png"> |
+ | <p style="padding-left:5em">pFadBA_34_luxR_plux_RFP experimental fitting of the model.</p> | ||
+ | <img style="width: 30%; display: block; margin: auto" src="https://static.igem.org/mediawiki/2018/c/cf/T--UPF_CRG_Barcelona--04mM.png"> | ||
+ | <img style="width: 30%; display: block; margin: auto" src="https://static.igem.org/mediawiki/2018/c/c4/T--UPF_CRG_Barcelona--1mM.png"> | ||
<p class="subapart1">Conclusions:</p> | <p class="subapart1">Conclusions:</p> | ||
Line 118: | Line 137: | ||
behaviour. | behaviour. | ||
</p> | </p> | ||
+ | |||
+ | <p class="subapart1">Parameters:</p> | ||
+ | <table class="fluid-table"> | ||
+ | <tr> | ||
+ | <thead> | ||
+ | <th><span style="font-weight:bold">Parameter</span></th> | ||
+ | <th><span style="font-weight:bold">Value</span></th> | ||
+ | <th>Reference</th> | ||
+ | </thead> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$Kd_1$<br></td> | ||
+ | <td data-label="Value" >1200 $\mu$M</td> | ||
+ | <td data-label="Reference">Data Fitted</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$Kd_2$</td> | ||
+ | <td data-label="Value" >0.13 $\mu$M</td> | ||
+ | <td data-label="Reference">iGEM Taida 2012</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$\alpha_{luxR}$</td> | ||
+ | <td data-label="Value" >0.000325 $\mu$M/min</td> | ||
+ | <td data-label="Reference">Data Fitted</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$\gamma_{luxR}$</td> | ||
+ | <td data-label="Value" >0.0231 $min^{-1}$</td> | ||
+ | <td data-label="Reference">iGEM Zurich ETH 2014</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$n_1$</td> | ||
+ | <td data-label="Value" >3</td> | ||
+ | <td data-label="Reference">iGEM Taida 2012</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$n_2$</td> | ||
+ | <td data-label="Value" >2</td> | ||
+ | <td data-label="Reference">iGEM Taida 2012</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$\alpha_{RFP}$</td> | ||
+ | <td data-label="Value" >3.5</td> | ||
+ | <td data-label="Reference">Data Fitted</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$\gamma_{RFP}$</td> | ||
+ | <td data-label="Value" >0.0041 $min{-1}$</td> | ||
+ | <td data-label="Reference">Data Fitted</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$Kd_3$<br></td> | ||
+ | <td data-label="Value" >0.0001892 nM</td> | ||
+ | <td data-label="Reference">Literature</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$Kd_4$</td> | ||
+ | <td data-label="Value" >0.015 nM</td> | ||
+ | <td data-label="Reference">Literature</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$n_3$</td> | ||
+ | <td data-label="Value" >2</td> | ||
+ | <td data-label="Reference">Literature</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">$n_4$</td> | ||
+ | <td data-label="Value" >1</td> | ||
+ | <td data-label="Reference">Literature</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td data-label="Parameter">[FadR]</td> | ||
+ | <td data-label="Value" >0.16 $\mu$M</td> | ||
+ | <td data-label="Reference">Assumption</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | |||
+ | |||
<div class="spacer"></div> | <div class="spacer"></div> |
Latest revision as of 03:59, 18 October 2018
Single cell model
Motivation:
The last decades have been marked by a high technological development in the field of computation. This has allowed synthetic biology to use modeling tools to predict behaviors and thus make more efficient the construction of genetic circuits.
We have used these tools to model two of our constructs: pFadBA_34_RFP and pFadBA_34_luxR_plux_32_RFP.
Regarding the model for pFadBA_34_RFP:
To analyze the response of pFadBA promoter to palmitic acid (PA) presence, we have used the RFP reporter protein to obtain fluorescence measurements. The issue of those experiments was that the PA has effects in the fluorescence over a concentration of 1mM, which generates a lot of noise. Because of this, we could not obtain experimentally the complete transfer function of this construct.
Nevertheless, we have made parameter optimization with weighted non-linear regression functions in order to obtain the curve that fitted the best our less noisy data. With this, we were able to estimate the whole transfer function of our construct.
Regarding the model for pFadba_34_luxR_plux_32_RFP:
In this case, even if we had the same problem that with pFadBA, since luxR/plux amplifies considerably the signal, we could obtain the whole curve with acceptable noise levels. Here, we also adapted our parameters to obtain a model that reflects the best the behaviour of the construct.
pFadBA model was created to complement and refine the experimental data, in contrast to pFadba_34_luxR_plux_32_RFP whose model was created for a different aim.The pFadba_34_luxR_plux_32_RFP system response is dependent on an eternal inducer: lactone. Having a model that predicts how much lactone is necessary for the modulation we want will facilitate its use in the future.
The ODE for the pFadBA_34_RFP and pFadba_34_luxR_plux_32_RFP systems are given below:
$$ X_{FadR}= [FadR]\left(\frac{Kd_2^{n_1}}{[FA]^{n_1}+Kd_1^{n_1}}\right) $$ $$ \frac{d[RFP]}{dt}=\alpha_{RFP} [RFP] \left(\frac{Kd_5^{n_5}}{Kd_5^{n_5}+X_{FadR}^{n_5}}\right) - \gamma_{RFP} [RFP] $$ $$ \frac{d[luxR]}{dt}=\alpha_{luxR} [luxR] \left(\frac{Kd_2^{n_2}}{Kd_2^{n_2}+X_{FadR}^{n_2}}\right) - \gamma_{luxR} [luxR] $$ $$ luxR_{AHL} = [luxR] \left(\frac{[AHL]^{n_3}}{[AHL]^{n_3}+Kd_3^{n_3}}\right) $$ $$ \frac{d[RFP]}{dt} = \alpha_{RFP} \left(\frac{[luxR_{AHL}]^{n_4}}{Kd_4^{n_4}+[luxR_{AHL}]^{n_4}}\right) - \gamma_{RFP} [RFP] $$Both sets of experimental data were convoluted in order to obtain a general deterministic model for pFadBA promoter.
pfadBA_34_RFP experimental fitting to the model.
pFadBA_34_luxR_plux_RFP experimental fitting of the model.
Conclusions:
We can say that with these models we have obtained a consistent way to predict and analyze the pFadBA promoter behaviour.
Parameters:
Parameter | Value | Reference |
---|---|---|
$Kd_1$ |
1200 $\mu$M | Data Fitted |
$Kd_2$ | 0.13 $\mu$M | iGEM Taida 2012 |
$\alpha_{luxR}$ | 0.000325 $\mu$M/min | Data Fitted |
$\gamma_{luxR}$ | 0.0231 $min^{-1}$ | iGEM Zurich ETH 2014 |
$n_1$ | 3 | iGEM Taida 2012 |
$n_2$ | 2 | iGEM Taida 2012 |
$\alpha_{RFP}$ | 3.5 | Data Fitted |
$\gamma_{RFP}$ | 0.0041 $min{-1}$ | Data Fitted |
$Kd_3$ |
0.0001892 nM | Literature |
$Kd_4$ | 0.015 nM | Literature |
$n_3$ | 2 | Literature |
$n_4$ | 1 | Literature |
[FadR] | 0.16 $\mu$M | Assumption |