Team:Edinburgh/Modelling/Kappa/Analysis

From 2012.igem.org

(Difference between revisions)
m
 
(8 intermediate revisions not shown)
Line 28: Line 28:
border:1px solid #ccc;
border:1px solid #ccc;
}
}
-
#fig2, #fig3,
+
#fig1, #fig2,
-
#fig4, #fig5,
+
#fig3, #fig4,
-
#fig6, #fig7,
+
#fig5, #fig6,
-
#fig8, #fig9{
+
#fig7, #fig8{
width:100%;
width:100%;
}
}
Line 44: Line 44:
<div class="text">
<div class="text">
<p class="h1">
<p class="h1">
-
Analysis of the Model
+
Analysis of the Sub-models
-
</p>
+
-
<p class="normal-text">
+
-
<br /><br />
+
-
<!--img src=""><br />
+
-
Fig 1. Contact map of the model of the whole electron transfer made with <a href="http://rulebender.cs.pitt.edu/wordpress/">RuleBender-1.1.415</a-->.
+
</p>
</p>
<p class="h2">
<p class="h2">
Line 57: Line 52:
The file 1_TCA.ka models the electron transfer between the TCA cycle and the quinol QH<sub>2</sub>. This process starts with the absorption of glucose by the TCA cycle and passes through Complex I and Complex II before it reaches the quinol.
The file 1_TCA.ka models the electron transfer between the TCA cycle and the quinol QH<sub>2</sub>. This process starts with the absorption of glucose by the TCA cycle and passes through Complex I and Complex II before it reaches the quinol.
<br /><br />
<br /><br />
-
<img id="fig8" src="https://static.igem.org/mediawiki/2012/6/67/Kappa-fig07.JPG"><br />
+
<img id="fig1" src="https://static.igem.org/mediawiki/2012/6/67/Kappa-fig07.JPG"><br />
-
The above graph shows the amounts of glucose, Complex I and Complex II as the electron transfer between the TCA cycle and the quinol goes on.
+
<b>Figure 1.</b> This graph shows the amounts of Complex I and Complex II as the electron transfer between the TCA cycle and the quinol goes on.
<br /><br />
<br /><br />
-
<img id="fig9" src="https://static.igem.org/mediawiki/2012/b/ba/Kappa-fig08.JPG">
+
The amount of Complex I with NADH together with the amount of Complex I without NADH gives us the total amount of Complex I. This is why the two corresponding graphs shown above are “complementary”.
 +
<br />
 +
The same can be said about the amount of Complex II with succinate and the amount of Complex II without succinate.
<br /><br />
<br /><br />
-
This is an idealized plot of the TCA cycle given optimum conditions. The actual model doesn’t work as cleanly as this, because adding the realistic external constraint seriously affects all involved parameters and also injects a great deal of randomness. Reality generally has a nasty habit of ruining clear and concise approaches. However we include this because it presents a good, if idealised, overview of the working of the TCA cycle.
+
When working on the TCA cycle, we started with an “ideal” model. What we mean is that there was no restriction on the total amount of agents or tokens of any kind. However, this is never the case in reality where the electron transfer process takes place in the limits of a cell.
 +
<br /><br />
 +
Even so, it is useful to consider such an ideal model in order to study the interaction between the different elements of the system when there are no restrictions.
 +
<br /><br />
 +
<img id="fig2" src="https://static.igem.org/mediawiki/2012/b/ba/Kappa-fig08.JPG">
 +
<br />
 +
<b>Figure 2.</b> Idealised electron transfer system
 +
<br /><br />
 +
Figure 2 presents such an idealised model. In it we can again see the “complimentary” property of the tuples (Complex I with NADH, Complex I without NADH) and (Complex II with succinate, Complex II without succinate).
</p>
</p>
<p class="h2">
<p class="h2">
Line 71: Line 76:
</p>
</p>
<p class="normal-text">
<p class="normal-text">
-
In this file the passage of the electrons through the periplasm is modelled.
+
In this file we model the passage of the electrons through the periplasm.
<br /><br />
<br /><br />
-
First the passage through the inner membrane occurs with the reaction between NapC and mtrA, then mtrA moves through the periplasm itself, eventually reaching an mtrB attached to the outer membrane to which it transfers the electrons it’s carrying. Additionally, when moving through the periplasm, mtrA can encounter soluble Iron molecules which will also accept electrons.
+
First, the electron transfer through the inner membrane occurs with the reaction between NapC and mtrA. Then mtrA moves through the periplasm, eventually reaching an mtrB attached to the outer membrane to which it transfers the electrons it is carrying. Additionally, when moving through the periplasm, mtrA can encounter soluble Iron molecules which will also accept electrons.In this file the passage of the electrons through the periplasm is modelled.
<br /><br />
<br /><br />
-
<img id="fig2" src="https://static.igem.org/mediawiki/2012/b/b4/Kappa-fig02.JPG">
+
<img id="fig3" src="https://static.igem.org/mediawiki/2012/b/b4/Kappa-fig02.JPG">
 +
<br />
 +
<b>Figure 3.</b>
<br /><br />
<br /><br />
-
Here we can see that, "mtrCAB_electrons_from_SolubleIron" is significantly larger than "mtrCAB_electrons_from_mtrB". As it turns out, the bulk of the produced electric current is via soluble Iron present in the periplasm. Noticing this we made a second model where MtrBC is absent and therefore we can have more mtrA agents. The "mtrA_electrons_from_SolubleIron" is even larger, as expected, but though the number of mtrAs has almost doubled(from 2100 to 4000), the growth in electron rate is ~1.4 It’s clear that the extra mtrAs suffer from diminishing returns.
+
As shown in Figure 3, "mtrCAB_electrons_from_SolubleIron" is significantly larger than "mtrCAB_electrons_from_mtrB". As it turns out, the bulk of the produced electric current is via soluble Iron present in the periplasm.
 +
<br /><br />
 +
Noticing this we made a second model where MtrBC is absent and therefore we can have more mtrA agents. The "mtrA_electrons_from_SolubleIron" is even larger, as expected, but though the number of mtrAs has almost doubled (from 2100 to 4000), the growth in electron rate is ~1.4 It's clear that the extra mtrAs suffer from diminishing returns.
</p>
</p>
<p class="h3">
<p class="h3">
Line 83: Line 92:
</p>
</p>
<p class="normal-text">
<p class="normal-text">
-
This file models the transfer of iron from the outer membrane of the cell to it’s final acceptor, insoluble iron. mtrC is always connected to a mtrB, so for the purpose of this model we consider them a single agent.
+
This file models the transfer of iron from the outer membrane of the cell to it's final acceptor, insoluble iron. mtrC is always connected to a mtrB, so for the purpose of this model we consider them a single agent.
<br /><br />
<br /><br />
-
If a particle of insoluble iron comes close enough to the mtrC then the electron can be directly transferred. Intercellular Flavin can also pick up the electrons from mtrC and then separately connect to insoluble iron particles. Since the number of mtrCs is quite small the Flavin transfer mechanism is important.
+
If a particle of insoluble iron comes close enough to the mtrC then the electron can be directly transferred. Inter-cellular Flavin can also pick up the electrons from mtrC and then separately connect to insoluble iron particles. Since the number of mtrCs is quite small the Flavin transfer mechanism is important.
<br /><br />
<br /><br />
-
<img id="fig3" src="https://static.igem.org/mediawiki/2012/1/1b/Kappa-fig03.JPG">
+
<img id="fig4" src="https://static.igem.org/mediawiki/2012/1/1b/Kappa-fig03.JPG">
 +
<br />
 +
<b>Figure 4.</b>
<br /><br />
<br /><br />
-
Here we can see the amount of electrons transferred directly to insoluble iron vs. the amount transferred via Flavin. We also test what happens when the insoluble iron particle count is significantly increased by using smaller particles. Specifically, we look at micrometer sized particles and nanometer sized particles. As expected, the rates are higher with the nanoParticles, but considering there are more than 1000x more particles of nanometer size the difference isn’t as dramatic as one might assume. The small number of mtrC agents is likely the limiting step here.
+
Figure 4 shows the amount of electrons transferred directly to insoluble iron vs. the amount transferred via Flavin.
-
<img id="fig4" src="https://static.igem.org/mediawiki/2012/6/61/Kappa-fig04.JPG">
+
<br /><br />
<br /><br />
-
Here we can see that similar to the number of insoluble iron particles, the number of Flavins can also have a big effect on the electron transfer rate. We show the direct mtrC-InsolubleIron rate as well as flavin-InsolubleIron rates for various different quantities of added Flavin.
+
We also test what happens when the insoluble iron particle count is significantly increased by using smaller particles. Specifically, we look at micrometer sized particles and nanometer sized particles. As expected, the rates are higher with the nanoParticles, but considering there are more than 1000x more particles of nanometer size the difference isn't as dramatic as one might assume. The small number of mtrC agents is likely the limiting step here.
 +
<br /><br />
 +
<img id="fig5" src="https://static.igem.org/mediawiki/2012/6/61/Kappa-fig04.JPG">
 +
<br />
 +
<b>Figure  5.</b>
 +
<br /><br />
 +
The number of Flavins can also have a big effect on the electron transfer rate, as seen in figure 5. We show the direct mtrC-InsolubleIron rate as well as flavin-InsolubleIron rates for various different quantities of added Flavin.
</p>
</p>
<p class="h3">
<p class="h3">
Line 102: Line 118:
We look at the model both with the previous values and with 10X more NapC, which brings its behaviour more in line with that expected from <i>Shewanella</i>.  
We look at the model both with the previous values and with 10X more NapC, which brings its behaviour more in line with that expected from <i>Shewanella</i>.  
<br /><br />
<br /><br />
-
<img id="fig5" src="https://static.igem.org/mediawiki/2012/c/c9/Kappa-fig05.JPG">
+
<img id="fig6" src="https://static.igem.org/mediawiki/2012/c/c9/Kappa-fig05.JPG">
 +
<br />
 +
<b>Figure 6.</b>
<br /><br />
<br /><br />
Here we see the familiar relationship between amount of electrons transferred directly through mtrC and those transferred through Flavin. We can also see that increasing the NapC to 10X the normal quantity scales up all rates by roughly the same amount.
Here we see the familiar relationship between amount of electrons transferred directly through mtrC and those transferred through Flavin. We can also see that increasing the NapC to 10X the normal quantity scales up all rates by roughly the same amount.
-
<img id="fig6" src="https://static.igem.org/mediawiki/2012/6/6c/Kappa-fig06.JPG">
+
<img id="fig7" src="https://static.igem.org/mediawiki/2012/6/6c/Kappa-fig06.JPG">
 +
<br />
 +
<b>Figure 7.</b>
<br /><br />
<br /><br />
-
The effect of the added NapC is more pronounced in the periplasm. Here originally the rate between electrons given to soluble iron and those given to mtrB is ~9, but with 10X NapC it becomes almost ~20, which is what would be expected from Shewanella. </p>
+
The effect of the added NapC is more pronounced in the periplasm. Here originally the rate between electrons given to soluble iron and those given to mtrB is ~9, but with 10X NapC it becomes almost ~20, which is what would be expected from <i>Shewanella</i>.
-
<p class="h2">
+
-
Conclusion and Future work
+
-
</p>
+
-
<p class="normal-text">
+
-
This model is a close approximation of the underlying biological system. It incorporates a lot of experimental work and has been verified against yet more research results, showing very promising results.
+
<br /><br />
<br /><br />
-
The next step is to use this model to try to better understand the biological underpinnings of the system. Specifically, we plan to do a series of sensitivity analysis type tests and see how varying the different parameters affects the results. In this fashion we can gain insight into which biological subsystems offer the most promise, thus hopefully reducing the work that needs to be undertaken in the lab.
 
<br /><br />
<br /><br />
-
Furthermore additions to the model trying to boost the electron output or perhaps replicate some other complex behaviour could first be tried in the model. This again will work as a preliminary step which will support or cast doubt upon the biological validity of the proposed addition and will, at the very least, raise some interesting questions about the working of any new proposed system.
+
<span class="intense-emphasis">So far we saw how isolated sub-systems work on their own. Next why consider the electron transfer system as a whole.</span>
 +
<br /><br />
</p>
</p>
<p style="border-top:1px solid #ccc;" class="normal-text">
<p style="border-top:1px solid #ccc;" class="normal-text">
-
<span class="intense-amphasis">* All graphs are aggregated over 50 runs of the model. They show the mean and the standard deviation of the curve. Many thanks to Donal Stewart for providing the tool which made aggregating all the hundreds of text files painless and actually fun!</span>
+
<span class="subtle-emphasis">* All graphs are aggregated over 50 runs of the model. They show the mean and the standard deviation of the curve. Many thanks to Donal Stewart for providing the tool which made aggregating all the hundreds of text files painless and actually fun!</span>
 +
</p>
 +
<p class="normal-text" style="text-align:center">
 +
<br /><br />
 +
<a href="https://2012.igem.org/Team:Edinburgh/Modelling/Kappa/Structure"><span class="intense-emphasis">&lt;&lt;Prev</span></a><span style="color:white;">__</span>3/4</span><span style="color:white;">__</span><a href="https://2012.igem.org/Team:Edinburgh/Modelling/Kappa/Analysis-All"><span class="intense-emphasis">Next&gt;&gt;</span></a>
 +
<br /><br />
</p>
</p>
</div><!-- /text -->
</div><!-- /text -->

Latest revision as of 20:03, 26 October 2012

Analysis of the Sub-models

The TCA Cycle

The file 1_TCA.ka models the electron transfer between the TCA cycle and the quinol QH2. This process starts with the absorption of glucose by the TCA cycle and passes through Complex I and Complex II before it reaches the quinol.


Figure 1. This graph shows the amounts of Complex I and Complex II as the electron transfer between the TCA cycle and the quinol goes on.

The amount of Complex I with NADH together with the amount of Complex I without NADH gives us the total amount of Complex I. This is why the two corresponding graphs shown above are “complementary”.
The same can be said about the amount of Complex II with succinate and the amount of Complex II without succinate.

When working on the TCA cycle, we started with an “ideal” model. What we mean is that there was no restriction on the total amount of agents or tokens of any kind. However, this is never the case in reality where the electron transfer process takes place in the limits of a cell.

Even so, it is useful to consider such an ideal model in order to study the interaction between the different elements of the system when there are no restrictions.


Figure 2. Idealised electron transfer system

Figure 2 presents such an idealised model. In it we can again see the “complimentary” property of the tuples (Complex I with NADH, Complex I without NADH) and (Complex II with succinate, Complex II without succinate).

The Electron Transport Chain

File 3_MtrABC.ka:

In this file we model the passage of the electrons through the periplasm.

First, the electron transfer through the inner membrane occurs with the reaction between NapC and mtrA. Then mtrA moves through the periplasm, eventually reaching an mtrB attached to the outer membrane to which it transfers the electrons it is carrying. Additionally, when moving through the periplasm, mtrA can encounter soluble Iron molecules which will also accept electrons.In this file the passage of the electrons through the periplasm is modelled.


Figure 3.

As shown in Figure 3, "mtrCAB_electrons_from_SolubleIron" is significantly larger than "mtrCAB_electrons_from_mtrB". As it turns out, the bulk of the produced electric current is via soluble Iron present in the periplasm.

Noticing this we made a second model where MtrBC is absent and therefore we can have more mtrA agents. The "mtrA_electrons_from_SolubleIron" is even larger, as expected, but though the number of mtrAs has almost doubled (from 2100 to 4000), the growth in electron rate is ~1.4 It's clear that the extra mtrAs suffer from diminishing returns.

File 4_UFe.ka:

This file models the transfer of iron from the outer membrane of the cell to it's final acceptor, insoluble iron. mtrC is always connected to a mtrB, so for the purpose of this model we consider them a single agent.

If a particle of insoluble iron comes close enough to the mtrC then the electron can be directly transferred. Inter-cellular Flavin can also pick up the electrons from mtrC and then separately connect to insoluble iron particles. Since the number of mtrCs is quite small the Flavin transfer mechanism is important.


Figure 4.

Figure 4 shows the amount of electrons transferred directly to insoluble iron vs. the amount transferred via Flavin.

We also test what happens when the insoluble iron particle count is significantly increased by using smaller particles. Specifically, we look at micrometer sized particles and nanometer sized particles. As expected, the rates are higher with the nanoParticles, but considering there are more than 1000x more particles of nanometer size the difference isn't as dramatic as one might assume. The small number of mtrC agents is likely the limiting step here.


Figure 5.

The number of Flavins can also have a big effect on the electron transfer rate, as seen in figure 5. We show the direct mtrC-InsolubleIron rate as well as flavin-InsolubleIron rates for various different quantities of added Flavin.

File 5.ka

This model is a concatenation of models 3 and 4 taking care that the mtrA-mtrBC reaction is accurate and also rewriting some other aspects so that it will be compatible with models 1 & 2.

We look at the model both with the previous values and with 10X more NapC, which brings its behaviour more in line with that expected from Shewanella.


Figure 6.

Here we see the familiar relationship between amount of electrons transferred directly through mtrC and those transferred through Flavin. We can also see that increasing the NapC to 10X the normal quantity scales up all rates by roughly the same amount.
Figure 7.

The effect of the added NapC is more pronounced in the periplasm. Here originally the rate between electrons given to soluble iron and those given to mtrB is ~9, but with 10X NapC it becomes almost ~20, which is what would be expected from Shewanella.



So far we saw how isolated sub-systems work on their own. Next why consider the electron transfer system as a whole.

* All graphs are aggregated over 50 runs of the model. They show the mean and the standard deviation of the curve. Many thanks to Donal Stewart for providing the tool which made aggregating all the hundreds of text files painless and actually fun!



<<Prev__3/4__Next>>