Team:TUDelft/Modeling

From 2009.igem.org

(Difference between revisions)
(Contact between colonies)
 
(10 intermediate revisions not shown)
Line 1: Line 1:
-
{{Template:TUDelftiGEM2009}}
+
{{Template:TUDelftiGEM2009_menu_Model}}
-
<br>
+
-
='''Delay Device'''=
+
-
[[Image:Negative Feedforward.jpg]]<br>
+
-
Figure 1: Negative Feedforward model<br>
+
-
[[Image:Riboregulator_delay.jpg]]<br>
+
-
Figure 2: Riboregulator model<br>
+
-
==Mathematical Modeling==
+
-
The Negative Feedforward model shown above can be expressed with the following set of ODEs (ordinary differential equations):<br>
+
-
[[Image:negative_feedforward_ODEs_2.jpg]]<br><br>
+
-
The Riboregulator model shown above can be expressed with the following set of ODEs:<br>
+
-
[[Image:riboregulator_ODEs.jpg]]<br><br>
+
-
The symbols in this system of equations are found in the table below:<br>
+
-
[[Image:symbols_table.jpg]]<br><br>
+
-
We would like to give credit to the [https://2008.igem.org/Team:Bologna/Modeling Bologna '08 team] and [https://2008.igem.org/Team:UNIPV-Pavia/Modeling 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'''==
+
-
#Cells are distributed randomly on the agar when introduced to it
+
-
#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)
+
-
#Donors, recipients andf transconjugants have identical colony growth rates, specific growth rates and cell yields on solidified LB
+
-
#Conjugation occurs such that all recipients become transconjugants after a certain conjugation time
+
-
#Plasmid loss is negligable
+
-
#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:<br>
+
-
F(r)=1-exp(-λ*r^2), where λ is the intensity (the expected number of bacteria per unit area).
+
-
=='''Growth of colonies'''==
+
='''Modeling'''=
-
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.<br>
+
-
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).<br>
+
-
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'''==
+
Please select one of the two modeling sections: [[Team:TUDelft/Modeling_Cascade |Transcriptional Cascade Modeling]] or [[Team:TUDelft/Modeling_Conjugation |Conjugation Modeling]].
-
Two colonies with radius ''r'' will touch if their centres are less then ''2r'' apart (i.e. overlapping colonies). Therefore, we can express the probability that the colonies will not be interacting by the function<br>
+
-
p(Not touching)=p(radius1+radius2>2r)=p1=exp(-4*r^2*π*λ)<br>
+
-
When this is considered as a function of time, the function would be<br>
+
-
p1(t)=exp(-4*π*λ*(r_0*exp(g_r*t))^2)=exp(-4*π*λ*(r_0)^2*exp(2*g_r*t))<br
+
-
When we approximate the distribution of colony clumps, it depends on the parameter ''p1'', so the probability that a colony chosen at random from the plate is part of a clump containing ''n'' colonies is<br>
+
-
p_n=p1(1-p1)^(n-1)=exp(-4*r^2*π*λ)(1-exp(-4*r^2*π*λ))^(n-1)<br>
+
-
This function in terms of time would be expressed as<br>
+
-
p_n(t)=exp(-4*π*λ*(r_0)^2*exp(2*g_r*t))(1-exp(-4*π*λ*(r_0)^2*exp(2*g_r*t)))^(n-1)<br>
+
-
In this case, the intensity λ=(N_d+N_r)/A, where A is the area of the medium, and N_d and N_r are the number of donors and recipients at t=0, respectively.
+
-
=='''Conjugation modeling'''==
+
<html>
-
 
+
<table border=0>
-
=='''Final cell number'''==
+
<tr>
 +
<td>
 +
<a border=0 href=https://2009.igem.org/Team:TUDelft/Modeling_Cascade><img src=https://static.igem.org/mediawiki/2009/f/f9/TUDtranscriptmodelani.gif></a>
 +
</td>
 +
<td>
 +
<img src=https://static.igem.org/mediawiki/2009/f/f1/TUD1pxspacer.png width=40px>
 +
</td>
 +
<td>
 +
<a border=0 href=https://2009.igem.org/Team:TUDelft/Modeling_Conjugation><img src=https://static.igem.org/mediawiki/2009/2/21/TUDconjmodelani.gif></a>
 +
</td>
 +
</tr>
 +
</table>
 +
</html>
{{Template:TUDelftiGEM2009_end}}
{{Template:TUDelftiGEM2009_end}}

Latest revision as of 16:43, 17 October 2009

Modeling

Please select one of the two modeling sections: Transcriptional Cascade Modeling or Conjugation Modeling.