HSV-2 and Substance Abuse amongst Adolescents : Insights through Mathematical Modelling

Herpes simplex virus infection is mostly spread and occurs more commonly among substance abusing adolescents as compared to the nonsubstance abusing. In this paper, amathematical model for the spread ofHSV-2within 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 onmaximising the susceptible nonsubstance abusing adolescents, whileminimising 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.


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 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.

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  0 of these adolescents are assumed to be nonsubstance abusers and the complementary fraction  1 = 1 −  0 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 ∈ (0, 1) 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  ∈ (0, 1). 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  = 0, then condom use is not effective;  = 1 corresponds to completely effective condom use, while 0 <  < 1 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 V 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  ≥ 0 and are bounded in G.
The solutions   ,   ,   , and   of model system (3) are positive for all  ≥ 0. For model system (3), the region G is positively invariant and all solutions starting in G approach enter or stay in G.
Since () ≥   () and () ≥   (), then implies that () is bounded and all solutions starting in G approach enter or stay in G.
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 G where existence, uniqueness, and continuation results for system (3) hold.

Model Analysis
3.1.Disease-Free Equilibrium and Its Stability Analysis.Model system (3) has the disease-free equilibrium:

Journal of Applied Mathematics
The linear stability of E 0 is governed by the basic reproductive number  0 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 V =  =  =  1 = 0 and  0 = 1.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  =  =  0 = 0 and  1 = 1.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  = 0. 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   ≤ 1, 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.

Global Stability of the Disease-Free Equilibrium
Theorem 3. The disease-free equilibrium E 0 of system ( 3) is globally asymptotically stable (GAS) if   ≤ 1 and unstable if   > 1.
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  ≥ 0 (full proof for  ≥ 0: 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   ≤ 1.

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.

Nonsubstance Abusers (Only) Endemic
Equilibrium.We set  1 =  =  = 0 and  0 = 1, 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 G  .Summing the two equation in system ( 22) above, we get   = Λ/.Furthermore, from the second equation of system (22), with   as defined in (13).

The local stability of E *
is given by the Jacobian evaluated at the point: From ( 22) we note that (1 − )(1 − ) *  / *  = , which is obtained from the second equation of system (22) with the right-hand side set equaling zero.Thus, (E *  ) can be written as Hence the characteristic polynomial of the linearised system is given by Thus, E *  is locally asymptotically stable for   > 1.We summarise the result in Theorem 5.
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 G  whenever   > 1.
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 1/  , it follows that Thus, by Dulac's criterion, there are no periodic orbits in G  .Since G  is positively invariant and the endemic equilibrium exists whenever   > 1, then it follows from the Poincaré-Bendixson Theorem [51] that all solutions of the limiting system originating in G remain in G  , for all  ≥ 0.
Further, the absence of periodic orbits in G  implies that the endemic equilibrium of the HSV-2 only model is globally asymptotically stable.

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 E *  makes biological sense whenever   > 1.Furthermore, using the analysis done in Section 3.2.1, the stability of E *  can be established.We now discuss the coexistence of nonsubstance abusers and substance abusers in the community.

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 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 E 0 since   > 1, and, thus, the equilibrium E *  exists.
Proof.Uniform persistence of system (3) implies that there exists a constant  > 0 such that any solution of model ( 3) starts in ∘ G, and the interior of G satisfies Define the following Korobeinikov and Maini [52] type of Lyapunov functional: which is continuous for all (  ,   ,   ,   ) > 0 and satisfies (see Korobeinikov and Wake [53] for more details).Consequently, the interior equilibrium E *  is the only extremum and the global minimum of the function  ∈ R 4 + .Also, (  ,   ,   ,   ) > 0 and   (  ,   ,   ,   ) = 0 only at E *  .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,   (  ,   ,   ,   ) ≤ 0 with equality only at E *  .The only invariant set in ∘ G is the set consisting of the equilibrium E *  .Thus, all solutions of model system (3) which intersect with ∘ G limit to an invariant set, the singleton {E *  }.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   > 1, then the interior equilibrium E *  is globally asymptotically stablein the interior of G.
Thus, system (3) does not have a limit cycle in the interior of G. Furthermore, the absence of periodic orbits in G implies that E *  is globally asymptotically stable whenever   > 1.

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 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  1 with  and , varying from 0.0 to 1.0.
(  ).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   > 1 and the disease will persist.However, if the rate of becoming a substance abuser is low, then   > 1 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.1. Parameters with positive PRCCs will increase   when they are increased, whereas parameters with negative PRCCs will decrease   when they are increased.
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  > .

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  1 ,  2 , and  3 represent the weight constants of the susceptible substance abusers, the infectious nonsubstance abusers, and the infectious substance abusers, while the parameters  1 and  2 are the weights and cost of condom use and educational campaigns for the controls  1 and  2 , respectively.The terms  2 1 and  2 2 reflect the effectiveness of condom use and educational campaigns in the control of the disease.The values  1 =  2 = 1 represent maximum effectiveness of condom use and educational campaigns.We therefore seek an optimal control  * 1 and  * 2 such that where 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 Λ = , V = 0, 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.

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  1 and  2 in the admissible control set along with each state equation being satisfied is not empty.
(2) The control set U 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  1 and  2 with coefficients depending on time and the state variables.
(4) The integrand of the functional is concave on U and is bounded above by An existence result in Lukes (1982) [58] (Theorem 9.2.1) 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 U. Furthermore,  1 ,  2 > 0 and  1 ,  2 ,  3 ,  1 ,  2 > 1, hence satisfying hence the optimal control exists, since the states are bounded.
Theorem 9. Given optimal controls  * 1 and  * 2 and solutions of the corresponding state system (39), there exist adjoint variables   ,  = 1, . . ., 4 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  1 () and  2 () and set to zero.Thus, we have the following. 1 () * is subject of formulae To determine the explicit expression for the control without  1 () and  2 (), a standard optimality technique is utilised.The specific characterisation of the optimal controls  1 () * and  2 () * .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  * 1 , and  * 2 for  1 (), and  2 (), 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.

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  1 = 0.005,  2 = 0.005,  1 = 0.06,  2 = 0.05, and  3 = 0.05.The assumed initial conditions for the differential equations are   0 = 0.4,   0 = 0.3,   0 = 0.2, and   0 = 0.1.
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  1 and  2 .The condom use control  1 is at the upper bound  1 for approximately 20 years and has a steady drop until it reaches the lower bound  1 = 0. Treatment control,  2 , is at the upper bound for approximately half a year and has a sharp drop until it reaches the lower bound  2 = 0.
These results suggest that more effort should be devoted to condom use ( 1 ), as compared to educational campaigns ( 2 ).

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.

Figure 1 :
Figure 1: Structure of the model.

Figure 2 :
Figure 2: Effects of varying the effective contact rate () on the reproduction numbers.

Figure 3 :
Figure 3: (a)  The relationship between the reproduction number (  ), rate of becoming a substance abuser (), and rate of quitting substance abuse ().(b) Numerical results of model system(3) showing the effect in the change of substance abusing habits.Parameter values used are in Table1with  and , varying from 0.0 to 1.0.

Figure 4 :
Figure 4: Partial rank correlation coefficients showing the effect of parameter variations on   using ranges in Table1.Parameters with positive PRCCs will increase   when they are increased, whereas parameters with negative PRCCs will decrease   when they are increased.

Figure 5 :
Figure 5: Monte Carlo simulations of 1000 sample values for three illustrative parameters (, , and ) chosen via Latin hypercube sampling.

Figure 6 :
Figure 6: Simulations of model system (3) showing the effects of an increase on condom use efficacy on the HSV-2 cases (  +   ) over a period of 60 years.The rest of the parameters are fixed on their baseline values, with Figure 8(a) having  >  and Figure 8(b) having  > .

Figure 7 :
Figure 7: Simulations of model system (3) showing the effects of increasing the proportion of recruited substance abusers ( 1 ), varying from 0.0 to 1.0 with a step size of 0.25, for the 2 cases where (a)  >  and (b)  > .

Table 1 :
Model parameters and their interpretations.Denote the right-hand side of model system (3) by Φ and Ψ for   and   , respectively, and choose a Dulac function as (  ,   ) = 1/    .Then we have Time series plots showing the effects of optimal control on the susceptible nonsubstance abusing adolescents   , susceptible substance abusing adolescents   , infectious nonsubstance abusing adolescents   , and the infectious substance abusing   , over a period of 20 years.