Abstract

Chronic myeloid leukemia (CML) is a blood cancer characterized by dysregulated production of maturing myeloid cells driven by the product of the Philadelphia chromosome, the BCR-ABL1 tyrosine kinase. Tyrosine kinase inhibitors (TKI) have proved effective in treating CML but there is still a cohort of patients who do not respond to TKI therapy even in the absence of mutations in the BCR-ABL1 kinase domain that mediate drug resistance. To discover novel strategies to improve TKI therapy in CML, we developed a nonlinear mathematical model of CML hematopoiesis that incorporates feedback control and lineage branching. Cell-cell interactions were constrained using an automated model selection method together with previous observations and new in vivo data from a chimeric BCR-ABL1 transgenic mouse model of CML. The resulting quantitative model captures the dynamics of normal and CML cells at various stages of the disease and exhibits variable responses to TKI treatment, consistent with those of CML patients. The model predicts that an increase in the proportion of CML stem cells in the bone marrow would decrease the tendency of the disease to respond to TKI therapy, in concordance with clinical data and confirmed experimentally in mice. The model further suggests that, under our assumed similarities between normal and leukemic cells, a key predictor of refractory response to TKI treatment is an increased maximum probability of self-renewal of normal hematopoietic stem cells. We use these insights to develop a clinical prognostic criterion to predict the efficacy of TKI treatment and to design strategies to improve treatment response. The model predicts that stimulating the differentiation of leukemic stem cells while applying TKI therapy can significantly improve treatment outcomes.

Data availability

Modelling code and parameter set data are available in a Github repository. Patient data is unavailable publicly as it could be used to potentially identify the patients. Deidentified raw patient transcript data will be made available to qualified researchers (academic or industry) upon request to Dr. Van Etten at vanetten@hs.uci.edu.

The following data sets were generated

Article and author information

Author details

  1. Jonathan Rodriguez

    Graduate Program in Mathematical, Computational and Systems Biology, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6414-0526
  2. Abdon Iniguez

    Graduate Program in Mathematical, Computational and Systems Biology, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Nilamani Jena

    Department of Medicine, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Prasanthi Tata

    Department of Medicine, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Joan Liu

    Department of Medicine, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Arthur D Lander

    Center for Complex Biological Systems, University of California, Irvine, Irvine, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4380-5525
  7. John S Lowengrub

    Department of Mathematics, University of California, Irvine, Irvine, United States
    For correspondence
    lowengrb@math.uci.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1759-0900
  8. Richard A Van Etten

    Chao Family Comprehensive Cancer Center, University of California, Irvine, Irvine, United States
    For correspondence
    vanetten@hs.uci.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institutes of Health (1U54CA217378-01A1)

  • Arthur D Lander
  • John S Lowengrub
  • Richard A Van Etten

National Institutes of Health (P30CA062203)

  • Arthur D Lander
  • John S Lowengrub
  • Richard A Van Etten

National Institutes of Health (R01 CA090576)

  • Richard A Van Etten

National Science Foundation (DMS-1763272)

  • Arthur D Lander
  • John S Lowengrub

National Science Foundation (DMS-1936833)

  • John S Lowengrub

National Science Foundation (DMS-1953410)

  • John S Lowengrub

National Science Foundation (GRFP 16-588)

  • Abdon Iniguez

Simons Foundation (594598QN)

  • Arthur D Lander
  • John S Lowengrub

National Institute of General Medical Sciences (GM136624)

  • Jonathan Rodriguez

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocol (AUP-19-159) of the University of California, Irvine.

Copyright

© 2023, Rodriguez et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 551
    views
  • 123
    downloads
  • 3
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Jonathan Rodriguez
  2. Abdon Iniguez
  3. Nilamani Jena
  4. Prasanthi Tata
  5. Joan Liu
  6. Arthur D Lander
  7. John S Lowengrub
  8. Richard A Van Etten
(2023)
Predictive nonlinear modeling of malignant myelopoiesis and tyrosine kinase inhibitor therapy
eLife 12:e84149.
https://doi.org/10.7554/eLife.84149

Share this article

https://doi.org/10.7554/eLife.84149

Further reading

    1. Cancer Biology
    Ismail M Meraz, Mourad Majidi ... Jack A Roth
    Research Article

    Expression of NPRL2/TUSC4, a tumor-suppressor gene, is reduced in many cancers including NSCLC. Restoration of NPRL2 induces DNA damage, apoptosis, and cell-cycle arrest. We investigated NPRL2 antitumor immune responses in aPD1R/KRAS/STK11mt NSCLC in humanized-mice. Humanized-mice were generated by transplanting fresh human cord blood-derived CD34 stem cells into sub-lethally irradiated NSG mice. Lung-metastases were developed from KRAS/STK11mt/aPD1R A549 cells and treated with NPRL2 w/wo pembrolizumab. NPRL2-treatment reduced lung metastases significantly, whereas pembrolizumab was ineffective. Antitumor effect was greater in humanized than non-humanized-mice. NPRL2 + pembrolizumab was not synergistic in KRAS/STK11mt/aPD1R tumors but was synergistic in KRASwt/aPD1S H1299. NPRL2 also showed a significant antitumor effect on KRASmt/aPD1R LLC2 syngeneic-tumors. The antitumor effect was correlated with increased infiltration of human cytotoxic-T, HLA-DR+DC, CD11c+DC, and downregulation of myeloid and regulatory-T cells in TME. Antitumor effect was abolished upon in-vivo depletion of CD8-T, macrophages, and CD4-T cells whereas remained unaffected upon NK-cell depletion. A distinctive protein-expression profile was found after NPRL2 treatment. IFNγ, CD8b, and TBX21 associated with T-cell functions were significantly increased, whereas FOXP3, TGFB1/B2, and IL-10RA were strongly inhibited by NPRL2. A list of T-cell co-inhibitory molecules was also downregulated. Restoration of NPRL2 exhibited significantly slower tumor growth in humanized-mice, which was associated with increased presence of human cytotoxic-T, and DC and decreased percentage of Treg, MDSC, and TAM in TME. NPRL2-stable cells showed a substantial increase in colony-formation inhibition and heightened sensitivity to carboplatin. Stable-expression of NPRL2 resulted in the downregulation of MAPK and AKT-mTOR signaling. Taken-together, NPRL2 gene-therapy induces antitumor activity on KRAS/STK11mt/aPD1R tumors through DC-mediated antigen-presentation and cytotoxic immune-cell activation.

    1. Biochemistry and Chemical Biology
    2. Cancer Biology
    Flavie Coquel, Sing-Zong Ho ... Philippe Pasero
    Research Article

    Cancer cells display high levels of oncogene-induced replication stress (RS) and rely on DNA damage checkpoint for viability. This feature is exploited by cancer therapies to either increase RS to unbearable levels or inhibit checkpoint kinases involved in the DNA damage response. Thus far, treatments that combine these two strategies have shown promise but also have severe adverse effects. To identify novel, better-tolerated anticancer combinations, we screened a collection of plant extracts and found two natural compounds from the plant, Psoralea corylifolia, that synergistically inhibit cancer cell proliferation. Bakuchiol inhibited DNA replication and activated the checkpoint kinase CHK1 by targeting DNA polymerases. Isobavachalcone interfered with DNA double-strand break repair by inhibiting the checkpoint kinase CHK2 and DNA end resection. The combination of bakuchiol and isobavachalcone synergistically inhibited cancer cell proliferation in vitro. Importantly, it also prevented tumor development in xenografted NOD/SCID mice. The synergistic effect of inhibiting DNA replication and CHK2 signaling identifies a vulnerability of cancer cells that might be exploited by using clinically approved inhibitors in novel combination therapies.