Cancer dynamics for identical twin brothers
- Ali Ghaffari^{1} and
- Mostafa Khazaee^{2}Email author
https://doi.org/10.1186/1742-4682-9-4
© Ghaffari and Khazaee; licensee BioMed Central Ltd. 2012
Received: 12 May 2011
Accepted: 6 February 2012
Published: 6 February 2012
Abstract
In this paper, a new mathematical model is developed to represent the interaction between healthy and cancer cells in the human body, focusing on the role of environmental factors and quality of life in the progression of cancer. We have investigated the dynamic effects of inputs on cancer growth, and provide an explanation of how cancer has variable behavior patterns throughout the lives of different patients. The behavior of the system with input and its trajectory patterns are investigated using trajectory patterns and stability analysis. The analysis suggests that a proper treatment method should change the dynamics of the cancer instead of only reducing the population of cancer cells and treatment burden.
Keywords
Background
- (i)
To reduce and control the tumor mass such that a specified volume is obtained at the end of the treatment [2].
- (ii)
To lessen the treatment burden of patients. This method considers some constraints on the treatment policy [3].
- (iii)
To evaluate the number of injected cells that affect the equilibrium points of the immune system and thus may ultimately be dangerous [4].
Previous studies have investigated the effects of therapeutic inputs, which are considered to have direct effects on the system states [5]. However, the behavior of cancer changes as the disease progresses [6]. External stresses that represent destructive inputs, such as environmental and quality of life factors, can cause disability in the DNA repair genes [7]. They can also interfere with, and alter, the functions of regulatory growth signals (TGF-α), growth-inhibiting signals (TGF-β), and apoptosis (TP53) [6].
We are interested in analyzing how inputs alter the dynamics of the human body, and whether this is the main factor involved in the occurrence of cancer. This paper deals with a comparison between the responses of human body cells in a set of twin identical brothers who live under two different conditions, or inputs. In the next section, a modified model for the human body cells is presented. In the third section, the equilibrium points of the system and its linear approximation are calculated. The fourth section investigates the stability analysis of the system, using stability theorems and trajectory patterns. The fifth and final section discusses the results and summarizes the conclusions.
Model
We have considered a set of identical twin brothers who were born with the same genetic structures. It is assumed that, after the birth, they reside and grow up in two different locations, and under different environmental conditions. One of them, say A, lives in a polluted environment and under external stresses, but brother B does not.
In Eq. 1, x and y are the healthy and cancerous cell concentrations, respectively. Because the state variables are physiologically possible, their values are nonnegative, i.e., x ≥ 0 and y ≥ 0. The coefficients a_{1} and a_{2} represent the growth rates of the healthy and cancer cells, respectively. The growth of both healthy tissue and tumor decelerates as the concentrations of both the healthy and tumor tissues approach the carrying capacities K_{1} and K_{2}, respectively [9]. The effect of the immune system is to kill the mutated and cancer cells in proportional rates d_{1} and d_{2}. The immune system agents force the cancer cells to suicide through apoptosis [10]. The coefficient c represents the proportion of healthy cells whose genome has been disordered by external stresses. These cells initiate the neoplastic transformation, and are added to the tumor cells [11].
The tumor competes with healthy tissue for resources, such as blood, nutrients, and space, so the organ "feels" the tumor [12]. Moreover, these same cancer cells compete with each other. The competition coefficients between different cells are b_{1}, b_{2}, and g.
The effects of the input on the dynamics of the system are introduced by the registration of input history in the system coefficients. The destructive inputs, in addition to increasing cell mutation, predispose the cancer cells to intensification of the growth rate, and reduce the death rate by conferring the ability to evade the intracellular controlling mechanism and the immune system.
Cancer progression represents a macro-evolutionary process where karyotype change or genome replacement plays the key dominant role [13]. In the present work, the transformation rate of healthy cells to cancer cells is assumed to be proportional to the input magnitude u. The biotransformation coefficient saturates at a definite limit K_{3}, which is related to the biological limits of body organs and the accumulation of external effects. The other variant parameters have similar formulations, but with different values and rates. When the destructive inputs terminate, their effects remain in the body and may or may not be compensated by the therapeutic inputs or by the body recovery.
Parameter estimations
Parameter | Unit | Estimated value | Source | |||||||
---|---|---|---|---|---|---|---|---|---|---|
Birth | Youth | Middle age | Old | After Old | ||||||
A, B | A | B | A | B | A | B | B | |||
a _{1} | week ^{-1} | 5 | 3.5 | 5 | 4.9 | Estimated | ||||
a _{2} | week ^{-1} | 0.12 | 3.816 | 1.152 | 6 | 1.872 | 6.7 | 1.934 | 3.531 | [4] |
K _{ 1 } | Cell | 80 × 1000 | [4] | |||||||
K _{ 2 } | Cell | 90 × 1000 | Estimated | |||||||
d _{1} | week ^{-1} | 0.1 | Estimated | |||||||
d _{2} | week ^{-1} | 2.5 | 1 | 2.26 | 0.05 | 2.09 | 0.03 | 1.94 | 1.167 | Estimated |
C | week ^{-1} | 0.001 | 0.032 | 0.010 | 0.05 | 0.016 | 0.05 | 0.019 | 0.028 | Estimated |
b _{1} | Cell ^{-1} | 0.06 | 0.191 | 0.063 | 0.3 | 0.107 | 0.6 | 0.115 | 0.206 | [9] |
b _{2} | week ^{-1} | 0.4 | 0.320 | 0.381 | 0.2 | 0.342 | 0.11 | 0.307 | 0.272 | [3] |
G | week ^{-1} | 0.002 | 0.064 | 0.002 | 0.1 | 0.005 | 0.1 | 0.007 | 0.016 | Estimated |
c _{1} | week ^{-1} | 0.001 | [4] | |||||||
K _{3} | week ^{-1} | 0.05 | Estimated |
Some important properties of the model are presented here, as described in the following two statements.
Statement 1. Let $\left[\begin{array}{c}\hfill x\left(t\right)\hfill \\ \hfill y\left(t\right)\hfill \end{array}\right]$be the solution of Eq. 1 Then the nonnegative orthant is invariant, i.e., $\left[\begin{array}{c}\hfill x\left(t\right)\hfill \\ \hfill y\left(t\right)\hfill \end{array}\right]\ge 0$.
- i)
Notice that input u has no direct effect on the variables x and y, it can be concluded that if $\left[\begin{array}{c}\hfill x\left(0\right)\hfill \\ \hfill y\left(0\right)\hfill \end{array}\right]=\left[\begin{array}{c}\hfill 0\hfill \\ \hfill 0\hfill \end{array}\right]$, then $\left[\begin{array}{c}\hfill x\left(t\right)\hfill \\ \hfill y\left(t\right)\hfill \end{array}\right]=\left[\begin{array}{c}\hfill 0\hfill \\ \hfill 0\hfill \end{array}\right]$ for all t ≥ 0.
- ii)
If the solution of Eq.1 approaches the horizontal axis from nonnegative orthant then, ${\left.\u1e8b\right|}_{x=0}=0$, ${\left.\u1e8f\right|}_{x=0}=\left[{a}_{2}\left(1-\frac{y}{{K}_{2}}\right)-\left({d}_{2}+g\right)\right]y\ge \phantom{\rule{0.3em}{0ex}}0$ and y does not decrease, so its value becomes negative.
- iii)
If the solution of Eq.1 approaches the vertical axis from nonnegative orthant then,${\left.\u1e8f\right|}_{y=0}=0$, ${\left.\u1e8b\right|}_{y=0}=\left[{a}_{1}\left(1-\frac{x}{{K}_{1}}\right)-\left({d}_{1}+c\right)\right]x\ge 0$ and x does not decrease, so its value becomes negative.
Hence, there is no solution that exits the first orthant. □
Proof:
Finally, with respect to Eqs. 6 and 8, the exponent convergence region is $\mathit{A}$. □
Equilibrium points
The higher order terms are neglected around the origin (x*, y*) = (0, 0) and the last term C is equal to zero at the equilibrium points.
Stability analysis
We study the stability of the equilibrium points of Eq. 1 in this section. The results of the analysis are stated as follows;
i. If the equilibrium point 3 is located in the first orthant and the equilibrium point 2 is not, then the state variables of Eq.1 will converge to equilibrium point 3 (the healthy state).
ii. If the equilibrium points 2 and 3 are located in the first orthant, then the state variables of Eq.1 will converge to one of these two points.
iii. If the equilibrium point 2 is located in the first orthant and the equilibrium point 3 is not, then the state variables of Eq.1 will converge to the equilibrium point 2 (the cancer state).
Proof:
As shown in Table 1, the values of a_{1} are larger than the sum of d_{1} and c. Also a_{2} is larger than the sum of d_{2} and g. Thus the eigenvalues of equilibrium point 1 are always positive and the origin is an unstable node.
In Eq.10 y is positive, then β_{2}/α_{2} < 0. Also, we notice from Table 1 that if the value of β_{2}/α_{ 2 } is larger than β_{1}, then the equilibrium point 2, if it exists, is a stable node.
3) At the equilibrium point 3, x = α_{6} + α_{7}, y = α(α_{6} + α_{7}) + β_{1}.
This means that the coefficient matrix A is negative definite at this equilibrium point and it is a stable node.
4) At the equilibrium point 4, x = α_{6} - α_{7}, y = α_{1}(α_{6} - α_{6}) + β_{1}.
Therefore, if this equilibrium point exists, then it is an unstable saddle point. □
Numerical analysis
In this section, the trajectory patterns of the dynamic systems for the twin brothers are evaluated. The life spans of both brothers are divided into three stages and the dynamic behaviors of the systems are analyzed before and after the presence of inputs. The first stage is from the birth to 850 weeks, the second stage is from 850 to 1950 weeks, and the third stage is from 1950 weeks to 2900 weeks.
The following figures illustrate the trajectory patterns of the dynamic systems representing twin brothers A and B at the beginning of the three stages of their life.
Figure 2-a shows that for brother A, the equilibrium points 2 and 4 are brought to the positive orthant by receiving the destructive input for a limited time, and the trajectories converge to the healthy equilibrium point 3 from a large portion of the state space. However, the behavior of the system for brother B, as shown in Figure 2-b, is similar to Figure 1, with the positions of the equilibrium points 2 and 4 showing only minor changes.
The body organ of B is changed a little, but the system does not pose any unstable behavior in the first orthant, as is seen in Figure 3-b. The small variation is due to aging and the probable effect of limited destructive inputs.
Finally, a situation is analyzed in which B(the healthy brother) encounters destructive inputs from 2900 to 3200 weeks.
Discussion and conclusion
The stable dynamics for the newborn brothers at birth means that they do not have cancer. If the equilibrium points of the system are in the nonnegative orthant then the system will settle in the equilibrium points 2 or 3. The history of the body organ's dynamics during the life of brother B indicates his healthy condition over his life (Figures 1, 2-b, 3-b, and 4-b). This brother is not afflicted by cancer since he has not faced the destructive inputs. The variation of the system dynamics where the destructive inputs are sufficiently effective, demonstrate cancer formation and its progress in the body of brother A(Figures 1, 2-a, 3-a, and 4-a). The stable behavior for A at birth is gradually converted to less stable conditions. The attraction domain of equilibrium point 2, indicating the cancer cells, becomes larger than equilibrium 3, indicating healthy cells. Thus, the cancer incidence is increased. In such a system, if the number of cancer cells reaches a specified limit, then the patient will die. In our model, brother A reaches the stage in which no treatment is able to cure the patient completely and he dies.
This research shows that focusing on the "indelible changes of the system dynamics" is the best way to describe the cancer formation process. The other major conclusions are as follows:
a) The changes in the dynamics of the system occur gradually in the body of the patient. Thus, the duration of effective destructive inputs is directly related to the probability of the cancer occurrence. This can be seen in Figures 2-a and 3-a.
b) As shown in Figures 2-b, 3-b, and 4-b, the behavior of the system for the healthy brother B, who is not subjected to the destructive inputs, also changes as he becomes older.
c) The effects of environmental and life quality factors are the main causes for the onset of cancer, in bodies where there are no dominant hereditary genetic disorders.
d) The sensitivity of the old healthy brother B to the destructive inputs is greater than that of the young affected brother A, and the system dynamics changes more rapidly for him.
Our next objective is to find a proper therapeutic input that can move equilibrium point 2 out of the nonnegative space. This treatment method would guarantee the impossibility of tumor recurrence. The undesirable changes of the system should be modified by using the corrective inputs, and the treatment not restricted to the system outputs. Furthermore, in future work, the parameters, especially the input, should relate to the qualitative conditions of a patient's life and to the clinical data.
Declarations
Acknowledgements
The authors wish to acknowledge Dr. T. Sadeghi in the Medical Science University of Tehran for her useful comments; and Eng. E.A. Shavazi in the Chemical Engineering Faculty of Newcastle University in Australia for his assistance.
Authors’ Affiliations
References
- Bellomo N, Li NK, Maini PK: On the foundations of cancer modeling: selected topics, speculations, and perspectives. Mathematical Models and Methods in Applied Sciences. 2008, 18 (4): 593-646. 10.1142/S0218202508002796.View ArticleGoogle Scholar
- Ribba B, Colin T, Schnell S: A multiscale mathematical model of cancer, and its use in analyzing irradiation therapies. Theoretical Biology and Medical Modelling. 2006, 3 (7):Google Scholar
- Ghaffari A, Naserifar N: Optimal therapeutic protocols in cancer immunotherapy. Computers in Biology and Medicine. 2010, 40 (3): 261-270. 10.1016/j.compbiomed.2009.12.001.View ArticlePubMedGoogle Scholar
- Fister KR, Panetta JC: Optimal control applied to cell-cycle-specific cancer chemotherapy. Siamese Journal of Applied Mathematics. 2000, 60 (3): 1059-1072. 10.1137/S0036139998338509.View ArticleGoogle Scholar
- De Pillis LG, Fister KR, Gu W, Collins C, Daub M, Gross D, Moore J, Preskill B: Seeking bang-bang solutions of mixed immune-chemotherapy of tumors. Electronic Journal of Differential Equations. 2007, 171: 1-24.Google Scholar
- Kumar V, Cotran RS, Robbins SL: Robbins Basic Pathology. 2003, Philadelphia: W.B. Saunders Co, 7Google Scholar
- Klein CA, Hölzel D: Systemic cancer progression and tumor dormancy. Cell Cycle. 2006, 5 (16): 1788-1798. 10.4161/cc.5.16.3097.View ArticlePubMedGoogle Scholar
- Freedman HI, Pinho STR: Stability criteria for the cure state in a cancer model with radiation. Nonlinear Analysis: Real World Applications. 2009, 10: 2709-2715. 10.1016/j.nonrwa.2008.08.001.View ArticleGoogle Scholar
- Sachs RK, Hlatky LR, Hahnfeldt P: Simple ODE models of tumor growth and anti-angiogenic or radiation treatment. Mathematical and Computer Modeling. 2001, 33: 1297-1305. 10.1016/S0895-7177(00)00316-2.View ArticleGoogle Scholar
- Wei X, Guo C: Global existence for a mathematical model of the immune response to cancer. Nonlinear Analysis: Real World Applications. 2010, doi:10.1016/j.nonrwa.2010.02.017Google Scholar
- Hirata Y, Bruchovsky N, Aihara K: Development of a mathematical model that predicts the outcome of hormone therapy for prostate cancer. Journal of Theoretical Biology. 2010, 264: 517-527. 10.1016/j.jtbi.2010.02.027.View ArticlePubMedGoogle Scholar
- Kuang Y, Nagy JD, Elser JJ: Biological stoichiometry of tumor dynamics: mathematical models and analysis. Dynamics of Continuous, Discrete and Impulsive Systems. 2004, 221-240. B4Google Scholar
- Heng HHQ, Stevens JB, Bremer SW, Ye KJ, Liu G, Ye CJ: The Evolutionary Mechanism of Cancer. Journal of Cellular Biochemistry. 2010, 109 (6): 1072-1084.PubMedGoogle Scholar
- Slotine JE, Li WL: Applied Nonlinear Control. 1991, Englewood Cliffs, New Jersey: Prentice-Hall IncGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.