Computational and Mathematical Methods in Medicine
Volume 2020 (2020), Article ID 5984095, 15 pages
https://doi.org/10.1155/2020/5984095

Research Article

Application of Optimal Control to Influenza Pneumonia Coinfection with Antiviral Resistance

Caroline W. Kanyiri,¹ Livingstone Luboobi,² and Mark Kimathi³

¹Department of Mathematics, Pan African University Institute of Basic Sciences, Technology and Innovation, P.O. Box 62000-00200, Nairobi, Kenya
²Institute of Mathematical Sciences, Strathmore University, P.O. Box 59857-00200, Nairobi, Kenya
³Department of Mathematics, Machakos University, P.O. Box 139-90100, Machakos, Kenya

Correspondence should be addressed to Caroline W. Kanyiri; kanyiricarolyne2@gmail.com

Received 20 November 2019; Revised 1 February 2020; Accepted 13 February 2020; Published 11 March 2020

Academic Editor: Chung-Min Liao

Copyright © 2020 Caroline W. Kanyiri 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.

Influenza and pneumonia independently lead to high morbidity and mortality annually among the human population globally; however, a glaring fact is that influenza pneumonia coinfection is more vicious and it is a threat to public health. Emergence of antiviral resistance is a major impediment in the control of the coinfection. In this paper, a deterministic mathematical model illustrating the transmission dynamics of influenza pneumonia coinfection is formulated having incorporated antiviral resistance. Optimal control theory is then applied to investigate optimal strategies for controlling the coinfection using prevalence reduction and treatment as the system control variables. Pontryagin’s maximum principle is used to characterize the optimal control. The derived optimality system is solved numerically using the Runge–Kutta-based forward-backward sweep method. Simulation results reveal that implementation of prevention measures is sufficient to eradicate influenza pneumonia coinfection from a given population. The prevention measures could be social distancing, vaccination, curbing mutation and reassortment, and curbing interspecies movement of the influenza virus.

1. Introduction

Clinical evidence points out that infection with a particular combination of pathogens results in an aggravated infection with more severe clinical outcome compared with infection with either pathogen alone [1]. This is specially true for influenza virus and bacterium Streptococcus pneumoniae [2–4]. Influenza and pneumonia each contributes greatly to the global burden of morbidity and leads to high death toll, typically over a relatively short period of time [5–8]. Streptococcus pneumoniae, Haemophilus influenzae, and Staphylococcus aureus are the most common causes of pneumonia, the chief being Streptococcus pneumoniae [9, 10]. Coinfection resulting from influenza virus and Streptococcus pneumoniae further increases morbidity and mortality and is a major public health concern. These two pathogens rank among the chief pathogens affecting humans, and their ability to work together presents a major threat to world health [11]. Coinfection greatly impairs the host’s immune system, increases antibacterial therapy intolerance, and can be detrimental to the diagnosis of the disease [12]. According to [13], it can be difficult to identify influenza patients experiencing bacterial coinfections due to symptom overlap of influenza and bacterial infections. In [14], it is indicated that a strong index of suspicion and additional diagnostic testing may be required for diagnosis and treatment of the infections.

The morbidity, mortality, and economic burden resulting from the lethal synergism that exists between influenza virus and pneumococcus are of major concern globally. The catastrophic 1918 influenza pandemic is an extreme example of the impact that results from this cooperative interaction [11]. Lung tissue samples examined from those who died during this pandemic revealed that the majority of deaths were as a result of secondary bacterial pneumonia. Data from the subsequent 1957, 1968, and 2009 influenza pandemics are consistent with these findings [15, 16]. In addition, during seasonal influenza outbreaks, coinfections resulting from influenza and Streptococcus pneumoniae have been associated with high morbidity and mortality rates [17–19]. According to [11], influenza virus alters the lungs in such a way that predisposes them to invasion by pneumococcus rendering a mild influenza infection severe or even fatal. This could be through several ways such as epithelial damage, changes in airway function, upregulation and exposure of receptors, dampening of the immune response, or amplification of inflammation. Several studies have been carried out to investigate the time course of susceptibility to Streptococcus pneumoniae after influenza virus infection. Results revealed that on average, individuals developed coinfection within 6 days after influenza virus infection [20–23].

Emergence of drug resistance, which has become a global concern, complicates influenza pneumonia coinfection even more. Drug resistance refers to the ability of disease-causing agents to resist the effects of drugs, thereby making the conventional treatment procedure ineffective. This leads to persistence of infections in the body, hence increasing the risk of spread to other individuals [24, 25]. The evolution of drug resistance is accelerated by overuse and misuse of antimicrobials, inappropriate use of antimicrobials, subtherapeutic dosing, and patient noncompliance with the recommended course of treatment [26]. There are two classes of antiviral drugs that are approved to treat influenza infections; these are M2 ion-channel inhibitors and neuraminidase (NA) inhibitors. However, due to antiviral drug resistance in influenza virus, neuraminidase (NA) inhibitors are the only class of antiinfluenza drugs currently in use as most of the circulating influenza viruses have acquired resistance to M2 ion-channel inhibitors [27, 28]. Moreover, many circulating influenza viruses have also acquired resistance to neuraminidase (NA) inhibitors [28, 29] raising an alarm in the health sector. Drug resistance continues to threaten effective prevention and treatment of influenza pneumonia coinfection. In addition, the cost of health care for patients with resistant infections is much higher than care for patients with nonresistant infections especially due to longer duration of illness.

Strategies such as vaccination, isolation, and treatment among others are necessary in order to curb the spread of various infectious diseases. However, if they are not administered at the right time and in the right amount, curtailing the spread of the infectious diseases remains a difficult task. The application of optimal control is therefore very vital since it is a necessary tool in making decisions of the viable control strategies to be employed in eradicating diseases. Optimal control theory has been applied in the study of influenza, for instance, in [30–34], and in the study of pneumonia, for instance, in [35]. Given that influenza pneumonia coinfection is more disastrous than either of the single infections alone, this paper seeks to investigate optimal control strategies for the influenza pneumonia coinfection and in the emergence of antiviral resistance.

2. Model Formulation

The model presented in this paper has the total population subdivided into eight compartments. These are susceptible (S), infected with wild-type influenza strain (), infected with resistant influenza strain (), infected with both wild-type influenza strain and pneumonia (), infected with both resistant influenza strain and pneumonia (), recovered from influenza (), recovered from pneumonia (), and recovered from both influenza and pneumonia (). In the model, individuals are first infected with influenza virus and then contract bacterial pneumonia. An assumption is made that there is no primary bacterial pneumonia infection. Individuals enter the population via immigration at the rate of π, and all recruited individuals are assumed to be susceptible. The susceptible can either become infected with wild-type influenza strain or resistant influenza strain at the rates of and , respectively, where , while , where , in which parameter b refers to the rate of developing drug resistance while parameters and refer to the transmission rate of wild-type strain and resistant strain, respectively. Parameters and are modification parameters accounting for the relative infectiousness of individuals in and classes in comparison with those in and , respectively. Those infected with influenza wild-type strain can recover at a rate of while those infected with influenza resistant strain recover at the rate of . The model incorporates development of antiviral resistance; hence, individuals infected with influenza wild-type strain progress to the class at a rate of b. Those infected with influenza wild-type strain and those infected with influenza resistant strain can contract secondary bacterial pneumonia at the rate of with the force of infection of pneumonia being . Parameter refers to the transmission rate of pneumonia. Individuals in and classes also suffer disease-induced death at the rates of and , respectively. Individuals infected with both the wild-type influenza strain and pneumonia and those infected with both resistant strain and pneumonia can recover from either influenza alone at rates of and , respectively, or pneumonia alone at rates of and , respectively, or recover from both influenza and pneumonia at rates of and , respectively. This means that the class contains some individuals from and classes that have only recovered from influenza but may still have pneumonia. class is however considered noninfectious given that bacterial pneumonia is considered weakly infectious after antimicrobials are administered. On the other hand, the class contains some individuals from and classes that have only recovered from pneumonia but may still have influenza. Given that the infectious period of influenza is about seven days after symptoms onset, is considered noninfectious. Individuals in and suffer disease-induced death at the rate of and , respectively. Those infected with both wild-type influenza strain and pneumonia could also develop antiviral resistance and progress to class at a rate of b. Individuals who have recovered from influenza, pneumonia, and both influenza and pneumonia lose their immunity and become susceptible again at rates of , , and , respectively. Individuals in all the epidemiological compartments suffer natural death at the rate of μ. Figure 1 shows the population flow between the different compartments.

Figure 1: Schematic diagram showing population flow between different epidemiological classes for influenza pneumonia coinfection.

2.1. Model Equations

Given the dynamics described in Figure 1, the following system of nonlinear ordinary differential equations, with nonnegative initial conditions, describes the dynamics of influenza pneumonia coinfection:

where

We assume that all the model parameters are positive and the initial conditions of model system (1) are given by

Table 1 gives the description of the various parameters used in the model.

Table 1: Description of parameters used.

Parameters
Description

π
Recruitment rate

Transmission rate of wild-type influenza strain

Transmission rate of resistant influenza strain

Transmission rate of pneumonia

Recovery rate of influenza for individuals in class

Recovery rate of influenza for individuals in class

Recovery rate of influenza for individuals in class

Recovery rate of influenza for individuals in class

Recovery rate of pneumonia for individuals in class

Recovery rate of pneumonia for individuals in class

Recovery rate of both influenza and pneumonia for individuals in class

Recovery rate of both influenza and pneumonia for individuals in class

, ,
Rate of losing immunity for influenza, pneumonia, and influenza and pneumonia, respectively

b
Rate of developing antiviral resistance

, , ,
Disease-induced death rates in ,,, and classes, respectively

μ
Natural death rate

3. Coinfection Model with Controls

In order to identify optimal control strategies that minimize the number of infected individuals and the cost of implementing the controls, a mathematical optimal control problem is formulated and analysed. Influenza pneumonia coinfection model (1) is extended to include time-dependent control measures. Let , represent the time-dependent controls whereby controls , and relate to prevalence reduction of wild-type influenza strain, resistant influenza strain, and pneumonia, respectively. The infection prevalence reduction could be through social distancing, vaccination, curbing mutation and reassortment, and curbing interspecies movement of the influenza virus. Controls and relate to treatment of the wild-type and resistant influenza strains, respectively.

Time is specified and is given by , where T is the final time. Given that there is a limitation on the maximum rate of treatment and prevalence reduction controls, bounded Lebesgue measurable control set is introduced and defined as

To identify the required level of effort to control the infection, an objective functional to be minimized is given by

subject to the differential system:

The coefficients , and represent the costs associated with minimizing the infected population. On the other hand, the expression represents costs associated with controls . Quadratic expressions of the cost of controls are considered because costs follow a nonlinear representation especially at high intervention levels. The objective of minimizing the infected population and the cost of controls can be achieved through proper implementation of the controls over a time interval given by . Therefore, we seek an optimal control set such that

The necessary conditions for the existence of an optimal solution come from Pontryagin’s Maximum Principle [36]. This principle converts (5)–(6) into a problem of minimizing pointwise Hamiltonian H, with respect to , and is obtained as

where are the corresponding adjoint or costate variables to be determined by applying Pontryagin’s Maximum Principle as in Theorem 1.

Theorem 1. Given optimal control set and the corresponding solutions of system (6) that minimizes over , there exist adjoint variables , such that

with transversality condition for .
Evaluating (9) leads to the following adjoint system:

Following the results in [37, 38], the existence of optimal control is stated and proved using Theorem 2.

Theorem 2. There exists optimal controls which minimizes over the region satisfying the optimality condition:

where

Proof. The Hamiltonian is minimized with respect to the controls , and at the optimal control functions. This is done by differentiating the Hamiltonian function with respect to each of the control variables on the set ; that is,

The following set of optimality conditions is thus obtained:

at , and , respectively. Solving for , and and using the bounds for the controls in , that is,

in compact notation yields,

Next, the optimality system is obtained as

with for and .

4. Numerical Simulation

In this section, the optimal solution of optimality system (17) is investigated. Given that there are initial conditions for the state variables and terminal conditions for the adjoints, the optimality system is a two-point boundary value problem with separated boundary conditions at times and . The forward-backward sweep method described in [39, 40] is hence used for the numerical solution of the optimal control problem. The method is named based on how the algorithm solves the problem’s state and adjoint equations. The state variables are solved using the forward difference scheme while the adjoint variables are solved using the backward difference scheme. The solution iterative scheme makes an initial guess of the controls and using that guess solves the state system forward in time using a 4th-order Runge–Kutta scheme. Next, using the stored values of the controls and the solution of the state system, the adjoint system is solved backward in time using a 4th-order Runge–Kutta scheme. The RK4 method has to be adapted to account for solving backward in time. The controls are then updated using a convex combination of the previous controls and the values obtained using the characterizations. The updated controls are then used to repeat the solution of the state and adjoint systems. This process is repeated with values in the current iteration being tested for convergence against a user provided tolerance and depending on that, the algorithm either starts the process over again using the updated control or the algorithm terminates. The final approximations for the control, state, and adjoint systems are considered to be the solution to the optimal control problem.

4.1. Simulation Results

The numerical simulations were carried out using the MATLAB software and the parameter values in Table 2. There is however paucity of published estimates of morbidity and mortality rates resulting from influenza pneumonia coinfection in the middle- and low-income countries. On the other hand, in high-income, temperate countries where influenza surveillance has been done for years, these rates are well documented [47]; hence the data used in this research study are from these countries. The population under consideration is children 5 years and adults 60 years. Studies show that influenza pneumonia coinfection affects people of all ages; however, the morbidity and mortality rates of the high risk group which includes the immunocompromised, children 5 years, and adults 60 years are disproportionately high [9, 47, 48].

Table 2: Description and values of the different parameters used.

Parameter
Description
Value
Reference.

π
Recruitment rate
0.0381
Assumed

Transmission rate of wild-type influenza strain
0.0102 day
Assumed

Transmission rate of resistant influenza strain
0.00026 day
Assumed

Transmission rate of pneumonia
0.000162 day
Reference [41]

Recovery rate of influenza for individuals in class
0.07143 day
Reference [42]

Recovery rate of influenza for individuals in class
0.0333 day
Assumed

Recovery rate of influenza for individuals in class
0.04762 day
Reference [43].

Recovery rate of influenza for individuals in class
0.0222 day
Assumed

Recovery rate of pneumonia for individuals in class
0.033 day
Reference [41]

Recovery rate of pneumonia for individuals in class
0.033 day
Reference [41]

Recovery rate of both influenza and pneumonia for individuals in class
0.0166 day
Assumed

Recovery rate of both influenza and pneumonia for individuals in class
0.0166 day
Assumed

Rate of losing immunity for influenza
0.00833 day
Reference [44]

Rate of losing immunity for pneumonia
0.00833 day
Assumed

Rate of losing immunity for influenza pneumonia coinfection
0.00833 day
Assumed

b
Rate of developing antiviral resistance
0.0118
Assumed

Wild-type influenza strain-induced death rate
0.01
Reference [45]

Resistant influenza strain-induced death rate
0.021
Assumed

class disease-induced death rate
0.05
Assumed

class disease-induced death rates
0.05
Assumed

Average human lifespan
days
Reference [46]

Most of the parameter values have a range as indicated in the references given in Table 2. Given that this research study focuses on a deterministic mathematical model to give an indication of the likely dynamics of influenza pneumonia coinfection, baseline parameter values were used for the simulations.

4.1.1. Control by Treatment Only

Figures 2 and 3 show the effect of using treatment efforts only in an attempt to curb the spread of influenza pneumonia coinfection.

Figure 2: Individuals coinfected with wild-type influenza and pneumonia.

Figure 3: Individuals coinfected with resistant influenza and pneumonia.

From Figure 2, it can be observed that with and without the treatment efforts, the wild-type influenza and pneumonia coinfection persists in the population.

Without any controls, it can be observed from Figure 3 that the number of individuals coinfected with resistant influenza and pneumonia increases initially and then decreases slowly but not to zero. With the treatment control strategies in place, the number of the infected individuals is lower but not significant and does not decrease to zero by day 30. This shows that use of treatment alone as a control strategy is not effective in curbing the spread of influenza pneumonia coinfection. As explained in details in [4, 49–51], treatment is often rendered ineffective because it is not always possible to administer the drugs at the right time and the diagnosis of the coinfection can be challenging because of timing of sample collection and false negative results when viruses replicate in the lower respiratory tract. In addition, treatment often aggravates development of resistance.

4.1.2. Control by Prevention Measures Only

Simulations are done when there is no control strategy in place and when there are controls involving prevention of wild-type influenza strain, prevention of influenza resistant strain, and prevention of pneumonia. Figures 4 and 5 show the results.

Figure 4: Individuals coinfected with wild-type influenza and pneumonia.

Figure 5: Individuals coinfected with resistant influenza and pneumonia.

It can be observed from Figure 4 that when preventive efforts are implemented as control strategies, the number of individuals coinfected with wild-type influenza and pneumonia decreases to zero.

Similarly, from Figure 5, it can be observed that with preventive efforts as control strategies, the number of individuals coinfected with resistant influenza strain and pneumonia decreases drastically right from the beginning to zero by day 30.

The prevention measures help to reduce the transmission of the coinfection. Comparing Figures 2 and 3 with Figures 4 and 5, it can be observed that the preventive control strategies are more effective in curbing the spread of the coinfection compared to treatment control strategies.

4.1.3. Control with Prevention and Treatment of Influenza

Simulations are carried out to investigate the effect of implementing control strategies involving the prevention and treatment of influenza. Figures 6 and 7 show the results.

Figure 6: Individuals coinfected with wild-type influenza and pneumonia.

Figure 7: Individuals coinfected with resistant influenza and pneumonia.

It can be observed from Figure 6 that with the prevention and treatment of influenza as control strategies, the number of individuals coinfected with wild-type influenza and pneumonia decreases right from the beginning and it is at zero by about day 25.

Figure 7 shows that with the implementation of these control strategies, the number of individuals coinfected with resistant influenza and pneumonia significantly decreases and it is inconsequential at day 30 implying that it is unlikely for the coinfection to persist in the population.

The prevention and treatment of influenza as control strategies aid in reducing the transmission and in treatment of those who are already infected; however, the treatment poses a danger of development of drug resistance. Therefore, caution should be taken during drug administration.

4.1.4. Control with All Strategies

When all the control strategies are applied, the number of infected individuals decreases as shown in Figures 8 and 9. From Figure 8, it can be observed that with all the control strategies in place, the number of individuals coinfected with wild-type influenza and pneumonia decreases to zero by day 30.

Figure 8: Individuals coinfected with wild-type influenza and pneumonia.

Figure 9: Individuals coinfected with resistant influenza and pneumonia.

It can also be observed from Figure 9 that when all the control strategies are applied, the number of individuals coinfected with resistant influenza and pneumonia drastically decreases, and by about day 28, the number is already at zero.

5. Conclusion

As observed from Figures 2–9, different intervention mechanisms produce different results. Of great importance to public health is the control strategy that will help eradicate the diseases. Given that availability of resources is always a factor to consider when implementing control strategies, it is paramount to have a strategy with maximum benefit. From the results of this study, with the implementation of prevention measures only, which include and not limited to social distancing, vaccination, curbing mutation and reassortment, and curbing interspecies movement of the influenza virus, influenza pneumonia coinfection can be eradicated from a given population. Compared to the other control strategies investigated and discussed in Section 4.1, these preventive control strategies are more effective in curbing the spread of influenza pneumonia coinfection; hence, the public health sector and other stake holders could apply them to eradicate the coinfection within a given population.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Supplementary Materials

The supplementary material is the source code (MATLAB) that was used to do the numerical simulations. (Supplementary Materials)

References

K. R. Short, M. N. Habets, P. W. Hermans, and D. A. Diavatopoulos, “Interactions between streptococcus pneumoniae and influenza virus: a mutually beneficial relationship?” Future Microbiology, vol. 7, no. 5, pp. 609–624, 2012.

J. M. Rudd, H. K. Ashar, V. T. Chow, and N. Teluguakula, “Lethal synergism between influenza and streptococcus pneumoniae,” Journal of Infectious Pulmonary Diseases, vol. 2, no. 2, 2016.

A. M. Smith, F. R. Adler, R. M. Ribeiro et al., “Kinetics of coinfection with influenza a virus and streptococcus pneumoniae,” PLoS Pathogens, vol. 9, no. 3, Article ID e1003238, 2013.

L. S. Cauley and A. T. Vella, “Why is co-infection with influenza virus and bacteria so difficult to control?” Discovery Medicine, vol. 19, no. 102, p. 33, 2015.

WHO, “WHO—influenza,” 2003, http://www.who.int/mediacentre/factsheets/2003/fs211/en/.

CDC, “Seasonal flu death estimate increases worldwide,” 2017, https://www.cdc.gov/media/releases/2017/p1213-flu-death-estimate.html.

K. Thomas, “Global burden of pneumonia,” International Journal of Infectious Diseases, vol. 45, no. 1, 2016.

UNICEF, “Statistics by topic/pneumonia,” 2018, https://data.unicef.org/topic/child-health/pneumonia/.

D. E. Morris, D. W. Cleary, and S. C. Clarke, “Secondary bacterial infections associated with influenza pandemics,” Frontiers in Microbiology, vol. 8, p. 1041, 2017.

R. K. Gupta, R. George, and J. S. Nguyen-Van-Tam, “Bacterial pneumonia and pandemic influenza planning,” Emerging Infectious Diseases, vol. 14, no. 8, pp. 1187–1192, 2008.

J. A. McCullers, “Insights into the interaction between influenza virus and pneumococcus,” Clinical Microbiology Reviews, vol. 19, no. 3, pp. 571–582, 2006.

X.-X. Li and X.-N. Zhou, “Co-infection of tuberculosis and parasitic diseases in humans: a systematic review,” Parasites & vectors, vol. 6, no. 1, p. 79, 2013.

E. Y. Klein, B. Monteforte, A. Gupta et al., “The frequency of influenza and bacterial coinfection: a systematic review and meta-analysis,” Influenza and Other Respiratory Viruses, vol. 10, no. 5, pp. 394–403, 2016.

A. Arora, U. Krishnaswamy, R. Moideen, and M. Padmaja, “Tubercular and bacterial coinfection: a case series,” Lung India: Official Organ of Indian Chest Society, vol. 32, no. 2, p. 172, 2015.

D. M. Morens, J. K. Taubenberger, and A. S. Fauci, “Predominant role of bacterial pneumonia as a cause of death in pandemic influenza: implications for pandemic influenza preparedness,” The Journal of Infectious Diseases, vol. 198, no. 7, pp. 962–970, 2008.

D. M. Weinberger, L. Simonsen, R. Jordan, C. Steiner, M. Miller, and C. Viboud, “Impact of the 2009 influenza pandemic on pneumococcal pneumonia hospitalizations in the United States,” The Journal of Infectious Diseases, vol. 205, no. 3, pp. 458–465, 2011.

K. Grabowska, L. Högberg, P. Penttinen, Å. Svensson, and K. Ekdahl, “Occurrence of invasive pneumococcal disease and number of excess cases due to influenza,” BMC Infectious Diseases, vol. 6, no. 1, p. 58, 2006.

I. C. Michelow, K. Olsen, J. Lozano et al., “Epidemiology and clinical characteristics of community-acquired pneumonia in hospitalized children,” Pediatrics, vol. 113, no. 4, pp. 701–707, 2004.

S. Shrestha, B. Foxman, J. Berus et al., “The role of influenza in the epidemiology of pneumonia,” Scientific Reports, vol. 5, p. 15314, 2015.

D. S. Chertow and M. J. Memoli, “Bacterial coinfection in influenza,” Jama, vol. 309, no. 3, pp. 275–282, 2013.

L. A. McNamee and A. G. Harmsen, “Both influenza-induced neutrophil dysfunction and neutrophil-independent mechanisms contribute to increased susceptibility to a secondary streptococcus pneumoniae infection,” Infection and Immunity, vol. 74, no. 12, pp. 6707–6721, 2006.

K. F. van der Sluijs, L. J. R. van Elden, M. Nijhuis et al., “Il-10 is an important mediator of the enhanced susceptibility to pneumococcal pneumonia after influenza infection,” The Journal of Immunology, vol. 172, no. 12, pp. 7603–7609, 2004.

S. Shrestha, B. Foxman, S. Dawid et al., “Time and dose-dependent risk of pneumococcal pneumonia following influenza: a model for within-host interaction between influenza and streptococcus pneumoniae,” Journal of The Royal Society Interface, vol. 10, no. 86, p. 20130233, 2013.

CDC, “about antimicrobial resistance—antibiotic/antimicrobial resistance,” 2017, https://www.cdc.gov/drugresistance/about.html.

V. Singh, ““Antimicrobial resistance,” microbial Pathogens and Strategies for combating them,” Science, Technology and Education, vol. 1, pp. 291–296, 2013.

Lumen, “Drug resistance. lumen-Microbiology,” 2017, https://courses.lumenlearning.com/microbiology/chapter/drug-resistance/.

R. G. Webster and E. A. Govorkova, “Continuing challenges in influenza,” Annals of the New York Academy of Sciences, vol. 1323, no. 1, pp. 115–139, 2014.

M. Hussain, H. D. Galvin, T. Y. Haw, A. N. Nutsford, and M. Husain, “Drug resistance in influenza a virus: the epidemiology and management,” Infection and Drug Resistance, vol. 10, p. 121, 2017.

A. S. Monto, “Implications of antiviral resistance of influenza viruses,” Clinical Infectious Diseases, vol. 48, no. 4, pp. 397–399, 2009.

J. Lee, J. Kim, and H.-D. Kwon, “Optimal control of an influenza model with seasonal forcing and age-dependent transmission rates,” Journal of Theoretical Biology, vol. 317, pp. 310–320, 2013.

M. Jaberi-Douraki and S. M. Moghadas, “Optimal control of vaccination dynamics during an influenza epidemic,” Mathematical Biosciences and Engineering, vol. 11, no. 5, pp. 1045–1063, 2014.

A. S. Kadi, M. S. P. Hiremath, and S. R. Gani, “Optimal control of an influenza model with multiple control strategies,” International Journal on Recent Trends in Life Science & Mathematics, vol. 3, no. 12, pp. 1–7, 2016.

M. Elhia, O. Balatif, J. Bouyaghroumni, E. Labriji, and M. Rachik, “Optimal control applied to the spread of influenza a (h1n1),” Applied Mathematical Sciences, vol. 6, no. 82, pp. 4057–4065, 2012.

L. Wang, Z. Liu, D. Xu, and X. Zhang, “Global dynamics and optimal control of an influenza model with vaccination, media coverage and treatment,” International Journal of Biomathematics, vol. 10, no. 5, p. 1750068, 2017.

G. T. Tilahun, O. D. Makinde, and D. Malonza, “Modelling and optimal control of pneumonia disease with cost-effective strategies,” Journal of Biological Dynamics, vol. 11, no. sup2, pp. 400–426, 2017.

L. Pontryagin, V. Boltyanskii, R. V. Gamkrelidze, and E. F. Mischenko, The Mathematical Theory of Optimal Processes, Interscience, Saint-Nom-la-Bretèche, France, 1962.

D. L. Lukes, Differential Equations: Classical to Controlled, Elsevier, Amsterdam, Netherlands, 1982.

W. H. Fleming and R. W. Rishel, Deterministic and Stochastic Optimal Control, vol. 1, Springer Science & Business Media, Berlin, Germany, 2012.

G. R. Rose, “Numerical methods for solving optimal control problems,” University of Tennessee, Knoxville, TN, USA, 2015, M.Sc. thesis.

O. Balatif, M. Rachik, and J. Bouyaghroumni, “Application of optimal control theory to an seir model with immigration of infectives,” International Journal of Computer Science Issues, vol. 10, no. 2, pp. 230–236, 2013.

Y.-W. Chien, B. R. Levin, and K. P. Klugman, “The anticipated severity of a 1918-like influenza pandemic in contemporary populations: the contribution of antibacterial interventions,” PLoS One, vol. 7, no. 1, Article ID e29219, pp. 1–10, 2012.

CDC, “Clinical signs and symptoms of influenza–health professionals,” 2016, https://www.cdc.gov/flu/professionals/acip/clinical.htm.

B. M. Davis, A. E. Aiello, S. Dawid, P. Rohani, S. Shrestha, and B. Foxman, “Influenza and community-acquired pneumonia interactions: the impact of order and time of infection on population patterns,” American Journal of Epidemiology, vol. 175, no. 5, pp. 363–367, 2012.

N. D. Lambert, I. G. Ovsyannikova, V. S. Pankratz, R. M. Jacobson, and G. A. Poland, “Understanding the immune response to seasonal influenza vaccination in older adults: a systems biology approach,” Expert Review of Vaccines, vol. 11, no. 8, pp. 985–994, 2012.

M. Imran, T. Malik, A. R. Ansari, and A. Khan, “Mathematical analysis of swine influenza epidemic model with optimal control,” Japan Journal of Industrial and Applied Mathematics, vol. 33, no. 1, pp. 269–296, 2016.

WHO, Global Health Observatory (Gho) Data–Life Expectancy, WHO, Geneva, Switzerland, 2016.

W. Clara, J. Armero, D. Rodriguez et al., “Estimated incidence of influenza-virus-associated severe pneumonia in children in El Salvador, 2008–2010,” Bulletin of the World Health Organization, vol. 90, pp. 756–763, 2012.

M. L. Metersky, R. G. Masterton, H. Lode, T. M. File Jr, and T. Babinchak, “Epidemiology, microbiology, and treatment considerations for bacterial pneumonia complicating influenza,” International Journal of Infectious Diseases, vol. 16, no. 5, pp. e321–e331, 2012.

K. Cawcutt and A. C. Kalil, “Pneumonia with bacterial and viral coinfection,” Current Opinion in Critical Care, vol. 23, no. 5, pp. 385–390, 2017.

A. Moscona, “Neuraminidase inhibitors for influenza,” New England Journal of Medicine, vol. 353, no. 13, pp. 1363–1373, 2005.

J. A. McCullers, “Preventing and treating secondary bacterial infections with antiviral agents,” Antiviral Therapy, vol. 16, no. 2, pp. 123–135, 2011.