Team:Paris Bettencourt/Modeling/Diffusion
From 2011.igem.org
Line 91: | Line 91: | ||
<br/> | <br/> | ||
<h3>Maya modeling</h3> | <h3>Maya modeling</h3> | ||
- | <p>this model is mapped in Maya for a user friendly aspect.</p> | + | <p>this model is mapped in Maya for a user friendly aspect. |
+ | This graphic reprsentaton is just for an idea of diffusion and it doesn't use the diffusion equation used by the java program</p> | ||
<iframe width="425" height="349" src="http://www.youtube.com/embed/jrX1Wa-CgCc?hl=fr&fs=1" frameborder="0" allowfullscreen></iframe> | <iframe width="425" height="349" src="http://www.youtube.com/embed/jrX1Wa-CgCc?hl=fr&fs=1" frameborder="0" allowfullscreen></iframe> | ||
<h2>Results</h2> | <h2>Results</h2> | ||
- | <p>results extracted from java output</p> | + | <p>results extracted from java output for 100 molecules and 10 row. the average time for one molecule is calculated from 10 molecule passing fromm cell1 to cell2 </p> |
+ | <table> | ||
+ | <tr> | ||
+ | <td>Molecules name</td> | ||
+ | <td>T7</td> | ||
+ | <td>tRNA</td> | ||
+ | <td>Insulin</td> | ||
+ | <td>GFP</td> | ||
+ | <td>Glucose</td> | ||
+ | </tr><tr> | ||
+ | <td>number of row</td> | ||
+ | <td>37870.95</td> | ||
+ | <td>30980.83</td> | ||
+ | <td></td> | ||
+ | <td></td> | ||
+ | <td></td> | ||
+ | </tr><tr> | ||
+ | <td>realTime(s)</td> | ||
+ | <td>0.021087233522727282</td> | ||
+ | <td>0.4520549124513621</td> | ||
+ | <td></td> | ||
+ | <td></td> | ||
+ | <td></td> | ||
+ | </tr><tr> | ||
+ | </tr><tr> | ||
+ | <td>result file</td> | ||
+ | <td><a href="">T7</a></td> | ||
+ | <td><a href="https://static.igem.org/mediawiki/2011/6/6b/Paris2011TRNA.txt">tRNA</a></td> | ||
+ | <td><a href="">Insulin</a></td> | ||
+ | <td><a href="">GFP</a></td> | ||
+ | <td><a href="">Glucose</a></td> | ||
+ | |||
+ | </tr> | ||
+ | </table> | ||
<code align="center"> | <code align="center"> | ||
====================================== | ====================================== |
Revision as of 15:55, 21 September 2011
Passive diffusion model
Introduction
Ours first calculations of the diffusion of the molecule inside a cell shows that the diffusion inside a cell is a very fast process (see the bottom of this page). It takes from 10 seconds to one minute for a component of the cell to have the average movement in the order of magnitude of the size of the cell. We wanted to explore further about the speed the passive diffusion can be achieved for a molecule to pass with passive diffusion, to see if the author is right or not in his statement.
On the other hand, we try to propose several ideas about the origin of the motion if this motion is indeed faster than the diffusion (see toward assisted diffusion). We also propose some experiments to verify these models.
The designs had been devised so that we can try to measure the speed of the diffusion through the tubes (see the Design & concepts page). We aimed to show if the diffusion is indeed a passive diffusion, and the author is wrong, or, if we have to think about more sophisticated models.
About the passive diffusion model
Earlier, we did some maths to calculate the speed of the diffusion inside a cell. But we were in spherical coordinates, that is to say a mathematics friendly conditions.Here, we have to deal with a more complex geometry, so we had to be helped by a computer.
The essence of the model remains the same than the stochastic motion approach, but in this space we introduce boundary limits with the shape of two cells linked by a tube. We observe the passage of the molecule through the tube and we calculate the time taken for a molecule (like a transcription factor) to leave the emitter cell, diffuse though the tube and meet a target in the receiver cell.
Model description
First we design a simple representation of the model with 2 boxes for cells and a tube between them. Each cell is designed like a 3D matrix(M*M*M) witch is proportional to the studied molecule. | |
We define the matrix size using our hypotesis about diffusion. So we create the cell matrix using this formula : We want to model each particle with a pixel, so each occupation site need to be 1px large. We calculate N, and with this other formula, we obtain the model size M of our cell |
For example : for the glucose, if we want to have a size of 1*1*1 the matrix representing the cell need to be 1000*1000*1000 and the nanotube which is 1/10 large of the cell will be 100*100*200.
We define the size of cells for each molecule of the hypothesis array :
We develop a java software to simulate this model. each simulation is done on 100 molecule in 10 row. All the 100 molecules of the simulation will start in one cell and move randomly until at least 10 molecule pass throw the nanotube to the other cell.
We use a synchronous model so all the 100 molecules move simultaneously. at each row (step of execution), all the molecule move randomly from there site to one of the 26 other possible position. |
this simulation take account of :
- Movement of particle : each particle have a random movement at each row.
- collision with cell membrane, nanotube and self collision. We have 2 model of collision :
- Random restart model
- Wait and see model
Random restart
In this collision model, when a particle collides with an other object it is reset to a random position without other particle on. This model is stastistically correct because of the random definition of the particle movement.
Wait and see
In this collision model, when a particle collide with an other object it stay in place, waiting the next row, so the molecule lost one movement, but it stays a credible model because of the random definition of the movement. This model have some problem due to the time lost by doing nothing with the particle. This model is actually only theoric, the java program taking more than 2 hours to calculate the movement.
How do we calculate the real time ?
A row correspond to the movement from one site to an other, so in reality it take
We obtain the real time of diffusion of 10 molecules through the nanotube with :
finally we obtain the time for diffusion from a cell to an other with :
For more accuracy, we do 10 simulation for each type of molecule, and we take the average time of those 10.
Maya modeling
this model is mapped in Maya for a user friendly aspect. This graphic reprsentaton is just for an idea of diffusion and it doesn't use the diffusion equation used by the java program
Results
results extracted from java output for 100 molecules and 10 row. the average time for one molecule is calculated from 10 molecule passing fromm cell1 to cell2
Molecules name | T7 | tRNA | Insulin | GFP | Glucose |
number of row | 37870.95 | 30980.83 | |||
realTime(s) | 0.021087233522727282 | 0.4520549124513621 | |||
result file | T7 | tRNA | Insulin | GFP | Glucose |
======================================
======================================
Characteristic size of particles :
--------------------------------------
Glucose : 1.0E-9
T7 : 1.0500000000000001E-8
tRNA : 1.5000000000000002E-8
Insulin : 5.0E-9
GFP : 6.000000000000001E-9
======================================
Diffusion coefficient of particles :
--------------------------------------
Glucose : 6.0E-10
T7 : 3.3E-11
tRNA : 2.5699999999999998E-12
Insulin : 1.5E-10
GFP : 2.7000000000000004E-11
======================================
======================================
you choose : Glucose
======================================
======================================
simulation 0
row : 0 time : 742140
row : 0 Realtime : 2.0615000000000005E-4
--------------------------------------
simulation 1
row : 1 time : 437512
row : 1 Realtime : 1.215312777777778E-4
--------------------------------------
simulation 2
row : 2 time : 726945
row : 2 Realtime : 2.0192927777777783E-4
--------------------------------------
simulation 3
row : 3 time : 1012357
row : 3 Realtime : 2.812104444444445E-4
--------------------------------------
simulation 4
row : 4 time : 329636
row : 4 Realtime : 9.156569444444446E-5
--------------------------------------
simulation 5
row : 5 time : 468676
row : 5 Realtime : 1.301878055555556E-4
--------------------------------------
simulation 6
row : 6 time : 626588
row : 6 Realtime : 1.7405227777777778E-4
--------------------------------------
simulation 7
row : 7 time : 376689
row : 7 Realtime : 1.0463602777777779E-4
--------------------------------------
simulation 8
row : 8 time : 486871
row : 8 Realtime : 1.3524216666666668E-4
--------------------------------------
simulation 9
row : 9 time : 738944
row : 9 Realtime : 2.0526238888888892E-4
--------------------------------------
======================================
======================================
average time for 10 row : 594636.25
--------------------------------------
average Realtime for 10 row : 1.6517673611111114E-4
======================================
======================================
Conclusion
Actually, because of the really short difference between passive end active diffusion, we can' say what is the nanotube type of fusion. the difference between those two model can be interpreted as noise in simulations.. Maybe a more precise method would provide better distinction between those two hypothesis.