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

\(M_{thymus}\) the maximum number of immunological cells

\(N_{thymus}\) the number of immunological cells

\( N_h\) high resistant immune cells number

\( N_l\) low resistant immune cells number

\( P_{thymus}\) the enviroment pressure for immune cell growth

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

\( f_{ht}\) the hiv-caused the number of immunological cells decreasing factor

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

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

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

\(f_m\) the mutation factor

\(f_{hr}\) the high resistance factor

\(f_{lr}\) the low resistance factor

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

\(N_{hiv}\) the number of hiv virus

\(vg_{th}\) the generation rate by the number of immunological cells

\(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 the number of immunological cells

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

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

\(rand\) random number between 0 and 1

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 system and other environment effects, we can roughly divide the immunological cells into two sides. One is with high resistance, the other is with low resistance.

English letters express variables, while Greek letters express parameters.

HIV infects vital cells in the human immune system such as helper T cells (specifically CD4+ T cells), macrophages, and so on. 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 equation (1) as follow.

\begin{equation} \frac{{\rm d} N^{(1)}_{thymus}}{{\rm d} t}=-\alpha_1 N_{hiv} \end{equation} Some kinds of immunological cells can replicate, so the increasing rate of immunological cells is influenced by the number of immunological cells. So we will get equation (2) as follow. \begin{equation} \frac{{\rm d} N^{(2)}_{thymus}}{{\rm d} t}=\alpha_2 N_h +\alpha_3 N_l \end{equation} T cells as hosts can boost the increasing of HIV, so we can get equation (3). \begin{equation} \frac{{\rm d} N_{hiv}}{{\rm d} t}=\beta_1 f_m(\beta_2 N_l+\beta_3 N_h) \end{equation} High resistance ones can decrease the number of HIV virus, which implies equation (4). \begin{equation} \frac{{\rm d} N_{hiv}}{{\rm d} t}=\epsilon(1-f_m)N_{hiv}N_h \end{equation} Of course, high resistance ones increase due to T cells' increasing. \begin{equation} \frac{{\rm d} N^{(1)}_{h}}{{\rm d} t}=\delta N_{thymus} \end{equation} HIV virus caused the high ones increasing slower than normal situation. \begin{equation} \frac{{\rm d} N^{(2)}_{h}}{{\rm d} t}=\gamma_1(1-\gamma_2 N_{hiv}) \end{equation} If we combine equation (5) and equation (6), we can get equation (7) as follows. \begin{equation} \frac{{\rm d} N_{h}}{{\rm d} t}=\delta N_{thymus}+\gamma_1(1-\gamma_2 N_{hiv}) \end{equation} Immunity decay rate for high ones due to mutation, which we can express by equation (9). \begin{equation} \frac{{\rm d} vd_{td}}{{\rm d} t}=\eta N_h N_{hiv}e^{\varepsilon f_{hr}} \end{equation} The number of high ones increase due to low ones that enhance the immune capacity. \begin{equation} \frac{{\rm d} N_h}{{\rm d} t}=\zeta N_l e^{ f_{lr}} \end{equation} Because low ones increase due to the number of immunological cells, we can get mathematic relation between the two variables. \begin{equation} \frac{{\rm d} N_l}{{\rm d} t}=\iota N_{thymus} \end{equation} HIV virus caused the high ones increasing slower than normal. \begin{equation} \frac{{\rm d} N_h}{{\rm d} t}=1-\kappa N_{hiv} \end{equation} Cells will die after HIV's attack. death formula for low ones due to HIV virus is as follow. \begin{equation} \frac{{\rm d} N_l}{{\rm d} t}=\lambda N_l N_{hiv} e^{\mu f_{lr}} \end{equation}

The values of parameters will depend on whether there is CRISPR/CAS9 system in the body.

\(N_{thymus}\) changes with the change of time \(t \) \begin{equation} N_{thymus}=\int_{0}^{t}[vg_{tg}M_{thymus}(1-\frac{P_{thymus}N_{thymus}}{M_{thymus}})+f_{ht}\alpha_1 N_{hiv}+f_{hlh}\alpha_2 N_h +\alpha_3 N_l]{{\rm d} t} \end{equation} \(N_{hiv}\) changes with the change of time \(t \) \begin{equation} N_{hiv}=\int_{0}^{t}[vg_{sh}N_{hiv}(1-\frac {N_{hiv}}{M_{hiv}}P_{hiv})\beta_1 f_m(\beta_2 N_l+\beta_3 N_h))-\epsilon(1-f_m)N_{hiv}N_h]{{\rm d} t} \end{equation} \(N_h\) changes with the change of time \(t \) \begin{equation} N_h=\int_{0}^{t}[vg_{th}\delta N_{thymus}\gamma_1(1-\gamma_2 N_{hiv})-vd_{td}\eta N_h N_{hiv}e^{\varepsilon f_hr}+vg_{lh}\zeta N_l e^{ f_{lr}} {{\rm d} t} \end{equation} \(N_l\) changes with the change of time \(t \) \begin{equation} N_l=\int_{0}^{t}[vg_{tl}\iota N_{thymus}(1-\kappa N_{hiv})-vd_{hl}\lambda N_l N_{hiv} e^{\mu f_{lr}}]{{\rm d} t} \end{equation}

\begin{equation} f_m=\xi_1(\frac{N_h}{N_h+N_l}+\xi_2rand) \end{equation}

\begin{equation} \end{equation}

Matlab analysis

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. Cunningham, A. L., Donaghy, H., Harman, A. N., Kim, M., & Turville, S. G. (2010). Manipulation of dendritic cell function by viruses. Current opinion in microbiology, 13(4), 524-529.