Exposure driven dose escalation design with overdose control: Concept and first real life experience in an oncology phase I trial

Background

Methods

Results

Conclusion

•

Formally leverage pharmacokinetic data and modeling in dose escalation studies.

•

Suitable for molecules with potential non-linear pharmacokinetic.

•

Smarter dose escalation and more informative recommended phase 2 dose.

Step-by-step comparison of ordinary differential equation and agent-based approaches to pharmacokinetic-pharmacodynamic models

Mathematical models in oncology aid in the design of drugs and understanding of their mechanisms of action by simulation of drug biodistribution, drug effects, and interaction between tumor and healthy cells. The traditional approach in pharmacometrics is to develop and validate ordinary differential equation models to quantify trends at the population level. In this approach, time‐course of biological measurements is modeled continuously, assuming a homogenous population. Another approach, agent‐based models, focuses on the behavior and fate of biological entities at the individual level, which subsequently could be summarized to reflect the population level. Heterogeneous cell populations and discrete events are simulated, and spatial distribution can be incorporated. In this tutorial, an agent‐based model is presented and compared to an ordinary differential equation model for a tumor efficacy model inhibiting the pERK pathway. We highlight strengths, weaknesses, and opportunities of each approach.

Association of circulating protein biomarkers with clinical outcomes of durvalumab in head and neck squamous cell carcinoma

The potential for durvalumab, a programmed cell death ligand-1 (PD-L1)-blocking monoclonal antibody, to treat head and neck squamous cell carcinoma (HNSCC) is being evaluated in multiple clinical trials. We assessed circulating proteins at baseline to identify potential biomarkers and to understand pathways related to clinical outcomes for durvalumab. Prior to treatment, 66 serum proteins were measured using multiplex immunoassays for 158 durvalumab-treated HNSCC patients in the phase II HAWK and CONDOR trials as a discovery dataset and 209 durvalumab-treated HNSCC patients in the phase III EAGLE trial as a validation dataset. Multivariate Cox modeling of HAWK and CONDOR datasets established that higher baseline concentrations of interleukin-6 (IL-6), C-reactive protein, S100 calcium-binding protein A12, and angiopoietin-2 (ANGPT2) were associated with shorter overall survival (OS), while higher concentrations of osteocalcin correlated with longer OS after durvalumab treatment (p < .05). All five proteins remained significantly correlated with OS after adjusting for baseline clinical factors, with consistent results across clinical efficacy endpoints based on univariate correlation analyses. The validation dataset from the EAGLE trial confirmed the independent association of IL-6 and osteocalcin with OS, and preserved directional trends for the other biomarkers identified in the discovery dataset. Our results demonstrate the important role of immunosuppressive proteins in the resistance of HNSCC to durvalumab treatment. Osteocalcin showed a positive correlation with clinical outcomes, which remains to be further investigated.

Early predictions of response and survival from a tumor dynamics model in patients with recurrent, metastatic head and neck squamous cell carcinoma treated with immunotherapy

We developed and evaluated a method for making early predictions of best overall response (BOR) and overall survival at 6 months (OS6) in patients with cancer treated with immunotherapy. This method combines machine learning with modeling of longitudinal tumor size data. We applied our method to data from durvalumab‐exposed patients with recurrent/metastatic head and neck cancer. A fivefold cross‐validation was used for model selection. Independent trial data, with various degrees of data truncation, were used for model validation. Mean classification error rates (90% confidence intervals [CIs]) from cross‐validation were 5.99% (90% CI 2.98%–7.50%) for BOR and 19.8% (90% CI 15.8%–39.3%) for OS6. During model validation, the area under the receiver operating characteristic curves was preserved for BOR (0.97, 0.97, and 0.94) and OS6 (0.85, 0.84, and 0.82) at 24, 18, and 12 weeks, respectively. These results suggest our method predicts trial outcomes accurately from early data and could be used to aid drug development.

Population Pharmacokinetics of Durvalumab in Cancer Patients and Association With Longitudinal Biomarkers of Disease Status

The objectives of this analysis were to develop a population pharmacokinetics (PK) model of durvalumab, an anti-PD-L1 antibody, and quantify the impact of baseline and time-varying patient/disease characteristics on PK. Pooled data from two studies (1,409 patients providing 7,407 PK samples) were analyzed with nonlinear mixed effects modeling. Durvalumab PK was best described by a two-compartment model with both linear and nonlinear clearances. Three candidate models were evaluated: a time-invariant clearance (CL) model, an empirical time-varying CL model, and a semimechanistic time-varying CL model incorporating longitudinal covariates related to disease status (tumor shrinkage and albumin). The data supported a slight decrease in durvalumab clearance with time and suggested that it may be associated with a decrease in nonspecific protein catabolic rate among cancer patients who benefit from therapy. No covariates were clinically relevant, indicating no need for dose adjustment. Simulations indicated similar overall PK exposures following weight-based and flat-dosing regimens.

Population Modeling of Tumor Kinetics and Overall Survival to Identify Prognostic and Predictive Biomarkers of Efficacy for Durvalumab in Patients With Urothelial Carcinoma

Durvalumab is an anti‐PD‐L1 monoclonal antibody approved for patients with locally advanced or metastatic urothelial carcinoma (UC) that has progressed after platinum‐containing chemotherapy. A population tumor kinetic model, coupled with dropout and survival models, was developed to describe longitudinal tumor size data and predict overall survival in UC patients treated with durvalumab (NCT01693562) and to identify prognostic and predictive biomarkers of clinical outcomes. Model‐based covariate analysis identified liver metastasis as the most influential factor for tumor growth and immune‐cell PD‐L1 expression and baseline tumor burden as predictive factors for tumor killing. Tumor or immune‐cell PD‐L1 expression, liver metastasis, baseline hemoglobin, and albumin levels were identified as significant covariates for overall survival. These model simulations provided further insights into the impact of PD‐L1 cutoff values on treatment outcomes. The modeling framework can be a useful tool to guide patient selection and enrichment strategies for immunotherapies across various cancer indications.

Dose-Exposure-Response Relationship of the Investigational Anti-Interleukin-13 Monoclonal Antibody Tralokinumab in Patients With Severe, Uncontrolled Asthma

Interleukin (IL)-13 is involved in the pathogenesis of some types of asthma. Tralokinumab is a human immunoglobulin G₄ monoclonal antibody that specifically binds to IL-13. Two placebo-controlled phase II studies (phase IIa, NCT00873860 and phase IIb, NCT01402986) have been conducted in which tralokinumab was administered subcutaneously. This investigation aimed to characterize tralokinumab's dose-exposure-response (forced expiratory volume in 1 s (FEV₁ )) relationship in patients with asthma and to predict the most appropriate dose for phase III. An integrated population pharmacokinetic-pharmacodynamic (PK/PD) modeling analysis was required for phase III dose selection, due to differing phase II patient populations, designs, and regimens. Analysis of combined datasets enabled the identification of tralokinumab's dose-exposure-FEV₁ response relationship in patients with asthma. Near-maximal FEV₁ increase was predicted at a dose of 300 mg SC once every 2 weeks (Q2W). This dose was chosen for tralokinumab in the phase III clinical development program for treatment of severe, uncontrolled asthma.

Pharmacokinetics of tralokinumab in adolescents with asthma: implications for future dosing

AIMS

Tralokinumab, an investigational human immunoglobulin G4 monoclonal antibody, potently and specifically neutralizes interleukin-13, a central mediator of asthma. Tralokinumab has shown improvements in clinical endpoints in adults with uncontrolled asthma. The present study explored the pharmacokinetics (PK) and safety of a single tralokinumab dose, and utilized a population PK modelling and simulation approach to evaluate the optimal dosing strategy for adolescents.

METHODS

Adolescent subjects with asthma, using daily controller medication, received a single subcutaneous dose of tralokinumab 300 mg. Safety, immunogenicity and PK data were collected during a 57-day follow-up. A population PK model was developed using data from the present study and prior studies in adults. Simulations were performed to evaluate dose adjustment requirements for adolescents.

RESULTS

Twenty adolescents (12-17 years) were enrolled; all completed the study. No clinically relevant safety findings or antidrug antibodies were detected. PK parameters were similar to those observed in adults. PK modelling showed that body weight was a minor predictor of tralokinumab PK; after incorporating body weight into the PK model, a 15% (nonparametric 95% confidence interval 5%, 26%) lower clearance was found in adolescents compared with adults [173 (151, 209) vs. 204 (191, 229) ml day(-1)]. Simulations showed no therapeutically relevant differences in exposures between adolescent and adult populations, and similar PK profiles for weight-based (4 mg kg(-1)) and fixed (300 mg) fortnightly subcutaneous doses of tralokinumab.

CONCLUSION

Single-dose administration of tralokinumab 300 mg in adolescents was well tolerated, with a PK profile similar to that in adults. Exposure predictions suggest that dose adjustment is not required for adolescents.

Pharmacodynamic modeling of cough responses to capsaicin inhalation calls into question the utility of the C5 end point

BACKGROUND

Inhaled capsaicin elicits cough reproducibly in human subjects and is widely used in the study of cough and antitussive therapies. However, the traditional end points C2 and C5 (the concentrations of capsaicin inducing at least 2 or 5 coughs, respectively) display extensive overlap between health and disease and therefore might not best reflect clinically relevant mechanisms.

OBJECTIVES

We sought to investigate capsaicin dose responses in different disease groups.

METHODS

Two novel capsaicin cough challenges were compared in patients with chronic cough (CC; n = 20), asthmatic patients (n = 18), and healthy volunteers (HVs; n = 20). Increasing doubling doses of capsaicin (0.48-1000 μmol/L, 4 inhalations per dose) were administered in challenge 1, whereas the order of the doses was randomized in challenge 2. A nonlinear mixed-effects model compared dose-response parameters by disease group and sex. Parameters were also correlated with objective cough frequency.

RESULTS

The model classified subjects based on maximum cough response evoked by any concentration of capsaicin (Emax) and the capsaicin dose inducing half-maximal response (ED50). HVs and asthmatic patients were not statistically different for either parameter and therefore combined for analysis (mean ED50, 38.6 μmol/L [relative SE, 28%]; mean Emax, 4.5 coughs [relative SE, 11%]). Compared with HVs/asthmatic patients, patients with CC had lower ED50 values (14.7 μmol/L [relative SE, 28%], P = .008) and higher Emax values (8.6 coughs [relative SE, 11%], P < .0001). Emax values highly correlated with 24-hour cough frequency (r = 0.71, P < .001) and were 37% higher in female compared with male subjects, regardless of disease group (P < .001).

CONCLUSIONS

Nonlinear mixed-effects modeling demonstrates that maximal capsaicin cough responses better discriminate health from disease and predict spontaneous cough frequency and therefore provide important insights into the mechanisms underlying CC.

Evaluation of the nonparametric estimation method in NONMEM VI: application to real data

The aim of the study was to evaluate the nonparametric estimation methods available in NONMEM VI in comparison with the parametric first-order method (FO) and the first-order conditional estimation method (FOCE) when applied to real datasets. Four methods for estimating model parameters and parameter distributions (FO, FOCE, nonparametric preceded by FO (FO-NONP) and nonparametric preceded by FOCE (FOCE-NONP)) were compared for 25 models previously developed using real data and a parametric method. Numerical predictive checks were used to test the appropriateness of each model. Up to 1000 new datasets were simulated from each model and with each method to construct 90% and 50% prediction intervals. The mean absolute error and the mean error of the different outcomes investigated were computed as indicators of imprecision and bias respectively and formal statistical tests were performed. Overall, less imprecision and less bias were observed with nonparametric methods than with parametric methods. Across the 25 models, t-tests revealed that imprecision and bias were significantly lower (P < 0.05) with FOCE-NONP than with FOCE for half of the NPC outcomes investigated. Improvements were even more pronounced with FO-NONP in comparison with FO. In conclusion, when applied to real datasets and evaluated by numerical predictive checks, the nonparametric estimation methods in NONMEM VI performed better than the corresponding parametric methods (FO or FOCE).