Team:Colombia/Deterministic

From 2014.igem.org

(Difference between revisions)
Line 119: Line 119:
When CAI is present in the media TetR production turns off, which can be seen as a small change in TetR concentration; even tough the change is small it is enough to decrease the dimer concentration by 10 times. Now is time to see if this decrease is enough to turn on the production of the signal.  
When CAI is present in the media TetR production turns off, which can be seen as a small change in TetR concentration; even tough the change is small it is enough to decrease the dimer concentration by 10 times. Now is time to see if this decrease is enough to turn on the production of the signal.  
 +
The first model we try was given by the equation:
 +
 +
[[File: EqJ1.png|center|680px]]
[[File: Colombia_Test6.png|center|400px]]
[[File: Colombia_Test6.png|center|400px]]

Revision as of 01:36, 30 September 2014

Wheeltz - CSS3 Navigational Wheel Menu

  • Home
  • iGEM
  • Facebook
  • Twitter



Deterministic model

Finding the mean behaviour of the system...

The first step is to develop a deterministic model based in differential equations. With the ordinary differential equations we can describe quantitatively the system at molecular level.

Using the law of mass action is possible to model most of the metabolic networks. For each substance of the system, we express the changes of concentration trough time. We consider processes related with the export and import of a molecule into the cell, the production by a gene or a chemical reaction. All these terms depend on reaction kinetics that could be expressed as a simple multiplication or a complex expression like a Hill equation.

Signal transduction →

The receptor CqsS acts as a kinase when there is no auto inducer (CAI). It phosphorylates LuxU, which transfers its phosphate to LuxO. And in our case LuxO phosphorylated activates Pqrr promoter.

Once the inducer is in the media the receptor changes to a phosphase (CqsSp) mode and the flow is reversed. LuxU is unphosphorylated and it removes the phosphate from LuxO. In consequence Pqrr promoter is repressed.

Colombia Rxn1.png
Colombia Rxn2.png
Colombia Rxn3.png
Colombia Rxn4.png
Colombia Rxn5.png


Taking this process into account and given the fact the phosphorylation process can been described as an enzymatic process; we used the Mass Action Law and the Michaelis-Menten Kinetics in order to obtain the mathematical expressions that describe the signal transduction:

EqA.png


EqB.png

Equations a and b describe the behavior of the receptor acting as a kinase (CqsS) and as a phosphase (CqsSP). The kinase is produced at a constant rate alpha CS, it is involved in Reaction 1 and has a decay rate of Gamma CS. The phosphase is also involved in Reaction 1 and has its own decay rate.


EqC.png


EqD.png


EqE.png


EqF.png


To describe the phosphorylation cascade we based our expressions on the model that Liu et al. did in 2012: LuxU and LuxO are involved in the enzymatic Reactions 2-5, both are produced by a constitute promoter at a rate a_u and a_o respectively and each molecule has its own decay rate, taking this into account Equations c-f were made.


Response →

LuxO controls Pqrr promoter. Downstream of this promoter there is the repressor TetR of the tetracycline promoter (pTet), which in our design is downstream the repressor and, controls de production of the colored response and the activator of the promoter TetA.


EqG.png


EqH.png


LuxO controls the production of TetR and it can be expressed with a Hill type equation where LuxO phosphorylated is the activator. The promoter has a small constitutive production rate alpha tr and TetR has a decay rate of gamma tr; also TetR is consumed when forms a dimer in order to repress pTet.

EqI.png


EqJ.png


Equations i and j describe the production of the response and the pTet activator. As in all the preview equation, both proteins have a constitutive production and a decay rate. In our design there is a promoter (pTet) with a unique operator sequence (TetO) where the repressor (TetR2) and the activator (tTA) bind (TetR2 and tTA have the same binding domain). For that reason it was necessary to propose a mathematical expression in which the repressor and the activator compete for the binding site.


Testing our equations


Before we established the equations above, we tried different versions of the systems. In order to get the desired response for our S. coli we needed to find the right circuit. In this section we will be showing how we did it, step by step.

The first we did, was to establish the desired response. We want the system to react with a low and long lasting concentration of CAI. The Figure below shows the impulse we used to test the system: two pulses with different concentrations and different durations.

Colombia Test1.png

We knew that the equations of the CqsS protein behaviour are right because they were previously reported by Lelu et al (2012). So we needed to check that behaviour, the figure below shows the answer of CqsS with CAI.

Colombia Test2.png

As soon as the CAI molecule is present the CqsS concentration is instantly decreased, and CqsS* (which stands for the CAI bound with the Cqss protein) increases. After half an hour the concentration of CqsS starts to decrease because it is consume by LuxU,just the way it’s supposed to.

The next step is to check what happens with LuxU and its phosphorylated version now that we know that Cqss is modeled properly we want to check the phosphorylation/ phosphate cascade. As with CqsS this molecule was already modelled by Lelu et al (2012).

Colombia Test3.png

Just as expected whenever when CAI is present (50h-100h and 250h) LuxU is greater than LuxU* , due to the phosphase function of Cqss*. At the beginning we can see peak in LuxU concentration because of the Le Chatelier principle the reaction tends to go forward; after 3 hours the LuxU stabilises itself thanks to the phosphorylation reaction.

Now let’s check LuxO and its phosphorylated version. Even tough this part of the cascade was modelled by Lelu et al, we did no take into account some of the regulation routes they did, like de miRNA, because our system does not have them. We wanted to check the respond of the most important molecule from the cascade

Colombia Test4.png

The first time we modelled the system, LuxO response was a slower than we want. So we changed the hill constant coefficient, proposed by Lelu et al in order to make the response faster and last longer. This change compensates the lack of the negative feedback caused by miRNA that is present in V. cholera but not in S. coli.

With this modification we can see that the systems takes 20 hours to react and after 30h the production of LuxO is completely on and LuxO* is almost zero, which is desirable, because we do not want any LuxO phosphorylated to turn on the Pqrr4 promoter.

So far the model is working great, and all this part correspond to the parts registered by Peking 2011 Team. But we are really interested in the response signal and we want to see if S. coli would be able to shows us when V. cholerae is in the media. Given that when CAI is in the media de Pqrr4 promoter turns off, we need an NOT logical gate to invert the transduction, so we used the pTet and TetR system.

Colombia Test5.png

When CAI is present in the media TetR production turns off, which can be seen as a small change in TetR concentration; even tough the change is small it is enough to decrease the dimer concentration by 10 times. Now is time to see if this decrease is enough to turn on the production of the signal.

The first model we try was given by the equation:

EqJ1.png
Colombia Test6.png

This is a really disappointing graph, but there’s plenty of tweaking we can do to increase the response, since it is obviously not going to be enough to be recognized. Luckily pTet also has an activator molecule TetA so we added it to the model as if it was also downstream of the pTet and not only should it increase the production of the chomoprotein but also itself

Colombia Test7.png

Now that we can see that the activator is actually being produced lets check how the response goes with this feedback.

Colombia Test8.png

And now the response is up, now this might be enough or it might not, there’s still much tweaking we can do to increase the response signal but if this is enough we can use this. When the part is ready at the lab we will know.