The role of interleukin-10 in immune response to hepatitis B virus and liver cancer co-existence dynamics ========================================================================================================= * Paul Chataa * Samuel M. Naandam * Francis T. Djankpa ## Abstract Experimental evidence confirms that interleukin-10 plays a critical role in clearing acute hepatitis B virus infection. This paper aims to develops a mathematical model to explore the dynamics of how the immune system responds to hepatitis B virus (HBV) and coexisting liver cancer within the liver cell population. Unlike previous models; we categorize liver cells into various stages of infection. We determine the invasion probability for transmission dynamics, specifically the basic reproduction number, ℝ, for populations of uninfected macrophages with and without cancer cells. Stability analyses of virus-free and virus equilibrium states are provided, along with numerical simulations to validate analytical findings. The impact of different branches of the immune response on model dynamics is assessed. Simulations predict the time at which T helper-1 cells surpass cytotoxic T cells (switching time), correlating positively with the proliferation rate of interleukin-10 (*ρ*3). Further numerical simulations demonstrate that interleukin-10 contributes to HBV persistence by inhibiting the immune response, thereby allowing the virus to evade immune surveillance and establish chronic infection through the suppression of cytotoxic T lymphocytes (CTLs), which are essential for clearing infected cells. Keywords * Macrophages * Invariant Region * Routh-Hurwitz criterion ## 1 Introduction A viral infection that targets the liver cells called hepatitis B can result in both acute and chronic illnesses. According to [1], this condition kills 750, 000 people per year, roughly 300, 000 of which result in liver cirrhosis and hepatocellular carcinoma [2]. A third of the world’s population is thought to be infected with the hepatitis B virus, making it not only one of the oldest but also one of the most dangerous viral dangers to human health [3]. According to the World Health Organisation (WHO) estimates, 1.5 million new cases of hepatitis B are reported annually, leaving 296 million people living with the virus. WHO estimates again indicated that hepatitis B caused 820000 in fatalities in 2019, with the majority being caused by cirrhosis and hepatocellular carcinoma, the predominant form of liver cancer. Every year, millions of individuals worldwide perish from liver cancer, and current trends suggest that millions more people will do so in the future. Reports from Europe and America shows that the hepatitis B infection is quite low (i.e., less than 1%). With a chronic infection incidence of 5 to 10 per 100 adults in Asia and Africa, the disease is still a significant burden in these regions. The hepatitis B virus, a member of the hepadnaviridae family of viruses, is what causes hepatitis B disease. Once a person is infected, it is particularly challenging to get rid of the hepatitis B virus due to its largely double-stranded structure [2]. The incubation periods for the hepatitis B virus range from 30 to 180 days. The virus may be detected within 30–60 days of infection, and if it is transferred during infancy or childhood, it may persist and result in chronic hepatitis B. The virus is frequently vertically transferred from mother to child, causing recurrent infection and, in the majority of instances (approximately 90%), chronic infection [4, 5]. The second potential route for virus spread is through horizontal adult-to-adult sexual contact, intravenous drug use, unhygienic behavior, and blood transfusion. Only 5-10% of adults have persistent infections as a result of this sort of transmission, according to [4] and [5]. When the body generates sufficient immune responses against the infection, these adults recover. Activation of robust and diversified CD4 (T helper 1) and CD8 (cytotoxic T lymphocytes) T-cells, production of protective, neutralizing antibodies against HBV surface antigen (HBsAg), and expression of antiviral cytokines in the liver, such as gamma interferon (type-2 interferons) and tumor necrosis factor alpha [6, 5], and the creation of refractory cells, which are immune to reinfection [7, 8] are examples of such responses. However, newborns who are not immunized and immunosuppressed adults tend to move to chronic HBV infection stage [9]. These people have humoral and cellular immune systems that are weak and ineffective, which leads to ongoing viral replication and HBV surface antigenemia [2, 10]. The relative contributions of the immune system’s many components are poorly understood, particularly the functions of anti-inflammatory cytokines like interleukin-10 in the development and progression of infection. The World Health Organization has set a goal of eliminating liver cancer and the hepatitis B virus (HBV) by 2030 [11, 12]. Treatment of sick people as well as immunization of newborns and susceptible adults with the hepatitis B vaccine are now used as control techniques to lessen the spread of the hepatitis B virus. Campaigns aimed at educating people about how alcohol and smoking can prevent liver cancer and help slow the spread of the illness is another way to control the spread of these two diseases. These preventative measures are meant to stop acute hepatitis B virus infection and prevent people with liver problems from developing liver cirrhosis again. However, the incidence of hepatitis B and liver cancer is alarmingly high, with millions of cases and deaths each year. Although there has been progress toward eliminating HBV, the disease still poses a serious threat to public health. 316 million persons worldwide were living with chronic HBV infection in 2019 [13]. This represents a prevalence of chronic HBV infection of 4.1%. The likelihood of clearing liver cancer is decreased when HBV infection coexist with liver cancer. Therefore, it is essential for the research community to conduct scientific study on hepatitis B and liver cancer in both clinical and theoretical forms. The development of effective hepatitis B and liver cancer therapies is a vast area of medical study that will influence our comprehension of tumor-immune dynamics. Although the majority of people (about 90% in adults) assemble a successful and defensive cell-driven mechanism that prevents them from developing chronic hepatitis B disease, which will eventually result in liver cancer, hepatitis B is still a leading cause of death globally (about 3 million deaths per year), according to estimates [2]. Numerous researchers have examined different facets of HBV and liver cancer dynamics and the immune response during infection using mathematical models. In order to investigate acute HBV infection and the significance of time lag in effector cell activation and expansion, [3, 14] expanded a standard model of immune response to incorporate the time delay in recruit naive T cells. Subsequently, they also investigated the function of pre-existing or vaccine-induced antibodies in containing the HBV infection [1]. Instead of using a mass action to account for a finite liver size and susceptibility to HBV infection, [15] employed a typical incidence function in their study of HBV transmission dynamics. A time-delayed version of the model put forth in [15] has been created by [16]. [17] model was able to more accurately describe the existing data and produce more realistic results for the basic reproduction number by using a standard incidence and a logistic growth for the hepatocyte population. [18] have examined potential inadequacies in the synchronization of distinct branches of the adaptive immune response, particularly the CTLs and antibodies, in the context of HBV infection. Yet, current modeling studies have not been able to aid in the development of successful HBV total eradication and drug treatments. A reason for this may lie in the fact that almost all models of HBV and liver cancer ignore the interleukin-10 and the differentiation of naive T cells aspects of infection which we have determined to be quite important for designing HBV and liver cancer control strategies. Moreover, it is impossible to ignore the critical part the immune system plays in the dynamics of the hepatitis B virus and liver cancer infection. Interleukin-10 (*I*10) is a cytokine that modulates both innate and adaptive immunity, primarily by exerting anti-inflammatory effects. So, formulating a new model to examine how interleukin-10 and other components of the immune system regulates both cell-mediated (most especially the production of T helper-1 and cytotoxic T cells) and innate immunity proved beneficial. Thus, we concentrated on limiting the concentration of the virus and the damage to the liver. We achieved this goal using different strategies which include interferon immunity (i.e type-1 and type-2) by removing the substrate that the virus needs for reproduction (i.e., the healthy cells), cellular immunity (T cells, NK cells, effector B cells, interleukin-10 cells) by removing the source of new viruses (i.e., the infected cells), and adaptive immunity (HBV-antibodies) by lowering the effectivity concentration of the virus. Another innovative part of our model is how it predict when T helper-1 (*T*1) cells will outnumber cytotoxic T cells (*T*2), a period we call the “switching time”. This period predict when acute hepatitis B infection transitions to chronic infection and eventually to liver cancer and the conditions that led to these transitions. The paper is organized as follows; Section 2 is devoted to the mathematical formulation of the model. The model basic properties was presented in Section 3. In Section 4 we present the model analysis. Local asymptotic stability of the virus-free equilibrium and global asymptotic stability of virus-free analysis was presented in Section 5. In Section 6, we study the numerical results and sensitivity analysis of the proposed models and present the results in the form of plots, and a discussion of the results is presented in section 7. ## 2 Materials and methods ### 2.1 Basic Model Formulation A compartmental mathematical model for the immune response dynamics to hepatitis B virus and liver cancer infection was constructed in order to comprehend the numerous changes that occur in the immune system’s response to the coexistence dynamics of these two diseases. The model is based on information obtained from new scientific explanations of the basic features of infections with the hepatitis B virus and infections that cause liver cancer, as well as on theoretical and experimental studies conducted by other researchers [5, 19, 14, 12]. The outcome of a hepatitis B virus infection and the rate at which acute infections lead to liver cancer are significantly influenced by the roles of antibodies, cell-mediated immune responses, and innate immune system responses as well as the cytokines. The total number of human liver cells, or hepatocytes, was categorized into two populations: uninfected macrophages without cancer *M*, and those with cancer *M*1. Both uninfected macrophages with and without cancer (i.e *M* and *M*1) are triggered from the blood stream to the site of infection because the hepatitis B virus is present in the body cells. Both of the uninfected macrophages then phagocytose the viruses, which causes them to become infected at time *t*. These infected macrophages are identified by *I**M*. It is assumed that both uninfected macrophages are created at constant rates ![Graphic][1] and ![Graphic][2] and die at constant rates ![Graphic][3] and ![Graphic][4] per cell. Hepatitis B virus particles infect them at a rate proportional to the product of *M*, *M*1, and *V* (i.e., *β*1*V M* and *β*2*V M*1) with constant proportionality rates of *β*1 and *β*2, respectively. Infected macrophages pass away at a steady rate of ![Graphic][5] per cell. In the class of infected macrophages, some of the individual macrophages develop cirrhosis, pick up liver cancer, and proceed at a rate of *ζ* to the liver cancer macrophage population. Through contaminated blood transfusions and having sex with an infected individual, liver cancer macrophages can get hepatitis. Hepatitis B virions *V* are recruited at a rate of Λ*V* whereas infected macrophages create additional free virions at a rate of *ω*. The virus degenerates at a constant rate of *µ**V* in each cell. Figure 1 shows interactions between all cell types in the human liver. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F1) Figure 1: Show the schematic diagram of immune response to HBV and liver cancer co-existence. Yellow rect-angles indicate host liver cells (uninfected macrophages with and without cancer, and infected macrophages) population, red rectangles denote T cell’s population, blue rectangle represent antibodies population, pink rectangle represent effector B cells population, green rectangles show cytokines (type-1 and type-2 interferon) population, the violet rectangle is the innate (*N**K* cells) population, and purple rectangle indicates virus particles (virions) population. Also, thick lines with arrowheads represent contacts, thick lines with bar heads represent clearance (destroy) while dash lines represent interactions. Invading infections are specially recognized by antibodies, which bind to them and render them ineffective. Sub-viral particles (SVPs), which can occur in quantities up to a thousand to one million times higher than the infectious virions, are produced in excess by HBV-infected cells during infections. These SVPs may have an impact on how the immune system of the host responds to the HBV infection. According to [1, 20], the sub-viral particles serve as immune system spies that tempt antiviral antibodies away from binding to the hepatitis B virus. They may also promote tolerance during neonatal infection, delaying the development of neutralizing antibodies and allowing the hepatitis B virus to avoid antibody detection. The creation of macrophages, T cells, and cytokines triggers a unique immune response that frequently mediates this issue. We assume that in the absence of infection, HBsAg-specific antibodies *A* are created over time rate Λ*A* and degrade at a per capita rate *µ**A*. Nevertheless, throughout an infection, antibodies are created at a rate *ψ* that is proportionate to the viral load and kill virions at a rate of *δ*4 [21]. Macrophages T cells are created at a rate of ![Graphic][6] and perish at a constant rate of ![Graphic][7] in each cell. Following internalization by the infected macrophages, the virus continues to spread by feigning an acute stage (i.e., the stage at which virus reproduction is extremely small) for a predetermined amount of time in order to evade intracellular killing systems. The virus population in the host liver cells grows as a result of the virus’s ongoing intracellular replication in an infected macrophages cell. Naive T cells from the bone marrow’s thymus are induced, triggered, and activated as a result of virus replication inside infected macrophages up to a point. The pro-inflammatory cytokines actions of the virus cause the naive T cells to differentiate into either T helper 1; *T*1 or T helper 2; *T*2 cells. The number of infected macrophages per unit area and the amount of virus present, respectively, determine the cellular *T*1 and *T*2 immunological responses to HBV infection. While T helper 1 cells *T*1 differentiate at a rate *ρ*1 and die at a per capita constant rate ![Graphic][8] with each cell specialized in carrying out a specific task, T helper 2 cells become cytotoxic T lymphocytes *T*2 at a differentiation rate *ρ*2 and die at a per capita constant rate ![Graphic][9]. While T helper 1 cells; *T*1 activate the infected macrophages cells to generate additional cytokines at a rate of *π*, cytotoxic T cells destroy their targeted virus cells primarily by releasing cytotoxic granules to kill the virus cell at a rate of *δ*3. The naive T cells undergo additional differentiation to form regulatory T cells, which control the generation of *T*1 and *T*2 cells. It has been discovered that the production of regulatory T cells, which have been found to perform a lytic activity for the production of interleukin-10; *I*10, is delayed as a result of the virus’s ability to appear temporarily inactive during the acute stage in the infected macrophages [22]. Interleukin-10 is a cytokine that predominantly has anti-inflammatory effects and affects both *T*1 and *T*2 adaptive immunity. During HBV infections, it also coordinates the innate immune system and antibodies. Naive T cells create *I*10 at a rate of *ρ*3 and it dies at a constant rate of ![Graphic][10] in each cell. It alternately suppresses and supports the development of *T*1 and *T*2 cells. Antigen-presenting cells (APCs) cannot produce pro-inflammatory cytokines like *I*12, and they are also prevented from up-regulating molecules involved in antigen presentation and lymphocyte activation (i.e. the minus sign stands for down regulation of cytokines), which is how *I*10 suppresses T helper 1 cells. Although these cytokines primarily have suppressive effects on the immune system, they also have certain immune-stimulating properties, such as increasing the generation of cytotoxic T lymphocytes also known as T helper 2 cells. The plus sign denotes cytokine up-regulation. Another method by which the immune system can combat and stop the progression of HBV infection into liver cancer is through effector B cells, also referred to as plasma cells. These effector cells releases antibodies, heals infected cells, and primes T cells (including cytotoxic T cells and T helper 1 cells), which trigger cell-mediated reactions. Effector cell levels are maintained at a certain homeostatic level after viral clearance as a result of long-lived plasma and memory B cells [6]. When viral peptide MHC class I molecules, which include viral proteins, are displayed on the cell surface, effector B lymphocytes begin to target those cells as infected cells. In the absence of infection, we assumed that effector B cells proliferate at a rate of Λ*E* and perish at a constant rate of *µ**E* per cell. We model these two impacts; killing and curing of infected cells in one reaction, where effector B cells kill and cure infected macrophages at a rate *δ*2. This is because it has been demonstrated that cured cells lose their resistance to productive infection at a slow rate up to the order of 105 per day, and effector cells also kill infected macrophage cells. When there is an infection, the population of effector cells increases by a factor of *σI**M* *E*, where *σ* is the maximum proliferation rate. This happens in a contaminated cell density-dependent manner which depends on the density of an antigen. The hepatitis B virus is highly adaptable and has developed ways of stopping MHC molecules from getting to the cell surface to display viral peptides to avoid being detected by T cells. If this happens, the T cell may not know that there is a virus inside the infected cell. However, another immune cell specializing in killing cells with a reduced number of MHC class I molecules on their surface is the natural killer cell or *N**K* cell for short. When the *N**K* cells finds a cell displaying fewer than normal MHC molecules it releases toxic substances, in a similar way to cytotoxic T cells to kill the viral-infected cell. Similar to effector B cells, in the absence of infection, natural killer *N**K* cells, which are the first line of defense against non-self pathogens are considered to grow at a rate ![Graphic][11], and die at a per capita constant rate ![Graphic][12]. When examining the immune response to HBV and liver cancer infection co-dynamics, the importance of cytokines in immune system dynamics cannot be understated. Interferon, a class of tiny proteins produced and released by virally infected macrophage cells, is crucial for immunological defense against hepatitis B virus infection. By directly impeding the ability of the virus to replicate within an infected macrophages cell, interferon stops viral replication. We refer to the Type-1 interferon (*IFN**α/β*) as *F*1, which is created by infected macrophages cells at a rate of *ε* and destroyed at a per-cell rate of ![Graphic][13] [2, 8]. Type-2 interferon *IFN**γ*, on the other hand, is denoted by *F*2 which is also produced by natural killer cells *N**K* [6, 9, 10] at rate *η* and they are lost at a per capita rate ![Graphic][14]. *IFN**γ* generates protein-10 that can both energize and enroll *N**K* cells [10], whereas *IFN**α/β* are able to activate the natural killer *N**K* cells during infection [23]. Therefore, it is assumed that the cumulative action of interferon on initiating *N**K* cells happens at *θ*1*N**K**F*1 + *θ*2*N**K**F*2; where *θ*1 represents the rate at which Type-1 interferon activates *N**K* cells and *θ*2 represents the rate at which Type-2 interferon activates *N**K* cells respectively. In addition to increasing the number of new *N**K* cells produced, *IFN**α/β* and *IFN**γ* boost the the *N**K* cells’ pathogenicity and effector B cells, respectively [24]. Hence, we assumed that natural killer cells and effector B cells destroy infected macrophages cells at *δ*1 (1 + χ1*F*1) *I**M* *N**K* and *δ*2 (1 + χ 2*F*2) *I**M* *E*; where χ1 represent the rate at which Type-1 interferon increase the pathogenicity of *N**K* cells and χ2 represent the rate at which Type-2 interferon increase the pathogenicity of effector B cells respectively. We additionally assume the following in addition to the presumptions listed above: * It is assumed that a portion of the liver cells have acquired malignancy from external sources such as excessive alcohol intake, smoking etc. We now introduce HBV infection and studied the new co-existence dynamics. * Time delays in the replication of individual cell components are not taken into account. * The populations of macrophages cells and virus are assumed to be uniformly distributed over the system at all times. The general process of the model is described in Figure 1 ### 2.2 Model Equations The following set of differential equations represents the immune response to hepatitis B virus and liver cancer co-existence model, based on the assumptions stated above. Clonal selection theory, mass-action kinetics, and the balances between cell and molecular population proliferation and mortality underpin interactions in the model system. ![Formula][15] The system described in reference 1 is not dimensionless. The term “non-dimensionless” implies that, in the absence of knowledge regarding the actual values of the variable, we assume these values to be proportional. ## 3 Basic Model Properties In this section, we present the positivity and boundedness of solutions of the model system. ### 3.1 Positivity of Solutions To ensure the mathematical and epidemiological relevance of system 1, it is essential to ensure that solutions originating from positive initial conditions remain positive for all *t >* 0. Therefore, it is necessary to establish that each state variable is non-negative, as negative cell population densities are illogical. This outcome is achieved through the lemma provided below. Lemma 1. If *T*(0) *>* 0, *T*1(0) *>* 0, *T*2(0) *>* 0, *M*(0) *>* 0, *M*1(0) *>* 0, *I**M* (0) *>* 0, *V* (0) *>* 0, *A*(0) *>* 0, *E*(0) *>* 0, *N**K*(0) *>* 0, *I*10(0) *>* 0, *F*1(0) *>* 0 and *F*2(0) *>* 0, then the solutions *T*(*t*) *>* 0, *T*1(*t*) *>* 0, *T*2(*t*) *>* 0, *M*(*t*) *>* 0, *M*1(*t*) *>* 0, *I**M* (*t*) *>* 0, *V* (*t*) *>* 0, *A*(*t*) *>* 0, *E*(*t*) *>* 0, *N**K*(*t*) *>* 0, *I*10(*t*) *>* 0, *F*1(*t*) *>* 0 and *F*2(*t*) *>* 0 of model 1 are positive for all *t* ≥ 0. *Proof*. We define ![Formula][16] this implies that ![Formula][17] Considering the first equation in system 1, we have ![Formula][18] It follows from (2) that ![Formula][19] Since ![Graphic][20]. Hence, ![Formula][21] Seperating variables and integrating both sides of (3) from 0 to *τ* gives ![Formula][22] Similarly the remaining equations gives ![Formula][23] Clearly, we have shown that *T*(*t*) *>* 0, *M*(*t*) *>* 0, *M*1(*t*) *>* 0, *V* (*t*) *>* 0, *A*(*t*) *>* 0, *E*(*t*) *>* 0, *N**K*(*t*) *>* 0, and *T*1(*t*) ≥ 0, *T*2(*t*) ≥ 0, *I**M* (*t*) ≥ 0, *I*10(*t*) ≥ 0, *F*1(*t*) ≥ 0, *F*2(*t*) ≥ 0. Therefore, the solutions *T*(*t*), *T*1(*t*), *T*2(*t*), *M*(*t*), *M*1(*t*) *>* 0, *I**M* (*t*), *V* (*t*), *A*(*t*), *E*(*t*), *N**K*(*t*), *I*10(*t*), *F*1(*t*) and *F*2(*t*) of the system 1 remains positive for all *t* ≥ 0. □ ### 3.2 Boundedness of Solutions The model population variables and parameters changes in system 1 were studied keenly under certain conditions. Considering the variables and parameters in the model system to be positive entities for all time values *t* ≥ 0, we show that all possible solutions are uniformly bounded. Lemma 2. The solutions of the system with any non-negative initial conditions are bounded for all *t* ≥ 0 in the biologically feasible region defines by the set ![Formula][24] *Proof*. Considering the various T cells (i.e *T*, *T*1 and *T*2) and summing them to obtain the total T cells population, we have ![Formula][25] It follows from (4) that ![Formula][26] Integrating (5) from 0 to *t*, we obtain ![Formula][27] Taking the limit of *T* at *t* → ∞ such that ![Formula][28] Hence, the T cells population is bounded. That is ![Formula][29] Similarly, the remaining populations give the following ![Formula][30] We can therefore conclude that for *t >* 0, any solution of the system is bounded in the region Ω. Thus, it is feasible to consider the dynamics and flow of the T cells, macrophages, virus, antibodies, effector B cells, natural killers and cytokines populations as described by the model system within the invariant region Ω. ## 4 Model Analysis In this section, we determine the virus-free and virus-persistence equilibria of the model system. ### 4.1 Virus-free Equilibrium and Reproduction Number It is believed that some liver cells have developed cancer due to external factors like excessive alcohol consumption and smoking. We have now introduced an HBV infection and examined the resulting co-existence dynamics. Hence, at the virus-free equilibrium state, no hepatitis B virus is present which implies that there are no infected macrophages, since both uninfected macrophages (with and without cancer) would have no virus to swallow-up to cause them to be infected. That is *V* ⋆ = 0 which immediately implies ![Graphic][31] and ![Graphic][32]. The lack of infected macrophages indicates that the hepatitis B disease is not being transmitted, which in turn means that naive T cells will not perform any lytic activity. Therefore, at the virus-free equilibrium, the populations of all cells involved in the immune interactions are given by ![Formula][33] Next, we use the Next Generation Matrix technique described by [25] and use in [26] to calculate the reproduction number of the model system. The number of further HBV and liver cancer infections brought on by a single infected macrophage placed into a population of liver cells that are completely susceptible to infection is known as the reproduction number (ℝ), for an infection. To put it another way, when all liver cells are free of infection, ℝ represents the amount of secondary infections that one infected macrophage cell produces. The model system’s equations describing the generation of fresh infected macrophages and modifications in their states are the first things we look at. The source of these equations is given by ![Formula][34] The set of these equations in system 6 is called infected subsystem. Linearizing the infected subsystem about the virus-free steady state is the initial stage, according to the existing rule. When we set ![Formula][35] where *T* denote the transpose, then the infected subsection can be written in the form: ![Formula][36] where ![Formula][37] Taking the Jacobian of 𝔽 (𝕏) at the virus-free equilibrium state, we have ![Formula][38] Also taking the Jacobian of 𝕍(𝕏) at the virus-free equilibrium, we obtain ![Formula][39] It follows from (7) that ![Formula][40] where ![Formula][41] Thus, ![Formula][42] Thus, the next generation matrix which represents the number of secondary infections produced by one infected macrophages cell when all liver cells are uninfected is given as ![Formula][43] Therefore, ![Formula][44] where ![Formula][45] ## 5 Model Stability Analysis The local and global stability analysis of the virus-free equilibrium point, *ε*, as well as the virus-persistence equilibria were investigated in this section. ### 5.1 Local Stability of Virus-free Equilibrium By calculating the eigenvalues of the linearized Jacobian matrix at the virus-free equilibrium, we are able to ascertain the local stability of the equilibrium when there is a slight disruption. That is when a small number of infected cells are injected into the population of susceptible liver cells, it is possible to completely eliminate the virus infection from the liver cells due to the stability of virus-free equilibrium when ℝ *<* 1. Theorem 1. The virus-free equilibrium of system 1 is locally asymptotically stable if ℝ *<* 1 and unstable if ℝ *>* 1. *Proof*. We investigate the proof of Theorem 1 by linearization approach. We make used of the Jacobian matrix associated with system 1 at the virus-free equilibrium and obtain ![Formula][46] where ![Graphic][47], *κ*5 = (*ψ* + *µ**V*), *κ*6 = *ψ, κ*7 = *µ**A*, ![Graphic][48]![Graphic][49] The eigenvaules from the jacobian matrix ![Graphic][50] are obtained to be ![Graphic][51], *λ*6 = −*µ**A*, *λ*7 = −*µ**E*, ![Graphic][52] with a polynomial equation ![Formula][53] Since all the coefficients of the characteristics polynomial in (8) are positive when ℝ *<* 1, by the Routh-Hurwitz criterion the solutions to the characteristic polynomial have negative real parts. Therefore all the eigenvalues of the Jacobian matrix ![Graphic][54] have negative real part when ℝ *<* 1. Hence, we can conclude based on Routh-Hurwitz criterion that the virus-free equilibrium ![Graphic][55] of system 1 is locally asymptotically stable. ### 5.2 Global Stability of Virus-free Equilibrium Theorem 2. From system 1, the virus-free equilibrium ![Formula][56] is globally asymptotically stable if ℝ *<* 1 and conditions (*G*1) and (*G*2) are satisfied. We apply the approach of [27], to prove the global stability of the virus-free equilibrium. Theorem 3. If a model system can be written in the form: ![Formula][57] where *X* ∈ ℝ*m* denotes the number of uninfected macrophages cells and *I* ∈ ℝ*n* denotes the number of infected macrophages cells including latent, acute and exposed cells. *U* (*X*⋆, 0) denotes the virus-free equilibrium of the system. Then the conditions (*G*1) and (*G*2) must be satisfied to guarantee local asymptotic stability. *G*1 : ![Graphic][58], *X*⋆ is globally asymptotically stable. *G*2 : *G*(*X, I*) = *AI* − *Ĝ* (*X*, 0) ≥ 0 for (*X, I*) ∈ Ω, where *A* = *D**i**G*(*X*⋆, 0) is a Metzler matrix (the off diagonal elements of *A* are non-negative) and Ω is the region where the model makes biological sense and mathematically well posed. Then the fixed point *U* = (*X*⋆, 0) is globally asymptotically stable equilibrium of the immune response to hepatitis B virus and liver cancer infection model 1 provided ℝ *<* 1. *Proof*. From the model system in 1, we have ![Formula][59] and *I* ∈ ℝ2 = (![Graphic][60], *V* ⋆). Hence, for condition (*G*1), we have ![Formula][61] and ![Formula][62] It follows that ![Formula][63] where ![Graphic][64], *ν*4 = *µ**A*, ![Graphic][65]. The eigenvaules from the matrix *F* (*X*, 0) are obtained to be ![Graphic][66], *λ*6 = −*µ**A*, *λ*7 = −*µ**E*, ![Graphic][67] Since all the eigenvalues are real and negative, it follows that *X*⋆ is always globally asymptotically stable. Also, applying Theorem 3 to the immune response model system 1 gives ![Formula][68] Hence, ![Formula][69] Therefore, ![Formula][70] So, *A* is a Metzler matrix with non-negative off-diagonal elements. We observed that ![Formula][71] because ![Graphic][72]. That is in the absence of infection, the proportional product of the macrophages with and without cancer at the virus-free equilibrium are obtained to be ![Graphic][73] and ![Graphic][74] respectively while *V M*1 and *V M* represent the proportion to the product of macrophages with and without cancer in the population. Therefore, the virus-free equilibrium *ε* is globally asymptotically stable. ### 5.3 Virus-persistence Equilibrium (VPE) and Stability Analysis In this section, we determine the virus-persistence equilibrium points by solving system 1 simultaneously for the state variables ![Graphic][75], *V* ⋆, *A*⋆, *E*⋆, ![Graphic][76] and ![Graphic][77] respectively. The virus-persistence equilibrium points are the steady state solutions where the hepatitis B virus and liver cancer infection cannot be totally eradicated but remains to invade the total cells populations. It is assumed that a portion of the liver cells have acquired liver cancer from external sources such as excessive alcohol intake, smoking etc. We now introduce HBV infection and studied the new co-existence dynamics. For the hepatitis B virus infection to remain in the host liver cells, the virus compartment does not turn to zero. As a result of this, the system depends on all the state variables. At the virus present state, uninfected macrophages become infected and this triggers activation of infected macrophages by cytokines. So, at the virus-persistence equilibrium the following equations are satisfied: ![Formula][78] Due to the high dimensionality of system 1, it is prudent to look at a significant number of possible boundary steady states. Hence, in order to systematically find and analyze these steady state, we begin with boundary steady states characterized by the presence hepatitis B virus particles (i.e *V ≠* 0) with different combinations of ![Graphic][79] or ![Graphic][80], *E*⋆ = 0 or *E*⋆ *≠* 0 and ![Graphic][81]or ![Graphic][82]. Case I: For the first scenario where *V* ⋆ *≠* 0, ![Graphic][83], *E*⋆ *≠* 0, ![Graphic][84], we obtain a boundary steady state of the form: ![Formula][85] where ![Graphic][86] and ![Graphic][87]. This equilibrium state ![Graphic][88] exists when ![Graphic][89], and ![Graphic][90]![Graphic][91]. Using the the parameter values in Table 1, we compute the components of ![Graphic][92] and obtain ![Formula][93] View this table: [Table 1:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/T1) Table 1: Parameters of model and their description Linearized system 1 about the equilibrium state ![Graphic][94], we obtain a Jacobian matrix of the form: ![Formula][95] where ![Formula][96] The coexisting equilibrium point ![Graphic][97] exhibits eigenvalues ![Formula][98] corresponding to an unstable inward spiral. The production of antibodies, effector B cells, naive T cells, and natural killers renders the virus persistence equilibrium state irrelevant. This stabilization promotes stability in the virus-free equilibrium, facilitating viral clearance. Case II: For the second scenario where *V ≠* 0, ![Graphic][99], *E*⋆ = 0,![Graphic][100], we obtain a boundary steady state of the form: ![Formula][101] Considering the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][102] The eigenvalues of ![Graphic][103] are: ![Formula][104] This indicates that ![Graphic][105] is a saddle-type critical point. Therefore, around this point, the system exhibits stable behavior. This equilibrium is considered a chronic equilibrium because there are no T cells, natural killer cells, effector B cells, or cytokines that can clear the virus. Thus, the presence of antibodies alone is not sufficient to clear the viral infection. Case III: For the third scenario where *V* ⋆ *≠* 0, ![Graphic][106], *E*⋆ = 0, ![Graphic][107], we obtain a boundary steady state of the form: ![Formula][108] Considering the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][109] The eigenvalues of ![Graphic][110] are: ![Formula][111] Regarding the equilibrium ![Graphic][112], the eigenvalues indicate that it is a saddle point, and the system is stable around ![Graphic][113]. As T cells, effector B cells, and other cytokines such as interferon alpha, beta, and gamma are absent at this equilibrium, there is no chance of viral clearance. This implies that the presence of antibodies and natural killer cells alone is not sufficient to destabilize the virus-persistence equilibrium required for viral clearance. Case IV: For the fourth scenario where *V* ⋆ *≠* 0, ![Graphic][114], *E*⋆ *≠* 0, ![Graphic][115], we obtain a boundary steady state of the form: ![Formula][116] Considering the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][117] The eigenvalues of ![Graphic][118] are: ![Formula][119] Corresponding to the equilibrium point ![Graphic][120], the eigenvalues indicate it is a stable inward spiral, and the system remains stable around this point. In this equilibrium state, the absence of T cells, natural killer cells, and other cytokines such as interferon alpha, beta, and gamma means there is no possibility of viral clearance. This indicates that the presence of antibodies and effector B cells alone cannot disrupt the virus-persistence equilibrium needed for viral clearance. Case V: For the fifth scenario where *V* ⋆ *≠* 0, ![Graphic][121], *E*⋆ = 0, ![Graphic][122], we obtain a boundary steady state of the form: ![Formula][123] Given the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][124] As corresponding to the equilibrium ![Graphic][125], the eigenvalues are ![Formula][126] The equilibrium ![Graphic][127] is an unstable saddle point. Biologically, boosting the production of T cells, aided by antibodies, is sufficient to disrupt the virus persistence equilibrium and stabilize the virus-free equilibrium. Consequently, the virus is cleared. Case VI: For the sixth scenario where *V* ⋆≠ 0, ![Graphic][128], *E*⋆≠ 0, ![Graphic][129], we obtain a boundary steady state of the form: ![Formula][130] The equilibrium state ![Graphic][131] will exists provided that *δ*2 *> σ*. Given the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][132] The eigenvalues of ![Graphic][133] are ![Formula][134] Referring to the eigenvalues of ![Graphic][135], this equilibrium is an unstable node. The presence of naive T cells and effector B cells, along with antibody support, is adequate to destabilize the virus persistence equilibrium and thereby facilitate virus clearance. Case VII: For the seventh scenario where *V* ⋆ ≠ 0, ![Graphic][136], *E*⋆ = 0, ![Graphic][137], we obtain a boundary steady state of the form: ![Formula][138] The equilibrium ![Graphic][139] will exists if ![Graphic][140]. Given the initial conditions for the state variables and parameter values in Table 1, we have ![Formula][141] The eigenvalues of the Jacobian matrix ![Graphic][142] are ![Formula][143] The eigenvalues of ![Graphic][144] indicate that this boundary steady state is an unstable spiral node. The presence of naive T cells and natural killer cells, along with antibody support, is sufficient to destabilize the virus persistence equilibrium and consequently lead to virus clearance. Case VIII: For the eighth scenario where *V* ⋆ *≠* 0, ![Graphic][145], *E*⋆ *≠* 0, ![Graphic][146], we obtain a boundary steady state of the form: ![Formula][147] Where ![Graphic][148], and ![Graphic][149]. This equilib-rium ![Graphic][150] will exists if ![Graphic][151] and ![Graphic][152]. Given the initial conditions for the state variables and parameter values in Table 1, we obtain ![Formula][153] The equilibrium point ![Graphic][154] has eigenvalues ![Formula][155] indicating that this boundary steady state is a stable saddle node. The absence of naive T cells results in a stable virus persistence equilibrium. Despite the production of antibodies and effector B cells, they are insufficient to disrupt the virus persistence equilibrium. This underscores the crucial role that T cells play in clearing hepatitis B virus infection. ## 6 Numerical Simulations and Sensitivity Analysis In this section, we present the numerical simulations and sensitivity analysis of the proposed model. We utilize values provided in Table 1 for the immune system response to the co-existence dynamics of HBV and liver cancer. To evaluate the usefulness and functionality of system (1) numerically, we use a set of reasonably approximated parameter values derived from published HBV article [1]. The following initial conditions: *T*(0) = 5 *×* 105, *T*1(0) = 0, *T*2(0) = 0, *M*(0) = 4 *×* 105, *M*1(0) = 4 *×* 105, *I**M* (0) = 0, *V* (0) = 300, *A*(0) = 0, *E*(0) = 0, *N**k*(0) = 0, *I*10(0) = 0, *F*1(0) = 0, and *F*2(0) = 0 and the parameters used in the numerical simulations are in Table 1. ### 6.1 Numerical Results at Virus-free Equilibrium State To showcase the various types of dynamical behaviors that model 1 can exhibit in different parameter regimes, we solve this system numerically using the baseline parameter values provided in Table 1. The results are displayed in Figures 2, 3, and 4. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F2) Figure 2: Simulation results showing the behaviour of the state variables at the virus-free equilibrium state. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F3) Figure 3: Simulation results showing the behaviour of the state variables at the virus-free equilibrium state. ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F4) Figure 4: Simulation results showing the behaviour of the state variables at the virus-free equilibrium state. Figures 2, 3, and 4 depict the immune response dynamics when there is no hepatitis B virus in the liver cell population. In this scenario, the initial viral growth leads to an increase in the numbers of antibodies (*A*), effector B cells (*E*), natural killer cells (*N**K*), cytotoxic T cells (*T*2), and both types of interferons. This results in the successful clearance of the HBV infection, after which type-1 interferons are also eliminated, and the system reaches a stable virus-free steady state. The absence of infected macrophages suggests that hepatitis B virus is not being transmitted, which consequently means that naive T cells will not engage in any lytic activity. Thus, the the T helper-1, cytotoxic T cells, antibodies, effector B cell, natural killer, interleukin-10 and interferon-gamma populations turns to zero as illustrated in Figures 2, 3 and 4 respectively. ### 6.2 Numerical Results at Virus Persistence Equilibrium State Figures 5, 6, and 7 illustrate the behavior profiles of the state variables of model 1 after infection with a few hepatitis B viruses, at the virus persistence equilibrium. From Figure 5(a), we observe the dynamics of naive T cells being fully produced with a total density of approximately 5 *×* 105. However, when HBV is introduced into the liver cells, the dynamics of naive T cells change. The presence of HBV, whose dynamics are depicted in Figure 6(a), causes a decrease in the total density of naive T cells. This decrease is due to the differentiation process they undergo to produce T-helper 1 (*T*1), cytotoxic T cells (*T*2), and interleukin-10 (*I*10) cells. ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F5) Figure 5: Simulation results showing the behaviour of the state variables at the virus persistence equilibrium state. ![Figure 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F6) Figure 6: Simulation results showing the behaviour of the state variables at the virus persistence equilibrium state. ![Figure 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F7) Figure 7: Simulation results showing the behaviour of the state variables at the virus persistence equilibrium state. We also observed that HBV and naive T cells have a negative correlation, meaning that a higher and more robust naive T cell density leads to a decrease in the HBV load. Consequently, with a constantly increasing HBV load, the rate at which naive T cells decrease accelerates, accompanied by successful production of both pro-inflammatory and anti-inflammatory cells, as shown in Figures 6(e), 6(f), and 7. From Figures 5(b) and 5(c), we examine the dynamics of T helper 1 and cytotoxic T cells, respectively. We notice an initial increase in both cell types due to naive T cell differentiation. However, we observed that cytotoxic T cells decrease immediately after they are produced within the shortest possible time due to the presence of HBV, which constantly inhibits their production. T helper 1 cells, on the other hand, decrease as a result of lytic activity caused by the replication rate of hepatitis B virus inside macrophages, impacting the overall rate of disease progression. Figures 5(d) and 5(e) show that uninfected macrophages without cancer cells and uninfected macrophages with cancer cells start with an initial density of approximately 4 *×* 105 and 7925, respectively, which decreases over time. Their decreasing nature is attributed to the increasing HBV load. As the replication rate of hepatitis B virus increases, more macrophages are required at the infection site to engulf the virus, as illustrated in Figure 5(f). This activity decreases the density of the uninfected macrophage population. Consequently, the number of infected macrophages increases as more uninfected macrophages engulf the virus and become infected. However, the rise in infected macrophages stabilizes once they become activated, as shown in Figure 5(f). Furthermore, the density of infected macrophages, which positively correlates with interferon gamma, interferon alpha and beta cytokines, antibodies, effector B cells, and natural killer cells, is illustrated in Figures 6(b), 6(c), and 6(d). As the concentration of interferon gamma increases, as shown in Figure 6(f), it triggers the activation of more infected macrophages, which in turn increases the population of uninfected macrophages. Thus, Figures 5, 6, and 7 depict the dynamics when the endemic steady state is feasible and stable. It can be observed that the initial viral growth is suppressed by the combined effects of different branches of the immune system. ### 6.3 Switching Time The results of HBV infection differ significantly among individuals, highlighting the importance of host genetic factors in determining susceptibility to HBV persistence and the progression of liver damage to cirrhosis and HCC. An effective antiviral response, primarily mediated by *CD*4+ and *CD*8+ T-cells, natural killer cells, and monocytes, can lead to immune-controlled HBV replication (functional cure). Conversely, in children and adults with weakened immune systems, active viral replication may persist. In its primary function, IL-10 acts as an immunosuppressive cytokine by inhibiting T-cell proliferation and the functions of antigen-presenting cells (APCs), as well as by modulating the synthesis of cytokines and chemokines [39]. Recently, a subset of IL-10 producing B-cells, known as regulatory B-cells (Bregs), has been shown to regulate HBV-specific *CD*8+ T-cell immunity [40, 41]. Down-regulation of IL-10 restores the function of exhausted HBV-specific *CD*8+ T-cells [40]. However, IL-10 may significantly impact the antiviral immune response, as it inhibits the production of pro-inflammatory cytokines such as IFN-*γ*, TNF-*α*, IL-1*β*, and IL-6. In chronic HBV infection, both the number of regulatory T-cells and the levels of inhibitory interleukin-10 (IL-10) and transforming growth factor beta (TGF-*β*) increase, leading to HBV-specific *CD*8+ T-cell exhaustion and making viral eradication from the liver impossible [42]. Thus, IL-10 expression is elevated during several chronic viral infections, serving as a viral strategy to down-regulate the host immune response and allow viral persistence in the host [43, 44, 45, 46, 47]. Over the past decade, numerous researchers have attempted to use mathematical models to predict when acute hepatitis B infection transitions to chronic infection and eventually to liver cancer; a period they have termed the “switching time.” We use system 1 to predict when T helper-1 (*T*1) cells will outnumber cytotoxic T cells (*T*2). The goal is to identify the conditions that lead an individual to progress from the acute stage to the chronic stage as the disease progresses. The period during which T helper-1 (*T*1) cells surpass cytotoxic T cells (*T*2) is considered crucial for an individual’s progression from the acute infection to the chronic stage. In the initial phase of the virus within infected macrophages, the virus tends to remain dormant while continuing to replicate. However, within the infected macrophages, there is a battle between the immune response and the hepatitis B virus. For the virus infection to reach the chronic stage, there must be a corresponding increase in viral load. Additionally, due to the complex nature of the hepatitis B virus, a higher viral load poses a significant threat to the host, potentially resulting in chronic disease. Therefore, the concept of switching time aims to identify the minimum time and factors that allow T helper-1 (*T*1) cells to quickly overcome cytotoxic T cells (*T*2), enabling infected macrophages to turn into acute macrophages for the immune response to take effect. Logically, a shorter switching time should correlate with a lower hepatitis B viral load. In our model, we observed that the parameters *ρ*1, *ρ*2, and *ρ*3 significantly impact the manipulation of switching time and viral load. These parameters are involved in the activities of interferon-alpha and beta, as well as interferon-gamma, and the production of pro-inflammatory cytokines such as IFN-*γ*, TNF-*α*, IL-1*β*, and IL-6. These factors play a significant role in the progression of the disease from the acute stage to chronic and liver cancer stages. Keeping all the other parameter values fixed and varying the value of *ρ*3, we obtained the following switching time: From Figures 8(a) to 8(d), we observed the various switching times for the model. Initially, the density of cytotoxic T cells increases in the first few days, after which it begins to decrease and is surpassed by the gradually increasing density of T helper-1 cells. For a parameter value of *ρ*3 = 2.03 *×* 10−7, we noticed that it takes over 100 days (approximately 112 days) for T helper-1 cells to surpass cytotoxic T cells. However, when we increase the parameter value of *ρ*3 from *ρ*3 = 2.03 *×* 10−7 to *ρ*3 = 2.03 *×* 10−6, *ρ*3 = 2.03 *×* 10−5, and *ρ*3 = 2.03 *×* 10−4, the switching time increases to 266, 457, and 600 days, respectively. This seems unexpected since one would usually expect that producing more interleukin-10 cytokines for the same viral load would help eliminate the infection. However, these findings support the argument that in HBV infection, ongoing antigen presentation by infected cells and exposure to high antigen levels are associated with *CD*8+ T cell exhaustion [48]. Biologically, the continuous production of interleukin-10 cytokines impairs cytotoxic T cells (*T*2), thus increasing the switching time. A switching time of approximately 112 days corresponds to the initial period after exposure to a core hepatitis B virus (i.e., the acute stage of HBV infection), which is the appropriate time for the immune system to respond. To further investigate and gain more understanding of the dynamics of HBV load as the switching time increases, we vary the parameter *ω* in relation to its switching time value and plot the hepatitis B virus compartment to observe the changes. From Figure 9, we observed that an increase in switching time corresponds to an increase in the viral load. For an *ω* value of 20.000000203, the density of the hepatitis B virus population is depicted by the blue line in the graph. When the value of *ω* is increased to 20.00000203, 20.0000203, and 20.000203, the density of the hepatitis B virus population is represented by the green, red, and magma lines, respectively. Generally, a lower switching time correlates with a lower hepatitis B viral load. Thus, from the graph, we noticed that a decrease in the value of *ω* from 20.000000203 to 20.0000000203, 20.00000000203, and 20.000000000203, which represents a lower switching time, corresponds to a lower hepatitis B viral load as illustrated by the black, cyan, and yellow dashed lines. This result is not surprising since we did not incorporate a time delay into the dynamics of our naive T cells. The lack of a time delay in the production of naive T cells allows for their continuous exposure to high levels of antigens, resulting in their impairment. Hence, the model predicts that a delay in the differentiation of the naive T cells into T helper 1 cells influences the determination of the switching time. ![Figure 8:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F8) Figure 8: The simulation graph for switching time. ![Figure 9:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F9) Figure 9: Graph illustrating the numerical solution that depicts the variation in HBV load with different switching times. ### 6.4 Sensitivity Analysis Due to uncertainties in the collection of data and estimated parameter values, the reproduction number ℝ is usually affected. We carry out sensitivity analysis to determine the relative importance of each epidemic parameter for transmission and control of both hepatitis B virus and liver cancer diseases and changes in the structure of the model. Latin Hypercube Sampling (LHS) and Partial Rank Correlation Coefficients (PRCC) are employed to identify model parameters that exert the greatest influence on the model, with the reproduction number (ℝ) acting as the dependent function. The main aim of this analysis is to ascertain the effects of parameters on model results. Parameters that are highly sensitive should be estimated with greater precision, as small changes in these parameters can lead to significant variations in the results [49, 50, 51]. Conversely, parameters that are not sensitive require less effort to estimate, since small changes in these parameters do not lead to large variations in the quantity of interest [50]. Parameters with PRCC values greater than positive 0.50 are considered to be highly positively correlated with the dependent function, while those with values less than negative 0.50 are considered to be highly negatively correlated with the dependent function [49, 50, 51]. The parameters included in the PRCC analysis are the interaction rate between *V* and *M* (*β*1), the interaction rate between *V* and *M*1 (*β*2), the supply rate of uninfected macrophages without cancer ![Graphic][156], the supply rate of uninfected macrophages with cancer ![Graphic][157], the rate at which infected macrophages develop into liver cancer (*ζ*), the replication rate of new infectious virions (*ω*), the decay rate of uninfected macrophages without cancer ![Graphic][158], the decay rate of uninfected macrophages with cancer ![Graphic][159], the virus decay rate (*µ**V*), the decay rate of infected macrophages ![Graphic][160], the virus lytic effect rate to release antibodies (*ψ*), and the lytic effect rate of infected macrophages to release interferon alpha and beta cytokines (*ε*). A PRCC analysis was conducted for five different periods; however, the parameters exhibited the same effect on the dependent function across all five periods. Therefore, we chose to present one plot, as shown in Figure 10. The results indicate that the twelve parameters that most significantly impact the response function (ℝ) are *β*1, *β*2, ![Graphic][161], *ζ*, ![Graphic][162], *µ**V*, ![Graphic][163], *ψ*, and *ε*. Base on the PRCC values, the parameters *β*1, *β*2, ![Graphic][164], and *ω* positively influence (ℝ), meaning that an increase (decrease) in these parameters will increase (decrease) (ℝ). In contrast, the parameters *ζ*, ![Graphic][165], *µ**V*, ![Graphic][166], *ψ*, and *ε* negatively influence (ℝ), and an increase in these parameters will decrease (ℝ). The results from PRCC analysis are summarized in Table 2. View this table: [Table 2:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/T2) Table 2: Parameters and their relationship with ℝ ![Figure 10:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/07/14/2024.07.14.24310388/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2024/07/14/2024.07.14.24310388/F10) Figure 10: Partial rank correlation coefficients (PRCCs) showing the impact of twelve model parameters on the reproduction number (ℝ) of the model. Parameter values used are as given in Table 1. ## 7 Discussion of results In this paper, we proposed and examined a new model for HBV infection and its coexistence with liver cancer, with a specific emphasis on the interactions between different branches of the immune system. This model includes the innate immune response, represented by natural killer cells (NK), and the adaptive immune response, represented by effector B cells and naive T cells, which differentiate into T helper-1 cells, cytotoxic T cells, interleukin-10, antibodies, and various cytokines. During infection, cytokines play a crucial role in recruiting both innate and adaptive immune components, boosting their efficacy, and enabling the non-cytolytic clearance of infected cells. The stability analysis of the steady states demonstrated how different parameters influence the dynamics of the immune response. Some findings, like the destabilization of the virus persistence equilibrium due to the presence of T cells, were intuitively obvious, while others were quite surprising. Naturally, increasing the number of NK cells, the rate at which antibodies clear free virus, the rate at which effector B cells eliminate infected cells, the rate at which cytotoxic T lymphocytes clear infected cells, and the rate at which IFN-*γ* inhibits viral production all contribute to more effective infection clearance, stabilizing the disease-free steady-state solutions. However, we noticed that for a very small or a very large rate of free virus clearance by antibodies, NK cells, and effector B cells, the stability of the endemic steady state is unaffected by how quickly the new antibodies, NK cells, and effector B cells are produced. But whenever the production of naive T cells increases, the virus-persistence equilibrium is destabilized. This result is very astonishing, as one would generally presume that a higher production rate of antibodies, NK cells, and effector B cells would lead to the clearance of the infection rather than the stabilization of a chronic state. This indicates that it is not the distinct production rates of antibodies, NK cells, and effector B cells that lead to viral clearance, but rather the ratio between these rates that determines whether the system maintains a chronic infection, exhibits periodic oscillations, or achieves viral clearance. Interleukin-10 (IL-10) is a cytokine with anti-inflammatory properties that plays a complex role in the immune response to hepatitis B virus (HBV) infection. IL-10 is primarily recognized for its ability to suppress the immune response. It is well known that IL-10 inhibits the production of pro-inflammatory cytokines such as IL-1, IL-6, TNF-*α*, and interferon-gamma (IFN-*γ*). By regulating these cytokines, IL-10 helps to decrease the overall inflammatory environment in the liver during HBV infection. IL-10 helps maintain a balance between necessary immune activity to control the virus and the need to limit excessive inflammation that can cause liver damage. This balance is crucial because excessive immune responses can lead to conditions like fulminant hepatitis, while inadequate responses can result in chronic infection. During HBV infection, IL-10 can suppress the activity of various immune cells, including T cells, macrophages, and dendritic cells. This suppression can reduce liver inflammation and damage, potentially preventing severe liver injury. Despite the significant role interleukin-10 plays in HBV infection, it can potentially contribute to viral persistence by dampening the immune system’s ability to clear the virus. By inhibiting the immune response, IL-10 may aid in the persistence of HBV. The suppression of cytotoxic T lymphocytes (CTLs), which are crucial for clearing infected cells, can allow the virus to evade the immune system and establish chronic infection as illustrated in Figure 9. Elevated levels of IL-10 have been associated with poor responses to antiviral therapies in HBV-infected individuals. High IL-10 levels can inhibit the effectiveness of the immune response needed for successful treatment. We observed various switching times for the model. Initially, the density of cytotoxic T cells increases in the first few days, after which it begins to decrease and is surpassed by the gradually increasing density of T helper-1 cells. However, when we increased the production rate of interleukin-10, the switching time increased. This seems unexpected since one would usually expect that producing more interleukin-10 cytokines for the same viral load would help eliminate the infection. Thus, these findings support the argument that in HBV infection, ongoing antigen presentation by infected cells and exposure to high antigen levels are associated with *CD*8+ T cell exhaustion. Biologically, the continuous production of interleukin-10 cytokines impairs cytotoxic T cells (*T*2), thus increasing the switching time. A switching time of approximately 112 days corresponds to the initial period after exposure to the core hepatitis B virus (i.e., the acute stage of HBV infection), which is the appropriate time for the immune system to respond. However, further production of interleukin-10 beyond this period will lead to immune impairment. The model further showed that under specific conditions (see Figures 8 and 9), the switching time can be adjusted. As illustrated in Figure 9, we observed that a reduction in the value of *ω*, indicating a shorter switching time, corresponds to a lower hepatitis B viral load, as shown by the black, cyan, and yellow dashed lines. This outcome is expected since a time delay was not included in the dynamics of our naive T cells. Consequently, the model suggests that a delay in the differentiation of naive T cells into T helper 1 cells affects the determination of the switching time. Therefore, it is reasonable to state that a positive lifestyle correlates positively with immune system functioning and is crucial in determining whether an acute infection progresses to the chronic stage. The implications of the model align with our expectations; logically, a longer switching time results in a higher viral load and vice versa. Minimizing the side effects of IL-10 in HBV infection requires carefully managing the immune response to ensure effective viral clearance while preventing excessive immune suppression that could lead to chronic infection or liver damage. Therefore, we recommend developing drugs that selectively inhibit IL-10 activity in specific cell types or tissues to help minimize systemic side effects. Additionally, using antibodies or small molecules to block IL-10 receptors on specific immune cells might reduce the suppression of antiviral responses while preserving the beneficial anti-inflammatory effects in the liver. Combining IL-10 modulation with potent antiviral drugs can lower the viral load, reducing the need for immune suppression. Furthermore, pairing IL-10 inhibitors with other immune modulators that enhance T cell activity or other antiviral immune responses can help balance immune activation and suppression. We have used arbitrary parameter values for simulating the model; it would be better if clinical data sets were used to fit the model. Time delay in the immune response to HBV infection also plays a crucial role in immune response dynamics. Thus, future research can extend this work by incorporating the time-delay factor in the immune response. Another refinement of the model could include different immunotherapy control techniques, such as antibody therapies currently in clinical trials, CAR-T cell therapy, molecular therapy, and exhaustion therapy. These therapies are extensively studied in [12, 21, 52, 53, 54] to find optimal treatment strategies for HBV and liver cancer management. ## Data Availability All data produced in the present work are contained in the manuscript * Received July 14, 2024. * Revision received July 14, 2024. * Accepted July 14, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Ciupe, S.M., Ribeiro, R.M., Perelson, A.S.: Antibody responses during hepatitis b viral infection. PLoS computational biology 10(7), 1003730 (2014) 2. [2].Guidotti, L.G., Chisari, F.V.: Noncytolytic control of viral infections by the innate and adaptive immuneresponse. Annual review of immunology 19(1), 65–91 (2001) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.immunol.19.1.65&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11244031&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168234600004&link_type=ISI) 3. [3].Ciupe, S.M., Ribeiro, R.M., Nelson, P.W., Dusheiko, G., Perelson, A.S.: The role of cells refractory to productive infection in acute hepatitis b viral dynamics. Proceedings of the National Academy of Sciences 104(12), 5050–5055 (2007) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA0LzEyLzUwNTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNC8yMDI0LjA3LjE0LjI0MzEwMzg4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 4. [4].Babiker, Z.O.E., Hogan, C., Ustianowski, A., Wilkins, E.: Does interferon-sparing tenofovir disoproxil fumarate-based therapy have a role in the management of severe acute hepatitis delta superinfection? Journal of medical microbiology 61(12), 1780–1783 (2012) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/jmm.0.046649-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22956751&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 5. [5].Nowak, M.A., Bonhoeffer, S., Hill, A.M., Boehme, R., Thomas, H.C., McDade, H.: Viral dynamics in hepatitis b virus infection. Proceedings of the National Academy of Sciences 93(9), 4398–4402 (1996) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo5OiI5My85LzQzOTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNC8yMDI0LjA3LjE0LjI0MzEwMzg4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 6. [6].DeVico, A.L., Gallo, R.C.: Control of hiv-1 infection by soluble factors of the immune response. Nature Reviews Microbiology 2(5), 401–413 (2004) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrmicro878&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15100693&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221589700014&link_type=ISI) 7. [7].Su, Y., Wen, Y., Min, L.: Analysis of a hbv infection model with alt. In: 2012 IEEE 6th International Conference on Systems Biology (ISB), pp. 97–100 (2012). IEEE 8. [8].Busca, A., Kumar, A.: Innate immune responses in hepatitis b virus (hbv) infection. Virology journal 11, 1–8 (2014) 9. [9].Guidotti, L.G., Guilhot, S., Chisari, F.V.: Interleukin-2 and alpha/beta interferon down-regulate hepatitis b virus gene expression in vivo by tumor necrosis factor-dependent and-independent pathways. Journal of Virology 68(3), 1265–1270 (1994) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjk6IjY4LzMvMTI2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzE0LzIwMjQuMDcuMTQuMjQzMTAzODguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 10. [10].Herbein, G., O’brien, W.A.: Tumor necrosis factor (tnf)–α and tnf receptors in viral pathogenesis. Proceedings of the Society for Experimental Biology and Medicine: Minireviews 223(3), 241–257 (2000) 11. [11].Vos, T., Barber, R.M., Bell, B., Bertozzi-Villa, A., Biryukov, S., Bolliger, I., Charlson, F., Davis, A., Degenhardt, L., Dicker, D., et al: Global, regional, and national incidence, prevalence, and years lived with disability for 301 acute and chronic diseases and injuries in 188 countries, 1990–2013: a systematic analysis for the global burden of disease study 2013. The lancet 386(9995), 743–800 (2015) 12. [12].Fatehi, F., Bingham, R.J., Stockley, P.G., Twarock, R.: An age-structured model of hepatitis b viral infection highlights the potential of different therapeutic strategies. Scientific Reports 12(1), 1252 (2022) 13. [13].Sheena, B.S., Hiebert, L., Han, H., Ippolito, H., Abbasi-Kangevari, M., Abbasi-Kangevari, Z., Abbastabar, H., Abdoli, A., Ali, H.A., Adane, M.M., et al: Global, regional, and national burden of hepatitis b, 1990– 2019: a systematic analysis for the global burden of disease study 2019. The lancet Gastroenterology & hepatology 7(9), 796–829 (2022) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-1253(22)00124-8&link_type=DOI) 14. [14].Ciupe, S.M., Ribeiro, R.M., Nelson, P.W., Perelson, A.S.: Modeling the mechanisms of acute hepatitis b virus infection. Journal of theoretical biology 247(1), 23–35 (2007) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2007.02.017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17428501&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247490000003&link_type=ISI) 15. [15].Min, L., Su, Y., Kuang, Y.: Mathematical analysis of a basic virus infection model with application to hbv infection. The Rocky Mountain Journal of Mathematics, 1573–1585 (2008) 16. [16].Gourley, S.A., Kuang, Y., Nagy, J.D.: Dynamics of a delay differential equation model of hepatitis b virus infection. Journal of Biological Dynamics 2(2), 140–153 (2008) 17. [17].Hews, S., Eikenberry, S., Nagy, J.D., Kuang, Y.: Rich dynamics of a hepatitis b viral infection model with logistic hepatocyte growth. Journal of Mathematical Biology 60, 573–590 (2010) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00285-009-0278-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19533136&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 18. [18].Yousfi, N., Hattaf, K., Tridane, A.: Modeling the adaptive immune response in hbv infection. Journal of mathematical biology 63, 933–957 (2011) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00285-010-0397-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21234570&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 19. [19].Perelson, A.S., Neumann, A.U., Markowitz, M., Leonard, J.M., Ho, D.D.: Hiv-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science 271(5255), 1582–1586 (1996) [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyNzEvNTI1NS8xNTgyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTQvMjAyNC4wNy4xNC4yNDMxMDM4OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. [20].Hu, J., Liu, K.: Complete and incomplete hepatitis b virus particles: formation, function, and application. Viruses 9(3), 56 (2017) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/v9030056&link_type=DOI) 21. [21].Chenar, F.F., Kyrychko, Y., Blyuss, K.: Mathematical model of immune response to hepatitis b. Journal of theoretical biology 447, 98–110 (2018) 22. [22].Hassuneh, M.R., Nagarkatti, M., Nagarkatti, P.S.: Role of interleukin-10 in the regulation of tumorigenicity of a t cell lymphoma. Leukemia & lymphoma 54(4), 827–834 (2013) 23. [23].Pawelek, K.A., Huynh, G.T., Quinlivan, M., Cullinane, A., Rong, L., Perelson, A.S.: Modeling withinhost dynamics of influenza virus infection including immune responses. PLoS computational biology 8(6), 1002588 (2012) 24. [24].Abbas, A., Lichtman, A., Pillai, S.: Cellular and Molecular Immunology E-book. Elsevier Health Sciences, ??? (2014) 25. [25].Diekmann, O., Heesterbeek, J.A.P.: Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation vol. 5. John Wiley & Sons, ??? (2000) 26. [26].Chataa, P., Nyabadza, F., Naandam, S.M.: Mathematical modelling of the transmission dynamics of hepatitis b virus in the presence of imperfect vaccination. J. Math. Comput. Sci. 11(6), 6949–6979 (2021) 27. [27].Castillo-Chavez, C., Song, B.: Dynamical models of tuberculosis and their applications. Mathematical Bio-sciences and Engineering 1(2), 361 (2004) 28. [28].Opoku, N.K.-D.O., Mazandu, G.K.: Modelling the human immune response dynamics during progression from mycobacterium latent infection to disease. Applied Mathematical Modelling 80, 217–237 (2020) 29. [29].Goyal, A., Ribeiro, R.M., Perelson, A.S.: The role of infected cell proliferation in the clearance of acute hbv infection in humans. Viruses 9(11), 350 (2017) 30. [30].Murray, J.M., Goyal, A.: In silico single cell dynamics of hepatitis b virus infection and clearance. Journal of Theoretical Biology 366, 91–102 (2015) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2014.11.020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25476731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 31. [31].Wigginton, J.E., Kirschner, D.: A model to predict cell-mediated immune regulatory mechanisms during human infection with mycobacterium tuberculosis. The Journal of Immunology 166(3), 1951–1967 (2001) 32. [32].Agarwal, M., Bhadauria, A.: Correlation between hepatitis and cancer: A mathematical model 33. [33].Guang, Y., Yuzhong, L., Hui, L.: Establishment of an analysis model based on measurement of hepatitis b viral infection serum markers. BMC Infectious Diseases 19, 1–5 (2019) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-019-4277-8&link_type=DOI) 34. [34].Zhang, Y., Wallace, D.L., De Lara, C.M., Ghattas, H., Asquith, B., Worth, A., Griffin, G.E., Taylor, G.P., Tough, D.F., Beverley, P.C., et al: In vivo kinetics of human natural killer cells: the effects of ageing and acute and chronic viral infection. Immunology 121(2), 258–265 (2007) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2567.2007.02573.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17346281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245988600013&link_type=ISI) 35. [35].Jamieson, A.M., Isnard, P., Dorfman, J.R., Coles, M.C., Raulet, D.H.: Turnover and proliferation of nk cells in steady state and lymphopenic conditions. The Journal of Immunology 172(2), 864–870 (2004) 36. [36].Guedj, J., Dahari, H., Rong, L., Sansone, N.D., Nettles, R.E., Cotler, S.J., Layden, T.J., Uprichard, S.L., Perelson, A.S.: Modeling shows that the ns5a inhibitor daclatasvir has two modes of action and yields a shorter estimate of the hepatitis c virus half-life. Proceedings of the National Academy of Sciences 110(10), 3991–3996 (2013) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEwLzEwLzM5OTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNC8yMDI0LjA3LjE0LjI0MzEwMzg4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 37. [37].Kitagawa, K., Nakaoka, S., Asai, Y., Watashi, K., Iwami, S.: A pde multiscale model of hepatitis c virus infection can be transformed to a system of odes. Journal of theoretical biology 448, 80–85 (2018) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2018.04.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29634960&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 38. [38].Schlub, T.E., Sun, J.C., Walton, S.M., Robbins, S.H., Pinto, A.K., Munks, M.W., Hill, A.B., Brossay, L., Oxenius, A., Davenport, M.P.: Comparing the kinetics of nk cells, cd4, and cd8 t cells in murine cytomegalovirus infection. The journal of Immunology 187(3), 1385–1392 (2011) 39. [39].Saraiva, M., O’garra, A.: The regulation of il-10 production by immune cells. Nature reviews immunology 10(3), 170–181 (2010) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri2711&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20154735&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274975600002&link_type=ISI) 40. [40].Das, A., Ellis, G., Pallant, C., Lopes, A.R., Khanna, P., Peppa, D., Chen, A., Blair, P., Dusheiko, G., Gill, U., et al: Il-10–producing regulatory b cells in the pathogenesis of chronic hepatitis b virus infection. The Journal of Immunology 189(8), 3925–3935 (2012) 41. [41].Liu, Y., Cheng, L.-s., Wu, S.-d., Wang, S.-q., Li, L., She, W.-m., Li, J., Wang, J.-y., Jiang, W.: Il-10-producing regulatory b-cells suppressed effector t-cells but enhanced regulatory t-cells in chronic hbv infection. Clinical Science 130(11), 907–919 (2016) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToicHBjbGluc2NpIjtzOjU6InJlc2lkIjtzOjEwOiIxMzAvMTEvOTA3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTQvMjAyNC4wNy4xNC4yNDMxMDM4OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 42. [42].Peeridogaheh, H., Meshkat, Z., Habibzadeh, S., Arzanlou, M., Shahi, J.M., Rostami, S., Gerayli, S., Teimourpour, R.: Current concepts on immunopathogenesis of hepatitis b virus infection. Virus research 245, 29–43 (2018) 43. [43].Hyodo, N., Nakamura, I., Imawari, M.: Hepatitis b core antigen stimulates interleukin-10 secretion by both t cells and monocytes from peripheral blood of patients with chronic hepatitis b virus infection. Clinical & Experimental Immunology 135(3), 462–466 (2004) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2249.2003.02376.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15008979&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) 44. [44].Ohga, S., Nomura, A., Takada, H., Tanaka, T., Furuno, K., Takahata, Y., Kinukawa, N., Fukushima, N., Imai, S., Hara, T.: Dominant expression of interleukin-10 and transforming growth factor-β genes in activated t-cells of chronic active epstein–barr virus infection. Journal of medical virology 74(3), 449–458 (2004) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.20197&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15368517&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224009800012&link_type=ISI) 45. [45].Brooks, D.G., Trifilo, M.J., Edelmann, K.H., Teyton, L., McGavern, D.B., Oldstone, M.B.: Interleukin-10 determines viral clearance or persistence in vivo. Nature medicine 12(11), 1301–1309 (2006) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1492&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17041596&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241844900028&link_type=ISI) 46. [46].Kaplan, D.E., Ikeda, F., Li, Y., Nakamoto, N., Ganesan, S., Valiga, M.E., Nunes, F.A., Reddy, K.R., Chang, K.-M.: Peripheral virus-specific t-cell interleukin-10 responses develop early in acute hepatitis c infection and become dominant in chronic hepatitis. Journal of hepatology 48(6), 903–913 (2008) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2008.01.030&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18384906&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F14%2F2024.07.14.24310388.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256683100004&link_type=ISI) 47. [47].Brockman, M.A., Kwon, D.S., Tighe, D.P., Pavlik, D.F., Rosato, P.C., Sela, J., Porichis, F., Le Gall, S., Waring, M.T., Moss, K., et al: Il-10 is up-regulated in multiple cell types during viremic hiv infection and reversibly inhibits virus-specific t cells. Blood, The Journal of the American Society of Hematology 114(2), 346–356 (2009) 48. [48].Revill, P.A., Chisari, F.V., Block, J.M., Dandri, M., Gehring, A.J., Guo, H., Hu, J., Kramvis, A., Lampertico, P., Janssen, H.L., et al: A global scientific strategy to cure hepatitis b. The lancet Gastroenterology & hepatology 4(7), 545–558 (2019) 49. [49].Cariboni, J., Gatelli, D., Liska, R., Saltelli, A.: The role of sensitivity analysis in ecological modelling. Ecological modelling 203(1-2), 167–182 (2007) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ecolmodel.2005.10.045&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246020300017&link_type=ISI) 50. [50].Blower, S.M., Dowlatabadi, H.: Sensitivity and uncertainty analysis of complex models of disease transmission: an hiv model, as an example. International Statistical Review/Revue Internationale de Statistique, 229–243 (1994) 51. [51].Iboi, E.A., Gumel, A.B.: Mathematical assessment of the role of dengvaxia vaccine on the transmission dynamics of dengue serotypes. Mathematical biosciences 304, 25–47 (2018) 52. [52].Abernathy, Z., Abernathy, K., Stevens, J., et al: A mathematical model for tumor growth and treatment using virotherapy. AIMS Math 5(5), 4136–4150 (2020) 53. [53].Cho, H., Levy, D.: The impact of competition between cancer cells and healthy cells on optimal drug delivery. Mathematical Modelling of Natural Phenomena 15, 42 (2020) 54. [54].Bukkuri, A.: Optimal control analysis of combined chemotherapy-immunotherapy treatment regimens in a pkpd cancer evolution model. Biomath 9(1), 2002137 (2020) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: /embed/inline-graphic-6.gif [7]: /embed/inline-graphic-7.gif [8]: /embed/inline-graphic-8.gif [9]: /embed/inline-graphic-9.gif [10]: /embed/inline-graphic-10.gif [11]: /embed/inline-graphic-11.gif [12]: /embed/inline-graphic-12.gif [13]: /embed/inline-graphic-13.gif [14]: /embed/inline-graphic-14.gif [15]: /embed/graphic-2.gif [16]: /embed/graphic-3.gif [17]: /embed/graphic-4.gif [18]: /embed/graphic-5.gif [19]: /embed/graphic-6.gif [20]: /embed/inline-graphic-15.gif [21]: /embed/graphic-7.gif [22]: /embed/graphic-8.gif [23]: /embed/graphic-9.gif [24]: /embed/graphic-10.gif [25]: /embed/graphic-11.gif [26]: /embed/graphic-12.gif [27]: /embed/graphic-13.gif [28]: /embed/graphic-14.gif [29]: /embed/graphic-15.gif [30]: /embed/graphic-16.gif [31]: /embed/inline-graphic-16.gif [32]: /embed/inline-graphic-17.gif [33]: /embed/graphic-17.gif [34]: /embed/graphic-18.gif [35]: /embed/graphic-19.gif [36]: /embed/graphic-20.gif [37]: /embed/graphic-21.gif [38]: /embed/graphic-22.gif [39]: /embed/graphic-23.gif [40]: /embed/graphic-24.gif [41]: /embed/graphic-25.gif [42]: /embed/graphic-26.gif [43]: /embed/graphic-27.gif [44]: /embed/graphic-28.gif [45]: /embed/graphic-29.gif [46]: /embed/graphic-30.gif [47]: /embed/inline-graphic-18.gif [48]: /embed/inline-graphic-19.gif [49]: /embed/inline-graphic-20.gif [50]: /embed/inline-graphic-21.gif [51]: /embed/inline-graphic-22.gif [52]: /embed/inline-graphic-23.gif [53]: /embed/graphic-31.gif [54]: /embed/inline-graphic-24.gif [55]: /embed/inline-graphic-25.gif [56]: /embed/graphic-32.gif [57]: /embed/graphic-33.gif [58]: /embed/inline-graphic-26.gif [59]: /embed/graphic-34.gif [60]: /embed/inline-graphic-27.gif [61]: /embed/graphic-35.gif [62]: /embed/graphic-36.gif [63]: /embed/graphic-37.gif [64]: /embed/inline-graphic-28.gif [65]: /embed/inline-graphic-29.gif [66]: /embed/inline-graphic-30.gif [67]: /embed/inline-graphic-31.gif [68]: /embed/graphic-38.gif [69]: /embed/graphic-39.gif [70]: /embed/graphic-40.gif [71]: /embed/graphic-41.gif [72]: /embed/inline-graphic-32.gif [73]: /embed/inline-graphic-33.gif [74]: /embed/inline-graphic-34.gif [75]: /embed/inline-graphic-35.gif [76]: /embed/inline-graphic-36.gif [77]: /embed/inline-graphic-37.gif [78]: /embed/graphic-42.gif [79]: /embed/inline-graphic-38.gif [80]: /embed/inline-graphic-39.gif [81]: /embed/inline-graphic-40.gif [82]: /embed/inline-graphic-41.gif [83]: /embed/inline-graphic-42.gif [84]: /embed/inline-graphic-43.gif [85]: /embed/graphic-43.gif [86]: /embed/inline-graphic-44.gif [87]: /embed/inline-graphic-45.gif [88]: /embed/inline-graphic-46.gif [89]: /embed/inline-graphic-47.gif [90]: /embed/inline-graphic-48.gif [91]: /embed/inline-graphic-49.gif [92]: /embed/inline-graphic-50.gif [93]: /embed/graphic-44.gif [94]: /embed/inline-graphic-51.gif [95]: /embed/graphic-46.gif [96]: /embed/graphic-47.gif [97]: /embed/inline-graphic-52.gif [98]: /embed/graphic-48.gif [99]: /embed/inline-graphic-53.gif [100]: /embed/inline-graphic-54.gif [101]: /embed/graphic-49.gif [102]: /embed/graphic-50.gif [103]: /embed/inline-graphic-55.gif [104]: /embed/graphic-51.gif [105]: /embed/inline-graphic-56.gif [106]: /embed/inline-graphic-57.gif [107]: /embed/inline-graphic-58.gif [108]: /embed/graphic-52.gif [109]: /embed/graphic-53.gif [110]: /embed/inline-graphic-59.gif [111]: /embed/graphic-54.gif [112]: /embed/inline-graphic-60.gif [113]: /embed/inline-graphic-61.gif [114]: /embed/inline-graphic-62.gif [115]: /embed/inline-graphic-63.gif [116]: /embed/graphic-55.gif [117]: /embed/graphic-56.gif [118]: /embed/inline-graphic-64.gif [119]: /embed/graphic-57.gif [120]: /embed/inline-graphic-65.gif [121]: /embed/inline-graphic-66.gif [122]: /embed/inline-graphic-67.gif [123]: /embed/graphic-58.gif [124]: /embed/graphic-59.gif [125]: /embed/inline-graphic-68.gif [126]: /embed/graphic-60.gif [127]: /embed/inline-graphic-69.gif [128]: /embed/inline-graphic-70.gif [129]: /embed/inline-graphic-71.gif [130]: /embed/graphic-61.gif [131]: /embed/inline-graphic-72.gif [132]: /embed/graphic-62.gif [133]: /embed/inline-graphic-73.gif [134]: /embed/graphic-63.gif [135]: /embed/inline-graphic-74.gif [136]: /embed/inline-graphic-75.gif [137]: /embed/inline-graphic-76.gif [138]: /embed/graphic-64.gif [139]: /embed/inline-graphic-77.gif [140]: /embed/inline-graphic-78.gif [141]: /embed/graphic-65.gif [142]: /embed/inline-graphic-79.gif [143]: /embed/graphic-66.gif [144]: /embed/inline-graphic-80.gif [145]: /embed/inline-graphic-81.gif [146]: /embed/inline-graphic-82.gif [147]: /embed/graphic-67.gif [148]: /embed/inline-graphic-83.gif [149]: /embed/inline-graphic-84.gif [150]: /embed/inline-graphic-85.gif [151]: /embed/inline-graphic-86.gif [152]: /embed/inline-graphic-87.gif [153]: /embed/graphic-68.gif [154]: /embed/inline-graphic-88.gif [155]: /embed/graphic-69.gif [156]: /embed/inline-graphic-89.gif [157]: /embed/inline-graphic-90.gif [158]: /embed/inline-graphic-91.gif [159]: /embed/inline-graphic-92.gif [160]: /embed/inline-graphic-93.gif [161]: /embed/inline-graphic-94.gif [162]: /embed/inline-graphic-95.gif [163]: /embed/inline-graphic-96.gif [164]: /embed/inline-graphic-97.gif [165]: /embed/inline-graphic-98.gif [166]: /embed/inline-graphic-99.gif