Team:Tec-Monterrey/ITESM14 modeling.html

From 2014.igem.org

<body>

Quorum Sensing: Modeling

<section id="column">

We built a mathematical model for the quorum sensing regulations on the basis of the Hill equation and made some slight modifications of the simple model exposed by Brian Ingalls in Mathematical Modelling in Systems Biology: An Introduction, when referring to the work of Sally James about the mathematical modelling the quorum sensing regulations of the V. fischeri (Ingalls, 2012).

First of all we define the different parameters we will be using:

<col style="width:20%"> <col style="width:80%">

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

</tr>

         </table>

</center>


We start our analysis with the active LuxR. Notice that the amount of unbound LuxR in a given moment is given by RT-2R*. this is because each LuxR-AHL complex requires an homodimer to be generated. The reaction in which this complex is involved, can be written as:


<figure> <a href="ITESM14_Ecuacion1.png" data-lightbox="Equation 1" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="ITESM14_Ecuacion1.png"></a> </figure>

Using simple chemical kinetics and adding a degradation/dissolution term we get the following differential equation:


<figure> <a href="ITESM14_Ecuacion2.png" data-lightbox="Equation 2" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="ITESM14_Ecuacion2.png"></a> </figure>

For the LuxI we use the Hill equation treated by Brian Ingalls (Ingalls, 2012), with a Hill coefficient equal 1, a basal expression level of a0 and the respective degradation/dissolution term:


<figure> <a href="ITESM14_Ecuacion3.png" data-lightbox="Equation 3" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="ITESM14_Ecuacion3.png"></a> </figure>

The interactions for the AHL are more complicate than the others. It is involved in the active LuxR's reaction. It is related with the LuxI concentration, adding a linear term k0I(t) from the LuxI concentration. Moreover, we must consider that the degradation of a LuxR-AHL molecule implies the creation of 2 AHL molecules, adding a term 2k2R*(t).Finally, there's also an interchange in the membrane, between the AHL inside and outside the bacteria. How fast this interchange occurs depends strongly on the membrane structure and the concentrations inside and outside the cell, which tends to created an equilibrium. Together, these interactions give us the following differential equation for the AHL inside the cell:


<figure> <a href="ITESM14_Ecuacion4.png" data-lightbox="Equation 4" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="ITESM14_Ecuacion4.png"></a> </figure>

By the same principle we develop the differential equation for the AHL outside the cell, considering that we should multiply the membrane-interchange term by the size of the population, assuming that the autoinducer production by each cell is uniform for all the bacterias. We also add a diffusion term for the autoinducer that dissipates away from the population:


<figure> <a href=" ITESM14_Ecuacion5.png" data-lightbox="Equation 5" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion5.png"></a> </figure>

For the bacteriophage production, we will consider only the production of the slowest part of the operon corresponding to the phague's parts, this will be also described with another Hill equation. We must include a diffusion parameter that will describe the amount of bacteriophagues that cross the membrane to the outside and will take into account the time required to ensemble the bacteriophague after the production of its components.


<figure> <a href=" ITESM14_Ecuacion6.png" data-lightbox="Equation 6" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion6.png"></a> </figure>

Finally we write an expression for the amount of phague outside the cells. Here we assume again an uniform production for all the bacterias and add also a diffusion term for the amount of phague that dissipates away from the population and a binding parameter rb for the amount of phague that binds to the cancer cells, this last parameter depends on the population of cancer cells:


<figure> <a href=" ITESM14_Ecuacion7.png" data-lightbox="Equation 7" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion7.png"></a> </figure>

Here rbBext(t) is the binding rate of the bacteriophage. If we want to predict the amount of cancer cells that are bound to a bacteriophague, the treated cells, we have to consider an stochastic model, taking into account that the amount of phagues bound to the cell will follows a binomial distribution (Smith H. and Trevino R., 2009):


<figure> <a href=" ITESM14_Ecuacion8.png" data-lightbox="Equation 8" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion8.png"></a> </figure>

Where pi is the probability for having i phagues bound to a given cell, p is the probability for a binding site to be occupied and N is the number of binding sites in the cell. From this we can have an expected value of Np and a standard deviation of √Np(1-p). More over we can estimate a maximum Emax and a minimum Emin for the number of phages bound per cancer cell. Taking this consideration together with the previous differential equations, we can write an expression for rate of the minimum and maximum possible values of treated cell.


<figure> <a href=" ITESM14_Ecuacion9.png" data-lightbox="Equation 9" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion9.png"></a> </figure>
<figure> <a href=" ITESM14_Ecuacion10.png" data-lightbox="Equation 10" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src=" ITESM14_Ecuacion10.png"></a> </figure>


</section>

</body>

</div>

A Average concentration of the AHL inside the E. Coli
Aext Average concentration of the AHL outside the E. Coli
KM The activation coefficient for the corresponding promoter
RT Total concentration of LuxR monomers
R*</sub></th> Concentration of active LuxR compound by a LuxR homodimer \newline bound to 2 molecules of AHL
R</th> Concentration of unbound LuxR
p</th> Population parameter
d</th> Diffusion parameter for the autoinductor outside the cell
Bin</th> Bacteriophage concentration inside the E. Coli
Bext</th> Bacteriophage concentration outside the E. Coli
l</th> Diffusion parameter for the bacteriophage outside the cell
k0</th> Rate constant for the AHL generation due to the LuxI concentration
k1</th> Rate constant for the LuxR activation reaction
k2</th> Degradation/dissolution parameter for the LuxR-AHL
b</th> Degradation/dissolution parameter for the LuxI
r</th> Rate of AHL diffusion across the membrane for each cell
s</th> Rate of bacteriophague diffusion across the membrane for each cell
rb</th> Binding rate for the bacteriophages that binds to the cancer cells