Team:SydneyUni Australia/Modelling Model

From 2013.igem.org

(Difference between revisions)
(Created page with "{{Team:SydneyUni_Australia/Style}} {{Team:SydneyUni_Australia/Header}} __NOTOC__ =='''Running the Model'''== The model was run using MATLAB through the ode solver ODE45. Th...")
 
(18 intermediate revisions not shown)
Line 1: Line 1:
{{Team:SydneyUni_Australia/Style}}
{{Team:SydneyUni_Australia/Style}}
{{Team:SydneyUni_Australia/Header}}
{{Team:SydneyUni_Australia/Header}}
-
 
__NOTOC__
__NOTOC__
-
=='''Running the Model'''==
+
<html>
 +
<h1><b> Table of Contents <b></h1>
 +
<br><br>
 +
<h3>
-
The model was run using MATLAB through the ode solver ODE45.
 
 +
<a href="#scheme">1 - Schematic of Engineered Metabolic pathway</a><br><br>
-
The modelling for both the pathays had initial conditions of:
+
<a href="#ODE">2 - ODE Model</a><br><br>
-
<center>
+
-
{| class="wikitable"
+
-
!Constant !! Value !! Comment
+
-
|-
+
-
|K<sub>M A</sub> || 0.530 mM || From literature [1] 
+
-
|-
+
-
|k<sub>cat A</sub> || 3.3 s<sup>-1</sup> || From literature [1] 
+
-
|-
+
-
|K<sub>M B</sub> || 0.940 mM || From literature [2] 
+
<a href="#info">3 - General Information regarding the Enzymes Involved in the Metabolic Pathway</a><br><br>
-
|-
+
-
|k<sub>cat B</sub> || 0.0871 s<sup>-1</sup> || From literature [2] 
+
<br><br>
-
|-
+
-
|K<sub>M C</sub> || 7.2 mM || From literature [3] 
+
</h3>
-
|-
+
-
|k<sub>cat c</sub> || 89.8 s<sup>-1</sup>|| From literature [3] 
+
<a name="scheme"></a>
-
|-
+
</html>
-
|K<sub>M D</sub> || 0.160 mM || From literature [4] 
+
==A Schematic of the Engineered Metabolic Pathway:==
-
|-
+
<center>
-
 
+
<img src="https://static.igem.org/mediawiki/2013/3/3d/Pathway_RRR.png" width="900"></html>
-
|k<sub>cat D</sub> || 0.600 s<sup>-1</sup> || From literature [4] 
+
-
|-
+
-
 
+
-
|K<sub>M E</sub> || 20 mM || From literature [5] 
+
-
|-
+
-
 
+
-
|k<sub>cat E</sub> || 25.4 s<sup>-1</sup>|| From literature [5] 
+
-
|-
+
-
 
+
-
|β, γ, δ & ε|| 0 mM || Not naturally present in cells. 
+
-
|-
+
-
 
+
-
|α<sub>out</sub> || 1 mM || Initial concentration of DCA in solution (arbitrary) 
+
-
|-
+
-
 
+
-
|α<sub>in</sub> || 0.001 mM || Initial concentration of DCA in cell 
+
-
|-
+
-
 
+
-
|A, B, C, D, E || 25.55 mM || Estimation from literature [11], as described in "principles".
+
-
|-
+
-
 
+
-
|  || 2 x 10<sup>8</sup> cells / mL || Initial cell concentration which allows appropriate growth 
+
-
|-
+
-
|}
+
</center>
</center>
-
Physical conditions: the model assumes that the cells are present in minimal media prior to DCA exposure and that the DCA is instantaneously and evenly present at time = 0 min in a solution of homogeneously mixed cells
 
-
Many assumptions have been made in the construction of the model and are outlined in the 'principles' section.
 
 +
<div align="justify">
 +
The above figure is a simplified schematic of the two metabolic pathways which are considered. The symbols in squares signify the intracellular concentrations of the associated metabolite (Greek letters) or enzyme (capital English letters).  These symbols will be used throughout the analysis. </div>
-
By using the constants summarised in the previous section the flux, J, took the value (alongside the bacterial surface area, S):
 
-
[[File: my spoon is too big.png|270px]]
+
<html>
 +
<div id="ODE"></div><a name="ODE"></a>
 +
</html>
 +
==ODE Model==
 +
Shown below are the two systems of ODEs we used to model each pathway.
 +
The non-monooxygenase pathway (with dhlA (A) and adh1b1*2  (B)).
 +
[[File:Igem ode 9111.png|center]]
-
=='''The Non-monooxygenase Pathway'''==
+
The monooxygenase pathway (with p450 (C)).
-
'''ODE overview:'''
+
[[File:Igem ode 922.png|center]]
-
[[File:Igem ode 111.png|center]]
+
Generally, each line represents how a metabolic intermediate changes over time; this is dictated by the relative rate at which the intermediate is formed vs. the rate at which it is removed.  
-
The summary of the ODE (explained and justified in the previous section).
+
The two above systems of ODEs were used to model two things:
 +
 +
'''1.''' The rate at which extracellular DCA is removed from solution (given the number of cells in solution). This is considering the rate at which DCA crosses the plasma membrane (from solution to cytoplasm) of all cells in solution.
 +
 +
'''2.''' The rate at which the intracellular concentrations of the metabolic intermediates change over time within a single cell. This considers the relative activities of the enzymes of the metabolic pathway.  
-
'''Raw MATLAB code:'''
 
-
<nowiki>
+
===='''Notation'''====
-
function dy = nop450(t,y)
+
* The symbols ''A, B, C, D'' & ''E'' represent the intracellular enzyme concentrations of 1,2-Dichloroethane Dechlorinase, Alcohol Dehydrogenase, Cytochrome p450, Chloroacetaldehyde Dehydrogenase and Haloacetate Dehydrogenase respectively. They have units of mM.
-
+
-
dy=zeros(7,1);
+
-
+
-
dy(1)= -y(7)*(6*10^(-12))*0.0463067*(y(1)-y(2));
+
-
dy(2)= ((6*10^(-12))*0.0463067*(y(1)-y(2)))-3.3*25.55*(y(2)/(0.53+y(2)));
+
-
dy(3)= 3.3*25.55*(y(2)/(0.53+y(2)))-0.0871*25.55*(y(3)/(0.94+y(3)));
+
-
dy(4)= .0871*25.55*(y(3)/(0.94+y(3)))- 0.6*25.55*(y(4)/(0.16+y(4)));
+
-
dy(5)= 0.6*25.55*(y(4)/(0.16+y(4))) - 25.4*25.55*(y(5)/(20+y(5)));
+
-
+
-
  if y(6) >  2*10^(-10)
+
-
      dy(6)= 25.4*25.55*(y(5)/(20+y(5))) -1.5789*10^(-10)
+
-
  else
+
-
      dy(6) = 25.4*25.55*(y(5)/(20+y(5)))
+
-
  end
+
-
 
+
-
  if y(6) > 0.0005
+
-
        if y(7) > 1.6*10^11
+
-
            dy(7)=0
+
-
        else
+
-
            dy(7) = 5*10^6
+
-
        end
+
-
  else
+
-
        dy(7) = 0
+
-
  end
+
-
+
-
end
+
 +
* The symbols'' α<sub>in</sub>, β, γ, δ'' and ''ε'' represent the intracellular concentrations of the metabolic intermediates DCA, 2-chloroethanol, chloroacetaldehyde, chloroacetate and glycolate respectively. They have units mM. ''α<sub>out</sub>'' represents the extracellular concentration of DCA, i.e. the concentration of DCA in solution. It has units mM.
-
</nowiki>
+
* The symbol Ψ represents the concentration of cells in solution and has units cells/mL.
 +
* The function ''J(α<sub>out</sub>,α<sub>in</sub>)'' represents the rate which extracellular DCA moves across the plasma membrane and into the single cell (the flux rate). It is a function of the extra and intracellular concentrations of DCA and has a value of surface area per time m<sup>2</sup>s<sup>-1</sup>. The value ''S'' represents the average surface area of the cellular membrane of <i>E. coli</i>. By multiplying ''S'' by ''J'' one achieves the total amount of DCA flowing into a single cell per unit time. The total amount of DCA flowing from solution is achieved by multiplying Ψ by ''S'' and ''J''.
 +
<html><a name="info"></a></html>
 +
==General Information regarding the Enzymes Involved in the Metabolic Pathway ==
 +
<center>
 +
{| class="wikitable"
 +
!Enzyme !! Gene !! Symbol !! Constants !! Substrate !! Product !! Ref
 +
|-
 +
|1,2-Dichloroethane Dechlorinase || dhlA || A || K<sub>M A</sub> = 0.53 mM, k<sub>cat A</sub> = 3.3 s<sup>-1</sup> || DCA (α<sub>in</sub>)|| 2-chloroethanol (β)|| [1]
 +
|-
 +
|Alcohol Dehydrogenase || Adhlb 1/2* || B || K<sub>M B</sub> = 0.94 mM, k<sub>cat B</sub> = 0.0871 s<sup>-1</sup>  || 2-chloroethanol (β) || chloroacetaldehyde (γ) || [2]
 +
|-
 +
|p450 || p450 || C || K<sub>M C</sub> = 0.120 mM, k<sub>cat C</sub> = 0.0113 s<sup>-1</sup>  || DCA (α<sub>in</sub>)|| chloroacetaldehyde (γ) || [3]
 +
|-
 +
|Chloroacetaldehyde Dehydrogenase || aldA || D || K<sub>M D</sub> = 0.06 mM, k<sub>cat D</sub> = 0.60 s<sup>-1</sup>  || chloroacetaldehyde (γ)|| chloroacetate (δ) || [4]
 +
|-
 +
|Haloacetate Dehydrogenase || dhlB || E || K<sub>M E</sub> = 20 mM, k<sub>cat E</sub> = 25.4 s<sup>-1</sup>  || chloroacetate (δ) || glycolate (ε)|| [5]
 +
|}
 +
</center>
-
===='''MATLAB output:'''====
 
-
 
-
 
-
'''The rate at which DCA is removed in solution:'''
 
-
 
-
[[File: DCAcominin.png|950px]]
 
-
 
-
Graph 1: The blue line represents how the concentration of DCA in solution (extracellular) decreases due to the action of our DCA degraders. The red line represents the intracellular concentions of the metabolites (disregard in this graph). One can see that an initial concentration of 2E8 cells/mL completely removes the DCA with a concentration 1mM in (roughly) 150mins.
 
-
 
-
 
-
 
-
 
-
'''The Rate at which the intracellular concentration of the metabolites change over time:'''
 
-
 
-
[[File: intermediates1111.png|950px]]
 
-
 
-
Graph 2: Each line represents the concentration of each of the metabolites . This graph is simply a rescaling of graph 1. Note: the glycolate won’t accumulate in the cell – it is metabolised – the model had glycolate as the final product. It is used to show that the presence of glycolate can attribute to cell growth.
 
-
 
-
 
-
 
-
 
-
 
-
'''Rescaling the graph once again: the rate at which the metabolic intermediate change over time.'''
 
-
 
-
[[File: intermediates2222.png|950px]]
 
-
 
-
Graph 3: Each line represents the concentration of each of the metabolites . This graph is simply a rescaling of graph 1 and 2.
 
-
 
-
 
-
 
-
 
-
'''The rate at which the cells grow over time:'''
 
-
 
-
[[File: cellscellscells.png|950px]]
 
-
 
-
Graph 4: the blue line represents the linear increase of cells due to the presence of intracellular glycolate.
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
=='''The Non-monooxygenase Pathway'''==
 
-
 
-
 
-
 
-
'''ODE overview''':
 
-
 
-
[[File: Igem_ode_22.png|450px]]
 
-
 
-
 
-
 
-
'''Raw MATLAB code:'''
 
-
 
-
<nowiki>
 
-
function dy = p450(t,y)
 
-
 
-
dy=zeros(6,1);
 
-
 
-
dy(1)= -y(6)*(6*10^(-12))*0.0463067*(y(1)-y(2));
 
-
dy(2)= ((6*10^(-12))*0.0463067*(y(1)-y(2))) - 0.0113*25.55*(y(2)/(.12+y(2)));
 
-
dy(3) = 0.0113*25.55*(y(2)/(.12+y(2))) - 0.6*25.55*(y(3)/(0.16+y(3)));
 
-
dy(4)= 0.6*25.55*(y(3)/(0.16+y(3))) - 25.4*25.55*(y(4)/(20+y(4)));
 
-
 
-
 
-
  if y(5) >  2*10^(-10)
 
-
      dy(5)= 25.4*25.55*(y(4)/(20+y(4))) -1.5789*10^(-10)
 
-
  else
 
-
      dy(5) = 25.4*25.55*(y(4)/(20+y(4)))
 
-
  end
 
-
 
 
-
  if y(5) > 0.0005
 
-
        if y(6) > 1.6*10^11
 
-
            dy(6)=0
 
-
        else
 
-
            dy(6) = 5*10^6
 
-
        end
 
-
  else
 
-
        dy(6) = 0
 
-
  end
 
-
 
-
end
 
-
 
-
</nowiki>
 
-
 
-
 
-
 
-
'''MATLAB output:'''
 
-
 
-
 
-
 
-
'''The rate at which DCA is removed in solution:'''
 
-
 
-
[[File: DCAcominin1.png|950px]]
 
-
 
-
Graph 5: The blue line represents how the concentration of DCA in solution (extracellular) decreases due to the action of our DCA degraders.
 
-
 
-
 
-
 
-
 
-
 
-
'''The Rate at which the intracellular concentration of the metabolites change over time:'''
 
-
 
-
[[File: intermediates111.png|950px]]
 
-
 
-
Graph 6: Each line represents the intracellular concentration of each of the metabolites over time within a single cell. The colour associated with each metabolite is depicted in the figure legend. This graph is simply a rescaling of graph 5.
 
-
 
-
 
-
 
-
 
-
 
-
'''Rescaling the graph once again: the rate at which the metabolic intermediate change over time.'''
 
-
 
-
[[File: thismessagewasbroughttoyoubuygoontheproudsponserofthismodel.png|950px]]
 
-
 
-
Graph 7: Again, each line represents the intracellular concentration of each of the metabolites over time within a single cell. This graph is simply a rescaling of graph 5 and 6.
 
-
 
-
 
-
 
-
 
-
 
-
'''The rate at which the cells grow over time:'''
 
-
 
-
[[File: cellscellscelllls.png|950px]]
 
-
 
-
Graph 8: The blue line represents the linear increase of the total number of cells due to the presence of intracellular glycolate.
 
-
 
-
 
-
 
-
 
-
 
-
'''Rescaling Shows a lag effect on cell growth:'''
 
-
 
-
[[File: cellscellscelllllllllllls.png|950px]]
 
-
 
-
Graph 9: The blue line represents the increase of cells due to the presence of intracellular glycolate. This graph is a rescaling of graph 8 to show the inital lag in cellualr growth.
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
== '''Conclusions:''' ==
 
-
 
-
From Graphs 1 and 5, one can see that 1mM of DCA is removed from solution within roughly 50 minutes when the DCA degrading cells are at a concentration of 2E8 cells/mL.
 
-
 
-
From Graphs 4, 8 and 9, it is evident that bacterial growth occurs. This growth is due to the production of glycolate, and by comparing graphs 6 and 9, one can see that bacterial growth correlates with glycolate accumulation.
 
-
 
-
The cytotoxic metabolic intermediate chloroactealdehyde doesn't accumulate to a significant concentration in any of the pathways and is consistently at a negligibly small concentration. From Graphs 3 and 7 one can see that chloroacetaldehyde reaches a maximum concentration of roughly 0.2 mM in both pathways. Chloroacetaldehyde is seen to be metabolised very quickly; this concentration maximum is very short lived where it peaks at roughly 0.03 seconds and returns back to 0 mM by 0.5 seconds. It is expected that chloroacetaldehyde toxicity will not be a problem in our engineered cells.
 
-
 
-
It is also possible to conclude that the pathways remove DCA at the same rate (through comparing graphs 1 and 5).
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
 
-
== '''References:''' ==
 
-
 
-
[1] Krooshof, G. H., Ridder, I. S., Tepper, A. W., Vos, G. J., Rozeboom, H. J., Kalk, K. H., Dijkstra, B. W. & Janssen, D. B. (1998). Kinetic Analysis and X-ray Structure of Haloalkane Dehalogenase with a Modified Halide-Binding Site. <i>Biochemistry</i>, <b>37</b>(43), 15013-15023.
 
-
 
-
[2]  Janecki, D. J., Bemis, K. G., Tegeler, T. J., Sanghani, P. C., Zhai, L., Hurley, T. D., Bosron, W. F. & Wang, M. (2007). A multiple reaction monitoring method for absolute quantification of the human liver alcohol dehydrogenase ADH1C1 isoenzyme. <i>Analytical Biochemistry</i>, <b>369</b>(1), 18-26.
 
-
 
-
[3] Pandey, A. V. & Flück, C. E. (2013). NADPH P450 oxidoreductase: Structure, function, and pathology of diseases. <i>Pharmacology & Therapeutics</i>, <b>138</b>(2), 229-254.
 
-
 
-
[4]  van der Ploeg, J., Shmidt, M. P., Landa, A. S., & Janssen, D. B. (1994). Identification of Chloroacetaldehyde Dehydrogenase Involved in 1,2-Dichloroethane Degradation. <i>Applied Environmental Microbiology</i> <b>60</b>(5), 1599-1605.
 
-
 
-
[5]  van der Ploeg, J., van Hall, G. & Janssen, D. B. (1991) Characterization of the haloacid dehalogenase from Xanthobacter autotrophicus GJ10 and sequencing of the dhlB gene. <i>Journal of Bacteriology</i>, <b>173</b>(24), 7925-33.
 
-
 
-
[6] Sinensky, M. I. (1974). Homeoviscous Adaption – A Homeostatic Process that Regulates the Viscosity of Membrane Lipids in <i>Escherichia coli</i>. <i>Proceedings from the National Academy of Science of the United States of America</i>, <b>71</b>(2), 522-525.
 
-
 
-
[7] CyberCell Database
 
-
 
-
[8]
 
-
 
-
[9]
 
-
 
-
[10] http://www.dtsc.ca.gov/AssessingRisk/Upload/12dca.pdf
 
-
[11] Ishihama, Y., Schmidt, T., Rappsilber, J., Mann, M., Hartl F. U., Kerner, M. J. & Frishman, D. (2008) Protein abundance profiling of the Escherichia coli cytosol. BMC Genomics, <b>9</b>:102
 
-
[12] Lord, J. M. (1972) Glycolate oxidoreductase in Escherichia coli. Biochemica et Biophysica Acta <b>267</b>:2, 227-327.
 
{{Team:SydneyUni_Australia/Footer}}
{{Team:SydneyUni_Australia/Footer}}

Latest revision as of 01:54, 29 October 2013

SydneyUniversity Top Banner.jpg SydneyUniversity Bottom Banner.jpg


Table of Contents



1 - Schematic of Engineered Metabolic pathway

2 - ODE Model

3 - General Information regarding the Enzymes Involved in the Metabolic Pathway



A Schematic of the Engineered Metabolic Pathway:

<img src="Pathway_RRR.png" width="900"></html>


The above figure is a simplified schematic of the two metabolic pathways which are considered. The symbols in squares signify the intracellular concentrations of the associated metabolite (Greek letters) or enzyme (capital English letters). These symbols will be used throughout the analysis.


ODE Model

Shown below are the two systems of ODEs we used to model each pathway.


The non-monooxygenase pathway (with dhlA (A) and adh1b1*2 (B)).

Igem ode 9111.png


The monooxygenase pathway (with p450 (C)).

Igem ode 922.png

Generally, each line represents how a metabolic intermediate changes over time; this is dictated by the relative rate at which the intermediate is formed vs. the rate at which it is removed.

The two above systems of ODEs were used to model two things:

1. The rate at which extracellular DCA is removed from solution (given the number of cells in solution). This is considering the rate at which DCA crosses the plasma membrane (from solution to cytoplasm) of all cells in solution.

2. The rate at which the intracellular concentrations of the metabolic intermediates change over time within a single cell. This considers the relative activities of the enzymes of the metabolic pathway.


Notation

  • The symbols A, B, C, D & E represent the intracellular enzyme concentrations of 1,2-Dichloroethane Dechlorinase, Alcohol Dehydrogenase, Cytochrome p450, Chloroacetaldehyde Dehydrogenase and Haloacetate Dehydrogenase respectively. They have units of mM.
  • The symbols αin, β, γ, δ and ε represent the intracellular concentrations of the metabolic intermediates DCA, 2-chloroethanol, chloroacetaldehyde, chloroacetate and glycolate respectively. They have units mM. αout represents the extracellular concentration of DCA, i.e. the concentration of DCA in solution. It has units mM.
  • The symbol Ψ represents the concentration of cells in solution and has units cells/mL.
  • The function J(αoutin) represents the rate which extracellular DCA moves across the plasma membrane and into the single cell (the flux rate). It is a function of the extra and intracellular concentrations of DCA and has a value of surface area per time m2s-1. The value S represents the average surface area of the cellular membrane of E. coli. By multiplying S by J one achieves the total amount of DCA flowing into a single cell per unit time. The total amount of DCA flowing from solution is achieved by multiplying Ψ by S and J.


General Information regarding the Enzymes Involved in the Metabolic Pathway

Enzyme Gene Symbol Constants Substrate Product Ref
1,2-Dichloroethane Dechlorinase dhlA A KM A = 0.53 mM, kcat A = 3.3 s-1 DCA (αin) 2-chloroethanol (β) [1]
Alcohol Dehydrogenase Adhlb 1/2* B KM B = 0.94 mM, kcat B = 0.0871 s-1 2-chloroethanol (β) chloroacetaldehyde (γ) [2]
p450 p450 C KM C = 0.120 mM, kcat C = 0.0113 s-1 DCA (αin) chloroacetaldehyde (γ) [3]
Chloroacetaldehyde Dehydrogenase aldA D KM D = 0.06 mM, kcat D = 0.60 s-1 chloroacetaldehyde (γ) chloroacetate (δ) [4]
Haloacetate Dehydrogenase dhlB E KM E = 20 mM, kcat E = 25.4 s-1 chloroacetate (δ) glycolate (ε) [5]



With thanks to: