Next Article in Journal / Special Issue
Solution of Optimal Harvesting Problem by Finite Difference Approximations of Size-Structured Population Model
Previous Article in Journal
Impact of Thermal Radiation and Heat Source/Sink on Eyring–Powell Fluid Flow over an Unsteady Oscillatory Porous Stretching Surface
Previous Article in Special Issue
Novel Spreadsheet Direct Method for Optimal Control Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control Analysis of a Mathematical Model for Breast Cancer

by
Segun Isaac Oke
*,
Maba Boniface Matadi
and
Sibusiso Southwell Xulu
Department of Mathematical Sciences, University of Zululand, Private Bag X1001, KwaDlangezwa 3886, South Africa
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2018, 23(2), 21; https://doi.org/10.3390/mca23020021
Submission received: 26 March 2018 / Revised: 19 April 2018 / Accepted: 20 April 2018 / Published: 24 April 2018
(This article belongs to the Special Issue Optimization in Control Applications)

Abstract

:
In this paper, a mathematical model of breast cancer governed by a system of ordinary differential equations in the presence of chemotherapy treatment and ketogenic diet is discussed. Several comprehensive mathematical analyses were carried out using a variety of analytical methods to study the stability of the breast cancer model. Also, sufficient conditions on parameter values to ensure cancer persistence in the absence of anti-cancer drugs, ketogenic diet, and cancer emission when anti-cancer drugs, immune-booster, and ketogenic diet are included were established. Furthermore, optimal control theory is applied to discover the optimal drug adjustment as an input control of the system therapies in order to minimize the number of cancerous cells by considering different controlled combinations of administering the chemotherapy agent and ketogenic diet using the popular Pontryagin’s maximum principle. Numerical simulations are presented to validate our theoretical results.

1. Introduction

Cancer is a generic name that refers to a group of diseases in which normal cells divide uncontrollably, that is, grow more rapidly than normal cells, and may eventually spread to other parts of the body by a process called metastasis [1]. According to the National Cancer Registry [2], cancer kills more people than tuberculosis (TB), AIDs and malaria combined. Statistics show that cancer related deaths amounted to about 8.2 million in 2010. The mortality rate from cancer is projected to continue to rise, with an estimated 13 million deaths by 2030 [3]. The most common types of cancer include: breast cancer, prostate cancer, brain cancer, lung cancer and skin cancer among others.
According to the [3] report, breast cancer is the most common invasive cancer in females worldwide. The formation of breast cancer can occur in the inner lining of the milk ducts, known as ductal carcinoma, or in the lobules of the breast, known as lobular carcinoma [4]. Breast cancer is one of the most widely recognized obstructive diseases in females around the world. The disease has presently been named as the most dangerous cancer in women [3]. However, little is known on the causes of the ailment. There are three major breast cancer risk factors namely hormonal imbalance (estrogen), genetic (family history), and environmental (poor diet, alcohol consumption, smoking, exposure to toxin, etc.) [5]. Surgery, chemotherapy, radiation therapy, hormonal therapy, hyperthermia, targeted therapy and ketogenic diet [5,6] amongst other therapeutics are used to inhibit tumor growth or kill the tumor cells in the body. However, each treatment has side effects attributed to it, for example, hair loss, vomiting, nausea and fatigue. Adverse effects occur as a result of chemotherapy, which is not able to differentiate between normal cells and tumor cells, consequently killing both of them [3].
Several dietary components and supplements have been examined as possible cancer prevention agents. Until recently, a few studies, such as [6,7,8], investigated diet as a possible adjuvant to cancer treatment, which includes a ketogenic diet. A ketogenic diet consists of high edible fat with moderate or low protein content and very low carbohydrates, which forces the body to burn fat instead of glucose for adenosine triphosphate (ATP) synthesis [6,9].
It is well-known that a mathematical model is a capable device used to investigate the spread of non-infectious diseases and to provide important insights into disease behaviors and control [10,11]. Over the years, it has become an important tool in comprehending the dynamics of diseases and in decision making processes regarding a medical intervention program for controlling breast cancer in many nations [12]. For instance, [13] explored the role of mathematical modeling on the optimal delivery of a combination therapy for tumors and to improve on the delivery of anti-tumor drugs.
Old and recent studies such as [4,12,14,15,16] amongst others have shown that mathematical modeling is a widely used tool for resolving questions on public health. For instance, it was used during the time of Bernoulli (on modeling the dynamics of Smallpox) in 1760 [17]. Kermack and McKendrick [14,15] and some other recent studies by [4,12,13,18,19,20,21] show that mathematical modeling is useful in solving the problem of epidemiology. However, these studies reveal that much has not been done in terms of the mathematical modeling of a nutritional diet (ketogenic diet) as a control or therapy on tumor cells. Hence, we improved the model in [19] for this paper by incorporating time dependent control parameters (use of ketogenic diet, immune booster, and anti-cancer drugs) based on the assumption that there is an interaction between normal cells and tumor cells that is due to a mutation in DNA as a result of excess estrogen in the body system [4,19,22].
Furthermore, we analyzed and applied an optimal control to the improved model to determine the possible impacts of ketogenic-diet use and anti-cancer drugs as a treatment on tumor cells. We carried out a rigorous qualitative optimal control analysis of the resulting model and found the necessary conditions for optimal control of the disease using Pontryagin’s maximum principle [11,23,24,25] in order to determine the optimal strategies for controlling the metastatic of the tumor cells.
This paper is organized as follows: In Section 2, four compartment models of ODEs to study the dynamics of breast cancer are developed. In Section 3, the existence of equilibria, their stabilities and basic reproductive numbers are discussed. In Section 4, an uncertainty and sensitivity analysis to check the most sensitive parameters in the model are discussed. In Section 5, an optimal control problem according to the model is proposed and an optimal solution is proffered. Numerical simulations are illustrated by implementing the forward and backward finite difference scheme in Section 6, while concluding remarks are provided in Section 7.

2. Model Formulation

Based on the existing model in [19], we developed a model by assuming the logistic (Verhulst) growth of a cell population and basic competition between normal cells and tumor cells. We considered the immune cells compartment to comprise Natural Killer cells (NK) and CD8+ T-cells as in [19] and we used a similar equation to model the immune response dynamic by introducing immune booster (ketone bodies) and anti-cancer drug efficacy.
We adapted an estrogen equation as presented in a model by [26]. Pinho and his co-workers in [26] considered that when a chemotherapy agent is continuously infused into the body and engulfed by different cell populations, natural death can occur. Excess estrogen was used in a similar way and assumed to be saturated daily through birth control (constant source rate) (1 − k). This was introduced to serve as anti-cancer drug efficacy (e.g., Tamoxifen) in order to bind estrogen receptors positive and to reduce excess estrogen from promoting cancer growth [27].
In this study, a model that splits the entire population of cells of the human breast tissues at any given period of time P ( t ) was reflected upon. Hence, normal cells compartment, represented by N ( t ) in the form of epithelial cells that constitute the breast tissue is described. The cells are assumed to develop and die normally as they have unaltered DNA that control all cell activities. It was suggested that the normal cells and tumor cells compete for nutrients and other resources in a small volume, which is the competition model used by [28]. Normal cells are represented by
d N d t = N α 1 μ 1 N 2 ϕ 1 N T ( 1 k ) λ 1 N E .
The first term represents the logistic growth rate α 1 of the normal cells, which are breast tissues that are made-up of epithelial cells. The second term represents the natural death rate of normal cells. ϕ 1 represents the rate at which normal cells inhibit due to an alteration in DNA that is responsible for cancer cells having an uncontrolled cycle that normal cells do not have [19]. The final term describes the gene transactivation that can be a contributing growth factor responsible for the estrogen stimulation of breast cancer, which can result in damage of DNA. Thus, there will be a reduction in the population of normal cells N ( t ) = being transformed into tumor cells by λ 1 N E where λ 1 represents the tumor formation rate resulting from DNA mutation caused by the presence of excess estrogen [4]. However, (1 − k) represents the effectiveness of anti-cancer drugs (Tamoxifen).
The tumor cells compartment can be denoted by T ( t ) in the form of an abnormal mass of tissue. Tumors are classic signs of inflammation, and can be benign or malignant (cancerous). Their names usually reflect the kind of tissue from where they arise, for example in breast or brain cancer, among others. There are about 51 breast cancer cell lines that mirror the 145 primary breast tumors [29]. These can be classified into two major branches: the Luminal, which has estrogen receptors ( E S R 1 + v e ) , and the Basal-like, which has no estrogen receptors ( E S R 1 v e ) . A homogeneous luminal type of cancer cells in the form of MDAMB361, MCF-7, BT474, T47D and ZR75 of the cell lines [19] are then assumed to be
d T d t = T α 2 d μ 2 T 2 γ 2 M T μ 5 T + ( 1 k ) λ 1 N E .
The first term of the equation is a limited growth term for tumor cells that depends on the rate of parameter d (ketogenic diets). Although, if d = 0 , tumor cells are automatically eradicated, but any DNA mutation that is caused by excess estrogen will repopulate the tumor cells again λ 1 N E . The induced death rate μ 5 is as a result of tumor starvation of nutrients, glucose and so on from the body system during the ketogenic diet, which alters nutrition. We assumed that γ 2 is the rate at which tumor cells are being removed due to the effectiveness of immune response.
The immune response compartment is represented by M ( t ) in the form of natural killer (NK) cells and CD8+ T cells. Their growth may be stimulated by the presence of the tumor and they can destroy tumor cells through the kinetics process. We also assumed that the presence of a detectable tumor in a body system does not necessarily imply that the tumor has completely escaped active immunosurveillance. However, a tumor is immunogenic. It is possible that the immune response may not be sufficient on its own to completely combat the rapid growth of the tumor cells population and their eventual development into a tumor.
d M d t = s β + ρ M T ω + T γ 3 M T μ 3 M ( ( 1 k ) λ 3 M E g + E )
The constant source parameter s denotes the source rate of immune response fully infused in the body daily. We introduced immune booster β (a supplement such as ketone bodies) to assist immune response whenever tumor cells overpower immune cells in order to activate the immune response and fight the cancer cells. The next term is a nonlinear growth term for immune response where ρ the rate of immune response is and ω is the immune cell threshold [12]. We denoted γ 3 as the rate at which immune response is inactivated upon interacting with tumor cells while μ 3 represents the immune cells natural death rate as a result of necrosis. The final term explains a limited rate at which estrogen suppresses immune cells activation where λ 3 is the rate of immune suppression and g is the estrogen threshold [19].
Finally, we considered estrogen compartment denoted by E ( t ) . Estrogen is a female steroid hormone that is produced by the ovaries in lesser amounts, and by the adrenal cortex, placenta and male testes. Estrogen helps to control and guide sexual development, including the physical changes associated with puberty [11,30]. However, an increase in estrogen levels can lead to the growth of the tumor cells. It also serves as a mitogen by triggering cell division in breast tissue [30]. Estrogen acts as a carcinogen by directly damaging DNA, forcing healthy epithelial cells to have a higher likelihood of malignant conversion [5,30].
d E d t = ( 1 k ) ϵ μ 4 E
The process of constantly replenishing excess estrogen is denoted by ϵ . We assumed that the majority of cancer cells are estrogen-receptor positive and only a small proportion of epithelial cells are estrogen-receptor positive, which can only be blocked by the anti-cancer drug ( 1 k ) Tamoxifen. μ 4 is the rate at which estrogen is being washed out from the body system. Thus, system (5) is our modified model.
d N d t = N ( α 1 μ 1 N ϕ 1 T ) ( 1 k ) ( λ 1 N E ) d T d t = T ( α 2 d μ 2 T ) γ 2 M T μ 5 T + ( 1 k ) ( λ 1 N E ) d M d t = s β + ρ M T ω + T γ 3 M T μ 3 M ( ( 1 k ) λ 3 M E g + E ) d E d t = ( 1 k ) ϵ μ 4 E

3. Model Analysis

3.1. Boundedness and Positivity of Solutions

The system of Equation (5) has an initial condition by
N ( 0 ) = N 0 0 ,   T ( 0 ) = T 0 0 ,   M ( 0 ) = M 0 0 ,   and   E ( 0 ) = E 0 0
since our model is to investigate cellular populations, therefore all the variables and parameters of the model are non-negative. Based on the biological finding, the system of Equation (5) will be studied in the following region such as:
Δ = { ( N , T , M , E ) + 4 }
The following theorem assures that the system of Equation (5) is well-posed such that solutions with non-negative initial conditions remain non-negative for all 0 < t < , and therefore makes the variable biologically meaningful. Hence, we have the following result:
Theorem 1.
The region Δ + 4 is positively invariant with respect to the system of Equation (5) and non-negative solution exists for all time 0 < t < .
Proof: 
Let Δ = Δ c + 4 with Δ = { ( N , T , M , E ) + 4 : N α 1 μ 1 } , then the solutions (N (t), T(t), M(t), E(t)) of system (5) are positive t 0 . It is obvious from the first compartment of system (5) that
d N d t N ( t ) α 1 μ 1 N 2 ( t ) .
Solving with Bernoulli method and taking N ( 0 ) = N 0 , we have
N ( t ) α 1 μ 1 + k α 1 e α 1 t
with
k = α 1 N 0 μ 1 N 0 α 1
N 0 = α 1 μ 1 + k α 1 .
Then,
N ( t ) α 1 μ 1 + ( α 1 N 0 μ 1 N 0 ) e α 1 t
N ( t ) α 1 μ 1   a s   t
hence, N ( t ) > 0 , t > 0 and if and only if ( 1 k ) 0 [31].
Consequently, it can be shown that T ( t ) > 0 ,   M ( t ) > 0 ,   a n d   E ( t ) > 0   t > 0 . This completes the proof. □

3.2. The Equilibrium Points of System (5)

The steady states occur by setting the left hand side (LHS) of system (5) to zero, i.e.,
d N d t = d T d t = d M d t = d E d t = 0
The model system admits six steady states in which there are four dead equilibria, one tumor-free equilibrium point and one co-existing equilibrium point P = ( N * , T * , M * , E * ) where N * , T * , M * , E * represent the tumor-free equilibrium values for the normal cells, tumor cells, immune cells and estrogen hormone respectively. We have N * > 0 , M * > 0 , E * > 0 since cell populations are non-negative and real. Therefore, all parameters s, β , g , μ 1 , μ 3 , μ 4 , ϵ , λ 3 , k , α 1 , and λ 1 are positive.
Tumor-Free equilibrium point
P 0 = ( α 1 μ 4 ( 1 k ) 2 λ 1 ϵ μ 1 μ 4 , 0 , s β ( g μ 4 + ( 1 k ) ϵ ) μ 3 ( g μ 4 + ( 1 k ) ϵ ) + ( 1 k ) 2 λ 3 ϵ , ( 1 k ) ϵ μ 4 )
Type 1 Dead equilibrium point
P d 1 = ( 0 , 0 , s β ( g μ 4 + ( 1 k ) ϵ ) μ 3 ( g μ 4 + ( 1 k ) ϵ ) + ( 1 k ) 2 λ 3 ϵ ) , ( 1 k ) ϵ μ 4 )
Type 2 Dead equilibrium point
P d 2 = ( 0 , d α 2 γ 2 m 1 * μ 5 μ 2 , m 1 * , ( 1 k ) ϵ μ 4 )
Type 3 Dead equilibrium point
P d 3 = ( 0 , d α 2 γ 2 m 2 * μ 5 μ 2 , m 2 * , ( 1 k ) ϵ μ 4 )
Type 4 Dead equilibrium point
P d 4 = ( 0 , d α 2 γ 2 m 3 * μ 5 μ 2 , m 3 * , ( 1 k ) ϵ μ 4 )
Co-existing equilibrium point
P e = ( N 4 * , T 4 * , M 4 * , E 4 * )

3.3. The Reproductive Number and Tumor-Free Equilibrium Point

In this section, we mainly analyzed the stability behaviors of system (5) by means of eigenvalues. We apply Hartman–Grobman Theorem which states that in the neighborhood of a hyperbolic equilibrium point, a nonlinear dynamical system is topologically equivalent to its linearization [32].
Theorem 2.
The tumor-free equilibrium (TFE) point P 0 of system (5) is locally asymptotically stable if R 0 < 1 , otherwise unstable.
Proof. 
Linearizing system (5) around TFE P 0 , we obtained the following Jacobian matrix J ( P 0 )
J = ( 2 μ 1 λ 1 ( 1 k ) 2 ϵ α 1 μ 1 μ 4 ( 1 k ) 2 λ 1 μ 1 ϵ μ 1 μ 4 ( 1 k ) 2 λ 1 Φ 1 ϵ Φ 1 α 1 μ 4 μ 1 μ 4 0 B 6 ( 1 k ) 2 λ 1 ϵ μ 4 B 3 0 B 6 0 B 4 B 5 B 7 0 0 0 μ 4 )
J ( P 0 ) = ( B 0 B 2 0 B 6 B 1 B 3 0 B 6 0 B 4 B 5 B 7 0 0 0 μ 4 )
| J ( P 0 ) | = | B 0 δ B 2 0 B 6 B 1 B 3 δ 0 B 6 0 B 4 B 5 δ B 7 0 0 0 μ 4 δ | = 0
Then the characteristic equation at P 0 of the linearized system of the model (5) is given below.
Obviously, there exists two negative characteristic roots
δ 1 = μ 4 ,   δ 2 = B 5
However, we only need to consider
δ 2 ( B 0 + B 3 ) δ + B 0 B 3 B 1 B 2 = 0
δ 2 ( B 0 + B 3 ) δ + B 0 B 3 ( 1 B 1 B 2 B 0 B 3 ) = 0
from (6), we have basic reproduction number
R 0 = B 1 B 2 B 0 B 3
δ 2 ( B 0 + B 3 ) δ + B 0 B 3 ( 1 R 0 ) = 0 ,
where
B 0 = 2 μ 1 λ 1 ( 1 k ) 2 ϵ α 1 μ 1 μ 4 ( 1 k ) 2 λ 1 μ 1 ϵ μ 1 μ 4 , B 1 = ( 1 k ) 2 λ 1 ϵ μ 4 , B 2 = ( 1 k ) 2 λ 1 Φ 1 ϵ Φ 1 α 1 μ 4 μ 1 μ 4
B 3 = α 2 μ 3 d ψ * + α 2 λ 3 d ϵ ( 1 k ) 2 γ 2 s β ψ * μ 5 μ 3 ψ * μ 5 λ 3 ϵ ( 1 k ) 2 μ 3 ψ * + ( 1 k ) 2 λ 3 ϵ
where
ψ * = ( g μ 4 + ( 1 k ) ϵ )
B 4 = s β ( g μ 4 + ( 1 k ) ϵ ( ρ γ 3 ω ) ω μ 3 ( g μ 4 + ( 1 k ) ϵ )
B 5 = μ 3 μ 4 g + μ 3 ( 1 k ) ϵ + λ 3 ( 1 k ) 2 ϵ g μ 4 + ( 1 k ) ϵ , B 6 = ( ( 1 k ) λ 1 α 1 μ 4 ( 1 k ) 3 λ 1 2 ϵ μ 1 μ 4 )
B 6 = ( 1 k ) λ 1 α 1 μ 4 ( 1 k ) 3 λ 1 2 ϵ μ 1 μ 4 , B 7 = ( μ 4 2 ( 1 k ) λ 3 g s β ( g μ 4 + ( 1 k ) ϵ ) ( g μ 4 + ( 1 k ) ϵ ) 2 [ μ 3 ( g μ 4 + ( 1 k ) ϵ ) + ( 1 k ) 2 λ 3 ϵ ] )
Here, we can apply the Routh-Hurwitz criterion namely,
( i )   T r ( A ) < 0    ( ii )   D e t ( A ) > 0
provided
a 0 = 1 > 0 ,   a 1 = ( B 0 + B 3 ) < 0 ,   B 0 B 3 ( 1 R 0 ) > 0   if   R 0 < 1 , B 0 B 3 > B 1 B 2 ,   B 1 > 0 ,   B 2 > 0 ,   B 3 < 0 ,   B 0 < 0
Since the Routh–Hurwitz criterion holds, all the eigenvalues are negative, i.e., δ 3 < 0 and δ 4 < 0 . Therefore, the TFE point of system (5) is locally asymptotically stable if (7) R 0 < 1 otherwise unstable. □
The epidemiological implication of the above result is that the tumor cells that are governed by system (5) can be eliminated from the population (normal cells or breast tissues) whenever an influx by tumor cells into the normal cells is small, such that R 0 < 1 . Therefore, the existence of a tumor-free equilibrium in this case depends on the estrogen level.
Theorem 3.
The Type 1 Dead equilibrium point P d 1 of system (5) is locally asymptotically stable if
( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1 ,
otherwise unstable.
Proof. 
Linearizing system (5) around the Type 1 Dead free equilibrium point P d 1 , we obtained the following Jacobian matrix J ( P d 1 )
J = ( α 1 ( 1 k ) λ 1 E 0 * 0 0 0 ( 1 k ) λ 1 E 0 * d α 2 γ 2 M 0 * μ 5 0 0 0 ρ ω M 0 * γ 3 M 0 * ω 2 ω 2 ( μ 3 ( g + E 0 * ) + ( 1 k ) λ 3 E 0 * ( g + E 0 * ) ) λ 3 g M 0 * ( 1 k ) ( g + E 0 * ) 2 0 0 0 μ 4 )
J ( P d 1 ) = ( C 0 0 0 0 C 1 C 2 0 0 0 C 3 C 4 C 5 0 0 0 μ 4 )
| J ( P 1 ) | = | C 0 δ 0 0 0 C 1 C 2 δ 0 0 0 C 3 C 4 δ C 5 0 0 0 μ 4 δ | = 0
Clearly, two eigenvalues of the system (5) at P d 1 are negative and real
δ 1 = μ 4
and
δ 2 = C 4 ( μ 3 ( g μ 4 + ( 1 k ) ϵ ) + ( 1 k ) 2 λ 3 ϵ g μ 4 + ( 1 k ) ϵ )
while the remaining two eigenvalues are obtained from 2 × 2 matrix.
A = ( C 0 0 C 1 C 2 )
Applying the Routh-Hurwitz criterion stated above; we have
( i )   T r ( A ) = C 0 + C 2 ( α 1 μ 4 ( 1 k ) 2 λ 1 ϵ μ 4 + ( μ 3 α 2 d γ 2 s β μ 5 μ 3 ) A * + d α 2 λ 3 ( 1 k ) 2 ϵ ( 1 k ) 2 λ 3 μ 5 ϵ μ 3 A * + ( 1 k ) 2 λ 3 ϵ ) < 0 i f   α 1 ( 1 ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 0 , ( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1
Therefore, T r ( A ) < 0
( ii )   D e t ( A ) = C 0 C 2 ( ( α 1 μ 4 ( 1 k ) 2 λ 1 ϵ μ 4 ) ( ( μ 3 α 2 d γ 2 s β μ 5 μ 3 ) A * + d α 2 λ 3 ( 1 k ) 2 ϵ ( 1 k ) 2 λ 3 μ 5 ϵ μ 3 A * + ( 1 k ) 2 λ 3 ϵ ) ) > 0 i f   α 1 ( 1 ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 0   p r o v i d e d   ( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1
and
{ ( α 1 ( 1 ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) ) ( ( ( μ 3 α 2 d γ 2 s β μ 5 μ 3 ) A * + d α 2 λ 3 ( 1 k ) 2 ϵ ( 1 k ) 2 λ 3 μ 5 ϵ μ 3 A * + ( 1 k ) 2 λ 3 ϵ ) ) } > 0
implies that D e t ( A ) > 0 . Thus, the remaining eigenvalues δ 3   a n d   δ 4 are negative and real since R-H criterion has been satisfied. Hence, the type 1 Dead equilibrium point P d 1 of the system (5) is locally asymptotically stable if ( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1 . □
Epidemiologically it is implied that the net growth of the tumor cells must be more than the immune cells values in order to have the tumor cells overpower the normal cells as the reactivation of the immune cells is due to the estrogen effects that are greater than the reactivation of the immune cells due to the tumor effect. However, ketogenic diet is inactive at the type 1 Dead equilibrium point.
Theorem 4.
The Type 2 Dead equilibrium point P d 2 of system (5) is locally asymptotically stable if
( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1
ω > A * μ 2 ( μ 2 C * ρ γ 3 A * C * + μ 2 μ 3 C * + ( 1 k ) 2 λ 3 ϵ 1 )
otherwise unstable.
Proof. 
We linearized system (5) around the Type 2 Dead free equilibrium point P d 2 , we obtained the following Jacobian matrix J ( P d 2 ) at P d 2 = ( 0 , d α 2 γ 2 m 1 * μ 5 μ 2 , m 1 * , ( 1 k ) ϵ μ 4 )
J = ( α 1 ( 1 k ) λ 1 E 1 * 0 0 0 ( 1 k ) λ 1 E 1 * Q 2 γ 2 T 1 * 0 0 Q 3 Q 5 λ 3 g M 1 * ( 1 k ) ( g + E 1 * ) 2 0 0 0 μ 4 )
where
Q 2 = ( d α 2 2 μ 2 T 1 * γ 2 M 1 * μ 5 ) , Q 3 = ( ρ ω μ 2 2 M 1 * γ 3 M 1 * ( ω μ 2 + d α 2 γ 2 M 1 * μ 5 ) 2 ( ω μ 2 + d α 2 γ 2 M 1 * μ 5 ) 2 ) , Q 5 = ( ρ T 1 * ω + T 1 * γ 3 T 1 * μ 3 ( 1 k ) λ 3 E 1 * g + E 1 * )
| J ( P 2 ) | = ( Q 0 0 0 0 Q 1 Q 2 Q 4 0 0 Q 3 Q 5 Q 6 0 0 0 μ 4 )
| J ( P 2 ) | = | Q 0 δ 0 0 0 Q 1 Q 2 δ Q 4 0 0 Q 3 Q 5 δ Q 6 0 0 0 μ 4 δ | = 0
Clearly, one of the eigenvalues of the system (5) at | J ( P 2 ) | is negative and real, i.e., δ 1 = μ 4 . However, the remaining can be analyzed by simple calculation.
( Q 0 δ ) ( Q 2 δ ) ( Q 5 δ ) = 0 Q 5 = δ 2 ,   Q 2 = δ 3 ,   Q 0 = δ 4
where
Q 0 = α 1 μ 4 ( 1 k ) 2 λ 1 ϵ μ 4 , Q 2 = γ 2 M 1 * α 2 d + μ 5
Q 5 = ( A * ρ μ 2 2 A * C * γ 3 ) ( d α 2 γ 2 M 1 * μ 5 ) C * μ 2 ( μ 3 A * + ( 1 k ) 2 λ 3 ϵ ) μ 2 A * C *
where A * = ( g μ 4 ( 1 k ) ϵ ) and C * = ( ω μ 2 + d α 2 γ 2 M 1 * μ 5 ) .
It follows the following conditions
(i) Q 0 < 0 if, α 1 1 , 0 k < 1 and ( ( 1 k ) 2 λ 1 ϵ α 1 μ 4 ) > 1 ;
(ii) Q 5 < 0 provided A * > 0 , 0 k 1 and ω > A * μ 2 ( μ 2 C * ρ γ 3 A * C * + μ 2 μ 3 C * + ( 1 k ) 2 λ 3 ϵ 1 ) . □

3.4. Co-Existing Equilibrium Points

Theorem 5.
The co-existing equilibrium point P e of system (5) is stable if the following Routh–Hurwitz criterion is satisfied,
T r a c e ( A ) = ( V 0 + V 3 + V 6 μ 4 ) < 0 D e t ( A ) = ( μ 4 ( V 0 V 6 V 3 + V 0 V 4 V 5 + V 1 V 2 V 6 ) ) > 0 ,
otherwise unstable.
Proof. 
We analyzed and linearized system (5) around the co-existing equilibrium point P e , we obtained the following Jacobian matrix J ( P e ) at P e = ( N 4 * , T 4 * , M 4 * , E 4 * ) where N 4 * , T 4 * , M 4 *   &   E 4 * represent the coexisting equilibrium values for normal cells, tumor cells, immune cells, and estrogen levels respectively.
A co-existing equilibrium state exists when all cells populations would have survived the competition.
N 4 * = 2 ( 1 k ) 4 λ 1 4 μ 1 μ 4 ϵ 2 + ϕ 1 α 1 2 μ 4 2 μ 1 2 ( 1 k ) 2 μ 1 μ 4 2 α 1 λ 1 ϕ 1 ϵ 2 α 1 ϕ 1 2 μ 1 μ 4 3 2 ( 1 k ) 2 α 1 μ 1 μ 4 2 λ 1 ϵ 2 ϕ 1 α 1 μ 1 2 μ 4 3 2 ( 1 k ) 2 μ 1 2 μ 4 2 λ 1 ϕ 1 ϵ
T 4 * = α 1 2 μ 1 μ 4 2 + 2 α 1 μ 1 μ 4 2 ϕ 1 2 ϕ 1 α 1 μ 1 μ 4 2 2 ( 1 k ) 2 μ 1 μ 4 λ 1 ϕ 1 ϵ
M 4 * = G * 2 Z * ( 1 k ) 2 λ 1 ϵ + ( α 1 2 α 2 μ 1 μ 4 3 d + 2 α 1 α 2 μ 1 μ 4 3 ϕ 1 d μ 4 3 μ 5 α 1 2 2 μ 4 3 μ 1 μ 5 α 1 ϕ 1 ) G * μ 3 α 1 4 μ 1 2 μ 4 5 4 α 1 2 μ 1 2 μ 2 μ 4 5 ϕ 1 4 ϕ 1 2 α 1 2 μ 1 2 μ 4 5 μ 2 G * 2 Q * μ 4
E 4 * = ( 1 k ) ϵ μ 4
where
G * 2 = 2 ϕ 1 α 1 μ 1 μ 4 2 2 ( 1 k ) 2 μ 1 μ 4 λ 1 ϕ 1 ϵ
Z * = 2 ( 1 k ) 4 λ 1 2 μ 1 μ 4 ϵ 2 + ϕ 1 α 1 2 μ 4 2 μ 1 2 ( 1 k ) 2 μ 1 μ 4 2 α 1 λ 1 ϕ 1 ϵ 2 α 1 ϕ 1 2 μ 1 μ 4 3 2 ( 1 k ) 2 α 1 μ 1 μ 4 2 λ 1 ϵ 2 ϕ 1 α 1 μ 1 2 μ 4 3 2 ( 1 k ) 2 μ 1 2 μ 4 2 λ 1 ϕ 1 ϵ
Q * = α 1 2 μ 1 μ 4 2 γ 2 2 α 1 μ 1 μ 4 2 ϕ 1 γ 2 2 ϕ 1 α 1 μ 1 μ 4 2 2 ( 1 k ) 2 μ 1 μ 4 λ 1 ϕ 1 ϵ
J = ( ( α 1 2 μ 1 N 4 * ( 1 k ) λ 1 E 4 * ) N 4 * ϕ 0 V 7 ( 1 k ) λ 1 E 4 * ( d α 2 2 μ 2 T 4 * γ 2 M 4 * μ 5 ) γ 2 T 4 * V 7 0 V 4 V 6 λ 3 g M 4 * ( 1 k ) ( g + E 4 * ) 2 0 0 0 μ 4 )
A = ( V 0 V 2 0 V 7 V 1 V 3 V 5 V 7 0 V 4 V 6 V 8 0 0 0 μ 4 )
| A | = | V 0 V 2 0 V 7 V 1 V 3 V 5 V 7 0 V 4 V 6 V 8 0 0 0 μ 4 | = 0
We need to show that T r a c e ( A ) < 0 , that is
T r ( A ) = ( V 0 + V 3 + V 6 μ 4 ) < 0
= α 1 ( 1 A 0 ) 2 μ 1 N 4 * + d α 2 ( 1 μ 5 ) μ 4 + T 4 * ( ρ γ 3 ( ω T 4 * ) ) ω + T 4 * μ 3 ( 1 k ) 4 λ 3 ϵ ( g μ 4 + ( 1 k ) ϵ
Thus,
T r ( A ) < 0 ,   if   A 0 > 1 ,   μ 5 > 1 ,   ρ < γ 3 ( ω + T 4 * )   with   A 0 = ( 1 k ) 2 λ 3 ϵ α 1 μ 4
To show that,
| A | = ( μ 4 ( V 0 V 3 V 6 + V 0 V 4 V 5 + V 1 V 2 V 6 ) ) > 0
Let ζ 1 = μ 4 V 0 V 3 V 6 , ζ 2 = μ 4 V 0 V 4 V 5 , ζ 3 = μ 4 V 1 V 2 V 6
ζ 1 = ( α 1 ( 1 A 0 ) 2 μ 1 N 4 * ) ( d α 2 ( 1 μ 5 ) 2 μ 2 T 4 * γ 2 μ 4 * ) ) ( T 4 * ( ρ γ 3 ( ω T 4 * ) ) ω + T 4 * μ 3 ( 1 k ) 4 λ 3 ϵ ( g μ 4 + ( 1 k ) ϵ ) .
This implies that, ζ 1 > 0 is a positive, if A 0 > 1 , μ 5 > 1 , ρ < γ 3 ( ω + T 4 * ) with A 0 = ( 1 k ) 2 λ 3 ϵ α 1 μ 4
ζ 2 = ( α 1 ( 1 A 0 ) 2 μ 1 N 4 * ) ( M 4 * ( ω + T 4 * ) 2 ( ρ ω γ 3 ( ω + T 4 * ) ) ) ( γ 2 T 4 * )
This implies that, ζ 2 > 0 is a positive, if A 0 > 1 , μ 5 > 1 , ρ ω < γ 3 ( ω + T 4 * ) 2 , with A 0 = ( 1 k ) 2 λ 3 ϵ α 1 μ 4
ζ 3 = μ 4 ( A 0 ) ( ϕ 1 N 4 * ) ( T 4 * ( ρ γ 3 ( ω T 4 * ) ) ω + T 4 * μ 3 ( 1 k ) 4 λ 3 ϵ ( g μ 4 + ( 1 k ) ϵ )
This implies that ζ 3 < 0 is a negative and by Routh-Hurwitz criterion the system cannot be stable. Thus the co-existing equilibrium point is always unstable if the cells coexist where
V 0 = α 1 μ 4 2 μ 1 μ 4 N 4 * ( 1 k ) 2 λ 1 ϵ μ 4 , V 1 = ( 1 k ) 2 λ 1 ϵ μ 4 , V 2 = ϕ 1 N 4 * ,
V 3 = ( d α 2 2 μ 2 T 4 * γ 2 M 4 * μ 5 ) , V 4 = ρ M 4 * ω γ 3 M 4 * ( ω + T 4 * ) 2 ( ω + T 4 * ) 2 , V 5 = γ 2 T 4 *
V 6 = ρ T 4 * ( g μ 4 + ( 1 k ) ϵ ) γ 3 T 4 * ( ω + T 4 * ) ( g μ 4 + ( 1 k ) ϵ ) μ 3 ( ω + T 4 * ) ( g μ 4 + ( 1 k ) ϵ ) ( 1 k ) 2 ( ω + T 4 * ) λ 3 ϵ ( ω + T 4 * ) ( g μ 4 + ( 1 k ) ϵ )
V 7 = ( 1 k ) λ 1 N 4 * , V 7 = ( 1 k ) λ 1 N 4 * , V 8 = λ 3 μ 4 2 g M 4 * ( 1 k ) ( g μ 4 + ( 1 k ) ϵ ) 2

4. Uncertainty and Sensitivity Analysis

In this section, we explore the dependence of the model solutions on the parameter values. We are able to figure-out a feasible range of parameter values and determine the most critical parameters in the model. We employed a similar method, which is discussed in detail by [20,33], using Latin Hypercube Sampling (LHS) for studying the uncertainty analysis and the Partial Rank Correlation Coefficient (PRCC) for analyzing the sensitivity analysis indexes of the parameters. LHS/PRCC was ran and analyzed with a sample size of 100. The choice of this sample size is due to the fact that PRCC produces accurate results for a lower sample size compared to other technique, such as eFAST [33].
Uncertainty and sensitivity analysis were performed on all non-dimensional system parameters in the system (5) with the aim of determining the most sensitive parameters to the model. The parameter baseline values in Table 1 were varied in the range of 25%. Figure 1 displays a bar graph of PRCCs plotted against the homogeneous parameter value with tumor compartment as the baseline dependent variable. The parameters that are significantly positively correlated with tumor cells, at P < 0.05 level of significance, are α 1 , g while μ 1 , γ 3 , and ω are significantly negatively correlated. An increase in the production of normal cells α 1 , leads to higher numbers of normal cells, thus the higher the α 1 , the higher the normal cells. While Figure 2 displays a bar graph of PRCCs plotted against the homogeneous parameter value with tumor compartment as the baseline dependent variable. The most sensitive parameters are shown to be P v a l u e s of s , γ 2 , μ 3 and ρ are less than 0.01.

5. Analysis of Optimal Control

In this section, we formulated a corresponding optimal control problem for the model in the system (5) considering ketogenic-diet and anti-cancer drugs as control interventions to minimize the breast cancer and tumor burden at final time. The units of cells were normalized in order for the carrying capacity of normal cells to be kept above threshold of 0 t t f [34,35,36]. On the other hand, the aim is to reduce the tumor-size which indicates the degree of the disease in the body and it requires the application of as much anti-cancer drugs as much as possible. However, it also minimized the systemic cost, which is based on the quantities of anti-cancer drugs, since large drug concentrations can be harmful and cause toxic side effects. In brief, the drug doses were minimized because the smaller the dose, the better. Then, we formulated the objective functional J 1
J 1 ( u 1 , u 2 ) = 0 T f ( A 1 T ( t ) + A 2 E ( t ) + 1 2 A 3 u 1 2 ( t ) + 1 2 A 4 u 2 2 ( t ) ) d t
System equations (5) is subject to:
d N d t = N α 1 μ 1 N 2 ϕ 1 T N ( 1 u 1 ( t ) ) ( λ 1 N E ) d T d t = ( 1 u 2 ( t ) ) T α 2 μ 2 T 2 γ 2 M T μ 5 T + ( 1 u 1 ( t ) ) ( λ 1 N E ) d M d t = s β + ρ M T ω + T γ 3 M T μ 3 M ( ( 1 u 1 ( t ) ) λ 3 M E g + E ) d E d t = ( 1 u 1 ( t ) ) ϵ μ 4 E
J 1 involves a quadratic control. In [37,38,39,40,41], it was established that quadratic control in the treatment terms has the added benefit of keeping the tumor in check both when it is small or large in size. The authors further explained that the quadratic control allows for a weaker treatment to minimize the toxic side-effects while permitting the system to maintain a low tumor size.
Furthermore, for us to address the tumor-to-therapy trade-off, we established the existence of an optimal control; by following the approach in [37,41,42], which required an analysis of the super-solutions (that is, the upper bounds on solutions) of the system (5). As soon as we were able to show that the system is bounded, we established the existence of an optimal control using a result from [43]. In addition, we proved that there exists an optimal control that minimizes the objective functional; using the established approach of [37,38,39,40,44]. We use the fact that super-solutions N ¯ , T ¯ , M ¯ , E ¯ of
d N ¯ d t = N α 1 , d T ¯ d t = T α 2 ( 1 u 2 ) , d M ¯ d t = s β + ρ M T ω + T , d E ¯ d t = 1
are bounded on a finite time interval. Since the sub-solutions are zero, the result obtained shows that our system is bounded. Since we have a bounded system, our next task was to establish the existence of the optimal control using a result from [43].

Existence of an Optimal Control

Theorem 6.
Given the objective functional in (8), where U = { u i * ( t ) ,   L e b e s g u e m e a s u r e : 0 u i * ( t ) 1 , t [ 0 , t ] } subject to system (9) with N ( 0 ) = N 0 , T ( 0 ) = T 0 , M ( 0 ) = M 0 , and E ( 0 ) = E 0 , then there exists an optimal control u i * ¯ such that m i n u i * ¯ ( t ) [ 0 , 1 ] J 1 ( u i * ¯ ) = J 1 ( u i * ( t ) ) if the following conditions holds:
  • f is not empty;
  • The admissible control set U is closed and convex;
  • Each right hand side of the state system is continuous, is bounded above by the sum of the bounded control and the state, and can be written as a linear function of u i * ( t ) ¯ with coefficients depending on time and the state;
  • The integrand of J 1 ( u i * ¯ ) is convex on U and is bounded below by c 2 + c 1 u ¯ 2 with c 1 > 0 .
Proof. 
Since the system (9) has bounded coefficients and the solutions are bounded on the finite time interval, we can apply the result of [45], to obtain the existence of the solution of the system (9). Furthermore, we note that U is closed and convex by definition. For the third conditions, the right hand side of the system (9) must be continuous. The right hand side is continuous since the denominators of all fractions from the right hand side of the system consists solely of positive entities. We let ϕ ( t , X ) be right hand side of the system (9) except for the terms of u i * ¯ and define.
| f ( t , X , u i * ) | = ϕ ( t , X ) + ( 0 λ 1 N E 0 u 1 ) , w i t h X = ( N T M E )
using the boundedness of the solutions (10), we have
| f ( t , X , u i * ) | | ( α 1 0 0 0 0 α 2 ( 1 u 2 ) 0 0 0 0 ρ 0 0 0 0 0 ) ( N T M E ) | + | ( 0 ( 1 u 1 ) λ 1 N E s β u 1 ϵ ) | c 1 ( | X | + | u i * | )
where c 1 depends on the coefficients of the system. For the fourth condition, we need to show
J ( t , T , E , ( 1 P i ) u i + P i V i ) ( 1 P i ) J ( t , T , E , u i ) + P i ( t , T , E , V i )
we analyze the difference of
J ( t , T , E , ( 1 P i ) u i + P i V i ) [ ( 1 P i ) J ( t , T , E , u i ) + P i ( t , T , E , V i ) ]
= T ( t ) + E ( t ) + ϵ 2 ( u i 2 2 P i u i 2 + P i 2 u i 2 + P i 2 V i 2 2 P i 2 V i 2 u i 2 + 2 P i V i u i ) ( T ( t ) + E ( t ) + ϵ 2 u i 2 ϵ 2 P i u i 2 + ϵ 2 P i V i 2 )
= ϵ 2 ( P i 2 P i ) ( u i V i ) 2
since, P i ( 0 , 1 ) implies ( P i 2 P i ) < 0 and ( u i V i ) 2 > 0 but ( P i 2 P i ) < 0 , which implies ϵ 2 ( P i 2 P i ) ( u i V i ) 2 is negative. This implies that,
J ( t , T , E , ( 1 P i ) u i + P i V i ) ( 1 P i ) J ( t , T , E , u i ) + P i ( t , T , E , V i )
Lastly,
T ( t ) + E ( t ) + ϵ 2 u i 2 ( t ) ϵ 2 u i 2 ( t ) c + ϵ 2 u i 2 ( t )
which gives c + ϵ 2 u i 2 ( t ) as the lower bound. With the existence of the optimal control established, we now characterized the optimal control using the Pontryagin’s maximum principle [11]. The constants A 1 , A 2 , A 3 and A 4 are a measure of the relative cost of the interventions over [ 0 , T ] . The optimal control problem is that of finding optimal functions ( u 1 * ( t ) , u 2 * ( t ) ) such that
J 1 ( u 1 * ( t ) , u 2 * ( t ) ) = min Ω J 1 ( u 1 ( t ) , u 2 ( t ) )
where
Ω = { u 1 ( t )   &   u 2 ( t ) : 0 u 1 ( t ) u 1 m a x ,0 u 2 ( t ) u 2 m a x , t [ 0 , T f ] }
Three different control strategies are explored. This approach can be used to test various options. However, we only looked at the following three alternatives:
  • Strategy 1: Anti-cancer drug treatment control on tumor cells (control u 1 ( t ) only);
  • Strategy 2: Ketogenic diet control on excess estrogen and tumor cells (control u 2 ( t ) only);
  • Strategy 3: Anti-cancer drug and ketogenic diet treatment combined control on tumor cells growth and excess estrogen (controls u 1 ( t ) and u 2 ( t ) ).
Thus, strategies (1–3) use the objective functional (8). We assumed that there are practical limitations on the maximum rate at which the anti-cancer treatment may be applied in a given time period. We defined the positive constant u m a x accordingly. We also define the set Ω of admissible controls to be all Lebesgue measurable functions that take on values in the control set [13,46,47] u = [ 0 , u m a x ] almost everywhere on [ 0 , T ] . We sought an optimal control u * Ω in (11) [13]. In order to find the optimal solutions, we first traced the Lagrangian and Hamiltonian for the optimal control problem (8) and (9). The Lagrangian of the optimal control problem is given by:
L ( N , T , M , E , u 1 , u 2 ) = A 1 T ( t ) + A 2 E ( t ) + 1 2 A 3 u 1 2 ( t ) + 1 2 A 4 u 2 2 ( t )
For the purpose of the necessary conditions for optimal control functions with the help of Pontryagin’s maximum principle [11]. We define the Hamiltonian, H for the control problem of the system (8) and (9)
H = L ( N , T , M , E , u 1 , u 2 ) + θ 1 N + θ 2 T + θ 3 M + θ 4 E
where L is the Lagrangian function (12),
H = ( A 1 T ( t ) + A 2 E ( t ) + 1 2 A 3 u 1 2 ( t ) + 1 2 A 4 u 2 2 ( t ) + θ 1 ( N α 1 μ 1 N 2 ϕ 1 T N ( 1 u 1 ( t ) ) ( λ 1 N E ) ) + θ 2 ( ( 1 u 2 ( t ) ) T α 2 μ 2 T 2 γ 2 M T μ 5 T + ( 1 u 1 ( t ) ) ( λ 1 N E ) ) + θ 3 ( s β + ρ M T ω + T γ 3 M T μ 3 M ( ( 1 u 1 ( t ) ) λ 3 M E g + E ) ) + θ 4 ( ( 1 u 1 ( t ) ) ϵ μ 4 E )
where θ 1 , θ 2 , θ 3 , θ 4 are the adjoints variables for the states N, T, M, E. However, with the help of Pontryagin’s Maximum Principle, we obtained a minimized Hamiltonian that minimizes the objective function or cost functional. We applied Pontryagin’s Maximum Principle [11], to characterize the optimal control pair u 1 *   &   u 2 * in the following result.
Theorem 7.
Given optimal control variables u 1 *   &   u 2 * and N * , T * , M *   &   E * are corresponding optimal state variables of the control system (8) and (9). Then there exists the adjoint variable θ i = ( θ 1 , θ 2 , θ 3 , θ 4 ) + 4 that satisfies the following equations.
d θ 1 d t = 2 θ 1 μ 1 N + ϕ 1 θ 1 T + ( θ 1 + θ 2 ) ( 1 u 2 ( t ) ) λ 1 E α 1 θ 1 d θ 2 d t = A 1 + θ 1 ϕ 1 N + θ 2 ( 2 T μ 2 + γ 2 M + μ 5 α 2 ( 1 u 2 ) ) + θ 3 ( γ 3 M ρ ω M ( ω + T ) 2 ) d θ 3 d t = θ 2 γ 2 T ρ θ 3 T + γ 3 θ 3 T + μ 3 θ 3 + θ 3 ( ( 1 u 1 ) λ 1 E g + E ) d θ 4 d t = A 2 + ( θ 1 θ 2 ) ( 1 u 1 ) λ 1 N θ 3 ( ( 1 u 1 ) λ 3 M g ( g + E ) 2 ) θ 4 μ 4
with transversality conditions
θ 1 ( T f ) = θ 2 ( T f ) = θ 3 ( T f ) = θ 4 ( T f ) = 0
The corresponding optimal controls u 1 *   &   u 2 * are given as,
u 1 * = m i n { m a x { 0 , 1 A 3 ( θ 2 λ 1 N * E * + θ 3 ϵ θ 1 λ 1 N * E * θ 3 λ 3 M * E * g + E * ) } ,1 }
and
u 2 * = m i n { m a x { 0 , 1 A 4 ( θ 2 α 2 T * ) } ,1 }
Proof. 
Let u 1 *   &   u 2 * be the given optimal control functions and N * , T * , M *   &   E * be the corresponding optimal state variables of the system (9) that minimize the cost functional or objective (8). Then by Pontryagin’s maximum principle [11], there exists adjoint variables (14) θ 1 , θ 2 , θ 3   &   θ 4 which satisfy the following equations
d θ 1 d t = H N , d θ 2 d t = H T , d θ 3 d t = H M , d θ 4 d t = H E
with transversality conditions
θ 1 ( T f ) = θ 2 ( T f ) = θ 3 ( T f ) = θ 4 ( T f ) = 0
where H is the Hamiltonian and defined as
H ( N , T , M , E , u 1 , u 2 , θ ) = L ( N , T , M , E , u 1 , u 2 ) + θ 1 N + θ 2 T + θ 3 M + θ 4 E
H = ( A 1 T ( t ) + A 2 E ( t ) + 1 2 A 3 u 1 2 ( t ) + 1 2 A 4 u 2 2 ( t ) + θ 1 ( N α 1 μ 1 N 2 ϕ 1 T N ( 1 u 1 ( t ) ) ( λ 1 N E ) ) + θ 2 ( ( 1 u 2 ( t ) ) T α 2 μ 2 T 2 γ 2 M T μ 5 T + ( 1 u 1 ( t ) ) ( λ 1 N E ) ) + θ 3 ( s β + ρ M T ω + T γ 3 M T μ 3 M ( ( 1 u 1 ( t ) ) λ 3 M E g + E ) ) + θ 4 ( ( 1 u 1 ( t ) ) ϵ μ 4 E )
from the optimality condition, we have
H u 1 = 0 ,   a t   u 1 = u 1 *   a n d   H u 2 = 0 ,   a t   u 2 = u 2 *
which implies that
0 = H u 1 = A 3 u 1 + θ 1 λ 1 N E θ 2 λ 1 N E + θ 3 λ 3 M E g + E θ 4 ϵ
0 = H u 1 = A 4 u 2 θ 2 α 2 T
Hence, we obtain (see [10])
u 1 * = 1 A 3 { θ 1 λ 1 N E + θ 4 ϵ θ 1 λ 1 N E θ 3 λ 3 M E g + E }
u 2 * = 1 A 4 { θ 2 α 2 T }
Thus we have (17) and (18).
By standard control arguments involving the bounds on the controls, we conclude that (15) and (16) can be written in this form
u 1 * = ( 0 i f   1 A 3 ( θ 1 λ 1 N E + θ 4 ϵ θ 1 λ 1 N E θ 3 λ 3 M E g + E ) < 0 1 A 3 ( θ 1 λ 1 N E + θ 4 ϵ θ 1 λ 1 N E θ 3 λ 3 M E g + E ) i f   0 1 A 3 ( θ 1 λ 1 N E + θ 4 ϵ θ 1 λ 1 N E θ 3 λ 3 M E g + E ) 1 1 i f   1 A 3 ( θ 1 λ 1 N E + θ 4 ϵ θ 1 λ 1 N E θ 3 λ 3 M E g + E ) > 1
and
u 2 * = ( 0 i f   1 A 4 ( θ 2 α 2 T * ) < 0 1 A 4 ( θ 2 α 2 T * ) i f   0 1 A 4 ( θ 2 α 2 T * ) 1 1 i f   1 A 4 ( θ 2 α 2 T * ) > 1
However, we discuss the numerical solution of the optimality system and the corresponding results of varying the optimal controls u 1   &   u 2 the parameter choices, and the interpretations from various cases.

6. Numerical Simulations and Discussion

A picture of the dynamical behavior of breast cancer cells in the presence of normal cells, tumor cells, immune cells, and estrogen is given by the numerical simulations of the model (5). The optimal control is acquired by solving the optimality system of four ordinary differential equations from the state variables and the adjoint system. An iterative scheme is used to solve the optimality system. All the numerical simulations were executed in MAPLE 18. We employed the forward-backward scheme method, beginning with an initial guess for optimal controls and solved the optimal state system forward in time and after that solved the adjoint state system backward in forward using the finite difference scheme in MAPLE. The two controls were then updated by using a convex combination of the previous controls as well as the characterization (17) and (18). The entire process was repeated until the values of the unknown at the previous iterations were closed to the one at the current iteration [39,41]. Key parameters are also noted in stabilizing the model in system (5), for example: ketogenic diet, anti-cancer, and immune booster. The initial values of variables are N(0) = 2000, T(0) = 800, M(0) = 500, E(0) = 20 and s = 1.3 × 10 4 adopted from [12]. All parameter values used for the numerical simulation are stated in Table 1 above.
Figure 3, indicates that the introduction of a ketogenic diet results in a reduction of activities of cancer cells and we also note that too much of a ketogenic diet will result in ketoacidosis. Ketoacidosis is the combination of ketosis and acidosis. Ketosis is the accumulation of substances called ketone bodies and acidosis is the increased acidity of the blood which can cause frequent urination (Polyuria), poor appetite, and a loss of consciousness. Therefore, our ketogenic diet’s parameter rate is best at d = 0.6 and it can complement the activity of the anti-cancer drug (Tamoxifen). Figure 4, shows the impact of anti-cancer drugs in reducing the production of excess estrogen in the system, but when there is less production of estrogen there will not be a rapid activation of the growth factor that expresses breast normal cells. However, the rapid production of estrogen results in abnormal breast cells expression, which will lead to breast cancer. Figure 5 shows the obvious effectiveness of anti-cancer drugs on tumor cells when there is no supply of nutrient or glucose to cancer cells.
Furthermore, Figure 6 illustrates that the red line β = 0 shows that during cancer formation the activities of both innate and adaptive reduces drastically, which is due to the expression of other proteins apart from those proteins that are responsible for the activation of the immune response, such as an immune booster introduced to the system, which reactivates the activities of the immune response towards the cancer cells.
The presence of abnormal estrogen level without anti-cancer drugs or a ketogenic diet will lead the system into critical condition and became unstable as shown in Figure 7. However, the system became stable as we introduced treatments, such as chemotherapy and the ketogenic diet as represented in Figure 8. In addition, Figure 9, indicates that there is DNA damage at λ 1 = 0 , which occurs naturally as a result of metabolic or hydrolytic processes. It is as a result of the Tumor Suppressor Gene (TSG), which is able to control the activity of DNA gene repair successfully. On the other hand, at λ 1 = 0.2 , 0.4 , 0.6 showed that TSG (such as BRCA 1, BRCA 2, P53) compromised the pathway that leads cells to grow uncontrollably and later form a tumor or it leads to accelerated aging.
However, the mathematical analysis of the model produced six equilibrium points. All the points have epidemiological implications in relation to explaining the dynamics of breast cancer growth. P 0 represents the situation where there is tumor-free equilibrium, that is when only tumor cell population has died off due to competition with other cells. P d 1 represents Type 1 dead equilibrium point where both normal cells and tumor cells die-off as a result of breast tissue removal through mastectomy surgery or death. This is because overtime the cancer cells which are depending on estrogen to develop into independent cells that grow regardless of estrogen receptors. P d 2 could be described by Type 2 dead equilibrium point where normal cells were only forced to extinction leaving the tumor cells surviving. P d 3 represent Type 3 dead equilibrium point which means immune system is weak and it cannot fight the tumor cells which eventually overpower normal cells and forced it to extinction. P d 4 show that Type 4 dead equilibrium point where ketogenic diet is not effective, immune booster is not active which lead to tumor cell over-compete normal cells as a result of infusion of excess estrogen to the body system.
We categorised this as ”dead” because biologically there is no recovery of damaged normal cells since they have died off of the cell population. It could be as a result of anti-cancer drug that destroy red blood cells which affected normal cells.

Effects of Control on the System (9)

By numerical simulation, optimal single control of anti-cancer drugs measure u 1 and ketogenic-diet optimal control measure u 2 are shown in Figure 10a,b respectively; where (red dots line) represented tumor cells and (solid green line) represented normal cells. Figure 10c is the use of combination of two control therapies which have significant impact on the increase of normal cells population against time. However, all the strategies are effectively restrain the tumor growth, they cannot totally eliminate a large tumor in 100 days. In Figure 11, optimal control using anti-cancer drugs and ketogenic diet as we optimized the system (54) with the objective function J for breast cancer model. It was observed that the combination of the two controls resulted in appreciable decreases in the number of tumor cells population in the presence of control (solid green line) while (dots red line) in the case of uncontrolled. However, tumor growth is driven to a very low but non-zero level.
Furthermore, it was noticed from Figure 12, that the level of estrogen was reduced drastically in the presence of controls (solid green line) against the constant increase level of estrogen (dots red line) in uncontrolled cases. However, anti-cancer drugs (for example Tamoxifen) blocks estrogen receptors on breast cells, that is, it stops estrogen from connecting to the cancer cells while tamoxifen also acts like an anti-estrogen in breast cells; it acts like an estrogen in other tissues like the uterus and the bones [48]. In addition, ketosis also regulating hormonal imbalance [8,27]. On the other hand, Figure 13, shows the effect of immune response with and without controls. Immune response can help to fight cancer cells while immune system recognize cancer cells as abnormal and kill them. However, this may not be enough to eliminate cancer cells from the body.

7. Conclusions

A four-dimensional compartmental deterministic model was designed and used to monitor the dynamics of breast cancer. The existing model in [19] was extended to incorporate treatments, ketogenic diet, and an immune booster. The system (5) was rigorously analyzed to gain insight into their dynamical behaviors. The study shows the following:
  • The conditions of stability of the tumor-free equilibrium (TFE) was established and the system is only local asymptotically stable (LAS) if a certain threshold quantity, known as the reproductive number, is less than unity ( R 0 < 1 ). It implies that the number of tumor cells in the body will be brought to zero if proper treatments and a ketogenic diet that can force make the threshold to a value less than unity are monitored.
  • An individual has the chance of developing breast cancer depending on the level of the immune system (s), the efficacy of the anti-cancer drug (k) and the rate at which the ketogenic diet (d) is being taken to fight tumor cells. We also found out that the presence of excess estrogen in system makes it unstable, as depicted in Figure 7. This implies that any additional estrogen quantity introduced into the body through the birth control, and hormone replacement therapy (HRT) enhances the rate of tumor formation. Thus, the development of breast cancer is certain.
  • The transition from normal cells class to tumor cells class plays a crucial role in breast cancer dynamics ( λ 1 ) . More tumor is formed if the DNA is damaged or altered as a result of excess estrogen, which reduces the number of normal cells being produced by red blood cells.
Furthermore, the results show that tumor cell formation depend on the level of excess estrogen introduced into the body system. It must be noted that the ability to resist changes in structure and amount of estrogen released during natural biological processes is dependent on an individual’s DNA. Such biological processes include: premenopausal and menopause stages. Other risk factors may also be incorporated in the model for future work, which might generate different results.
However, the focus of this study has been identifying the advantages that come with the process of breast cancer relief policies that combined anti-cancer drugs and ketogenic diet procedures to knit the circumstances of unlimited and limited resources. The effort to moderate the effect of breast cancer on the body can be fruitful, especially if our basic reproductive number R0 is properly analyzed. In addition, moderation is conceivable if the planning of intercessions is sufficiently quick and if the arrangement includes the utilization of more than one therapy procedure. No therapy (ketogenic diet and anti-cancer drug) is possible, unless minimal resources are accessible.

8. Further Research

Breast cancer is a health challenge disease, especially among women world-wide. This study explored the use of a quadratic control law to formulate the optimal control problem for the objective function. Hence, the authors hope to conduct further research into the application of a switching function and to investigate the side-effects of anti-cancer drugs by employing a linear control law to formulate the optimal control problem for further study.

Author Contributions

The authors contributed towards mathematical model formulation; Segun Isaac Oke analyzed the model both analytical and numerical simulation.

Acknowledgments

The corresponding author appreciate National Research Foundation (NRF), South Africa for the grant towards my PhD; Grant Number: 109824. The authors also acknowledges the support of Research Office of University of Zululand for providing the funds for the publication. The authors are grateful to Adeniyi Michael (LASPOTECH, Nigeria) and Chinaza Uleanya (Unizulu) for their useful comments in the preparation of the manuscript. The authors are grateful to the anonymous Reviewers and the Handling Editor for their constructive comments, which have enhanced the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Evans, C.W. The invasion and metastatic behaviour of malignant cells. In The Metastatic Cell: Behavior and Biochemistry; Chapman and Hall: London, UK, 1991; pp. 137–214. [Google Scholar]
  2. National Cancer Registry of South Africa. Available online: http://www.cansa.org.za/files/2015/10/NCR_Final_2010_tables1.pdf (accessed on 24 April 2018).
  3. World Health Organization. Global Action Plan for the Prevention and Control on NCDs; World Health Organization: Geneva, Switzerland, 2014. [Google Scholar]
  4. Abernathy, K.; Abernathy, Z.; Baxter, A.; Stevens, M. Global Dynamics of a Breast Cancer Competition Model. Differ. Equ. Dyn. Syst. 2017, 3, 1–15. [Google Scholar] [CrossRef]
  5. Patel, M.I.; Nagl, S. The Role of Model Integration in Complex Systems Modelling: An Example from Cancer Biology; Springer: Berlin, Germany, 2010. [Google Scholar]
  6. Allen, B.G.; Bhatia, S.K.; Anderson, C.M.; Eichenberger-Gilmore, J.M.; Sibenaller, Z.A.; Mapuskar, K.A.; Schoenfeld, J.D.; Buatti, J.M.; Spitz, D.R.; Fath, M.A. Ketogenic diets as an adjuvant cancer therapy: History and potential mechanism. Redox Biol. 2014, 2, 963–970. [Google Scholar] [CrossRef] [PubMed]
  7. Gilbert, D.L.; Pyzik, P.L.; Freeman, J.M. The ketogenic diet: Seizure control correlates better with serum β-hydroxybutyrate than with urine ketones. J. Child Neurol. 2000, 15, 787–790. [Google Scholar] [CrossRef] [PubMed]
  8. Westman, E.C.; Yancy, W.S.; Mavropoulos, J.C.; Marquart, M.; McDuffie, J.R. The effect of a low-carbohydrate, ketogenic diet versus a low-glycemic index diet on glycemic control in type 2 diabetes mellitus. Nutr. Metab. 2008, 5, 36. [Google Scholar] [CrossRef] [PubMed]
  9. Kareva, I.; Berezovskaya, F. Cancer immunoediting: A process driven by metabolic competition as a predator–prey–shared resource type model. J. Theor. Biol. 2015, 380, 463–472. [Google Scholar] [CrossRef] [PubMed]
  10. Lenhart, S.; Workman, J.T. Optimal Control Applied to Biological Models; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  11. Pontryagin, L.S. Mathematical Theory of Optimal Processes; CRC Press: Boca Raton, FL, USA, 1987. [Google Scholar]
  12. De Pillis, L.G.; Radunskaya, A. A mathematical tumor model with immune resistance and drug therapy: An optimal control approach. Comput. Math. Methods Med. 2001, 3, 79–100. [Google Scholar] [CrossRef]
  13. D’Onofrio, A.; Ledzewicz, U.; Maurer, H.; Schättler, H. On optimal delivery of combination therapy for tumors. Math. Biosci. 2009, 222, 13–26. [Google Scholar] [CrossRef] [PubMed]
  14. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. In Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, London, UK, 1 August 1927; The Royal Society: London, UK, 1927; Volume 115, pp. 700–721. [Google Scholar]
  15. Kermack, W.O.; McKendrick, A.G. Contributions to the mathematical theory of epidemics: V. Analysis of experimental epidemics of mouse-typhoid; a bacterial disease conferring incomplete immunity. Epidemiol. Infect. 1939, 39, 271–288. [Google Scholar] [CrossRef]
  16. Cai, P.; Tang, J.S.; Li, Z.B. Analysis and controlling of Hopf Bifurcation for chaotic Van der Pol-Duffing system. Math. Comput. Appl. 2014, 19, 184–193. [Google Scholar] [CrossRef]
  17. Bernoulli, D. Essai d’une nouvelle analyse de la mortalité causée par la petite vérole et des avantages de l’inoculation pour la prévenir. Hist. Acad. R. Sci. (Paris) Mém. Math. Phys. Mém. 1760, 1, 1–45. [Google Scholar]
  18. De Pillis, L.G.; Radunskaya, A.E.; Wiseman, C.L. A validated mathematical model of cell-mediated immune response to tumor growth. Cancer Res. 2005, 65, 7950–7958. [Google Scholar] [CrossRef] [PubMed]
  19. Mufudza, C.; Walter, S.; Chiyaka, E.T. Assessing the effects of estrogen on the dynamics of breast cancer. Comput. Math. Methods Med. 2012, 2012, 473572. [Google Scholar] [CrossRef] [PubMed]
  20. Malinzi, J.; Eladdadi, A.; Sibanda, P. Modelling the spatiotemporal dynamics of chemovirotherapy cancer treatment. J. Biol. Dyn. 2017, 11, 244–274. [Google Scholar] [CrossRef] [PubMed]
  21. Swierniak, A.; Krzeslak, M.; Student, S.; Rzeszowska-Wolny, J. Development of a population of cancer cells: Observation and modeling by a mixed spatial evolutionary games approach. J. Theor. Biol. 2016, 405, 94–103. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, C.H.; Motohashi, T.; Abdel-Rahman, H.A.; Flickinger, G.L.; Mikhail, G. Free and protein-bound plasma estradiol-17β during the menstrual cycle. J. Clin. Endocrinol. Metab. 1976, 43, 436–445. [Google Scholar] [CrossRef] [PubMed]
  23. Kumar, A.; Srivastava, P.K. Vaccination and treatment as control interventions in an infectious disease model with their cost optimization. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 334–343. [Google Scholar] [CrossRef]
  24. Kimmel, M.; Swierniak, A. Control theory approach to cancer chemotherapy: Benefiting from phase dependence and overcoming drug resistance. In Tutorials in Mathematical Biosciences III; Springer: Berlin/Heidelberg, Germany, 2006; pp. 185–221. [Google Scholar]
  25. Buonomo, B. Modeling itns usage: Optimal promotion programs versus pure voluntary adoptions. Mathematics 2015, 3, 1241–1254. [Google Scholar] [CrossRef]
  26. Pinho, S.T.R.D.; Freedman, H.I.; Nani, F. A chemotherapy model for the treatment of cancer with metastasis. Math. Comput. Model. 2002, 36, 773–803. [Google Scholar] [CrossRef]
  27. American Cancer Society. Breast Cancer; American Cancer Society: Atlanta, GA, USA, 2013. [Google Scholar]
  28. Gatenby, R.A. The potential role of transformation-induced metabolic changes in tumor-host interaction. Cancer Res. 1995, 1766, 4151–4156. [Google Scholar]
  29. Neves, A.A.; Brindle, K.M. Assessing responses to cancer therapy using molecular imaging. Biochim. Biophys. Acta 2006, 1766, 242–261. [Google Scholar] [CrossRef] [PubMed]
  30. Valayannopoulos, V.; Bajolle, F.; Arnoux, J.B.; Dubois, S.; Sannier, N.; Baussan, C.; Petit, F.; Labrune, P.; Rabier, D.; Ottolenghi, C.; et al. Successful treatment of severe cardiomyopathy in glycogen storage disease type III With d, l-3-hydroxybutyrate, ketogenic and high-protein diet. Pediatr. Res. 2011, 70, 638–641. [Google Scholar] [CrossRef] [PubMed]
  31. De Leenheer, P.; Aeyels, D. Stability properties of equilibria of classes of cooperative systems. IEEE Trans. Autom. Control 2001, 46, 1996–2001. [Google Scholar] [CrossRef]
  32. Perko, L. Differential Equations and Dynamical Systems; Springer: Berlin, Germany, 2013; Volume 7. [Google Scholar]
  33. Marino, S.; Hogue, I.B.; Ray, C.J.; Kirschner, D.E. A methodology for performing global uncertainty and sensitivity analysis in systems biology. J. Theor. Biol. 2008, 254, 178–196. [Google Scholar] [CrossRef] [PubMed]
  34. Acar, E.; Aplak, H.S. A Model Proposal for a Multi-Objective and Multi-Criteria Vehicle Assignment Problem: An Application for a Security Organization. Math. Comput. Appl. 2016, 21, 39. [Google Scholar] [CrossRef]
  35. Madhi, M.; Mohamed, N. An Initial Condition Optimization Approach for Improving the Prediction Precision of a GM (1,1) Model. Math. Comput. Appl. 2017, 22, 21. [Google Scholar] [CrossRef]
  36. Zhuang, K. Spatiotemporal Dynamics of a Delayed and Diffusive Viral Infection Model with Logistic Growth. Math. Comput. Appl. 2017, 22, 7. [Google Scholar] [CrossRef]
  37. De Pillis, L.G.; Gu, W.; Fister, K.R.; Head, T.A.; Maples, K.; Murugan, A.; Neal, T.; Yoshida, K. Chemotherapy for tumors: An analysis of the dynamics and a study of quadratic and linear optimal controls. Math. Biosci. 2007, 209, 292–315. [Google Scholar] [CrossRef] [PubMed]
  38. De Pillis, L.G.; Fister, K.R.; Gu, W.; Head, T.; Maples, K.; Neal, T.; Murugan, A.; Kozai, K. Optimal control of mixed immunotherapy and chemotherapy of tumors. J. Biol. Syst. 2008, 16, 51–80. [Google Scholar] [CrossRef]
  39. Kirschner, D.; Lenhart, S.; Serbin, S. Optimal control of the chemotherapy of HIV. J. Math. Biol. 1997, 35, 775–792. [Google Scholar] [CrossRef] [PubMed]
  40. Swan, G.W. Role of optimal control theory in cancer chemotherapy. Math. Biosci. 1990, 101, 237–284. [Google Scholar] [CrossRef]
  41. Ratajczyk, E.; Ledzewicz, U.; Schättler, H. Optimal Control for a Mathematical Model of Glioma Treatment with Oncolytic Therapy and TNF-α Inhibitors. J. Optim. Theory Appl. 2018, 456–477. [Google Scholar] [CrossRef]
  42. Di Liddo, A. Optimal Control and Treatment of Infectious Diseases. The Case of Huge Treatment Costs. Mathematics 2016, 4, 21. [Google Scholar] [CrossRef]
  43. Fleming, W.H.; Rishel, R.W. Deterministic and Stochastic Optimal Control; Springer: Berlin, Germany, 2012; Volume 1. [Google Scholar]
  44. Ghaddar, C.K. Novel Spreadsheet Direct Method for Optimal Control Problems. Math. Comput. Appl. 2018, 23, 6. [Google Scholar] [CrossRef]
  45. Lukes, D.L. Differential Equations: Classical to Controlled; Elsevier: New York, NY, USA, 1982. [Google Scholar]
  46. Schattler, H.; Ledzewicz, U. Optimal Control for Mathematical Models of Cancer Therapies; Springer: New York, NY, USA, 2015. [Google Scholar]
  47. Otieno, G.; Koske, J.K.; Mutiso, J.M. Cost effectiveness analysis of optimal malaria control strategies in kenya. Mathematics 2016, 4, 14. [Google Scholar] [CrossRef]
  48. Davies, C.; Pan, H.; Godwin, J.; Gray, R.; Arriagada, R.; Raina, V.; Abraham, M.; Medeiros Alencar, V.H.; Badran, A.; Bonfill, X.; et al. Long-term effects of continuing adjuvant tamoxifen to 10 years versus stopping at 5 years after diagnosis of oestrogen receptor-positive breast cancer: ATLAS, a randomised trial. Lancet 2013, 381, 805–816, Erratum in 2013, 381, 804. [Google Scholar] [CrossRef]
Figure 1. PRCCs of homogeneous model parameters with the tumor cells as the baseline variable. All parameter values were varied in 25% of their baseline values in Table 1. The most sensitive parameters are shown to be P v a l u e s of α 1 , g , μ 1 , γ 3 and ω are less than 0.01.
Figure 1. PRCCs of homogeneous model parameters with the tumor cells as the baseline variable. All parameter values were varied in 25% of their baseline values in Table 1. The most sensitive parameters are shown to be P v a l u e s of α 1 , g , μ 1 , γ 3 and ω are less than 0.01.
Mca 23 00021 g001
Figure 2. PRCCs of homogeneous model parameters with the tumor cells as the baseline variable. All parameter values were varied in 25% of their baseline values in Table 1. The most sensitive parameters are shown to be P v a l u e s of s , γ 2 , μ 3 and ρ are less than 0.01.
Figure 2. PRCCs of homogeneous model parameters with the tumor cells as the baseline variable. All parameter values were varied in 25% of their baseline values in Table 1. The most sensitive parameters are shown to be P v a l u e s of s , γ 2 , μ 3 and ρ are less than 0.01.
Mca 23 00021 g002
Figure 3. The variation of proportion of Tumor cell population for different values of d with other parameters fixed.
Figure 3. The variation of proportion of Tumor cell population for different values of d with other parameters fixed.
Mca 23 00021 g003
Figure 4. The variation of proportion of Estrogen level population for different values of k with other parameters fixed.
Figure 4. The variation of proportion of Estrogen level population for different values of k with other parameters fixed.
Mca 23 00021 g004
Figure 5. The variation of proportion of Tumor cell population for different values of k with other parameters fixed.
Figure 5. The variation of proportion of Tumor cell population for different values of k with other parameters fixed.
Mca 23 00021 g005
Figure 6. The variation of proportion of Immune booster population for different values of β with other parameters fixed.
Figure 6. The variation of proportion of Immune booster population for different values of β with other parameters fixed.
Mca 23 00021 g006
Figure 7. The variation of Total cells population depicted as locally asymptotically unstable.
Figure 7. The variation of Total cells population depicted as locally asymptotically unstable.
Mca 23 00021 g007
Figure 8. The variation of Total cells population depicted as locally asymptotically stable.
Figure 8. The variation of Total cells population depicted as locally asymptotically stable.
Mca 23 00021 g008
Figure 9. The variation of proportion of Normal cell population for different values of λ 1 with other parameters fixed.
Figure 9. The variation of proportion of Normal cell population for different values of λ 1 with other parameters fixed.
Mca 23 00021 g009
Figure 10. Simulation result of the model (9), showing normal cell population against time with and without control.
Figure 10. Simulation result of the model (9), showing normal cell population against time with and without control.
Mca 23 00021 g010aMca 23 00021 g010b
Figure 11. Simulation result of the model (9), showing tumor cell population against time with and without control.
Figure 11. Simulation result of the model (9), showing tumor cell population against time with and without control.
Mca 23 00021 g011
Figure 12. Simulation result of the model (9), showing estrogen level against time with and without control.
Figure 12. Simulation result of the model (9), showing estrogen level against time with and without control.
Mca 23 00021 g012
Figure 13. Simulation result of the model (9), showing immune response against time with and without control.
Figure 13. Simulation result of the model (9), showing immune response against time with and without control.
Mca 23 00021 g013aMca 23 00021 g013b
Table 1. Description of parameters in the model.
Table 1. Description of parameters in the model.
Parameter SymbolValueUnitRefs
Per capita growth rate of normal cells α 1 0.70day−1[12]
Per capita growth rate of tumor cells α 2 0.514day−1[5]
Natural death rate of normal cells μ 1 0.00003day−1Assumed
Natural death rate of tumor cells μ 2 0.01day−1[7]
Rate of inhibition of normal cells ϕ 1 6 × 10−8day−1[1]
Tumor cells death rate due to immune response γ 2 3 × 10−6day−1[12]
Interaction coefficient rate with immune response γ 3 1 × 10−7day−1[5]
Source rate of immune cells s 1.3 × 104day−1[12]
Source rate of estrogen ϵ 1.3 × 104day−1est
Immune threshold rate ω 3 × 105day−1[5]
Immune response rate ρ 0.20day−1[13]
Natural death rate of immune cells μ 3 0.29day−1[5]
Efficacy of anti-cancer drug k 0–1day−1Assumed
Supplement for immune booster β 0.01day−1est
Tumor formation rate as a result of DNA damage by excess estrogen λ 1 0.20(Pg/mL)−1day−1est
Immune suppression rate due to excess estrogen λ 3 0.002day−1est
Assume constant of value of decay factor g 0.1day−1est
Natural death rate of estrogen μ 4 0.97day−1[19]
Death rate due to ketogenic diet μ 5 2.0day−1est
Constant rate of ketogenic diet d 0.5day−1est

Share and Cite

MDPI and ACS Style

Oke, S.I.; Matadi, M.B.; Xulu, S.S. Optimal Control Analysis of a Mathematical Model for Breast Cancer. Math. Comput. Appl. 2018, 23, 21. https://doi.org/10.3390/mca23020021

AMA Style

Oke SI, Matadi MB, Xulu SS. Optimal Control Analysis of a Mathematical Model for Breast Cancer. Mathematical and Computational Applications. 2018; 23(2):21. https://doi.org/10.3390/mca23020021

Chicago/Turabian Style

Oke, Segun Isaac, Maba Boniface Matadi, and Sibusiso Southwell Xulu. 2018. "Optimal Control Analysis of a Mathematical Model for Breast Cancer" Mathematical and Computational Applications 23, no. 2: 21. https://doi.org/10.3390/mca23020021

APA Style

Oke, S. I., Matadi, M. B., & Xulu, S. S. (2018). Optimal Control Analysis of a Mathematical Model for Breast Cancer. Mathematical and Computational Applications, 23(2), 21. https://doi.org/10.3390/mca23020021

Article Metrics

Back to TopTop