Difference between revisions of "Template:Groningen/Model"

(Created page with "<html> <head> <link rel="stylesheet" type="text/css" href="https://2018.igem.org/Template:Groningen/CSS&action=raw&ctype=text/css" /> <link rel="stylesheet" type="text/css"...")
 
 
(58 intermediate revisions by 5 users not shown)
Line 9: Line 9:
 
<script src="https://code.jquery.com/jquery-3.2.1.min.js" integrity="sha256-hwg4gsxgFZhOsEEamdOYGBf13FyQuiTwlAQgxVSNgt4=" crossorigin="anonymous"></script>
 
<script src="https://code.jquery.com/jquery-3.2.1.min.js" integrity="sha256-hwg4gsxgFZhOsEEamdOYGBf13FyQuiTwlAQgxVSNgt4=" crossorigin="anonymous"></script>
 
<script type="text/javascript" src="https://2018.igem.org/Template:Groningen/Javascript&action=raw&ctype=text/javascript"></script>
 
<script type="text/javascript" src="https://2018.igem.org/Template:Groningen/Javascript&action=raw&ctype=text/javascript"></script>
 +
 +
<style>
 +
 +
.row_m {
 +
  display: flex;
 +
  flex-wrap: wrap;
 +
  margin: auto;
 +
  width: 98%;
 +
 +
}
 +
 +
/* Create four equal column_ms that sits next to each other */
 +
.column_m {
 +
  flex: 31%;
 +
  max-width: 31%;
 +
 +
}
 +
.col_box {
 +
  padding: 4px;
 +
  background-color: #00973930;
 +
  border-left: 6px solid;
 +
  border-color: #00973970;
 +
  margin: 8px;
 +
}
 +
.column_m img {
 +
  margin-top: 8px;
 +
  vertical-align: middle;
 +
  max-width: 80%;
 +
}
 +
 +
/* Responsive layout - makes a two column_m-layout instead of four column_ms */
 +
@media screen and (max-width: 800px) {
 +
  .column_m {
 +
    flex: 50%;
 +
    max-width: 50%;
 +
  }
 +
}
 +
 +
/* Responsive layout - makes the two column_ms stack on top of each other instead of next to each other */
 +
@media screen and (max-width: 600px) {
 +
  .column_m {
 +
    flex: 100%;
 +
    max-width: 100%;
 +
  }
 +
}
 +
 +
</style>
 
</head>
 
</head>
 
<body>
 
<body>
Line 15: Line 62:
  
 
     <div class="igem_2018_team_column_wrapper">
 
     <div class="igem_2018_team_column_wrapper">
         <div class="column two_thirds_size">
+
         <div class="column">
             <h1> Work in Progress </h1>
+
         
<h4> </h4>
+
             <img src="https://static.igem.org/mediawiki/2018/a/a5/T--Groningen--banner_modeling.png" class="responsive-img">
    <p>   </p>
+
<h1>Overview</h1>
        </div>
+
<p>
<div class="clear"></div>
+
The IGEM team Groningen has invested a lot of effort into developing <u>sophisticated models</u> that simulates all parts of our project. In our quest for producing styrene from the polysaccharide cellulose, the first step is to get our enzymes to the place they need to go; cellulose. As the cellulose binding domain of our mini-cellulosome is responsible for this task, we characterized its cellulose binding properties by creating a <u>cutting edge</u> coarse grained molecular dynamics simulation and running it on our <b>6652 core supercomputer cluster</b> peregrine. The simulation shows the cellulose binding domain as an affinity for cellulose several <b>orders of magnitude higher</b> than the enzymes alone and draw novel insights from this. However by restraining the enzymes together in a scaffold protein, the added rigidity might prove detrimental to enzyme activity. We used an <u>advanced mathematical model</u> to work out the <b>complex system of differential equations</b> that describe this restrained situation, and compared the results to the solubilized enzymes. Luckily, the model shows that restraining the enzymes only impacts their performance negligibly. Finally, we once more harnessed the <b>supercomputing power</b> at our disposal to simulate our synthetic styrene production pathway in the metabolism of <i>S. cerevisiae</i> using a <u>flux based model</u>. We confirmed that yeast is indeed capable of simultaneous growth and high theoretical styrene production. Most strikingly however, we discovered several <b>important metabolic engineering targets</b>, some of which are corroborated by empirical evidence, while others are <b>entirely novel discoveries</b>. Overall all our models have provided us with <u>key insights</u> to aid us in reaching our goal: <b>a sustainable future</b>.
 +
</p>
  
<div class="clear"></div>
+
<div class="clear extra_space"></div>
 +
<div class="line_divider"></div>
 +
<div class="clear_extra_space"></div>
  
 +
<img src="https://static.igem.org/mediawiki/2018/f/f2/T--Groningen--TAPDATASSyesweretiredPARTTWO.png" width="98%">
  
<div class="column full_size">
+
<div class="clear extra_space"></div>
<h1> Modeling</h1>
+
<div class="line_divider"></div>
 +
<div class="clear_extra_space"></div>
 +
<br>
 +
<div class="row_m" align="middle">
 +
  <div class="column_m col_box">
 +
  <a href="https://2018.igem.org/Team:Groningen/Model/Molecular_Dynamics">
 +
  <h5>Proximity by affinity</h5>
 +
  <img src="https://static.igem.org/mediawiki/2018/9/96/T--Groningen--MD_thumb.png">
 +
  </a>
 +
  <br>
 +
  <p style="text-align: justify; padding: 8%">Our cutting edge molecular dynamics simulation uses the in-house developed MARTINI coarse grained force field to effectively model the thermodynamic properties of a reductionist view of our system. </p>
 +
  </div>
 +
    <div class="column_m col_box">
 +
    <a href="https://2018.igem.org/Team:Groningen/Model/Mathematical_Modeling">
 +
    <h5>Cellulose degradation</h5>
 +
   
 +
    <img src="https://static.igem.org/mediawiki/2018/0/03/T--Groningen--math_thumb.png">
 +
    </a>
 +
    <p style="text-align: justify; padding: 8%">We were able to model the kinetic behaviour of a reduced enzyme scaffold. This system was described by a system of differential equations that describe the behaviour of the enzyme complex.</p>
 +
    </div>
 +
<div class="column_m col_box">
 +
<a href="https://2018.igem.org/Team:Groningen/Model/Flux_Based_Analysis">
 +
<h5>Optimizing styrene production</h5>
 +
<img src="https://static.igem.org/mediawiki/2018/4/4e/T--Groningen--fba_thumb.png">
 +
</a>
 +
<p style="text-align: justify; padding: 8%">By using flux based analysis, the complex network of reactions in the metabolism of <i>S. cerevisiae</i> could be modeled. We used this approach to simulate styrene production and to find metabolic engineering targets.</p>
 +
</div>
 +
</div>
  
<p>Mathematical models and computer simulations provide a great way to describe the function and operation of BioBrick Parts and Devices. Synthetic Biology is an engineering discipline, and part of engineering is simulation and modeling to determine the behavior of your design before you build it. Designing and simulating can be iterated many times in a computer before moving to the lab. This award is for teams who build a model of their system and use it to inform system design or simulate expected behavior in conjunction with experiments in the wetlab.</p>
 
  
</div>
+
<!--
<div class="clear"></div>
+
    <p>Our integrated styrene production pathway in Saccharomyces cerevisiae relies on several key properties to make the process the most efficient in ideal conditions. Our enzyme scaffold which supposedly enhances enzyme activity by bringing the enzymes in closer proximity to each other, and closer to cellulose by a cellulose binding domain. Furthermore our yeast strain needs several metabolic optimizations to maximize styrene production levels. </p>
  
<div class="column full_size">
+
    <p>To verify whether an enzyme scaffold indeed enhances enzyme activity, a <a href="https://2018.igem.org/Team:Groningen/Model/Mathematical_Modeling">mathematical  model</a> was created to describe catalyzed degradation of cellulose in the presence of our enzymes. The model, based on work by Levine et al. [1], includes an exoglucanase and an endoglucanase, both in bound together and separately in solution.  
<h3> Gold Medal Criterion #3</h3>
+
  The model provides a justification for the use of scaffolding and insight into how to choose enzymes for a scaffold</p>
<p>
+
Convince the judges that your project's design and/or implementation is based on insight you have gained from modeling. This could be either a new model you develop or the implementation of a model from a previous team. You must thoroughly document your model's contribution to your project on your team's wiki, including assumptions, relevant data, model results, and a clear explanation of your model that anyone can understand.  
+
<br><br>
+
The model should impact your project design in a meaningful way. Modeling may include, but is not limited to, deterministic, exploratory, molecular dynamic, and stochastic models. Teams may also explore the physical modeling of a single component within a system or utilize mathematical modeling for predicting function of a more complex device.
+
</p>
+
  
<p>
+
    <p>As it is likely that substrate availability is higher when the scaffold is strongly bound to cellulose, we explored the binding characteristics of the cellulose binding domain that is attached to our minicellulosome scaffold using <a href="https://2018.igem.org/Team:Groningen/Model/Molecular_Dynamics">coarse grained molecular dynamics</a>. The cellulose binding domain (CBD) of Cip3A from Clostridium thermocellum (PDB 1NBC by Tormo et al. [2]) was modeled in the presence of a single cellulose fiber using coarse grained molecular dynamics simulations [3]. The simulations affirm the hypothesized binding mechanism of CBD to cellulose, and show that CBD from C. thermocellum has a very high affinity for cellulose fibers. From these simulations, the potential of mean force of the binding interaction could be computed, from which the binding free energy could be derived. The results show that our CBD binds even stronger to cellulose than some other cellulose binding domains known from literature. These results justify including the CBD in the enzyme scaffold as it will likely increase the activity of the complex by substrate proximity.</p>
Please see the <a href="https://2018.igem.org/Judging/Medals"> 2018
+
Medals Page</a> for more information.  
+
</p>
+
</div>
+
  
<div class="column two_thirds_size">
+
    <p>Finally <a href="https://2018.igem.org/Team:Groningen/Model/Flux_Based_Analysis">flux balance modeling</a> was used to model the metabolism of S. cerevisiae after introduction of our heterologous pathway enabling styrene production. Yeast GEM 7.6 [4] was used to model the metabolic network of yeast. Furthermore, the OptForce algorithm [5], implemented in the COBRA toolbox [6], was used to find compatible sets of up- or down-regulated genes.
<h3>Best Model Special Prize</h3>
+
    When optimizing for biomass, it was shown that 65% of maximum styrene production is still compatible with normal growth of yeast, justifying our choice for yeast as our host organism.</p>
 +
    <p>To speed up simulations, our models were run using the Groningen University’s state of the art 4000+ core Peregrine cluster.</p>
  
<p>
 
To compete for the <a href="https://2018.igem.org/Judging/Awards">Best Model prize</a>, please describe your work on this page  and also fill out the description on the <a href="https://2018.igem.org/Judging/Judging_Form">judging form</a>. Please note you can compete for both the gold medal criterion #3 and the best model prize with this page.
 
<br><br>
 
You must also delete the message box on the top of this page to be eligible for the Best Model Prize.
 
</p>
 
  
</div>
+
<h4>References</h4>
 +
    <p>[1] Levine, S. E., Fox, J. M., Blanch, H. W., & Clark, D. S. (2010). A mechanistic model of the enzymatic hydrolysis of cellulose. Biotechnology and Bioengineering. <a href="https://doi.org/10.1002/bit.22789" target="_blank">https://doi.org/10.1002/bit.22789</a></p>
  
 +
    <p>[2] Tormo, J., Lamed, R., Chirino, A. J., Morag, E., Bayer, E. A., Shoham, Y., & Steitz, T. A. (1996). Crystal structure of a bacterial family-III cellulose-binding domain: a general mechanism for attachment to cellulose. The EMBO Journal, 15(21), 5739–5751. Retrieved from <a href="http://www.ncbi.nlm.nih.gov/pubmed/8918451" target="_blank">http://www.ncbi.nlm.nih.gov/pubmed/8918451</a></p>
  
<div class="column third_size">
+
    <p>[3] Marrink, S. J., Risselada, H. J., Yefimov, S., Tieleman, D. P., & Vries, A. H. de. (2007). The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. <a href="https://doi.org/10.1021/JP071097F" target="_blank">https://doi.org/10.1021/JP071097F</a></p>
<div class="highlight decoration_A_full">
+
<h3> Inspiration </h3>
+
<p>
+
Here are a few examples from previous teams:
+
</p>
+
<ul>
+
<li><a href="https://2016.igem.org/Team:Manchester/Model">2016 Manchester</a></li>
+
<li><a href="https://2016.igem.org/Team:TU_Delft/Model">2016 TU Delft</li>
+
<li><a href="https://2014.igem.org/Team:ETH_Zurich/modeling/overview">2014 ETH Zurich</a></li>
+
<li><a href="https://2014.igem.org/Team:Waterloo/Math_Book">2014 Waterloo</a></li>
+
</ul>
+
</div>
+
</div>
+
  
        <div class="column third_size">
+
    <p>[4] Sánchez, B. J., Zhang, C., Nilsson, A., Lahtvee, P., Kerkhoven, E. J., & Nielsen, J. (2017). Improving the phenotype predictions of a yeast genome‐scale metabolic model by incorporating enzymatic constraints. Molecular Systems Biology. <a href="https://doi.org/10.15252/msb.20167411" target="_blank">https://doi.org/10.15252/msb.20167411</a></p>
            <div class="clear extra_space"></div>
+
 
            <img src="https://static.igem.org/mediawiki/2018/4/48/T--Groningen--Earth-fossil-fuels.png">
+
    <p>[5] Ranganathan, S., Suthers, P. F., & Maranas, C. D. (2010). OptForce: An optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Computational Biology. <a href="https://doi.org/10.1371/journal.pcbi.1000744" target="_blank">https://doi.org/10.1371/journal.pcbi.1000744</a></p>
 +
 
 +
    <p>[6] Heirendt, L., Arreckx, S., Pfau, T., Mendoza, S. N., Richelle, A., Heinken, A., … Fleming, R. M. T. (2017). Creation and analysis of biochemical constraint-based models: the COBRA Toolbox v3.0. ArXiv. <a href="https://doi.org/10.1038/protex.2011.234" target="_blank">https://doi.org/10.1038/protex.2011.234</a></p>
 +
    -->
 +
        </div>
 +
<div class="clear"></div>
 +
       
 +
        </div>
 
         </div>
 
         </div>
            <div class="clear extra_space"></div>
 
            <div class="line_divider"></div>
 
            <div class="clear_extra_space"></div>
 
   
 
    </div>
 
</div>
 
  
 
<!-- EINDE CONTENT --------------------------------------------------->
 
<!-- EINDE CONTENT --------------------------------------------------->

Latest revision as of 18:30, 7 December 2018

Overview

The IGEM team Groningen has invested a lot of effort into developing sophisticated models that simulates all parts of our project. In our quest for producing styrene from the polysaccharide cellulose, the first step is to get our enzymes to the place they need to go; cellulose. As the cellulose binding domain of our mini-cellulosome is responsible for this task, we characterized its cellulose binding properties by creating a cutting edge coarse grained molecular dynamics simulation and running it on our 6652 core supercomputer cluster peregrine. The simulation shows the cellulose binding domain as an affinity for cellulose several orders of magnitude higher than the enzymes alone and draw novel insights from this. However by restraining the enzymes together in a scaffold protein, the added rigidity might prove detrimental to enzyme activity. We used an advanced mathematical model to work out the complex system of differential equations that describe this restrained situation, and compared the results to the solubilized enzymes. Luckily, the model shows that restraining the enzymes only impacts their performance negligibly. Finally, we once more harnessed the supercomputing power at our disposal to simulate our synthetic styrene production pathway in the metabolism of S. cerevisiae using a flux based model. We confirmed that yeast is indeed capable of simultaneous growth and high theoretical styrene production. Most strikingly however, we discovered several important metabolic engineering targets, some of which are corroborated by empirical evidence, while others are entirely novel discoveries. Overall all our models have provided us with key insights to aid us in reaching our goal: a sustainable future.


Proximity by affinity

Our cutting edge molecular dynamics simulation uses the in-house developed MARTINI coarse grained force field to effectively model the thermodynamic properties of a reductionist view of our system.

Cellulose degradation

We were able to model the kinetic behaviour of a reduced enzyme scaffold. This system was described by a system of differential equations that describe the behaviour of the enzyme complex.

Optimizing styrene production

By using flux based analysis, the complex network of reactions in the metabolism of S. cerevisiae could be modeled. We used this approach to simulate styrene production and to find metabolic engineering targets.