Team:KU Leuven/Journal/Modelling

From 2013.igem.org

(Difference between revisions)
 
(34 intermediate revisions not shown)
Line 2: Line 2:
{{:Team:KU Leuven/Template:Style}}
{{:Team:KU Leuven/Template:Style}}
{{:Team:KU Leuven/Template:Menu}}
{{:Team:KU Leuven/Template:Menu}}
 +
<html>
 +
<body>
 +
<div id="container">
 +
<div class="container">
 +
<!--LOGOS-->
-
<html>
+
<br/>
-
<body>
+
<div class="row-fluid">
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">
 +
  <div class="visible-desktop span3 bg-yellow">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/Modelling">
 +
    <i class="livicon activeicon" data-name="laptop" data-onparent="true" data-color="white"></i>
 +
  </div>
 +
  <div class="span7 icon-text">
 +
    <h3>Modelling</h3> </a>
 +
    <p>You are here!</p>
 +
  </div>
 +
  </div>
 +
</div>
 +
 
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">
 +
  <div class="farnesene visible-desktop span3 bg-grey1">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/Ecology">
 +
    <i class="livicon activeicon" data-name="tree" data-onparent="true" data-color="white"></i>
 +
    </div>
 +
  <div class="span7 icon-text">
 +
    <h3>Ecology</h3>
 +
    </a>
 +
    <p>Plant and insect experiments</p>
 +
  </div>
 +
  </div>
 +
</div>
 +
 
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">
 +
  <div class="MeS visible-desktop span3 bg-grey1">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/EBF/wetlab">
 +
    <i class="livicon activeicon" data-name="bell" data-onparent="true" data-color="white"></i>
 +
    </div>
 +
  <div class="span7 icon-text">
 +
    <h3>E-β-farnesene</h3> </a>
 +
    <p>Work about our EBF-biobrick</p>
 +
  </div>
 +
  </div>
 +
</div>
 +
   
 +
</div>
 +
 
 +
<div class="row-fluid">
 +
 
 +
 
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">   
 +
    <div class="oscillator visible-desktop span3 bg-grey1">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/FFL/wetlab">
 +
    <i class="livicon activeicon" data-name="gears" data-onparent="true" data-color="white"></i>
 +
    </div>
 +
  <div class="span7 icon-text">
 +
    <h3>Feed forward loop</h3> </a>
 +
    <p>Implementation of the oscillator <i>in vivo</i></p>
 +
  </div>
 +
  </div>
 +
</div>
 +
 
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">
 +
  <div class="ecology visible-desktop span3 bg-grey1">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/MeS/wetlab">
 +
    <i class="livicon activeicon" data-name="lab" data-onparent="true" data-color="white"></i>
 +
  </div>
 +
  <div class="span7 icon-text">
 +
    <h3>Methyl salicylate</h3>
 +
    </a>
 +
    <p>Work about our MeS-biobrick</p>
 +
  </div>
 +
  </div>
 +
</div>
 +
 
 +
<div class="span4 icon white">
 +
  <div class="row-fluid">
 +
  <div class="parts visible-desktop span3 bg-grey1">
 +
    <a href="https://2013.igem.org/Team:KU_Leuven/Journal/MeS/qPCR">
 +
  <i class="livicon activeicon" data-name="spinner-six" data-onparent="true" data-color="white"></i>
 +
  </div>
 +
  <div class="span7 icon-text">
 +
    <h3>qPCR</h3> </a>
 +
    <p>Finding wetlab data for the models</p>
 +
  </div>
 +
  </div>
 +
</div>
 +
 +
</div>
 +
 
 +
<!--TITEL-->
 +
 
 +
<div id="header" class="row-fluid">
 +
<div class="span12">
 +
  <h3 class="bg-yellow">Modelling Journal</h3>
 +
</div>
 +
</div>
 +
 
 +
<!-- TEKST -->
 +
<div class="row-fluid">
 +
<div class="span12 white">
 +
  <p align="justify">This is the journal of the modelling team, they <a href="https://2013.igem.org/Team:KU_Leuven/Project/Oscillator/Description"> designed</a> and <a href="https://2013.igem.org/Team:KU_Leuven/Project/Oscillator/Modelling">analysed</a> an <a href ="https://2013.igem.org/Team:KU_Leuven/Project/Oscillator"> oscillator</a>, delved into Cobra, <a href="https://2013.igem.org/Team:KU_Leuven/Project/Glucosemodel/MeS/Modelling">modelled the production system of methyl salicylate</a> and even found the time to start a <a href="https://2013.igem.org/Team:KU_Leuven/Project/Ecological/Modelling">model</a> about how our products will be dispersed on a field over time.</p>
 +
</div>
 +
</div>
<div class="row-fluid">
<div class="row-fluid">
<div class="span12 timeline">
<div class="span12 timeline">
<ul class="cbp_tmtimeline">
<ul class="cbp_tmtimeline">
-
 
  <li>
  <li>
Line 19: Line 124:
     <h3>Week 1: The modelling team is on fire!</h3>
     <h3>Week 1: The modelling team is on fire!</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">The modelling team is ready to rock Matlab after watching some tutorials and webinars about modelling and <a href="http://www.mathworks.nl/products/simbiology/webinars.html">Symbiology</a><br/>
+
     <p align="justify">The modelling team is ready to rock MATLAB after watching some tutorials and webinars about modelling and <a href="http://www.mathworks.nl/products/simbiology/webinars.html">Symbiology</a><br/>
-
We are excited to get our modelling  team going and therefore we had an appointment with professor <a href="http://www.kuleuven.be/optec/people/74-Prof-Kristel-Bernaerts">Kristel Bernaerts</a>. It was very useful,as she gave some good ideas on how to get started. We will have to design a way of expressing the enzyme(s) in a cyclic manner and a possible way to achieve this would be by expressing it during a short pulse, stimulated by the presence of a signal from other cells.<br/>
+
We are excited to get our modelling  team going and therefore we had an appointment with professor <a href="http://www.kuleuven.be/optec/people/74-Prof-Kristel-Bernaerts">Kristel Bernaerts</a>. It was very useful, as she gave some good ideas on how to get started. We will have to design a way of expressing the enzyme(s) in a cyclic manner and a possible way to achieve this would be by expressing it during a short pulse, stimulated by the presence of a signal from other cells.<br/>
And to tackle this problem we also started brainstorming about some possible networks with oscillating behaviour. Here we have to keep in mind that the bacterial colonies have to get and stay synchronized and this might be achieved by a rapid (protein-protein interaction) feedback which is proportional to the phase difference.
And to tackle this problem we also started brainstorming about some possible networks with oscillating behaviour. Here we have to keep in mind that the bacterial colonies have to get and stay synchronized and this might be achieved by a rapid (protein-protein interaction) feedback which is proportional to the phase difference.
In an attempt to collaborate with the wetlab team we discussed our main focus for the upcoming weeks: figuring out an oscillating construct and simulating the behaviour of the methyl salicylate BioBrick.<br/>
In an attempt to collaborate with the wetlab team we discussed our main focus for the upcoming weeks: figuring out an oscillating construct and simulating the behaviour of the methyl salicylate BioBrick.<br/>
-
And since organization comes after thought Bert and Sander started  working on simulating the Mortier Oscillator (as Bert likes to call it) in SimBiology, while Tina and Tomas started modelling the network to produce methyl salicylate by looking at BSMT1(the enzyme that catalyzes the reaction from salicylate to methyl salicylate) and at pchA and pchB which catalyze the reactions from chorismate to isochorismate to salicylate.<br/>
+
And since organization comes after thought Bert and Sander started  working on simulating the Mortier Oscillator (as Bert likes to call it) in SimBiology, while Tina and Tomas started modelling the network to produce methyl salicylate by looking at BSMT1 (the enzyme that catalyzes the reaction from salicylate to methyl salicylate) and at pchA and pchB which catalyze the reactions from chorismate to isochorismate to salicylate.<br/>
Later this week Tina had a look on how to use the COBRA toolbox for our purpose and we would have to check the constraints on reaction rates and add our new reactions.<br/>
Later this week Tina had a look on how to use the COBRA toolbox for our purpose and we would have to check the constraints on reaction rates and add our new reactions.<br/>
-
Bert contacted professor <a href="http://www.kuleuven.be/wieiswie/en/person/u0015385">Johan Suykens</a> and professor <a href="http://www.kuleuven.be/wieiswie/nl/person/00005764">Jan Degrève</a> to get an idea of how to analyse the MO. He gave the advice towards a bifurcation analysis, on which Bert will have a closer look.<br/>Tomas on the other hand started on the development of another oscillator and contacted professor href="http://www.kuleuven.be/wieiswie/en/person/u0006989">Dirk Roose</a> of mathematical engineering, who suggested the use of MatCont.Furthermore we are also considering a strategy of a knock-out/knock-down of genes responsible for chorismate mutase and chorismate lyase to optimize methyl salicylate production.<br/>
+
Bert contacted professor <a href="http://www.kuleuven.be/wieiswie/en/person/u0015385">Johan Suykens</a> and professor <a href="http://www.kuleuven.be/wieiswie/nl/person/00005764">Jan Degrève</a> to get an idea of how to analyse the MO. He gave the advice towards a bifurcation analysis, on which Bert will have a closer look.<br/>Tomas on the other hand started on the development of another oscillator and contacted professor <a href="http://www.kuleuven.be/wieiswie/en/person/u0006989">Dirk Roose</a> of mathematical engineering, who suggested the use of MatCont. Furthermore we are also considering a strategy of a knock-out/knock-down of genes responsible for chorismate mutase and chorismate lyase to optimize methyl salicylate production.<br/>
</p> </div>
</p> </div>
   </div>
   </div>
Line 37: Line 142:
     <div class="content">
     <div class="content">
     <p align="justify">This week the AutoBioCAD software is keeping us busy! <br/>
     <p align="justify">This week the AutoBioCAD software is keeping us busy! <br/>
-
Besides this we started investigating the literature of ecology and looking into a better model for our quorum sensing system using PDEs. We chose aphid and crop as our model species and now we’ll try to model aphid reproduction together with the influence of bèta-farnesene and methyl salicylate on life cycle and movement. We’ll have to make an estimation of the damage/loss of crops due to a certain population of aphids.<br/>
+
Besides this we started investigating the literature of ecology and looking into a better model for our quorum sensing system using PDEs. We chose aphid and crop as our model species and now we’ll try to model aphid reproduction together with the influence of EBF and methyl salicylate on life cycle and movement. We’ll have to make an estimation of the damage/loss of crops due to a certain population of aphids.<br/>
Meanwhile Bert added spatial heterogeneity to his Mortier Oscillator to see what happens on the population level.<br/>
Meanwhile Bert added spatial heterogeneity to his Mortier Oscillator to see what happens on the population level.<br/>
(For some more mathematical information we would like to refer you to : Murray - Mathematical Biology)<br/>
(For some more mathematical information we would like to refer you to : Murray - Mathematical Biology)<br/>
-
We also had a meeting with Jonas Demeulemeester this week, who explained us the methods the modelling team of 2008 used. It was a lot of help to us and now we are off for a new start! Sander started to complete our model in SimBiology, integrating the DNA, RNA, protein and metabolite level. Tina will look further into the application of FBA and COBRA in our project. Tomas will continue modelling the diffusion of the pheromones and their effect on the insects and was finally able to solve the 2D PDE in Matlab, he refined the PDE and boundary conditions for the diffusion. A more realistic model for the wind speed was now used.<br/>
+
We also had a meeting with Jonas Demeulemeester this week, who explained us the methods <a href="https://2008.igem.org/Team:KULeuven/Model/Overview"> the modelling team of 2008 </a> used. It was a lot of help to us and now we are off to a new start! Sander started to complete our model in SimBiology, integrating the DNA, RNA, protein and metabolite level. Tina will look further into the application of FBA and COBRA in our project. Tomas will continue modelling the diffusion of the pheromones and their effect on the insects and was finally able to solve the 2D PDE in MATLAB, he refined the PDE and boundary conditions for the diffusion. A more realistic model for the wind speed was now used.<br/>
We also had a small meeting with our wetlabteam colleagues to look for a possible integration of our modelling approaches in the wetlab approach.<br/>
We also had a small meeting with our wetlabteam colleagues to look for a possible integration of our modelling approaches in the wetlab approach.<br/>
And in between all this we started an extensive literature study in the hope to find the right kinetic paramaters for our model in SimBiology.<br/>
And in between all this we started an extensive literature study in the hope to find the right kinetic paramaters for our model in SimBiology.<br/>
-
We also discussed the possibility and usefulness of a q-PCR experiment and decided that Sander will perform the experiment next week at the CMPG (centre for microbiological and plant genetics)<br/>
+
We also discussed the possibility and usefulness of a qPCR experiment and decided that Sander will perform the experiment next week at the <a href="http://www.biw.kuleuven.be/dtp/cmpg/">CMPG</a> (centre for microbiological and plant genetics)<br/>
</p> </div>
</p> </div>
Line 55: Line 160:
     <h3>Week 3: Meeting professionals and digging into the literature</h3>
     <h3>Week 3: Meeting professionals and digging into the literature</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">This week we continued the literature search for optimal transcription and degradation rates.
+
     <p align="justify">This week we continued the literature search for optimal transcription and degradation rates  
and the exploration of the COBRA toolbox. <br/>
and the exploration of the COBRA toolbox. <br/>
We also looked at the possible implementation of COBRA in the general approach of our project.
We also looked at the possible implementation of COBRA in the general approach of our project.
-
Tina tried to implement our model in COBRA by adding the different parameters for the LB medium, found in the literature as a list of extracellular metabolites available in the LB medium for uptake but there are still some questions left unanswered regarding the COBRA toolbox, therefore we contacted Professor Bernaerts again,but unfortunately she is out of office until the 5th of august. <br/>
+
Tina tried to implement our model in COBRA by adding the different parameters for the LB medium, found in the literature as a list of extracellular metabolites available in the LB medium for uptake but there are still some questions left unanswered regarding the COBRA toolbox, therefore we contacted Professor <a href="http://www.kuleuven.be/optec/people/74-Prof-Kristel-Bernaerts">Kristel Bernaerts</a> again, but unfortunately she is out of office until the 5th of august. <br/>
Therefore we had a look at what different approaches we could take while putting the COBRA part on a low profile for a while until her return.<br/>
Therefore we had a look at what different approaches we could take while putting the COBRA part on a low profile for a while until her return.<br/>
-
       We looked into motiontracking for which we contacted a person skilled in the art, Tim Tambuyzer. To gather some information on the nCounter Technology we contacted Nucleomics Core and scheduled a meeting with dr. Allemeersch, the ncounter technology was not useful for our project though since it would be rather an overkill, finally we also had a meeting with Prof. Aertsen of the centre for Food and Microbial Technology about implementation of networks in vivo.
+
       We looked into motiontracking for which we contacted a person skilled in the art, <a href="http://www.kuleuven.be/wieiswie/nl/person/u0073220">Tim Tambuyzer</a> . To gather some information on the nCounter Technology we contacted Nucleomics Core and scheduled a meeting with dr.<a href="http://www.kuleuven.be/wieiswie/nl/person/u0038989">Joke Allemeersch</a>, the nCounter technology was not useful for our project though, since it would be rather an overkill, finally we also had a meeting with professor <a href="http://www.kuleuven.be/wieiswie/en/person/u0033548">Abram Aertsen</a> of the centre for Food and Microbial Technology about implementation of networks <i>in vivo</i>.
</p> </div>
</p> </div>
   </div>
   </div>
Line 70: Line 175:
   <div class="cbp_tmicon cbp_tmicon-phone"><i class="livicon activeicon" data-name="dashboard" data-size="23" data-onparent="true" data-color="white"></i></div>
   <div class="cbp_tmicon cbp_tmicon-phone"><i class="livicon activeicon" data-name="dashboard" data-size="23" data-onparent="true" data-color="white"></i></div>
   <div class="cbp_tmlabel">
   <div class="cbp_tmlabel">
-
     <h3>Week 4: Computational superspeed and oldschool Matlabcoding</h3>
+
     <h3>Week 4: Computational superspeed and oldschool MATLAB coding</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">Tomas is back on track to finetune his model for diffusion, while Bert is thinking about the correct way to verify the model. He is writing a matlab code in order to verify the oscillator.<br/>
+
     <p align="justify">Tomas is back on track to finetune his model for diffusion, while Bert is thinking about the correct way to verify the model. He is writing a MATLAB code in order to verify the oscillator.<br/>
-
Tina on the other hand is looking deeper into GapFilling of our model in the COBRA toolbox while Tomas got an account to the HPC cluster and 1000 credits. As it Seems waiting for a MATLAB job to be done is something from the past!<br/>
+
Tina on the other hand is looking deeper into GapFilling of our model in the COBRA toolbox while Tomas got an account to the HPC cluster and 1000 credits. As it seems waiting for a MATLAB job to be done is something from the past!<br/>
-
Later this week Bert and Tina had a meeting with Tim Tambuyzer, to talk about the options of videotracking and with prof. dr. ir. Dirk Roose to talk about possible analysis methods for the oscillator model. Meanwhile Sander was shedding blood, sweat and tears on the q-PCR!<br/>
+
Later this week Bert and Tina had a meeting with <a href="http://www.kuleuven.be/wieiswie/nl/person/u0073220">Tim Tambuyzer</a>, to talk about the options of videotracking and with professor <a href="http://www.kuleuven.be/wieiswie/en/person/u0006989">Dirk Roose</a> to talk about possible analysis methods for the oscillator model. Meanwhile Sander was shedding blood, sweat and tears on the qPCR!<br/>
-
Besides this we already started rocking some wiki texts.The pages are getting filled!
+
Besides this we already started rocking some wiki texts. The pages are getting filled!
</p> </div>
</p> </div>
Line 87: Line 192:
     <h3>Week 5: Awfully long wiki text alert!</h3>
     <h3>Week 5: Awfully long wiki text alert!</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">Our days have started to look very similiar now, so besides the usual stuff we finally had a meeting with prof. Bernaerts again to tackle some COBRA problems and she helped us back on track!<br/>
+
     <p align="justify">Our days have started to look very similiar now, so besides the usual stuff we finally had a meeting with professor <a href="http://www.kuleuven.be/optec/people/74-Prof-Kristel-Bernaerts">Kristel Bernaerts</a>  again to tackle some COBRA problems and she helped us back on track!<br/>
Bert started writing a second text regarding the Oscillator analysis while Tina is becoming super best friends with COBRA and Sander is preparing for the presentation of tomorrow.<br/>
Bert started writing a second text regarding the Oscillator analysis while Tina is becoming super best friends with COBRA and Sander is preparing for the presentation of tomorrow.<br/>
-
Later this week Sander had a meeting with prof. Engelborghs to get some clarification about enzyme kinetics for the parameter model and Tomas is back on track for his ecological model and installed COMSOL for this purpose. Bert is also still writing awfully long wiki texts by the way.<br/>
+
While Tomas is back on track for his ecological model and installed COMSOL for this purpose, Bert is writing awfully long wiki texts.<br/>
Line 103: Line 208:
     <h3>Week 6: No obligatory days off for the modelling team!</h3>
     <h3>Week 6: No obligatory days off for the modelling team!</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">On the 18th of august, the KU Leuven was closed because of a national holiday. But this couldn’t stop us from working! We got interviewed by the ethics team! We had some philosophical reflections on ourselves and our project with our Philosophy students Michiel and Sam.<br/>
+
     <p align="justify">On the 15th of august, the KU Leuven was closed because of a national holiday. But this couldn’t stop us from working! We got interviewed by the ethics team! We had some philosophical reflections on ourselves and our project with our philosophy students Michiel and Sam.<br/>
-
When we got back to modelling, Tomas took a look at the supercomputer framework in order to parallelize the calculations needed to analyze the oscillator. Bert is looking in the literature for good parameter ranges.
+
When we got back to modelling, Tomas took a look at the supercomputer framework in order to parallelise the calculations needed to analyse the oscillator. Bert is looking in the literature for good parameter ranges.
Sander continued on the parameters model.<br/>
Sander continued on the parameters model.<br/>
Line 115: Line 220:
   <div class="cbp_tmicon cbp_tmicon-phone"><i class="livicon activeicon" data-name="laptop" data-size="23" data-onparent="true" data-color="white"></i></div>
   <div class="cbp_tmicon cbp_tmicon-phone"><i class="livicon activeicon" data-name="laptop" data-size="23" data-onparent="true" data-color="white"></i></div>
   <div class="cbp_tmlabel">
   <div class="cbp_tmlabel">
-
     <h3>Week 7: Wikitexts,modelling,wikitexts...</h3>
+
     <h3>Week 7: Wikitexts, modelling, wikitexts...</h3>
     <div class="content">
     <div class="content">
     <p align="justify">COMSOL modelling seems to work out and realistic parameters have been found.<br/>
     <p align="justify">COMSOL modelling seems to work out and realistic parameters have been found.<br/>
Tina started writing a wikitext regarding the COBRA analyses while Sander found out that doing a qPCR on plasmids is hard as it’s not easy to get rid of plasmid DNA contamination, so we’re working on a solution.<br/>
Tina started writing a wikitext regarding the COBRA analyses while Sander found out that doing a qPCR on plasmids is hard as it’s not easy to get rid of plasmid DNA contamination, so we’re working on a solution.<br/>
-
The wiki text about Ecology and the matlab code for the biggest part of the verification of the oscillator are near completion.
+
The wiki text about Ecology and the MATLAB code for the biggest part of the verification of the oscillator are near completion.
</p> </div>
</p> </div>
Line 131: Line 236:
     <h3>Week 8: Tomas and Tina traded Bert and Sander for some insects (the difference is insignificant)</h3>
     <h3>Week 8: Tomas and Tina traded Bert and Sander for some insects (the difference is insignificant)</h3>
     <div class="content">
     <div class="content">
-
     <p align="justify">Tomas and Tina are spending a week at PCFruit to perform the insect experiments while Sander and Bert continue the good old fashioned modelling work! The finish line is in sight! The future will be dedicated to writing wiki texts.
+
     <p align="justify">Tomas and Tina are spending a week at pcfruit to perform the insect experiments while Sander and Bert continue the good old fashioned modelling work! The finish line is in sight! The future will be dedicated to writing wiki texts.
</p> </div>
</p> </div>
   </div>
   </div>

Latest revision as of 03:36, 29 October 2013

iGem

Secret garden

Congratulations! You've found our secret garden! Follow the instructions below and win a great prize at the World jamboree!


  • A video shows that two of our team members are having great fun at our favourite company. Do you know the name of the second member that appears in the video?
  • For one of our models we had to do very extensive computations. To prevent our own computers from overheating and to keep the temperature in our iGEM room at a normal level, we used a supercomputer. Which centre maintains this supercomputer? (Dutch abbreviation)
  • We organised a symposium with a debate, some seminars and 2 iGEM project presentations. An iGEM team came all the way from the Netherlands to present their project. What is the name of their city?

Now put all of these in this URL:https://2013.igem.org/Team:KU_Leuven/(firstname)(abbreviation)(city), (loose the brackets and put everything in lowercase) and follow the very last instruction to get your special jamboree prize!

tree ladybugcartoon

Modelling

You are here!

Ecology

Plant and insect experiments

E-β-farnesene

Work about our EBF-biobrick

Feed forward loop

Implementation of the oscillator in vivo

Methyl salicylate

Work about our MeS-biobrick

qPCR

Finding wetlab data for the models

This is the journal of the modelling team, they designed and analysed an oscillator, delved into Cobra, modelled the production system of methyl salicylate and even found the time to start a model about how our products will be dispersed on a field over time.

  • Week 1: The modelling team is on fire!

    The modelling team is ready to rock MATLAB after watching some tutorials and webinars about modelling and Symbiology
    We are excited to get our modelling team going and therefore we had an appointment with professor Kristel Bernaerts. It was very useful, as she gave some good ideas on how to get started. We will have to design a way of expressing the enzyme(s) in a cyclic manner and a possible way to achieve this would be by expressing it during a short pulse, stimulated by the presence of a signal from other cells.
    And to tackle this problem we also started brainstorming about some possible networks with oscillating behaviour. Here we have to keep in mind that the bacterial colonies have to get and stay synchronized and this might be achieved by a rapid (protein-protein interaction) feedback which is proportional to the phase difference. In an attempt to collaborate with the wetlab team we discussed our main focus for the upcoming weeks: figuring out an oscillating construct and simulating the behaviour of the methyl salicylate BioBrick.
    And since organization comes after thought Bert and Sander started working on simulating the Mortier Oscillator (as Bert likes to call it) in SimBiology, while Tina and Tomas started modelling the network to produce methyl salicylate by looking at BSMT1 (the enzyme that catalyzes the reaction from salicylate to methyl salicylate) and at pchA and pchB which catalyze the reactions from chorismate to isochorismate to salicylate.
    Later this week Tina had a look on how to use the COBRA toolbox for our purpose and we would have to check the constraints on reaction rates and add our new reactions.
    Bert contacted professor Johan Suykens and professor Jan Degrève to get an idea of how to analyse the MO. He gave the advice towards a bifurcation analysis, on which Bert will have a closer look.
    Tomas on the other hand started on the development of another oscillator and contacted professor Dirk Roose of mathematical engineering, who suggested the use of MatCont. Furthermore we are also considering a strategy of a knock-out/knock-down of genes responsible for chorismate mutase and chorismate lyase to optimize methyl salicylate production.

  • Week 2: From the Mortier Oscillator to COBRA and back

    This week the AutoBioCAD software is keeping us busy!
    Besides this we started investigating the literature of ecology and looking into a better model for our quorum sensing system using PDEs. We chose aphid and crop as our model species and now we’ll try to model aphid reproduction together with the influence of EBF and methyl salicylate on life cycle and movement. We’ll have to make an estimation of the damage/loss of crops due to a certain population of aphids.
    Meanwhile Bert added spatial heterogeneity to his Mortier Oscillator to see what happens on the population level.
    (For some more mathematical information we would like to refer you to : Murray - Mathematical Biology)
    We also had a meeting with Jonas Demeulemeester this week, who explained us the methods the modelling team of 2008 used. It was a lot of help to us and now we are off to a new start! Sander started to complete our model in SimBiology, integrating the DNA, RNA, protein and metabolite level. Tina will look further into the application of FBA and COBRA in our project. Tomas will continue modelling the diffusion of the pheromones and their effect on the insects and was finally able to solve the 2D PDE in MATLAB, he refined the PDE and boundary conditions for the diffusion. A more realistic model for the wind speed was now used.
    We also had a small meeting with our wetlabteam colleagues to look for a possible integration of our modelling approaches in the wetlab approach.
    And in between all this we started an extensive literature study in the hope to find the right kinetic paramaters for our model in SimBiology.
    We also discussed the possibility and usefulness of a qPCR experiment and decided that Sander will perform the experiment next week at the CMPG (centre for microbiological and plant genetics)

  • Week 3: Meeting professionals and digging into the literature

    This week we continued the literature search for optimal transcription and degradation rates and the exploration of the COBRA toolbox.
    We also looked at the possible implementation of COBRA in the general approach of our project. Tina tried to implement our model in COBRA by adding the different parameters for the LB medium, found in the literature as a list of extracellular metabolites available in the LB medium for uptake but there are still some questions left unanswered regarding the COBRA toolbox, therefore we contacted Professor Kristel Bernaerts again, but unfortunately she is out of office until the 5th of august.
    Therefore we had a look at what different approaches we could take while putting the COBRA part on a low profile for a while until her return.
    We looked into motiontracking for which we contacted a person skilled in the art, Tim Tambuyzer . To gather some information on the nCounter Technology we contacted Nucleomics Core and scheduled a meeting with dr.Joke Allemeersch, the nCounter technology was not useful for our project though, since it would be rather an overkill, finally we also had a meeting with professor Abram Aertsen of the centre for Food and Microbial Technology about implementation of networks in vivo.

  • Week 4: Computational superspeed and oldschool MATLAB coding

    Tomas is back on track to finetune his model for diffusion, while Bert is thinking about the correct way to verify the model. He is writing a MATLAB code in order to verify the oscillator.
    Tina on the other hand is looking deeper into GapFilling of our model in the COBRA toolbox while Tomas got an account to the HPC cluster and 1000 credits. As it seems waiting for a MATLAB job to be done is something from the past!
    Later this week Bert and Tina had a meeting with Tim Tambuyzer, to talk about the options of videotracking and with professor Dirk Roose to talk about possible analysis methods for the oscillator model. Meanwhile Sander was shedding blood, sweat and tears on the qPCR!
    Besides this we already started rocking some wiki texts. The pages are getting filled!

  • Week 5: Awfully long wiki text alert!

    Our days have started to look very similiar now, so besides the usual stuff we finally had a meeting with professor Kristel Bernaerts again to tackle some COBRA problems and she helped us back on track!
    Bert started writing a second text regarding the Oscillator analysis while Tina is becoming super best friends with COBRA and Sander is preparing for the presentation of tomorrow.
    While Tomas is back on track for his ecological model and installed COMSOL for this purpose, Bert is writing awfully long wiki texts.

  • Week 6: No obligatory days off for the modelling team!

    On the 15th of august, the KU Leuven was closed because of a national holiday. But this couldn’t stop us from working! We got interviewed by the ethics team! We had some philosophical reflections on ourselves and our project with our philosophy students Michiel and Sam.
    When we got back to modelling, Tomas took a look at the supercomputer framework in order to parallelise the calculations needed to analyse the oscillator. Bert is looking in the literature for good parameter ranges. Sander continued on the parameters model.

  • Week 7: Wikitexts, modelling, wikitexts...

    COMSOL modelling seems to work out and realistic parameters have been found.
    Tina started writing a wikitext regarding the COBRA analyses while Sander found out that doing a qPCR on plasmids is hard as it’s not easy to get rid of plasmid DNA contamination, so we’re working on a solution.
    The wiki text about Ecology and the MATLAB code for the biggest part of the verification of the oscillator are near completion.

  • Week 8: Tomas and Tina traded Bert and Sander for some insects (the difference is insignificant)

    Tomas and Tina are spending a week at pcfruit to perform the insect experiments while Sander and Bert continue the good old fashioned modelling work! The finish line is in sight! The future will be dedicated to writing wiki texts.