Model-Informed Clinical Development of Once-Every-6-Month Injection of Paliperidone Palmitate in Patients with Schizophrenia: A Pharmacometric Bridging Approach (Part I).

BACKGROUND AND OBJECTIVE: A model-informed drug development (MIDD) approach was implemented for paliperidone palmitate (PP) 6-month (PP6M) clinical development, using pharmacokinetics and pharmacokinetic/pharmacodynamic model-based simulations. METHODS: PP6M pharmacokinetics were simulated by extending the PP 3-month (PP3M) pharmacokinetic model to account for increased injection volume, and hence dose. Contribution of the MIDD approach to the design of the pivotal PP6M phase-3 study (PP6M/PP3M noninferiority study, NCT03345342) investigating schizophrenia relapse rates was twofold: (1) PP6M dose selection, and (2) hypothesis generation that lower trough concentrations (Ctrough) associated with PP6M, relative to PP3M, were not associated with lower efficacy, which was to be evaluated in the phase-3 study. Moreover, accompanied by an intense sampling scheme to adequately characterize paliperidone pharmacokinetics and to elucidate the potential relationship between concentration and safety/efficacy, the bridging strategy eliminated the need for additional phase-1/phase-2 clinical studies. RESULTS: Using a MIDD bridging strategy, PP6M doses were selected that, compared with PP3M, were expected to have a similar range of exposures and a noninferior relapse rate and safety profile. Clinical data from PP6M/PP3M noninferiority study confirmed that PP6M, compared with PP3M, had a similar range of exposures (T'jollyn et al. in Eur J Drug Metab Pharmacokinet 2024), as well as a noninferior relapse rate and safety profile (this manuscript). CONCLUSIONS: Consistency of the MIDD approach with observed clinical outcomes confirmed the hypothesis that lower Ctrough did not lead to increased relapse rates at the doses administered. Although higher paliperidone peak concentrations are achieved with corresponding doses of PP6M relative to PP3M in the phase-3 clinical study, types and incidences of treatment-related adverse events were comparable between PP6M and PP3M groups and no new safety concerns emerged for PP6M (Najarian et al. in Int J Neuropsychopharmacol 25(3):238-251, 2022).

Model-Informed Clinical Development of 6-Monthly Injection of Paliperidone Palmitate in Patients with Schizophrenia: Dosing Strategies Guided by Population Pharmacokinetic Modeling and Simulation (Part II).

BACKGROUND AND OBJECTIVE: Paliperidone palmitate 6-month (PP6M) intramuscular (IM) injection is the longest-acting treatment available for patients with schizophrenia. A population pharmacokinetic (popPK) modeling and simulation approach was deployed to inform dosing strategies for PP6M. METHODS: The extensive analysis database included 15,932 paliperidone samples from 700 patients receiving gluteal paliperidone palmitate 3-month (PP3M) or PP6M injections in the double-blind phase of a phase-3 noninferiority study (NCT03345342). Exposure parameters for paliperidone appeared to increase dose-proportionally within each dosing schedule (PP3M/PP6M). The range of paliperidone exposures after IM administration of PP6M overlaps with that of corresponding doses of oral paliperidone extended release, PP 1-month (PP1M), and PP3M. Model-based simulations were performed to investigate paliperidone exposures in different PP6M dosing scenarios and relevant subpopulations. RESULTS: A dosing window of ≤ 2 weeks earlier and ≤ 3 weeks later than the target 6-month interval for maintenance treatment with PP6M dosing maintains paliperidone exposures at levels that are not expected to substantially impact its safety and efficacy. For missed-dose scenarios, tailored re-initiation regimens are proposed that should be applied before resuming PP6M maintenance dosing. Regarding subpopulations, PP6M 700 mg eq. is the highest dose recommended in mild renal-impairment patients; the paliperidone pharmacokinetics after PP6M administration is not affected by sex, body mass index, or age in a clinically meaningful way. CONCLUSION: Paliperidone concentration-time profiles after PP6M and PP3M dosing were adequately described by the popPK model. Model-based simulation results provide guidance for clinicians on initiating PP6M therapy, transitioning between paliperidone formulations, the dosing windows to use for maintenance dosing, and managing missed PP6M doses.

Model-based meta-analysis to quantify the effects of short interfering RNA therapeutics on hepatitis B surface antigen turnover in hepatitis B-infected mice.

The objective of this study was to compare the efficacy of short interfering RNA therapeutics (siRNAs) in reducing hepatitis B surface antigen (HBsAg) levels in hepatitis B-infected (HBV) mice across multiple siRNA therapeutic classes using model-based meta-analysis (MBMA) techniques. Literature data from 10 studies in HBV-infected mice were pooled, including 13 siRNAs, formulated as liposomal nanoparticles (LNPs) or conjugated to either cholesterol (chol) or N-acetylgalactosamine (GalNAc). Time course of the baseline- and placebo-corrected mean HBsAg profiles were modeled using kinetics of drug effect (KPD) model coupled to an indirect response model (IRM) within a longitudinal non-linear mixed-effects MBMA framework. Single and multiple dose simulations were performed exploring the role of dosing regimens across evaluated siRNA classes. The HBsAg degradation rate (0.72 day-1) was consistent across siRNAs but exhibited a large between-study variability of 31.4% (CV%). The siRNA biophase half-life was dependent on the siRNA class and was highest for GalNAc-siRNAs (21.06 days) and lowest for chol-siRNAs (2.89 days). ID50 estimates were compound-specific and were lowest for chol-siRNAs and highest for GalNAc-siRNAs. Multiple dose simulations suggest GalNAc-siRNAs may require between 4 and 7 times less frequent dosing at higher absolute dose levels compared to LNP-siRNAs and chol-siRNAs, respectively, to reach equipotent HBsAg-lowering effects in HBV mice. In conclusion, non-clinical HBsAg concentration-time data after siRNA administration can be described using the presented KPD-IRM MBMA framework. This framework allows to quantitatively compare the effects of siRNAs on the HBsAg time course and inform dose and regimen selection across siRNA classes. These results may support siRNA development, optimize preclinical study designs, and inform data analysis methodology of future anti-HBV siRNAs; and ultimately, support siRNA model-informed drug development (MIDD) strategies.

Mechanistic modeling projections of antibody persistence after homologous booster regimens of COVID-19 vaccine Ad26.COV2.S in humans.

Mechanistic model-based simulations can be deployed to project the persistence of humoral immune response following vaccination. We used this approach to project the antibody persistence through 24 months from the data pooled across five clinical trials in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-seronegative participants following vaccination with Ad26.COV2.S (5 × 1010 viral particles), given either as a single-dose or a homologous booster regimen at an interval of 2, 3, or 6 months. Antibody persistence was quantified as the percentage of participants with detectable anti-spike binding and wild-type virus neutralizing antibodies. The projected overall 24-month persistence after single-dose Ad26.COV2.S was 70.5% for binding antibodies and 55.2% for neutralizing antibodies, and increased after any homologous booster regimen to greater than or equal to 89.9% for binding and greater than or equal to 80.0% for neutralizing antibodies. The estimated model parameters quantifying the rates of antibody production attributed to short-lived and long-lived plasma cells decreased with increasing age, whereas the rate of antibody production mediated by long-lived plasma cells was higher in women relative to men. Accordingly, a more pronounced waning of antibody responses was predicted in men aged greater than or equal to 60 years and was markedly attenuated following any homologous boosting regimen. The findings suggest that homologous boosting might be a viable strategy for maintaining protective effects of Ad26.COV2.S for up to 24 months following prime vaccination. The estimation of mechanistic modeling parameters identified the long-lived plasma cell pathway as a key contributor mediating antibody persistence following single-dose and homologous booster vaccination with Ad26.COV2.S in different subgroups of recipients stratified by age and sex.

Mechanistic Model Describing the Time Course of Humoral Immunity Following Ad26.COV2.S Vaccination in Non-Human Primates.

Mechanistic modeling can be used to describe the time course of vaccine-induced humoral immunity and to identify key biologic drivers in antibody production. We used a six-compartment mechanistic model to describe a 20-week time course of humoral immune responses in 56 non-human primates (NHPs) elicited by vaccination with Ad26.COV2.S according to either a single-dose regimen (1 × 1011 or 5 × 1010 viral particles [vp]) or a two-dose homologous regimen (5 × 1010 vp) given in an interval of 4 or 8 weeks. Humoral immune responses were quantified by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike-specific binding antibody concentrations as determined by spike protein-enzyme-linked immunosorbent assay. The mechanistic model adequately described the central tendency and variability of binding antibody concentrations through 20 weeks in all vaccination arms. The estimation of mechanistic modeling parameters revealed greater contribution of the antibody production mediated by short-lived cells as compared with long-lived cells in driving the peak response, especially post second dose when a more rapid peak response was observed. The antibody production mediated by long-lived cells was identified as relevant for generating the first peak and for contributing to the long-term time course of sustained antibody concentrations in all vaccination arms. The findings contribute evidence on the key biologic components responsible for the observed time course of vaccine-induced humoral immunity in NHPs and constitute a step toward defining immune biomarkers of protection against SARS-CoV-2 that might translate across species. SIGNIFICANCE STATEMENT: We demonstrate the adequacy of a mechanistic modeling approach describing the time course of binding antibody concentrations in non-human primates (NHPs) elicited by different dose levels and regimens of Ad26.COV2.S. The findings are relevant for informing the mechanism-based accounts of vaccine-induced humoral immunity in NHPs and translational research efforts aimed at identifying immune biomarkers of protection against SARS-CoV-2 infection.

Pharmacodynamic-Mediated Drug Disposition (PDMDD) Model of Daratumumab Monotherapy in Patients with Multiple Myeloma.

BACKGROUND AND OBJECTIVE: We aimed to quantify the daratumumab concentration- and CD38 dynamics-dependent pharmacokinetics using a pharmacodynamic mediated disposition model (PDMDD) in patients with multiple myeloma (MMY) following daratumumab IV or SC monotherapy. Daratumumab, a human IgG monoclonal antibody targeting CD38 with a direct on-tumor and immunomodulatory mechanism of action, has been approved to treat patients with multiple myeloma (MM). METHODS: In total, 7788 daratumumab plasma samples from 850 patients with diagnosis of MMY were used. The serum concentration-time data of daratumumab were analysed using nonlinear mixed-effects modeling with NONMEM®. The PDMDD with quasi steady-state approximation (QSS) was compared to the previously developed Michaelis-Menten (MM) approximation with respect to the parameter estimates, the goodness-of-fit plots and prediction-corrected visual predictive check, as well as model-based simulations. The effect of patients' covariates on daratumumab pharmacokinetics was also investigated. RESULTS: The QSS approximation characterized the concentration- and CD38 dynamics-dependency of daratumumab pharmacokinetics within the doses ranging from 0.1 to 24 mg/kg after IV administration and 1200 and 1800 mg after SC administration in patients with MMY, mechanistically describing the binding of daratumumab and CD38, the internalization of the daratumumab-CD38 complex and the CD38 turnover. Compared to the previously developed MM approximation, the MM approximation with the non-constant total target and dose-correction provided substantial improvement of the model fit, but it was still not as good as the QSS approximation. The effect of the previously identified covariates as well as the newly identified covariate (baseline M protein) on daratumumab pharmacokinetics was confirmed, but the magnitude of the effect was deemed not clinically relevant. CONCLUSIONS: Accounting for the CD38 turnover and its binding capacity to daratumumab, the QSS approximation provided a mechanistic interpretation of daratumumab PK parameters and consequently well described the concentration- and CD38 dynamics-dependency of daratumumab pharmacokinetics. CLINICAL STUDIES INCLUDED IN THE ANALYSIS WERE REGISTERED WITH THE NCT NUMBER BELOW AT HTTP://WWW. CLINICALTRIALS: GOV : MMY1002 (ClinicalTrials.gov: NCT02116569), MMY1003 (NCT02852837), MMY1004 (NCT02519452), MMY1008 (NCT03242889), GEN501 (NCT00574288), MMY2002 (NCT01985126), MMY3012 (NCT03277105).

Pharmacokinetic-Pharmacodynamic Modeling of the Ponesimod Effect on Heart Rate in Patients With Multiple Sclerosis.

The purpose of this study was to characterize the ponesimod effect on the heart rate (HR) in patients with multiple sclerosis (MS). A previous pharmacokinetic (PK) and pharmacodynamic model developed in healthy participants was updated using data from phase II and III trials conducted in patients with MS. Clinically relevant covariates were assessed. Simulations were conducted to evaluate the impact of the lack of adherence to ponesimod treatment and provide guidance in cases of treatment re-initiation. The maximal effect parameter of the PK/HR model was lower in patients with MS (23.5% decrease) compared with healthy volunteers (43.2%). The effect of patient covariates on PK/HR was similar to those identified in healthy participants and not clinically relevant in patients with MS. The population PK/HR model well characterized the effect of ponesimod on the time course of HR in patients with MS. After 2 weeks of treatment with 10 mg or higher doses, the model indicated full tolerance development. After repeated dosing at 20 mg, tolerance was maintained > 60% of the steady-state tolerance for up to 4 days after the last dose. Re-initiating with gradual uptitration is recommended if drug discontinuation lasts ≥ 4 days. This managed the negative chronotropic effects of ponesimod. No bradycardia events were observed within the first 2 weeks of treatment in patients with relapsing MS with a baseline HR > 55 bpm. This justifies the recommendation included in the human prescription drug labeling to monitor HR after the first ponesimod dose in these patients.

An exposure-response analysis of ponesimod clinical efficacy in a randomized phase III study in patients with relapsing multiple sclerosis.

The efficacy of ponesimod and teriflunomide for the treatment of relapsing multiple sclerosis (MS) was compared in a randomized phase III trial. This study explores the exposure-response (E-R) relationships of efficacy end points (annualized relapse rate [ARR] and combined unique active lesions [CUALs]) of ponesimod observed in this trial. The E-R relationships were described using nonlinear mixed effects models for count data. The effect of baseline covariates (demography and prognostic factors) was also explored. Ponesimod 20 mg reduced ARR (primary end point) by 30.5% (95% confidence interval [CI]: 9.8% to 46.4%) and the number of CUALs by 56% (95% CI: 46% to 64%) between baseline and week 108 compared to teriflunomide 14 mg. The E-R analyses indicated a significant relationship between ARR and CUAL. In turn, CUAL was significantly related to ponesimod systemic exposure. Based on these relationships, the predicted reduction of ARR was relatively flat in the range of ponesimod systemic exposure achieved with the 20 mg clinical dose: the expected ARR decrease ranged from 28% (95% CI: 11% to 42%) at the 5th percentile of ponesimod exposure to 34% (95% CI: 19% to 47%) at the 95th percentile. No significant baseline covariates affected the ponesimod effects and, consequently, dosage adjustments are not warranted by these analyses. Although significant relationships were found between ARR and CUAL and between ponesimod exposure and CUAL, these analyses were supportive of the use of a flat 20 mg maintenance dose for ponesimod in adult patients with MS.

Plasma and Liver Pharmacokinetics of the N-Acetylgalactosamine Short Interfering RNA JNJ-73763989 in Recombinant Adeno-Associated-Hepatitis B Virus-Infected Mice.

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 (>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.

Exposure-response analyses of erdafitinib in patients with locally advanced or metastatic urothelial carcinoma.

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's effect on serum phosphate justifies its pharmacodynamically guided dosing in patients with cancer.

A population pharmacokinetic (PK)-pharmacodynamic (PD) model was developed using data from 345 patients with cancer. The population PK-PD model evaluated the effect of erdafitinib total and free plasma concentrations on serum phosphate concentrations after once-daily oral continuous (0.5-12 mg) and intermittent (10-12 mg for 7 days on/7 days off) dosing, and investigated the potential covariates affecting erdafitinib-related changes in serum phosphate levels. Phosphate is used as a biomarker for erdafitinib's efficacy and safety: increases in serum phosphate were observed after dosing with erdafitinib, which were associated with fibroblast growth factor receptor target engagement via inhibition of renal fibroblast growth factor 23-mediated signaling. PK-PD model-based simulations were performed to assess the approved PD-guided dosing algorithm of erdafitinib (8 mg once-daily continuous dosing, with up-titration to 9 mg based on phosphate levels [<5.5 mg/dl] and tolerability at 14-21 days of treatment). The serum phosphate concentrations increased after the first dose and reached near maximal level after 14 days of continuous treatment. Serum phosphate increased with erdafitinib free drug concentrations: doubling the free concentration resulted in a 1.8-fold increase in drug-related phosphate changes. Dose adjustment after at least 14 days of dosing was supported by achievement of >95% maximal serum phosphate concentration. The peak-to-trough fluctuation within a dosing interval was limited for serum phosphate concentrations (5.68-5.65 mg/dl on Day 14), supporting phosphate monitoring at any time relative to dosing. Baseline phosphate was higher in women, otherwise, none of the investigated covariate-parameter relationships were considered clinically relevant. Simulations suggest that the starting dose of 8-mg with up-titration to 9-mg on Days 14-21 maximized the number of patients within the target serum phosphate concentrations (5.5-7 mg/dl) while limiting the number of treatment interruptions. The findings from the PK-PD model provided a detailed understanding of the erdafitinib concentration-related phosphate changes over time, which supports erdafitinib's dosing algorithm.

Bioequivalence and food effect of a fixed-dose combination of macitentan and tadalafil: Adaptive design in the COVID-19 pandemic.

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 (Cmax ) 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.

A quantitative systems pharmacology model for acute viral hepatitis B.

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.

Multiscale model of hepatitis C virus dynamics in plasma and liver following combination therapy.

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.

Delay differential equations based models in NONMEM.

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.

A Population Pharmacokinetic Model of Macitentan and Its Active Metabolite Aprocitentan in Healthy Volunteers and Patients with Pulmonary Arterial Hypertension.

BACKGROUND: 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. OBJECTIVE: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Population Pharmacokinetic Analysis of Darunavir and Tenofovir Alafenamide in HIV-1-Infected Patients on the Darunavir/Cobicistat/Emtricitabine/Tenofovir Alafenamide Single-Tablet Regimen (AMBER and EMERALD Studies).

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.

Effect of Ponesimod Exposure on Total Lymphocyte Dynamics in Patients with Multiple Sclerosis.

OBJECTIVE: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

An Exposure-Response Analysis of the Clinical Efficacy of Ponesimod in a Randomized Phase II Study in Patients with Multiple Sclerosis.

BACKGROUND AND OBJECTIVE: 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. METHODS: 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). RESULTS: 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. CONCLUSIONS: 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. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01006265, registration date 1 November, 2009.

Population Pharmacokinetics of Esketamine Nasal Spray and its Metabolite Noresketamine in Healthy Subjects and Patients with Treatment-Resistant Depression.

BACKGROUND: Esketamine nasal spray is approved for treatment-resistant depression. OBJECTIVE: The objective of this study was to characterize the pharmacokinetics of esketamine and noresketamine in healthy subjects and patients with treatment-resistant depression. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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).