Abstract

Herpes simplex virus infection is mostly spread and occurs more commonly among substance abusing adolescents as compared to the nonsubstance abusing. In this paper, a mathematical model for the spread of HSV-2 within a community with substance abusing adolescents is developed and analysed. The impacts of condom use and educational campaigns are examined. The study suggests that condom use is highly effective among adolescents, when we have more of them quitting than becoming substance abusers. Measures such as educational campaigns can be put in place to try and reduce adolescents from becoming substance abusers. Further, we applied optimal control theory to the proposed model. The controls represent condom use and educational campaigns. The objective is based on maximising the susceptible nonsubstance abusing adolescents, while minimising the susceptible substance abusing adolescents, the infectious nonsubstance abusing adolescents, and the infectious substance abusing adolescents. We used Pontrygin’s maximum principle to characterise the optimal levels of the two controls. The resulting optimality system is solved numerically. Overall, the application of the optimal control theory suggests that more effort should be devoted to condom use as compared to educational campaigns.

1. Introduction

Herpes simplex virus type 2 (HSV-2), a double-stranded DNA virus, is a sexually transmitted disease that inflicts severe public health burden globally. The most common symptoms of HSV-2 usually appear within 1-2 weeks after sexual exposure to the virus. The first signs are a tingling sensation in the affected areas (genitalia, buttocks, and thighs) and group of small red bumps that develop into blisters. Herpes simplex virus-2 (HSV-2) is a highly prevalent STD that causes severe public health burden globally, with the highest prevalence in the Sub-Saharan Africa and some Asian countries. There are 1.5 million new herpes simplex virus-2 (HSV-2) infections among 15–19 year old females in Sub-Saharan African every year [1]. In 2003, for instance, up to 536 million people aged 15–49 years were living with HSV-2 globally and 23.6 million new infections were recorded in that year [2]. In the United States of America (USA), annual health costs for STIs has reached US $17 Billion [3], with HSV-2 chewing up $541 million, making it the third most costly STI after HIV-1 and human papillomavirus (HPV) [4].

Adolescents are the range of people who are 10–19 years old [5], making up approximately 20% of the world’s population and about 85% living in developing countries [6]. During the period of adolescence, many individuals begin to engage in risky sexual behaviours, hence making them the age group at greatest risk for nearly all STIs [7]. The reasons for this trend are many, including cognitive development, physiological susceptibility, peer pressure, logistic issues, and specific sexual behaviours [7]. Adolescents aged 15–19 years account for approximately 3 million cases, meaning one out of four sexually active teenagers reports an STI every year [8, 9]. Adolescents are more vulnerable to STIs because of lack of education, yet, in some schools, sex education in schools begins late when adolescents have already been initiated into sexual intercourse [10]. Biologically, adolescent girls are at higher risk of contracting STIs because, in adolescent girls, the cervix and vagina undergo dramatic histological changes due to estrogen exposure [11]. Adolescent girls are also the ones at greater risk of HSV-2 infections [12, 13]. Some studies have found out that being young at age really increased the likelihood of acquiring HSV-2 infection [1416], implying that we have more adolescents being infected. Prevalence estimates for adolescents are available through community-based studies. One study in Zimbabwe among adolescent females aged 15–19 years reported that 12% of participants tested positive for HSV-2 [17]. A study in western Kenya found HSV-2 prevalence to be 9% for females and 4% for males aged 13-14 years and 28% females and 17% for males aged 15–19 years [18]. It was shown that sexual activity in Namibia begins at a very young age, often as low as 10 years [10], and also some 61% percent of twelfth graders in United States of America (USA) have had sexual intercourse at least once [19]. Many adolescents experiment with substances, parents, teachers, and policy makers want to know the consequences of adolescent substance abuse. Since the abuse of substances affects their thinking capacity, most of them who are substance abusers have higher rates of HSV-2 infection [15, 20]. Those adolescents who are infected with HSV-2 type 2 (HSV-2) infection may be at greater risk for transmission and acquisition of human immunodeficiency virus (HIV) [21]. Interventions to control HSV-2 infection could prevent as many as 50% to 60% of new HIV infections [22]. Furthermore, pregnant adolescents infected with genital HSV (particularly those with a primary infection) can transmit infection to the neonate, which can lead to serious complications for the neonate, such as neurologic problems and even death [23, 24]. Since currently there is no vaccination for HSV-2 [25], current intervention options for controlling HSV-2 are limited to patient education and promotion of condom use or treatment with oral medication [26]. Some preventive efforts should be made to target children before initiation of sexual activity, since prevention is always better than cure.

A number of mathematical models have looked into mathematical modelling of HSV-2 transmission dynamics focusing on a number of different issues (see [3642], to mention just a few, but none has looked into the effects of substance abuse on the transmission dynamics of HSV-2 among adolescents). In contrast to other STIs such as chlamydia, gonorrhea, syphilis, and trichomoniasis, infection with HSV-2 is lifelong, and, once established, there is currently no treatment to eliminate [43, 44]. The other problem with HSV-2 infection in adolescents is that it is somewhat difficult to determine given that many infections are asymptomatic or go unrecognised; hence it differs with most of other STIs in general. To the best of our knowledge as authors, no mathematical model has looked into the transmission dynamics of HSV-2 and substance abuse with a focus on adolescents. It is against this background that our study finds relevance and motivation, by formulating a mathematical model to investigate the impact of substance abuse among adolescents on the transmission dynamics of HSV-2. The model incorporates some key epidemiological features such as the impact of condom use and educational campaigns among the adolescents. The main objective in this study is to forecast future trends in the incidence of HSV-2 epidemic and also to quantify the association between substance abuse and HSV-2 epidemic amongst adolescents in the community.

The paper is structured as follows. The HSV-2 transmission model is formulated in the next section. Analytic results of the model system are presented in Section 3. Simulation results and projection profiles of HSV-2 are presented in Section 4. In Section 5 optimal control theory has been applied to the model formulation in Section 2. Summary and concluding remarks round up the paper.

2. Model Formulation

The total sexually active population at time , denoted by , is subdivided into mutually exclusive compartments, namely, susceptible adolescents who are nonsubstance abusers , susceptible adolescents who are substance abusers , adolescents infected with HSV-2 who are nonsubstance abusers , and adolescents infected with HSV-2 who are substance abusers . Although there are many causes of substance abuse amongst adolescents here we dwell on peer pressure as the main one. The total population is given by The susceptible is increased by a constant inflow into the population at rate . A fraction of these adolescents are assumed to be nonsubstance abusers and the complementary fraction being substance abusers. This is reasonable because the underlying population is made of adolescents only into these two classes (influence from some peers is probably the main drive behind the adolescents abusing substances and hence their recruitment into the substance abusers class). It has been found that adolescents start abusing substances for four main reasons: to improve their mood, to receive social rewards, to reduce negative feelings, and to avoid social rejection [45], and these make them fall victims to peer influence easily. Individuals in different human subgroups suffer from natural death, at rate . Assuming homogeneous mixing of the population, the susceptible individuals acquire HSV-2 infection at rate , with is a modification parameter accounting for reduced number a nonsubstance abuser would infect as compared to a substance abuser [15, 20], perhaps because adolescents who engage in substance abuse also engage in other risky behaviours. Substance abuse diminishes judgement, leading to the engagement in the high risk behaviours [46], and the same explanation also applies for the modification parameter . is the probability of being infected from a sexual partner and is the rate at which an individual acquires sexual partners per unit time and is the effective contact rate for HSV-2 infection (contact sufficient to result in HSV-2 infection). The role of condom use is represented by . If , then condom use is not effective; corresponds to completely effective condom use, while implies that condom use would be effective to some degree. Substance abusers experience related diseases like liver cirrhosis, pancreatitis, and heart disease, among other substance abuse induced diseases and die at a rate because of these diseases [28]. Individuals are assumed to become substance abusers mostly by peer pressure, maybe by choice, or any other reasons at rate , and likewise individuals may cease being substance abusers due to educational campaigns (mostly), health reasons, pressure at work or school, poverty, or any other reason at rate . The model flow diagram is depicted in Figure 1.

From the descriptions and assumptions on the dynamics of the epidemic made above, the following are the model equations:

2.1. Positivity and Boundedness of Solutions

Model system (3) describes human population and therefore it is necessary to prove that all the variables , , , and are nonnegative for all time. Solutions of the model system (3) with positive initial data remain positive for all time and are bounded in .

Theorem 1. Let , , , and . The solutions , , , and of model system (3) are positive for all . For model system (3), the region is positively invariant and all solutions starting in approach enter or stay in .

Proof. Under the given conditions, it is easy to prove that each solution component of model system (3) is positive; otherwise assume by contradiction that there exists a first time, and , , for  or there exists a,, and , , for  or there exists a,, and , , for  or there exists a,, and , , for . In the first case we have, which is a contradiction and consequently remains positive. In the second case we have which is a contradiction; therefore is positive. In the third case we have which is a contradiction, therefore is positive. In the final case we have which is a contradiction; that is, remains positive. Thus, in all cases , , , and remain positive.
Since and , then implies that is bounded and all solutions starting in approach enter or stay in .

For system (3), the first octant in the state space is positively invariant and attracting; that is, solutions that start in this octant where all the variables are nonnegative stay there. Thus, system (3) will be analysed in a suitable region which is positively invariant and attracting. It is sufficient to consider solutions in where existence, uniqueness, and continuation results for system (3) hold.

3. Model Analysis

3.1. Disease-Free Equilibrium and Its Stability Analysis

Model system (3) has the disease-free equilibrium: The linear stability of is governed by the basic reproductive number which is defined as the spectral radius of an irreducible or primitive nonnegative matrix [47]. Biologically, the basic reproductive number can be interpreted as the expected number of secondary infections produced by a single infectious individual during his/her infectious period when introduced in a completely susceptible population. Using the notation in [48], the nonnegative matrix and the nonsingular matrix for the new infection terms and the remaining transfer terms are, respectively, given (at DFE) by We now consider different possibilities in detail.

Case a (there are no substance abusers in the community). We set and . Then the spectral radius is given by This reproductive rate sometimes referred to as the back of the napkin [49] is simply the ratio of the per capita rate of infection and the average lifetime of an individual in class . It is defined as the number of secondary HSV-2 cases produced by a single infected individual during his/her entire infectious period in a totally naive (susceptible) population in the absence of substance abusers.

Case b (the entire population is made up of substance abusers). We set and . Then the spectral radius is given by biologically, measures the average number of new infections generated by a single HSV-2 infective who is also a substance abuser during his/her entire infectious period when he/she is introduced to a susceptible population of substance abusers.

Case c (there is no condom use). We set . Then the spectral radius is given by where biologically, measures the average number of new secondary HSV-2 infections generated by a single HSV-2 infective during his/her infectious period whilst in a community of nonsubstance abusers and substance abusers, in the absence of condom use.

Case d (the general case). Existence of nonsubstance abusers and substance abusers. Then the spectral radius is given by where biologically, measures the average number of new HSV-2 infections generated by a single HSV-2 infective during his/her infectious period within a community, in the presence of condom use, nonsubstance abusers, and substance abusers. Theorem 2 follows from van den Driessche and Watmough [48].

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

Using a proof based on the comparison theorem, we can even show the global stability of the disease-free equilibrium in the case that the disease-free equilibrium is less than unity.

3.1.1. Global Stability of the Disease-Free Equilibrium

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

Proof. The proof is based on using a comparison theorem [50]. Note that the equations of the infected components in system (3) can be written as where and are as defined earlier in (11) and (12), respectively, and also where Since (full proof for see Appendix A) (for all time ) in , it follows that Using the fact that the eigenvalues of the matrix all have negative real parts, it follows that the linearized differential inequality system (21) is stable whenever . Consequently, as . Thus, by a comparison theorem [50] as , and evaluating system (3) at gives and for . Hence, the disease-free equilibrium (DFE) is globally asymptotically stable for .

3.2. Endemic Equilibrium and Stability Analysis

Model system (3) has three possible endemic equilibria: the (HSV-2 only) substance abuse-free endemic equilibria with a population of nonsubstance abusers only, the endemic equilibrium when the whole population is made up of substance abusers only and lastly the equilibrium where nonsubstance abusers and substance abusers co-exist, herein referred to as the interior equilibrium point.

3.2.1. Nonsubstance Abusers (Only) Endemic Equilibrium

We set and , so that there are no substance abusers in the community. Thus, model system (3) reduces to For system (22), it can be shown that the region is invariant and attracting. Thus, the dynamics of nonsubstance abusing (only) adolescents would be considered in .

Summing the two equation in system (22) above, we get . Furthermore, from the second equation of system (22), or with as defined in (13).

Since , hence , and thus Thus, system (22) has an endemic equilibrium which makes biological sense only when . This leads to Theorem 4 below.

Theorem 4. The substance abuse-free endemic equilibrium exists whenever .

The local stability of is given by the Jacobian evaluated at the point: From (22) we note that , which is obtained from the second equation of system (22) with the right-hand side set equaling zero. Thus, can be written as Hence the characteristic polynomial of the linearised system is given by Thus, is locally asymptotically stable for . We summarise the result in Theorem 5.

Theorem 5. The endemic equilibrium is locally asymptotically stable whenever .

We now study the global stability of system (22) using the Poincaré-Bendixon Theorem [51]. We claim the following result.

Theorem 6. The endemic equilibrium of the HSV-2 only model system (22) is globally asymptotically stable in whenever .

Proof. From model system (22), we note that . Using the fact that , it follows that , and, substituting this result into (22), we obtain Using Dulac’s multiplier , it follows that Thus, by Dulac’s criterion, there are no periodic orbits in . Since is positively invariant and the endemic equilibrium exists whenever , then it follows from the Poincaré-Bendixson Theorem [51] that all solutions of the limiting system originating in remain in , for  all . Further, the absence of periodic orbits in implies that the endemic equilibrium of the HSV-2 only model is globally asymptotically stable.

3.2.2. Substance Abusers (Only) Endemic Equilibrium

This occurs when the whole community consists of substance abusers only. Using the similar analysis as in Section 3.2.1, it can be easily shown that the endemic equilibrium is with as defined in (14), and it follows that the endemic equilibrium makes biological sense whenever . Furthermore, using the analysis done in Section 3.2.1, the stability of can be established. We now discuss the coexistence of nonsubstance abusers and substance abusers in the community.

3.2.3. Interior Endemic Equilibrium

This endemic equilibrium refers to the case where both the nonsubstance abusers and the substance abusers coexist. This state is denoted by with We now study the global stability of the interior equilibrium of system (3) using the Poincaré-Bendinxon Theorem (see Appendix B). The permanence of the disease destabilises the disease-free equilibrium since , and, thus, the equilibrium exists.

Lemma 7. Model system (3) is uniformly persistent on .

Proof. Uniform persistence of system (3) implies that there exists a constant such that any solution of model (3) starts in , and the interior of satisfies Define the following Korobeinikov and Maini [52] type of Lyapunov functional: which is continuous for all and satisfies (see Korobeinikov and Wake [53] for more details). Consequently, the interior equilibrium is the only extremum and the global minimum of the function . Also, and only at . Then, substituting our interior endemic equilibrium identities into the time derivative of along the solution path of the model system (3), we have The function, , can be shown to be nonpositive using Barbalat Lemma [54] or the approach in McCluskey [55]. Hence, with equality only at . The only invariant set in is the set consisting of the equilibrium . Thus, all solutions of model system (3) which intersect with limit to an invariant set, the singleton . Therefore, from the Lyapunov-Lasalle invariance principle, model system (3) is uniformly persistent.

From the epidemiological point of view, this result means that the disease persists at endemic state, which is not good news from the public health viewpoint.

Theorem 8. If , then the interior equilibrium is globally asymptotically stable in the interior of .

Proof. Denote the right-hand side of model system (3) by and for and , respectively, and choose a Dulac function as . Then we have Thus, system (3) does not have a limit cycle in the interior of . Furthermore, the absence of periodic orbits in implies that is globally asymptotically stable whenever .

4. Numerical Results

In order to illustrate the results of the foregoing analysis, we have simulated model (3) using the parameters in Table 1. Unfortunately, the scarcity of the data on HSV-2 and substance abuse in correlation with a focus on adolescents limits our ability to calibrate; nevertheless, we assume some of the parameters in the realistic range for illustrative purpose. These parsimonious assumptions reflect the lack of information currently available on HSV-2 and substance abuse with a focus on adolescents. Reliable data on the risk of transmission of HSV-2 among adolescents would enhance our understanding and aid in the possible intervention strategies to be implemented.

Figure 2 shows the effects of varying the effective contact rate on the reproduction number. Results on Figure 2 suggest that the increase or the existence of substance abusers may increase the reproduction number. Further analysis of Figure 2 reveals that the combined use of condoms and educational campaigns as HSV-2 intervention strategies among substance abusing adolescents in a community have an impact on reducing the magnitude of HSV-2 cases, as shown by the lowest graph for .

Figure 3(a) depicts that the reproduction number can be reduced by low levels of adolescents becoming substance abusers and high levels of adolescents quitting the substance abusive habits. Numerical results in Figure 3(b) suggest that more adolescents should desist from substance abuse, since the trend shows a decrease in the reproduction number with increasing . Trend line shows the effects of increasing which results in an increase of the reproductive number thereby increasing the prevalence of HSV-2 in community. The point depicts the balance in transfer rates from being a substance abuser to a nonsubstance abuser or from a nonsubstance abuser to a substance abuser , and this suggests that changes in the substance abusive habit will not have an impact on the reproduction number. Point also suggests that, besides substance abuse among adolescents, other factors are influencing the spread of HSV-2 amongst them in the community. Measures such as educational campaigns and counselling may be introduced so as to reduce while increasing .

In many epidemiological models, the magnitude of the reproductive number is associated with the level of infection. The same is true in model system (3). Sensitivity analysis assesses the amount and type of change inherent in the model as captured by the terms that define the reproductive number . If is very sensitive to a particular parameter, then a perturbation of the conditions that connect the dynamics to such a parameter may prove useful in identifying policies or intervention strategies that reduce epidemic prevalence.

Partial rank correlation coefficients (PRCCs) were calculated to estimate the degree of correlation between values of and the ten model parameters across 1000 random draws from the empirical distribution of and its associated parameters.

Figure 4 illustrates the PRCCs using as the output variable. PRCCs show that the rate of becoming substance abusers by the adolescents has the most impact in the spread of HSV-2 among adolescents. Effective contact rate is found to support the spread of HSV-2 as noted by its positive PRCC. However, condom use efficacy is shown to have a greater impact on reducing .

Since the rate of becoming a substance abuser , condom use efficacy , and the effective contact rate have significant effects on , we examine their dependence on in more detail. We used Latin hypercube sampling and Monte Carlo simulations to run 1000 simulations, where all parameters were simultaneously drawn from across their ranges.

Figure 5 illustrates the effect that varying three sample parameters will have on . In Figure 5(a), if the rate of becoming a substance abuser is sufficiently high, then and the disease will persist. However, if the rate of becoming a substance abuser is low, then and the disease can be controlled. Figure 5(b) illustrates the effect of condom use efficacy on controlling HSV-2. The result demonstrates that an increase in results in a decrease in the reproduction number. Figure 5(b) demonstrates that if condom use efficacy is more than 70% of the time, then the reproduction number is less than unity, and then the disease will be controlled.

Numerical results on Figure 8 are in agreement with the earlier findings, which have shown that condom use has a significant effect on HSV-2 cases among adolescents. However, Figure 8(b) shows that, for more adolescents leaving their substance abuse habit than becoming substance abusers, condom use is more effective. Figure 8(b) also shows that more adolescents should desist from substance abuse to reduce the HSV-2 cases. It is also worth noting that, as the condom use efficacy is increasing, the cumulative HSV-2 cases also are reduced.

In general, simulations in Figure 7 suggest that an increase in recruitment of substance abusing adolescents into the community will increase the prevalence of HSV-2 cases. It is worth noting that the prevalence is much higher when we have more adolescents becoming substance abusers than quitting the substance abusing habit, that is, for .

5. Optimal Condom Use and Educational Campaigns

In this section our goal is to solve the following problem: given initial population sizes of all the four subgroups, , , , and , find the best strategy in terms of combined efforts of condom use and educational campaigns that would minimise the cumulative HSV-2 cases while at the same time minimising the cost of condom use and educational campaigns. Naturally, there are various ways of expressing such a goal mathematically. In this section, for a fixed , we consider the following objective functional: subject to The parameters , , and represent the weight constants of the susceptible substance abusers, the infectious nonsubstance abusers, and the infectious substance abusers, while the parameters and are the weights and cost of condom use and educational campaigns for the controls and , respectively. The terms and reflect the effectiveness of condom use and educational campaigns in the control of the disease. The values represent maximum effectiveness of condom use and educational campaigns. We therefore seek an optimal control and such that where is measurable: , , .

The basic framework of this problem is to characterize the optimal control and prove the existence of the optimal control and uniqueness of the optimality system.

In our analysis, we assume , , as applied in [56]. Thus the total population is constant. We can also treat the nonconstant population case by these same techniques, but we choose to present the constant population case here.

5.1. Existence of an Optimal Control

The existence of an optimal control is proved by a result from Fleming and Rishel (1975) [57]. The boundedness of solutions of system (39) for a finite interval is used to prove the existence of an optimal control. To determine existence of an optimal control to our problem, we use a result from Fleming and Rishel (1975) (Theorem 4.1 pages 68-69), where the following properties are satisfied.(1)The class of all initial conditions with an optimal control set and in the admissible control set along with each state equation being satisfied is not empty.(2)The control set is convex and closed.(3)The right-hand side of the state is continuous, is bounded above by a sum of the bounded control and the state, and can be written as a linear function of each control in the optimal control set and with coefficients depending on time and the state variables.(4)The integrand of the functional is concave on and is bounded above by , where .

An existence result in Lukes (1982) [58] (Theorem ) for the system of (39) for bounded coefficients is used to give condition 1. The control set is closed and convex by definition. The right-hand side of the state system equation (39) satisfies condition 3 since the state solutions are a priori bounded. The integrand in the objective functional, , is Lebesgue integrable and concave on . Furthermore, and , hence satisfying hence the optimal control exists, since the states are bounded.

5.2. Characterisation

Since there exists an optimal control for maximising the functional (42) subject to (39), we use Pontryagin’s maximum principle to derive the necessary conditions for this optimal control. The Lagrangian is defined as where , , , and are penalty multipliers satisfying , , and , at the optimal point and .

Theorem 9. Given optimal controls and and solutions of the corresponding state system (39), there exist adjoint variables , satisfying

Proof. The form of the adjoint equation and transversality conditions are standard results from Pontryagin’s maximum principle [57, 59]; therefore, solutions to the adjoint system exist and are bounded. To determine the interior maximum of our Lagrangian, we take the partial derivative with respect to and and set to zero. Thus, we have the following.
is subject of formulae
is subject of formulae
To determine the explicit expression for the control without and , a standard optimality technique is utilised. The specific characterisation of the optimal controls and . Consider The optimality system consists of the state system coupled with the adjoint system with the initial conditions, the transversality conditions and the characterisation of the optimal control. Substituting , and for , and , in (43) gives the optimality system. The state system and adjoint system have finite upper bounds. These bounds are needed in the uniqueness proof of the optimality system. Due to a priori boundedness of the state and adjoint functions and the resulting Lipschitz structure of the ODEs, we obtain the uniqueness of the optimal control for small [60]. The uniqueness of the optimal control follows from the uniqueness of the optimality system.

5.3. Numerical Simulations

The optimality system is solved using an iterative method with Runge-Kutta fourth-order scheme. Starting with a guess for the adjoint variables, the state equations are solved forward in time. Then these state values are used to solve the adjoint equations backward in time, and the iterations continue until convergence. The simulations were carried out using parameter values in Table 1 and the following are the values , , , , and . The assumed initial conditions for the differential equations are , , , and .

Figure 8(a) illustrates the population of the susceptible nonsubstance abusing adolescents in the presence and absence of the controls over a period of 20 years. For the period 0–3 years, in the absence of the control, the susceptible nonsubstance abusing adolescents have a sharp increase, and for the same period when there is no control the population of the nonsubstance abusing adolescents decreases sharply. Further, we note that for the case when there is no control, for the period 3–20 years, the population of the susceptible nonsubstance abusing adolescents increases slightly.

Figure 8(b) illustrates the population of the susceptible substance abusers adolescents in the presence and absence of the controls over a period of 20 years. For the period 0–3 years, in the presence of the controls, the susceptible substance abusing adolescents have a sharp decrease, and for the same period when there are controls the population of the susceptible substance abusers increases sharply. It is worth noting that, after 3 years, the population for the substance abusers with control remains constant while the population of the susceptible substance abusers without control continues increasing gradually.

Figure 8(c) highlights the impact of controls on the population of infectious nonsubstance abusing adolescents. For the period 0–2 years, in the absence of the control, cumulative HSV-2 cases for the nonsubstance abusing adolescents decrease rapidly, before starting to increase moderately for the remaining years under consideration. The cumulative HSV-2 cases for the nonsubstance abusing adolescents decrease sharply for the whole 20 years under consideration. It is worth noting that, for the infectious nonsubstance abusing adolescents, the controls yield positive results after a period of 2 years; hence the controls have an impact in controlling cumulative HSV-2 cases, within a community.

Figure 8(d) highlights the impact of controls on the population of infectious substance abusing adolescents. For the whole period under investigation, in the absence of the control, cumulative infectious substance abusing HSV-2 cases increase rapidly, and, for the similar period, when there is a control the population of the exposed individuals decreases sharply. Results on Figure 8(d) clearly suggest that the presence of controls has an impact on reducing cumulative HSV-2 cases within a society with substance abusing adolescents.

Figure 9 represents the controls and . The condom use control is at the upper bound for approximately 20 years and has a steady drop until it reaches the lower bound . Treatment control, , is at the upper bound for approximately half a year and has a sharp drop until it reaches the lower bound .

These results suggest that more effort should be devoted to condom use , as compared to educational campaigns .

6. Discussion

Adolescents who abuse substances have higher rates of HSV-2 infection, since adolescents who engage in substance abuse also engage in risky behaviours; hence most of them do not use condoms. In general adolescents are not fully educated on HSV-2. A mathematical model for investigating the effects of substance abuse amongst adolescents on the transmission dynamics of HSV-2 in the community is formulated and analysed. We computed and compared the reproduction numbers. Results from the analysis of the reproduction numbers suggest that condom use and educational campaigns have a great impact in the reduction of HSV-2 cases in the community. The PRCCs were calculated to estimate the correlation between values of the reproduction number and the epidemiological parameters. It has been seen that an increase in condom use and a reduction on the number of adolescents becoming substance abusers have the greatest influence on reducing the magnitude of the reproduction number than any other epidemiological parameters (Figure 6). This result is further supported by numerical simulations which show that condom use is highly effective when we are having more adolescents quitting substance abusing habit than those becoming ones, whereas it is less effective when we have more adolescents becoming substance abusers than quitting. Substance abuse amongst adolescents increases disease transmission, and prevalence of disease increases with increased rate of becoming a substance abuser. Thus, in the event when we have more individuals becoming heavy substance abusers, there is urgent need for intervention strategies such as counselling and educational campaigns in order to curtail the HSV-2 spread. Numerical simulations also show and suggest that an increase in recruitment of substance abusers in the community will increase the prevalence of HSV-2 cases.

The optimal control results show how a cost effective combination of aforementioned HSV-2 intervention strategies (condom use and educational campaigns) amongst adolescents may influence cumulative cases over a period of 20 years. Overall, optimal control theory results suggest that more effort should be devoted to condom use compared to educational campaigns.

Appendices

A. Positivity of

To show that , it is sufficient to show that To do this we prove it by contradiction. Assume the following statements are true: Adding the inequalities (i) and (ii) in (A.2), we have that which is not true, hence a contradiction (since a number cannot be less than itself).

Equations (A.2) and (A.5) should justify that . In order to prove this we shall use the direct proof; for example, if Adding items (i) and (ii) one gets .

Using the above argument if and only if Adding items (i) and (ii) one gets Thus, from the above argument .

B. Poincaré-Bendixson’s Negative Criterion

Consider the nonlinear autonomous system where . If , let be the solution of (B.1) which satisfies . Bendixon’s negative criterion: when , a sufficient condition for the nonexistence of nonconstant periodic solutions of (B.1) is that, for each , div.

Theorem B.1. Suppose that one of the inequalities , holds for all . Then the system (B.1) has nonconstant periodic solutions.

For a detailed discussion of the Bendixon’s negative criterion, we refer the reader to [51, 61].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the handling editor and anonymous reviewers for their comments and suggestions which improved the paper.