Based on experiments, we have designed sulforaphane production system, curcumin production system, and temperature-induced suicide system, as well as feasibility verification. However, quantified analysis is also needed, which can be realized by modelling. In order to analyze the dynamics properties of myrosinase and DCS-CURS, the inhibition effect to cancer cells of sulforaphane and curcumin, and the efficiency of temperature-induced suicide switch, we use Michaelis-Menten equation, Arrhenius equation, Hill equation, and Logistic equation to describe the system. By fitting, we derive dynamic properties of myrosinase and DCS-CURS in 37℃ and pH7.4. Numerical simulation via MATLAB shows that the suicide switch and inhibition system work well, whereas the inhibition rate is slow and not thorough.
Dynamic models
Introduction
This section is organized as follows.
-
Fitting of dynamic properties of myrosinase and DCS-CURS using Michaelis-Menten equations.
-
Description of TcI42 suicide switch using Hill equation and Arrhenius equation. This part also compute the generating velocity of sRRz.
-
Description of population growth vua Logistics equation. The lysis of engineering bacteria by sRRz and the inhibition of cancer cell growth by sulforaphane and curcumin are also considered.
-
Summarization of overall system ODEs.
Activity of enzymes
The Michaelis-Menten equation, which is widely used to describe the velocity of enzymatic reactions, is
where is the maximum velocity for given enzyme concentration, is the concentration of substrate, is Michaelis-Menten constant, which is equal to the substrate concentration when the velocity reaches . According to experimental results, the dynamic properties for myrosinase and DCS-CURS in 37℃ and pH7.4 is shown below.
Table 1: Dynamic properties of myrosinase and DCS-CURS in 37℃ and pH7.4
Emzyme | ||
---|---|---|
Myrosinase | 15.51 | 84.15 |
DCS-CURS | 0.2862 | 71.2 |
TcI42 and the production of sRRz
Temperature-inducible promoters are a type of sequences, which can adject transcriptional activity in a particular temperature range. They are widely used in bioengineering and synthetic biology to regulate temperature-sensitive gene expression. Temperature-inducible promoters work by changing the activity of promoters through changes in temperature, thereby regulating the transcription level of genes. At low temperatures, they are inhibited and the transcriptional activity of genes is low. At high temperatures, they are activated and the transcriptional activity of the gene increases. This temperature sensitivity is achieved by the interaction of specific regulatory elements and temperature-sensing factors in the promoter sequence.
For the activity of TcI42 versus temperature, as the principle of TcI42 can be simplified into the association of ligand[1], we can consider the dissociation reaction: , where , . and represent the concentrations of ligand, TcI42, and associated TcI42, respectively, and represents Hill coefficient. According to Arrhenius equation, the velocity of positive and reversed reaction can be described as and , where and represent the activation energy of positive and reversed reactions, respectively. When the reaction reaches equilibrium, according to , the apparent dissociation constant is
Let , the enthalpy change , then the apparent dissociation constant is . If the ratio of associated TcI42 is , then according to Hill equation,
As and do not change when temperature changes, we have
where and are two constants. Denote as the temperature when half of TcI42 is associated. As , and for temperature-induced switch, is usually a real number which is much smaller than 1, then
In suicide switch, the transcription and translation process can be described by
\frac{d[mRNA]}{dt}=\hat\alpha_1+\hat\alpha_2\theta-\hat\delta_1[mRNA], \\ \frac{d[sRRz]}{dt}=\hat\alpha_3[mRNA]-\hat\delta_2[P].where represents the transcription rate when all TcI42 is not associated. By assuming that the concentration of mRNA stablizes instantly, we have
For zero initial value of , integration result is
where and .
Growth
In an environment where the space is finite, the growth rate is restricted by population. For an environment, if the growth rate without restriction is , environmental capacity is , then the change of population follows Logistic model.
For engineering bacteria, temperature-induced suicide system can synthesize sRRz, the concentration of which affects the lysis rate of bacteria. That is, the number of engineering bacteria satisfies
By substituting the result from previous section, we have
where and . For cancer cells, if we use to represent the population, and to represent the concentration of sulforaphane and curcumin, respectively, then
where and represent inhibition intensity of sulforaphane and curcumin, respectively.
Production of sulforaphane and curcumin
We assume that the concentrations of enzymes reach stability instantly, then for sulforaphane and curcumin, there is
where and represent the concentrations of substates synthesizing sulforaphane and curcumin (the synthesis of curcumin is simplified to a one-substrate-one-enzyme reaction), respectively, and represent the maximum speed and Michaelis-Menten constants when the number of engineering bacteria reaches environmental capacity, respectively, and and represent the degradation rate of sulforaphane and curcumin, respectively.
Summary
By summarizing the above results, we obtain that the entire process satisfies the following differential equations:
Numerical simulation
We set the concentration of substrates synthesizing sulforaphane and curcumin to the corresponding Michaelis-Menten constants in order to synthesize sulforaphane and curcumin at half of the maximum velocity. Initially, we set the number of engineering bacteria to a low level, and the number of cancer cells to half of the environmental capacity. For the environment, we set the temperature to 37℃ for the first 7 hours, and after 7 hours the temperature is set to 42℃. We use MATLAB to simulate. The parameter list and the simulation result is shown below.
Table 2: Parameter list
Symbol | Meaning | Unit | Value |
---|---|---|---|
Number of engineering bacteria | None | Initial 10 | |
Number of cancer cells | None | Initial 50000 | |
Concentration of sulforaphane | M | Initial 0 | |
Concentration of curcumin | M | Initial 0 | |
Growth rate of engineering bacteria | min | 0.0083 | |
Growth rate of cancer cells | min | 0.01 | |
Environmental capacity of engineering bacteria | None | 10000 | |
Environmental capacity of cancer cells | None | 100000 | |
Maximum lysis rate when TcI42 is incativated | min | 0.0017 | |
Maximum lysis rate caused by activation of TcI42 | min | 0.1 | |
Activation constant | None | 300 | |
Temperature when half of TcI42 is associated | K | 315 | |
Environmental temperature | K | 300$\rightarrow$315 | |
Degradation rate of sRRz | min | 0.5 | |
Degradation rate of myrosinase | min | 0.0033 | |
Degradation rate of DCS-CURS | min | 0.0008 | |
Time | min | Variable | |
Maximum velocity of myrosinase when the number of engineering bacteria reaches environmental capacity | M/min | 15.51 | |
Maximum velocity of DCS-CURS when the number of engineering bacteria reaches environmental capacity | M/min | 0.2862 | |
Michaelis-Menten constant of myrosinase | M | 84.15 | |
Michaelis-Menten constant of DCS-CURS | M | 71.2 | |
Concentration of substrate synthesizing sulforaphane | M | 84.15 | |
Concentration of substrate synthesizing curcumin | M | 71.2 | |
Inhibition intensity of sulforaphane | (Mmin) | 0.0025 | |
Inhibition intensity of curcumin | (Mmin) | 0.0025 |
![](https://static.igem.wiki/teams/4986/wiki/model/result.jpg)
When the temperature increases to 42℃, the suicide switch turns on, leading to fast death of engineering bacteria. As the degradation of sulforaphane and curcumin needs time, the number of cancer cells will still decrease for some time after increase of temperature. However, after 14h, the concentration of sulforaphane and curcumin will be not low enough to inhibit the growth of cancer cells due to degradation. The simulation result demonstates that the suicide switch and inhibition system work well, whereas the inhibition rate is slow and not thorough.
References
[1] Abedi, M. H., Yao, M. S., Mittelstein, D. R., Bar-Zion, A., Swift, M. B., Lee-Gosselin, A., Barturen-Larrea, P., Buss, M. T., & Shapiro, M. G. (2022). Ultrasound-controllable engineered bacteria for cancer immunotherapy. Nature communications, 13(1), 1585.