Team:Colombia/Modeling/Diff

From 2012.igem.org

(Difference between revisions)
(Differential equations for Pestsbuster)
(RALSTONIA)
 
(48 intermediate revisions not shown)
Line 9: Line 9:
<p align="justify">
<p align="justify">
-
Before experimentation with the designed bacteria and plants, we can check if this design works through a mathematical model of it. This model is based on the use of differential equations. This type of equations are deterministic because it describes the behavior for each of the substances in the synthetic circuit. However, they do not take into account the probabilities involved in each of the events described.   
+
Before starting to work at the lab with the designed bacteria and plants, we can check if this design works through a mathematical model. The first step is to develop a deterministic model based in differential equations. This type of equations are deterministic because they describe the behavior for each of the substances in the synthetic circuit over time. However, they do not take into account the probabilities involved in each of the events described, the population interactions and the noise of the system.   
-
Based on mass conservation:  
+
All of these equations are based on the law of mass conservation:  
-
[[File:Mass balance.jpg|thumb|center|500x300px]]
+
[[File:Mass balance1.png|center|530px]]
-
 
+
-
 
+
-
In a biological system, the accumulation stands for the change in time, the input and output are related with the processes of exportation and importation of a molecule into the cell, the generation is related to the production by a gene or a chemical reaction, and the consumption is related to the chemical reaction too. All these terms depend on the reaction kinetics that could be a simple multiplication or a complex expression like a Hill equation.
+
 +
In a biological system, the accumulation stands for the concentration changes over time, the input and output are related with the processes of export and import of a molecule into the cell, the generation is related to the production by a gene or a chemical reaction, and the consumption is related to chemicals reactions too. All these terms depend on reaction kinetics that could be expressed as a simple multiplication or a complex expression like a Hill equation.
Line 24: Line 22:
'''Example:'''
'''Example:'''
-
[[File:difeq.png|thumb|center|500x300px|Example process]]
+
[[File:difeq.png|thumb|center|530px| Figure 1. Biological system interacting]]
-
The molecule W promotes the production of A, which is used in a reaction with B for producing C and for the promotion of D.   
+
The molecule W activates de promoter for the production of A, which is used in a reaction with B for producing C and for the activation of the D's promoter.   
-
Based on the mass balance for substance A, it is possible to see that there is not exportation or importation. Besides, there is not input or output term. The generation term is the production by the gene promoted by W, and it can be expressed with a Hill kinetic:
+
To describe the changes over time of the protein A it is possible to see that there is not export or import; as a consequence, there is not input or output term. The generation term is the production by the gene promoted by W, and it can be expressed with a Hill type equation:
-
[[File:hill.png|center|100x100px]]
+
[[File:hil1l.png|center|80x80px]]
-
The consumption term is divided in two: first the reaction with B, that can be expressed as a first order reaction with a kinetic constant m and the promotion of the protein D that has a Hill kinetic:
+
The consumption term is divided in two: first the reaction with B, that can be expressed as a first order reaction with a kinetic constant m and the activation of the promoter of D; however this activation is an event that occurs in a short time scale so it can be removed
-
[[File:consumption.png|center|200x100px]]
+
[[File:consumption12.png|center|80x80px]]
Writing the complete mass balance equation for substance A:
Writing the complete mass balance equation for substance A:
-
[[File:eqq.png|center|400x100px]]
+
[[File:eqq1.png|center|280x100px]]
-
This procedure was made for each of substances in the designed circuit. As shown above, Hill equations were used. Once all parameters in Hill equation are known, these can be solved using a specialized software like Matlab. This software allows us to solve the differential equations with a Runge - Kutta algorithm or using one of its tools like ODE45.  
+
This procedure was made for each of substances in the designed circuit. As shown above, Hill type equations were used to describe the production by the activation of a promoter and the reactions terms followed the law of mass action. Once all the equations are describe and the constants are set we can solve the differential equation with specialized software like MATLAB and its tools such as ODE.
   
   
-
== '''Differential equations for Pestsbuster''' ==
+
== '''Differential equations for Pest-busters''' ==
-
==='''''RUST'''''===
+
==='''''FUNGAL INFECTION'''''===
-
Below it is shown how differential equations for Pest-busters were made. However, before describing the equations it is necessary to know all the symbols that will be use in the documents and the simulation. The following table contains all the substances involved in the processes and the constants required for the simulation.
+
Below it is shown how differential equations for Pest-busters were made. However, before describing the modeled system ,it is necessary to know all the symbols that will be use in the documents and the simulation. The following table contains all the substances involved in the processes and the constants required for the simulation.
-
[[File:tablerust.png|center|500x1000px]]
+
<html>
 +
<br>
 +
</br>
 +
</html>
-
'''''CBP:'''''  The Chitin binding protein participates in two major processes. First, it binds to the sensor. Second, when the chitin is in the cell, the CBP binds to it.
+
[[File:tablerust.png|thumb|center|500x1000px|Table 1. Symbols used in the document and simulation ]]
 +
 
 +
<html>
 +
<br>
 +
</br>
 +
</html>
 +
 
 +
 
 +
'''''CBP:'''''  The Chitin binding protein participates in two major processes. First, it binds to the sensor. then, when the chitin is in the cell, the CBP binds to it.
[[File:reactionCBP.png|center|200x150px]]
[[File:reactionCBP.png|center|200x150px]]
Line 64: Line 73:
[[File:eq12.png|center|400x280px]]
[[File:eq12.png|center|400x280px]]
-
'''''ChiP:''''' The chitioporin is the one that allows the entrance of the chitin monomers into the cell. Its production is regulated by the sensor (S), this regulation is expressed as a hill equation. The differential expression also includes the basal production and the rate of destruction; it is assumed that all the chitoporin produced goes to the membrane.
+
'''''ChiP:''''' The chitioporin is the protein that allows the entrance of the chitin monomers into the cell. Its production is regulated by the sensor (S), this regulation is expressed as a hill type equation. The differential expression also includes the basal production and the rate of destruction; it is assumed that all the chitoporin produced goes to the membrane.
[[File:eq3.png|center|280x170px]]
[[File:eq3.png|center|280x170px]]
-
'''''ChiA : chitinase'''''
+
'''''ChiA : The Chitinase'''''
-
The chitinase inside the cell depends of its basal production and the induction made by the sensor protein (S). The exportation of this enzyme is also involved in the equation.  
+
inside the cell depends of its basal production and the induction made by the sensor protein (S) as a positive feedback. The exportation of this enzyme is also involved in the equation.  
[[File:eq4.png|center|380x200px]]
[[File:eq4.png|center|380x200px]]
-
The chitinase outside the cell depends of the rate of exportation of the chitinase inside, the diffusion in the environment and the reaction with the dimers of chitin.  
+
The chitinase outside the cell depends of the export rate of the chitinase inside, the diffusion in the environment and the reaction with the dimers of chitin.  
[[File:eq5.png|center|400x280px]]
[[File:eq5.png|center|400x280px]]
-
'''''Chitin:''''' This equation measures the monomers of chitin inside the cell.  This depends of a function of importation that depends on the quantity of chitoporins and the creation of chitin monomers outside the cell. This monomer is consumed in the reaction made with the complex C.S.  
+
'''''Chitin:''''' This equation measures the chitin oligomers inside the cell.  These depends of a function of importation that depends on the quantity of chitoporins and the creation of chitin monomers outside the cell. This oligomers are consumed in the reaction made with the complex C.S.  
[[File:eq6.png|center|400x200px]]
[[File:eq6.png|center|400x200px]]
-
'''''Sensor:''''' This protein is one of the most important in the detection system, its function is to activate the promoters of P,A and I. Once the chitin enters to the cell S frees from the complex C.S and starts activating promoters.
+
'''''Sensor:''''' This protein is one of the most important in the detection system, its function is to activate the promoters of P,A and I. Once the chitin enters to the cell S is is freed from the complex C.S and starts activating promoters.
-
This protein has two equations, the first one is and expression that shows the total of S in the cell. And the last one describes de total of the S protein free in the cell.
+
This protein has two equations, the first one is an expression that shows the total of S in the cell, and suppose that it remains constant. The second one shows the terms involved in the formation of the complex C.S and the liberation of S after the chitin enters into the cell. It also includes terms associated to binding of the S protein with promoters regions and its reverse reaction, the unbinding of the promoter site. This expression was not taken into account, because the activation of the promoters are insignificant in the scale of time we are using.   
-
 
+
-
The second one shows the terms involved in the formation of the complex C.S and the liberation of S after the chitin enters into the cell. It also includes terms associated to the binding of the S protein with promoters regions and its reverse reaction, the unbinding of the promoter site. This expression wasn't taken into account, because its change it insignificant compared with the all system.   
+
[[File:eq7.png|center|210x200px]]
[[File:eq7.png|center|210x200px]]
-
'''''LuxI:''''' The I protein that is inside the cell depends on the basal production and normal degradation rate, also of its promotion by the S protein, the exportation and importation from the cell and its reaction to make the complex with R
+
'''''LuxI:''''' The I protein inside the cell depends on the basal production and normal degradation rate.  Also, its promotion by the S protein, export and import from the cell and the reaction with R to make the I.R complex.
[[File:eq9.png|center|530x280px]]
[[File:eq9.png|center|530x280px]]
-
The I protein outside the cell only depends on the rates of exportation and importation from the cells and the diffusion into the environment.  
+
The I protein outside the cell only depends on the rates of exportation and importation of the cells and the diffusion in the environment.  
[[File:eq10.png|center|400x280px]]
[[File:eq10.png|center|400x280px]]
-
The complex I.R activates one of the promoters of CI and HB and its creation depends of the concentration of I and R .
+
The complex I.R creation depends of the concentration of I and R. It activates the production of HipB, CI and Salicylic acid.  
-
[[File:eq11.png|center|400x280px]]
+
[[File:eq11b.png|center|300x280px]]
-
'''''LuxR:''''' The R protein is only involved in the creation of the complex L.R.
+
'''''LuxR:''''' The R protein is only involved in the creation of the complex L.R. It is downstream the I protein so its production also depends on the activation of the promoter by S
[[File:eq12r.png|center|350x280px]]  
[[File:eq12r.png|center|350x280px]]  
-
'''''CI and HipB: ''''' The same promoters produce these two proteins downstream. Its creation depends on the activation of the promoters, the basal production and rate of destruction.  Even though the CI promoter has three boxes it was found in the literature a paper where it can be modeled as a simple box and the parameter are already known .
+
'''''CI and HipB: ''''' The same promoters produce these two proteins downstream. Its creation depends on the activation of the promoters, the basal production and rate of destruction.  Although, CI promoter has three boxes, it was found in the literature that can be modeled as a simple box and the experimental parameters are known. HipB forms dimers and interact with HipA7 to "awake" the cell
[[File:eq1314.png|center|400x280px]]  
[[File:eq1314.png|center|400x280px]]  
Line 111: Line 118:
[[File:eq14.png|center|550x280px]]  
[[File:eq14.png|center|550x280px]]  
-
'''''HipA7 and Salicylic acid:''''' These two proteins are activated by CI. Its concentration on the cell depend of the basal production, the rate of production when it is activated and the rate of destruction. Additionally the Salicylic acid is exported outside the cell for the plant signaling, and the toxin HA reacts with the antitoxin to awake the cell.  
+
'''''HipA7 ''''' This protein is activated by CI. Its concentration on the cell depends of the basal production, the rate of production when it is activated and the rate of destruction. The toxin HA reacts with the antitoxin to awake the cell.  
[[File:eq1516.png|center|430x280px]]
[[File:eq1516.png|center|430x280px]]
-
[[File:eq16.png|center|380x280px]]
+
''' Salicylic Acid: The desired response'''
 +
 
 +
The salicylic acid is produced when the IR and the CI promoters are activated. It also is exported outside the cell because it works as a phytohormone that alerts the plant.
 +
 
 +
[[File:eq16b.png|center|450x280px]]
Line 122: Line 133:
==='''''RALSTONIA'''''===
==='''''RALSTONIA'''''===
-
In this model the only thing that changes is the detection system. Now we want to sense 3-OH-PAME instead of Chitin. All the system involving CBP, chitinase and chitoporin is replaced with the sensor pchS and the transcription activator pchA, so the differences in the differential equations are in the first part only. The following table presents all the new substances involved in the process:
+
For this mathematical model the only thing that changes is the detection system. Now we want to sense 3-OH-PAME instead of Chitin. All the system involving CBP, chitinase and chitoporin is replaced with the sensor pchS and the transcription activator pchA. Thus, differences in the differential equations are in the first part only. The following table presents all the new substances involved in the process:
-
[[File:tableralstonia.png|center|380x350px]]
+
[[File:tableralstonia.png|thumb|center|380x350px|Table 2. Substances involved in Ralstonia detection]]
-
'''''phcS:'''''  This molecule is the sensor of the OH. When the OH is in the environment the sensor is phosphorylated and phosphorylates the phcR-phcA complex. This set of equations includes the basal production and normal degradation of the sensor and the phosphorylation processes.
+
'''''phcS:'''''  This molecule is the sensor of the OH. When the OH is in the environment, sensor is phosphorylated and phosphorylates the phcR-phcA complex. This set of equations includes the basal production and normal degradation of the sensor and the phosphorylation processes.
[[File:eq1ral.png|center|280x180px]]
[[File:eq1ral.png|center|280x180px]]
-
'''''phcR:''''' This complex is normally in the cell with a rate of basal production and destruction. It is broken when the Sf phosphorylates it and the phcA is freed.
+
'''''phcR:''''' This complex is normally in the cell with a rate of basal production and destruction. It is broken when the Sf phosphorylates it and the phcA is released.
[[File:eq2ral.png|center|230x200px]]
[[File:eq2ral.png|center|230x200px]]
-
'''''phcA:'''''  This protein is the activator, it is freed after the phosphorylation of R and it binds to a promoter, in this case the CI and HipB promoter.
+
'''''phcA:'''''  This protein is the activator and its released after the phosphorylation of R.  It binds to a promoter, in this case the CI and HipB promoter.
-
[[File:eq3ral.png|center|320x280px]]
+
[[File:eq3ralb.png|center|230x280px]]
-
The equations of the second part are the same as in the rust model (Equations 9-16), the only difference is that the promoter of the LuxI-LuxR system is now induced by A:
+
The equations of the second part are the same as in the rust model (Equations 9-16). The difference is that promoter of the LuxI-LuxR system is now induced by A:
[[File:eqchangral.png|center|510x280px]]
[[File:eqchangral.png|center|510x280px]]

Latest revision as of 01:57, 27 October 2012

Contents

Team Colombia @ 2012 iGEM

Template:Https://2012.igem.org/User:Tabima

Differential equations

Before starting to work at the lab with the designed bacteria and plants, we can check if this design works through a mathematical model. The first step is to develop a deterministic model based in differential equations. This type of equations are deterministic because they describe the behavior for each of the substances in the synthetic circuit over time. However, they do not take into account the probabilities involved in each of the events described, the population interactions and the noise of the system. All of these equations are based on the law of mass conservation:

Mass balance1.png

In a biological system, the accumulation stands for the concentration changes over time, the input and output are related with the processes of export and import of a molecule into the cell, the generation is related to the production by a gene or a chemical reaction, and the consumption is related to chemicals reactions too. All these terms depend on reaction kinetics that could be expressed as a simple multiplication or a complex expression like a Hill equation.


Below it is shown an example to understand the previous ideas:

Example:

Figure 1. Biological system interacting


The molecule W activates de promoter for the production of A, which is used in a reaction with B for producing C and for the activation of the D's promoter.

To describe the changes over time of the protein A it is possible to see that there is not export or import; as a consequence, there is not input or output term. The generation term is the production by the gene promoted by W, and it can be expressed with a Hill type equation:


Hil1l.png

The consumption term is divided in two: first the reaction with B, that can be expressed as a first order reaction with a kinetic constant m and the activation of the promoter of D; however this activation is an event that occurs in a short time scale so it can be removed


Consumption12.png


Writing the complete mass balance equation for substance A:

Eqq1.png


This procedure was made for each of substances in the designed circuit. As shown above, Hill type equations were used to describe the production by the activation of a promoter and the reactions terms followed the law of mass action. Once all the equations are describe and the constants are set we can solve the differential equation with specialized software like MATLAB and its tools such as ODE.


Differential equations for Pest-busters

FUNGAL INFECTION

Below it is shown how differential equations for Pest-busters were made. However, before describing the modeled system ,it is necessary to know all the symbols that will be use in the documents and the simulation. The following table contains all the substances involved in the processes and the constants required for the simulation.



Table 1. Symbols used in the document and simulation




CBP: The Chitin binding protein participates in two major processes. First, it binds to the sensor. then, when the chitin is in the cell, the CBP binds to it.

ReactionCBP.png

This set of equations involves de change of CBS free in the cell and the change of the complex C.S. With the production, the destruction rate and the reactions describe above the equations are:

Eq12.png

ChiP: The chitioporin is the protein that allows the entrance of the chitin monomers into the cell. Its production is regulated by the sensor (S), this regulation is expressed as a hill type equation. The differential expression also includes the basal production and the rate of destruction; it is assumed that all the chitoporin produced goes to the membrane.

Eq3.png

ChiA : The Chitinase inside the cell depends of its basal production and the induction made by the sensor protein (S) as a positive feedback. The exportation of this enzyme is also involved in the equation.

Eq4.png

The chitinase outside the cell depends of the export rate of the chitinase inside, the diffusion in the environment and the reaction with the dimers of chitin.

Eq5.png

Chitin: This equation measures the chitin oligomers inside the cell. These depends of a function of importation that depends on the quantity of chitoporins and the creation of chitin monomers outside the cell. This oligomers are consumed in the reaction made with the complex C.S.

Eq6.png

Sensor: This protein is one of the most important in the detection system, its function is to activate the promoters of P,A and I. Once the chitin enters to the cell S is is freed from the complex C.S and starts activating promoters.

This protein has two equations, the first one is an expression that shows the total of S in the cell, and suppose that it remains constant. The second one shows the terms involved in the formation of the complex C.S and the liberation of S after the chitin enters into the cell. It also includes terms associated to binding of the S protein with promoters regions and its reverse reaction, the unbinding of the promoter site. This expression was not taken into account, because the activation of the promoters are insignificant in the scale of time we are using.

Eq7.png

LuxI: The I protein inside the cell depends on the basal production and normal degradation rate. Also, its promotion by the S protein, export and import from the cell and the reaction with R to make the I.R complex.

Eq9.png

The I protein outside the cell only depends on the rates of exportation and importation of the cells and the diffusion in the environment.

Eq10.png

The complex I.R creation depends of the concentration of I and R. It activates the production of HipB, CI and Salicylic acid.

Eq11b.png

LuxR: The R protein is only involved in the creation of the complex L.R. It is downstream the I protein so its production also depends on the activation of the promoter by S

Eq12r.png

CI and HipB: The same promoters produce these two proteins downstream. Its creation depends on the activation of the promoters, the basal production and rate of destruction. Although, CI promoter has three boxes, it was found in the literature that can be modeled as a simple box and the experimental parameters are known. HipB forms dimers and interact with HipA7 to "awake" the cell

Eq1314.png
Eq14.png

HipA7 This protein is activated by CI. Its concentration on the cell depends of the basal production, the rate of production when it is activated and the rate of destruction. The toxin HA reacts with the antitoxin to awake the cell.

Eq1516.png

Salicylic Acid: The desired response

The salicylic acid is produced when the IR and the CI promoters are activated. It also is exported outside the cell because it works as a phytohormone that alerts the plant.

Eq16b.png



RALSTONIA

For this mathematical model the only thing that changes is the detection system. Now we want to sense 3-OH-PAME instead of Chitin. All the system involving CBP, chitinase and chitoporin is replaced with the sensor pchS and the transcription activator pchA. Thus, differences in the differential equations are in the first part only. The following table presents all the new substances involved in the process:

Table 2. Substances involved in Ralstonia detection


phcS: This molecule is the sensor of the OH. When the OH is in the environment, sensor is phosphorylated and phosphorylates the phcR-phcA complex. This set of equations includes the basal production and normal degradation of the sensor and the phosphorylation processes.

Eq1ral.png

phcR: This complex is normally in the cell with a rate of basal production and destruction. It is broken when the Sf phosphorylates it and the phcA is released.

Eq2ral.png


phcA: This protein is the activator and its released after the phosphorylation of R. It binds to a promoter, in this case the CI and HipB promoter.

Eq3ralb.png


The equations of the second part are the same as in the rust model (Equations 9-16). The difference is that promoter of the LuxI-LuxR system is now induced by A:

Eqchangral.png
Eqchangra2l.png