Modeling optimal treatment strategies in a heterogeneous mixing model

BackgroundMany mathematical models assume random or homogeneous mixing for various infectious diseases. Homogeneous mixing can be generalized to mathematical models with multi-patches or age structure by incorporating contact matrices to capture the dynamics of the heterogeneously mixing populations. Contact or mixing patterns are difficult to measure in many infectious diseases including influenza. Mixing patterns are considered to be one of the critical factors for infectious disease modeling.MethodsA two-group influenza model is considered to evaluate the impact of heterogeneous mixing on the influenza transmission dynamics. Heterogeneous mixing between two groups with two different activity levels includes proportionate mixing, preferred mixing and like-with-like mixing. Furthermore, the optimal control problem is formulated in this two-group influenza model to identify the group-specific optimal treatment strategies at a minimal cost. We investigate group-specific optimal treatment strategies under various mixing scenarios.ResultsThe characteristics of the two-group influenza dynamics have been investigated in terms of the basic reproduction number and the final epidemic size under various mixing scenarios. As the mixing patterns become proportionate mixing, the basic reproduction number becomes smaller; however, the final epidemic size becomes larger. This is due to the fact that the number of infected people increases only slightly in the higher activity level group, while the number of infected people increases more significantly in the lower activity level group. Our results indicate that more intensive treatment of both groups at the early stage is the most effective treatment regardless of the mixing scenario. However, proportionate mixing requires more treated cases for all combinations of different group activity levels and group population sizes.ConclusionsMixing patterns can play a critical role in the effectiveness of optimal treatments. As the mixing becomes more like-with-like mixing, treating the higher activity group in the population is almost as effective as treating the entire populations since it reduces the number of disease cases effectively but only requires similar treatments. The gain becomes more pronounced as the basic reproduction number increases. This can be a critical issue which must be considered for future pandemic influenza interventions, especially when there are limited resources available.

[1]  Frank Diederich,et al.  Mathematical Epidemiology Of Infectious Diseases Model Building Analysis And Interpretation , 2016 .

[2]  R. Antia,et al.  How do pathogen evolution and host heterogeneity interact in disease emergence? , 2006, Proceedings of the Royal Society B: Biological Sciences.

[3]  M. Kretzschmar,et al.  Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. , 2006, American journal of epidemiology.

[4]  Julien Arino,et al.  A model for influenza with vaccination and antiviral treatment. , 2008, Journal of theoretical biology.

[5]  A. Nizam,et al.  Containing pandemic influenza with antiviral agents. , 2004, American journal of epidemiology.

[6]  H. Kwon,et al.  Optimal control of an influenza model with seasonal forcing and age-dependent transmission rates. , 2013, Journal of theoretical biology.

[7]  A. Nold Heterogeneity in disease-transmission modeling , 1980 .

[8]  J. Watmough,et al.  Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. , 2002, Mathematical biosciences.

[9]  Herbert W. Hethcote,et al.  Mixing patterns between age groups in social networks , 2007, Soc. Networks.

[10]  Benjamin J. Cowling,et al.  Assortativity and the Probability of Epidemic Extinction: A Case Study of Pandemic Influenza A (H1N1-2009) , 2010, Interdisciplinary perspectives on infectious diseases.

[11]  M. Safan,et al.  Title Pros and cons of estimating the reproduction number from earlyepidemic growth rate of influenza A ( H 1 N 1 ) 2009 , 2010 .

[12]  John T. Workman,et al.  Optimal Control Applied to Biological Models , 2007 .

[13]  C. Castillo-Chavez,et al.  A note on the use of influenza vaccination strategies when supply is limited. , 2011, Mathematical biosciences and engineering : MBE.

[14]  M. L. Chambers The Mathematical Theory of Optimal Processes , 1965 .

[15]  L. S. Pontryagin,et al.  Mathematical Theory of Optimal Processes , 1962 .

[16]  Christopher A. Gilligan,et al.  Optimal control of epidemics in metapopulations , 2009, Journal of The Royal Society Interface.

[17]  Herbert W. Hethcote,et al.  Epidemiological models for heterogeneous populations: proportionate mixing, parameter estimation, and immunization programs , 1987 .

[18]  L. Meyers,et al.  When individual behaviour matters: homogeneous and network models in epidemiology , 2007, Journal of The Royal Society Interface.

[19]  J. Janssen,et al.  Deterministic and Stochastic Optimal Control , 2013 .

[20]  O Diekmann,et al.  The construction of next-generation matrices for compartmental epidemic models , 2010, Journal of The Royal Society Interface.

[21]  G. Chowell,et al.  Optimal control for pandemic influenza: the role of limited antiviral treatment and isolation. , 2010, Journal of theoretical biology.

[22]  A. Nizam,et al.  Containing Pandemic Influenza at the Source , 2005, Science.

[23]  P. Kaye Infectious diseases of humans: Dynamics and control , 1993 .

[24]  Julien Arino,et al.  A final size relation for epidemic models. , 2007, Mathematical biosciences and engineering : MBE.

[25]  R. Mikolajczyk,et al.  Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases , 2008, PLoS medicine.

[26]  Jia Li,et al.  Behavior Changes in SIS STD Models with Selective Mixing , 1997, SIAM J. Appl. Math..

[27]  Steve Leach,et al.  Potential Impact of Antiviral Drug Use during Influenza Pandemic , 2005, Emerging infectious diseases.

[28]  Lisa Sattenspiel,et al.  Modeling and analyzing HIV transmission: the effect of contact patterns , 1988 .

[29]  Joel C. Miller,et al.  Supplementary Text S1 , 2014 .

[30]  C. Castillo-Chavez,et al.  A note on the use of optimal control on a discrete time model of influenza dynamics. , 2011, Mathematical biosciences and engineering : MBE.

[31]  O. Diekmann,et al.  Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation , 2000 .

[32]  Daniel Hlubinka,et al.  Kermack-McKendrick epidemic model revisited , 2007, Kybernetika.

[33]  M. Safan,et al.  Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009 , 2010, Theoretical Biology and Medical Modelling.

[34]  V. Colizza,et al.  Metapopulation epidemic models with heterogeneous mixing and travel behaviour , 2014, Theoretical Biology and Medical Modelling.

[35]  Fred Brauer,et al.  Epidemic Models with Heterogeneous Mixing and Treatment , 2008, Bulletin of mathematical biology.

[36]  Karan P. Singh,et al.  Theoretical Biology and Medical Modelling , 2007 .

[37]  Erik M. Volz,et al.  Correction: Effects of Heterogeneous and Clustered Contact Patterns on Infectious Disease Dynamics , 2011, PLoS Computational Biology.

[38]  F. Brauer Heterogeneous mixing in epidemic models , 2012 .

[39]  D. Cummings,et al.  Strategies for containing an emerging influenza pandemic in Southeast Asia , 2005, Nature.

[40]  G. Chowell,et al.  Modeling Optimal Age-Specific Vaccination Strategies Against Pandemic Influenza , 2011, Bulletin of Mathematical Biology.

[41]  F. Brauer,et al.  Simple models for containment of a pandemic , 2006, Journal of The Royal Society Interface.

[42]  D. Earn,et al.  Generality of the Final Size Formula for an Epidemic of a Newly Invading Infectious Disease , 2006, Bulletin of mathematical biology.