� � “main” — 2014/9/1 — 11:35 — page 211 — #1 � � � � � � Tema Tendências em Matemática Aplicada e Computacional, 15, N. 2 (2014), 211-221 © 2014 Sociedade Brasileira de Matemática Aplicada e Computacional www.scielo.br/tema doi: 10.5540/tema.2014.015.02.0211 A Prelude to the Fractional Calculus Applied to Tumor Dynamic* N. VARALTA1, A.V. GOMES1 and R.F. CAMARGO2** Received on December 9, 2013 / Accepted on June 14, 2014 ABSTRACT. In order to refine the solution given by the classical logistic equation and extend its range of applications in the study of tumor dynamics, we propose and solve a generalization of this equation, using the so-called Fractional Calculus, i.e., we replace the ordinary derivative of order 1, in one version of the usual equation, by a non-integer derivative of order 0 < α ≤ 1, and recover the classical solution as a particular case. Finally, we analyze the applicability of this model to describe the growth of cancer tumors. Keywords: biomathematics, fractional calculus, logistics equation, dynamics of cancer tumor. 1 INTRODUCTION The art of getting a differential equation whose solution describes well the reality, brings a great difficulty, generally the closer we are to perfectly describe a real problem the harder and complex the equations involved are, in Albert Einstein words [2]“One thing I have learned in a long life: all our science, measured against reality, is primitive and childlike – and yet it is the most precious thing we have”. In this sense, the so-called Fractional Calculus, which is the branch of mathematics that deals with the study of integrals and derivatives of non-integer orders, plays an outstanding role. There are countless problems that, when described in terms of differential equations of non integer order, provide a more accurate description of reality. Specially phenomena that posses time de- pendence, since the fractional derivatives excellently describe effects of memory and hereditary properties3[2, 19]. Among the different ways to solve a fractional differential equation (FDE), or still an initial value problem, we mention the so-called method of integral transforms, i.e., we transform the *Paper presented at CMAC-SE. **Corresponding author: Rubens de Figueiredo Camargo. 1Department of Biostatistics, IBB, UNESP – Universidade Estadual Paulista, 18618-970 Botucatu, SP, Brazil. E-mails: najla@ibb.unesp.br; ariannevellasco@gmail.com 2Department of Mathematics, FC, UNESP – Universidade Estadual Paulista, 17033-360 Bauru, SP, Brazil. E-mail: rubens@fc.unesp.br 3Indeed, fractional derivative is not a local operator, since it is defined in terms of the fractional integral, as it will be seen in equation (2.10). � � “main” — 2014/9/1 — 11:35 — page 212 — #2 � � � � � � 212 A PRELUDE TO THE FRACTIONAL CALCULUS APPLIED TO TUMOR DYNAMIC initial problem into another one (called transformed problem and normally easier to be solved), solve this new problem and recover the original solution applying the inverse integral transform [2, 3, 19]. Usually the solution of a FDE is expressed in terms of a parameter corresponding to the order of the derivative, and the solution of the corresponding integer order equation is recovered for a particular value of this parameter. In many cases, it is not a whole order of the derivative that makes the solution of equation closer to reality [2, 4, 5]. Moreover, the logistic equation was published in 1838 by Pierre François Verhulst to model the growth of world population and was based on the assessment of population statistics available, complementing the theory of exponential growth of Thomas Robert Malthus. It can be applied to models with time dependence and has a wide application area since the inhibiting factors are taken into consideration. Furthermore, the logistic equation proved applicable in a series of probabilistic events and related to chaos theory and industrial and business dynamics [7, 25]. Recently, the logistic equation has been applied to describe the growth of populations both in the laboratory natural habitat, limiting the growth by influence of competition, mortality and fertility factors. However, the logistic model does not fit very well in cases where there are more complex relationships acting as interactions within food webs or several features that occur in most cases in nature dependencies[5, 6, 7, 8, 9, 22]. Aiming to generalize the logistic equation and give a better description for some of those events, El-Sayed, at all [5], studied the corresponding fractional-order logistic by a numerical analyze. The purpose of this paper is to present the analytic solution of the fractional logistic equation, in the inverse form, and analyze the applicability of this fractional equation to improve the de- scription of the dynamic of cancer tumor and compare this model with some classical models presented in the literature [9, 18], and is organized as follows. In Section 1 we make a brief introduction to our work. In Section 2 we present some functions inherent to fractional calculus and the definition of the fractional integral of Riemann-Liouville and Caputo’s fractional deriva- tive. In Section 3 we present the Classical Logistic Equation and its corresponding solution and also its generalization via Fractional Calculus. Besides that, we recover the classical solution as a particular case. In Section 4 we present some classical models presented in literature that de- scribe the growth of tumors and compare them with the solution of fractional logistic equation. Finally in Section 5 we present our concluding remarks. 2 FRACTIONAL CALCULUS In this section we present the definition of some important special functions related to fractional calculus and their basic properties, that will allow us to propose and solve the fractional version of the classical logistic equation. 2.1 Gamma Function We denote by �(z) the Gamma function, this function can be defined by the improper integral [19]: �(z) = ∫ ∞ 0 e−t t z−1dt, (2.1) Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 213 — #3 � � � � � � VARALTA, VELLASCO and CAMARGO 213 with Re(z) > 0, the integral is convergent. Solving this integral by parts and a change of variable, it follows that: �(z + 1) = z�(z). Once that �(1) = 1, it is easy to see that, for n ∈ N, �(n + 1) = n!, i.e., the Gamma function is a generalization of factorial4. 2.2 Gel’fand-Shilov function definition Let n and ν be, respectively, a natural and non-integer number with Re(ν) > 0, we define the Gel’fand-Shilov function as φn(t) = ⎧⎨ ⎩ t n−1 (n − 1)! se t ≥ 0 0 se t < 0. and φν(t) = ⎧⎨ ⎩ tν−1 �(ν) se t ≥ 0 0 se t < 0. (2.2) Thus, the Laplace transform of φν(t) is given by: �[φν(t)] = ∫ ∞ 0 e−st tν−1 �(ν) dt = 1 �(ν) ∫ ∞ 0 e−a (a s )ν−1 da s = s−ν, (2.3) in which, the second equality is due to the change of variable (st = a) and the last equality is due to the definition of gamma function. 2.3 Mittag-Leffler functions The classical Mittag-Leffler function is a complex function depending on a complex parameter and was defined and studied by Mittag-Lefller [14, 15, 16], in the year 1903, as Eα(z) = ∞∑ k=0 zk �(αk + 1) , Re(α) > 0. (2.4) This function is a direct generalization of the exponential function since for α = 1 we have the exponential function. The function defined by Eα,β (z) = ∞∑ k=0 zk �(αk + β) Re(α), Re(β) > 0, (2.5) gives a generalization of Eα(z) since Eα,1(z) = Eα(z). This generalization was studied by Wiman [26] in 1905. The Laplace transform of the function tβ−1Eα,β(atα), is given by: � [ tβ−1Eα,β(±atα) ] = sα−β sα ∓ a . (2.6) 4As we will see, the generalization of factorial, by the gamma function, is the key for the generalization of the usual calculus by the fractional calculus. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 214 — #4 � � � � � � 214 A PRELUDE TO THE FRACTIONAL CALCULUS APPLIED TO TUMOR DYNAMIC One important relation between the Mittag-Lffler functions with one and two parameters is given by Eα,α+1(z) = ∞∑ n=0 (z)n �(αn + α + 1) = 1 z ∞∑ n=0 (z)n+1 �[α(n + 1) + 1] = 1 z [Eα(z) − 1]. (2.7) 2.4 Fractional integral Here we first define the integral operator of order n and then, using the generalization of factorial by Gamma function, we define the fractional integral of Riemann-Liouville. Definition. Let n ∈ N and f (t) : R→ R an integrable function. We define the Integral Operator, I , of order 1 and n respectively as: I f (t) = ∫ t 0 f (t1) dt1 and I n f (t) = ∫ t 0 ∫ t1 0 · · · ∫ tn−1 0 f (tn) dtn dtn−1 . . . dt1. Theorem. For f (t) : R→ R integrable, we have that [2]: I n f (t) = φn(t) ∗ f (t) = ∫ t 0 φn(t − τ) f (τ ) dτ = ∫ t 0 (t − τ)n−1 (n − 1)! f (τ ) dτ, (2.8) where ∗ denotes the Laplace convolution product and φn(t) is the Gel’fand-Shilov function. Definition. Let f (t) be an integrable function and ν ∈ C, such that Re(ν) > 0, the Riemann- Liouville fractional integral of order ν of f (t), denoted by I ν f (t), is defined as5: I ν f (t) = φν(t) ∗ f (t) = ∫ t 0 (t − τ)ν−1 �(ν) f (τ ) dτ (2.9) 2.5 Caputo’s Fractional Derivative Caputo’s fractional derivative is based on the fact that derivative is the left inverse operator of integral and on Lagrange’s exponential law. According to Caputo, the fractional derivative is the fractional integral of an integer order derivative such that exponential law makes sense [10]. Thus, let f (t) be a differentiable function, m ∈ N and α ∈ N such that m − 1 < Re(α) < m, we have: Dα f (t) = I m−α Dm f (t) = φm−α(t) ∗ Dm f (t). (2.10) 5Note that this definition comes from equation (2.8) and from the fact that Gel’fand-Shillov function is well defined for non integer numbers. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 215 — #5 � � � � � � VARALTA, VELLASCO and CAMARGO 215 2.6 Laplace Transform From the definition (2.10) and from the Laplace convolution theorem we may write: �[Dα f (t)] = �[�m−α(t)] �[Dm f (t)] = sα−m �[Dm f (t)]. (2.11) 3 FRACTIONAL LOGISTIC EQUATION In this section, we present the classic logistic equation and its solution and also their respective generalized using fractional calculus. 3.1 The Classic Logistic Equation In 1838, Verhulst [25] published logistic equation defined as: d dt N(t) = kN(t) ( 1 − N(t) r ) . (3.1) where N(t) is the number of individuals in time t , k is the intrinsic growth rate and r is carrying capacity. Without loss of generality, taking r = 1, this equation can be rewritten in the form: d dt N(t) = kN(t)[1 − N(t)]. (3.2) Note that this is a non-linear and separable equation. Though solving this equation is quite simple we cannot say the same about its fractional version. Using the variable change v(t) = 1/N(t) in the previous equation we obtain the following Linear and separable equation: dv(t) dt = k[1 − v(t)], (3.3) whose solution is given by v(t) = 1 + 1 c e−kt ⇒ v(0) = 1 + 1 c · Since N(t) = v(t)−1, we can write 1/c = 1/N(0) − 1, as a result; N = 1 e−kt c−1 + 1 ⇒ N = 1 1 + [ 1 N(0) − 1 ] e−kt . (3.4) Note that 0 < N(0) < 1 and limt→∞ N(t) = 1 (i.e. the carrying capacity). 3.2 Fractional Logistic Equation (Inverse form) We will use the results presented above to propose a generalization, via fractional calculus, for the logistic equation in the linear form given by equation (3.3). It is very important to note that the fractional version with order 0 < α ≤ 1 of equation (3.1) was numerically solved by El-Sayed, at all [5], however, we will present the analytic solution of equation (3.3). Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 216 — #6 � � � � � � 216 A PRELUDE TO THE FRACTIONAL CALCULUS APPLIED TO TUMOR DYNAMIC The fractional version of equation (3.3) with order 0 < α ≤ 1 is given by dαv(t) dtα = Dαv(t) = k[1 − v(t)]. (3.5) Applying Laplace transform to both sides, following: � [ Dαv(t) ] = k� [1 − v(t)] . Using equation (2.11), with m = 1 (since 0 < α ≤ 1) we have: sα F(s) − sα−1v(0) = k [ 1 s − F(s) ] Then F(s) = k [ s−1 sα + k ] + v(0) [ sα−1 sα + k ] , in which F(s) = �[v(t)]. Thus we have: v(t) = �−1� [v(t)] = k�−1 [ s−1 sα + k ] + v(0)�−1 [ sα−1 sα + k ] . From this, the result of equation (2.6), we have: v(t) = tαkEα,α+1(−ktα) + v(0)Eα(−ktα). (3.6) Then, using equation (2.7) with z = −ktα v(t) = 1 + Eα(−ktα)[v(0) − 1]. (3.7) Once we v(t) = 1/N(t), we obtain: N(t) = 1 1 + [ 1 N(0) − 1 ] Eα(−ktα) · (3.8) Note that, limα→1 N(t) = 1 1+ [ 1 N(0) −1 ] e−kt , i.e., the solution of integer order is a particular case of the fractional solution. Below is the graph corresponding to the solution of equation (3.5), taking N(0) = 0.2 carrying capacity r = 1, for different values of α, we have: Since limt→∞ Eα(−ktα) = 0 for all values 0 < α ≤ 1 we have: lim t→∞ N(t) = lim t→∞ 1 1 + [ 1 N(0) − 1 ] Eα(−ktα) = 1, that is, the all the values considered converge to the value of the carrying capacity. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 217 — #7 � � � � � � VARALTA, VELLASCO and CAMARGO 217 0 �0.2 �0.4 �0.6 �0.8 �1 �1 2 4 6 8 10 0.2 0.4 0.6 0.8 1 t N(t) α α α α α k Figure 1: The solution of equation (3.5). 4 TUMOR DYNAMIC With the artifice of scientific research and the use of mathematics in modeling, description and even the prediction of a given physical process we can better understand some factors in the biological context. However, such use is recent, particularly in the study of cancer tumors one of the reasons is that we still finding difficulty in modeling satisfactorily the tumor’s growth. The models found are currently based on computer simulations, differential equations and curve fitting, but even the most sophisticated models are unable to predict how the size of the host interferes in tumor’s growth. There are several models of tumor growth in the literature, thus, it is very important to know which one best fits with experimental data. It is of paramount importance to understand the assumptions and the consequences of such models because often these models supports more complex models of tumor growth [5, 6, 7, 8, 9, 22]. It can be considered a volume of tissue as a “cellular community” comprises “species” of mes- enchymal and epithelial cells in dynamic equilibrium with the environment and with each other [8]. A small number of tumor cells produced by successive non-lethal mutations begin to inter- act with the community of normal cells and do not recognize, thus triggering the acquisition of space and the vital resources of the existing community. In the case of tumor dynamic saturation in the growth of various types of tumors is not well modeled by the exponential model. For this reason, this model applies only to avascular tumors, i.e., when angiogenesis has not occurred, which under Kerbel [11], have around 1 to 2 mm diameter [21, 22]. Indeed, tumor cells compete for oxygen and vital resources that is the reason why the logistic model fits well in several cases[12, 25]. In the literature, one realizes that tumor growth does not follow a universal law [20], so we can mention another two of the most used models6: 6There are several other models that describe the dynamic of cancer tumor that will not be considered here [8], due to the fact that they are very similar to those that we considered. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 218 — #8 � � � � � � 218 A PRELUDE TO THE FRACTIONAL CALCULUS APPLIED TO TUMOR DYNAMIC i) The so called Generalized Logistic model is given by [8] d N(t) dt = r θ N(t) ( 1 − N(t) k )θ . This model describes well the growth of breast cancer [23]. ii) The Gompertz model [8], d N(t) dt = −r N(t) ln ( N(t) k ) , best fits the volumetric growth in vivo [13]; Below, we present some graphs of models in the literature that describe the increase in tu- mor mass/volume over time, from the simple Exponential growth model from one parameter to the more sophisticated models such as Gompertz, Generalized Logistic and Fractional Logistic Equation. 0 500 1000 1500 2000 2500 3000 0 100 200 300 400 500 600 700 800 900 1000 Tumor growth models Time (days) N um be r of tu m or al c el ls ( bi lli on s) Logistic Exponential Generalised Logistic Gompertz Fractional Logistics order 0,97 Fractional Logistics order 0,99 Figure 2: Tumor growth in humans according to the Exponential Models, Logistic, Gompertz and Logistic Fractional. In the chart above tumor growth in humans we consider the initial population N(0) = 4x109, a carrying capacity cells r = 1012, k = 10−2/day and for the generalized logistic model θ = 2 with scale 1 : 1.000.000.000. One disadvantage of all the usual models is that they go to the carrying support faster than what is expected [5, 6, 7, 8, 9, 22]. This is one of the advantages of the fractional logistic model, i.e., we note that as the order of the derivative decreases the amount of convergence to the support becomes slower. This slow convergence to the carrying support is consistent with the growth of some types of cancer tumor [8] making this highly relevant for the study of tumor dynamic equa- tion, since this model include competition between tumor cells for vital resources and predicts that the maximum size of a tumor takes longer to be achieved. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 219 — #9 � � � � � � VARALTA, VELLASCO and CAMARGO 219 5 CONCLUDING REMARKS AND FUTURE WORK The Fractional Calculus is nearly as old as the entire order one, but just in end of last century this valuable tool became evident, especially in modeling phenomena that possess time dependence, since fractional derivatives excellently describe memory effects when encountered derived from their entire order. In this work, we had the purpose of emphasizing the importance of Fractional Calculus to gener- alize and refine the solution of differential equations of integer order, these being of paramount importance in describing various phenomena in biological field. With this, we propose a general- ization via Fractional Calculus for classical logistic equation and show that the fractional model may provides a better description than the usual on, in the study of cancer tumor growth. One important aspect is to analyze if the solutions of the fractional versions of equations (3.1) and (3.3) are the same. Despite having scientific value the search for an answer to the previous question, is not indispensable for the purpose of the article, since equations (3.1) and (3.3) are equivalents, to consider the fractional growth we may to choose one of them. We are working on an answer for this question and we intend to publish it in a future work. A natural continuation of this work is to study the system presented by Rodrigues [21, 22] in fractional version, since we know the behavior of tumor free from chemotherapy we will use this knowledge to analyze the increase or decrease in tumor cancer on the action of chemotherapeutic agents. RESUMO. Com o intuito de refinar a solução da clássica equação logı́stica e ampliar seu campo de aplicação para o estudo da dinâmica tumoral, propomos e resolvemos a gene- ralização fracionária desta equação, utilizando o assim chamado cálculo fracionário, i.e., substituı́mos a derivada ordinária da equação usual por uma derivada fracionária de ordem 0 < α ≤ 1, e recuperamos a solução clássica como um caso particular. Por fim, analisamos a aplicabilidade deste modelo para descrever o crescimento de tumores de câncer. Palavras-chave: biomatemática, cálculo fracionário, equação logı́stica, dinâmica tumoral de câncer. REFERENCES [1] R.F. Camargo, E.C. Oliveira & Ary O. Chiacchio. Teorema de Adição para as Funções de Mittag- Leffler. TEMA, 10(1) (2009), 1–8. [2] R.F. Camargo. Cálculo Fracionário e Aplicações. Tese de Doutorado, IMECC, UNICAMP (2009). [3] R.F. Camargo, Ary O. Chiacchio & E. Capelas de Oliveira. Differentiation to Fractional Orders and the Fractional Telegraph Equation. J. Math. Phys., 49, 033505, (2008). [4] R.F. Camargo, Ary O. Chiacchio & E. Capelas de Oliveira. On anomalous diffusion and the fractional generalized Langevin equation for a harmonic oscillator. J. Math. Phys., 50, 123518, (2009). Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 220 — #10 � � � � � � 220 A PRELUDE TO THE FRACTIONAL CALCULUS APPLIED TO TUMOR DYNAMIC [5] A.M.A El-Sayed, A.E.M. El-Mesiry & H.A.A El-Saka. One the fractional-order logistic equation. Applied Mathematics Letters, 20 (2007), 817–823. [6] G.H. Erjaee, M. Shahbazi & A. Erjaee. Dynamical analysis of mathematical model presented by fractional differential equations, describing the interaction between leukemic cancer cells, T cells and drug treatment with a drug optimal control. Open Access Scientific Reports, 1 (2012), 1–8. [7] U. Forys & A. Marciniak-Czochra. Logistic equations in tumour growth modelling. International Journal of Applied Mathematics and Computer Science, 13 (2003), 317–325. [8] R.A. Gatenby & T.L. Vincent. Application of quantitative models from population biology and evo- lutionary game theory to tumor therapeutic strategies. American Association for Cancer Research, 2 (2013), 919–927. [9] P. Gerlee. The model muddle: in search of tumour growth laws. Cancer Research, 73 (2013), 2407– 2411. [10] R. Gorenflo & F. Mainardi. Fractional Calculus: Integral and Differential Equations of Fractional Order. CISM Lectures Notes, (2000), 223–276. [11] R.S. Kerbel. Tumor angiogenesis: past, present and the near future. Carcinogenesis, 21 (2000), 505– 515. [12] A. Lotka. Meeting on the problem o forecasting city populations with special reference to New York city. Journal of the American Statistical Association, 20 (1925). [13] S. Michelson, A.S. Glicksman & J.T. Leith. Growth in solid heterogeneous human colon adenocarci- nomas: comparison of simple logistical models. Cell Prolif, 20 (1987), 343–355. [14] G.M. Mittag-Leffler. Une generalisation de l. integrale de Laplace-Abel. Comptes Rendus de l’Academie des Sciences Serie II, 137 (1903), 537–539. [15] G.M. Mittag-Leffler. Sur la nouvelle fonction Eα(x). Comptes Rendus de l’Academie des Sciences, 137 (1903), 554–558. [16] G.M. Mittag-Leffler. Sur la representation analytiqie d’une fonction monogene (cinquieme note). Acta Mathematica, 29(1) (1905), 101–181. [17] E.C. Oliveira. Funções Especiais com Aplicações. Editora Livrariada Fı́sica, São Paulo (2005). [18] C. Phipps. Combination of Chemotherapy and Antiangiogenic Therapies: A Mathematical Modelling Approach. Thesis the Master Degree, University of Waterloo, Canada (2009). [19] I. Podlubny. Fractional Differential Equation – An Introduction to Fractional Derivates, Fractional Differential Equations, to Methods os their Solution and some of their Applications. Academic Press, San Diego, 198 (1999). [20] M. Retsky. Letter to the editor. J. Theor. Biol., 229 (2004), p. 289. [21] D.S. Rodrigues. Modelagem matemática em câncer: dinâmica angiogência e quimioterapia anti- neoplásica. Dissertação de Mestrado, UNESP, IBB (2011). [22] D.S. Rodrigues, P.F.A. Mancera & S.T.R. Pinho. Modelagem Matemática em Câncer e quimioterapia: uma introdução. Notas em Matemática Aplicada, SBMAC, São Carlos – SP, Brasil, Volume 58, e- ISSN 2236–5915, (2011). [23] J.S. Spratt, J.S. Meyer & J.A. Spratt. Rates of growth of human neoplasms: part II. J. Surg. Oncol., 61 (1996), 68–73. Tend. Mat. Apl. Comput., 15, N. 2 (2014) � � “main” — 2014/9/1 — 11:35 — page 221 — #11 � � � � � � VARALTA, VELLASCO and CAMARGO 221 [24] V.G. Vaidya & F.J. Alexandro-Jr. Evaluation of some mathematical models for tumor growth. Int. J. Bio-med. Comp., 13 (1982), 19–35. [25] P.F. Verhulst. Notice sur la loi que la population poursuit dans son accroissement. Correspondance mathématique et physique, 10 (1838), 113–121. [26] A. Wiman. Über den fundamentalsatz in der theorie der funktionen Eα(z). Acta Math., 29 (1905), 191–201. Tend. Mat. Apl. Comput., 15, N. 2 (2014)