Team:KULeuven/Modelling

From 2009.igem.org


Introduction

As an introduction to modelling we made a short presentation. This presentation tells about the following:

  • some definitions and the role of modelling
  • black and white box modelling
  • the role of ordinary differential equations
  • modelling applied to iGEM

this presentation is mainly based on the presentation of last year and the wiki of ETH Zürich 2007.

The Full Model

Controller system.png

The full model consists of four parts: the blue light receptor, the comparator, the vanillin production and the vanillin receptor similar to the real world essencia coli bacteria. The blue light sensor and the vanillin receptor are acting as the inputs of the system. The intensity of blue light gives the wanted vanillin concentration. The vanillin receptor gives the actual concentration of vanillin outside the cell. Comparing these two inputs gives us a control system for the vanillin production.

Modelling steps

The scientific procedure for the study of a physical system can be divided in the following 3 steps

  • Parameterization of the system: discover the minimal set of parameters that completely define the system.
  • Foward modelling: define the physical laws that, given the values of the parameters of the system, determine the value of the observable parameters.
  • Inverse modelling: given observed parameters, infer the actual parameters that produced the observed data.

Most of the time this is an iterative approach, the last step will contain an indication of how good the model fits the observed data. Sometimes it will be necessary to adjust the proposed model.

Step 1 and 2: Modelling and simulation

By a model, we mean an abstraction of some real system that can be used to obtain predictions and formulate control strategies. In order to be useful, a model must necessarily incorporate elements of two conflicting attributes:

  • realism
  • simplicity

On the one hand, the model should serve as a reasonably close approximation to the real system, on the other hand, the model must not be overly complex so as to preclude its understanding and manipulation. All that's required is a high correlation between predictions and real-life performance.

We created a single cell model of interacting biochemical reacting to describe the behaviour of our 'Miss Blue Vanilla' using simbiology (Mathworks) graphical representation. This reactions can sufficiently accurate be described by Ordinary Differential Equations (ODEs) like Mass-Action laws, Hill Kinetic laws and so on. These equations are solved by a deterministic solver available in simbiology.

Step 3: Inverse modelling

Although some parameters can be obtained by mining the literature in a search of relevant information, some reactions are still not fully described or exactly known by science. This requires another method for obtaining the needed information, inverse modelling also known as parameter estimation.

While the first two steps are mainly deductive, this step is inductive. The inverse problem consists of using the actual result of some measurements to infer the value of the parameters that characterize the system. While the forward problem has (in deterministic physics) a unique solution, the inverse problem does not. The most general theory is obtained when using a probabilistic point of view, where the a priori information on the model parameters is represented by a probability distribution over the 'model space'. Because an exhaustive search through the model space is computationally very demanding, more intelligent Monte Carlo techniques will be used, we used a Metropolis-Hastings sampling algorithm. Notice that if one wants to resolve all the parameters in the model, a vast number of experimental data has to be obtained.

Most of the time we are only interested in the model that best fits the observed data. This best model can be obtained by the solution of a (non linear) optimization problem. This is the main method we will use to estimate the unknown parameters is the model, although the probabilistic approach is more general it suffers form major computational overhead.


References

"Inverse problem theory and methods for model parameter estimation", Albert Tarantola
"Modern simulation and modeling", Reuven Y. Rubinstein, Benjamin Melamed