Team:Newcastle/Modelling/SporulationTuning

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:

Equation 1

Equation 2

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]

Equation 4

Equation 5

Equation 6

Equation 7

Equation 8

Equation 9

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

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

Equation 10

Equation 11

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.





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.

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

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

The CDS lengths was found on the 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:
 * Sporulation Tuning Model (SBML and COPASI)

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

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.

Figure 1.2

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

Figure 1.3

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.

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.

Figure 2.2

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



Figure 2.3

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