From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth
Article Information
Abstract
This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm$^3$ to 157~cm$^3$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction.
Graphical Abstract
Keywords
Data Availability Statement
Funding
Conflicts of Interest
AI Use Statement
Ethical Approval and Consent to Participate
References
- American Cancer Society. (2025). Cancer facts & figures 2025. American Cancer Society. Retrieved from https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/2025-cancer-facts-figures.html
[Google Scholar] - National Cancer Institute. (n.d.). Cancer stat facts: Leukemia. Surveillance, Epidemiology, and End Results Program. Retrieved from https://seer.cancer.gov/statfacts/html/leuks.html
[Google Scholar] - Mackey, M. C., & Glass, L. (1977). Oscillation and chaos in physiological control systems. Science, 197(4300), 287-289.
[CrossRef] [Google Scholar] - Berezansky, L., Braverman, E., & Idels, L. (2013). Mackey-Glass model of hematopoiesis with monotone feedback revisited. Applied Mathematics and Computation, 219(9), 4892--4907.
[CrossRef] [Google Scholar] - Berezansky, L., Braverman, E., & Idels, L. (2013). Mackey--Glass model of hematopoiesis with non-monotone feedback: Stability, oscillation and control. Applied Mathematics and Computation, 219(11), 6268--6283.
[CrossRef] [Google Scholar] - Adimy, M., Crauste, F., & Ruan, S. (2005). A mathematical study of the hematopoiesis process with applications to chronic myelogenous leukemia. SIAM journal on applied mathematics, 65(4), 1328-1352.
[CrossRef] [Google Scholar] - Pujo-Menjouet, L. (2016). Blood cell dynamics: half of a century of modelling. Mathematical Modelling of Natural Phenomena, 11(1), 92-115.
[CrossRef] [Google Scholar] - Majdoul, S., Colson-Durand, L., To, N. H., & Belkacemi, Y. (2017). Adaptive radiotherapy for an uncommon chloroma. Case Reports in Oncology, 9(3), 593-598.
[CrossRef] [Google Scholar] - Shampine, L. F., & Thompson, S. (2001). Solving ddes in MATLAB. Applied Numerical Mathematics, 37(4), 441-458.
[CrossRef] [Google Scholar] - Crauste, F. (2009). Delay model of hematopoietic stem cell dynamics: asymptotic stability and stability switch. Mathematical Modelling of Natural Phenomena, 4(2), 28-47.
[CrossRef] [Google Scholar] - Halanay, A., Cândea, D., & Rădulescu, I. R. (2014). Existence and Stability of Limit Cycles in a Two-delays Model of Hematopoiesis Including Asymmetric Division. Mathematical Modelling of Natural Phenomena, 9(1), 58-78.
[CrossRef] [Google Scholar] - Badralexi, I., & Halanay, A. (2017, January). Stability and oscillations in a CML model. In AIP Conference Proceedings (Vol. 1798, No. 1, p. 020011). AIP Publishing LLC.
[CrossRef] [Google Scholar] - Cândea, D., Halanay, A., Rădulescu, R., & Tălmaci, R. (2017, January). Parameter estimation and sensitivity analysis for a mathematical model with time delays of leukemia. In AIP Conference Proceedings (Vol. 1798, No. 1, p. 020034). AIP Publishing LLC.
[CrossRef] [Google Scholar] - Parajdi, L. G., Precup, R., Bonci, E. A., & Tomuleasa, C. (2020). A mathematical model of the transition from normal hematopoiesis to the chronic and accelerated-acute stages in myeloid leukemia. Mathematics, 8(3), 376.
[CrossRef] [Google Scholar] - Bakst, R. L., Dabaja, B. S., Specht, L. K., & Yahalom, J. (2018). Use of radiation in extramedullary leukemia/chloroma: guidelines from the International Lymphoma Radiation Oncology Group. International Journal of Radiation Oncology* Biology* Physics, 102(2), 314-319.
[CrossRef] [Google Scholar] - Loscocco, G. G., & Vannucchi, A. M. (2023). Myeloid sarcoma: more and less than a distinct entity. Annals of Hematology, 102(8), 1973-1984.
[CrossRef] [Google Scholar] - Adimy, M., Crauste, F., & Ruan, S. (2006). Periodic oscillations in leukopoiesis models with two delays. Journal of Theoretical Biology, 242(2), 288-299.
[CrossRef] [Google Scholar] - Colijn, C., & Mackey, M. C. (2005). A mathematical model of hematopoiesis—I. Periodic chronic myelogenous leukemia. Journal of Theoretical Biology, 237(2), 117-132.
[CrossRef] [Google Scholar] - Safarishahrbijari, A., & Gaffari, A. (2013). Parameter identification of hematopoiesis mathematical model–periodic chronic myelogenous leukemia. Contemporary Oncology/Współczesna Onkologia, 17(1), 73-77.
[CrossRef] [Google Scholar] - Foley, C., & Mackey, M. C. (2009). Dynamic hematological disease: a review. Journal of mathematical biology, 58(1), 285-322.
[CrossRef] [Google Scholar]
Cite This Article
TY - JOUR AU - Farias, Matheus dos Santos PY - 2026 DA - 2026/09/20 TI - From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth JO - Journal of Numerical Simulations in Physics and Mathematics T2 - Journal of Numerical Simulations in Physics and Mathematics JF - Journal of Numerical Simulations in Physics and Mathematics VL - 2 IS - 2 SP - 130 EP - 140 DO - 10.62762/JNSPM.2026.604346 UR - https://www.icck.org/article/abs/JNSPM.2026.604346 KW - Mackey-Glass equation KW - leukemia growth KW - delay differential equations KW - mathematical oncology KW - radiotherapy KW - numerical simulation AB - This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm$^3$ to 157~cm$^3$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction. SN - 3068-9082 PB - Institute of Central Computation and Knowledge LA - English ER -
@article{Farias2026From,
author = {Matheus dos Santos Farias},
title = {From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth},
journal = {Journal of Numerical Simulations in Physics and Mathematics},
year = {2026},
volume = {2},
number = {2},
pages = {130-140},
doi = {10.62762/JNSPM.2026.604346},
url = {https://www.icck.org/article/abs/JNSPM.2026.604346},
abstract = {This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm\$^3\$ to 157~cm\$^3\$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction.},
keywords = {Mackey-Glass equation, leukemia growth, delay differential equations, mathematical oncology, radiotherapy, numerical simulation},
issn = {3068-9082},
publisher = {Institute of Central Computation and Knowledge}
}
Article Metrics
Publisher's Note
ICCK stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and Permissions
Copyright © 2026 by the Author(s). Published by Institute of Central Computation and Knowledge. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
Portico