Team:TUDelft/Modeling

From 2009.igem.org

Revision as of 10:32, 24 July 2009 by Orr87 (Talk | contribs)


Delay Device

Negative Feedforward.jpg
Figure 1: Negative Feedforward model
Riboregulator delay.jpg
Figure 2: Riboregulator model

Mathematical Modeling

The Negative Feedforward model shown above can be expressed with the following set of ODEs (ordinary differential equations):
Negative feedforward ODEs 2.jpg

The Riboregulator model shown above can be expressed with the following set of ODEs:
Riboregulator ODEs.jpg

The symbols in this system of equations are found in the table below:
Symbols table.jpg

We would like to give credit to the Bologna '08 team and UNIPV-Pavia '08 team, who's project was useful in determining our own equations.

Conjugation

Using the paper "[http://www3.interscience.wiley.com/journal/118837812/abstract?CRETRY=1&SRETRY=0 A model for bacterial conjugal gene transfer on solid surfaces]", we established the assumptions and parameters that would be needed for the modeling of the conjugation.

Assumptions

  1. Cells are distributed randomly on the agar when introduced to it
  2. Number of cells grows until one or more of the medium components are exhausted, unless the initial nutrient concentration is lager than the saturation constant for growth (i.e. the ability of the colony to expand)
  3. Donors, recipients andf transconjugants have identical colony growth rates, specific growth rates and cell yields on solidified LB
  4. Conjugation occurs such that all recipients become transconjugants after a certain conjugation time
  5. Plasmid loss is negligable
  6. Cells only move on the surface through colony expansion

Parameters

  • Surface area of media (A)
  • Initial colony radius (r_0)
  • Specific growth rate (g_n)
  • Colony radial growth rate (g_r)
  • Maximum numbers of cells sustained by system (N_max)
  • Initial number of donors (N_d) and recipients (N_r)

All these parameters can be determined experimentally, or by literary research into previous conjugation experiments.

Spatial distribution of bacteria

The bacteria are initially distributed at random on the medium, so a distribution function to find the closest cell r would be:
F(r)=1-e^(-λ*r^2), where λ is the intensity (the expected number of bacteria per unit area).

Growth of colonies

With g_r being the exponential increasing rate of the colony radius at a certain time period, we can express the radius of the colony at time t as rad(t)=r_0*exp(g_r*t), where r_0 is the initial radius.
With g_n being the exponential increasing rate of the cell number within a colony at a certain time period, we can express the number of bacteria in a colony at time t as N(t)=exp(g_n*t).
Finally, there is a pre-given maximum yield of cells N_max, and growth stops when the total number of cells reaches this number.

Contact between colonies

Conjugation modeling

Final cell number