Team:ETHZ Basel/Modeling

From 2010.igem.org

Revision as of 14:46, 4 August 2010 by Georgerosenberger (Talk | contribs)

E-Lemming @ ETHZ – Modeling Overview

ETHZ’s project, E-Lemming, aims to modify the chemotaxis property of E.coli such that, instead of response to a chemical attractant/repellent, the bacterium responds to a light stimulus (phototaxis). Furthermore, this light sensitivity is used to control E.coli’s movement by deciding, at any given time, which type of motion will our ‘Lemming – in – disguise’ adopt (tumbling or straight run). This leads to a controllable E.coli, which can follow any pre – defined spatial path/maze/game, as a result of the combination of tumbling and going straight. The bacterium colony is imaged and, by image processing, the position of a single tracked cell is inferred. By activating a light switch, the user decides whether the bacterium should continue running or should change direction.

In the theory world, the steps we are following in mindlessly driving E.coli to our pre - defined target are the following:

• deterministic (ODE) & stochastic models of the chemotaxis pathway (documented from the literature)
• model of the movement of E.coli (built on the information of the pathway derived from the molecular models)
• control algorithms ( built on the user’s desire to play around with E.coli)
• image tracking & image processing algorithms
• java applications/movies of the E.Lemming (the fun part)

Model coupling

Coupled models for the simulation of the whole process and their interfaces.


Questions

Team:ETHZ_Basel/Modeling/Questions

A - Molecular Models

The network topology

Many bacteria posses a signal transduction network that makes it possible to sense the changes in concentration of a certain chemical attractant/repellent in the extracellular environment and direct the movement towards the attractant and away from the repellent. This property is what is known as ‘chemotaxis’. There are basically two types of movement that the bacterium can employ: 'tumbling', which means change of direction and occurs when the bacterium doesn’t sense an increase of attractant concentration in the extracellular environment anymore and 'running straight', which occurs when the attractant concentration is increasing. The running straight stands for “correct direction, ok, keep going” (towards the attractant/away from repellent) while the tumbling stands for “wrong direction, give it another shot”. In the case of E.coli, these two types of movements correspond to different rotation directions of its 4 flagellar motors (clockwise: tumbling; counter-clockwise: runs).
The direction of rotation of the motors can be determined by the ratio of concentrations of certain proteins inside the cell. That is, there are some quantitative indicators from which one can infer whether the bacterium will run or it will go straight, as a response to changes in input concentration. These indicators (concentration levels for certain proteins) are the result of a well – studied signal transduction network, consisting of membrane receptor proteins and intracellular proteins (Che).
Our models were constructed based on the following network topology:

<no worries, picture to be inserted>


CheW and CheA are localized right next to the membrane, inside the cell. They form a complex together with the receptor clusters on the membrane (which sense the changes in input concentration in the extracellular environment). This complex can be successively methylated (on different methylation states), a constant process done by CheR. In the literature, there are usually 5 methylation sites considered (m = 0,1,2,3,4; m=0 stands for no methylation and m=4 stands for the highest methylation level).
The key process that influences the outcome of the network (tumbling/running) is the autophosphorylation of the protein CheA, which is favoured either by the increase in the methylation state of the membrane complex, done by CheR, or by decreasing attractant concentration.
With decreased attractant (increased repellent) concentration, CheA is autophosphorylated. The phosphoryl groups are being used to produce either CheYp or CheBp. CheYp diffuses through the cytoplasm to the motors, generating clockwise spins, therefore tumbling. CheBp is responsible for de-methylation of the receptors, which reduces the chances of CheA being autophosphorylated, therefore the system is returning to its previous state.
With increased attractant (decreased repellent) concentration, the rate of CheA autophosphorylation is reduced, meaning that there are less phosphate groups available for CheYp and CheBp. Therefore the level of CheYp decreases and the level of CheY increases, which leads to straight runs. Then, the level of CheBp decreases, increasing the chances of CheA being autophosphorylated (since the receptors are methylated by CheR), therefore returning the system again to its previous state.

A perfect model of the chemotaxis pathway should reconstruct a system that exhibits the following characteristics:
adaptation: the ability to respond to changes in the external environment and return the intracellular protein phosphorylation levels to their pre-stimulus levels.
sensitivity: even small changes in the concentration of the attractant/repellent trigger chemotactic responses from bacteria.
gain: the ability to amplify the received signal from the extracellular environment, so as to modulate the intracellular signaling cascade.
robustness: perturbations of the network’s parameters should not influence the behavior of the system.

Even though the chemotaxis system has been long and in depth studied and it is well documented in the literature (from a modeling point of view), it is to be noted that none of the existing models, accounting for experimental findings, can satisfy all of the above requirments. So, we decided to implement several models to answer our questions.

Tasks for our models

1. Estimate CheYp for no tumbling: we already know from the literature that the tumbling frequency of the bacterium is a function of the concentration of CheY phosphorylated. Therefore, our first task was to determine the concentration of CheYp (defined as the ‘Threshold’) necessary to get a bias (= time not tumbling/total time) of 1 or close enough to 1 so that we could consider that the bacterium exhibits no tumbling (for a given input).
2. Construct the light input: Given the value of CheYp from 1), the second task would be to construct this light input which induces the effect of ‘no tumbling’ (or ‘running straight’). This is to be done by fusing PIF3 (for example) with any of the Che proteins (CheB, CheR, CheZ). The best light input signal is the one for which the bias is close enough to 1 for the highest amount of time (in other words, CheYp stays the most beyond the Threshold).
3. Repeat the previous tasks for different values of the input concentrations (different input values trigger different effects on the bias of the system).

ODE models of the chemotaxis pathway

Spiro (1997)

This model basically predicts that for all Che proteins (CheR, B, Y and Z) would make a good light input. For CheR and CheY the concentration of CheYp drops more than Delta (initial value - Threshold), for CheB and CheZ it increases more than Delta. For all Che proteins the concentrations stay like forever below/above the Threshold, until we deactivate them with far-red light. The best results was obtained for assuming a high ligand concentration (saturation, so that the methylation level of the receptors is high). For CheY and CheZ the reaction times were (as expected), much faster than for CheB and CheR. However, for these proteins the reaction times were still fast enough.

Mello (2003)

This model is already implemented compatible to the Spiro1997 model, but not yet documented.

Rao et al. (2004)

This model predicts quite a similar behavior of the Che proteins, as demonstrated by the previous model. A notable difference is that the response of CheYp level to the change of the level of the other Che proteins is robust to the fluctuations of the ligand (attractant) concentration. In the previous model (Spiro), the most desirable results were obtained only for higher concentrations of the ligand.

A detailed analysis of the simulation results obtained from this model can be found here.


Barkai & Leibler (1997)

Modeling of the light-input, light-output system

The light-input, light-output system PIF3-PhyB will be used to control spatial localization of key-proteins from chemotaxis. A model will be implemented according to a recent publication. In a next step, this model will be integrated with the existing chemotaxis models to predict photo-activation / -deactivation by the experimental setup.

Status: in development (subgoals reached)


Stochastic model

The stochastic features of the chemotaxis pathway seem to be quite strongly expressed, so we also looked into stochastic simulation models. The goal was to get a more realistic view over the chemotaxis process and complete the big picture obtained by implementing the above deterministic ODE models.

Status: in development.

B - Movement model

Stochastic model

C - Imaging

In order to control (or direct) the cell, we need to detect the current direction of its movement and the location, so that it could be communicated to the control logic that decides if the cell should be allowed to swim straight or make a turn.

The algorithm we implement will have to be fast enough to work on a live video stream so that there is sufficiently less processing delay for the control logic to work in real time. This is the main challenge of the cell tracking algorithm.

More information about the development of cell tracking and its current status can be found in the Cell Tracking page.

D - Controller