Computational and Mathematical Methods in Medicine
Volume 2017 (2017), Article ID 3573082, 9 pages
https://doi.org/10.1155/2017/3573082

Research Article

Dynamics of Immune Checkpoints, Immune System, and BCG in the Treatment of Superficial Bladder Cancer

Farouk Tijjani Saad, Evren Hincal, and Bilgen Kaymakamzade

Department of Mathematics, Near East University, North Nicosia, Northern Cyprus, Mersin 10, Turkey

Correspondence should be addressed to Farouk Tijjani Saad; farouksaaad@yahoo.co.uk

Received 8 June 2017; Revised 29 August 2017; Accepted 26 September 2017; Published 26 October 2017

Academic Editor: Nadia A. Chuzhanova

Copyright © 2017 Farouk Tijjani Saad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper aims to study the dynamics of immune suppressors/checkpoints, immune system, and BCG in the treatment of superficial bladder cancer. Programmed cell death protein-1 (PD-1), cytotoxic T-lymphocyte-associated antigen 4 (CTLA4), and transforming growth factor-beta (TGF-β) are some of the examples of immune suppressors/checkpoints. They are responsible for deactivating the immune system and enhancing immunological tolerance. Moreover, they categorically downregulate and suppress the immune system by preventing and blocking the activation of T-cells, which in turn decreases autoimmunity and enhances self-tolerance. In cancer immunotherapy, the immune checkpoints/suppressors prevent and block the immune cells from attacking, spreading, and killing the cancer cells, which leads to cancer growth and development. We formulate a mathematical model that studies three possible dynamics of the treatment and establish the effects of the immune checkpoints on the immune system and the treatment at large. Although the effect cannot be seen explicitly in the analysis of the model, we show it by numerical simulations.

1. Introduction

Cancer is a class of diseases characterized by out-of-control cell growth which affects and damages the DNA. Cancer prevalence is increasing in many countries [1]. Many treatment options of cancer exist, which include surgery, immunotherapy, chemotherapy, radiotherapy, vaccine therapy, and hormonal therapy [1, 2]. The mode and type of treatment depend on the type, location, and grade of the cancer and the patient’s body. The bladder is a hollow organ in the lower abdomen which collects urine produced by the kidneys. Bladder cancer is a growth of malignant cells initiating in the urinary bladder. It is common, with around 38,000 men and 15,000 women diagnosed every year in the United States. Approximately 400,000 new cases are diagnosed and about 150,000 die directly from the disease every year across the globe [3, 4].

The bladder wall is lined with cells called transitional and squamous cells. The most common type of bladder cancer is urothelial carcinoma or transitional cell carcinoma (TCC). It mostly originates from the transitional cells and further progresses and grows superficially on the inner surface of the bladder; as a result, it invades the bladder wall and vessels, dispersing into the neighboring organs as well as forming distant metastases [5–7].

One of the most effective ways of treating bladder cancer is immunotherapy. This is the process of stimulating, activating, and triggering the immune system to spread, locate, and kill cancer cells [8].

Intravesical Bacillus Calmette-Guerin (BCG) is an attenuated nonpathogenic strain of Mycobacterium bovis that was initially used as a vaccine against tuberculosis. The attenuation was reached via manipulation of the bacillus by serial growth on a culture medium. As a result, the genes causing virulence will be lost and inoculated into humans [9, 10]. It is undoubtedly the most efficient and successful immunotherapy of cancer [9]. BCG therapy is used for various types of cancers, including acute lymphoblastic leukemia and melanoma. The first report of successful use of BCG to treat patients with bladder cancer was in 1976 by Morales et al. They obtained the efficacy of BCG therapy and established it as the pillar for the treatment of non-muscle-invasive bladder cancer after transurethral resection [5, 11].

Intravesical BCG is a type of immunotherapy that is also used to treat superficial bladder cancer [12, 13]. It is usually applied after local surgery to prevent tumor recurrence. It is given in 6-weekly intravesical instillation of bacteria, which has been proven to be superior to chemotherapy in reducing recurrence rates of the tumor [12–14]. When the BCG is instilled and processed into the bladder, it creates an inflammatory environment which in turn stimulates an immune response, resulting in attacking the cancer cells. Therefore, many researchers believed that BCG reduces tumor progression and stated that the primary role of BCG treatment is to stimulate, trigger, and activate the immune effector cells in order to attack the cancer cells. In spite of the fact that BCG instillation is regarded as the “gold standard” treatment, it has many side effects which include hematuria, pain, dysuria, and fever, to mention a few [7–14].

Immune checkpoints are negative regulators of the immune system which play important roles in maintaining self-tolerance, preventing autoimmunity, and protecting tissues from immune collateral damage. These immune checkpoints are often hijacked by tumors to restrain the ability of the immune system to mount an effective antitumor response. The tumors neutralize some immune checkpoint pathways in order to maintain immune resistance, particularly against T-cells. The T-cells are specific tumor antigens. Examples of the aforementioned checkpoints are PD-1 and CTLA4 [15–17].

Programmed cell death protein-1 (PD-1) is a protein that is encoded by the PDCD1 gene in humans. It is a cell surface receptor which belongs to the immunoglobulin superfamily and is expressed on T-cells and pro-B-cells. PD-1 binds two ligands, PD-L1 and PD-L2. The PD-1 acts as an immune checkpoint, which plays an important role in downregulating the immune system by preventing the activation of the T-cells. Hence, it decreases autoimmunity and encourages self-tolerance [18, 19]. The immune system is directly affected by the activities of PD-1 in the sense that it suppresses, blocks, and deactivates the immune cells from spreading, fighting, and attacking the cancer cells. Therefore, PD-1 aids in growth, development, and progression of the cancer. In conclusion, it disrupts and affects immunotherapy [20–24].

Transforming growth factor-beta 1 (TGF-β1) is a regulatory cytokine which suppresses immune function in cancers and in chronic viral infections. It inhibits the activation of the T-cells and subdues their proliferation. Hence, cancer cells take advantage of this immune checkpoint pathway as a way to escape and evade detection. This leads to the inhibition of antitumor immune response, resulting in cancer growth and development [25, 26].

Mathematical modeling and simulation helps in predicting treatments’ outcome, as well as describing the behavior and complex dynamics involved. Bunimovich-Mendrazitsky et al. (2007, 2008, and 2011) modeled mathematically the use of BCG in noninvasive bladder cancer, where their study identified fixed points and conditions for stability of the dynamical system [6, 8, 14]. In 2016, Bunimovich-Mendrazitsky developed a new mathematical model for combined BCG and IL-2 bladder cancer treatment which introduces the effect of TAA T-cells. Furthermore, Starkov utilized a mathematical approach for bladder cancer treatment model in the derivation of ultimate upper and lower bounds. He also presented tumor clearance conditions for BCG treatment of bladder cancer [13].

In this research, we formulate a mathematical model to study the dynamics of immune checkpoints/suppressors, immune system, and the BCG immunotherapy of bladder cancer. Moreover, we highlight the effects of immune checkpoints/suppressors on the immune system and the treatment numerically.

This paper is organized as follows. Section 1 is the introduction. Section 2 is the formulation and presentation of our model. We give the stability analysis and numerical simulations in Sections 3 and 4, respectively. In the final section, we state our conclusions and discussions.

2. Formulation of the Model

The model consists of a system of four nonlinear differential equations, which characterize the dynamics of the interaction between cancer cells , different arms of the immune system regarded as effector cells , the BCG , and all categories of immune suppressors/checkpoints as .

2.1. Dynamics of Cancer Cells

The dynamics of cancer cells is given by

Here, we assume that, in the absence of the immune system, the cancer cells grow exponentially with growth rate . The second term shows the elimination of cancer cells by the effector cells at the rate , while is the immunosuppressive factor by the immune checkpoints/suppressors, which interrupts the activities of the effector cells, with being an inhibitory parameter.

2.2. Dynamics of the Effector Cells

The dynamics of the effector cells is given by

The first term here gives the recruitment of effector cells at the rate which is directly proportional to the population of cancer cells (i.e., occurring due to the direct presence of cancer cells). shows the activation of effector cells by BCG at the rate . is the antigenicity of cancer cells which triggers an immune response in the host. It is believed that the immune checkpoints will distort both the recruitment and the activation of effector cells; hence, is the immunosuppressive response which puts a limitation on the recruitment level and interrupts the activation of effector cells, with here being an inhibitory parameter. The next term gives the elimination of effector cells by the cancer cells at the rate , and the last term describes the degradation of effector cells at the rate .

2.3. Dynamics of BCG

The dynamics of BCG is given by

The first term is the constant rate of introduction of BCG into the bladder, the second term describes the elimination of BCG by effector cells at the rate , and the third term gives the decay of BCG at the rate .

2.4. Dynamics of Immune Suppressors/Checkpoints

The dynamics of the immune checkpoints is given by

The first term gives the source of immune checkpoints at a constant rate , and the second term is the degradation of the immune checkpoints at the rate .

Finally, the interactions of the cancer cells, effector cells, BCG, and immune checkpoints together lead to the following nonlinear ordinary differential equations:

with initial conditions , , , and .

2.5. Invariance of Positive Orthant

We show that the system is positively invariant.

From the system, assume , , , and .

From , the solution is given by This implies given that Also, from , if , then Therefore, since Moreover, if , then implying since Using , if , then since .

Also, if and , then which implies given that Now consider given that .

This implies that if Hence, the positive orthant is invariant and , , , and .

3. Equilibrium and Stability Analysis

3.1. Model without Treatment

We first analyze our model in the absence of treatment :

The equilibrium points of the model are obtained by equating the equations in (6) to zero and solving simultaneously for , , , and . They are as follows:

From the invariance of the positive orthant, we concentrate only on the nonnegative equilibria assuming all initial conditions are positive.

As a result, the equilibrium point will not be considered. Moreover, exists only if the following condition is satisfied:

The Jacobian matrix obtained from (6) is given by

3.2. Stability Analysis of Equilibria of Model (6)

3.2.1. Immune Checkpoints Equilibrium:

The Jacobian matrix evaluated at yields

The eigenvalues of are

Since one of the eigenvalues is always positive, then is an unstable saddle point. Clinically, is referred to as the death equilibrium.

3.2.2. BCG-Free Equilibrium:

Assume exists; that is, ; then, substituting in yields the following eigenvalues:

Two of the eigenvalues have a real part equal to zero, which signifies neutral stability. Therefore, the equilibrium point is neutrally stable.

Conclusively, in the absence of treatment, none of the equilibrium points was found to be stable.

3.3. Model without Immune Checkpoints

Now, we analyze the model without any suppression on the immune system by the immune checkpoints. The model is given by

The equilibrium points are as follows:

The equilibrium point exists only if . This means that the cancer cells will disappear if the constant rate of introduction of BCG and activation rate of BCG are bigger than the degradation rates of both the effector cells and the BCG.

The equilibrium point also exists if

(i),(ii).

From model (13), we have the following Jacobian matrix:

3.4. Stability Analysis of Equilibria of Model (13)

3.4.1. BCG Equilibrium:

The eigenvalues of evaluated at are

The eigenvalue is always positive and the rest are negative. Therefore, the equilibrium point is an unstable saddle point.

3.4.2. Cancer-Free Equilibrium:

Assume the equilibrium point exists; then, substituting in will give the following matrix:

The eigenvalues of are

Now, if

(i) and are complex roots, then is a stable fixed point if ;(ii) and are real roots, then is a stable fixed point if and .

But since we already assume that the equilibrium point exists, then , and we can conclude that is a stable fixed point if .

This means that the effector cells activated by BCG will eradicate/destroy the cancer cells, if the constant rate of introduction of BCG, recruitment rate of effector cells, and the activation rate of effector cells by BCG are bigger than or can overcome the cancer growth rate, the rate of elimination of BCG by effector cells, and the degradation rates of effector cells and BCG altogether. Therefore, to eliminate the cancer, we increase the rate of introduction of BCG, rate of recruitment of effector cells, and activation rate of effector cells by BCG and at the same time decrease the rate of elimination of BCG by effector cells, degradation rates of both effector cells and BCG, and the cancer growth rate.

3.5. Model with Treatment and Immune Checkpoints

We now consider the dynamics of cancer cells, effector cells BCG, and immune checkpoints (see (5)).

The equilibrium points of model (5) are as follows:

Also, exists if

(i) and ;(ii) and

From model (5), we obtain the following Jacobian matrix:

3.6. Stability Analysis of Equilibria of Model (5)

3.6.1. BCG and Immune Checkpoints Equilibrium:

The eigenvalues of evaluated at are

Since one of the eigenvalues is always positive, then is an unstable saddle point.

3.6.2. Tumor-Free Equilibrium:

Assume this equilibrium point exists; then, the eigenvalues of evaluated at are as follows:

The equilibrium point is a stable fixed point if

However, condition is already true; then, is a stable fixed point if

3.6.3. Interior Equilibrium:

The eigenvalues of the Jacobian matrix are very long, complicated, and difficult to analyze. Therefore, we use numerical simulations to show the stability of the equilibrium point .

4. Numerical Illustrations

In this section, the numerical simulations of the three models will be shown. The aim here is to show the effect of immune checkpoints on the effector cells. We use MATLAB version 2016b to plot the graphs with initial populations of the compartments involved taken to be equal. Other parameters used in the numerical simulations are given in Table 1.

Table 1: List of all parameters used in numerical simulations.

Parameter
Interpretation
(units)
Estimated value
Reference

Tumor growth rate
= day⁻¹

Shochat et al., 1999

Rate of elimination of cancer cells by effector cells
cell day⁻¹

Kuznetsov et al., 1994

Inhibitory parameter

Not found

Recruitment rate of effector cells
= day⁻¹

Sud D. et al., 2006

Activation rate of effector cells by the BCG
cells⁻¹ day⁻¹

Wigginton and Kirschner, 2001

Internal production of immune checkpoints

Sandip Banerjee et al., 2015

Elimination rate of effector cells by cancer cells
cells⁻¹ day⁻¹

Kuznetsov et al., 1994

Degradation rate of effector cells
= day⁻¹

Kuznetsov et al., 1994

Rate of BCG decay
= day⁻¹

Archuleta et al., 2002

Bioeffective concentration of BCG
c.f.u./day

Cheng et al., 2004

Destruction of BCG by effector cells
cells⁻¹ day⁻¹

Wigginton and Kirschner, 2001

Degradation rate of immune checkpoints
= day⁻¹

Sandip Banerjee et al., 2015

We first plot the graph of model (6) to illustrate what happens in the absence of treatment. As expected, the cancer cells develop with the help of suppression on the effector cells by the immune checkpoints, hence dominating the effector cells and resulting in the growth and maturation of the cancer. Therefore, the numerical simulations of model (6) support this notion as shown in Figure 1.

Figure 1: Model (6) (without treatment): cancer cells grow exponentially, overcoming the effector cells , with the help of immune checkpoints .

Next, we show the behavior of model (13) (i.e., without the immune checkpoints). Here, we will see how the effector cells attack and kill the cancer cells as a result of the stimulation/activation by the BCG. Unlike in Figure 1, Figure 2 shows how the growth of the cancer cells is restricted and eventually leads to their extinction by the effector cells.

Figure 2: Model (13) (without immune suppressors): the effector cells overcome the development of cancer cells as a result of the stimulation and activation by the BCG .

The general model will now be considered. Despite stimulation and activation of the effector cells by the BCG, the immune suppressors block and deactivate their function; hence, this leads to the reduction of autoimmunity of the effector cells. Therefore, the cancer develops and grows exponentially as shown in Figure 3.

Figure 3: Model (5): despite the stimulation of effector cells by the BCG , the immune checkpoints block and deactivate the activities of the effector cells, thereby leading to the development and progression of cancer cells .

Therefore, comparing Figures 2 and 3, we will notice the effect of immune checkpoints on the effector cells. In Figure 2, the effector cells in the absence of immune suppressors fight the cancer cells, resulting in stopping their development and progression, while Figure 3 shows the progression and development of the cancer cells as a result of the presence of immune suppressors.

5. Conclusion and Discussion

In this paper, we used a system of four nonlinear ordinary differential equations to model the dynamics of cancer cells, effector cells, BCG, and immune checkpoints/suppressors in the immunotherapy of bladder cancer. We derived three possible dynamics from our model. Firstly, the model was analyzed in the absence of treatment and we studied the stability analysis of the equilibria involved. Figure 1 shows how the cancer progressed in the absence of treatment and presence of immune checkpoints/suppressors.

Secondly, we study the model without the immune checkpoints/suppressors. Conditions for stability of the equilibria involved were also given. In the absence of immune checkpoints/suppressors, the activated-effector cells have unlimited freedom to roam about and detect the cancer cells; as a result, they kill them and stop the cancer from progressing. This was shown in Figure 2.

Thirdly, we considered the dynamics of the model with treatment and the immune checkpoints/suppressors. Conditions for stability of the equilibrium points were given, and Figure 3 shows how the cancer cells grow and develop despite the application of the treatment (BCG). This is believed to be as a result of the blockage and suppression that the effector cells suffered by the immune checkpoints.

Therefore, the figures used in this paper assist in showing the effect of immune checkpoints/suppressors on the effector cells and the treatment at large. To avoid cancer progression and advancement, there is a need for action to block or limit the production of the immune checkpoints. This will take the brakes off the immune system and thereby allow it to mount a stronger and more effective attack against cancer cells.

Nivolumab is a drug recently approved by the FDA to be used alone or with other drugs to treat cancer. It is a fully human immunoglobulin (Ig) G4 monoclonal antibody directed against the negative immunoregulatory human cell surface receptor programmed cell death protein-1 (PD-1) with immune checkpoint inhibitory and antineoplastic activities. Nivolumab binds to and blocks the activation or production of immune checkpoints like PD-1. This results in the activation of T-cells and cell-mediated responses against cancer cells. So, the primary role of nivolumab is to block the immune checkpoints from suppressing the immune systems. Hence, this helps in allowing the immune cells to rise against cancer cells without any interference [16, 17].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

A. Böhle and S. Brandau, “Immune mechanisms in bacillus Calmette-Guerin immunotherapy for superficial bladder cancer,” The Journal of Urology, vol. 170, no. 3, pp. 964–969, 2003.

D. Kirschner and J. C. Panetta, “Modeling immunotherapy of the tumo—immune interaction,” Journal of Mathematical Biology, vol. 37, no. 3, pp. 235–252, 1998.

R. Kapoor, V. Vijjan, and P. Singh, “Bacillus Calmette-Guerin in the management of superficial bladder cancer,” Indian Journal of Urology, vol. 24, no. 1, pp. 72–76, 2008.

S. Bunimovich-Mendrazitsky, S. Halachmi, and N. Kronik, “Improving BCG immunotherapy for bladder cancer by adding interleukin-2 (IL-2): A mathematical model,” Mathematical Medicine and Biology. A Journal of the IMA, vol. 33, no. 2, pp. 159–188, 2016.

O. Fuge, N. Vasdev, P. Allchorne, and J. S. Green, “Immunotherapy for bladder cancer,” Research and Reports in Urology, vol. 7, pp. 65–79, 2015.

S. Bunimovich-Mendrazitsky, E. Shochat, and L. Stone, “Mathematical model of {BCG} immunotherapy in superficial bladder cancer,” Bulletin of Mathematical Biology, vol. 69, no. 6, pp. 1847–1870, 2007.

K. Kawai, J. Miyazaki, A. Joraku, H. Nishiyama, and H. Akaza, “Bacillus Calmette-Guerin (BCG) immunotherapy for bladder cancer: Current understanding and perspectives on engineered BCG vaccine,” Cancer Science, vol. 104, no. 1, pp. 22–27, 2013.

S. Bunimovich-Mendrazitsky, H. Byrne, and L. Stone, “Mathematical model of pulsed immunotherapy for superficial bladder cancer,” Bulletin of Mathematical Biology, vol. 70, no. 7, pp. 2055–2076, 2008.

E. J. Askeland, M. R. Newton, M. A. O’Donnell, and Y. Luo, “Bladder Cancer Immunotherapy: BCG and Beyond,” Advances in Urology, vol. 2012, pp. 1–13, 2012.

G. Redelman-Sidi, M. S. Glickman, and B. H. Bochner, “The mechanism of action of BCG therapy for bladder cancer-a current perspective,” Nature Reviews Urology, vol. 11, no. 3, pp. 153–162, 2014.

P. Andius and S. Holmäng, “Bacillus Calmette-Guérin therapy in stage Ta/T1 bladder cancer: Prognostic factors for time to recurrence and progression,” BJU International, vol. 93, no. 7, pp. 980–984, 2004.

D. L. Lamm, D. E. Thor, S. C. Harris, J. A. Reyna, V. D. Stogdill, and H. M. Radwin, “Bacillus Calmette-guerin Immunotherapy of Superficial Bladder Cancer,” The Journal of Urology, vol. 124, no. 1, pp. 38–42, 1980.

K. E. Starkov and S. Bunimovich-Mendrazitsky, “Dynamical properties and tumor clearance conditions for a nine-dimensional model of bladder cancer immunotherapy,” Mathematical Biosciences and Engineering, vol. 13, no. 5, pp. 1059–1075, 2016.

S. Bunimovich-Mendrazitsky, J. Claude Gluckman, and J. Chaskalovic, “A mathematical model of combined Bacillus Calmette-Guerin (BCG) and interleukin-2 (IL-2) immunotherapy of superficial bladder cancer,” Journal of Theoretical Biology, vol. 277, pp. 27–40, 2011.

M. Thibult, E. Mamessier, J. Gertner-Dardenne et al., “PD-1 is a novel regulator of human B-cell activation,” International Immunology, vol. 25, no. 2, pp. 129–137, 2013.

D. M. Pardoll, “The blockade of immune checkpoints in cancer immunotherapy,” Nature Reviews Cancer, vol. 12, no. 4, pp. 252–264, 2012.

M. A. Postow, M. K. Callahan, and J. D. Wolchok, “Immune checkpoint blockade in cancer therapy,” Journal of Clinical Oncology, vol. 33, no. 17, pp. 1974–1982, 2015.

F. Paillard, “Commentary: Immunosuppression Mediated by Tumor Cells: A Challenge for Immunotherapeutic Approaches,” Human Gene Therapy, vol. 11, no. 5, pp. 657-658, 2000.

K. C. Ohaegbulam, A. Assal, E. Lazar-Molnar, Y. Yao, and X. Zang, “Human cancer immunotherapy with antibodies to the PD-1 and PD-L1 pathway,” Trends in Molecular Medicine, vol. 21, no. 1, pp. 24–33, 2015.

M. E. Keir, M. J. Butte, G. J. Freeman, and A. H. Sharpe, “PD-1 and its ligands in tolerance and immunity,” Annual Review of Immunology, vol. 26, pp. 677–704, 2008.

K. A. Hofmeyer, H. Jeon, and X. Zang, “The PD-1/PD-L1 (B7-H1) pathway in chronic infection-induced cytotoxic T lymphocyte exhaustion,” Journal of Biomedicine and Biotechnology, vol. 2011, Article ID 451694, 2011.

Y. Latchman, C. R. Wood, T. Chernova et al., “PD-L2 is a second ligand for PD-1 and inhibits T cell activation,” Nature Immunology, vol. 2, no. 3, pp. 261–268, 2001.

Y. Wang, Y. Cao, X. Yang et al., “Effect of TLR4 and B7-H1 on Immune Escape of Urothelial Bladder Cancer and its Clinical Significance,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 3, pp. 1321–1326, 2014.

Y. Iwai, M. Ishida, Y. Tanaka, T. Okazaki, T. Honjo, and N. Minato, “Involvement of PD-L1 on tumor cells in the escape from host immune system and tumor immunotherapy by PD-L1 blockade,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 99, no. 19, pp. 12293–12297, 2002.

A. Yoshimura and G. Muto, “TGF-beta function in immune suppression,” Current Topics in Microbiology and Immunology, vol. 350, pp. 127–147, 2011.

S. Banerjee, S. Khajanchi, S. Chaudhuri, and M. M. Alonso, “A Mathematical Model to Elucidate Brain Tumor Abrogation by Immunotherapy with T11 Target Structure,” PLoS ONE, 2015.