Difference between revisions of "Team:Waterloo/Modeling"
(way more plant defense results) |
|||
| Line 44: | Line 44: | ||
<section id="antiviral" title="Plant Defense"> | <section id="antiviral" title="Plant Defense"> | ||
<h2>CRISPR Plant Defense</h2> | <h2>CRISPR Plant Defense</h2> | ||
| − | + | <p> | |
| + | Mathematical models are not limited by the same time constraints as growing plants, so they were used to validate the CRISPR Plant Defense aspect of our project. The team looked at the antiviral effects of CRISPR/Cas9 targeting on three scales: CaMV genomes, plant cells and plant leaves. A background primer on Cauliflower Mosaic Virus (CaMV) genetics, replication and spread can be found on the <a href="https://2015.igem.org/Team:Waterloo/Modeling/CaMV_Biology">CaMV Biology page</a>.</p> | ||
<div class="row"> | <div class="row"> | ||
<div class="col-sm-4"> | <div class="col-sm-4"> | ||
| Line 67: | Line 68: | ||
<p>On the <strong>scale of individual genomes</strong>, we used our <a href="https://2015.igem.org/Team:Waterloo/Modeling/Cas9_Dynamics">dynamics model of CRISPR/Cas9</a> to model viral genomes becoming non-functional over time as frameshift mutations were introduced. The results of running 1000 simulations using our <a href="">sgRNA target design</a> to deactivate the CaMV P6 gene predict the following pattern of P6 functionality over time. | <p>On the <strong>scale of individual genomes</strong>, we used our <a href="https://2015.igem.org/Team:Waterloo/Modeling/Cas9_Dynamics">dynamics model of CRISPR/Cas9</a> to model viral genomes becoming non-functional over time as frameshift mutations were introduced. The results of running 1000 simulations using our <a href="">sgRNA target design</a> to deactivate the CaMV P6 gene predict the following pattern of P6 functionality over time. | ||
| + | |||
<figure> | <figure> | ||
| − | <img src="/wiki/images/ | + | <img src="/wiki/images/a/a4/Waterloo_p6_fit.png" alt="Exponential decay of functional viral genomes to non-functional genomes" /> |
| − | <figcaption> | + | <figcaption>CRISPR effect parameter (t1/2) is derived from an exponential fit to the fraction of CaMV genomes active at each timestep across 1000 stochastic genome simulations.</figcaption> |
</figure> | </figure> | ||
| + | |||
<p>The exponential decay fit shown in the graph above was then passed to a model on the scale of <strong>plant cells</strong>. We modelled <a href="https://2015.igem.org/Team:Waterloo/Modeling/CaMV_Replication">CaMV Replication</a> in a single infected cell using ordinary differential equations. The ODE replication model allowed us to see the effect that the CRISPR/Cas9 decay of P6 would have on <strong>virion production</strong>. We also examined the interaction bewteen intracellular plant defenses (RNA interference) and the CRISPR/Cas9 defense and conducted a sensitivity analysis on our parameters.</p> | <p>The exponential decay fit shown in the graph above was then passed to a model on the scale of <strong>plant cells</strong>. We modelled <a href="https://2015.igem.org/Team:Waterloo/Modeling/CaMV_Replication">CaMV Replication</a> in a single infected cell using ordinary differential equations. The ODE replication model allowed us to see the effect that the CRISPR/Cas9 decay of P6 would have on <strong>virion production</strong>. We also examined the interaction bewteen intracellular plant defenses (RNA interference) and the CRISPR/Cas9 defense and conducted a sensitivity analysis on our parameters.</p> | ||
<figure> | <figure> | ||
| − | <img src="/wiki/images/ | + | <img src="/wiki/images/3/38/Waterloo_CRISPRCRISPieRtimeseries.png" alt="Time series of functional and non-functional virions in the cell in the presence of CRISPR and CRISPieR systems."/> |
| − | <figcaption> | + | <figcaption>Replication model shows total CaMV virions within a cell over time with and without a CRISPR/CRISPieR defense system. Non-functional virion production is lowered because fewer functional proteins that protect against RNAi are produced.</figcaption> |
</figure> | </figure> | ||
| + | |||
<p>The ODE model was then scaled up to the level of <strong>plant leaves</strong>. We modeled the spread of infection between plant cells using an agent-based framework, which is described in detail on the <a href="https://2015.igem.org/Team:Waterloo/Modeling/Intercellular_Spread">Intercellular Viral Spread</a> page. The agent-based model allowed us to explore the effect of intercellular defense signalling (which leads plant cells to resist to infection and become suseptible to apoptosis).</p> | <p>The ODE model was then scaled up to the level of <strong>plant leaves</strong>. We modeled the spread of infection between plant cells using an agent-based framework, which is described in detail on the <a href="https://2015.igem.org/Team:Waterloo/Modeling/Intercellular_Spread">Intercellular Viral Spread</a> page. The agent-based model allowed us to explore the effect of intercellular defense signalling (which leads plant cells to resist to infection and become suseptible to apoptosis).</p> | ||
| − | |||
<p>Overall, our multi-scale model of CRISPR Plant Defense allowed us to investigate the extent to which we could slow CaMV spread using CRISPR/Cas9 and its interaction with host defense systems such as RNA interference and intercellular defense signalling.</p> | <p>Overall, our multi-scale model of CRISPR Plant Defense allowed us to investigate the extent to which we could slow CaMV spread using CRISPR/Cas9 and its interaction with host defense systems such as RNA interference and intercellular defense signalling.</p> | ||
| + | |||
| + | </p> | ||
| + | |||
| + | </section> | ||
</div> | </div> | ||
</html> | </html> | ||
{{Waterloo_Footer}} | {{Waterloo_Footer}} | ||
Latest revision as of 01:33, 21 November 2015
Modeling
Mathematical modeling is a core part of Waterloo iGEM: we have nearly as many team members typing furiously away in our dry lab as we do wrangling transformations in our wet lab. This year, we created models in Python, MATLAB and NetLogo that examine the feasibility of our design and provide tools for assessing future designs. The code for each model is available on our GitHub and details on the formulation of each model may be found in the pages linked below.
Cas9 Frameshift Dynamics
PAM Structural Bioinformatics
A key limitation of the current use of S. pyogenes Cas9 for gene editing is the inability to target sequences which are not adjacent to a NGG PAM sequence. To overcome this issue, we developed a pipeline to probe Cas9 mutants with altered PAM specificity using simulation tools in PyRosetta and Python. A detailed description of the analysis and software used in this work can be found on the Modeling Engineered Cas9 PAM Flexibility.
The simulation managed to successfully predict the PAM specificity of wild type spCAs9. However, validation against known Cas9 mutants was less successful, predicting a mix of correct and incorrect, though biochemically similar, PAM sequences.
CRISPR Plant Defense
Mathematical models are not limited by the same time constraints as growing plants, so they were used to validate the CRISPR Plant Defense aspect of our project. The team looked at the antiviral effects of CRISPR/Cas9 targeting on three scales: CaMV genomes, plant cells and plant leaves. A background primer on Cauliflower Mosaic Virus (CaMV) genetics, replication and spread can be found on the CaMV Biology page.
On the scale of individual genomes, we used our dynamics model of CRISPR/Cas9 to model viral genomes becoming non-functional over time as frameshift mutations were introduced. The results of running 1000 simulations using our sgRNA target design to deactivate the CaMV P6 gene predict the following pattern of P6 functionality over time.
The exponential decay fit shown in the graph above was then passed to a model on the scale of plant cells. We modelled CaMV Replication in a single infected cell using ordinary differential equations. The ODE replication model allowed us to see the effect that the CRISPR/Cas9 decay of P6 would have on virion production. We also examined the interaction bewteen intracellular plant defenses (RNA interference) and the CRISPR/Cas9 defense and conducted a sensitivity analysis on our parameters.
The ODE model was then scaled up to the level of plant leaves. We modeled the spread of infection between plant cells using an agent-based framework, which is described in detail on the Intercellular Viral Spread page. The agent-based model allowed us to explore the effect of intercellular defense signalling (which leads plant cells to resist to infection and become suseptible to apoptosis).
Overall, our multi-scale model of CRISPR Plant Defense allowed us to investigate the extent to which we could slow CaMV spread using CRISPR/Cas9 and its interaction with host defense systems such as RNA interference and intercellular defense signalling.