Model-Based Optimal AML Consolidation Treatment

Objective: Neutropenia is an adverse event commonly arising during intensive chemotherapy of acute myeloid leukemia (AML). It is often associated with infectious complications. Mathematical modeling, simulation, and optimization of the treatment process would be a valuable tool to support clinical decision making, potentially resulting in less severe side effects and deeper remissions. However, until now, there has been no validated mathematical model available to simulate the effect of chemotherapy treatment on white blood cell (WBC) counts and leukemic cells simultaneously. Methods: We developed a population pharmacokinetic/pharmacodynamic (PK/PD) model combining a myelosuppression model considering endogenous granulocyte-colony stimulating factor (G-CSF), a PK model for cytarabine (Ara-C), a subcutaneous absorption model for exogenous G-CSF, and a two-compartment model for leukemic blasts. This model was fitted to data of 44 AML patients during consolidation therapy with a novel Ara-C plus G-CSF schedule from a phase II controlled clinical trial. Additionally, we were able to optimize treatment schedules with respect to disease progression, WBC nadirs, and the amount of Ara-C and G-CSF. Results: The developed PK/PD model provided good prediction accuracies and an interpretation of the interaction between WBCs, G-CSF, and blasts. For 14 patients (those with available bone marrow blast counts), we achieved a median 4.2-fold higher WBC count at nadir, which is the most critical time during consolidation therapy. The simulation results showed that relative bone marrow blast counts remained below the clinically important threshold of 5%, with a median of 60% reduction in Ara-C. Conclusion: These in silico findings demonstrate the benefits of optimized treatment schedules for AML patients. Significance: Until 2017, no new drug had been approved for the treatment of AML, fostering the optimal use of currently available drugs.

[1]  Mats O. Karlsson,et al.  Characterization of Endogenous G-CSF and the Inverse Correlation to Chemotherapy-Induced Neutropenia in Patients with Breast Cancer Using Population Modeling , 2014, Pharmaceutical Research.

[2]  Sebastian Sager,et al.  Optimal Control for Cancer Chemotherapy ODE Models: Potential of Optimal Schedules and Choice of Objective Function , 2010 .

[3]  Bob Löwenberg,et al.  Diagnosis and management of AML in adults: 2017 ELN recommendations from an international expert panel. , 2017, Blood.

[4]  A. Marciniak-Czochra,et al.  Mathematical Modeling of Leukemogenesis and Cancer Stem Cell Dynamics , 2012 .

[5]  R. Mayer,et al.  Intensive postremission chemotherapy in adults with acute myeloid leukemia. Cancer and Leukemia Group B. , 1994, The New England journal of medicine.

[6]  Markus Scholz,et al.  Model-based optimization of G-CSF treatment during cytotoxic chemotherapy , 2017, Journal of Cancer Research and Clinical Oncology.

[7]  Verena I Gaidzik,et al.  Measurable Residual Disease Monitoring in Acute Myeloid Leukemia with t(8;21)(q22;q22.1): Results of the AML Study Group. , 2019, Blood.

[8]  Sebastian Sager,et al.  Mathematical models for cytarabine-derived myelosuppression in acute myeloid leukaemia , 2019, PloS one.

[9]  Leonid Kagan,et al.  Pharmacokinetic Modeling of the Subcutaneous Absorption of Therapeutic Proteins , 2014, Drug Metabolism and Disposition.

[10]  Gisela Schwab,et al.  Pharmacokinetic/Pharmacodynamic Modeling of Pegfilgrastim in Healthy Subjects , 2006, Journal of clinical pharmacology.

[11]  Moritz Diehl,et al.  CasADi: a software framework for nonlinear optimization and optimal control , 2018, Mathematical Programming Computation.

[12]  S. I. Rubinow,et al.  A mathematical model of the chemotherapeutic treatment of acute myeloblastic leukemia. , 1976, Biophysical journal.

[13]  Efstratios N. Pistikopoulos,et al.  Chemotherapy Drug Scheduling for the Induction Treatment of Patients With Acute Myeloid Leukemia , 2014, IEEE Transactions on Biomedical Engineering.

[14]  J. Dick,et al.  Acute myeloid leukemia originates from a hierarchy of leukemic stem cell classes that differ in self-renewal capacity , 2004, Nature Immunology.

[15]  Eleni Pefani,et al.  Modelling and optimisation based drug delivery systems for the treatment of Acute Myeloid Leukemia (AML) , 2013 .

[16]  Neil D. Evans,et al.  Structural identifiability for mathematical pharmacology: models of myelosuppression , 2018, Journal of Pharmacokinetics and Pharmacodynamics.

[17]  Vered Rom-Kedar,et al.  Novel Strategies for Granulocyte Colony-Stimulating Factor Treatment of Severe Prolonged Neutropenia Suggested by Mathematical Modeling , 2008, Clinical Cancer Research.

[18]  Yu Wang,et al.  MRD-directed risk stratification treatment may improve outcomes of t(8;21) AML in the first complete remission: results from the AML05 multicenter trial. , 2013, Blood.

[19]  J. Cannon,et al.  Pharmacokinetics of 1-beta-D-arabinofuranosylcytosine (ARA-C) deamination in several species. , 1973, Biochemical pharmacology.

[20]  Giancarlo Mauri,et al.  Modeling cell proliferation in human acute myeloid leukemia xenografts , 2019, Bioinform..

[21]  Michael C Mackey,et al.  Mathematical model for G-CSF administration after chemotherapy. , 2009, Journal of theoretical biology.

[22]  A V Savkin,et al.  Optimal chemotherapy regimens: influence of tumours on normal cells and several toxicity constraints. , 2001, IMA journal of mathematics applied in medicine and biology.

[23]  M Marusić,et al.  Relationship between differing volumes of bone marrow aspirates and their cellular composition. , 1990, Bone marrow transplantation.

[24]  Wei Zhang,et al.  Exploring the Antitumor Mechanism of High-Dose Cytarabine through the Metabolic Perturbations of Ribonucleotide and Deoxyribonucleotide in Human Promyelocytic Leukemia HL-60 Cells , 2017, Molecules.

[25]  A. Ferreri,et al.  Clinical implications, safety, efficacy of Recombinant Human Granulocyte Colony-Stimulating Factors and Pegylated Equivalent , 2013, Epidemiology, Biostatistics, and Public Health.

[26]  Torsten Haferlach,et al.  Minimal/measurable residual disease in AML: a consensus document from the European LeukemiaNet MRD Working Party. , 2018, Blood.

[27]  Wolfgang Kern,et al.  High antileukemic activity of sequential high dose cytosine arabinoside and mitoxantrone in patients with refractory acute leukemias , 1997, Cancer.

[28]  Michael Stadler,et al.  Measurable residual disease monitoring by NGS before allogeneic hematopoietic cell transplantation in AML. , 2018, Blood.

[29]  Julio R. Banga,et al.  SUPPLEMENTARY FILE: Full observability and estimation of unknown inputs, states, and parameters of nonlinear biological models , 2019 .

[30]  Bob Löwenberg,et al.  Criteria for defining a complete remission in acute myeloid leukaemia revisited. An analysis of patients treated in HOVON‐SAKK co‐operative group studies , 2005, British journal of haematology.

[31]  R. Schlenk,et al.  Is there justification for 4 cycles of consolidation therapy in AML? , 2016, Best practice & research. Clinical haematology.

[32]  A Benner,et al.  Condensed versus standard schedule of high-dose cytarabine consolidation therapy with pegfilgrastim growth factor support in acute myeloid leukemia , 2016, Blood Cancer Journal.

[33]  Lorenz T. Biegler,et al.  Nonlinear Waves in Integrable and Nonintegrable Systems , 2018 .

[34]  Mats O Karlsson,et al.  Model of chemotherapy-induced myelosuppression with parameter consistency across drugs. , 2002, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[35]  E. Afenya,et al.  Acute leukemia and chemotherapy: a modeling viewpoint. , 1996, Mathematical biosciences.

[36]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[37]  G Mele,et al.  Fludarabine, cytarabine, and G‐CSF (FLAG) for the treatment of poor risk acute myeloid leukemia , 1998, American journal of hematology.

[38]  Emilio Fernández-Varón,et al.  Granulocyte and granulocyte macrophage colony-stimulating factors as therapy in human and veterinary medicine. , 2007, Veterinary journal.

[39]  J. Gabrilove,et al.  Filgrastim (r-metHuG-CSF): the first 10 years. , 1996, Blood.

[40]  T. Inoue,et al.  Pharmacokinetics and adverse events following 5-day repeated administration of lenograstim, a recombinant human granulocyte colony-stimulating factor, in healthy subjects , 2000, Bone Marrow Transplantation.

[41]  Mats O Karlsson,et al.  Population Pharmacokinetic-Pharmacodynamic Model for Neutropenia with Patient Subgroup Identification: Comparison across Anticancer Drugs , 2006, Clinical Cancer Research.

[42]  J. Dick,et al.  Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell , 1997, Nature Medicine.

[43]  M. Oken,et al.  Maintenance chemotherapy prolongs remission duration in adult acute nonlymphocytic leukemia. , 1988, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[44]  Andrew C. Hooker,et al.  A Population Pharmacokinetic-Pharmacodynamic Model of Pegfilgrastim , 2018, The AAPS Journal.

[45]  Markus Scholz,et al.  Pharmacokinetic and -dynamic modelling of G-CSF derivatives in humans , 2012, Theoretical Biology and Medical Modelling.

[46]  Wojciech Krzyzanski,et al.  Population Modeling of Filgrastim PK‐PD in Healthy Adults Following Intravenous and Subcutaneous Administrations , 2010, Journal of clinical pharmacology.

[47]  E. Yukawa,et al.  Estimation of rhG-CSF absorption kinetics after subcutaneous administration using a modified Wagner-Nelson method with a nonlinear elimination model. , 2001, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[48]  A Vacca,et al.  Comparison between filgrastim and lenograstim plus chemotherapy for mobilization of PBPCs , 2010, Bone Marrow Transplantation.

[49]  A. Benner,et al.  Randomized phase-II trial evaluating induction therapy with idarubicin and etoposide plus sequential or concurrent azacitidine and maintenance therapy with azacitidine , 2019, Leukemia.

[50]  S. Higuchi,et al.  Pharmacokinetic and Pharmacodynamic Analysis of Subcutaneous Recombinant Human Granulocyte Colony Stimulating Factor (Lenograstim) Administration , 1999, Journal of clinical pharmacology.

[51]  H. Preisler,et al.  Relationship between plasma Ara-C and intracellular Ara-CTP pools under conditions of continuous infusion and high-dose Ara-C treatment. , 1982, Medical and pediatric oncology.

[52]  Roman Jaksik,et al.  Predicting Minimal Residual Disease in Acute Myeloid Leukemia through Stochastic Modeling of Clonality , 2019, bioRxiv.

[53]  Jan Liliemark,et al.  Relationship of 1-β-d-Arabinofuranosylcytosine in Plasma to 1-β-d-Arabinofuranosylcytosine 5′-Triphosphate Levels in Leukemic Cells during Treatment with High-Dose 1-β-d-Arabinofuranosylcytosine , 1985 .

[54]  L. Biegler An overview of simultaneous strategies for dynamic optimization , 2007 .

[55]  Zinnia P. Parra-Guillen,et al.  Semimechanistic Bone Marrow Exhaustion Pharmacokinetic/Pharmacodynamic Model for Chemotherapy-Induced Cumulative Neutropenia , 2017, The Journal of Pharmacology and Experimental Therapeutics.

[56]  Doraiswami Ramkrishna,et al.  Optimal Chemotherapy for Leukemia: A Model-Based Strategy for Individualized Treatment , 2014, PloS one.

[57]  HT Swan,et al.  A Colour Atlas of Thymus and Lymph Node Histopathology with Ultrastructure , 1982 .

[58]  Eva Balsa-Canto,et al.  Bioinformatics Applications Note Systems Biology Genssi: a Software Toolbox for Structural Identifiability Analysis of Biological Models , 2022 .

[59]  P. Becker,et al.  Controversies Regarding Use of Myeloid Growth Factors in Leukemia. , 2017, Journal of the National Comprehensive Cancer Network : JNCCN.

[60]  Bing Wang,et al.  Population Pharmacokinetic–Pharmacodynamic Modeling of Filgrastim (r-metHuG-CSF) in Healthy Volunteers , 2001, Journal of Pharmacokinetics and Pharmacodynamics.

[61]  Bing‐Bing Yang,et al.  Pharmacokinetics and Pharmacodynamics of Pegfilgrastim , 2011, Clinical pharmacokinetics.

[62]  M. Manz,et al.  Quantification and three-dimensional microanatomical organization of the bone marrow. , 2017, Blood advances.

[63]  Philippe Jacqmin,et al.  Pharmacokinetic–pharmacodynamic modelling of neutrophil response to G‐CSF in healthy subjects and patients with chemotherapy‐induced neutropenia , 2018, British journal of clinical pharmacology.

[64]  B. Djulbegovic,et al.  Mathematical Model of Acute Myeloblastic Leukaemia: an Investigation of the Relevant Kinetic Parameters , 1985, Cell and tissue kinetics.

[65]  Leonard Leibovici,et al.  Subcutaneous versus intravenous granulocyte colony stimulating factor for the treatment of neutropenia in hospitalized hemato‐oncological patients: Randomized controlled trial , 2014, American journal of hematology.

[66]  Jordi Esteve,et al.  Treatment with G-CSF reduces acute myeloid leukemia blast viability in the presence of bone marrow stroma , 2015, Cancer Cell International.

[67]  Thomas Stiehl,et al.  Mathematical modeling of the impact of cytokine response of acute myeloid leukemia cells on patient prognosis , 2018, Scientific Reports.

[68]  Wojciech Krzyzanski,et al.  Population Pharmacokinetic Modelling of Filgrastim in Healthy Adults following Intravenous and Subcutaneous Administrations , 2009, Clinical pharmacokinetics.

[69]  Stephen P. Boyd,et al.  Extending Scope of Robust Optimization: Comprehensive Robust Counterparts of Uncertain Problems , 2006, Math. Program..

[70]  R Ohno,et al.  G-CSF in the treatment of acute myeloid leukemia: is it safe? , 1993, Leukemia & lymphoma.

[71]  M. Aapro,et al.  Management of febrile neutropaenia: ESMO Clinical Practice Guidelines. , 2010, Annals of oncology : official journal of the European Society for Medical Oncology.

[72]  Arijit Chakravarty,et al.  Dose Schedule Optimization and the Pharmacokinetic Driver of Neutropenia , 2014, PloS one.