Bioavailability of three novel oral, sustained-release pellets, relative to an immediate-release tablet containing 500 mg flucytosine: A randomized, open-label, crossover study in healthy volunteers

The opportunistic fungal infection cryptococcal meningoencephalitis is a major cause of death among people living with HIV in sub-Saharan Africa. We report pharmacokinetic (PK) and safety data from a randomized, four-period crossover phase I trial of three sustained-release (SR) oral pellet formulations of 5-flucytosine conducted in South Africa. These formulations were developed to require less frequent administration, to provide a convenient alternative to the current immediate release (IR) formulation, A. Formulations B, C, and D were designed to release 5-flucytosine as a percentage of the nominal dose in vitro. We assessed their safety and PK profiles in a single dose (1 × 3000 mg at 0 h), relative to commercial IR tablets (Ancotil 500 mg tablets; 3 × 500 mg at 0 h and 3 × 500 mg at 6 h) in healthy, fasted participants. Forty-two healthy participants were included. All treatments were well-tolerated. The primary PK parameters, maximum observed plasma concentration (Cmax ) and area under the concentration-time profiles, were significantly lower for the SR formulations than for the IR tablets, and the geometric mean ratios fell outside the conventional bioequivalence limits. The median maximum time to Cmax was delayed for the SR pellets. Physiologically-based PK modeling indicated a twice-daily 6400 mg dose of SR formulation D in fasted condition would be optimal for further clinical development. This regimen is predicted to result in a rapid steady-state plasma exposure with effective and safe trough plasma concentration and Cmax values, within the therapeutic boundaries relative to plasma exposure after four times per day administration of IR tablets (PACTR202201760181404).

Virtual patients, digital twins and causal disease models: paving the ground for in silico clinical trials

Computational models are being explored to simulate in silico the efficacy and safety of drug candidates and medical devices. Disease models that are based on patients' profiling data are being produced to represent interactomes of genes or proteins and to infer causality in the pathophysiology {AuQ: Edit OK?}, which makes it possible to mimic the impact of drugs on relevant targets. Virtual patients designed from medical records as well as digital twins were generated to simulate specific organs and to predict treatment efficacy at the individual patient level {AuQ: Edit OK?}. As the acceptance of digital evidence by regulators grows, predictive artificial intelligence (AI)-based models will support the design of confirmatory trials in humans and will accelerate the development of efficient drugs and medical devices.

How to Successfully Generate an Alternative Approach to a Thorough QT Study: GLPG1972 as an Example

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.

Mechanistic Modeling of the Interplay Between Host Immune System, IL-7 and UCART19 Allogeneic CAR-T Cells in Adult B-cell Acute Lymphoblastic Leukemia

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.

Modeling restoration of gefitinib efficacy by co-administration of MET inhibitors in an EGFR inhibitor-resistant NSCLC xenograft model: A tumor-in-host DEB-based approach

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.

Reduced physiologically-based pharmacokinetic model of dabigatran etexilate-dabigatran and its application for prediction of intestinal P-gp-mediated drug-drug interactions

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/AUC_control (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/AUC_control = 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.

Modeling Approach to Predict the Impact of Inflammation on the Pharmacokinetics of CYP2C19 and CYP3A4 Substrates

PURPOSE

For decades, inflammation has been considered a cause of pharmacokinetic variability, mainly in relation to the inhibitory effect of pro-inflammatory cytokines on the expression level and activity of cytochrome P450 (CYP). In vitro and clinical studies have shown that two major CYPs, CYP2C19 and CYP3A4, are both impaired. The objective of the present study was to quantify the impact of the inflammatory response on the activity of both CYPs in order to predict the pharmacokinetic profile of their substrates according to systemic C-reactive protein (CRP).

METHODS

The relationships between CRP concentration and both CYPs activities were estimated and validated using clinical data first on midazolam then on voriconazole. Finally, clinical data on omeprazole were used to validate the findings. For each substrate, a physiologically based pharmacokinetics model was built using a bottom-up approach, and the relationships between CRP level and CYP activities were estimated by a top-down approach. After incorporating the respective relationships, we compared the predictions and observed drug concentrations.

RESULTS

Changes in pharmacokinetic profiles and parameters induced by inflammation seem to be captured accurately by the models.

CONCLUSIONS

These findings suggest that the pharmacokinetics of CYP2C19 and CYP3A4 substrates can be predicted depending on the CRP concentration.

Quantitative Systems Pharmacology Approaches for Immuno-Oncology: Adding Virtual Patients to the Development Paradigm

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.

Application of Item Response Theory to Model Disease Progression and Agomelatine Effect in Patients with Major Depressive Disorder

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.

Prediction of Human Nonlinear Pharmacokinetics of a New Bcl-2 Inhibitor Using PBPK Modeling and Interspecies Extrapolation Strategy

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 a new selective Bcl-2 (B-cell lymphoma 2) inhibitor developed by Servier Laboratories and used to restore apoptosis functions in cancer patients. The aim of this work was to develop a translational approach using physiologically based (PB) pharmacokinetic (PK) modeling for interspecies extrapolation to anticipate the nonlinear PK behavior of this new compound in patients. A PBPK mouse model was first built using a hybrid approach, defining scaling factors (determined from in vitro data) to correct in vitro clearance parameters and predicted Kp (partition coefficient) values. The qualification of the hybrid model using these empirically determined scaling factors was satisfactorily completed with rat and dog data, allowing extrapolation of the PBPK model to humans. Human PBPK simulations were then compared with clinical trial data from a phase 1 trial in which the drug was given orally and daily to cancer patients. Human PBPK predictions were within the 95% prediction interval for the eight dose levels, taking into account both the nonlinear dose and time dependencies occurring in S 55746 kinetics. Thus, the proposed PK interspecies extrapolation strategy, based on preclinical and in vitro information and physiologic assumptions, could be a useful tool for predicting human plasma concentrations at the early stage of drug development.

Correction to: Model-Based Adaptive Optimal Design (MBAOD) Improves Combination Dose Finding Designs: an Example in Oncology

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.

Integrated Pharmacokinetic/Pharmacodynamic Model of a Bispecific CD3xCD123 DART Molecule in Nonhuman Primates: Evaluation of Activity and Impact of Immunogenicity

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.

Getting Innovative Therapies Faster to Patients at the Right Dose: Impact of Quantitative Pharmacology Towards First Registration and Expanding Therapeutic Use

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.

Prediction of renal transporter-mediated drug-drug interactions for a drug which is an OAT substrate and inhibitor using PBPK modelling

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 C_max and IC₅₀, therefore this can be considered a potentially valuable tool within drug development.

Combined Analysis of Phase I and Phase II Data to Enhance the Power of Pharmacogenetic Tests

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.

Good Practices in Model-Informed Drug Discovery and Development: Practice, Application, and Documentation

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 R&D efficiency; 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.

Improvement of Parameter Estimations in Tumor Growth Inhibition Models on Xenografted Animals: a Novel Method to Handle the Interval Censoring Caused by Measurement of Smaller Tumors

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%).

Pharmacometrics Markup Language (PharmML): Opening New Perspectives for Model Exchange in Drug Development

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.

Comparison of Nonlinear Mixed Effects Models and Noncompartmental Approaches in Detecting Pharmacogenetic Covariates

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 pharmacokinetics (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 (H1), six 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.

Influence of covariance between random effects in design for nonlinear mixed-effect models with an illustration in pediatric pharmacokinetics

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.

Pharmacokinetic/Pharmacodynamic modeling of abexinostat-induced thrombocytopenia across different patient populations: application for the determination of the maximum tolerated doses in both lymphoma and solid tumour patients

BACKGROUND

In the clinical development of oncology drugs, the recommended dose is usually determined using a 3 + 3 dose-escalation study design. However, this phase I design does not always adequately describe dose-toxicity relationships.

METHODS

125 patients, with either solid tumours or lymphoma, were included in the study and 1217 platelet counts were available over three treatment cycles. The data was used to build a population pharmacokinetic/pharmacodynamic (PKPD) model using a sequential modeling approach. Model-derived Recommended Doses (MDRD) of abexinostat (a Histone Deacetylase Inhibitor) were determined from simulations of different administration schedules, and the higher bound for the probability of reaching these MDRD with a 3 + 3 design were obtained.

RESULTS

The PKPD model developed adequately described platelet kinetics in both patient populations with the inclusion of two platelet baseline counts and a disease progression component for patients with lymphoma. Simulation results demonstrated that abexinostat administration during the first 4 days of each week in a 3-week cycle led to a higher MDRD compared to the other administration schedules tested, with a maximum probability of 40 % of reaching these MDRDs using a 3 + 3 design.

CONCLUSIONS

The PKPD model was able to predict thrombocytopenia following abexinostat administration in both patient populations. A model-based approach to determine the recommended dose in phase I trials is preferable due to the imprecision of the 3 + 3 design.

Development of a physiologically-based pharmacokinetic model of the rat central nervous system

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.

Influence of the duration of intravenous drug administration on tumor uptake

Enhancing tumor uptake of anticancer drugs is an important therapeutic goal, because insufficient drug accumulation is now considered to be an important reason for unresponsiveness or resistance to antitumor therapy. Based on a mechanistic tumor uptake model describing tumor exposure to molecules of different molecular size after bolus administration, we have investigated the influence of the duration of intravenous administration on tumor uptake. The model integrates empirical relationships between molecular size and drug disposition (capillary permeability, interstitial diffusivity, available volume fraction, and plasma clearance), together with a compartmental pharmacokinetics model and a drug/target binding model. Numerical simulations were performed using this model for protracted intravenous drug infusion, a common mode of administration of anticancer drugs. The impact of mode of administration on tumor uptake is described for a large range of molecules of different molecular size. Evaluation was performed not only for the maximal drug concentration achieved in the tumor, but also for the dynamic profile of drug concentration. It is shown that despite a lower maximal uptake for a given dose, infusion allows for a prolonged exposure of tumor tissues to both small- and large-sized molecules. Moreover, infusion may allow higher doses to be administered by reducing Cmax-linked toxicity, thereby achieving a similar maximal uptake compared to bolus administration.

Optimal sampling times for a drug and its metabolite using SIMCYP(®) simulations as prior information

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.

Application of hematological toxicity modeling in clinical development of abexinostat (S-78454, PCI-24781), a new histone deacetylase inhibitor

PURPOSE

A population pharmacokinetic/pharmacodynamic (PK/PD) model was developed to describe the thrombocytopenia (dose-limiting toxicity) of abexinostat, a new histone deacetylase inhibitor. An optimal administration schedule of the drug was determined using a simulation-based approach.

METHODS

Early PK and PK/PD data were analysed using a sequential population modeling approach (NONMEM 7), allowing for the description of a PK profile and platelet-count decrease after abexinostat administration with various administration schedules. Simulations of platelet count with several administration schedules over 3-week treatment cycles (ASC) and over a day (ASD) were computed to define the optimal schedule that limits the depth of thrombocytopenia.

RESULTS

An intermediate PK/PD model accurately described the data. The administration of abexinostat during the first 4 days of each week in a 3-week cycle resulted in fewer adverse events (with no influence of ASD on platelet count profiles), and corresponded to the optimal treatment schedule. This administration schedule was clinically evaluated in a phase I clinical trial and allowed for the definition of a new maximum tolerated dose (MTD), leading to a nearly 30% higher dose-intensity than that of another previously tested schedule. Lastly, a final model was built using all of the available data.

CONCLUSIONS

The final model, characterizing the dose-effect and the dose-toxicity relationships, provides a useful modeling tool for clinical drug development.

Abstract A192: PK/PD modeling-based optimization of administration schedule for the histone deacetylase inhibitor (HDACi) S78454/PCI-24781 in phase I

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.

Development of a complex parent-metabolite joint population pharmacokinetic model

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.

Predictions of metabolic drug-drug interactions using physiologically based modelling: Two cytochrome P450 3A4 substrates coadministered with ketoconazole or verapamil

Nowadays, evaluation of potential risk of metabolic drug-drug interactions (mDDIs) is of high importance within the pharmaceutical industry, in order to improve safety and reduce the attrition rate of new drugs. Accurate and early prediction of mDDIs has become essential for drug research and development, and in vitro experiments designed to evaluate potential mDDIs are systematically included in the drug development plan prior to clinical assessment. The aim of this study was to illustrate the value and limitations of the classical and new approaches available to predict risks of DDIs in the research and development processes. The interaction of cytochrome P450 (CYP) 3A4 inhibitors (ketoconazole and verapamil) with midazolam was predicted using the inhibitor concentration/inhibition constant ([I]/K(i)) approach, the static approach with added variability (Simcyp(R)), and whole-body physiologically based pharmacokinetic (WB-PBPK) modelling (acslXtreme(R)). Then an in-house reference drug was used to challenge the different approaches based on the midazolam experience. Predicted values (pharmacokinetic parameters, the area under the plasma concentration-time curve [AUC] ratio and plasma concentrations) were compared with observed values obtained after intravenous and oral administration in order to assess the accuracy of the prediction methods. With the [I]/K(i) approach, the interaction risk was always overpredicted for the midazolam substrate, regardless of its route of administration and the coadministered inhibitor. However, the predictions were always satisfactory (within 2-fold) for the reference drug. For the Simcyp(R) calculations, two of the three interaction results for midazolam were overpredicted, both when midazolam was given orally, whereas the prediction obtained when midazolam was administered intravenously was satisfactory. For the reference drug, all predictions could be considered satisfactory. For the WB-PBPK approach, all predictions were satisfactory, regardless of the substrate, route of administration, dose and coadministered inhibitor. DDI risk predictions are performed throughout the research and development processes and are now fully integrated into decision-making processes. The regulatory approach is useful to provide alerts, even at a very early stage of drug development. The 'steady state' approach in Simcyp(R) improves the prediction by using physiological knowledge and mechanistic assumptions. The DDI predictions are very useful, as they provide a range of AUC ratios that include individuals at the extremes of the population, in addition to the 'average tendency'. Finally, the WB-PBPK approach improves the predictions by simulating the concentration-time profiles and calculating the related pharmacokinetic parameters, taking into account the time of administration of each drug - but it requires a good understanding of the absorption, distribution, metabolism and excretion properties of the compound.

Pharmacogenetics and population pharmacokinetics: impact of the design on three tests using the SAEM algorithm

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 (<10%). For each design, the corrected power is analogous for the three tests. Among the three designs with a total of 160 observations, the design N = 80/n = 2 optimized with PFIM provides both the lowest standard error on the effect coefficients and the best power for the Wald test and the LRT while a high shrinkage decreases the power of the ANOVA. In conclusion, a correction method should be used for model-based tests in pharmacogenetic studies with reduced sample size and/or sparse sampling and, for the same amount of samples, some designs have better power than others.

Drug-drug interaction predictions with PBPK models and optimal multiresponse sampling time designs: application to midazolam and a phase I compound. Part 2: clinical trial results

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.

Norfloxacin blood-brain barrier transport in rats is not affected by probenecid coadministration

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.

Optimal Blood Sampling Time Windows for Parameter Estimation Using a Population Approach: Design of a Phase II Clinical Trial

The objective of this paper is to determine optimal blood sampling time windows for the estimation of pharmacokinetic (PK) parameters by a population approach within the clinical constraints. A population PK model was developed to describe a reference phase II PK dataset. Using this model and the parameter estimates, D-optimal sampling times were determined by optimising the determinant of the population Fisher information matrix (PFIM) using PFIM_ _M 1.2 and the modified Fedorov exchange algorithm. Optimal sampling time windows were then determined by allowing the D-optimal windows design to result in a specified level of efficiency when compared to the fixed-times D-optimal design. The best results were obtained when K(a) and IIV on K(a) were fixed. Windows were determined using this approach assuming 90% level of efficiency and uniform sample distribution. Four optimal sampling time windows were determined as follow: at trough between 22 h and new drug administration; between 2 and 4 h after dose for all patients; and for 1/3 of the patients only 2 sampling time windows between 4 and 10 h after dose, equal to [4 h-5 h 05] and [9 h 10-10 h]. This work permitted the determination of an optimal design, with suitable sampling time windows which was then evaluated by simulations. The sampling time windows will be used to define the sampling schedule in a prospective phase II study.

Pharmacokinetic modeling of free amoxicillin concentrations in rat muscle extracellular fluids determined by microdialysis

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 (Rmodel) was estimated to 0.80 by the model, which is close to the mean value obtained by noncompartmental data analysis (Rarea= 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.

Simultaneous central nervous system distribution and pharmacokinetic-pharmacodynamic modelling of the electroencephalogram effect of norfloxacin administered at a convulsant dose in rats

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.

Dose ranging pharmacokinetics and brain distribution of norfloxacin using microdialysis in rats

The purpose of this study was to investigate the effect of dose on norfloxacin pharmacokinetics and distribution into the brain extracellular fluid (ECF), in freely moving rats. Unbound concentrations of norfloxacin in hippocampus were determined by microdialysis after an i.v. bolus dose of 12.5, 25, 50, 100, or 150 mg/kg in rats. In vivo recovery of norfloxacin was determined by retrodialysis by calibrator. Among three fluoroquinolones (enoxacin, pefloxacin, and ciprofloxacin) selected as potential calibrators, ciprofloxacin was selected as the best one. Maximum ECF brain norfloxacin concentrations are rapidly obtained but the ECFbrain/plasma areas under curves (AUC) ratios are low and independent of dose with a mean value of 8.2 +/- 5.8%. By contrast, norfloxacin systemic pharmacokinetics was nonlinear, with total plasma clearance decreasing significantly from 23.0 +/- 3.4 to 14.4 +/- 3.8 mL/min/kg when dose increased from 12.5 to 150 mg/kg.