Application of Brody Growth Function to Describe Dynamics of Breast Cancer Cells
Abdulsamad Engida Sado, Purnachandra Rao Koya
School of Mathematical and Statistical Sciences, Hawassa University, Hawassa, Ethiopia
Abdulsamad Engida Sado, Purnachandra Rao Koya. Application of Brody Growth Function to Describe Dynamics of Breast Cancer Cells. American Journal of Applied Mathematics. Vol. 3, No. 3, 2015, pp. 138-145. doi: 10.11648/j.ajam.20150303.20
Abstract: In this paper we have constructed a mathematical model using Brody function and applied to describe the dynamics of breast cancer. To construct the mathematical model we considered that the linear cancer network technique describes the growth of estrogen receptor positive breast cancers. Model validity is verified using simulation study and mathematical analysis. It is verified that the hormone therapy is a technique to treat endocrine receptor positive breast cancers. Hormone therapy is considered as a treatment and used to block the estrogens receptors from the cancer and health cells. Important observations are made from the simulation study and physical interpretations are drawn and presented lucidly in the paper.
Cancer has been a major cause, among all the diseases, for human death throughout the world. An estimation of 14.10 million cancer cases of all types were found around the world in the year 2012, of those cases 7.40 million cases were found among men and 6.70 million cases in women. This number is expected to increase steadily and reach up to 24 million by the year 2035. Breast cancer cases in women occupy the second most in the list of common cancer types with nearly 1.7 million such new cases found in 2012. Breast cancer cases among women diagnosed in 2012 is about 11.9% of all types of cancer cases found in both men and women except non-melanoma skin cancer. Among women breast cancer alone occupy more than 25% of the cumulative number of all types of the cases identified during2012worldwide.So far the developed countries have been experiencing about 55 % of all the globally identified cancer cases but the cancer occurrence rates have started to rise rapidly among the developing countries due to globalization and other influencing factors [13-14].
In Ethiopia breast cancer is the second most occurring cancer among women which is next to the cervical cancer. It is estimated that around ten thousand Ethiopian women have been affected with breast cancer every year. However, thousands of cases go unreported because the rural women often approach traditional healers for treatment. A very few cases seek treatment from the government health organizations where only these cases are recorded [3,6,9].
The cancer decease is a result of the uncontrolled growth of abnormal cells in human body. Uncontrolled growth of abnormal cells is a result of a set of mutations in the healthy cells. Cells are building blocks of the body and cancer starts from normal cells. Cells in the body born and die naturally. In order to balance natural death of the cells and to maintain cell population equilibrium state, the normal cells in human body divide themselves and grow in number. Cancer occurs when unbounded and fast growth of cells takes place in the body. As a whole, cancer is a collection of more than a hundred diseases or reasons in human body [2,7].
A variety of cancer diseases can occur in human body. In general, name to the cancer disease is given with the name of the body part where the disease occurs viz., Brain cancer, Lung cancer, Bone cancer, blood cancer and Breast cancer are a few in the list. Any party of the human body can be affected by cancer disease [2-3,5,10].Breast cancer is the most common cancer affecting women worldwide. The breast cancer in women can occur in many forms including Endocrine receptor positive for estrogen, Endocrine receptor negative for estrogen, Endocrine receptor positive for progesterone, Endocrine receptor negative for progesterone, HER2 positive and HER2 negative. The most common type of breast cancer in women is Endocrine receptor positive breast cancer. The rise of Estrogen hormone in women human body to the higher levels above the normal level causes an increase in the rate of cell growth and that leads to mutation in healthy cells that resulting in the birth of cancer cells.
The common signs of breast cancer include one or more of the following: (i) a hard compact and unshapely mass develop in the breast or armpit. This lump is painless and it occurs only on one side of the breast, (ii) the breast size increases and / or its shape changes. Indentation (concave shaped cuts on an edge) of breast occurs. Veins (pipes in the body through which blood circulates) grow. Skin erosion (be destroyed gradually) take place, (iii) various changes will appear in the skin. Skin becomes hard. Small hollows (dimpling) will occur in the skin. Swellings (bumps) occur. The color of the skin will change to redness or to orange, and (iv) various changes will appear in the nipple. Nipple goes back into the breast (retraction).Unusual liquids are produced and discharged around the nipple area. Breaking of skin causing red spots (rashes) will occur around the nipple area.
The treatment choices for breast cancer in women include radiation therapy, hormone therapy, targeted therapy, chemotherapy and surgery. Hormone therapy is mainly used to work against Hormone receptive positive breast cancer. Hormone therapy helps to lower the risk of reoccurring of breast cancer after surgery. Hormone therapy is applied to treat breast cancer very frequently. Recall that the high amount of estrogen causes breast cancer. Hormone therapy reduces estrogen levels and helps breast cancer to be under control. Hormone therapy effectively prevents the effect of estrogen on the breast cancer cells [2,5,10].
The technique of anti-estrogens has been the first among all these therapies. This treatment method was very classical and adopted to fight against the breast cancer. Also for a long time the drug containing tamoxifen has been the standard front-line treatment to fight against the breast cancer that occurs after menopausal. Symptoms of menopause or temporary menopause may cause due to the hormonal treatments. Hormone therapy is effectively used in pre-menopausal and post-menopausal women and it works also well in the remission stage of the breast cancer [2,7,10].
In Section 2, we have introduced the Brody growth function as a special case of Koya – Goshu generalized growth function. In Section 3, mathematical modeling to represent the dynamics of breast cancer cells is constructed using Brody growth function. Model variables and parameters are defined. In Section 4, simulation study with various sets of parametric values is taken up and important observations are interpreted. In Section 5, mathematical analysis of the model is considered. Equilibrium points are identified and their positivity and stability are discussed. The paper ends with concluding remarks in Section 6.
2. Brody Model
Brody growth model is the special case of Koya –Goshu generalized biological growth function where the latter is defined as
Here the symbol
Represents a notation and is derived from the four parameters
and . The eight parameters involved in the definition for the function are defined as follows: (i)where is the upper asymptote of , (ii) is the lower asymptote of, (iii), (iv)is the growth rate parameter, (v) is the time shift parameter and is a constant, (vi) is time scale parameter and is a constant, and(vii) and a are the shape parameters. The Koya-Goshu growth model is an eight parameter function. The commonly used growth models such as Brody, Von Bertalanffy, Richards, Weibull, Monomolecular, Mitscherlich, Gompertz, Logistic and Generalized Logistic functions all are special cases of Koya-Goshu growth. The model is a more general solution of the rate-state equation represented by the first ordinary differential equation and the detailed description is given in [1, 4]. Brody function is a special case of Koya – Goshu growth function with the parameters taking the values. and , However, Brody function (Brody, 1945) describing biological growths is defined by
Here in (1) the parameters in their original notations represent: (i) , (ii) is upper asymptote of, (iii) denotes growth at the initial point of time , and (iv) represents the growth parameter [1, 4]. In fact the Brody function given in (1) is a solution of the rate-state equation where the relative rate function is given by . Thus, the rate-state equation for Brody function can be expressed as , see .
3. Brody Function to Represent Breast Cancer Cell Growth
The most common type of breast cancer in women forms when estrogen levels in them are increased intensively. A sequence of reasons contribute for the disease formation and the flow contains the following: Increased levels of estrogen leads to increase in the rate of cellular, latter leads to mutation in health cells which in turn resulting in endocrine receptor positive breast cancer
The idea of Eric Werner [8, 10] can be followed as a framework, in order to construct a mathematical model to represent and describe (i) growth of breast cancer cells and (ii) the effect of Hormone therapy treatment on breast cancer. We use Brody function to construct a set of two differential equations for modeling the growth of cancer stem cells and cancer tumor cells. Eric
It has been proposed a new design to represent cancer growth. This new design suggests that cancer develops as a result of mutations in the normal cells. The network of mutations describes the behavior of cancer stem cells. The stem cells divided and produce tumor cells depending on the instructions of the corresponding network. It is proposed to consider linear cancer networks in which the cancer stem cells can either be divided into cells of the same type or produce tumor cells. Once the model of linear cancer network is constructed, it can be extended to include the effect of estrogen presence and that of the Hormone therapy treatment [8, 10].
The number of cancer stem cells is represented by while the number of cancer tumor cells is represented by respectively at any point of time t. Further, to compute even before hand the effect of hormone therapy on both cancerous and healthy cells, a third differential equation that describes the growth rate of healthy cells within the breast tissue is required to be constructed. The number of healthy cells at any point of time is represented by and we will use Brody function to model too.
In our present study we consider Brody function as where the constant parameter representing the caring capacity or upper asymptotic value or maximum value. The value of the constant can be computed using the formula where the constant represents lower asymptotic value or minimum value. Also the analogues rate-state equation representing Brody function can be expressed as . We now apply the Brody function in its form as described here to construct a system of differential equations that represent the cancer cell and health cell of breast cancer. The system of differential equations representing cancer stem cells, cancer tumor cells and health cells can be constructed using Brody function as follows:
Now, the model (2) is required to be supported by another differential equation that describes the dynamics of estrogens in the human body. Note that it is required to bring estrogen level to natural equilibrium state so as to control the mutation of the normal cells or control formation of new cancer cells and to reduce already formed cancer cells. Hence, reforming or updating of model (2) is required for the following two reasons. Firstly, the new model has to be improved so as to accommodate a new equation representing the dynamics of estrogen in the body. Secondly, in the new model the differential equations representing cancer stem cells, cancer tumor cells and health cells are to be improved so as to include the effect of estrogens on them.
That is, in support to the linear cancer network given in (2) the new model has to be built on the following assumptions made by Eric Werner [8, 10]:
(i) The cancer cells are divided into two types of cells viz., (i) cancer stems cell and (ii) cancer tumor cell
(ii) The computation between cancers stems cells and cancer tumor cells is negligible.
(iii) All the three types of cells viz., cancer stem cells, cancer tumor cells and health cells are Estrogen dependent.
(iv) The tumor cells represented by B are terminal. Terminal means, they do not divide on their own. The cells die with a natural death rate of .
Now using the Brody growth function, a system of four ordinary differential equations can be constructed so as to accommodate the cited assumptions
The term in (3a), (3b) and (3c) accounts for how estrogen helps the growth of stem, tumor and health cells respectively. The parameter represents the carrying capacity or normal amount of Estrogen. Further (i) If , then the cells are dividing with normal rates, (ii) if, then the cell are dividing with slower rates and (iii) If , then the cells are dividing with faster rates. Notations of the variables and parameters used in the model (3) are given in table 1.
|Dying rates of cancer, tumors and health cells respectively due to the lack of sufficient amount of estrogen|
|Amount of estrogen available to bind estrogen receptors in the cell|
|The rate of division of cancer stem cells|
|Normal amount of estrogen needed for all cells to divide normally|
|Carrying capacities of cancer stem cells, tumor cells and health cells respectively|
|Death rate of cancer tumor cell|
|Normal rate of healthy cells near the tumor cells|
|Binding rate of estrogen with cells viz., cancer and health cells|
|Effectiveness (0 to 1) of the Tamoxifen in blocking cells from receiving estrogen|
|Dose of Tamoxifen|
4. Simulation Study
|Parameter||Parametric values used in|
|Figure 1||Figure 2||Figure 3||Figure 4|
We now perform numerical simulations of the model (3) with various parametric values using MATLAB. The parametric values used in these simulations are given in Table 2
4.1. Simulation Without Treatment
We first simulate the model without treatment to visualize normal interactions between cancerous cells, healthy cells, and available estrogen.
Figure 1 shows that before the treatment (i) the number of cancer stem cells represented by blue curve increases with time from zero to the maximum carrying capacity following Brody growth function, (ii) the number of cancer tumor cells represented by red curve increases with time from zero to the maximum carrying capacity following Brody growth function, (iii) the number of health cells is normal (remains constant) and (iv) the amount of estrogen remains constant at the maximum level.
4.2. Simulation with Treatment of Effectiveness
Recall that giving tamoxifen drug to a cancer patient is a kind of hormone therapy and the hormone therapy is a kind of treatment given to a cancer patient. The drug Tamoxifen is used for hormone therapy treatment. We now simulate the model (3) to show the effect of introducing treatment on all the three types the cells and on estrogen.
Figure 2. Brody growth model for breast cancer cells after treatment with effectiveness of .
Figure 2 shows the affect of hormone therapy on all the three types of cells and on estrogen. Tamoxifin is given to meet an effectiveness of . (i) The number of cancer stem cells represented by blue curve decreases with time from maximum carrying capacity 20 to zero following Brody growth function, (ii) the number of cancer tumor cells represented by red curve decreases with time from the maximum carrying capacity 40 to zero following Brody growth function, (iii) the number of health cells represented by green curve decreases with time from the maximum carrying capacity 60 to 55 following Brody growth function, and (iv) the amount of estrogen represented by black curve decreases with time from the maximum carrying capacity 100 to 43 following Brody growth function.
4.3. Simulation with Treatment of Effectiveness
Figure 3 shows that the (i) amount of estrogen, (ii) growth of health cells and (iii) cancer cells coexist in the body of cancer patient when the treatment with an effectiveness of is given. Since the sizes of both types of cancer cells viz., stem and tumor, decreased it may be possible to perform surgical treatment so as to nullyfy the cancer.
Further figure 3 shows the affect of hormone therapy on all the three types of cells and on estrogen. Tamoxifin is given to meet an effectiveness of . (i) The number of cancer stem cells represented by blue curve decreases with time from maximum carrying capacity 20 to 10 following Brody growth function, (ii) the number of cancer tumor cells represented by red curve decreases with time from the maximum carrying capacity 40 to 3 following Brody growth function, (iii) the number of health cells represented by green curve decreases with time from the maximum carrying capacity 60 to 57 following Brody growth function, and (iv) the amount of estrogen represented by black curve decreases with time from the maximum carrying capacity 100 to 72 following Brody growth function.
Figure 3. Brody growth model for breast cancer cells after treatment with effectiveness of.
4.4. Simulation with Treatment of Effectiveness
Figure 4. Brody growth model for breast cancer cells after treatment with effectiveness of.
Further figure 4 shows the affect of hormone therapy on all the three types of cells and on estrogen. Tamoxifin is given to meet an effectiveness of . (i) The number of cancer stem cells represented by blue curve decreases with time from maximum carrying capacity 20 to zero following Brody growth function, (ii) the number of cancer tumor cells represented by red curve decreases with time from the maximum carrying capacity 40 to zero following Brody growth function, (iii) the number of health cells represented by green curve decreases with time from the maximum carrying capacity 60 to 45 following Brody growth function, and (iv) the amount of estrogen represented by black curve decreases with time from the maximum carrying capacity 100 to 72 following Brody growth function.
Further in figure 4 it can be observed that the sizes of both cancer cells viz., tumor and stem, fall off to zero. However the fall of tumor cells is much faster than that of stem cells.
The observations made from the simulation study are tabulated below:
|Without treatment (Figure 1)||Estrogen is Maximum Cancer cells grow to maximum Health cells remain constant||Patient is suffering with breast cancer and treatment is required.|
|Estrogen decreases to normal level Cancer cells die out Health cells decrease||The treatment is effective and the patient is free of breast cancer disease.|
|Health cells, cancer cells and estrogen exist together||The treatment with harmone therpy is not sufficient. Surgical treatment is required.|
|Tumor cells die Stem cells decrease||The treatment with harmone therpy is not sufficient. Some additional or adjuvant treatment is required.|
5. Analysis of the Breast Cancer Model
In this section we perform analysis of the breast cancer model designed using Brody function. The analysis includes (i) dimensionless form of the model, (ii) identifying equilibrium points, (iii) stability analysis of the equilibrium points, and (iv) positivity of equilibrium points. This analysis is required and of helpful to manipulate the parametric values in simulation study so as to obtain biologically meaning full results or conclusions.
5.1. Dimension Less Form
Dimensionless form of a mathematical model helps (i) to decrease the number of parameters involved in the model and (ii) to make the variables dimensionless, that is, the units of the variables are removed. Note that the minimization of the number of parameters and eliminating the units of variables will make the mathematical analysis simple. Hence, we construct the dimensionless form the model (3) in this section as follows:
We now introduce four dimensionless variables and defined by,, and respectively. Also we introduce nine parameters defined by ,, , , , and respectively. In view of introduction of these new variables and parameters the dimensionless form of the mathematical model (3) takes the form as follows:
In (4), the over bar of the variables on the left hand side of each equation represent a differentiation with respect to.
From the dimension less form of Brody model given in (4) we now identify the equilibrium points and undertake the stability analysis of these equilibrium points as follows:
5.2. Equilibrium Points
A point in real line is said to be an equilibrium point of one dimensional ordinary differential equation if the point satisfies the condition, Equilibrium point can also be termed as stationary point or also as fixed point. Now in order to obtain the stationary point of the model (4) let us set the right hand sides of each of the four equations in (4) simultaneously equal to 0 and solve, so as to get the following:
The dimensionless model (4) has the equilibrium point where the expressions for and are described in (5).We now take up the analysis of the stability of this equilibrium point by (i) constructing the Jacobian matrix and (ii) finding the eigenvalues of the Jacobian matrix and (iii) finding the signs of these Eigen values at the equilibrium point.
5.3. Positivity of Equilibrium Point
In order to show that the model (4) is well defined and is biological meaning full, we restrict that all of the four population sizes can approach a non-negative steady state. Hence, we establish the condition on the parameters as . These non-negativity conditions are required to be satisfied for the equilibrium points. To guarantee the for all value of the parameters , ,,,,, found in the dimension less form is positive hence the mathematical model of Brody growth model for breast cancer is well defined both mathematical and biological.
5.4. Stability Analysis of the Equilibrium Points
The Jacobian matrix is used to analyze the local stability of equilibrium points where the Jacobian matrix for the model (4) is given by
and the Jacobian matrix at the equilibrium point
takes the form
Where we have used , and to represent, respectively, the following:
The four Eigen values of the Jacobian matrices evaluated at the equilibrium point are given by
Based on the algebraic sign of the Eigen values of the Jacobian matrix, we can determine the nature of the equilibrium points. Since all the Eigen values are less than zero (negative) the dimensionless system(4) described using Brody growth function is stable at the equilibrium point.
Mathematical modeling has been an ever increasing tool to explore complex systems and unsolved problems in biology so as to make the study of their behavior simple. In this paper, we have applied Brody growth function on the mathematical modeling of growth of breast cancer cells and response of hormone therapy on cancer disease. The Brody growth function is a special case of Koya-Goshu generalized model and the former is a better model to address the breast cancer cell growth. In this study we have considered the linear network growth of cancer cells and also the growth of all cells depend on estrogens.
The maximum levels of estrogen in the human body results in (i) making the abnormal growth of health cells by protecting them from the natural death and (ii) occurrence of estrogen positive respiratory breast cancer and its rapid growth (see Figure 1). When hormone therapy is given to a breast cancer patient it is observed that (i) the estrogen level is reduced rapidely to the normal size and (ii) the tumor and stem cells die together. Hence, the hormone therapy tretements is effectiveon breast cancer (see Figure 2). Even though the health cells, cancer stem and tumor cells and estrogen present in the body of breast cancer patient, if the cancer cell size is reduced to a minimum level surgery can be performed to cure the cancer. The simulation study supports such a reduction of cancer cells as displayed in Figure 3.Thesimulation study in figure 4 show that all cancer tumor cells are died,cancer stem cells are reduced in size but not compeltly removed, health cells are reduced to the normal size, and estrogen is rapidely decreased to normal size. Thus, in this case, an adjuvant treatment i.e., an additive that enhances the effectiveness of medical treatment, is required.
Further, the simulation study supports that the application of brody growthfunction to breast cancer model is biological meaining full and matematical well defined. The model constructed and studied in this paper is local stable at the equlibruim point.