Research Article | Open Access
Modeling the Treatment of Glioblastoma Multiforme and Cancer Stem Cells with Ordinary Differential Equations
Despite improvements in cancer therapy and treatments, tumor recurrence is a common event in cancer patients. One explanation of recurrence is that cancer therapy focuses on treatment of tumor cells and does not eradicate cancer stem cells (CSCs). CSCs are postulated to behave similar to normal stem cells in that their role is to maintain homeostasis. That is, when the population of tumor cells is reduced or depleted by treatment, CSCs will repopulate the tumor, causing recurrence. In this paper, we study the application of the CSC Hypothesis to the treatment of glioblastoma multiforme by immunotherapy. We extend the work of Kogan et al. (2008) to incorporate the dynamics of CSCs, prove the existence of a recurrence state, and provide an analysis of possible cancerous states and their dependence on treatment levels.
Dynamical systems continue to play an important role in understanding cancer dynamics [1–4]. A recent development in cancer dynamics is the cancer stem cell hypothesis. The cancer stem cell hypothesis states that malignant tumors are initiated and maintained by a population of tumor cells that share similar biologic properties to normal adult stem cells . With evidence mounting in support of the cancer stem cell hypothesis, recent work has been devoted to the inclusion of cancer stem cells (CSCs) in current cancer models [6–9].
Cancer stem cells are a specialized type of cancer cell that are believed to be responsible for populating tumors. CSCs have a very small population in comparison to normal cancer cells because of their specialized function. While tumor cells are only able to undergo a limited number of divisions, CSCs are able to repopulate a depleted tumor, even if there are only a few CSCs left . Once the number of CSCs begins to drop, usually due to treatment, they cease creating cancer cells and focus on repopulating themselves. The small population of CSCs is hard to detect and therefore treatment is often stopped before all CSCs have been eradicated, which leads to recurrence of cancer . It is clear that, for treatment to be effective, we must focus our efforts on eliminating both tumor cells and CSCs.
The cancer stem cell hypothesis has been biologically verified for many solid tumors, including brain cancer . Glioblastoma is a type of brain tumor which forms in the cerebral hemisphere of the brain, often in the frontal and temporal lobe. These types of tumors are highly malignant, forming from normal brain cells, astrocytes, or star-shaped glial cells which support nerve cells. These cells can grow rapidly due to ample amounts of blood available in the brain. Immunotherapy is a cancer treatment which stimulates the immune system to work harder to attack cancer cells. The therapy uses additive components, such as man-made proteins, vaccines, or white blood cells, to further attack cancer cells. Immunotherapy is essential to treating multiforme glioblastomas because of their sensitive location, the brain, which is too delicate to be treated by chemotherapy or surgery [10, 13].
2. Presentation of the Model
In this paper, we extend the previous work of Kogan et al.  to include cancer stem cells in modeling the treatment of glioblastoma multiforme with immunotherapy. We present an abstract model that can be adapted to fit various biological assumptions. We analyze the stability of the model both with and without treatment and derive sufficient conditions on treatment to ensure a globally asymptotically stable cure state. We conclude with an example illustrating the transition from coexistence of cancer cells to eradication of cancer cells with various treatment levels. Much of this model is based on experimental results obtained by Kruse et al. .
The following system models the dynamics of tumor cells (), cancer stem cells (), alloreactive cytotoxic-T-lymphocytes (), TGF- (), IFN- (), and major histocompatibility complex classes I and II ( and , resp.):
To understand the formation of the system above, we discuss the biological interpretations of each equation in the model:
The first term on the right hand side (RHS) of the equation represents differentiated tumor cells produced by the CSCs without immune intervention, where is the rate at which CSCs produce tumor cells and is the number of nonstem tumor cells (TCs) currently present. The second term stands for normal tumor growth, the cells produced by regular reproduction of nonstem tumor cells. Both the first and second terms use classical logistic growth (note that the carrying capacities for and are distinct). The third term on the RHS represents tumor elimination by CTL in proportion to both and CTLs are white blood cells responsible for attacking tumor cells, in this case. The third term also introduces the effects of MHC class I receptors () and TGF- (), which is assumed to be a major immunosuppressive factor for CTL activity . Consider
The first term on the RHS stands for the rate of stem cell growth without immune intervention. As before, this follows a logistical growth model with a carrying capacity, the maximal tumor cell burden. The second term, like the first of the previous equation, stands for differentiated tumor cells produced by CSCs. The third term is almost identical to the third term of the previous equation, except that the functions and represent the interaction of CSCs and CTLs with regard to TGF- and MHC class I (as opposed to TCs in the above equation) . Consider
The first summand of the RHS stands for CTL recruitment from the blood system. The recruitment function is positively affected by MHC class II () and the number of TCs () and CSCs (). The cytokine TGF- suppresses the proliferation and activation of T-lymphocytes , as well as leukocyte migration across the brain-blood boundary (BBB) ; these are collectively represented by the function . We assume a constant death rate for , represented by . The term describes the rate of infusion of primed CTLs directly to the tumor site. is set equal to 0 in absence of immunotherapy [14, 16]. Consider
The above two equations describe the dynamics of TGF- and IFN-, respectively. In the first equation, the first term represents the natural basal level in the CNS (central nervous system). This includes TGF- produced by the tumor, which is assumed to be proportional to the tumor’s size. The degradation of TGF- is assumed to be constant with the rate and is represented by the second term.
In the second equation, the first term on the RHS is a linear production of IFN-, . We assume the only source of IFN- is CTL. The second term is the natural degradation of IFN- with constant rate . Consider
The above two equations represent the dynamics of MHC classes I and II, respectively. For the first equation, the first term on the RHS is the basal rate of receptor expression per tumor cell. This includes the stimulation by IFN- of expression on the surface of a glioblastoma cell. The second term is the natural degradation of with constant rate .
3. Preliminary Results
Throughout the paper, we will assume that system (1) is subject to nonnegative initial conditions. In addition, the functions , , , , , , , , , , , , , , , and are all functions with nonnegative values. Here, we use to denote the space of continuously differentiable functions. These assumptions imply that the nonnegative orthant is invariant under (1) and there exists a unique solution to (1) subject to initial conditions. To ensure solutions to (1) stay bounded over time, we need additional assumptions. We make the following mathematical assumptions modified from  to account for cancer stem cells (A1):(1) and are at most linear;(2) is increasing;(3) and are decreasing and bounded below; and are increasing and bounded above;(4) and are decreasing and bounded below;(5) is increasing and bounded above; is decreasing and bounded below;(6) is nonnegative and bounded above;(7), , and are increasing;(8) is increasing and bounded above;(9) is increasing and bounded above; is decreasing and bounded below.
We will use the following substitutions to simplify our equations:
These substitutions give us the following system equivalent to (1):
Also as in  and included here for the reader’s benefit, we make the following biological assumptions on our system (A2):(1) and are decreasing on [0, ], [0, ], respectively, where and are their respective carrying capacities. Furthermore, and . Also, and .(2) (when the CSC population is small, CSCs repopulate themselves, but as their population grows, they focus on populating the TC population) .(3) and are decreasing, , and , , for some (TGF- decreases the efficacy of tumors killed by CTLs up to some limit).(4) are increasing, , and , (MHC class I receptors are necessary for the action of CTLs and increase their efficiency up to some limit).(5) (CTLs are less efficient when attacking CSCs) .(6) and are decreasing, , and , (large tumor mass hampers the access of CTLs to tumor cells and reduces their kill rate).(7) is increasing from 0 to , , and (the total number of MHC class II receptors on all tumor cells and stem cells determines the recruitment of CTLs; the rate of CTL recruitment is limited and its growth decreases to zero).(8) is decreasing from 1 to some bound greater than 0 (TGF- reduces recruitment of CTLs).(9), , where (TGF- and IFN- are secreted by the cancerous cells and CTLs, resp., at constant rate; there is base level secretion of TGF-).(10) and (there is a constant basic production of MHC class II receptors at the cell surface, while IFN- increases this production up to some level).(11) and is decreasing to 0 (TGF- decreases MHC class II production to 0).(12), is increasing to , , and (IFN- is necessary to induce production of MHC class II receptors and increases up to some level with increase declining to zero).(13) (treatment is assumed to be constant with respect to time).
Theorem 1. Under the (A2) assumptions, system (1) is dissipative on .
Proof. Let , for all values of and , , where , , and These values are well defined based on the (A2) assumptions.
In order to show that our system is dissipative, we need to construct a such that where is the RHS of system (1). Let . Then, since and , for some on Similarly, for some on Thus, where and
Since our inequalities , hold for all values of , in this space, the system is dissipative everywhere on , and by a theorem from Robinson , we get the following corollary.
Corollary 2. System (1) has a compact global attractor on .
4. Stability Analysis
In the following section, we present an analysis of three potential steady states: tumor elimination (where CSC and TC populations are eradicated), recurrence (where the TC population is eradicated, but the CSC population persists), and coexistence (where CSC and TC populations persist). In each case, we discuss sufficiency conditions on the treatment term which will allow for a globally asymptotically stable cure state (tumor elimination).
4.1. Semitrivial Tumor Elimination
We begin our analysis with tumor elimination. Setting and , we find the equilibrium values: Substituting this equilibrium point into the Jacobian matrixyields eigenvalues , , , , , , and . So as long as we choose large enough so that we are guaranteed to have a locally asymptotically stable cure state.
We now show that, under necessary condition (13), is locally asymptotically stable. Let , where meets condition (13). Then, for some initial conditions, there exists such that, for , for arbitrarily small . If we also let and , for arbitrarily small, we have that from some starting moment, and therefore . By our assumptions (A1), and , so if we let be large enough so that , then we have that A parallel argument works to show that, for large enough, Thus, by increasing the treatment value , we are able to show that, for some set of initial conditions, the CSC population and TC population decay exponentially to 0, leading to tumor elimination.
4.2. Persistence of Tumor
We now wish to consider steady states in which some subset of cancer cell populations persist. Let our hypothetical equilibrium point be , where Then We wish to show existence of , , and . To do so, we must solve the system Defining the auxiliary function we see that, for every , .
In addition, taking the derivative with respect to , we get From assumptions (A2), we have that is decreasing, so there is exactly one positive for which for any given . Thus, such exists, but to further solve system (17), we need more information about the arbitrary functions present.
To move forward in the analysis of system (8), we will need to make the following simplifying assumptions (A3):(1)The dynamics of TGF- are much faster than those of the other system components.(2)The inflow of CTLs is constant.
With these assumptions, we can assume is determined by its steady-state and we get the simplified system:where Note that if , these equations will simply be referred to in terms of and vice versa.
With simplifying assumptions (A3), we are able to study the possible dynamics of persistence of cancer, recurrence, and coexistence, in the following subsections.
4.3. Recurrence State Stability
We know that, for system (20), the equilibrium points for , , and must be , , and . In this section, we wish to study the recurrence steady state, so we will set the TC population and observe the consequences for the CSC population steady state: will be a steady state when or For further analysis, we denote where . Note that (23) now becomes where is increasing in at least linearly and and are both decreasing. We recall from assumptions (A2) that and . We define for which we know for any .
We also know that . Since is an increasing function, its inverse exists and we can define . Notice that, for , . For , , and . Therefore, for , has at least one solution and, in general, since the function must cross the -axis an odd number of times, there are an odd number of solutions, .
Without loss of generality, we can assume , , and are at their respective steady states since these variables will all converge to their steady-state values exponentially. We can also assume that, for large enough , (20) is arbitrarily well approximated by
For , the equilibrium point is unstable, since any values of less than 0 will yield a negative value for , and any values of between 0 and will yield a positive value for . Thus our first stable equilibrium point is at .
In contrast, if is large enough so that , our equilibrium point is locally stable, since, for , . There could, however, still exist positive solutions , where is even. If these solutions are organized in nondecreasing order, is locally stable for even and unstable for odd , since between an odd and even root and between an even and odd root. Note that if is the only equilibrium point, it is globally asymptotically stable.
We now show that if we increase our treatment term , we can guarantee the existence of a globally asymptotically stable cure state. Let be the maximum of and choose such that Notice that this is possible since is increasing and continuous. Then for , for all , and so is now a globally asymptotically stable equilibrium point.
4.4. Coexistence State Stability
Still working under simplifying assumptions (A3), we conclude our stability analysis with consideration of a coexistence steady state. As above, for large values of , we can expect , , and to be at steady state, and so we can reduce our system to the two equations where , , , and . Define , where and .
Proposition 3. For , system (28) has a locally stable coexistence steady state.
Proof. We begin by showing the existence of such a steady state. Without loss of generality, suppose . Then for values of , and for all values of Therefore, there exists a value such that for any In general, there exist an odd number of solutions such that for any , .
Likewise, by our choice of , and for all values of . Therefore, for , there exists a value such that for all and, in general, there exist an odd number of solutions such that for any , .
Thus, we are able to deduce the existence of an equilibrium point in for system (28). Moreover, notice that while is an equilibrium point of system (28), it is unstable since for , , and for , In fact, if we denote and , equilibrium points , , , are locally unstable when or is even and locally stable when and are odd.
Remark 4. Note that, in the absence of treatment, . Therefore, in the case where the tumor is left untreated, cancer persists.
We now show that, by increasing the treatment term , we will achieve a globally asymptotically stable cure steady state. For , let us consider the system Define as the maximum of and choose such that . Similarly, let be the maximum of and choose such that Let .
Proposition 5. For , is a globally asymptotically stable equilibrium point of system (29).
Proof. For ,for since is decreasing based on assumptions (A2). Thus, for all Similarly, for , for all Therefore, is the only equilibrium point for system (29) and is globally asymptotically stable.
We conclude our analysis by noting that for all
Corollary 6. For , is a globally asymptotically stable equilibrium solution to system (28).
In this section, we present an example to illustrate the theory presented above. This model is a modification of the biologically verified system presented in Kronik et al.  to address the CSC Hypothesis as presented in this paper. Modifications include the incorporation of CSCs and amending the CTL equation to satisfy assumptions (A3). Consider the following system: subject to the initial conditions , , , , , , and . Parameter values are given by Table 1; calculations for these parameter values can be found in [8, 10, 20].
The equilibrium values for , , , , and are
Following calculations from Section 4.4, we get and . This gives
Thus, in the case of no treatment, cancer will persist (see Figure 1).
When we increase treatment past , we are able to find initial conditions that will allow a locally stable cure or recurrence state. This can be seen by adding a treatment value of and increasing our initial amount of CTLs to (see Figure 2).
When we maximize on the interval , we calculate Similarly, when we maximize on the interval , we calculate Since and , we have
Taking the maximum of these two values, we find (see Figure 3).
In this paper we extend a previous model for the treatment of glioblastoma multiforme with immunotherapy by accounting for the existence of cancer stem cells that can lead to the recurrence of cancer when not treated to completion. We prove existence of a coexistence steady state (one where both tumor and cancer stem cells survive treatment), a recurrence steady state (one where cancer stem cells survive treatment, but tumor cells do not; hence, upon discontinuation of treatment, the tumor would be repopulated by the surviving cancer stem cells), and a cure state (one where both tumor and cancer stem cells are eradicated by treatment). Furthermore, we categorize the stability of the previously men