Difference between revisions of "Team:TECHNION-ISRAEL/Assay Model"

 
(29 intermediate revisions by 5 users not shown)
Line 149: Line 149:
 
 
 
 
<h4 id="Purpose">The purpose of our model: </h4>
+
<h4 id="Purpose">The purpose of our model </h4>
 
<p>
 
<p>
 
We set out to model the interaction between our engineered model cell lines (HPC-7 /HEK293 cells presenting antigen) and an immature B cell model (WEHI-231 naturally presenting IgM). We chose to model this interaction because data and protocols for this kind of assay are scarce or non-existent. It was our goal to produce valuable information, which would be incorporated into our wet-lab experiment, in an attempt to define more concretely the necessary conditions and concentrations for our assay.  
 
We set out to model the interaction between our engineered model cell lines (HPC-7 /HEK293 cells presenting antigen) and an immature B cell model (WEHI-231 naturally presenting IgM). We chose to model this interaction because data and protocols for this kind of assay are scarce or non-existent. It was our goal to produce valuable information, which would be incorporated into our wet-lab experiment, in an attempt to define more concretely the necessary conditions and concentrations for our assay.  
 
</p>
 
</p>
 
<br>
 
<br>
<h4 id="Does">What our model does: </h4>
+
<h4 id="Does">What our model does </h4>
 
<p>
 
<p>
Our model simulates receptor induced apoptosis in heterogeneous mammalian cell co-culture. Given initial conditions by the user, our model predicts the rate at which the cells within the system will die. While seemingly a simple task, due to the variable nature of cells and receptor presentation combined with non-binary signaling, the system is very complex. Determining the optimal initial concentrations of each cell line as well as the necessary transfection efficiency would require tens of experiments and weeks of work. Using our model, researchers may narrow down the ranges of the aforementioned parameters saving valuable time and money.
+
Our model simulates receptor induced apoptosis in heterogeneous mammalian cell co-culture. Given initial conditions by the user, our model predicts the rate at which the cells within the system will die. While seemingly a simple task, due to the variable nature of cells and receptor presentation combined with non-binary signaling, the system is very complex. Determining the optimal initial concentrations of each cell line as well as the necessary transfection efficiency would require tens of experiments and weeks of work. Using our model, researchers may narrow down the ranges of the aforementioned parameters, saving valuable time and money.
 
</p>
 
</p>
 
 
 
<br>
 
<br>
 
<h4 id="Assumptions">
 
<h4 id="Assumptions">
Assumptions:
+
Assumptions
 
</h4>
 
</h4>
 
 
Line 180: Line 180:
 
<div class = "row">
 
<div class = "row">
 
<div class="col-md-offset-2 col-md-8">
 
<div class="col-md-offset-2 col-md-8">
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/7/79/T--TECHNION-ISRAEL--our_model_reactions_new.png" alt = "" style= "width: 100%; margin:auto;">
+
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/7/7b/%D7%9E%D7%A9%D7%95%D7%95%D7%90%D7%95%D7%AA1.png" alt = "" style= "width: 100%; margin:auto;">
 
<br>
 
<br>
<p style="text-align: center;"><strong> Eq. 1</strong>: Chemical equation describing system  </p>
+
<p style="text-align: center;"><strong> Eq. 1</strong>: Chemical equations describing the system  </p>
 
 
 
<br>
 
<br>
Line 274: Line 274:
 
<tr>
 
<tr>
 
                             <td><span style="font-size: 27px;">ᵞ</span><span style="font-size: 11px;">W</span></td>
 
                             <td><span style="font-size: 27px;">ᵞ</span><span style="font-size: 11px;">W</span></td>
                             <td>1</td>
+
                             <td>0.01</td>
 
</tr>
 
</tr>
  
Line 310: Line 310:
 
</sup>.   
 
</sup>.   
 
<br>
 
<br>
The initial concentrations were determined based on a co-culture experiment that we conducted. In this experiment we aimed to test the viability of HEK cells that were not transfected with the assay plasmid, and thus don't induce apoptosis, and WEHI cells, in co-culture. Different ratios of [WEHI]:[HEK/HPC] cells were tested and the optimal ratio was found to be 6:4 for [WEHI]:[HEK/HPC] respectively.  
+
The initial concentrations were determined based on a <a target="_blank" href="https://2017.igem.org/Team:TECHNION-ISRAEL/Notebook"> co-culture experiment</a> that we conducted. In this experiment we aimed to test the viability of HEK cells that were not transfected with the assay plasmid, and thus don't induce apoptosis, and WEHI cells, in co-culture. Different ratios of [WEHI]:[HEK/HPC] cells were tested and the optimal ratio was found to be 6:4 for [WEHI]:[HEK/HPC] respectively.  
 
</p>
 
</p>
 
<br>
 
<br>
<h4 id="Results">Results and Analysis: </h4>
+
<h4 id="Results">Results and Analysis </h4>
 
<br>
 
<br>
 
</div>
 
</div>
Line 319: Line 319:
 
<br>
 
<br>
 
<div class = "row">
 
<div class = "row">
<div class="col-md-offset-2 col-md-8">
+
<div class="col-md-offset-1 col-md-10">
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/d/d2/T--TECHNION-ISRAEL--our_model_fig_1.jpg" alt = "" style= "width: 100%; margin:auto;">
+
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/3/36/%D7%9E%D7%A6%D7%92%D7%AA.png" alt = "" style= "width: 100%; margin:auto;">
 
<br>
 
<br>
<p style="text-align: center;"><strong> Figure 1: </strong> Constants as listed in Chart 1  </p>
+
<p style="text-align: center;"><strong> Figure 1: </strong> The model simulation results - when using the constants as listed in Chart 1  </p>
 
</div>
 
</div>
 
</div>
 
</div>
Line 329: Line 329:
 
<div class =  "col-md-12" >
 
<div class =  "col-md-12" >
 
<p>
 
<p>
As can be seen in figure 1, our model is in line with our expectations. We can see that the WEHI cells undergo apoptosis as a result of co-culture with the antigen displaying HEK/HPC cells, and in direct response to the increase of H-W complexes.
+
As can be seen in <strong>figure 1</strong>, our model is in line with our expectations. We can see that the WEHI cells undergo apoptosis as a result of co-culture with the antigen displaying HEK/HPC cells, and in direct response to the increase of H-W complexes.
 
<br>
 
<br>
 
In order to further characterize our system, we tested our model using different constants and cell ratios. We discovered that setting the value of  
 
In order to further characterize our system, we tested our model using different constants and cell ratios. We discovered that setting the value of  
 
K<span style="font-size: 11px;">a</span>  equal to
 
K<span style="font-size: 11px;">a</span>  equal to
K<span style="font-size: 11px;">H</span>  yields a second behavioral regime in the model:
+
K<span style="font-size: 11px;">H</span>  yields a second behavioral regime in the model <strong>(figure 2)</strong>:
  
 
</p>
 
</p>
Line 343: Line 343:
 
<br>
 
<br>
 
<div class = "row">
 
<div class = "row">
<div class="col-md-offset-2 col-md-8">
+
<div class="col-md-offset-1 col-md-10">
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/9/9c/T--TECHNION-ISRAEL--our_model_fig_2_new.jpg" alt = "" style= "width: 100%; margin:auto;">
+
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/b/b5/%D7%A0%D7%AA%D7%95%D7%9F_2.png" alt = "" style= "width: 100%; margin:auto;">
 
<br>
 
<br>
<p style="text-align: center;"><strong> Figure 2: </strong>K<span style="font-size: 11px;">a</span>=K<span style="font-size: 11px;">H</span>  </p>
+
<p style="text-align: center;"><strong> Figure 2:</strong> The model simulation results when K<span style="font-size: 11px;">a</span>=K<span style="font-size: 11px;">H</span>  </p>
 
</div>
 
</div>
 
</div>
 
</div>
Line 360: Line 360:
  
 
<p>
 
<p>
We also found  that our model is robust and shows the same behavior for a wide range of initial concentration ratios-  a decrease in WEHI cells to near zero. For example, setting the ratio of [WEHI]:[HEK/HPC-7] to 4:6 yielded the following graph:
+
We also found  that our model is robust and shows the same behavior for a wide range of initial concentration ratios-  a decrease in WEHI cells to near zero. For example, setting the ratio of [WEHI]:[HEK/HPC-7] to 4:6 yielded the following graph <strong>(figure 3)</strong>:
 
</p>
 
</p>
 
 
Line 368: Line 368:
 
<br>
 
<br>
 
<div class = "row">
 
<div class = "row">
<div class="col-md-offset-2 col-md-8">
+
<div class="col-md-offset-1 col-md-10">
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/9/94/T--TECHNION-ISRAEL--our_model_fig_3.jpg" alt = "" style= "width: 100%; margin:auto;">
+
<img class="no-rad" src="https://static.igem.org/mediawiki/2017/4/41/%D7%A0%D7%AA%D7%95%D7%9F_3.png" alt = "" style= "width: 100%; margin:auto;">
 
<br>
 
<br>
<p style="text-align: center;"><strong> Figure 3:</strong> [WEHI]:[HEK/HPC] ratio of 4:6  </p>
+
<p style="text-align: center;"><strong> Figure 3:</strong> The model simulation results for the initial concentrations [WEHI]:[HEK/HPC] ratio of 4:6  </p>
 
</div>
 
</div>
 
</div>
 
</div>
Line 380: Line 380:
 
<div class =  "col-md-12" >
 
<div class =  "col-md-12" >
 
<p>
 
<p>
As the model predicted excellent results using a [WEHI]:[HEK] ratio of 6:4, we conducted the tolerance assay experiment using this ratio. Our RESULTS established our model’s assumptions and demonstrated a qualitatively equivalent outcome.
+
As the model predicted excellent results using a [WEHI]:[HEK] ratio of 6:4, we conducted the  
 +
tolerance assay experiment using this ratio. Our  
 +
<a target="_blank" href="https://2017.igem.org/Team:TECHNION-ISRAEL/Results#Proofofconcept">Results</a> established our model’s assumptions and demonstrated a qualitatively equivalent outcome.
 
</p>
 
</p>
 
 
Line 392: Line 394:
 
<div class= "references">
 
<div class= "references">
 
<ol>
 
<ol>
<li id="ref1"> STEVEN M. RUDICH, KENNETH H. ROUX, ROBERTJ. WINCHESTER, AND PATRICIA K. A. MONGINI, ANTI-IgM-MEDIATED B CELL SIGNALING, <i>The Rockefeller University press</i>, Volume 168, 247-266, July 1988</li>
+
<li id="ref1"> Rudich, S. M., et al. "Anti-IgM-mediated B cell signaling. Molecular analysis of ligand binding requisites for human B cell clonal expansion and tolerance." <i>Journal of Experimental Medicine</i> 168.1 (1988): 247-266.</li>
<li id="ref2">] Rodolfo Molina-Peña and Mario Moisés Álvarez, A Simple Mathematical Model Based on the Cancer Stem Cell Hypothesis Suggests Kinetic Commonalities in Solid Tumor Growth, <i>PLoS One.</i> 7(2): e26233.2012; </li>
+
<li id="ref2"> Molina-Peña, Rodolfo, and Mario Moisés Álvarez. "A simple mathematical model based on the cancer stem cell hypothesis suggests kinetic commonalities in solid tumor growth." <i>PloS one</i> 7.2 (2012): e26233.  </li>
 
</ol>
 
</ol>
 
</div>
 
</div>

Latest revision as of 13:59, 12 December 2017

navbar

Our Model

Tolerance Assay Model



The purpose of our model

We set out to model the interaction between our engineered model cell lines (HPC-7 /HEK293 cells presenting antigen) and an immature B cell model (WEHI-231 naturally presenting IgM). We chose to model this interaction because data and protocols for this kind of assay are scarce or non-existent. It was our goal to produce valuable information, which would be incorporated into our wet-lab experiment, in an attempt to define more concretely the necessary conditions and concentrations for our assay.


What our model does

Our model simulates receptor induced apoptosis in heterogeneous mammalian cell co-culture. Given initial conditions by the user, our model predicts the rate at which the cells within the system will die. While seemingly a simple task, due to the variable nature of cells and receptor presentation combined with non-binary signaling, the system is very complex. Determining the optimal initial concentrations of each cell line as well as the necessary transfection efficiency would require tens of experiments and weeks of work. Using our model, researchers may narrow down the ranges of the aforementioned parameters, saving valuable time and money.


Assumptions

  1. The growth rate of HEK293 cells is similar to the growth rate of WEHI cells
  2. WEHI-231 cells undergo apoptosis as a result of interaction with antigen presenting HEK/HPC cells
  3. Kapop > Ka

In order to accurately model our system, we designed a kinetic model based on the following equations:




Eq. 1: Chemical equations describing the system




Eq. 2: Set of three ODEs describing the interactions within the system


Where:

  • H represents our model cell lines HEK293/HPC-7
  • W represents the WEHI cells
  • H-W represents HEK cells bound to WEHI cells
  • KH is the proliferation rate constant of HEK/HPC cells
  • H is the death rate constant of HEK/HPC cells
  • KW is the proliferation rate constant of WEHI cells
  • W is the death rate constant of WEHI cells
  • Ka is the binding rate constant for combination of HEK and WEHI to form H-W
  • Kd is the rate constant for disassociation of H-W to form HEK and WEHI
  • Kapop is the rate constant for the apoptosis of WEHI cells as a result of binding HEK cells

The equations were solved by MATLAB and the simulation was done using the following values for constants and initial concentrations:


Chart 1 – values of constants and initial concentrations


The ratios between the constants were defined based on the ratios found in relevant literature [1] [2] .
The initial concentrations were determined based on a co-culture experiment that we conducted. In this experiment we aimed to test the viability of HEK cells that were not transfected with the assay plasmid, and thus don't induce apoptosis, and WEHI cells, in co-culture. Different ratios of [WEHI]:[HEK/HPC] cells were tested and the optimal ratio was found to be 6:4 for [WEHI]:[HEK/HPC] respectively.


Results and Analysis




Figure 1: The model simulation results - when using the constants as listed in Chart 1


As can be seen in figure 1, our model is in line with our expectations. We can see that the WEHI cells undergo apoptosis as a result of co-culture with the antigen displaying HEK/HPC cells, and in direct response to the increase of H-W complexes.
In order to further characterize our system, we tested our model using different constants and cell ratios. We discovered that setting the value of Ka equal to KH yields a second behavioral regime in the model (figure 2):




Figure 2: The model simulation results when Ka=KH


For this value of Ka, and for a wide range of smaller values, the same behavior was observed. We can conclude that when Ka is higher than KH, at least 50% of HEK cells must bind WEHI cells in order to reduce the amount of WEHI cells to a negligible level.


We also found that our model is robust and shows the same behavior for a wide range of initial concentration ratios- a decrease in WEHI cells to near zero. For example, setting the ratio of [WEHI]:[HEK/HPC-7] to 4:6 yielded the following graph (figure 3):




Figure 3: The model simulation results for the initial concentrations [WEHI]:[HEK/HPC] ratio of 4:6


As the model predicted excellent results using a [WEHI]:[HEK] ratio of 6:4, we conducted the tolerance assay experiment using this ratio. Our Results established our model’s assumptions and demonstrated a qualitatively equivalent outcome.




  1. Rudich, S. M., et al. "Anti-IgM-mediated B cell signaling. Molecular analysis of ligand binding requisites for human B cell clonal expansion and tolerance." Journal of Experimental Medicine 168.1 (1988): 247-266.
  2. Molina-Peña, Rodolfo, and Mario Moisés Álvarez. "A simple mathematical model based on the cancer stem cell hypothesis suggests kinetic commonalities in solid tumor growth." PloS one 7.2 (2012): e26233.
My First Website