Difference between revisions of "Team:TU Delft/Modeling"

Line 164: Line 164:
 
<figcaption>Fig. 2: csgA-GFPmut3 construct (left), csgA construct (right).</figcaption>
 
<figcaption>Fig. 2: csgA-GFPmut3 construct (left), csgA construct (right).</figcaption>
 
</figure>
 
</figure>
 
 
  
 
<p class="lead">As csgA and GFP are under the same promoter, the fluorescence of GFP can be related to the production of csgA. It is presumed that the promoter activity is independent of the addition of the GFPmut3 gene. </p>
 
<p class="lead">As csgA and GFP are under the same promoter, the fluorescence of GFP can be related to the production of csgA. It is presumed that the promoter activity is independent of the addition of the GFPmut3 gene. </p>
  
 
<p class="lead">As already mentioned in the section “Kinetic experiment” under “Project”, both the fluorescence and OD600 were measured. It is important to also measure the cell density, as we are interested in the fluorescence per cell per second. To calculate how much fluorescence corresponds to a certain amount of GFP (in nanogram), a calibration curve has been constructed with wtGFP. As the ratio of emission maxima for GFPmut3 relative to wtGFP is 21 at an excitation of 488nm (Brendan et. al), the correct mass amount of GFP mut3 could be determined. </p>
 
<p class="lead">As already mentioned in the section “Kinetic experiment” under “Project”, both the fluorescence and OD600 were measured. It is important to also measure the cell density, as we are interested in the fluorescence per cell per second. To calculate how much fluorescence corresponds to a certain amount of GFP (in nanogram), a calibration curve has been constructed with wtGFP. As the ratio of emission maxima for GFPmut3 relative to wtGFP is 21 at an excitation of 488nm (Brendan et. al), the correct mass amount of GFP mut3 could be determined. </p>
 +
 +
<figure><img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/f/f3/Modelling_pic_8.png" style="width:100%; background-size: cover;" alt="Generic placeholder image"><figcaption>Fig. 3: OD600 for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose (for both the csgA-GFPmut3 and csgA constructs). </figcaption></figure>
 +
 +
<p class="lead">In Fig. 3, the OD600 is plotted for each of the samples in the kinetic experiment. As can be seen from Fig. 3, there is some difference in OD600 between the samples. Therefore, the fluorescence signal normalized by the number of cells as shown in Fig. 4:</p>
 +
 +
<figure><img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/f/f3/Modelling_pic_8.png" style="width:100%; background-size: cover;" alt="Generic placeholder image"><figcaption>Fig. 3: OD600 for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose (for both the csgA-GFPmut3 and csgA constructs). </figcaption></figure>
 +
 +
<p class="lead">In Fig. 3, the OD600 is plotted for each of the samples in the kinetic experiment. As can be seen from Fig. 3, there is some difference in OD600 between the samples. Therefore, the fluorescence signal normalized by the number of cells as shown in Fig. 4:</p>
 +
 +
<figure><img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/d/d3/Modelling_pic_9_.png" style="width:100%; background-size: cover;" alt="Generic placeholder image"><figcaption>Fig. 4: Fluorescence signal normalized by the number of cells for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose with the csgA and csgA-GFPmut3 construct. The error bars are included for all experiments. </figcaption></figure>
 +
 +
<p class="lead">As can be seen from Fig. 4, GFPmut3 only appears if the culture has been induced by rhamnose. Furthermore, it can be seen that a higher induction level of rhamnose leads to an increase in GFPmut3 and thus fluorescence. Finally, as the fluorescence signal is normalized by the cell density, one can make statements about the activity of the rhamnose promoter. The promoter seems to not be active right after induction, but more after 3 or 4 hours.  This is in accordance with data from literature (Wegerer et. al), in which a low amount of fluorescence with a rhamnose promoter was observed after 2 hours of induction. </p>
 +
 +
<p class="lead">The calibration line of fluorescence versus mass amount GFPmut3 is given in Fig. 5. </p>
 +
 +
<p class="lead"> corresponding function of the GFPmut3 calibration line with <i>mass<sub>GFP</sub></i> in ng is:</p>
 +
<img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/7/79/Modelling_eq_1.PNG" style="width:100%; background-size: cover;" alt="Generic placeholder image">
 +
 +
<p class="lead">With the calibration line, which was made with the exact same settings as the kinetic experiment, the mass amount of GFPmut3 per cell can be calculated. With Eq. 1, the fluorescence signal was converted to mass amount of GFPmut3 for rhamnose induction of 0.2% (w/v) and 0.5% (w/v). </p>
 +
 +
<p class="lead">After this conversion, the units are still in nanogram per OD600. As the desired unit is molecules per cells, the following dimension analysis has been performed:</p>
 +
 +
<img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/d/d8/Modelling_eq_2.PNG" style="width:100%; background-size: cover;" alt="Generic placeholder image">
 +
 +
<p class="lead">The result of this unit conversion is shown in Fig. 6:</p>
 +
 +
<figure><img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/a/ae/Modelling_pic_11.png" style="width:100%; background-size: cover;" alt="Generic placeholder image"><figcaption>Fig. 6: Calibration line of the fluorescence signal of GFPmut3 versus its amount in molecules / cell.</figcaption></figure>
 +
 +
<p class="lead">The GFPmut3 steady state concentrations are thus:</p>
 +
<img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/c/cd/Modelling_eq_3.PNG" style="width:100%; background-size: cover;" alt="Generic placeholder image">
 +
 +
<p class="lead">As the data is now converted to molecules per cell, it is possible to construct a model with the purpose of predicting the same steady state levels as in the performed experiments. The only unknown remaining in this model is the promoter activity, which can be varied to correctly fit the data.  </p>
 +
 +
<p class="lead">The model that was made is a mathematical model of ordinary differential equations (ODE’s) describing the GFP mRNA concentration (Eq. 2), immature GFP protein concentration (Eq. 3) and mature GFP protein concentration (Eq. 4) in Table 1. All concentrations are modeled per cell per second. The immature protein is the nonfluorescent GFP, the mature protein is the fluorescent protein (which will be fitted to the data). </p>
 +
 +
<img class="featurette-image img-responsive center-block" src="https://static.igem.org/mediawiki/2015/4/46/Modelling_table_1.PNG" style="width:100%; background-size: cover;" alt="Generic placeholder image">
  
 
   <a href="#accordion1gene" class="btn btn-info" role="button" onclick="$('html,body').animate({scrollTop: $('#accordion1gene').offset().top}, 'slow');" style="text-decoration:none; color:white; float:right;">Back to Top</a>
 
   <a href="#accordion1gene" class="btn btn-info" role="button" onclick="$('html,body').animate({scrollTop: $('#accordion1gene').offset().top}, 'slow');" style="text-decoration:none; color:white; float:right;">Back to Top</a>

Revision as of 11:54, 17 September 2015

Modeling

..

Overview

Our main goal is to make a reproducible biofilm. The strength of the biofilm is determined by the degree of intercellular connectivity through curli fibers. With modeling, it is possible to determine which factors have a strong influence on the intercellular connectivity. For instance, one could argue that a higher csgB nucleator production would lead to more curli fibers and therefore an improved connectivity. But if the csgA production would be limiting, solely short curli chains would be formed possibly preventing the cells to even interconnect.

Generic placeholder image
Fig. 1: Illustration depicting the cells interconnecting by curli fibers. The csgA protein is secreted by the cells to the extracellular space. A curli can be formed when the csgA reacts with the nucleator protein csgB, located at the outer membrane of the cells.

The modeling can be divided in 4 sections. In section 1, the csgA production rate, intracellular csgA concentration and csgB membrane concentration are determined. As our csgA production is controlled under the induction of rhamnose, all rates and concentrations are calculated for two levels of induction (0.2% (w/v) and 0.5% (w/v)) rhamnose. In section 2, the rates and concentrations are used in a grid model which is able to make an estimate on the characteristic time of curli formation. In section 3 the characteristic time for curli formation is used to predict the strength of the biofilm. Finally, in section 4, an application in MATLAB is presented able to calculate the printing time of a certain figure or shape with our Biolinker printer. To go to either of these sections, click on one of the buttons below or select your section of interest in the modeling submenu!

Generic placeholder image Generic placeholder image Generic placeholder image Generic placeholder image

Gene Modeling

Subtitle or summary goes here. Should be short - two or three sentences.

Text1 Back to Top
Text2 Back to Top
text3 Back to Top

csgA production rate

The csgA production rate is defined as the production of csgA per second per cell. To be able to calculate this number, the following unknown parameters need to be characterized:

I. Activity of the promoter (with different levels of rhamnose induction).

II. Export rate of csgA to the extracellular space.

If these parameters are known, a model for the csgA production can be formulated and the export of csgA to the extracellular space can be calculated.

The activity of the promoter, denoted as PoPs, can be calculated by a kinetic experiment with GFP. For this experiment, we made a construct with GFPmut3 (modified GFP variant with stronger fluorescent signal) and csgA under the same rhamnose promoter. The constructs used for this experiment are depicted in Fig. 2.

Generic placeholder image
Generic placeholder image
Fig. 2: csgA-GFPmut3 construct (left), csgA construct (right).

As csgA and GFP are under the same promoter, the fluorescence of GFP can be related to the production of csgA. It is presumed that the promoter activity is independent of the addition of the GFPmut3 gene.

As already mentioned in the section “Kinetic experiment” under “Project”, both the fluorescence and OD600 were measured. It is important to also measure the cell density, as we are interested in the fluorescence per cell per second. To calculate how much fluorescence corresponds to a certain amount of GFP (in nanogram), a calibration curve has been constructed with wtGFP. As the ratio of emission maxima for GFPmut3 relative to wtGFP is 21 at an excitation of 488nm (Brendan et. al), the correct mass amount of GFP mut3 could be determined.

Generic placeholder image
Fig. 3: OD600 for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose (for both the csgA-GFPmut3 and csgA constructs).

In Fig. 3, the OD600 is plotted for each of the samples in the kinetic experiment. As can be seen from Fig. 3, there is some difference in OD600 between the samples. Therefore, the fluorescence signal normalized by the number of cells as shown in Fig. 4:

Generic placeholder image
Fig. 3: OD600 for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose (for both the csgA-GFPmut3 and csgA constructs).

In Fig. 3, the OD600 is plotted for each of the samples in the kinetic experiment. As can be seen from Fig. 3, there is some difference in OD600 between the samples. Therefore, the fluorescence signal normalized by the number of cells as shown in Fig. 4:

Generic placeholder image
Fig. 4: Fluorescence signal normalized by the number of cells for 0% (w/v), 0.2% (w/v) and 0.5% (w/v) rhamnose with the csgA and csgA-GFPmut3 construct. The error bars are included for all experiments.

As can be seen from Fig. 4, GFPmut3 only appears if the culture has been induced by rhamnose. Furthermore, it can be seen that a higher induction level of rhamnose leads to an increase in GFPmut3 and thus fluorescence. Finally, as the fluorescence signal is normalized by the cell density, one can make statements about the activity of the rhamnose promoter. The promoter seems to not be active right after induction, but more after 3 or 4 hours. This is in accordance with data from literature (Wegerer et. al), in which a low amount of fluorescence with a rhamnose promoter was observed after 2 hours of induction.

The calibration line of fluorescence versus mass amount GFPmut3 is given in Fig. 5.

corresponding function of the GFPmut3 calibration line with massGFP in ng is:

Generic placeholder image

With the calibration line, which was made with the exact same settings as the kinetic experiment, the mass amount of GFPmut3 per cell can be calculated. With Eq. 1, the fluorescence signal was converted to mass amount of GFPmut3 for rhamnose induction of 0.2% (w/v) and 0.5% (w/v).

After this conversion, the units are still in nanogram per OD600. As the desired unit is molecules per cells, the following dimension analysis has been performed:

Generic placeholder image

The result of this unit conversion is shown in Fig. 6:

Generic placeholder image
Fig. 6: Calibration line of the fluorescence signal of GFPmut3 versus its amount in molecules / cell.

The GFPmut3 steady state concentrations are thus:

Generic placeholder image

As the data is now converted to molecules per cell, it is possible to construct a model with the purpose of predicting the same steady state levels as in the performed experiments. The only unknown remaining in this model is the promoter activity, which can be varied to correctly fit the data.

The model that was made is a mathematical model of ordinary differential equations (ODE’s) describing the GFP mRNA concentration (Eq. 2), immature GFP protein concentration (Eq. 3) and mature GFP protein concentration (Eq. 4) in Table 1. All concentrations are modeled per cell per second. The immature protein is the nonfluorescent GFP, the mature protein is the fluorescent protein (which will be fitted to the data).

Generic placeholder image Back to Top

Curli Formation

Subtitle or summary goes here. Should be short - two or three sentences.

Biofilm Formation

Subtitle or summary goes here. Should be short - two or three sentences.

Policy and Practice

Subtitle or summary goes here. Should be short - two or three sentences.