# Team:Tokyo Tech/Project/Artificial Cooperation System/modeling

### From 2010.igem.org

(→Abstract) |
|||

(16 intermediate revisions not shown) | |||

Line 38: | Line 38: | ||

<div id="tf_SubWrapper"> | <div id="tf_SubWrapper"> | ||

<font size="5"><b>3-4 modeling</b></font> | <font size="5"><b>3-4 modeling</b></font> | ||

- | |||

==Abstract== | ==Abstract== | ||

- | In order to confirm | + | In order to confirm the feasibility of the Artificial Cooperation System, we simulated the system under typical four experimental conditions. First, when the system worked and any antibiotics did not exist, cell A and cell B in the system showed no difference of growth (Fig.3-4-1). Second, when the system did not work and chloramphenicol existed, the number of cell A decreased while the number of cell B increased (Fig.3-4-2). Third, when the system worked and chloramphenicol existed, the number of cell A increased after temporal decline (Fig.3-4-3), which means it was rescued by Artificial Cooperation System. Fourth, in the contrast, when the system worked and Kanamycin existed, the number of cell B increased after temporal decline (Fig.3-4-4). From the above results, the feasibility of the Artificial Cooperation System was confirmed. |

+ | |||

+ | [[IMAGE: tokyotech_no_anti_sim_result.jpg|300px|left|thumb|Fig3-4-1. When the system worked and any antibiotics did not exist (work done by Shohei Kitano)]] | ||

+ | [[IMAGE: tokyotech_no_system_sim_result.jpg|300px|left|thumb|Fig3-4-2. when the system did not work and chloramphenicol existed (work done by Shohei Kitano)]] | ||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | [[IMAGE: tokyotech_cm_sim_result.jpg|300px|left|thumb|Fig3-4-3 when the system worked and chloramphenicol existed (work done by Shohei Kitano)]] | ||

+ | [[IMAGE: tokyotech_kan_sim_result.jpg|300px|left|thumb|Fig3-4-4 when the system worked and Kanamycin existed (work done by Shohei Kitano) ]] | ||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

+ | |||

- | |||

- | |||

Line 51: | Line 92: | ||

- | |||

- | |||

==Model development== | ==Model development== | ||

- | In order to | + | In order to simulate the Artificial Cooperation System, we developed an ordinary differential equation model. The state variables and parameters used in the model are described in detail in [[Team:Tokyo_Tech/Project/Artificial_Cooperation_System/modeling|Table1]]. |

- | [[IMAGE:Tokyotech_S1-S8.jpg| | + | [[IMAGE:Tokyotech_S1-S8.jpg|600px]] |

- | + | In the Artificial Cooperation System, the major kinetic events are: cell population growth; resistance gene expression under the control of quorum sensing activation promoters; AHLs synthesis and degradation; I protein synthesis under the control of quorum sensing repressive promoters. These kinetic events are contained in the equations. Following sentences describe how the equations are developed. | |

'''I, Cell Population Growth''' | '''I, Cell Population Growth''' | ||

- | The growth | + | The growth rate of each type of cell without effect of antibiotics can be described by the following logistic growth equation: |

- | |||

- | + | [[IMAGE:Tokyotech S9.jpg|400px]], | |

- | + | where N denotes the number of cells, Nm denotes the maximum allowed number of cells due to nutrient limitation, and μ is a constant coefficient. The maximum growth rate is proportional to the number of cells and the growth rate is inhibited by the number of cells due to nutrient limitation. | |

- | + | When two types of cell coexist, the growth rates of the both type of cell are described in the following equations: | |

- | [[IMAGE:Tokyotech | + | [[IMAGE:Tokyotech S10S11.jpg|400px]] |

- | The | + | The maximum growth rate of each type of cell is proportional to the number of their type of cells and the growths are inhibited by the total number of the cells in the system (1-(N<sub>A</sub>+ N<sub>B</sub>)/ N<sub>m</sub>) |

- | + | The growth inhibition by antibiotics is represented by the following equations: | |

- | + | [[IMAGE:Tokyotech S12S13.jpg|400px]] | |

- | [[IMAGE:Tokyotech S1S2.jpg| | + | where Cm and Kan denote chloramphenicol and kanamycin concentration respectively. The rate of decreasing the number of cells is proportional to the product of multiplying antibiotics concentration and the number of cell. Here, &gamma is the rate constant. |

+ | |||

+ | |||

+ | The antibiotic resistance enzyme releases the antibiotic growth inhibition dependent on the enzyme concentration. Then, S12 and S13 are replaced by the following equations: | ||

+ | |||

+ | [[IMAGE:Tokyotech S14S15.jpg|400px]] | ||

+ | |||

+ | Where CmR and KanR are the concentration of the chloramphenicol and kanamycin resistance enzyme, respectively, and η is constant for specific gene. Here, if (Cm-η<sub>Cm</sub>CmR) or (Kan-η<sub>Kan</sub>KanR) is negative, these terms are set to zero during simulation. As a consequence, the dynamics of the number of cells are described by equation S1 and S2． | ||

+ | |||

+ | [[IMAGE:Tokyotech S1S2.jpg|400px]] | ||

'''II, Resistance gene expression under the control of AHLs''' | '''II, Resistance gene expression under the control of AHLs''' | ||

- | In our system, the expression level of resistance gene is | + | In our system, the expression level of resistance gene is regulated by R protein-AHL complex. R proteins, LuxR and LasR, are assumed to be synthesized sufficiently so that the concentration of the complex depends on only the AHL concentration. Thus the synthesis rates of resistance enzyme are described by Hill function dependent on the AHL concentration (equation S16 and S17). |

- | [[IMAGE:Tokyotech S16S17.jpg| | + | [[IMAGE:Tokyotech S16S17.jpg|400px]] |

- | Together with | + | Together with the leaky production of the resistance enzymes and the degradation of resistance enzymes, equation S3 and S4 are employed to simulate the expression of resistance genes. |

- | |||

- | + | ||

+ | [[IMAGE:Tokyotech S3S4.jpg|400px]] | ||

+ | |||

+ | where leakyCmR and leakyKanR denote the leaky production of chloramphenicol and kanamycin resistance enzyme respectively, k<sub>CmR</sub> and k<sub>KanR</sub> denote the maximum synthesis rate of each enzyme, respectively. n<sub>1</sub> and n<sub>2</sub> are Hill coefficient. m<sub>Cm</sub> and m<sub>KanR</sub> denote the AHL concentration required for producing half level of each enzyme, respectively. | ||

'''III, AHL synthesis and degradation''' | '''III, AHL synthesis and degradation''' | ||

- | |||

- | + | AHL is enzymatically synthesized by I protein from some substrates. For the simplicity, we assumed the amount of substrates is sufficient so that the AHL synthesis rate are estimated to be proportional to the cognate I protein. The degradation rate of AHL is assumed to decay with first-order kinetic. The dynamical behaviors of AHLs concentration are described by the following equations (S5 and S7): | |

- | k<sub>3OC6HSL</sub>, k<sub>3OC12HSL</sub>, d<sub>3OC6HSL</sub> and d<sub>3OC12HSL</sub> | + | [[IMAGE:Tokyotech S5S7.jpg|400px]] |

+ | |||

+ | where, k<sub>3OC6HSL</sub>, k<sub>3OC12HSL</sub>, d<sub>3OC6HSL</sub> and d<sub>3OC12HSL</sub> denote constant coefficients. | ||

'''IV, I protein production under the control of AHL dependent promoter''' | '''IV, I protein production under the control of AHL dependent promoter''' | ||

- | The | + | The synthesis of I proteins can be generated by Hill function according to the concentration of AHL like the expression of resistance gene. Though the resistance genes are activated by AHLs, the production of I proteins is repressed by AHLs. So, we employ the following Hill function (equation S18 and S19) |

- | |||

- | k<sub>LuxI</sub> and k<sub>LasI</sub> | + | [[IMAGE:Tokyotech S18S19.jpg|400px]] |

+ | |||

+ | Where, k<sub>LuxI</sub> and k<sub>LasI</sub> denote maximum production rate of LuxI and LasI. When the concentration of AHL is 0nM. n<sub>3</sub> and n<sub>4</sub> are Hill coefficients. m<sub>LuxI</sub> and m<sub>LasI</sub>denote the AHL concentration required for producing half level of I proteins. Therefore dynamics equations of I proteins production are listed in equation S6 and S8. | ||

<!--S6 S8--> | <!--S6 S8--> | ||

- | [[IMAGE:Tokyotech S6S8.jpg| | + | [[IMAGE:Tokyotech S6S8.jpg|400px]] |

- | d<sub>LuxI</sub>and d<sub>LasI</sub> | + | Where d<sub>LuxI</sub>and d<sub>LasI</sub> denote constant coefficient. |

==Prameters== | ==Prameters== | ||

- | + | The parameters used in equation S1 to S8 are listed in Table1. Many parameter values are taken from references. However, quantitative information is not sufficient for some parameters; therefore we use educated guesses that are biologically feasible. | |

[[IMAGE:tokyotech_parameter_table.jpg|550px|left|thumb|Table1 Variables and Parameters]] | [[IMAGE:tokyotech_parameter_table.jpg|550px|left|thumb|Table1 Variables and Parameters]] | ||

+ | |||

+ | ==References== | ||

+ | 1, Bo Hu, Jin Du, Rui-yang Zou, Ying-jin Yuan, An Environment-Sensitive Synthetic Microbial Ecosystem, PLoS ONE, volume 5 issue 5, 2010<br> | ||

+ | 2, Balagadde FK, Song H, Ozaki J, Collins CH, Barnet M, et al. (2008) A synthetic Escherichia coli predator-prey ecosystem. Mol Syst Biol 4: 187.<br> | ||

+ | 3, Gardner TS, Cantor CR, Collins JJ (2000) Construction of a genetic toggle switch in Escherichia coli. Nature 403: 339–342.<br> | ||

+ | 4, Brenner K, Karig DK, Weiss R, Arnold FH (2007) Engineered bidirectional | ||

+ | communication mediates a consensus in a microbial biofilm consortium. Proc | ||

+ | Natl Acad Sci USA 104: 17300–17304.<br> | ||

+ | |||

+ | |||

+ | |||

## Latest revision as of 03:58, 28 October 2010

**3-4 modeling**

## Contents |

## Abstract

In order to confirm the feasibility of the Artificial Cooperation System, we simulated the system under typical four experimental conditions. First, when the system worked and any antibiotics did not exist, cell A and cell B in the system showed no difference of growth (Fig.3-4-1). Second, when the system did not work and chloramphenicol existed, the number of cell A decreased while the number of cell B increased (Fig.3-4-2). Third, when the system worked and chloramphenicol existed, the number of cell A increased after temporal decline (Fig.3-4-3), which means it was rescued by Artificial Cooperation System. Fourth, in the contrast, when the system worked and Kanamycin existed, the number of cell B increased after temporal decline (Fig.3-4-4). From the above results, the feasibility of the Artificial Cooperation System was confirmed.

## Model development

In order to simulate the Artificial Cooperation System, we developed an ordinary differential equation model. The state variables and parameters used in the model are described in detail in **Table1**.

In the Artificial Cooperation System, the major kinetic events are: cell population growth; resistance gene expression under the control of quorum sensing activation promoters; AHLs synthesis and degradation; I protein synthesis under the control of quorum sensing repressive promoters. These kinetic events are contained in the equations. Following sentences describe how the equations are developed.

**I, Cell Population Growth**
The growth rate of each type of cell without effect of antibiotics can be described by the following logistic growth equation:

where N denotes the number of cells, Nm denotes the maximum allowed number of cells due to nutrient limitation, and μ is a constant coefficient. The maximum growth rate is proportional to the number of cells and the growth rate is inhibited by the number of cells due to nutrient limitation.

When two types of cell coexist, the growth rates of the both type of cell are described in the following equations:

The maximum growth rate of each type of cell is proportional to the number of their type of cells and the growths are inhibited by the total number of the cells in the system (1-(N_{A}+ N_{B})/ N_{m})

The growth inhibition by antibiotics is represented by the following equations:

where Cm and Kan denote chloramphenicol and kanamycin concentration respectively. The rate of decreasing the number of cells is proportional to the product of multiplying antibiotics concentration and the number of cell. Here, &gamma is the rate constant.

The antibiotic resistance enzyme releases the antibiotic growth inhibition dependent on the enzyme concentration. Then, S12 and S13 are replaced by the following equations:

Where CmR and KanR are the concentration of the chloramphenicol and kanamycin resistance enzyme, respectively, and η is constant for specific gene. Here, if (Cm-η_{Cm}CmR) or (Kan-η_{Kan}KanR) is negative, these terms are set to zero during simulation. As a consequence, the dynamics of the number of cells are described by equation S1 and S2．

**II, Resistance gene expression under the control of AHLs**
In our system, the expression level of resistance gene is regulated by R protein-AHL complex. R proteins, LuxR and LasR, are assumed to be synthesized sufficiently so that the concentration of the complex depends on only the AHL concentration. Thus the synthesis rates of resistance enzyme are described by Hill function dependent on the AHL concentration (equation S16 and S17).

Together with the leaky production of the resistance enzymes and the degradation of resistance enzymes, equation S3 and S4 are employed to simulate the expression of resistance genes.

where leakyCmR and leakyKanR denote the leaky production of chloramphenicol and kanamycin resistance enzyme respectively, k_{CmR} and k_{KanR} denote the maximum synthesis rate of each enzyme, respectively. n_{1} and n_{2} are Hill coefficient. m_{Cm} and m_{KanR} denote the AHL concentration required for producing half level of each enzyme, respectively.

**III, AHL synthesis and degradation**

AHL is enzymatically synthesized by I protein from some substrates. For the simplicity, we assumed the amount of substrates is sufficient so that the AHL synthesis rate are estimated to be proportional to the cognate I protein. The degradation rate of AHL is assumed to decay with first-order kinetic. The dynamical behaviors of AHLs concentration are described by the following equations (S5 and S7):

where, k_{3OC6HSL}, k_{3OC12HSL}, d_{3OC6HSL} and d_{3OC12HSL} denote constant coefficients.

**IV, I protein production under the control of AHL dependent promoter**
The synthesis of I proteins can be generated by Hill function according to the concentration of AHL like the expression of resistance gene. Though the resistance genes are activated by AHLs, the production of I proteins is repressed by AHLs. So, we employ the following Hill function (equation S18 and S19)

Where, k_{LuxI} and k_{LasI} denote maximum production rate of LuxI and LasI. When the concentration of AHL is 0nM. n_{3} and n_{4} are Hill coefficients. m_{LuxI} and m_{LasI}denote the AHL concentration required for producing half level of I proteins. Therefore dynamics equations of I proteins production are listed in equation S6 and S8.

Where d_{LuxI}and d_{LasI} denote constant coefficient.

## Prameters

The parameters used in equation S1 to S8 are listed in Table1. Many parameter values are taken from references. However, quantitative information is not sufficient for some parameters; therefore we use educated guesses that are biologically feasible.

## References

1, Bo Hu, Jin Du, Rui-yang Zou, Ying-jin Yuan, An Environment-Sensitive Synthetic Microbial Ecosystem, PLoS ONE, volume 5 issue 5, 2010

2, Balagadde FK, Song H, Ozaki J, Collins CH, Barnet M, et al. (2008) A synthetic Escherichia coli predator-prey ecosystem. Mol Syst Biol 4: 187.

3, Gardner TS, Cantor CR, Collins JJ (2000) Construction of a genetic toggle switch in Escherichia coli. Nature 403: 339–342.

4, Brenner K, Karig DK, Weiss R, Arnold FH (2007) Engineered bidirectional
communication mediates a consensus in a microbial biofilm consortium. Proc
Natl Acad Sci USA 104: 17300–17304.