Research Article  Open Access
Bahatdin Daşbaşı, "Stability Analysis of Mathematical Model including PathogenSpecific Immune System Response with FractionalOrder Differential Equations", Computational and Mathematical Methods in Medicine, vol. 2018, Article ID 7930603, 10 pages, 2018. https://doi.org/10.1155/2018/7930603
Stability Analysis of Mathematical Model including PathogenSpecific Immune System Response with FractionalOrder Differential Equations
Abstract
In this study, the mathematical model examined the dynamics between pathogen and specific immune system cells (memory T cells) for diseases such as chronic infection and cancer in which nonspecific immune system cells are inadequate to destroy the pathogen and has been suggested by using a system of the fractionalorder differential equation with multiorders. Qualitative analysis of the proposed model reveals the equilibrium points giving important ideas about the proliferation of the pathogen and memory T cells. According to the results of this analysis, the possible scenarios are as follows: the absence of both pathogen and memory T cells, only the existence of pathogen, and the existence of both pathogen and memory T cells. The qualitative analysis of the proposed model has expressed the persistent situations of the disease where the memory T cells either do not be able to respond to the pathogen or continue to exist with the diseasecausing pathogen in the host. Results of this analysis are supported by numerical simulations. In the simulations, the timedependent size of the tumor population under the pressure of the memory T cells was tried to be estimated.
1. Introduction
For three centuries, the theory of fractional derivatives was developed as a pure theoretical field of mathematics, useful only for mathematicians. But, the use of fractionalorders differential and integral operators in mathematical models has become increasingly common of late years. Therefore, various forms of fractionalorder differential equations are suggested for standard models. In this sense, the fractionalorder calculus plays an important role in physics [1], thermodynamics [2], viscoelasticity [3], electrical circuits theory [4], fractances [5], mechatronics systems [6], signal processing [7], chemical mixing [8], chaos theory [9], engineering [10], biological system [11], and other applications [12]. Also, a large number of literatures on the application of fractionalorder differential equations (FODEs) in nonlinear dynamics have been improved. Especially, when the biological applications of FODEs have considered, it is a rich source for mathematical ideas [13].
The mathematical modeling of diseases in biological applications is a subject discussed in the literature. Such models are considered under two main headings as by modeling the size of the spread of infected individuals in a population () and modeling the population size of the pathogens such as the tumor in an individual as it is here. The word tumor simply refers to a mass. This is a general term that can refer to benign (generally harmless) or malignant (cancerous) growths. Many types of tumors are considered to be a major factor in many fatal diseases in human history. Fundamentally, it is said that this disease is a complex process for both tumor and host. Although different treatment strategies are proposed for tumors, the first and foremost role in disease progression belongs to the immune system of the individual (or host) [14]. The immune system is stated as a system of biological structures and processes in an organism that protects the body from the possible hazardous organism by recognizing and responding to antigens. In more detail, the immune system cells such as T cells are generally described in terms of two different types. These are the effector and the memory of T cells. The ordinary behaviour of the immune system is generally an acute infection, controlled initially by effector T cells (aspecific response or the innate immune system response), later by memory T cells (specific response or the adaptive immune response), until complete clearance of the pathogen. The adaptive immune response is much slower to respond to threats and infections than the innate immune response, which is primed and ready to fight at all times [15]. Especially, T cells are a crucial component of the adaptive immune response against malignancies. Antigenexperienced T cells specific for tumor antigens can be recovered from the blood, lymphoid organs, and tumors of both cancer patients and tumorbearing mice [16]. Concordantly, the reactions of different hosts in case of the same disease may be different because of the immune system response given by host, and so, the disease progression varies from person to person. Within this context, dynamics of relevances between immune systems cells (T cells) and tumor are significant to find out the nature of the disease. The problem is to try to obtain the known biological features without making the mathematics too complicated.
The basic of a most useful explanation of fractional calculus is memory concept. If the output of a system at each time depends only on the input at time , then such systems are said to be memoryless systems. On the contrary, if the system has to remember previous values of the input in order to determine the current value of the output, then such systems are said memory systems [17, 18]. Accordingly, the behaviour of most biological systems has memory or aftereffects. The modeling of these systems by FODEs has more advantages than classical integerorder modeling, in which such effects are neglected. Also, FODEs are, at least, as stable as their integer order counterpart [11]. In the process of modeling reallife situations, the created mathematical models by using the fractionalorder differential operations allow to display the some extra cases regarding the stability region of the equilibrium point of the mathematical model caused by parameters such as derivative orders. For this reason, the mathematical models formed by FODEs are more realistic and feasible [19]. Additionally, stability analysis of equilibrium points for mathematical models consisting of FODEs with multi‐orders and its systems is more general than those of the sameorders too.
Although there are many studies that examined the dynamics between tumor and immune system response, the proposed model in this study differs from them in terms of both mathematical structure such as the use of Holling function type2 (functional and numerical responses) in the model consisting of the FODE system with multiorders and examination of qualitative analysis of the proposed model. In this sense, it was tried to bring a different perspective from the previous studies.
In this study, a FODE model with multiorders considering the basic mechanisms of tumor and the memory T cells having functional and numerical responses, respectively, has been constructed, and so, the qualitative analysis of the proposed model was performed. The reason for using the Holling function type2 is to show the limit cycle behaviour of system [20]. The certain conditions dependent on the development of the tumor population under the pressure of memory T cells was obtained. In this respect, all of the possible scenarios related to the tumor size were tried to be explained as parameterdependent. Additionally, numerical analysis of the model was given as to be compatible with the qualitative analysis.
2. Preliminaries and Definitions
In here, the main definitions and properties of fractional derivative operators have been expressed. Also, the FODE systems with multiorders have been introduced, and the properties such as stability and existence of the equilibrium points of such systems are given.
2.1. Fractional Differential Operators
There are various definitions of a fractional derivative with the order . The definitions of Riemann–Liouville and Caputo are used most widely. The Caputo sense was used in this study. Taking into account the definition of Caputo sense, the fractional derivative of the function is identified asfor , , [21].
2.2. The FODE System with Multiorders
Let us consider that is the time parameter. We have assumed that the system of FODE with multiorders is given as the following equation:where the variable , the initial conditions by , the functions by and for , and the derivative orders by .
Also, when it is considered as , indicates thorder fractional derivative in the Caputo sense. In this sense, it is . The multiorders can be mathematically any real or complex vector. In this study, the real case was only taken into account. Throughout the paper, we restrict to a rational number in the interval [22].
Remark 1. From (2), we have assumed thatwhere the independent variable is not clearly seen in the function . The equilibrium point of (3) is the point obtained from the equations .
Remark 2. For each equilibrium point of the autonomous system in (3), the eigenvalues obtain from the following equation:where is the Jacobian matrix evaluated at the equilibrium point and is the smallest of the common multiples of the denominators of the rational numbers [23].
Theorem 1. For each equilibrium point of system (3), is locally asymptotically stable (LAS), if the eigenvalues obtained from (4) satisfy Routh–Hurwitz Stability Criteria or the inequalities for . Here, has been defined in Remark 2 [22]. Because the 2dimensional of system (3) is used in the proposed model in this study, the stability analysis of such systems are described in detail below.
Remark 3. Let us assume that the autonomous system of FODE with multiorders is as following:with the nonnegative initial conditionswhere the derivative orders and are rational numbers in the interval . The equilibrium point of system (5) is the point obtained from the equations for . To evaluate locally asymptotically stability (LAS) of equilibrium point, the Jacobian matrix, , is used. Considering Remark 2, the eigenvalues for are obtained from the following equation:Therefore, the characteristic equation for eigenvalues isIf the eigenvalues for satisfy Routh–Hurwitz stability criteria or the conditions,then the equilibrium point is the LAS point for system (5).
For the system of FODE with multiorders and , the stability region is as shown in Figure 1 (where and are the real and imaginary parts of the eigenvalues, respectively, and ). By Figure 1, we openly see that the stability region of the equilibrium point of the FODE with multiorders is greater than the stability regions of the integerorder case and the same fractionalorder case [24].
Remark 4. Let in system (5). In this case, we have the systemwith the nonnegative initial conditions . From the equations for , we have presumed that the equilibrium point of system (10) is . If the eigenvalues and obtained from the equationprovide the conditionsthen the equilibrium point is the LAS point for system (10).
Conditions expressed in (12) can be detailed as the followings. Characteristic equation of (11) is the following generalized polynomial:When both the conditions (12) and the polynomial (13) are considered together, the conditions for LAS of the equilibrium point are either Routh–Hurwitz conditions [25, 26]:or
3. Model Formulation
The proposed model is particularly well suited for describing diseases such as chronic infection and cancer in which the nonspecific immune system cells are inadequate to destroy the pathogen. Consequently, it has been proposed, and another extension of the models in [20, 27–33] has been analyzed.
It has been identified mathematically as a pathogen load, specifically tumor population, and level of memory T cells, namely, the adaptive immune response or specific response, in an individual. In this sense, we have assumed that the population densities of pathogen and memory T cells at time are denoted by and , respectively. Additionally, the memory T cells predate the tumor cells by a Holling function type2.
Under the assumptions aforementioned, we have proposed the following system of FODE with multiorders and :where for indicates thorder fractional derivatives in the Caputo sense, it is and , and the parameters have the following properties:
In addition that, system (16) has to be finished with positive initial conditions and . The parameters used in the model are defined as follows.
It is presumed that the pathogen follows a logistic growth rule with the carrying capacity and the growth rate . The memory T cells proliferate proportionally to the pathogen load by the Holling function type2. Since the pathogen capture rate of memory T cells is assumed to be proportional to the per capita growth rate of memory T cells, the constant represents the maximum growth rate for memory T cells and the constant is the pathogen population size at which the growth rate of memory T cells in half of its maximum. These situations are very suitable for the growth of memory T cells especially in case of chronic infection or tumor. Memory T cells have per capita natural death rate . Moreover, the pathogen die due to the action of the memory T cells, and we have presumed that the effect of these cells on pathogen is modeled using a saturating response, , subject to a maximum killing rate and the level of memory T cells required for the half maximum effect, .
Remark 5. Rate of replication of the immune system cells is higher than its death rate, at least every time the pathogen load is very high [20]. In this case, we have by (16). Therefore, the following inequality is obtained:by this limit.
Proposition 1. System (16) provides the followings. The freedisease equilibrium point and the equilibrium point , where only the pathogen exists and always exist. In addition to and , there exists a third equilibrium point as forwhen .
Proof. The steady states of the model (16) are again the intersection of null clines and in (16). We have accepted that the solutions of theses equations consist of the pairs . Then, we have the following system:From the first equation of (20), it is or . Let , and then . Therefore, the system (16) has the freedisease equilibrium point . On the contrary, let , that is, . If the value is rewritten in the second equation of system (20), then we have found the equilibrium points and . Considering (19), if is rewritten, the point is obtained. always exists due to (17). Let us consider . is positive due to (17) and (18). On the contrary, is positive due to (17), whenTherefore, we have a positive equilibrium point where and are in (19).
In Table 1, biological existence conditions of equilibrium points of system (16) are showed.

Proposition 2. In system (16), let us consider derivative orders aswhere the smallest common multiple of and is . System (16) satisfies the following:(a) is a unstable point(b) is LAS, when . Also, if , then this point is a unstable point(c) where and defined in (19) is LAS, when all roots for found from the equationsatisfy Routh–Hurwitz stability criteria or the condition . In here, it is
Proof. For the stability analysis of the equilibrium points, the functions in system (16) are assigned asIn this respect, the Jacobian matrix evaluated at each equilibrium point showed in Table 1 isfor .(a)For , the Jacobian matrix in (26) isNow, we will investigate whether the inequality (9) has been achieved. By (27), the eigenvalues have been found from the following determinant:In this respect, the characteristic equation obtained from (28) isThus, we have and . is real positive due to (17). Moreover, it is obtained asFor the eigenvalues in (30), it is for . These eigenvalues are positive real number on the right side of the complex plane, and so, it is . Since the stability condition is not supplied, the equilibrium point is a unstable point for system (16).(b)From (24), the Jacobian matrix related to isFrom the equation , the characteristic equation of eigenvalues isTherefore, it is and . These equations are examined as the following:
(1) is a negative real number due to (17). By De Moivre’s formula, we have , such that . Considering (9), the stability condition for is due to . This condition has been always provided since in (16).(2)On the other hand, we have considered the equation . If , then for is obtained from by De Moivre formulas. In this respect, the stability condition is by . This condition has been always provided from (16). Additionally, If , then the eigenvalues are positive real number due to for . In this sense, we have . The stability condition is not provided due to . Therefore, the equilibrium point is an unstable point for system (16).Consequently, if the inequalityis provided, then the equilibrium point is LAS, and if , then this point is an unstable point for system (16). When (33) is rearranged, the stability condition of is(c)Jacobian matrix evaluated at isthat is,where and are defined in (24).Let us give more details for and . Accordingly, it issince the value , the pathogen size may take, is less than or equal to its carrying capacity . In addition, it isdue to (17), and the components of the equilibrium point in are positive.
From (36), we have the characteristical equation as follows:To be LAS of , it should be that all roots for found from the (39) satisfy the inequalities in (9).
Proposition is proved.
Corollary 1. Equation (39) can be examined in more detail as shown below. This equation can be rewritten by De Moivre formulas such thatwhere , angle , and . By (40), (39) transforms toand so,By arranging (42), there is the following system:and so,From the first equation in system (44), we have foundBy substituting (45) in the second equation in (44), it is foundConsequently, if the angles for obtained from (46) satisfy Routh–Hurwitz stability criteria or the condition (9) , then is LAS.
For equilibria of system (16), the conditions found for LAS and biological existence are summarized in Table 2.
Corollary 2. is an unstable point, when exists biologically. Therefore, these equilibrium points cannot be stable when together. Similarly, is biologically meaningless, when is LAS. These circumstances appeared are also seen in Table 2.
Corollary 3. Let us consider the special case of for . In this case, we have Remark 4. The characteristical equation obtained from isThe conditions for LAS of the equilibrium point are either Routh–Hurwitz stability conditions.due to (37) and (38), or the conditionsIn addition that, by considering equation (46), there may be found a similar condition for stability of this point. In this sense, it isby formulas of half angle. Thereby, the LAS condition of iswhere . Corollary 3 is summarized in Table 3.
4. Numerical Analysis for Model
The proposed model in this study summarizes the general dynamics of the pathogenimmune system. The immune system cells described herein are memory T cells specially produced by the host against the pathogen. Hence, the proposed model is suitable for modeling of diseases such as chronic infections or tumors in which the nonspecific immune system cells of the host at the beginning of the disease have failed to destroy the pathogen.
The conditions found in Table 2 have been supported by numerical studies shown below. In this section, the timedependent sizes of the tumor and memory T cells for cancer tried to be estimated by giving the different values to the parameters in the proposed model. The reason for this is to be able to obtain different scenarios and to better demonstrate the results of qualitative analysis. The values of parameters used in system (16) are shown in Table 4.

The values calculated by Table 4 are given in Table 5.
