^{1}

^{1}

^{2}

^{1}

^{*}

A mathematical modeling of tumor therapy with oncolytic viruses is discussed. The model consists of two coupled, deterministic differential equations allowing for cell reproduction and death, and cell infection. The model is one of the conceptual mathematical models of tumor growth that treat a tumor as a dynamic society of interacting cells. In this paper, we obtain an approximate analytical expression of uninfected and infected cell population by solving the non-linear equations using Homotopy analysis method (HAM). Furthermore, the results are compared with the numerical simulation of the problem using Matlab program. The obtained results are valid for the whole solution domain.

Oncolytic viruses are viruses that infect and kill cancer cells but not normal cells [1-4]. Oncolytic virus therapy originated early in the last century upon the observation of occasional tumor regressions in cancer patients suffering from virus infections or those receiving vaccinations. Many types of oncolytic viruses have been studied as candidate therapeutic agents including adenoviruses, herpes viruses, reoviruses, paramyxoviruses, retroviruses, and others [2,4]. Probably, the best-characterized oncolytic virus, that has drawn a lot of attention, is ONYX- 015, an attenuated adenovirus that selectively infects tumor cells with a defect in the p53 gene [

The oncolytic effect has several possible mechanisms that yield complex results. The first such mechanism is the result of viral replication within the cell and rupture out of the cell [11,12]. The third mechanism involves virus infection of cancer cells that induces antitumoral immunity. Surviving mice acquired a resistance to rechallenge with tumor cells [

Several mathematical models that describe the evolution of tumors under viral injection were recently developed. Wodarz [13,16] presented a mathematical model that describes interaction between two types of tumor cells (the cells that are infected by the virus and the cells that are not infected by the virus) and the immune system. However, to the best of our knowledge, till date no general analytical expressions for the mathematical modeling of two populations of cells namely uninfected tumor cells and infected cells [

The model, which considers two types of cells growing in the logistic fashion, has the following form [

The equation can be solved subject to the following initial conditions:

where X is the size of the uninfected cell population; Y is the size of the infected cell population; and are the maximum per capita growth rates of uninfected and infected cells; K is the carrying capacity; b is the transmission coefficient; and a is the rate of infected cell killing by virus. We introduce the following set of dimensionless variables as follows [

The governing non-linear differential equations (Equations (1) and (2)) expressed in the following non-dimensionless format:

An appropriate set of boundary condition are given by:

, (8)

The Homotopy analysis method (HAM) is a powerful and easy-to-use analytic tool for nonlinear problems [21- 23]. It contains the auxiliary parameter h, which provides us with a simple way to adjust and control the convergence region of solution series. Furthermore, the obtained result is of high accuracy. Solving the Equations (6) and (7) using HAM (see Appendix A) and simultaneous equation method (see Appendix B and C), the steady state and transient contributions to the model are given by:

where

Similarly we can obtain as follows:

where

Here and represent a time independent steady state term and and denote the time dependent transient component. The steady state term will be important at long times as In contrast the transient term will be of important at short times as

The non-linear equations (Equations (6) and (7)) for the boundary conditions (Equation (8)) are solved by numerically. The function ode45 in Scilab software is used to solve two-point boundary value problems (BVPs) for ordinary differential equations. The Matlab program is also given in Appendix C. The numerical results are also compared with the obtained analytical expressions (Equations (9) and (12)) for all values of parameters, , , and.

Equations (9) and (12) represent the simple approximate analytical expressions of the uninfected and infected tumor cells for all values of parameters, , , and. The two types of tumor cell growing using Equations (9) and (12) are represented in Figures 1-4. In Figures 1 and 2, the dimensionless uninfected tumor cells reach the constant value when for some fixed value of and different values of and. The dimensionless infected tumor cells reaches the

steady state value when.

In Figures 3 and 4, the dimensionless infected tumor cells reach the constant value when for some fixed value of and different values of and. The dimensionless infected tumor cells reaches the steady state value when. Figures 5 and 6 give us the confirmation for the above discussion in three-dimensional graphs also.

In this work, the coupled system of time dependent differential equations for the two types of cells growing has been solved analytically using the Homotopy Analysis Method. Approximate analytical expressions for uninfected and infected cell population are derived for all values of parameters. Furthermore, on the basis of the outcome of this work, it is possible to calculate the approximate rate of the tumor cells growth. The extension of the procedure to other systems that include interaction between tumor cells and anticancer agents seem possible.

This work was supported by the Council of Scientific and Industrial Research (CSIR No.: 01(2442)/10/EMR-II), Government of India. The authors also thank Mr. M.S. Meenakshisundaram, Secretary, The Madura College Board, Dr. R. Murali, The Principal, and Prof. S. Thigarajan, HOD, Department of Mathematics, The Madura College, Madurai, Tamilnadu, India for their constant encouragement. The authors S. Usha is very thankful to the Manonmaniam Sundaranar University, Tirunelveli for allowing the research work.

Consider the following nonlinear differential equation

where N is a nonlinear operator, t denotes an independent variable, is an unknown function. For simplicity, we ignore all boundary or initial conditions, which can be treated in the similar way. By means of generalizing the conventional Homotopy method, Liao constructed the so-called zero-order deformation equation as:

where is the embedding parameter, is a nonzero auxiliary parameter, is an auxiliary function, L is an auxiliary linear operator, is an initial guess of, is an unknown function. It is important, that one has great freedom to choose auxiliary unknowns in HAM. Obviously, when and, it holds:

and (A3)

respectively. Thus, as p increases from 0 to 1, the solution varies from the initial guess to the solution. Expanding in Taylor series with respect to p, we have:

where

If the auxiliary linear operator, the initial guess, the auxiliary parameter h, and the auxiliary function are so properly chosen, the series (A4) converges at p = 1 then we have:

Define the vector

Differentiating Equation (A2) for m times with respect to the embedding parameter p, and then setting and finally dividing them by m!, we will have the socalled mth-order deformation equation as:

where

and

Applying on both side of Equation (A8), we get

In this way, it is easily to obtain for at order, we have

When, we get an accurate approximation of the original Equation (A1). For the convergence of the above method we refer the reader to Liao. If Equation (A1) admits unique solution, then this method will produce the unique solution. If Equation (A1) does not possess unique solution, the HAM will give a solution among many other (possible) solutions.

For the case of steady-state condition, the Equations (6) and (7) becomes as follows:

Solving the above Equations (B1) and (B2), we get

and

Thus we can obtain and as in the text (Equations (10) and (13)).

The given differential equations for the non-steady state condition are of the form as:

For the transient part, the initial conditions are redefined as

In order to solve the Equations (C1) and (C2) by means of the HAM, we first construct the Zeroth-order deformation equation by taking,

We have the solution series as

and

where

Substituting Equations (C7) and (C8) into Equations (C5) and (C6) and comparing the coefficient of like powers of p, we get

and so on.

The initial conditions are redefined as

and

for (C16)

Solving the Equations (C10) and (C10) by using the boundary conditions given in Equations (C14) and (C15), we get

and

Substituting the Equations (C17) and (C18) in Equations (C12) and (C12) and by using the boundary conditions given in Equation (C16), we get

and

Adding the Equations (C17) and (C19) and the Equations (C18) and (C20), we get the Equations (9) and (12) as in the text.

function main1 options= odeset('RelTol',1e-6,'Stats','on');

Xo = [10; 2];

tspan = [0,10];

tic

[t,X] = ode45(@TestFunction,tspan,Xo,options);

toc figure hold on plot(t, X(:,1),'blue')

plot(t, X(:,2),'green')

return function [dx_dt]= TestFunction(t,x)

a=16;

b=3;

r=10;

dx_dt(1) =x(1)*(1-(x(1)+x(2)))-(b*x(1)*x(2))/(x(1)+x(2));

dx_dt(2)=r*x(2)*(1-(x(1)+x(2)))+(b*x(1)*x(2))/(x(1)+x(2))-a*x(2);

dx_dt = dx_dt';

Symbol Meaning

Size of the uninfected cell population

Size of the infected cell population

Rate of infected cell killing by the virus

Transmission coefficient

Maximum per capita growth rates of uninfected cells

Maximum per capita growth rates of infected cells

Carrying capacity

Time

Size of the dimensionless uninfected cell population

Size of the dimensionless infected cell population

Dimensionless time