Research Article  Open Access
Chunmei Zeng, Shaojuan Ma, "Dynamic Analysis of a TumorImmune System under Allee Effect", Mathematical Problems in Engineering, vol. 2020, Article ID 4892938, 11 pages, 2020. https://doi.org/10.1155/2020/4892938
Dynamic Analysis of a TumorImmune System under Allee Effect
Abstract
In this paper, we develop a definite tumorimmune model considering Allee effect. The deterministic model is studied qualitatively by mathematical analysis method, including the positivity, boundness, and local stability of the solution. In addition, we explore the effect of random factors on the transition of the tumorimmune system from a stable coexistence equilibrium point to a stable tumorfree equilibrium point. Based on the method of stochastic averaging, we obtain the expressions of the steadystate probability density and the mean firstpassage time. And we find that the Allee effect has the greatest impact on the number of cells in the system when the Allee threshold value is within a certain range; the intensity of random factors could affect the likelihood of the system crossing from the coexistence equilibrium to the tumorfree equilibrium.
1. Introduction
With the continuous improvement of living standards, individuals pay more attention to their own health. As one of the most harmful diseases to human health and life, malignant tumor has aroused wide concern in the world. At the same time, it has also become a research hotspot in academia. Thompson [1] discussed strategies to target YAP/TAZ activity in cancer and suggested that YAP and TAZ are important drivers of solid tumor growth, metastasis, and resistance to therapy. Du et al. [2] represented a firsttime report of NSCLCintrinsic PD1 expression and a potential mechanism by which PD1 blockade may promote cancer growth. Studies on ascorbic acid, which chemosensitize colorectal cancer cells and synergistically inhibit tumor growth, were presented in [3]. Chaoul et al. [4] found that the lack of MHC class II molecules was beneficial to the infiltration of Tcells into tumors, thus enhancing the control of tumor growth. Villalobos et al. [5] demonstrated the anticancer effect of SOCS1 in prostate cancer and identified its potential inhibitory mechanism. Volz et al. [6] studied a new family of potent TLR9 agonists, EnanDIM, and it can inhibit tumor growth.
Although there are some achievements on finding factors which can effectively inhibit or control tumor growth. What we should not ignore is that the human body itself has a mechanism that can resist tumor cells—the immune system. It can not only protect the body from various internal and external factors against the normal immune system but also defend the body from tumor interference so that the normal operation of the body tissues. This function of the immune system has urged mathematicians and biologists exert to the research on its dynamic properties. Adriana et al. [7] discussed the regulation of angiogenesis by innate immune cells in the tumor microenvironment, specific features, and roles of major players. Key steps for visualizing the complex interaction between malignant tumors and the immune system were proposed and discussed in [8]. Jin et al. [9] summarized the current understanding of autophagy and its regulation in cancer. Medzhitov and Janeway [10] introduced the role of innate immune mechanism in the development of organisms and illustrated that it is the first line of defense of human immunity. The dynamics between tumor and immune cells in microenvironment can be referred to [11–14]. Khajanchi [15] described how tumor cells evolve and survive after encounter with the immune system. Terry et al. [16] presented new insights on the role of EMT in tumorimmune escape.
At present, the research on tumorimmune system has obtained above results, and Katrin et al. [17] showed that tumors have the Allee effect. According to the study of literature [18], the Allee effect is a correlation between any component of individual fitness and the number or density of the same individual. Since the Allee effect has been shown to change the optimal control decisions and control costs and estimate the potential risks in ecology, Katrin et al. [17] believed that it was also crucial for tumor growth control. However, so far, it has been neglected in tumor growth and persistence. Therefore, under the excitation of the above research, we attempt to add the Allee effect on the basis of Pillis [19] model, study its dynamic behavior, and explore the influence of Allee effect on tumor cells.
The remainder of the paper is organized as follows. In Section 2, we develop a definite tumorimmune model with the Allee effect and analyze the property of the deterministic model, including the positivity, boundness, and local stability. Section 3 is devoted to discuss the tumorimmune model with the Allee effect, which is disturbed by white noise. We carry out numerical simulations in Section 4. In Section 5, we take a discussion and draw a conclusion.
2. The Deterministic TumorImmune Model with Allee Effect
2.1. Establishment of Deterministic Model
In 2003, Pillis and Radunskaya [20, 21] proposed a mathematical model to describe the relationship between the tumor and immune system:
Yet, for the tumor itself, it has its own optimal density in the process of growth and reproduction, whether it is too dense or too thin will have an inhibitory effect on its growth. Based on this, we study the relationship between tumor cells, normal cells, and immune cells under the background of the Allee effect. The mathematical model is as follows:where represents the number of immune cells, is the amount of tumor cells, can be thought of the quantity of normal cells, and represent the resource constant and mortality of effector immune cells, respectively, and is Allee threshold and satisfies the condition of . Since the presence of tumor cells would activate effector immune cells, this term is expressed as by Michaelis–Menten expression. Competition between effector immune cells and tumor cells would lead to the death of two types of cells, which are expressed as and . The maximum growth rate and environmental carrying capacity of normal cells are represented by and , respectively. Similarly, the tumor cells are expressed by and . stand for the cell death rate caused by the competition between normal cells and tumor cells, respectively.
In this study, we assume that the tumor is homogeneous. And according to the physical meaning of each variable, the initial conditions of model (2) are
2.2. Positiveness and Boundedness of Solutions
Theorem 1. For any greater than 0, the solution of system (2) is positive and bounded under the initial conditions (3), when .
Proof. (1)Positiveness: integrating both sides of formula (2), we have the following relations: The initial condition is . When , it can be obtained from the above equation:(2)Boundedness: considering a function when according to equation (2), we can obtain that Then, the above equation can be written as follows: where and , where is the maximum value obtained by within the range of and is a positive number selected which is greater than the difference between and . We then have When , . This shows that, under initial conditions (3), the solution of system (2) is bounded.
2.3. Classification of Equilibrium Points
According to equation (2), we get the following three kinds of equilibrium points:(1)Tumorfree equilibrium point: .(2)Coexisting equilibrium point: , where satisfies the equation and the coefficients of equation (9) are(3)Abnormal equilibrium point: , where is determined by the equation and the coefficients of above equation are
2.4. Stability of Equilibrium Points
In this section, we analyze the local stability of the equilibrium points. The calculation shows that the number of normal cells at the abnormal equilibrium point is zero, which is not consistent with the practical significance of biology, so we only discuss the local asymptotic stability of tumorfree equilibrium point and coexistence equilibrium point.
The Jacobian matrix of system (2) at the equilibrium point is as follows:
2.4.1. Stability of TumorFree Equilibrium Point
The eigenvalues of Jacobian matrix at are
And by calculation, the equilibrium solution is locally asymptotically stable if and only if , while is the saddle point with index two when .
2.4.2. Stability of Coexistence Equilibrium Point
When the equilibrium point is , the characteristic equation of matrix (13) at this point iswhere satisfies equation (9), and
According to Routh–Hurwitz criterion, we know that
The determinant of Hurwitz can be calculated according to the above expressions:
Thus, when are satisfied at the same time, coexistence balance is locally asymptotically stable.
By equation (9), we know that coexistence equilibrium expression is the root of quartic equation, namely, there are at most four coexistence equilibrium points. Considering the meaningful number of equilibrium and its stability is related to parameter selection, and in order to make the results more intuitive, we carry out numerical simulation according to effective range of variables. As a result, the number and stability of the coexistence equilibrium point are shown in Figure 1. In this numerical simulation, the horizontal axis is value from 0 to 2, while the vertical axis is value from −2 to 1.
In order to reveal the role of the Allee effect in the system, we only change the value of . By referring to [20, 21], we make ; then, the time history diagrams of three cells in different are as follows.
Observing Figure 2, we know when , the number of the three types of cells do not change significantly. When belongs to , the cell change trend is the most significant. This indicates that when the Allee threshold of tumor cells is within the range of −1.2 to −0.5, the Allee effect has the greatest influence on the number of cells in the system, and this also proves that the Allee effect has an impact on the tumorimmune system.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
3. Stochastic TumorImmune Model with Allee Effect
3.1. Establishment of Stochastic Model
Considering that the growth of tumor cells can be affected by various uncertainties, such as vibration and temperature, and these can be simulated with noise, therefore, based on the deterministic model (2), we propose the following stochastic model:
Among them, , , and are three random processes, which represent the number of immune cells, normal cells, and tumor cells, respectively. is the independent Gaussian white noise with zero mean intensity , which acts on the natural growth rate of normal cells and the nature growth rate of tumor , and we assume that both damping and excitation are small. The specific meanings of other parameters are the same as those in model (2).
3.2. The Stochastic Averaging Method
In order to facilitate the application of the stochastic average method, we first simplified model (19), setting
Then, model (19) can be written as follows:
The above equation (21) is converted into the following stochastic differential equation by formula:where and are independent Wiener process and and are the Wong–Zakai amendments. Defining a random process,
Among them, . , , and are positive, . Then, the stochastic differential equation of random process is
Using the random average principle of StratorovichKhasminskii, the average stochastic differential equation of is
Whence,and according to the formula,
We have the following expressions:where
Then, with the help of equations (26)–(30), the drift coefficient and diffusion coefficient can be simplified as follows:
3.2.1. Stationary Probability Density
Base on equation (23), we obtain that
Substitute (33) and (34) into (23); then, it can be written as follows:
Due to equations,thus, we obtain that
When the boundary is nonsingular; when it is the second kind of singular boundary. The diffusion index the drift index . Therefore, the nontrivial stationary probability density exists. And its steadystate probability density function satisfies the following FPK equation:
Whence,
Solving the FPK equation, we then havewhere is a normalization constant, and
Hence, equation (40) can be simplified as follows:
Suppose , the conditional probability density of random process is ; then, the joint probability density of and can be expressed as follows:
As a result, we obtain the stationary probability density of and :
Thus, we get the steadystate probability densities of and , respectively. They are
3.2.2. Mean FirstPassage Time
When people are suffering from tumor, we hope that the number of tumor cells go to zero. However, according to the calculation in the above section, we know that the steadystate probability density of tumor cells is present, which indicates that tumor cells will not be extinct under natural state, but the effects of random factors may change this. Therefore, it is practical to study the average firstpassage time from the coexistence equilibrium point to the tumorfree equilibrium point, and the Pontryagin equation is satisfied thatwhere
From the above studies, we know that is nonsingular and is natural. According to the constant variation method, can be written as follows:where is the integral constant, and
4. Numerical Simulation
Refer to the experimental data of article [20–29] (Table 1) and the analysis of stability of equilibria in the previous section, and we obtain that system (19) of the tumorfree equilibrium is stable, and the corresponding critical state is denoted by . Calculate based on Table 1 and reference to the Figure 1, we know the result of coexisting equilibrium belongs to the red area when the value is −0.8, −0.9, and −1.0. Namely, system (19) has only one stable coexistence equilibrium, and the corresponding critical state is denoted by . As a result, is nonsingular, is natural, and the expression of the firstpassage time from to satisfies equation (46). In order to investigate the effect of the Allee threshold and noise intensity on the firstpassage time, we numerically simulate the above results and obtain Figures 3 and 4.

(a)
(b)
(c)
(a)
(b)
(c)
As can be seen from Figures 3 and 4, the time of first passage is limited only when the random intensity is within a certain range. This indicates that the existence of random factors does not cause the firstpassage to occur, and only when the random intensity is within a certain range will the system cross from the coexistence equilibrium to the tumorfree equilibrium.
5. Conclusions
In this paper, we analyze the property of the deterministic tumorimmune model under the Allee effect, including the positivity and boundness. The results show that, in the condition of , the solution of system (2) is positive and bounded. At the same time, we know that the Allee effect affects the growth and reproduction of cells in the tumorimmune system according to Figure 2. And the influence is most significant when the Allee threshold of tumor cells is in the range of −1.2 to −0.5.
And in the meantime, we are devoted to explore the effect of random factors on the transition of the tumorimmune system from the stable coexistence equilibrium to the stable tumorfree equilibrium. Based on the method of stochastic averaging, we obtain the expressions of the steadystate probability density and the mean firstpassage time. Through the numerical simulation about the expression of firstpassage time, we know that the existence of random factors will not make the firstpassage must happen, in other words, the random strength only within a certain range; the system will pass from the stable coexistence equilibrium point to the tumorfree equilibrium point. Our conclusions provide a possible method to change the number of tumor cells, which is expected to be helpful for the treatment of tumors.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the National Natural Science Foundation (11772002), Ningxia Higher Education Firstclass Discipline Construction Funding Project (NXYLXK2017B09), Major Special Project of North Minzu University (ZDZX201902), Open Project of the Key Laboratory of Intelligent Information and Big Data Processing of NingXia Province (2019KLBD008), “Light of the West” Talent Training Program of Chinese Academy of Sciences (XAB2016AW04), and Graduate Innovation Project of North Minzu University (YCX19130).
References
 B. J. Thompson, “YAP/TAZ: drivers of tumor growth, metastasis, and resistance to therapy,” BioEssays, vol. 42, no. 5, Article ID 1900162, 2020. View at: Publisher Site  Google Scholar
 S. Du, N. Mccall, K. Park et al., “Blockade of tumorexpressed PD1 promotes lung cancer growth,” OncoImmunology, vol. 7, no. 4, Article ID e1408747, 2018. View at: Publisher Site  Google Scholar
 A. S. Pires, R. C. Marque, A. M. Abrantes et al., “Ascorbic acid chemosensitizes colorectal cancer cells and synergistically inhibits tumor growth,” Frontiers in Physiology, vol. 9, p. 911, 2018. View at: Publisher Site  Google Scholar
 N. Chaoul, A. Tang, B. Desrues et al., “Lack of MHC class II molecules favors CD8+ Tcell infiltration into tumors associated with an increased control of tumor growth,” Oncoimmunology, vol. 7, no. 3, Article ID e1404213, 2017. View at: Publisher Site  Google Scholar
 H. A. Villalobos, D. Bobbala, R. Kandhi et al., “SOCS1 inhibits migration and invasion of prostate cancer cells, attenuates tumor growth and modulates the tumor stroma,” Prostate Cancer and Prostatic Diseases, vol. 20, no. 1, pp. 36–47, 2017. View at: Publisher Site  Google Scholar
 B. Volz, K. Kapp, D. Oswald et al., “A new family of potent TLR9 agonists, EnanDIM, which inhibit tumor growth in various syngeneic murine models,” European Journal of Cancer, vol. 92, no. 1, pp. S5–S6, 2018. View at: Publisher Site  Google Scholar
 A. Adriana, B. Antonino, D. M. Noonan et al., “Contribution to tumor angiogenesis from innate immune cells within the tumor microenvironment: implications for immunotherapy,” Frontiers in Immunology, vol. 9, p. 527, 2018. View at: Publisher Site  Google Scholar
 M. Preusser, A. S. Berghoff, C. Thallinger et al., “Cancer immune cycle: a video introduction to the interaction between cancer and the immune system,” ESMO Open, vol. 1, no. 3, Article ID e000056, 2016. View at: Publisher Site  Google Scholar
 Y. Jin, Y. Hong, and C. Hong, “Molecular interactions of autophagy with the immune system and cancer,” International Journal of Molecular Sciences, vol. 18, no. 8, pp. 1694–1703, 2017. View at: Publisher Site  Google Scholar
 R. Medzhitov and C. A. Janeway, “Innate immunity: the virtues of a nonclonal system of recognition,” Cell, vol. 91, no. 3, pp. 295–298, 1997. View at: Publisher Site  Google Scholar
 J. Wang, D. Li, H. Cang et al., “Crosstalk between cancer and immune cells: role of tumorassociated macrophages in the tumor microenvironment,” Cancer Medicine, vol. 8, no. 10, pp. 1–13, 2019. View at: Publisher Site  Google Scholar
 H. Simon, H. L. Christian, and K. Simone, “MiRNAs: dynamic regulators of immune cell functions in inflammation and cancer,” Cancer Letters, vol. 431, pp. 11–21, 2018. View at: Publisher Site  Google Scholar
 E. Naghmehossadat, F. L. Tamara, and H. K. Phillip, “Dynamic adaptation of tumor immune response with nivolumab demonstrated by 18FFDG PET/CT,” Clinical Nuclear Medicine, vol. 43, no. 2, pp. 114–116, 2018. View at: Publisher Site  Google Scholar
 L. Pang, S. Liu, X. Zhang et al., “Mathematical modeling and dynamic analysis of antitumor immune response,” Journal of Applied Mathematics and Computing, vol. 62, no. 12, pp. 473–488, 2020. View at: Publisher Site  Google Scholar
 S. Khajanchi, “Chaotic dynamics of a delayed tumorimmune interaction model,” International Journal of Biomathematics, vol. 13, no. 2, pp. 1039–1091, 2020. View at: Publisher Site  Google Scholar
 S. Terry, P. Savagner, S. OrtizCuaran et al., “New insights into the role of EMT in tumor immune escape,” Molecular Oncology, vol. 11, no. 7, pp. 824–846, 2017. View at: Publisher Site  Google Scholar
 B. Katrin, H. Hatzikirou, V. B. Anja et al., “An emerging Allee effect is critical for tumor initiation and persistence,” PLoS Computational Biology, vol. 11, no. 9, Article ID e1004366, 2015. View at: Publisher Site  Google Scholar
 A. Konstorum, T. Hillen, and J. Lowengrub, “Feedback regulation in a cancer stem cell model can cause an Allee effect,” Bulletin of Mathematical Biology, vol. 78, no. 4, pp. 754–785, 2016. View at: Publisher Site  Google Scholar
 P. A. Stephens, W. J. Sutherland, and R. P. Freckleton, “What is the Allee effect?” Oikos, vol. 87, no. 1, pp. 185–190, 1999. View at: Publisher Site  Google Scholar
 L. G. D. Pillis and A. Radunskaya, “The dynamics of an optimally controlled tumor model: a case study,” Mathematical and Computer Modelling, vol. 37, no. 11, pp. 1221–1244, 2003. View at: Publisher Site  Google Scholar
 L. G. D. Pillis and A. Radunskaya, “A mathematical tumor model with immune resistance and drug therapy: an optimal control approach,” Journal of Theoretical Medicine, vol. 3, no. 2, pp. 79–100, 2001. 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
 J. C. Panetta, “A mathematical model of periodically pulsed chemotherapy: tumor recurrence and metastasis in a competitive environment,” Bulletin of Mathematical Biology, vol. 58, no. 3, pp. 425–447, 1996. View at: Publisher Site  Google Scholar
 S. Michelson and J. T. Leith, “Host response in tumor growth and progression,” Invasion & Metastasis, vol. 16, no. 45, pp. 235–246, 1996. View at: Google Scholar
 S. Kusama, J. S. Spratt, W. L. Donegan, F. R. Watson, and C. Cunningham, “The gross rates of growth of human mammary carcinoma,” Cancer, vol. 30, no. 2, pp. 594–599, 1972. View at: Publisher Site  Google Scholar
 V. C. Arnerl, S. O. Emdin, B. Lundgren et al., “Breast carcinoma growth rate described by mammographic doubling time and Sphase fraction. Correlations to clinical and histopathologic factors in a screened population,” Cancer, vol. 70, no. 7, pp. 1928–1934, 1992. View at: Google Scholar
 N. Gunduz, B. Fisher, and E. A. Saffer, “Effect of surgical removal on the growth and kinetics of residual tumor,” Cancer Research, vol. 39, no. 10, pp. 3861–3865, 1979. View at: Google Scholar
 J. D. Ford, “Organizational change as shifting conversations,” Journal of Organizational Change Management, vol. 12, no. 6, pp. 480–500, 1999. View at: Publisher Site  Google Scholar
 L. G. D. Pillis, W. Gu, K. R. Fister et al., “Chemotherapy for tumors: an analysis of the dynamics and a study of quadratic and linear optimal controls,” Mathematical Biosciences, vol. 209, no. 1, pp. 292–315, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Chunmei Zeng and Shaojuan Ma. 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.