Difference between revisions of "Team:Waterloo/Modeling/CaMV Replication"

m
(Added sensitivity analysis explanation)
Line 343: Line 343:
  
 
         <h3>Sensitivity Analysis: Effect of parameter values on steady-state</h3>
 
         <h3>Sensitivity Analysis: Effect of parameter values on steady-state</h3>
         <p>Exploration of parameter values important, especially when we were unable to find characterization of many of the processes integral to CaMV replication</p>
+
         <p>Exploration of parameter values is crucial to determining how robust our system is, especially for parameters which we were unable to find values for. Below is a table showing the relative sensitivity of our system to each parameter at the value listed in the parameter table above. Note we omitted $f_u$, $k_p$, $k_a$, $\delta_3$, $\delta_4$, and $\delta_5$ since for each of those parameters the relative sensitivity was less than 0.01.</p>
 +
 
 +
        <p>We note that virion concentration at 16 hours is increased the most by increasing the transcription rate of 35S RNA. This is intuitive to us since this leads to more 35S RNA, which means more pregenomes to be packaged and more P3, P4, and P5 available to produce virions. We also see a significant increase from increasing the maximum number of genomes in the nucleus, transcription rates for P3, P4, P5, and P6, and the transcription rate of 19S RNA.</p>
 +
 
 +
        <p>Virion concentration was also very sensitive to $k_g$ and $k_c$, our rates of P6 deactivation by CRISPR/Cas9. An increase in $k_g$ and $k_c$ effectiveness leads to a relatively larger decrease in virion concentrations. This is exciting since our application is focused around the implementation of the CRISPR/Cas9 system. Other parameters that will cause a large relative decrease in virion concentration are the maximum RNAi effectiveness ($L$) and the half-saturation constant for P6 activation.</p>
  
 
                 <h4>Relative Sensitivity at 16 hours</h4>
 
                 <h4>Relative Sensitivity at 16 hours</h4>

Revision as of 23:56, 18 September 2015

CaMV Replication

Intracellular replication of CaMV may be strongly affected by the addition of CRISPR/Cas9 and in turn may influence viral spread. Therefore, it's crucial for us to accurately model the replication process to predict viral production over time. In this way, our model of replication will connect our work on CRISPR/Cas9 Dynamics with our agent-based modelling of Intercellular Viral Spread. The first goal of the model is to demonstrate the reduction of viral load over time due to CRISPR/Cas9 modification of ORFVI. This reduction in viral load will hopefully allow for plant defenses to effectively contain the virus before the onset of systemic infection.

Stylized genome
Plant Cells

This leads to our second goal which is to integrate the intracellular replication and intercellular spread models to fully understand impact of CRISPR/Cas9. Altogether, we are attempting to demonstrate the feasibility of our anti-viral system and use our findings to direct the project design.

Model Formation

There are several core processes involved in CaMV replication within an infected host cell. The details of CaMV replication are discussed on the CaMV Biology page and shown in the figure. In brief, however, we determined that in addition to typical processes such as transcription, constitutive translation, and degradation we must also account for the following processes in our model:

  • Virions reinfecting the nucleus
  • Repair of pdsDNA
  • Splicing of 35S RNA
  • RNA interference
  • Packaging of 35S RNA
  • P6 transactivating production of P1-P5
  • Anchoring of P3 to the viral capsid
  • Virions exiting the cell

Additionally, since our project uses CRISPR to produce a targeted mutation in the P6 gene, we track both wild-type and mutated versions of genomes, pregenomes, and virions.

CaMV Replication
i

Network and Equations

The replication pathway of CaMV, simplified to contain only the processes of interest, is shown in the figure below. Gapped CaMV genomes ($g_{gap}$) are closed ($g_{cc}$) and modified by CRISPR/Cas9 ($g_{m.gap}$, $g_{m.cc}$). The modified genomes produce 35S transcripts $35S_m$, but do not produce P6, which is needed to activate translation of the $35S$ transcript into $P3$, $P4$ and $P5$. Proteins $P1$ and $P2$ are disregarded because they mostly play a role in aphid transmission, which we did not model.

CaMV Replication Network
CaMV replication pathway, showing main steps of CRISPR/Cas9 modification, transcription, translation and viral packaging.

The proteins produced from either 35S transcript ($35S$ or $35S_m$) may be used to package virions, though virions produced from $35S_m$, while they are packaged and transported to other cells, will not be able to initiate infection due to their deactivated P6 gene. It should be noted that transport in and out of the nucleus is not modelled explicitly, though certain species only reside in the nucleus and others only in the cytoplasm.

The principles of mass-action kinetics, reviewed by Ingalls , were used to extract a series of ordinary differential equations from this network. The construction of the network and extraction of equations required a number of assumptions, which are detailed below.

CaMV DNA

CaMV DNA is present in both gapped and ungapped forms in the nucleus and may be converted to "modified" form by CRISPR/Cas9, which removes its P6 activity. There is assumed to be a carrying capacity for viruses within the nucleus because of limited cellular resources and the nucleus may be reinfected by virions packaged in the cytoplasm.

$$\frac{d d_g}{dt} = k_v V (d_{max} - d_{total}) - \alpha_c d_g - k_g d_g - \gamma_d d_g$$ $$\frac{d d_{cc}}{dt} = \alpha_c d_g - k_c d_{cc} - \gamma_d d_{cc}$$ $$\frac{d d_{m.g}}{dt} = k_v V_m (d_{max} - d_{total}) - \alpha_c d_{m.g} + k_g d_g - \gamma_d d_{m.g}$$ $$\frac{d d_{m.cc}}{dt} = \alpha_c d_{m.g} + k_c d_{cc} - \gamma_d d_{m.cc}$$

RNA

We consider three transcripts, all assumed to be of importance only in the cytoplasm. The 19S mRNA transcript, the 35S transcript and the modified 35S transcript. These transcripts are degraded both passively, according to the $\gamma_{19S}$ and $\gamma_{35S}$ parameters, and actively by the plant RNAi defenses, captured in the $\gamma_{r}$ term. The 35S transcripts are also packaged into virions by a P4/P5 complex present on the inclusion bodies.

$$\frac{d r_{19S}}{dt} = \alpha_{19S} d_{cc} - (\gamma_{19S}+\gamma_{r}) r_{19S}$$ $$\frac{d r_{35S}}{dt} = \alpha_{35S} d_{cc} - k_p p_4 p_5 f_u r_{35S} - (\gamma_{35S}+\gamma_{r}) r_{35S}$$ $$\frac{d r_{35Sm}}{dt} = \alpha_{35S} d_{m.cc} - k_p p_4 p_5 f_u r_{35Sm} - (\gamma_{35S}+\gamma_{r}) r_{35Sm}$$

Protein

The concentrations of four proteins are modelled. Translation of P3, P4 and P5 from the either of the 35S transcripts occurs because of activation by P6, while P6 is translated from the 19S transcript. All transcripts proteins are degraded and P3, P4 and P5 participate in viral packaging.

$$\frac{d p_3}{dt} = \beta_3 \left( \frac{p_6}{p_6+K_6} \right) (r_{35S} + r_{35Sm}) - k_a p_3 (V_i+V_{im}) - \delta_3 p_3$$ $$\frac{d p_4}{dt} = \beta_4 \left( \frac{p_6}{p_6+K_6} \right) (r_{35S} + r_{35Sm}) - k_p p_4 p_5 f_u (r_{35S} + r_{35Sm}) - \delta_4 p_4$$ $$\frac{d p_5}{dt} = \beta_5 \left( \frac{p_6}{p_6+K_6} \right) (r_{35S} + r_{35Sm}) - k_p p_4 p_5 f_u (r_{35S} + r_{35Sm}) - \delta_5 p_5$$ $$\frac{d p_6}{dt} = \beta_6 r_{19S} - \delta_6 p_6$$

Virions

The concentrations of four proteins are modelled. Translation of P3, P4 and P5 from the either of the 35S transcripts occurs because of activation by P6, while P6 is translated from the 19S transcript. All transcripts proteins are degraded and P3, P4 and P5 participate in viral packaging. Due to the complicated nature and relative lack of characterization of viral packaging, we decided in our model to combine the packaging reaction into one term, dependent on P4, P5, and unspliced 35S RNA concentration.

$$\frac{d V_i}{dt} = k_p p_4 p_5 f_u r_{35S} - k_a p_3 V_i$$ $$\frac{d V}{dt} = k_a p_3 V_i - k_v V (d_{max} - d_{total}) - v_e V - \delta_v V$$ $$\frac{d V_{im}}{dt} = k_p p_4 p_5 f_u r_{35Sm} - k_a p_3 V_{im}$$ $$\frac{d V_m}{dt} = k_a p_3 V_{im} - k_v V_m (d_{max} - d_{total}) - v_e V_m - \delta_v V_m$$

Algebraic Equations

Not all equations pertinent to the dynamics are differential: in particular, the carrying capacity of genomes in the nucleus and effect of RNAi are modelled using algebraic equations.

$$d_{total} = d_g + d_{cc} + d_{m.g} + d_{m.cc}$$ $$\gamma_r = \frac{L}{1+e^{k (p_6-x_0)}}$$ $$x_0 = \frac{1}{2} p_6 ^{ss}$$ $$p_6 ^{ss} = \frac{\beta_6}{\delta_6} \frac{\alpha_{19}}{\gamma_{19}} d_{max}$$

  • Cell concentrations are continuous
  • Molecules in the cytosplasm are well-mixed
  • No outside infection
    • This ODE model only tracks replication within one cell, it cannot track multiple cells. This is handled by the viral spread model instead
  • Limited number of genomes in nucleus
    • "Another pool of viral genomes, in the order of 10-100 copies of minichromosomes comprising supercoiled circular viral DNA and host histones, accumulates in the nucleus."
  • Rate of repair of gapped DNA follows law of mass-action
  • Rate of P6 gene mutation is proportional to number of wild-type genomes
  • DNA, RNA, proteins, and complete virions degrade, incomplete virions do not
    • We assume P3 binds to virions rapidly enough for incomplete virion degradation to be negligible
  • We are only targeting 19S RNA
  • RNA production follows mass-action
  • Concentration of spliced/unspliced RNA is at rapid equilibrium
  • Spliced and unspliced 35S RNA degrade at the same rate
  • Encapsidation and reverse transcription occur simultaneously
    • Information for CaMV is lacking, but for HBV (another pararetrovirus) reverse transcription is initiated during encapsidation and encapsidation is initiated by RT . The complexity of this packaging process is not fully captured in this model
  • Only unspliced RNA is packaged
    • Spliced RNA lacks P1 and P2 and so, although able to replicate within the cell, virions which package spliced RNA will be unable to propagate
  • P1 and P2 do not affect replication dynamics
    • P1 is primarily involved in cell-to-cell movement and P2 with host-to-host movement. Neither have a significant impact on the replication process
  • P6 is only translated from 19S RNA, translation rate is proportional to cccDNA copy number
  • P3, P4, and P5 are translated from 35S RNA, translation is activated by P6
    • This is very well established
  • All P6 is incorporated in inclusion bodies
  • Mass-action anchoring of P3 to virions
  • All P4 is instantaneously spliced
    • As discussed above the subtleties of encapsidation are not captured by this model
  • Virions may reinfect nucleus
    • This is very well established
  • Virions leave the cell at a constant rate

Model Parameters

Finding and implementing accurate parameter values is crucial when trying to draw conclusions from a model. Our model contains nearly thirty parameters, each of which needed to be found to exactly reproduce the viral replication process in silico. Model dynamics may be drastically affected by our choice in parameters and this must be accounted for in our analysis. However, as with many models in systems biology our network suffers from the "parameter problem" described in chapter 1 of the book by Gunawardena which arises from the complexity of biological processes and the difficulty in measuring parameter values.

A number of parameter values are the same as those used in an HBV model developed by Nakabayashi & Sasaki . Parameters that have a reference in the "Source" column are ones we were able to find or derive accurate values for from the literature. Unfortunately, there are a few parameter values that could not be found (e.g. the half-saturation constant of P6 activation) and so we explicitly mention the absence of data in the table. In order to compensate for the lack of accurate parameters, we analyze our model over a wide range of parameters to determine how behaviour will be affected.

Lastly, for simplicity we normalized each concentration by the volume of the cell (volume=1.414x10$^{-13}$m$^3$).

Symbol Value Units Description Source
$k_v$ 0.1 min$^{-1}$ Rate at which virions produced by the cell reinfect the nucleus. No source found, parameter is further analyzed in the results section
$d_{max}$ 100 molecules/volume Maximum concentration of viral genomes in the nucleus. "Another pool of viral genomes, in the order of 10-100 copies of minichromosomes comprising supercoiled circular viral DNA and host histones, accumulates in the nucleus"
$\alpha_c$ 0.1 min$^{-1}$ Rate at which gaps are repaired in gapped DNA to form cccDNA. HBV model by Nakabayashi & Sasaki
$k_g$ 0.01 min$^{-1}$ Rate at which the P6 gene on the gapped DNA is modified. No source found, parameter is further analyzed in the results section
$k_c$ 0.01 min$^{-1}$ Rate at which the P6 gene on cccDNA is modified. No source found, parameter is further analyzed in the results section
$\gamma_d$ 0.001 min$^{-1}$ DNA Degradation rate. HBV model by Nakabayashi & Sasaki
$\alpha_{19S}$ 0.01 min$^{-1}$ Transcription rate of 19S RNA. HBV model by Nakabayashi & Sasaki
$\alpha_{35S}$ 0.09 min$^{-1}$ Transcription rate of 35S RNA. HBV model by Nakabayashi & Sasaki
$\gamma_{19S}$, $\gamma_{35S}$ 0.001 min$^{-1}$ Degradation rate of 19S RNA, 35S RNA. HBV model by Nakabayashi & Sasaki
$f_u$ 0.3 unitless Fraction of unspliced 35S RNA in the cell (assumed at equilibrium). Froissart et al. 2004 <\cite> and Bouton et al. 2015 <\cite>
$\beta_3$, $\beta_4$, $\beta_5$, $\beta_6$ 0.1 min$^{-1}$ Translation rate of P3, P4, P5, P6. HBV model by Nakabayashi & Sasaki
$K_6$ 5000 molecules/volume Half-saturation constant for transactivation of P1-P5 production. No source found, parameter is further analyzed in the results section
$\delta_3$, $\delta_4$, $\delta_5$, $\delta_6$ 0.001 min$^{-1}$ Degradation rate of P3, P4, P5, P6. HBV model by Nakabayashi & Sasaki
$k_p$ 0.1 molecules$^{-2}$ min$^{-1}$ Packaging rate. No source found, parameter is further analyzed in the results section
$k_a$ 1 molecules$^{-1}$min$^{-1}$ (converted) Rate of P3 anchoring to virions. No source found, parameter is further analyzed in the results section
$v_e$ 0.1 min$^{-1}$ Rate at which virions exit the cell. No source found, parameter is further analyzed in the results section
$\delta_v$ 0.001 min$^{-1}$ Rate of virion degradation. HBV model by Nakabayashi & Sasaki
$L$ 0.005 min$^{-1}$ Maximum RNAi effectiveness. Estimated to be on the order of $\gamma_{35S}$, but slightly larger since RNAi is more effective than degradation
$k$ 10$^{-4}$ molecules$^{-1}$ Slope of logistic function describing RNAi effectiveness No source found, parameter is further analyzed in the results section.

Results

Time Course of CaMV Infection in a Single Cell

The effects of RNAi and CRISPR/Cas9 the time course of an infection by CaMV in a single cell were aexplored.

Time course x 4

Investigating Final Virion Concentration

Heatmaps

Sensitivity Analysis: Effect of parameter values on steady-state

Exploration of parameter values is crucial to determining how robust our system is, especially for parameters which we were unable to find values for. Below is a table showing the relative sensitivity of our system to each parameter at the value listed in the parameter table above. Note we omitted $f_u$, $k_p$, $k_a$, $\delta_3$, $\delta_4$, and $\delta_5$ since for each of those parameters the relative sensitivity was less than 0.01.

We note that virion concentration at 16 hours is increased the most by increasing the transcription rate of 35S RNA. This is intuitive to us since this leads to more 35S RNA, which means more pregenomes to be packaged and more P3, P4, and P5 available to produce virions. We also see a significant increase from increasing the maximum number of genomes in the nucleus, transcription rates for P3, P4, P5, and P6, and the transcription rate of 19S RNA.

Virion concentration was also very sensitive to $k_g$ and $k_c$, our rates of P6 deactivation by CRISPR/Cas9. An increase in $k_g$ and $k_c$ effectiveness leads to a relatively larger decrease in virion concentrations. This is exciting since our application is focused around the implementation of the CRISPR/Cas9 system. Other parameters that will cause a large relative decrease in virion concentration are the maximum RNAi effectiveness ($L$) and the half-saturation constant for P6 activation.

Relative Sensitivity at 16 hours

Parameter Relative Sensitivity
$\alpha_{35}$ $1.48$
$d_{max}$ $1.26$
$\beta_3, \beta_4, \beta_5, \beta_6$ $1.06$
$\beta_6$ $1.03$
$\alpha_{19}$ $1.03$
$\alpha_c$ $0.20$
$k_v$ $0.01$
$...$ $...$
$\delta_v$ $-0.01$
$k$ $-0.05$
$\gamma_d$ $-0.09$
$\gamma_{19}$ $-0.14$
$\gamma_{35}$ $-0.18$
$\delta_6$ $-0.34$
$v_e$ $-0.99$
$K_6$ $-1.03$
$L$ $-1.55$
$k_g, k_c$ $-1.74$

Phase Plane Analysis: Effect of initial conditions on steady state

Discussion

References

    Top