<p class="description"><strong>Conclusions:</strong>From this simulation, we can conclude that our model simulates the functioning of our circuit as expected and intended. We use the parameters that we obtained from fitting of our experimental data and tuning them with regards to the intended application context. The above results verify the functioning of the Tumor Sensing circuit and thus validate CATE's behavior in a tumor.</p>
−
<p class="description"></p>
+
</section>
</section>
Revision as of 22:44, 1 November 2017
Modelling the Behavior of CATE inside Tumor
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 our 3D 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. We use the parameters obtained by fitting our experimental data and also finally integrating our self-designedhybrid promoter.
Geometry
Our model structure consists of 2 domains - Tumor and Layer.
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.
Figure 1: Geometry of the tumor and bacteria colony (green area: colonized by E. coli Nissle)
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.
Figure 2: Geometry of the COMSOL Model - 2D Axisymmetric; Coordinate system: (r,z); axes in mm; 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.
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, and bacterial layer: production and intracellular degradation.
Expand the sections below for more details about the domain-wise reaction rates for each species (AHL, LuxI and Azurin) and equations used.
Disclaimer
The concentration of LuxI, LuxR and LuxR-AHL modelled above represents the intracellular concentration i.e. the concentration at a point had a bacteria been there.
However, since AHL diffuses everywhere rapidly and freely through the cell-memberane, a rapid equilibrium between the intercellular and extracellular concentration of AHL is reached, and thus we multiply the LuxI controlled AHL production rate by dcell, since the layer has all the bacteria colonizing the tumor.
Similarly, since we are interested in the Azurin concentration after lysis, we have a dcell multiplication factor in its reaction rate, as after lysis Azurin from all the cells diffused out of the layer. This is also then used to find the effective killing area to estimate the percentage of volume of the tumor treated.
Growth Model
An exponential growth model is used, as shown below.
Transport of Diluted Species physics of COMSOL was used to model the diffusion of AHL and Azurin. The equations below represent the diffusion models in both the Tumor and Layer domains. For simulating the effect of lysis, the diffusion of Azurin is triggered when temperature reaches 42°C, while AHL diffuses all the time. The value of the diffusion coefficients of AHL, D and Azurin, DAzu are given in the parameter list below.
The parameters that were used in our 3D 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.
Expand the details section to read more about the different parameters and their values used.
The steady state cell density for growth, dcell,ss is chosen to be 0.05 in colonized tumor and 0.0005 in healthy tissue, based on the conclusions derived in the System Specifications. Also, [Lac] is chosen to be 1 mM for a healthy tissue and 5 mM for a tumor, and the initial cell density dcell,0 is taken to be a very small non-zero value.
Parameters fitted to our experimental data and tuned with respect to the context of the intended application
To test our model, we simulated it for a given set of parameters that describe the conditions of tumor colonization as close as possible to a reality. This is given bsy the relative criteria for Azurin production. The 3D model enabled us to use exact partial differential equations to model diffusion of the different species and helped us bridge the gap between experiments and reality, as it was not feasible to conduct experiments for tumor colonization in the lab. This also helped us to validate the assumptions used in the in-vitro model.
Our 3D model helped us to simulate the three main phases of the CATE treatment: Growth, AND-Gate Switching (Environment sensing) and finally Lysis & Azurin Diffusion. For details see CATE in Action.
Phase 1: Growth
The growth is modelled using the exponential growth rate \[\frac{\mathrm{d} d_{\text{cell}}}{\mathrm{d} t} = \frac{1}{\tau} d_{\text{cell}} \left ( 1-\frac{d_{\text{cell}}}{d_{\text{cell,ss}}} \right )\], where τ is the doubling time, that we obtained by fitting E. coli Nissle growth curves from our lab experiments. Read more in details for the equations.
Phase 2: AND-Gate Tumor Sensing Switch
The AND gate Switch senses the environment based on Quorum sensing (dcell) and Lactatae concentration ([Lac]), as explained in detail in the description of the Tumor Sensing circuit. Based on the different combinations of dcell and [Lac], as mentioned in the system specifications, the AND gate switches 'ON' or 'OFF'. For detailed equations see the equation details.
ON state:
Rapid and high-fold increase in [LuxI] or [Azu]
OFF state:
Slow and negligible-fold increase in LuxI] or [Azu]
Phase 3: Lysis and Azurin Diffusion
The effect of lysis is simulated by a temperature controlled trigger of the diffusion of azurin that is produced in the cells. The temperature is increased from 0 to 42°C as a step function and when the temperature reaches 42°C, the production of Azurin and AHL stops, since the cells are lysed, and diffusion of Azurin begins which depletes Azurin out of the cell into the tumor, effectively inducing apoptosis of the tumor cells.
We simulated our model to test the working of the Tumor Sensing circuit for the case of tumor colonization i.e. High dcell AND High [Lac]. The three phases, as described above were simulated and the results for the case of bacteria colonization of tumor are shown in Figures 3, 4 and 5. In Figure 6, dcell shows the growth of the cell density inside the layer and PLux-Lac shows the switching based on the environmental conditions of dcell and lactate.
During the growth phase, our sensing circuit is OFF (visible in Figure 3), and is triggered to turn ON once the desired cell density for quorum sensing is reached at around 40 hr, as shown in Figure 6 and 2. Then once steady state is reached, the temperature step triggers the cell-lysis and stops the production of Azurin and AHL. Finally as is visible in Figure 5, all the Azurin diffuses out of the layer very rapidly, thus completing the last phase of the treatment.
Figure 3: Growth phase: 0 hr to 35 hr; Colonization of the tumor by bacteria.Figure 4: AND-Gate Switch Sensing phase: 35 hr to 69 hr; The rapid increase in [Azu] shows the switch ON of the AND-gate tumor sensor.Figure 5: Cell Lysis and Azurin Diffusion phase: 69 hr to 100 hrFigure 6: Normalized concentration of AHL and Azurin, Cell density (as a ratio of its steady state value), PLux-Lac and Temperature (as a ratio of its steady state value 42°C) probe plots at a point inside the layer of bacterial colonization in the tumor.
Conclusions:From this simulation, we can conclude that our model simulates the functioning of our circuit as expected and intended. We use the parameters that we obtained from fitting of our experimental data and tuning them with regards to the intended application context. The above results verify the functioning of the Tumor Sensing circuit and thus validate CATE's behavior in a tumor.
AND-Gate Tumor Sensor Characterization
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.verify the behavior of our tumor-sensing circuit in more real-life conditions pertaining to the intended application context of a solid spherical tumor.To test our sensing circuit with the intended application in mind, we use the relative criteria set in the system specifications. Since our AND-Gate switch has 2 inputs for environment sensing viz. dcell and [Lac], there are 4 possible binary combinations that CATE can encounter in real-life scenario, viz.:
Tumor Colonization
High dcell AND High [Lac]
Healthy tissue Colonization
High dcell AND Low [Lac]
Tumor NOT colonized
Low dcell AND High [Lac]
Healthy tissue NOT colonized
Low dcell AND Low [Lac]
The simulation results between 35 hr and 75 hr for all the 4 possible scenarios that CATE can encounter, are shown below.
Figure 6: Case: Tumor colonized - High dcell AND High [Lac]; AND-gate Switch is turned ON fully.Figure 7: Case: Healthy tissue colonized - High dcell AND Low [Lac]; AND-gate Switch is turned ON partially.Figure 8: Case: Tumor NOT colonized - Low dcell AND High [Lac]; AND-gate Switch is turned OFF fully.Figure 9: Case: Healthy tissue NOT colonized - Low dcell AND Low [Lac]; AND-gate Switch is turned OFF fully.
As is clear from the Figures 6-9, the AND-gate switch is fully turned ON for the case of tumor colonization with steady state Azurin concentration reaching 75000 nM (as shown in Figure 10). For the case of colonization of healthy tissue, the AND-gate switch is partially turned ON with steady state Azurin concentration at about 22000 nM. This is about 3.5 times lower than when the switch is fully ON. In the case of tumor and healthy tissue not colonized the steady state Azurin concentrations are around 20 nM and 17 nM, respectively, which represents a fully turned OFF AND-gate switch.
Figure 10: [Azu] vs Time with linear scale; Shows the behavior of our tumor sensing circuit in the different dcell and [Lac] conditions.
Killing Area
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.
Using the COMSOL data, the maximum Azurin concentration over time was found for every point in the system, and the results are as shown in Figure 15.
Figure 15: Maximum [Azu] (nM) over time at all points.Figure 16: Killing Area: Points where maximum [Azu] over time exceeds 20 µM.
From Figure 15, the effective killing area can be estimated by choosing the points in the (r,z) plane that have [Azu]max|t > 20 µM. Figure 16 shows the maximum effective killing area is at z = 0, between 9.1 and 10.1 mm. This amounts to approximately 3.5 % of tumor being treated by volume. So we can conclude, that a more potent cytotoxic agent that Azurin is required to achieve significant tumor treatment using CATE.
Bacterial colonization patterns
To further push the limits of our system, we tested the environment sensing AND-gate switching functioning in different colonization patterns. Although any possible colonization pattern can be simulated, we simulated the logically closest patterns (compared to a spherical shell layer), that we could think of, and introduced heterogenity through partitions.
Case 1: Homogeneous distribution in a Single spherical shell layer in tumour
already simulated as shown in previous results
Case 2: Heterogeneous distribution (Partitions) in a Single spherical shell layer in tumour
another possible feasible scenario
Case 3: Heterogeneous distribution (Partitions) in Double spherical shell layer in tumour
introducing some more heterogenity in the system
Case 4: Homogeneous distribution throughout Healthy tissue
with dcell = 0.0005, as in this geometry the AND-gate is more easily activated as comparaed to other geometries, since difusion can not occur within the inner region where bacteria colonize, and thus this colonized inner region acts as a huge source, with no diffusion between the cells, as shown in the Figure 15.
The simulation results between 35 hr and 75 hr for all the 4 possible scenarios, are shown below.
Figure 11: Case 1: Tumor colonized in a single spherical shell-shaped layer; AND-gate Switch is turned ON fully.Figure 12: Case 2: Tumor is colonized in partitions in the shape of a single shell-shaped layer; AND-gate Switch is turned ON fully.Figure 13: Case 3: Tumor is colonized in alternative paritions in two spherical shell-shaped layers; AND-gate Switch is turned ON fully.Figure 14: Case 4: Homogeneous colonization (dcell) of a healthy tissue; AND-gate switch turns ON but only partially.
The above figures show that our Tumor sensing circuit functions properly in the different colonization patterns simulated. According to Figure 14, in case of colonization of a healthy tissue, the switch is only partially triggered (slower rise and lesser final steady state [Azu]), when compared to Case 1 in Figure 12, that has a fully ON switch indicated by the rapid rise and approximately 3.5 times higher-fold increase in [Azu]
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 of our Model
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.
Jordi Garcia-Ojalvo, Michael B. Elowitz, and Steven H. Strogatz Modeling a synthetic multicellular clock: Repressilators coupled by quorum sensing PNAS 2004 101 (30) 10955-10960 | doi:10.1073/pnas.0307095101
A.B. Goryachev, D.J. Toh T.Lee, Systems analysis of a quorum sensing network: Design constraints imposed by the functional requirements, network topology and kinetic constants Biosystems, Volume 83, Issues 2–3, February–March 2006, Pages 178-187 | doi:10.1016/j.biosystems.2005.04.006
A synthetic multicellular system for programmed pattern formation Subhayu Basu, Yoram Gerchman, Cynthia H. Collins, Frances H. Arnold & Ron WeissNature 434, 1130-1134 (28 April 2005) | doi:10.1038/nature03461