Search

Chimeric antigen receptor (CAR)-T cell therapies have shown tremendous results against various hematologic cancers. Prior to cell infusion, a host preconditioning regimen is required to achieve lymphodepletion and improve CAR-T cell pharmacokinetic exposure, leading to greater chances of therapeutic success. To better understand and quantify the impact of the preconditioning regimen, we built a population-based mechanistic pharmacokinetic-pharmacodynamic model describing the complex interplay between lymphodepletion, host immune system, homeostatic cytokines, and pharmacokinetics of UCART19, an allogeneic product developed against CD19+ B cells. Data were collected from a phase I clinical trial in adult relapsed/refractory B-cell acute lymphoblastic leukemia and revealed three different UCART19 temporal patterns: (i) expansion and persistence, (ii) transient expansion with subsequent rapid decline, and (iii) absence of observed expansion. On the basis of translational assumptions, the final model was able to capture this variability through the incorporation of IL-7 kinetics, which are thought to be increased owing to lymphodepletion, and through an elimination of UCART19 by host T cells, which is specific to the allogeneic context. Simulations from the final model recapitulated UCART19 expansion rates in the clinical trial, confirmed the need for alemtuzumab to observe UCART19 expansion (along with fludarabine cyclophosphamide), quantified the importance of allogeneic elimination, and suggested a high impact of multipotent memory T-cell subpopulations on UCART19 expansion and persistence. In addition to supporting the role of host cytokines and lymphocytes in CAR-T cell therapy, such a model could help optimizing the preconditioning regimens in future clinical trials. Significance: A mathematical mechanistic pharmacokinetic/pharmacodynamic model supports and captures quantitatively the beneficial impact of lymphodepleting patients before the infusion of an allogeneic CAR-T cell product. Mediation through IL-7 increase and host T lymphocytes decrease is underlined, and the model can be further used to optimize CAR-T cell therapies lymphodepletion regimen.

Chimeric antigen receptor (CAR)-T cell therapies have shown tremendous results against various hematologic cancers. Prior to cell infusion, a host preconditioning regimen is required to achieve lymphodepletion and improve CAR-T cell pharmacokinetic exposure, leading to greater chances of therapeutic success. To better understand and quantify the impact of the preconditioning regimen, we built a population-based mechanistic pharmacokinetic-pharmacodynamic model describing the complex interplay between lymphodepletion, host immune system, homeostatic cytokines, and pharmacokinetics of UCART19, an allogeneic product developed against CD19+ B cells. Data were collected from a phase I clinical trial in adult relapsed/refractory B-cell acute lymphoblastic leukemia and revealed three different UCART19 temporal patterns: (i) expansion and persistence, (ii) transient expansion with subsequent rapid decline, and (iii) absence of observed expansion. On the basis of translational assumptions, the final model was able to capture this variability through the incorporation of IL-7 kinetics, which are thought to be increased owing to lymphodepletion, and through an elimination of UCART19 by host T cells, which is specific to the allogeneic context. Simulations from the final model recapitulated UCART19 expansion rates in the clinical trial, confirmed the need for alemtuzumab to observe UCART19 expansion (along with fludarabine cyclophosphamide), quantified the importance of allogeneic elimination, and suggested a high impact of multipotent memory T-cell subpopulations on UCART19 expansion and persistence. In addition to supporting the role of host cytokines and lymphocytes in CAR-T cell therapy, such a model could help optimizing the preconditioning regimens in future clinical trials.Significance:A mathematical mechanistic pharmacokinetic/pharmacodynamic model supports and captures quantitatively the beneficial impact of lymphodepleting patients before the infusion of an allogeneic CAR-T cell product. Mediation through IL-7 increase and host T lymphocytes decrease is underlined, and the model can be further used to optimize CAR-T cell therapies lymphodepletion regimen.

Supplementary Data

Purpose: Flotetuzumab (MGD006 or S80880) is a bispecific molecule that recognizes CD3 and CD123 membrane proteins, redirecting T cells to kill CD123-expressing cells for the treatment of acute myeloid leukemia. In this study, we developed a mathematical model to characterize MGD006 exposure–response relationships and to assess the impact of its immunogenicity in cynomolgus monkeys.Experimental Design: Thirty-two animals received multiple escalating doses (100-300-600-1,000 ng/kg/day) via intravenous infusion continuously 4 days a week. The model reflects sequential binding of MGD006 to CD3 and CD123 receptors. Formation of the MGD006/CD3 complex was connected to total T cells undergoing trafficking, whereas the formation of the trimolecular complex results in T-cell activation and clonal expansion. Activated T cells were used to drive the peripheral depletion of CD123-positive cells. Anti-drug antibody development was linked to MGD006 disposition as an elimination pathway. Model validation was tested by predicting the activity of MGD006 in eight monkeys receiving continuous 7-day infusions.Results: MGD006 disposition and total T-cell and CD123-positive cell profiles were well characterized. Anti-drug antibody development led to the suppression of T-cell trafficking but did not systematically abolish CD123-positive cell depletion. Target cell depletion could persist after drug elimination owing to the self-proliferation of activated T cells generated during the first cycles. The model was externally validated with the 7-day infusion dosing schedule.Conclusions: A translational model was developed for MGD006 that features T-cell activation and expansion as a key driver of pharmacologic activity and provides a mechanistic quantitative platform to inform dosing strategies in ongoing clinical studies. Clin Cancer Res; 24(11); 2631–41. ©2018 AACR.

Supplementary figures and tables Figure S1 PK/PD experimental design. This figure details the dose and sampling strategies for the PK/PD endpoints for all animal groups. Figure S2 Detailed structural PK/PD model of MGD006 administrated in cynomolgus monkeys. This figure supplements Fig. 1. Figure S3 Categorization of anti-drug antibody (ADA) development. This figure describes the method used to categorize the ADA effect in three different groups. Table S1 ADA development groups. This table summarizes the different identified groups and the number of animals in each group. Table S2 Parameter definitions of PK/ADA, total T-cell, and CD123+-cell models. Figure S4 /S5 /S6 Time-course of observed and model-fitted plasma MGD006 concentrations, total T-cell counts, and CD123+-cell counts following continuous 4-day infusion in cynomolgus monkeys. This figure supplements Fig. 3. Figure S7 Goodness-of-fit plots for PK, total T-cell, and CD123+-cell models: individual predictions vs. observations and individual weighted residuals vs. time. Figure S8 Prediction-corrected visual predictive checks (PC-VPC) of PK, total T-cell, and CD123+-cell models. Figure S9 /S10 /S11 External validation of the PK, total T-cell, and CD123+-cell models by VPC with data collected in animals following continous 7-day infusions. These figures supplement Fig. 4.

During development of a drug, the requirement of evaluating the proarrhythmic risk and delayed repolarization needs to be fulfilled. Would it be possible to create an alternative to a thorough QT (TQT) study or is there a need to perform a dedicated TQT study? How is an alternative approach generated, what information is available, and which instructions are considered missing today to generate such an approach? This tutorial describes the considerations and path followed to create an early and feasible alternative to a TQT study using experience-based insights from a successful application to the US Food and Drug Administration for GLPG1972, an ADAMTS-5 inhibitor, and discusses the approach used in light of the current guidelines and literature.

MET receptor tyrosine kinase inhibitors (TKIs) can restore sensitivity to gefitinib, a TKI targeting epidermal growth factor receptor (EGFR), and promote apoptosis in non‐small cell lung cancer (NSCLC) models resistant to gefitinib treatment in vitro and in vivo. Several novel MET inhibitors are currently under study in different phases of development. In this work, a novel tumor‐in‐host modeling approach, based on the Dynamic Energy Budget (DEB) theory, was proposed and successfully applied to the context of poly‐targeted combination therapies. The population DEB‐based tumor growth inhibition (TGI) model well‐described the effect of gefitinib and of two MET inhibitors, capmatinib and S49076, on both tumor growth and host body weight when administered alone or in combination in an NSCLC mice model involving the gefitinib‐resistant tumor line HCC827ER1. The introduction of a synergistic effect in the combination DEB‐TGI model allowed to capture gefitinib anticancer activity enhanced by the co‐administered MET inhibitor, providing also a quantitative evaluation of the synergistic drug interaction. The model‐based comparison of the two MET inhibitors highlighted that S49076 exhibited a greater anticancer effect as well as a greater ability in restoring sensitivity to gefitinib than the competitor capmatinib. In summary, the DEB‐based tumor‐in‐host framework proposed here can be applied to routine combination xenograft experiments, providing an assessment of drug interactions and contributing to rank investigated compounds and to select the optimal combinations, based on both tumor and host body weight dynamics. Thus, the combination tumor‐in‐host DEB‐TGI model can be considered a useful tool in the preclinical development and a significant advance toward better characterization of combination therapies.

Clinical assessment of drug-drug interactions (DDIs) in children is not a common practice in drug development. Therefore, physiologically-based pharmacokinetic (PBPK) modeling can be beneficial for informing drug labeling. Using ivabradine and its metabolite (both cytochrome P450 3A4 enzyme (CYP3A4) substrates), the objectives were (i) to scale ivabradine-metabolite adult PBPK/PD to pediatrics, (ii) to predict the DDIs with a strong CYP3A4 inhibitor, and (iii) to compare the sensitivity of children to DDIs using two CYP3A4 hepatic ontogeny functions: Salem and Upreti. A scaled parent-metabolite PBPK/PD model from adults to children satisfactorily predicted pharmacokinetics (PK) and pharmacodynamics (PD) in 74 children (0.5-18 years) regardless of CYP3A4 hepatic ontogeny function applied. However, using the Salem ontogeny, mean predicted parent and metabolite area under the concentration-time curve over 12 hours (AUC12h ) and heart rate change from baseline were 2-fold, 1.5-fold, and 1.4-fold higher in young children (0.5-3 years old) compared with Upreti ontogeny, respectively. Despite these differences, choice of appropriate hepatic CYP3A4 ontogeny was challenging due to sparse PK and PD data. Different sensitivity to ivabradine-ketoconazole DDIs was simulated in young children relative to adults depending on the choice of hepatic CYP3A4 ontogeny. Predicted ivabradine and metabolite AUCDDI /AUCcontrol were 2-fold lower in the youngest children (0.5-1 year old) compared with adults (Salem function). In contrast, the Upreti function predicted comparable ivabradine DDIs across all age groups, although predicted metabolite AUCDDI/ AUCcontrol was 1.3-fold higher between the youngest children and adults. In the case of PD, differences in predicted DDIs were minor across age groups and between both functions. Current work highlights the importance of careful consideration of hepatic CYP3A4 ontogeny function and implications on labeling recommendations in the pediatric population.

Drug development in oncology commonly exploits the tools of molecular biology to gain therapeutic benefit through reprograming of cellular responses. In immuno-oncology (IO) the aim is to direct the patient's own immune system to fight cancer. After remarkable successes of antibodies targeting PD1/PD-L1 and CTLA4 receptors in targeted patient populations, the focus of further development has shifted toward combination therapies. However, the current drug-development approach of exploiting a vast number of possible combination targets and dosing regimens has proven to be challenging and is arguably inefficient. In particular, the unprecedented number of clinical trials testing different combinations may no longer be sustainable by the population of available patients. Further development in IO requires a step change in selection and validation of candidate therapies to decrease development attrition rate and limit the number of clinical trials. Quantitative systems pharmacology (QSP) proposes to tackle this challenge through mechanistic modeling and simulation. Compounds' pharmacokinetics, target binding, and mechanisms of action as well as existing knowledge on the underlying tumor and immune system biology are described by quantitative, dynamic models aiming to predict clinical results for novel combinations. Here, we review the current QSP approaches, the legacy of mathematical models available to quantitative clinical pharmacologists describing interaction between tumor and immune system, and the recent development of IO QSP platform models. We argue that QSP and virtual patients can be integrated as a new tool in existing IO drug development approaches to increase the efficiency and effectiveness of the search for novel combination therapies.

The Bcl-2 inhibitor S 55746 ((S)-N-(4-hydroxyphenyl)-3-(6-(3-(morpholinomethyl)-1,2,3,4-tetrahydroisoquinoline-2-carbonyl)benzo[d][1,3]dioxol-5-yl)-N-phenyl-5,6,7,8-tetrahydroindolizine-1-carboxamide) is able to restore apoptosis functions impaired by tumorigenesis in mice. Data from pharmacokinetic (PK), biomarker, and tumor growth studies in a xenograft mouse model were considered for population modeling. The aim of the modeling exercise was to link the kinetics of the drug to the biomarker and tumor-size time profiles to better understand its dose-effect relationship. The PK, caspase kinetics, and tumor dynamics were successfully characterized by the proposed pharmacokinetic-pharmacodynamic model. The nonlinear plasma PK was best described by a two-compartment disposition model with both saturable absorption and elimination. Caspase was activated above the effective drug-concentration threshold (CTHRE ), at which near-maximal activity was reached. Increasing the dose did not increase the activation but better sustained it. Tumor growth followed a biphasic pattern, with caspase having an all-or-none inhibiting effect, consistent with the bistability property of the caspase pathway. For tumor eradication, the CTHRE in plasma was 2876 ng ml-1, and the relative caspase activity threshold (CaspTHRE) was 46.5. There was a strong relationship between the time spent above these thresholds and tumor growth inhibition. Tumor growth was inhibited by 50% when CaspTHRE was exceeded 13.8% of the time and when CTHRE was exceeded 8.1% of the time per dosing. This semimechanistic approach, based on experimental mice data and in vitro parameters, provides an interesting tool to quantify or simulate antitumor effects and, eventually, to plan phase 1 studies.

International audience; INTRODUCTION:In this paper, we studied the effect over time of agomelatine, an antidepressant drug administered in patient with major depressive disorder, through item response theory (IRT), taking into account a strong placebo effect and missing not at random. We also assessed the informativeness of the HAMD-17 scale's item.MATERIALS AND METHODS:The data includes five phase III clinical trials sponsored by Servier Institute, totalling 1549 patients followed during a maximum of 1 year. At each observation, individual scores for the 17 items of the HAMD scale were recorded. The probability for each score was modelled with IRT. A non-linear mixed effects model was used to describe the evolution of the disease and was coupled with a time to event model to predict dropout. Clinical trial simulations were then used to compare placebo and active treatment. Informativeness of each item was evaluated using the Fisher information theory.RESULTS:The best model combined an IRT model, a longitudinal model for underlying depression which describes the remission and then a possible relapse, and a hazard model for dropout depending on the evolution from baseline. The drug effect was best modelled as an effect on the remission and the relapse phases. The median predicted drop in HAMD between baseline and 6 weeks was 8.8 (90% PI, 8.3-9.2) when on placebo and 13.1 (90% PI, 12.8-13.4) when treated. Nine items were found to be the most informative.CONCLUSION:The IRT framework allowed to characterise the evolution of depression with time and estimate the effect of agomelatine, as well as the link between symptoms and disease.

Quantitative pharmacology (QP) applications in translational medicine, drug-development, and therapeutic use were crowd-sourced by the ASCPT Impact and Influence initiative. Highlighted QP case studies demonstrated faster access to innovative therapies for patients through 1) rational dose selection for pivotal trials; 2) reduced trial-burden for vulnerable populations; or 3) simplified posology. Critical success factors were proactive stakeholder engagement, alignment on the value of model-informed approaches, and utilizing foundational clinical pharmacology understanding of the therapy.

Background Dabigatran etexilate (DABE) has been suggested as a clinical probe for intestinal P-glycoprotein (P-gp)-mediated drug-drug interaction (DDI) studies and, as an alternative to digoxin. Clinical DDI data with various P-gp inhibitors demonstrated a dose-dependent inhibition of P-gp with DABE. The aims of this study were to develop a joint DABE (prodrug)-dabigatran reduced physiologically-based-pharmacokinetic (PBPK) model and to evaluate its ability to predict differences in P-gp DDI magnitude between a microdose and a therapeutic dose of DABE. Methods A joint DABE-dabigatran PBPK model was developed with a mechanistic intestinal model accounting for the regional P-gp distribution in the gastrointestinal tract. Model input parameters were estimated using DABE and dabigatran pharmacokinetic (PK) clinical data obtained after administration of DABE alone or with a strong P-gp inhibitor, itraconazole, and over a wide range of DABE doses (from 375 µg to 400 mg). Subsequently, the model was used to predict extent of DDI with additional P-gp inhibitors and with different DABE doses. Results The reduced DABE-dabigatran PBPK model successfully described plasma concentrations of both prodrug and metabolite following administration of DABE at different dose levels and when co-administered with itraconazole. The model was able to capture the dose dependency in P-gp mediated DDI. Predicted magnitude of itraconazole P-gp DDI was higher at the microdose (predicted vs. observed median fold-increase in AUC+inh/AUCcontrol (min-max) = 5.88 (4.29–7.93) vs. 6.92 (4.96–9.66) ) compared to the therapeutic dose (predicted median fold-increase in AUC+inh/AUCcontrol = 3.48 (2.37–4.84) ). In addition, the reduced DABE-dabigatran PBPK model predicted successfully the extent of DDI with verapamil and clarithromycin as P-gp inhibitors. Model-based simulations of dose staggering predicted the maximum inhibition of P-gp when DABE microdose was concomitantly administered with itraconazole solution; simulations also highlighted dosing intervals required to minimise the DDI risk depending on the DABE dose administered (microdose vs. therapeutic). Conclusions This study provides a modelling framework for the evaluation of P-gp inhibitory potential of new molecular entities using DABE as a clinical probe. Simulations of dose staggering and regional differences in the extent of intestinal P-gp inhibition for DABE microdose and therapeutic dose provide model-based guidance for design of prospective clinical P-gp DDI studies.

Recent work on Model Informed Drug Discovery and Development (MID3) has noted the need for clarity in model description used in quantitative disciplines such as pharmacology and statistics. 1-3 Cur ...

This document was developed to enable greater consistency in the practice, application, and documentation of Model-Informed Drug Discovery and Development (MID3) across the pharmaceutical industry. A collection of "good practice" recommendations are assembled here in order to minimize the heterogeneity in both the quality and content of MID3 implementation and documentation. The three major objectives of this white paper are to: i) inform company decision makers how the strategic integration of MID3 can benefit RD ii) provide MID3 analysts with sufficient material to enhance the planning, rigor, and consistency of the application of MID3; and iii) provide regulatory authorities with substrate to develop MID3 related and/or MID3 enabled guidelines.

We show through a simulation study how the joint analysis of data from phase I and phase II studies enhances the power of pharmacogenetic tests in pharmacokinetic (PK) studies. PK profiles were simulated under different designs along with 176 genetic markers. The null scenarios assumed no genetic effect, while under the alternative scenarios, drug clearance was associated with six genetic markers randomly sampled in each simulated dataset. We compared penalized regression Lasso and stepwise procedures to detect the associations between empirical Bayes estimates of clearance, estimated by nonlinear mixed effects models, and genetic variants. Combining data from phase I and phase II studies, even if sparse, increases the power to identify the associations between genetics and PK due to the larger sample size. Design optimization brings a further improvement, and we highlight a direct relationship between η-shrinkage and loss of genetic signal.

The middle initial in the fifth author's name is incorrect in the original article. "Lena F. Friberg" should be "Lena E. Friberg". The original article was corrected.

Design of phase 1 combination therapy trials is complex compared to single therapy trials. In this work, model-based adaptive optimal design (MBAOD) was exemplified and evaluated for a combination of paclitaxel and a hypothetical new compound in a phase 1 study to determine the best dosing regimen for a phase 2 trial. Neutropenia was assumed as the main toxicity and the dose optimization process targeted a 33% probability of grade 4 neutropenia and maximal efficacy (based on preclinical studies) by changing the dose amount of both drugs and the dosing schedule for the new drug. Different starting conditions (e.g., initial dose), search paths (e.g., maximal change in dose intensity per step), and stopping criteria (e.g., "3 + 3 rule") were explored. The MBAOD approach was successfully implemented allowing the possibility of flexible designs with the modification of doses and dosing schedule throughout the trial. The 3 + 3 rule was shown to be highly conservative (selection of a dosing regimen with at least 90% of the possible maximal efficacy in less than 21% of the cases) but also safer (selection of a toxic design in less than 2% of the cases). Without the 3 + 3 rule, better performance was observed (>67% of selected designs were associated with at least 90% of possible maximal efficacy) while the proportion of DLTs per trial was similar. Overall, MBAOD is a promising tool in the context of dose finding studies of combination treatments and was showed to be flexible enough to be associated with requirements imposed by clinical protocols. Correction in: AAPS JOURNAL, Volume:20, Issue:3, Article Number:55, DOI:10.1208/s12248-018-0218-5

International audience; Purpose - Normalised prediction distribution errors (npde) are used to graphically and statistically evaluate mixed-effect models for continuous responses. In this study, our aim was to extend npde to time-to-event (TTE) models and evaluate their performance.Methods - Let V denote a dataset with censored TTE observations. The null hypothesis (H) is that observations in V can be described by model M. We extended npde to TTE models using imputations to take into account censoring. We then evaluated their performance in terms of type I error and power to detect model misspecifications for TTE data by means of a simulation study with different sample sizes.Results - Type I error was found to be close to the expected 5% significance level for all sample sizes tested. The npde were able to detect misspecifications in the baseline hazard as well as in the link between the longitudinal variable and the survival function. The ability to detect model misspecifications increased as the difference in the shape of the survival function became more apparent. As expected, the power also increased as the sample size increased. Imputing the censored events tended to decrease the percentage of rejections.Conclusions - We have shown that npde can be readily extended to TTE data and that they perform well with an adequate type I error.

The lack of a common exchange format for mathematical models in pharmacometrics has been a long-standing problem. Such a format has the potential to increase productivity and analysis quality, simplify the handling of complex workflows, ensure reproducibility of research, and facilitate the reuse of existing model resources. Pharmacometrics Markup Language (PharmML), currently under development by the Drug Disease Model Resources (DDMoRe) consortium, is intended to become an exchange standard in pharmacometrics by providing means to encode models, trial designs, and modeling steps.

Purpose: Flotetuzumab (MGD006 or S80880) is a bispecific molecule that recognizes CD3 and CD123 membrane proteins, redirecting T cells to kill CD123-expressing cells for the treatment of acute myeloid leukemia. In this study, we developed a mathematical model to characterize MGD006 exposure–response relationships and to assess the impact of its immunogenicity in cynomolgus monkeys. Experimental Design: Thirty-two animals received multiple escalating doses (100-300-600-1,000 ng/kg/day) via intravenous infusion continuously 4 days a week. The model reflects sequential binding of MGD006 to CD3 and CD123 receptors. Formation of the MGD006/CD3 complex was connected to total T cells undergoing trafficking, whereas the formation of the trimolecular complex results in T-cell activation and clonal expansion. Activated T cells were used to drive the peripheral depletion of CD123-positive cells. Anti-drug antibody development was linked to MGD006 disposition as an elimination pathway. Model validation was tested by predicting the activity of MGD006 in eight monkeys receiving continuous 7-day infusions. Results: MGD006 disposition and total T-cell and CD123-positive cell profiles were well characterized. Anti-drug antibody development led to the suppression of T-cell trafficking but did not systematically abolish CD123-positive cell depletion. Target cell depletion could persist after drug elimination owing to the self-proliferation of activated T cells generated during the first cycles. The model was externally validated with the 7-day infusion dosing schedule. Conclusions: A translational model was developed for MGD006 that features T-cell activation and expansion as a key driver of pharmacologic activity and provides a mechanistic quantitative platform to inform dosing strategies in ongoing clinical studies. Clin Cancer Res; 24(11); 2631–41. ©2018 AACR.

A PBPK modelling approach was used to predict organic anion transporter (OAT) mediated drug-drug interactions involving S44121, a substrate and an inhibitor of OAT1 and OAT3. Model predictions were then compared to the results of a clinical DDI study which was carried out to investigate the interaction of S44121 with probenecid, tenofovir and ciprofloxacin. PBPK models were developed and qualified using existing clinical data, and inhibition constants were determined in vitro. The model predictions for S44121 as an OAT inhibitor were similar to the results obtained from the clinical DDI study, with no interaction observed for tenofovir or ciprofloxacin in the presence of S44121. An observed AUC ratio of 2.2 was obtained for S44121 in the presence of probenecid, which was slightly higher than the model predicted AUC ratio of 1.6. A DDI study in the monkey was also carried out for the interaction between S44121 and probenecid, since the monkey has previously been reported to be a good preclinical model for OAT-mediated DDI. However, this study highlighted a species difference in the major route of S44121 elimination between monkey (mainly hepatic metabolism) and human (mainly renal excretion of unchanged drug), rendering a comparison between the two DDI studies difficult. Overall, for S44121 the PBPK modelling approach gave a better prediction of the extent of DDI than the static predictions based on inhibitor Cmax and IC50, therefore this can be considered a potentially valuable tool within drug development.

International audience; Nonlinear mixed-effect models are used increasingly during drug development. For design, an alternative to simulations is based on the Fisher information matrix. Its expression was derived using a first-order approach, was then extended to include covariance and implemented into the R function PFIM. The impact of covariance on standard errors, amount of information, and optimal designs was studied. It was also shown how standard errors can be predicted analytically within the framework of rich individual data without the model. The results were illustrated by applying this extension to the design of a pharmacokinetic study of a drug in pediatric development.

Central nervous system (CNS) drug disposition is dictated by a drug’s physicochemical properties and its ability to permeate physiological barriers. The blood–brain barrier (BBB), blood-cerebrospinal fluid barrier and centrally located drug transporter proteins influence drug disposition within the central nervous system. Attainment of adequate brain-to-plasma and cerebrospinal fluid-to-plasma partitioning is important in determining the efficacy of centrally acting therapeutics. We have developed a physiologically-based pharmacokinetic model of the rat CNS which incorporates brain interstitial fluid (ISF), choroidal epithelial and total cerebrospinal fluid (CSF) compartments and accurately predicts CNS pharmacokinetics. The model yielded reasonable predictions of unbound brain-to-plasma partition ratio (Kpuu,brain) and CSF:plasma ratio (CSF:Plasmau) using a series of in vitro permeability and unbound fraction parameters. When using in vitro permeability data obtained from L-mdr1a cells to estimate rat in vivo permeability, the model successfully predicted, to within 4-fold, Kpuu,brain and CSF:Plasmau for 81.5% of compounds simulated. The model presented allows for simultaneous simulation and analysis of both brain biophase and CSF to accurately predict CNS pharmacokinetics from preclinical drug parameters routinely available during discovery and development pathways.

The purpose of this study was to explore the impact of censoring due to animal sacrifice on parameter estimates and tumor volume calculated from two diameters in larger tumors during tumor growth experiments in preclinical studies. The type of measurement error that can be expected was also investigated. Different scenarios were challenged using the stochastic simulation and estimation process. One thousand datasets were simulated under the design of a typical tumor growth study in xenografted mice, and then, eight approaches were used for parameter estimation with the simulated datasets. The distribution of estimates and simulation-based diagnostics were computed for comparison. The different approaches were robust regarding the choice of residual error and gave equivalent results. However, by not considering missing data induced by sacrificing the animal, parameter estimates were biased and led to false inferences in terms of compound potency; the threshold concentration for tumor eradication when ignoring censoring was 581 ng.ml(-1), but the true value was 240 ng.ml(-1).

The purpose of this study was to explore the interval censoring induced by caliper measurements on smaller tumors during tumor growth experiments in preclinical studies and to show its impact on parameter estimations. A new approach, the so-called interval-M3 method, is proposed to specifically handle this type of data. Thereby, the interval-M3 method was challenged with different methods (including classical methods for handling below quantification limit values) using Stochastic Simulation and Estimation process to take into account the censoring. In this way, 1000 datasets were simulated under the design of a typical of tumor growth study in xenografted mice, and then, each method was used for parameter estimation on the simulated datasets. Relative bias and relative root mean square error (relative RMSE) were consequently computed for comparison purpose. By not considering the censoring, parameter estimations appeared to be biased and particularly the cytotoxic effect parameter, k 2 , which is the parameter of interest to characterize the efficacy of a compound in oncology. The best performance was noted with the interval-M3 method which properly takes into account the interval censoring induced by caliper measurement, giving overall unbiased estimations for all parameters and especially for the antitumor effect parameter (relative bias = 0.49%, and relative RMSE = 4.06%).

International audience; Genetic data is now collected in many clinical trials, especially in population pharmacokinetic studies. There is no consensus on methods to test the association between pharmacokinetics and genetic covariates. We performed a simulation study inspired by real clinical trials, using the PK of a compound under development having a nonlinear bioavailability along with genotypes for 176 Single Nucleotide Polymorphisms (SNPs). Scenarios included 78 subjects extensively sampled (16 observations per subject) to simulate a phase I study, or 384 subjects with the same rich design. Under the alternative hypothesis (H 1), 6 SNPs were drawn randomly to affect the log-clearance under an additive linear model. For each scenario 200 PK data sets were simulated under the null hypothesis (no gene effect) and H1. We compared 16 combinations of four association tests, a stepwise procedure and three penalised regressions (ridge regression, Lasso, HyperLasso), applied to four pharmacokinetic phenotypes, two observed concentrations, area under the curve estimated by noncompartmental analysis and model-based clearance. The different combinations were compared in terms of true and false positives and probability to detect the genetic effects. In presence of nonlinearity and/or variability in bioavailability, model-based phenotype allowed a higher probability to detect the SNPs than other phenotypes. In a realistic setting with a limited number of subjects, all methods showed a low ability to detect genetic effects. Ridge regression had the best probability to detect SNPs, but also a higher number of false positives. No association test showed a much higher power than the others.

The aim of the present study was to investigate amoxicillin (AMX) distribution in muscle interstitial fluid by microdialysis in healthy, awake rats. Microdialysis probes were inserted into the jugular vein and hind leg muscle. Probe recoveries in each rat were determined by retrodialysis with cefadroxil. AMX was administered as a bolus dose of 50 mg · kg −1 , and microdialysis samples were collected during 180 min. Concentrations of unbound drug in blood and muscle were analyzed simultaneously by a population approach. Simulations were conducted using a hybrid, physiologically based pharmacokinetic model to investigate the potential impact of tissue blood flow on muscle AMX distribution. A two-compartment pharmacokinetic model described adequately the unbound amoxicillin concentration-time profiles in blood and muscle. Muscle AMX distribution equilibrium was rapidly achieved. Consequently, the best results were obtained by considering concentrations in muscle as part of the central compartment. The ratio of the concentration of unbound drug in muscle to that in blood ( R model ) was estimated to 0.80 by the model, which is close to the mean value obtained by noncompartmental data analysis ( R area = 0.86 ± 0.29). Simulations conducted with a hybrid, physiologically based pharmacokinetic model suggest that a muscle blood flow reduction of 30% to 50%, such as could be encountered in critical care patients, has virtually no effect on muscle AMX concentration profiles. In conclusion, this study has clearly demonstrated that AMX distributes rapidly and extensively within muscle interstitial fluid, consistent with theory, and that altered muscle blood flow seems unlikely to have a major effect on these distribution characteristics.

A previously developed pharmacokinetic-pharmacodynamic (PK-PD) modeling approach was used to investigate the epileptogenic activity of norfloxacin as a representative antibiotic with concentration-dependent antimicrobial activity. Rats received an intravenous infusion of norfloxacin at a rate of 5 mg kg of body weight −1 min −1 over 30 min. Blood samples were collected for drug assay, and an electroencephalogram (EEG) was recorded during infusion and postinfusion. An important delay was observed between concentrations of norfloxacin in plasma and the EEG effect. Indirect effect models failed to describe these data, which were successfully fitted by using an effect compartment model with a spline function to describe the relationship between effect and concentration at the effect site, as previously observed with imipenem. The robustness of the PK-PD model was then assessed by keeping the dose constant but increasing the duration of infusion to 120 and 240 min. Although this was accompanied by PK modifications, PD parameters did not vary significantly, and the PK-PD model still applied. In conclusion, the successful PK-PD modeling of the norfloxacin EEG effect in rats should be considered to predict and reduce the epileptogenic risk associated with this antibiotic as a representative fluoroquinolone (E. Fuseau and L. B. Sheiner, Clin. Pharmacol. Ther. 35:733-741, 1984).

International audience; Nonlinear mixed effects models (NLMEM) are used in pharmacokinetics to analyse concentrations of patients during drug development, particularly for pediatric studies. Approaches based on the Fisher information matrix can be used to optimize their design. Local design needs some a priori parameter values which might be difficult to guess. Therefore two-stage adaptive designs are useful to provide some flexibility. We implemented in the R function PFIM the Fisher matrix for two-stage designs in NLMEM. We evaluated, with simulations, the impact of one-stage and two-stage designs on the precision of parameter estimation when the true and a priori parameters are different.

The effect of probenecid (PRO) on norfloxacin (NOR) blood-brain barrier transport was investigated with rats by microdialysis. Maximum brain drug concentrations were rapidly attained, and the brain penetration factor was close to 5% in the absence and presence of PRO. In conclusion, PRO has no effect on NOR blood-brain barrier transport.

International audience; Methods and software tools for optimal design in nonlinear mixed effect models, based on the Fisher information matrix, have been developed for a decade. Academic groups regularly proposed new versions. Present tools do not incorporate adaptive designs for these models although prior information is needed and adaptive designs are increasingly used in drug development. We conducted a study among drug companies of the Drug and Disease Model Resources consortium to identify current practices and expectations.

International audience; BACKGROUND: Since 2007, it is mandatory for the pharmaceutical companies to submit a Paediatric Investigation Plan to the Paediatric Committee at the European Medicines Agency for any drug in development in adults, and it often leads to the need to conduct a pharmacokinetic study in children. Pharmacokinetic studies in children raise ethical and methodological issues. Because of limitation of sampling times, appropriate methods, such as the population approach, are necessary for analysis of the pharmacokinetic data. The choice of the pharmacokinetic sampling design has an important impact on the precision of population parameter estimates. Approaches for design evaluation and optimization based on the evaluation of the Fisher information matrix (M(F)) have been proposed and are now implemented in several software packages, such as PFIM in R. OBJECTIVES: The objectives of this work were to (1) develop a joint population pharmacokinetic model to describe the pharmacokinetic characteristics of a drug S and its active metabolite in children after intravenous drug administration from simulated plasma concentration-time data produced using physiologically based pharmacokinetic (PBPK) predictions; (2) optimize the pharmacokinetic sampling times for an upcoming clinical study using a multi-response design approach, considering clinical constraints; and (3) evaluate the resulting design taking data below the lower limit of quantification (BLQ) into account. METHODS: Plasma concentration-time profiles were simulated in children using a PBPK model previously developed with the software SIMCYP(®) for the parent drug and its active metabolite. Data were analysed using non-linear mixed-effect models with the software NONMEM(®), using a joint model for the parent drug and its metabolite. The population pharmacokinetic design, for the future study in 82 children from 2 to 18 years old, each receiving a single dose of the drug, was then optimized using PFIM, assuming identical times for parent and metabolite concentration measurements and considering clinical constraints. Design evaluation was based on the relative standard errors (RSEs) of the parameters of interest. In the final evaluation of the proposed design, an approach was used to assess the possible effect of BLQ concentrations on the design efficiency. This approach consists of rescaling the M(F), using, at each sampling time, the probability of observing a concentration BLQ computed from Monte-Carlo simulations. RESULTS: A joint pharmacokinetic model with three compartments for the parent drug and one for its active metabolite, with random effects on four parameters, was used to fit the simulated PBPK concentration-time data. A combined error model best described the residual variability. Parameters and dose were expressed per kilogram of bodyweight. Reaching a compromise between PFIM results and clinical constraints, the optimal design was composed of four samples at 0.1, 1.8, 5 and 10 h after drug injection. This design predicted RSE lower than 30 % for the four parameters of interest. For this design, rescaling M(F) for BLQ data had very little influence on predicted RSE. CONCLUSION: PFIM was a useful tool to find an optimal sampling design in children, considering clinical constraints. Even if it was not forecasted initially by the investigators, this approach showed that it was really necessary to include a late sampling time for all children. Moreover, we described an approach to evaluate designs assuming expected proportions of BLQ data are omitted.

Pharmacokinetic-pharmacodynamic modeling of the electroencephalogram (EEG) effect was combined with intracerebral glutamate determinations using microdialysis for rats receiving norfloxacin intravenously (150 mg/kg of body weight). The EEG effect (accompanied by tremors and seizures) was consistently observed without glutamate level modifications. Therefore, norfloxacin-inducted seizures are not triggered by intracerebral glutamate level enhancement.

Nonlinear mixed effects models allow investigating individual differences in drug concentration profiles (pharmacokinetics) and responses. Pharmacogenetics focuses on the genetic component of this variability. Two tests often used to detect a gene effect on a pharmacokinetic parameter are (1) the Wald test, assessing whether estimates for the gene effect are significantly different from 0 and (2) the likelihood ratio test comparing models with and without the genetic effect. Because those asymptotic tests show inflated type I error on small sample size and/or with unevenly distributed genotypes, we develop two alternatives and evaluate them by means of a simulation study. First, we assess the performance of the permutation test using the Wald and the likelihood ratio statistics. Second, for the Wald test we propose the use of the F-distribution with four different values for the denominator degrees of freedom. We also explore the influence of the estimation algorithm using both the first-order conditional estimation with interaction linearization-based algorithm and the stochastic approximation expectation maximization algorithm. We apply these methods to the analysis of the pharmacogenetics of indinavir in HIV patients recruited in the COPHAR2-ANRS 111 trial. Results of the simulation study show that the permutation test seems appropriate but at the cost of an additional computational burden. One of the four F-distribution-based approaches provides a correct type I error estimate for the Wald test and should be further investigated.

International audience; This study aimed to develop a joint population pharmacokinetic model for an antipsychotic agent in development (S33138) and its active metabolite (S35424) produced by reversible metabolism. Because such a model leads to identifiability problems and numerical difficulties, the model building was performed using the FOCE-I and the Stochastic Approximation Expectation Maximization (SAEM) estimation algorithms in NONMEM and MONOLIX, respectively. Four different structural models were compared based on Bayesian information criteria. Models were first written as ordinary differential equations systems and then in closed form (CF) to facilitate further analyses. The impact of polymorphisms on genes coding for the CYP2C19 and CYP2D6 enzymes, respectively involved in the parent drug and the metabolite elimination were investigated using permutation Wald test. The parent drug and metabolite plasma concentrations of 101 patients were analyzed on two occasions after 4 and 8 weeks of treatment at 1, 3, 6, and 24 h following daily oral administration. All configurations led to a two compartment model with back-transformation of the metabolite into the parent drug and a first-pass effect. The elimination clearance of the metabolite through other processes than back-transformation was decreased by 35% [9-53%] in CYP2D6 poor metabolizer. Permutation tests were performed to ensure the robustness of the analysis, using SAEM and CF. In conclusion, we developed a complex joint pharmacokinetic model adequately predicting the impact of CYP2D6 polymorphisms on the parent drug and its metabolite concentrations through the back-transformation mechanism.

The original publication is available at www.springerlink.com; International audience; Pharmacogenetics is now widely investigated and health institutions acknowledge its place in clinical pharmacokinetics. Our objective is to assess through a simulation study, the impact of design on the statistical performances of three different tests used for analysis of pharmacogenetic information with nonlinear mixed effects models: (i) an ANOVA to test the relationship between the empirical Bayes estimates of the model parameter of interest and the genetic covariate, (ii) a global Wald test to assess whether estimates for the gene effect are significant, and (iii) a likelihood ratio test (LRT) between the model with and without the genetic covariate. We use the stochastic EM algorithm (SAEM) implemented in MONOLIX 2.1 software. The simulation setting is inspired from a real pharmacokinetic study. We investigate four designs with N the number of subjects and n the number of samples per subject: (i) N = 40/n = 4, similar to the original study, (ii) N = 80/n = 2 sorted in 4 groups, a design optimized using the PFIM software, (iii) a combined design, N = 20/n = 4 plus N = 80 with only a trough concentration and (iv) N = 200/n = 4, to approach asymptotic conditions. We find that the ANOVA has a correct type I error estimate regardless of design, however the sparser design was optimized. The type I error of the Wald test and LRT are moderatly inflated in the designs far from the asymptotic (

International audience; PURPOSE: To compare results of population PK analyses obtained with a full empirical design (FD) and an optimal sparse design (MD) in a Drug-Drug Interaction (DDI) study aiming to evaluate the potential CYP3A4 inhibitory effect of a drug in development, SX, on a reference substrate, midazolam (MDZ). Secondary aim was to evaluate the interaction of SX on MDZ in the in vivo study. Methods To compare designs, real data were analysed by population PK modelling technique using either FD or MD with NONMEM FOCEI for SX and with NONMEM FOCEI and MONOLIX SAEM for MDZ. When applicable a Wald test was performed to compare model parameter estimates, such as apparent clearance (CL/F), across designs. To conclude on the potential interaction of SX on MDZ PK, a Student paired test was applied to compare the individual PK parameters (i.e. log(AUC) and log(C(max))) obtained either by a non-compartmental approach (NCA) using FD or from empirical Bayes estimates (EBE) obtained after fitting the model separately on each treatment group using either FD or MD. RESULTS: For SX, whatever the design, CL/F was well estimated and no statistical differences were found between CL/F estimated values obtained with FD (CL/F = 8.2 l/h) and MD (CL/F = 8.2 l/h). For MDZ, only MONOLIX was able to estimate CL/F and to provide its standard error of estimation with MD. With MONOLIX, whatever the design and the administration setting, MDZ CL/F was well estimated and there were no statistical differences between CL/F estimated values obtained with FD (72 l/h and 40 l/h for MDZ alone and for MDZ with SX, respectively) and MD (77 l/h and 45 l/h for MDZ alone and for MDZ with SX, respectively). Whatever the approach, NCA or population PK modelling, and for the latter approach, whatever the design, MD or FD, comparison tests showed that there was a statistical difference (P < 0.0001) between individual MDZ log(AUC) obtained after MDZ administration alone and co-administered with SX. Regarding C(max), there was a statistical difference (P < 0.05) between individual MDZ log(C(max)) obtained under the 2 administration settings in all cases, except with the sparse design with MONOLIX. However, the effect on C(max) was small. Finally, SX was shown to be a moderate CYP3A4 inhibitor, which at therapeutic doses increased MDZ exposure by a factor of 2 in average and almost did not affect the C(max). CONCLUSION: The optimal sparse design enabled the estimation of CL/F of a CYP3A4 substrate and inhibitor when co-administered together and to show the interaction leading to the same conclusion as the full empirical design.

The objective of this study was to investigate the contribution of norfloxacin blood–brain barrier (BBB) transport to its delayed electroencephalogram (EEG) effect in rats. Norfloxacin was injected as a bolus dose of 150 mg kg−1. Blood samples were collected for total norfloxacin plasma concentration measurements. The corresponding unbound levels were determined in brain extracellular fluid (ECF) using microdialysis. Quantitative EEG recording was conducted during 9 h post-dose. Brain ECF norfloxacin concentrations were much lower than plasma levels (AUC ratio=9.7±2.8%) but peaked very early, and concentration versus time profiles were parallel in both biological fluids. The best pharmacokinetic (PK) modelling was obtained by considering that ECF concentrations were part of the central compartment, with a proportionality factor. The peak of EEG effect was delayed and the effect versus plasma concentration curves exhibited a dramatic hysteresis. A PK–pharmacodynamic (PD) effect compartment model with a spline function to describe the relationship between effect and concentration at the effect site successfully described the data. Comparisons of PK–PD parameters estimated from plasma and ECF concentrations show that most of the delayed norfloxacin EEG effect is not due to BBB transport, but also that PD parameters derived from plasma data must be carefully interpreted when drug distribution at the effect site is restricted, as may often be the case for centrally acting drugs. Keywords: Microdialysis, electroencephalogram, seizures, PK–PD modelling, fluoroquinolones, rat Introduction Pharmacokinetic–pharmacodynamic (PK–PD) modelling has recently been used to characterise the electroencephalographic (EEG) effect of imipenem (Dupuis et al., 2001) and norfloxacin (Chenel et al., 2003) in rats. A dramatic time delay was observed between the EEG effect versus time and plasma concentration versus time curves. Experimental data were fitted using PK–PD models, including an effect compartment, and the robustness of these models was assessed in various experimental settings. However, these models may lack physiological significance. Owing to their physicochemical characteristics, the distribution of these two antibiotics in the brain is presumably low, which could possibly explain or at least contribute to the delayed effect. In contrast, according to Hammarlund-Udenaes et al. (1997), when the drug influx rate to the brain is slower than the efflux, drug concentrations within the brain are expected to be lower than in the plasma, but with an early peak. Furthermore, in a recent study using intracerebral microdialysis in rats, we have observed that the ECF brain concentrations of norfloxacin are actually much lower than the corresponding plasma levels, but also that they peak early after drug administration as a bolus dose (Marchand et al., 2003). Therefore, the delayed EEG effect of norfloxacin should not be explained by its central nervous system (CNS) distribution characteristics. The above studies (Chenel et al., 2003; Marchand et al., 2003) were conducted in different rats with different experimental settings (e.g. different doses). Thus, the objective of the current study was to appropriately characterise the contribution of the CNS distribution process to the overall delayed effect. This study was based on simultaneous CNS distribution using microdialysis and PK–PD modelling of the EEG effect in the same rat treated with a convulsant dose of norfloxacin.

S78454 (also called PCI-24781) is a novel, orally bioavailable, hydroxamate-based pan-HDACi currently being tested in clinical trials in the US and EU. Two single agent Phase I studies have been completed in patients with advanced solid tumors. The first study (PCYC-0402) tested 4 different administration schedules with 4-week (wk) cycles: 1) TID and 2) BID 5 days/wk 3 wks/4; 3) BID on days 1–7 & 15–21; 4) BID on days 1–5 & 15–19. The planned dose levels were 30, 45, 60, 75 and 90 mg/m2. The second study (CL1-002) tested an additional schedule: 5) BID 14 days on/7 days off in 3-week cycles. The maximum tolerated dose (MTD), defined by dose-limiting toxicity (DLT) occurring in at least 2/3 or 2/6 patients so that further dose-escalation is not undertaken, was established at 75 mg/m2 BID and the recommended dose at 60 mg/m2 BID. The DLT, consistently observed across all these schedules, was reversible thrombocytopenia, common to all HDACi and directly related to their mechanism of action. To address this issue, a pharmacokinetic/pharmacodynamic (PK/PD) model was used to predict the optimal administration schedule that could allow higher doses with minimal thrombocytopenia. Platelet-time profiles from the first study (PCYC-0402) were analyzed using a mixed effect model (Sheiner et al, 1980). This approach allowed the estimation of population mean parameters and their variability by the simultaneous fit of all data records, despite their heterogeneity. This PK/PD model consisted in a previously developed PK model of S78454 used as an input into a semi-physiological platelet time-course model (Friberg et al, 2002). This model was made up of a self-renewing progenitor compartment connected to transit compartments representing maturation process, and leading to a circulating compartment. A feedback mechanism mimicked the increase in proliferation rate when platelet count falls below baseline. The model was parameterized with baseline, feedback intensity, mean maturation time of platelets, and drug concentration/effect relationship. It adequately fitted observed platelet counts after different administration schedules and showed good prediction of adverse events related to haematological toxicity. Several administration schedules were simulated using this model. The results showed that a 4 days on/ 3 days off schedule was associated with the smallest platelet decrease. Accordingly, the protocol for study CL1–002 was amended. After reaching MTD in schedule 5, subsequent cohorts received S78454 on a revised schedule of 4 days on/ 3 days off, starting at one dose level below MTD (60 mg/m2 BID). The dose escalation followed the same design as the initial schedule with the same pre-defined dose levels. As expected, the dose escalation continued for two more dose levels beyond the MTD in schedule 5. At 105 mg/m2 BID, 2/3 patients included experienced DLTs of grade 4 thrombocytopenia, and therefore the MTD in the revised schedule was established at 105 mg/m2 BID. In conclusion, early understanding of toxicities and PK determination allowed us to build a PK/PD model of thrombocytopenia, which predicted the optimal administration schedule. This optimized schedule is planned to be used in future trials with S78454/PCI-24781. This also provides a successful example of modeling early in development that can be applied to other therapeutics. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2011 Nov 12-16; San Francisco, CA. Philadelphia (PA): AACR; Mol Cancer Ther 2011;10(11 Suppl):Abstract nr A192.