Paris/Cell auto 2
From 2007.igem.org
Line 20: | Line 20: | ||
BactS is a Bact which produce DAP and can grow<br> | BactS is a Bact which produce DAP and can grow<br> | ||
BactG is a Bact which consume DAP and can divide or differentiate <br><br><br> | BactG is a Bact which consume DAP and can divide or differentiate <br><br><br> | ||
+ | |||
+ | ='''Case 2'''= | ||
+ | |||
+ | =='''''We produce this set of rules'''''== | ||
+ | |||
+ | <br> | ||
+ | Mechanic forces<br> | ||
+ | * We create a spring between the center of each Bact, then we compute the forces related to this spring and we update the position of the cells (adding noise to it)<br> | ||
+ | For bactS<br> | ||
+ | *if random < Probability of death then | ||
+ | :::BactS=null<br> | ||
+ | :else | ||
+ | ::if random < probability to grow & size < max cell size then | ||
+ | :::BactS=BactS+{new size=size+delta} | ||
+ | ::else | ||
+ | :::nothing<br> | ||
+ | *Produce DAP | ||
+ | <br> | ||
+ | For BactG<br> | ||
+ | *DAP'=DAP - self consumed DAP - diffused DAP | ||
+ | *if enough DAP then | ||
+ | ::if random< probability of differentiation then | ||
+ | :::BactG=BactS | ||
+ | ::else | ||
+ | :::BactG= BactG+{DAP'} | ||
+ | :else | ||
+ | ::if size > max size then | ||
+ | :::if probability to divide > random & DAP'> minimal needed to divide then | ||
+ | ::::BactG = 2 BactG with minimal size | ||
+ | :::else BactG= BactG +{DAP=DAP'} | ||
+ | ::else | ||
+ | :::if random < probability to grow then | ||
+ | ::::BactG = BactG + {new size= size + delta} | ||
+ | :::else | ||
+ | ::::nothing | ||
+ | <br><br><br> | ||
+ | |||
+ | =='''''Initial state'''''== | ||
+ | <br> | ||
+ | 4 BactS and a BactG in the middle | ||
+ | <br><br><br> | ||
+ | |||
+ | =='''''Parameters'''''== | ||
+ | <br> | ||
+ | We have 8 parameters and we can had noise for each of them.<br> | ||
+ | '''Mechanic'''<br> | ||
+ | *DT time step<br> | ||
+ | *K constant of the spring<br> | ||
+ | *Mu variation of position<br> | ||
+ | *R0_Gm minimal size of a BactG (after division) | ||
+ | *R0_G maximal size of a BactG (before division) | ||
+ | *R0_S maximal size of BactS | ||
+ | <br> | ||
+ | '''In Bact'''<br> | ||
+ | *Diff diffusion constant<br> | ||
+ | <br> | ||
+ | '''In BactS:'''<br> | ||
+ | *Diffp probability of differentiation<br> | ||
+ | *DEPOT production of DAP<br> | ||
+ | *DeathSP probability of death <br> | ||
+ | *CroitS probability of growth <br> | ||
+ | <br> | ||
+ | '''In BactG:'''<br> | ||
+ | *CONS Dap consumed<br> | ||
+ | *DivG probability of division<br> | ||
+ | *CroitG probability of growth<br> | ||
+ | <br><br><br> | ||
+ | |||
+ | =='''''Output'''''== | ||
+ | <br> | ||
+ | The output is two videos showing the evolution of the organism<br> | ||
+ | <center> | ||
+ | <html> | ||
+ | <embed style="width:350px; height:300px;" id="VideoPlayback" type="application/x-shockwave-flash" src="http://video.google.com/googleplayer.swf?docid=2051019221671269340&hl=en" flashvars=""; > </embed> | ||
+ | <embed style="width:350px; height:300px;" id="VideoPlayback" type="application/x-shockwave-flash" src="http://video.google.com/googleplayer.swf?docid=2422231061053866723&hl=en=en =en" flashvars=""> </embed></html></center> | ||
+ | |||
+ | *The first video show a first comportment | ||
+ | :Red : BactG | ||
+ | :Green : BactS | ||
+ | :dark<->light Bact : low<->high DAP | ||
+ | :The number of bacteria (G or S) increase with the time | ||
+ | *The second video show a second comportment | ||
+ | :Red BactG | ||
+ | :Green BactS | ||
+ | :dark<->light Bact : low<->high DAP | ||
+ | :The number of cells is constant and maintains itself | ||
+ | |||
+ | After playing with the parameters we can isolate 4 kinds of comportment.<br> | ||
+ | 2 of them are not really interested the system doesn't evolve or collapse (all the bacteria become S type).<br> | ||
+ | The other 2 comportment show that our system can lead to an evolving organism developing itself and colonizing the environment or it can stay stable like a tissue or an organ. | ||
+ | |||
+ | |||
='''Case 1'''= | ='''Case 1'''= |
Revision as of 14:03, 22 October 2007
Contents |
Spatial simulation
We try with this work to characterize the diffusion of the DAP and the effect on the cells. We have a growing culture with germinal cells and somatic cells.
We want to see if we can have different kinds of evolution for our cells
We make some hypothesis:
We work with a evolving population ( death for BactS and division BactG).
- Case 1 : The differentiation is DAP dependent, it's append when the cell as enough DAP to evolve but not enough to divide.
- Case 2 : The differentiation has a constant rate, it will be the same rate for each division cycle
The DAP is made in bacteria S, the production rate is the difference between the total production and self consummation
We consider a global variable DAP (no internal/external DAP)
The DAP is consumed in bacteria G
All the cells grow
We have 3 entities in our model
Bact it has a concentration internal of DAP and a radius. It's a cell in our automaton
BactS is a Bact which produce DAP and can grow
BactG is a Bact which consume DAP and can divide or differentiate
Case 2
We produce this set of rules
Mechanic forces
- We create a spring between the center of each Bact, then we compute the forces related to this spring and we update the position of the cells (adding noise to it)
For bactS
- if random < Probability of death then
- BactS=null
- BactS=null
- else
- if random < probability to grow & size < max cell size then
- BactS=BactS+{new size=size+delta}
- else
- nothing
- nothing
- if random < probability to grow & size < max cell size then
- Produce DAP
For BactG
- DAP'=DAP - self consumed DAP - diffused DAP
- if enough DAP then
- if random< probability of differentiation then
- BactG=BactS
- else
- BactG= BactG+{DAP'}
- if random< probability of differentiation then
- else
- if size > max size then
- if probability to divide > random & DAP'> minimal needed to divide then
- BactG = 2 BactG with minimal size
- else BactG= BactG +{DAP=DAP'}
- if probability to divide > random & DAP'> minimal needed to divide then
- else
- if random < probability to grow then
- BactG = BactG + {new size= size + delta}
- else
- nothing
- if random < probability to grow then
- if size > max size then
Initial state
4 BactS and a BactG in the middle
Parameters
We have 8 parameters and we can had noise for each of them.
Mechanic
- DT time step
- K constant of the spring
- Mu variation of position
- R0_Gm minimal size of a BactG (after division)
- R0_G maximal size of a BactG (before division)
- R0_S maximal size of BactS
In Bact
- Diff diffusion constant
In BactS:
- Diffp probability of differentiation
- DEPOT production of DAP
- DeathSP probability of death
- CroitS probability of growth
In BactG:
- CONS Dap consumed
- DivG probability of division
- CroitG probability of growth
Output
The output is two videos showing the evolution of the organism
- The first video show a first comportment
- Red : BactG
- Green : BactS
- dark<->light Bact : low<->high DAP
- The number of bacteria (G or S) increase with the time
- The second video show a second comportment
- Red BactG
- Green BactS
- dark<->light Bact : low<->high DAP
- The number of cells is constant and maintains itself
After playing with the parameters we can isolate 4 kinds of comportment.
2 of them are not really interested the system doesn't evolve or collapse (all the bacteria become S type).
The other 2 comportment show that our system can lead to an evolving organism developing itself and colonizing the environment or it can stay stable like a tissue or an organ.
Case 1
We produce this set of rules
Mechanic forces
- We create a spring between the center of each Bact, then we compute the forces related to this spring and we update the position of the cells (adding noise to it)
For bactS
- if random < Probability of death then
- BactS=null
- BactS=null
- else
- if random < probability to grow & size < max cell size then
- BactS=BactS+{new size=size+delta}
- else
- nothing
- nothing
- if random < probability to grow & size < max cell size then
- Produce DAP
For BactG
- DAP'=DAP - self consumed DAP - diffused DAP
- if enough DAP then
- if random< probability of differentiation then
- BactG=BactS
- else
- BactG= BactG+{DAP'}
- if random< probability of differentiation then
- else
- if size > max size then
- if probability to divide > random & DAP'> minimal needed to divide then
- BactG = 2 BactG with minimal size
- else BactG= BactG +{DAP=DAP'}
- if probability to divide > random & DAP'> minimal needed to divide then
- else
- if random < probability to grow then
- BactG = BactG + {new size= size + delta}
- else
- nothing
- if random < probability to grow then
- if size > max size then
Initial state
4 BactS and a BactG in the middle
Parameters
We have 8 parameters and we can had noise for each of them.
Mechanic
- DT time step
- K constant of the spring
- Mu variation of position
- R0_Gm minimal size of a BactG (after division)
- R0_G maximal size of a BactG (before division)
- R0_S maximal size of BactS
In Bact
- Diff diffusion constant
In BactS:
- Diffp probability of differentiation
- DEPOT production of DAP
- DeathSP probability of death
- CroitS probability of growth
In BactG:
- CONS Dap consumed
- DivG probability of division
- CroitG probability of growth
Output
The output is two videos showing the evolution of the organism
- The first video show a first comportment
- Red : BactG
- Green : BactS
- dark<->light Bact : low<->high DAP
- The number of bacteria (G or S) increase with the time
- The second video show a second comportment
- Red BactG
- Green BactS
- dark<->light Bact : low<->high DAP
- The number of cells is constant and maintains itself
After playing with the parameters we can isolate 4 kinds of comportment.
2 of them are not really interested the system doesn't evolve or collapse (all the bacteria become S type).
The other 2 comportment show that our system can lead to an evolving organism developing itself and colonizing the environment or it can stay stable like a tissue or an organ.