Team:UNITS Trieste/Project

From 2011.igem.org

(Difference between revisions)
 
Line 111: Line 111:
                     <div class="four_col">
                     <div class="four_col">
                         <p>
                         <p>
-
                            Our goals were 1) To test if our genetic design works as planned; 2) To get preliminary ideas on the main parameters of the two-compartment reactor to be built for the lab experiments. In vivo tests would be clearly too costly and time consuming for this purpose.
+
Our goals were 1) To test if our genetic design works as planned; 2) To get preliminary ideas on the main parameters of the two-compartment reactor to be built for the lab experiments. In vivo tests would be clearly too costly and time consuming for this purpose.
                         </p>
                         </p>
                     </div>
                     </div>
Line 117: Line 117:
                     <div class="four_col">
                     <div class="four_col">
                         <p>
                         <p>
-
                            Principle We designed a simple model in which QS bacteria have two states. In the i) ground state, glucose intake and signal production is at a low level. If the signal concentration reaches a threshold, the cells enter into an ii) active state characterized by higher glucose intake, higher signal production and the production of glucose from cellobiose. (Figure 1).
+
Principle we designed a simple model in which QS bacteria have two states. In the i) ground state, glucose intake and signal production is at a low level. If the signal concentration reaches a threshold, the cells enter into an ii) active state characterized by higher glucose intake, higher signal production and the production of glucose from cellobiose. (Figure 1).
                         </p>
                         </p>
                         <p>
                         <p>
-
                            The cells are in a perfectly mixed, closed environment. At the beginning of the simulation, a given number of ground state cells are placed into the medium, and the simulation proceeds in discrete time steps. At every step, the cells take up nutrients and carry out the signal and/or glucose production depending on their state of activation. When the cells accumulate a certain amount of energy ("glucose equivalents"), they divide. This program is repeated at every time step. As a result, the cells produce a growth curve quite similar to that seen in liquid cultures.
+
The cells are in a perfectly mixed, closed environment. At the beginning of the simulation, a given number of ground state cells are placed into the medium, and the simulation proceeds in discrete time steps. At every step, the cells take up nutrients and carry out the signal and/or glucose production depending on their state of activation. When the cells accumulate a certain amount of energy ("glucose equivalents"), they divide. This program is repeated at every time step. As a result, the cells produce a growth curve quite similar to that seen in liquid cultures.
                         </p>
                         </p>
                     </div>
                     </div>
Line 134: Line 134:
                     <div class="four_col">
                     <div class="four_col">
                         <p>
                         <p>
-
                            This is a very simple model. In which the cells do not move and the signal is not diffusing within the compartment, since the medium is perfectly mixed. At each time step we can record the number of cells, the concentration of the solutes, etc. This is a so-called agent-based model since each cell-agent executes its own program that depends on its state of activation and stored energy.
+
This is a very simple model, in which the cells do not move and the signal is not diffusing within the compartment, since the medium is perfectly mixed. At each time step we can record the number of cells, the concentration of the solutes, etc. This is a so-called agent-based model since each cell-agent executes its own program that depends on its state of activation and stored energy.
                         </p>
                         </p>
                         <p>
                         <p>
-
                            Two-compartment reactor. In this setup, two different kinds of cells are put into two equal, perfectly mixed compartments that are separated by a large pore-size semi-permeable membrane. The cells themselves can not pas through the membrane, but the dissolved materials can freely move between the compartments. As compared to the previous one-comaprtment model, the only difference is that the solutes need to equilibrate between the two compartments at every time step. If we assume perfect equilibration between the compartments, the concentration of a solute would be the average of those mesured in the two compartments at each time step. Here we employ a simple trick, we introduce a virtual diffusion coefficient that regulates the exchange between the two compartments. The diffusion coefficient D was defined in such a manner that its values be between zero and 1.0. If it is zero, the two compartments do not communicate. It if is 1.0, the compartments are fully equilibrated, i.e. the equilibrated concentration wiill be the average of the concentrations within the two compartments at every time-step.
+
Two-compartment reactor. In this setup, two different kinds of cells are put into two equal, perfectly mixed compartments that are separated by a large pore-size semi-permeable membrane. The cells themselves can’t pass through the membrane, but the dissolved materials can freely move between the compartments. As compared to the previous one-compartment model, the only difference is that the solutes need to equilibrate between the two compartments at every time step. If we assume perfect equilibration between the compartments, the concentration of a solute would be the average of those measured in the two compartments at each time step. Here we employ a simple trick, we introduce a virtual diffusion coefficient that regulates the exchange between the two compartments.<br/>
 +
The diffusion coefficient D was defined in such a manner that its values be between zero and 1.0. If it is zero, the two compartments do not communicate. It if is 1.0, the compartments are fully equilibrated, i.e. the equilibrated concentration will be the average of the concentrations within the two compartments at every time-step.                          
                         </p>
                         </p>
                         <p>
                         <p>
-
                            The model was implemented using an existing Matlab code written by S. Netotea and A. Kerenyi for modelling the swarming of quorum sensing bacteria on agar plates (Netotea et al, 2009, Venturi et al, 2010) [1][2], kindly provided to us by the authors. This code had to be slightly simplified, as described in Appendix 1.  
+
The model was implemented using an existing Matlab code written by S. Netotea and A. Kerenyi for modeling the swarming of quorum sensing bacteria on agar plates (Netotea et al, 2009, Venturi et al, 2010) [1][2], kindly provided to us by the authors. This code had to be slightly simplified, as described in Appendix 1.                            
                         </p>
                         </p>
                     </div>
                     </div>
Line 146: Line 147:
                     <div class="four_col">
                     <div class="four_col">
                         <p>
                         <p>
-
                            We carried out simulation runs at a large number of parameter settings, changing the number of starting population, the ratio of the populations, the starting concentrations of the various solutes, etc. We valuated the results in a qualitatie way, i.e. groth vs. no growth, slower or faster growth. As a comparison, we used a wt model that responds to its own signal.
+
We carried out simulation runs at a large number of parameter settings, changing the number of starting population, the ratio of the populations, the starting concentrations of the various solutes, etc. We evaluated the results in a qualitative way, i.e. growth vs. no growth, slower or faster growth. As a comparison, we used a wt model that responds to its own signal.
                         </p>
                         </p>
                     </div>
                     </div>
Line 160: Line 161:
                         <p class="descr-img">
                         <p class="descr-img">
                             <b>Figure 2b.</b>
                             <b>Figure 2b.</b>
-
                            Carbon source transformation - from cellobiose to glucose - and consumption.
+
Carbon source transformation - from cellobiose to glucose - and consumption.
                         </p>
                         </p>
<img src="https://static.igem.org/mediawiki/2011/d/d4/Units-model_2b2.png" width="600"/>
<img src="https://static.igem.org/mediawiki/2011/d/d4/Units-model_2b2.png" width="600"/>
Line 169: Line 170:
                     </div>
                     </div>
                     <div class="four_col">
                     <div class="four_col">
-
                    Then we extended the experiments to the species designd to depend on each other.The main finding were as follows:
+
Then we extended the experiments to the species designed to depend on each other. The main finding were as follows:
                     </br>
                     </br>
-
                     <span class="list">- Mutually dependent bacterial cells that respond only to the signal of the other species, can grow within the same compartment, same as the wild type. The ratio of the species is 50:50 %. This is a stable equilibrium that can be reached even if one of the populations starts with a single cell.</span>
+
                     <span class="list">-&nbsp;&nbsp; Mutually dependent bacterial cells that respond only to the signal of the other species, can grow within the same compartment, same as the wild type. The ratio of the species is 50:50 %. This is a stable equilibrium that can be reached even if one of the populations starts with a single cell.</span>
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
Line 177: Line 178:
                     <p class="descr-img">
                     <p class="descr-img">
                         <b>Figure 3.</b>
                         <b>Figure 3.</b>
-
                        When the two bacterial populations are in the same compartment (c = 1) their growth is synchronized.
+
                      When the two bacterial populations are in the same compartment (c = 1) their growth is synchronized.
                     </p>
                     </p>
                     <div class="clear">
                     <div class="clear">
Line 183: Line 184:
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
-
                     <span class="list">- When put into separate but communicating compartments, mutually dependent cell populations will grow, but the speed of growth will depend on the intensity of material exchange between the compartments. If there is no exchange, mutually dependent cells will not grow.</span>
+
                     <span class="list">-&nbsp;&nbsp; When put into separate but communicating compartments, mutually dependent cell populations will grow, but the speed of growth will depend on the intensity of material exchange between the compartments. If there is no exchange, mutually dependent cells will not grow.</span>
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
Line 196: Line 197:
                     <p class="descr-img">
                     <p class="descr-img">
                         <b>Figure 4b.</b>
                         <b>Figure 4b.</b>
-
                        In this last condition (Figure 3) we can observe the usual pattern in the carbon source transformation and consumption. The concentration of the signals molecules (OC8 and OC 12 AHL) follow the cells growth patter for both spicies.
+
In this last condition (Figure 3) we can observe the usual pattern in the carbon source transformation and consumption. The concentration of the signals molecules (OC8 and OC 12 AHL) follow the cells growth pattern for both species.
                     </p>
                     </p>
<img src="https://static.igem.org/mediawiki/2011/e/e7/Units-model_4b2.png" width="600"/>
<img src="https://static.igem.org/mediawiki/2011/e/e7/Units-model_4b2.png" width="600"/>
Line 203: Line 204:
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
-
                     <span class="list">- If we put a ternary community into the two compartments, in such a way that lactamase producing partner is in one compartment, and the mutually dependent bacteria are in the other compartment, the system will start only if we jump start it by adding glucose and signal. </span>
+
                     <span class="list">-&nbsp;&nbsp; If we put a ternary community into the two compartments, in such a way that lactamase producing partner is in one compartment, and the mutually dependent bacteria are in the other compartment, the system will start only if we jump start it by adding glucose and signal.</span>
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
Line 209: Line 210:
                     <p class="descr-img">
                     <p class="descr-img">
                         <b>Figure 5a.</b>
                         <b>Figure 5a.</b>
-
                         Cell growth of the ternary system without and with the "jump start" .
+
                         Cell growth of the ternary system without and with the "jump start".
                     </p>
                     </p>
<img src="https://static.igem.org/mediawiki/2011/5/5e/Units-model_5a2.png" width="600"/>
<img src="https://static.igem.org/mediawiki/2011/5/5e/Units-model_5a2.png" width="600"/>
Line 217: Line 218:
                     <p class="descr-img">
                     <p class="descr-img">
                         <b>Figure 5b.</b>
                         <b>Figure 5b.</b>
-
                        The bacterial signaling also switch-on the production of beta-lactamase by HeLa cells that allows the bacterial growth .
+
                      The bacterial signaling also switch-on the production of beta-lactamase by HeLa cells that allows the bacterial growth.
                     </p>
                     </p>
                     <div class="clear">
                     <div class="clear">
Line 224: Line 225:
                 <span class="centered"><b>Conclusion of the simulations</b></span>
                 <span class="centered"><b>Conclusion of the simulations</b></span>
                 <div class="four_col">
                 <div class="four_col">
-
                     We used a highly simplified qualitative model in order to have preliminary insights into the logics of our system, in hope of obtaining indications regarding the physical setup of the reaction chamber to be used in the lab experiments. The model shows that the system is in princple "viable", howeer the growth can be limited by the material exchange between the compartments. Also, the system may need to be jump-started by adding a certain amount of glucose and signal. Naturally, these predictions are qualitative and have to be checked by experiment.
+
                     We used a highly simplified qualitative model in order to have preliminary insights into the logics of our system, in hope of obtaining indications regarding the physical setup of the reaction chamber to be used in the lab experiments. The model shows that the system is in princple "viable", howeer the growth can be limited by the material exchange between the compartments. Also, the system may need to be jump-started by adding a certain amount of glucose and signal. Naturally, these predictions are qualitative and have to be checked by experiment.
                 </div>
                 </div>
                 <div class="four_col">
                 <div class="four_col">
Line 240: Line 241:
                     <span class=""><b>Appendix 1.</b></span>
                     <span class=""><b>Appendix 1.</b></span>
                     <p>
                     <p>
-
                        The QS modelling program of S. Netotea and A. Kerenyi (Netotea et al, 2009, Venturi et al, 2010)is designed to model the swarming of QS bacteria. In this modell i) bacteria are modeled as individuals freely moving on a 2D plane, i.e. in an open environment; ii) bacteria have an individual program that activates them depending on the threshold concentration of a) signals, b) public goods, according to the known rules of QS; and iii) nutrients and solutes freely diffuse on a 2D plane which is discretized into square zones (Netotea et al, 2009). This modell was instrumental in showing that QS regulation is sufficient for a modell popultion to show density-dependent activtion, tracking of exdternal signals, co-swarming of species and community collapse.
+
The QS modeling program of S. Netotea and A. Kerenyi (Netotea et al, 2009, Venturi et al, 2010) is designed to model the swarming of QS bacteria. In this modell i) bacteria are modeled as individuals freely moving on a 2D plane, i.e. in an open environment; ii) bacteria have an individual program that activates them depending on the threshold concentration of a) signals, b) public goods, according to the known rules of QS; and iii) nutrients and solutes freely diffuse on a 2D plane which is discretized into square zones (Netotea et al, 2009). This model was instrumental in showing that QS regulation is sufficient for a model population to show density-dependent activation, tracking of external signals, co-swarming of species and community collapse.
                     </p>
                     </p>
                     <p>
                     <p>
-
                        Our system is simplified, since due to the perfect mixing within a reactor compartment, a) cells do not move by themselves; b) solutes do not diffuse (concentration is uniform throughout the compartment). So, in order to model growth within a closed compartment, the diffusion and the movement part of the original program had to be simply switched off.
+
Our system is simplified, since due to the perfect mixing within a reactor compartment, a) cells do not move by themselves; b) solutes do not diffuse (concentration is uniform throughout the compartment). So, in order to model growth within a closed compartment, the diffusion and the movement part of the original program had to be simply switched off.
                     </p>
                     </p>
                     <p>
                     <p>
-
                        In order to simulate growth in a two-compartment system, one can model the growth of bacteria in such a way, that solutes accumulated separately, i.e. the same program is executed for the two separate components. If we now want to modell the passage of solutes between the compartments, we can use a diffusion-like concentration equilibration at each time step. This setup corresponds to two compartments separated by a semipermeable membrane. The concentration of a solute can be calculated as follows:
+
 
 +
In order to simulate growth in a two-compartment system, one can model the growth of bacteria in such a way, that solutes accumulated separately, i.e. the same program is executed for the two separate components. If we now want to model the passage of solutes between the compartments, we can use a diffusion-like concentration equilibration at each time step. This setup corresponds to two compartments separated by a semi-permeable membrane. The concentration of a solute can be calculated as follows:
                     </p>
                     </p>
                     <span class="centered">c1(t), equilibrated = c1(t) – 0.5 * D* [(c1(t) - c2(t)]</span>
                     <span class="centered">c1(t), equilibrated = c1(t) – 0.5 * D* [(c1(t) - c2(t)]</span>
                     <br/>
                     <br/>
                     <p>
                     <p>
-
                        where c1 and c2 are the concentrations in compartment 1 and 2, respectively, and D is the virtual diffusion coefficient, with a value between zero and 1.0. It is easy to see, that at D=0, the concentration remains the same, at m=1, the concentration will be the average of c1 and c2.
+
where c1 and c2 are the concentrations in compartment 1 and 2, respectively, and D is the virtual diffusion coefficient, with a value between zero and 1.0. It is easy to see, that at D=0, the concentration remains the same, at m=1, the concentration will be the average of c1 and c2.
                     </p>
                     </p>
                     <p>
                     <p>
-
                        Based on the above design, the Matlab code was modified by A. Kerenyi and put to our disposal for modelling as a *.exe file that runs under Windows. The results are cell counts, signal and food concentrations as a function of the time steps These results were visualized by Excell.
+
Based on the above design, the Matlab code was modified by A. Kerenyi and put to our disposal for modeling as a *.exe file that runs under Windows. The results are cell counts, signal and food concentrations as a function of the time steps These results were visualized by Excel.
                     </p>
                     </p>
                 </div>
                 </div>

Latest revision as of 14:47, 28 October 2011

SYNBIOME OVERVIEW

An important challenge in the near future will be the optimization of bioreactors for the production of complex molecules, and the employment of bacterial systems towards this goal is gaining wide prominence. We believe that using different cell communities in these systems would improve their ability to synthesize biological molecules through cooperation. The aim of our project is thus to combine cell types belonging to different kingdoms using synthetic biology. The innovation and challenge will be to obtain a stable communication between different kingdoms and establish mutualism among them. This interdependence will be obtained through metabolic and signaling pathways in which the survival and/or growth of each cell type depends on the others.
Our project, SynBiome, is based on a three-element system: two different bacterial strains and one eukaryotic cell type that communicate through quorum sensing (QS) signalling molecules.
In order to achieve the goal of constructing this synthetic community, both the bacterial cells and the eukaryotic cell will be engineered with a genetic circuit under the regulation of N-acyl homoserine lactone (AHL) QS signals. More specifically, the two bacterial strains will respectively produce 3-oxo-dodecanoyl-HL (OC12) and 3-oxo-octanoyl-HL (OC8), which in turn will induce the production of the bacterial enzyme cellobiosidase and a eukaryotic secreted beta-lactamase (sBLA), whose expression will be guaranteed by a eukaryotic trans-activator sensible to bacterial OC8.

Our set up will ensure interdependence among the three cell types, having the bacteria surviving in an ampicillin-containing culture medium and all the cells benefitting from the free glucose made available by the cellobiosidase.
Importantly, this genetic circuit will be designed in such a way that the system can be adapted to different bacterial species and eukaryotic cell types.

Previous Next
data_frame0
SYNBIOME... an innovative synthetic biome in which bacteria and eukaryotic cells need each other to survive!

MODELING

The goal

Our goals were 1) To test if our genetic design works as planned; 2) To get preliminary ideas on the main parameters of the two-compartment reactor to be built for the lab experiments. In vivo tests would be clearly too costly and time consuming for this purpose.

The model

Principle we designed a simple model in which QS bacteria have two states. In the i) ground state, glucose intake and signal production is at a low level. If the signal concentration reaches a threshold, the cells enter into an ii) active state characterized by higher glucose intake, higher signal production and the production of glucose from cellobiose. (Figure 1).

The cells are in a perfectly mixed, closed environment. At the beginning of the simulation, a given number of ground state cells are placed into the medium, and the simulation proceeds in discrete time steps. At every step, the cells take up nutrients and carry out the signal and/or glucose production depending on their state of activation. When the cells accumulate a certain amount of energy ("glucose equivalents"), they divide. This program is repeated at every time step. As a result, the cells produce a growth curve quite similar to that seen in liquid cultures.

Figure 1. The two states of the bacterial model.

This is a very simple model, in which the cells do not move and the signal is not diffusing within the compartment, since the medium is perfectly mixed. At each time step we can record the number of cells, the concentration of the solutes, etc. This is a so-called agent-based model since each cell-agent executes its own program that depends on its state of activation and stored energy.

Two-compartment reactor. In this setup, two different kinds of cells are put into two equal, perfectly mixed compartments that are separated by a large pore-size semi-permeable membrane. The cells themselves can’t pass through the membrane, but the dissolved materials can freely move between the compartments. As compared to the previous one-compartment model, the only difference is that the solutes need to equilibrate between the two compartments at every time step. If we assume perfect equilibration between the compartments, the concentration of a solute would be the average of those measured in the two compartments at each time step. Here we employ a simple trick, we introduce a virtual diffusion coefficient that regulates the exchange between the two compartments.
The diffusion coefficient D was defined in such a manner that its values be between zero and 1.0. If it is zero, the two compartments do not communicate. It if is 1.0, the compartments are fully equilibrated, i.e. the equilibrated concentration will be the average of the concentrations within the two compartments at every time-step.

The model was implemented using an existing Matlab code written by S. Netotea and A. Kerenyi for modeling the swarming of quorum sensing bacteria on agar plates (Netotea et al, 2009, Venturi et al, 2010) [1][2], kindly provided to us by the authors. This code had to be slightly simplified, as described in Appendix 1.

Simulation results

We carried out simulation runs at a large number of parameter settings, changing the number of starting population, the ratio of the populations, the starting concentrations of the various solutes, etc. We evaluated the results in a qualitative way, i.e. growth vs. no growth, slower or faster growth. As a comparison, we used a wt model that responds to its own signal.

Figure 2a. Growth of wild type cells that produce a signal and respond to it by producing glucose from cellobiose.

Figure 2b. Carbon source transformation - from cellobiose to glucose - and consumption.

Then we extended the experiments to the species designed to depend on each other. The main finding were as follows:
-   Mutually dependent bacterial cells that respond only to the signal of the other species, can grow within the same compartment, same as the wild type. The ratio of the species is 50:50 %. This is a stable equilibrium that can be reached even if one of the populations starts with a single cell.

Figure 3. When the two bacterial populations are in the same compartment (c = 1) their growth is synchronized.

-   When put into separate but communicating compartments, mutually dependent cell populations will grow, but the speed of growth will depend on the intensity of material exchange between the compartments. If there is no exchange, mutually dependent cells will not grow.

Figure 4a. Growth of the two bacterial populations in separated compartments. c = mixing coefficient.

Figure 4b. In this last condition (Figure 3) we can observe the usual pattern in the carbon source transformation and consumption. The concentration of the signals molecules (OC8 and OC 12 AHL) follow the cells growth pattern for both species.

-   If we put a ternary community into the two compartments, in such a way that lactamase producing partner is in one compartment, and the mutually dependent bacteria are in the other compartment, the system will start only if we jump start it by adding glucose and signal.

Figure 5a. Cell growth of the ternary system without and with the "jump start".

Figure 5b. The bacterial signaling also switch-on the production of beta-lactamase by HeLa cells that allows the bacterial growth.

Conclusion of the simulations
We used a highly simplified qualitative model in order to have preliminary insights into the logics of our system, in hope of obtaining indications regarding the physical setup of the reaction chamber to be used in the lab experiments. The model shows that the system is in princple "viable", howeer the growth can be limited by the material exchange between the compartments. Also, the system may need to be jump-started by adding a certain amount of glucose and signal. Naturally, these predictions are qualitative and have to be checked by experiment.

References

[1] Netotea S, Bertani I, Steindler L, Kerényi A, Venturi V, Pongor S. (2009) .A simple model for the early events of quorum sensing in Pseudomonas aeruginosa: modeling bacterial swarming as the movement of an "activation zone". Biology Direct, 4:6 . PDF

[2] Venturi V., Bertani I., Kerényi Á., Netotea S. and Pongor S. (2010) Co-Swarming and Local Collapse: Quorum Sensing Conveys Resilience to Bacterial Communities by Localizing Cheater Mutants in Pseudomonas aeruginosa, PlosOne, 5,4, PDF


Appendix 1.

The QS modeling program of S. Netotea and A. Kerenyi (Netotea et al, 2009, Venturi et al, 2010) is designed to model the swarming of QS bacteria. In this modell i) bacteria are modeled as individuals freely moving on a 2D plane, i.e. in an open environment; ii) bacteria have an individual program that activates them depending on the threshold concentration of a) signals, b) public goods, according to the known rules of QS; and iii) nutrients and solutes freely diffuse on a 2D plane which is discretized into square zones (Netotea et al, 2009). This model was instrumental in showing that QS regulation is sufficient for a model population to show density-dependent activation, tracking of external signals, co-swarming of species and community collapse.

Our system is simplified, since due to the perfect mixing within a reactor compartment, a) cells do not move by themselves; b) solutes do not diffuse (concentration is uniform throughout the compartment). So, in order to model growth within a closed compartment, the diffusion and the movement part of the original program had to be simply switched off.

In order to simulate growth in a two-compartment system, one can model the growth of bacteria in such a way, that solutes accumulated separately, i.e. the same program is executed for the two separate components. If we now want to model the passage of solutes between the compartments, we can use a diffusion-like concentration equilibration at each time step. This setup corresponds to two compartments separated by a semi-permeable membrane. The concentration of a solute can be calculated as follows:

c1(t), equilibrated = c1(t) – 0.5 * D* [(c1(t) - c2(t)]

where c1 and c2 are the concentrations in compartment 1 and 2, respectively, and D is the virtual diffusion coefficient, with a value between zero and 1.0. It is easy to see, that at D=0, the concentration remains the same, at m=1, the concentration will be the average of c1 and c2.

Based on the above design, the Matlab code was modified by A. Kerenyi and put to our disposal for modeling as a *.exe file that runs under Windows. The results are cell counts, signal and food concentrations as a function of the time steps These results were visualized by Excel.