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

Front Immunol. 2024 May 7:15:1384509. doi: 10.3389/fimmu.2024.1384509. eCollection 2024.

Abstract

Introduction: The 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 methods: We 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).

Results: The 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).

Discussion: Systematic 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: JAK2 V617F; blood cancer; mathematical modelling; myeloproliferative neoplasms (MPN); ordinary differential equations; parameter estimation; ruxolitinib; stem cells.

MeSH terms

  • Humans
  • Janus Kinase 2* / antagonists & inhibitors
  • Janus Kinase 2* / genetics
  • Models, Theoretical
  • Myeloproliferative Disorders* / drug therapy
  • Myeloproliferative Disorders* / genetics
  • Neoplastic Stem Cells / drug effects
  • Nitriles*
  • Protein Kinase Inhibitors / pharmacology
  • Protein Kinase Inhibitors / therapeutic use
  • Pyrazoles* / pharmacology
  • Pyrazoles* / therapeutic use
  • Pyrimidines* / therapeutic use

Substances

  • Nitriles
  • ruxolitinib
  • Pyrazoles
  • Pyrimidines
  • Janus Kinase 2
  • JAK2 protein, human
  • Protein Kinase Inhibitors

Grants and funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the fellowship “Personalized prediction of blood cancer progression using clinical data and mathematical modelling” (R335-2019-2020) from the Lundbeck Foundation (PI: TS).