Team:Imperial College London/Modelling

From 2010.igem.org

(Difference between revisions)
(Undo revision 208542 by Anita (Talk))
 
(23 intermediate revisions not shown)
Line 5: Line 5:
|-
|-
-
|In the process of designing our construct two major questions arose which could be answered by computer modelling:
+
|During the design of our construct three major questions arose which could be answered by computer modelling:
 +
 
<ol>
<ol>
-
<li><b>Output Amplification Model</b><br/>We came up with an idea of using the amplification of a colour output to make it show within minutes after the stimulus has been added. The question that arose was whether amplification will actually perform better than simple production in the cellular environment. Furthermore, we had trouble deciding whether we should design the amplification module to consist of 1,2 or even more amplification steps. These issues seemed to be difficult enough to employ modelling.</li>
+
 
-
<li><b>Signalling Module Model</b><br/>We decided to use the ComCDE signalling pathway from ''S.pneumoniae'' and so  questions arose on whether it would work appropriately in ''B.subtilis''. We modelled this system to make sure that the signalling pathway would be working as anticipated.
+
{| style="width:800px;background:#f5f5f5;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;"
-
</li>
+
|1. <b>Detection Model</b>
-
<li><b>Surface Protein Model</b><br/>We came up with a novel idea of detecting organisms that we do not have a specific receptor for. In our particular example, the protease of Schistosoma was meant to cleave a protein displayed on the bacteria's cell wall. The cleaved peptide was supposed to be recognized by the receptor which would activate the colour expression. This solution raised questions about the risk of false positive or whether there are any chances for ComD receptors to be activated in the diluted environment. Modelling of this module would answer these questions.</li>
+
|style="font-family: helvetica, arial, sans-serif;font-size:1.2em;color:#ea8828;" align="right"|[[Team:Imperial_College_London/Modelling/Protein_Display/Objectives|Objectives]] | [[Team:Imperial_College_London/Modelling/Protein_Display/Detailed_Description|Description]] | [[Team:Imperial_College_London/Modelling/Protein_Display/Parameters_and_Constants|Constants]] | [[Team:Imperial_College_London/Modelling/Protein_Display/Results_and_Conclusion|Results]] | [[Team:Imperial_College_London/Modelling/Protein_Display/Download_MatLab_Files|MATLAB Code]]
 +
|}
 +
We came up with a novel idea of detecting organisms that we do not have a specific receptor for. In our particular example, the protease of Schistosoma is meant to cleave a protein displayed on the bacteria's cell wall (<i>the pink rectangle on the diagram below</i>). The cleaved peptide is recognized by the receptor which would then activate the colour expression. This solution raised questions about the risk of false positive or whether there are any chances for ComD receptors to be activated in the diluted environment. Modelling of this module would enable us to answer these questions.<br/><br/>
 +
{| style="width:800px;background:#f5f5f5;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;"
 +
|2. <b>Signaling Model</b>
 +
|style="font-family: helvetica, arial, sans-serif;font-size:1.2em;color:#ea8828;" align="right"|[[Team:Imperial_College_London/Modelling/Signalling/Objectives|Objectives]] | [[Team:Imperial_College_London/Modelling/Signalling/Detailed_Description|Description]] | [[Team:Imperial_College_London/Modelling/Signalling/Parameters_and_Constants|Constants]] | [[Team:Imperial_College_London/Modelling/Signalling/Results_and_Conclusion|Results]] | [[Team:Imperial_College_London/Modelling/Signalling/Download_MatLab_Files|MATLAB Code]]
 +
|}
 +
We decided to use the ComCDE signaling pathway from ''S.pneumoniae'' and so  questions arose on whether it would work appropriately in ''B.subtilis''. We modelled this system to make sure that the signalling pathway would be working as anticipated (<i>the green rectangle on the diagram below</i>).<br/><br/>
 +
{| style="width:800px;background:#f5f5f5;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;"
 +
|3. <b>Fast Response Model</b>
 +
|style="font-family: helvetica, arial, sans-serif;font-size:1.2em;color:#ea8828;" align="right"|[[Team:Imperial_College_London/Modelling/Output/Objectives|Objectives]] | [[Team:Imperial_College_London/Modelling/Output/Detailed_Description|Description]] | [[Team:Imperial_College_London/Modelling/Output/Parameters_and_Constants|Constants]] | [[Team:Imperial_College_London/Modelling/Output/Results_and_Conclusion|Results]] | [[Team:Imperial_College_London/Modelling/Output/Download_MatLab_Files|MATLAB Code]]
 +
|}
 +
We came up with the idea of using the amplification of a colour output, which would show within minutes of the stimulus being added (<i>the blue rectangle on the diagram below</i>). The question that arose was whether amplification will actually perform better than simple production (i.e. transcription and translation) in the cellular environment. Furthermore, we had difficulties deciding whether we should design the amplification module to consist of 1,2 or even more amplification steps. These issues were quite difficult to answer, so we decided to employ modelling.<br/>
</ol>
</ol>
 +
<div ALIGN=CENTER>
 +
{| style="width:404px;background:#e7e7e7;text-align:center;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5";
 +
|-
 +
|[[Image:IC_Overview.jpg|400px]]
 +
|-
 +
|Each of the shadowed rectangles corresponds to the model of a particular module of the system. Pink corresponds to the model of the detection mechanism, green corresponds to the model of the signal transduction and blue corresponds to the model of the fast response module.
 +
|}
 +
</div>
|}
|}
Line 19: Line 40:
|-
|-
-
|<b>Output Amplification Model</b><br/>
+
|<b>Detection Model</b><br/>
<ol>
<ol>
-
  <li>It was shown that amplified systems easily outperform the simple production system (control)
+
  <li>We determined initial TEV protease concentrations which would result in the optimal activation of the receptor. This optimal activation would happen within 1.5 minutes after elastases come into contact with our cell.</li>
-
<li>It was concluded that there is no advantage of 3-step amplification over 2-step amplification. Therefore, the design of a 3-step amplifier was abandoned.</li>
+
-
<li>The results concerning the 2-step amplification module were not conclusive. It could not be firmly decided whether 2-step amplification is going to perform better than 1-step amplification. This is because several of the parameters that 2- and 1-step amplifiers are sensitive to could not be determined with certainty. 2 parameters have been recognised as crucial and decisive. Those are protein production rates and catalytic constants of the enzymes.</li>
+
-
<li>Hence, the conditions for effective amplification were determined.</li>
+
</ol>
</ol>
-
 
<div ALIGN=CENTER>
<div ALIGN=CENTER>
  {| style="background:#e7e7e7;text-align:center;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5";
  {| style="background:#e7e7e7;text-align:center;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5";
  |-
  |-
-
  |[[Image:Output_model.png|450px]]
+
  |[[Image:IC_AIP_threshold.png|600px]]
  |-
  |-
-
  |Concentration of coloured compound.
+
  |Notice log-log scale.
  |}
  |}
</div>
</div>
 +
|-
 +
|<b>Signaling Model</b><br/>
 +
Even though our model of the signaling module is more simplistic than the real life situation, it provided very important results. We were able to determine under which conditions the signaling pathway would be working and could obtain the major constraints of our system. These constraints are that the necessary concentrations for ComD and AIP are reached before signal transduction is started.
-
<b>Signalling Module Model</b><br/>
+
In this model, we will treat phosphorylated ComE (referred to as ComE*) as our "output". Phosphorylated ComE is a transcription factor, which will bind to the DNA to enable transcription.  
-
Even though our model of the signalling module is more simplistic than the real life situation, it provided very important results. We were able to determine under which conditions the signalling pathway would be working and could obtain the major constraints of our system. These constraints are that the necessary concentrations for ComD and AIP are reached before signal transduction is started.  
+
<div ALIGN=CENTER>
<div ALIGN=CENTER>
Line 44: Line 63:
  |[[Image:IC_Signalling_Results1.png|450px]]
  |[[Image:IC_Signalling_Results1.png|450px]]
  |-
  |-
-
  |Graph showing how [ComE]<html><sub>final</sub> eventually reaches the value 5&times;10<sup>-11</sup>M</html>.  
+
  |Graph showing how the concentration of the transcription
 +
|-
 +
|factor (phosphorylated ComE) eventually reaches the value 5&times;10<sup>-11</sup>M.  
|}
|}
</div>
</div>
-
 
+
|-
-
<b>Surface Protein Model</b><br/>
+
|<b>Fast Response Model</b><br/>
<ol>
<ol>
-
  <li>Initial TEV protease concentrations we determined for the optimal activation of the receptor within 1.5 minutes after elastases would have come into contact with our cell.</li>
+
  <li>It was shown that our amplification systems easily outperform the simple production system. Our system can respond within several minutes rather than several hours, which is the case for many currently used reporters.
 +
<li>It was concluded that there is no advantage of 3-step amplification over 2-step amplification. Therefore, the design of a 3-step amplifier was abandoned.</li>
 +
<li>The results concerning the 2-step amplification module were not totally conclusive. It could not be firmly decided whether 2-step amplification is going to perform better than 1-step amplification. This is because several of the parameters that 2- and 1-step amplifiers are sensitive to could not be determined with certainty. Two parameters have been recognised as crucial and decisive. These are the protein production rates and the catalytic constants of the enzymes.</li>
 +
<li>Hence, the conditions for effective amplification were determined (more details are revealed in the Fast Response Model [[Team:Imperial_College_London/Modelling/Output/Results_and_Conclusion|Results]] page).</li>
</ol>
</ol>
 +
<div ALIGN=CENTER>
<div ALIGN=CENTER>
-
  {| style="background:#e7e7e7;text-align:center;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5";
+
  {| style="width:454px;background:#e7e7e7;text-align:center;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5";
  |-
  |-
-
  |[[Image:IC_AIP_threshold.png|450px]]
+
  |[[Image:Output_model.png|450px]]
  |-
  |-
-
  |Graph showing when threshold AIP concentration is reached <br/>(for different initial TEV concentrations). Notice log-log scale.
+
  |Graph showing how our 2-step amplification system outperforms 1-step amplification. Note that time=0 corresponds to when the first transcription happens upon binding of the transcription factor to DNA.
  |}
  |}
</div>
</div>
Line 63: Line 88:
{| style="width:900px;background:#f5f5f5;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;margin-top:5px;" cellspacing="20"
{| style="width:900px;background:#f5f5f5;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;margin-top:5px;" cellspacing="20"
-
|style="font-family: helvetica, arial, sans-serif;font-size:2em;color:#ea8828;"|Quick overview of models
+
|style="font-family: helvetica, arial, sans-serif;font-size:2em;color:#ea8828;"|Quick overview of the models
|-
|-
-
|<b>Output Amplification Model</b><br/>
+
|<b>Detection Model</b><br/>
Goals:
Goals:
-
<ol><p>This model was mainly developed in order to determine whether simple production is better than 1-, 2- or 3-step amplification.</p>Furthermore, an estimation of the speed of the response was desirable.
+
<ol><p>The aim of this model is to determine the concentration of Schistosoma elastase or TEV protease that should be added to the bacteria in order to trigger a response.</p>
 +
<p>It was also attempted to model how long it takes for the protease or elastase to cleave the required amount of peptides to activate receptors.</p>
</ol>
</ol>
 +
Elements of the system:
Elements of the system:
-
<ol><li>Dioxygenase (<i>blue on the diagrams below</i>) is an enzyme that acts on catechol to produce a yellow output. In most of our models dioxygenase was treated as an output because it was found that active dioxygenase acting on catechol produces the coloured output within a split second.</li>
+
<ol><li>The surface protein consists of a cell wall binding domain, a linker and AIP (Auto Inducing Peptide)</li>
-
<li>GFP-Dioxygenase fusion protein (<i>GFP is shown green on the diagrams</i>). Dioxygenase joined by the linker to GFP was assumed to be inactive.</li>
+
<li>Schistosoma elastase (this is the enzyme released by the parasite) cleaves AIP from the cell wall binding domain at the linker site. In the laboratory, we used TEV protease as we could not obtain the Schistosoma elastase.</li>
-
<li>TEV protease (<i>pink on the diagrams below</i>) has the ability to cleave the GFP-Dioxygenase fusion protein, hence, it activates dioxygenase</li>
+
<li>The ComD receptor is activated (i.e. AIP concentration is high enough).</li>
-
<li>Split TEV protease (<i>purple on the diagrams below</i>) is an inactive split form of TEV mounted on coiled coils. It can be activated again by coiled coils being cleaved by another active TEV.</li>
+
</ol>
</ol>
-
 
-
<div ALIGN=CENTER>
 
-
{| style="background:#e7e7e7;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5" width=704px;
 
-
|-
 
-
|[[Image:Simple_Production.png|700px]]
 
-
|-
 
-
|Simple production upon activation of arbitrary colour output by transcription and translation indicated by the <i>blue arrow</i>.
 
-
|-
 
-
|[[Image:1-step_amplification.png|700px]]
 
-
|-
 
-
| Dioxygenase (C230) is simply produced. Upon activation at time t=0, it acts on catechol (cat.) to produce yellow output - muconic acid. Catechol is not shown to be produced by cell as it is added by person at arbitrary time.
 
-
|-
 
-
|[[Image:2-step_amplification.png|700px]]
 
-
|-
 
-
| The species that are shown in front of vertical line which indicates beginning of experiment mean that they have been accumulated beforehand in the cell. TEV protease activates inactive dioxygenase which acts on catechol to produce colour. 
 
-
|-
 
-
|[[Image:3-step_amplification.png|700px]]
 
-
|-
 
-
| This diagram introduces inactive split TEV protease attached to a coiled-coil as the third amplification step. Both inactive compounds have active site for TEV to activate tehm which results in multiple possibilities of action.
 
-
|}
 
-
</div>
 
Major assumptions:
Major assumptions:
<ol>
<ol>
<li>The chemical and enzymatic reactions are modelled according to the Law of Mass Action.</li>
<li>The chemical and enzymatic reactions are modelled according to the Law of Mass Action.</li>
-
<li>Our model assumes that the modelled system is inert within the bacterial body or that reactions with other species within the bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.</li>
+
<li>Our model assumes that the modelled system is inert within the bacterial cell or that reactions with other species within the bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.</li>
 +
<li>Due to our carefully chosen cell concentrations, the diffusion of free AIPs could be neglected. However, this restricts the model to the considered cell concentrations only.</li>
 +
<li>The threshold for receptor activation was defined by one specific value as opposed to considering intermediate states between fully "off" and "on".</li>
</ol>
</ol>
-
 
+
|-
-
<b>Signalling Modelling Model</b><br/>
+
|<b>Signaling Model</b><br/>
Goals:
Goals:
-
The aim of this model is to determine under which conditions the signalling transduction will happen in our bacteria.  
+
The aim of this model is to determine under which conditions the signaling transduction will happen in our bacteria.  
Elements of the system:
Elements of the system:
<ol>
<ol>
-
<li>The signalling model consists of ComD and ComE which are expressed by our bacteria, as well as AIP and Phosphate. These species are all assumed to be present in the cell at a sufficient concentration.</li>
+
<li>The signaling model consists of ComD and ComE, which are expressed by our bacteria, as well as AIP and Phosphate. These species are all assumed to be present in the cell at a sufficient concentration.</li>
<li>The ComD receptor is activated by the AIP. This triggeres phosphorylation of the ComD receptor. The Phosphate group of the ComD receptor then binds to ComE. The phosphorylated ComE binds to the DNA and acts as a transcription factor.</li>
<li>The ComD receptor is activated by the AIP. This triggeres phosphorylation of the ComD receptor. The Phosphate group of the ComD receptor then binds to ComE. The phosphorylated ComE binds to the DNA and acts as a transcription factor.</li>
</ol>
</ol>
Line 120: Line 127:
<li>Phosphorylation of the ComD receptor is modelled as an enzymatic reaction, neglecting the formation of an intermediate complex.</li>
<li>Phosphorylation of the ComD receptor is modelled as an enzymatic reaction, neglecting the formation of an intermediate complex.</li>
</ol>
</ol>
-
 
+
|-
-
<b>Surface Protein Model</b><br/>
+
|<b>Fast Response Model</b><br/>
Goals:
Goals:
-
<ol><p>The aim of this model is to determine the concentration of Schistosoma elastase or TEV protease that should be added to the bacteria in order to trigger a response. This would allow us to correlate the required concentration for the activation with the concentration of Schistosoma elastase in the lake.</p>
+
<ol><p>This model was mainly developed in order to determine whether simple production is better than 1-, 2- or 3-step amplification.</p>Furthermore, an estimation of the speed of the response was desirable.
-
<p>It was also attempted to model how long it takes for the protease or elastase to cleave the required amount of peptides.</p>
+
</ol>
</ol>
-
 
Elements of the system:
Elements of the system:
-
<ol><li>The surface protein consists of a cell wall binding domain, linker, AIP (Auto Inducing Peptide)</li>
+
<ol><li>Dioxygenase (<i>blue on the diagrams below</i>) is an enzyme that acts on catechol to produce a yellow output. In most of our models, dioxygenase was treated as an output because it was found that active dioxygenase acting on catechol produces the coloured output within a split second.</li>
-
<li>Schistosoma elastase (this is the enzyme released by the parasite) cleaves AIP from the cell wall binding domain at the linker site. In the laboratory we used TEV protease as we could not obtain the Schistosoma elastase.</li>
+
<li>GFP-Dioxygenase fusion protein (<i>GFP is shown green on the diagrams</i>). Dioxygenase, which is joined to GFP  by the linker, was assumed to be inactive.</li>
-
<li>The ComD receptor is activated (i.e. AIP concentration is high enough).</li>
+
<li>TEV protease (<i>pink on the diagrams below</i>) has the ability to cleave the GFP-Dioxygenase fusion protein, hence, it activates dioxygenase</li>
 +
<li>Split TEV protease (<i>purple on the diagrams below</i>) is an inactive split form of TEV mounted on coiled coils. It can be activated again by coiled coils being cleaved by another active TEV.</li>
</ol>
</ol>
 +
 +
<div ALIGN=CENTER>
 +
{| style="background:#e7e7e7;text-align:justify;font-family: helvetica, arial, sans-serif;color:#555555;margin- top:5px;padding: 2px;" cellspacing="5" width=704px;
 +
|-
 +
|[[Image:Simple_Production.png|700px]]
 +
|-
 +
|Simple production upon activation of an arbitrary colour output by transcription and translation indicated by the <i>blue arrow</i>.
 +
|-
 +
|[[Image:1-step_amplification.png|700px]]
 +
|-
 +
| Dioxygenase (C230) is simply produced. Upon activation at time t=0, it acts on catechol (cat.) to produce the yellow output - muconic acid. Catechol is not shown to be produced by the cell as it is added manually at an arbitrary time.
 +
|-
 +
|[[Image:2-step_amplification.png|700px]]
 +
|-
 +
| The species that are shown in front of the vertical line (which indicates the beginning of the experiment) indicate that they have been accumulated beforehand in the cell. TEV protease activates inactive dioxygenase which acts on catechol to produce colour. 
 +
|-
 +
|[[Image:3-step_amplification.png|700px]]
 +
|-
 +
| This diagram introduces inactive split TEV protease attached to a coiled-coil as the third amplification step. Both inactive compounds have active sites for TEV to activate them, which results in multiple possibilities of action.
 +
|}
 +
</div>
Major assumptions:
Major assumptions:
<ol>
<ol>
<li>The chemical and enzymatic reactions are modelled according to the Law of Mass Action.</li>
<li>The chemical and enzymatic reactions are modelled according to the Law of Mass Action.</li>
-
<li>Our model assumes that the modelled system is inert within the bacterial body or that reactions with other species within bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.</li>
+
<li>Our model assumes that the modelled system is inert within the bacterial cell or that reactions with other species within the bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.</li>
-
<li>Due to our carefully chosen cell concentrations, the diffusion of free AIPs could be neglected. However, this restricts the model to the considered cell concentrations only.</li>
+
-
<li>The threshold for receptor activation was defined by one specific value as opposed to considering intermediate states between fully "off" and "on".</li>
+
</ol>
</ol>
-
 
-
 
|}
|}

Latest revision as of 03:59, 28 October 2010

Modelling Overview | Detection Model | Signaling Model | Fast Response Model | Interactions
A major part of the project consisted of modelling each module. This enabled us to decide which ideas we should implement. Look at the Fast Response page for a great example of how modelling has made a major impact on our design!
Introduction to modelling
During the design of our construct three major questions arose which could be answered by computer modelling:


    1. Detection Model Objectives | Description | Constants | Results | MATLAB Code

    We came up with a novel idea of detecting organisms that we do not have a specific receptor for. In our particular example, the protease of Schistosoma is meant to cleave a protein displayed on the bacteria's cell wall (the pink rectangle on the diagram below). The cleaved peptide is recognized by the receptor which would then activate the colour expression. This solution raised questions about the risk of false positive or whether there are any chances for ComD receptors to be activated in the diluted environment. Modelling of this module would enable us to answer these questions.

    2. Signaling Model Objectives | Description | Constants | Results | MATLAB Code

    We decided to use the ComCDE signaling pathway from S.pneumoniae and so questions arose on whether it would work appropriately in B.subtilis. We modelled this system to make sure that the signalling pathway would be working as anticipated (the green rectangle on the diagram below).

    3. Fast Response Model Objectives | Description | Constants | Results | MATLAB Code

    We came up with the idea of using the amplification of a colour output, which would show within minutes of the stimulus being added (the blue rectangle on the diagram below). The question that arose was whether amplification will actually perform better than simple production (i.e. transcription and translation) in the cellular environment. Furthermore, we had difficulties deciding whether we should design the amplification module to consist of 1,2 or even more amplification steps. These issues were quite difficult to answer, so we decided to employ modelling.

IC Overview.jpg
Each of the shadowed rectangles corresponds to the model of a particular module of the system. Pink corresponds to the model of the detection mechanism, green corresponds to the model of the signal transduction and blue corresponds to the model of the fast response module.
Results & Conclusions
Detection Model
  1. We determined initial TEV protease concentrations which would result in the optimal activation of the receptor. This optimal activation would happen within 1.5 minutes after elastases come into contact with our cell.
IC AIP threshold.png
Notice log-log scale.
Signaling Model

Even though our model of the signaling module is more simplistic than the real life situation, it provided very important results. We were able to determine under which conditions the signaling pathway would be working and could obtain the major constraints of our system. These constraints are that the necessary concentrations for ComD and AIP are reached before signal transduction is started.

In this model, we will treat phosphorylated ComE (referred to as ComE*) as our "output". Phosphorylated ComE is a transcription factor, which will bind to the DNA to enable transcription.

IC Signalling Results1.png
Graph showing how the concentration of the transcription
factor (phosphorylated ComE) eventually reaches the value 5×10-11M.
Fast Response Model
  1. It was shown that our amplification systems easily outperform the simple production system. Our system can respond within several minutes rather than several hours, which is the case for many currently used reporters.
  2. It was concluded that there is no advantage of 3-step amplification over 2-step amplification. Therefore, the design of a 3-step amplifier was abandoned.
  3. The results concerning the 2-step amplification module were not totally conclusive. It could not be firmly decided whether 2-step amplification is going to perform better than 1-step amplification. This is because several of the parameters that 2- and 1-step amplifiers are sensitive to could not be determined with certainty. Two parameters have been recognised as crucial and decisive. These are the protein production rates and the catalytic constants of the enzymes.
  4. Hence, the conditions for effective amplification were determined (more details are revealed in the Fast Response Model Results page).
Output model.png
Graph showing how our 2-step amplification system outperforms 1-step amplification. Note that time=0 corresponds to when the first transcription happens upon binding of the transcription factor to DNA.
Quick overview of the models
Detection Model

Goals:

    The aim of this model is to determine the concentration of Schistosoma elastase or TEV protease that should be added to the bacteria in order to trigger a response.

    It was also attempted to model how long it takes for the protease or elastase to cleave the required amount of peptides to activate receptors.

Elements of the system:

  1. The surface protein consists of a cell wall binding domain, a linker and AIP (Auto Inducing Peptide)
  2. Schistosoma elastase (this is the enzyme released by the parasite) cleaves AIP from the cell wall binding domain at the linker site. In the laboratory, we used TEV protease as we could not obtain the Schistosoma elastase.
  3. The ComD receptor is activated (i.e. AIP concentration is high enough).

Major assumptions:

  1. The chemical and enzymatic reactions are modelled according to the Law of Mass Action.
  2. Our model assumes that the modelled system is inert within the bacterial cell or that reactions with other species within the bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.
  3. Due to our carefully chosen cell concentrations, the diffusion of free AIPs could be neglected. However, this restricts the model to the considered cell concentrations only.
  4. The threshold for receptor activation was defined by one specific value as opposed to considering intermediate states between fully "off" and "on".
Signaling Model

Goals: The aim of this model is to determine under which conditions the signaling transduction will happen in our bacteria.

Elements of the system:

  1. The signaling model consists of ComD and ComE, which are expressed by our bacteria, as well as AIP and Phosphate. These species are all assumed to be present in the cell at a sufficient concentration.
  2. The ComD receptor is activated by the AIP. This triggeres phosphorylation of the ComD receptor. The Phosphate group of the ComD receptor then binds to ComE. The phosphorylated ComE binds to the DNA and acts as a transcription factor.

Major assumptions:

  1. ComD and ComE are present in the cell/cell wall at a high concentration. ComD and ComE are both in steady-state, so the production and degradation constants are negligible.
  2. AIP and Phosphate are present inside/outside the cell at a high concentration. The degradation rates for these two species are negligible.
  3. Phosphorylation of the ComD receptor is modelled as an enzymatic reaction, neglecting the formation of an intermediate complex.
Fast Response Model

Goals:

    This model was mainly developed in order to determine whether simple production is better than 1-, 2- or 3-step amplification.

    Furthermore, an estimation of the speed of the response was desirable.

Elements of the system:

  1. Dioxygenase (blue on the diagrams below) is an enzyme that acts on catechol to produce a yellow output. In most of our models, dioxygenase was treated as an output because it was found that active dioxygenase acting on catechol produces the coloured output within a split second.
  2. GFP-Dioxygenase fusion protein (GFP is shown green on the diagrams). Dioxygenase, which is joined to GFP by the linker, was assumed to be inactive.
  3. TEV protease (pink on the diagrams below) has the ability to cleave the GFP-Dioxygenase fusion protein, hence, it activates dioxygenase
  4. Split TEV protease (purple on the diagrams below) is an inactive split form of TEV mounted on coiled coils. It can be activated again by coiled coils being cleaved by another active TEV.
Simple Production.png
Simple production upon activation of an arbitrary colour output by transcription and translation indicated by the blue arrow.
1-step amplification.png
Dioxygenase (C230) is simply produced. Upon activation at time t=0, it acts on catechol (cat.) to produce the yellow output - muconic acid. Catechol is not shown to be produced by the cell as it is added manually at an arbitrary time.
2-step amplification.png
The species that are shown in front of the vertical line (which indicates the beginning of the experiment) indicate that they have been accumulated beforehand in the cell. TEV protease activates inactive dioxygenase which acts on catechol to produce colour.
3-step amplification.png
This diagram introduces inactive split TEV protease attached to a coiled-coil as the third amplification step. Both inactive compounds have active sites for TEV to activate them, which results in multiple possibilities of action.

Major assumptions:

  1. The chemical and enzymatic reactions are modelled according to the Law of Mass Action.
  2. Our model assumes that the modelled system is inert within the bacterial cell or that reactions with other species within the bacterium is negligible. For example, the TEV protease is not supposed to cleave other molecules due to its specifity.