Team:BCCS-Bristol/BSim/Case studies/Repressilators
From 2009.igem.org
Line 1: | Line 1: | ||
{{:Team:BCCS-Bristol/Header}} | {{:Team:BCCS-Bristol/Header}} | ||
- | + | <html><h1>Repressilators coupled by quorum sensing - an agent based approach</h1></html> | |
- | + | {| class="panel" align="center" width="90%" | |
- | + | |width="33%" padding="3px"| <html><center> | |
+ | <img src="http://dl.getdropbox.com/u/1944619/Repressilators_Title_A1.png" border="0"> | ||
</center></html> | </center></html> | ||
- | + | |width="33%" padding="3px"| <html><center> | |
- | + | <img src="http://dl.getdropbox.com/u/1944619/Repressilators_Title_B1.png" border="0"> | |
+ | </center></html> | ||
+ | |width="33%" padding="3px"| <html><center><img src="http://dl.getdropbox.com/u/1944619/Repressilators_Title_C1.png"></center></html> | ||
+ | |} | ||
==Background== | ==Background== | ||
The emergence of new GRNs such as the repressilator [[#ref_elowitz|[1]]], which may be affected by factors on the population rather than the individual level, has sparked a new interest in modelling GRNs across a bacterial population. Synthetic clocks such as the repressilator may help to provide us with a deeper understanding of oscillatory behaviour in natural systems. Theoretical modelling of such systems across a population is an important step towards better understanding of natural oscillators such as the circadian clock. | The emergence of new GRNs such as the repressilator [[#ref_elowitz|[1]]], which may be affected by factors on the population rather than the individual level, has sparked a new interest in modelling GRNs across a bacterial population. Synthetic clocks such as the repressilator may help to provide us with a deeper understanding of oscillatory behaviour in natural systems. Theoretical modelling of such systems across a population is an important step towards better understanding of natural oscillators such as the circadian clock. | ||
+ | |||
+ | ===What is a repressilator?=== | ||
+ | |||
+ | ===A mathematical description of the repressilator=== | ||
+ | |||
+ | ===Previous work=== | ||
+ | |||
+ | ==Why we used BSim== | ||
Recent mathematical modelling approaches in systems biology tend to model a gene regulatory network in a single cell, and agent based models are considered in a separate context. However, some GRNs such as the repressilator can be coupled across a population of bacteria. In the case of the repressilator the GRNs are coupled by an autoinducer chemical which is free to diffuse in and out of the cell. Previous approaches to modelling this problem (see for example [[#ref_garcia|[2]]]) have all assumed that the chemical is well-mixed across the population. In reality external chemical concentrations will vary across a large space, therefore it is important to consider the effects of a nonuniform chemical field on network dynamics. | Recent mathematical modelling approaches in systems biology tend to model a gene regulatory network in a single cell, and agent based models are considered in a separate context. However, some GRNs such as the repressilator can be coupled across a population of bacteria. In the case of the repressilator the GRNs are coupled by an autoinducer chemical which is free to diffuse in and out of the cell. Previous approaches to modelling this problem (see for example [[#ref_garcia|[2]]]) have all assumed that the chemical is well-mixed across the population. In reality external chemical concentrations will vary across a large space, therefore it is important to consider the effects of a nonuniform chemical field on network dynamics. | ||
Line 17: | Line 29: | ||
In [[#ref_garcia|[2]]] it was shown that communication between cells via quorum sensing can result in population level synchronisation given a large enough cellular density. By bringing together an agent based approach with the standard ordinary differential equation methods used for modelling GRNs we hope to be able to extensively study spatial factors affecting the behaviour of the repressilator in a population. | In [[#ref_garcia|[2]]] it was shown that communication between cells via quorum sensing can result in population level synchronisation given a large enough cellular density. By bringing together an agent based approach with the standard ordinary differential equation methods used for modelling GRNs we hope to be able to extensively study spatial factors affecting the behaviour of the repressilator in a population. | ||
+ | ==BSim coupled repressilators== | ||
+ | |||
+ | ===Overview=== | ||
+ | |||
+ | ===Results=== | ||
+ | |||
+ | <html><center> | ||
+ | <object width="425" height="344"><param name="movie" value="http://www.youtube.com/v/V1aaygajIcM&hl=en&fs=1"></param><param name="allowFullScreen" value="true"></param><param name="allowscriptaccess" value="always"></param><embed src="http://www.youtube.com/v/V1aaygajIcM&hl=en&fs=1" type="application/x-shockwave-flash" allowscriptaccess="always" allowfullscreen="true" width="425" height="344"></embed></object> | ||
+ | </center></html> | ||
+ | |||
+ | The video shows 200 bacteria swimming in a 100x100x100 micron volume. Each bacterium has a system of ODE's inside which model the essential dynamics of a repressilator. In this case, the individual repressilators are coupled via a diffusing autoinducer signal (in this case AHL). The colour of a bacterium represents the level of lacI mRNA in that bacterium (the lacI gene is one of the genes present in the repressilator); the bacterium will change from yellow to red as the internal level of lacI mRNA increases. In the example shown here, all 200 of the individual repressilators are initialised with random conditions, but quickly synchronise due to the effect of the AHL communication. | ||
+ | |||
+ | ====Frequency==== | ||
+ | |||
+ | ====Phase locking==== | ||
+ | |||
+ | ====Synchronisation transition==== | ||
+ | |||
+ | ==Further work== | ||
- | |||
Revision as of 19:52, 21 October 2009
iGEM 2009
Repressilators coupled by quorum sensing - an agent based approach
| | |
Contents |
Background
The emergence of new GRNs such as the repressilator [1], which may be affected by factors on the population rather than the individual level, has sparked a new interest in modelling GRNs across a bacterial population. Synthetic clocks such as the repressilator may help to provide us with a deeper understanding of oscillatory behaviour in natural systems. Theoretical modelling of such systems across a population is an important step towards better understanding of natural oscillators such as the circadian clock.
What is a repressilator?
A mathematical description of the repressilator
Previous work
Why we used BSim
Recent mathematical modelling approaches in systems biology tend to model a gene regulatory network in a single cell, and agent based models are considered in a separate context. However, some GRNs such as the repressilator can be coupled across a population of bacteria. In the case of the repressilator the GRNs are coupled by an autoinducer chemical which is free to diffuse in and out of the cell. Previous approaches to modelling this problem (see for example [2]) have all assumed that the chemical is well-mixed across the population. In reality external chemical concentrations will vary across a large space, therefore it is important to consider the effects of a nonuniform chemical field on network dynamics.
In [2] it was shown that communication between cells via quorum sensing can result in population level synchronisation given a large enough cellular density. By bringing together an agent based approach with the standard ordinary differential equation methods used for modelling GRNs we hope to be able to extensively study spatial factors affecting the behaviour of the repressilator in a population.
BSim coupled repressilators
Overview
Results
The video shows 200 bacteria swimming in a 100x100x100 micron volume. Each bacterium has a system of ODE's inside which model the essential dynamics of a repressilator. In this case, the individual repressilators are coupled via a diffusing autoinducer signal (in this case AHL). The colour of a bacterium represents the level of lacI mRNA in that bacterium (the lacI gene is one of the genes present in the repressilator); the bacterium will change from yellow to red as the internal level of lacI mRNA increases. In the example shown here, all 200 of the individual repressilators are initialised with random conditions, but quickly synchronise due to the effect of the AHL communication.
Frequency
Phase locking
Synchronisation transition
Further work
References
- [1] Michael B. Elowitz & Stanislas Leibler - A synthetic oscillatory network of transcriptional regulators | [http://dx.doi.org/doi:10.1073/pnas.0307095101 doi:10.1073/pnas.0307095101]
- [2] J. Garcia-Ojalvo, Michael B. Elowitz, Steven H. Strogatz - Modeling a synthetic multicellular clock: Repressilators coupled by quorum sensing | [http://dx.doi.org/doi:10.1038/35002125 doi:10.1038/35002125]