Stability and Optimal Control of Tuberculosis Spread with an Imperfect Vaccine in the Case of Co-Infection with HIV

This paper focuses on the study and control of a non-linear mathematical epidemic model ( vih SS VELI ) based on a system of ordinary differential equation modeling the spread of tuberculosis infectious with HIV/AIDS coinfection. Existence of both disease free equilibrium and endemic equilibrium is discussed. Reproduction number 0  is determined. Using Lyapunov-Lasalle methods, we analyze the stability of epidemic system around the equilibriums (disease free and endemic equilibrium). The global asymptotic stability of the disease free equilibrium whenever 1 vac <  is proved, where 0  is the reproduction number. We prove also that when 0  is less than one, tuberculosis can be eradicated. Numerical simulations are conducted to approve analytic results. To achieve control of the disease, seeking to reduce the infectious group by the minimum vaccine coverage, a control problem is formulated. The Pontryagin’s maximum principle is used to characterize the optimal control. The optimality system is derived and solved numerically using the Runge Kutta fourth procedure.


Introduction
Tuberculosis is one of the top 10 causes of death in the world [1]. In Cameroon, tuberculosis and HIV infection account for 25% of the morbid burden. TB is caused by bacteria (Mycobacterium tuberculosis) that most often affects the lungs. Tuberculosis is curable and preventable. TB is spread from person to person through the air. When people with lung TB cough, sneeze or spit, they propel the TB germs into the air. A person needs to inhale only a few of these germs to become infected. About one-third of the world population has latent TB, which means people have been infected by TB bacteria but are not (yet) ill with the disease and cannot transmit the disease. People infected with TB bacteria have a 10% lifetime risk of falling ill with TB. However, persons with compromised immune systems, such as people living with HIV, malnutrition or diabetes, or people who use tobacco, have a much higher risk of falling ill. When a person develops active TB disease, the symptoms (such as cough, fever, night sweats, or weight loss) may be mild for many months. This can lead to delays in seeking care, and results in transmission of the bacteria to others. People with active TB can infect 10 other people through close contact over the course of a year. The population of Cameroon is estimated in 2017 to 24 million inhabitants.
In the same year, 24,905 cases of tuberculosis were reported. The medical coverage is estimated at 52%. The death rate due to the tuberculose is approximatively 0.29% [2]. In 2018, the national budget allocated to the control of the tuberculosis is 10 million American dollars, and distributed as follows: 14% of internal source, 27% of external source and 59% of deficit.
The theme "End the global TB" has been decided by the world health organization and it covers the period 2016-2035 and the overall goal is to end the global TB epidemic. As part of the necessary multidisciplinary research approach, mathematical models have been extensively used to provide a framework for understanding tuberculosis transmission dynamics and control strategies of the infection spread in the host population [3] [4] [5] [6] [7]. Arthur et al. [6] study social and cultural factors in the successful control of tuberculosis. Optimal control theory has been applied to some biological models (VIH/AIDS treatment, Cancer, Fish harvestime) [8]- [14] and transmission model. Some mathematical models have been used to control tuberculosis. In [15], optimal control theory is applied to investigate optimal strategies for controlling the spread of tuberculosis using treatment of infected individuals with TB as the system control variables. In [16], the tuberculosis control is formulated and solved as an optimal control theory problem, indicating how a control term on the chemoprophylaxis should be introduced in the population to reduce the number of individuals with active TB. The papers [17] consider the optimal control of tuberculosis through education, diagnosis campaign and chemoprophylaxis of latently infected. In [6], it talks about the Impact of an effective multidrug-resistant tuberculosis control programme in the setting of an immature HIV epidemic.
This paper deals with the stability analysis of an , vih SVS ELI transmission model and uses optimal control technique to find and evaluate the impact of a mass vaccination schedule in the spread of TB/VIH coinfection. Individuals are Open Journal of Modelling and Simulation classified as one of susceptible (S) V vaccinated, AIDS patients ( vih S ), earlier latent (E) late latent (L) infectious or tuberculous (I), but allow that susceptible individuals may be given an imperfect vaccine that reduces their susceptibility to the disease, the V-compartment of vaccinated individuals is considered [18] as a susceptible compartment. AIDS patients are considered as another type of susceptible with hight level of susceptibility due to the fact that their immune system are compromised.
Since V and vih S are considered as the susceptible compartments, thus we are dealing with a differential susceptibility system with bilinear mass action as in Hyman and Li [19]. However, we include one-way flow between these two compartments, to denote the prevalence of AIDS in community. For analysis and stability issues;We use Lyapunov-LaSalle methods, we fully resolve the global dynamics of the model for the full parameter space. We demonstrate that the model exhibits threshold behavior with a globally stable disease-free equilibrium if the basic reproduction number is less than unity and a globally stable endemic equilibrium if the basic reproduction number is greater than unity. Our goal is to minimise number of TB infectious persons and the overall cost of the vaccine during a fixed period.
For numerical simulation, the programs used in this paper are designed so that no knowledge of MATLAB is required. For the control problem, there is a user-friendly interface that will guide you through. We have two different MATLAB programs, plotTB.m and codeTB.m. Here, .m is the extension given to all files intended for use in MATLAB.The file codeTB.m is the Runge-Kutta based, forward-backward sweep solver. It takes as input the values of the various parameters in the problem and outputs the solution to the optimality system. The ?le plotTB.m is the user friendly interface. It will ask you to enter the values of the parameters one by one, compile codeTB.m with these values, and plot the resulting solutions. All the files must be in the directory that MATLAB treats as the home directory. This is usually the Work directory. This paper is organized as follows.
In next section, model is described, in Section 3 we investigate stability analysis for the ( vih SS VELI ) epidemic model in this section We followed the methods of Nkamba, Leontine Nkague et al. 2019 [20] A control system for the optimality and its existence,and the optimal control are derived in Section 4. In Section 5, utilizing the representation of the optimal control, we describe a numerical solution of the optimality system consisting of the original state system, the adjoint system, and their boundary conditions. In Section 6, we describe, in detail, a real application of our optimal control theory. Finally, we conclude in Section 7.

vih SS VELI , Preliminaries
When first infected with TB bacteria, a person typically goes through a latent, asymptomatic and non-infectious period during which the body's immune Open Journal of Modelling and Simulation system fights the TB bacteria. There are two distinct stages of the latent TB infection. During the first two years, the risk of developing active disease is much higher, whereas during the later stage, the progression to active disease is much slower.
Compartmental modeling is used among epidemiologists to simulate disease dynamics. These models treat each disease state as a different compartment that contains a homogeneous population of individuals. Using a compartmental approach, the total host population can be partitioned into seven compartments:  µ , E µ , L µ and I µ respectively. Once infected, individuals progress through the early latent stage with an average rate ω . A fraction ;0 1 p p < ≤ ; of these individuals progress directly to the active TB stage, and the remaining 1 p − fraction progresses to the late latent stage. Once there, the rate of progression to active disease is at a lower rate ν . The recruitment makes respectively into the susceptible class, the vaccinated class, the VIH/patient class with the constant rate 1 π 2 π and 3 π . α is the vaccination coverage rate. 2 π is the recruitment of vaccinated a few day after birth, so we suppose that immunity is passed during the birth. 3 π is VIH/AIDS vertical transmission recruitment, it's means some peoples born with VIH/AIDS infection.
The dynamical transfer among the seven compartments is depicted in the transfer diagram ( Figure 1).

Description of Variables and Parameters
All parameters described in Table 1 are assumed to be positive.

vih SS VELT Model
Our model consists of the following system of ordinary differential equations: Cost of immune deficiency rate

A Compact Positively Invariant Absorbing Set
In order that the model be well-posed, it is necessary that the state variables

Existence of an Disease Free Equilibrium (DFE)
It is easy to check that model 1 always has the disease-free equilibrium In order to assume that vaccinated people don't produced infected more than susceptible people we should have Additionally, an endemic equilibrium

Basic Reproduction Ratio
Using the method of James Watmouth and all the next generation matrix [21] Open Journal of Modelling and Simulation [22], the basic reproduction number 0  is giving by We will see in the Section 3.4 theorem 1 that, when 0  is less than unity, infection can disappear in the population. Numerical simulations will confirm our results.

Stability of Disease Free Equilibrium (DFE)
In this section, we show that the disease-free equilibrium 0 P is globally asymptotically stable with respect to if 0 1 R ≤ ; and 0 P is unstable if 0 1 R > : , then the disease-free equilibrium is globally asymptotically stable.
Proof. Consider a Lyapunov function, Direct calculation leads to µ ω µ ν µ τ β θ β θ β ω µ ν νω µ ω µ ν µ τ = + + + Because , , vih S V S ∈ Γ we have:  parameter. Namely, if 0 1 R < ; all solutions in the feasible region converge to the disease-free equilibrium 0 P ; and the TB will die out from the population irrespective of the initial conditions. If 0 1 R > ; 0 P is unstable, it could exist and endemic equilibrium and the system is uniformly persistent, TB epidemic will always become endemic.

Optimal Control of System 1
In this section, an optimal control is formulated and it examined to study properties of optimal control strategies.
Let us pose Hamiltonian of our control problem is The adjoint equations and transversality conditions can be obtained by using From relations 14 we have: From relations 16 we have: Taking account the optimality conditions 14 induce by the bounds conditions of control u we have

The Numerical Algorithm
The numerical algorithm presented below is a classical Rung -Kutta 4 method. We discretize the interval [t0,T = tf] at the points 0 ti t ih = Next, we define the state and adjoint variables respectivily ( ) The control u in terms of nodal points ( )  Table 2 for the description of parameters and their based line or range value.

Optimal Strategies of Mass Vaccination
Some key parameters like the effective contact rate β medical coverage rate of VIH/Patients σ , HIV vertical transmission rate 3 π and HIV prevalence ε have a great impact in the spread of TB infectious. we are going to simulate five different scenarios and observed the mass vaccination optimal strategy induced respectively by the low effective contact rate, the hight effective contact rate, the absence of HIV medical coverage, the hight HIV medical coverage and the combination (hight HIV medical coverage, low HIV vertical transmission rate, and the low HIV prevalence rate) which assure the eradication of TB infectious. The first one is tuberculosesida.m the main file and the second one is val t uberculosesida.m in order to plot some figures.
In matlab code, the parameters are named as followed. pi1 = π 1 ; pi2 = π 2 pi3 = π 1 ; beta = β varepsilon = ε; muvih = μ vih ; sigma = σ; theta1 = θ 1 ; theta2 = θ 2 p; muS = μ S ; muV = μ V ; muE = μ E ; omega = ω; muL = μ L ; nu = ν; muI = μ I ; tau = τ A; alpha = α; S0 = S 0 ; V0 = S 0 ; Svih0 = S 0 ; E0 = E 0 ; L0 = L 0 ; I0 = I 0 . With the precedents values of parameters we obtain a Reproduction number ( 0 R ) bigger than unity. In this case Figure 3, the system will reach definitively and endemic equilibrium. Controlling the system is not necessary until the trajectory of TB patients approach the endemic equilibrium and then optimal control appears in step 110 and goes until step 200 in order to have a better endemic equilibrium. Here in Figure 4, the effective contact rate β is higher than the precedent case, optimal control is more aggressive and appears earlier in the step 70, and is stopped before step 200; the consequence is that the trajectory of vaccinated  reached more later abscise axis and TB patients leave abscise axis later also, this means that spread of epidemic is retarded. we could conclude that more earlier system is controlling more later epidemic occurs in the population. In this case 5 4 cd T is nul, that's suppose HIV patients are not cured. In Figure 5, we can see that optimal control appears quickly at years fourth, in the same time hen TB patients are leaving the abscise axis and epidemic occurs. We remark also that optimal control disappears after few months. We can conclude that:the lack of HIV medical coverage accelerate the propagation of TB infectious, and when the velocity of epidemic is hight optimal control appears quickly and disappears quickly also. Also mass vaccination is useless when we have lack of HIV medical coverage.

Optimal Control Induced by a VIH Medical Coverage:
When VIH patients are cured, optimal control appears later after 40 years and stay more than five years Figure 6 here Optimal control is better than the precedent case 5 more than 200 TB patients at equilibrium. In Figure 6, we have above 150 TB patients at equilibrium.   When HIV medical coverage is hight, control optimal appears more later step 80 and stay until step 200. In Figure 7, optimal control is more efficient because at endemic equilibrium we have 100 TB patients lesser than the precedent case where he have 150 TB patients in endemic equilibrium. Here we have reproduction number 0 R less than unity, and the system 1 reach a disease free equilibrium (Figure 8). The number of TB patients decreases and reach abscise axis at step 50. Nowadays optimal control appears at step 70 and the vaccination rate coverage decreases progressively until 0.5% at step 100.   ) assures the eradication of TB infectious in community. We can conclude that, well controlling the spread of VIH infectious (hight medical coverage of HIV patients, voluntary HIV screening campaign and HIV awareness campaigns) has positive effects in the propagation of TB infectious. Mass vaccination is not necessary when we have at least a good percentage of peoples who are vaccinated at birth.

Eradication of Desease with Hight VIH
All the scenarios are resumed in the following Figure 9.

Discussions and Conclusion
The goals of this paper were to study the overall and asymptotic stability of the system around the point of equilibrium on one hand and, to use optimal control techniques to find mass vaccination strategies for each situation and assess impact on the second hand. We simulated the spread of tuberculosis/HIV coinfection and mass vaccination schedule. The database used was essentially made of the World Health organisation report 2017 and the data collected at the Centre Jamot Hospital where is housed the screening center of tuberculosis. The Open Journal of Modelling and Simulation approach of optimal control used in various cases with adapted tools, leads to some lessons learnt. The higher the effective contact rate will be, the earlier we should start the mass vaccination that could economically expensive ( Figure 3 and Figure 4). The better the medical coverage will be, the the later we should engage mass vaccination as we reach the equilibriums (Figures 5-7). We eradicate the TB propagation by ensuring a better medical coverage, by reducing the HIV transmission from the mother to the child and by reducing the prevalence of HIV within the population (Figure 8).
In the area of HIV and TB co-infection, to reach the target of eradication of the TB propagation, we need to control the HIV propagation and make an emphasis on the immunization against the TB infection, the early screening and treatment of TB patients.