We propose a fractional order model for human papillomavirus (HPV) dynamics, including the effects of vaccination and public health education on developing cervical cancer. First, we discuss the general structure of Caputo fractional derivatives and integrals. Next, we define the fractional HPV model using Caputo derivatives. The model equilibrium quantities, with their stability, are discussed based on the magnitude of the reproduction number. We compute and simulate numerical solutions of the presented fractional model using the Adams-Bashforth-Moulton scheme. Meanwhile, real data sourced from reports from the World Health Organization is used to establish the parameters and compute the basic reproduction number. We present figures of state variables for different fractional orders and the classical integer order. The impacts of vaccination and public health education are discussed through numerical simulations. From the results, we observe that an increase in both vaccination rates and public health education increases the quality of life, and thus, reduces disease burden and suffering in communities. The results also confirm that modeling HPV transmission dynamics using fractional derivatives includes history effects in the model, making the model further insightful and appropriate for studying HPV dynamics.
Citation: Simphiwe M. Simelane, Justin B. Munyakazi, Phumlani G. Dlamini, Oluwaseun F. Egbelowo. Projections of human papillomavirus vaccination and its impact on cervical cancer using the Caputo fractional derivative[J]. Mathematical Biosciences and Engineering, 2023, 20(7): 11605-11626. doi: 10.3934/mbe.2023515
We propose a fractional order model for human papillomavirus (HPV) dynamics, including the effects of vaccination and public health education on developing cervical cancer. First, we discuss the general structure of Caputo fractional derivatives and integrals. Next, we define the fractional HPV model using Caputo derivatives. The model equilibrium quantities, with their stability, are discussed based on the magnitude of the reproduction number. We compute and simulate numerical solutions of the presented fractional model using the Adams-Bashforth-Moulton scheme. Meanwhile, real data sourced from reports from the World Health Organization is used to establish the parameters and compute the basic reproduction number. We present figures of state variables for different fractional orders and the classical integer order. The impacts of vaccination and public health education are discussed through numerical simulations. From the results, we observe that an increase in both vaccination rates and public health education increases the quality of life, and thus, reduces disease burden and suffering in communities. The results also confirm that modeling HPV transmission dynamics using fractional derivatives includes history effects in the model, making the model further insightful and appropriate for studying HPV dynamics.
[1] | CDC fact sheet. Centers for Disease Control and Prevention, Genital HPV infection. Available from: http://www.cdc.gov/std/HPV/STDFact-HPV.htm. |
[2] | S. L. Lee, A. M. Tameru, A mathematical model of Human Papillomavirus (HPV) in the United States and its impact on Cervical Cancer, J. Cancer, 3 (2012), 262–268. https://doi.org/10.7150/jca.4161 doi: 10.7150/jca.4161 |
[3] | National Cancer Institute, HPV and Cancer. Available from: https://www.cancer.gov/about-cancer/causes-prevention/risk/infectious-agents/hpv-and-cancer. |
[4] | World Health Organization, Cervical Cancer Fact Sheet, Access date: 02 June 2022. Available from: https://www.who.int/news-room/fact-sheets/detail/cervical-cancer. |
[5] | K. P. Braaten, M. R. Laufer, Human Papillomavirus (HPV), HPV-related disease, and the HPV vaccine, Rev. Obstet. Gynecol., 1 (2008), 2–10. |
[6] | Centers for Disease Control and Prevention, Human Papillomavirus (HPV) Vaccination: What Everyone Should Know, Accesse date: 02 June 2022. https://www.cdc.gov/vaccines/vpd/hpv/public/index.html. |
[7] | U. K. Nwajeri, A. B. Panle, A. Omame, M. C. Obi, C. P. Onyenegecha, On the fractional order model for HPV and Syphilis using non–singular kernel, Results Phys., 37 (2022), 105463. https://doi.org/10.1016/j.rinp.2022.105463 doi: 10.1016/j.rinp.2022.105463 |
[8] | I. Podlubny, Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and Some of their Applications, Academic Press, San Diego, CA, 1999. |
[9] | D. Baleanu, M. Hassan Abadi, A. Jajarmi, K. Zarghami Vahid, J. J. Nieto, A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects, Alexandria Eng. J., 61 (2022), 4779–4791. https://doi.org/10.1016/j.aej.2021.10.030 doi: 10.1016/j.aej.2021.10.030 |
[10] | P. A. Naik, J. Zu, K. M. Owolabi, Global dynamics of a fractional order model for the transmission of HIV epidemic with optimal control, Chaos, Solitons Fractals, 138 (2020), 109826. https://doi.org/10.1016/j.chaos.2020.109826 doi: 10.1016/j.chaos.2020.109826 |
[11] | A. Atangana, E. Alabaraoye, Solving a system of fractional partial differential equations arising in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equations, Adv. Differ. Equations, 94 (2013). https://doi.org/10.1186/1687-1847-2013-94 doi: 10.1186/1687-1847-2013-94 |
[12] | R. Scherer, S. L. Kalla, Y. Tang, J. Huang, The Grunwald-Letnikov method for fractional differential equations, Comput. Math. Appl., 62 (2011), 902–917. https://doi.org/10.1016/j.camwa.2011.03.054 doi: 10.1016/j.camwa.2011.03.054 |
[13] | U. K. Nwajeri, A. Omame, C. P. Onyenegecha, Analysis of a fractional order model for HPV and CT co-infection, Results Phys., 28 (2021), 104643. https://doi.org/10.1016/j.rinp.2021.104643 doi: 10.1016/j.rinp.2021.104643 |
[14] | K. M. Owolabi, A. Atangana, Numerical Methods for Fractional Differentiation, Springer Singapore, 2019. https://doi.org/10.1007/978-981-15-0098-5 |
[15] | U. K. Nwajeri, A. B. Panle, A. Omame, M. C. Obi, C. P. Onyenegecha, On the fractional order model for HPV and Syphilis using non–singular kernel, Results Phys., 37 (2022), 105463. https://doi.org/10.1016/j.rinp.2022.105463 doi: 10.1016/j.rinp.2022.105463 |
[16] | K. Diethelm, N. J. Ford, Analysis of fractional differential equations, J. Math. Anal. Appl., 265 (2002), 229–248. https://doi.org/10.1006/jmaa.2000.7194 doi: 10.1006/jmaa.2000.7194 |
[17] | S. M. Simelane, P. G. Dlamini, A fractional order differential equation model for Hepatitis B virus with saturated incidence, Results Phys., 24 (2021), 104114. https://doi.org/10.1016/j.rinp.2021.104114 doi: 10.1016/j.rinp.2021.104114 |
[18] | Z. M. Odibat, N. T. Shawagfeh, Generalized Taylor's formula, Appl. Math. Comput., 186 (2007), 286–293. https://doi.org/10.1016/j.amc.2006.07.102 doi: 10.1016/j.amc.2006.07.102 |
[19] | N. Ozalp, E. Demirci, A fractional order SEIR model with vertical transmission, Math. Comput. Model., 54 (2011), 1–6. https://doi.org/10.1016/j.mcm.2010.12.051 doi: 10.1016/j.mcm.2010.12.051 |
[20] | W. Lin, Global existence theory and chaos control of fractional differential equations, J. Math. Anal. Appl., 332 (2007), 709–726. https://doi.org/10.1016/j.jmaa.2006.10.040 doi: 10.1016/j.jmaa.2006.10.040 |
[21] | O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324 |
[22] | J. Jia, J. Xiao, Stability analysis of a disease resistance seirs model with nonlinear incidence rate, Adv. Differ. Equations, 75 (2018). https://doi.org/10.1186/s13662-018-1494-1 |
[23] | P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6 |
[24] | J. L. S. Allen, Introduction to Mathematical Biology, Pearson Education, Inc, Upper Saddle River, 2007. |
[25] | MATLAB. 9.7.0.1190202 (R2019b). Natick, Massachusetts: The MathWorks Inc, 2018. |
[26] | Wolfram Research, Inc, Mathematica, Version 9.0, Champaign, IL, 2012. |
[27] | Indian Council of Medical Research, Consensus Document for Management of Cancer Cervix, 2016. |
[28] | J. Ferlay, M. Ervik, F. Lam, M. Colombet, L. Mery, M. Pineros, et al., Global Cancer Observatory: Cancer Today. Lyon, France: International Agency for Research on Cancer, 2020. Available from: https://gco.iarc.fr/today. |
[29] | E. Amponsah-Dacosta, N. Blose, V. V. Nkwinika, V. Chepkurui, Human papillomavirus vaccination in South Africa: Programmatic challenges and opportunities for integration with other adolescent health services, Front. Public Health, 10 (2022). https://doi.org/10.3389/fpubh.2022.799984 doi: 10.3389/fpubh.2022.799984 |
[30] | United Nations, Department of Economic and Social Affairs, Population Division, World Population Prospects 2019, Online Edition. Rev. 1, Access date: 02 January 2021, Available from: https://population.un.org/wpp/Download/Standard/Population/. |
[31] | World Health Organization, Global strategy to accelerate the elimination of cervical cancer as a public health problem, Access date: 05 March 2021, Available from: https://www.who.int/publications/i/item/9789240014107. |
[32] | N. Gupta, A. S. Chauhan, S. Prinja, A. K. Pandey, Impact of COVID-19 on outcomes for patients with cervical cancer in India, JCO Global Oncol., 7 (2021), 716–725. https://doi.org/10.1200/GO.20.00654 doi: 10.1200/GO.20.00654 |
[33] | L. Bruni, G. Albero, B. Serrano, M. Mena, J. J. Collado, D. Gómez, et al., ICO/IARC Information Centre on HPV and Cancer (HPV Information Centre), Human Papillomavirus and Related Diseases in India, Summary Report, Access date: 18 June 2022. |