BioMed Research International

Volume 2016 (2016), Article ID 5621313, 8 pages

http://dx.doi.org/10.1155/2016/5621313

## Optimal Control Model of Tumor Treatment with Oncolytic Virus and MEK Inhibitor

School of Mathematics and Physics, University of Science and Technology Beijing, Beijing, China

Received 21 September 2016; Accepted 27 November 2016

Academic Editor: Tun-Wen Pai

Copyright © 2016 Yongmei Su et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Tumors are a serious threat to human health. The oncolytic virus is a kind of tumor killer virus which can infect and lyse cancer cells and spread through the tumor, while leaving normal cells largely unharmed. Mathematical models can help us to understand the tumor-virus dynamics and find better treatment strategies. This paper gives a new mathematical model of tumor therapy with oncolytic virus and MEK inhibitor. Stable analysis was given. Because mitogen-activated protein kinase (MEK) can not only lead to greater oncolytic virus infection into cancer cells, but also limit the replication of the virus, in order to provide the best dosage of MEK inhibitors and balance the positive and negative effect of the inhibitors, we put forward an optimal control problem of the inhibitor. The optimal strategies are given by theory and simulation.

#### 1. Introduction

Tumors are a serious threat to human health, and chemotherapy and radiotherapy may not only kill cancer cells, but also damage human body normal cells at the same time [1]. The oncolytic virus is a kind of tumor killer virus which can infect and lyse cancer cells and spread through the tumor, while leaving normal cells largely unharmed [2]. When oncolytic viruses are inoculated into a cancer patient or directly injected into a tumor, these viruses will spread throughout the tumor and infect tumor cells. The viruses can be replicated in the infected tumor cells. When an infected tumor cell is lysed, it can burst out a mass of new oncolytic viruses. Then, these new viruses can infect much more neighboring tumor cells [3].

Experiments using oncolytic viruses such as adenovirus, CN706 [4], and ONYX-15 [5] in animal tumors show that these viruses are nontoxic and infect tumor cells specifically. Now, treatment of cancer with oncolytic virus has been clinically tested [6–8]. This treatment of cancer with oncolytic viruses has been explored by clinicians [9–11].

In recent years, in order to understand the cancer-virus dynamics and find better treatment strategies, some mathematical models have been set up [12–19]. Tian proposed a mathematical model to describe the development of a growing tumor and an oncolytic virus population as follows [18]:where variables stand for the population of uninfected cells, infected tumor cells, and oncolytic viruses, respectively. The coefficient represents the infection of the virus. The tumor growth is modeled by logistic growth, and is the maximal tumor size. is the per capita tumor growth rate. means the lysis rate of the infected tumor cells. represents the burst size of new viruses coming out from the lysis of an infected tumor cell. represents the death rate of the virus.

It was shown that when the threshold , the equilibrium solution is globally asymptotically stable [18], indicating that the oncolytic virus therapy finally has no effect. Obviously, the smaller the value of , the more easily holds. Since represents the total number of tumor cells, smaller tumors may be more resistant to the treatment by oncolytic virus than large ones, which should be a contradiction. In [19], by replacing with , we proposed the model The meanings of variables and parameters , and are the same as those in model (1), and is positive and sufficiently small. The threshold obtained by our model is , which is almost independent of when is sufficiently small.

On the other hand, all the above papers did not consider coxsackie-adenovirus receptor (CAR). In fact, CAR is a main receptor when oncolytic viruses enter into tumor cells [20–22]. The successful entry of viruses into cancer cells is related to the presence of CAR. When oncolytic viruses infect the tumor cells, firstly, they combine with the CAR and are absorbed into the cells.

Mitogen-activated protein kinase (MEK) inhibitors have been shown to promote CAR expression and could increase oncolytic viruses infection into tumor cells. But MEK inhibitors may also limit the replication of viruses [23–25], which will affect the treatment by oncolytic virus. With the function of MEK, [25] gave a model:The variables , and have the same meanings as those in model (2); represents the average expression level of CAR on the surface of the cells. The intensity of MEK inhibitor application is captured in the parameter . If , there is no MEK inhibitor application, and the CAR expression level will gradually decline. If , the MEK inhibitor has the maximum possible effect. The model assumes that exponential growth can be slowed down by the inhibitor with expression . CAR grow at the rate of and become extinct at the rate of .

Based on models (2) and (3), we establish the following mathematical model:The variables , and have the same meanings as those in model (3). The parameters , and are the same as those of (1). The parameter has the same meaning as that in model (3). All the parameters are strictly positive.

Since the use of MEK inhibitors not only results in enhanced oncolytic virus entry into the tumor cells, but also renders infected cells temporarily unable to produce viruses, the maximum dosage of MEK use may not result in the best treatment effect, so the optimal control-based schedules of MEK inhibitor application should be studied. The optimal MEK inhibitor application strategy can increase the efficacy of this treatment in an economical fashion. So, first, in this paper, we let the control variable be a constant; a stability analysis of our model is conducted, and then the optimal control strategy is discussed; we also compare the optimal control with constant control by simulation.

#### 2. Materials and Methods

##### 2.1. Stability Analysis

System (4) always has two equilibrium points:If is sufficiently small, when , the third steady state exists in whichHere, It should be noted that is equivalent to to ensure that when is sufficiently small.

The Jacobi matrix at point isObviously, is one eigenvalue of which means is unstable. The unstable result of seems consistent with the biological meaning that, without viruses and infected tumor cells, the tumor will grow from an initial small value around .

As for equilibrium point , we have the following theorem.

Theorem 1. *When , is locally asymptotically stable. When , is unstable.*

*Proof. *At the equilibrium point , the Jacobi matrix iswhere . The eigenvalues of are ,Here, in which .

When , we havewhich ensure that and are negative, so is locally asymptotically stable.

Similarly, when , is positive, and is unstable.

Actually, we can prove that the equilibrium solution is globally asymptotically stable when . But we need to show the boundness of system (4). From the first two equations, we obtainBy the comparison principle, we can obtain .

From the third equation of (4), we can have It is easily shown that .

Similarly, fromwe can get that holds.

Theorem 2. *When , is globally asymptotically stable.*

*Proof. *Consider the Lyapunov function ; the derivative along a solution is given by Since , we have , which implies and because , therefore, When , we can have .

Let ; it is clear that . Let be the largest positively invariant subset of the set ; by the third equation of system (4), we can know that , so By LaSalle invariance principal [26], we knowSo, the limit equation of system (4) is Therefore, when . So, is globally attractive; note that can also ensure the local asymptotical stability of , so we can know that of system (4) is globally asymptotically stable when .

Although we can prove the global asymptotical stability of* E*_{1}, we would not want this to happen, because the global asymptotical stability means the therapy does not have any effect. When holds, the coexistent steady state exists, but it is difficult to give the stable analysis of* E*_{2}, so we just give some simulations about it.

We choose , , , , , , , , and .

The initial condition is , where the unit of each is cells.

We choose , respectively, and all hold; the simulation results are shown in Figures 1 and 2.

The simulation results show that oncolytic virus therapy may keep the tumor stable at some level as shown in Figure 1 or keep oscillating at a certain range as shown in Figure 2. Since the cured equilibrium is always unstable, just from our model, we could not give the condition that ensures the tumor can be cured by oncolytic virus therapy, but if we choose appropriate* u* which satisfiesthe simulation shows that oncolytic virus therapy can prevent the tumor from getting worse and worse. Some other therapy methods should be combined to cure the tumor.