Revisiting dosing regimen using PK/PD modeling: the MODEL1 phase I/II trial of docetaxel plus epirubicin in metastatic breast cancer patients

Dynamical behavior-based approach for the evaluation of treatment efficacy: The case of immuno-oncology

Sophistication of mathematical models in the pharmacological context reflects the progress being made in understanding physiological, pharmacological, and disease relationships. This progress has illustrated once more the need for advanced quantitative tools able to efficiently extract information from these models. While dynamical systems theory has a long history in the analysis of systems biology models, as emphasized under the dynamical disease concept by Mackey and Glass [Science 197, 287-289 (1977)], its adoption in pharmacometrics is only at the beginning [Chae, Transl. Clin. Pharmacol. 28, 109 (2020)]. Using a quantitative systems pharmacology model of tumor immune dynamics as a case study [Kosinsky et al., J. Immunother. Cancer 6, 17 (2018)], we here adopt a dynamical systems analysis to describe, in an exhaustive way, six different statuses that refer to the response of the system to therapy, in the presence or absence of a tumor-free attractor. To evaluate the therapy success, we introduce the concept of TBA, related to the Time to enter the tumor-free Basin of Attraction, and corresponding to the earliest time at which the therapy can be stopped without jeopardizing its efficacy. TBA can determine the optimal time to stop drug administration and consequently quantify the reduction in drug exposure.

Treatment of evolving cancers will require dynamic decision support

Cancer research has traditionally focused on developing new agents, but an underexplored question is that of the dose and frequency of existing drugs. Based on the modus operandi established in the early days of chemotherapies, most drugs are administered according to predetermined schedules that seek to deliver the maximum tolerated dose and are only adjusted for toxicity. However, we believe that the complex, evolving nature of cancer requires a more dynamic and personalized approach. Chronicling the milestones of the field, we show that the impact of schedule choice crucially depends on processes driving treatment response and failure. As such, cancer heterogeneity and evolution dictate that a one-size-fits-all solution is unlikely-instead, each patient should be mapped to the strategy that best matches their current disease characteristics and treatment objectives (i.e. their 'tumorscape'). To achieve this level of personalization, we need mathematical modeling. In this perspective, we propose a five-step 'Adaptive Dosing Adjusted for Personalized Tumorscapes (ADAPT)' paradigm to integrate data and understanding across scales and derive dynamic and personalized schedules. We conclude with promising examples of model-guided schedule personalization and a call to action to address key outstanding challenges surrounding data collection, model development, and integration.

Adaptive therapy for ovarian cancer: An integrated approach to PARP inhibitor scheduling

Toxicity and emerging drug resistance are important challenges in PARP inhibitor (PARPi) treatment of ovarian cancer. Recent research has shown that evolutionary-inspired treatment algorithms which adapt treatment to the tumor's treatment response (adaptive therapy) can help to mitigate both. Here, we present a first step in developing an adaptive therapy protocol for PARPi treatment by combining mathematical modelling and wet-lab experiments to characterize the cell population dynamics under different PARPi schedules. Using data from in vitro Incucyte Zoom time-lapse microscopy experiments and a step-wise model selection process we derive a calibrated and validated ordinary differential equation model, which we then use to test different plausible adaptive treatment schedules. Our model can accurately predict the in vitro treatment dynamics, even to new schedules, and suggests that treatment modifications need to be carefully timed, or one risks losing control over tumour growth, even in the absence of any resistance. This is because our model predicts that multiple rounds of cell division are required for cells to acquire sufficient DNA damage to induce apoptosis. As a result, adaptive therapy algorithms that modulate treatment but never completely withdraw it are predicted to perform better in this setting than strategies based on treatment interruptions. Pilot experiments in vivo confirm this conclusion. Overall, this study contributes to a better understanding of the impact of scheduling on treatment outcome for PARPis and showcases some of the challenges involved in developing adaptive therapies for new treatment settings.

Mathematical modeling of regulatory networks of intracellular processes - Aims and selected methods

Regulatory networks structure and signaling pathways dynamics are uncovered in time- and resource consuming experimental work. However, it is increasingly supported by modeling, analytical and computational techniques as well as discrete mathematics and artificial intelligence applied to to extract knowledge from existing databases. This review is focused on mathematical modeling used to analyze dynamics and robustness of these networks. This paper presents a review of selected modeling methods that facilitate advances in molecular biology.

Reinforcement learning as an innovative model-based approach: Examples from precision dosing, digital health and computational psychiatry

Model-based approaches are instrumental for successful drug development and use. Anchored within pharmacological principles, through mathematical modeling they contribute to the quantification of drug response variability and enables precision dosing. Reinforcement learning (RL)-a set of computational methods addressing optimization problems as a continuous learning process-shows relevance for precision dosing with high flexibility for dosing rule adaptation and for coping with high dimensional efficacy and/or safety markers, constituting a relevant approach to take advantage of data from digital health technologies. RL can also support contributions to the successful development of digital health applications, recognized as key players of the future healthcare systems, in particular for reducing the burden of non-communicable diseases to society. RL is also pivotal in computational psychiatry-a way to characterize mental dysfunctions in terms of aberrant brain computations-and represents an innovative modeling approach forpsychiatric indications such as depression or substance abuse disorders for which digital therapeutics are foreseen as promising modalities.

Artificial Intelligence and Mechanistic Modeling for Clinical Decision Making in Oncology

The amount of "big" data generated in clinical oncology, whether from molecular, imaging, pharmacological, or biological origin, brings novel challenges. To mine efficiently this source of information, mathematical models able to produce predictive algorithms and simulations are required, with applications for diagnosis, prognosis, drug development, or prediction of the response to therapy. Such mathematical and computational constructs can be subdivided into two broad classes: biologically agnostic, statistical models using artificial intelligence techniques, and physiologically based, mechanistic models. In this review, recent advances in the applications of such methods in clinical oncology are outlined. These include machine learning applied to big data (omics, imaging, or electronic health records), pharmacometrics and quantitative systems pharmacology, as well as tumor kinetics and metastasis modeling. Focus is set on studies with high potential of clinical translation, and particular attention is given to cancer immunotherapy. Perspectives are given in terms of combinations of the two approaches: "mechanistic learning."

The timing of cyclic cytotoxic chemotherapy can worsen neutropenia and neutrophilia

Despite recent advances in immunotherapies, cytotoxic chemotherapy continues to be a first-line treatment option for the majority of cancers. Unfortunately, a common side effect in patients undergoing chemotherapy treatment is neutropenia. To mitigate the risk of neutropenia and febrile neutropenia, prophylactic treatment with granulocyte-colony stimulating factor (G-CSF) is administered. Extensive pharmacokinetic/pharmacodynamic modelling of myelosuppression during chemotherapy has suggested avenues for therapy optimization to mitigate this neutropenia. However, the issue of resonance, whereby neutrophil oscillations are induced by the periodic administration of cytotoxic chemotherapy and the coadministration of G-CSF, potentially aggravating a patient's neutropenic/neutrophilic status, is not well-characterized in the clinical literature. Here, through analysis of neutrophil data from young acute lymphoblastic leukaemia patients, we find that resonance is occurring during cyclic chemotherapy treatment in 26% of these patients. Motivated by these data and our previous modelling studies on adult lymphoma patients, we examined resonance during treatment with or without G-CSF. Using our quantitative systems pharmacology model of granulopoiesis, we show that the timing of cyclic chemotherapy can worsen neutropenia or neutrophilia, and suggest clinically-actionable schedules to reduce the resonant effect. We emphasize that delaying supportive G-CSF therapy to 6-7 days after chemotherapy can mitigate myelosuppressive effects. This study therefore highlights the importance of quantitative systems pharmacology for the clinical practice for developing rational therapeutic strategies.

Optimizing adaptive cancer therapy: dynamic programming and evolutionary game theory

Recent clinical trials have shown that adaptive drug therapies can be more efficient than a standard cancer treatment based on a continuous use of maximum tolerated doses (MTD). The adaptive therapy paradigm is not based on a preset schedule; instead, the doses are administered based on the current state of tumour. But the adaptive treatment policies examined so far have been largely ad hoc. We propose a method for systematically optimizing adaptive policies based on an evolutionary game theory model of cancer dynamics. Given a set of treatment objectives, we use the framework of dynamic programming to find the optimal treatment strategies. In particular, we optimize the total drug usage and time to recovery by solving a Hamilton-Jacobi-Bellman equation. We compare MTD-based treatment strategy with optimal adaptive treatment policies and show that the latter can significantly decrease the total amount of drugs prescribed while also increasing the fraction of initial tumour states from which the recovery is possible. We conclude that the use of optimal control theory to improve adaptive policies is a promising concept in cancer treatment and should be integrated into clinical trial design.

Optimal dynamic regimens with artificial intelligence: The case of temozolomide

We determine an optimal protocol for temozolomide using population variability and dynamic optimization techniques inspired by artificial intelligence. We use a Pharmacokinetics/Pharmacodynamics (PK/PD) model based on Faivre and coauthors (Faivre, et al., 2013) for the pharmacokinetics of temozolomide, as well as the pharmacodynamics of its efficacy. For toxicity, which is measured by the nadir of the normalized absolute neutrophil count, we formalize the myelosuppression effect of temozolomide with the physiological model of Panetta and coauthors (Panetta, et al., 2003). We apply the model to a population with variability as given in Panetta and coauthors (Panetta, et al., 2003). Our optimization algorithm is a variant in the class of Monte-Carlo tree search algorithms. We do not impose periodicity constraint on our solution. We set the objective of tumor size minimization while not allowing more severe toxicity levels than the standard Maximum Tolerated Dose (MTD) regimen. The protocol we propose achieves higher efficacy in the sense that –compared to the usual MTD regimen– it divides the tumor size by approximately 7.66 after 336 days –the 95% confidence interval being [7.36–7.97]. The toxicity is similar to MTD. Overall, our protocol, obtained with a very flexible method, gives significant results for the present case of temozolomide and calls for further research mixing operational research or artificial intelligence and clinical research in oncology.

Modeling therapeutic response to radioiodine in metastatic thyroid cancer: a proof-of-concept study for individualized medicine

Purpose

Radioiodine therapy (RAI) has traditionally been used as treatment for metastatic thyroid cancer, based on its ability to concentrate iodine. Propositions to maximize tumor response with minimizing toxicity, must recognize the infinite possibilities of empirical tests. Therefore, an approach of this study was to build a mathematical model describing tumor growth with the kinetics of thyroglobulin (Tg) concentrations over time, following RAI for metastatic thyroid cancer.

Experimental Design

Data from 50 patients with metastatic papillary thyroid carcinoma treated within eight French institutions, followed over 3 years after initial RAI treatments, were included in the model. A semi-mechanistic mathematical model that describes the tumor growth under RAI treatment was designed.

Results

Our model was able to separate patients who responded to RAI from those who did not, concordant with the physicians' determination of therapeutic response. The estimated tumor doubling-time (Td was found to be the most informative parameter for the distinction between responders and non-responders. The model was also able to reclassify particular patients in early treatment stages.

Conclusions

The results of the model present classification criteria that could indicate whether patients will respond or not to RAI treatment, and provide the opportunity to perform personalized management plans.

Mathematical optimisation of the cisplatin plus etoposide combination for managing extensive-stage small-cell lung cancer patients

Background:

Small-cell lung cancer (SCLC) represents one of the most aggressive forms of lung cancer. Despite the fair sensitivity of SCLC to chemotherapy and radiotherapy, the current standard treatment regimens have modest survival rates and are associated with potential life-threatening adverse events. Therefore, research into new optimised regimens that increase drug efficacy while respecting toxicity constraints is of primary importance.

Methods:

A PK/PD model for the combination of cisplatin and etoposide to treat extensive-stage SCLC patients was generated. The model takes into consideration both the efficacy of the drugs and their haematological toxicity. Using optimisation techniques, the model can be used to propose new regimens.

Results:

Three new regimens with varying timing for combining cisplatin and etoposide have been generated that respect haematological toxicity constraints and achieve better or similar tumour regression. The proposed regimens are: (1) Protocol OP1: etoposide 80 mg m⁻² over 1 h D1, followed by a long infusion 12 h later (over 3 days) of 160 mg m⁻² plus cisplatin 80 mg m⁻² over 1 h D1, D1–D1 21 days; (2) Protocol OP2: etoposide 80 mg m⁻² over 1 h D1, followed by a long infusion 12 h later (over 4 days) of 300 mg m⁻² plus cisplatin 100 mg m⁻² over 1 h D1, D1–D1 21 days; and (3) Protocol OP3: etoposide 40 mg m⁻² over 1 h, followed by a long infusion 6 h later (3 days) of 105 mg m⁻² plus cisplatin 50 mg m⁻² over 1 h, D1–D1 14 days.

Conclusions:

Mathematical modelling can help optimise the design of new cisplatin plus etoposide regimens for managing extensive-stage SCLC patients.