Line 11: | Line 11: | ||
<h1 class="headline">Modelling the Behavior of CATE inside Tumor</h1> | <h1 class="headline">Modelling the Behavior of CATE inside Tumor</h1> | ||
− | + | <figure class="fig-nonfloat"> | |
<img src="https://static.igem.org/mediawiki/2017/3/3e/T--ETH_Zurich--InVivo_Head.png" alt="FIXME"> | <img src="https://static.igem.org/mediawiki/2017/3/3e/T--ETH_Zurich--InVivo_Head.png" alt="FIXME"> | ||
− | </figure | + | </figure> |
<section class="first"> | <section class="first"> | ||
<!--<h1>What</h1>--> | <!--<h1>What</h1>--> | ||
− | <p><em>We developed a model to gauge the behavior of our | + | <p><em>We developed a model to gauge the behavior of our circuit in the real life conditions of solid tumor colonization.</em></p> |
+ | <p>This page gives details about the 3D model we developed to simulate the behavior of CATE inside the tumor and healthy tissue with different colonization patterns. This section presents details about the following: </p> | ||
+ | <ul> | ||
+ | <li> | ||
+ | The Geometry used in the COMSOL model.<br> | ||
+ | <p class="description"></p> | ||
+ | </li> | ||
+ | <li> | ||
+ | The domain-wise Partial Differential Equations used for modelling the Growth, AND-gate Tumor sensing Switch and the Lysis & Diffusion effects.<br> | ||
+ | <p class="description"></p> | ||
+ | </li> | ||
+ | <li> | ||
+ | Values of the Parameters used<br> | ||
+ | <p class="description"></p> | ||
+ | </li> | ||
+ | <li> | ||
+ | Simulation Results for the AND-Gate tumor sensing switch in the different environmental conditions to validate its functioning in context of the intended application.<br> | ||
+ | <p class="description"></p> | ||
+ | </li> | ||
+ | </ul> | ||
+ | <p>The 3D model presented here is used to model the behavior in conditions as close as possible to real-life scenario of tumor colonization by including the partial differential equations for diffusion of AHL and Azuin, which also helps validate the assumptions and simplifications used in the <a href="https://2017.igem.org/Team:ETH_Zurich/Model/Environment_Sensing/model"><em>in-vitro</em> model</a>. Since in experiments, there is no diffusion, our 3D model helps us model scenarios that CATE will encounter, as close as possible to reality. We are also able to test the behavior of our <a href="https://2017.igem.org/Team:ETH_Zurich/Experiments/Tumor_Sensor">tumor sensor</a>.</p> | ||
</section> | </section> | ||
<section> | <section> | ||
− | <h1> | + | <h1>Geometry</h1> |
− | + | <p class="description">As mentioned in <a href="https://2017.igem.org/Team:ETH_Zurich/Model/Environment_Sensing/system_specifications">system specifications</a>, the tumour has been chosen as a solid sphere of radius 20 mm and the bactierial colonization pattern as a homogenous distribution in a spherical shell-shaped 0.5 mm thick layer in the tumour at a distance of 10 mm from the centre of the tumor, as shown in Figure 1.</p> | |
− | + | ||
− | + | ||
− | <p class="description">As mentioned in <a href="https://2017.igem.org/Team:ETH_Zurich/Model/Environment_Sensing/system_specifications">system specifications</a>, the tumour has been chosen as a solid sphere of radius | + | |
<figure class="fig-nonfloat" style="max-width: 600px;"> | <figure class="fig-nonfloat" style="max-width: 600px;"> | ||
Line 32: | Line 49: | ||
<figcaption>Figure 1: Geometry of the tumor and bacteria colony (green area: colonized by <span class="bacterium">E. coli</span> Nissle)</figcaption> | <figcaption>Figure 1: Geometry of the tumor and bacteria colony (green area: colonized by <span class="bacterium">E. coli</span> Nissle)</figcaption> | ||
</figure> | </figure> | ||
+ | |||
+ | <div class="multi-summary"> | ||
+ | <details> | ||
+ | <summary>Details about the geometry</summary> | ||
+ | <p class="description">Due to the spherical symmetry of the system, a 2D axisymmetric COMSOL model was used as shown in Figure 2 - a semicircle of radius 20 mm represents the tumor and the 0.5 mm thick layer at a distance of 10 mm from the center of the tumor represents the bacteria colonization pattern as explained in the <a href="https://2017.igem.org/Team:ETH_Zurich/Model/Environment_Sensing/system_specifications">system specifications</a>. COMSOL then sweeps the semi-circle between 0° and 360 ° to simulate the entire 3D problem. The symmetry helps decrease the computational time and space requirements, without having to apply simplifications based on assumptions.</p> | ||
+ | <figure class="fig-nonfloat"> | ||
+ | <img src="https://static.igem.org/mediawiki/2017/a/a5/T--ETH_Zurich--Geometry.png" alt="FIXME"> | ||
+ | <figcaption>Figure 2: Geometry of the COMSOL Model - 2D Axisymmetric; <em>Coordinate system: (r,z); axes in mm</em>; Tumor is a solid sphere of radius 20 mm located at (0,0), Bacterial colonization is in a spherical layer of thickness 0.5 mm at a distance of 10 mm from the center of the tumor.</figcaption> | ||
+ | </figure> | ||
+ | </details> | ||
+ | </div> | ||
+ | |||
+ | <p class="description emphasize"><em>Thus our model structure consists of 2 domains - <strong>Tumor</strong> and <strong>Layer</strong>. In these two domains, the reaction rates of the different species is different as is explained by the equations below.</em> </p> | ||
+ | |||
+ | </section> | ||
<h2>Equations <button><a href="https://2017.igem.org/Team:ETH_Zurich/Model/In_Vivo_Detailed#InVivo_Equations" class="more">Details</a></button> </h2> | <h2>Equations <button><a href="https://2017.igem.org/Team:ETH_Zurich/Model/In_Vivo_Detailed#InVivo_Equations" class="more">Details</a></button> </h2> | ||
Line 45: | Line 77: | ||
<section> | <section> | ||
− | <h1> | + | <h1>Results</h1> |
<!--<p>The COMSOL model helped us to extend our MATLAB model in the following ways:</p>--> | <!--<p>The COMSOL model helped us to extend our MATLAB model in the following ways:</p>--> | ||
<h2>We could simulate for a geometry of the system closer to the real-life tumor conditions <button><a href="https://2017.igem.org/Team:ETH_Zurich/Model/In_Vivo_Detailed#InVivo_Treatment" class="more">Simulating the Treatment</a></button></h2> | <h2>We could simulate for a geometry of the system closer to the real-life tumor conditions <button><a href="https://2017.igem.org/Team:ETH_Zurich/Model/In_Vivo_Detailed#InVivo_Treatment" class="more">Simulating the Treatment</a></button></h2> |
Revision as of 20:44, 1 November 2017
![CATE logo](https://static.igem.org/mediawiki/2017/d/d6/T--ETH_Zurich--LogoMinimal.png)
Modelling the Behavior of CATE inside Tumor
![FIXME](https://static.igem.org/mediawiki/2017/3/3e/T--ETH_Zurich--InVivo_Head.png)
We developed a model to gauge the behavior of our circuit in the real life conditions of solid tumor colonization.
This page gives details about the 3D model we developed to simulate the behavior of CATE inside the tumor and healthy tissue with different colonization patterns. This section presents details about the following:
-
The Geometry used in the COMSOL model.
-
The domain-wise Partial Differential Equations used for modelling the Growth, AND-gate Tumor sensing Switch and the Lysis & Diffusion effects.
-
Values of the Parameters used
-
Simulation Results for the AND-Gate tumor sensing switch in the different environmental conditions to validate its functioning in context of the intended application.
The 3D model presented here is used to model the behavior in conditions as close as possible to real-life scenario of tumor colonization by including the partial differential equations for diffusion of AHL and Azuin, which also helps validate the assumptions and simplifications used in the in-vitro model. Since in experiments, there is no diffusion, our 3D model helps us model scenarios that CATE will encounter, as close as possible to reality. We are also able to test the behavior of our tumor sensor.
Geometry
As mentioned in system specifications, the tumour has been chosen as a solid sphere of radius 20 mm and the bactierial colonization pattern as a homogenous distribution in a spherical shell-shaped 0.5 mm thick layer in the tumour at a distance of 10 mm from the centre of the tumor, as shown in Figure 1.
![Geometry of tumor and bacterial colony](https://static.igem.org/mediawiki/2017/3/36/T--ETH_Zurich--tumor_geometry.png)
Details about the geometry
Due to the spherical symmetry of the system, a 2D axisymmetric COMSOL model was used as shown in Figure 2 - a semicircle of radius 20 mm represents the tumor and the 0.5 mm thick layer at a distance of 10 mm from the center of the tumor represents the bacteria colonization pattern as explained in the system specifications. COMSOL then sweeps the semi-circle between 0° and 360 ° to simulate the entire 3D problem. The symmetry helps decrease the computational time and space requirements, without having to apply simplifications based on assumptions.
![FIXME](https://static.igem.org/mediawiki/2017/a/a5/T--ETH_Zurich--Geometry.png)
Thus our model structure consists of 2 domains - Tumor and Layer. In these two domains, the reaction rates of the different species is different as is explained by the equations below.
Equations
Transport of Diluted Species physics was used in COMSOL to integrate diffusion into our model. The partial differential equation for diffusion of a species C with reaction source rate RC is \[\frac{\partial \text{[C]}}{\partial t} + \nabla \cdot (-D_{\text{C}} \nabla \text{[C]})= R_{\text{C}}\]. The reaction rates of the species depends on the domain – tumor (no production and only extracellular degradation) or bacterial layer (production and intracellular degradation). Read here about the details of the domain-wise reaction rates for each species (AHL, LuxI and Azurin) and equations used.
Parameters
The parameters that were used in the COMSOL model were obtained partly from literature, partly from characterizations of previous iGEM teams and finally the most important ones were estimated by fitting our experimental data and tuning the fitted-results in the context of the intended applciation, as explained in detail by the Functional Parameter Search. Check out details about the model to read more about the different parameter values used.
For more details about the model go to the detailed description and Functional Parameter Search.
Results
We could simulate for a geometry of the system closer to the real-life tumor conditions
Since it was not practically feasible to conduct experiments of bacterial colonization inside tumors, we simulated the bacterial colonization in a thin spherical layer inside a solid tumor considering the simplifications and assumptions as mentioned in the system specifications. This helped us to test our tumour sensing AND-Gate switch functionality in all the possible real-life scenarios that CATE might encounter in context of the intended application.
Exact diffusion physics of AHL was included witout any simplifications
Our MATLAB model uses a simplified AHL diffusion model with the assumption of negligible degradation inside the layer and and not taking into consideration the diffusion of AHL far from the source. Extending the diffusion physics ordinary differential equations into partial differential equations using the COMSOL model helped us gauge and verify the behavior of our tumor-sensing circuit in more real-life conditions pertaining to the intended application context of a solid spherical tumor. Using the results obtained from our simulations, we could check the behavior of the AND Gate Switching in different conditions of dcell and lactate.
Diffusion physics of Azurin was included to simulate the effect of lysis
To simulate the effect of lysis, our COMSOL model stops the production of Azurin and starts its diffusion when temperature reaches 42°C. This simulates the effect of increase in temperature with FUS to cause cell lysis. Using data obtained from such a simulation, we could also find the temporal-maximum concentrations of Azurin at each point in the tumor, effectively helping us to estimate the killing area and the time-scale of the treatment.
Simulation of different colonization patterns
Using our model, we also tried a few other colonization patterns to show our system works as expected inside a tumor while stays dormant in healthy tissue. We simulated the following patterns:
-
Homogeneous distribution in a Single spherical-shell-shaped layer in Tumor
-
Heterogeneous distribution in a Single spherical-shell-shaped layer in Tumor
-
Heterogeneous distribution in Double spherical-shell-shaped layer in Tumor
-
Homogeneous distribution in Healthy tissue
Final Conclusions
Our model finally helped us implement a comprehensive in-silico test to prove that our system already exhibits an excellent performance for the clinical application it has been designed for. This helped us to take into account the fitting of the parameters using our experimental data and integrate our own hybrid promoter tuned in reference to the intended application context, to verify the functioning of our environment sensing function.
Limitations
Our model has some limitations. We do not model protein E production and cell lysis caused by it. Instead lysis is just simulated in effect as the end of production of AHL and Azurin and start of diffusion of Azurin. Moreover, a step signal is used as a trigger for the lysis. As mentioned in the parameters description, Azurin production is taken to be 10 times proportional to LuxI production. Also, killing mechanism of Azurin has not been modelled since that was not necessary to demonstrate the working of our project CATE in the scope of iGEM.
Tools used
- COMSOL Multiphysics 5.2a by COMSOL Inc.
- MATLAB R2016b by MathWorks