Team:NTNU Trondheim/Journal Modeling

From 2011.igem.org

(Difference between revisions)
(Friday 1/7)
(Journal for Modeling Group)
 
(55 intermediate revisions not shown)
Line 2: Line 2:
 +
= Journal for Modeling Group =
-
=== Journal for Modeling Group ===
+
{{:Team:NTNU_Trondheim/Journal_Modeling/Calendar}}
 +
You can find an overview of the modeling [[Team:NTNU_Trondheim/Modeling|here]].
-
{{:Team:NTNU_Trondheim/Journal_Modeling/Calendar}}
 
-
You can find an overview of the modeling[[Team:NTNU_Trondheim/Modeling| here]].
 
== Wednesday 22/6 ==
== Wednesday 22/6 ==
-
Introductory lectures in systems biology was held by Marius Eidsaa. Different approaches were discussed and it was decided to work with two different models: Bayesian Networks and systems of Ordinary Differential Equations, and possibly extend it to Stochastic Differential Equations.
+
 
 +
Introductory lectures in systems biology was held by [[Team:NTNU_Trondheim/Team#Marius_Eidsaa|Marius Eidsaa]]. Different approaches were discussed and it was decided to work with two different models: Bayesian Networks and systems of Ordinary Differential Equations, and possibly extend it to Stochastic Differential Equations.
 +
 
==Friday 24/6==
==Friday 24/6==
 +
Dizzy was installed and we started to look at how to simulate Michealis Menten reactions.
Dizzy was installed and we started to look at how to simulate Michealis Menten reactions.
 +
==Monday 27/6==
==Monday 27/6==
 +
We continued where we ended last week, with having lectures and to get aquinted with different software.
We continued where we ended last week, with having lectures and to get aquinted with different software.
 +
== Thursday 30/6 ==
== Thursday 30/6 ==
 +
A program which solves ordinary differential equations numerically using the Runge Kutta method of order 4 and 5
A program which solves ordinary differential equations numerically using the Runge Kutta method of order 4 and 5
-
was made in Matlab. Using equation (*) and (**) (see the )  and having the level of ppGpp fixed the reponse of the system could be found.
+
was made in Matlab. Using expressions based on proteins inhibiting mRNA and having the level of ppGpp fixed the reponse of the system could be found.
 +
 
 +
Tried to find suitable parameters, but found this to be hard. Contacted the wetlab group and got help to find the  processes which cause inhibation. This can be used for modeling in Dizzy.  
The Bayesian model started to take form, and is ready to be implemented in C++.
The Bayesian model started to take form, and is ready to be implemented in C++.
 +
==Friday 1/7 ==
==Friday 1/7 ==
-
[[File:PpGpp.png|thumb|The level of LacI decreases when ppGpp is introduced and leads to an increase of mCherry]]
 
-
The program was tested to see if it gave reasonable results. Minor errors were corrected and the first results were obtained, due to a lack of accurate parameters the units on the axis are arbitrary. Here is the Matlab code:
 
 +
[[File:PpGpp.png|thumb|The level of LacI decreases when ppGpp is introduced and leads to an increase of mCherry]]
 +
The Matlab program was tested to see if it gave reasonable results. Minor errors were corrected and the first results were obtained, due to a lack of accurate parameters the units on the axis are arbitrary. Here is the Matlab code:
tspan=[0;20]; % defines the start and stop <br />
tspan=[0;20]; % defines the start and stop <br />
Line 59: Line 69:
The resulting plot is shown in this [[Media:PpGpp.png|figure]].
The resulting plot is shown in this [[Media:PpGpp.png|figure]].
 +
== Monday 4/7 ==
== Monday 4/7 ==
 +
[[File:MCherry.png|thumb|Response of mCherry to different values of ppGpp]]
[[File:MCherry.png|thumb|Response of mCherry to different values of ppGpp]]
The program was modified to make it possible to have several values of ppGpp in one run making it possible to see the response to different values of ppGpp, see [[Media:MCherry.png|figure]].
The program was modified to make it possible to have several values of ppGpp in one run making it possible to see the response to different values of ppGpp, see [[Media:MCherry.png|figure]].
The theory for the Bayesian model seems to be in place, see [[Media:Bayesian.png|figure]] [[File:Bayesian.png|thumb]], some of the implementation has already been done in C++.
The theory for the Bayesian model seems to be in place, see [[Media:Bayesian.png|figure]] [[File:Bayesian.png|thumb]], some of the implementation has already been done in C++.
 +
 +
 +
== Tuesday 5/7 ==
 +
 +
Started to implement the model in Dizzy and looking at the relationship between this model and the model implemented in Matlab. The model in Dizzy involves many more processes, if these are assumed to be fast the species involved can be assumed to be in equilibrium and be replaced with a steady state value. The current model in Dizzy involves 10 processes, while the Matlab model only have four.
 +
== Thursday 7/7 ==
== Thursday 7/7 ==
 +
[[File:step.png|thumb]]
[[File:step.png|thumb]]
-
A new program was made in Matlab to treat the concentration of ppGpp as an input signal. A possible feedback effect was also included as production of mCherry might lead to increase in ppGpp. The concentration of ppGpp was split into two parts, one external controlled and one caused by feedback , see [[Media:step.png|figure]].
+
A new program was made in Matlab to treat the concentration of ppGpp as an input signal. A possible feedback effect was also included as production of mCherry might lead to increase in ppGpp. The concentration of ppGpp was split into two parts, one external controlled and one caused by feedback, see [[Media:step.png|figure]].
 +
 
== Friday 8/7 ==
== Friday 8/7 ==
-
We was given an introductory lecture in Biology and Biotechnology by  [[Team:NTNU_Trondheim/Team#Gunvor_Røkke|Gunnvor ]], to get an understanding of the work beeing done at the lab.
+
 
 +
Derived expressions for the steady state concentrations and compared it to the results obtained in Matlab.
 +
 
 +
We were given an introductory lecture in molecular biology and biotechnology by  [[Team:NTNU_Trondheim/Team#Gunvor_Røkke|Gunvor]], to get an understanding of the work beeing done at the lab.
 +
 
== Monday 11/7 ==
== Monday 11/7 ==
 +
After a team meeting we was given a Master Thesis by Sigurd Johansen describing a toogle switch. We started to expand the ODE model to include more reactions than we have done so far. Both to be aware of possible limitations of the current model and to get a clearer understanding of the fundamental processes involved.
After a team meeting we was given a Master Thesis by Sigurd Johansen describing a toogle switch. We started to expand the ODE model to include more reactions than we have done so far. Both to be aware of possible limitations of the current model and to get a clearer understanding of the fundamental processes involved.
Started to investigate how to model and solve stochastic differential equations, and how to derive them from physical properties.
Started to investigate how to model and solve stochastic differential equations, and how to derive them from physical properties.
 +
 +
 +
==Wednesday 13/7==
 +
 +
The Master Thesis includes a lot more step then currently included in our model and from these step a set of four equations are obtained. Processes currently not considered in our model includes: Several ways of inhibation (i.e both one or two ppGpp kolecules might inactivate the rrnB promoter) dimerization of the Lac repressor, several steps during transcription.
 +
== Thursday 14/7 ==
== Thursday 14/7 ==
 +
We contacted the lab group to clarify the way ppGpp inhibits the rrnB promotor and was informed that it affects the RNAp and not the promotor in itself and can therefore  also affect pLac.
We contacted the lab group to clarify the way ppGpp inhibits the rrnB promotor and was informed that it affects the RNAp and not the promotor in itself and can therefore  also affect pLac.
The theory behind stochastic differential equation was investigated further.
The theory behind stochastic differential equation was investigated further.
 +
== Wednesday 20/7 ==
== Wednesday 20/7 ==
 +
We contacted Autodesk to explore the possibility of using Maya and the mMaya toolkit to visiualize the system.
We contacted Autodesk to explore the possibility of using Maya and the mMaya toolkit to visiualize the system.
We also visited the lab to see how gelelectrophoresis was performed.
We also visited the lab to see how gelelectrophoresis was performed.
 +
==Thursday 21/7 ==
==Thursday 21/7 ==
 +
We recieved mail from Autodesk, installed mMaya and explored the possibility of creating animations.
We recieved mail from Autodesk, installed mMaya and explored the possibility of creating animations.
 +
== Friday 22/7 ==
== Friday 22/7 ==
 +
A simple animation of a green bacteria (or something resembling that) turning red was made in Maya.
A simple animation of a green bacteria (or something resembling that) turning red was made in Maya.
 +
 +
[[Media:bacteria_anim.ogg|Simple animation]]
 +
==Wednesday 27/7==
==Wednesday 27/7==
 +
Meeting with prof. Almaas.
Meeting with prof. Almaas.
   
   
-
We decided to  drop the bayesian model [[Team:NTNU_Trondheim/Modeling#Bayesian_Hierarchy|see here]], and focus on stochastic models.
+
We decided to  drop the bayesian model ([[Team:NTNU_Trondheim/Modeling#Bayesian_Hierarchy|see here]]), and focus on stochastic models.
 +
 
==Thursday 28/7 ==
==Thursday 28/7 ==
 +
Dizzi was installed on one of the calculations server at NTNU, Dept. of Mathematics.  
Dizzi was installed on one of the calculations server at NTNU, Dept. of Mathematics.  
An electronic mail was send to see if another iGEM team needed help (requirement for gold medal).  
An electronic mail was send to see if another iGEM team needed help (requirement for gold medal).  
 +
==Friday 29/7 ==
==Friday 29/7 ==
 +
Looked into regression models for stress response, using non-linear regression with heteroskedasticity (from mass spectrometry??). Think that Kernel smoothing estimation might be the best approach.       
Looked into regression models for stress response, using non-linear regression with heteroskedasticity (from mass spectrometry??). Think that Kernel smoothing estimation might be the best approach.       
 +
==Saturday 30/7 ==
==Saturday 30/7 ==
 +
Kernel smoothing might be a bad idea, found out how to do smoothing B - splines with unequal weighting in R, which is just what we need. The weight should be inverse proporsjonal with the variance at each point.
Kernel smoothing might be a bad idea, found out how to do smoothing B - splines with unequal weighting in R, which is just what we need. The weight should be inverse proporsjonal with the variance at each point.
 +
== Monday 1/8 ==
 +
 +
Sozi was installed in the computer lab, and we looked into the presentation and poster. How to make that one great? Using Inkscape..
 +
Looked at some of the poster/presentation from last year, and found out what NOT to do, but not what to do.
 +
 +
Searched for reasonable parameters for the dizzy simulation.
 +
 +
 +
== Tuesday 2/8 ==
 +
 +
Found out how to use loops in Dizzy, which makes it easy to vary parameters and start condition. Experimented with storing data from dizzy and made a script in matlab which plots the data.
 +
 +
 +
== Wednesday 3/8 ==
 +
 +
[[File:rnap_var1.png|thumb|&gamma; denotes how strong ppGpp affects mCherry, if &gamma;=1.0 it is not affected at all and if &gamma;=0 it is affected just as strong as LacI ]]
 +
Began to take advantage of the features developed yesterday, started to study what effect a possible inhibation of mCherry by ppGpp could have, see [[Media:rnap_var1.png|figure]].
 +
 +
Made a short presentation in sozi, to show how it should be done in Amsterdam
 +
 +
 +
== Thursday 4/8 ==
 +
 +
Met with prof. Eivind Almaas, to discuss further progress, started to extend the tau-leap program with complex prameters.
 +
 +
 +
== Friday 5/8 ==
 +
 +
Found a reasonable way of simulating ppGpp in the cell as a "restricted" random walk, using a mixture of a uniform and a exponential random variable.
 +
 +
[[File:randomWalk_ppGpp.png|thumb|Shows how ppGpp might change over time under stressful conditions, negative values on y-axis is irrelevant]]
 +
 +
Looked into how to plot the output from the C++ program on tau leap, decide that python might be best (but never used it before)
 +
 +
 +
== Saturday 6/8 ==
 +
 +
Writing on the model page... and some basic python trail; "Hello World!"
 +
 +
 +
== Monday 8/8 ==
 +
 +
Threw the old poster overboard, and started on a new one.
 +
 +
 +
== Tuesday 9/8 ==
 +
 +
Used Dizzi on the calculation server, produced weird results, no idea why. 
 +
 +
 +
== Wednesday 10/8 ==
 +
 +
Continued to work on the new poster
 +
== Thursday 11/8 ==
 +
Tried to learn python, and make a script that plots the output from the tau - leap program. It doesn't work at all, who says python is easy?!?
 +
== Friday 12/8 ==
 +
Parameter testing..
{{:Team:NTNU_Trondheim/NTNU_footer}}
{{:Team:NTNU_Trondheim/NTNU_footer}}

Latest revision as of 15:24, 20 September 2011



Journal for Modeling Group

June
Mon Tue Wed Thu Fri Sat Sun
1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30
July
Mon Tue Wed Thu Fri Sat Sun
1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30 31
August
Mon Tue Wed Thu Fri Sat Sun
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31

You can find an overview of the modeling here.


Wednesday 22/6

Introductory lectures in systems biology was held by Marius Eidsaa. Different approaches were discussed and it was decided to work with two different models: Bayesian Networks and systems of Ordinary Differential Equations, and possibly extend it to Stochastic Differential Equations.


Friday 24/6

Dizzy was installed and we started to look at how to simulate Michealis Menten reactions.


Monday 27/6

We continued where we ended last week, with having lectures and to get aquinted with different software.


Thursday 30/6

A program which solves ordinary differential equations numerically using the Runge Kutta method of order 4 and 5 was made in Matlab. Using expressions based on proteins inhibiting mRNA and having the level of ppGpp fixed the reponse of the system could be found.

Tried to find suitable parameters, but found this to be hard. Contacted the wetlab group and got help to find the processes which cause inhibation. This can be used for modeling in Dizzy.

The Bayesian model started to take form, and is ready to be implemented in C++.


Friday 1/7

The level of LacI decreases when ppGpp is introduced and leads to an increase of mCherry

The Matlab program was tested to see if it gave reasonable results. Minor errors were corrected and the first results were obtained, due to a lack of accurate parameters the units on the axis are arbitrary. Here is the Matlab code:

tspan=[0;20]; % defines the start and stop
x=10.0; %the initial conentration of ppGpp
%the initial concentrations of LacI mRNA, LacI, mCherry mRNA and mCherry
y0=[20;20;0;0];
[T Y]=ode45(@F,tspan,y0);
plot(T,Y(:,2),T,Y(:,4),'r')
ppGpp=num2str(x);
title(['ppGpp=',ppGpp])
legend('LacI','mCherry','Location','East')
xlabel('t(a.u)')
ylabel('Concentration (a.u)')

Where F is a function defined by:

function Y = F(t,y)
Y=zeros(4,1);
x=10.0;
k1=[5e-4;5e-4];
k_prime=[20.0;20.0];
d1=[1;1];
n=[1;1];
k2=[20;20];
d2=[20;20];
Y(1)=k1(1)+k_prime(1)/(1+x^n(1))-d1(1)*y(1);
Y(2)=k2(1)*y(1)-d2(1)*y(2);
Y(3)=k1(2)+k_prime(2)/(1+y(2)^n(2))-d1(2)*y(3);
Y(4)=k2(2)*y(3)-d2(2)*y(4);
end

The resulting plot is shown in this figure.


Monday 4/7

Response of mCherry to different values of ppGpp

The program was modified to make it possible to have several values of ppGpp in one run making it possible to see the response to different values of ppGpp, see figure.

The theory for the Bayesian model seems to be in place, see figure
Bayesian.png
, some of the implementation has already been done in C++.


Tuesday 5/7

Started to implement the model in Dizzy and looking at the relationship between this model and the model implemented in Matlab. The model in Dizzy involves many more processes, if these are assumed to be fast the species involved can be assumed to be in equilibrium and be replaced with a steady state value. The current model in Dizzy involves 10 processes, while the Matlab model only have four.


Thursday 7/7

Step.png

A new program was made in Matlab to treat the concentration of ppGpp as an input signal. A possible feedback effect was also included as production of mCherry might lead to increase in ppGpp. The concentration of ppGpp was split into two parts, one external controlled and one caused by feedback, see figure.


Friday 8/7

Derived expressions for the steady state concentrations and compared it to the results obtained in Matlab.

We were given an introductory lecture in molecular biology and biotechnology by Gunvor, to get an understanding of the work beeing done at the lab.


Monday 11/7

After a team meeting we was given a Master Thesis by Sigurd Johansen describing a toogle switch. We started to expand the ODE model to include more reactions than we have done so far. Both to be aware of possible limitations of the current model and to get a clearer understanding of the fundamental processes involved.

Started to investigate how to model and solve stochastic differential equations, and how to derive them from physical properties.


Wednesday 13/7

The Master Thesis includes a lot more step then currently included in our model and from these step a set of four equations are obtained. Processes currently not considered in our model includes: Several ways of inhibation (i.e both one or two ppGpp kolecules might inactivate the rrnB promoter) dimerization of the Lac repressor, several steps during transcription.


Thursday 14/7

We contacted the lab group to clarify the way ppGpp inhibits the rrnB promotor and was informed that it affects the RNAp and not the promotor in itself and can therefore also affect pLac.

The theory behind stochastic differential equation was investigated further.


Wednesday 20/7

We contacted Autodesk to explore the possibility of using Maya and the mMaya toolkit to visiualize the system. We also visited the lab to see how gelelectrophoresis was performed.


Thursday 21/7

We recieved mail from Autodesk, installed mMaya and explored the possibility of creating animations.


Friday 22/7

A simple animation of a green bacteria (or something resembling that) turning red was made in Maya.

Simple animation


Wednesday 27/7

Meeting with prof. Almaas.

We decided to drop the bayesian model (see here), and focus on stochastic models.


Thursday 28/7

Dizzi was installed on one of the calculations server at NTNU, Dept. of Mathematics.

An electronic mail was send to see if another iGEM team needed help (requirement for gold medal).


Friday 29/7

Looked into regression models for stress response, using non-linear regression with heteroskedasticity (from mass spectrometry??). Think that Kernel smoothing estimation might be the best approach.


Saturday 30/7

Kernel smoothing might be a bad idea, found out how to do smoothing B - splines with unequal weighting in R, which is just what we need. The weight should be inverse proporsjonal with the variance at each point.


Monday 1/8

Sozi was installed in the computer lab, and we looked into the presentation and poster. How to make that one great? Using Inkscape.. Looked at some of the poster/presentation from last year, and found out what NOT to do, but not what to do.

Searched for reasonable parameters for the dizzy simulation.


Tuesday 2/8

Found out how to use loops in Dizzy, which makes it easy to vary parameters and start condition. Experimented with storing data from dizzy and made a script in matlab which plots the data.


Wednesday 3/8

γ denotes how strong ppGpp affects mCherry, if γ=1.0 it is not affected at all and if γ=0 it is affected just as strong as LacI

Began to take advantage of the features developed yesterday, started to study what effect a possible inhibation of mCherry by ppGpp could have, see figure.

Made a short presentation in sozi, to show how it should be done in Amsterdam


Thursday 4/8

Met with prof. Eivind Almaas, to discuss further progress, started to extend the tau-leap program with complex prameters.


Friday 5/8

Found a reasonable way of simulating ppGpp in the cell as a "restricted" random walk, using a mixture of a uniform and a exponential random variable.

Shows how ppGpp might change over time under stressful conditions, negative values on y-axis is irrelevant

Looked into how to plot the output from the C++ program on tau leap, decide that python might be best (but never used it before)


Saturday 6/8

Writing on the model page... and some basic python trail; "Hello World!"


Monday 8/8

Threw the old poster overboard, and started on a new one.


Tuesday 9/8

Used Dizzi on the calculation server, produced weird results, no idea why.


Wednesday 10/8

Continued to work on the new poster


Thursday 11/8

Tried to learn python, and make a script that plots the output from the tau - leap program. It doesn't work at all, who says python is easy?!?


Friday 12/8

Parameter testing..