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.