Team:TU Delft-Leiden/Modeling

From 2014.igem.org

(Difference between revisions)
 
(97 intermediate revisions not shown)
Line 11: Line 11:
<h2> Modeling Overview</h2>
<h2> Modeling Overview</h2>
 +
<p>We developed models for each of the three different modules of our project: the <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Curli">conductive curli module</a>, the <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/EET">extracellular electron transport (EET) module</a> and the <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine">landmine detection module</a>. <br>
 +
For the conductive curli module, we wanted to know if a conductive path between two electrodes of a chip filled with curli growing <i> E. coli </i> arise at a certain point in time. We also wanted to make quantitative predictions about the resistance between the two electrodes of our system in time. <br>
 +
For the EET module, our goal was to investigate the carbon metabolism providing the electrons for the EET module. Also, we want the EET pathway used by the cells in order to have a measurable electrical signal for our biosensor, see the <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Gadget">gadget section</a> of our wiki. Furthermore, in our modeling of the assembly of the EET complex, we wanted to predict how many EET complexes are formed under different initial conditions. We focused, in addition to the assembly mechanism, also on the apparent reduced cell viability.<br>
 +
For the landmine module, we tried to find a model which would be able to reproduce the response curves of both the landmine promoters, as found in [1]. <br>
 +
For the EET and landmine modules, we used deterministic modeling. For the curli module, we used a stochastic modeling approach, and considered the system at the gene, cell and colony level. At the colony levvel, we employed <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Techniques#PercolationTheory">percolation theory</a> in order to predict if a conductive path between the two electrodes arise at a certain point in time and to predict at which time this happens. Our application of percolation theory to describe the formation of a conductive biological network represents a novel approach that has not been used in the literature before.
 +
</p>
-
<p> We modeled all three different modules our project consists of, namely the landmine module, the Extracellular Electron Transport (EET) module and the conductive curli module. In order to achieve this, we had to use all kinds of different modeling methods.</p>
+
<br>
 +
 
 +
<p>
 +
We used Matlab for most of the calculations; the scripts we made can be found in the <a href="/Team:TU_Delft-Leiden/Modeling/CodeRepository">Code Repository</a>. We had great interactions with the Life Science and Microfluidics departments, which for the conductive curli module can be read <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/curli/integration">here</a>, for the EET module can be read <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/EET/integration">here</a> and for the landmine detection module can be read <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/landmine/integration">here</a>.  
 +
</p>
Line 21: Line 31:
<ul>
<ul>
      
      
-
     <li>
+
      
-
         
+
       
-
          <a href="/Team:TU_Delft-Leiden/Modeling/Landmine">
+
        <a href="/Team:TU_Delft-Leiden/Modeling/Curli">
-
          <p>Landmine Module</p>
+
        <p>Curli Module</p>
-
          </a>
+
        </a>
-
         
+
              <ul>
-
          <ul>
+
                    <li>
-
              <li>
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Curli/Gene">
-
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine#simplemodel">
+
                     <p>Gene Level Modeling</p>
-
                     <p>Simple Binding Model </p>
+
                     </a>
                     </a>
-
              </li>
+
                       
-
              <li>
+
                    </li>
-
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine#coopbinding">
+
 
-
                     <p>Cooperative Binding Model</p>
+
                    <li>
 +
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Curli/Cell">
 +
                     <p>Cell Level Modeling</p>
                     </a>
                     </a>
-
              </li>
+
                       
 +
                    </li>
-
          </ul>
+
                    <li>
-
    </li>
+
                    <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Curli/Colony">
-
+
                    <p>Colony Level Modeling</p>
-
    <li>
+
                    </a>
 +
                     
 +
                              <li>
 +
                              <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Curli/Reflection">
 +
                              <p>Critical reflection on our model</p>
 +
                              </a>
 +
                              </li>
 +
 
 +
                    </li>
 +
              </ul>
 +
 
 +
 
-
           <a href="/Team:TU_Delft-Leiden/Modeling/EET">
+
           <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/EET">
           <p>EET Module</p>
           <p>EET Module</p>
           </a>
           </a>
Line 53: Line 76:
               <li>
               <li>
                    
                    
-
                     <a href="/Team:TU_Delft-Leiden/Modeling/EET/FBA">
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/EET/FBA">
                     <p> Flux Balance Analysis of the EET Module</p>
                     <p> Flux Balance Analysis of the EET Module</p>
                     </a>
                     </a>
Line 62: Line 85:
               <li>
               <li>
-
                     <a href="/Team:TU_Delft-Leiden/Modeling/EET/Deterministic">  
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/EET/Deterministic">  
                     <p> Deterministic Model of EET Complex Assembly</p>
                     <p> Deterministic Model of EET Complex Assembly</p>
                     </a>
                     </a>
Line 73: Line 96:
         </ul>
         </ul>
-
  </li>
+
 
-
 
+
          <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine">
-
  <li>
+
          <p>Landmine Module</p>
-
       
+
          </a>
-
        <a href="/Team:TU_Delft-Leiden/Modeling/Curli">
+
         
-
        <p>Curli Module</p>
+
          <ul>
-
        </a>
+
-
              <ul>
+
-
                    <li>
+
-
                     <a href="/Team:TU_Delft-Leiden/Modeling/Curli#Gene Level">
+
              <li>
-
                     <p>Gene Level Modeling</p>
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine#simplemodel">
 +
                     <p>Simple Binding Model </p>
                     </a>
                     </a>
-
                       
+
              </li>
-
                    </li>
+
              <li>
-
 
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine#coopbinding">
-
                    <li>
+
                     <p>Cooperative Binding Model</p>
-
                     <a href="/Team:TU_Delft-Leiden/Modeling/Curli#Cell Level">
+
-
                     <p>Cell Level Modeling</p>
+
                     </a>
                     </a>
-
                       
+
              </li>
-
                    </li>
+
    <li>
-
 
+
                     <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Landmine#experimentaldata">
-
                    <li>
+
                     <p>Fitting to Experimental Data</p>
-
                     <a href="/Team:TU_Delft-Leiden/Modeling/Curli#Colony Level">
+
-
                     <p>Colony Level Modeling</p>
+
                     </a>
                     </a>
-
                     
+
               </li>
-
                    </li>
+
-
               </ul>
+
-
 
+
-
  </li>
+
-
  <li>
+
          </ul>
 +
 
 +
         
 +
          <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/interactions">
 +
          <p>Interaction with Life Science and Microfluidics</p>
 +
          </a>
 +
         
 +
          <ul>
 +
 
 +
              <li>  
 +
                    <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/curli/integration">
 +
                    <p>Curli Module </p>
 +
                    </a>
 +
              </li>
 +
              <li>
 +
                    <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/EET/integration">
 +
                    <p>EET Module </p>
 +
                    </a>
 +
              </li>
 +
                <li>
 +
                    <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Project/Life_science/landmine/integration">
 +
                    <p>Landmine Module </p>
 +
                    </a>
 +
              </li>
 +
             
 +
          </ul>
 +
           <a href="/Team:TU_Delft-Leiden/Modeling/Techniques">
           <a href="/Team:TU_Delft-Leiden/Modeling/Techniques">
Line 136: Line 175:
           </ul>     
           </ul>     
 +
 +
   
 +
         
 +
          <a href="/Team:TU_Delft-Leiden/Modeling/CodeRepository">
 +
          <p>Code Repository</p>
 +
          </a>
 +
         
-
    </li>
 
</ul>   
</ul>   
Line 143: Line 188:
</div>
</div>
-
<h3>Landmine module </h3>
+
<h3> References </h3>
-
<p>
 
-
An important part of our iGEM project is a promoter sensitive to DNT/TNT. We will use two promoters that are sensitive to DNT/TNT, namely <i>ybiJ</i> and <i>ybiFB2A1</i>, in our project. Of these promoters, not much is known other than the fact that they have a DNT/TNT-dependent response curve . Our goal was to find a model which would be able to reproduce the response curves of both promoters. To achieve this, we constructed two different models, both using <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Techniques#DeterministicTheory">Deterministic Modeling Methods</a>. One model is based on a simple binding model of DNT to the promoter, the other is based on cooperative binding of DNT to the promoter.
 
-
When based on the simple binding model, fits of promoter activation with respect to DNT concentration to the experimental data of [1] did not yield good results. However, when the fits were based on the cooperative binding model, we were able to match the experimental data in [1] really well, see figure 1.
 
-
</p>
 
-
 
-
<figure>
 
-
<img src="https://static.igem.org/mediawiki/2014/5/59/TUDelft_2014_model_fit_tnt_promoters_coop.png" width="100%" width="100%" height="100%">
 
-
<figcaption>
 
-
Figure 1: Fits of the promoter activation model described by cooperative promoter activation to the data of [1]. The left panel shows the fit for the jbiJ promoter, the right panel the fit for the yqjFB2A1 promoter. For comparison, also the fits described by the simple binding model are displayed.
 
-
</figcaption>
 
-
</figure>
 
-
 
-
<br>
 
-
 
-
<h3>EET Module </h3>
 
-
 
-
<p>
 
-
In the wet lab, we integrated the Extracellular Electron Transport (EET) module of <i>S. oneidensis</i> into <i>E. coli</i> <font color="red">[reference]</font>. For the modeling of the EET module, we wanted at first to gain insight in the consequences of the integration of the EET module into <i>E. coli</i>. To achieve this, we simulated the cell metabolism of <i>E. coli</i> including the EET module using the <a href="https://2014.igem.org/Team:TU_Delft-Leiden/Modeling/Techniques#FBATheory">Flux Balance Analysis (FBA) method</a>. <br>
 
-
Our goal was to investigate the carbon metabolism providing the electrons for the EET module <font color="red">[reference]</font>. Also, we want the EET pathway used by the cells in order to have a measurable electrical signal for our biosensor <font color="red">[reference]</font>. <br>
 
-
From the FBA method, we conclude that in aerobic conditions the cell does not use the EET pathway, but oxygen gets reduced instead, as it is a stronger oxidizing agent. However, in anaerobic growth the cell does use the EET pathway to export electrons out of the cell. When the cell is grown on glucose, the growth rate will be higher than when the cell is grown on D-lactate. <br>
 
-
We also conclude that in an experimental setting the EET pathway has a higher chance of being used when the cells are grown on D-lactate as the EET pathway is necessary in order for the cells to grow, while when grown on glucose and the EET pathway is turned off (represented by 0 \(mmol \ (gDW)^{-1} \ hr^{-1}\) (per gram dry weight per hour) maximum EET flux) growth is still possible, see figure 2. From Flux Variability Analysis (FVA) we conclude that for maximum growth for each specific combination of carbon source uptake flux and maximum EET flux, only one possible EET flux is possible for both growth on glucose and growth on D-lactate, namely the EET flux shown in figure 2.
 
-
</p>
 
-
 
-
<figure>
 
-
<img src="https://static.igem.org/mediawiki/2014/9/9b/TUDelft_2014_Phase_plane%28anaerobe%29_opt_growth_v2.png" width="100%" height="100%">
 
-
<figcaption>
 
-
Figure 2: Above: Phenotypic phase planes for growth rate, related to the maximum EET flux and carbon uptake flux, under anaerobic conditions, optimized for growth. The left panel display growth on glucose, the right panel growth on D-lactate. Green means a low growth rate, yellow means a high growth rate. Regions indicated by I correspond to no growth, regions indicated by II correspond to carbon source-limited growth, regions indicated by III correspond to carbon source-limited and maximum EET flux-limited growth. Below: Phenotypic phase planes for EET flux, related to the maximum EET flux and carbon uptake flux, under anaerobic conditions, optimized for growth. The left panel display growth on glucose, the right panel growth on lactate. Green means a low EET flux, yellow means a high EET flux. Regions indicated by I correspond to no EET flux (and no growth), regions indicated by II correspond to carbon source-limited EET flux, and regions indicated by III correspond to maximum EET flux-limited EET flux.
 
-
</figcaption>
 
-
</figure>
 
-
 
-
<br>
 
-
 
-
<p>
 
-
As performing the FBA method while maximizing the EET flux yielded no growth, we wondered if there are pathways possible that would yield growth. So, we performed FVA, the results can be found in figure 3. From this figure, we conclude that when maximizing the EET flux, there are pathways possible that yield growth. Note that the EET flux is 2.8 times higher and 2 times higher in comparison to maximizing for growth rate, for glucose and D-lactate as a carbon source, respectively.
 
-
</p>
 
-
 
-
<figure>
 
-
<img src="https://static.igem.org/mediawiki/2014/0/04/TUDelft_2014_Growth_diff_phase_plane%28anaerobe%29_opt_EET_v2.png" width="100%" height="100%">
 
-
<figcaption>
 
-
Figure 3: Above: Phenotypic phase plane for EET flux, related to the maximum EET flux and carbon uptake flux, under anaerobic conditions, optimized for EET flux. The left panel display growth on glucose, the right panel growth on D-lactate. Green means a low EET flux, yellow means a high EET flux. Regions indicated by I correspond to no EET flux, regions indicated by II correspond to carbon source-limited EET flux, and regions indicated by III correspond to maximum EET flux-limited EET flux. Below: Phenotypic phase planes for growth rate, related to the maximum EET flux and carbon uptake flux, under anaerobic conditions, optimized for EET flux. The left panel display growth on glucose, the right panel growth on D-lactate. Both panels give the difference between the optimized maximum value and the optimized minimum value of growth. Green means a low growth rate, yellow means a high growth rate. Regions indicated by I correspond to no growth, regions indicated by II correspond to carbon source-limited EET flux, and regions indicated by III correspond to carbon source-limited and maximum EET flux-limited growth.
 
-
</figcaption>
 
-
</figure>
 
-
 
-
<br>
 
-
 
-
<p>
 
-
From figures 2 and 3 we conclude that there are different regions in which the cell can operate. In an experimental setting, it can be investigated in which region the cell actually operates and if it maximizes its growth rate or its EET flux. To be able to do this, the experimental observed pathway has to be compared to the possible pathways when maximizing the EET flux and to the pathway when maximizing the growth rate. From these regions, it can be deduced if the experimentally observed EET flux and growth rate are carbon source limited or limited by the maximum possible EET flux. <br>
 
-
Finally, we investigated an extended model of <i>E. coli</i> metabolism. This model contains, in contrast to the previously used core model, L-lactate as a metabolite. Using the extended model, we found that for glucose and D-lactate as carbon sources, the maximized growth rate agreed quite well to the previous analysis's using the core model. Using L-lactate as a carbon source, we conclude that a steady state solution in which <i>E. coli</i> can grow on L-lactate and use the EET pathway is not possible. A possible way to obtain information about the EET flux when the cells are not in steady state as observed by Goldbeck <i>et al.</i> [2], would be by the use of dynamic flux balance analysis (dFBA), which can also model the dynamics of a system before it reaches steady state [3].
 
-
</p>
 
-
 
-
<h3>References </h3>
 
<p>
<p>
[1] S. Yagur-Kroll, S. Belkin <i>et al.</i>, “<i>Escherichia Coli</i> bioreporters for the detection of 2,4-dinitrotoluene and 2,4,6-trinitrotoluene”, Appl. Microbiol. Biotechnol. 98, 885-895, 2014.  
[1] S. Yagur-Kroll, S. Belkin <i>et al.</i>, “<i>Escherichia Coli</i> bioreporters for the detection of 2,4-dinitrotoluene and 2,4,6-trinitrotoluene”, Appl. Microbiol. Biotechnol. 98, 885-895, 2014.  
</p>
</p>
-
<p>
 
-
[2] C.P. Goldbeck, H.M. Jensen <i>et al.</i>, “Tuning Promoter Strengths for Improved Synthesis and Function of Electron Conduits in <i>Escherichia coli</i>”, ACS Synth. Biol. 2, 150-159, 2013.
 
-
</p>
 
-
 
-
<p>
 
-
[3] R. Mahadevan, J.S. Edwards & F.J. Doyle, “Dynamic Flux Balanace Analysis of Diauxic Growth in <i>Escherichia coli</i>”, Biophys. J. 83, 1331-1340, 2002.
 
-
</p>
 
</div>
</div>

Latest revision as of 23:44, 17 October 2014


Modeling Overview

We developed models for each of the three different modules of our project: the conductive curli module, the extracellular electron transport (EET) module and the landmine detection module.
For the conductive curli module, we wanted to know if a conductive path between two electrodes of a chip filled with curli growing E. coli arise at a certain point in time. We also wanted to make quantitative predictions about the resistance between the two electrodes of our system in time.
For the EET module, our goal was to investigate the carbon metabolism providing the electrons for the EET module. Also, we want the EET pathway used by the cells in order to have a measurable electrical signal for our biosensor, see the gadget section of our wiki. Furthermore, in our modeling of the assembly of the EET complex, we wanted to predict how many EET complexes are formed under different initial conditions. We focused, in addition to the assembly mechanism, also on the apparent reduced cell viability.
For the landmine module, we tried to find a model which would be able to reproduce the response curves of both the landmine promoters, as found in [1].
For the EET and landmine modules, we used deterministic modeling. For the curli module, we used a stochastic modeling approach, and considered the system at the gene, cell and colony level. At the colony levvel, we employed percolation theory in order to predict if a conductive path between the two electrodes arise at a certain point in time and to predict at which time this happens. Our application of percolation theory to describe the formation of a conductive biological network represents a novel approach that has not been used in the literature before.


We used Matlab for most of the calculations; the scripts we made can be found in the Code Repository. We had great interactions with the Life Science and Microfluidics departments, which for the conductive curli module can be read here, for the EET module can be read here and for the landmine detection module can be read here.

Contents

References

[1] S. Yagur-Kroll, S. Belkin et al., “Escherichia Coli bioreporters for the detection of 2,4-dinitrotoluene and 2,4,6-trinitrotoluene”, Appl. Microbiol. Biotechnol. 98, 885-895, 2014.

Top
facebook twitter

Retrieved from "http://2014.igem.org/Team:TU_Delft-Leiden/Modeling"