Frontiers in Immunology (May 2024)

Mathematical modelling of stem and progenitor cell dynamics during ruxolitinib treatment of patients with myeloproliferative neoplasms

  • Tobias Idor Boklund,
  • Jordan Snyder,
  • Johanne Gudmand-Hoeyer,
  • Morten Kranker Larsen,
  • Trine Alma Knudsen,
  • Christina Schjellerup Eickhardt-Dalbøge,
  • Vibe Skov,
  • Lasse Kjær,
  • Hans C. Hasselbalch,
  • Morten Andersen,
  • Johnny T. Ottesen,
  • Thomas Stiehl,
  • Thomas Stiehl

DOI
https://doi.org/10.3389/fimmu.2024.1384509
Journal volume & issue
Vol. 15

Abstract

Read online

IntroductionThe Philadelphia chromosome-negative myeloproliferative neoplasms are a group of slowly progressing haematological malignancies primarily characterised by an overproduction of myeloid blood cells. Patients are treated with various drugs, including the JAK1/2 inhibitor ruxolitinib. Mathematical modelling can help propose and test hypotheses of how the treatment works.Materials and methodsWe present an extension of the Cancitis model, which describes the development of myeloproliferative neoplasms and their interactions with inflammation, that explicitly models progenitor cells and can account for treatment with ruxolitinib through effects on the malignant stem cell response to cytokine signalling and the death rate of malignant progenitor cells. The model has been fitted to individual patients’ data for the JAK2 V617F variant allele frequency from the COMFORT-II and RESPONSE studies for patients who had substantial reductions (20 percentage points or 90% of the baseline value) in their JAK2 V617F variant allele frequency (n = 24 in total).ResultsThe model fits very well to the patient data with an average root mean square error of 0.0249 (2.49%) when allowing ruxolitinib treatment to affect both malignant stem and progenitor cells. This average root mean square error is much lower than if allowing ruxolitinib treatment to affect only malignant stem or only malignant progenitor cells (average root mean square errors of 0.138 (13.8%) and 0.0874 (8.74%), respectively).DiscussionSystematic simulation studies and fitting of the model to the patient data suggest that an initial reduction of the malignant cell burden followed by a monotonic increase can be recapitulated by the model assuming that ruxolitinib affects only the death rate of malignant progenitor cells. For patients exhibiting a long-term reduction of the malignant cells, the model predicts that ruxolitinib also affects stem cell parameters, such as the malignant stem cells’ response to cytokine signalling.

Keywords