Human metapneumovirus (HMPV) is a pathogen that causes severe respiratory infections, especially in children, the elderly, and immunocompromised people. Deterministic epidemiological models cannot capture memory effects or stochastic variability of the environment and thus limit the formulation of effective control measures. In this paper, we developed a stochastic optimal control problem, based on a fractional-order stochastic dynamic, to promote vaccination and treatment solutions to the disease burden. In this study, we designed a modified fractional-order stochastic SEIR model with a Brownian motion to consider long-range time memory effects and stochastic environmental variation. We verified the existence, uniqueness, and positivity of the model, as well as identified disease-free equilibrium, the basic reproduction number $ R_0 $, and endemic equilibrium. Local and global stability of both states were examined, and a bifurcation analysis disclosed a forward transcritical bifurcation at the point $ R_0 = 1 $. Sensitivity analysis highlighted the transmission rate $\beta$ (+1.000), progression rate $\sigma$ (+0.4545), and recovery rate $\gamma$ (-0.6818) as critical parameters affecting $R_0$. We proposed an optimal control analysis that uncovered time-dependent optimal vaccination and treatment techniques. Under the optimal treatment-focused strategy, decreasing the fractional order $ \alpha $ from 1.0 to 0.7 reduced peak infections by $ 36.1\% $ (from 244 to 156 cases) and delayed the epidemic peak by $ 57.8\% $ (from 26.8 to 42.3 days). A strong linear correlation ($R^2 = 0.96$) disclosed that lower fractional orders $\alpha$ impose earlier but less intensive interventions. This adaptive control framework offers the public health authorities a tool for designing an outbreak response that is resilient to real-life uncertainty.
Citation: Md Hossain Islam, Md. Nur Alam, Noor Muhammad, Xinsong Yang, Md. Jakir Hossen. Mathematical modeling of human metapneumovirus transmission dynamics with optimal control analysis[J]. AIMS Mathematics, 2026, 11(7): 21570-21606. doi: 10.3934/math.2026874
Human metapneumovirus (HMPV) is a pathogen that causes severe respiratory infections, especially in children, the elderly, and immunocompromised people. Deterministic epidemiological models cannot capture memory effects or stochastic variability of the environment and thus limit the formulation of effective control measures. In this paper, we developed a stochastic optimal control problem, based on a fractional-order stochastic dynamic, to promote vaccination and treatment solutions to the disease burden. In this study, we designed a modified fractional-order stochastic SEIR model with a Brownian motion to consider long-range time memory effects and stochastic environmental variation. We verified the existence, uniqueness, and positivity of the model, as well as identified disease-free equilibrium, the basic reproduction number $ R_0 $, and endemic equilibrium. Local and global stability of both states were examined, and a bifurcation analysis disclosed a forward transcritical bifurcation at the point $ R_0 = 1 $. Sensitivity analysis highlighted the transmission rate $\beta$ (+1.000), progression rate $\sigma$ (+0.4545), and recovery rate $\gamma$ (-0.6818) as critical parameters affecting $R_0$. We proposed an optimal control analysis that uncovered time-dependent optimal vaccination and treatment techniques. Under the optimal treatment-focused strategy, decreasing the fractional order $ \alpha $ from 1.0 to 0.7 reduced peak infections by $ 36.1\% $ (from 244 to 156 cases) and delayed the epidemic peak by $ 57.8\% $ (from 26.8 to 42.3 days). A strong linear correlation ($R^2 = 0.96$) disclosed that lower fractional orders $\alpha$ impose earlier but less intensive interventions. This adaptive control framework offers the public health authorities a tool for designing an outbreak response that is resilient to real-life uncertainty.
| [1] |
B. G. Van den Hoogen, J. C. de Jong, J. Groen, T. Kuiken, R. de Groot, R. A. M. Fouchier, et al., A newly discovered human pneumovirus isolated from young children with respiratory tract disease, Nat. Med., 7 (2001), 719–724. https://doi.org/10.1038/89098 doi: 10.1038/89098
|
| [2] |
S. Panda, N. K. Mohakud, L. Pena, S. Kumar, Human metapneumovirus: Review of an important respiratory pathogen, Int. J. Infect. Dis., 25 (2014), 45–52. https://doi.org/10.1016/j.ijid.2014.03.1394 doi: 10.1016/j.ijid.2014.03.1394
|
| [3] |
K. Mohammadi, S. Faramarzi, S. Yaribash, Z. Valizadeh, E. Rajabi, M. Ghavam, et al. Human metapneumovirus (hMPV) in 2025: Emerging trends and insights from community and hospital-based respiratory panel analyses-a comprehensive review, Virol. J., 22 (2025), 150. https://doi.org/10.1186/s12985-025-02782-y doi: 10.1186/s12985-025-02782-y
|
| [4] |
W. Mahikul, L. J. White, K. Poovorawan, N. Soonthornworasiri, P. Sukontamarn, P. Chanthavilay, et al., Modeling household dynamics on respiratory syncytial virus (RSV), PLoS One, 14 (2019). https://doi.org/10.1371/journal.pone.0219323 doi: 10.1371/journal.pone.0219323
|
| [5] |
A. Sobanjo-ter Meulen, A. V. Gutierrez, O. Ruiz, J. Eeuwijk, H. Vroling, N. Kanesa-Thasan, The burden of human metapneumovirus (hMPV) disease in older and high-risk adults in developed countries: A systematic literature review, Infect. Dis. Ther., 14 (2025), 1917–1933. https://doi.org/10.1007/s40121-025-01187-2 doi: 10.1007/s40121-025-01187-2
|
| [6] |
W. Cao, J. Xu, F. Zhu, B. Cao, Addressing the unmet needs: Optimized vaccine design for human metapneumovirus, Innov. Med., 3 (2025), 100112. https://doi.org/10.59717/j.xinn-med.2025.100112 doi: 10.59717/j.xinn-med.2025.100112
|
| [7] | N. Muhammad, M. N. Alam, Z. Shiqing, Intelligent-based neural networks and optimal control of fractional order Ebola virus dynamics, arXiv preprint, 2025. |
| [8] |
H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000), 599–653. https://doi.org/10.1137/S0036144500371907 doi: 10.1137/S0036144500371907
|
| [9] |
N. Muhammad, M. N. Alam, Z. Shiqing, Modeling and analysis of dynamic waveforms in nonlinear fractional models of fifth order, Comput. Methods Diffe., 2025. https://doi.org/10.22034/cmde.2025.68014.3265 doi: 10.22034/cmde.2025.68014.3265
|
| [10] |
W. O. Kermack, A. G. McKendrick, A contribution to the mathematical theory of epidemics, P. Roy. Soc. A, 115 (1927), 700–721. https://doi.org/10.1098/rspa.1927.0118 doi: 10.1098/rspa.1927.0118
|
| [11] | F. Brauer, P. Driessche, J. H. Wu, Mathematical epidemiology, Springer, 2008. https://doi.org/10.1007/978-3-540-78911-6 |
| [12] |
P. Van den Driessche, J. Watmough, Reproduction numbers and subthreshold 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
|
| [13] | L. J. S. Allen, An introduction to stochastic epidemic models, in Mathematical Epidemiology, Springer, 2008, 81–130. |
| [14] |
M. Wang, P. van den Driessche, L. L. E. Cowen, J. Ma, Distributions of prevalence and daily new cases in a stochastic linear SEIR model, Math. Biosci., 388 (2025), 109508. https://doi.org/10.1016/j.mbs.2025.109508 doi: 10.1016/j.mbs.2025.109508
|
| [15] |
D. Temfack, J. Wyse, Sequential Monte Carlo squared for online inference in stochastic epidemic models, Epidemics, 52 (2025). https://doi.org/10.1016/j.epidem.2025.100847 doi: 10.1016/j.epidem.2025.100847
|
| [16] |
W. Xu, H. Liu, C. Qin, Dynamics of nonlinear stochastic SEIR infectious disease model with isolation and latency period, Symmetry, 17 (2025), 155. https://doi.org/10.3390/sym17020155 doi: 10.3390/sym17020155
|
| [17] |
Y. Sabbar, A. Zeb, N. Gul, D. Kiouach, S. P. Rajasekar, N. Ullah, et al., Stationary distribution of an SIR epidemic model with three correlated Brownian motions and general Lévy measure, AIMS Math., 8 (2023), 1329–1344. https://doi.org/10.3934/math.2023066 doi: 10.3934/math.2023066
|
| [18] |
A. Gray, D. Greenhalgh, L. Hu, X. Mao, J. Pan, A stochastic differential equation SIS epidemic model, SIAM J. Appl. Math., 71 (2011), 876–902. https://doi.org/10.1137/10081856X doi: 10.1137/10081856X
|
| [19] |
J. Du, C. Qin, Y. Hui, Optimal control and analysis of a stochastic SEIR epidemic model with nonlinear incidence and treatment, AIMS Math., 9 (2024), 33532–33550. https://doi.org/10.3934/math.20241600 doi: 10.3934/math.20241600
|
| [20] |
H. M. Wanjala, M. O. Okongo, J. O. Ochwach, Mathematical modelling with optimal control of infectious diseases with vaccination, Comput. Methods Diffe., 14 (2026), 616–632. https://doi.org/10.22034/cmde.2024.62350.2743 doi: 10.22034/cmde.2024.62350.2743
|
| [21] | K. K. Ali, K. Raslsn, A. F. Koura, M. A. Shaalan, Numerical treatment and optimal control of the hepatitis B virus spatio-temporal model, Comput. Methods Diffe., 14 (2026), 95–109. |
| [22] |
M. S. Hasan, M. N. Alam, M. Fayz-Al-Asad, N. Muhammad, C. Tunç, B-spline curve theory: An overview and applications in real life, Nonlinear Eng., 13 (2024), 20240054. https://doi.org/10.1515/nleng-2024-0054 doi: 10.1515/nleng-2024-0054
|
| [23] |
M. Qasim, H. Ali, A. Farooq, M. Kamran, H. Ahmad, F. A. Awwad, et al., Intelligent neural framework for modeling the lifestyle-induced remission in the type 2 diabetes, Chaos Soliton. Fract., 205 (2026), 117841. https://doi.org/10.1016/j.chaos.2025.117841 doi: 10.1016/j.chaos.2025.117841
|
| [24] |
L. Tang, H. Wang, X. Zhao, N. Xu, L. Li, Adaptive fixed-time bipartite containment control for saturated nonlinear multi-agent systems based on optimized backstepping technique, Math. Methods Appl. Sci., 205 (2026), 117841. https://doi.org/10.1016/j.chaos.2025.117841 doi: 10.1016/j.chaos.2025.117841
|
| [25] |
C. Chu, Y. He, A unified neural event-triggered control approach of high-order switched uncertain systems with time-varying state constraints, Robot. Intell. Automat., 46 (2026), 290–302. https://doi.org/10.1108/RIA-09-2025-0295 doi: 10.1108/RIA-09-2025-0295
|
| [26] |
S. N. K. A. Khan, M. Y. Misro, Hybrid B-spline collocation method with particle swarm optimization for solving linear differential problems, AIMS Math., 10 (2025), 5399–5420. https://doi.org/10.3934/math.2025249 doi: 10.3934/math.2025249
|
| [27] | A. Chertock, A. S. Iskhakov, S. Janajra, A. Kurganov, Spline-based stochastic collocation methods for uncertainty quantification in nonlinear hyperbolic PDEs, in European Conference on Numerical Mathematics and Advanced Applications, Springer, Cham, 2023,239–248. https://doi.org/10.1007/978-3-031-86173-4_24 |
| [28] |
J. C. Pedjeu, G. S. Ladde, Stochastic fractional differential equations: Modeling, method and analysis, Chaos Soliton. Fract., 45 (2012), 279–293. https://doi.org/10.1016/j.chaos.2011.12.009 doi: 10.1016/j.chaos.2011.12.009
|
| [29] |
H. Zhang, G. Chen, X. Wu, Y. Zhao, Y. Wang, Dynamic analysis of a fractional-order SEAIR model for influenza transmission with optimal control and stochastic stability, AIMS Math., 10 (2025), 20157–20198. https://doi.org/10.3934/math.2025901 doi: 10.3934/math.2025901
|
| [30] | Z. U. A. Zafar, S. Ijaz, C. Tunc, Numerical analysis of the SEIR epidemic model with fractional order, Comput. Methods Diffe., 14 (2026), 392–410. |
| [31] |
F. Evirgen, E. Uçar, S. Uçar, N. Özdemir, Modelling influenza a disease dynamics under Caputo-Fabrizio fractional derivative with distinct contact rates, Math. Model. Numer. Simul. Appl., 3 (2023), 58–73. https://doi.org/10.53391/mmnsa.1274004 doi: 10.53391/mmnsa.1274004
|
| [32] |
B. G. van den Hoogen, J. C. de Jong, J. Groen, T. Kuiken, R. de Groot, R. A. M. Fouchier, et al., A newly discovered human pneumovirus isolated from young children with respiratory tract disease, Nat. Med., 7 (2001), 719–724. https://doi.org/10.1038/89098 doi: 10.1038/89098
|
| [33] |
N. Muhammad, M. N. Alam, Z. Shiqing, Mathematical modeling and analysis of human metapneumovirus transmission dynamics using neural network intelligence and optimal control, Comput. Biol. Chem., 123 (2026), 109008. https://doi.org/10.1016/j.compbiolchem.2026.109008 doi: 10.1016/j.compbiolchem.2026.109008
|
| [34] |
S. Ghosh, P. J. Birrell, D. De Angelis, An approximate diffusion process for environmental stochasticity in infectious disease transmission modelling, PLoS Comput. Biol., 19 (2023), e1011088. https://doi.org/10.1371/journal.pcbi.1011088 doi: 10.1371/journal.pcbi.1011088
|
| [35] |
E. Tornatore, S. M. Buccellato, P. Vetro, Stability of a stochastic SIR system, Physica A, 354 (2005), 111–126. https://doi.org/10.1016/j.physa.2005.02.057 doi: 10.1016/j.physa.2005.02.057
|
| [36] | D. Matignon, Stability results for fractional differential equations with applications to control processing, in Computational Engineering in Systems Applications, Vol. 2, IMACS, IEEE-SMC, Lille, France, 1996,963–968. |
| [37] |
Y. Li, Y. Chen, I. Podlubny, Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability, Comput. Math. Appl., 59 (2010), 1810–1821. https://doi.org/10.1016/j.camwa.2009.08.019 doi: 10.1016/j.camwa.2009.08.019
|
| [38] |
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
|
| [39] | K. Diethelm, The analysis of fractional differential equations: An application-oriented exposition using differential operators of Caputo type, Springer-Verlag, 2010. https://doi.org/10.1007/978-3-642-14574-2 |
| [40] |
I. Area, J. Losada, J. J. Nieto, A note on the fractional logistic equation, Physica A, 444 (2015), 182–187. https://doi.org/10.1016/j.physa.2015.10.037 doi: 10.1016/j.physa.2015.10.037
|
| [41] |
E. Ahmed, A. M. A. El-Sayed, H. A. El-Saka, On some Routh-Hurwitz conditions for fractional order differential equations and their applications in Lorenz, Rössler, Chua and Chen systems, Phys. Lett. A, 358 (2006), 1–4. https://doi.org/10.1016/j.physleta.2006.04.087 doi: 10.1016/j.physleta.2006.04.087
|
| [42] |
N. Aguila-Camacho, M. A. Duarte-Mermoud, J. A. Gallegos, Lyapunov functions for fractional order systems, Commun. Nonlinear Sci., 19 (2014), 2951–2957. https://doi.org/10.1016/j.cnsns.2014.01.022 doi: 10.1016/j.cnsns.2014.01.022
|
| [43] | J. Carr, Applications of centre manifold theory, Springer Science & Business Media, 35 (2012). |
| [44] |
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
|
| [45] |
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
|
| [46] | A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, 204 (2006). |
| [47] |
B. Mezerdi, S. Bahlali, Necessary conditions for optimality in relaxed stochastic control problems, Stochastics, 73 (2002), 201–218. https://doi.org/10.1080/1045112021000025925 doi: 10.1080/1045112021000025925
|
| [48] | Y. S. Mishura, Stochastic calculus for fractional Brownian motion and related processes, Springer Science & Business Media, 2008. https://doi.org/10.1007/978-3-540-75873-0 |
| [49] |
C. Sorgentone, E. Pellegrino, F. Pitolli, A spline-based framework for solving the space-time fractional convection-diffusion problem, Appl. Math. Lett., 161 (2025), 109370. https://doi.org/10.1016/j.aml.2024.109370 doi: 10.1016/j.aml.2024.109370
|
| [50] | S. Malge, R. K. Lodhi, Quintic B-spline method for numerical solution of second-order singularly perturbed delay differential equations, Comput. Methods Diffe., 14 (2026), 36–49. |
| [51] |
M. S. Arif, K. Abodayeh, Y. Nawaz, Precision in disease dynamics: Finite difference solutions for stochastic epidemics with treatment cure and partial immunity, Part. Differ. Equ. Appl. Math., 9 (2024), 100660. https://doi.org/10.1016/j.padiff.2024.100660 doi: 10.1016/j.padiff.2024.100660
|
Figures(7) / Tables(3)
Md Hossain Islam, Md. Nur Alam, Noor Muhammad, Xinsong Yang, Md. Jakir Hossen. Mathematical modeling of human metapneumovirus transmission dynamics with optimal control analysis[J]. AIMS Mathematics, 2026, 11(7): 21570-21606. doi: 10.3934/math.2026874
DownLoad: