Randomized Phase 3 Trial of the Hypoxia Modifier Nimorazole Added to Radiation Therapy With Benefit Assessed in Hypoxic Head and Neck Cancers Determined Using a Gene Signature (NIMRAD).

PURPOSE: Tumor hypoxia is an adverse prognostic factor in head and neck squamous cell carcinoma (HNSCC). We assessed whether patients with hypoxic HNSCC benefited from the addition of nimorazole to definitive intensity modulated radiation therapy (IMRT). METHODS AND MATERIALS: NIMRAD was a phase 3, multicenter, placebo-controlled, double-anonymized trial of patients with HNSCC unsuitable for concurrent platinum chemotherapy or cetuximab with definitive IMRT (NCT01950689). Patients were randomized 1:1 to receive IMRT (65 Gy in 30 fractions over 6 weeks) plus nimorazole (1.2 g/m2 daily, before IMRT) or placebo. The primary endpoint was freedom from locoregional progression (FFLRP) in patients with hypoxic tumors, defined as greater than or equal to the median tumor hypoxia score of the first 50 patients analyzed (≥0.079), using a validated 26-gene signature. The planned sample size was 340 patients, allowing for signature generation in 85% and an assumed hazard ratio (HR) of 0.50 for nimorazole effectiveness in the hypoxic group and requiring 66 locoregional failures to have 80% power in a 2-tail log-rank test at the 5% significance level. RESULTS: Three hundred thirty-eight patients were randomized by 19 centers in the United Kingdom from May 2014 to May 2019, with a median follow-up of 3.1 years (95% CI, 2.9-3.4). Hypoxia scores were available for 286 (85%). The median patient age was 73 years (range, 44-88; IQR, 70-76). There were 36 (25.9%) locoregional failures in the hypoxic group, in which nimorazole + IMRT did not improve FFLRP (adjusted HR, 0.72; 95% CI, 0.36-1.44; P = .35) or overall survival (adjusted HR, 0.96; 95% CI, 0.53-1.72; P = .88) compared with placebo + IMRT. Similarly, nimorazole + IMRT did not improve FFLRP or overall survival in the whole population. In total (N = 338), 73% of patients allocated nimorazole adhered to the drug for ≥50% of IMRT fractions. Nimorazole + IMRT caused more acute nausea compared with placebo + IMRT (Common Terminology Criteria for Adverse Events version 4.0 G1+2: 56.6% vs 42.4%, G3: 10.1% vs 5.3%, respectively; P < .05). CONCLUSIONS: Addition of the hypoxia modifier nimorazole to IMRT for locally advanced HNSCC in older and less fit patients did not improve locoregional control or survival.

Radius additivity score: a novel combination index for tumour growth inhibition in fixed-dose xenograft studies.

The effect of combination therapies in many cancers has often been shown to be superior to that of monotherapies. This success is commonly attributed to drug synergies. Combinations of two (or more) drugs in xenograft tumor growth inhibition (TGI) studies are typically designed at fixed doses for each compound. The available methods for assessing synergy in such study designs are based on combination indices (CI) and model-based analyses. The former methods are suitable for screening exercises but are difficult to verify in in vivo studies, while the latter incorporate drug synergy in semi-mechanistic frameworks describing disease progression and drug action but are unsuitable for screening. In the current study, we proposed the empirical radius additivity (Rad-add) score, a novel CI for synergy detection in fixed-dose xenograft TGI combination studies. The Rad-add score approximates model-based analysis performed using the semi-mechanistic constant-radius growth TGI model. The Rad-add score was compared with response additivity, defined as the addition of the two response values, and the bliss independence model in combination studies derived from the Novartis PDX dataset. The results showed that the bliss independence and response additivity models predicted synergistic interactions with high and low probabilities, respectively. The Rad-add score predicted synergistic probabilities that appeared to be between those predicted with response additivity and the Bliss model. We believe that the Rad-add score is particularly suitable for assessing synergy in the context of xenograft combination TGI studies, as it combines the advantages of CI approaches suitable for screening exercises with those of semi-mechanistic TGI models based on a mechanistic understanding of tumor growth.

Novel Thienoduocarmycin-Trastuzumab ADC Demonstrates Strong Antitumor Efficacy with Favorable Safety Profile in Preclinical Studies.

New antibodies-drug conjugate (ADC) payloads overcoming chemoresistance and killing also poorly proliferating tumors at well-tolerated doses are much desired. Duocarmycins are a well-known class of highly potent cytotoxic agents, with DNA minor groove-binding and alkylation properties, active also in chemoresistant tumors. Although different duocarmycin derivatives have been used during the years as payloads for ADC production, unfavorable physicochemical properties impaired the production of ADCs with optimal features. Optimization of the toxin to balance reactivity and stability features and best linker selection allowed us to develop the novel duocarmycin-like payload-linker NMS-P945 suitable for conjugation to mAbs with reproducible drug-antibody ratio (DAR) >3.5. When conjugated to trastuzumab, it generated an ADC with good internalization properties, ability to induce bystander effect and immunogenic cell death. Moreover, it showed strong target-driven activity in cells and cytotoxic activity superior to trastuzumab deruxtecan tested, in parallel, in cell lines with HER2 expression. High in vivo efficacy with cured mice at well-tolerated doses in HER2-driven models was also observed. A developed pharmacokinetic/pharmacodynamic (PK/PD) model based on efficacy in mice and cynomolgus monkey PK data, predicted tumor regression in patients upon administration of 2 doses of trastuzumab-NMS-P945-ADC at 0.5 mg/kg. Thus, considering the superior physicochemical features for ADC production and preclinical results obtained with the model trastuzumab ADC, including bystander effect, immunogenic cell death and activity in chemoresistant tumors, NMS-P945 represents a highly effective, innovative payload for the creation of novel, next-generation ADCs.

A two-stages global sensitivity analysis by using the δ sensitivity index in presence of correlated inputs: application on a tumor growth inhibition model based on the dynamic energy budget theory.

Global sensitivity analysis (GSA) evaluates the impact of variability and/or uncertainty of the model parameters on given model outputs. GSA is useful for assessing the quality of Pharmacometric model inference. Indeed, model parameters can be affected by high (estimation) uncertainty due to the sparsity of data. Independence between model parameters is a common assumption of GSA methods. However, ignoring (known) correlations between parameters may alter model predictions and, then, GSA results. To address this issue, a novel two-stages GSA technique based on the δ index, which is well-defined also in presence of correlated parameters, is here proposed. In the first step, statistical dependencies are neglected to identify parameters exerting causal effects. Correlations are introduced in the second step to consider the real distribution of the model output and investigate also the 'indirect' effects due to the correlation structure. The proposed two-stages GSA strategy was applied, as case study, to a preclinical tumor-in-host-growth inhibition model based on the Dynamic Energy Budget theory. The aim is to evaluate the impact of the model parameter estimate uncertainty (including correlations) on key model-derived metrics: the drug threshold concentration for tumor eradication, the tumor volume doubling time and a new index evaluating the drug efficacy-toxicity trade-off. This approach allowed to rank parameters according to their impact on the output, discerning whether a parameter mainly exerts a causal or 'indirect' effect. Thus, it was possible to identify uncertainties that should be necessarily reduced to obtain robust predictions for the outputs of interest.

mvLognCorrEst: an R package for sampling from multivariate lognormal distributions and estimating correlations from uncomplete correlation matrix.

BACKGROUND AND OBJECTIVE: Pharmacometrics (PMX) is a quantitative discipline which supports decision-making processes in all stages of drug development. PMX leverages Modeling and Simulations (M&S), which represents a powerful tool to characterize and predict the behavior and the effect of a drug. M&S-based methods, such as Sensitivity Analysis (SA) and Global Sensitivity Analysis (GSA), are gaining interest in PMX as they allow the evaluation of model-informed inference quality. Simulations should be correctly designed to obtain reliable results. Neglecting correlations between model parameters can significantly alter the results of simulations. However, the introduction of a correlation structure between model parameters can cause some issues. Sampling from a multivariate lognormal distribution, which is the typically distribution assumed for PMX model parameters, is not straightforward when a correlation structure is introduced. Indeed, correlations need to respect some constraints which depend by the CVs (i.e., coefficients of variation) of lognormal variables. In addition, when correlation matrices have some unspecified values, they should be properly fixed preserving the positive semi-definiteness of the correlation structure. In this paper, we present mvLognCorrEst, an R package developed to address these issues. METHODS: The proposed sampling strategy was based on reconducting the extraction from the multivariate lognormal distribution of interest to the underlying Normal distribution. However, with high lognormal CVs, a positive semi-definite Normal covariance matrix cannot be obtained due to the violation of some theoretical constraints. In these cases, the Normal covariance matrix was approximated to its nearest positive definite matrix using Frobenius norm as matrix distance. For the estimation of unknown correlations terms, the graph theory was used to represent the correlation structure as weighed undirected graph. Plausible value ranges for the unspecified correlations were derived considering the paths between variables. Then, their estimation was performed by solving a constrained optimization problem. RESULTS: Package functions are presented and applied on a real case study, that is the GSA of a PMX model that has been recently developed to support preclinical oncological studies. CONCLUSIONS: mvLognCorrEst package is an R tool to support simulation-based analysis for which sampling from multivariate lognormal distributions with correlated variables and/or estimation of partially defined correlation matrix are required.

Use of In Vivo Imaging and Physiologically-Based Kinetic Modelling to Predict Hepatic Transporter Mediated Drug-Drug Interactions in Rats.

Gadoxetate, a magnetic resonance imaging (MRI) contrast agent, is a substrate of organic-anion-transporting polypeptide 1B1 and multidrug resistance-associated protein 2. Six drugs, with varying degrees of transporter inhibition, were used to assess gadoxetate dynamic contrast enhanced MRI biomarkers for transporter inhibition in rats. Prospective prediction of changes in gadoxetate systemic and liver AUC (AUCR), resulting from transporter modulation, were performed by physiologically-based pharmacokinetic (PBPK) modelling. A tracer-kinetic model was used to estimate rate constants for hepatic uptake (khe), and biliary excretion (kbh). The observed median fold-decreases in gadoxetate liver AUC were 3.8- and 1.5-fold for ciclosporin and rifampicin, respectively. Ketoconazole unexpectedly decreased systemic and liver gadoxetate AUCs; the remaining drugs investigated (asunaprevir, bosentan, and pioglitazone) caused marginal changes. Ciclosporin decreased gadoxetate khe and kbh by 3.78 and 0.09 mL/min/mL, while decreases for rifampicin were 7.20 and 0.07 mL/min/mL, respectively. The relative decrease in khe (e.g., 96% for ciclosporin) was similar to PBPK-predicted inhibition of uptake (97-98%). PBPK modelling correctly predicted changes in gadoxetate systemic AUCR, whereas underprediction of decreases in liver AUCs was evident. The current study illustrates the modelling framework and integration of liver imaging data, PBPK, and tracer-kinetic models for prospective quantification of hepatic transporter-mediated DDI in humans.

Physiologically Based Pharmacokinetic Modeling of Transporter-Mediated Hepatic Disposition of Imaging Biomarker Gadoxetate in Rats.

Physiologically based pharmacokinetic (PBPK) models are increasingly used in drug development to simulate changes in both systemic and tissue exposures that arise as a result of changes in enzyme and/or transporter activity. Verification of these model-based simulations of tissue exposure is challenging in the case of transporter-mediated drug-drug interactions (tDDI), in particular as these may lead to differential effects on substrate exposure in plasma and tissues/organs of interest. Gadoxetate, a promising magnetic resonance imaging (MRI) contrast agent, is a substrate of organic-anion-transporting polypeptide 1B1 (OATP1B1) and multidrug resistance-associated protein 2 (MRP2). In this study, we developed a gadoxetate PBPK model and explored the use of liver-imaging data to achieve and refine in vitro-in vivo extrapolation (IVIVE) of gadoxetate hepatic transporter kinetic data. In addition, PBPK modeling was used to investigate gadoxetate hepatic tDDI with rifampicin i.v. 10 mg/kg. In vivo dynamic contrast-enhanced (DCE) MRI data of gadoxetate in rat blood, spleen, and liver were used in this analysis. Gadoxetate in vitro uptake kinetic data were generated in plated rat hepatocytes. Mean (%CV) in vitro hepatocyte uptake unbound Michaelis-Menten constant (Km,u) of gadoxetate was 106 µM (17%) (n = 4 rats), and active saturable uptake accounted for 94% of total uptake into hepatocytes. PBPK-IVIVE of these data (bottom-up approach) captured reasonably systemic exposure, but underestimated the in vivo gadoxetate DCE-MRI profiles and elimination from the liver. Therefore, in vivo rat DCE-MRI liver data were subsequently used to refine gadoxetate transporter kinetic parameters in the PBPK model (top-down approach). Active uptake into the hepatocytes refined by the liver-imaging data was one order of magnitude higher than the one predicted by the IVIVE approach. Finally, the PBPK model was fitted to the gadoxetate DCE-MRI data (blood, spleen, and liver) obtained with and without coadministered rifampicin. Rifampicin was estimated to inhibit active uptake transport of gadoxetate into the liver by 96%. The current analysis highlighted the importance of gadoxetate liver data for PBPK model refinement, which was not feasible when using the blood data alone, as is common in PBPK modeling applications. The results of our study demonstrate the utility of organ-imaging data in evaluating and refining PBPK transporter IVIVE to support the subsequent model use for quantitative evaluation of hepatic tDDI.

A latent variable approach to account for correlated inputs in global sensitivity analysis.

In drug development decision-making is often supported through model-based methods, such as physiologically-based pharmacokinetics (PBPK). Global sensitivity analysis (GSA) is gaining use for quality assessment of model-informed inference. However, the inclusion and interpretation of correlated factors in GSA has proven an issue. Here we developed and evaluated a latent variable approach for dealing with correlated factors in GSA. An approach was developed that describes the correlation between two model inputs through the causal relationship of three independent factors: the latent variable and the unique variances of the two correlated parameters. The latent variable approach was applied to a set of algebraic models and a case from PBPK. Then, this method was compared to Sobol's GSA assuming no correlations, Sobol's GSA with groups and the Kucherenko approach. For the latent variable approach, GSA was performed with Sobol's method. By using the latent variable approach, it is possible to devise a unique and easy interpretation of the sensitivity indices while maintaining the correlation between the factors. Compared methods either consider the parameters independent, group the dependent variables into one unique factor or present difficulties in the interpretation of the sensitivity indices. In situations where GSA is called upon to support model-informed decision-making, the latent variable approach offers a practical method, in terms of ease of implementation and interpretability, for applying GSA to models with correlated inputs that does not violate the independence assumption. Prerequisites and limitations of the approach are discussed.

Inter-compound and Intra-compound Global Sensitivity Analysis of a Physiological Model for Pulmonary Absorption of Inhaled Compounds.

In recent years, global sensitivity analysis (GSA) has gained interest in physiologically based pharmacokinetics (PBPK) modelling and simulation from pharmaceutical industry, regulatory authorities, and academia. With the case study of an in-house PBPK model for inhaled compounds in rats, the aim of this work is to show how GSA can contribute in PBPK model development and daily use. We identified two types of GSA that differ in the aims and, thus, in the parameter variability: inter-compound and intra-compound GSA. The inter-compound GSA aims to understand which are the parameters that mostly influence the variability of the metrics of interest in the whole space of the drugs' properties, and thus, it is useful during the model development. On the other hand, the intra-compound GSA aims to highlight how much the uncertainty associated with the parameters of a given drug impacts the uncertainty in the model prediction and so, it is useful during routine PBPK use. In this work, inter-compound GSA highlighted that dissolution- and formulation-related parameters were mostly important for the prediction of the fraction absorbed, while the permeability is the most important parameter for lung AUC and MRT. Intra-compound GSA highlighted that, for all the considered compounds, the permeability was one of the most important parameters for lung AUC, MRT and plasma MRT, while the extraction ratio and the dose for the plasma AUC. GSA is a crucial instrument for the quality assessment of model-based inference; for this reason, we suggest its use during both PBPK model development and use.

Mechanistic Multiscale Pharmacokinetic Model for the Anticancer Drug 2',2'-difluorodeoxycytidine (Gemcitabine) in Pancreatic Cancer.

The aim of this work is to build a mechanistic multiscale pharmacokinetic model for the anticancer drug 2',2'-difluorodeoxycytidine (gemcitabine, dFdC), able to describe the concentrations of dFdC metabolites in the pancreatic tumor tissue in dependence of physiological and genetic patient characteristics, and, more in general, to explore the capabilities and limitations of this kind of modeling strategy. A mechanistic model characterizing dFdC metabolic pathway (metabolic network) was developed using in vitro literature data from two pancreatic cancer cell lines. The network was able to describe the time course of extracellular and intracellular dFdC metabolites concentrations. Moreover, a physiologically-based pharmacokinetic model was developed to describe clinical dFdC profiles by using enzymatic and physiological information available in the literature. This model was then coupled with the metabolic network to describe the dFdC active metabolite profile in the pancreatic tumor tissue. Finally, global sensitivity analysis was performed to identify the parameters that mainly drive the interindividual variability for the area under the curve (AUC) of dFdC in plasma and of its active metabolite (dFdCTP) in tumor tissue. From this analysis, cytidine deaminase (CDA) concentration was identified as the main driver of plasma dFdC AUC interindividual variability, whereas CDA and deoxycytidine kinase concentration mainly explained the tumor dFdCTP AUC variability. However, the lack of in vitro and in vivo information needed to characterize key model parameters hampers the development of this kind of mechanistic approach. Further studies to better characterize pancreatic cell lines and patient enzymes polymorphisms are encouraged to refine and validate the current model.

Accounting for inter-correlation between enzyme abundance: a simulation study to assess implications on global sensitivity analysis within physiologically-based pharmacokinetics.

Physiologically based pharmacokinetic (PBPK) models often include several sets of correlated parameters, such as organ volumes and blood flows. Because of recent advances in proteomics, it has been demonstrated that correlations are also present between abundances of drug-metabolising enzymes in the liver. As the focus of population PBPK has shifted the emphasis from the average individual to theoretically conceivable extremes, reliable estimation of the extreme cases has become paramount. We performed a simulation study to assess the impact of the correlation between the abundances of two enzymes on the pharmacokinetics of drugs that are substrate of both, under assumptions of presence or lack of such correlations. We considered three semi-physiological models representing the cases of: (1) intravenously administered drugs metabolised by two enzymes expressed in the liver; (2) orally administered drugs metabolised by CYP3A4 expressed in the liver and gut wall; (3) intravenously administered drugs that are substrates of CYP3A4 and OATP1B1 in the liver. Finally, the impact of considering or ignoring correlation between enzymatic abundances on global sensitivity analysis (GSA) was investigated using variance based GSA on a reduced PBPK model for repaglinide, substrate of CYP3A4 and CYP2C8. Implementing such correlations can increase the confidence interval for population pharmacokinetic parameters (e.g., AUC, bioavailability) and impact the GSA results. Ignoring these correlations could lead to the generation of implausible parameters combinations and to an incorrect estimation of pharmacokinetic related parameters. Thus, known correlations should always be considered in building population PBPK models.

Variance based global sensitivity analysis of physiologically based pharmacokinetic absorption models for BCS I-IV drugs.

Regulatory agencies have a strong interest in sensitivity analysis for the evaluation of physiologically-based pharmacokinetic (PBPK) models used in pharmaceutical research and drug development and regulatory submissions. One of the applications of PBPK is the prediction of fraction absorbed and bioavailability for drugs following oral administration. In this context, we performed a variance based global sensitivity analysis (GSA) on in-house PBPK models for drug absorption, with the aim of identifying key parameters that influence the predictions of the fraction absorbed and the bioavailability for neutral, acidic and basic compounds. This analysis was done for four different classes of drugs, defined according to the Biopharmaceutics Classification System, differentiating compounds by permeability and solubility. For class I compounds (highly permeable, highly soluble), the parameters that mainly influence the fraction absorbed are related to the formulation properties, for class II compounds (highly permeable, lowly soluble) to the dissolution process, for class III (lowly permeable, highly soluble) to both absorption process and formulation properties and for class IV (lowly permeable, lowly soluble) to both absorption and dissolution processes. Considering the bioavailability, the results are similar to those for the fraction absorbed, with the addition that parameters related to gut wall and liver clearance influence as well the predictions. This work aimed to give a demonstration of the GSA methodology and highlight its importance in improving our understanding of PBPK absorption models and in guiding the choice of parameters that can safely be assumed, estimated or require data generation to allow informed model prediction.