Abstract

To begin with, in this paper, single immunotherapy, single chemotherapy, and mixed treatment are discussed, and sufficient conditions under which tumor cells will be eliminated ultimately are obtained. We analyze the impacts of the least effective concentration and the half-life of the drug on therapeutic results and then find that increasing the least effective concentration or extending the half-life of the drug can achieve better therapeutic effects. In addition, since most types of tumors are resistant to common chemotherapy drugs, we consider the impact of drug resistance on therapeutic results and propose a new mathematical model to explain the cause of the chemotherapeutic failure using single drug. Based on this, in the end, we explore the therapeutic effects of two-drug combination chemotherapy, as well as mixed immunotherapy with combination chemotherapy. Numerical simulations indicate that combination chemotherapy is very effective in controlling tumor growth. In comparison, mixed immunotherapy with combination chemotherapy can achieve a better treatment effect.

1. Introduction

Immunotherapies are becoming a crucial component in the multipronged approaches which are developed to treat certain cancer [1]. By strengthening the antitumor function of the immune system, immunotherapy can enhance the body’s own natural ability to combat cancer. In past decades, advances in cancer immunology have been increasingly translated into clinical testing of immune-based approaches to cancer treatment, including monoclonal antibody treatment and adoptive transfer of the cytotoxic T lymphocytes (ACT) [2, 3]. In detail, this technique [4] need identify autologous or allogeneic lymphocytes with antitumor activity, which are then infused into cancer patients, often along with appropriate growth factors (such as IL-2) to stimulate their survival. It is indispensable to identify only a small number of antitumor cells with the appropriate properties which can then be expanded to large numbers ex vivo for treatment [5]. In vitro tests can identify the exact populations and effector functions required for cancer regression, which can then be selected for expansion [6]. The cells that are activated in the laboratory can wipe off endogenous inhibitory factors and thus can be induced to exhibit the required antitumor effector functions.

Recently, clinical data have indicated that there is a potential benefit in making use of the power of the immune system in conjunction with traditional chemotherapy [7]. As a conventional treatment, chemotherapy has become a part of treatment plan of most cancer patients. Chemotherapy aims at shrinking primary tumors, slowing the tumor growth, and killing cancer cells that may have spread (metastasized) to other parts of the body from the original, primary tumor. Currently, more than 50 kinds of chemotherapy drugs are available to treat cancer and many more are being tested for their ability to destroy cancer cells [8]. Although chemotherapy is one of the principal modes of treatment for cancer patients, one of the limitations of chemotherapy is that it also kills the normal fast dividing cells, which causes serious side effects in patients.

The immune response to a tumor is usually cell mediated with cytotoxic T lymphocytes (CTL) cells and natural killer (NK) cells playing a dominant role. Mathematical modelling has become an important and useful tool in studying the interactions between the immune system and a growing tumor (Bianca et al. [9, 10], Pappalardo et al. [1113]). Bell (1973) [14] proposed a model consisting of a system of two equations based on the classic predator-prey interaction. Kuznetsov et al., in 1992 [15] and 1994 [16], presented a mathematical model of CTL cells response to the growth of immunogenic tumor, which exhibits a number of phenomena that are observed in vivo, including immunostimulation, “sneaking through,” and “dormant state” of the tumor. Moreover, the parameters of the target model were estimated by using the experimental data of chimeric mice. de Pillis, in [1, 1719], analyzed the interactions among tumor cells and various immune effector cells and applied the numerical calculations to discuss the treatment effects of different therapeutic regimens. Due to considering the treatment processes that are subject to short-term perturbations, the model with impulsive treatments conforms better to the practice than the continuous models mentioned above. Hence, Borges et al. (2014) [20] introduced continuous and pulsed chemotherapy to investigate the treatment effect of tumor with the help of numerical calculations.

In this paper, our purpose is to provide a useful reference and guidance for experimental workers and scientists of human cancer research by designing treatment protocols of chimeric mice with pulsed chemotherapy and immunotherapy. Hence, we introduce pulsed immunotherapy and chemotherapy into the mathematical model proposed by Kuznetsov, and the model readswhere denotes the concentration of CTL cells with antitumor activity in the tumor site, represents the number of tumor cells, and is the blood drug concentration. is the therapeutic period, is the infusion dose of CTL cells with antitumor activity every time, and denotes an increment of the blood drug concentration due to delivering drug at time . , , and denote the right limits of , , and at time , respectively. The descriptions and estimated values of all remaining parameters (which were estimated by using the experimental data of chimeric mice [16]) are listed in Table 1. For convenience, we suppose that immunotherapy and chemotherapy are executed at almost the same time and use to denote the therapeutic regimen, where , , and represent the dosage of immunotherapy, an increment of the blood drug concentration caused by chemotherapy and therapeutic period, respectively.

Although drug targeted therapy is yielding promising results in the treatment of some specific cancers, drug resistance caused mainly by mutation plays a critical role of the chemotherapy failure [21]. Hence, incorporating drug resistance into our model can help us to find ways to eliminate the problem. The paper is organized in the following manner. In Section 2, we will investigate therapeutic effects of immunotherapy, chemotherapy, and mixed treatment and design the corresponding therapeutical schedules. In Section 3, the efficacy of cancer chemotherapy often becomes severely limited as cancer cells become resistant to chemotherapy drugs. Hence, we will introduce the drug resistance into system (1) to account for the failure of chemotherapy and develop some new therapeutic regimens so as to achieve the goal of clinical cure. This paper ends with a brief conclusion and discussion in Section 4.

2. Investigation of Therapeutic Regimens

In this section, we discuss the effects of single immunotherapy, single chemotherapy, and mixed immunotherapy with chemotherapy and provide the corresponding therapeutic regimens. First of all, we discuss the single immunotherapy.

2.1. Single Immunotherapy

Suppose that system (1) only involves immunotherapy (i.e., and ), which is equivalent to the following system:Firstly, we give some basic properties about the following subsystem of (2):Clearly,is a positive periodic solution of system (3), where Since the solution of system (3) with any initial value iswe have the following lemma.

Lemma 1. System (3) has a positive periodic solution . And, for every solution of (3), it follows that . Furtherly, system (2) has a tumor-free periodic solution .

Next, we discuss the local stability of the tumor-free periodic solution .

Theorem 2. Let be any solution of system (2); then is locally asymptotically stable provided that

Proof. The local stability of the periodic solution can be determined by the small amplitude perturbations of the solution. Define where are small perturbations. The linearized equations of system (2) are given byLet be fundamental matrix of (9) and is the identity matrix; then must satisfyand . Hence, the fundamental solution matrix is where and the exact form of is not required in this analysis.
When , the linearization of the resetting impulsive conditions of (2) becomesThus, the monodromy matrix of (2) is Let and be eigenvalues of matrix , and then Therefore, all eigenvalues of , namely, and , have absolute values less than one if and only if (7) holds. According to Floquet theory of impulsive differential equations, the tumor-free periodic solution is locally asymptotically stable.

Remark 3. The tumor-free periodic solution is unstable if .

Theorem 4. A supercritical bifurcation occurs at in the sense that there is such that for all there is a stable positive nontrivial periodic solution of (2) with period .

Detailed proof is similar to Theorem  1 in [22].

From (7), we know that when the infusion dose of CTL cells with antitumor activity every time is not less than the threshold value , then the tumor will be eliminated eventually; otherwise the number of tumor cells will present periodic oscillation. Hence, we consider a tumor with size and take therapeutic period days. From (7) and Table 1, we calculate and achieve a critical value CTL cells. Here, we take (more than ) and (less than ), respectively. The variations of the number of tumor cells with time are depicted in Figures 1(a) and 1(b). Figure 1(b) indicates that since , the number of tumor cells exhibits periodic oscillation. However, from Figure 1(a), we can see that although the tumor cells are wiped out eventually, this process takes too long time of about 37 years (13000 days). In order to shorten the time of curing a tumor, we increase the infusion dosage of CTL cells. Supposing the time of curing a cancer is a year, by numerical calculations, we achieve a new threshold value CTL cells. Hence, we obtain a single immunotherapy regimen:The corresponding change of the number of tumor cells is shown in Figure 1(c) when is carried out.

Furthermore, by Theorem 2, we derive a condition which is given by where straight line is defined as “critical boundary of treatment regimens.” The region which lies above is named as “acceptable region” (i.e., the region of successful treatment), and the region which sits below is intituled as “rejected region” (i.e., the region of failed treatment), which is shown in Figure 1(d).

2.2. Single Chemotherapy

We consider the case of single chemotherapy (i.e., for system (1)). Since the third equation in system (1) is independent of the variables and , we consider the following subsystem:Clearly, system (18) has a positive periodic solution:Since the solution of (18) with initial value is we have .

Furthermore, from (19), we have

The dose-delivery schedule of chemotherapy drugs can determinate their efficacy in killing cancer cells and degree of toxicity to the patients [23]. Besides, conventional chemotherapy drugs often have a therapeutic window which is defined as a range of a drug’s serum concentration at which a desired effect occurs, below which there is little effect, and above which toxicity occurs [24]. Hence, we denote the least effective concentration (LEC) as and the maximum tolerated concentration (MTC) as , respectively. Thus, we have . Without loss of generality, we suppose and . Hence, we can getwhich is equivalent toFurthermore, we know that Adriamycin is a drug used in cancer chemotherapy and is commonly used in the treatment of a wide range of cancers, including hematological malignancies (blood cancers, like leukaemia and lymphoma), many types of carcinoma (solid tumors), and soft tissue sarcomas. In our research, we select Adriamycin as a chemotherapy drug and follow the dosage suggested by the manufactures of the Adriamycin drug [19]. The suggested procedure entails a single dose of 60–75  once every 21 days [25]. Further, we approximate mice to have surface area of and have body volume of and use the dosing as every 21 days. Supposing that Adriamycin is distributed uniformly in every tissues, we have Since Adriamycin has a half-life of about 48 hours, from the first equation of (18), we getand then have

Further, from (23), we can obtain

According to dosage suggested by the manufactures of the Adriamycin drug, we can get a single chemotherapeutic regimen:Taking as a chemotherapy regimen, the dynamic behavior of tumor cells is depicted in Figure 2(a) which indicates that the number of tumor cells exhibits periodic oscillation. In other words, single chemotherapeutic regimen is not adequate to wipe out a tumor eventually. Since high concentration of chemotherapy drug can kill more tumor cells, we can try to elevate the average blood drug concentration of Adriamycin by increasing the least effective concentration of therapeutic window so that chemotherapy can obtain a better treatment result. In order to facilitate the following discussion, we denote new least effective concentration, new infusion dosage of Adriamycin, and new chemotherapy period by , , and , respectively. From (23), we haveFor different values of , from (29), we achieve different chemotherapy regimens which are listed in Table 2.

From Table 2, we clearly see that higher least effective concentration requires smaller infusion dosage of chemotherapy drug and shorter therapeutic period. By numeric calculations, we obtain a critical value of the least effective concentration , which implies that if , then the corresponding chemotherapy regimen is successful; otherwise chemotherapy regime is failing. Since and are more than and and are smaller than , performing treatment regimens and can effectively control tumor growth, but executing therapeutic regimes and makes the number of tumor cells present periodical oscillation (see Figure 2(b)).

Supposing that the time of curing a cancer is still a year, we can obtain that the least effective concentration is not lower than . Hence, a single chemotherapy treatment regimen is given by (see Figure 2(b))

Further, we consider that the half-life of chemotherapy drug impacts on therapeutic result. Equation (25) indicates that extending the half-life of drug makes decay rate of chemotherapy drug reduce so that the average concentration of drug goes up. Thus, supposing that the half-life of Adriamycin is , , and , from (25), we obtain that the decay rates of it () are , , and , respectively. In addition, we still take (i.e., recommended dosage suggested by the manufactures of the Adriamycin drug) as chemotherapy regimen. The variations of the number of tumor cells with time are depicted in Figure 3(a) when the half-lives of Adriamcin are , , and , respectively. The changes of the concentration of drugs with time are shown in Figure 3(b) when the half-lives of them are and , respectively. Figures 3(a) and 3(b) show that the half-life of drug is longer and the average concentration of drug is higher so that the effect of chemotherapy is better.

2.3. Mixed Immunotherapy with Chemotherapy

According to recommended dosage suggested by the manufactures of the Adriamycin drug, we know that the treatment regimen is not sufficient to wipe out a tumor with size alone. In addition, increasing the dosage of Adriamycin may cause grievous damage to cancer patients. However, in contrast with chemotherapy, the immunotherapy can differentiate between normal and malignant cells and, thus, decrease the damage to normal cells or tissues due to chemotherapy. Hence, we consider mixed immunotherapy with chemotherapy and aim at designing a mixed treatment schedule which can prevent the tumor growth more effectively. In detail, we still take suggested dosage of the Adriamycin drug and try to seek out a critical dosage of CTL cells , so that the tumor will die out eventually when . With the help of computer, we get CTL cells. As a result, we achieve a mixed treatment regimen:The dynamic behavior of tumor cells is shown in Figure 4(a) for therapeutic regimen . Figure 4(a) shows that although taking treatment regimen can eliminate the tumor eventually, this process also takes too long time (up to 100 years). Obviously, is not a good treatment means to eliminate a tumor with size . In order to shorten the time of curing a cancer to a year, we need add infusion dosage of CTL cells. By numerical simulations, we achieve a new critical value CTL cells. In other words, as long as infusion dosage of CTL cells is not less than , then a tumor with size can be cured within a year (see Figure 4(b)). Thus, a new mixed treatment regimen is given byFurther, we also find that mixed immunotherapy and chemotherapy can make the number of tumor cells always in a lower level, which indicates that mixed treatment is evidently superior to single chemotherapy or single immunotherapy.

3. Considering the Cases of Drug Resistance

Although chemotherapies are effective treatment for metastatic tumors, the ability of cancer cells to become resistant to chemotherapy drugs remains a significant impediment to successful chemotherapy [26]. Drug resistance results from a variety of factors including individual variations in patients and somatic cell genetic differences in tumors, even those from the same tissue of origin. The most common reason for acquisition of resistance to a broad range of anticancer drugs is expression of one or more energy-dependent transporters that detect and eject anticancer drugs from cells, but other mechanisms of resistance including insensitivity to drug-induced apoptosis and induction of drug-detoxifying mechanisms probably play an important role in acquired anticancer drug resistance [27]. Hence, in this section, we will investigate the impact of drug resistance on therapeutic results and start a discussion about chemotherapy with a single drug at first.

3.1. Chemotherapy with a Single Drug

We introduce drug resistance into model (1) to explain the cause of chemotherapy failure with single drug. New model is given bywhere denotes the number of tumor cells which are sensitive to Adriamycin. presents the number of tumor cells which are resistant to Adriamycin. is the conversion rate of tumor cells from being sensitive to being resistant to Adriamycin. Usually is very small since cancer cells mutate at a rate of about in every cells; that is, [28].

As mentioned before, taking single chemotherapy regimen can wipe out a tumor with size within a year without consideration of drug resistance. Now, we still take as a treatment regimen. By numerical simulations, the variations of the numbers of sensitive tumor cells and resistant tumor cells are shown in Figures 5(a) and 5(b), respectively. From Figure 5(a), we know that the number of sensitive tumor cells to Adriamycin reaches quickly. However, Figure 5(b) indicates that the number of resistant tumor cells to Adriamycin is stable at a fixed value (about ). All these findings show that the sensitive tumor cells will completely convert into resistant tumor cells as chemotherapy is executed, which implies that the drug resistance is the main cause of chemotherapy failure. In order to solve the problem of drug resistance, next, we will investigate the treatment effect of combination chemotherapy with two drugs.

3.2. Combination Chemotherapy with Two Drugs

As analyzed above, we know that the designs of cancer chemotherapeutic regimens have become increasingly sophisticated, and a single chemotherapy drug is very difficult to cure a tumor. The use of multiple therapeutic agents in combination has become the primary strategy to treat drug resistant cancers. This approach is called combination chemotherapy which provides a higher chance of destroying cancer cells. As a consequence, chemotherapy with two or more cytotoxic drugs that kill tumor cells by one or more mechanisms will be considered. To make the discussion easier, we assume that two drugs are administered to treat a tumor with drug resistance.

For convenience, we denote two chemotherapy drugs by and which are not toxic for the same normal organ. The concentrations of them are denoted by and . and are the decay rates of them, respectively. The tumor cells population is divided into three subpopulations which are sensitive to drug and drug , sensitive to drug but resistant to drug , and sensitive to drug but resistant to drug , respectively. The numbers of them are denoted by , , and . Model (33) can be modified in the following way to account for combination chemotherapy: where and are the fractional immune cells killed by drugs and . and are the fractional tumor cells killed by drugs and . and are the conversion rates of tumor cells from being sensitive to being resistant to drugs and . As mentioned above, and represent increments of the blood drug concentrations caused by combination chemotherapy every time. Analogously, we use to denote the combination chemotherapy regimen.

To facilitate discussion, we suppose that dynamical features and therapeutic effects of drugs and are the same as Adriamycin chemotherapy drug. Hence, we have , , and . The conversion rates of tumor cells from being sensitive to being resistant and are still taken as the values . The meanings and estimated values of all remaining parameters are the same as model (1). Next, in order to investigate the impact of the difference of two-drug dosages on therapeutic results, we consider two different modes, one case is , and the other is .

Case 1 (). To make tumor cells die out completely within a year and the concentrations of drugs be in a therapeutic window simultaneously, by calculating numerically, we achieve a new critical value of the least effective concentration, . That is, the tumor cells will be eliminated within a year if the least effective concentration is not less than . Hence, by (23), we obtain a combination chemotherapy regimen:Comparing the therapeutic regimen with (see (30)), it is not difficult to find that curing a tumor with drug resistance is required to administrate more dosage of chemotherapy drug, which indicates that drug resistance is one of the difficult questions of tumor treatment.

When shortening therapeutic period to 10 days and requiring to cure a tumor within a year, by calculating, then we get another combination chemotherapy regimen:Executing treatment regimen , the variation of the number of tumor cells with time is shown in Figure 6(a), which shows that the number of tumor cells is always in a lower level and reaches 0 eventually. In the following, we will discuss another case.

Case 2 (). Here, we take delivering dosages of drugs and as and every 10 days, respectively. Thus, we obtain another combination treatment regimen:The dynamic behavior of tumor cells for treatment regimen is exhibited in Figure 6(b). However, Figure 6(b) shows that the number of tumor cells is in a lower level at the early stage but later increases gradually. In other words, taking the therapeutic regimen easily causes the tumor recurrence.

Cases 1 and 2 indicate that therapeutic effect of is obviously superior to that of although the total amount of delivering drug of is the same as that of . Hence, a topic of how to determine the delivering dosage of every chemotherapy drug deserves deep exploration. The analysis result indicates that combination chemotherapy with multidrugs is very effective in controlling tumor growth, but larger dosage of chemotherapy drugs is destined to bring greater harm to normal tissues. As a result, in order to diminish the damage of chemotherapy to normal tissues, we will combine immunotherapy with combination chemotherapy.

3.3. Mixed Immunotherapy with Combination Chemotherapy

We substitute the seventh equation of model (34) by to analyze a new mixed immunotherapy with combination chemotherapy model. For convenience, we denote new therapeutic regimen as a quad , where the meanings of , , and are the same as above. Since two chemotherapy drugs and are all toxic to normal cells, we assume delivering dosage is of recommend dosage of Adriamycin. Thus, we have . Hence, we obtain a combination chemotherapy regimen: When treatment regimen is performed, the variation of the number of tumor cells with time is exhibited in Figure 7(a) which shows that the number of tumor cells presents periodical oscillation. In other words, executing can not cure a tumor successfully. Based on this, we introduce immunotherapy into treatment regimen . By numerical calculation, we achieve a threshold value CTL cells. In detail, tumor will be eliminated within a year if infusion dosage of CTL cells is not less than . Thus, a mixed immunotherapy with combination chemotherapy regimen is given byThe dynamic behavior of tumor cells is shown in Figure 7(b) when the mixed treatment regimen is carried out. From Figure 7(b), we can obviously see that a tumor with size can be eliminated very quickly.

4. Conclusion

In this paper, we investigate the therapeutic effects of single immunotherapy, single chemotherapy, and mixed treatment and provide corresponding therapeutic regimens. For single immunotherapy, we derive a condition under which tumor cells will be eliminated ultimately, but this process may take too long time. Thus, we suppose that the time of curing a tumor is a year and then achieve a single immunotherapy regimen . For single chemotherapy, we select Adriamycin as chemotherapy drug. By exploring the impacts of the least effective concentration and the drug half-life on therapeutic results, we draw a conclusion that increasing the least effective concentration and extending the half-life of the drug can make average drug concentration maintain a higher level so that the effect of chemotherapy is better. Using numerical calculations, we obtain a threshold value ; that is, single chemotherapy can eliminate a tumor within a year if the least effective concentration is not lower than . As a result, we obtain an ideal single chemotherapeutic regimen . Further, in order to decrease the damage of chemotherapy drug to normal tissues, we consider mixed immunotherapy with chemotherapy. Taking recommended dosage of Adriamycin and supposing the time of curing a tumor with size as a year, by numerical calculations, we achieve a mixed immunotherapy with immunotherapy regimen . In addition, we find that mixed treatment can make the number of tumor cells be always in a lower level.

Since most tumors are resistant to chemotherapy drugs, we consider influence of drug resistance on therapeutic results and improve a new mathematical model. By analyzing the target model, we explain the cause of chemotherapeutic failure. Further, we consider the effect of combination chemotherapy with two drugs and form a combination chemotherapy regimen which can get rid of a tumor within a year. Finally, in order to cure a tumor more effectively, we consider mixed immunotherapy with combination chemotherapy to treat a tumor and establish an ideal mixed treatment regimen which can make tumor cells be always in a lower level and be also wiped out completely in a year.

In a word, combination chemotherapy is very effective in controlling tumor growth, and further mixed immunotherapy with combination chemotherapy can obtain a better treatment effect. But, with tumor cells becoming resistant to many structurally and mechanistically unrelated drugs, the efficacy of chemotherapy of tumor often becomes severely limited. Hence, the problem of how to combine reasonably those treatment modes and design an optimal mixed therapeutic regimen deserves deep research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publishing of this paper.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (no. 11371164), NSFC-Talent Training Fund of Henan (no. U1304104), and Key Scientific Research Project of Higher Education Institutions of Henan Province (16A110005).