Dynamical behavior of an HIV infection model with treatment failure, adherence, and drug resistance

  • Maria Batool orcid

    Institute of Numerical Sciences, Kohat University of Science and Technology, Kohat 26000, Pakistan
  • Israr Ali Khan orcid

    Institute of Numerical Sciences, Kohat University of Science and Technology, Kohat 26000, Pakistan
  • Rashid Jan orcid

    Department of Mathematics, College of Science, Qassim University, Buraydah 51452, Saudi Arabia
  • Muhammad Shoaib Arif orcid

    Department of Mathematics and Sciences, College of Sciences and Humanities, Prince Sultan University, Riyadh 11586, Saudi Arabia
  • Ateeq Ur Rehman

    Department of Mathematics and Sciences, College of Sciences and Humanities, Prince Sultan University, Riyadh 11586, Saudi Arabia
Article ID: 4725
DOI: https://doi.org/10.59400/adecp4725
Keywords: HIV dynamics; nonlinear systems; antiretroviral therapy adherence; treatment failure; acquired drug resistance; stability analysis; optimal control

Abstract

Imperfect adherence to antiretroviral therapy (ART), virological failure, and the emergence of drug-resistant human immunodeficiency virus (HIV) strains continue to limit the long-term effectiveness of HIV control programmes. Existing deterministic models generally treat these mechanisms separately, so the pathway linking imperfect adherence to acquired resistance is seldom represented explicitly. We formulate a seven-compartment HIV transmission model in which stage-structured progression, ART initiation, adherence-mediated treatment failure, and acquired drug resistance appear as sequential, mechanistically linked transitions within a single tractable system. The model is shown to be well-posed: solutions remain non-negative and bounded, and a biologically feasible region is positively invariant. The basic reproduction number ℛ0 is obtained in closed form by the next-generation matrix method and separates into progression-driven and acquired immunodeficiency syndrome (AIDS)-mediated contributions, which makes the weight of each transmission pathway explicit. The disease-free equilibrium is locally and globally asymptotically stable when ℛ0 <  1, whereas a unique endemic equilibrium exists whenever ℛ0 > 1. Local asymptotic stability of the endemic equilibrium is conditional on the associated Routh–Hurwitz criteria and is verified numerically for the baseline parameter set through direct computation of the Jacobian eigenvalues. Normalised sensitivity indices rank the transmission rate first, followed by parameters governing the treatment-failure and drug-resistant compartments. An optimal control problem combining prevention, treatment uptake, and adherence support is characterised through Pontryagin’s Maximum Principle and solved by a forward–backward sweep; the combined strategy lowers the infected burden relative to the uncontrolled case over the horizon considered.

Published
2026-09-04
How to Cite
Batool, M., Khan, I. A., Jan, R., Arif, M. S., & Rehman, A. U. (2026). Dynamical behavior of an HIV infection model with treatment failure, adherence, and drug resistance. Advances in Differential Equations and Control Processes, 33(3). https://doi.org/10.59400/adecp4725
Issue
Vol. 33 No. 3 (2026): In Progress
Section
Article

Copyright (c) 2026 Maria Batool, Israr Ali Khan, Rashid Jan, Muhammad Shoaib Arif, Ateeq Ur Rehman

Creative Commons License

This work is licensed under a Creative Commons Attribution 4.0 International License.

References

[1]UNAIDS. Global HIV & AIDS Statistics - Fact Sheet 2025. UNAIDS; 2025. Available online:

[2]https://www.unaids.org/en/resources/fact-sheet

[3]World Health Organization. HIV Data and Statistics. World Health Organization; 2025. Available online:

[4]https://www.who.int/teams/global-hiv-hepatitis-and-stis-programmes/hiv/strategic-information/hiv-data-and-statistics

[5]Granich R, Gupta S, Suthar AB, et al. Antiretroviral Therapy in Prevention of HIV and TB: Update on Current Research

[6]Efforts. Current HIV Research. 2011; 9(6): 446–469.

[7]Naik PA, Yavuz M, Zu J. The role of prostitution on HIV transmission with memory: A modeling approach. Alexandria

[8]Engineering Journal. 2020; 59(4): 2513–2531.

[9]Nowak MA, May RM. Virus Dynamics: Mathematical Principles of Immunology and Virology. Oxford University

[10]Press; 2000.

[11]Kirschner D, Webb GF. A model for treatment strategy in the chemotherapy of AIDS. Bulletin of Mathematical Biology.

[12]; 58(2): 367–390.

[13]Blower SM, Aschenbach AN, Gershengorn HB, et al. Predicting the unpredictable: Transmission of drug-resistant

[14]HIV. Nature Medicine. 2001; 7(9): 1016–1020.

[15]Perelson AS, Nelson PW. Mathematical Analysis of HIV-1 Dynamics in Vivo. SIAM Review. 1999; 41(1): 3–44.

[16]Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford University Press; 1991.

[17]Hethcote HW. The Mathematics of Infectious Diseases. SIAM Review. 2000; 42(4): 599–653.

[18]Jan R, Boulaaras S, Khan IA, et al. Nonlocal and nonsingular kernel-based fractional dynamics of infectious disease

[19]with waning immunity. Journal of Mathematics and Computer Science. 2024; 38(1): 98–110.

[20]Kongson J, Thaiprayoon C, Neamvonk A, et al. Investigation of fractal-fractional HIV infection by evaluating

[21]the drug therapy effect in the Atangana–Baleanu sense. Mathematical Biosciences and Engineering. 2022; 19(11): 10762–10808. doi: 10.3934/mbe.2022504

[22]Khan H, Gómez‐Aguilar JF, Alkhazzan A, et al. A fractional order HIV‐TB coinfection model with nonsingular

[23]Mittag‐Leffler Law. Mathematical Methods in the Applied Sciences. 2020; 43(6): 3786–3806. doi: 10.1002/mma.6155

[24]Haq IU, Ali A, Ullah A, et al. Mathematical analysis of fractional order co-abuse infection model using power law

[25]type kernel. Journal of Mathematics and Computer Science. 2025; 36: 352–370. doi: 10.22436/jmcs.036.04.01

[26]Naeem M, Kuang Y, Li J. Stability of linear delay differential systems via Laplace transformation. Mathematical

[27]Biosciences. 2026; 400: 109778. doi: 10.1016/j.mbs.2026.109778

[28]Liu X, Wang H, Hu Z, et al. Global stability of an HIV pathogenesis model with cure rate. Nonlinear Analysis: Real

[29]World Applications. 2011; 12(6): 2947–2961. doi: 10.1016/j.nonrwa.2011.04.016

[30]Huo HF, Feng LX. Global stability for an HIV/AIDS epidemic model with different latent stages and treatment. Applied

[31]Mathematical Modelling. 2013; 37(3): 1480–1489. doi: 10.1016/j.apm.2012.04.013

[32]Xue L, Sun Y, Ren X, et al. Modelling the transmission dynamics and optimal control strategies for HIV infection in

[33]China. Journal of Biological Dynamics. 2023; 17(1): 2174275.

[34]Abimbade SF, Chuma FM, Olaniyi S, et al. Lyapunov Stability and Optimization of a Fractional‐Order HIV/AIDS

[35]Model Coupling Vertical Transmission With Saturated Treatment. International Journal of Mathematics and

[36]Mathematical Sciences. 2026; 2026(1): 3059683.

[37]Kumar N, Kashif M, Chauhan TS, et al. Modeling and Control of HIV/AIDS Epidemics: A Stochastic and Optimal

[38]Control Perspective. Journal of Nonlinear Mathematical Physics. 2025; 32(1): 80.

[39]Elmojtaba IM. An optimal control analysis of HIV-Visceral Leishmaniasis co-infection model. Partial Differential

[40]Equations in Applied Mathematics. 2025; 14: 101216.

[41]Ahmad I, Ahmad I, Ali N. Impact of Awareness and Protective Measures on the Transmission Dynamics of HIV/AIDS:

[42]A Mathematical Analysis. International Journal of Applied and Computational Mathematics. 2025; 11(3): 102.

[43]Diekmann O, Heesterbeek JAP, Roberts MG. The construction of next-generation matrices for compartmental

[44]epidemic models. Journal of The Royal Society Interface. 2010; 7(47): 873–885.

[45]Van Den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental

[46]models of disease transmission. Mathematical Biosciences. 2002; 180(1–2): 29–48.

[47]Castillo-Chavez C, Song B. Dynamical Models of Tuberculosis and Their Applications. Mathematical Biosciences

[48]and Engineering. 2004; 1(2): 361–404.

[49]Castillo-Chavez C, Feng Z, Huang W. On the computation of ℛ0 and its role on global stability. In: Mathematical

[50]Approaches for Emerging and Reemerging Infectious Diseases: An Introduction. Springer; 2002. pp. 229–250.

[51]Omame A, Isah ME, Abbas M, et al. A fractional order model for Dual Variants of COVID-19 and HIV co-infection

[52]via Atangana-Baleanu derivative. Alexandria Engineering Journal. 2022; 61(12): 9715–9731.

[53]Teklu SW, Guya TT, Kotola BS, et al. Analyses of an age structure HIV/AIDS compartmental model with optimal

[54]control theory. Scientific Reports. 2025; 15(1): 5491.

[55]Chitnis N, Hyman JM, Cushing JM. Determining Important Parameters in the Spread of Malaria Through the

[56]Sensitivity Analysis of a Mathematical Model. Bulletin of Mathematical Biology. 2008; 70(5): 1272–1296.

[57]Pontryagin LS, Boltyanskii VG, Gamkrelidze RV, et al. The Mathematical Theory of Optimal Processes. Wiley; 1962.

[58]Lenhart S, Workman JT. Optimal Control Applied to Biological Models. CRC Press; 2007.

[59]Hussein AA, Mebrate B. A Fractional Order Model for HIV/AIDS With Treatment and Optimal Control Using

[60]Caputo Derivative. International Journal of Mathematics and Mathematical Sciences. 2025; 2025(1): 9342227. doi:

[61]1155/ijmm/9342227

[62]Peter OJ, Teklu SW, Balogun GB, et al. HIV/AIDS compartmental model analysis with drug resistance treatment,

[63]vertical transmission, and optimal control theory. Scientific Reports. 2025; 15(1): 44156.

[64]Marsudi M, Trisilowati T, Musafir RR. Bifurcation and Optimal Control Analysis of an HIV/AIDS Model with

[65]Saturated Incidence Rate. Mathematics. 2025; 13(13): 2149.

[66]Poonia A, Chakrabarty SP. Strategic Control of Drug-Resistant HIV: Multi-Strain Modeling with Diagnosis,

[67]Adherence, and Treatment Switching. Bulletin of Mathematical Biology. 2025; 87(12): 173.

[68]Brauer F. Compartmental models in epidemiology. In: Mathematical Epidemiology. Springer; 2008. pp. 19–79.

[69]Perko L. Differential Equations and Dynamical Systems, 3rd ed. Springer; 2001. doi: 10.1007/978-1-4613-0003-8

[70]Khalil HK. Nonlinear Systems, 3rd ed. Prentice Hall; 2002. ISBN: 978-0-13-067389-3.

[71]Shuai Z, Van Den Driessche P. Global Stability of Infectious Disease Models Using Lyapunov Functions. SIAM Journal

[72]on Applied Mathematics. 2013; 73(4): 1513–1532.

[73]Fleming WH, Rishel RW. Deterministic and Stochastic Optimal Control. Springer-Verlag; 1975.