Team:TU-Eindhoven/Modeling

From 2013.igem.org

...no models can be 'correct'. The role of any model is to provide
a framework for viewing known facts and to suggest experiments.

Suresh Moolgavkar

Modeling Summary

The main objective of the modeling component of our project is to provide insight to the rest of team in what concerns the effect of altering the necessary conditions for the proper operation of each component of our device. Keeping this in mind we established a series of questions that if answered thoroughly would benefit our project. We categorized these questions in four main topics:

TU-Eindhoven Images eColiModelsDiagram.jpg
eColiModelsFigure General Models of Processes Happening Inside the E. Coli.

Protein Selection

  1. Which proteins that can be expressed in bacteria are best suited as CEST based markers?


Gene expression

  1. How does the oxygen concentration affect the concentration of the transcription activator (FNR) in its active form?
  2. How does introducing decoy sites influence the gene expression rate?


Tumor Imaging In Vivo

  1. How do bacteria distribute within the human body?
  2. How long does it take for the bacteria to reach the target area?
  3. How do the bacteria react to the activation of the killing mechanism?
  4. What is the bacteria elimination rate?


Based on the previous question we developed seven models that can be divided into two categories, processes occurring inside the bacteria and the interaction of the bacteria with its surroundings. presents a schematic diagram of the developed models related to the first category, while depicts the model that accounts for the interaction between the bacteria and its surroundings when introduced into the human body for tumor imaging. Both figures will be further explained in the next paragraphs. The approaches used for the models range from molecular dynamics simulations to deterministic models and stochastic simulations as well as pharmacokinetics.


As previously mentioned, aims to schematically depict all the models that are related to the diverse components of our device, e.g the proteins that will be produced by the device, the FNR promoter and decoy sites which are related to the gene expression, and the CEST protein production directly related to the CEST MRI contrast agent.

TU-Eindhoven Images PBPKModelDiagram.jpg
PBPKFigure Model That Depicts The Administration Of The CEST MRI Contrast Producing Bacteria Into The Patient's Bloodstream.


To define which proteins are most suitable for the production of the MRI CEST contrast agent by the bacteria we used molecular dynamics simulations. By using this approach we were able to determine the list of proteins that would be generated by our device, based on the the accessibility of the various exchangeable hydrogen atoms of the Arginine and Lysine amino acids of the proteins being considered. In what concerns the FNR promotor and decoy sites, both deterministic models and stochastic simulations were used to predict the concentration of transcription factor and influence of introducing decoy sites. Aside from the individual models made for these two components of the device, a merged deterministic model was created to analyze the complete gene expression process. A sensitivity analysis was used on the merged model to determine which elements of the model have the greatest influence on the outcome. Similarly, a PLE analysis was utilized on the decoy sites model.


Furthermore, a deterministic model was developed for analyzing the protein production and degradation. Using this information it was possible to predict the average concentration of the CEST proteins produced by our device. These results where then compared to the experimental data.


is related to the proposed applications of our contrast agent. This figure presents the complete process of injecting the bacteria into the bloodstream in order to reach the tumor region and start producing the MRI CEST contrast agent. When the contrast agent is no longer needed a prodrug is injected into the bloodstream to activate the killing mechanism that will induce apoptosis in the bacteria. The whole process is modeled using a pharmacokinetics approach. By using this approach we aim at predicting the distribution of the bacteria in the body and model the kill switch that should be implemented when further developing this application.

The source code of all models can be found here.