Team:Edinburgh/modelling(reallifemodelling)

From 2009.igem.org

Revision as of 19:00, 14 October 2009 by Dzonni (Talk | contribs)

Modelling - Real Life Modelling
Personal note

Since our system designed to detect landmines it was thought that we should include a visual representation of how the system would operate in real life. It is known that TNT leaks out of the landmine into the soil and we predict that the bacteria will move toward the TNT due to the chemotaxis response of the ribose binding protein which has previously computationally designed TNT receptor1
TNT and Nitrite diffusion has been modelled using a finite difference method which approximates the 2 dimensional time dependent diffusion equation below



where α is the diffusion coefficient, to a numerical solution which can be solved . By splitting the area into a grid M2 by N with a space step of h=1/M (in both the x and y directions) and time step size of k=T/N. This gives a grid of the form



which was set up as a matrix in the program.

An initial concentration of TNT in the field was then set to 107800µg/cm3 and the diffusion of TNT out from this site was calculated from the following equation,


and



The diffusion coefficient for TNT is 1.18 *10-6 cm2 s-1 4 and the diffusion of Nitrite is 1.9*10-5 cm2 s-1 .

This system was coded in both MATLAB and JAVA. MATLAB was used to produce graphs and JAVA was used for the chemotaxis modelling which is described below.

The TNT diffusion graph after a year is shown below



Figure 1 shows the diffusion of TNT over a small area, only 10m by 10m, so that the uniform distribution of the TNT concentration can be shown. However a better representation of what would happen in the field can be shown when the size of the grid is much larger, 100m by 100m as shown below


It is clear from Figure 2 that TNT is only concentrated around the landmines location rather than across the whole field and hence the bacteria could be used to determine the approximate location of the landmine due to the light emit on detection of the TNT. The bacteria also produces enzymes, nitroreducases, which breaks down TNT into various products including Nitrites. Our system is designed to emit light when these Nitrites are detected. For the model it has been assumed that a very low conversion of TNT to Nitrites will be achieved. The diffusion graph for Nitrite is very similar to that for TNT.
As previously stated, it is expected that the bacteria will move toward the landmine due to chemotaxis of the TNT receptor. This was incorporated into the model by assuming the bacteria will move toward the highest concentration of TNT, essentially swimming up the concentration gradient. This was programmed by incorporating a 3rd matrix where the grid matches up with both the TNT and the Nitrite diffusion graphs. The bacteria is then placed in the matrix using a random number generator to get its grid location and this location is set to 1 while the rest of the grid is set to zero.

The model then calculates the concentration of TNT at its grid position and its 4 surrounding grid positions. It then stores the current concentration as the maximum concentration and then compares each of the surrounding locations concentrations to this maximum. If any of the concentrations are greater than the current locations concentration then it set this concentration as the maximum and will then compare subsequent concentrations to this new maximum. It also indexes the location of the maximum concentration which is used to determine where the bacteria should move to. The moving of the bacteria is done by a series of IF statements as shown below

if (index==1){
Bact[xbact][ybact]=1;
}

if (index==2){
Bact[xbact-1][ybact]=1;
xbact=xbact-1;
}

if (index==3){
Bact[xbact+1][ybact]=1;
xbact=xbact+1;
}

if (index==4){
Bact[xbact][ybact-1]=1;
ybact=ybact-1;
}

if (index==5){
Bact[xbact][ybact+1]=1;
ybact=ybact+1;
}



Depending on the value of the index, depend on where the bacteria moves to by setting the new location in the grid for the bacteria to 1 and by reseting the location of the bacteria.

The random motion of chemotaxis has not been included in the model as over the distances that system operates it was thought that this would be unimportant and in general it would travel in a straight line.

This model is shown in an applet where the background is coloured green to represent a field, the landmine is coloured black and the bacteria is coloured yellow.
The model could be developed further by introducing more bacteria which would all move towards the highest concentration of TNT. Further to this the model could incorporate the concentration of Nitrite and dependant on the concentration of TNT and Nitrite the bacteria could be shown to represent the colour that would experienced in field, ie there would be a yellow and blue circle where the bacteria are. Further to this, the model could be linked with the gene expression modelling so that the colour of the light emitted could be dependent on the level of each protein expressed at each location.

It has also been suggested that the model could be made into a game, similar to minesweeper, where you use the movement of the bacteria to determine where the landmine is located.
display ('***********************************************************')
%Finite difference model for the diffusion of TNT through soil

xsize=100;
ysize=100;
MaxSteps=8760;
TimeStep=0.01;
SpaceStep=1;
diffusivityC=1.18*(10-6)*3600/100/100;
%calculates s and s2 for TNT based on the diffusion coefficient of TNT
sC=(diffusivityC*TimeStep/SpaceStep/SpaceStep);
s2C=(1-4*sC);
diffusivityN=1.18*(10-6)*3600/100/100;
%Calculates s and s2 for Nitrites based on the diffusion coefficient of
%Nitrites
sN=(diffusivityN*TimeStep/SpaceStep/SpaceStep);
s2N=(1-4*sC);


% Sets up 2 matricies for the calculation
C=zeros (xsize,ysize);
CNew = zeros (xsize,ysize);
N= zeros (xsize,ysize);
NNew= zeros (xsize,ysize);


%Set up the concentration and position of the landmine
MineX=5;
MineY=5;
C(MineX,MineY)=107800;
N (MineX,MineY)=1;

React=0.99;
Produced=1/React;

for k=1:MaxSteps,
for i=2:xsize-1,
for j=2 :ysize-1,
%Calculates concentrations due to finite difference equation
CNew(i,j)= React*(sC*(C((i-1),j)+ C((i+1),j)+C(i,(j-1))+C(i,(j+1)))+s2C*C(i,j));
NNew(i,j)=Produced*(sN*(N((i-1),j)+ N((i+1),j)+N(i,(j-1))+N(i,(j+1)))+s2N*N(i,j));

end
end

C=CNew;
N=NNew;
display (C)
display(N)
>br /> surfc (C)

refreshdata

end
Edinburgh University iGEM Team 2009