Team:Imperial College London/Modelling/Output
From 2010.igem.org
(Adding the buttons) |
(Output amplification with major math problems) |
||
Line 23: | Line 23: | ||
#wrapper { | #wrapper { | ||
width: 800px; | width: 800px; | ||
- | margin-left: | + | margin-left: 20px; |
margin-right: auto; | margin-right: auto; | ||
} | } | ||
Line 65: | Line 65: | ||
<div class="accordionButton">Abandoned Initial Attempts</div> | <div class="accordionButton">Abandoned Initial Attempts</div> | ||
- | <div class="accordionContent"></div> | + | <div class="accordionContent"> |
+ | </html> | ||
+ | =Model based on Michaelis Menten Kinetics (Weeks 4 and 5)= | ||
+ | ==Motivation== | ||
+ | We came up with a simple concept of output amplification, which is enhanced by using enzymes. It is beneficial for us to model the behaviour of our design so that we will be able to answer the following questions. | ||
+ | |||
+ | # How beneficial is the use of amplification? (Compare speed of response of transcription (and translation) with 1- or 2-step amplification) | ||
+ | # How many amplification steps are beneficial to have? Will further adding of amplification steps introduce too many time delays? | ||
+ | # Is it better to use TEV all or HIV1? | ||
+ | |||
+ | Modelling should allows us to make a decision on which design is the most efficient one. | ||
+ | |||
+ | ==First Model== | ||
+ | ===HIV1=== | ||
+ | {| class="wikitable" style="text-align: center; width: 100%; height: 170px;" | ||
+ | | [[Image:Slide2.JPG|300px|thumb|center|alt=A|At each stage of amplification a distinct protease is being used ]] | ||
+ | | align="left"| | ||
+ | '''Equations''' | ||
+ | |||
+ | *<math>\dot{m}=k_{to} - d_{to}m\,\!</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_h} = k_hm - d_hp_h</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_t} = k_tp_h - d_tp_t</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_g} = k_gp_t - d_gp_g</math> | ||
+ | | align="left"| | ||
+ | '''Parameters''' | ||
+ | *<math>k_{to}\mbox{...transcription rate of HIV1}</math> | ||
+ | *<math>d_{to}\mbox{...degradation rate of mRNA coding for HIV1}</math> | ||
+ | *<math>k_h\mbox{...translation rate of HIV1}</math> | ||
+ | *<math>d_h\mbox{...degradation rate of HIV1}</math> | ||
+ | *<math>k_t\mbox{...production rate of TEV by HIV1}</math> | ||
+ | *<math>d_t\mbox{...degradation rate of TEV}</math> | ||
+ | *<math>k_g\mbox{...production rate of GFP by TEB}</math> | ||
+ | *<math>d_g\mbox{...degradation rate of GFP}</math> | ||
+ | |} | ||
+ | |||
+ | ===TEV=== | ||
+ | {| class="wikitable" style="text-align: center; width: 100%; height: 170px;" | ||
+ | |[[Image:TEV.jpg|300px|thumb|center|alt=A|TEV is used at both stages of amplification]] | ||
+ | |align="left"| | ||
+ | '''Equations''' | ||
+ | |||
+ | |||
+ | *<math>\dot{m} = k_{to} - d_{to}m</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_t} = k_tm - d_tp_t</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_{ts}} = k_{ts}p_t - d_{ts}p_{ts}</math> | ||
+ | |||
+ | |||
+ | *<math>\dot{p_g} = k_{g1}p_t + k_{g2}p_{ts} - d_gp_g</math> | ||
+ | |align="left"| | ||
+ | '''Parameters''' | ||
+ | *<math>k_{to}\mbox{...rate of transcription by TEV}</math> | ||
+ | *<math>d_{to}\mbox{...degradation rate of mRNA coding for TEV}</math> | ||
+ | *<math>k_t\mbox{...rate of translation of TEV}</math> | ||
+ | *<math>d_t\mbox{...degradation rate of TEV}</math> | ||
+ | *<math>k_{ts}\mbox{...rate of production (fusion) of split TEV}</math> | ||
+ | *<math>d_{ts}\mbox{...degradation rate of split TEV}</math> | ||
+ | *<math>k_{g1}\mbox{...rate of production of GFP by full TEV}</math> | ||
+ | *<math>k_{g2}\mbox{...rate of production of GFP by split TEV}</math> | ||
+ | *<math>d_g\mbox{...degradation rate of GFP}</math> | ||
+ | |} | ||
+ | |||
+ | ==Improved Model which accounts for enzyme reactions (28/07/2010)== | ||
+ | ===TEV=== | ||
+ | {| class="wikitable" style="text-align: center; width: 100%; height: 170px;" | ||
+ | |[[Image:TEV.jpg|500px|thumb|center|alt=A|TEV is used at both stages of amplification]] | ||
+ | |align="left"| | ||
+ | '''Equations''' | ||
+ | * 1. Production of TEV from transcription | ||
+ | <math>\dot{p_t} = s_t - d_tp_t</math> | ||
+ | |||
+ | <math>s_t = \dfrac{k_tk_{to}}{d_{to}}</math> | ||
+ | |||
+ | |||
+ | * 2. Production of split TEV from transcription | ||
+ | <math>\dot{p_{st}} = s_{st} - d_{st}p_{st}</math> | ||
+ | |||
+ | |||
+ | * 3. Production of split GFP from transcription | ||
+ | <math>\dot{p_{sg}} = s_{sg} - d_{sg}p_{sg}</math> | ||
+ | |||
+ | |||
+ | * 4. Production of fused split TEV catalysed by TEV (1) | ||
+ | <math>\dot{p_{ts}} = \dfrac{V_{max,t}[p_{st}]}{K_{m,ts} + [p_{st}]} - d_{ts}p_{ts}</math> | ||
+ | |||
+ | |||
+ | * 5. Production of GFP catalysed by TEV (1) and fused split TEV (4) | ||
+ | <math>\dot{p_g} = \dfrac{V_{max,tg}[p_{sg}]}{K_{m,tg} + [p_{sg}]} + \dfrac{V_{max,tsg}[p_{sg}]}{K_{m,tsg} + [p_{sg}]} - d_gp_g</math> | ||
+ | |} | ||
+ | |||
+ | ===Implementation in Matlab=== | ||
+ | The Matlab code for the different stages of amplification and diagrams can be found [http://www.openwetware.org/wiki/Image:Modelling.docx here]. | ||
+ | |||
+ | ===Kinetic constants=== | ||
+ | |||
+ | {| border="1" | ||
+ | ! | ||
+ | ! GFP | ||
+ | ! TEV | ||
+ | ! split TEV | ||
+ | ! split GFP | ||
+ | |- | ||
+ | |'''<math>Km</math> and <math>k_{cat}</math>''' | ||
+ | | - | ||
+ | |<math>K_m = 0.061</math>; <math>k_{cat} = 0.16</math>; [http://peds.oxfordjournals.org/cgi/reprint/14/12/993] | ||
+ | | 40% of value for TEV | ||
+ | | - | ||
+ | |- | ||
+ | | '''Half-life or degradation rate''' | ||
+ | | Half-life in B.sub approximately 1.5 hours | ||
+ | | ? | ||
+ | | ? | ||
+ | | Half-life shorter than GFP | ||
+ | |- | ||
+ | | '''Production rate in B.sub''' | ||
+ | | ? | ||
+ | | ? | ||
+ | | ? | ||
+ | | ? | ||
+ | |} | ||
+ | |||
+ | ==Conclusion== | ||
+ | |||
+ | We were not able to obtain all the necessary constants. Hence, we decided to make educated guesses about possible relative values between the constants as well as varying them and observing the change in output. | ||
+ | |||
+ | As the result, we concluded that the amplification happens at each amplification level proposed. The magnitude of amplification varies depending on the constants. There is not much difference between using TEV or HIV1. | ||
+ | |||
+ | ==References== | ||
+ | #Kapust R. et al (2001) Tobacco etch virus protease: mechanism of autolysis and rational design of stable mutants with wild-type catalytic proficiency. Protein Engineering. [Online] 14(12), 993-1000. Available from: http://peds.oxfordjournals.org/cgi/reprint/14/12/993 [Accessed 28th July 2010] | ||
+ | <html> | ||
+ | |||
+ | </div> | ||
<div class="accordionButton">Detailed Description of the Model</div> | <div class="accordionButton">Detailed Description of the Model</div> | ||
- | <div class="accordionContent"></div> | + | <div class="accordionContent"> |
+ | </html> | ||
+ | =Model based on Law of Mass Action (Weeks 6 - 8)= | ||
+ | ==Michaelis Menten kinetics does not apply== | ||
+ | We cannot use Michaelis-Menten kinetics because of its preliminary assumptions, which our system does not fulfil. | ||
+ | These assumptions are: | ||
+ | *'''<math>V_{max}</math> is proportional to the overall concentration of the enzyme.''' | ||
+ | But we are producing enzyme, so <math>V_{max}</math> will change! Therefore, the conservation <math>E_0 = E + E_S</math> does not hold for our system. | ||
+ | |||
+ | *'''Substrate >> Enzyme.''' | ||
+ | Since we are producing both substrate and enzyme, we have roughly the same amount of substrate and enzyme. | ||
+ | |||
+ | *'''Enzyme affinity to substrate has to be high.''' | ||
+ | |||
+ | Therefore, the model above is not representative of the enzymatic reaction. As we cannot use the Michaelis-Menten model we will have to solve from first principle (which just means writing down all of the biochemical equations and solving for these in Matlab). | ||
+ | |||
+ | ==Change of output== | ||
+ | Instead of GFP, it is now dioxygenase acting on catechol (activating it into colourful form). Catechol will be added to the bacteria manually (i.e. the bacteria will not produce Catechol). Hence, in our models dioxynase is going to be treated as an output as this enzyme is the only activator of catechol in our system. This means that the change of catechol into its colourful form is dependent on dioxygenase concentration. | ||
+ | |||
+ | == Models== | ||
+ | ===Model preA: Simple Production of Dioxygenase=== | ||
+ | [[Image:Simple_production.JPG|400px|thumb|center|alt=A|Simple production of Dioxygenase]] | ||
+ | This model includes transcription and translation of the dioxygenase. | ||
+ | It does not involve any amplification steps. It is our control model against which we will be comparing the results of other models. | ||
+ | *[http://www.openwetware.org/wiki/Image:Wiki_Model_preA.txt Here is the Matlab code for Model pre-A.] | ||
+ | |||
+ | |||
+ | ===Model A: Activation of Dioxygenase by TEV enzyme=== | ||
+ | [[Image:1-step_amplification.JPG|400px|thumb|center|alt=A|1-step amplification]] | ||
+ | The reaction can be rewritten as: | ||
+ | |||
+ | <math>{TEV} + {split Dioxygenase} \rightleftarrows {TsD} \rightarrow {TEV} + {Dioxygenase}</math> | ||
+ | |||
+ | This is a simple enzymatic reaction, where TEV is the enzyme, Dioxygenase the product and split Dioxygenase the substrate. | ||
+ | Choosing <math>k_1, k_2, k_3</math> as reaction constants, the reaction can be rewritten in these four sub-equations: | ||
+ | |||
+ | #<math>\dot{T} = -k_1[T][sD] + (k_2+k_3)[TsD] + s_T - d_T[T]</math> | ||
+ | #<math>\dot{sD}= -k_1[T][sD] + k_2[TsD] + s_{sD} - d_{sD}[sD]</math> | ||
+ | #<math>\dot{TsD} = k_1[T][sD] - (k_2+k_3)[TsD] - d_{TsD}[TsD]</math> | ||
+ | #<math>\dot{D} = k_3[TsD] - d_D[D]</math> | ||
+ | |||
+ | These four equations were implemented in Matlab, using a built-in function (ode45) which solves ordinary differential equations. | ||
+ | *[http://www.openwetware.org/wiki/Image:Wiki_Model_A.txt Here is the Matlab code for Model A.] | ||
+ | |||
+ | |||
+ | '''Implementation in TinkerCell''' | ||
+ | |||
+ | Another approach to model the amplification module would be to implement it in a program such as TinkerCell (or CellDesigner). This would be useful to check whether the Matlab model works. | ||
+ | [[Image:Tinkercell_Model_A.PNG|450px|thumb|center|alt=A|Network implemented in TinkerCell]] | ||
+ | |||
+ | |||
+ | ===Model B: Activation of Dioxygenase by TEV or activated split TEV enzyme=== | ||
+ | [[Image:2-step_amplification.JPG|400px|thumb|center|alt=A|2-step amplification]] | ||
+ | |||
+ | This version includes the following features: | ||
+ | *2 amplification steps (TEV and split TEV) | ||
+ | *Split TEV is specified to have a and b parts | ||
+ | *TEVa is forbidden to interact with TEVa (though in reality there could be some affinity between the two). Same for the interaction between Tevb and Tevb | ||
+ | *Both TEV and TEVs are allowed to activate dioxygenase | ||
+ | *Dioxygenase is assumed to be active as a monomer | ||
+ | *Activate split TEV (TEVs) is not allowed to activate sTEVa or sTEVb (this kind of interaction is accounted for in the next model version) | ||
+ | *This model does not include any specific terms for time delays | ||
+ | *[http://www.openwetware.org/wiki/Image:Wiki_Model_B.txt Here is the Matlab code for Model B.] | ||
+ | |||
+ | |||
+ | ===Model C: Further improvement=== | ||
+ | This model has not been implemented because of the results from Models A and B. | ||
+ | |||
+ | This version adds the following features: | ||
+ | * activated split TEV (TEVs) is allowed to activate not only sD but sTEVa and sTEVb | ||
+ | |||
+ | ==Results== | ||
+ | The first results that were obtained seemed to be flawed since they indicated negative concentrations would be obtained from the amplification step. In particular, for concentrations smaller than <math>10^{-4} mol/dm^3</math> the results were inconclusive since they were oscillating around zero. | ||
+ | We realised that this could be due to the ode-solver that we were using (ode45 in Matlab). | ||
+ | |||
+ | Trying to correct this problem with the ode-solver, the following precuations were implemented: | ||
+ | *NonNegative function in Matlab preventing solver from reaching negative values - still some marginally negative values show | ||
+ | *Scaling - all the values were scaled up by a factor of <math>10^6</math> as working on small numbers could be problematic for Matlab. Once the result is generated by the solver the resulting matrix is scaled back down by <math>10^6</math>. | ||
+ | |||
+ | When we entered the real production and degradation rates into our model, we once again obtained nagetive values. This was due to our set of differential equations being stiff. Since ode45 cannot solve stiff differential equations, we had to switch to using ode15s - an ode-solver designed to handle stiff equations. | ||
+ | |||
+ | |||
+ | ===Model pre-A=== | ||
+ | This is the result of the simulation of simple production of Dioxygenase. It can be seen that the concentration will tend towards a final value of approximately <math>8*10^{-6} mol/dm^3</math>. This final value is dependent on the production rate (which has been estimated for all of the models). | ||
+ | [[Image:Model_prea.bmp|450px|thumb|center|alt=A|Results of the Matlab simulation of Model preA]] | ||
+ | |||
+ | ===Model A=== | ||
+ | *'''Initial Concentration''' | ||
+ | The initial concentration of split Dioxygenase, <math>c_0</math>, determines whether the system is amplifying. The minimum concentration for any amplification to happen is <math>10^{-5} mol/dm^3</math>. If the initial concentration of split Dioxygenase is higher, then the final concentration of Dioxygenase will be higher as well (see graphs below). '''Note that the obtained threshold value is higher than the maximum value that can be generated in the cell according to Model pre-A.''' | ||
+ | [[Image:Comparison_a%2C_prea.bmp|450px|thumb|center|alt=A|Comparison between Model pre-A and Model A. Initial concentration of split Dioxygenase: <math>10^{-5} mol/dm^3</math>]] | ||
+ | |||
+ | *'''Changing <math>K_m</math>:''' | ||
+ | <math>K_m</math> is indirectly proportional to the "final concentration" (which is the concentration at the end of the simulation), i.e. the bigger the value of <math>K_m</math>, the smaller the "final concentration" will be. | ||
+ | Different <math>K_m</math> values determine how quickly the amplification will take place. | ||
+ | |||
+ | (Also, it was found that the absolute value of <math>k_1</math> and <math>k_2</math> entered into Matlab does not change the outcome as long as the ratio between them (<math>K_m</math>≈<math>k_2/k_1</math>) is kept constant. This is important when simulating (in case entering very high values for <math>k_1</math> and <math>k_2</math> slows down the simulation). | ||
+ | |||
+ | *'''Changing <math>k_{cat}</math>''' | ||
+ | If <math>k_{cat} = k_3</math> is increased, the model predicts that the dioxygenase concentration will rise quicker and to higher values. This indicates that <math>k_3</math> is the slowest step in the enzymatic reaction. | ||
+ | |||
+ | *'''Changing production rate''' | ||
+ | At the moment, our biggest source of error could be the production rate, which we could not obtain from literature. Hence, we had to estimate the value of the production rate (see variables). We hope to be able to take a measurement of this value in the lab as it has a big effect on model's behaviour. | ||
+ | |||
+ | |||
+ | '''Sensitivity of Model A (20/08/2010)''' | ||
+ | |||
+ | We want to determine how our system reacts if different parameters are changed. This is to find out which parameters our system is very sensitive to. | ||
+ | {| class="wikitable" style="text-align: center; width: 80%; height: 170px;" border="1" | ||
+ | |- | ||
+ | ! Parameter !! Sensitivity | ||
+ | |- | ||
+ | | Initial concentration of split Dioxygenase || Change of one order of magnitude in the initial concentration, <math>c_0</math>, gives change of one order of magnitude in the output concentration (range: <math>1>c_0>10^{-5}</math>). Sensitivity is lost for extremely high or low values. | ||
+ | |- | ||
+ | | <math>K_m</math> ||Change of one order of magnitude results in change of output concentration by one order of magnitude (<math>0.01<K_m<100</math>). At values smaller than 0.01, the sensitivity is affected. For higher values than 100 the sensitivity is slightly higher than the change of order of magnitude. | ||
+ | |- | ||
+ | | <math>k_{cat}</math> || <math>k_{cat}</math> is proportional to dioxygenase production (1-to-1 sensitivity for all values) for an initial concentration of <math>0.01 mol/dm^3</math>. For very high initial concentrations, the system is very sensitive to changes in <math>k_{cat}</math>. | ||
+ | |- | ||
+ | | Production rate of TEV || 1-1 sensitivity for most values. At some point the system’s response is limited by the initial concentration of sD, so for very high TEV production rates not much change is observed. | ||
+ | |- | ||
+ | | Production rate of split Dioxygenase || Not much influence on 1-step amplification. However, the value seems to be crucial for simple production of Dioxygenase (1-1 order of magnitude sensitivity). | ||
+ | |- | ||
+ | | Degradation rates || Sensitive within the relevant range. Not very sensitive for values smaller than <math>10^{-6}</math>. For high degradation rates (1>degradation rate>0.01): unexplainable behaviour. | ||
+ | |} | ||
+ | |||
+ | Hence, the system is sensitive to most of the constants (given a particular range of values). The most crucial one, however, seems to be the initial concentration of split Dioxygenase. | ||
+ | |||
+ | ===Model B=== | ||
+ | *'''Initial Concentration''' | ||
+ | The initial concentration of split Dioxygenase, <math>c_0</math>, determines whether the system is amplifying. | ||
+ | |||
+ | If the initial concentration is changed, the observed behaviour is similar to the one from Model A. If the initial concentration of split Dioxygenase is increased, then the final concentration of Dioxygenase will increase as well (see graphs below). | ||
+ | *'''Model A vs. B''' | ||
+ | Running both models with the same initial conditions (<math>c_0=10^{-5} mol/dm^3</math>), it has been noted that Model B does not generate a siginificant amplfication over Model A. Hence, it would be more sensible to integrate a one step amplification module into our system. | ||
+ | |||
+ | [[Image:Comparison_prea%2C_a%2C_b.bmp|450px|thumb|center|alt=A|Comparison between Models pre-A, A and B]] | ||
+ | |||
+ | ===Colour response model (25/08/2010)=== | ||
+ | Initially, dioxygenase was being treated as "output" because the last reaction (catechol and dioxygenase react to form a coloured output) is common to all models. We concluded that 2-step amplification (Model B) presented little improvement over 1-step amplification (Model A). However, we realized that it is very important to add this last reaction, since it is an enzymatic reaction, hence an amplification step by itself. Adding it to all models means that Model A becomes a 2-step amplifier and Model B becomes a 3-step amplifier. Having previously drawn the conclusion that 2-step amplification is not much better than 1-step amplification, we decided to model this last step to deduce whether our initial conclusion was valid. If it was true, this would mean that our construct is not innovative at all. | ||
+ | |||
+ | The important information about the last amplification step is that the coloured compound (i.e. the product of the last enzymatic reaction) is toxic to the cells. | ||
+ | It is suggested that product of Catechol destroys the cell membrane by inhibiting lipid peroxidation. It causes significant changes in the structure and functioning of membrane components (e.g. disruption of membrane potential, removal of lipids and proteins, loss of magnesium and calcium ions). These effects cause the loss of membrane functions, leading to cell death. | ||
+ | |||
+ | Since the product of Catechol acts on the cell membrane, it might not affect our enzymatic reaction immediately. In our simulation, we will try to model immediate cell death as well as neglecting the effect that the coloured output has on the cell. Comparing these two models will show if there are significant differences in the results. | ||
+ | |||
+ | |||
+ | Important Articles: | ||
+ | *Toxicity of three phenolic compounds and their mixtures on the gram-positive bacteria Bacillus subtilis in the aquatic environment[http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V78-4XXNV4Y-4-7&_cdi=5836&_user=217827&_pii=S0048969709011668&_orig=search&_coverDate=02%2F01%2F2010&_sk=995919994&view=c&wchp=dGLbVzz-zSkzk&_valck=1&md5=f8dca6227c29db659ddbeb588ad115e7&ie=/sdarticle.pdf] | ||
+ | *Toxic Effects of Catechol and 4-Chlorobenzoate Stresseson Bacterial Cells[http://www.msk.or.kr/jsp/downloadPDF1.jsp?fileName=393-15.pdf] | ||
+ | |||
+ | ===Initial conclusions(26/08/2010)=== | ||
+ | Despite our model not working entirely correctly (some negative concentrations were obtained), it was possible to deduce several points. | ||
+ | *The images presented below show cathecol being added at 3 different points in time. '''Cross-section''' refers to a point in time at which concentration of dioxygenase in the amplified systems crosses the concentration of the non-amplified system. From these graphs it can be seen that the output amplification is only visible after the cross-section has been reached. Note that these simulations were run for 1M solutions of catechol (which is quite high). This allows to see the differences between various amplification models easily. | ||
+ | {| | ||
+ | |- | ||
+ | || [[Image:Cathecol_model_before_cross.png|350px|thumb|center|alt=A|Concentration of coloured compound for catechol being added before the cross-section is reached ]] | ||
+ | |||
+ | || [[Image:Cathecol_model_at_cross.png|350px|thumb|center|alt=A|Concentration of coloured compound for catechol being added when the cross-section is reached ]] | ||
+ | |||
+ | || [[Image:Cathecol_model_after_cross.png|350px|thumb|center|alt=A|Concentration of coloured compound for catechol being added after the cross-section has been reached ]] | ||
+ | |||
+ | |} | ||
+ | *We noticed that amplification is only benefical if the initial concentration of catechol is quite high (>0.01M). For smaller concentrations of catechol, the dioxygenase conetration in different systems do not seem to be crucial for the speed of response (no difference between all 3 models). Hence, for small catechol concentrations the amplified systems are redundant (dioxygenase is overproduced) as concnetration of dioxygenase from simple production seems to be high enough to convert catechol almost instantenously. Amplification models become useful when there is a lot of substrate present to act on (ie. high concentration of catechol). Therefore, we need to determine the threshold value of coloured output for visibility. It will be a crucial factor in deciding whether the amplifiers designed by our team obtain the response faster than by simple production. | ||
+ | |||
+ | ===Problems with the simulation=== | ||
+ | Once we have implemented the colour change into the models, we noticed that there are some inconclusive results. After adding catechol, some concentrations were reaching negative values. We checked our equations and constants but could not find a mistake. Hence, we concluded that there must be something wrong with the way that MatLab evaluates the equation or deals with the numbers. The problem seems to originate from the very rapid concentration change of catechol which disrupts the whole system. | ||
+ | |||
+ | '''Prospective solutions''' | ||
+ | *'''Implementation in TinkerCell (31/08/2010)''' | ||
+ | We hoped that TinkerCell imposes non-negative conditions on its solutions. Hence, we implemented the whole amplification model (including coloured output) in TinkerCell. However, we realized that TinkerCell does not deal well with very high or low numbers (For example, values higher than <math>10^5</math> are not acceptable - this is important since our rate constants (<math>k_1</math>) are usually bigger than <math>10^5</math>. Also, the low degradation rates (<math>10^{-9}</math>) result in a zero output line). | ||
+ | However, TinkerCell can still be used for testing that our Matlab programs behave the way we anticipated (by using default values of 1), as well as producing nice diagrams of our system. | ||
+ | *'''Varying ODE solver options in MatLab (31/08/2010)''' | ||
+ | We had a close look at the ODE solver options in MatLab. However, we were already using the one that produced the most reasonable results (ode15s). We found that decreasing Relative and Aboslute tolerances (to values as small as <math>10^{-15}</math>) significantly improved the simulation. However, this is not an ultimate solution as in the simulations negative numbers still appear (order of <math>10^{-15}</math>). We decided that such small negative concentrations were acceptable. We also decided that the point of interest lies between the first 100 to 150 seconds after adding catechol, while concentrations hit the negatve values at much higher time values. | ||
+ | |||
+ | The images below show the influence of the relative and absolute tolerance values on the model. Note that it was important to allow the ode-olver to adjust the time step automatically, as big time steps (1 second) were generating wrong answers for the catechol model. Adjusting time steps manually to very small values was not efficient (the whole simulation does not require very high definition simluation). | ||
+ | {| | ||
+ | |- | ||
+ | || [[Image:High_tolerance_inaccurate.png|420px|thumb|center|alt=A|Graphs representing conctrations of catechol, dioxygenase, catechol-dioxygenase complex and colur compound when adding catechol at t=2000s. Note the negative values that concentration of catechol hits.]] | ||
+ | |||
+ | || [[Image:Low_tolerance_accurate.png|600px|thumb|center|alt=A|Here the Relative and the Absolute Tolerances are chosen to be <math>10^{-18}</math> and <math>10^{-12}</math>, respectively. When lookin closely it is clear that catechol and the catechol-dioxygenase complex reach slightly negative values (eg. <math>10^{-120}</math>]] | ||
+ | |} | ||
+ | *'''Using SimBiology to model (31/08/2010)''' | ||
+ | We hoped that SimBiology could be more suited for our modelling than using ode-solvers, so we implemented our models into SimBiology. This package offers an interactive user interface similar to Tinker Cell, but uses MatLab to simulate. Initially, we confirmed that our simple production model (Model PreA) and 1-step amplification model (Model A) implemented in SimBiology generated exactly the same results as our ODE equation based models. The interface allowed us to have clearer control over paramters. It also allowed modelling special events, for example, adding catechol at certain point in time. Previously we had to split simulation into to two parts. | ||
+ | |||
+ | ==References== | ||
+ | #Chen, H. et al (2009) Toxicity of three phenolic compounds and their mixtures on the gram-positive bacteria Bacillus subtilis in the aquatic environment. Science of the Total Environment. [Online] 408(2010), 1043-1049. Available from: http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V78-4XXNV4Y-4-7&_cdi=5836&_user=217827&_pii=S0048969709011668&_orig=search&_coverDate=02%2F01%2F2010&_sk=995919994&view=c&wchp=dGLbVzz-zSkzk&_valck=1&md5=f8dca6227c29db659ddbeb588ad115e7&ie=/sdarticle.pdf [Accessed 6th September 2010] | ||
+ | #Park, S., Ko, Y. & Kim, C. (2001) Toxic Effects of Catechol and 4-Chlorobenzoate Stresses on Bacterial Cells. The Journal of Microbiology. [Online] 39(3), 206-212. Available from: http://www.msk.or.kr/jsp/downloadPDF1.jsp?fileName=393-15.pdf [Accessed 6th September 2010] | ||
+ | #Habibi Nazhad, B. (2008) E. coli Statistics. [Online] Available from: http://gchelpdesk.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi [Accessed 6th September 2010] | ||
+ | <html> | ||
+ | |||
+ | </div> | ||
<div class="accordionButton">Paramters & Constants</div> | <div class="accordionButton">Paramters & Constants</div> | ||
- | <div class="accordionContent"></div> | + | <div class="accordionContent"> |
+ | |||
+ | <body style="background-color:FFFFCC"> | ||
+ | <h1>Constants for the Output Amplification Model</h1> | ||
+ | </html> | ||
+ | |||
+ | <html> | ||
+ | <table width="800px" border="0"> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFFF66;height:50px;width:100;text-align:center"><b>Type of Constant</b> | ||
+ | </td> | ||
+ | <td style="background-color:#FFFF99;height:50px;width:700;text-align:center"><b>Derivation of Value</b> | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>TEV Enzyme Dynamics</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">Enzymatic Reaction: E+S <var>↔</var> ES <var>→</var> E+P | ||
+ | <br /> | ||
+ | Let | ||
+ | <ul> | ||
+ | <li>k<sub>1</sub> = rate constant for E+S <var>→</var> ES | ||
+ | <li>k<sub>2</sub> = rate constant for E+S <var>←</var> ES</math> | ||
+ | <li>k<sub>cat</sub> = rate constant for ES <var>→</var> E+P | ||
+ | </ul> | ||
+ | We know that K<sub>m</sub> = (k<sub>cat</sub> + k<sub>2</sub>)/k<sub>1</sub> | ||
+ | |||
+ | Assuming that k<sub>cat</sub> << k<sub>2</sub> << k<sub>1</sub>, we can rewrite K<sub>m</sub> <var>≈</var> k<sub>2</sub>/k<sub>1</sub> | ||
+ | <br /> | ||
+ | From this paper <a href="http://peds.oxfordjournals.org/cgi/reprint/14/12/993">[1]</a> the constants for TEV can be found: | ||
+ | <br /> | ||
+ | For example, for wildtype TEV: K<sub>m</sub> = 0.061<var>±</var>0.010mM and k<sub>cat</sub> = 0.16<var>±</var>0.01s<sup>-1</sup> | ||
+ | <br /> | ||
+ | These values correspond with our assumption that k<sub>cat</sub> = 0.1 s<sup>-1</sup> and K<sub>m</sub> = 0.01 mM. | ||
+ | <br /> | ||
+ | Hence, we can estimate the following orders of magnitude for the rate constants: | ||
+ | <br /> | ||
+ | k<sub>1</sub> = 10<sup>8</sup>M<sup>-1</sup>s<sup>-1</sup> | ||
+ | <br /> | ||
+ | k<sub>2</sub> = 10<sup>3</sup>s<sup>-1</sup> | ||
+ | <br /> | ||
+ | Using these values should be a good approximation for our model. | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:130px;text-align:center;"><b>Degradation rate (common for all)</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:130px;text-align:center;">Assumption: To be approximated by cell division (dilution of media) as none of the proteins are involved in any active degradation pathways | ||
+ | Growth rate, gr (divisions/h): 0.53 <var>≤</var> gr <var>≤</var> 2.18 <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC235685/pdf/jbacter00326-0019.pdf">[2]</a> | ||
+ | <br /> | ||
+ | Hence on average, gr = 1.5 divisions per hour, which gives one division every 40mins | ||
+ | <br /> | ||
+ | To deduce degradation rate we use the following formula: | ||
+ | <br /> | ||
+ | <var>τ</var><sub>1/2</sub> = ln2/k, where <var>τ</var><sub>1/2</sub> = 0.667 hours and k = degradation rate | ||
+ | <br /> | ||
+ | k = ln2/<var>τ</var><sub>1/2</sub> = 0.000289s<sup>-1</sup> | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>Production rate (TEV and Dioxygenase)</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">We had difficulties finding values of the production rate in the literature and we hope to be able to perform experiments to obtain those values (for TEV protease and catechol 2,3-dioxygenase). Before any values can be obtained from the Lab, we suggest very simplistic approach for estimating production rates. | ||
+ | <br /> | ||
+ | We have found production rates for two arbitrary proteins in E.Coli. We want to get estimates of production rates by comparing the lengths of the proteins (number of amino-acids). | ||
+ | <br /> | ||
+ | As this approach is very vague, it is important to realise its limitations and inconsistencies: | ||
+ | <ul> | ||
+ | <li>Values are taken from E.Coli not B.sub.</li> | ||
+ | <li>The two production rates are of the same value for quite different amino-acid number which indicates that protein folding is limiting the production rates.</li> | ||
+ | </ul> | ||
+ | LacY production = 100 molecules/min<a href="http://bionumbers.hms.harvard.edu/bionumber.aspx?s=y&id=100738&ver=0&hlid=29205">[3]</a> (417 Amino Acids<a href="http://www.uniprot.org/uniprot/P02920">[4]</a>) | ||
+ | <br /> | ||
+ | LacZ production = 100 molecules/min<a href="http://bionumbers.hms.harvard.edu/bionumber.aspx?s=y&id=100737&ver=0&hlid=29206">[5]</a> (1024 AA<a href="http://www.uniprot.org/uniprot/P00722">[6]</a>) | ||
+ | <br /> | ||
+ | Average production ≈ 100molecules/min 720 AA | ||
+ | <br /> | ||
+ | This gives us: | ||
+ | TEV production ≈ 24 molecules/min = 0.40 molecules/s (3054 AA<a href="http://www.uniprot.org/uniprot/P04517">[7]</a>) | ||
+ | <br /> | ||
+ | As production rate needs to be expressed in concentration units per unit volume, the above number is converted to mols/s and divided by the cell volume: 2.3808x10<sup>-10</sup> mol/dm<sup>3</sup>/s | ||
+ | <br /> | ||
+ | C23D production ≈ 252 molecules/min = 4.2 molecules/s (285 AA<a href="http://www.uniprot.org/uniprot/P54721#section_x-ref">[8]</a>) → 2.4998x10<sup>-9</sup> mol/dm<sup>3</sup>/s | ||
+ | <br /> | ||
+ | We will treat these numbers as guiding us in terms of range of orders of magnitudes. We will try to run our models for variety of values and determine system’s limitations. | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>Kinetic Parameters of Dioxygenase</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">Initial velocity of the enzymatic reaction was investigated at pH 7.5 and 30 °C. | ||
+ | <br /> | ||
+ | Wild type (used for our simulations): K<sub>m</sub> = 10 <var>μ</var>M; k<sub>cat</sub> = 52s<sup>-1</sup> | ||
+ | <br /> | ||
+ | Mutated type: K<sub>m</sub> = 40 <var>μ</var>M; k<sub>cat</sub> = 192s<sup>−1</sup> | ||
+ | <br /> | ||
+ | Consequently, the ratio of K<sub>m</sub>/k<sub>cat</sub> of the mutant (K<sub>m</sub>/k<sub>cat</sub> = 4.8) is slightly lower than the ratio of the wild type (K<sub>m</sub>/k<sub>cat</sub> = 5.2), indicating that the mutation has little effect on the catalytic efficiency <a href="http://www.springerlink.com/content/e3718758m5052214/fulltext.pdf">[9]</a>. | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>Dimensions of B.sub cell</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">Dimensions of B.sub (cylinder/rod shape) in rich media: | ||
+ | <br /> | ||
+ | diameter: d = 0.87<var>μ</var>m; length: l = 4.7<var>μ</var>m | ||
+ | <br /> | ||
+ | This gives: Volume= <var>π</var>d<sup>2</sup>l/4 = 2.793999<var>μ</var>m<sup>3</sup> <var>≈</var> 2.79x10<sup>-15</sup> dm<sup>3</sup> | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>Production Rate of split TEV</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">Assuming that both parts of split TEV are half the size of the whole TEV (3054/2=1527 AA). | ||
+ | <br /> | ||
+ | The length of the coil is 90 AA. | ||
+ | <br /> | ||
+ | The whole construct is then: 1617 AA | ||
+ | <br /> | ||
+ | Therefore, split TEV production rate ≈ 1.2606x10<sup>-10</sup> mol/dm<sup>3</sup>/s | ||
+ | </td> | ||
+ | </tr> | ||
+ | |||
+ | <tr> | ||
+ | <td style="background-color:#FFCC66;height:100px;text-align:center;"><b>Relevant concentrations of Catechol</b> | ||
+ | </td> | ||
+ | <td style="background-color:#e7e7e7;height:100px;text-align:center;">We have catechol in the lab in powder form so we are only limited by it's solubility. | ||
+ | <br /> | ||
+ | For a concentration of 0.1 M with built up levels of dioxygenase the colour change happens within seconds. | ||
+ | <br /> | ||
+ | We will run our models for 0.1M ± several orders of magnitude to determine the smallest catechol concentration that will give a significant difference between the simple production response and the amplified response. | ||
+ | </td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | </html> | ||
+ | |||
+ | <html> | ||
+ | <h2>References</h2> | ||
+ | <ol> | ||
+ | <li>Kapust, R. et al (2001) Tobacco etch virus protease: mechanism of autolysis and rational design of stable mutants with wild-type catalytic proficiency. Protein Engineering. [Online] 14(12), 993-1000. Available from: http://peds.oxfordjournals.org/content/14/12/993.full.pdf+html [Accessed 20th August 2010]</li> | ||
+ | <li>Sargent, M. (1975) Control of Cell Length in Bacillus subtilis. Journal of Bacteriology. [Online] 123(1), 7-19. Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC235685/pdf/jbacter00326-0019.pdf [Accessed 20th August 2010]</li> | ||
+ | <li>Milo, R., Jorgensen, P. & Springer, M. (2007) BioNumbers. [Online] Available from: http://bionumbers.hms.harvard.edu/bionumber.aspx?s=y&id=100738&ver=0&hlid=29205 [Accesed 25th August 2010]</li> | ||
+ | <li>UniProt Consortium (2002-2010) UniProt. [Online] Available from: http://www.uniprot.org/uniprot/P02920 [Accessed 24th August 2010]</li> | ||
+ | <li>Milo, R., Jorgensen, P. & Springer, M. (2007) BioNumbers. [Online] Available from: http://bionumbers.hms.harvard.edu/bionumber.aspx?s=y&id=100737&ver=0&hlid=29206 [Accesed 25th August 2010]</li> | ||
+ | <li>UniProt Consortium (2002-2010) UniProt. [Online] Available from: http://www.uniprot.org/uniprot/P00722 [Accessed 24th August 2010]</li> | ||
+ | <li>UniProt Consortium (2002-2010) UniProt. [Online] Available from: http://www.uniprot.org/uniprot/P04517 [Accessed 24th August 2010]</li> | ||
+ | <li>UniProt Consortium (2002-2010) UniProt. [Online] Available from: http://www.uniprot.org/uniprot/P54721#section_x-ref [Accessed 24th August 2010]</li> | ||
+ | <li>Wei, J. et al (2009) Rational Design of Catechol-2, 3-dioxygenase for Improving the Enzyme Characteristics. Appl Biochem Biotechnol. [Online] 162, 116-126. Available from: http://www.springerlink.com/content/e3718758m5052214/fulltext.pdf [Accessed 25th August 2010]</li> | ||
+ | </ol> | ||
+ | </body> | ||
+ | |||
+ | </div> | ||
<div class="accordionButton">Results & Conclusions</div> | <div class="accordionButton">Results & Conclusions</div> | ||
- | <div class="accordionContent"></div> | + | <div class="accordionContent"> |
+ | </html> | ||
+ | ===Final Conclusions (31/08/2010)=== | ||
+ | *'''Changing time when catechol is added''' | ||
+ | If Catechol is added before t= 1000s, then the coloured output will reach its threshold value faster by simple production. If Catechol is added when t>1000s, then the coloured output will increase (marginally) faster through the amplification step in Model A. There does not seem to be a significant difference between the two models (Model preA and Model A). These observations are true for intial concentration of dioxygenase equal to <math>10^{-5}mol/dm^3</math>. However, we noticed that if the initial concentration is raised to <math>10^{-4}mol/dm^3</math>, then Model A can be more beneficial than Model preA after only 100 seconds. | ||
+ | |||
+ | Hence, the question arises whether the concentration of protein in the cell can be as high as <math>10^{-4}mol/dm^3</math>. Our simple production model predicted that the concentration of protein could not reach such a high value. However, we decided to research more on ribosomal concentrations in bacteria to determine whether it is possible to establish such a high concentration in the cell. | ||
+ | |||
+ | On the website E.coli Statistics [http://gchelpdesk.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi] it is stated that number of ribosomal proteins per cell is 900,000. | ||
+ | In a cellular volume of order of <math>1\mu m^3 = 10^{-15}dm^3=10^{-15}L</math>, the above number of ribosomes converts to <math>1.5*10^{-3}mol/L</math>. This means that a concentration of <math>10^{-4}mol/dm^3</math> is not completely out of scale. | ||
+ | |||
+ | *'''Changing concentration of catechol added''' | ||
+ | There seem to be 3 regions of catechol concentration that influence the system in different ways. These regions are: c>1M, 1M>c>0.01M, 0.01M>c. The boundaries of these regions tend to vary depending on the choice of other initial conditions. The values given above apply to boundary conditions that are currently considered to be physiologically relevant. Varying the initial concentration of catechol within the highest region does not result in any change of colour output response (It is possible that all enzymes are occupied and the solution is over saturated with catechol). In the middle region the catechol concentration influences the amplfication. Amplification decreased when the concentration tends towards 0.01M. When this region is entered, there is no difference in output production by the two models. | ||
+ | |||
+ | *'''Cell death''' | ||
+ | The coloured product of catechol kills cells by destroying the cell membrane. | ||
+ | However, we do not know how quickly the cells will die. | ||
+ | Therefore, we examined two different cases: immediate cell death and negligible cell death (i.e. cells death is negligible because it takes too long) | ||
+ | |||
+ | Running the simluation in Matlab (not Simbiology!), our conclusions are: | ||
+ | #Immediate cell death slows down production of coloured output. Depending on the threshold concentration this can delay the detectable response by a few minutes. | ||
+ | #If Catechol is added before t=1000s, then cell death slows down the response considerably. | ||
+ | #In case of cells being modelled as alive, the difference between the amplified and the simple production model is smaller than it is in case of cell death. | ||
+ | |||
+ | Since it appears that the time of cell death is important, we decided to discuss this issue with Wolf and Harriet. Referring to this paper [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V78-4XXNV4Y-4-7&_cdi=5836&_user=217827&_pii=S0048969709011668&_orig=search&_coverDate=02%2F01%2F2010&_sk=995919994&view=c&wchp=dGLbVzz-zSkzk&_valck=1&md5=f8dca6227c29db659ddbeb588ad115e7&ie=/sdarticle.pdf (1)] we decided that cell death induced by catechool is a very slow process (we estimate that it will take a few hours) in comparison to the time scale that we are interested in (several seconds to minutes). | ||
+ | {| | ||
+ | |- | ||
+ | |[[Image:Alive_cells.png|500px|thumb|center|alt=A|Colour response model for cells being kept alive]] | ||
+ | | | ||
+ | [[Image:Dead_cells.png|500px|thumb|center|alt=A|Colour response model with the same parameters for cells being killed instantaneously by catechol ]] | ||
+ | |} | ||
+ | |||
+ | <html> | ||
+ | </div> | ||
</body> | </body> | ||
</html> | </html> | ||
|} | |} |
Revision as of 08:36, 13 October 2010
Output Amplification Model | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Objectives
We came up with a simple concept of output amplification, which is enhanced by using enzymes. It is beneficial for us to model the behaviour of our design so that we will be able to answer the following questions.
Abandoned Initial Attempts
Model based on Michaelis Menten Kinetics (Weeks 4 and 5)MotivationWe came up with a simple concept of output amplification, which is enhanced by using enzymes. It is beneficial for us to model the behaviour of our design so that we will be able to answer the following questions.
Modelling should allows us to make a decision on which design is the most efficient one. First ModelHIV1
TEV
Improved Model which accounts for enzyme reactions (28/07/2010)TEV
Implementation in MatlabThe Matlab code for the different stages of amplification and diagrams can be found [http://www.openwetware.org/wiki/Image:Modelling.docx here]. Kinetic constants
ConclusionWe were not able to obtain all the necessary constants. Hence, we decided to make educated guesses about possible relative values between the constants as well as varying them and observing the change in output. As the result, we concluded that the amplification happens at each amplification level proposed. The magnitude of amplification varies depending on the constants. There is not much difference between using TEV or HIV1. References
Detailed Description of the Model
Model based on Law of Mass Action (Weeks 6 - 8)Michaelis Menten kinetics does not applyWe cannot use Michaelis-Menten kinetics because of its preliminary assumptions, which our system does not fulfil. These assumptions are:
But we are producing enzyme, so <math>V_{max}</math> will change! Therefore, the conservation <math>E_0 = E + E_S</math> does not hold for our system.
Since we are producing both substrate and enzyme, we have roughly the same amount of substrate and enzyme.
Therefore, the model above is not representative of the enzymatic reaction. As we cannot use the Michaelis-Menten model we will have to solve from first principle (which just means writing down all of the biochemical equations and solving for these in Matlab). Change of outputInstead of GFP, it is now dioxygenase acting on catechol (activating it into colourful form). Catechol will be added to the bacteria manually (i.e. the bacteria will not produce Catechol). Hence, in our models dioxynase is going to be treated as an output as this enzyme is the only activator of catechol in our system. This means that the change of catechol into its colourful form is dependent on dioxygenase concentration. ModelsModel preA: Simple Production of DioxygenaseFile:Simple production.JPG Simple production of Dioxygenase This model includes transcription and translation of the dioxygenase. It does not involve any amplification steps. It is our control model against which we will be comparing the results of other models.
Model A: Activation of Dioxygenase by TEV enzymeFile:1-step amplification.JPG 1-step amplification The reaction can be rewritten as: <math>{TEV} + {split Dioxygenase} \rightleftarrows {TsD} \rightarrow {TEV} + {Dioxygenase}</math> This is a simple enzymatic reaction, where TEV is the enzyme, Dioxygenase the product and split Dioxygenase the substrate. Choosing <math>k_1, k_2, k_3</math> as reaction constants, the reaction can be rewritten in these four sub-equations:
These four equations were implemented in Matlab, using a built-in function (ode45) which solves ordinary differential equations.
Another approach to model the amplification module would be to implement it in a program such as TinkerCell (or CellDesigner). This would be useful to check whether the Matlab model works. File:Tinkercell Model A.PNG Network implemented in TinkerCell
Model B: Activation of Dioxygenase by TEV or activated split TEV enzymeFile:2-step amplification.JPG 2-step amplification This version includes the following features:
Model C: Further improvementThis model has not been implemented because of the results from Models A and B. This version adds the following features:
ResultsThe first results that were obtained seemed to be flawed since they indicated negative concentrations would be obtained from the amplification step. In particular, for concentrations smaller than <math>10^{-4} mol/dm^3</math> the results were inconclusive since they were oscillating around zero. We realised that this could be due to the ode-solver that we were using (ode45 in Matlab). Trying to correct this problem with the ode-solver, the following precuations were implemented:
When we entered the real production and degradation rates into our model, we once again obtained nagetive values. This was due to our set of differential equations being stiff. Since ode45 cannot solve stiff differential equations, we had to switch to using ode15s - an ode-solver designed to handle stiff equations.
Model pre-AThis is the result of the simulation of simple production of Dioxygenase. It can be seen that the concentration will tend towards a final value of approximately <math>8*10^{-6} mol/dm^3</math>. This final value is dependent on the production rate (which has been estimated for all of the models). File:Model prea.bmp Results of the Matlab simulation of Model preA Model A
The initial concentration of split Dioxygenase, <math>c_0</math>, determines whether the system is amplifying. The minimum concentration for any amplification to happen is <math>10^{-5} mol/dm^3</math>. If the initial concentration of split Dioxygenase is higher, then the final concentration of Dioxygenase will be higher as well (see graphs below). Note that the obtained threshold value is higher than the maximum value that can be generated in the cell according to Model pre-A. File:Comparison a, prea.bmp Comparison between Model pre-A and Model A. Initial concentration of split Dioxygenase: <math>10^{-5} mol/dm^3</math>
<math>K_m</math> is indirectly proportional to the "final concentration" (which is the concentration at the end of the simulation), i.e. the bigger the value of <math>K_m</math>, the smaller the "final concentration" will be. Different <math>K_m</math> values determine how quickly the amplification will take place. (Also, it was found that the absolute value of <math>k_1</math> and <math>k_2</math> entered into Matlab does not change the outcome as long as the ratio between them (<math>K_m</math>≈<math>k_2/k_1</math>) is kept constant. This is important when simulating (in case entering very high values for <math>k_1</math> and <math>k_2</math> slows down the simulation).
If <math>k_{cat} = k_3</math> is increased, the model predicts that the dioxygenase concentration will rise quicker and to higher values. This indicates that <math>k_3</math> is the slowest step in the enzymatic reaction.
At the moment, our biggest source of error could be the production rate, which we could not obtain from literature. Hence, we had to estimate the value of the production rate (see variables). We hope to be able to take a measurement of this value in the lab as it has a big effect on model's behaviour.
We want to determine how our system reacts if different parameters are changed. This is to find out which parameters our system is very sensitive to.
Hence, the system is sensitive to most of the constants (given a particular range of values). The most crucial one, however, seems to be the initial concentration of split Dioxygenase. Model B
The initial concentration of split Dioxygenase, <math>c_0</math>, determines whether the system is amplifying. If the initial concentration is changed, the observed behaviour is similar to the one from Model A. If the initial concentration of split Dioxygenase is increased, then the final concentration of Dioxygenase will increase as well (see graphs below).
Running both models with the same initial conditions (<math>c_0=10^{-5} mol/dm^3</math>), it has been noted that Model B does not generate a siginificant amplfication over Model A. Hence, it would be more sensible to integrate a one step amplification module into our system. File:Comparison prea, a, b.bmp Comparison between Models pre-A, A and B Colour response model (25/08/2010)Initially, dioxygenase was being treated as "output" because the last reaction (catechol and dioxygenase react to form a coloured output) is common to all models. We concluded that 2-step amplification (Model B) presented little improvement over 1-step amplification (Model A). However, we realized that it is very important to add this last reaction, since it is an enzymatic reaction, hence an amplification step by itself. Adding it to all models means that Model A becomes a 2-step amplifier and Model B becomes a 3-step amplifier. Having previously drawn the conclusion that 2-step amplification is not much better than 1-step amplification, we decided to model this last step to deduce whether our initial conclusion was valid. If it was true, this would mean that our construct is not innovative at all. The important information about the last amplification step is that the coloured compound (i.e. the product of the last enzymatic reaction) is toxic to the cells. It is suggested that product of Catechol destroys the cell membrane by inhibiting lipid peroxidation. It causes significant changes in the structure and functioning of membrane components (e.g. disruption of membrane potential, removal of lipids and proteins, loss of magnesium and calcium ions). These effects cause the loss of membrane functions, leading to cell death. Since the product of Catechol acts on the cell membrane, it might not affect our enzymatic reaction immediately. In our simulation, we will try to model immediate cell death as well as neglecting the effect that the coloured output has on the cell. Comparing these two models will show if there are significant differences in the results.
Initial conclusions(26/08/2010)Despite our model not working entirely correctly (some negative concentrations were obtained), it was possible to deduce several points.
Problems with the simulationOnce we have implemented the colour change into the models, we noticed that there are some inconclusive results. After adding catechol, some concentrations were reaching negative values. We checked our equations and constants but could not find a mistake. Hence, we concluded that there must be something wrong with the way that MatLab evaluates the equation or deals with the numbers. The problem seems to originate from the very rapid concentration change of catechol which disrupts the whole system. Prospective solutions
We hoped that TinkerCell imposes non-negative conditions on its solutions. Hence, we implemented the whole amplification model (including coloured output) in TinkerCell. However, we realized that TinkerCell does not deal well with very high or low numbers (For example, values higher than <math>10^5</math> are not acceptable - this is important since our rate constants (<math>k_1</math>) are usually bigger than <math>10^5</math>. Also, the low degradation rates (<math>10^{-9}</math>) result in a zero output line). However, TinkerCell can still be used for testing that our Matlab programs behave the way we anticipated (by using default values of 1), as well as producing nice diagrams of our system.
We had a close look at the ODE solver options in MatLab. However, we were already using the one that produced the most reasonable results (ode15s). We found that decreasing Relative and Aboslute tolerances (to values as small as <math>10^{-15}</math>) significantly improved the simulation. However, this is not an ultimate solution as in the simulations negative numbers still appear (order of <math>10^{-15}</math>). We decided that such small negative concentrations were acceptable. We also decided that the point of interest lies between the first 100 to 150 seconds after adding catechol, while concentrations hit the negatve values at much higher time values. The images below show the influence of the relative and absolute tolerance values on the model. Note that it was important to allow the ode-olver to adjust the time step automatically, as big time steps (1 second) were generating wrong answers for the catechol model. Adjusting time steps manually to very small values was not efficient (the whole simulation does not require very high definition simluation).
We hoped that SimBiology could be more suited for our modelling than using ode-solvers, so we implemented our models into SimBiology. This package offers an interactive user interface similar to Tinker Cell, but uses MatLab to simulate. Initially, we confirmed that our simple production model (Model PreA) and 1-step amplification model (Model A) implemented in SimBiology generated exactly the same results as our ODE equation based models. The interface allowed us to have clearer control over paramters. It also allowed modelling special events, for example, adding catechol at certain point in time. Previously we had to split simulation into to two parts. References
Paramters & Constants
Constants for the Output Amplification Model
References
Results & Conclusions
Final Conclusions (31/08/2010)
If Catechol is added before t= 1000s, then the coloured output will reach its threshold value faster by simple production. If Catechol is added when t>1000s, then the coloured output will increase (marginally) faster through the amplification step in Model A. There does not seem to be a significant difference between the two models (Model preA and Model A). These observations are true for intial concentration of dioxygenase equal to <math>10^{-5}mol/dm^3</math>. However, we noticed that if the initial concentration is raised to <math>10^{-4}mol/dm^3</math>, then Model A can be more beneficial than Model preA after only 100 seconds. Hence, the question arises whether the concentration of protein in the cell can be as high as <math>10^{-4}mol/dm^3</math>. Our simple production model predicted that the concentration of protein could not reach such a high value. However, we decided to research more on ribosomal concentrations in bacteria to determine whether it is possible to establish such a high concentration in the cell. On the website E.coli Statistics [http://gchelpdesk.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi] it is stated that number of ribosomal proteins per cell is 900,000. In a cellular volume of order of <math>1\mu m^3 = 10^{-15}dm^3=10^{-15}L</math>, the above number of ribosomes converts to <math>1.5*10^{-3}mol/L</math>. This means that a concentration of <math>10^{-4}mol/dm^3</math> is not completely out of scale.
There seem to be 3 regions of catechol concentration that influence the system in different ways. These regions are: c>1M, 1M>c>0.01M, 0.01M>c. The boundaries of these regions tend to vary depending on the choice of other initial conditions. The values given above apply to boundary conditions that are currently considered to be physiologically relevant. Varying the initial concentration of catechol within the highest region does not result in any change of colour output response (It is possible that all enzymes are occupied and the solution is over saturated with catechol). In the middle region the catechol concentration influences the amplfication. Amplification decreased when the concentration tends towards 0.01M. When this region is entered, there is no difference in output production by the two models.
The coloured product of catechol kills cells by destroying the cell membrane. However, we do not know how quickly the cells will die. Therefore, we examined two different cases: immediate cell death and negligible cell death (i.e. cells death is negligible because it takes too long) Running the simluation in Matlab (not Simbiology!), our conclusions are:
Since it appears that the time of cell death is important, we decided to discuss this issue with Wolf and Harriet. Referring to this paper [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V78-4XXNV4Y-4-7&_cdi=5836&_user=217827&_pii=S0048969709011668&_orig=search&_coverDate=02%2F01%2F2010&_sk=995919994&view=c&wchp=dGLbVzz-zSkzk&_valck=1&md5=f8dca6227c29db659ddbeb588ad115e7&ie=/sdarticle.pdf (1)] we decided that cell death induced by catechool is a very slow process (we estimate that it will take a few hours) in comparison to the time scale that we are interested in (several seconds to minutes).
|