Team:Newcastle/Modelling/SporulationTuning

From 2009.igem.org


Sporulation Tuning Model

The expression of KinA has been modelled as seen above, therefore we can now proceed further into the Sporulation Tuning Model, which is built from the KinA Expression Model with COPASI.

To proceed with the modelling of our Sporulation Tuning Model, we have decided that in response to an unidentified stimuli, where KinA autophosphorylates and then donates its phosphate groups to the response regulator Spo0F, the unidentified stimuli will be termed as 'sporulation signal'.[5]

Also, to simplify the modelling, the phosphatases that provide access for negative signals to influence the cell’s decision whether to sporulate or to continue vegetative growth,[4] such as Spo0E[1], YisI, and YnzD,[5] are not taken into consideration while building this model.

The following equations describe the model:

TeamNewcastleSporeTuneEqn1.png
Equation 1


TeamNewcastleSporeTuneEqn2.png
Equation 2


TeamNewcastleSporeTuneEqn3.png
Equation 3


As Spo0F lacks an output domain and is incapable of activating transcription; it serves only as an intermediary in the phosphorelay. The phosphotransferase Spo0B transfers the phosphate from Spo0F~P to Spo0A.[4]

TeamNewcastleSporeTuneEqn4.png
Equation 4


TeamNewcastleSporeTuneEqn5.png
Equation 5


TeamNewcastleSporeTuneEqn6.png
Equation 6


TeamNewcastleSporeTuneEqn7.png
Equation 7


TeamNewcastleSporeTuneEqn8.png
Equation 8


TeamNewcastleSporeTuneEqn9.png
Equation 9


With respect to Equations 1 to 9, as seen above, the corresponding fluxes are as follows:

TeamNewcastleSporeTuneFluxes2.png
Fluxes


The phophorylated proteins Spo0FP, Spo0BP and Spo0AP will go through degradation, and the following equations, Equation 10, 11 and 12 describes degradation.

TeamNewcastleSporeTuneEqn10.png
Equation 10


TeamNewcastleSporeTuneEqn11.png
Equation 11


TeamNewcastleSporeTuneEqn12.png
Equation 12


As parameters such as transcription and translation rates could not be found in literature, the following equations were used to calculate the transcription and translation rates whenever possible.

TeamNewcastleSporeTuneTranscriptionRate.png


TeamNewcastleSporeTuneTranslationRate.png


Using the transcription and translation rate equations as seen above, the rates were calculated for the components Spo0A, Spo0B, Spo0F and KinA, and the results are reported in the table below.

TeamNewcastleSporeTuneTranscriptionTranslateRate.png
Table 1: Transcription and Translation Rates


For components where information on the CDS length could not be found, we used the value 0.1 for both the translation and transcription rate.

Also, the values for the various rate constants could not be found in literature either. Therefore, the following values were used:

kforward = 0.001

kreverse = 0.05

Other parameters and formulas which were used in this model are as follows:

Protein Degradation Rate = 0.0012

mRNA Degradation Rate = 0.0058


TeamNewcastleSporeTuneFunctions.png
Table 2: Functions used in the KinA Expression and Sporulation Tuning Models


The CDS lengths was found on the [http://www.ncbi.nlm.nih.gov/Genbank/ GenBank] using the accession number, AL009126.

These parameters and functions were also used in the KinA Expression Model.

The Sporulation Tuning Model is available for download as follows:

Results

At a constant IPTG concentration of 5000nmol/fl, we vary the concentration of the sporulation signal, [SS] and the following graphs show the results.

At [SS] = 1000nmol/fl

TeamNewcastleSporeTunePic1.png
Figure 1.1


Figure 1.1 shows that at an IPTG concentration of 5000nmol/fl and sporulation signal concentration of 1000nmol/fl, the Spo0F concentration gets significantly higher than the other components.

In order to take a better look at the behaviour of the rest of the components, the following graphs were produced.


TeamNewcastleSporeTunePic2.png
Figure 1.2


From Figure 1.2, the concentration of Spo0AP seems to be at a healthy concentration.


TeamNewcastleSporeTunePic3.png
Figure 1.3


TeamNewcastleSporeTunePic4.png
Figure 1.4


At [SS] = 3000nmol/fl

The concentration of the sporulation signal, [SS] is increased to 3000nmol/fl in the following graphs, and the results were investigated.

TeamNewcastleSporeTunePic5.png
Figure 2.1


With a stronger sporulation signal, the concentration of Spo0FP and Spo0BP increased significantly, while the concentration of Spo0F dropped significantly from approximately 9000nmol/fl to 2300nmol/fl, comparing Figure 1.1 to 2.1.


TeamNewcastleSporeTunePic6.png
Figure 2.2


The concentration of LacI* remained approximately the same, so as Spo0AP.


TeamNewcastleSporeTunePic7.png


Figure 2.3


TeamNewcastleSporeTunePic8.png
Figure 2.4


Comparing Figure 1.3 and 2.4, the concentration of Spo0A falls from approximately 100nmol/fl to 50nmol/fl when [SS] is increased.


Referring to the graphs obtained from the model, it seems like increasing the sporulation signal concentration, does not increase the Spo0AP concentration. However, the concentrations of Spo0BP and SpoFP did increase as the sporulation signal concentration increased.


Other Models

See KinA Expression Model
See Sin Operon Model





News

Events

Social Net

  • Newcastle iGEM Twitter
  • [http://www.facebook.com/home.php#/group.php?gid=131709337641 Newcastle on Facebook]
  • [http://www.youtube.com/user/newcastle2009igem Newcastle Youtube Channel]