Team:British Columbia/Modeling

From 2013.igem.org

(Difference between revisions)
Line 38: Line 38:
We use the Poisson distribution to statistically determine the expected infected bacterial population.
We use the Poisson distribution to statistically determine the expected infected bacterial population.
\begin{align}
\begin{align}
-
P(n) = \frac{m^n e^{-m}}{n!}
+
P(n) = \frac{MOI^n e^{-MOI}}{n!}
\end{align}
\end{align}
-
If one or more than one virus collides with a cell, the cell will be infected. Thus, we calculate the fraction of cells to not be infected.
+
Where $n$ is the number of bacteriophages attacking a cell. If one or more than one virus collides with a cell, the cell will be infected. Therefore, we calculate the fraction of cells to not be infected.
\begin{align}
\begin{align}
 +
P(n > 0) = 1 - P(0) = 1 - e^{-MOI}
\end{align}
\end{align}

Revision as of 22:00, 31 August 2013

iGEM Home

Reaction Kinetics: - Modelling the production of cinnamaldehyde, vanillin and caffeine production


Population dynamics: - Modelling the population of a bacteria in a batch reactor with phage dependency

\begin{align} \frac{dX}{dt} = \frac{dX_u}{dt} + \frac{dX_i}{dt} + \delta \end{align} Where $\frac{dX}{dt}$ is rate of change in bacteria population over time, $\frac{dX_u}{dt}$ and $\frac{dX_i}{dt}$ are the rates of change in the uninfected and infected populations respectively. $\Gamma$ is defined as a stepwise function. \begin{align} \delta &= - E(X_i), \ \ \ during \ phage \ caused \ lysis\\ &= 0, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ otherwise \end{align} Where $E(X_i)$ is the statistically expected infected bacterial populations. We define the rate of change in population in () and (). \begin{align} &\frac{dX_u}{dt} = \big(\mu_u + k_{d_u}\big)X_u\\ &\frac{dX_i}{dt} = \big(\mu_i + k_{d_i}\big)X_i \end{align} Where $\mu_u$ and $\mu_i$ are the growth rate of uninfected and infected cells respectively. And $k_{d_u}$, $k_{d_i}$ are the decay rates of the uninfected and infected bacteria. To create a more complete model, we have chosen to model the bacteriophage population as well. \begin{align} P = X_i \beta ^{\frac{\tau}{LT}} \end{align} However, due to the nature of the phages in our experiment (in that they have a very low or non-existent phage secretion rate), we can redefine our phage population as a stepwise function. \begin{align} \Delta P &= E(X_i)\beta, \ \ \ \ for \ lysis\\ &= 0, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ otherwise \end{align}

The multiplicity of infection (MOI) is defined as the ratio of bacteriophages to bacterial cells. \begin{align} MOI = \frac{P}{X} \end{align} We use the Poisson distribution to statistically determine the expected infected bacterial population. \begin{align} P(n) = \frac{MOI^n e^{-MOI}}{n!} \end{align} Where $n$ is the number of bacteriophages attacking a cell. If one or more than one virus collides with a cell, the cell will be infected. Therefore, we calculate the fraction of cells to not be infected. \begin{align} P(n > 0) = 1 - P(0) = 1 - e^{-MOI} \end{align}