Team:Tec-Monterrey/ITESM14 modeling.html

From 2014.igem.org

Revision as of 06:27, 15 October 2014 by Ejzardainc (Talk | contribs)

<body> <section id="column">

Quorum Sensing: Modeling

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%">
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></td> Concentration of active LuxR compound by a LuxR homodimer \newline bound to 2 molecules of AHL</td>

</tr>

R</td> Concentration of unbound LuxR </td>

</tr>

p</td> Population parameter</td>

</tr>

d</td> Diffusion parameter for the autoinductor outside the cell </td>

</tr>

Bin</td> Bacteriophage concentration inside the E. Coli </td>

</tr>

Bext</td> Bacteriophage concentration outside the E. Coli </td>

</tr>

l</td> Diffusion parameter for the bacteriophage outside the cell</td>

</tr>

k0</td> Rate constant for the AHL generation due to the LuxI concentration </td>

</tr>

k1</td> Rate constant for the LuxR activation reaction </td>

</tr>

k2</td> Degradation/dissolution parameter for the LuxR-AHL</td>

</tr>

b</td> Degradation/dissolution parameter for the LuxI </td>

</tr>

r</td> Rate of AHL diffusion across the membrane for each cell </td>

</tr>

s</td> Rate of bacteriophague diffusion across the membrane for each cell </td>

</tr>

rb</td> Binding rate for the bacteriophages that binds to the cancer cells </td>

</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:

...

<script type="text/javascript"> katex.render("2A+2R \\leftarrow\\rightarrow R^*", eq1); </script>

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

...

<script type="text/javascript"> katex.render("\\frac{dR^*}{dt}=k_1\\left(A(t)\\right)^2\\left(R_T-2R^*(t)\\right)^2-k_2R^*(t)", eq2); </script>

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:

...

<script type="text/javascript"> katex.render("\\frac{dI}{dt}= a_0+\\frac{aR^*(t)}{K_M+R^*(t)}-bI(t)", eq3); </script>

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:

...

<script type="text/javascript"> katex.render("\\frac{dA_{ext}}{dt}= pr\\left(A(t)-A_{ext}(t)\\right)-dA_{ext}(t)", eq4); </script>

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.

...

<script type="text/javascript"> katex.render("\\frac{dB_{in}}{dt}= \\frac{bR^*(t)}{K_M+R^*(t)}-sB_{in}(t)", eq5); </script>

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:

...

<script type="text/javascript"> katex.render("\\frac{dB_{ext}}{dt}= (sp-l-r_b)B_{ext}(t)", eq6); </script>

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):

...

<script type="text/javascript"> katex.render("p_i=\\binom{N}{i}p^i(1-p)^{N-i}", eq7); </script>

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.

...

<script type="text/javascript"> katex.render("\\frac{dCells_{treated_{min}}}{dt}= \\frac{r_bB_{ext}(t)}{E_{max}}", eq8); </script>

...

<script type="text/javascript"> katex.render("\\frac{dCells_{treated_{max}}}{dt}= \\frac{r_bB_{ext}(t)}{E_{min}}", eq9); </script>


</section>

</body>

</div>