Team:Paris Liliane Bettencourt/Project/Population counter/model
From 2010.igem.org
(Difference between revisions)
Line 55: | Line 55: | ||
We can develop the information we have in a system of equations:<br /> | We can develop the information we have in a system of equations:<br /> | ||
<ul> | <ul> | ||
- | <li> In a first approach, we will consider that the dynamics of production of LuxI is represented by a step function: [LuxI] is 0 before the recombination, and 1 after. If the time between the two states cannot be neglected compared to the average time between two recombinations, we may have to reconsider this simple model. (which is probably the case since it is roughly the expected time of a cell recombinating eg close to 20 minutes, way larger than the diffusion significant times)</li> | + | <li> <h5>LuxI</h5> In a first approach, we will consider that the dynamics of production of LuxI is represented by a step function: [LuxI] is 0 before the recombination, and 1 after. If the time between the two states cannot be neglected compared to the average time between two recombinations, we may have to reconsider this simple model. (which is probably the case since it is roughly the expected time of a cell recombinating eg close to 20 minutes, way larger than the diffusion significant times)</li> |
<li> Let's call <em>R</em> the constant concentration of LuxR in a cell (how <em>"constant"</em> is it?). Some of it is binded with AHL and we called it LuxR*, so this gives:<br /> | <li> Let's call <em>R</em> the constant concentration of LuxR in a cell (how <em>"constant"</em> is it?). Some of it is binded with AHL and we called it LuxR*, so this gives:<br /> | ||
Line 74: | Line 74: | ||
</li> | </li> | ||
- | <li>Following the paper "A synchronized quorum of genetic clocks" by Hasty and al., we define the intracellular concentration of AHL, [AHL<sub>i</sub>], and the external one, [AHL<sub>e</sub>].<br /> | + | <li> <h5>AHL</h5> Following the paper "A synchronized quorum of genetic clocks" by Hasty and al., we define the intracellular concentration of AHL, [AHL<sub>i</sub>], and the external one, [AHL<sub>e</sub>].<br /> |
The concentration of AHL in the room [AHL<sub>e</sub>]decreases as the microfluidic device evacuates it, which gives birth to a diffusion term in the equation caracterized by the kinetic constant kdiff. In parallel, the internal and external AHL concentrations change as AHL is produced by the cells and can go throuh the cell membrane (kinetic constant kI and kmembrane respectively). For now, we shall assume the diffusion of AHL inside the device is instantaneous. <img src="https://static.igem.org/mediawiki/2010/2/22/Eq_AHLi.jpg" alt="Eq_AHLi" title="Equation for the intra-cellular AHL" align="center" height="50"/><br /> | The concentration of AHL in the room [AHL<sub>e</sub>]decreases as the microfluidic device evacuates it, which gives birth to a diffusion term in the equation caracterized by the kinetic constant kdiff. In parallel, the internal and external AHL concentrations change as AHL is produced by the cells and can go throuh the cell membrane (kinetic constant kI and kmembrane respectively). For now, we shall assume the diffusion of AHL inside the device is instantaneous. <img src="https://static.igem.org/mediawiki/2010/2/22/Eq_AHLi.jpg" alt="Eq_AHLi" title="Equation for the intra-cellular AHL" align="center" height="50"/><br /> | ||
Line 94: | Line 94: | ||
<li> Arabinose is modeled simply in the model with a diffusion equation: <br/> | <li> Arabinose is modeled simply in the model with a diffusion equation: <br/> | ||
<center> | <center> | ||
- | <img src="https://static.igem.org/mediawiki/2010/f/ff/Eq_arabinose.jpg" alt="arabinose" title="Diffusion equation, arabinose" align="center" height=" | + | <img src="https://static.igem.org/mediawiki/2010/f/ff/Eq_arabinose.jpg" alt="arabinose" title="Diffusion equation, arabinose" align="center" height="53"/></center><br /> |
with the boundary condition at the border between the device and the tunnel corresponding to a positive value if a pulse is being sent, and 0 otherwise. Thus we get the profile of arabinose independently of enerything else, except the pulse itself. | with the boundary condition at the border between the device and the tunnel corresponding to a positive value if a pulse is being sent, and 0 otherwise. Thus we get the profile of arabinose independently of enerything else, except the pulse itself. | ||
</li> | </li> |
Revision as of 20:51, 27 October 2010
Introduction
We aim at setting up a biological device counting random events occuring in a microfluidic device (see figure 1). These events are recombinations in the cells of the device, stimulated by pulses of arabinose coming from the tunnel. Once a certain number of pulses have been triggered, we expect our device to start emitting green fluorescence. First, we shall describe the dynamics of the population of bacteria in our device as for instance its recombinations, which further lead to the rise of the concentration of produced AHL. AHL is a quorum sensing molecule that, once a certain concentration threshold is reached, triggers in our device the production of GFP. It then alerts us of the end of counting. To count how many recombinations (so how many events) are needed to get the quorum sensing response, we shall determine the role and assess the value of a few critical parameters in our model.The process
Here we will give the bigger picture, discussing how every component of the process reacts and interacts with each other. First, we send a pulse of arabinose through the tunnel. The bacteria are then expected to react in a way that is described in the next section. Those bacteria that recombinated start producing LuxI, which in turns promotes the production of AHL. This AHL will play two roles:- it can bind with LuxR in the cell. LuxR is a protein produced so that its concentration inside the cell will be considered as constant. We call it LuxRf as long as it is free from AHL, LuxR* otherwise;
- it can cross the cell membrane to bind with LuxR of other bacteria in the device.
![Scheme of the microfluidic device Mf device](https://static.igem.org/mediawiki/2010/4/40/Mf_device.jpg)
Figure 1: Scheme of the microfluidic device
Model
We can develop the information we have in a system of equations:-
LuxI
In a first approach, we will consider that the dynamics of production of LuxI is represented by a step function: [LuxI] is 0 before the recombination, and 1 after. If the time between the two states cannot be neglected compared to the average time between two recombinations, we may have to reconsider this simple model. (which is probably the case since it is roughly the expected time of a cell recombinating eg close to 20 minutes, way larger than the diffusion significant times) - Let's call R the constant concentration of LuxR in a cell (how "constant" is it?). Some of it is binded with AHL and we called it LuxR*, so this gives:
The complexation between LuxR and AHL gives a chemical equilibrium determined by the constant Kreac:
so at the equilibrium:
which gives:
-
AHL
Following the paper "A synchronized quorum of genetic clocks" by Hasty and al., we define the intracellular concentration of AHL, [AHLi], and the external one, [AHLe].
The concentration of AHL in the room [AHLe]decreases as the microfluidic device evacuates it, which gives birth to a diffusion term in the equation caracterized by the kinetic constant kdiff. In parallel, the internal and external AHL concentrations change as AHL is produced by the cells and can go throuh the cell membrane (kinetic constant kI and kmembrane respectively). For now, we shall assume the diffusion of AHL inside the device is instantaneous.
The first term says that if the cell produces LuxI, then it produces AHL at a certain rate. The second term concerns the exchange between the cell and the medium, and is proportional to the difference of concentration.
where kdiff is the constant of diffusion, related to parameters of the fluid and of the device, kmembrane is the same as previously, but is now compensated by d, the cell density, so that the exchanges between cells and the medium is balanced. A note on diffusion further.
Here we consider only the case where the concentration of AHL is homogeneous, ie the concentration inside and outside the cell are the same. This is justified since the caracteristic time of diffusion of AHL between the outside and the inside of a cell is short compared kI. The contrary would lead to a discussion we started in paragraph \ref{modif_model}. We end up with this equation:
The synthesis of LuxI.
Above, we introduced directly the concentration of LuxI without giving more details. However, we need to model the production of LuxI to calibrate the counter. This means modelling how the bacteria evolve in the device and how they react to arabinose.- Arabinose is modeled simply in the model with a diffusion equation:
with the boundary condition at the border between the device and the tunnel corresponding to a positive value if a pulse is being sent, and 0 otherwise. Thus we get the profile of arabinose independently of enerything else, except the pulse itself. - The device may have several sizes, however one device will not ever host more than 2.105 cells, which is approximately the ration between the volume of a device and the volume of a cell. Moreover, considering both populations, switched and unswitched, we will consider the bacteria are in the exponential phase, and that they can be evacuated through the path to the tunnel. Therefore, we need a model for the behaviour of the bacteria (movement and birth/death) that allows us to deal with a few cells \footnote{at the beginning, when we only have a few recombinated cells.} up to more than 105 cells. An individually-centered model would be computationnaly inefficient, while a continuous model would perform wrong computations for small numbers of bacteria.
I have implemented something relying on a binomial distribution: among n people, we pick a random number of them to move, or give birth; this random number being ruled by a probability p proportional to the gradient of concentration of bacteria. Discretising time and space, one can draw at every step time a random binomial number\footnote{X being non-negative if the first argument in the binomial law below is non-negative, negative otherwise}:
where...
Algorithm and simulations
In this paragraph, we use several constants suggested in the paper of Hasty et al. (see references), or sometimes some approximations. We shall mention these numerical values along this paragraph.Solving the diffusion equation using an implicit scheme
mainly refers to NRCStochastic Gillespie algorithm
else
![Our simulation tool, stage 1 soft1](https://static.igem.org/mediawiki/2010/e/e5/Pic_soft1.jpg)
![Our simulation tool, stage 2 soft2](https://static.igem.org/mediawiki/2010/b/be/Pic_soft2.jpg)
For now, we use a simple spatial discretisation scheme, assuming lateral diffusion is instantaneous\footnote{this may make simulations inaccurate and thus yields a piece of the program to be improved}, so that diffusion is only considered along one axis (the one from the tunnel to the top of the device (see the first figure).
![Number of switched cells, LuxI and AHL against time AHL_chart](https://static.igem.org/mediawiki/2010/3/31/Chart_AHL.jpg)
![Number of switched cells, LuxI and AHL against time AHL_chart3](https://static.igem.org/mediawiki/2010/5/5d/Chart_AHL3.jpg)
Discuss the existence of two phases, the critical parameters being the diffusion coefficient of AHL and the time between two pulses.
Although we have been able to count six pulses thanks to the previous graph of AHL, it is not always possible to do so. For instance, if we shorten the time between two pulses, keeping the same diffusion coefficient, then it is no longer possible to distinguish pulses:
![AHL against time AHL_chart2](https://static.igem.org/mediawiki/2010/1/14/Chart_AHL2.jpg)
References
- Hasty...- Numerical recipes in C, ??, CUP
- Gillespie...
-