rb</th>
| Binding rate for the bacteriophages that binds to the cancer cells |
</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="
" data-lightbox="Equation 1" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></a>
</figure>
Using simple chemical kinetics and adding a degradation/dissolution term we get the following differential equation:
<figure>
<a href="
" data-lightbox="Equation 2" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 3" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 4" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 5" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 6" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 7" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 8" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></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="
" data-lightbox="Equation 9" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></a>
</figure>
<figure>
<a href="
" data-lightbox="Equation 10" data-title=""><img class="img img-responsive" style="margin:0px auto;display:block width:20%;" src="
"></a>
</figure>
</section>
</body>
</div>