Journal of Tumor Research

Journal of Tumor Research
Open Access

ISSN: 2684-1258

+44 1223 790975

Research Article - (2016) Volume 2, Issue 1

Dynamics of Tumor-Immune System with Fractional-Order

Fathalla A Rihan1*, Adel Hashish2, Fatma Al-Maskari3, Mohamud Sheek-Hussein3, Elsayed Ahmed4, Muhammad B Riaz5 and Radoune Yafia6
1Department of Mathematical Sciences, College of Science, UAE University, 15551, Al-Ain, UAE
2Department of Physics, College of Science, UAE University, 15551, Al-Ain, UAE
3Zayed bin Sultan Al Nahyan Center for Health Sciences, College of Medicine and Health Sciences, UAE University, Al-Ain, UAE
4Department of Mathematics, Faculty of Science, Mansoura University, Mansoura, 35516, Egypt
5Department of Mathematics, University of Management and Technology, Lahore, Pakistan
6Ibn Zohr University, 638, Ouarzazate, Morocco
*Corresponding Author: Fathalla A Rihan, Department of Mathematical Sciences, College of Science, UAE University, 15551, Al-Ain, UAE, Tel: +971-3-7673333 Email:

Abstract

Most of biological systems have long-range temporal memory. Modeling of such systems by fractional-order (or arbitrary-order) models provides the systems with long-time memory and gains them extra degrees of freedom. Herein, we suggest a simple fractional-order model to describe the dynamics of tumor-immune interactions. Two effector cells are considered, in the model, with a Holling function response of type-III. The model is extended to include treatment terms which represent an external source of the effectors cells by ACI and an external input of IL- 2. Asymptotic stabilities of tumour-free steady state and persistent- tumour steady state are studied. The threshold parameter R0 (average number of newly infected cells produced by a single councerous cell) is deduced. The numerical simulations show that the fractional-order derivative enriches the dynamics of the system and increases the complexity of the observed behaviours, which confirms that the fractional-order may play the role of memory in the system.

Keywords: Cancer; Fractional-order; Numerical simulation; Stability; Steady states; Tumor-immune system

Introduction

Tumors are a family of high-mortality diseases, exhibiting a derangement of cellular proliferation which often lead to uncontrolled cell growth [1,2]. Research efforts are being devoted to understand the interaction between the tumour cells and the immune system [3-6]. Mathematical models, using ordinary differential equations with integer-order, have been proven valuable in understanding the dynamics of tumour-immune system and how host immune cells and cancerous cells evolve and interact; See e.g. [7-12]. However, modeling of biological systems by fractional-order differential equations has more advantages than classical integer-order mathematical modeling, in which such affects the memory are neglected. Accordingly, the subject of fractional calculus (that is, calculus of integral and derivatives of arbitrary order) has gained popularity and importance, mainly due to its demonstrated applications in system biology [13,14] and other fields of sciences [15-18]. The Fractional-Order Differential Equations (FODEs) models are more consistent with the biological phenomena than those of integer-orders. This is due to the fact that fractionalorder derivatives the description of the memory and hereditary properties inherent in the processes [19]. It has been deduced in [13] that the membranes of the organism have fractional-order electrical conductance, which are classified under the non-integer order models.

Fractional Models of Tumor-Immune System

Immune system is considered as one of the most fascinating schemes in terms of biology and mathematics. Immune system is multi-functional with several metabolic pathways; therefore most effector cells perform more than one function. Plus, each function of the immune system is typically done by more than one effector, which makes it more complex [20]. Differential equations, with integerorders, have long been used in modeling cancer phenomena [21-23], but fractional- order differential equations have short history in modeling such systems with memory [24]. Herein, we use FODEs in modeling tumor-immune interactions, which are naturally related to systems with memory.

Assume that model of cancer-immune system includes two immune effectors: E1(t) and E2(t) (such as cytotoxic T-cells, and natural killer cells), interacting with the cancer cells, T(t) with tumour’s functional response of Holling Type-III [25]. The model takes the form

image

The parameter d1 represents natural decay rate of the effect cells E1(t). k1 and k2 are half saturation parameters, a is the growth rate of tumour cells. r1 and r2 are the reduction rate of tumour cells due to presence of the effector cells E1(t) and E2(t) respectively. d2 represents natural decay rate of the effect cells E2(t). All of these parameters are supposed to be positive constants. The interaction terms in the first and third equations of model (1) satisfy the cross reactivity property of the immune system. It has been assumed that (d1k1/(1−d1))<<(d2k2/(1− d2)), to avoid the non-biological interior solution where both immune effectors coexist. For more details about model (1), we refer to [16].

Here, we have modified model (1) to include external sources of effector cells and immune stimulation effects by treatment Interleukin-2 (IL-2). Assume that three populations of the activated immune-system cells, E(t); the tumor cells, T(t); and the concentration of IL-2 in the single tumor-site compartment, IL(t) (Figure 1). To ease the analysis, consider a classic bilinear model that includes Holling Type-I function and external effector cells s1 and input of IL-2, s2. The interactions of the three populations are then governed by the fractional-order differential model:

tumor-diagnostics-reports-tumour-immune

Figure 1: Key-players in tumour-immune interactions described in model (2).

image

with initial conditions: image . The first equation describes the rate of change in the effector cells population. The parameter p1 represents the antigenicity rate of the tumor (immune response to the appearance of the tumor), p3 is the cooperation rate of effector cells with Interleukin-2 parameter, and s1 represents the external source of the effector cells, with rate of death p2. The second equation shows the rate of change of the tumour cells which follows a logistic growth (a type of limiting growth) in the absence of immune response. The parameter p4 incorporates growth rate of tumor cells. The maximal carrying capacity of the biological environment for tumor cell is p−15 . Whereas, p6 is the rate of tumor cells. The third equation gives the rate of change for the concentration of IL-2. Its source is the effector cells, which are stimulated by interaction with the tumour. The parameter p7 is the competition rate between the effector cells and the tumor cells. The external input of IL-2 into the system is s2 and the loss-rate parameter of IL-2 cells is p8.

In the absence of immunotherapy with IL-2, we have

image

To minimize sensitivity (or robustness) of the model, we nondimensionalize the bilinear system (3) by taking the following rescaling:

image

substitution into (3), we have

image

(Here t is replaced by t*.) The analytical stability region of fractionalorder system (4) is given in Figure 2.

We then study the stability of the steady states of models (1) and (4).

Equilibria and local stability of model (1)

The equilibrium points of system (1) are:

image

Here ε0 is the naive first equilibrium, ε1 is the memory equilibrium and the ε2 is endemic according to the value of the tumor size. Stability analysis shows that the naive state is unstable.

However, the memory state is locally asymptotically stable if

image

While the endemic state is locally asymptotically stable if 0 R > 1 and 0 < d2 < 1. There is bifurcation at R0=1. The stability of the memory state depends on the value of one parameter namely the immune effector death rate.

Equilibria and local stability of model (4)

The steady states of the reduced model (4) are again the intersection of the null-clines image If y=0, the tumorfree equilibrium is at image This steady state always exist, since image From the analysis, it is easy to prove that the tumorfree equilibrium image of the model (4) is asymptotically stable if threshold parameter (the minimum tumor-clearance parameter) image and unstable ifimage

However, if y≠0, the steady states are obtained by solving

image

In this case, we have two endemic equilibria, image

image

with image The Jacobian matrix of system (4) at the endemic equilibrium ε1 is

image

tumor-diagnostics-reports-fractional-order

Figure 2: Stability region of the fractional-order system (4) is enlarged when 0<α ≤ 1, where λ is the root of the characteristic equation.

Proposition 1: Assume that the endemic equilibrium ε1 exists and has nonnegative coordinates.

If image thenimage is unstable.

Proof Since

image

image

then inequality image is true if

image

Therefore, when image we haveimage and hence both sides of the inequality are positive. Therefore if the equilibrium point image exists and has nonnegative coordinates, then imageand the point (image) is unstable wheneverimage

Similarly, we arrive at the following Proportion.

Proposition 2: If the point image exists and has nonnegative coordinates, then it is asymptotically stable.

Proposition 3: The presence of a fractional differential order in a differential equation can lead to a notable increase in the complexity of the observed behaviour, and the solution is continuously depends on all the previous states (Figure 1).

Results and Discussion

The numerical technique discussed in [16] was used to numerically simulate the qualitative behaviours of the fractional-order models. The numerical technique is based on Caputo sense for fractional derivative (Appendix A) and implicit Euler’s approximation, for the resulting systems, with step-size h=0.05 and 0.5<α ≤ 1 and parameters values given in the captions of the figures.

Figure 2 provides indication that the fractional differential order enlarges the stability region of the solution when 0<α ≤ 1. Figure 3 displays the numerical simulation of (1). The left banner shows the memory equilibrium point ε1 ; While right banner shows an endemic equilibrium ε2 with sustained oscillations. Figures 4 and 5 confirm that the fractional derivative damps the oscillation behaviour for model (1).

In this paper, we presented two fractional-order models for tumorimmune interactions. In the first model, we provide a Holling type- III function and cross reactivity in fractional-order model for tumorimmune interactions. Two immune effectors have been considered because of the property of multi-functional and multi-pathways of the immune system. We obtained memory states, using the fractionalorder, whose stability depends on the value of one parameter namely the immune effector death rate. However, in the second model, we extend the model to include external treatments then reduce the system into a prey-predator model. The models possess non-negative solutions, as desired in any population dynamics. We deduced the threshold parameters Rand image respectively for each model. These parameters represent the minimum tumor-clearance parameter or minimum infection free. It has been seen from the numerical simulations that the fractional-order enriches the dynamics of the system and enlarges the stability regions of the solutions. Although the fractional-order model (4) is simple, it displays up to three steady states. Figure 6 shows different types of steady states: Tumour-free steady state image, dormancy-, medium and high persistent-tumour steady states (right)image for the model (4).

tumor-diagnostics-reports-numerical-simulations

Figure 3: Shows the numerical simulations of model (1) when α=1.0 and a=r1=r2=1; d1=0.3, d2=0.7, k1=0.3, k2=0.7 (left banners) where the system converges to steady state ε1 ; and when d1=0.7, d2=0.3 (right banner), where the system converges to steady state ε2 . We note sustained oscillations.

tumor-diagnostics-reports-system-converges

Figure 4: Shows the numerical simulations of model (1) when α=0.95 and a=r1=r2=1; d1=0.3, d2=0.7, k1=0.3, k2=0.7 (left banners) where the system converges to steady state ε1 ; and when d1=0.7, d2=0.3 (right banner), where the system converges to steady state ε2 . The fractional derivative damps the oscillation behaviour.

tumor-diagnostics-reports-system-converges

Figure 5: Shows the numerical simulations of model (1) when α=0.75 and a=r1=r2=1; d1=0.3, d2=0.7, k1=0.3, k2=0.7 (left banners) where the system converges to steady state ε1 ; and when d1=0.7, d2=0.3 (right banner), where the system converges to steady state ε2 . The fractional derivative damps the oscillation behaviour.

As a result, in the endemic steady states, Figure 7 shows that the fractional-order derivative kills the oscillation behaviour. Figure 8 displays the numerical simulations for the model (4) with different values of the fractional order and parameter values given in the caption. The tumour-free image is locally asymptotically stable asimage .From the graphs, it can be seen that FODEs have rich dynamics and are better descriptors of biological systems than traditional integerorder models. The equilibria for infection-free and endemic fractionalorder cases are the same as integer-order counterparts.

Conclusion

In this paper, the authors conclude that the fractional-order derivatives in the models provide an excellent instrument for the description of memory and hereditary properties of inter and intra cells. It is possible that the tumour may result in either equilibrium with (dormancy) or escape from the immune system. The fractionalorder differential equations are, at least, stable as their integer-order counterparts. The presence of a fractional differential order in the differential equation can lead to a notable increase in the complexity of the observed behaviour, as the solution continuously depends on all the previous states of the solutions. The analysis can be extended to include more components of immune response and control variables. It can also be extended to describe the dynamics of hepatitis B and C virus infections.

Appendix A: Preliminaries

The subject of fractional calculus deals with the investigations of derivatives and integrals of any arbitrary real or complex order, which unify and extend the notions of integer-order derivative and n-fold integral. Here, we provide some definitions of fractional-order integration and fractional-order differentiation [26]. There are several approaches to the generalization of the notion of differentiation to fractional-orders e.g. Riemann-Liouville, Caputo and Generalized Functions approach. Let image be the class of Lebesgue integrable functions on [a, b], image.

tumor-diagnostics-reports-manifold

Figure 6: Shows (left) tumour-free steady state image and low-, medium and high persistent-tumour steady states (right) image,for the model (4). The stable steady states are [•], unstable steady states are [◦], stable manifold is [−−−] and initial conditions are [+], with external treatment doses σ>0 and particular parameters.

tumor-diagnostics-reports-fractional-order

Figure 7: Shows numerical simulations of model (4) when s=0.1181, ω=0.1184, θ=0.1747, r=0.636, b=0.002. _α12=1, 0.9, 0.7. Note that image is locally asymptoticallyimage. The fractional-order derivative kills the oscillation behaviour.

tumor-diagnostics-reports-asymptotically-stable

Figure 8: Shows the numerical simulations of model (4) when σ=0.1211, ω=0.1184, θ=0.3746, a=1.635, b=0.002.α12=1, 0.9, 0.7.image is locally asymptotically stable as image

Definition 1 The fractional-integral (or the RiemannLiouville integral) of order image of the functionimage is defined by

image

The fractional derivative of order aε(n−1,n) of f(t) is defined by two nonequivalent approaches.

(i) Riemann-Liouville fractional derivative in which we take fractional integral of order (n− α) and then take nth derivative,

image

(ii) Caputo fractional derivative in which we take nth derivative and then take a fractional integral of order (n− α)

image

From the definition, we notice that the time-fractional derivative of a function f(t) at t=tn involves an integration and calculating timefractional derivative that requires all the past history, i.e. all the values of f(t) from t=0 to t=tn. In this paper, we have adopted Caputo’s definition which has the advantage of dealing properly with initial value problems. The following Remark addresses some of the main properties of the fractional derivatives and integrals [12,17].

Remark 1

image

image

image

(iv) If f(t) is absolutely continuous on image

image

image

image

If (vi) holds, and image then f(t) is nondecreasing for each image then f(t) is non-increasing for eachimage

Acknowledgement

The work is supported by Center-Based ZCHS-12-2014 and UPAR Projects, UAE University. The authors are grateful to the reviewers and Editor Shonisani C Tshidino for their valuable suggestions and constructive comments which improved the manuscript.

References

  1. Bellomo N, Bellouquid A, Delitala M (2004) Mathematical topics on the modeling of multicellular systems in competition between tumor and immune cells. Math Models Methods Appl Sci 14: 1683-1733.
  2. Preziosi L (2003) Cancer modeling and simulation. Chapman & Hal, CRC Press.
  3. Chaplain MA Modelling aspects of cancer growth: Insight from mathematical and numerical analysis and computational simulation. In: Multiscale Problems in the Life Sciences 1940: 147-200.
  4. Kirschner D, Tsygvintsev A (2009) On the global dynamics of a model for tumor immunotherapy. J Math Biosci Eng 6 3 573–583.
  5. Nagy JD (2005) The ecology and evolutionary biology of cancer: a review of mathematical models of necrosis and tumor cells diversity. Math Biosci Eng 2: 381-418
  6. Roose T, Chapman SJ, Maini PK (2007) Mathematical models of avascular tumor growth. SIAM Rev 49: 179–208.
  7. Bellomo N, Li NK, Maini PK (2008) On the foundations of cancer modeling: selected topics, speculations and perspectives. Math Mod Methods Appl Sci 18: 593-646.
  8. Bellomo N, Preziosi L (2000) Modelling and mathematical problems related to tumor evolution and its interactions with the immune system. Math Comput Model 32: 413-452.
  9. Byrne HM, Alarcon T, Owen MR, Webb SD, Maini PK, et al. (2006) Modeling aspects of cancer dynamics: A review. Philos Trans R Soc A 364: 1563–1578.
  10. Castiglione F, Piccoli B (2007) Cancer immunotherapy, mathematical modeling and optimal control. J TheorBiol 247: 723-732.
  11. Rihan FA, Abdelrahman DH, Al-Maskari F, Ibrahim F, Abdeen MA (2014) Delay differential model for tumour-immune response with chemo–immunotherapy and optimal control. Comput Math Methods Med 2014: 982978-982993.
  12. Rihan FA, Rahman DA, Lakshmanan S, Alkhajeh AS (2014) A time delay model of tumour-immune system interactions: Global dynamics, parameter estimation, sensitivity analysis. Appl Math Comput 232: 606-623.
  13. Cole KS (1993) Electric conductance of biological systems. Cold Spring HarbSymp Quant Biol 1: 107-116.
  14. Rihan FA, Lakshmanan S, Hashish AH, Rakkiyappan R, Ahmed E (2015) Fractional order delayed predator-prey systems with Holling type-II functional response. Nonlinear Dynam 1: 777-789.
  15. Ferdri, Y (2012) Some applications of fractional order calculus to design digital filters for biomedical signal processing. J Mech Med Biol 12: 1-13.
  16. Hilfer R (2000) Applications of fractional calculus in physics. World Scientific, River Edge, NJ, USA.
  17. Yuste SB, Acedo L, Lindenberg K (2004) Sub diffusion-limited A+B → C reaction- sub diffusion process. Physical Review E 69: 036-126.
  18. Zaslavsky GM (2002) Chaos, fractional kinetics, and anomalous transport. Physics Reports 371: 461-580.
  19. Machado JT (2010) Entropy analysis of integer and fractional dynamical systems. Nonlinear Dyn. 62: 371–378.
  20. Ahmed E, Hashish A, Rihan FA (2012) On fractional order cancer model. J Fractional CalcAppl 3: 1–6.
  21. Bellomo N, Bellouquid A, Nieto J, Soler J (2010) Multiscale biological tissue models and fluxlimited chemotaxis from binary mixtures of multicellular growing systems. Math. Models Methods Appl Sci 20: 1179-1207.
  22. Gökdoǧan A, Yildirim A, Merdan M (2011) Solving a fractional order model of HIV infection of CD4+ T cells. math. Comput Modeling 54: 2132-2138.
  23. Kirschner D, Panetta J (1998) Modeling immunotherapy of the tumor immune interaction. J Math Biol 37: 235-252.
  24. Rihan FA (2013) Numerical modeling of fractional-order biological systems. In Abstract and Applied Analysis 2013: 1-11.
  25. Apreutesei N, Dimitriu G (2010) On a prey-predator reaction-diffusion system with Holling type III functional response. J ComputApplMath 235: 366-379.
Citation: Rihan FA, Hashish A, Al-Maskari F, Hussein MS, Ahmed E, et al. (2016) Dynamics of Tumor-Immune System with Fractional-Order. J Tumor Res 2: 109.

Copyright: ©2016 Rihan FA, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Top