Search

JNJ-73763989 is an N-acetylgalactosamine conjugated short interfering RNA combination product consisting of two triggers in clinical development for chronic hepatitis B virus (HBV) infection treatment that induces a selective degradation of all HBV mRNA transcripts. Our aim is to characterize the plasma and liver pharmacokinetics (PK) of JNJ-73763989 after intravenous and subcutaneous administration in recombinant adeno-associated (rAAV) HBV infected mice. Forty-two male rAAV-HBV infected C57Bl/6 mice received JNJ-73763989 doses of 10 mg/kg i.v. or 1, 3 and 10 mg/kg s.c. Plasma and liver concentrations were analyzed simultaneously using nonlinear mixed-effects modeling with the NONMEM 7.4. A population PK model consisting of a two-compartment disposition model with transporter-mediated drug disposition, including internalization to the liver compartment, linear elimination from plasma and liver, and first-order absorption following subcutaneous administration, was suitable to describe both plasma and liver PK. After subcutaneous dosing, absolute bioavailability was complete and flip-flop kinetics were observed. JNJ-73763989 distributes from plasma to liver via transporter-mediated liver internalization in less than 24 hours, with sustained (gt;42 days) liver exposure. The saturation of transporter-mediated liver internalization was hypothesized to be due to asialoglycoprotein receptor saturation. Increasing the dose decreased the relative liver uptake efficiency in mice for intravenously and, to a lesser extent, subcutaneously administered JNJ-73763989. Lower dose levels administered subcutaneously in mice can maximize the proportion of the dose reaching the liver. SIGNIFICANCE STATEMENT: Pharmacokinetic modeling of JNJ-73763989 liver and plasma concentration-time data in mice indicated that the proportion of JNJ-73763989 reaching the liver may be increased by administering lower subcutaneous doses compared to higher intravenous doses. Model-based simulations can be applied to optimize the dose and regimen combination.

Purpose: Rilotumumab is an investigational, fully human monoclonal antibody to hepatocyte growth factor. In a randomized phase II study, trends toward improved survival were observed with rilotumumab (7.5 or 15 mg/kg) plus epirubicin, cisplatin, and capecitabine (ECX) versus placebo plus ECX in gastric/gastroesophageal junction (GEJ) cancer patients, especially in MET-positive patients. Here, we quantitatively characterized the longitudinal exposure–response [tumor growth (TG) and overall survival (OS)] relationship for rilotumumab.Experimental Design: Rilotumumab concentrations, tumor sizes, and survival time from the phase II study were pooled to develop a longitudinal exposure versus TG model and parametric OS model that explored predictive/prognostic/treatment effects (MET expression, rilotumumab exposure, relative tumor size). Model evaluation included visual predictive checks, nonparametric bootstrap, and normalized prediction distribution errors. Simulations were undertaken to predict the relationship between rilotumumab dose and OS.Results: Rilotumumab exhibited linear time-independent pharmacokinetics not affected by MET expression. The TG model adequately described tumor size across arms. A Weibull distribution best described OS. Rilotumumab exposure and change in tumor size from baseline at week 24 were predictive of OS. MET-positive patients showed shorter survival and responded better to rilotumumab than MET-negative patients. Simulations predicted a median (95% confidence interval) HR of 0.38 (0.18–0.60) in MET-positive patients treated with 15 mg/kg rilotumumab Q3W.Conclusions: Rilotumumab plus ECX demonstrated concentration-dependent effects on OS, influenced by MET expression, and tumor size in gastric/GEJ cancer patients. These findings support the phase II testing of rilotumumab 15 mg/kg every 3 weeks in MET-positive gastric/GEJ cancer (RILOMET-1; NCT01697072). Clin Cancer Res; 21(11); 2453–61. ©2015 AACR.

Supplementary Figure 1. Goodness-of-Fit Plots for the Pharmacokinetic Model

Supplementary Figure 2. Goodness-of-Fit Plots for the Tumor Growth Model

Background Exposure–response analyses were conducted to explore the relationship between selected efficacy and safety endpoints and serum phosphate (PO4) concentrations, a potential biomarker of efficacy and safety, in locally advanced or metastatic urothelial carcinoma patients with FGFR alterations treated with erdafitinib. Methods Data from two dosing regimens of erdafitinib in a phase 2 study (NCT02365597), 6 and 8-mg/day with provision for pharmacodynamically guided titration per serum PO4 levels, were analyzed using Cox proportional hazard or logistic regression models. Efficacy endpoints were overall survival (OS), progression-free survival (PFS), and objective response rate (ORR). Safety endpoints were adverse events typical for FGFR inhibitors. Results Exposure-efficacy analyses on 156 patients (6-mg = 68; 8-mg = 88) showed that patients with higher serum PO4 levels within the first 6 weeks showed better OS (hazard ratio 0.57 [95% CI 0.46–0.72] per mg/dL of PO4; p = 0.01), PFS (hazard ratio 0.80 [0.67–0.94] per mg/dL of PO4; p = 0.01), and ORR (odds ratio 1.38 [1.02–1.86] per mg/dL of PO4; p = 0.04). Exposure-safety analyses on 177 patients (6-mg = 78; 8-mg = 99) showed that the incidence of selected adverse events associated with on-target off-tumor effects significantly rose with higher PO4. Conclusions The exploratory relationship between serum PO4 levels and efficacy/safety outcomes supported the use of pharmacodynamically guided dose titration to optimize erdafitinib’s therapeutic benefit/risk ratio. Clinical trial registration number NCT02365597.

Erdafitinib (JNJ-42756493, BALVERSA) is a tyrosine kinase inhibitor indicated for the treatment of advanced urothelial carcinoma. In this work, a translational model-based approach to inform the choice of the doses in phase 1 trials is illustrated. A pharmacokinetic (PK) model was developed to describe the time course of erdafitinib plasma concentrations in mice and rats. Data from multiple xenograft studies in mice and rats were analyzed using the Simeoni tumor growth inhibition (TGI) model. The model parameters were used to derive a range of erdafitinib exposures that might inform the choice of the doses in oncology phase 1 trials. Conversion of exposures to doses was based on preliminary PK assessments from the first-in human (FIH) study. A one-compartment PK disposition model, with linear absorption and dose-dependent clearance, adequately described the PK data in both mice and rats via an allometric scaling approach. The TGI model was able to describe tumor growth dynamics, providing quantitative measurements of erdafitinib antitumor potency in mice and rats. Based on these estimates, ranges of efficacious unbound concentration were identified for erdafitinib in mice (0.642–5.364 μg/L) and rats (0.782–2.565 μg/L). Based on the FIH data, it was possible to transpose exposures into doses and doses of above 4 mg/day provided erdafitinib exposures associated with significant TGI in animals. The findings were in agreement with the results of the FIH trial, in which the first hints of clinical activities were observed at 6 mg. The successful modeling exercise of erdafitinib preclinical data showed how translational PK-PD modeling might be a tool to help to inform the choice of the doses in FIH studies.

Graphical abstract, Highlights • Mechanistic model characterizing acute immune response and HBV system interactions. • Key role of the cellular and regulatory response triggering hepatitis B chronicity. • Modelling framework to easily incorporate and explore additional biological mechanisms., Hepatitis B liver infection is caused by hepatitis B virus (HBV) and represents a major global disease problem when it becomes chronic, as is the case for 80–90% of vertical or early life infections. However, in the vast majority (>95%) of adult exposures, the infected individuals are capable of mounting an effective immune response leading to infection resolution. A good understanding of HBV dynamics and the interaction between the virus and immune system during acute infection represents an essential step to characterize and understand the key biological processes involved in disease resolution, which may help to identify potential interventions to prevent chronic hepatitis B. In this work, a quantitative systems pharmacology model for acute hepatitis B characterizing viral dynamics and the main components of the innate, adaptive, and tolerant immune response has been successfully developed. To do so, information from multiple sources and across different organization levels has been integrated in a common mechanistic framework. The final model adequately describes the chronology and plausibility of an HBV-triggered immune response, as well as clinical data from acute patients reported in the literature. Given the holistic nature of the framework, the model can be used to illustrate the relevance of the different immune pathways and biological processes to ultimate response, observing the negligible contribution of the innate response and the key contribution of the cellular response on viral clearance. More specifically, moderate reductions of the proliferation of activated cytotoxic CD8+ lymphocytes or increased immunoregulatory effects can drive the system towards chronicity.

This work explores the application of a physiologically structured population (PSP) framework in modeling hepatitis C virus (HCV) kinetics. To do so, a model was developed for the viral RNA load in plasma and liver as observed in 15 patients treated with a combination therapy of pegylated interferon, ribavirin, and telaprevir. By including both intracellular and extracellular processes of the HCV lifecycle, the model provided a description of the treatment effect on the intracellular HCV lifecycle. Combining PSP models with a nonlinear mixed effects approach in a single model permits a natural inclusion of the direct‐acting antiviral effect on intracellular processes, which can then be integrated with the viral kinetics within the host while accounting for the interindividual variability between patients. This should allow an exploration of the treatment effect within the entire chronic HCV‐infected population.

Ponesimod is a sphingosphine-1-phosphate receptor modulator being developed for the treatment of multiple sclerosis. The effects of disease-modifying treatments on magnetic resonance imaging (MRI) lesions in relapsing multiple sclerosis accurately predict effects on clinical relapses, therefore MRI lesion counts are generally accepted efficacy endpoints in phase II clinical studies of multiple sclerosis disease-modifying treatments. Here, we characterize the effect of ponesimod systemic exposure on the cumulative number of T1 gadolinium-enhancing (Gd+) lesions and the annualized relapse rate in a phase IIb study. This study assessed the cumulative number of new Gd+ lesions on T1-weighted MRI scans (primary endpoint) at weeks 12, 16, 20, and 24 and the annualized relapse rate (secondary endpoint). The effect of the demographic and prognostic covariates of sex, age, weight, T1 Gd+ lesions at baseline, and Expanded Disability Status Scale score at baseline were explored. Analyses were performed using NONMEM, Version 7.3.0 (ICON plc). An increase in ponesimod exposure led to a statistically significant decrease in the cumulative T1 Gd+ lesions on MRI from week 12 to 24 of treatment. Increasing the ponesimod daily dose beyond 20 mg did not provide significant additional benefits. Sex, age, T1 Gd+ lesions at baseline, and Expanded Disability Status Scale score at baseline were associated with a higher number of new cumulative T1 Gd+ from week 12 to 24 of treatment. This analysis shows a relationship between ponesimod exposure and the cumulative number of new T1 Gd+ lesions. Sex, age, T1 Gd+ lesions at baseline, and Expanded Disability Status Score at baseline were not found to be importantly associated with the magnitude of ponesimod effect, and consequently, there is no indication from these analyses that dosage adjustments based on the explored covariates are warranted. ClinicalTrials.gov Identifier: NCT01006265, registration date 1 November, 2009.

Purpose: To evaluate the relationship between exposure of apalutamide and its active metabolite, N-desmethyl-apalutamide, and selected clinical efficacy and safety parameters in men with high-risk nonmetastatic castration-resistant prostate cancer. Patients and Methods: An exploratory exposure–response analysis was undertaken using data from the 1,207 patients (806 apalutamide and 401 placebo) enrolled in the SPARTAN study, including those who had undergone dose reductions and dose interruptions. Univariate and multivariate Cox regression models evaluated the relationships between apalutamide and N-desmethyl-apalutamide exposure, expressed as area under the concentration–time curve at steady state, and metastasis-free survival (MFS). Univariate and multivariate logistic regression models assessed the relationship between apalutamide and N-desmethyl-apalutamide exposure and common treatment-emergent adverse events including fatigue, fall, skin rash, weight loss, and arthralgia. Results: A total of 21% of patients in the apalutamide arm experienced dose reductions diminishing the average daily dose to 209 mg instead of 240 mg. Within the relatively narrow exposure range, no statistically significant relationship was found between MFS and apalutamide and N-desmethyl-apalutamide exposure. Within apalutamide-treated subjects, skin rash and weight loss had a statistically significant association with higher apalutamide exposure. Conclusions: The use of apalutamide at the recommended dose of 240 mg once daily provided a similar delay in metastases across the SPARTAN patient population, regardless of exposure. The exploratory exposure–safety analysis supports dose reductions in patients experiencing adverse events.

Acute exposure to high doses of radiation leads to severe myelosuppression, but few treatments are currently available to treat hematopoietic syndrome of acute radiation syndrome. Granulocyte colony stimulating factors (e.g., filgrastim) stimulate proliferation of neutrophil precursors and enhance mature neutrophil function. Owing to ethical constraints on conducting clinical research in lethally irradiated humans, we developed a model‐based strategy to integrate preclinical experience in irradiated nonhuman primates (NHPs) and other clinical myelosuppressive conditions to inform filgrastim dosing to treat hematopoietic syndrome of acute radiation syndrome. Models predicting neutrophil counts and overall survival based on drug exposures were calibrated and scaled from NHPs to adult and pediatric human subjects. Several scenarios were examined investigating variations in filgrastim doses, dose frequency, treatment initiation, and duration, as well as the effect of age and radiation dose rate. Model‐based simulations and established safety profiles supported that a subcutaneous filgrastim dose of 10 µg/kg once daily provides a significant survival benefit (50%) over placebo in both adults and children, provided that the treatment is initiated within 1–14 days after radiation exposure and lasts 2–3 weeks. For treatment durations of longer than 3 weeks, filgrastim treatment is not expected to provide significantly greater benefit. This survival benefit is expected to hold for the wide range of radiation doses and dose rates (0.01–1,000 Gy/hours) examined.

Apalutamide plus androgen-deprivation therapy (ADT) has been approved for treatment of patients with metastatic castration-sensitive prostate cancer (mCSPC) based on data from phase 3 TITAN study. This analysis was conducted to describe pharmacokinetics of apalutamide and N-desmethyl-apalutamide and explore relationships between apalutamide exposure and selected clinical efficacy and safety observations.1052 patients were randomized to apalutamide + ADT (n = 525) or placebo + ADT (n = 527). A previously developed population pharmacokinetic model was applied. Cox regression analysis investigated the relationships between apalutamide exposure and overall survival (OS; n = 1004) and radiographic progression-free survival (rPFS; n = 1003). Logistic regression analysis assessed the relationships between apalutamide exposure and selected clinically relevant adverse events (n = 1051).Apalutamide + ADT treatment was efficacious in extending rPFS and OS versus placebo + ADT. Within a relatively narrow apalutamide exposure range (coefficient of variation: 22%), no statistical association was detected between rPFS, OS and apalutamide exposure quartiles. Incidence of skin rash and pruritus increased significantly with increasing apalutamide exposure.Differences in apalutamide exposure were not associated with clinically relevant differences in rPFS or OS in patients with mCSPC. Patients with increased apalutamide exposure are more likely to develop skin rash and pruritus. Dose reductions may improve these adverse events, based on an individual risk-benefit approach.

A population pharmacokinetic (PK) model was developed using data pooled from 6 clinical studies (3 in healthy volunteers and 3 in cancer patients) to characterize total and free plasma concentrations of erdafitinib following single- and multiple-dose administration, to understand clinically relevant covariates, and to quantify the inter- and intraindividual variability in erdafitinib PK. An open, linear, 3-compartment disposition model with first-order absorption and a lag time was used to describe the PK profile of total and free erdafitinib plasma concentrations. The PK of erdafitinib were linear and time independent. After oral administration, erdafitinib was rapidly absorbed, with a time to maximum concentration between 2 and 4 hours. In patients, erdafitinib total apparent oral clearance was 0.200 L/h (median free fraction = 0.24%), and the effective terminal half-life of total drug was 76.4 hours. Interindividual variability in PK parameters was moderate for oral clearance and central volume of distribution, and large for absorption rate and peripheral volume of distribution. Sex and renal function were significant covariates on free oral clearance, while weight, sex, and α1 -acid-glycoprotein were significant on oral central volume of distribution. Age, race, and mild hepatic impairment were not significant covariates of erdafitinib exposure. Given that the magnitude of the covariate effects were within 25% of reference values and that the recommended dosing regimen of erdafitinib comprises individual dose up-titrations and reductions based on presence or absence of toxicities, the clinical relevance of the investigated covariates is expected to be limited, and no dose adjustments are warranted.

Background Macitentan is a clinically approved endothelin receptor antagonist for the treatment of pulmonary arterial hypertension (PAH). Increasing use of combination drug therapy in PAH means that it is important to recognize potential drug–drug interactions (DDIs) that could affect the efficacy and safety of macitentan in patients with PAH. Objective Two Phase 1 studies were conducted to investigate the effect of macitentan at steady-state on the pharmacokinetics of the breast cancer resistance protein (BCRP) substrates, rosuvastatin and riociguat in healthy male subjects. Another objective was to determine the safety and tolerability of concomitant administration of rosuvastatin or riociguat with macitentan. Methods Healthy male subjects received a single oral dose of rosuvastatin 10 mg (n = 20) or riociguat 1 mg (n = 20) on Day 1 (reference treatment). A loading oral dose of macitentan 30 mg was administered on Day 5 followed by macitentan 10 mg once-daily from Day 6 to Day 15 (riociguat study) or Day 6 to Day 16 (rosuvastatin study). A concomitant oral dose of rosuvastatin 10 mg or riociguat 1 mg was administered on Day 10 (test treatment). Pharmacokinetics were evaluated for 96 h after treatment on Day 1 and for 144 h (riociguat study) or 168 h (rosuvastatin study) after treatment on Day 10. To compare the reference and test treatments, the geometric mean ratio was calculated for the maximum plasma concentration (Cmax), the area under the plasma concentration-time curve (AUC) from zero (pre-dose) to time of the last measured concentration above the limit of quantification (AUC0–t), the AUC from zero to infinity (AUC0–∞) and the terminal elimination half-life (t½) of rosuvastatin, riociguat and riociguat’s metabolite, M1. The difference in the time to reach maximum plasma concentration (tmax) was determined by the Wilcoxon test. Trough levels of macitentan and its metabolite, ACT-132577, were measured and safety was monitored throughout. Results Ninety percent confidence intervals of the geometric mean ratios were within the bioequivalence criteria of 0.80–1.25. There was no significant difference between test and reference tmax. Rosuvastatin or riociguat did not affect the steady-state concentrations of macitentan and ACT-132577. The adverse event profile was consistent with the known safety profiles of the drugs. Conclusions Macitentan 10 mg did not affect the pharmacokinetics of BCRP substrates, rosuvastatin or riociguat in healthy male subjects. EudraCT numbers: 2017–003095–31 and 2017–003502–41.

Tobacco use is a major health concern. To assist smoking cessation, nicotine replacement therapy (NRT) is used to reduce nicotine craving. We quantitatively described the relationship between nicotine pharmacokinetics (PKs) from NRTs and momentary craving, linking two different pharmacodynamic (PD) scales for measuring craving. The dataset comprised retrospective data from 17 clinical studies and included 1,077 adult smokers with 39,802 craving observations from four formulations: lozenge, gum, mouth spray, and patch. A PK/PD model was developed that linked individual predicted nicotine concentrations with the categorical and visual analogue PD scales through a joint bounded integer model. A maximum effect model, accounting for acute tolerance development, successfully related nicotine concentrations to momentary craving. Results showed that all formulations were similarly effective in reducing craving, albeit with a fourfold lower potency for the patch. Women were found to have a higher maximal effect of nicotine to reduce craving, compared with men.

Apalutamide is a next-generation androgen receptor inhibitor approved for treatment of subjects with high-risk, non-metastatic, castration-resistant prostate cancer (NM-CRPC). The objective of this study was to characterize the population pharmacokinetics of apalutamide and its metabolite N-desmethyl-apalutamide in healthy male and castration-resistant prostate cancer subjects. Plasma concentration data for apalutamide and N-desmethyl-apalutamide from 1092 subjects (seven clinical studies) receiving oral apalutamide (30–480 mg) once daily were pooled for a population pharmacokinetic analysis using a non-linear mixed-effect modelling approach. The impact of clinically relevant covariates was also assessed. Apalutamide absorption was rapid, and the apparent steady-state volume of distribution was large (276 L), reflecting a wide body distribution. Apalutamide was eliminated slowly, with its apparent clearance increasing from 1.31 L/h after the first dose to 2.04 L/h at steady state. No evidence of time-dependent disposition was observed for N-desmethyl-apalutamide, which was also widely distributed and slowly cleared (1.5 L/h). After 4 weeks of treatment, more than 95% of steady-state exposure of apalutamide and N-desmethyl-apalutamide was reached. At a dose of apalutamide 240 mg/day, apalutamide and N-desmethyl-apalutamide exposure exhibited 5.3- and 85.2-fold accumulation in plasma, respectively. Inter-individual variability in apalutamide apparent clearance is low (

Macitentan and its active metabolite, aprocitentan, are non-peptide, potent, dual endothelin receptor antagonists. Macitentan is approved for the treatment of pulmonary arterial hypertension in adults, at a dose of 10 mg/day. The objective of this study was to develop a comprehensive population model to describe the pharmacokinetics of macitentan and aprocitentan in healthy adults and adult subjects with pulmonary arterial hypertension. Pharmacokinetic data of 452 subjects in nine studies, after single and repeated doses (dose range 0.2–600 mg), were pooled for a non-linear mixed-effects analysis and the assessment of covariates, i.e., body weight, age, sex, race, renal and hepatic impairment, health status (healthy volunteers vs patients with pulmonary arterial hypertension), and formulation (capsules vs tablets) on pharmacokinetic parameters. The final model was an open one-compartment disposition model, with linear elimination for macitentan and linear formation and elimination for aprocitentan. A semi-mechanistic absorption model described the dose dependency and multiple peaks observed for macitentan. For a female patient with pulmonary arterial hypertension after oral administration at 10 mg, macitentan reached a maximum concentration after 9 h and, following daily dosing, reached steady state after 3 days with a twofold accumulation factor. The apparent volume of distribution was 34 L and clearance was 1.39 L/h. Aprocitentan reached maximum concentration after 51 h and steady state after 9 days, with a 12.5-fold accumulation factor. Body weight, sex, race, renal impairment, health status, and formulation were statistically significant covariates on pharmacokinetic parameters. The comprehensive population pharmacokinetic model adequately described the pharmacokinetics of macitentan and aprocitentan across different dose concentrations, regimens, and formulations. Several covariates significantly influenced the pharmacokinetics of macitentan and aprocitentan, but none was considered clinically relevant.

The single-tablet regimen darunavir/cobicistat/emtricitabine/tenofovir alafenamide (D/C/F/TAF) 800/150/200/10 mg has undergone phase III studies AMBER (NCT02431247) and EMERALD (NCT02269917) in HIV-infected patients. An existing population pharmacokinetic (PopPK) model for cobicistat-boosted darunavir (DRV) was updated to describe DRV PK in AMBER and EMERALD. For TAF, a PopPK model was developed using richly sampled phase I/II data and updated with sparsely sampled AMBER data. Individual exposure metrics for DRV and TAF in patients receiving D/C/F/TAF were derived (AMBER, n=356; EMERALD, n=750). The DRV PopPK model is a two-compartment model with sequential zero-order, first-order input. TAF PK is described by a one-compartment model with dual parallel input for absorption (slow and fast pathway). DRV covariates were α1-acid-glycoprotein and body weight. TAF covariates were lean body weight and α1-acid-glycoprotein. DRV and TAF PK were unaffected by age, race, or gender. Estimated DRV mean (SD) C0h and AUC24h, respectively, were 1899 (759) ng/mL and 87,909 (20,232) ng*h/mL in AMBER; 1813 (859) ng/mL and 85,972 (22,413) ng*h/mL in EMERALD. Estimated TAF mean (SD) AUC24h was 132 (41) ng*h/mL. These PK parameters were in line with historical data. No apparent relationships of DRV or TAF exposure with efficacy (virologic response) or safety (metabolic, cardiac, liver, gastrointestinal, skin, bone, renal, pancreas, lipid events) parameters were seen. Additionally, our findings demonstrate that in patients with low plasma concentrations, there is no risk of decreased virologic response or virologic rebound. This supports the use of a once-daily, single-tablet regimen of D/C/F/TAF 800/150/200/10 mg for the treatment of HIV-1-infected subjects.

The COVID‐19 pandemic has forced clinical studies to accommodate imposed limitations. In this study, the bioequivalence part could not be conducted as planned. Thus, the aim was to demonstrate bioequivalence, using an adaptive study design, of tadalafil in fixed‐dose combination (FDC) tablets of macitentan/tadalafil with single macitentan and tadalafil (Canadian‐sourced) tablets and assess the effect of food on FDC tablets in healthy subjects. This Phase 1, single‐center, open‐label, single‐dose, two‐part, two‐period, randomized, crossover study enrolled 62 subjects. Tadalafil bioequivalence as part of FDC of macitentan/tadalafil (10/40 mg) with single‐component tablets of macitentan (10 mg) and tadalafil (40 mg) was determined by pharmacokinetic (PK) assessment under fasted conditions. The effect of food on FDC was evaluated under fed and fasted conditions. Fasted 90% confidence intervals (CIs) for geometric mean ratios (GMRs) were within bioequivalence limits for tadalafil and macitentan. Fed and fasted 90% CIs for area under the curve (AUC) GMR were within bioequivalence limits. However, 90% CIs for maximum plasma concentration (C max) GMR for macitentan and tadalafil were outside bioequivalence limits. One FDC‐treated subject experienced a serious adverse event of transient ischemic attack (bioequivalence part). To address pandemic‐imposed limitations, an adaptive study design was implemented to demonstrate that the FDC tablet was bioequivalent to the free combination of macitentan and tadalafil (Canadian‐sourced). No clinically significant differences in PK were determined between fed and fasted conditions; the FDC formulation could be taken irrespective of meals. The FDC formulation under fasted and fed conditions was well tolerated with no clinically relevant differences in safety profiles between the treatment groups. NCT Number: NCT04235270., The macitentan 10 mg/tadalafil 40 mg Fixed Dose Combination is bioequivalent to the free combination. No clinically significant differences in PK were determined between fed and fasted conditions of the FDC formulation.

The aim of this study was to characterize the relationship between ponesimod plasma concentrations and the temporal evolution of lymphocyte counts in multiple sclerosis (MS) patients. Population pharmacokinetic (PK) and PK/pharmacodynamic (PD) models were developed using data from phase I, II, and III trials, and the impact of clinically relevant covariates on PK and PD parameters was assessed. Simulations were conducted to evaluate the maximal lymphocyte count reduction after ponesimod treatment, and the time required for total lymphocyte counts to return to normal values after treatment interruption. In MS patients, ponesimod PK were characterized by a low mean apparent plasma clearance (5.52 L/h) and a moderate mean apparent volume of distribution at steady state (239 L). The model developed indicated that none of the evaluated covariates (age, sex, formulation, food, body weight, clinical condition, and renal impairment) had a clinically relevant impact on the PK/PD parameters. In MS patients, total lymphocyte counts were characterized by a maximum reduction of 88.0% and a half maximal inhibitory concentration (IC50) of 54.9 ng/mL. Simulations indicated that in patients with normal hepatic function treated with ponesimod 20 mg daily, total lymphocyte counts were reduced to 41% of baseline at trough. After stopping treatment, lymphocyte counts were restored to normal levels within one week. The population PK/PD model well-characterized the PK of ponesimod and the time course of total lymphocyte counts in MS patients. Additionally, none of the evaluated covariates had a clinically relevant impact. This should be taken into consideration when assessing the risk of infection, administration of live-attenuated vaccines, and concomitant use of immunosuppressants.

Objectives To characterize the population pharmacokinetics of the rilpivirine long-acting (LA) formulation after intramuscular administration. Methods Rich and sparse rilpivirine plasma concentration data were obtained from seven clinical studies. In total, 18 261 rilpivirine samples were collected from 986 subjects (131 healthy subjects from Phase I studies and 855 people living with HIV from Phase IIb/III studies). Doses ranged from 300 to 1200 mg, as single-dose or multiple-dose regimens (every 4 or 8 weeks). In Phase III studies, an initiation injection of 900 mg followed by continuation injections of 600 mg every 4 weeks was used. Non-linear mixed-effects modelling was performed using NONMEM® software. Results A one-compartment model with linear elimination and two parallel absorption pathways (fast and slow) with sequential zero-first-order processes adequately captured rilpivirine flip-flop pharmacokinetics after intramuscular administration of the LA formulation. The estimated apparent elimination half-life of rilpivirine LA was 200 days. None of the evaluated covariates (age, body weight, BMI, sex, race, health status and needle length) had a clinically relevant impact on rilpivirine pharmacokinetics. Conclusions The population pharmacokinetic model suitably describes the time course and associated variability of rilpivirine plasma concentrations after rilpivirine LA intramuscular administration. The monthly regimen consists of an oral lead-in period (rilpivirine 25 mg tablets once daily for 4 weeks), followed by an initiation injection of 900 mg rilpivirine LA, then 600 mg rilpivirine LA continuation injections monthly. The absence of a clinically relevant effect of covariates on rilpivirine pharmacokinetics suggests that rilpivirine LA dose adjustments for specific subgroups are not warranted.

Delay differential equations (DDEs) are commonly used in pharmacometric models to describe delays present in pharmacokinetic and pharmacodynamic data analysis. Several DDE solvers have been implemented in NONMEM 7.5 for the first time. Two of them are based on algorithms already applied elsewhere, while others are extensions of existing ordinary differential equations (ODEs) solvers. The purpose of this tutorial is to introduce basic concepts underlying DDE based models and to show how they can be developed using NONMEM. The examples include previously published DDE models such as logistic growth, tumor growth inhibition, indirect response with precursor pool, rheumatoid arthritis, and erythropoiesis-stimulating agents. We evaluated the accuracy of NONMEM DDE solvers, their ability to handle stiff problems, and their performance in parameter estimation using both first-order conditional estimation (FOCE) and the expectation–maximization (EM) method. NONMEM control streams and excerpts from datasets are provided for all discussed examples. All DDE solvers provide accurate and precise solutions with the number of significant digits controlled by the error tolerance parameters. For estimation of population parameters, the EM method is more stable than FOCE regardless of the DDE solver. published

Esketamine nasal spray is approved for treatment-resistant depression. The objective of this study was to characterize the pharmacokinetics of esketamine and noresketamine in healthy subjects and patients with treatment-resistant depression. Esketamine and noresketamine were measured in > 9000 plasma samples collected from 820 individuals who received esketamine by the intranasal, intravenous, and oral routes. An open linear model for esketamine (three compartments) and noresketamine (two compartments) that included a hepato-portal compartment was developed using NONMEM® VII. The effects of covariates on esketamine pharmacokinetics and a model evaluation were performed using conventional methods. The fraction of a 28-mg intranasal dose absorbed through the nasal cavity (FRn) is 54% (100% of this fraction is completely absorbed); the remaining 46% is swallowed and undergoes intestinal and first-pass metabolism and 18.6% of the swallowed dose reaches the systemic circulation. The absolute bioavailability of 56 and 84 mg of intranasal esketamine is 54 and 51%, respectively. Esketamine volume at steady state and clearance were 752 L and 114 L/h, respectively. Noresketamine volume at steady state and apparent clearance were 185 L and 38 L/h, respectively. Relative to non-Asian subjects, Asian subjects showed a 64.0 and 19.4% decrease in the esketamine elimination rate constant and noresketamine apparent clearance, respectively. Japanese subjects exhibited a 34% increase in FRn vs other races. Hepatic blood flow decreased by 21.9 L/h for each decade in age in subjects aged > 60 years. These changes resulted in esketamine and noresketamine maximum concentration and area under the concentration–time curve after 24 h post-dose values that were up to 36% higher than those observed in other races or in younger adult subjects. Esketamine and noresketamine pharmacokinetics was successfully characterized in healthy subjects and patients with treatment-resistant depression. The model quantified esketamine absolute nasal and oral bioavailability, its hepatic flow-limited clearance and biotransformation to the major metabolite noresketamine, and the influence of intrinsic and extrinsic factors on esketamine pharmacokinetics. Clinical trials registration numbers of the studies included in the analysis: ESKETINTRD1001 (NCT01780259), ESKETINTRD1002 (NCT01980303), ESKETINTRD1003 (NCT02129088), ESKETINTRD1008 (NCT02846519), ESKETINTRD1009 (NCT02343289), ESKETINTRD1010 (NCT02568176), ESKETINTRD1012 (NCT02345148), 54135419TRD1015 (NCT02682225), ESKETINTRD2003 (NCT01998958), ESKETINSUI2001 (NCT02133001), ESKETINTRD3001 (NCT02417064), ESKETINTRD3002 (NCT02418585), and ESKETINTRD3005 (NCT02422186).

Purpose To model absolute neutrophil count (ANC) suppression in response to acute radiation (AR) exposure and evaluate ANC time course as a predictor of overall survival (OS) in response to AR exposure with or without treatment with granulocyte colony-stimulating factor in nonhuman primates. Methods Source data were obtained from two pivotal studies conducted in rhesus macaques exposed to 750 cGy of whole body irradiation on day 0 that received either placebo, daily filgrastim, or pegfilgrastim (days 1 and 8 after irradiation). Animals were observed for 60 days with ANC measured every 1 to 2 days. The population model of ANC response to AR and the link between observed ANC time course and OS consisted of three submodels characterizing injury due to radiation, granulopoiesis, and a time-to-event model of OS. Results The ANC response model accurately described the effects of AR exposure on the duration of neutropenia. ANC was a valid surrogate for survival because it explained 76% (95% CI, 41%–97%) and 73.2% (95% CI, 38.7%–99.9%) of the treatment effect for filgrastim and pegfilgrastim, respectively. Conclusion The current model linking radiation injury to neutropenia and ANC time course to OS can be used as a basis for translating these effects to humans.

Aim The objective of the present study was to use pharmacokinetic-pharmacodynamic modelling to characterize the effects of chemotherapy on the granulopoietic system and to predict the absolute neutrophil counts (ANCs) for patients with chemotherapy-induced neutropenia treated with filgrastim and pegfilgrastim. Methods Data were extracted from 10 phase I-III studies conducted in 110 healthy adults, and 618 adult and 52 paediatric patients on chemotherapy following administration of filgrastim or pegfilgrastim. The structural model accounted for ANC dynamics and the effects of filgrastim and pegfilgrastim, chemotherapy and corticosteroids. The impact of neutrophils on drug disposition was based on a drug receptor-binding model that assumed quasi-equilibrium and stimulation of the production and maturation of neutrophils upon treatment. The chemotherapy and corticosteroid effects were represented by kinetic-pharmacodynamic-type models, where chemotherapy stimulated elimination of neutrophil precursors at the mitotic stage, and corticosteroids stimulated neutrophil production. Results The systemic half-lives of filgrastim (2.6 h) and pegfilgrastim (10.1 h) were as expected. The effective half-life of chemotherapy was 9.6 h, with a 2-day killing effect. The rate of receptor elimination from mitotic compartments exhibited extreme interindividual variability (% coefficient of variation >200), suggesting marked differences in sensitivity to chemotherapy effects on ANCs. The stimulatory effects of pegfilgrastim were significantly greater than those of filgrastim. Model qualification confirmed the predictive capability of this model. Conclusion This qualified model simulates the time course of ANC in the absence or presence of chemotherapy and predicts nadir, time to nadir and time of recovery from different grades of neutropenia upon treatment with filgrastim and pegfilgrastim.

This study aimed to develop a method for implementing dose correction in a Michaelis-Menten (M-M) approximation of a target-mediated drug disposition (TMDD) model with multiple intravenous (IV) bolus administrations. We derived the formula of a correction factor (Fcorr) for each dose in a multiple IV bolus dosing regimens for M-M model. Fcorr depends on the residual free drug amount prior IV bolus dosing event and dose amount. We conducted a stochastic simulation and estimation (SSE) exercise based on therapeutic antibody PK parameters to evaluate the effect of Fcorr on parameter estimation. Previously published clinical PK data of recombinant human erythropoietin (rHuEPO) from four clinical trials in healthy subjects receiving multiple IV bolus doses were analyzed by both M-M model with and without dose correction (MMC and MMNC) as well as the rapid-binding/quasi-steady-state (RB/QSS) TMDD models. Our results showed that MMNC introduced bias to fixed-effect parameter estimates and overestimated random-effect variables. Compared with MMC, MMNC was not able to adequately characterize the nonlinearity in the PK data of antibody and rHuEPO. The MMC-based simulation demonstrated that thricely weekly 10 IU/kg rHuEPO dosing regimen yielded Fcorr = 0.5. This result suggested that the lower-than-expected exposure for rHuEPO at low dose is due to target binding. An M-M approximation of the TMDD model should include a dose correction to avoid model misfitting and potential bias in the estimated PK parameters.

AIMS Selexipag is a prostacyclin receptor agonist approved for the treatment of pulmonary arterial hypertension. Cytochrome P450 (CYP) 2C8 is involved in the metabolism of selexipag and its active metabolite, ACT-333679. This study evaluated the interaction of selexipag and clopidogrel, a CYP2C8 inhibitor. METHODS The study had a 2-treatment, 1-sequence, crossover design. Pharmacokinetics (PK) and CYP2C8 genotype were assessed in healthy male subjects administered selexipag (200 μg twice daily [b.i.d.]) alone or with clopidogrel (300 mg single dose or 75 mg once daily [o.d.]). PK modelling and simulation were conducted to support dosing recommendations. RESULTS Clopidogrel had a comparatively small effect on selexipag (

To characterize the effect of erdafitinib on electrocardiogram (ECG) parameters and the relationship between erdafitinib plasma concentrations and QTc interval changes in patients with advanced or refractory solid tumors. Triplicate ECGs and continuous 12-lead Holter data were collected in the dose escalation part (Part 1) of the first-in-human study, with doses ranging from 0.5 to 12 mg. Triplicate ECG monitoring continued in Parts 2–4 where 2 dose regimens selected from Part 1 were expanded in prespecified tumor types. Analyses of ECG data included central tendency analyses, identification of categorical outliers and morphological assessment. A concentration–QTc analysis was conducted using a linear mixed-effect model based on extracted time matching Holter data. Central tendency, categorical outlier, and ECG morphologic analyses from 187 patients revealed no clinically significant effect of erdafitinib on heart rate, atrioventricular conduction or cardiac depolarization (PR and QRS), and no effect on cardiac repolarization (QTc). Concentration–QTc analysis from 62 patients indicated that the slopes of relationship between total and free erdafitinib plasma concentrations and QTcI (mean exponent of 0.395) were estimated as − 0.00269 ms/(ng/mL) and − 1.138 ms/(ng/mL), respectively. The predicted change in QTcI at the observed geometric mean of total and free concentration at the highest therapeutic erdafitinib dose (9 mg daily) was

The liver is a well-known immunotolerogenic environment, which provides the adequate setting for liver infectious pathogens persistence such as the hepatitis B virus (HBV). Consequently, HBV infection can derive in the development of chronic disease in a proportion of the patients. If this situation persists in time, chronic hepatitis B (CHB) would end in cirrhosis, hepatocellular carcinoma and eventually, the death of the patient. It is thought that this immunotolerogenic environment is the result of complex interactions between different elements of the immune system and the viral biology. Therefore, the purpose of this work is to unravel the mechanisms implied in the development of CHB and to design a tool able to help in the study of adequate therapies. Firstly, a conceptual framework with the main components of the immune system and viral dynamics was constructed providing an overall insight on the pathways and interactions implied in this disease. Secondly, a review of the literature was performed in a modular fashion: (i) viral dynamics, (ii) innate immune response, (iii) humoral and (iv) cellular adaptive immune responses and (v) tolerogenic aspects. Finally, the information collected was integrated into a single topological representation that could serve as the plan for the systems pharmacology model architecture. This representation can be considered as the previous unavoidable step to the construction of a quantitative model that could assist in biomarker and target identification, drug design and development, dosing optimization and disease progression analysis.

Non-adherence to antipsychotic medication is a primary factor in disease relapse in schizophrenic patients. We sought to evaluate if plasma concentrations of the antipsychotic risperidone can be used as a predictor of treatment adherence and to identify the optimal plasma concentration threshold to reliably distinguish between adherent and non-adherent patients. A population pharmacokinetic model was used to simulate plasma risperidone steady-state trough concentrations in 1000 virtual patients, where 60% of the patients were 100% adherent to their medication, while 40% of the patients were non-adherent to their medication. The probability of adherence was assessed by receiver operating characteristic (ROC) analysis on Ctrough. The area under the ROC curve (AUCROC) was used to identify the optimal Ctrough threshold. Single vs multiple Ctrough at steady state was also evaluated. After a single risperidone Ctrough measurement, the AUCROC (95% CI) was estimated to be 0.71 (0.69–0.72) and the optimal Ctrough threshold accounting for the lowest number of adherent and non-adherent misclassifications was estimated to be 11.9 ng/mL. After multiple Ctrough measurements, the AUCROC (95% CI) increased up to 0.85 (0.84–0.87) for three Ctrough measurements. The optimal probability threshold to reliably discriminate between adherent and non-adherent patients was estimated to be 0.51. Using this model which is reflective of typical adherence to antipsychotic medication, we found that three consecutive steady-state Ctrough measurements are needed for an accurate and precise diagnostic test to discriminate between patients who are adherent or non-adherent to treatment.

Etelcalcetide is a novel calcimimetic in development for the treatment of secondary hyperparathyroidism (SHPT). A population pharmacokinetic/pharmacodynamic (PK/PD) model was developed relating etelcalcetide exposures to markers of efficacy (parathyroid hormone [PTH]) and safety (calcium) using data from three clinical studies. The semimechanistic model was developed that included allosteric activation pharmacology and understanding of calcium homeostasis. The temporal profiles for all biomarkers were well described by the model. The cooperativity constant was 4.94, confirming allosteric activation mechanism. Subjects with more severe disease (higher PTH baseline) were predicted to experience less pronounced reduction in PTH (percentage change from baseline), but more reduction in calcium (Ca; percentage change from baseline). There was no evidence that dose adjustment by any covariate was needed. Model-based simulations provided quantitative support to several elements of dosing, such as starting dose, monitoring, and titration timing for registration trials.

Hepcidin (H25) is a hormone peptide synthesized by the liver that binds to ferroportin and blocks iron export. In this study, H25 was inhibited by administration of single and multiple doses of an anti-H25 monoclonal antibody Ab 12B9m in cynomolgus monkeys. The objective of this analysis was to develop a pharmacodynamic model describing the role of H25 in regulating iron homeostasis and the impact of hepcidin inhibition by Ab 12B9m. Total serum H25 and Ab 12B9m were determined in each animal. Corresponding measurements of serum iron and hemoglobin (Hb) were obtained. The PD model consisted of iron pools in serum (FeS), reticuloendothelial macrophages (FeM), hemoglobin (FeHb), and liver (FeL). The iron was assumed to be transported between the FeS, FeHb, and FeM unidirectionally at rates k S, k Hb, and k M. H25 serum concentrations were described by the previously developed PK model with the parameters fixed at their estimates. The serum iron and Hb data were fitted simultaneously. The corresponding estimates of the rate constants were k S/Fe0 = 0.113 h−1, k M = 0.00191 h−1, and k Hb = 0.00817 h−1. The model-based IC50 value for the H25 inhibitory effect on ferroportin activity was 0.398 nM. The PD model predicted a negligible effect of Ab 12B9m on Hb levels for the tested doses. The presented PD model adequately described the serum iron time courses following single and multiple doses of Ab 12B9m. Ab 12B9m-induced inhibition of H25 resulted in a temporal increase in serum and liver iron and a decrease in the iron stored in reticuloendothelial macrophages.

To characterize the trabectedin population pharmacokinetics in children and adolescent patients with cancer and compare it with the trabectedin pharmacokinetics in adults. Plasma concentrations from ten adolescent and three children with cancer (age range 4.0–17.0 years) treated with trabectedin at doses ranging from 1.1 to 1.7 mg/m2, administered as a 24-h continuous intravenous infusion every 3 weeks, were available for the analysis. An external model evaluation was performed to verify whether a previously developed adult population pharmacokinetic model was predictive of the pediatric plasma concentrations of trabectedin. The maximum a posteriori estimation of the individual pharmacokinetic parameters for pediatric patients was conducted, after successful completion of the external evaluation step. The relationships between pharmacokinetic parameters and body size were evaluated. External evaluation methods showed no major differences between the adult population and children and adolescent patients of this study. The mean ± standard deviation (SD) of the individual estimated clearance and central volume of distribution in these children/adolescent patients was 36.4 ± 16.1 L/h and 13.2 ± 6.54 L, respectively. These values were similar to the typical values reported for adult patients—37.6 L/h and 13.9 L (for females) and 16.1 L (for males). The median area under the plasma concentration versus time curve (AUC) in children/adolescent patients was 55.1 µg h/L, while in the adult population the median AUC was 61.3 µg h/L, both administered a 1.5 mg/m2 dose regimen with mean (range) BSA for adults = 1.86 (0.90–2.80) vs children/adolescent patients = 1.49 (0.66–2.54). The adult population pharmacokinetic model adequately described the trabectedin plasma concentrations and its variability in the pediatric population of patients involved in this assessment that mostly comprised adolescents. The trabectedin systemic exposure achieved in this population was comparable (within 12%) to the exposure obtained in adult population when the same dose, expressed in mg/m2, was administered.

The aim of this study was to characterize the pharmacokinetic drug–drug interaction (DDI) between simeprevir (NS3/4A protease inhibitor) and odalasvir (NS5A inhibitor) after oral administration to support the design and dose selection of clinical studies with this combination for the treatment of chronic hepatitis C infection (HCV). Simeprevir and odalasvir plasma concentrations from 30 healthy subjects receiving these drugs in monotherapy as well as in combination were pooled and analyzed using a population pharmacokinetic modeling approach. Previous pharmacokinetic models developed to characterize the pharmacokinetics for each drug were used as starting point. The dual effect of simeprevir and odalasvir on their pharmacokinetic parameters was explored. Simulations were performed to assess the impact of the DDI on exposure parameters. In presence of odalasvir, the relative bioavailability of simeprevir increased by 26% and the apparent clearance was reduced following competitive inhibition depending on odalasvir plasma concentrations, with an inhibitory constant (Ki) estimated to be 1610 ng/mL. The apparent odalasvir clearance was reduced by simeprevir plasma concentrations following an Imax model, characterized by a maximum inhibitory effect of 46.7% and an IC50 of 257 ng/mL. Model-based simulations indicated that both Cmax and AUC24h increased for both drugs, when co-administered. The pharmacokinetic model adequately describes the time course of plasma concentrations and their variability when simeprevir and/or odalasvir were orally administered. This model can be used as a first step to predict the exposures of concomitant administration of simeprevir and odalasvir in HCV-infected subjects. Data from study AL355-602 (NCT02512562) were used for this analysis.

The aim of the current study was to characterize the time course of plasma concentrations of AL-335 and its main metabolites (ALS-022399 and ALS-022227) after oral administration in healthy and hepatitis C virus (HCV)-infected subjects, in monotherapy as well as in combination with simeprevir and/or odalasvir. AL-335, ALS-022399, and ALS-022227 plasma concentrations from subjects receiving 800 mg of AL-335 orally once daily (qd) as monotherapy or in combination were pooled and analyzed using a nonlinear mixed effect modeling approach. The typical values (between subject variability) of AL-335 and ALS-022399 apparent linear clearances were 3300 L/h (33.9%) and 1910 L/h (30.0%), respectively. ALS-022227 elimination was characterized as a nonlinear process, with typical values of Vmax,ALS-022227 and Km,ALS-022227 estimated to be 84,799 ng/h (14.9%) and 450.2 ng/mL, respectively. AL-335 and ALS-022399 plasma concentrations were increased more than 2-fold in presence of simeprevir and/or odalasvir, while the effect on ALS-022227 plasma concentrations was limited. The effect of simeprevir and/or odalasvir might be explained by their capacity to inhibit P-glycoprotein. Internal evaluation confirmed that the population pharmacokinetic model developed was deemed appropriate to describe the time course of AL-335, ALS-022399, and ALS-022227 plasma concentrations and their associated variability in both healthy and HCV-infected subjects, as well as the interaction effect of simeprevir and/or odalasvir over AL-335 and its metabolites in healthy subjects. This model can be used as a starting point to evaluate drug-drug interaction processes in HCV-infected patients and support the development of a direct-acting antiviral (DAA) combination.

Optimal dose selection in clinical trials is problematic when efficacious and toxic concentrations are close. A novel quantitative approach follows for optimizing dose titration in clinical trials. A system of pharmacokinetics (PK), pharmacodynamics, efficacy, and toxicity was simulated for scenarios characterized by varying degrees of different types of variability. Receiver operating characteristic (ROC) and clinical trial simulation (CTS) were used to optimize drug titration by maximizing efficacy/safety. The scenarios included were a low-variability base scenario, and high residual (20%), interoccasion (20%), interindividual (40%), and residual plus interindividual variability scenarios, and finally a shallow toxicity slope scenario. The percentage of subjects having toxicity was reduced by 87.4% to 93.5%, and those having efficacy was increased by 52.7% to 243%. Interindividual PK variability may have less impact on optimal cutoff values than other sources of variability. ROC/CTS methods for optimizing dose titration offer an individualized approach that leverages exposure-response relationships.

A population pharmacokinetic pharmacodynamic (PK/PD) model describing the effect of epoetin alfa on hemoglobin (Hb) response in hemodialysis patients was developed. Epoetin alfa pharmacokinetics was described using a linear 2-compartment model. PK parameter estimates were similar to previously reported values. A maturation-structured cytokinetic model consisting of 5 compartments linked in a catenary fashion by first-order cell transfer rates following a zero-order input process described the Hb time course. The PD model described 2 subpopulations, one whose Hb response reflected epoetin alfa dosing and a second whose response was unrelated to epoetin alfa dosing. Parameter estimates from the PK/PD model were physiologically reasonable and consistent with published reports. Numerical and visual predictive checks using data from 2 studies were performed. The PK and PD of epoetin alfa were well described by the model.

Purpose: Rilotumumab is an investigational, fully human monoclonal antibody to hepatocyte growth factor. In a randomized phase II study, trends toward improved survival were observed with rilotumumab (7.5 or 15 mg/kg) plus epirubicin, cisplatin, and capecitabine (ECX) versus placebo plus ECX in gastric/gastroesophageal junction (GEJ) cancer patients, especially in MET-positive patients. Here, we quantitatively characterized the longitudinal exposure–response [tumor growth (TG) and overall survival (OS)] relationship for rilotumumab. Experimental Design: Rilotumumab concentrations, tumor sizes, and survival time from the phase II study were pooled to develop a longitudinal exposure versus TG model and parametric OS model that explored predictive/prognostic/treatment effects (MET expression, rilotumumab exposure, relative tumor size). Model evaluation included visual predictive checks, nonparametric bootstrap, and normalized prediction distribution errors. Simulations were undertaken to predict the relationship between rilotumumab dose and OS. Results: Rilotumumab exhibited linear time-independent pharmacokinetics not affected by MET expression. The TG model adequately described tumor size across arms. A Weibull distribution best described OS. Rilotumumab exposure and change in tumor size from baseline at week 24 were predictive of OS. MET-positive patients showed shorter survival and responded better to rilotumumab than MET-negative patients. Simulations predicted a median (95% confidence interval) HR of 0.38 (0.18–0.60) in MET-positive patients treated with 15 mg/kg rilotumumab Q3W. Conclusions: Rilotumumab plus ECX demonstrated concentration-dependent effects on OS, influenced by MET expression, and tumor size in gastric/GEJ cancer patients. These findings support the phase II testing of rilotumumab 15 mg/kg every 3 weeks in MET-positive gastric/GEJ cancer (RILOMET-1; NCT01697072). Clin Cancer Res; 21(11); 2453–61. ©2015 AACR.

Pazopanib, a new oral angiogenesis inhibitor, has demonstrated clinical activity against multiple solid tumors and was approved for the treatment of patients with advanced renal cell carcinoma. As an exposure-response relationship has been observed for pazopanib, its therapeutic drug monitoring could be a valuable tool in clinical practice. Therefore, the aim of this study was to develop and validate a selective and precise high performance liquid chromatography-ultraviolet method for the measurement of pazopanib in plasma from patients with cancer.After liquid-liquid extraction with diethyl ether, pazopanib and gefitinib (internal standard) were separated using isocratic elution on an Ultrabase C18 column using a mobile phase consisting of a mixture in vol/vol proportion of 47:53 of ammonium acetate (pH, 7; 0.02 mol/L) and acetonitrile/methanol (70:30, vol/vol) pumped at a constant flow rate of 1 mL/min. Quantification was performed at 260 nm. Method validation was undertaken as per the guidelines for Bioanalytical Method Validation published by the Food and Drug Administration and European Medicines Agency.Calibration curves were linear over the range 0.5-100 mcg/mL. Interday and intraday coefficients of variations were less than 4.5%. The limit of detection and the lower limit of quantification were 0.2 and 0.5 mcg/mL, respectively. Recovery of pazopanib from plasma was80%.This is the first high performance liquid chromatography-ultraviolet method for pazopanib quantification that has been validated within a wide range of plasma concentrations and is a suitable method for therapeutic drug monitoring of pazopanib.

Denosumab (XGEVA®) is a recombinant, fully human IgG2 monoclonal antibody directed against the receptor activator of nuclear factor kappa-B ligand (RANKL) that prevents differentiation of osteoclast precursors into mature osteoclasts and acceleration of bone resorption, resulting in the inhibition of osteoclast activation. Denosumab is indicated for the prevention of skeletal-related events (SREs) in adult patients with bone metastases from solid tumors at the dose of 120 mg administered subcutaneously (SC) every 4 weeks. This review is focused on describing its target-mediated disposition and direct inhibitory effect on bone resorption, as well as the modeling and simulation techniques used to integrate the PKPD information collected during clinical development of denosumab. In addition, this review further discusses the clinical relevance of patient covariate effects on denosumab systemic exposure, target engagement and downstream pharmacodynamics biomarkers, and the rationale for dosing regimen selection for Phase 3 studies. Phase 3 clinical studies demonstrated that denosumab was superior to zoledronic acid in inhibiting bone resorption and, consequently, delaying the time to first SRE by a median of 8.2 months in patients with bone metastases from solid tumors. Thus, denosumab may be considered a better alternative treatment than zoledronic acid for the prevention of SRE in patients with bone metastases from solid tumors.

Objective The aim of this study was to understand the influence of disease and patient characteristics on exposure to daratumumab, an immunoglobulin Gκ (IgGκ) monoclonal antibody, and clinical outcomes in relapsed or refractory multiple myeloma (MM). Patients and Methods Baseline myeloma type, albumin levels, renal/hepatic function, age, sex, race, weight, Eastern Cooperative Oncology Group (ECOG) status, refractory status, and number of prior therapies were evaluated using data from two clinical studies—GEN501 (N = 104) and SIRIUS (N = 124). Results Daratumumab clearance was approximately 110% higher in IgG myeloma patients than non-IgG myeloma patients, leading to significantly lower exposure in IgG myeloma patients based on maximum trough serum concentrations (p

Our objective was to compare the time course of bone mineral density (BMD) changes at the lumbar spine (LS) and total hip (TH) in postmenopausal women during treatment with denosumab, bisphosphonates, selective estrogen receptor modulators, PTH, or calcitonin.Data were extracted from 142 randomized controlled trials for prevention or treatment of postmenopausal osteoporosis representing over 113,000 women.The percent change from baseline in BMD was analyzed using a nonlinear least-squares random-effects meta-regression analysis. The dose-response relationship of BMD changes was well characterized by a maximal effect (EMax) model with a different EMax for LS and TH for each drug class. The ratio of LS and TH BMD changes was significantly different across the different drug classes. The time course of BMD changes was well characterized by an exponential onset with a different rate for LS and TH for each drug class. The dose-response relationship for denosumab showed that the approved dosing regimen of 60 mg every 6 months resulted in maximal BMD changes.This exploratory analysis shows that 3 years of treatment with denosumab resulted in bigger changes in LS and TH BMD compared with 3 years of treatment with 10 mg/d oral alendronate, 5 mg/y i.v. zoledronic acid, 5 mg/d oral risedronate, 150 mg/mo oral ibandronate, 3 mg i.v. ibandronate every 3 months, 60 mg/d oral raloxifene, and 200 IU/d calcitonin. Treatment with PTH resulted in larger changes in LS BMD compared with denosumab; however, denosumab treatment provided larger changes in TH BMD.

To investigate the relationships between AMG 811 exposure, concentration changes in serum IFN-γ, and IFN-γ-induced protein 10 (CXCL10), and to identify important contributions of baseline covariates to these relationships. A mechanism based pharmacokinetic (PK)-pharmacodynamic (PD) model was developed. A target mediated disposition model was used to describe AMG 811 and target IFN-γ interaction. CXCL10 was predicted to be driven by estimated free IFN-γ levels. For an average systemic lupus erythematosus (SLE) subject, the linear clearance (CL) of AMG 811 was 0.176 L/day, and the central (Vc) and peripheral (Vp) volumes of distribution were 1.48 and 2.12 L, respectively. Body weight was found to correlate with CL, Vc, Vp, and inter compartment clearance (Q); and age was found to correlate with Vc. The relationship between estimated free serum IFN-γ concentration levels and serum CXCL10 in logarithmic scales was best described by a linear model with slope and intercept estimated to be 0.197 and -0.3, respectively. The largest observed reduction of serum CXCL10 concentration was achieved at the highest AMG 811 dose tested (180 mg SC). This model enables simulations of AMG 811 PK-PD profiles under various dosing regimens to support future clinical studies.

In pharmacokinetics/pharmacodynamics (PKPD) the measured response is often delayed relative to drug administration, individuals in a population have a certain lifespan until they maturate or the change of biomarkers does not immediately affects the primary endpoint. The classical approach in PKPD is to apply transit compartment models (TCM) based on ordinary differential equations to handle such delays. However, an alternative approach to deal with delays are delay differential equations (DDE). DDEs feature additional flexibility and properties, realize more complex dynamics and can complementary be used together with TCMs. We introduce several delay based PKPD models and investigate mathematical properties of general DDE based models, which serve as subunits in order to build larger PKPD models. Finally, we review current PKPD software with respect to the implementation of DDEs for PKPD analysis.

Denosumab has received approval in many countries and indications include treating women with postmenopausal osteoporosis (PMO) at increased or high risk for fracture and men at high risk for fracture receiving androgen deprivation therapy (ADT) for non-metastatic prostate cancer. Increases in total hip bone mineral density (BMD) with denosumab explained a large percentage of new vertebral fracture risk reduction in women with PMO; however, this effect has not been studied in men with prostate cancer receiving ADT. We compared the relationship between the time course of BMD changes and new vertebral fracture risk reduction with denosumab in women with PMO and men with prostate cancer. After adjusting for different baseline hazards, a significant and similar relationship between time course of total hip and lumbar spine BMD changes and new vertebral fracture risk was observed in both patient populations. Time course of total hip BMD changes with denosumab was the best predictor for changes in fracture risk and explained 88% of the new vertebral fracture risk reduction in women with PMO and 91% in men with prostate cancer. Therefore, total hip BMD is a useful surrogate to measure the clinical impact of denosumab on fracture risk reduction in both patient populations.