Computational and Mathematical Methods in Medicine
Volume 2018 (2018), Article ID 2657461, 21 pages
https://doi.org/10.1155/2018/2657461

Research Article

Mathematical Modelling of Bacterial Meningitis Transmission Dynamics with Control Measures

Joshua Kiddy K. Asamoah,1,2 Farai Nyabadza,³ Baba Seidu,⁴ Mehar Chand,⁵ and Hemen Dutta⁶

¹Department of Mathematics, Kwame Nkrumah University of Science and Technology, Kumasi, Ghana
²African Institute for Mathematical Sciences, Biriwa, Ghana
³Division of Mathematics, Stellenbosch University, Western Cape, South Africa
⁴Department of Mathematics, University for Development Studies, Navrongo, Ghana
⁵Department of Applied Sciences, Guru Kashi University, Bathinda, India
⁶Department of Mathematics, Gauhati University, Guwahati 781014, India

Correspondence should be addressed to Joshua Kiddy K. Asamoah; jasamoah@aims.edu.gh

Received 23 September 2017; Accepted 4 January 2018; Published 27 March 2018

Academic Editor: Konstantin Blyuss

Copyright © 2018 Joshua Kiddy K. Asamoah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Vaccination and treatment are the most effective ways of controlling the transmission of most infectious diseases. While vaccination helps susceptible individuals to build either a long-term immunity or short-term immunity, treatment reduces the number of disease-induced deaths and the number of infectious individuals in a community/nation. In this paper, a nonlinear deterministic model with time-dependent controls has been proposed to describe the dynamics of bacterial meningitis in a population. The model is shown to exhibit a unique globally asymptotically stable disease-free equilibrium , when the effective reproduction number , and a globally asymptotically stable endemic equilibrium , when ; and it exhibits a transcritical bifurcation at . Carriers have been shown (by Tornado plot) to have a higher chance of spreading the infection than those with clinical symptoms who will sometimes be bound to bed during the acute phase of the infection. In order to find the best strategy for minimizing the number of carriers and ill individuals and the cost of control implementation, an optimal control problem is set up by defining a Lagrangian function to be minimized subject to the proposed model. Numerical simulation of the optimal problem demonstrates that the best strategy to control bacterial meningitis is to combine vaccination with other interventions (such as treatment and public health education). Additionally, this research suggests that stakeholders should press hard for the production of existing/new vaccines and antibiotics and their disbursement to areas that are most affected by bacterial meningitis, especially Sub-Saharan Africa; furthermore, individuals who live in communities where the environment is relatively warm (hot/moisture) are advised to go for vaccination against bacterial meningitis.

1. Introduction

Meningitis is an inflammation of the meninges which are membranes that surround the spinal cord and the brain [1]. It is often caused by viruses, bacteria, and protozoa. Bacterial meningitis is common in children and young adults. This disease mostly spreads in communities/societies that live in close quarters (e.g., police staff, police cells, college students, military staff, and prisons) [2]. Bacterial meningitis is generally caused by germs such as Listeria monocytogenes, Streptococcus pneumoniae, Group B Streptococcus, Neisseria meningitidis, and Haemophilus influenzae, which spreads from one person to another [3]. This infection varies by age groups: Group B Streptococcus, Streptococcus pneumoniae, Listeria monocytogenes, and Escherichia coli are mostly found in newborn babies; Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae type b (Hib), and Group B Streptococcus are common in babies and children; Neisseria meningitidis and Streptococcus pneumoniae are predominant in teens and young adults; and Streptococcus pneumoniae, Neisseria meningitidis, Haemophilus influenzae type b (Hib), Group B Streptococcus, and Listeria monocytogenes are commonly found in older adults [3]. Bacterial meningitis is characterized by intense headache and fever, vomiting, sensitivity to light, and stiff neck, which result in convulsion, delirium, and death.

It is estimated that meningococcal meningitis causes over 10,000 deaths annually in Sub-Saharan Africa [4]. About 4,100 cases of bacterial meningitis occurred between 2003 and 2007 in the United States [3, 5]. Between 5% to 40% of children and 20% to 50% of adults with this condition die [6]. Infections from bacterial meningitis can cause permanent disabilities such as brain damage, hearing loss, and learning disabilities [3]. The illness of bacterial meningitis becomes worse when symptoms are not detected early enough; even with proper treatment, the individual could die [2].

Prevention of bacterial meningitis can be achieved through vaccination and/or preventing contact with infectious individuals. Vaccination is the most effective way of protecting children against certain types of bacterial meningitis [3]. Vaccines that can prevent meningitis include Haemophilus influenza type B (Hib), pneumococcal conjugate, and meningococcal vaccine [6, 7]. The conjugate meningitis A vaccine, MenAfrivac, is recommended to protect people in Sub-Saharan Africa against the most common type, serotype A [8]. In the United States, the primary means of preventing meningococcal meningitis is antimicrobial chemoprophylaxis [9]. Empirical therapy includes ceftriaxone or cefotaxime and vancomycin for Streptococcus pneumoniae [2]. There is a vaccine against meningococcal disease which is 85%–100% effective in preventing four kinds of bacteria (serogroups A, C, Y, and W-135) that cause about 70% of the disease in the United States [2].

Trotter and Ramsay [10] outlined some recommendations on the use of conjugate vaccines in Europe based on the experience with meningococcal C conjugate (MCC) vaccines. In areas with limited health infrastructure and resources, there are a number of antibiotics including penicillin, ampicillin, and chloramphenicol that can be used to treat the infection meningitis.

Mathematical models have been shown to help increase the understanding of the spread and control of infectious diseases. Martínez et al. [2] studied the spread of meningococcal meningitis with the use of a discrete mathematical model, based on cellular automata where the population was divided into five classes: susceptible, asymptomatic infected, infected with symptoms, carriers, recovered, and died classes. Broutin et al. [11] studied the dynamics of meningococcal meningitis in nine African countries by adopting some mathematical tools to time series analysis and wavelet method, the results of their studies suggest that “international cooperation in Public Health and cross disciplines studies are highly recommended to help in controlling this infectious disease.” Miller and Shahab [12] studied the cost effectiveness of immunisation strategies for the control of epidemic meningococcal meningitis. The research work in [13] gives a detailed description of the use of antibiotics for the prevention and treatment of meningitis infection. Irving et al. [14] used deterministic compartmental models to investigate how well simple model structures with seasonal forcing were able to qualitatively capture the patterns of meningitis infection. They demonstrated that the complex and irregular timing of epidemics could be caused by the interaction of temporary immunity conferred by carriage of the bacteria together with seasonal changes in the transmissibility of infection. Actually, there have been a significant number of studies of various types of Meningitis in Africa and Europe without the use of optimal control analysis (see [15–28]).

It is obvious that mathematical modelling has become crucial in investigating the epidemiological behaviour of meningitis. Furthermore, mathematical modelling helps to identify the risk factors for diseases, so as to find out why everyone does not have the same infection uniformly [29].

The application of optimal control in disease modelling gives valuable information on how to apply control measures. Through vaccination, treatment, public education, and so forth, many infectious diseases have been controlled [29]. Since the introduction of optimal control theory in disease modelling, there have been a considerable number of studies of infectious diseases using optimal control analysis (see [30–41]). With the significant influence of optimal control theory in disease modelling, this paper presents an optimal control model for bacterial meningitis in the presence of vaccination and treatment due to public health education. The model is qualitatively analyzed and numerically simulated in order to help give policy direction on how to control the spread of the disease.

The rest of the paper is organized as follows. Section 2 presents the model formulation and analysis. Section 3 presents the analysis of the optimal control problem, leading to the existence and characterization of the control measures. Section 4 contains the numerical simulations and discussion. Section 5 presents the conclusion of the study.

2. Model Formulation and Analysis

Adopting the epidemiological studies of a meningitis model as presented in [4], we consider four mutually exclusive compartments to indicate individuals with unique natures (i.e., susceptibles, , carriers, , ill individuals, , and recovered individuals, ) in relation to the disease. It is assumed that the susceptible compartment, , is populated through recruitment at the rate, (thus migration and/or birth rate), and is the rate of effective contact of carriers and/or infected (ill) individuals in the susceptible population. The carrier compartment consists of individuals that have the infection and do not show any clinical symptoms but contribute to the spread of the disease. When a susceptible individual is exposed to this infection, that individual can harbour the bacterium for weeks or even months [42]; but, in a normal circumstance, an individual develops symptoms of the infection within 3 to 7 days after exposure [3]. Carriers are assumed to develop clinical symptoms (i.e., move to ill individuals compartment, ) at rate . Ill individuals who are seriously infected are assumed to have no natural recovery except when given treatment on time. From an epidemiological perspective, individuals in the removed/recovered compartment, , do not attain permanent immunity. After vaccination, immunity develops within 7–10 days and remains effective for approximately 3–5 years [2]. Therefore, it is assumed that the immunity acquired from developing the diseases or carrying the bacteria or through vaccination is of the same intensity (they all lead to the recovered compartment from which people return to the susceptible compartment at a given unique rate ). Research indicates that carriers may recover naturally from the infection without treatment, and we denote such natural recovery rate as . Infected individuals are assumed to die from disease at rate . Since natural death is inevitable, is assumed to be the natural death rate of individuals in all the compartments. Vaccination and treatment due to public health education have been shown to be strategies of control of diseases. Therefore we introduced this two control measures in the model as and , respectively, here ; thus vaccination is comprised of both reactive vaccination and preventative vaccination. The effectiveness of both control measures in minimizing the disease is denoted by and , respectively. If , it signifies that vaccination and treatment have no effect on the model; if , it also signifies that vaccination and treatment are perfectly effective (i.e., 100% effectiveness) (see [43]). If , it signifies that both vaccine and treatment are imperfect [44]. In view of this, it is assumed that administering treatment to the ill individuals leads to recovery at rate . It is also assumed that the vaccinated individuals develop partial immunity at rate . From an epidemiological perspective treatment is not given to carriers in real life (since we do not know who carries the bacteria or not), but it is assumed that carriers of meningitis are just like people with the HIV virus who do not know their status unless they go for medical test; hence this paper seeks to encourage individuals to go for regular test of this bacterial disease; therefore it is assumed that a certain portion of carriers could be treated before any symptoms of the infection show up, which results in as shown in Figure 1. The vaccine is assumed to be imperfect and thus has a failure rate of . Therefore the force of new infections is given by

Equation (1) is often referred to as standard incidence rate of new infections, which is normalized by the total population . Table 1 gives a full description of parameters used in the model.

Table 1: Description of parameter values used in the model. “Due to lack of relevant data, most of the parameter values are assumed within realistic ranges for a typical scenario in a rural community/society for the purpose of illustration” [44], with most of the parameter values where appreciable per year.

Parameter
Description
Value/Range
Reference

Recruitment rate (migration and/or birth rate)
100–100000
Assumed

Effective contact rate

[28]

Loss of immunity
0.04–2
[14, 47]

Natural death rate

[14]

Disease-induced mortality
0.05–0.5
[48]

Rate of progression from to
0.1–0.52
[14]

Natural recovery rate
0.06–0.2
Assumed

Per capita infection rate by ill individuals
0.2–0.95
Assumed

Per capita infection rate by Carriers
0.2–0.85
Assumed

Vaccination rate of susceptible individuals against Meningitis

Assumed

Treatment rate for carriers without natural recovery and ill individuals

Assumed

Effectiveness of vaccination
0.85–1
[2]

Effectiveness of Treatment
0.1–0.9
Assumed

Figure 1: The flowchart diagram describing bacterial meningitis transmission dynamics within the population. The four circles represent the four compartments of individuals, the movement between the compartments is indicated by the continuous arrows, is a control measure (vaccination), and is the second control measure (treatment), with the consideration that both control measures lies in .

The set of differential equations and flow diagram corresponding to the bacterial meningitis dynamics and disease pathway with control terms is given in system (2) and Figure 1.

Adding the equations in model (2) gives the rate of change of total population as

Let

System (2) is well-posed with all solutions in remaining in if initial conditions are positive. It can easily be shown that if the initial conditions start outside the solutions tend to .

2.1. Equilibrium Points

To obtain the equilibrium points of system (2), we assume that the control measures are time-independent (see [45] for similar analysis).

2.1.1. Disease-Free Equilibrium

To obtain the disease-free equilibrium, ,  , and the right-hand-side of system (2) are set to zero. If susceptible individuals are assumed to receive vaccination against the disease at a constant rate, then the disease-free equilibrium will be given by

2.1.2. Effective Reproduction Number

Using the next generation matrix method [46], the effective reproduction number of the bacterial meningitis model with vaccination and treatment is obtained as

where

To determine how the two control measures impact on the reproduction number, we make a plot of , on the plane for arbitrary constant values of model parameters in Figure 2. From the figure, it is shown that decreases with increasing . So an increase in vaccination levels decreases the need for treatment. Figure 2 also shows the region in which the vaccination and treatment values should lie for the disease control.

Figure 2: Region where and in the parameter space, with the parameter values .

2.1.3. Endemic Equilibrium

System (2) can be shown to have a unique endemic equilibrium of the form , where

with and

Remark 1. (i) If , then system (2) will have only one equilibrium: the disease-free equilibrium.
(ii) If , then system (2) will have two equilibria: the disease-free equilibrium, , and the endemic equilibrium, .
(iii) The case is a critical threshold point where the disease-free equilibrium loses its local asymptotic stability. Thus gives the idea of transcritical bifurcation where the stability of system (2) moves between and [37].

2.2. Stability Analysis

In analyzing the local stability of the disease-free equilibrium, the Routh-Hurwitz criteria are used, and, for the global stability of the two equilibria, the direct Lyapunov technique is employed.

2.2.1. Local Stability Analysis of

Theorem 2. The disease-free equilibrium of system (2) is locally asymptotically stable if and unstable if .

Proof. Evaluating the Jacobian matrix of system (2) at the disease-free equilibrium gives

The characteristic polynomial of the Jacobian matrix is given by

where

The Routh-Hurwitz conditions [46] that guarantee that the eigenvalues of the characteristic polynomial in (10) have negative real parts are given by

These conditions are easily seen to be satisfied when . Thus, the disease-free equilibrium of system (2) is locally asymptotically stable when and unstable when . This completes the proof.

2.2.2. Global Stability of

Theorem 3. The disease-free equilibrium of system (2) is globally asymptotically stable if and unstable if .

Proof. Let , with positive constants, and , be a Lyapunov function defined as

Taking the time derivative of the Lyapunov function we obtain

Substituting ,  ,  , and in (2) into (14) gives

Since

this implies that

Equating the coefficient of in (17) to zero gives

Choosing   and   and plugging and into (17), we have

Additionally if and only if . Hence, the largest compact invariant set in is the singleton set . Therefore, from LaSalle’s invariance principle, we conclude that is globally asymptotically stable in if [37, 49].

2.2.3. Global Stability of

Theorem 4. The endemic equilibrium of system (2) is globally asymptotically stable whenever .

Proof. Suppose , and then the existence of the endemic equilibrium point is assured. Using the common quadratic Lyapunov function

as illustrated in [50], we consider the following candidate Lyapunov function:

The time derivative of in (21) is given by

Plugging the equations in system (2) into (22) yields

Now setting

we have

Further simplification gives

It has therefore been shown that is negative, and additionally at (i.e., if , , , and ), . It follows from LaSalle’s invariant principle [51] that all solutions of system (2) approach as if . Therefore, the endemic equilibrium is globally asymptotically stable in whenever [37, 49]. This completes the proof.

2.3. Sensitivity Analysis

Sensitivity analysis is used to determine the response of a model to variations in its parameter values. In the present case, the focus is given to determining how changes in the model parameters impact the effective reproduction number. This is done through the normalized forward-sensitivity index. We also use the Latin hypercube sampling and the partial rank correlation coefficients (PRCC) to plot scatter diagrams and Tornado plots to determine the relative importance of the parameters in for the disease transmission and prevalence (see also [52]).

Definition 5. The normalized forward-sensitivity index of to any parameter, say , as given in [46] can be defined as

The sensitivity indexes of with respect to its parameters are computed as follows:

Similarly, we can compute the sensitivity indexes of with respect to the remaining parameters in , in the same manner. Using the parameter values , , , , , , , , and , with and , the sensitivity indexes of are shown in Table 2.

Table 2: Sensitivity indexes of to the parameters in (6).

Parameter
Description
Sensitivity index

Effective contact rate
1

Per capita infection rate by Carriers
0.6154

Per capita infection rate by ill individuals
0.3846

Loss of immunity
0.6686

Rate of falling ill
0.0620

Natural death rate
0.1031

Disease-induced death rate
−0.1202

Natural recovery rate
−0.3226

Vaccination rate
−1.8280

Effectiveness of vaccination
−0.8280

Treatment rate
−0.5630

Effectiveness of treatment
−0.5230

The corresponding Tornado plots based on a random sample of points for the twelve parameters in are shown in Figure 3(a). The positive values in Table 2 show a promotion of the propagation of the disease. Therefore an increase in the values of , , , , and will have an increase in the spread of the disease. For example, indicates that increasing the effective contact rate by 10% increases the number of secondary infections by 10%. The negative values in Table 2 indicate a reduction in the effective reproduction number if the values of the corresponding parameters are increased. Thus, a reduction in the values of vaccination , treatment , and natural recovery will lead to an increase in the number of secondary infections in the population.

(a) The Tornado plots for the elven parameters in

(b) The effect of effective contact rate on

(c) The effect of vaccination on

(d) The effect of treatment on

Figure 3: Sensitivity plots.

Figure 3(a) shows the Tornado plots for the twelve parameters in . It can be seen that, in controlling the spread of bacterial meningitis in a population, more susceptible individuals should be given vaccination. Figure 3(a) also suggests that carriers are likely to have more contacts with the susceptible population than the ill individuals who will typically be bound to their beds during the acute phase of the disease. Therefore, the probability of ill individuals transmitting the infections to susceptibles may be lower than that of carriers who are able to mix well with others within the population. Figures 3(b), 3(c), and 3(d) show the regression plots of effective contact rate (), vaccination rate (), and treatment rate (, respectively. Figure 3(b) shows that transmission rate has a positive correlation in the spread of bacterial meningitis. Figure 3(c) shows that vaccination has a negative correlation in the spread of bacterial meningitis and hence vaccination increases the immunity of individuals against the meningitis infection, thereby reducing the spread of the infection. Figure 3(d) shows that individuals who receive treatment after being infected with bacterial meningitis have a higher chance of recovery and that reduces the spread of the infection and deaths due to bacterial meningitis.

3. Optimal Control Problem

Since the goal of this paper is to find the best ways to control the spread of meningitis, we define the following optimal control problem:

subject to model (2).

The admissible control set is Lebesgue measurable, which is defined by

Our objective is to find which minimizes the associated cost of the vaccination and the associated cost of the treatment over the specified time interval, as well as minimizing the number of infections at a terminal time (see also [37]). The coefficients and are constants that are introduced to maintain a balance in the size of and , respectively. and are the corresponding weights associated with the cost of vaccination and treatment , respectively. The higher bounds (maximum) attainable for the control measures and are and , respectively. We fix the control measures and to lie between and so that and . Therefore the attainment of and depends on the number of resources available [37]. These resources may include the human effort, material resources, cost of producing vaccine and disbursement, infrastructural resources, the number of health facilities in the community, and the number of hospital beds at the health facilities. The cost of hospitalization, medical test, diagnosis, drug cost, and so forth (see [38–40]) can be associated with treatment. The cost of vaccination may include the cost of the vaccine, the cost of production, the cost of disbursement, the vaccine storage cost, and other related overheads [37]. The severity of the side effects and overdoses of the vaccination and treatment is taken care of by squaring the control measures, and is the final time during the optimal simulation.

3.1. Existence of the Optimal Control

Model (2) can be written as

where

and is the times derivative of . System (31) is nonlinear with a bounded coefficient.

Setting and gives

Therefore

where

so that

The function is therefore uniformly Lipschitz continuous. From the definition of the control measures and and the constraint on the state variables, such that , , , and , we observe that a solution of system (31) exists [40, 53, 54]. From the objective functional and its associated constraints in model (2), we can find the optimal solution for our model. Firstly, we find the Lagrangian () and Hamiltonian () for the control problem [55]. The Lagrangian of the optimal problem is given by

Our focus is to find the minimal value of the Lagrangian function, which is done by a pointwise minimization of the Hamiltonian () defined as follows (using Pontryagin’s maximum principle):

where , are the adjoint variables associated with , , , and , defined by

Theorem 6. There exists an optimal control pair such that

Proof. We start our proof by considering the properties of the existence of the optimal control (see [56]). Following [57], the set of control measures with corresponding state variables are positive. The set is convex and closed by definition. Therefore, our optimal system is closed and bounded which ascertains the compactness required for the existence of the optimal control. Additionally, the integrand in the objective functional (29), , is convex on the control set . Furthermore, we can state that there exists a positive constant [58], and nonnegative numbers and such that the objective functional has a lower bound of so that

since the control measures and the state variables are bounded, this leads us to a compact proof of existence of the optimal control.

3.2. Characterization of the Optimal Control

We will apply Pontryagin’s maximum principle to the Hamiltonian function above to derive the necessary condition of optimality for our control problem.

Theorem 7. Let , , , and be optimal state solutions with corresponding optimal control variables and for the objective functional and its constraints in model (2) with . Then, there exist four adjoint variables , , , and that satisfy

with transversality conditions

Therefore, the optimal control pair is given by

Proof. We use the Hamiltonian function in (38) in order to obtain the adjoint relations and the transversality conditions. We set the state variables in the Hamiltonian function to , and , and differentiating the Hamiltonian with respect to , and , respectively, yields (42). Also, differentiating the Hamiltonian with respect to the control measures and in the interior of , we obtain the optimality conditions below:

Plugging and into (45) and solving the optimal control pair , we have

The two control measures which are bounded with lower bounds zero and upper bounds , where , give

Using (47), the optimal control measures are characterized as (44), completing the proof.

Therefore, our optimality system is given by

4. Numerical Simulations of the Optimal Control

The optimality system consisting of the state equations (2) and the adjoint equations (45) is solved using the forward-backward sweep scheme. The solution is started with an initial guess for the control measures and the final time set to days and later varied to days and days. The state system is solved forward in time while the costate system is solved backward in time. The current solutions of the state system together with the initial guess for the control measures and are used to solve the costate system. The controls and are updated using the characterizations in (47). With the current cosystem solution and updated controls, the state system is solved and the whole solution process continues until convergence is achieved. Since bacterial meningitis is an endemic disease especially in Sub-Sahara Africa, parameter values that make effective reproduction number more than unity are considered. Taking the associated costs on carriers and ill individuals as , , and , the graphical results of the model are as follows. Figures 4 and 5, show the dynamics of the infection in the presences of control measures and without control measures in a more localised form. Figure 5(f) shows the effect of varying the cost associated with the treatment control. It is seen that if the weight of treatment cost is continuously increased, the upper bound time reduces to about days as indicated by the red line in Figure 5(f). This is true because a higher cost of treatment causes the use of the treatment control to be less. Figure 6 shows the spreading rate of meningitis with and without control measures in a highly dense population, thus in colleges, prisons, cities, and so forth, which suggests that when carriers become more in a highly populated settings there is a likelihood of the infection spreading faster, thereby leading to a corresponding increase in the number of ill (sick) individuals; hence the trajectories in Figures 6(a) and 6(b) show that there is a significant difference in the number of carriers and infected individuals with and without control measures. Therefore, applying both vaccination and treatment has a higher rate of controlling bacterial meningitis than depending only on one control method. Figures 6(a) and 6(b) also show that as infections increase in the population, there is a need for introducing vaccination and treatment at the earlier stage so as to minimize the spread of the infection in the population. This further suggests that college students and prisoners should be given regular vaccination against meningitis, since bacterial meningitis is more likely to affect college students and prisoners than other people, due to the close proximate of beds in the dormitories and prison cells, especially in Sub-Saharan Africa. Figure 6(c) shows that the trajectory for treatment moves to an upper bound at about days and slowly decreases to a lower bound, which implies that, during the period before days, substantial amount of the control should be applied while administering the vaccination control in the susceptible population so as to reduce the number of carriers and ill individuals until the 50th day. Afterwards, less amount of treatment can be used since the number of carriers and ill individuals will be considerably reduced by the earlier investment in treatment within the period before days.

(a) The optimal effect on with vaccination control

(b) The optimal effect on with vaccination control

(c) The optimal effect on with an increase in time

(d) The optimal effect on with an increase in time

(e) Vaccination control profile

(f) Optimal control profile with an increase in time

Figure 4: (a), (b), (c), and (d) show that administering vaccination on susceptible individuals reduces the number of secondary infections in the carrier and the ill population. (e) and (f) show the optimal control profile of vaccination. Initial conditions , and and the parameter values were used in this simulation.

(a) The optimal effect on with treatment control

(b) The optimal effect on with treatment control

(c) The optimal effect on with an increase in time

(d) The optimal effect on with an increase in time

(e) Treatment control effect

(f) Optimal control plot of treatment with different cost weights

Figure 5: (a) and (b) show that giving treatment to carriers and the ill individuals reduces the number of secondary infections in the carrier and the ill population but the disease may not be eradicated if the use of only treatment is stuck to. Furthermore, we increased time so as to see whether the disease could be eradicated when using treatment only. We see in (c) and (d) that there is a possibility of recurrence of the disease using only treatment control strategy as indicated by the curved green line in (c) and the linear flow of the green line in (d). (e) shows that, to minimize bacterial meningitis outbreak in days, the treatment control should be held intensively for days at a constant rate. (f) shows the optimal control plot of treatment with different cost weights.

(a) The optimal effect on with both control measures

(b) The optimal effect on with both control measures

(c) Simulation diagram for control measures

(d) Immunity of recovered individuals

Figure 6: Further simulation trajectories. (d) shows the immunity levels when the two control measures are applied with the following parameter values: ,  , and .

Figures 7(a), 7(b), 8(c), and 8(d) show the impact of varying and . It can be seen that an increase in and without control measures has a corresponding decrease in the number of carriers and a considerable increase in the number of infected individuals. Therefore, applying the control measures gives a drastic decrease in the number of infections in both carriers and the ill compartment at a higher rate of and . Furthermore, in Figures 7(a) and 7(b), authors kept the natural recovery rate () fixed at 0.2 while assessing the effect of varying the rate of falling ill on carriers and on ill individuals, respectively. Hence, considering an increased rate from 0.2 to 0.7 of the rate of falling ill as shown in Figures 8(a) and 8(b) demonstrates how serious the infection is in the absence of vaccination or treatment, which indicates that ill individuals have a lower means of recovering from the disease, not even through a higher rate of natural defences of as indicated in Figure 8(b) without vaccination and treatment, which clearly shows that vaccination and treatment are essential in controlling meningitis irrespective of the age group; therefore we suggest that health authorities should increase the rate at which susceptible individuals get vaccination through public education and also encourage individuals with and without symptoms of the infection to visit health centers for a quick check-up and if the infection is detected an immediate treatment/vaccination should be given on time to avert the spreading of the infection to community members, relatives, or close friends.

(a) Optimal simulation effect in varying on

(b) Optimal simulation effect in varying on

Figure 7: Simulation results of the optimal control model with α = 0.2-0.5, on carriers and the ill compartment.

(a) Numerical simulation without control measures on

(b) Numerical simulation without control measures on

(c) Optimal simulation effect in varying on

(d) Optimal simulation effect in varying on

(e) Effectiveness of treatment on

(f) Effectiveness of treatment on

Figure 8: Simulation results of the optimal control model with η = 0.35-0.65. The solutions for the carriers and ill individuals with and give and , respectively. The case where and gives and , respectively, which indicates the existence of an unstable disease-free equilibrium and a stable endemic equilibrium. In (e) and (f), vaccination and treatment rates are set to and the effectiveness of the treatment control varied. (e) and (f) show that as the effectiveness of the treatment control increases, the number of ill individuals reduces faster as compared to the number of carriers.

4.1. Vaccination Control Only

Figures 4(a)–4(f) show simulation results of implementing only control strategy with set to zero.

4.2. Treatment Control Only

Considering treatment control, , by setting the vaccination control, , to zero, the simulation results are as shown in Figures 5(a)–5(f).

4.3. Applying Both Control Measures

Finally, both control measures are implemented using the following initial population size of , , , and and the results are plotted in Figures 6(a)–6(d).

Figures 9(a), 9(b), 10(a), and 10(b) show the interepidemic nature of some selected parameters in the model presented above. Figures 9(a) and 9(b) indicate that as individuals lose immunity and become susceptible again, there is a corresponding increase in the rate of infection, and this suggests that an individual who receives vaccination and/or treatment and remains immune for 3–5 years should immediately go for another vaccination against the infection so as to avert transmission of the infection to that individual. Figures 10(a) and 10(b) depict that the per capita infection rate by ill individuals, , and the per capita infection rate by carriers, , are interlinked and so have a corresponding influence in the spread of the infection. Figure 10(b) also indicates that as the infectivity rate of carriers and ill individuals goes beyond , there could be a high spread of the infection in the community/society should there be an outbreak, since the effective reproduction number will be getting closer to one with a possibility of going beyond a unit.

(a) The contour plot of and to

(b) The 3D plot of to and

Figure 9: The interepidemic period of the SCIRS model depending on parameters and .

(a) The contour plot of and to

(b) The 3D plot of to and

Figure 10: The interepidemic period of the SCIRS model depending on parameters and .

5. Conclusion

We presented a mathematical framework of vaccination and treatment on bacterial meningitis model. The investigation of the stability of the model shows that the disease-free equilibrium is locally and globally asymptotically stable and the endemic equilibrium depicts a global stability. Scatter plots and the Tornado plots of the twelve parameters in show that the effective contact rate has a major impact in transmitting the disease, followed by the infectivity potential of carriers . This supports the fact that asymptomatic carriers are likely to have more contacts in a community/nation when there is an outbreak of bacterial meningitis compared to ill individuals who will typically be bound to their beds during the acute phase of the disease, thereby lowering the probability of infecting susceptible individuals. Finally, the numerical simulation shows that the optimal (best) way of controlling the transmission of meningitis in Sub-Saharan Africa and the world at large is to encourage susceptible individuals to go for vaccination against meningitis at the health centers and also report any suspected symptoms of meningitis to health practitioners for early detection and immediate care (treatment).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

W. P. Howlett, Neurology in Africa: Clinical Skills and Neurological Disorders, Cambridge University Press, 2015.

M. J. F. Martínez, E. G. Merino, E. G. Sánchez, J. E. G. Sánchez, A. M. Del Rey, and G. R. Sánchez, “A mathematical model to study the meningococcal meningitis,” in Proceedings of the 13th Annual International Conference on Computational Science, ICCS 2013, pp. 2492–2495, Spain, June 2013.

CDC Centers for Disease Control, Bacterial Meningitis, 2017, http://www.who.int/gho/epidemic_diseases/meningitis/suspected_cases_deaths_text/en/.

K. B. Blyuss, “Mathematical modelling of the dynamics of meningococcal meningitis in Africa,” UK Success Stories in Industrial Mathematics, pp. 221–226, 2016.

M. C. Thigpen, C. G. Whitney, N. E. Messonnier et al., “Bacterial meningitis in the United States, 1998–2007,” The New England Journal of Medicine, vol. 364, no. 21, pp. 2016–2025, 2011.

CDC Centers for Disease Control, Non-polio enterovirus Infections, 2010, http://www.cdc.gov/ncidod/dvrd/revb/.

B. Shmaefsky and H. Babcock, MMeningitis, Infobase Publishing, 2010.

Worldbank, Meningitis vaccine provides hope to people in Ghana, 2012, http://www.who.int/features/2012/meningitis_ghana/en/.

CDC Centers for Disease Control, Control and Prevention of Meningococcal Disease: Recommendations of the Advisory Committee on Immunization Practices (ACIP), 2017, https://www.cdc.gov/mmwr/preview/mmwrhtml/00046263.htm.

C. L. Trotter and M. E. Ramsay, “Vaccination against meningococcal disease in Europe: review and recommendations for the use of conjugate vaccines,” FEMS Microbiology Reviews, vol. 31, no. 1, pp. 101–107, 2007.

H. Broutin, S. Philippon, G. Constantin de Magny, M.-F. Courel, B. Sultan, and J.-F. Guégan, “Comparative study of meningitis dynamics across nine African countries: A global perspective,” International Journal of Health Geographics, vol. 6, article no. 29, 2007.

M. A. Miller and C. K. Shahab, “Review of the cost effectiveness of immunisation strategies for the control of epidemic meningococcal meningitis,” PharmacoEconomics, vol. 23, no. 4, pp. 333–343, 2005.

National Collaborating Centre for Women’s and Children’s Health, Antibiotics for Early-Onset Neonatal Infection: Antibiotics for The Prevention And Treatment of Early-Onset Neonatal Infection, RCOG Press, 2012.

T. J. Irving, K. B. Blyuss, C. Colijn, and C. L. Trotter, “Modelling meningococcal meningitis in the African meningitis belt,” Epidemiology and Infection, vol. 140, no. 5, pp. 897–905, 2012.

B. Greenwood, “Priorities for research on meningococcal disease and the impact of serogroup A vaccination in the African meningitis belt,” Vaccine, vol. 31, no. 11, pp. 1453–1457, 2013.

A. C. Cohn and L. H. Harrison, “Meningococcal vaccines: Current issues and future strategies,” Drugs, vol. 73, no. 11, pp. 1147–1155, 2013.

R. Yogev and T. Tan, “Meningococcal disease: the advances and challenges of meningococcal disease prevention.,” Human Vaccines & Immunotherapeutics, vol. 7, no. 8, pp. 828–837, 2011.

S. B. Gordon, S. Kanyanda, A. L. Walsh et al., “Poor potential coverage for 7-valent pneumococcal conjugate vaccine, Malawi,” Emerging Infectious Diseases, vol. 9, no. 6, pp. 747–749, 2003.

A. J. Pollard, J. Ochnio, M. Ho, M. Callaghan, M. Bigham, and S. Dobson, “Disease Susceptibility to ST11 complex meningococci bearing serogroup C or W135 polysaccharide capsules, North America,” Emerging Infectious Diseases, vol. 10, no. 10, pp. 1812–1815, 2004.

W. P. Hanage, “Serotype-specific problems associated with pneumococcal conjugate vaccination,” Future Microbiology, vol. 3, no. 1, pp. 23–30, 2008.

G. B. Coulson, A. Von Gottberg, M. Du Plessis, A. M. Smith, L. De Gouveia, and K. P. Klugman, “Meningococcal disease in South Africa, 1999-2002,” Emerging Infectious Diseases, vol. 13, no. 2, pp. 273–281, 2007.

MC. Thomson, E. Firth, M. Jancloes et al., “A climate and health partnership to inform the prevention and control of meningoccocal meningitis in sub-Saharan Africa,” in the MERIT initiative. in. Climate Science for Serving Society, pp. 459–484, Springer, Climate Science for Serving Society, 2013.

J. E. Mueller, L. Sangaré, B.-M. Njanpop-Lafourcade et al., “Molecular characteristics and epidemiology of meningococcal carriage, Burkina Faso, 2003,” Emerging Infectious Diseases, vol. 13, no. 6, pp. 847–854, 2007.

V. Dukic, M. Hayden, A. A. Forgor, and et al., “The role of weather in meningitis outbreaks in Navrongo, Ghana: a generalized additive modeling approach,” Journal of Agricultural, Biological, and Environmental Statistics, vol. 17, no. 3, pp. 442–460, 2012.

A. Beresniak, E. Bertherat, W. Perea et al., “A Bayesian network approach to the study of historical epidemiological databases: Modelling meningitis outbreaks in the Niger,” Bulletin of the World Health Organization, vol. 90, no. 6, pp. 412–416, 2012.

S. Tartof, A. Cohn, F. Tarbangdo et al., “Identifying Optimal Vaccination Strategies for Serogroup A Neisseria meningitidis Conjugate Vaccine in the African Meningitis Belt,” PLoS ONE, vol. 8, no. 5, Article ID e63605, 2013.

L. Agier, H. Broutin, E. Bertherat et al., “Timely detection of bacterial meningitis epidemics at district level: A study in three countries of the African Meningitis Belt,” Transactions of the Royal Society of Tropical Medicine and Hygiene, vol. 107, no. 1, Article ID trs010, pp. 30–36, 2013.

K. Vereen, An SCIR Model of Meningococcal Meningitis, 2008.

H. W. Hethcote, “Three basic epidemiological models,” in Applied Mathematical Ecology, S. A. Levin, T. G. Hallam, and L. J. Gross, Eds., vol. 18 of Biomathematics, pp. 119–144, Springer, Berlin, Germany, 1989.

R. L. Miller Neilan, Optimal Control Applied to Population and Disease Models (Doctoral Dissertations), 2009.

A. O. Isere, J. E. Osemwenkhae, and D. Okuonghae, “Optimal control model for the outbreak of cholera in Nigeria,” African Journal of Mathematics and Computer Science Research, vol. 7, no. 2, pp. 24–30, 2014.

A. A. Lashari, Mathematical Modeling And Optimal Control Of A Vector Borne Disease, National University of Modern Languages, Islamabad, Pakistan, 2012.

A. Abdelrazec, S. Lenhart, and H. Zhu, “Transmission dynamics of West Nile virus in mosquitoes and corvids and non-corvids,” Journal of Mathematical Biology, vol. 68, no. 6, pp. 1553–1582, 2014.

A. Stashko, The effects of prevention and treatment interventions in a microeconomic model of HIV transmission, Duke University Durham, 2012.

B. Seidu and O. D. Makinde, “Optimal control of HIV/AIDS in the workplace in the presence of careless individuals,” Computational and Mathematical Methods in Medicine, vol. 2014, Article ID 831506, 2014.

S. D. Djiomba Njankou and F. Nyabadza, “An optimal control model for Ebola virus disease,” Journal of Biological Systems, vol. 24, no. 1, pp. 29–49, 2016.

T. T. Yusuf and F. Benyah, “Optimal control of vaccination and treatment for an SIR epidemiological model,” World Journal of Modelling and Simulation, vol. 8, no. 3, pp. 194–204, 2012.

H. Gaff and E. Schaefer, “Optimal control applied to vaccination and treatment strategies for various epidemiological models,” Mathematical Biosciences and Engineering, vol. 6, no. 3, pp. 469–492, 2009.

S. Nanda, H. Moore, and S. Lenhart, “Optimal control of treatment in a mathematical model of chronic myelogenous leukemia,” Mathematical Biosciences, vol. 210, no. 1, pp. 143–156, 2007.

G. Zaman, Y. Han Kang, and I. H. Jung, “Stability analysis and optimal vaccination of an SIR epidemic model,” BioSystems, vol. 93, no. 3, pp. 240–249, 2008.

J. K. K. Asamoah, F. T. Oduro, E. Bonyah, and B. Seidu, “Modelling of Rabies Transmission Dynamics Using Optimal Control Analysis,” Journal of Applied Mathematics, vol. 2017, Article ID 2451237, 2017.

U. Princeton, Emergency Guidelines for the Campus Community: Meningitis FAQ General Information, 2015, http://web.princeton.edu/sites/emergency/meningitis/FAQ-General.html.

C. Sun and Y.-H. Hsieh, “Global analysis of an SEIR model with varying population size and vaccination,” Applied Mathematical Modelling, vol. 34, no. 10, pp. 2685–2697, 2010.

S. M. Moghadas, “Modelling the effect of imperfect vaccines on disease epidemiology,” Discrete and Continuous Dynamical Systems - Series B, vol. 4, no. 4, pp. 999–1012, 2004.

D. Aldila, N. Rarasati, N. Nuraini, and E. Soewono, “Optimal control problem of treatment for obesity in a closed population,” International Journal of Mathematics and Mathematical Sciences, vol. 2014, Article ID 273037, 2014.

M. Martcheva, An Introduction to Mathematical Epidemiology, vol. 61 of Texts in Applied Mathematics, Springer, New York, NY, USA, 2015.

A. Karachaliou, A. J. K. Conlan, M.-P. Preziosi, and C. L. Trotter, “Modeling long-term vaccination strategies with MenAfriVac in the African Meningitis Belt,” Clinical Infectious Diseases, vol. 61, pp. S594–S600, 2015.

WHO World Health Organization, Global Health Observatory (GHO) data: Number of suspected meningitis cases and deaths reported, 2010, http://www.who.int/emergencies/diseases/en/.

S. D. Hove-Musekwa, F. Nyabadza, and H. Mambili-Mamboundou, “Modelling hospitalization, home-based care, and individual withdrawal for people living with HIV/AIDS in high prevalence settings,” Bulletin of Mathematical Biology, vol. 73, no. 12, pp. 2888–2915, 2011.

C. V. De Le and C. V. De León, “Constructions of Lyapunov functions for classics SIS, SIR and SIRS epidemic model with variable population size,” Foro-Red-Mat: Revista electrónica de contenido matemático, vol. 26, no. 5, p. 1, 2009.

J. LaSalle, “The stability of dynamical systems,” in Proceedings of the CBMS-NSF regional conference series in applied mathematics 25, Philadelphia, Pennsylvania, 1976.

S. Marino, I. B. Hogue, C. J. Ray, and D. E. Kirschner, “A methodology for performing global uncertainty and sensitivity analysis in systems biology,” Journal of Theoretical Biology, vol. 254, no. 1, pp. 178–196, 2008.

G. Birkhoff and G.-C. Rota, Ordinary Differential Equations, Ginn, 1982.

D. Iacoviello and N. Stasio, “Optimal control for SIRC epidemic outbreak,” Computer Methods and Programs in Biomedicine, vol. 110, no. 3, pp. 333–342, 2013.

H. Laarabi, M. Rachik, O. El Kahlaoui, and E. H. Labriji, “Optimal vaccination strategies of an SIR epidemic model with a saturated treatment,” Universal Journal of Applied Mathematics, vol. 1, no. 3, pp. 185–191, 2013.

S. C. Mpeshe, L. S. Luboobi, and Y. Nkansah-Gyekye, “Modeling the impact of climate change on the dynamics of rift valley fever,” Computational and Mathematical Methods in Medicine, Article ID 627586, Art. ID 627586, 12 pages, 2014.

D. L. Lukes, Differential Equations: Classical to Controlled, Mathematics in Science and Engineering, Academic Press, New York, NY. USA, 1982.

A. A. Lashari, K. Hattaf, G. Zaman, and X. Li, “Backward bifurcation and optimal control of a vector borne disease,” Applied Mathematics & Information Sciences, vol. 7, no. 1, pp. 301–309, 2013.