Mathematical Modelling for the Role of CD4<sup>+</sup>T Cells in Tumor-Immune Interactions

Makhlouf, Ahmed M.; El-Shennawy, Lamiaa; Elkaranshawy, Hesham A.

doi:https://doi.org/10.1155/2020/7187602

Computational and Mathematical Methods in Medicine

On this page

Abstract Introduction Materials and Methods Results and Discussion Conclusions Data Availability Disclosure Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 7187602 | https://doi.org/10.1155/2020/7187602

Mathematical Modelling for the Role of CD4⁺T Cells in Tumor-Immune Interactions

Ahmed M. Makhlouf,¹Lamiaa El-Shennawy,²and Hesham A. Elkaranshawy¹

Academic Editor: David Diller

Received04 Oct 2019

Revised17 Dec 2019

Accepted20 Jan 2020

Published19 Feb 2020

Abstract

Mathematical modelling has been used to study tumor-immune cell interaction. Some models were proposed to examine the effect of circulating lymphocytes, natural killer cells, and CD8⁺T cells, but they neglected the role of CD4⁺T cells. Other models were constructed to study the role of CD4⁺T cells but did not consider the role of other immune cells. In this study, we propose a mathematical model, in the form of a system of nonlinear ordinary differential equations, that predicts the interaction between tumor cells and natural killer cells, CD4⁺T cells, CD8⁺T cells, and circulating lymphocytes with or without immunotherapy and/or chemotherapy. This system is stiff, and the Runge–Kutta method failed to solve it. Consequently, the “Adams predictor-corrector” method is used. The results reveal that the patient’s immune system can overcome small tumors; however, if the tumor is large, adoptive therapy with CD4⁺T cells can be an alternative to both CD8⁺T cell therapy and cytokines in some cases. Moreover, CD4⁺T cell therapy could replace chemotherapy depending upon tumor size. Even if a combination of chemotherapy and immunotherapy is necessary, using CD4⁺T cell therapy can better reduce the dose of the associated chemotherapy compared to using combined CD8⁺T cells and cytokine therapy. Stability analysis is performed for the studied patients. It has been found that all equilibrium points are unstable, and a condition for preventing tumor recurrence after treatment has been deduced. Finally, a bifurcation analysis is performed to study the effect of varying system parameters on the stability, and bifurcation points are specified. New equilibrium points are created or demolished at some bifurcation points, and stability is changed at some others. Hence, for systems turning to be stable, tumors can be eradicated without the possibility of recurrence. The proposed mathematical model provides a valuable tool for designing patients’ treatment intervention strategies.

1. Introduction

Cancer is one of the leading causes of death worldwide. According to the World Health Organization (WHO), there were 8.8 million deaths in 2015 due to cancer [1]. The global cancer burden is expected to rise to nearly 21.4 million cases and 13.5 million deaths by 2030 [2]. Cancer treatment includes surgery, radiation [3], hormonal therapy [4], virotherapy [5], chemotherapy, and recently immunotherapy [6]. Chemotherapy is a well-known method for treatment of cancer. It is based on the administration of drugs that can kill tumor cells. These drugs do not only kill tumor cells but can also kill normal cells. Thus, many patients suffer from side effects of treatment as well as resistance to therapy and recurrence [7].

New approaches to treatment have been investigated, and immunotherapy has been recently approved for the treatment of many types of cancers [8, 9]. Immunotherapy is based on enhancing the effectiveness of the immune system to identify and kill tumor cells using two strategies: (1) passive immunotherapy, where effector components of the immune system are used to directly attack tumor cells. This strategy includes antibody-targeted therapy and genetically engineered T cells (e.g., chimeric antigen receptor [CAR]-T). (2) Active immunotherapy, where the activity of the immune system is enhanced. This includes cancer vaccines, cytokines, and adoptive cell therapy [10]. Cancer vaccines enhance cytotoxic T lymphocytes response to specific antigens produced by the tumor cells. Cytokines are proteins important for cell signalling and are produced by many cells, such as macrophages, B lymphocytes, and T lymphocytes. However, not all cytokines are approved for the treatment of cancer. The Food and Drug Administration (FDA) in USA approved only two cytokines, interleukin-2 (IL-2) and interferon alpha (IFN-α), since they demonstrated moderate clinical benefit [11]. IL-2 is most commonly used to trigger the immune system by activating T-cells and natural killer cells that identify and kill tumor cells. INF-α has a similar response rate to IL-2, but it does not achieve long-term patient survival compared to IL-2 [2, 10, 12, 13]. Adoptive cell therapy (ACT) involves ex vivo stimulation and expansion of tumor infiltrated T cells, then infusing the cells back into cancer patients [10]. The majority of ACT clinical trials and animal studies comprise CD4⁺T helper 1 cells and CD8⁺ cytotoxic T cells [14].

Mathematical modelling provides a powerful tool to describe and analyse many engineering and physical problems. It has also been used to describe some biological processes such as heart beats [15], diffusion of drugs [15], hepatitis C virus [16], management of HIV infection [17], treatment of diabetes [18], clearing the antibiotic resistant infection [19], aortic aneurysm formation [20], and tumor growth and cancer treatment. Studies in the field of cancer biology focus on developing suitable mathematical models involving quantitative approaches to understand various aspects of tumor growth and the response of cancer cells to clinical interventions. Jones et al. proposed a model describing tumor growth due to its internal pressure [15]. The model was in the form of a system of partial differential equations and included tumor size and internal pressure and nutrients' concentration without any treatment. Similar models, in the form of systems of partial differential equations, were also proposed by Tao et al. [21], Wei and Cui [22], Wise et al. [23], Frieboes et al. [24], Lee et al. [25], Zhang et al. [26], Knopoff et al. [27], de Pillis et al. [28], Villasana and Radunskaya [29], and Xu and Bai [30]. The last two models involved time delayed equations as well. A statistical model for cancer mortality has been studied by Ghosh and Samanta [31].

An important division of the mathematical models about cancer intends to investigate the effect of varying the concentrations of immunotherapy and chemotherapy and to explore the influence of various parameters. Valuable and comprehensive reviews can be found in [32–35]. Most of these models are in the form of a system of coupled nonlinear ordinary differential equations (see, for example, references [12, 36–39]). Along this line, de Pillis and colleagues developed a useful model for the interaction between the tumor and immune system [13, 40–43]. They obtained a numerical solution and performed both stability analysis and sensitivity analysis. They took into account the following assumptions: (a) a tumor grows logistically in absence of an immune response, (b) both natural killer and CD8⁺T cells are capable of killing tumor cells, (c) both natural killers and CD8⁺T cells respond to tumor cells by expanding and increasing cytolytic activity, (d) natural killers are normally present in the body, even when no tumor cells are present, (e) active tumor specific CD8⁺T cells are only present in large numbers when tumor cells are present, and (f) both CD8⁺T cells and natural killers become inactive after some number of encounters with tumor cells [40, 42]. Arabameri et al. [8] proposed a model for immunotherapy using mature dendritic cells. Mamat et al. [2] proposed a modified model that includes the effect of both IL-2 and INF-α. A model was proposed by Eftimie et al. [44] and was modified by Anderson et al. [45] and by Hu and Jang [46] to study the effect of CD4⁺T cells on the system. CD4⁺T cells are helper cells that activate CD8⁺T cells to achieve a regulated immune response and protect against tumors. However, CD4⁺T cells can play a direct role in killing tumor cells via secretion of cytokines in addition to its traditional helper role [45–49].

Intriguingly, previous studies modelled the effect of natural killer cells, circulating lymphocytes and CD8⁺T cells, and neglected the role of CD4⁺T cells, or studied the role of CD4⁺T cells, neglecting the role of other immune cells. In this study, a mathematical model is proposed for active immunotherapy that includes the effects of natural killer cells, circulating lymphocytes, CD8⁺T cells, CD4⁺T cells, and chemotherapy on the tumor. Both roles of CD4⁺T cells are considered in the model. The work of de Pillis et al. [13, 40–42], Eftimie et al. [44], Anderson et al. [45], and Hu and Jang [46] are the main references for the proposed model. The model is solved numerically using the “Adams predictor-corrector” method with Adams–Bashforth predictor and Adams–Moulton corrector. A stability analysis is performed for patients' data to determine if the system is stable. For a nonzero tumor case, the tumor is totally eradicated by the treatment intervention and the possibility of recurrence is to be studied. At the zero tumor case, the system must be stable to avoid recurrence and a condition for stability is deduced. It is recommended to work on modifying the patient’s body parameters, using cancer vaccines, for example, before applying the treatment intervention in order to satisfy the condition needed for stability [42]. Finally, bifurcation analysis is performed to study the effect of changing patient’s parameters on stability.

2. Materials and Methods

2.1. Mathematical Model

To illustrate the interaction between tumor cells and immune cells present in the microenvironment of the tumor without external effects, such as treatment, the following model is proposed:where is the natural killer cell population, is the CD8⁺T cell population, is the CD4⁺T cell population, is the circulating lymphocyte cell population not including natural killer cells and active CD8⁺T and CD4⁺T cells [50, 51], is the concentration of the IL-2 cytokine in the body, and is defined asand the parameters , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and are given in Table 1. The role of CD4⁺T cells is included in this model, and it depends mainly upon the work of Eftimie et al. [44], Anderson et al. [45], and Hu and Jang [46]. In this model, CD4⁺T is included in equation (4), and its effect is included in the terms and in equations (1) and (6) [46, 54]. In equation (1), represents tumor growth and is assumed to be logistic [40, 42, 46]. The term represents the natural killer-induced tumor death [40, 42]. The next term is the tumor lysis by CD8⁺T cells [40, 42]. Unlike CD8⁺T cells, CD4⁺T cells cannot kill tumor cells directly but through the cytokines that they produce. The role of cytokine in killing the tumor is represented by the term in equation (1) [46, 54]. In equation (2), represents the production of natural killer cells from circulating lymphocytes [40, 42]. This growth term is tied to the overall immune health levels as measured by the population of circulating lymphocytes [54–56]. represents natural killers’ turnover [51]. The term represents the recruitment of natural killer cells induced due to the existence of the tumor [42, 54, 57]. The next term is natural killer death by exhaustion of tumor-killing resources [40, 42, 51]. In equation (3), represents the inactivation of CD8⁺T cells consisting only of natural death rate since no CD8⁺T cells are assumed to be present in the absence of tumor cells [40, 42]. The term represents the recruitment in CD8⁺T cell population due to the presence of the tumor [40, 42, 58]. The next term is the CD8⁺T cell death by exhaustion of tumor-killing resources [40, 42, 51]. The term represents CD8⁺T cell stimulation by natural killer-lysed tumor cell debris, and the term represents the activation of native CD8⁺T cells in the general lymphocyte population [42, 51]. The term is the natural killer cell regulation of CD8⁺T cells. It describes regulation and suppression of CD8⁺T cell activity, which occurs when there are very high levels of activated CD8⁺T cells without responsiveness to cytokines present in the system [40, 42, 59, 60]. The activating effect of cytokines on CD8⁺T cells is represented in equation (3) by the term [40, 42, 61].

CD4⁺T cell population is represented by equation (4) which has a constant decaying rate . The other decaying term exists due to the interaction with the tumor. Activating effect of cytokines on CD4⁺T cells due to existence of the tumor is represented by the term in equation (4) [46, 54]. Equation (5) represents the circulating lymphocytes cell population which are assumed to be generated at a constant rate and death rate [40, 42, 54–56]. Equation (6) represents the concentration of cytokines which is decaying with a constant rate . The activating role of the CD4⁺T cells on cytokines due to the existence of the tumor is represented by the term in equation (6) [46, 54].

The previous model represents the interaction of patient’s body with the tumor only without any treatment. It can be generalized to introduce the effects of chemotherapy and/or immunotherapy; refer, for example, de Pillis and Radunskaya [40] and Hu and Jang [46]. The general model takes the following form:such that , , , and are the concentrations of chemotherapy, IL-2 cytokine therapy, CD8⁺T, and CD4⁺T adoptive immunotherapy, respectively. The effect of chemotherapy on immune cells is represented by the terms , , and . The three terms illustrate the reduction in number of each type of immune cell populations, while tumor cell killing due to chemotherapy is represented by and the concentration of the chemotherapeutic agent is represented by equation (14) [40, 42]. The parameters , , , , and are defined in Table 1.

2.2. Method of Solution

In numerical methods for solving differential equations, a problem may occur when small changes in the independent variable lead to large changes in the dependent variable. If this problem exists in a system of differential equations, the system is called a stiff system [62, 63]. The treatment doses , , , and change suddenly from zero to very large values within small intervals of , and systems (8)–(14) become a stiff system. We tried to solve the system using the Runge–Kutta method, but the solution was divergent in most cases due to the stiffness. The Adams–Bashforth method [62–65] is a powerful method for solving stiff problems and was originally proposed in 1883. For the following nonlinear initial value problem,

The approximate solution at is given bywhere is the step size and is the order of the method. The coefficients , , can be deduced using Taylor’s expansion of near the point with error of order [65]. This means that, at any point, we can find the solution numerically using any number of the previous points.

The Adams–Moulton method [62, 63, 65, 66] is a modification of the Adams–Bashforth method, which was proposed in 1926 by adding an implicit term to equation (16) to become

Coefficients can be found in the same way of Adams–Bashforth method.

Finally, a combined two-staged method is proposed by first computing the predictor stage using the Adams–Bashforth method [62, 65, 67]:and then the corrector stage using the Adams–Moulton formula:

This combined method is called the “Adams predictor-corrector method.” In this paper, we use it to solve the systems of nonlinear ordinary differential equations presented by equations (1)–(6) and (8)–(14).

3. Results and Discussion

The model proposed in this study is solved numerically using the Adams predictor-corrector method up to the 12^th order, and the solutions are interpolated. The method has nonzero stability ordinates for this order [68]. Values of the parameters , , , , ,, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and are given in Table 1 for humans. The solutions are obtained for three cases: without any treatment, with continuous treatment, and with pulsed treatment as detailed below.

3.1. First: A Case of No Treatment

In this case, we illustrate the interaction between tumor cells and the immune system without any treatment. In Figures 1(a) and 1(b), simulation was performed for first patient’s data where the initial natural killer cells , initial CD8⁺T cells , initial CD4⁺T cells , and initial circulating lymphocytes . The immune system could overcome tumor of size , as shown in Figure 1(a), but cannot overcome a tumor of size , as shown in Figure 1(b). The simulation is performed again with tumor size but with a stronger immune system such that , , , and , which means that initial CD8⁺T cells and initial circulating lymphocytes are increased. Figure 1(c) shows that tumor cells were destroyed in such case.

(a)

(b)

(c)

In this section, we demonstrate the role of different immune cells in abolishing the tumor. The results suggest that patient’s body can overcome cancer cells up to a certain size. In case of larger tumors, however, a stronger immune system is necessary to overcome cancer cells.

3.2. Second: A Case of Continuous Treatment

In this case, we simulated the model with continuous immunotherapy. This is performed using the first patient’s data, with initial natural killer cells , initial CD8⁺T cells , initial CD4⁺T cells , initial circulating lymphocytes , and initial tumor size . First, a continuous adoptive therapy is introduced using CD8⁺T cells with concentration and cytokine IL-2 with concentration . Figure 2(a) demonstrates that the tumor cannot be eliminated using this treatment. Second, a continuous adoptive therapy is introduced using CD4⁺T cells with concentration instead of both CD8⁺T cells and IL-2 cytokine therapies, and the tumor could be eliminated, as shown in Figure 2(b).

(a)

(b)

Treatment is applied using the second patient’s data with initial conditions: , , , , and tumor size . Figure 3(a) illustrates the results of using continuous CD8⁺T cell therapy together with cytokine IL-2 treatment , and Figure 3(b) shows the result when using only CD4⁺T cell therapy . In both cases, the tumor is not eliminated, but when using the three therapies together, , , and , the tumor is eliminated in about 10 days, as shown in Figure 3(c).

(a)

(b)

(c)

3.3. Third: A Case of Pulsed Treatment

In this case, the patient is administered therapy on certain days only. In Figures 4(b), 5(c), and 5(d), a tumor of initial size is assumed to exist in the first patient with initial conditions: , , , and . In each figure, we use a certain regimen of treatment. In Figure 4(a), five pulses of chemotherapy of intensity starting from the 6^th day with 5 days period can eliminate the tumor, as shown in Figure 4(b). An alternative to chemotherapy, 7 pulses of CD8⁺T cell treatment combined with IL-2 cytokine with 2 days period starting from day 6 were proposed, as illustrated in Figures 5(a) and 5(b). Figure 5(c) demonstrates that this combination is not enough to eradicate the tumor. Oscillations of CD4⁺T cell population and IL-2 cytokine concentration are experienced due to fast consumption of the treatment by the resisting tumor cells. These large changes in small time intervals explain the reason of the failure of the Runge–Kutta method in obtaining the numerical solution. When we use the same treatment schedule with additional 5 pulses of CD4⁺T cell treatment with 2 days period starting from day 6, the tumor is eliminated, as shown in Figure 5(d). These simulations suggest that CD4⁺T cells could replace chemotherapy and thus could spare patients the agony of chemotherapy-induced toxicities and side effects.

(a)

(b)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5

Simulation of pulsed treatment case for first patient’s data: (a) 7 pulses of CD8⁺T cell treatment = 10⁵ with 2 days period starting from day 6; (b) 7 pulses of IL-2 cytokine = 10⁵ with 2 days period starting from day 6; (c) pulsed CD8⁺T cells and IL-2 treatment with initial tumor size T(0) = 8 × 10⁵; (d) pulsed CD8⁺T cells and IL-2 combined with CD4⁺T cell treatments with initial tumor size T(0) = 8×10⁵; (e) pulsed chemotherapy with initial tumor size T(0) = 3 × 10⁶; (f) pulsed chemotherapy combined with CD4⁺T cell treatment with initial tumor size T(0) = 3 × 10⁶; (g) pulsed chemotherapy combined with CD8⁺T cells and IL-2 treatments with initial tumor size T(0) = 10⁷; (h) pulsed chemotherapy combined with CD4⁺T cell treatment with initial tumor size T(0) = 10⁷.

Pulsed treatment is applied to the first patient with initial conditions: , , , and but with a larger tumor size . In Figure 5(e), five pulses of chemotherapy of intensity is applied starting from the 6^th day with 5 days period. Chemotherapy alone is insufficient to overcome tumor growth. In Figure 5(f), the same chemotherapy drug is applied combined with 5 pulses of CD4⁺T cell therapy with 2 days period starting from day 6. Combined treatment could abolish the tumor, while chemotherapy alone was not able to do so. The results suggest that, depending on tumor size, CD4⁺T cell treatment can be an alternative to chemotherapy or is needed in combination with therapy.

De Pillis and Radunskaya [40] studied the case of a larger tumor of initial size . Pulsed treatment is applied to the first patient with initial conditions: , , and . They applied 9 pulses of chemotherapy of intensity starting from the 1^st day with 10 days period combined with 6 pulses of the IL-2 cytokine therapy from day 8 to day 11 and CD8⁺T cells treatment during days 7 and 8. Combined treatment could abolish the tumor in this case, as shown in Figure 5(g). As an alternative to this therapy, with the same initial conditions, however, we added the effect of CD4⁺T cells . We suggest 5 pulses of CD4⁺T cell therapy with 2 days period starting from day 6 combined with three pulses of chemotherapy of intensity applied starting from the 1^st day with 10 days period. Figure 5(h) illustrates that, in this case, the tumor can be also abolished. It is observed that we used a lower dose of the chemotherapeutic agent. Also, we used CD4⁺T cell therapy as an alternative to both CD8⁺T cell therapy and cytokine therapy.

3.4. Stability Analysis

In this section, stability analysis is performed for the system of equations (1) to (6). The main advantage of studying stability is to investigate the possibility of eradicating the tumor without the occurrence of recurrence. If an equilibrium point is stable, it will attract the system towards it, in its domain of attraction. However, if the equilibrium point is unstable, any small perturbation from equilibrium will cause the system to move away from this equilibrium point; i.e., upon terminating treatment, the tumor will escape immune surveillance unless every single tumor cell is killed. To avoid recurrence of tumor after treatment, the zero tumor case has to be stable. In this section, a condition for stability of the zero tumor case is deduced. The first step is equating the right-hand sides of these unforced autonomous equations to zero [69–72]. From equation (1), we can deduce that either or is at equilibrium. In the following part, the two equilibrium conditions are considered; first, the zero tumor equilibrium, , and then the nonzero tumor equilibrium, and .

3.4.1. Equilibrium and Stability for Zero Tumor

In this case, . Substituting in equations (4)–(6), the zero tumor equilibrium yields , , and . By substituting in equation (2) at equilibrium, it follows that . Finally, substituting in equation (3) at equilibrium, the result is , then , and the negative solution for is refused.

The zero tumor equilibrium point is

Linearizing system of equations (1) to (6) near the zero tumor equilibrium point, then

The Jacobian matrix [69–71] of this system is

The system has the eigenvalues , , , , , and . All the eigenvalues are always negative except the first eigenvalue which is negative only when . Then, the tumor-free state is a stable state if and unstable otherwise.

3.4.2. Equilibrium and Stability for Nonzero Tumor

In this case,then,

From equations (2) and (5) at equilibrium,

From equation (7) at equilibrium,

From equation (3) at equilibrium,

From equations (4) and (6) at equilibrium,

Then, from equation (28), the nonzero tumor equilibrium will be divided into two cases: the zero cytokine case at which and nonzero cytokine case at which and . In the nonzero tumor equilibrium with zero cytokine, substitution in equation (4) at equilibrium gives . Equilibrium points in this case can be found by solving equations (26) and (27) for and and then calculating other variables. Data for the two patients are used for numerical solution. Since equations (26) and (27) are highly nonlinear, the solution is obtained by plotting both of them in the plane and determining the intersections. In the nonzero tumor equilibrium with nonzero cytokine,

Rearranging,

The equilibrium points in this case can be found by solving equation (30) for . For each value of , equations (26) and (27) are solved for and and then substituted for other variables. Due to the high nonlinearity of the two equations, the solution is obtained by plotting both of them in the plane and determining the intersections. For each equilibrium point, a linearization method is used to write the Jacobian and to determine the eigenvalues. Consequently, the stability of the system is specified.

3.4.3. Results of Stability Analysis for Both the Zero Tumor and Nonzero Tumor

Table 2 illustrates the results for the parameter values given in Table 1. For more illustration, streamlines are plotted. We use the values of , , , and of each equilibrium point and allow and to vary with time according to equations (1) and (4). For the first patient, the streamlines are plotted for equilibrium point 1, , and , in Figure 6(a). One can notice that, for any state of the system in the vicinity of this unstable equilibrium point, the system moves away from that point. The streamlines are plotted for equilibrium point 2, and , in Figure 6(b). Also, for this unstable equilibrium point, any state of the system in the vicinity of this point, the system moves away from it. In the same manner, the streamlines for the second patient are obtained; however, to be concise, they are not included.

(a)

(b)

In this analysis, it is shown that all equilibrium points are unstable for both patients. This means that, for zero tumor case, the tumor can recur, and in nonzero tumor case, the size of the tumor is uncontrollable; i.e., it could increase without a limit. For the zero tumor, an analytical expression for the stability of the equilibrium point is obtained. Simply, if the condition is satisfied, then the zero tumor case is a stable state and the tumor cannot recur. If the condition is not satisfied, then the tumor-free state is an unstable state and the tumor is recurrent. Though for the nonzero tumor the stability analysis is analysed numerically, changing the system parameters could lead to stable equilibrium. This is tackled in the next section. Hence, it is recommended for any successful treatment to work on modifying the patient’s body parameters in order to satisfy the stability condition.

3.5. Bifurcation Analysis

A bifurcation analysis is performed to study the effect of changing patient’s parameters on the system stability [15, 40, 42, 73]. This change could lead to turning unstable point(s) into stable; hence, for nonzero tumor, the size is limited and can be reduced by surgery or radiation and then immunotherapy and/or chemotherapy to reach a stable zero tumor case. In the zero tumor case, it is required to change patient’s parameters, using cancer vaccines, for example, to satisfy the condition of stability stated in the previous section and prevent recurrence. Bifurcation diagram is plotted for parameters and for the first patient and the second patient. In the same manner, bifurcation analysis can be studied for other system parameters.

First, the value of parameter is changed from 0 to 10. Equilibrium points are monitored, and then bifurcation diagrams are plotted. In Figure 7(a), the bifurcation diagram for parameter with first patient’s data is illustrated. For , there is only one equilibrium point at which and the system is stable. The value is a bifurcation point because the system is turned into unstable after it. For , there is another equilibrium point due to the nonzero tumor case with zero cytokine. This point moves on curve 2 according to the value of , and it is always unstable. The boundaries of curve 2, and , are bifurcation points because this equilibrium point is created at and vanished at . For , there is only one equilibrium point at the zero tumor state and it is unstable.

(a)

(b)

(c)

(d)

In Figure 7(b), the bifurcation diagram for the parameter with second patient’s data is plotted. It shows that, for , there is only one equilibrium point at the zero tumor case and the system is stable. The value is a bifurcation point because the system is turned into unstable after it. For , there is another equilibrium point due to the nonzero tumor case with zero cytokine. This point moves on curve 2 according to the value of and is always unstable. The points and are bifurcation points. For the nonzero tumor with nonzero cytokine, there are two equilibrium points, illustrated by curves 3 and 4. The two equilibrium points created at the bifurcation points and , respectively. Curve 3 is always unstable and exists for . Curve 4 appears for , and it represents the case with small size tumor and is always stable. A healthy state might be one which maintains the system at this stable low-tumor size level. Parameter can be altered through treatments such as vaccination, and the tumor can be abolished/reduced by surgery or radiation and then immunotherapy and/or chemotherapy to reach the attraction zone of the tumor-free equilibrium state in the stable region of curve 1.

Bifurcation analysis is then performed for parameter . It is changed from 0 to . Figure 7(c) shows that, for first patient’s data, there are two equilibrium points for all values of . The equilibrium point moving on curve 2 for the nonzero tumor is always unstable. However, the equilibrium point moving on curve 1 for the zero tumor is unstable for and stable for . The point is a bifurcation point.

Figure 7(d) shows that, for second patient’s data, there are two equilibrium points for and three equilibrium points for . The equilibrium points moving on curves 2 and 3 are always unstable; however, the equilibrium point moving on curve 1 is unstable for and stable for . The point is a bifurcation point. Point at which curve 3 is created is a bifurcation point also. It can be deduced that changing the value of parameter leads to changing the stability of the zero tumor case only. For other equilibrium points, changing is not practically effective because other points are unstable. So, surgery is slight effective.

Bifurcation analysis performed in this section shows that stability of the tumor can be controlled by changing patient’s parameters. It is important to perform this analysis before deciding the treatment strategy and thus avoiding tumor recurrence.

4. Conclusions

In this study, a new mathematical model has been proposed in the form of a system of nonlinear ordinary differential equations to study the interaction between tumor cells and some immune cells. Effects of natural killer cells, circulating lymphocytes, CD8⁺T cells, and CD4⁺T cells on cancer cells have been investigated. Traditionally, CD4⁺T cells are considered helper cells that activate CD8⁺T cells. However, in the current study, the indirect role of CD4⁺T cells in killing tumor cells by secreting cytokines has been considered. Due to stiffness of the system of nonlinear ordinary differential equations, the model has been solved numerically using the method of Adams predictor-corrector with Adams–Bashforth predictor and Adams–Moulton corrector. The Adams method is usually used to solve a single nonlinear ordinary differential equation; however, in the current study, it is used to solve a system of nonlinear ordinary differential equations.

Three cases have been studied: first, the no-treatment case where interaction between immune cells and tumors with different sizes has been discussed. The results reveal that the ability of the patient’s body to overcome tumor cells depending upon the patient’s parameters, initial tumor size, and initial immune cells populations. For large tumors, the immune system alone fails to overcome the tumor and needs the intervention of treatment. Second, the case of continuous immunotherapy is considered where the role of CD4⁺T cells has been found to be very important and significant. In some cases, adoptive CD4⁺T cell therapy can be an alternative to combined CD8⁺T cells and cytokine therapy. In other cases, combined CD8⁺T cells and cytokine therapy is not enough to eliminate the tumor and introducing CD4⁺T cells is necessary. Third, the case of pulsed chemotherapy and immunotherapy is considered where CD4⁺T cell therapy has shown to be an alternative to chemotherapy in some cases. In other cases, specific pulses of chemotherapy alone could not be enough for treatment and immunotherapy could be introduced in order to completely eradicate the tumor. For large tumors, previous models abolished the cancer cells using intensive combined chemotherapy, CD8⁺T cells, and cytokine. In the current study, the same tumor has been abolished with a lower dose of chemotherapy combined with CD4⁺T cell therapy.

Stability analysis has been performed for the mathematical model. A condition for the stability of the zero tumor point has been deduced, and it depends on some immunological parameters of the patient. The stability analysis has answered the question of whether the tumor can recur after reaching size zero or not. If , then the zero tumor state will be a stable state and the tumor cannot recur. The results reveal that all equilibrium points are unstable for both the first and the second patients. This means that the size of the tumor is uncontrollable; i.e., it could increase without a limit. Also, after reaching the zero tumor equilibrium point using treatment, recurrence of the tumor can occur unless every single tumor cell is killed. We would like here to emphasize that the previously mentioned conclusions about the three studied cases would not be disturbed by the fact that the equilibrium points are unstable. For instance, it is still valid that CD4⁺T cells can replace chemotherapy or at least reduce its amount and, hence, spare patients the complications of chemotherapy-induced toxicities and side effects.

Finally, bifurcation analysis has been performed to show the effect of changing patient’s parameters on the stability. Bifurcation analysis has been considered for parameters and with the data for the first patient and the second patient. It has specified the values of and for the system to be stable, and hence, tumor cannot recur after reaching the zero tumor case. The analysis has also shown that, for the second patient, there is a nonzero small amount of tumor case which is stable. Bifurcation analysis provides a powerful tool to control the stability of the tumor by changing patient’s parameters.

The proposed model offers a powerful method for predicting the body response to immunotherapy and/or chemotherapy. Thus, the model provides a valuable tool for planning patients’ treatment intervention strategies.

Data Availability

Previously reported human data were used to support this study and are available in Table 1. These prior studies are cited at relevant places within Table 1.

Disclosure

The authors did not receive any funding for this research.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

References

WHO, “WHO,” April 2019, http://www.who.int/en/news-room/fact-sheets/detail/cancer.
View at: Google Scholar
M. Mamat, Subiyanto, and A. Kartono, “Mathematical model of cancer treatments using immunotherapy, chemotherapy and biochemotherapy,” Applied Mathematical Sciences, vol. 7, no. 5, pp. 247–261, 2013.
View at: Publisher Site | Google Scholar
R. Rockne, E. C. Alvord Jr., J. K. Rockhill, and K. R. Swanson, “A mathematical model for brain tumor response to radiation therapy,” Journal of Mathematical Biology, vol. 58, pp. 561–578, 2009.
View at: Publisher Site | Google Scholar
Y. Tao, Q. Guo, and K. Aihara, “A partial differential equation model and its reduction to an ordinary differential equation model for prostate tumor growth under intermittent hormone therapy,” Journal of Mathematical Biology, vol. 69, no. 4, pp. 817–838, 2014.
View at: Publisher Site | Google Scholar
J. Malinzi, “Mathematical analysis of a mathematical model of chemovirotherapy: effect of drug infusion method,” Computational and Mathematical Methods in Medicine, vol. 2019, Article ID 7576591, 16 pages, 2019.
View at: Publisher Site | Google Scholar
L. Pang, L. Shen, and Z. Zhao, “Mathematical modelling and analysis of the tumor treatment regimens with pulsed immunotherapy and chemotherapy,” Computational and Mathematical Methods in Medicine, vol. 2016, Article ID 6260474, 12 pages, 2016.
View at: Publisher Site | Google Scholar
T. P. Kingham, O. I. Alatise, V. Vanderpuye et al., “Treatment of cancer in Sub-Saharan Africa,” The Lancet Oncology, vol. 10, no. 4, pp. 158–167, 2013.
View at: Google Scholar
A. Arabameri, D. Asemani, and J. Hajati, “Mathematical modeling of in-vivo tumor-immune interactions for the cancer immunotherapy using matured dendritic cells,” Journal of Biological Systems, vol. 26, no. 1, pp. 1–22, 2018.
View at: Publisher Site | Google Scholar
A. Cesano and S. Warren, “Bringing the next generation of immuno-oncology biomarkers to the clinic,” Biomedicines, vol. 6, no. 14, pp. 14–24, 2018.
View at: Google Scholar
M. F. Sanmamed and L. Chen, “A paradigm shift in cancer immunotherapy: from enhancement to normalization,” Cell, vol. 175, no. 2, pp. 313–326, 2018.
View at: Publisher Site | Google Scholar
P. Berraondo, M. F. Sanmamed, M. C. Ochoa et al., “Cytokines in clinical cancer immunotherapy,” British Journal of Cancer, vol. 120, no. 1, pp. 6–15, 2019.
View at: Publisher Site | Google Scholar
S. Sharma and G. P. Samanta, “Dynamical behaviour of a tumor-immune system with chemotherapy and optimal control,” Journal of Nonlinear Dynamics, vol. 2013, Article ID 608598, 13 pages, 2013.
View at: Publisher Site | Google Scholar
L. G. de Pillis, K. R. Fister, W. Gu et al., “Optimal control of mixed immunotherapy and chemotherapy of tumors,” Journal of Biological Systems, vol. 16, no. 1, pp. 51–80, 2008.
View at: Publisher Site | Google Scholar
K. Li, B. Donaldson, V. Young et al., “Adoptive cell therapy with CD4+ T helper 1 cells and CD8+ cytotoxic T cells enhances complete rejection of an established tumour, leading to generation of endogenous memory responses to non-targeted tumour epitopes,” Clinical and Translational Immunology, vol. 6, no. 10, p. e160, 2017.
View at: Publisher Site | Google Scholar
D. S. Jones, M. J. Plank, and B. D. Sleeman, Differential Equations and Mathematical Biology, CRC Press, Boca Raton, FL, USA, 2nd edition, 2009.
H. A. Elkaranshawy, H. M. Ezzat, Y. Abouelseoud, and N. N. Ibrahim, “Innovative approximate analytical solution for standard model of viral dynamics: hepatitis C with direct-acting agents as an implemented case,” Mathematical Problems in Engineering, vol. 2019, Article ID 1454739, 8 pages, 2019.
View at: Publisher Site | Google Scholar
O. M. Ogunlaran and S. C. O. Noutchie, “Mathematical model for an effective management of HIV infection,” BioMed Research International, vol. 2016, Article ID 4217548, 6 pages, 2016.
View at: Publisher Site | Google Scholar
B. Shtylla, M. Gee, A. Do, S. Shabahang, L. Eldevik, and L. de Pillis, “A mathematical model for DC vaccine treatment of type I diabetes,” Frontiers in Physiology, vol. 10, 2019.
View at: Publisher Site | Google Scholar
P. A. Roberts, R. M. Huebinger, E. Keen, A.-M. Krachler, and S. Jabbari, “Mathematical model predicts anti-adhesion-antibiotic-debridement combination therapies can clear an antibiotic resistant infection,” PLoS Computational Biology, vol. 15, no. 7, Article ID e1007211, 2019.
View at: Publisher Site | Google Scholar
W. Hao, S. Gong, S. Wu et al., “A mathematical model of aortic aneurysm formation,” PLoS One, vol. 12, no. 2, Article ID e0170807, 2017.
View at: Publisher Site | Google Scholar
Y. Tao, Q. Guo, and K. Aihara, “A mathematical model of prostate tumor growth under hormone therapy with mutation inhibitor,” Journal of Nonlinear Science, vol. 20, no. 2, pp. 219–240, 2010.
View at: Publisher Site | Google Scholar
X. Wei and S. Cui, “Existence and uniqueness of global solutions for a mathematical model of antiangiogenesis in tumor growth,” Nonlinear Analysis: Real World Applications, vol. 9, no. 5, pp. 1827–1836, 2008.
View at: Publisher Site | Google Scholar
S. M. Wise, J. S. Lowengrub, H. B. Frieboes, and V. Cristini, “Three-dimensional multispecies nonlinear tumor growth-I,” Journal of Theoretical Biology, vol. 253, no. 3, pp. 524–543, 2008.
View at: Publisher Site | Google Scholar
H. B. Frieboes, F. Jin, Y.-L. Chuang, S. M. Wise, J. S. Lowengrub, and V. Cristini, “Three-dimensional multispecies nonlinear tumor growth-II: tumor invasion and angiogenesis,” Journal of Theoretical Biology, vol. 264, no. 4, pp. 1254–1278, 2010.
View at: Publisher Site | Google Scholar
H. G. Lee, Y. Kim, and J. Kim, “Mathematical model and its fast numerical method for the tumor growth,” Mathematical Biosciences and Engineering, vol. 12, no. 6, pp. 1173–1187, 2015.
View at: Publisher Site | Google Scholar
Y. Zhang, Y. He, and H. Chen, “Boundary element method for a free boundary problem modeling three dimensional tumor growth,” Computers and Mathematics with Applications, vol. 73, no. 7, pp. 1627–1641, 2017.
View at: Publisher Site | Google Scholar
D. A. Knopoff, D. R. Fernández, G. A. Torres, and C. V. Turner, “Adjoint method for A tumor growth PDE-constrained optimization problem,” Computers and Mathematics with Applications, vol. 66, no. 6, pp. 1104–1119, 2013.
View at: Publisher Site | Google Scholar
L. G. de Pillis, D. G. Mallet, and A. E. Radunskaya, “Spatial tumor-immune modeling,” Computational and Mathematical Methods in Medicine, vol. 7, no. 2-3, pp. 159–176, 2006.
View at: Publisher Site | Google Scholar
M. Villasana and A. Radunskaya, “A delay differential equation model for tumor growth,” Journal of Mathematical Biology, vol. 47, no. 3, pp. 270–294, 2003.
View at: Publisher Site | Google Scholar
S. Xu and M. Bai, “Stability of solutions to a mathematical model for necrotic tumor growth with time delays in proliferation,” Journal of Mathematical Analysis and Applications, vol. 421, no. 1, pp. 955–962, 2015.
View at: Publisher Site | Google Scholar
S. Ghosh and G. P. Samanta, “Model justification and stratification for confounding of Chlamydia trachomatis disease,” Letters in Biomathematics, pp. 1–13, 2019.
View at: Publisher Site | Google Scholar
R. Eftimie, J. L. Bramson, and D. J. D. Earn, “Interactions between the immune system and cancer: a brief review of non-spatial mathematical models,” Bulletin of Mathematical Biology, vol. 73, no. 1, pp. 2–32, 2011.
View at: Publisher Site | Google Scholar
P. M. Altrock, L. L. Liu, and F. Michor, “The mathematics of cancer: integrating quantitative models,” Nature Reviews Cancer, vol. 15, no. 12, pp. 730–745, 2015.
View at: Publisher Site | Google Scholar
R. Eftimie, J. J. Gillard, and D. A. Cantrell, “Mathematical models for immunology: current state of the art and future research directions,” Bulletin of Mathematical Biology, vol. 78, no. 10, pp. 2091–2134, 2016.
View at: Publisher Site | Google Scholar
K. P. Wilkie, “A review of mathematical models of cancer-immune interactions in the context of tumor dormancy,” in Systems Biology of Tumor Dormancy, H. Enderling, N. Almog, and L. Hlatky, Eds., vol. 734, Springer, New York, NY, USA, 2013.
View at: Google Scholar
M. Robertson-Tessi, A. El-Kareh, and A. Goriely, “A mathematical model of tumor-immune interactions,” Journal of Theoretical Biology, vol. 294, pp. 56–73, 2011.
View at: Google Scholar
K. J. Mahasa, R. Ouifki, A. Eladdadi, and L. d. Pillis, “Mathematical model of tumor-immune surveillance,” Journal of Theoretical Biology, vol. 404, pp. 312–330, 2016.
View at: Publisher Site | Google Scholar
G. P. Samanta, R. G. Aíza, and S. Sharma, “Analysis of a mathematical model of periodically pulsed chemotherapy treatment,” International Journal of Dynamics and Control, vol. 5, no. 3, pp. 842–857, 2017.
View at: Publisher Site | Google Scholar
S. Sharma and G. P. Samanta, “Analysis of the dynamics of a tumor-immune system with chemotherapy and immunotherapy and quadratic optimal control,” Differential Equations and Dynamical Systems, vol. 24, no. 2, pp. 149–171, 2016.
View at: Publisher Site | Google Scholar
L. G. de Pillis and A. E. Radunskaya, “Modeling tumor-immune dynamics,” in Mathematical Models of Tumor-Immune System Dynamics, Springer, New York, NY, USA, 2014.
View at: Google Scholar
L. dePillis, H. Savage, and A. E. Radunskaya, “Mathematical model of colorectal cancer with monoclonal antibody treatments,” British Journal of Medicine and Medical Research, vol. 4, no. 16, pp. 3101–3131, 2014.
View at: Publisher Site | Google Scholar
L. G. de Pillis, W. Gu, and A. E. Radunskaya, “Mixed immunotherapy and chemotherapy of tumors: modeling, applications and biological interpretations,” Journal of Theoretical Biology, vol. 238, no. 4, pp. 841–862, 2006.
View at: Publisher Site | Google Scholar
L. de Pillis, K. R. Fister, W. Gu et al., “Mathematical model creation for cancer chemo-immunotherapy,” Computational and Mathematical Methods in Medicine, vol. 10, no. 3, pp. 165–184, 2009.
View at: Publisher Site | Google Scholar
R. Eftimie, J. L. Bramson, and D. J. D. Earn, “Modeling anti-tumor Th1 and Th2 immunity in the rejection of melanoma,” Journal of Theoretical Biology, vol. 265, no. 3, pp. 467–480, 2010.
View at: Publisher Site | Google Scholar
L. Anderson, S. Jang, and J.-L. Yu, “Qualitative behavior of systems of tumor-CD4+-cytokine interactions with treatments,” Mathematical Methods in the Applied Sciences, vol. 38, no. 17, pp. 4330–4344, 2015.
View at: Publisher Site | Google Scholar
X. Hu and S. R.-J. Jang, “Dynamics of tumor-CD4+-cytokine-host cells interactions with treatments,” Applied Mathematics and Computation, vol. 321, pp. 700–720, 2018.
View at: Publisher Site | Google Scholar
J. Mattes, M. Hulett, W. Xie et al., “Immunotherapy of cytotoxic T cell-resistant tumors by T helper 2 cells,” The Journal of Experimental Medicine, vol. 197, no. 3, pp. 387–393, 2003.
View at: Publisher Site | Google Scholar
A. Perez-Diez, N. T. Joncker, K. Choi et al., “CD4 cells can be more efficient at tumor rejection than CD8 cells,” Blood, vol. 109, no. 12, pp. 5346–5354, 2007.
View at: Publisher Site | Google Scholar
S. Zhang, D. Bernard, W. I. Khan, M. H. Kaplan, J. L. Bramson, and Y. Wan, “CD4+ T-cell-mediated anti-tumor immunity can be uncoupled from autoimmunity via the STAT4/STAT6 signaling axis,” European Journal of Immunology, vol. 39, no. 5, pp. 1252–1259, 2009.
View at: Publisher Site | Google Scholar
A. Ghaffari, M. Nazari, and F. Arab, “Suboptimal mixed vaccine and chemotherapy in finite duration cancer treatment: state-dependent riccati equation control,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 37, no. 1, pp. 45–56, 2015.
View at: Publisher Site | Google Scholar
A. Ghaffari, B. Bahmaie, and M. Nazari, “A mixed radiotherapy and chemotherapy model for treatment of cancer with metastasis,” Mathematical Methods in the Applied Sciences, vol. 39, no. 15, pp. 4603–4617, 2016.
View at: Publisher Site | Google Scholar
M. E. Dudley, J. R. Wunderlich, P. F. Robbins et al., “Cancer regression and autoimmunity in patients after clonal repopulation with antitumor lymphocytes,” Science, vol. 298, no. 5594, pp. 850–854, 2002.
View at: Publisher Site | Google Scholar
J. C. Arciero, T. L. Jackson, and D. E. Kirschner, “A mathematical model of tumor-immune evasion and siRNA treatment,” Discrete and Continuous Dynamical Systems—Series B, vol. 4, no. 1, pp. 39–58, 2004.
View at: Google Scholar
D. Kirschner and J. C. Panetta, “Modeling immunotherapy of the tumor - immune interaction,” Journal of Mathematical Biology, vol. 37, no. 3, pp. 235–252, 1998.
View at: Publisher Site | Google Scholar
M. A. Cooper, T. A. Fehniger, and M. A. Caligiuri, “The biology of human natural killer-cell subsets,” Trends in Immunology, vol. 22, no. 11, pp. 633–640, 2001.
View at: Publisher Site | Google Scholar
R. Glas, L. Franksson, C. Une et al., “Recruitment and activation of natural killer (nk) cells in vivo determined by the target cell phenotype,” The Journal of Experimental Medicine, vol. 191, no. 1, pp. 129–138, 2000.
View at: Publisher Site | Google Scholar
V. A. Kuznetsov, I. A. Makalkin, M. A. Taylor, and A. S. Perelson, “Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis,” Bulletin of Mathematical Biology, vol. 56, no. 2, pp. 295–321, 1994.
View at: Publisher Site | Google Scholar
A. Huang, P. Golumbek, M. Ahmadzadeh, E. Jaffee, D. Pardoll, and H. Levitsky, “Role of bone marrow-derived cells in presenting MHC class I-restricted tumor antigens,” Science, vol. 264, no. 5161, pp. 961–965, 1994.
View at: Publisher Site | Google Scholar
S. M. Gilbertson, P. D. Shah, and D. A. Rowley, “NK cells suppress the generation of Lyt-2+ cytolytic T cells by suppressing or eliminating dendritic cells,” The Journal of Immunology, vol. 136, no. 10, pp. 3567–3571, 1986.
View at: Google Scholar
A. V. Gett, F. Sallusto, A. Lanzavecchia, and J. Geginat, “T cell fitness determined by signal strength,” Nature Immunology, vol. 4, no. 4, pp. 355–360, 2003.
View at: Publisher Site | Google Scholar
P. Unni and P. Seshaiyer, “Mathematical modeling, analysis, and simulation of tumor dynamics with drug interventions,” Computational and Mathematical Methods in Medicine, vol. 2019, Article ID 4079298, 13 pages, 2019.
View at: Publisher Site | Google Scholar
D. F. Griffiths and D. J. Higham, Numerical Methods for Ordinary Differential Equations: Initial Value Problems, Springer, Basel, Switzerland, 2010.
M. H. Holmes, Introduction to Numerical Methods in Differential Equations, Springer, Basel, Switzerland, 2010.
E. Hairer and G. Wanner, Solving Ordinary Differential Equations II: Stiff and Differential-A Lgebraic Problems, Springer, Basel, Switzerland, 2nd edition, 1996.
J. C. Butcher, Numerical Methods for Ordinary Differential Equations, Wiley, Hoboken, NJ, USA, 3rd edition, 2016.
I. N. Katz, M. A. Franklin, and A. Sen, “Optimally stable parallel predictors for adams-moulton correctors,” Computers & Mathematics with Applications, vol. 3, no. 3, pp. 217–233, 1977.
View at: Publisher Site | Google Scholar
J. Peinado, J. Ibáñez, E. Arias, and V. Hernández, “Adams-bashforth and adams-moulton methods for solving differential riccati equations,” Computers & Mathematics with Applications, vol. 60, no. 11, pp. 3032–3045, 2010.
View at: Publisher Site | Google Scholar
M. L. Ghrist, B. Fornberg, and J. A. Reeger, “Stability ordinates of Adams predictor-corrector methods,” BIT Numerical Mathematics, vol. 55, no. 3, pp. 733–750, 2015.
View at: Publisher Site | Google Scholar
S. Sastry, Nonlinear Systems: Analysis, Stability, and Control, Springer-Verlag, New York, NY, USA, 1999.
H. K. Khalil, Nonlinear Systems, Prentice Hall, Upper Saddle River, NJ, USA, 3rd edition, 2002.
D. W. Jordan and P. Smith, Nonlinear Ordinary Differential Equations, Oxford University Press, New York, NY, USA, 4th edition, 2007.
X. Hou and J. Wu, “Stability analysis of a multi-Layer tumor model with free boundary,” Computers and Mathematics with Applications, vol. 77, no. 1, pp. 199–208, 2019.
View at: Publisher Site | Google Scholar
S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry and Engineering, Perseus Books Publishing, USA, 1994.

Copyright

Copyright © 2020 Ahmed M. Makhlouf 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.

PDF Download Citation

Download other formats

Order printed copies

Views

2797

Downloads

2062

Citations

Computational and Mathematical Methods in Medicine

Mathematical Modelling for the Role of CD4+T Cells in Tumor-Immune Interactions

Abstract

1. Introduction

2. Materials and Methods

2.1. Mathematical Model

2.2. Method of Solution

3. Results and Discussion

3.1. First: A Case of No Treatment

3.2. Second: A Case of Continuous Treatment

3.3. Third: A Case of Pulsed Treatment

3.4. Stability Analysis

3.4.1. Equilibrium and Stability for Zero Tumor

3.4.2. Equilibrium and Stability for Nonzero Tumor

3.4.3. Results of Stability Analysis for Both the Zero Tumor and Nonzero Tumor

3.5. Bifurcation Analysis

4. Conclusions

Data Availability

Disclosure

Conflicts of Interest

References

Copyright

Mathematical Modelling for the Role of CD4⁺T Cells in Tumor-Immune Interactions