Difference between revisions of "Team:Munich/Model"

Line 119: Line 119:
 
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/8/87/T--Munich--ModellingPagePicture_kinetics_Cas13a_only.png" width="700">
 
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/8/87/T--Munich--ModellingPagePicture_kinetics_Cas13a_only.png" width="700">
 
<p>
 
<p>
Figure 1: Kinetics of Cas13a using 1 nM and 10 nM at different target concentrations.  
+
Figure 1: Kinetics of Cas13a using 1 nM Cas13a and 10 nM crRNA at different target concentrations.  
 
</p>
 
</p>
 
</div>
 
</div>
Line 132: Line 132:
 
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/b/b8/T--Munich--ModellingPagePicture_CascAID.png" width="600">
 
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/b/b8/T--Munich--ModellingPagePicture_CascAID.png" width="600">
 
<p>
 
<p>
Figure 2: Estimated detection limit determined for the Cas13a system using 20 nM concentrations of Cas13a and crRNA.  
+
Figure 2: Estimated detection limit determined for the Cas13a system using 1 nM Cas13a and 10 nM crRNA.  
 
</p>
 
</p>
 
</div>
 
</div>
  
<div class="captionPicture">
 
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/1/14/T--Munich--ModellingPagePicture_kinetics_Cas13a.png" width="600">
 
<p>
 
Figure 3: Kinetics of the  Cas13a system using 1 nM concentrations of Cas13a and  10 nM crRNA at different target concentrations using the reaction cascade.
 
</p>
 
</div>
 
  
 
</td>
 
</tr>
 
 
<tr><td colspan=6 align=center valign=center>
 
 
<h3>Estimating the Detection Limit of the Reaction Cascade</h3>
 
<h3>Estimating the Detection Limit of the Reaction Cascade</h3>
 
<p>   
 
<p>   
Line 166: Line 155:
 
</p>
 
</p>
 
</div>
 
</div>
 +
 +
 +
 +
<div class="captionPicture">
 +
<img alt="LightbringerReal" src="https://static.igem.org/mediawiki/2017/1/14/T--Munich--ModellingPagePicture_kinetics_Cas13a.png" width="600">
 +
<p>
 +
Figure 3: Kinetics of the  Cas13a systemusing 1 nM Cas13a and 10 nM crRNA at different target concentrations using the reaction cascade.
 +
</p>
 +
</div>
 +
 +
 +
</td>
 +
</tr>
 +
 +
<tr><td colspan=6 align=center valign=center>
  
  

Revision as of 03:27, 2 November 2017


Modelling

Modelling in Biosciences is a powerful tool that allows us to get a deeper understanding of our system. It guided and assisted the design of our detection system which helped us saving a lot of time by avoiding dead-end designs. We used simple models to simulate the kinetics of our enzyme cascade to develop intuition for the design of experiments. We then used our models to optimize the design of our reaction cascades for an improved detection limit and optimal lysis times. Our scripts can be found on our GitHub repository.

Detection Limit

One major concern when dealing with the problem of diagnostics on patients is extracting the sample with which detection can actually be performed. Since we wanted our method to be non-invasive, one concern that we needed to deal with is the concentration of pathogens and thus detectable RNA in the patients mucus or other non-invasive sample. First approximations from different papers already showed that virological samples show concentrations no higher than low pM and can even go as low as aM.

Our wetlab experiments indicated that the detection limit of the Cas13a RNase activity is in the range of 10 nM. Using our kinetic data, we estimated the rate constants for the different reactions to create a simple ODE model.
The chemical and differential equations for the model are shown below:

LightbringerReal

As shown in Figure 1, our simulations are able to reproduce the behavior observed experimentally.

LightbringerReal

Figure 1: Kinetics of Cas13a using 1 nM Cas13a and 10 nM crRNA at different target concentrations.

Next, we analysed the amount of readout RNA that was cleaved after 30 minutes for varying target concentrations. As shown in Figure 2, the curve follows a sigmoidal behavior and suggests a detection limit in the range of 10 nM. Due to this result, our initial design of applying the lysed and purified RNA sample directly on the detection paper strip had to be discarded. Since it is known from literature that Cas proteins show activity independent of their activation mechanism at high concentrations, we could not increase the concentration of Cas13a to improve the sensitivity. Instead, we explored amplification methods upstream in the detection process.

LightbringerReal

Figure 2: Estimated detection limit determined for the Cas13a system using 1 nM Cas13a and 10 nM crRNA.

Estimating the Detection Limit of the Reaction Cascade

In order to compare the detection limit of the Cas13a protein kinetics by itself with the detection limit implementing the reaction cascade of reverse transcription, recombinase polymerase amplification and in-vitro transcription, we created a model that takes exponential amplification of the target RNA into account. This was constructed in a way that amplification of the RNA signal takes in an exponential manner, illustrated in Figure. In order to assure that the amplification of the target RNA reaches an upper limit, we capped the amplification at 1000 nM. As visible in Figure XXX, the detection limit of the reaction cascade decreases the detection limit by approximately three orders of magnitude.

RT-RPA-TX

Figure 3: Schematic representation of the target RNA amplification during the estimation of the detection limit using the reaction cascade.

RT-RPA-TX_scheme

Figure 4: Scheme for the RT-RPA-Tx Amplification system.

LightbringerReal

Figure 3: Kinetics of the Cas13a systemusing 1 nM Cas13a and 10 nM crRNA at different target concentrations using the reaction cascade.

Lysis on Chip

We modelled the lysis process on chip to get an idea of how long lysis would need to take place in order to release enough RNA for downstream amplification. For this, we constructed a very simplistic model for bacterial cell lysis. In this, we estimated the rate constants for cell lysis by common colony PCR protocols which use a 10 minute lysis step at 95 °C for thermolysis. Thus, we considered a half-time of Bacteria of 2 minutes at 95 °C. This would result in a lysis efficiency of 96.875%. Starting from this estimation, we considered the rate constant of lysis and thus the half-time using Arrhenius equation as commonly done in the literature:

(1)
(2)

with rate constants k1 and k2 at temperature T1 and T2 and Boltzmann constant R.

where R is the gas constant and k1 and k2 are the rate constant at temperature T1 and T2 The activation energy difference E_A was fitted to a barrier that follows the common rule of thumb that lysis should increase twice in efficiency for every temperature increase of 10 °C. The model for lysis is shown derived in the following:

The full model can then be described by the coupled ordinary differential equations:

(3)
(4)

with klysis being the rate constant of bacterial lysis, kRNase the rate constant of RNA degradation, count of target RNA [targetRNA] and count of bacteria Baks(t). The solution to equation 3 is of course simply:

(5)

Plugging equation 5 into equation 4 gives

(6)

where ratio determines the copy number of a target RNA in a single cell. This differential equation has the form and thus the analytical solution:

Equation 7 + 8

(7)
(8)

with initial condition of

(10)

we get the final solution to the lysis equation:

(11)

The full model at different temperatures looks as follows:

Lysis_Temperature

Figure 3: Effect of lysis temperature on the lysis efficiency of bacterial cells and determination of the released concentration of target RNA from lysis assuming a ratio of 30 RNA molecules per cell.

References

  1. Moody, C., et al. (2016). “A mathematical model of recombinase polymerase amplification under continuously stirred conditions.” Biochemical Engineering Journal 112: 193-201.
  2. Moody, C., et al. (2016). “A mathematical model of recombinase polymerase amplification under continuously stirred conditions.” Biochemical Engineering Journal 112(Supplement C): 193-201.
  3. Li, L., et al. (2011). “Kinetics of hydrothermal inactivation of endotoxins.” Appl Environ Microbiol 77(8): 2640-2647.
  4. Valente, W., et al. (2009). “A Kinetic Study of In Vitro Lysis of Mycobacterium smegmatis.” Chem Eng Sci 64(9): 1944-1952.
  5. Fabritz, H. (2007). “Autoclaves Qualification & Validation.” Experts Meeting in Baden
  6. Marras, S. A., et al. (2004). “Real-time measurement of in vitro transcription.” Nucleic Acids Res 32(9): e72.
  7. Mafart, P., et al. (2002). “On calculating sterility in thermal preservation methods: application of the Weibull frequency distribution model.” Int J Food Microbiol 72(1-2): 107-113.
  8. Chiruta, J. (2000). “Thermat Sterilisation Kinetics of Bacteria as Influenced by Combined Temperature and pH in Continuous Processing of Liquid.” Thesis, The Universit of Adelaide Department of Chemical Engineering Faculty of Engineering.
  9. Rauhut, R. and G. Klug (1999). “mRNA degradation in bacteria.” FEMS Microbiol Rev 23(3): 353-370.
  10. Licciardello, J. J., & Nickerson, J. T. R. (1963). “Some Observations on Bacterial Thermal Death Time Curves.” Applied Microbiology, 11(6), 476–480.
  11. Deindoerfer, F. H. (1957). “Calculation of Heat Sterilization Times for Fermentation Media.” Applied Microbiology, 5(4), 221–228.