Team:Macquarie Australia/Project/Model

From 2014.igem.org

(Difference between revisions)
 
(23 intermediate revisions not shown)
Line 11: Line 11:
tr {
tr {
border: 1px solid;
border: 1px solid;
-
margin-right: auto;
+
text-align: left;
 +
margin: 5px;
}
}
td {
td {
border: 1px solid;
border: 1px solid;
-
margin-right: auto;
+
text-align: left;
 +
margin: 5px;
}
}
</style>
</style>
Line 37: Line 39:
<div class="cont-out">
<div class="cont-out">
-
<h3>Modelling</h3>
 
-
<p>Jumping the GUN: Modelling GUN4 binding to protoporphyrin IX in the heme-chlorophyll biosynthetic pathway.</p>
+
<img src="https://static.igem.org/mediawiki/2014/9/97/Modeling2.jpg" width=600/>
-
<p>Modelling biological and biochemical systems in silico is a central part of synthetic biology. The use of models to predict, test, and rationalise experimental observations is key in understanding the complex systems present in cells, even more so when those systems are being deliberately modified. </p>
+
<h3>Modelling - Jumping the GUN</h3>
 +
<br/>
-
<p>The issue when starting the modelling component of a project is often finding something to model, and finding something which will give useful results. When we set out to construct the chlorophyll pathway in our cells, we had a number of questions. How much chlorophyll are we going to end up making? How is the gene expression going to affect and result in that production? Is the chlorophyll going to kill the E. coli? Is the hydrogen production we’re proposing feasible? Questions from all sides of our project ended up having some component we could model, making the problem deciding what was the most important, central, and above all feasible part to construct mathematically.</p>
+
<h4>Modelling GUN4 binding to protoporphyrin IX in the heme-chlorophyll biosynthetic pathway</h4>
 +
 
 +
<p>Modelling biological and biochemical systems <i>in silico</i> is a central part of synthetic biology. The use of models to predict, test, and rationalise experimental observations is key in understanding the complex systems present in cells, even more so when those systems are being deliberately modified. </p>
 +
 
 +
<p>The issue when starting the modelling component of a project is often finding something to model, and finding something which will give useful results. When we set out to construct the chlorophyll pathway in our <i>E. coli</i> cells, we had a number of questions. How much chlorophyll are we going to end up making? How is the gene expression going to affect and result in that production? Is the chlorophyll going to kill the <i>E. coli</i>? Is the hydrogen production we’re proposing feasible? Questions from all sides of our project ended up having some component we could model, making the problem deciding what was the most important, central, and above all feasible part to construct mathematically.</p>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2014/3/32/Modelfig1.png" width=700 />
 +
<p><i><b>Figure 1</b> The first part of the Chlorophyll biosynthesis pathway. Our modelling examined the effect of co-factor GUN4 on the incorporation of Mg into Protoporphyrin IX by Mg-Chelastase (Operon 1).</i></p>
<h4>Why GUN4?</h4>
<h4>Why GUN4?</h4>
-
<p>Our pathway begins at protoporphyrin IX (PPIX), an intermediate in the synthesis of heme in E. coli. Here, we’ve constructed a branch point, and are siphoning off this molecule from this pathway onto the chlorophyll production line. The amount of PPIX which we redirect onto the chlorophyll pathway dictates almost everything about our project, including our end-product formation, whether the gene product enzymes will function and to what extent, and how much gene expression to promote in order to have a functioning biosynthesis pathway. This redirection is primarily dependent on the binding of the GUN4 protein to PPIX, which makes it inaccessible to the ferrochelatase enzyme which begins the heme pathway, and allows the magnesium chelatase to insert the magnesium ion to make magnesium protoporphyrin IX (Mg-PPIX). Since this is the “point of no return”, any amounts of PPIX which are metallated by magnesium chelatase are committed to the synthesis of chlorophyll, meaning all further questions about the project depend on this one branch step. Due to this importance, it is the system around this point which we have modelled.</p>
+
<p>Our pathway begins at protoporphyrin IX (PPIX), an intermediate in the synthesis of heme in <i>E. coli</i>. Here, we’ve constructed a branch point, and are siphoning off this molecule from this pathway onto the chlorophyll production line (Figure 1). The amount of PPIX which we redirect onto the chlorophyll pathway dictates almost everything about our project, including our end-product formation, whether the gene product enzymes will function and to what extent, and how much gene expression to promote in order to have a functioning biosynthesis pathway. This redirection is primarily dependent on the binding of the GUN4 protein to PPIX, which makes it inaccessible to the ferrochelatase enzyme which begins the heme pathway, and allows the magnesium chelatase to insert the magnesium ion yielding magnesium protoporphyrin IX (Mg-PPIX). Since this is the “point of no return”, any amounts of PPIX which are metalated by magnesium chelatase are committed to the synthesis of chlorophyll, meaning all further questions about the project depend on this one branch step. Due to this importance, it is the system around this point which we have modelled.</p>
<h4>Constructing the Equations</h4>
<h4>Constructing the Equations</h4>
<p>The first step in modelling our system was to find out more about its biochemistry, and let that knowledge dictate the mathematics. In the case of our system, there are two main catalytic equilibria involved, namely the ferrochelatase-catalysed production of heme, and the magnesium chelatase-catalysed production of Mg-PPIX. Linking these two pathways is the binding equilibrium of GUN4 with PPIX, with the detachment of GUN4 from Mg-PPIX the final step to be modelled. This product would then go on to react with the next enzyme of our pathway and continue to make chlorophyll, and is thus a quantity of interest.</p>
<p>The first step in modelling our system was to find out more about its biochemistry, and let that knowledge dictate the mathematics. In the case of our system, there are two main catalytic equilibria involved, namely the ferrochelatase-catalysed production of heme, and the magnesium chelatase-catalysed production of Mg-PPIX. Linking these two pathways is the binding equilibrium of GUN4 with PPIX, with the detachment of GUN4 from Mg-PPIX the final step to be modelled. This product would then go on to react with the next enzyme of our pathway and continue to make chlorophyll, and is thus a quantity of interest.</p>
-
<p>The reaction scheme in Figure 1 shows the main equations of the pathway which were included in the model. In addition to this, there is a feedback loop from protoheme to the production of PPIX; when the concentration of protoheme gets above a certain value, the synthesis of PPIX is turned off, while when the value drops below another value, the synthesis is started again. </p>
+
<img src="https://static.igem.org/mediawiki/2014/a/a1/Modelfig2.png" width=700 />
 +
<p><i><b>Figure 2 </b> Equations used (<i>top</i>) to model incorporation of Mg into Protoporphyrin by GUN4 (co-factor) (<i>bottom</i>). Abbreviations: Protoporphyrin (PPIX); Magnesium Protoporphyrin IX (PPIX); Ferrochelatase (Ferro); Magnesium Chelatase (MgChel).</i></p>
-
<p>Due to the insertion and expression of the GUN4 gene, there will also be some level of production of GUN4. However, this is likely to cease above a certain concentration due to transcription of the operon. These two factors are then taken into account in the ordinary differential equations (ODEs) which are presented in Figure 2, with the code shown in Table 1.</p>
+
<p>The reaction scheme in Figure 2 shows the main equations of the pathway which were included in the model. In addition to this, there is a feedback loop from protoheme to the production of PPIX; when the concentration of protoheme gets above a certain value, the synthesis of PPIX is turned off, while when the value drops below another value, the synthesis is started again. </p>
 +
 
 +
<p>Due to the insertion and expression of the GUN4 gene, there will also be some level of production of GUN4. However, this is likely to cease above a certain concentration due to transcription of the Operon 1. These two factors are then taken into account in the ordinary differential equations (ODEs) which are presented in Figure 2, with the code shown in Table 1.</p>
<p>This coupled set of 10 differential equations was then put into MATLAB, with constants and initial conditions as given and justified in Table 2.
<p>This coupled set of 10 differential equations was then put into MATLAB, with constants and initial conditions as given and justified in Table 2.
-
We are primarily concerned with the production of MgPPIX, as well as the kinetics of the system as a whole. Figure 3 presents all ten molecules and complexes involved and their cellular concentrations over a five hour time scale. There are several interesting long-term behaviours present in the plot. First, the product which is being produced the most is in fact the complexed MgPPIX-GUN4, rather than the detached MgPPIX product. In the full reaction scheme, the MgPPIX would be further modified by other enzymes and thus would cause the MgPPIX-GUN4 equilibrium to move towards the right, leading to a decrease in the accumulation of this product. Initially, the PPIX concentration is highly dependent on the feedback mechanism by the protoheme regulation; however, after approximately 90 minutes, this concentration settles to an equilibrium which is more than the initial concentration. This is expected behaviour, as two systems are now relying on this pool of PPIX substrate, leading the cell to produce a higher standing concentration to compensate. Over long periods, the production of MgPPIX-GUN4, PPIX-GUN4, and GUN4 becomes linear, due to the constant production of GUN4 in our in silico cell. The feedback modelling of protoheme seems to function as desired, with a very fast feedback loop creating a very constant protoheme concentration in the cell. </p>
+
We are primarily concerned with the production of MgPPIX, as well as the kinetics of the system as a whole. Figure 3 presents all ten molecules and complexes involved and their cellular concentrations over a five hour time scale. There are several interesting long-term behaviours present in the plot. First, the product which is being produced the most is in fact the complexed MgPPIX-GUN4, rather than the detached MgPPIX product. In the full reaction scheme, the MgPPIX would be further modified by other enzymes and thus would cause the MgPPIX-GUN4 equilibrium to move towards the right, leading to a decrease in the accumulation of this product. Initially, the PPIX concentration is highly dependent on the feedback mechanism by the protoheme regulation; however, after approximately 90 minutes, this concentration settles to an equilibrium which is more than the initial concentration. This is expected behaviour, as two systems are now relying on this pool of PPIX substrate, leading the cell to produce a higher standing concentration to compensate. Over long periods, the production of MgPPIX-GUN4, PPIX-GUN4, and GUN4 becomes linear, due to the constant production of GUN4 in our <i>in silico</i> cell. The feedback modelling of protoheme seems to function as desired, with a very fast feedback loop creating a very constant protoheme concentration in the cell. </p>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2014/2/20/Modelfig3.png" width=700 />
 +
<p><i><b>Figure 3 </b>Modelling of concentration of all 10 molecules/proteins/protein complexes produced in chlorophyll biosynthesis pathway over a 5 hour time period.</i></p>
<p>Figure 4 shows in more detail the concentrations of the initial PPIX substrate, as well as our final desired product of the system, the MgPPIX. It is important to note that the fluctuating behaviour around 150-200 minutes is likely due to the time sampling of MATLAB rather than a behaviour of the system. This shows that the PPIX concentration indeed settles to an equilibrium, while the MgPPIX concentration follows quasi-hyperbolic kinetics.</p>
<p>Figure 4 shows in more detail the concentrations of the initial PPIX substrate, as well as our final desired product of the system, the MgPPIX. It is important to note that the fluctuating behaviour around 150-200 minutes is likely due to the time sampling of MATLAB rather than a behaviour of the system. This shows that the PPIX concentration indeed settles to an equilibrium, while the MgPPIX concentration follows quasi-hyperbolic kinetics.</p>
 +
 +
<img src="https://static.igem.org/mediawiki/2014/c/cc/Modelfig4.png" width=700 />
 +
<p><i><b>Figure 4 </b>Modelling of initial concentration of PPIX substrate and final product (MgPPIX). Note: fluctuation likely due to sampling by MATLAB and not representative of biology of system.</i></p>
<h4>Linking to Experiment</h4>
<h4>Linking to Experiment</h4>
<p>Initially, there were high hopes to measure the concetration of MgPPIX as a tool to verify the model. However, as the project neared its end, the best point of comparison to functional assays was to measure the impact of GUN4 only. This meant removing the MgChel component of the pathway, and looking at the measurable comparison between PPIX and the complexed PPIX-GUN4. The results of the model over an hour is shown in Figure 5. This illustrates that the assumption of constant production rate of GUN4 is likely invalid, due to a lack of equilibrium concentration of PPIX and thus an unlikely result for the system logically. From functional assays, the difference in PPIX concentration before and after the addition of GUN4 could be measured, and the model adjusted accordingly.</p>
<p>Initially, there were high hopes to measure the concetration of MgPPIX as a tool to verify the model. However, as the project neared its end, the best point of comparison to functional assays was to measure the impact of GUN4 only. This meant removing the MgChel component of the pathway, and looking at the measurable comparison between PPIX and the complexed PPIX-GUN4. The results of the model over an hour is shown in Figure 5. This illustrates that the assumption of constant production rate of GUN4 is likely invalid, due to a lack of equilibrium concentration of PPIX and thus an unlikely result for the system logically. From functional assays, the difference in PPIX concentration before and after the addition of GUN4 could be measured, and the model adjusted accordingly.</p>
-
<table>
+
<img src="https://static.igem.org/mediawiki/2014/7/77/Modelfig5.png" width=700 />
-
<tr>
+
<p><i><b>Figure 5</b> Modelling the effect of GUN4 on first step of chlorophyll biosynthesis pathway. Modelling measured the concentration of Protoporphyrin IX (PPIX) produced be removing Mg Chelatase from pathway. Modelling assumed a constant production rate of GUN4.</i></p>
-
<td>Reagent</td>
+
 
-
<td>Code</td>
+
<table class="tableizer-table">
 +
<tr class="tableizer-firstrow">
 +
<th>Reagent</th>
 +
<th>Code</th>
</tr>
</tr>
<tr>
<tr>
<td>Protoporphyrin IX (PPIX)</td>
<td>Protoporphyrin IX (PPIX)</td>
-
<td>S1</td>
+
<td>S<sub>1</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>Ferrochelatase (Ferro)</td>
<td>Ferrochelatase (Ferro)</td>
-
<td>E1</td>
+
<td>E<sub>1</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>PPIX-Ferro</td>
<td>PPIX-Ferro</td>
-
<td>E1S1</td>
+
<td>E<sub>1</sub>S<sub>1</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>Protoheme</td>
<td>Protoheme</td>
-
<td>P1</td>
+
<td>P<sub>1</sub></td>
</tr>
</tr>
<tr>
<tr>
Line 90: Line 111:
<tr>
<tr>
<td>PPIX-GUN4</td>
<td>PPIX-GUN4</td>
-
<td>S¬2</td>
+
<td>S<sub>2</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>Magnesium chelatase (MgChel)</td>
<td>Magnesium chelatase (MgChel)</td>
-
<td>E2</td>
+
<td>E<sub>2</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>PPIX-GUN4-MgChel</td>
<td>PPIX-GUN4-MgChel</td>
-
<td>E2S2</td>
+
<td>E<sub>2</sub>S<sub>2</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>Magnesium protoporphyrin IX (MgPPIX)-GUN4</td>
<td>Magnesium protoporphyrin IX (MgPPIX)-GUN4</td>
-
<td>P2</td>
+
<td>P<sub>2</sub></td>
</tr>
</tr>
<tr>
<tr>
<td>MgPPIX</td>
<td>MgPPIX</td>
-
<td>Pf</td>
+
<td>P<sub>f</sub></td>
</tr>
</tr>
</table>
</table>
 +
<p><b>Table 1.</b> Acronyms of the components of the model.</p>
-
Table 2
+
<table class="tableizer-table">
 +
<tr class="tableizer-firstrow">
 +
<th>Parameter Symbol</th>
 +
<th>Meaning</th>
 +
<th>Value</th>
 +
<th>Reference / Comments</th>
 +
</tr>
 +
<tr>
 +
<td>K<sub>D1</sub></td>
 +
<td>Dissociation constant of ferrochelatase with PPIX</td>
 +
<td>110 nM/min </td>
 +
<td>Ref. 1</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>D2</sub></td>
 +
<td>Dissociation constant of GUN4 and PPIX</td>
 +
<td>11 nM/min</td>
 +
<td>Ref. 2</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>D3</sub></td>
 +
<td>Dissociation constant of magnesium chelatase and PPIX-GUN4</td>
 +
<td>10 nM/min</td>
 +
<td>Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>D4</sub></td>
 +
<td>Dissociation constant of MgPPIX and GUN4</td>
 +
<td>1.6 nM/min</td>
 +
<td>Ref. 2</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>cat,1</sub></td>
 +
<td>Turnover number for ferrochelatase for PPIX</td>
 +
<td>0.78 /min</td>
 +
<td>Ref. 1</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>cat,2</sub></td>
 +
<td>Turnover number for magnesium chelatase for PPIX-GUN4</td>
 +
<td>0.8 /min</td>
 +
<td>Ref. 3 & 4</td>
 +
</tr>
 +
<tr>
 +
<td>K<sub>f</sub></td>
 +
<td>Rate of formation of PPIX due to protoheme regulation</td>
 +
<td>20 nM/min</td>
 +
<td>Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).</td>
 +
</tr>
 +
<tr>
 +
<td>χ<sub>f</sub></td>
 +
<td>Switch triggering PPIX production when protoheme concentration drops below 50 nM</td>
 +
<td>1 (when [Protoheme] < 50 nM), 0 (when [Protoheme] > 50 nM)</td>
 +
<td>Thorough lit search yielded no switch value; estimated from Willows, R. 2014 (unpublished).</td>
 +
</tr>
 +
<tr>
 +
<td>β<sub>deg</sub></td>
 +
<td>Degradation rate/rate at which protoheme is removed by conversion to heme</td>
 +
<td>1 nM/min</td>
 +
<td>Thorough lit search yielded no value; since we haven’t included the formation of heme in our system, this compensates for this process; estimated from Willows, R. 2014 (unpublished).</td>
 +
</tr>
 +
<tr>
 +
<td>α<sub>G</sub></td>
 +
<td>Formation rate of GUN4 from gene</td>
 +
<td>1 nM/min</td>
 +
<td>No value available due to lack of characterisation in E. coli; estimated from Willows, R. 2014 (unpublished).</td>
 +
</tr>
 +
</table>
 +
<p><b>Table 2.</b> Constant values of the model components based on reported literature.</p>
-
Parameter Symbol Meaning Value Reference/Comments
+
<h3>References</h3>
-
KD1 Dissociation constant of ferrochelatase with PPIX 110 nM/min Ref. 1
+
<p><ul style="list-style-type: decimal;">
-
KD2 Dissociation constant of GUN4 and PPIX 11 nM/min Ref. 2
+
<li>Storm, P., Tibiletti, T., Hall, M., Funk, C. Refolding and enzyme kinetic studies on the ferrochelatase of the cyanobacterium Synechocystis sp. PCC 6803. PloS one. 2013. 8.</li>
-
KD3 Dissociation constant of magnesium chelatase and PPIX-GUN4 10 nM/min Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).
+
<li>Adhikari, N.D., Orler, R., Chory, J., Froehlich, J.E., Larkin, R.M. Porphyrins promote the association of GENOMES UNCOUPLED 4 and a Mg-chelatase subunit with chloroplast membranes. JBC. 2009. 284: 24783-96.</li>
-
KD4 Dissociation constant of MgPPIX and GUN4 1.6 nM/min Ref. 2
+
<li>Reid, J.D., Hunter, C.N. Magnesium-dependent ATPase activity and cooperativity of magnesium chelatase from Synechocystis sp. PCC6803. JBC. 2004. 279: 26893-9.</li>
-
kcat,1 Turnover number for ferrochelatase for PPIX 0.78 /min Ref. 1
+
-
kcat,2 Turnover number for magnesium chelatase for PPIX-GUN4 0.8 /min Ref. 3
+
<li>Adams, N.B., Marklew, C.J., Brindley, A.A., Hunter, C.N., Reid, J.D. Characterization of the magnesium chelatase from Thermosynechococcus elongatus. Biochemical Journal. 2014. 457: 163-70. </li>
-
kf Rate of formation of PPIX due to protoheme regulation 20 nM/min Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).
+
</ul></p>
-
χf Switch triggering PPIX production when protoheme concentration drops below 50 nM 1
+
<h3></h3>
-
(when [Protoheme] < 50 nM)
+
-
0
+
-
(when [Protoheme] > 50 nM) Thorough lit search yielded no switch value; estimated from Willows, R. 2014 (unpublished).
+
-
βdeg Degradation rate/rate at which protoheme is removed by conversion to heme 1 nM/min Thorough lit search yielded no value; since we haven’t included the formation of heme in our system, this compensates for this process; estimated from Willows, R. 2014 (unpublished).
+
-
αG Formation rate of GUN4 from gene 1 nM/min No value available due to lack of characterisation in E. coli; estimated from Willows, R. 2014 (unpublished).
+
-
References
 
-
1. Storm, P., Tibiletti, T., Hall, M., Funk, C. Refolding and enzyme kinetic studies on the ferrochelatase of the cyanobacterium Synechocystis sp. PCC 6803. PloS one. 2013. 8.
 
-
2. Adhikari, N.D., Orler, R., Chory, J., Froehlich, J.E., Larkin, R.M. Porphyrins promote the association of GENOMES UNCOUPLED 4 and a Mg-chelatase subunit with chloroplast membranes. JBC. 2009. 284: 24783-96.
 
-
3. Reid, J.D., Hunter, C.N. Magnesium-dependent ATPase activity and cooperativity of magnesium chelatase from Synechocystis sp. PCC6803. JBC. 2004. 279: 26893-9.
 
-
4. Castelfranco, P.A., Jones, O.T. Protoheme turnover and chlorophyll synthesis in greening barley tissue. Plant physiology. 1975. 55: 485-90.
 
-
5. Tanaka, R., Tanaka, A. Tetrapyrrole biosynthesis in higher plants. Annual review of plant biology. 2007. 58: 321-46.
 
-
6. McLean, S., Hunter, C.N. An enzyme-coupled continuous spectrophotometric assay for magnesium protoporphyrin IX methyltransferases. Analytical biochemistry. 2009. 394: 223-8.
 
-
7. Adams, N.B., Marklew, C.J., Brindley, A.A., Hunter, C.N., Reid, J.D. Characterization of the magnesium chelatase from Thermosynechococcus elongatus. Biochemical Journal. 2014. 457: 163-70.
 
</div>
</div>

Latest revision as of 03:02, 18 October 2014

Modelling - Jumping the GUN


Modelling GUN4 binding to protoporphyrin IX in the heme-chlorophyll biosynthetic pathway

Modelling biological and biochemical systems in silico is a central part of synthetic biology. The use of models to predict, test, and rationalise experimental observations is key in understanding the complex systems present in cells, even more so when those systems are being deliberately modified.

The issue when starting the modelling component of a project is often finding something to model, and finding something which will give useful results. When we set out to construct the chlorophyll pathway in our E. coli cells, we had a number of questions. How much chlorophyll are we going to end up making? How is the gene expression going to affect and result in that production? Is the chlorophyll going to kill the E. coli? Is the hydrogen production we’re proposing feasible? Questions from all sides of our project ended up having some component we could model, making the problem deciding what was the most important, central, and above all feasible part to construct mathematically.

Figure 1 The first part of the Chlorophyll biosynthesis pathway. Our modelling examined the effect of co-factor GUN4 on the incorporation of Mg into Protoporphyrin IX by Mg-Chelastase (Operon 1).

Why GUN4?

Our pathway begins at protoporphyrin IX (PPIX), an intermediate in the synthesis of heme in E. coli. Here, we’ve constructed a branch point, and are siphoning off this molecule from this pathway onto the chlorophyll production line (Figure 1). The amount of PPIX which we redirect onto the chlorophyll pathway dictates almost everything about our project, including our end-product formation, whether the gene product enzymes will function and to what extent, and how much gene expression to promote in order to have a functioning biosynthesis pathway. This redirection is primarily dependent on the binding of the GUN4 protein to PPIX, which makes it inaccessible to the ferrochelatase enzyme which begins the heme pathway, and allows the magnesium chelatase to insert the magnesium ion yielding magnesium protoporphyrin IX (Mg-PPIX). Since this is the “point of no return”, any amounts of PPIX which are metalated by magnesium chelatase are committed to the synthesis of chlorophyll, meaning all further questions about the project depend on this one branch step. Due to this importance, it is the system around this point which we have modelled.

Constructing the Equations

The first step in modelling our system was to find out more about its biochemistry, and let that knowledge dictate the mathematics. In the case of our system, there are two main catalytic equilibria involved, namely the ferrochelatase-catalysed production of heme, and the magnesium chelatase-catalysed production of Mg-PPIX. Linking these two pathways is the binding equilibrium of GUN4 with PPIX, with the detachment of GUN4 from Mg-PPIX the final step to be modelled. This product would then go on to react with the next enzyme of our pathway and continue to make chlorophyll, and is thus a quantity of interest.

Figure 2 Equations used (top) to model incorporation of Mg into Protoporphyrin by GUN4 (co-factor) (bottom). Abbreviations: Protoporphyrin (PPIX); Magnesium Protoporphyrin IX (PPIX); Ferrochelatase (Ferro); Magnesium Chelatase (MgChel).

The reaction scheme in Figure 2 shows the main equations of the pathway which were included in the model. In addition to this, there is a feedback loop from protoheme to the production of PPIX; when the concentration of protoheme gets above a certain value, the synthesis of PPIX is turned off, while when the value drops below another value, the synthesis is started again.

Due to the insertion and expression of the GUN4 gene, there will also be some level of production of GUN4. However, this is likely to cease above a certain concentration due to transcription of the Operon 1. These two factors are then taken into account in the ordinary differential equations (ODEs) which are presented in Figure 2, with the code shown in Table 1.

This coupled set of 10 differential equations was then put into MATLAB, with constants and initial conditions as given and justified in Table 2. We are primarily concerned with the production of MgPPIX, as well as the kinetics of the system as a whole. Figure 3 presents all ten molecules and complexes involved and their cellular concentrations over a five hour time scale. There are several interesting long-term behaviours present in the plot. First, the product which is being produced the most is in fact the complexed MgPPIX-GUN4, rather than the detached MgPPIX product. In the full reaction scheme, the MgPPIX would be further modified by other enzymes and thus would cause the MgPPIX-GUN4 equilibrium to move towards the right, leading to a decrease in the accumulation of this product. Initially, the PPIX concentration is highly dependent on the feedback mechanism by the protoheme regulation; however, after approximately 90 minutes, this concentration settles to an equilibrium which is more than the initial concentration. This is expected behaviour, as two systems are now relying on this pool of PPIX substrate, leading the cell to produce a higher standing concentration to compensate. Over long periods, the production of MgPPIX-GUN4, PPIX-GUN4, and GUN4 becomes linear, due to the constant production of GUN4 in our in silico cell. The feedback modelling of protoheme seems to function as desired, with a very fast feedback loop creating a very constant protoheme concentration in the cell.

Figure 3 Modelling of concentration of all 10 molecules/proteins/protein complexes produced in chlorophyll biosynthesis pathway over a 5 hour time period.

Figure 4 shows in more detail the concentrations of the initial PPIX substrate, as well as our final desired product of the system, the MgPPIX. It is important to note that the fluctuating behaviour around 150-200 minutes is likely due to the time sampling of MATLAB rather than a behaviour of the system. This shows that the PPIX concentration indeed settles to an equilibrium, while the MgPPIX concentration follows quasi-hyperbolic kinetics.

Figure 4 Modelling of initial concentration of PPIX substrate and final product (MgPPIX). Note: fluctuation likely due to sampling by MATLAB and not representative of biology of system.

Linking to Experiment

Initially, there were high hopes to measure the concetration of MgPPIX as a tool to verify the model. However, as the project neared its end, the best point of comparison to functional assays was to measure the impact of GUN4 only. This meant removing the MgChel component of the pathway, and looking at the measurable comparison between PPIX and the complexed PPIX-GUN4. The results of the model over an hour is shown in Figure 5. This illustrates that the assumption of constant production rate of GUN4 is likely invalid, due to a lack of equilibrium concentration of PPIX and thus an unlikely result for the system logically. From functional assays, the difference in PPIX concentration before and after the addition of GUN4 could be measured, and the model adjusted accordingly.

Figure 5 Modelling the effect of GUN4 on first step of chlorophyll biosynthesis pathway. Modelling measured the concentration of Protoporphyrin IX (PPIX) produced be removing Mg Chelatase from pathway. Modelling assumed a constant production rate of GUN4.

Reagent Code
Protoporphyrin IX (PPIX) S1
Ferrochelatase (Ferro) E1
PPIX-Ferro E1S1
Protoheme P1
GUN4 G
PPIX-GUN4 S2
Magnesium chelatase (MgChel) E2
PPIX-GUN4-MgChel E2S2
Magnesium protoporphyrin IX (MgPPIX)-GUN4 P2
MgPPIX Pf

Table 1. Acronyms of the components of the model.

Parameter Symbol Meaning Value Reference / Comments
KD1 Dissociation constant of ferrochelatase with PPIX 110 nM/min Ref. 1
KD2 Dissociation constant of GUN4 and PPIX 11 nM/min Ref. 2
KD3 Dissociation constant of magnesium chelatase and PPIX-GUN4 10 nM/min Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).
KD4 Dissociation constant of MgPPIX and GUN4 1.6 nM/min Ref. 2
Kcat,1 Turnover number for ferrochelatase for PPIX 0.78 /min Ref. 1
Kcat,2 Turnover number for magnesium chelatase for PPIX-GUN4 0.8 /min Ref. 3 & 4
Kf Rate of formation of PPIX due to protoheme regulation 20 nM/min Thorough lit search yielded no value; estimated from Willows, R. 2014 (unpublished).
χf Switch triggering PPIX production when protoheme concentration drops below 50 nM 1 (when [Protoheme] < 50 nM), 0 (when [Protoheme] > 50 nM) Thorough lit search yielded no switch value; estimated from Willows, R. 2014 (unpublished).
βdeg Degradation rate/rate at which protoheme is removed by conversion to heme 1 nM/min Thorough lit search yielded no value; since we haven’t included the formation of heme in our system, this compensates for this process; estimated from Willows, R. 2014 (unpublished).
αG Formation rate of GUN4 from gene 1 nM/min No value available due to lack of characterisation in E. coli; estimated from Willows, R. 2014 (unpublished).

Table 2. Constant values of the model components based on reported literature.

References

  • Storm, P., Tibiletti, T., Hall, M., Funk, C. Refolding and enzyme kinetic studies on the ferrochelatase of the cyanobacterium Synechocystis sp. PCC 6803. PloS one. 2013. 8.
  • Adhikari, N.D., Orler, R., Chory, J., Froehlich, J.E., Larkin, R.M. Porphyrins promote the association of GENOMES UNCOUPLED 4 and a Mg-chelatase subunit with chloroplast membranes. JBC. 2009. 284: 24783-96.
  • Reid, J.D., Hunter, C.N. Magnesium-dependent ATPase activity and cooperativity of magnesium chelatase from Synechocystis sp. PCC6803. JBC. 2004. 279: 26893-9.
  • Adams, N.B., Marklew, C.J., Brindley, A.A., Hunter, C.N., Reid, J.D. Characterization of the magnesium chelatase from Thermosynechococcus elongatus. Biochemical Journal. 2014. 457: 163-70.