Modeling

Introduction

The project is intended to treat retrovirus diseases. And we take AIDS as an example. In this project, we want to establish a more effective HIV-curing system with less side-effects by integrating CRISPR/Cas9 system into human hemopoietic stem cells, which aims to protect the helper T cells from virus infection.

So we want to know how effective the CRISPR/Cas9 system is for HIV. In this model, we discuss about dynamic changes of different cells in a person’s body using Matlab. And in this model, we can see how the system works. And this model will help us to forecast the future applications of our project.

In our model, we can see the influence the CRISPR/Cas system produces by the change of virus and cells. But there are too many parameters and that makes it impossible to get clear relation of these parameters easily. If the environment doesn't change, we can stimulate the situation in Matlab for a short time.

Derivation of the model

Theoretical analysis

First, stably transfect hematopoietic stem cells with the plasmid which encodes Cas9 protein. To make the question easier, we suppose that the efficiency is constant. As is shown in other part in this project, second, we can construct a plasmid which encodes the gRNA sequence which is a 20bp conserved regions in HIV DNA and has no conservation in human genome. When the human body is infected by AIDS, using non-viral DNA transfecting methods (in our project we choose A-B toxin based shuttle) to deliver the plasmid into human somatic cells (especially CD4+ cells) and the CRISPR/Cas system is activated. The Cas9 protein recognizes viral DNA with the help of gRNA and cut the target sequence, causing double strand break (DSB). The DSB is then repaired by the non-homologous end joining mechanism in the cell but with high error rate (in this model we suppose the rate doesn’t change), causing frame-shift mutation and gene knockout. The viral protein synthesis is stopped and further infection is impossible, but it provides potentiation for eradicating AIDS.

Symbol description

\(vg_{tg}\) the self-generation rate of stem cell

\( N_{thymus}\) the possible total thymus counting

\( N_h\) high resistant thymus number

\( N_l\) low resistant thymus number

\( P_{thymus}\) the enviroment pressure for thymus growth in M-model

\( f_{ht}\) the hiv-caused thymus counting decreasing factor

\( f_{hlh}\) the high and low immune cell induced increasing factor

\(vg_{sh}\) the self-generation rate of hiv virus

\(f_m\) the mutation factor

\( N_{hiv}\) the possible total hiv counting, can be set to be infinity

\(P_{hiv}\) Almost no self-pressure for hiv virus

\(vg_{th}\) the generation rate by thymus counting

\(f_{hr}\) the resistance factor, can be invicible if CRISPR is 100% effective

\(vd_{td}\) the decay(death) rate from high resistant T cell to low ones

\(vg_{lh}\) the gain rate from low resistant T cell to high ones

\(vg_{tl}\) generation rate from thymus counting

\(f_{lr}\) the resistance factor

\(vd_{hl}\) the death rate caused by hiv

\(vg_{gt}\) the gain rate by the immune decaying of high resistant T cells

Formula derivation

We suppose that a person was infected by HIV for some time. And the environment doesn't change in our model. Because the existence of CRISPR/CAS9 systems, we can roughly divide the immunological cells into two sides. One is with high resistance, the other is with low resistance.

HIV virus will attract T cells, so we can suppose that there is a linear relationship between the decreasing rate and the number of HIV and we will get eqution (1) as follow. \begin{equation} \frac{{\rm d} N^{(1)}_{thymus}}{{\rm d} t}=-\alpha_1 N_{hiv} \end{equation}

\begin{equation} \frac{{\rm d} N^{(2)}_{thymus}}{{\rm d} t}=\alpha_2 N_h +\alpha_3 N_l \end{equation}

\begin{equation} \frac{{\rm d} N_{thymus}}{{\rm d} t}=-\alpha_1 N_{hiv}+\alpha_2 N_h +\alpha_3 N_l \end{equation}

\begin{equation} \frac{{\rm d} N_{hiv}}{{\rm d} t}=-\beta_1 f_m(\beta_2 N_l+\beta_3 N_h) \end{equation}

\begin{equation} \frac{{\rm d} N_{hiv}}{{\rm d} t}=\epsilon(1-f_m)N_{hiv}N_h \end{equation}

\begin{equation} \frac{{\rm d} N_{h}}{{\rm d} t}=\delta N_{thymus} \end{equation}

\begin{equation} \frac{{\rm d} N_{h}}{{\rm d} t}=\gamma_1(1-\gamma_2 N_{hiv}) \end{equation}

\begin{equation} \frac{{\rm d} N_{h}}{{\rm d} t}=\delta N_{thymus}+\gamma_1(1-\gamma_2 N_{hiv}) \end{equation}

\begin{equation} \frac{{\rm d} vd_{td}}{{\rm d} t}=(\eta N_h N_{hiv}e^{\varepsilon f_{hr}}) \end{equation}

\begin{equation} \frac{{\rm d} N_h}{{\rm d} t}=\zeta N_l e^{ f_{lr}} \end{equation}

\begin{equation} \frac{{\rm d} N_l}{{\rm d} t}=\iota N_{thymus} \end{equation}

\begin{equation} \frac{{\rm d} N_h}{{\rm d} t}=1-\kappa N_{hiv} \end{equation}

\begin{equation} \frac{{\rm d} N_l}{{\rm d} t}=\lambda N_l N_{hiv} e^{\mu f_{lr}} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

\begin{equation} \end{equation}

Matlab analysis

this is a \(\left( \sum_{k=1}^n a_k b_k \right)^2 \leq \left( \sum_{k=1}^n a_k^2 \right) \left( \sum_{k=1}^n b_k^2 \right)\) .

However , \begin{equation} \left( \sum_{k=1}^n a_k b_k \right)^2 \leq \left( \sum_{k=1}^n a_k^2 \right) \left( \sum_{k=1}^n b_k^2 \right) \end{equation}

Reference

Berezovskaya, F. S., Wolf, Y. I., Koonin, E. V., & Karev, G. P. (2014). Pseudo-chaotic oscillations in CRISPR-virus coevolution predicted by bifurcation analysis. Biology direct, 9(1), 13.