Team:UNITN-Trento/Modeling

From 2013.igem.org

(Difference between revisions)
 
(39 intermediate revisions not shown)
Line 1: Line 1:
-
<!--OLD PAGE-->
+
<!--NEW PAGE-->
 +
{{:Team:UNITN-Trento/Templates/Default|<html><!--start content-->
-
<!-- *** What falls between these lines is the Alert Box!  You can remove it from your pages once you have read and understood the alert *** -->
+
<!--CSS-->
 +
<link rel="stylesheet" href="https://2013.igem.org/wiki/index.php?title=Team:UNITN-Trento/CSS/Project_pages&action=raw&ctype=text/css" type="text/css" />
 +
<link rel="stylesheet" href="https://2013.igem.org/wiki/index.php?title=Team:UNITN-Trento/CSS/Project/Modeling&action=raw&ctype=text/css" type="text/css" />
-
<html>
+
<!-- PAGE -->
-
<div id="box" style="width: 700px; margin-left: 137px; padding: 5px; border: 3px solid #000; background-color: #fe2b33;">
+
<div class="container">
-
<div id="template" style="text-align: center; font-weight: bold; font-size: large; color: #f6f6f6; padding: 5px;">
+
<div class="sheet">
-
This is a template page. READ THESE INSTRUCTIONS.
+
-
</div>
+
-
<div id="instructions" style="text-align: center; font-weight: normal; font-size: small; color: #f6f6f6; padding: 5px;">
+
-
You are provided with this team page template with which to start the iGEM season.  You may choose to personalize it to fit your team but keep the same "look." Or you may choose to take your team wiki to a different level and design your own wiki.  You can find some examples <a href="https://2008.igem.org/Help:Template/Examples">HERE</a>.
+
-
</div>
+
-
<div id="warning" style="text-align: center; font-weight: bold; font-size: small; color: #f6f6f6; padding: 5px;">
+
-
You <strong>MUST</strong>  have all of the pages listed in the menu below with the names specified.  PLEASE keep all of your pages within your teams namespace. 
+
-
</div>
+
-
</div>
+
-
</html>
+
-
<!-- *** End of the alert box *** -->
+
<span class="tn-title">Modeling</span>
 +
<p>
 +
<i>B. fruity</i> exploits bacterial ethylene to ripen fruit. Given ethylene dangerous characteristics (flammable, asphyxiating, explosive) it is even more important to maintain ethylene concentration under a certain threshold or, even better, inside a certain interval.
 +
</p>
 +
<p>
 +
For this purpose, we planned to model our light-regulated system and use the results to determine the required light induction pattern. Given the complexity of the system and the fact that many variables are not present in literature, we hypothesized intra-cellular concentrations of reporter protein and estimated the parameters from these hypothetic data. This might be used as a proof of concept.
 +
</p>
 +
<p>
 +
Of course a better approach would be to obtain experimental data on the intracellular concentration of the reporter protein. We are in the process of improving the model.
 +
</p>
 +
<span class="tn-subtitle">The model</span>
 +
<p>
 +
For simplicity, we did not model the entire circuit (that can be found in the <a href="https://2013.igem.org/Team:UNITN-Trento/Project/Datapage">datapage</a>). Instead, the model is based on the <a href="http://parts.igem.org/Part:BBa_K1065302">BBa_K1065302</a> part, that consists of YF1 and FixJ under R0010 followed by the amilGFP reporter gene under pFixK2. YF1 is the photo-sensor, FixJ is the response regulator (phosphorylated by YF1) and pFixK2 is the target of FixJ-P.
 +
</p>
-
{| style="color:#1b2c8a;background-color:#0c6;" cellpadding="3" cellspacing="1" border="1" bordercolor="#fff" width="62%" align="center"
+
<p style="margin-bottom: 0;">
-
!align="center"|[[Team:UNITN-Trento|Home]]
+
The main differences with the complete circuit are:
-
!align="center"|[[Team:UNITN-Trento/Team|Team]]
+
</p>
-
!align="center"|[https://igem.org/Team.cgi?year=2013&team_name=UNITN-Trento Official Team Profile]
+
<ol>
-
!align="center"|[[Team:UNITN-Trento/Project|Project]]
+
<li>the inverter, which adds complexity, is not included in the model;</li>
-
!align="center"|[[Team:UNITN-Trento/Parts|Parts Submitted to the Registry]]
+
<li>instead of EFE, the model contains the amilGFP reporter gene, which will simplify the retrieving of experimental data for parameter estimation.</li>
-
!align="center"|[[Team:UNITN-Trento/Modeling|Modeling]]
+
</ol>
-
!align="center"|[[Team:UNITN-Trento/Notebook|Notebook]]
+
-
!align="center"|[[Team:UNITN-Trento/Safety|Safety]]
+
-
!align="center"|[[Team:UNITN-Trento/Attributions|Attributions]]
+
-
|}
+
 +
<p style="margin-bottom: 0;">
 +
Also, a couple of assumptions/approximations have been arranged to simplify the model even more:
 +
</p>
 +
<ul>
 +
<li>the LOV domain is responsible for YF1 autodimerization. From literature it appears that YF1 has kinase net activity in the dark and phosphatase net activity in the light. It is not clear whether YF1 is able to phosphorylate FixJ only if dimerized, so the model does not include the dimerization reaction.</li>
 +
<li>in the light, YF1 has phosphatase net activity and so is responsible for FixJ-P de-phosphorylation. The model does not contain all the details of this reaction since they are not meaningful for our purpose.</li>
 +
</ul>
-
If you choose to include a '''Modeling''' page, please write about your modeling adventures here. This is not necessary but it may be a nice list to include.
+
<b>Figure 1</b> shows the modeled pathway. You can download our model (COPASI) <a href="https://static.igem.org/mediawiki/2013/8/84/Tn-2013-model.zip">here</a>.
 +
<div class="tn-doublephoto-wrap">
 +
<img class="photo" src="https://static.igem.org/mediawiki/2013/7/7a/Tn-2013-mod-Dark_modeling.jpg" />
 +
<img class="photo" src="https://static.igem.org/mediawiki/2013/9/95/Tn-2013-mod-LIGHT_MODELING.jpg" />
 +
</div>
 +
<span class="tn-caption"><b>Figure 1:</b> the image shows the model behavior in the light (right) and in the dark (left). Red arrows indicate degradation reactions. In the light YF1 acts both as kinase and phosphatase.</span>
-
<!--NEW PAGE-->
+
<p style="margin-bottom: 0;">
-
{{:Team:UNITN-Trento/Templates/Default|<html><!--start content-->
+
Some kinetic constants have been found in literature and used as center of an interval to estimate the variables of our model, other variables were estimated from hypothetical data within constraints of our choice:
 +
</p>
 +
<ul style="margin-left: 3em; margin-right: 4em;">
 +
<li>Starting DNA concentrations were estimated based on hypothetical data within the interval [1e-6;1e-2]</li>
 +
<li>FixJ phosphorylation and FixJ-P de-phosphorylation kinetics, k4 and k4b, were estimated based on hypothetical data within the interval [0;10] 1/s</li>
 +
<li>Transcription kinetic was taken from <a href="http://kirschner.med.harvard.edu/files/bionumbers/Kinetic%20parameters%20used%20in%20transcription%20factor%20DNA%20binding%20model.pdf">BIONUMBERS</a> (k = 0.33 1/s) and estimated for each reaction inside the interval [1e-4;1e-1] 1/s</li>
 +
<li>RNA degradation kinetic was taken from <a href="http://kirschner.med.harvard.edu/files/bionumbers/Kinetic%20parameters%20used%20in%20transcription%20factor%20DNA%20binding%20model.pdf">BIONUMBERS</a> (k = 0.011 1/s) and estimated for each reaction inside the interval [1e-3;1e-1] 1/s</li>
 +
<li>Protein degradation kinetic was taken from <a href="http://kirschner.med.harvard.edu/files/bionumbers/Kinetic%20parameters%20used%20in%20transcription%20factor%20DNA%20binding%20model.pdf">BIONUMBERS</a> (k = 0.00083 1/s) and estimated for each reaction inside the interval [1e-3;1e-1] 1/s</li>
 +
<li>Protein translation kinetic was taken from <a href="http://kirschner.med.harvard.edu/files/bionumbers/Kinetic%20parameters%20used%20in%20transcription%20factor%20DNA%20binding%20model.pdf">BIONUMBERS</a> (k = 22 1/s) and estimated within the interval [1e-2;1] 1/s</li>
 +
<li>The autophosphorylation kinetics in the light and in the dark, k3 and k3b, were taken from <span class="tn-ref">(M&ouml;glich A, <i>J MOL BIOL</i>, 2009, 385(5):1433-44)</span></li>
 +
</ul>
-
:D MODELING PAGE :D
+
<p>
 +
The result of a time-course simulation of 2e+05 seconds with a 8e+04 dark exposure window is shown in <b>Figure 2</b>.
 +
</p>
 +
 
 +
<img style="width: 600px;" src="https://static.igem.org/mediawiki/2013/6/65/Tn-2013-mod-Plot1-amilGFP-protein.png" />
 +
<span class="tn-caption"><b>Figure 2:</b> the image shows the protein (also in the phosphorylated form, x-P) concentration over time. AmilGFP is produced only during the window of dark exposure.</span>
 +
 
 +
<span class="tn-subtitle">The light induction pattern</span>
 +
<p>We added windows of "dark" and "light" exposure into the model. At first we simply added single windows of "dark" of various length, then we repeated the windows to make the enzyme production fluctuate into the desired interval.</p>
 +
 +
<img src="https://static.igem.org/mediawiki/2013/9/97/Tn-2013-mod-Interval-total.png" />
 +
<span class="tn-caption"><b>Figure 3.</b> Left panel shows the purpose of our model: keeping the concentration of the protein under pFixK2 inside a certain interval. Right panel shows windows of light/dark exposure. Windows of 4,8e+03 s of dark exposure are spaced out by windows of 1e+03 s of light exposure. Thus, certain patterns of light induction can be theoretically used to maintain EFE concentration (and consequently ethylene concentration) under control.</span>
 +
 
 +
<p>This represents a proof of concept, since we are now in the process of implementing the model with experimental data.</p>
 +
</div>
 +
 
 +
<div class="sheet-2">
 +
    <a href="https://2013.igem.org/Team:UNITN-Trento/Project/Bacillus">
 +
        <img class="tn-arr-prev" src="https://static.igem.org/mediawiki/2013/5/58/Tn-2013-arr-Modeling_prev.png" />
 +
    </a>
 +
    <a href="javascript:toTop('#tn-main-wrap-wrap');">
 +
        <img src="https://static.igem.org/mediawiki/2013/6/6e/Tn-2013-arr-AAA_TOP.png" />
 +
    </a>
 +
    <a href="https://2013.igem.org/Team:UNITN-Trento/Parts">
 +
        <img class="tn-arr-next" src="https://static.igem.org/mediawiki/2013/9/95/Tn-2013-arr-Modeling_next.png" />
 +
    </a>
 +
</div>
 +
 
 +
</div>
-
<!--end content--></html>}}
+
<!--end content--></html>|<html>https://static.igem.org/mediawiki/2013/a/a4/Tn-2013-headerbg-Sfondolm_or.jpg</html>|<html>https://static.igem.org/mediawiki/2013/1/10/Tn-2013-headerbg-Sfondolm.jpg</html>}}

Latest revision as of 23:32, 4 October 2013

Modeling

B. fruity exploits bacterial ethylene to ripen fruit. Given ethylene dangerous characteristics (flammable, asphyxiating, explosive) it is even more important to maintain ethylene concentration under a certain threshold or, even better, inside a certain interval.

For this purpose, we planned to model our light-regulated system and use the results to determine the required light induction pattern. Given the complexity of the system and the fact that many variables are not present in literature, we hypothesized intra-cellular concentrations of reporter protein and estimated the parameters from these hypothetic data. This might be used as a proof of concept.

Of course a better approach would be to obtain experimental data on the intracellular concentration of the reporter protein. We are in the process of improving the model.

The model

For simplicity, we did not model the entire circuit (that can be found in the datapage). Instead, the model is based on the BBa_K1065302 part, that consists of YF1 and FixJ under R0010 followed by the amilGFP reporter gene under pFixK2. YF1 is the photo-sensor, FixJ is the response regulator (phosphorylated by YF1) and pFixK2 is the target of FixJ-P.

The main differences with the complete circuit are:

  1. the inverter, which adds complexity, is not included in the model;
  2. instead of EFE, the model contains the amilGFP reporter gene, which will simplify the retrieving of experimental data for parameter estimation.

Also, a couple of assumptions/approximations have been arranged to simplify the model even more:

  • the LOV domain is responsible for YF1 autodimerization. From literature it appears that YF1 has kinase net activity in the dark and phosphatase net activity in the light. It is not clear whether YF1 is able to phosphorylate FixJ only if dimerized, so the model does not include the dimerization reaction.
  • in the light, YF1 has phosphatase net activity and so is responsible for FixJ-P de-phosphorylation. The model does not contain all the details of this reaction since they are not meaningful for our purpose.
Figure 1 shows the modeled pathway. You can download our model (COPASI) here.
Figure 1: the image shows the model behavior in the light (right) and in the dark (left). Red arrows indicate degradation reactions. In the light YF1 acts both as kinase and phosphatase.

Some kinetic constants have been found in literature and used as center of an interval to estimate the variables of our model, other variables were estimated from hypothetical data within constraints of our choice:

  • Starting DNA concentrations were estimated based on hypothetical data within the interval [1e-6;1e-2]
  • FixJ phosphorylation and FixJ-P de-phosphorylation kinetics, k4 and k4b, were estimated based on hypothetical data within the interval [0;10] 1/s
  • Transcription kinetic was taken from BIONUMBERS (k = 0.33 1/s) and estimated for each reaction inside the interval [1e-4;1e-1] 1/s
  • RNA degradation kinetic was taken from BIONUMBERS (k = 0.011 1/s) and estimated for each reaction inside the interval [1e-3;1e-1] 1/s
  • Protein degradation kinetic was taken from BIONUMBERS (k = 0.00083 1/s) and estimated for each reaction inside the interval [1e-3;1e-1] 1/s
  • Protein translation kinetic was taken from BIONUMBERS (k = 22 1/s) and estimated within the interval [1e-2;1] 1/s
  • The autophosphorylation kinetics in the light and in the dark, k3 and k3b, were taken from (Möglich A, J MOL BIOL, 2009, 385(5):1433-44)

The result of a time-course simulation of 2e+05 seconds with a 8e+04 dark exposure window is shown in Figure 2.

Figure 2: the image shows the protein (also in the phosphorylated form, x-P) concentration over time. AmilGFP is produced only during the window of dark exposure. The light induction pattern

We added windows of "dark" and "light" exposure into the model. At first we simply added single windows of "dark" of various length, then we repeated the windows to make the enzyme production fluctuate into the desired interval.

Figure 3. Left panel shows the purpose of our model: keeping the concentration of the protein under pFixK2 inside a certain interval. Right panel shows windows of light/dark exposure. Windows of 4,8e+03 s of dark exposure are spaced out by windows of 1e+03 s of light exposure. Thus, certain patterns of light induction can be theoretically used to maintain EFE concentration (and consequently ethylene concentration) under control.

This represents a proof of concept, since we are now in the process of implementing the model with experimental data.

[http://2013.igem.org/wiki/index.php?title=Team:UNITN-Trento/Modeling&action=edit Edit this page] | Main Page