Chronic kidney disease (CKD), affecting approximately 843 million individuals globally, substantially increases susceptibility to severe COVID-19 outcomes, while SARS-CoV-2 infection independently accelerates renal deterioration through cytokine storms and microvascular injury. Despite this clinically significant bidirectional interaction, rigorous mathematical characterization of their co-dynamics within a unified framework integrating stability theory, bifurcation analysis, and optimal control remains limited. We formulated a seven-compartment deterministic model incorporating CKD'S irreversibility, the enhanced blue susceptibility of CKD patients, and vaccine imperfection, calibrated against Indian epidemiological data. The basic reproduction number $ \mathcal{R}_0 $ was derived using the next-generation matrix method, and stability and bifurcation analysis were performed. Sensitivity analysis identified transmission and immunity waning as dominant drivers, while optimal control strategies significantly reduced co-infections and hospitalizations, demonstrating the effectiveness of coordinated intervention policies.
Citation: Mallela Ankamma Rao, Emad K Jaradat, Medisetty Padma Devi, Prasantha Bharathi Dhandapani, Carlos Martin-Barreiro, Mohannad Al-Hmoud. Mathematical modeling of COVID-19 and chronic kidney disease co-infection with vaccination and optimal control: a bifurcation and sensitivity analysis approach[J]. AIMS Mathematics, 2026, 11(6): 17239-17292. doi: 10.3934/math.2026707
Chronic kidney disease (CKD), affecting approximately 843 million individuals globally, substantially increases susceptibility to severe COVID-19 outcomes, while SARS-CoV-2 infection independently accelerates renal deterioration through cytokine storms and microvascular injury. Despite this clinically significant bidirectional interaction, rigorous mathematical characterization of their co-dynamics within a unified framework integrating stability theory, bifurcation analysis, and optimal control remains limited. We formulated a seven-compartment deterministic model incorporating CKD'S irreversibility, the enhanced blue susceptibility of CKD patients, and vaccine imperfection, calibrated against Indian epidemiological data. The basic reproduction number $ \mathcal{R}_0 $ was derived using the next-generation matrix method, and stability and bifurcation analysis were performed. Sensitivity analysis identified transmission and immunity waning as dominant drivers, while optimal control strategies significantly reduced co-infections and hospitalizations, demonstrating the effectiveness of coordinated intervention policies.
| [1] | World Health Organization, Tracking SARS-CoV-2 variants, 2021. Available from: https://www.who.int/activities/tracking-SARS-CoV-2-variants. |
| [2] |
W. T. Harvey, A. M. Carabelli, B. Jackson, R. K. Gupta, E. C. Thomson, E. M. Harrison, et al., SARS-CoV-2 variants, spike mutations and immune escape, Nat. Rev. Microbiol., 19 (2021), 409–424. https://doi.org/10.1038/s41579-021-00573-0 doi: 10.1038/s41579-021-00573-0
|
| [3] |
A. Telenti, E. B. Hodcroft, D. L. Robertson, The evolution and biology of SARS-CoV-2 variants, Cold Spring Harb. Perspect. Med., 12 (2022), a041390. https://doi.org/10.1101/cshperspect.a041390 doi: 10.1101/cshperspect.a041390
|
| [4] |
E. Callaway, Why a highly mutated coronavirus variant has scientists on alert, Nature, 620 (2023), 934. https://doi.org/10.1038/d41586-023-02656-9 doi: 10.1038/d41586-023-02656-9
|
| [5] |
P. B. Mark, L. K. Stafford, M. E. Grams, H. Aalruz, S. Abd ElHafeez, A. A. Abdelgalil, et al., Global, regional, and national burden of chronic kidney disease in adults, 1990–2023, and its attributable risk factors: a systematic analysis for the Global Burden of Disease Study 2023, Lancet, 406 (2025), 2461–2482. https://doi.org/10.1016/S0140-6736(25)01853-7 doi: 10.1016/S0140-6736(25)01853-7
|
| [6] |
B. Bikbov, C. A. Purcell, A. S. Levey, M. Smith, A. Abdoli, M. Abebe, et al., Global, regional, and national burden of chronic kidney disease, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017, Lancet, 395 (2020), 709–733. https://doi.org/10.1016/S0140-6736(20)30045-3 doi: 10.1016/S0140-6736(20)30045-3
|
| [7] |
G. Eknoyan, T. Hostetter, G. L. Bakris, L. Hebert, A. S. Levey, H. H. Parving, et al., Proteinuria and other markers of chronic kidney disease: a position statement of the national kidney foundation (NKF) and the national institute of diabetes and digestive and kidney diseases (NIDDK), Am. J. Kidney Dis., 42 (2003), 617–622. https://doi.org/10.1016/S0272-6386(03)00826-6 doi: 10.1016/S0272-6386(03)00826-6
|
| [8] |
R. Borg, N. Carlson, J. Søndergaard, F. Persson, The growing challenge of chronic kidney disease: an overview of current knowledge, Int. J. Nephrol., 2023 (2023), 9609266. https://doi.org/10.1155/2023/9609266 doi: 10.1155/2023/9609266
|
| [9] |
P. E. Stevens, S. B. Ahmed, J. J. Carrero, B. Foster, A. Francis, R. K. Hall, et al., KDIGO 2024 clinical practice guideline for the evaluation and management of chronic kidney disease, Kidney Int., 105 (2024), S117–S314. https://doi.org/10.1016/j.kint.2023.10.018 doi: 10.1016/j.kint.2023.10.018
|
| [10] |
C. W. Mende, Chronic kidney disease and SGLT2 inhibitors: a review of the evolving treatment landscape, Adv. Ther., 39 (2024), 148–164. https://doi.org/10.1007/s12325-021-01994-2 doi: 10.1007/s12325-021-01994-2
|
| [11] |
S. S. Jdiaa, R. Mansour, A. El Alayli, A. Gautam, P. Thomas, R. A. Mustafa, COVID-19 and chronic kidney disease: an updated overview of reviews, J. Nephrol., 35 (2022), 69–85. https://doi.org/10.1007/s40620-021-01206-8 doi: 10.1007/s40620-021-01206-8
|
| [12] | E. Y. M. Chung, S. C. Palmer, P. Natale, A. Krishnan, T. E. Cooper, V. M. Saglimbene, et al., Incidence and outcomes of COVID-19 in people with CKD: a systematic review and meta-analysis, Am. J. Kidney Dis., 78 (2021), 804–815. https://doi.org/10.1053/j.ajkd.2021.07.003 |
| [13] |
V. Ambalarajan, A. R. Mallela, V. Sivakumar, P. B. Dhandapani, V. Leiva, C. Martin-Barreiro, et al., A six-compartment model for COVID-19 with transmission dynamics and public health strategies, Sci. Rep., 14 (2024), 22226. https://doi.org/10.1038/s41598-024-72487-9 doi: 10.1038/s41598-024-72487-9
|
| [14] |
Y. Zhang, Y. Zhao, J. Wang, X. Zheng, D. Xu, J. Lv, et al., Long-term renal outcomes of patients with COVID-19: a meta-analysis of observational studies, J. Nephrol., 36 (2023), 2441–2456. https://doi.org/10.1007/s40620-023-01731-8 doi: 10.1007/s40620-023-01731-8
|
| [15] |
I. Chen, L. Chang, C. Ho, J. Wu, Y. Tsai, C. Lin, et al., Association between COVID-19 and the development of chronic kidney disease in patients without initial acute kidney injury, Sci. Rep., 15 (2025), 10924. https://doi.org/10.1038/s41598-025-96032-4 doi: 10.1038/s41598-025-96032-4
|
| [16] |
H. Rabb, Kidney diseases in the time of COVID-19: major challenges to patient care, J. Clin. Invest., 130 (2020), 2749–2751. https://doi.org/10.1172/JCI138871 doi: 10.1172/JCI138871
|
| [17] |
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
|
| [18] |
O. Diekmann, J. A. P. Heesterbeek, M. G. Roberts, The construction of next-generation matrices for compartmental epidemic models, J. R. Soc. Interface, 7 (2010), 873–885. https://doi.org/10.1098/rsif.2009.0386 doi: 10.1098/rsif.2009.0386
|
| [19] | Z. Liu, G. Wan, B. A. Prakash, M. S. Y. Lau, W. Jin, A review of graph neural networks in epidemic modeling, Proceedings of the 30th ACM SIGKDD Conference on Knowledge Discovery and Data Mining, 2024, 6577–6587. https://doi.org/10.1145/3637528.3671455 |
| [20] |
I. Takaidza, O. D. Makinde, O. K. Okosun, Computational modelling and optimal control of Ebola virus disease with non-linear incidence rate, J. Phys.: Conf. Ser., 818 (2017), 012003. https://doi.org/10.1088/1742-6596/818/1/012003 doi: 10.1088/1742-6596/818/1/012003
|
| [21] |
O. Kiseleva, S. Yakovlev, O. Prytomanova, O. Kuzenkov, Mathematical modeling of regional infectious disease dynamics based on extended compartmental models, Computation, 13 (2025), 187. https://doi.org/10.3390/computation13080187 doi: 10.3390/computation13080187
|
| [22] |
S. Marino, I. B. Hogue, C. J. Ray, D. E. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254 (2008), 178–196. https://doi.org/10.1016/j.jtbi.2008.04.011 doi: 10.1016/j.jtbi.2008.04.011
|
| [23] |
D. S. Mgonja, A. Hugo, A. Hassan, O. D. Makinde, Mathematical modeling and analysis of COVID-19 and typhoid fever co-dynamics with treatment, Sci. Rep., 15 (2025), 33206. https://doi.org/10.1038/s41598-024-82955-x doi: 10.1038/s41598-024-82955-x
|
| [24] |
J. O. Akanni, S. Ajao, S. F. Abimbade, Fatmawati, Dynamical analysis of COVID-19 and tuberculosis co-infection using mathematical modelling approach, Math. Model. Control, 4 (2024), 208–229. https://doi.org/10.3934/mmc.2024018 doi: 10.3934/mmc.2024018
|
| [25] |
K. G. Mekonen, L. L. Obsu, Mathematical modeling and analysis for the co-infection of COVID-19 and tuberculosis, Heliyon, 8 (2022), e11195. https://doi.org/10.1016/j.heliyon.2022.e11195 doi: 10.1016/j.heliyon.2022.e11195
|
| [26] |
M. Ankamma Rao, E. K. Jaradat, M. Padma Devi, P. B. Dhandapani, R. M. Nalule, M. Al-Hmoud, Modeling COVID-19 pneumonia and COVID-associated pulmonary aspergillosis: sensitivity analysis and optimal control, BMC Infect. Dis., 25 (2025), 1192. https://doi.org/10.1186/s12879-025-11606-x doi: 10.1186/s12879-025-11606-x
|
| [27] |
M. A. Hye, M. H. A. Biswas, M. F. Uddin, M. M. Rahman, A mathematical model for the transmission of co-infection with COVID-19 and kidney disease, Sci. Rep., 14 (2024), 5680. https://doi.org/10.1038/s41598-024-56399-2 doi: 10.1038/s41598-024-56399-2
|
| [28] |
D. Gibertoni, C. Reno, P. Rucci, M. P. Fantini, A. Buscaroli, G. Mosconi, et al., COVID-19 incidence and mortality in non-dialysis chronic kidney disease patients, PLoS ONE, 16 (2021), e0254525. https://doi.org/10.1371/journal.pone.0254525 doi: 10.1371/journal.pone.0254525
|
| [29] |
C. Castillo-Chavez, B. Song, Dynamical models of tuberculosis and their applications, Math. Biosci. Eng., 1 (2004), 361–404. https://doi.org/10.3934/mbe.2004.1.361 doi: 10.3934/mbe.2004.1.361
|
| [30] |
N. Chitnis, J. M. Hyman, J. M. Cushing, Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model, Bull. Math. Biol., 70 (2008), 1272–1296. https://doi.org/10.1007/s11538-008-9299-0 doi: 10.1007/s11538-008-9299-0
|
| [31] | L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, E. F. Mishchenko, The mathematical theory of optimal processes, Boca Raton, FL: CRC Press, 2018. https://doi.org/10.1201/9780203749319 |
| [32] | S. Lenhart, J. T. Workman, Optimal control applied to biological models, Boca Raton, FL: Chapman & Hall/CRC, 2007. https://doi.org/10.1201/9781420011418 |
| [33] |
K. Jain, A. Bhattacharjee, S. Krishnamurhty, Mathematical analysis of COVID-19 and TB co-infection dynamics with optimal control, Model. Earth Syst. Environ., 11 (2025), 65. https://doi.org/10.1007/s40808-024-02197-8 doi: 10.1007/s40808-024-02197-8
|
| [34] |
S. Kurmi, U. Chouhan, A multicompartment mathematical model to study the dynamic behaviour of COVID-19 using vaccination as control parameter, Nonlinear Dyn., 109 (2022), 2185–2201. https://doi.org/10.1007/s11071-022-07591-4 doi: 10.1007/s11071-022-07591-4
|
| [35] |
S. Bugalia, J. P. Tripathi, H. Wang, Estimating the time-dependent effective reproduction number and vaccination rate for COVID-19 in the USA and India, Math. Biosci. Eng., 20 (2023), 4673–4689. https://doi.org/10.3934/mbe.2023216 doi: 10.3934/mbe.2023216
|
| [36] |
R. K. Rai, S. Khajanchi, P. K. Tiwari, E. Venturino, A. K. Misra, Impact of social media advertisements on the transmission dynamics of COVID-19 pandemic in India, J. Appl. Math. Comput., 68 (2022), 19–44. https://doi.org/10.1007/s12190-021-01507-y doi: 10.1007/s12190-021-01507-y
|
| [37] |
S. J. Schrauben, H. Chen, E. Lin, C. Jepson, W. Yang, J. J. Scialla, et al., Hospitalizations among adults with chronic kidney disease in the United States: a cohort study, PLOS Med., 17 (2020), e1003470. https://doi.org/10.1371/journal.pmed.1003470 doi: 10.1371/journal.pmed.1003470
|
| [38] |
Y. Gao, Y. Chen, M. Liu, S. Shi, J. Tian, Impacts of immunosuppression and immunodeficiency on COVID-19: a systematic review and meta-analysis, J. Infect., 81 (2020), e93–e95. https://doi.org/10.1016/j.jinf.2020.05.017 doi: 10.1016/j.jinf.2020.05.017
|
| [39] |
V. P. Bajiya, S. Bugalia, J. P. Tripathi, Mathematical modeling of COVID-19: impact of non-pharmaceutical interventions in India, Chaos, 30 (2020), 113143. https://doi.org/10.1063/5.0021353 doi: 10.1063/5.0021353
|
| [40] |
F. Shahbazi, A. Doosti-Irani, A. Soltanian, J. Poorolajal, Global forecasting of chronic kidney disease mortality rates and numbers with the generalized additive model, BMC Nephrol., 25 (2024), 286. https://doi.org/10.1186/s12882-024-03720-w doi: 10.1186/s12882-024-03720-w
|
| [41] |
Y. Cheng, R. Luo, K. Wang, M. Zhang, Z. Wang, L. Dong, et al., Kidney disease is associated with in-hospital death of patients with COVID-19, Kidney Int., 97 (2020), 829–838. https://doi.org/10.1016/j.kint.2020.03.005 doi: 10.1016/j.kint.2020.03.005
|
| [42] | COVID-19 India, COVID-19 India dataset, 2024. Available from: https://data.covid19india.org/. |
| [43] | Institute for Health Metrics and Evaluation (IHME), GBD Results Tool, 2025. Available from: https://vizhub.healthdata.org/gbd-results/. |
| [44] | D. R. Powell, J. Fair, R. J. LeClaire, L. M. Moore, D. Thompson, Sensitivity analysis of an infectious disease model, ISDC 2005, 2005. Available from: https://proceedings.systemdynamics.org/2005/proceed/papers/LECLA330.pdf. |
| [45] |
J. Wu, R. Dhingra, M. Gambhir, J. V. Remais, Sensitivity analysis of infectious disease models: methods, advances and their application, J. R. Soc. Interface, 10 (2013), 20121018. https://doi.org/10.1098/rsif.2012.1018 doi: 10.1098/rsif.2012.1018
|
| [46] | M. Martcheva, An introduction to mathematical epidemiology, New York, NY: Springer, 2015. https://doi.org/10.1007/978-1-4899-7612-3 |
| [47] |
Z. Zhang, R. Gul, A. Zeb, Global sensitivity analysis of COVID-19 mathematical model, Alex. Eng. J., 60 (2020), 565–572. https://doi.org/10.1016/j.aej.2020.09.035 doi: 10.1016/j.aej.2020.09.035
|
| [48] |
V. Ambalarajan, A. R. Mallela, P. B. Dhandapani, V. Sivakumar, V. Leiva, C. Castro, Multi-strain COVID-19 dynamics with vaccination strategies: mathematical modeling and case study, Alex. Eng. J., 119 (2025), 665–684. https://doi.org/10.1016/j.aej.2025.01.105 doi: 10.1016/j.aej.2025.01.105
|
| [49] | P. Hartman, Ordinary differential equations (2nd ed.), Society for Industrial and Applied Mathematics, 2002. https://doi.org/10.1137/1.9780898719222 |
| [50] | J. P. LaSalle, The stability of dynamical systems, Philadelphia, PA: Society for Industrial and Applied Mathematics, 1976. |