Team:NTNU Trondheim/Modeling

From 2011.igem.org

Revision as of 15:15, 18 July 2011 by Lutro (Talk | contribs)




Modeling

To describe and understand the biological reactions and process, as the bacteria turns red under stress we developed multiple mathematical and statistical models. As the most basic model we used a system om ordinary differential equations (ODE), this is a fully deterministic model, describing the change in concentration for all molecules involved. The process for a bacteria to go from normal to glowing red is involved, and one can think of this as a stochastic process, that at each step, the process can ether succeed or fail, with a given probability. This gives rise to a Bayesian model.

Combining these two models will result in a system of stochastic differential equations (SDE), which will be solved using numerical algorithms. As a last model we will explore the relationship between variants of stress and fluoricene intensity using (non)-conventional regression.

Model Introduction

At the heart of the modeling lies biological consistency and data integration. The modeling will be focused on interpretation simplicity and data consistency. That is to develop models that can be easily interpreted by biologist and mathematicians, but the models should also strive to describe that which is observed at the laboratory.

The two main ways to model biological systems. One way is deterministic, the other is stochastic. In this project we will attempt to approach to problems in both ways. Using a deterministic model, with fixed parameters, and a stochastic model to integrate data more dynamically.

The Models

Four basic models where constructed, which will be described below.

Systems of ODE

This model is a fully deterministic model, describing the change of concentration of multiple molecules in a system over time. As time goes by, this will describe the evolution of the systems and the stability of the dynamics. At the basis of the model is a set of ODE's, describing the change in concentration over time in relation to the other substances present in the bacteria.


Using ode45() in matlab, solving the system numerically.

Bayesian Hierarchy

As shown in the introduction, the process before the bacteria turns red is a linear chain of events. This can be modeled as a stochastic process, and solved with Monte Carlo methods. This is a way of including stochastic elements and data integration. Setting the probability of a transition to be successful to be Beta distributed, with hyper parameters ("known constants"), αi and βi, the the probability of total success, that the start input will transcend and cause the output (observable), is the product of all the beta distributed variables, which itself is Beta distributed whit parameter α = ∑ αi-q+1 and β = ∑ βi-q+1, where q is the number of parameters.

Having done some previous test at the lab, we wish to extend this in a probabilistic context. To integrate this data and the model assumptions, we can get a better picture of what is happening in the cell, and how likely it is to succeed in each trail.

Stochastic Differential equations

This is AWESOME

Regression Models

Model Validation

No data yet.

References