Team:Purdue/Modeling

From 2008.igem.org

Revision as of 19:54, 24 October 2008 by Cbarcus (Talk | contribs)
PUlogo.jpg

Contents

Modeling Objectives

  1. Develop a mechanistics ODE model of the population to predict gene expression dynamics
  2. Question? .... Will the color production be fast enough to be useful to the user? Or will it be too late?
  3. What is the relationship between UV exposure and reporter gene expression?
  4. Can we construct a useful calibration curve of color as a function of UV?


Modeling References

  1. 1997. "Mathematical model of the SOS response regulation of an excision repair deficient mutant of Escherichia coli after UV light irradation". [http://www2.lib.purdue.edu:2184/science?_ob=ArticleURL&_udi=B6WMD-45KKS62-5S&_user=29441&_coverDate=05%2F21%2F1997&_rdoc=1&_fmt=&_orig=search&_sort=d&view=c&_version=1&_urlVersion=0&_userid=29441&md5=7297869ffc31d1edfcd20856301793a5]. Off-the-shelf mechanistic model of the SOS response. Utilized as the basic model for our system.
  2. 2005. "Response times and mechanisms of SOS induction by attaching promoters to GFP: "Precise Temporal Modulation in the Response of the SOS DNA Repair Network in Individual Bacteria" [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1151601]. Potential Validation Data Set. Model should predict similar dynamics.


Modeling the SOS response in a uvr- mutant (No nucleotide excision repair)

Assumptions

  1. The UV light intensity is constant and instantaneous.
  2. The bacteria are not undergoing any type of DNA repair at the time of UV exposure.
  3. Thymine dimer formation is the major DNA damage occurring.


The first figure below shows the response to a 5 J/m^2 UV irradiation. We see that the concentration of bound LexA drops considerably within the first four minutes. This also correlates with the concentration of activated RecA (RecA*) going up appreciably. After approximately 60 minutes, the concentration of RecA returns to a normal level. We therefore consider this the "stopping" point of SOS.

The problems with this model include:

  1. The model is based on an instantaneous irradiation of UV at 5 J/m^2.
  2. The model does not account for continuous UV exposure.
  3. The model does not account for any other proteins/genes that may be involved in SOS (ie. SulA).

The benefits of this model include:

  1. Giving us a mathematical, manipulateable model to mend for our purposes.
  2. Showing the general trend of how SOS behaves.
  3. Gives us a time frame for how our color reporters need to work to be feasible.


The second figure provides us with a close-up view of the increase in RecA* and decreases of LexA and RecA. We see that an artifact occurs around the 30 second mark. This artifact is related to the step-wise function that occurs in the model. Over the larger time frame, this artifact is negligible.

Since we have no quantitative data on the activation and deactivation of SOS, we must make an assumption as to when SOS truly starts to take effect. We will consider the intersection of LexA and RecA* to be the initial time when SOS can start repairing the damage accrued by UV irradiation.

All Calculations and the figures above were performed in MathCad, utilizing built in Runge-Kutta 4th order function, by Craig Barcus, utilizing the mathematical model presented by SV. Aksenov in 1997.


PurdueFullModel.jpg
PurdueIntersectionsModel.jpg


Modeling the cleavage of X-gal by Beta-galactosidase

Assumptions

  1. The bacteria have been on the X-gal plate sufficiently long for the X-gal to be in equilibrium with the surface and the cell. This means that the cell has the same concentration of X-gal within its cell wall as the surface outside.
  2. There is no diffusion limitation with this system. As soon as an enzyme is produced, it immediately can bind and cleave X-gal.
  3. The enzymes follow a Michaelis-Menten type enzymatic reaction. Km for this system is 0.2 mM and Vm was calculated to be 4.826 M/min. (Sharp et al. 1969).
  4. Kcat was calculated at 1.52*10^(-20) moles X-Gal / molecule beta-Gal*min. This is the rate of X-gal Cleava
  5. The initial concentration of X-gal in the cell is: 7.48*10^(-16) moles.
  6. Half-life of beta-gal is 60 minutes. After the half-life, the enzyme no longer functions properly and is recycled. (Bachmair et al. 1986).

Schemes for the different models

For the three graphs below, the following scheme was utilized in Microsoft Excel:

  1. The enzyme will only function for 60 minutes, therefore, any complete cleavage time over 60 minutes requires that the enzymes acting for the first 60 minutes of their life will not act after that.
  2. The cleavage of X-gal is as follows: Rate of X-gal Cleavage*number of enzyme molecules*time / Initial X-gal cleavage. This gives us a percentage which is easier to represent and analyze visually.
  3. There are three different rates of formation of the enzyme. Data could not be found for the rate of formation, so different values were utilized.
X gal cleave 40 min half life.jpg


Home The Team The Project Parts Submitted to the Registry Modeling Notebook Tools and References