Dose adjustment of venetoclax when co-administered with posaconazole: clinical drug-drug interaction predictions using a PBPK approach.

PURPOSE: Venetoclax, a targeted anticancer agent approved for the treatment of chronic lymphocytic leukemia and acute myeloid leukemia, is a substrate of cytochrome P450 (CYP) 3A enzyme (CYP3A4). Posaconazole, commonly used to prevent invasive fungal infections in neutropenic patients with hematological malignancies, potently inhibits CYP3A4. The purpose of this evaluation was to predict venetoclax exposures following co-administration of posaconazole at doses not previously studied clinically. METHODS: Two physiologically based pharmacokinetic (PBPK) models were developed for posaconazole based on published parameters, one for an oral suspension and another for delayed released tablets. Parameter optimization, guided by sensitivity analyses, was conducted such that the models could replicate clinical exposures of posaconazole and drug-drug interactions with sensitive CYP3A substrates including venetoclax. The clinically verified posaconazole PBPK models were then utilized to predict DDI with a previously published venetoclax PBPK model at clinically relevant dosing scenarios. RESULTS: The posaconazole PBPK models predicted posaconazole exposure and DDI related fold changes with acceptable prediction errors for both posaconazole formulations. The model predicted exposures of venetoclax, when co-administered with a 300 mg QD dose of delayed release tablets of posaconazole, were in concordance with observed data. Increasing the posaconazole dose to 500 mg QD increased venetoclax exposures by about 12% relative to 300 mg QD, which were still within the venetoclax safe exposure range. CONCLUSIONS: The posaconazole PBPK models were developed and clinically verified. Predictions using the robust PBPK model confirmed the venetoclax label recommendation of 70 mg in the presence of posaconazole at doses up to 500 mg QD.

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.

Mechanistic Modeling of Intra-Tumor Spatial Distribution of Antibody-Drug Conjugates: Insights into Dosing Strategies in Oncology.

Antibody drug conjugates (ADCs) provide targeted delivery of cytotoxic agents directly inside tumor cells. However, many ADCs targeting solid tumors have exhibited limited clinical efficacy, in part, due to insufficient penetration within tumors. To better understand the relationship between ADC tumor penetration and efficacy, previously applied Krogh cylinder models that explore tumor growth dynamics following ADC administration in preclinical species were expanded to a clinical framework by integrating clinical pharmacokinetics, tumor penetration, and tumor growth inhibition. The objective of this framework is to link ADC tumor penetration and distribution to clinical efficacy. The model was validated by comparing virtual patient population simulations to observed overall response rates from trastuzumab-DM1 treated patients with metastatic breast cancer. To capture clinical outcomes, we expanded upon previous Krogh cylinder models to include the additional mechanism of heterogeneous tumor growth inhibition spatially across the tumor. This expansion mechanistically captures clinical response rates by describing heterogeneous ADC binding and tumor cell killing; high binding and tumor cell death close to capillaries vs. low binding, and high tumor cell proliferation far from capillaries. Sensitivity analyses suggest that clinical efficacy could be optimized through dose fractionation, and that clinical efficacy is primarily dependent on the ADC-target affinity, payload potency, and tumor growth rate. This work offers a mechanistic basis to predict and optimize ADC clinical efficacy for solid tumors, allowing dosing strategy optimization to improve patient outcomes.

Informing Development of Bispecific Antibodies Using Physiologically Based Pharmacokinetic-Pharmacodynamic Models: Current Capabilities and Future Opportunities.

Antibody therapeutics continue to represent a significant portion of the biotherapeutic pipeline, with growing promise for bispecific antibodies (BsAbs). BsAbs can target 2 different antigens at the same time, such as simultaneously binding tumor-cell receptors and recruiting cytotoxic immune cells. This simultaneous engagement of 2 targets can be potentially advantageous, as it may overcome disadvantages posed by a monotherapy approach, like the development of resistance to treatment. Combination therapy approaches that modulate 2 targets simultaneously offer similar advantages, but BsAbs are more efficient to develop. Unlike combination approaches, BsAbs can facilitate spatial proximity of targets that may be necessary to induce the desired effect. Successful development of BsAbs requires understanding antibody formatting and optimizing activity for both targets prior to clinical trials. To realize maximal efficacy, special attention is required to fully define pharmacokinetic (PK)/pharmacodynamic (PD) relationships enabling selection of dose and regimen. The application of physiologically based pharmacokinetics (PBPK) has been evolving to inform the development of novel treatment modalities such as bispecifics owing to the increase in our understanding of pharmacology, utility of multiscale models, and emerging clinical data. In this review, we discuss components of PBPK models to describe the PK characteristics of BsAbs and expand the discussion to integration of PBPK and PD models to inform development of BsAbs. A framework that can be adopted to build PBPK-PD models to inform the development of BsAbs is also proposed. We conclude with examples that highlight the application of PBPK-PD and share perspectives on future opportunities for this emerging quantitative tool.

Quantitative Assessment of Elagolix Enzyme-Transporter Interplay and Drug-Drug Interactions Using Physiologically Based Pharmacokinetic Modeling.

INTRODUCTION: Elagolix is approved for the management of moderate-to-severe pain associated with endometriosis. The aim of this analysis was to develop a physiologically based pharmacokinetic (PBPK) model that describes the enzyme-transporter interplay involved in the disposition of elagolix and to predict the magnitude of drug-drug interaction (DDI) potential of elagolix as an inhibitor of P-glycoprotein (P-gp) and inducer of cytochrome P450 (CYP) 3A4. METHODS: A PBPK model (SimCYP® version 15.0.86.0) was developed using elagolix data from in vitro, clinical PK and DDI studies. Data from DDI studies were used to quantify contributions of the uptake transporter organic anion transporting polypeptide (OATP) 1B1 and CYP3A4 in the disposition of elagolix, and to quantitatively assess the perpetrator potential of elagolix as a CYP3A4 inducer and P-gp inhibitor. RESULTS: After accounting for the interplay between elagolix metabolism by CYP3A4 and uptake by OATP1B1, the model-predicted PK parameters of elagolix along with the DDI AUC∞ and Cmax ratios, were within 1.5-fold of the observed data. Based on model simulations, elagolix 200 mg administered twice daily is a moderate inducer of CYP3A4 (approximately 56% reduction in midazolam AUC∞). Simulations of elagolix 150 mg administered once daily with digoxin predicted an increase in digoxin Cmax and AUC∞ by 68% and 19%, respectively. CONCLUSIONS: A PBPK model of elagolix was developed, verified, and applied to characterize the disposition interplay between CYP3A4 and OATP1B1, and to predict the DDI potential of elagolix as a perpetrator under dosing conditions that were not tested clinically. PBPK model-based predictions were used to support labeling language for DDI recommendations of elagolix.

Accelerating Drug Development in Pediatric Oncology With the Clinical Pharmacology Storehouse.

Pediatric drug development is a challenging process due to the rarity of the population, the need to meet regulatory requirements across the globe, the associated uncertainty in extrapolating data from adults, the paucity of validated biomarkers, and the lack of systematic testing of drugs in pediatric patients. In oncology, pediatric drug development has additional challenges that have historically delayed availability of safe and effective medicines for children. In particular, the traditional approach to pediatric oncology drug development involves conducting phase 1 studies in children once the drug has been characterized and in some cases approved for use in adults. The objective of this article is to describe clinical pharmacology factors that influence pediatric oncology trial design and execution and to highlight efficient approaches for designing and expediting oncology drug development in children. The topics highlighted in this article include (1) study design considerations, (2) updated dosing approaches, (3) ways to overcome the significant biopharmaceutical challenges unique to the oncology pediatric population, and (4) use of data analysis strategies for extrapolating data from adults, with case studies. Finally, suggestions for ways to use clinical pharmacology approaches to accelerate pediatric oncology drug development are provided.

Population Pharmacokinetic and Pharmacodynamic Modeling of Etelcalcetide in Patients with Chronic Kidney Disease and Secondary Hyperparathyroidism Receiving Hemodialysis.

INTRODUCTION: Etelcalcetide is a novel calcimimetic that binds and activates calcium-sensing receptors (CaSRs) for the treatment of secondary hyperparathyroidism (SHPT). METHODS: To assess titrated dosing regimens, population pharmacokinetic (PK) and PK/pharmacodynamic (PKPD) modeling of etelcalcetide was performed using NONMEM 7.2. In this analysis, plasma etelcalcetide, serum parathyroid hormone (PTH) and calcium (Ca) concentration-time data were collected from five phase I, II, and III clinical trials following single or multiple intravenous doses of etelcalcetide ranging from 2.5 to 60 mg. A semi-mechanistic model was used to describe the relationship between etelcalcetide, PTH, and Ca. This model included the role of PTH in Ca regulation, the feedback of Ca onto PTH production via the CaSR, and the activity of etelcalcetide plasma levels in increasing the sensitivity of the CaSR to Ca via the cooperative binding model. The impact of relevant covariates was evaluated by stepwise forward/backward selection. Model evaluation was based on standard goodness-of-fit plots and prediction-corrected visual predictive checks (pcVPCs). Simulation was conducted to evaluate titrated dosing regimens. RESULTS AND DISCUSSION: The time courses of etelcalcetide, PTH, and Ca were well-described by the model. The clearance and central volume of distribution (Vc) of etelcalcetide were 0.472 L/h and 49.9 L, respectively, while estimates of the turnover half-lives of PTH and Ca were 0.36 and 23 h, respectively. The extent of interindividual variability in model parameters was low to moderate (6-67%), and no covariates were identified as significant predictors of PK and PD variability. pcVPCs confirmed the predictive ability of the model. CONCLUSIONS: The current analysis confirms the putative mechanism of action of etelcalcetide as an allosteric activator of CaSR. Simulations showed that dose titration of etelcalcetide, rather than fixed dose, is needed to effectively decrease the PTH level in patient populations.

SRM-based measurements of proprotein convertase subtilisin/kexin type 9 and lipoprotein(a) kinetics in nonhuman primate serum.

AIM: PCSK9 and Lp(a) have been identified as potential biomarkers for cardiovascular disease. The ability to measure protein turnover rates will provide insights into the dynamic properties of these proteins and lead to better understanding of their biological roles. We aimed to implement the stable isotope-labeled tracers ([2H3]-leucine) and develop a novel LC-selected reaction monitoring (SRM) mass spectrometry (MS) method to study the kinetics of PCSK9 and Lp(a). RESULTS: A sensitive method using immunoaffinity enrichment coupled with LC-SRM MS was developed to measure the production and degradation rates of PCSK9 and Lp(a) in naive nonhuman primate serum. Comparable results were obtained from two different routes of tracer administration. CONCLUSION: Immunoaffinity enrichment coupled with LC-SRM MS demonstrated success in in vivo kinetic measurements of proteins with relatively slow turnover rate (Lp[a]) or low abundance (PCSK9) in serum.

Population pharmacokinetic/pharmacodynamic model of subcutaneous adipose 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD1) activity after oral administration of AMG 221, a selective 11ß-HSD1 inhibitor.

Inhibition of 11ß-HSD1 is hypothesized to improve measures of insulin sensitivity and hepatic glucose output in patients with type II diabetes. AMG 221 is a potent, small molecule inhibitor of 11ß-HSD1. The objective of this analysis is to describe the pharmacokinetic/pharmacodynamic (PK/PD) relationship between AMG 221 and 11ß-HSD1 inhibition in ex vivo adipose tissue samples. Healthy, obese subjects were administered a single dose of 3, 30, or 100 mg of oral AMG 221 (n = 44) or placebo (n = 11). Serial blood samples were collected over 24 hours. Subcutaneous adipose tissue samples were collected by open biopsy. Population PK/PD analysis was conducted using NONMEM. The inhibitory effects (mean ± standard error of the estimate) of AMG 221 on 11ß-HSD1 activity were directly related to adipose concentrations with I(max) (the maximal inhibition of 11ß-HSD1 activity) and IC50 (the plasma AMG 221 concentration associated with 50% inhibition of enzyme activity) of 0.975 ± 0.003 and 1.19 ± 0.12 ng/mL, respectively. The estimated baseline 11ß-HSD1 enzyme activity was 755 ± 61 pmol/mg. An equilibration rate constant (k(eo)) of 0.220 ± 0.021 h⁻¹ described the delay between plasma and adipose tissue AMG 221 concentrations. AMG 221 potently blocked 11ß-HSD1 activity, producing sustained inhibition for the 24-hour study duration as measured in ex vivo adipose samples. Early characterization of concentration-response relationships can support rational selection of dose and regimen for future studies.

Prediction of exposure-response relationships to support first-in-human study design.

In drug development, phase 1 first-in-human studies represent a major milestone as the drug moves from preclinical discovery to clinical development activities. The safety of human subjects is paramount to the conduct of these studies and regulatory considerations guide activities. Forces of evolution on the pharmaceutical industry are re-shaping the first-in-human dose selection strategy. Namely, high attrition rates in part due to lack of efficacy have led to the re-organization of research and development organizations around the umbrella of translational research. Translational research strives to bring basic research advances into the clinic and support the reverse transfer of information to enhance compound selection strategies. Pharmacokinetic/pharmacodynamic (PK/PD) modeling holds a unique position in translational research by attempting to integrate diverse sets of information. PK/PD modeling has demonstrated utility in dose selection and trial design for later stages of drug development and is now being employed with greater prevalence in the translational research setting to manage risk (i.e., oncology and inflammation/immunology). Moving from empirical E (max) models to more mechanistic representations of the biological system, a higher fidelity of human predictions is expected. Strategies that have proven useful for PK predictions are being applied to PK/PD predictions. This review article examines examples of the application of PK/PD modeling in establishing target concentrations for supporting first-in-human study design.

Application of in vivo animal models to characterize the pharmacokinetic and pharmacodynamic properties of drug candidates in discovery settings.

A goal of preclinical discovery is the identification of drug candidates suitable for clinical testing. Successful integration of in vitro and in vivo experimental data sets can afford projections of human dose regimens anticipated to be safe and therapeutically beneficial. While in vitro experiments guide new chemical syntheses and are essential to understanding drug action and disposition, in vivo characterizations provide unique insight into complex biological systems that control concentrations at the site of action and pharmacologic response. Pharmacokinetic and pharmacodynamic (PK/PD) concepts underlying drug disposition and response provide a quantitative framework with which to identify potential clinical candidates. To improve throughput in earlier stages of drug discovery, in vivo pharmacokinetic study designs such as cassette dosing and sparse sampling schemes have been utilized. In later stages of discovery, pharmacokinetic studies using chemical inhibitors or surgical and genetic animal models are used to characterize the underlying determinants of drug disposition. In a complimentary fashion, modeling of in vivo pharmacodynamic effects may quantitatively link biomarkers to pharmacological response, validate in vitro to in vivo correlations and underwrite predictions of efficacious exposure targets. When applied to in vivo discovery data, PK/PD models have aided in understanding mechanisms of pharmacological response such as receptor theory in the central nervous system and cell turnover concepts in infectious disease and oncology. This review considers the role of in vivo testing toward understanding the pharmacokinetic and pharmacodynamic attributes of lead candidates in drug discovery.

Valproic acid uptake by bovine brain microvessel endothelial cells: role of active efflux transport.

The basis for low brain permeability of valproic acid (VPA) appears to be the result of efflux transport at the blood-brain barrier (BBB); however, the identity of the putative efflux transporter has not been investigated. The objective of our studies was to determine whether the multidrug resistance-associated protein (MRP) might be involved in efflux transport of VPA. Brain microvessel endothelial cells (BMEC) were isolated from cow brains and grown to confluence. MRP messenger RNA (mRNA) in BMEC were verified by reverse transcriptase-polymerase chain reaction (RT-PCR). Functional activity was demonstrated using the steady-state retention of calcein and MRP inhibitors, indomethacin (IND) and probenecid (PRB). Probenecid (0.50 mM) and indomethacin (10 microM) produced a 26 and 13% ( P<0.05 ) elevation in steady-state cellular VPA uptake following a 30-min-incubation with tracer 3H-VPA and 30 microM cold VPA. In contrast, at higher concentrations of probenecid (2 mM) and indomethacin (500 microM), an 11 and 31% reduction in VPA uptake was observed. The biphasic pattern of VPA uptake suggested concurrent inhibition of uptake and efflux transporters by the inhibitor with differing sensitivities, i.e. the efflux transporter being more susceptible to inhibition than the influx transporter. Similar results were obtained in the MRP overexpressing cell line A549. Overall, the results suggest that MRP(s) is(are) involved in the efflux transport of VPA, but do not preclude the possible contribution(s) of other organic anion transporters. The findings also adds to the growing evidence that up-regulation of active drug efflux transporters at the BBB may contribute to the development of drug resistance to antiepileptic drug therapy.

Association of high-dose cyclophosphamide, cisplatin, and carmustine pharmacokinetics with survival, toxicity, and dosing weight in patients with primary breast cancer.

This report investigates relationships between the pharmacokinetics and pharmacodynamics of high-dose alkylators used for the treatment of primary breast cancer. Eighty-five women with primary breast cancer involving >or=10 lymph nodes received four cycles of standard-dose chemotherapy followed by a high-dose regimen consisting of: cyclophosphamide (1875 mg/m(2) once daily x 3), cisplatin (165 mg/m(2) given over 72 h), carmustine (600 mg/m(2)), and stem cell transplantation. Dosages were attenuated in patients whose body weight exceeded their calculated ideal weight by >20%. Pharmacokinetics of the high-dose chemotherapeutic agents were evaluated in each patient by collection and analysis of serial blood samples. Area under the concentration time curve (AUC) for cyclophosphamide and carmustine was highly variable (>10-fold inter-patient range) with coefficients of variation > 50%, in contrast to cisplatin exposures (2-fold range; coefficient of variation 12%). The dosing method for overweight patients resulted in significantly lower systemic exposure to cisplatin (P = 0.035). The parent cyclophosphamide clearance on the 1st day of administration was significantly higher in patients who experienced acute cardiac toxicity (n = 5; P = 0.011), whereas carmustine disposition was not found to be different in those developing pulmonary toxicity (n = 25; P = 0.96). Kaplan-Meier analysis (median follow-up of 5.9 years) demonstrated that patients with lower cyclophosphamide AUC (faster parent drug clearance to potentially cytotoxic compounds) survived longer (P = 0.031). Inter-individual differences in the pharmacokinetic disposition of high-dose chemotherapy may explain variability in both response and toxicity. Prospective strategies, which attempt to individualize AUC, should be evaluated in this setting.