Team:Calgary/Modelling/Method

From 2009.igem.org

Revision as of 18:38, 20 October 2009 by Cmukerji (Talk | contribs)

University of Calgary

UNIVERSITY OF CALGARY



MODELLING INDEX
Overview

Membrane Computing Modelling
Differential Equation Modelling

A TOUR OF THE UNIVERSITY OF CALGARY iGEM TEAM


We've reached modelling, the fifth stop on our tour! We've looked in to two different methods of modelling our system: Differential Equation Based Modelling and Membrane Computing. Here, you can explore the similarities and differences, as well as the functions of each method. As well, you can find the results of our characterization of the signalling pathway. Once you're done, we'll move on to the Second Life component of the project HERE.


DIFFERENTIAL EQUATIONS MODELLING METHODS
The simbiology interface from Matlab was used to simulate the differential equations model. Chemical Kinetic equations were used to build the model for simulation.

Fig : The Reaction of Species A with B to produce C and D


Fig : The Chemical Kinetic Rate Equation
k is the kinetic rate constant. The size of k will determine the speed of the reaction. A smaller value of k will produce a slow reaction rate while a larger value of k will produce a fast reaction rate.
[A] is the amount of reactant A present.
The simulations were run for 50000 seconds . It was considered to be enough time for the system to reach equilibrium after disturbance. The Sundials Solver (how do they work ?)was used to run these simulations because for this model it was considered to produce optimal results. (how do simulations work?)

The Reactions
The system was represented by the following reactions. The reactions with double headed arrows have two rate constants(forward/ reverse rate constant). All reactions were assumed to be elementary reactions.

Parameter Rationale
Table: Initial Values of the Species in the System

Species Initial Value Rationale
AI-2 0 Initially the amount of AI-2 is constant at 0. After an equilibruim is established variable amounts of AI-2 are added at different simulations.
LuxPQ 10 The amount of LuxPQ varies depending on the simulation run.
AI-2:LuxPQ 0 This value is kept at time = 0 because the initial concentration of AI-2 is 0.
LuxU:p 2 ----
LuxU 1000 There is a lot of this species present in the cell in nature. To signify plenty we assign it a value 1000.
LuxO:p 2 Equal amounts of LuxO:p and LuxU:p was considered in the model because LuxU:p phosphorylates LuxO . The phosphorylation reaction is considered to be a fast reaction therefore there are equal amounts of the two proteins present.
LuxO 1000 There are equal amounts of LuxO and LuxU present since they are encoded within the same operon in the cell.
p 10.0658 We assume that there is enough p is the environment that it doesn’t become a limiting factor. For that reason we assign p as a constant value in simbiology. (It doesn’t really matter that the initial amount is presented as a comparatively small number in this case. )
sigma54 0.14183 This amount is kept at a constant value to ensure that this value does not become a limiting factor.
sigma54:LuxO:p:Pqrr4 0.63 There is only 1 copy of Pqrr4 present in each cell . since in the reaction equations Pqrr4 is shared between 2 other equations we decided to break the concentration of Pqrr4 between 3 species: sigma54:LuxO:p:Pqrr4 , Pqrr4 , sigma54:Pqrr4 . The initial values of the three species add up to one. The fractions of the Pqrr4 combination species are weighted differently . Since the Pqrr4 promotor stays on most of the time we decided the sigma54:LuxO:p:Pqrr4 complex should recieve the most weight. Pqrr4 is assumed to stay unbound from any complex for the least amount of time therefore Pqrr4 initial amount is the smallest.
Sigma54:Pqrr4 0.345
Pqrr4 0.025
GFP    
mRNA    
Table: The Kinetic Rate Constant Values
     
     
     

Retrieved from "http://2009.igem.org/Team:Calgary/Modelling/Method"