Population pharmacokinetic-pharmacodynamic modelling in oncology: a tool for predicting clinical response

Model-informed precision dosing: State of the art and future perspectives

Model-informed precision dosing (MIPD) stands as a significant development in personalized medicine to tailor drug dosing to individual patient characteristics. MIPD moves beyond traditional therapeutic drug monitoring (TDM) by integrating mathematical predictions of dosing, and considering patient-specific factors (patient characteristics, drug measurements) as well as different sources of variability. For this purpose, rigorous model qualification is required for the application of MIPD in patients. This review delves into new methods in model selection and validation, also highlighting the role of machine learning in improving MIPD, the utilization of biosensors for real-time monitoring, as well as the potential of models integrating biomarkers for efficacy or toxicity for precision dosing. The clinical evidence of TDM and MIPD is discussed for various medical fields including infection medicine, oncology, transplant medicine, and inflammatory bowel diseases, thereby underscoring the role of pharmacokinetics/pharmacodynamics and specific biomarkers. Further research, particularly randomized clinical trials, is warranted to corroborate the value of MIPD in enhancing patient outcomes and advancing personalized medicine.

Disentangling Anti-Tumor Response of Immunotherapy Combinations: A Physiologically Based Framework for V937 Oncolytic Virus and Pembrolizumab

Immuno-oncology (IO) is a growing strategy in cancer treatment. Oncolytic viruses (OVs) can selectively infect cancer cells and lead to direct and/or immune-dependent tumor lysis. This approach represents an opportunity to potentiate the efficacy of immune checkpoint inhibitors (ICI), such as pembrolizumab. Currently, there is a lack of comprehensive quantitative models for the aforementioned scenarios. In this work, we developed a mechanistic framework describing viral kinetics, viral dynamics, and tumor response after intratumoral (i.t.) or intravenous (i.v.) administration of V937 alone or in combination with pembrolizumab. The model accounts for tumor shrinkage, in both injected and non-injected lesions, induced by: viral-infected tumor cell death and activated CD8 cells. OV-infected tumor cells enhanced the expansion of CD8 cells, whereas pembrolizumab inhibits their exhaustion by competing with PD-L1 in their binding to PD-1. Circulating viral levels and treatment effects on tumor volume were adequately characterized in all the different scenarios. This mechanistic-based model has been developed by combining top-down and bottom-up approaches and provides individual estimates of viral and ICI responses. The robustness of the model is reflected by the description of the tumor size time profiles in a variety of clinical scenarios. Additionally, this platform allows us to investigate not only the contribution of processes related to the viral kinetics and dynamics on tumor response, but also the influence of its interaction with an ICI. Additionally, the model can be used to explore different scenarios aiming to optimize treatment combinations and support clinical development.

Quantification of biochemical PSA dynamics after radioligand therapy with [177Lu]Lu-PSMA-I&T using a population pharmacokinetic/pharmacodynamic model

BACKGROUND

There is an unmet need for prediction of treatment outcome or patient selection for [177Lu]Lu-PSMA therapy in patients with metastatic castration-resistant prostate cancer (mCRPC). Quantification of the tumor exposure-response relationship is pivotal for further treatment optimization. Therefore, a population pharmacokinetic (PK) model was developed for [177Lu]Lu-PSMA-I&T using SPECT/CT data and, subsequently, related to prostate-specific antigen (PSA) dynamics after therapy in patients with mCRPC using a pharmacokinetic/pharmacodynamic (PKPD) modelling approach.

METHODS

A population PK model was developed using quantitative SPECT/CT data (406 scans) of 76 patients who received multiple cycles [177Lu]Lu-PSMA-I&T (± 7.4 GBq with either two- or six-week interval). The PK model consisted of five compartments; central, salivary glands, kidneys, tumors and combined remaining tissues. Covariates (tumor volume, renal function and cycle number) were tested to explain inter-individual variability on uptake into organs and tumors. The final PK model was expanded with a PD compartment (sequential fitting approach) representing PSA dynamics during and after treatment. To explore the presence of a exposure-response relationship, individually estimated [177Lu]Lu-PSMA-I&T tumor concentrations were related to PSA changes over time.

RESULTS

The population PK model adequately described observed data in all compartments (based on visual inspection of goodness-of-fit plots) with adequate precision of parameters estimates (< 36.1% relative standard error (RSE)). A significant declining uptake in tumors (k14) during later cycles was identified (uptake decreased to 73%, 50% and 44% in cycle 2, 3 and 4-7, respectively, compared to cycle 1). Tumor growth (defined by PSA increase) was described with an exponential growth rate (0.000408 h-1 (14.2% RSE)). Therapy-induced PSA decrease was related to estimated tumor concentrations (MBq/L) using both a direct and delayed drug effect. The final model adequately captured individual PSA concentrations after treatment (based on goodness-of-fit plots). Simulation based on the final PKPD model showed no evident differences in response for the two different dosing regimens currently used.

CONCLUSIONS

Our population PK model accurately described observed [177Lu]Lu-PSMA-I&T uptake in salivary glands, kidneys and tumors and revealed a clear declining tumor uptake over treatment cycles. The PKPD model adequately captured individual PSA observations and identified population response rates for the two dosing regimens. Hence, a PKPD modelling approach can guide prediction of treatment response and thus identify patients in whom radioligand therapy is likely to fail.

New 2-alkoxycyanopyridine derivatives as inhibitors of EGFR, HER2, and DHFR: Synthesis, anticancer evaluation, and molecular modeling studies

New series of substituted 2-alkoxycyanopyridine derivatives were synthesized and evaluated for their in vitro and in vivo anticancer activities. Comparing the evaluated activities against cancer cell lines to the broad-spectrum anticancer doxorubicin, and the kinase inhibitor sorafenib, compounds 3a, 4b, 4c, 7a, and 8d demonstrated superior anticancer efficacy with elevated safety profiles and selectivity indices, particularly against MCF7 breast cancer. For exploration of their mechanism of action, assays for inhibition of EGFR, HER2 kinase, and DHFR were performed. The promising synthesized compounds exhibited potent dual kinase EGFR/HER2 inhibitory activity with IC50values of 0.248/0.156 μM for 4b and 0.138/0.092 μM for 4c. Additionally, with IC50 values of 0.138 and 0.193 M, respectively, 4b and 4c had the greatest DHFR inhibitory activity that was comparable to methotrexate. In the MCF7 breast cancer cell line, they caused arrest at the S phase of the cell cycle and exhibited apoptosis induction activity. With restored caspase-3 immunoexpression, the anti-breast cancer assay performed in vivo of 4b and 4c demonstrated a substantial decrease in tumor volume. Results from molecular modeling were in agreement with biological assays proving the importance of the 3-caynopyridine, two substituted phenyl rings attached to central pyridine ring, and propoxy side chain moieties for binding with the receptors. As 4c works by inhibiting both EGFR/HER2 kinase, DHFR enzymes, in addition to cellular apoptosis, it could be viewed as a model of compounds possessing a multi-targeting anticancer activity. Collectively, compounds 4b and 4c might represent prototypes for further development as anticancer molecules.

C-Reactive Protein as an Early Predictor of Efficacy in Advanced Non-Small-Cell Lung Cancer Patients: A Tumor Dynamics-Biomarker Modeling Framework

In oncology, longitudinal biomarkers reflecting the patient's status and disease evolution can offer reliable predictions of the patient's response to treatment and prognosis. By leveraging clinical data in patients with advanced non-small-cell lung cancer receiving first-line chemotherapy, we aimed to develop a framework combining anticancer drug exposure, tumor dynamics (RECIST criteria), and C-reactive protein (CRP) concentrations, using nonlinear mixed-effects models, to evaluate and quantify by means of parametric time-to-event models the significance of early longitudinal predictors of progression-free survival (PFS) and overall survival (OS). Tumor dynamics was characterized by a tumor size (TS) model accounting for anticancer drug exposure and development of drug resistance. CRP concentrations over time were characterized by a turnover model. An x-fold change in TS from baseline linearly affected CRP production. CRP concentration at treatment cycle 3 (day 42) and the difference between CRP concentration at treatment cycles 3 and 2 were the strongest predictors of PFS and OS. Measuring longitudinal CRP allows for the monitoring of inflammatory levels and, along with its reduction across treatment cycles, presents a promising prognostic marker. This framework could be applied to other treatment modalities such as immunotherapies or targeted therapies allowing the timely identification of patients at risk of early progression and/or short survival to spare them unnecessary toxicities and provide alternative treatment decisions.

Population pharmacokinetic modeling of treosulfan and rationale for dose recommendation in children treated for conditioning prior to allogeneic hematopoietic stem cell transplantation

Intravenously infused treosulfan was evaluated in adult and pediatric patients for conditioning regimen prior to allogeneic hematopoietic stem cell transplantation. A population pharmacokinetic (PK) model was initially developed on 116 adult and pediatric PK profiles from historical trials, to support treosulfan dose recommendations for children in 2 prospective trials. The aim was to assess and update the initial population PK model by inclusion of additional 83 pediatric PK profiles from these 2 trials. The final population PK model was 2-compartmental with dosing in the central compartment, linear elimination, and inter-compartmental clearance. Inter-individual variability was included on clearance (CL), central volume (V1), peripheral volume (V2), and inter-compartmental clearance (Q). The final model described an effect of the body surface area (BSA) on CL, V1, V2, and Q. The final model resulted in a modified dose recommendation for children and advises treosulfan doses of 10 g/m2, 12 g/m2, and 14 g/m2 for BSAs of <0.4 m2, ≥0.4 to <0.9 m2, and ≥0.9 m2, respectively. This simplified BSA-dependent dose recommendation was developed for children, ensuring a well comparable treosulfan exposure as a dose of 14 g/m2 in adults - irrespective of their age and without applying individual therapeutic drug monitoring.

External validation of a tumor growth inhibition-overall survival model in non-small-cell lung cancer based on atezolizumab studies using alectinib data

BACKGROUND

A modeling framework was previously developed to simulate overall survival (OS) using tumor growth inhibition (TGI) data from six randomized phase 2/3 atezolizumab monotherapy or combination studies in non-small-cell lung cancer (NSCLC). We aimed to externally validate this framework to simulate OS in patients with treatment-naive advanced anaplastic lymphoma kinase (ALK)-positive NSCLC in the alectinib ALEX study.

METHODS

TGI metrics were estimated from a biexponential model using longitudinal tumor size data from a Phase 3 study evaluating alectinib compared with crizotinib in patients with treatment-naive ALK-positive advanced NSCLC. Baseline prognostic factors and TGI metric estimates were used to predict OS.

RESULTS

286 patients were evaluable (at least baseline and one post-baseline tumor size measurements) out of 303 (94%) followed for up to 5 years (cut-off: 29 November 2019). The tumor growth rate estimate and baseline prognostic factors (inflammatory status, tumor burden, Eastern Cooperative Oncology Group performance status, race, line of therapy, and sex) were used to simulate OS in ALEX study. Observed survival distributions for alectinib and crizotinib were within model 95% prediction intervals (PI) for approximately 2 years. Predicted hazard ratio (HR) between alectinib and crizotinib was in agreement with the observed HR (predicted HR 0.612, 95% PI 0.480-0.770 vs. 0.625 observed HR).

CONCLUSION

The TGI-OS model based on unselected or PD-L1 selected NSCLC patients included in atezolizumab trials is externally validated to predict treatment effect (HR) in a biomarker-selected (ALK-positive) population included in alectinib ALEX trial suggesting that TGI-OS models may be treatment independent.

A framework for precision "dosing" of mental healthcare services: algorithm development and clinical pilot

BACKGROUND

One in five adults in the US experience mental illness and over half of these adults do not receive treatment. In addition to the access gap, few innovations have been reported for ensuring the right level of mental healthcare service is available at the right time for individual patients.

METHODS

Historical observational clinical data was leveraged from a virtual healthcare system. We conceptualize mental healthcare services themselves as therapeutic interventions and develop a prototype computational framework to estimate their potential longitudinal impacts on depressive symptom severity, which is then used to assess new treatment schedules and delivered to clinicians via a dashboard. We operationally define this process as "session dosing": 497 patients who started treatment with severe symptoms of depression between November 2020 and October 2021 were used for modeling. Subsequently, 22 mental health providers participated in a 5-week clinical quality improvement (QI) pilot, where they utilized the prototype dashboard in treatment planning with 126 patients.

RESULTS

The developed framework was able to resolve patient symptom fluctuations from their treatment schedules: 77% of the modeling dataset fit criteria for using the individual fits for subsequent clinical planning where five anecdotal profile types were identified that presented different clinical opportunities. Based on initial quality thresholds for model fits, 88% of those individuals were identified as adequate for session optimization planning using the developed dashboard, while 12% supported more thorough treatment planning (e.g. different treatment modalities). In the clinical pilot, 90% of clinicians reported using the dashboard a few times or more per member. Although most clinicians (67.5%) either rarely or never used the dashboard to change session types, numerous other discussions were enabled, and opportunities for automating session recommendations were identified.

CONCLUSIONS

It is possible to model and identify the extent to which mental healthcare services can resolve depressive symptom severity fluctuations. Implementation of one such prototype framework in a real-world clinic represents an advancement in mental healthcare treatment planning; however, investigations to assess which clinical endpoints are impacted by this technology, and the best way to incorporate such frameworks into clinical workflows, are needed and are actively being pursued.

Semiautomated Pipeline to Quantify Tumor Evolution From Real-World Positron Emission Tomography/Computed Tomography Imaging

PURPOSE

A semiautomated pipeline for the collection and curation of free-text and imaging real-world data (RWD) was developed to quantify cancer treatment outcomes in large-scale retrospective real-world studies. The objectives of this article are to illustrate the challenges of RWD extraction, to demonstrate approaches for quality assurance, and to showcase the potential of RWD for precision oncology.

METHODS

We collected data from patients with advanced melanoma receiving immune checkpoint inhibitors at the Lausanne University Hospital. Cohort selection relied on semantically annotated electronic health records and was validated using process mining. The selected imaging examinations were segmented using an automatic commercial software prototype. A postprocessing algorithm enabled longitudinal lesion identification across imaging time points and consensus malignancy status prediction. Resulting data quality was evaluated against expert-annotated ground-truth and clinical outcomes obtained from radiology reports.

RESULTS

The cohort included 108 patients with melanoma and 465 imaging examinations (median, 3; range, 1-15 per patient). Process mining was used to assess clinical data quality and revealed the diversity of care pathways encountered in a real-world setting. Longitudinal postprocessing greatly improved the consistency of image-derived data compared with single time point segmentation results (classification precision increased from 53% to 86%). Image-derived progression-free survival resulting from postprocessing was comparable with the manually curated clinical reference (median survival of 286 v 336 days, P = .89).

CONCLUSION

We presented a general pipeline for the collection and curation of text- and image-based RWD, together with specific strategies to improve reliability. We showed that the resulting disease progression measures match reference clinical assessments at the cohort level, indicating that this strategy has the potential to unlock large amounts of actionable retrospective real-world evidence from clinical records.

Bayesian forecasting of tumor size metrics and overall survival

The tumor size ratio (TSR), time-to-tumor growth (TTG), and tumor growth rate (kG) are frequently suggested as model-based predictors of overall survival (OS) for different types of tumors. When the tumor metrics are applied in forecasting of the outcome for individual patients at an early stage, the tumor data might be sparse resulting in imprecise prediction. This simulation study aimed to investigate how the tumor follow-up data and estimation approaches influence the accuracy in the tumor size metrics and the predicted hazard of death for individual patients. Longitudinal tumor size and OS data were simulated using tumor growth inhibition and Weibull distribution models, respectively. Based on the model and increasing measurement durations, the accuracy (defined as 80-125% of the simulated "true" value) in individual metrics and hazard was computed. TSR week 6 (TSRw6) accuracy was adequate for 91% of the patients when tumor size was measured up to 12 weeks. For TTG and kG metrics, the highest accuracy observed was lower (43 and 77%, respectively) and occurred later (42 and 60 weeks, respectively). The simultaneous (joint) and sequential estimation approaches resulted in similar accuracies, however, in general, the sequential approach where individual tumor size parameters are fixed, demonstrated inferior estimation properties. The TSRw6 and the model-predicted tumor time course (absolute or relative change) had better forecasting properties than TTG or kG. The population pharmacokinetic (PK) parameters and data approach performed similarly or better than the simultaneous approach and had a better accuracy in estimating individuals' hazard of death than the individual PK parameters method.

Impact of Cabozantinib Exposure on Proteinuria and Muscle Toxicity in Patients with Unresectable Hepatocellular Carcinoma

This prospective study investigated the impact of cabozantinib exposure on proteinuria and muscle toxicity, in a cohort of 14 Japanese patients with unresectable hepatocellular carcinoma (uHCC). We measured the trough concentration of cabozantinib (Ctrough) weekly for 6 weeks after starting treatment. Although the initial dose was less than 60 mg in most cases, dose interruption occurred in 79%, primarily because of proteinuria and/or malaise. The median and coefficient of variation of maximum Ctrough at 7−42 d were 929.0 ng/mL and 59.2%, respectively. The urinary protein-to-creatinine ratio (UPCR), serum creatine kinase, and serum aldolase values were all significantly elevated following treatment. Moreover, maximum changes in serum creatine kinase and aldolase were significantly associated with maximum Ctrough (r = 0.736, p < 0.01; r = 0.798, p < 0.001; respectively). Receiver operating characteristic (ROC) curve analysis showed that changes in serum creatine kinase ≥70.5 U/L and aldolase ≥6.1 U/L from baseline relatively accurately predicted inclusion in the high-maximum Ctrough (≥929.0 ng/mL) group, with an area under the ROC of 0.929 and 0.833, respectively. Measurement of serum creatine kinase and aldolase may increase the clinical usefulness of cabozantinib treatment for uHCC and help alleviate difficulties with dose adjustments.

Longitudinal and time-to-event modeling for prognostic implications of radical surgery in retroperitoneal sarcoma

Retroperitoneal sarcoma (RPS) is a rare malignancy which can be difficult to manage due to the variety of clinical behaviors. In this study, we aimed to develop a parametric modeling framework to quantify the relationship between postoperative dynamics of several biomarkers and overall/progression-free survival of RPS. One hundred seventy-four patients with RPS who received surgical resection with curative intent at the Peking University Cancer Hospital Sarcoma Center were retrospectively included. Potential prognostic factors were preliminarily identified. Longitudinal analyses of body mass index (BMI), serum total protein (TP), and white blood cells (WBCs) were performed using nonlinear mixed effects models. The impacts of time-varying and time-invariant predictors on survival were investigated by parametric time-to-event (TTE) models. The majority of patients experienced decline in BMI, recovery of TP, as well as transient elevation in WBC counts after surgery, which significantly correlated with survival. An indirect-response model incorporating surgery effect described the fluctuation in percentage BMI. The recovery of TP was captured by a modified Gompertz model, and a semimechanistic model was selected for WBCs. TTE models estimated that the daily cumulative average of predicted BMI and WBC, the seventh-day TP, as well as certain baseline variables, were significant predictors of survival. Model-based simulations were performed to examine the clinical significance of prognostic factors. The current work quantified the individual trajectories of prognostic biomarkers in response to surgery and predicted clinical outcomes, which would constitute an additional strategy for disease monitoring and intervention in postoperative RPS.

Longitudinal nonlinear mixed effects modeling of EGFR mutations in ctDNA as predictor of disease progression in treatment of EGFR-mutant non-small cell lung cancer

Correlations between increasing concentrations of circulating tumor DNA (ctDNA) in plasma and disease progression have been shown. A nonlinear mixed effects model to describe the dynamics of epidermal growth factor receptor (EGFR) ctDNA data from patients with non-small cell lung cancer (NSCLC) combined with a parametric survival model were developed to evaluate the ability of these modeling techniques to describe ctDNA data. Repeated ctDNA measurements on L858R, exon19del, and T790M mutants were available from 54 patients with EGFR mutated NSCLC treated with erlotinib or gefitinib. Different dynamic models were tested to describe the longitudinal ctDNA concentrations of the driver and resistance mutations. Subsequently, a parametric time-to-event model for progression-free survival (PFS) was developed. Predicted L858R, exon19del, and T790M concentrations were used to evaluate their value as predictor for disease progression. The ctDNA dynamics were best described by a model consisting of a zero-order increase and first-order elimination (19.7/day, 95% confidence interval [CI] 14.9-23.6/day) of ctDNA concentrations. In addition, time-dependent development of resistance (5.0 × 10^-4 , 95% CI 2.0 × 10^-4 -7.0 × 10^-4 /day) was included in the final model. Relative change in L858R and exon19del concentrations from baseline was identified as most significant predictor of disease progression (p = 0.001). The dynamic model for L858R, exon19del, and T790M concentrations in ctDNA and time-to-event model adequately described the observed concentrations and PFS data in our clinical cohort. In addition, it was shown that nonlinear mixed effects modeling is a valuable method for the analysis of longitudinal and heterogeneous biomarker datasets obtained from clinical practice.

Anti-cancer treatment schedule optimization based on tumor dynamics modelling incorporating evolving resistance

Quantitative characterization of evolving tumor resistance under targeted treatment could help identify novel treatment schedules, which may improve the outcome of anti-cancer treatment. In this study, a mathematical model which considers various clonal populations and evolving treatment resistance was developed. With parameter values fitted to the data or informed by literature data, the model could capture previously reported tumor burden dynamics and mutant KRAS levels in circulating tumor DNA (ctDNA) of patients with metastatic colorectal cancer treated with panitumumab. Treatment schedules, including a continuous schedule, intermittent schedules incorporating treatment holidays, and adaptive schedules guided by ctDNA measurements were evaluated using simulations. Compared with the continuous regimen, the simulated intermittent regimen which consisted of 8-week treatment and 4-week suspension prolonged median progression-free survival (PFS) of the simulated population from 36 to 44 weeks. The median time period in which the tumor size stayed below the baseline level (T_TS<TS0) was prolonged from 52 to 60 weeks. Extending the treatment holiday resulted in inferior outcomes. The simulated adaptive regimens showed to further prolong median PFS to 56–64 weeks and T_TS<TS0 to 114–132 weeks under different treatment designs. A prospective clinical study is required to validate the results and to confirm the added value of the suggested schedules.

Impact of genetic polymorphisms on the pharmacokinetics and pharmacodynamics of lenvatinib in patients with hepatocellular carcinoma

This prospective study examined the impact of genetic polymorphisms on the pharmacokinetics and clinical efficacy and safety of lenvatinib, a substrate of ATP-binding transporters, in a cohort of 48 Japanese patients with hepatocellular carcinoma (HCC). Pharmacokinetic studies were performed at the start of lenvatinib therapy (day 1) and on day 15. The coefficients of variation in AUC_0-24h of lenvatinib on days 1 and 15 were 44.0% and 52.4%, respectively. Although the ABCB1 3435C > A, 1236C > T, and 2677G>T/A polymorphisms did not influence pharmacokinetic parameters, the AUC_0-24h values on days 1 and 15 of the ABCG2 C/A or A/A group were approximately 1.1-fold and 1.4-fold that in the ABCG2 C/C group (P = 0.164 and 0.024). There were no significant differences in AUC_0-24h on days 1 and 15 between the responders (complete or partial response) and non-responders (stable or progressive disease). The AUC_0-24h on day 15 in those developing anorexia of any grade was significantly higher than that without such development (P = 0.017). In multivariate analysis, ABCG2 421C > A C/A or A/A was significantly associated with the development of anorexia (odds ratio 9.009, P = 0.009). ABCG2 421C > A polymorphism could affect exposure to lenvatinib and the development of anorexia.

Multistate model for pharmacometric analyses of overall survival in HER2-negative breast cancer patients treated with docetaxel

The aim of this study was to develop a multistate model for overall survival (OS) analysis, based on parametric hazard functions and combined with an investigation of predictors derived from a longitudinal tumor size model on the transition hazards. Different states - stable disease, tumor response, progression, second-line treatment, and death following docetaxel treatment initiation (stable state) in patients with HER2-negative breast cancer (n = 183) were used in model building. Past changes in tumor size prospectively predicts the probability of state changes. The hazard of death after progression was lower for subjects who had longer treatment response (i.e., longer time-to-progression). Young age increased the probability of receiving second-line treatment. The developed multistate model adequately described the transitions between different states and jointly the overall event and survival data. The multistate model allows for simultaneous estimation of transition rates along with their tumor model derived metrics. The metrics were evaluated in a prospective manner so not to cause immortal time bias. Investigation of predictors and characterization of the time to develop response, the duration of response, the progression-free survival, and the OS can be performed in a single multistate modeling exercise. This modeling approach can be applied to other cancer types and therapies to provide a better understanding of efficacy of drug and characterizing different states, thereby facilitating early clinical interventions to improve anticancer therapy.

Exposure-response analysis of efficacy and safety for pexidartinib in patients with tenosynovial giant cell tumor

This analysis was conducted to assess exposure–response relationships for efficacy and safety of pexidartinib in patients with tenosynovial giant cell tumor. Efficacy was assessed categorically by overall response rate (ORR) with Response Evaluation Criteria in Solid Tumors version 1.1 and longitudinally (changes in tumor size and volume). Safety included hepatic parameters (i.e., alanine aminotransferase [ALT], aspartate aminotransferase [AST], and total bilirubin). Average pexidartinib concentration (C_avg) was identified as the primary exposure parameter correlated with response. In categorical and longitudinal analyses, higher C_avg coincided with greater ORR and tumor size reduction, respectively, with smaller joint size having a greater impact. For safety, a significant relationship was observed between C_avg and incidence of ALT‐related and AST‐related adverse events (AEs). With increased exposure, an increase in efficacy was predicted with near maximum effect at 800 mg/day. Higher initial dose (1000 mg/day) during the first 2 weeks did not improve efficacy. Higher doses were associated with an increased risk of ALT‐related and AST‐related AEs. These results support the US Food and Drug Administration–approved dose (400 mg two times/day without initial loading dose).

The Role of Mathematical Models in Immuno-Oncology: Challenges and Future Perspectives

Immuno-oncology (IO) focuses on the ability of the immune system to detect and eliminate cancer cells. Since the approval of the first immune checkpoint inhibitor, immunotherapies have become a major player in oncology treatment and, in 2021, represented the highest number of approved drugs in the field. In spite of this, there is still a fraction of patients that do not respond to these therapies and develop resistance mechanisms. In this sense, mathematical models offer an opportunity to identify predictive biomarkers, optimal dosing schedules and rational combinations to maximize clinical response. This work aims to outline the main therapeutic targets in IO and to provide a description of the different mathematical approaches (top-down, middle-out, and bottom-up) integrating the cancer immunity cycle with immunotherapeutic agents in clinical scenarios. Among the different strategies, middle-out models, which combine both theoretical and evidence-based description of tumor growth and immunological cell-type dynamics, represent an optimal framework to evaluate new IO strategies.

Drug Exposure to Establish Pharmacokinetic-Response Relationships in Oncology

In the oncology field, understanding the relationship between the dose administered and the exerted effect is particularly important because of the narrow therapeutic index associated with anti-cancer drugs and the high interpatient variability. Therefore, in this review, we provide a critical perspective of the different methods of characterising treatment exposure in the oncology setting. The increasing number of modelling applications in oncology reflects the applicability and the impact of pharmacometrics on all phases of the drug development process and patient management as well. Pharmacometric modelling is a worthy component within the current paradigm of model-based drug development, but pharmacometric modelling techniques are also accessible for the clinician in the optimisation of current oncology therapies. Consequently, the application of population models in a hospital setting by generating close collaborations between physicians and pharmacometricians is highly recommended, providing a systematic means of developing and assessing model-based metrics as 'drivers' for various responses to treatments, which can then be evaluated as predictors for treatment success. Characterising the key determinants of variability in exposure is of particular importance for anticancer agents, as efficacy and toxicity are associated with exposure. We present the different strategies to describe and predict drug exposure that can be applied depending on the data available, with the objective of obtaining the most useful information in the patients' favour throughout the full drug cycle. Therefore, the objective of the present article is to review the different approaches used to characterise a patient's exposure to oncology drugs, which will result in a better understanding of the time course of the response and the magnitude of interpatient variability.

Tumor growth inhibition modeling of individual lesion dynamics and interorgan variability in HER2-negative breast cancer patients treated with docetaxel

Information on individual lesion dynamics and organ location are often ignored in pharmacometric modeling analyses of tumor response. Typically, the sum of their longest diameters is utilized. Herein, a tumor growth inhibition model was developed for describing the individual lesion time-course data from 183 patients with metastatic HER2-negative breast cancer receiving docetaxel. The interindividual variability (IIV), interlesion variability (ILV), and interorgan variability of parameters describing the lesion time-courses were evaluated. Additionally, a model describing the probability of new lesion appearance and a time-to-event model for overall survival (OS), were developed. Before treatment initiation, the lesions were largest in the soft tissues and smallest in the lungs, and associated with a significant IIV and ILV. The tumor growth rate was 2.6 times higher in the breasts and liver, compared with other metastatic sites. The docetaxel drug effect in the liver, breasts, and soft tissues was greater than or equal to 1.2 times higher compared with other organs. The time-course of the largest lesion, the presence of at least 3 liver lesions, and the time since study enrollment, increased the probability of new lesion appearance. New lesion appearance, along with the time to growth and time-course of the largest lesion at baseline, were identified as the best predictors of OS. This tumor modeling approach, incorporating individual lesion dynamics, provided a more complete understanding of heterogeneity in tumor growth and drug effect in different organs. Thus, there may be potential to tailor treatments based on lesion location, lesion size, and early lesion response to provide better clinical outcomes.

Identification of high-dimensional omics-derived predictors for tumor growth dynamics using machine learning and pharmacometric modeling

Pharmacometric modeling can capture tumor growth inhibition (TGI) dynamics and variability. These approaches do not usually consider covariates in high-dimensional settings, whereas high-dimensional molecular profiling technologies ("omics") are being increasingly considered for prediction of anticancer drug treatment response. Machine learning (ML) approaches have been applied to identify high-dimensional omics predictors for treatment outcome. Here, we aimed to combine TGI modeling and ML approaches for two distinct aims: omics-based prediction of tumor growth profiles and identification of pathways associated with treatment response and resistance. We propose a two-step approach combining ML using least absolute shrinkage and selection operator (LASSO) regression with pharmacometric modeling. We demonstrate our workflow using a previously published dataset consisting of 4706 tumor growth profiles of patient-derived xenograft (PDX) models treated with a variety of mono- and combination regimens. Pharmacometric TGI models were fit to the tumor growth profiles. The obtained empirical Bayes estimates-derived TGI parameter values were regressed using the LASSO on high-dimensional genomic copy number variation data, which contained over 20,000 variables. The predictive model was able to decrease median prediction error by 4% as compared with a model without any genomic information. A total of 74 pathways were identified as related to treatment response or resistance development by LASSO, of which part was verified by literature. In conclusion, we demonstrate how the combined use of ML and pharmacometric modeling can be used to gain pharmacological understanding in genomic factors driving variation in treatment response.

Characterizing Exposure-Response Relationship for Therapeutic Monoclonal Antibodies in Immuno-Oncology and Beyond: Challenges, Perspectives, and Prospects

Recent data from immuno-oncology clinical studies have shown the exposure-response (E-R) relationship for therapeutic monoclonal antibodies (mAbs) was often confounded by various factors due to the complex interplay of patient characteristics, disease, drug exposure, clearance, and treatment response and presented challenges in characterization and interpretation of E-R analysis. To tackle the challenges, exposure relationships for therapeutic mAbs in immuno-oncology and oncology are reviewed, and a general framework for an integrative understanding of E-R relationship is proposed. In this framework, baseline factors, drug exposure, and treatment response are envisioned to form an interconnected triangle, driving the E-R relationship and underlying three components that compose the apparent relationship: exposure-driven E-R, baseline-driven E-R, and response-driven E-R. Various strategies in data analysis and study design to decouple those components and mitigate the confounding effect are reviewed for their merits and limitations, and a potential roadmap for selection of these strategies is proposed. Specifically, exposure metrics based on a single-dose pharmacokinetic model can be used to mitigate response-driven E-R, while multivariable analysis and/or case control analysis of data obtained from multiple dose levels in a randomized study may be used to account for the baseline-driven E-R. In this context, the importance of collecting data from multiple dose levels, the role of prognostic factors and predictive factors, the potential utility of clearance at baseline and its change over time, and future directions are discussed.

Comparing Circulating Tumor Cell Counts with Dynamic Tumor Size Changes as Predictor of Overall Survival: A Quantitative Modeling Framework

PURPOSE

Quantitative relationships between treatment-induced changes in tumor size and circulating tumor cell (CTC) counts, and their links to overall survival (OS), are lacking. We present a population modeling framework identifying and quantifying such relationships, based on longitudinal data collected in patients with metastatic colorectal cancer (mCRC) to evaluate the value of tumor size and CTC counts as predictors of OS.

EXPERIMENTAL DESIGN

A pharmacometric approach (i.e., population pharmacodynamic modeling) was used to characterize the changes in tumor size and CTC count and evaluate them as predictors of OS in 451 patients with mCRC treated with chemotherapy and targeted therapy in a prospectively randomized phase III study (CAIRO2).

RESULTS

A tumor size model of tumor quiescence and drug resistance was used to characterize the tumor size time-course, and was, in addition to the total normalized dose (i.e., of all administered drugs) in a given cycle, related to the CTC counts through a negative binomial model (CTC model). Tumor size changes did not contribute additional predictive value when the mean CTC count was a predictor of OS. Treatment reduced the typical mean count from 1.43 to 0.477 (HR = 3.94). The modeling framework was applied to explore whether dose modifications (increased and reduced) would result in a CTC count below 1/7.5 mL after 1 to 2 weeks of treatment.

CONCLUSIONS

Time-varying CTC counts can be useful for early predicting OS in patients with mCRC, and may therefore have potential for model-based treatment individualization. Although tumor size was connected to CTC, its link to OS was weaker.

Mathematical Approach to Differentiate Spontaneous and Induced Evolution to Drug Resistance During Cancer Treatment

Purpose

Drug resistance is a major impediment to the success of cancer treatment. Resistance is typically thought to arise from random genetic mutations, after which mutated cells expand via Darwinian selection. However, recent experimental evidence suggests that progression to drug resistance need not occur randomly, but instead may be induced by the treatment itself via either genetic changes or epigenetic alterations. This relatively novel notion of resistance complicates the already challenging task of designing effective treatment protocols.

Materials and Methods

To better understand resistance, we have developed a mathematical modeling framework that incorporates both spontaneous and drug-induced resistance.

Results

Our model demonstrates that the ability of a drug to induce resistance can result in qualitatively different responses to the same drug dose and delivery schedule. We have also proven that the induction parameter in our model is theoretically identifiable and propose an in vitro protocol that could be used to determine a treatment’s propensity to induce resistance.

Pharmacogenetics in Model-Based Optimization of Bevacizumab Therapy for Metastatic Colorectal Cancer

Vascular endothelial growth factor A (VEGF-A) and intercellular adhesion molecule 1 (ICAM-1) are significant regulators of angiogenesis, an important biological process involved in carcinogenesis. Bevacizumab, an anti-VEGF monoclonal antibody (MAB), is approved for the treatment of metastatic Colorectal cancer (mCRC), however clinical outcomes are highly variable. In the present study, we developed a pharmacokinetic (PK), a simplified quasi-steady state (QSS) and a pharmacokinetic/pharmacodynamic (PK/PD) model to identify potential sources of variability. A total of 46 mCRC patients, who received bevacizumab in combination with chemotherapy were studied. VEGF-A (rs2010963, rs1570360, rs699947) and ICAM-1 (rs5498, rs1799969) genes’ polymorphisms, age, gender, weight, and dosing scheme were investigated as possible co-variates of the model’s parameters. Polymorphisms, trough, and peak levels of bevacizumab, and free VEGF-A were determined in whole blood and serum. Data were analyzed using nonlinear mixed-effects modeling. The two-compartment PK model showed that clearance (CL) was significantly lower in patients with mutant ICAM-1 rs1799969 (p < 0.0001), inter-compartmental clearance (Q) was significantly higher with mutant VEGF-A rs1570360 (p < 0.0001), and lower in patients with mutant VEGF-A rs699947 (p < 0.0001). The binding QSS model also showed that mutant ICAM-1 rs1799969 was associated with a lower CL (p = 0.0177). Mutant VEGF-A rs699947 was associated with a lower free VEGF-A levels, prior to the next dose (p = 0.000445). The above results were confirmed by the PK/PD model. Findings of the present study indicated that variants of the genes regulating angiogenesis might affect PK and PD characteristics of bevacizumab, possibly influencing the clinical outcomes.

Integrated Data Analysis of Six Clinical Studies Points Toward Model-Informed Precision Dosing of Tamoxifen

Introduction

At tamoxifen standard dosing, ∼20% of breast cancer patients do not reach proposed target endoxifen concentrations >5.97 ng/mL. Thus, better understanding the large interindividual variability in tamoxifen pharmacokinetics (PK) is crucial. By applying non-linear mixed-effects (NLME) modeling to a pooled ‘real-world’ clinical PK database, we aimed to (i) dissect several levels of variability and identify factors predictive for endoxifen exposure and (ii) assess different tamoxifen dosing strategies for their potential to increase the number of patients reaching target endoxifen concentrations.

Methods

Tamoxifen and endoxifen concentrations with genetic and demographic data of 468 breast cancer patients from six reported studies were used to develop a NLME parent-metabolite PK model. Different levels of variability on model parameters or measurements were investigated and the impact of covariates thereupon explored. The model was subsequently applied in a simulation-based comparison of three dosing strategies with increasing degree of dose individualization for a large virtual breast cancer population. Interindividual variability of endoxifen concentrations and the fraction of patients at risk for not reaching target concentrations were assessed for each dosing strategy.

Results and Conclusions

The integrated NLME model enabled to differentiate and quantify four levels of variability (interstudy, interindividual, interoccasion, and intraindividual). Strong influential factors, i.e., CYP2D6 activity score, drug–drug interactions with CYP3A and CYP2D6 inducers/inhibitors and age, were reliably identified, reducing interoccasion variability to <20% CV. Yet, unexplained interindividual variability in endoxifen formation remained large (47.2% CV). Hence, therapeutic drug monitoring seems promising for achieving endoxifen target concentrations. Three tamoxifen dosing strategies [standard dosing (20 mg QD), CYP2D6-guided dosing (20, 40, and 60 mg QD) and individual model-informed precision dosing (MIPD)] using three therapeutic drug monitoring samples (5–120 mg QD) were compared, leveraging the model. The proportion of patients at risk for not reaching target concentrations was 22.2% in standard dosing, 16.0% in CYP2D6-guided dosing and 7.19% in MIPD. While in CYP2D6-guided- and standard dosing interindividual variability in endoxifen concentrations was high (64.0% CV and 68.1% CV, respectively), it was considerably reduced in MIPD (24.0% CV). Hence, MIPD demonstrated to be the most promising strategy for achieving target endoxifen concentrations.

Model-Based Biomarker Selection for Dose Individualization of Tyrosine-Kinase Inhibitors

Tyrosine-kinase inhibitors (TKIs) demonstrate high inter-individual variability with respect to safety and efficacy and would therefore benefit from dose or schedule adjustments. This study investigated the efficacy, safety, and economical aspects of alternative dosing options for sunitinib in gastro-intestinal stromal tumors (GIST) and axitinib in metastatic renal cell carcinoma (mRCC). Dose individualization based on drug concentration, adverse effects, and sVEGFR-3 was explored using a modeling framework connecting pharmacokinetic and pharmacodynamic models, as well as overall survival. Model-based simulations were performed to investigate four different scenarios: (I) the predicted value of high-dose pulsatile schedules to improve clinical outcomes as compared to regular daily dosing, (II) the potential of biomarkers for dose individualizations, such as drug concentrations, toxicity measurements, and the biomarker sVEGFR-3, (III) the cost-effectiveness of biomarker-guided dose-individualizations, and (IV) model-based dosing approaches versus standard sample-based methods to guide dose adjustments in clinical practice. Simulations from the axitinib and sunitinib frameworks suggest that weekly or once every two weeks high-dosing result in lower overall survival in patients with mRCC and GIST, compared to continuous daily dosing. Moreover, sVEGFR-3 appears a safe and cost-effective biomarker to guide dose adjustments and improve overall survival (€36 784.- per QALY). Model-based estimations were for biomarkers in general found to correctly predict dose adjustments similar to or more accurately than single clinical measurements and might therefore guide dose adjustments. A simulation framework represents a rapid and resource saving method to explore various propositions for dose and schedule adjustments of TKIs, while accounting for complicating factors such as circulating biomarker dynamics and inter-or intra-individual variability.

Tumor Time-Course Predicts Overall Survival in Non-Small Cell Lung Cancer Patients Treated with Atezolizumab: Dependency on Follow-Up Time

The large heterogeneity in response to immune checkpoint inhibitors is driving the exploration of predictive biomarkers to identify patients who will respond to such treatment. We extended our previously suggested modeling framework of atezolizumab pharmacokinetics, IL18, and tumor size (TS) dynamics, to also include overall survival (OS). Baseline and model‐derived variables were explored as predictors of OS in 88 patients with non‐small cell lung cancer treated with atezolizumab. To investigate the impact of follow‐up length on the inclusion of predictors of OS, four different censoring strategies were applied. The time‐course of TS change was the most significant predictor in all scenarios, whereas IL18 was not significant. Identified predictors of OS were similar regardless of censoring strategy, although OS was underpredicted when patients were censored 5 months after last dose. The study demonstrated that the tumor‐time course‐OS relationship could be identified based on early phase I data.

In silico trials: Verification, validation and uncertainty quantification of predictive models used in the regulatory evaluation of biomedical products

•

Regulators now consider also evidences produced in silico.

•

We need accepted methods to evaluate the credibility of models.

•

In this paper we describe the use of the ASME V&V-40 technical standard.

•

We also discuss its application to various types of modelling methods.

Historically, the evidences of safety and efficacy that companies provide to regulatory agencies as support to the request for marketing authorization of a new medical product have been produced experimentally, either in vitro or in vivo. More recently, regulatory agencies started receiving and accepting evidences obtained in silico, i.e. through modelling and simulation. However, before any method (experimental or computational) can be acceptable for regulatory submission, the method itself must be considered “qualified” by the regulatory agency. This involves the assessment of the overall “credibility” that such a method has in providing specific evidence for a given regulatory procedure. In this paper, we describe a methodological framework for the credibility assessment of computational models built using mechanistic knowledge of physical and chemical phenomena, in addition to available biological and physiological knowledge; these are sometimes referred to as “biophysical” models. Using guiding examples, we explore the definition of the context of use, the risk analysis for the definition of the acceptability thresholds, and the various steps of a comprehensive verification, validation and uncertainty quantification process, to conclude with considerations on the credibility of a prediction for a specific context of use. While this paper does not provide a guideline for the formal qualification process, which only the regulatory agencies can provide, we expect it to help researchers to better appreciate the extent of scrutiny required, which should be considered early on in the development/use of any (new) in silico evidence.

Progress and Opportunities to Advance Clinical Cancer Therapeutics Using Tumor Dynamic Models

There is a need for new approaches and endpoints in oncology drug development, particularly with the advent of immunotherapies and the multiple drug combinations under investigation. Tumor dynamics modeling, a key component to oncology "model-informed drug development," has shown a growing number of applications and a broader adoption by drug developers and regulatory agencies in the past years to support drug development and approval in a variety of ways. Tumor dynamics modeling is also being investigated in personalized cancer therapy approaches. These models and applications are reviewed and discussed, as well as the limitations and issues open for further investigations. A close collaboration between stakeholders like clinical investigators, statisticians, and pharmacometricians is warranted to advance clinical cancer therapeutics.

In silico mouse study identifies tumour growth kinetics as biomarkers for the outcome of anti-angiogenic treatment

Angiogenesis is a crucial step in tumour progression, as this process allows tumours to recruit new blood vessels and obtain oxygen and nutrients to sustain growth. Therefore, inhibiting angiogenesis remains a viable strategy for cancer therapy. However, anti-angiogenic therapy has not proved to be effective in reducing tumour growth across a wide range of tumours, and no reliable predictive biomarkers have been found to determine the efficacy of anti-angiogenic treatment. Using our previously established computational model of tumour-bearing mice, we sought to determine whether tumour growth kinetic parameters could be used to predict the outcome of anti-angiogenic treatment. A model trained with datasets from six in vivo mice studies was used to generate a randomized in silico tumour-bearing mouse population. We analysed tumour growth in untreated mice (control) and mice treated with an anti-angiogenic agent and determined the Kaplan–Meier survival estimates based on simulated tumour volume data. We found that the ratio between two kinetic parameters, k₀ and k₁, which characterize the tumour's exponential and linear growth rates, as well as k₁ alone, can be used as prognostic biomarkers of the population survival outcome. Our work demonstrates a robust, quantitative approach for identifying tumour growth kinetic parameters as prognostic biomarkers and serves as a template that can be used to identify other biomarkers for anti-angiogenic treatment.

A Review of Mathematical Models for Tumor Dynamics and Treatment Resistance Evolution of Solid Tumors

Increasing knowledge of intertumor heterogeneity, intratumor heterogeneity, and cancer evolution has improved the understanding of anticancer treatment resistance. A better characterization of cancer evolution and subsequent use of this knowledge for personalized treatment would increase the chance to overcome cancer treatment resistance. Model‐based approaches may help achieve this goal. In this review, we comprehensively summarized mathematical models of tumor dynamics for solid tumors and of drug resistance evolution. Models displayed by ordinary differential equations, algebraic equations, and partial differential equations for characterizing tumor burden dynamics are introduced and discussed. As for tumor resistance evolution, stochastic and deterministic models are introduced and discussed. The results may facilitate a novel model‐based analysis on anticancer treatment response and the occurrence of resistance, which incorporates both tumor dynamics and resistance evolution. The opportunities of a model‐based approach as discussed in this review can be of great benefit for future optimizing and personalizing anticancer treatment.

Optimizing drug development in oncology by clinical trial simulation: Why and how?

In therapeutic research, the safety and efficacy of pharmaceutical products are necessarily tested on humans via clinical trials after an extensive and expensive preclinical development period. Methodologies such as computer modeling and clinical trial simulation (CTS) might represent a valuable option to reduce animal and human assays. The relevance of these methods is well recognized in pharmacokinetics and pharmacodynamics from the preclinical phase to postmarketing. However, they are barely used and are poorly regarded for drug approval, despite Food and Drug Administration and European Medicines Agency recommendations. The generalization of CTS could be greatly facilitated by the availability of software for modeling biological systems, by clinical trial studies and hospital databases. Data sharing and data merging raise legal, policy and technical issues that will need to be addressed. Development of future molecules will have to use CTS for faster development and thus enable better patient management. Drug activity modeling coupled with disease modeling, optimal use of medical data and increased computing speed should allow this leap forward. The realization of CTS requires not only bioinformatics tools to allow interconnection and global integration of all clinical data but also a universal legal framework to protect the privacy of every patient. While recognizing that CTS can never replace 'real-life' trials, they should be implemented in future drug development schemes to provide quantitative support for decision-making. This in silico medicine opens the way to the P4 medicine: predictive, preventive, personalized and participatory.

Neoadjuvant therapy for locally advanced gastric cancer patients. A population pharmacodynamic modeling

Background

Perioperative chemotherapy (CT) or neoadjuvant chemoradiotherapy (CRT) in patients with locally advanced gastric (GC) or gastroesophageal junction cancer (GEJC) has been shown to improve survival compared to an exclusive surgical approach. However, most patients retain a poor prognosis due to important relapse rates. Population pharmacokinetic-pharmacodynamic (PK/PD) modeling may allow identifying at risk-patients. We aimed to develop a mechanistic PK/PD model to characterize the relationship between the type of neoadjuvant therapy, histopathologic response and survival times in locally advanced GC and GEJC patients.

Methods

Patients with locally advanced GC and GEJC treated with neoadjuvant CT with or without preoperative CRT were analyzed. Clinical response was assessed by CT-scan and EUS. Pathologic response was defined as a reduction on pTNM stage compared to baseline cTNM. Metastasis development risk and overall survival (OS) were described using the population approach with NONMEM 7.3. Model evaluation was performed through predictive checks.

Results

A low correlation was observed between clinical and pathologic TNM stage for both T (R = 0.32) and N (R = 0.19) categories. A low correlation between clinical and pathologic response was noticed (R = -0.29). The OS model adequately described the observed survival rates. Disease recurrence, cTNM stage ≥3 and linitis plastica absence, were correlated to a higher risk of death.

Conclusion

Our model adequately described clinical response profiles, though pathologic response could not be predicted. Although the risk of disease recurrence and survival were linked, the identification of alternative approaches aimed to tailor therapeutic strategies to the individual patient risk warrants further research.

Use of Machine-Learning Algorithms in Intensified Preoperative Therapy of Pancreatic Cancer to Predict Individual Risk of Relapse

Background: Although surgical resection is the only potentially curative treatment for pancreatic cancer (PC), long-term outcomes of this treatment remain poor. The aim of this study is to describe the feasibility of a neoadjuvant treatment with induction polychemotherapy (IPCT) followed by chemoradiation (CRT) in resectable PC, and to develop a machine-learning algorithm to predict risk of relapse. Methods: Forty patients with resectable PC treated in our institution with IPCT (based on mFOLFOXIRI, GEMOX or GEMOXEL) followed by CRT (50 Gy and concurrent Capecitabine) were retrospectively analyzed. Additionally, clinical, pathological and analytical data were collected in order to perform a 2-year relapse-risk predictive population model using machine-learning techniques. Results: A R0 resection was achieved in 90% of the patients. After a median follow-up of 33.5 months, median progression-free survival (PFS) was 18 months and median overall survival (OS) was 39 months. The 3 and 5-year actuarial PFS were 43.8% and 32.3%, respectively. The 3 and 5-year actuarial OS were 51.5% and 34.8%, respectively. Forty-percent of grade 3-4 IPCT toxicity, and 29.7% of grade 3 CRT toxicity were reported. Considering the use of granulocyte colony-stimulating factors, the number of resected lymph nodes, the presence of perineural invasion and the surgical margin status, a logistic regression algorithm predicted the individual 2-year relapse-risk with an accuracy of 0.71 (95% confidence interval [CI] 0.56–0.84, p = 0.005). The model-predicted outcome matched 64% of the observed outcomes in an external dataset. Conclusion: An intensified multimodal neoadjuvant approach (IPCT + CRT) in resectable PC is feasible, with an encouraging long-term outcome. Machine-learning algorithms might be a useful tool to predict individual risk of relapse. A small sample size and therapy heterogeneity remain as potential limitations.

Population pharmacokinetics and individualized lobaplatin regimen for the treatment of Chinese small cell lung cancer in the elderly

BACKGROUND

Lobaplatin (LBP) is a third-generation platinum compound.

MATERIAL AND METHODS

This prospective study was performed in 7 institutions in 2014-2016. Elderly small cell lung cancer (SCLC) patients (≥65 years old) were divided into 2 groups to receive LBP regimens according to endogenous creatinine clearance rate (Ccr). LBP was administered at 30 and 20 mg/m in groups A (Ccr ≥ 80 ml/min) and B (60 ml/min ≤ Ccr < 80 ml/min), respectively. The primary endpoint was plasma LBP concentrations. Secondary endpoints were safety and efficacy parameters, including progression-free survival (PFS) and overall survival (OS).

RESULTS

One-hundred patients were enrolled. Median PFS and OS in groups A and B were 155 vs170 days and 306 vs 272 days, respectively. The rates of grade III/IV AEs in groups A and B were 60.8% (n = 31) and 51.0% (n = 25), respectively. In population pharmacokinetics, the area under the curve (AUC) value for group B was 39% lower than that of group A. With LBP administration based on body surface area (BSA), AUC differences between individuals were small.

CONCLUSION

With Ccr ≥ 60 ml/min, BSA based administration is necessary. Meanwhile, LBP-based regimens are reliable in treating elderly patients with SCLC.

Design and Conduct Considerations for First-in-Human Trials

A milestone step in translational science to transform basic scientific discoveries into therapeutic applications is the advancement of a drug candidate from preclinical studies to initial human testing. First-in-human (FIH) trials serve as the link to advance new promising drug candidates and are conducted primarily to determine the safe dose range for further clinical development. Cross-functional collaboration is essential to ensure efficient and successful FIH trials. The aim of this publication is to serve as a tutorial for conducting FIH trials for both small molecule and biological drug candidates with topics covering regulatory requirements, preclinical safety testing, study design considerations, safety monitoring, biomarker assessment, and global considerations. An emphasis is placed on FIH trial design considerations, including starting dose selection, study size and population, dose escalation scheme, and implementation of adaptive designs. In light of the recent revision of the European Medicines Agency (EMA) guideline on FIH trials to promote safety and mitigate risk, we also discuss new measures introduced in the guideline that impact FIH trial design.

Impact of relative dose intensity on bone marrow suppression induced by S-1: retrospective observational study

Background

S-1 (a combination of tegafur, gimeracil, and oteracil) is used to treat various cancers. Bone marrow suppression is a dose-limiting toxicity of S-1. The relationship between relative dose intensity (RDI) and bone marrow suppression has not been investigated. Hence, we aimed to elucidate the threshold for RDI to identify bone marrow suppression induced by S-1.

Methods

In this retrospective cohort study, patients who initiated S-1 treatment at Tokyo Women’s Medical University, Medical Center East between June 2015 and June 2017 were included. Bone marrow suppression induced by S-1 was assessed using Common Terminology Criteria for Adverse Events version 4.0. The relationships between grade 3 or higher bone marrow suppression induced by S-1 and RDIs (i.e., 70, 75, and 80%) were investigated using the multivariate Cox proportional hazard model.

Results

We identified 143 patients in this study. The median RDI was 78.8%. Bone marrow suppression induced by S-1 developed in 19 (13.3%) patients. The multivariate Cox proportional hazard model revealed that grade ≥ 2 lymphocytopenia was associated with bone marrow suppression induced by S-1 regardless of the threshold for RDI. In addition, RDI > 75% [hazard ratio (HR) = 1.71, p < 0.05] and RDI > 80% (HR = 1.65, p < 0.05) were associated with bone marrow suppression induced by S-1.

Conclusions

Reduced dose of S-1 still has the risk of developing bone marrow suppression. Clinicians should assess RDI to identify high risk patients with bone marrow suppression induced by S-1.

Exposure-response analysis and simulation of lenvatinib safety and efficacy in patients with radioiodine-refractory differentiated thyroid cancer

Purpose

Once-daily lenvatinib 24 mg is the approved dose for radioiodine-refractory differentiated thyroid cancer. In a phase 3 trial with lenvatinib, the starting dose of 24 mg was associated with a relatively high incidence of adverse events that required dose reductions. We used an exposure–response model to investigate the risk–benefit of different dosing regimens for lenvatinib.

Methods

A population pharmacokinetics/pharmacodynamics modeling analysis was used to simulate the potential benefit of lower starting doses to retain efficacy with improved safety. The seven lenvatinib regimens tested were: 24 mg; and 20 mg, 18 mg, and 14 mg, all with or without up-titration to 24 mg. Exposure–response models for efficacy and safety were created using a 24-week time course.

Results

The approved dose of lenvatinib at 24 mg, predicted the best efficacy. However, the lenvatinib dosing regimens of 14 mg with up-titration or 18 mg without up-titration potentially provides comparable efficacy (objective response rate at 24 weeks) and a better safety profile.

Conclusions

Treatment with lenvatinib at starting doses lower than the approved once-daily 24 mg dose could provide comparable antitumor efficacy and a similar or better safety profile. Based on the results from this modeling and simulation study, a comparator dose of lenvatinib 18 mg without up-titration was selected for evaluation in a clinical trial.

Electronic supplementary material

The online version of this article (10.1007/s00280-018-3687-4) contains supplementary material, which is available to authorized users.

A PK/PD Analysis of Circulating Biomarkers and Their Relationship to Tumor Response in Atezolizumab-Treated non-small Cell Lung Cancer Patients

To assess circulating biomarkers as predictors of antitumor response to atezolizumab (anti-programmed death-ligand 1 (PD-L1), Tecentriq) serum pharmacokinetic (PK) and 95 plasma biomarkers were analyzed in 88 patients with relapsed/refractory non-small cell lung cancer (NSCLC) receiving atezolizumab i.v. q3w (10-20 mg/kg) in the PCD4989g phase I clinical trial. Following exploratory analyses, two plasma biomarkers were chosen for further study and correlation with change in tumor size (the sum of the longest diameter) was assessed in a pharmacokinetic/pharmacodynamic (PK/PD) tumor modeling framework. When longitudinal kinetics of biomarkers and tumor size were modeled, tumor shrinkage was found to significantly correlate with area under the curve (AUC), baseline factors (metastatic sites, liver metastases, and smoking status), and relative change in interleukin (IL)-18 level from baseline at day 21 (RCFB_{IL -18,d21} ). Although AUC was a major predictor of tumor shrinkage, the effect was estimated to dissipate with an average half-life of 80 days, whereas RCFB_{IL -18,d21} seemed relevant to the duration of the response.

Assessing Similarity Among Individual Tumor Size Lesion Dynamics: The CICIL Methodology

Mathematical models of tumor dynamics generally omit information on individual target lesions (iTLs), and consider the most important variable to be the sum of tumor sizes (TS). However, differences in lesion dynamics might be predictive of tumor progression. To exploit this information, we have developed a novel and flexible approach for the non‐parametric analysis of iTLs, which integrates knowledge from signal processing and machine learning. We called this new methodology ClassIfication Clustering of Individual Lesions (CICIL). We used CICIL to assess similarities among the TS dynamics of 3,223 iTLs measured in 1,056 patients with metastatic colorectal cancer treated with cetuximab combined with irinotecan, in two phase II studies. We mainly observed similar dynamics among lesions within the same tumor site classification. In contrast, lesions in anatomic locations with different features showed different dynamics in about 35% of patients. The CICIL methodology has also been implemented in a user‐friendly and efficient Java‐based framework.

Predicting tumour growth and its impact on survival in gemcitabine-treated patients with advanced pancreatic cancer

The aim of this evaluation was to characterize the impact of the tumour size (TS) effects driven by the anticancer drug gemcitabine on overall survival (OS) in patients with advanced pancreatic cancer by building and validating a predictive semi-mechanistic joint TS-OS model. TS and OS data were obtained from one phase II and one phase III study where gemcitabine was administered (1000-1250 mg/kg over 30-60 min i.v infusion) as single agent to patients (n = 285) with advanced pancreatic cancer. Drug exposure, TS and OS were linked using the population approach with NONMEM 7.3. Pancreatic tumour progression was characterized by exponential growth (doubling time = 67 weeks), and tumour response to treatment was described as a function of the weekly area under the gemcitabine triphosphate concentration vs time curve (AUC), including treatment-related resistance development. The typical predicted percentage of tumour growth inhibition with respect to no treatment was 22.3% at the end of 6 chemotherapy cycles. Emerging resistance elicited a 57% decrease in drug effects during the 6th chemotherapy cycle. Predicted TS profile was identified as main prognostic factor of OS, with tumours responders' profiles improving median OS by 30 weeks compared to stable-disease TS profiles. Results of NCT00574275 trial were predicted using this modelling framework, thereby validating the approach as a prediction tool in clinical development. Our analyses show that despite the advanced stage of the disease in this patient population, the modelling framework herein can be used to predict the likelihood of treatment success using early clinical data.

Advanced Methods for Dose and Regimen Finding During Drug Development: Summary of the EMA/EFPIA Workshop on Dose Finding (London 4-5 December 2014)

Inadequate dose selection for confirmatory trials is currently still one of the most challenging issues in drug development, as illustrated by high rates of late‐stage attritions in clinical development and postmarketing commitments required by regulatory institutions. In an effort to shift the current paradigm in dose and regimen selection and highlight the availability and usefulness of well‐established and regulatory‐acceptable methods, the European Medicines Agency (EMA) in collaboration with the European Federation of Pharmaceutical Industries Association (EFPIA) hosted a multistakeholder workshop on dose finding (London 4–5 December 2014). Some methodologies that could constitute a toolkit for drug developers and regulators were presented. These methods are described in the present report: they include five advanced methods for data analysis (empirical regression models, pharmacometrics models, quantitative systems pharmacology models, MCP‐Mod, and model averaging) and three methods for study design optimization (Fisher information matrix (FIM)‐based methods, clinical trial simulations, and adaptive studies). Pairwise comparisons were also discussed during the workshop; however, mostly for historical reasons. This paper discusses the added value and limitations of these methods as well as challenges for their implementation. Some applications in different therapeutic areas are also summarized, in line with the discussions at the workshop. There was agreement at the workshop on the fact that selection of dose for phase III is an estimation problem and should not be addressed via hypothesis testing. Dose selection for phase III trials should be informed by well‐designed dose‐finding studies; however, the specific choice of method(s) will depend on several aspects and it is not possible to recommend a generalized decision tree. There are many valuable methods available, the methods are not mutually exclusive, and they should be used in conjunction to ensure a scientifically rigorous understanding of the dosing rationale.

Exposure-Response Analysis of Necitumumab Efficacy in Squamous Non-Small Cell Lung Cancer Patients

We sought to describe the exposure-response relationship of necitumumab efficacy in squamous non-small cell lung cancer patients and evaluate intrinsic and extrinsic patient descriptors that may guide dosing. SQUIRE was a phase III study comparing necitumumab in combination with gemcitabine and cisplatin vs. gemcitabine and cisplatin alone in 1,014 patients. An integrated model for tumor size dynamics and overall survival was developed, where reduction in tumor size results in a decrease in survival hazard. The change in tumor size was characterized using linear growth and first-order shrinkage. Overall survival was described using a combination of a Weibull function and Gompertz function for the hazard, with dynamic tumor size being a predictor for the hazard. Although body weight resulted in higher clearance and lower exposure, simulations showed that an 800 mg flat dose provided optimal response regardless of body weight.

Model-based prediction of myelosuppression and recovery based on frequent neutrophil monitoring

Purpose

To investigate whether a more frequent monitoring of the absolute neutrophil counts (ANC) during myelosuppressive chemotherapy, together with model-based predictions, can improve therapy management, compared to the limited clinical monitoring typically applied today.

Methods

Daily ANC in chemotherapy-treated cancer patients were simulated from a previously published population model describing docetaxel-induced myelosuppression. The simulated values were used to generate predictions of the individual ANC time-courses, given the myelosuppression model. The accuracy of the predicted ANC was evaluated under a range of conditions with reduced amount of ANC measurements.

Results

The predictions were most accurate when more data were available for generating the predictions and when making short forecasts. The inaccuracy of ANC predictions was highest around nadir, although a high sensitivity (≥90%) was demonstrated to forecast Grade 4 neutropenia before it occurred. The time for a patient to recover to baseline could be well forecasted 6 days (±1 day) before the typical value occurred on day 17.

Conclusions

Daily monitoring of the ANC, together with model-based predictions, could improve anticancer drug treatment by identifying patients at risk for severe neutropenia and predicting when the next cycle could be initiated.

Electronic supplementary material

The online version of this article (doi:10.1007/s00280-017-3366-x) contains supplementary material, which is available to authorized users.

A Pharmacometric Framework for Axitinib Exposure, Efficacy, and Safety in Metastatic Renal Cell Carcinoma Patients

The relationships between exposure, biomarkers (vascular endothelial growth factor (VEGF), soluble VEGF receptors (sVEGFR)-1, -2, -3, and soluble stem cell factor receptor (sKIT)), tumor sum of longest diameters (SLD), diastolic blood pressure (dBP), and overall survival (OS) were investigated in a modeling framework. The dataset included 64 metastatic renal cell carcinoma patients (mRCC) treated with oral axitinib. Biomarker timecourses were described by indirect response (IDR) models where axitinib inhibits sVEGFR-1, -2, and -3 production, and VEGF degradation. No effect was identified on sKIT. A tumor model using sVEGFR-3 dynamics as driver predicted SLD data well. An IDR model, with axitinib exposure stimulating the response, characterized dBP increase. In a time-to-event model the SLD timecourse predicted OS better than exposure, biomarker- or dBP-related metrics. This type of framework can be used to relate pharmacokinetics, efficacy, and safety to long-term clinical outcome in mRCC patients treated with VEGFR inhibitors. (ClinicalTrial.gov identifier NCT00569946.).