Computational and Mathematical Methods in Medicine

Volume 2017 (2017), Article ID 3573082, 9 pages

https://doi.org/10.1155/2017/3573082

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

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

Correspondence should be addressed to Farouk Tijjani Saad

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 byHere, 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 byThe 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 byThe 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 byThe 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 yieldsThe eigenvalues of areSince 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 byThe 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 areThe 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 areSince 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 ifHowever, 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.