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.

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.

Quantifying Antibody Persistence After a Single Dose of COVID-19 Vaccine Ad26.COV2.S in Humans Using a Mechanistic Modeling and Simulation Approach.

Understanding persistence of humoral immune responses elicited by vaccination against coronavirus disease 2019 (COVID-19) is critical for informing the duration of protection and appropriate booster timing. We developed a mechanistic model to characterize the time course of humoral immune responses in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-seronegative adults after primary vaccination with the Janssen COVID-19 vaccine, Ad26.COV2.S. The persistence of antibody responses was quantified through mechanistic modeling-based simulations. Two biomarkers of humoral immune responses were examined: SARS-CoV-2 neutralizing antibodies determined by wild-type virus neutralization assay (wtVNA) and spike protein-binding antibodies determined by indirect spike protein enzyme-linked immunosorbent assay (S-ELISA). The persistence of antibody responses was defined as the period of time during which wtVNA and S-ELISA titers remained above the lower limit of quantification. A total of 442 wtVNA and 1,185 S-ELISA titers from 82 and 220 participants, respectively, were analyzed following administration of a single dose of Ad26.COV2.S (5 × 1010 viral particles). The mechanistic model adequately described the time course of observed wtVNA and S-ELISA serum titers and its associated variability up to 8 months following vaccination. Mechanistic model-based simulations show that single-dose Ad26.COV2.S elicits durable but waning antibody responses up to 24 months following immunization. Of the estimated model parameters, the production rate of memory B cells was decreased in older adults relative to younger adults, and the antibody production rate mediated by long-lived plasma cells was increased in women relative to men. A steeper waning of antibody responses was predicted in men and in older adults.

Population Pharmacokinetics of Paliperidone Palmitate (Once-Monthly Formulation) in Japanese, Korean, and Taiwanese Patients With Schizophrenia.

The paliperidone pharmacokinetics after intramuscular administration of once-monthly paliperidone palmitate in Japanese patients were studied in 3 phase 1 studies and in 2 phase 3 studies performed in Japan, Korea, and Taiwan. These data (Japanese, n = 509; Korean, n = 31; Taiwanese, n = 47) were used to describe the paliperidone palmitate pharmacokinetics in Japanese, to compare with non-Japanese, and to validate the historical population pharmacokinetic (Pop-PK) model for paliperidone palmitate, developed using data from studies in patients with schizophrenia outside Japan. The final historical Pop-PK model, including all significant patient covariates of Japanese studies, was used to simulate paliperidone plasma concentration-time data using nonlinear mixed effects, followed by comparison with actual data. Visual predictive checks displayed considerable overlap between predicted and actual plasma concentrations; the majority of observations were within the 90% prediction interval. Japanese, Korean, and Taiwanese patients had comparable plasma concentrations. Covariate distributions demonstrated comparatively lower median body mass index in Japanese, Korean, and Taiwanese patients versus rest-of-world population. Prediction errors for the data set used for external validation were within cutoff values, confirming accuracy/precision of the model. Paliperidone pharmacokinetics were adequately predicted for Japanese studies using the historical Pop-PK model, confirming its robustness. Pharmacokinetics in Japanese, Korean, and Taiwanese patients with schizophrenia were comparable with rest-of-world population.

Systemic Exposure of Rituximab Increased by Ibrutinib: Pharmacokinetic Results and Modeling Based on the HELIOS Trial.

INTRODUCTION: In the HELIOS trial, bendamustine/rituximab (BR) plus ibrutinib (BR-I) improved disease outcomes versus BR plus placebo in previously treated chronic lymphocytic leukemia/small lymphocytic lymphoma. Here, we describe the pharmacokinetic (PK) observations, along with modeling to further explore the interaction between ibrutinib and rituximab. METHODS: 578 subjects were randomized to ibrutinib or placebo with BR (6 cycles). Ibrutinib PK samples and tumor measurements were obtained from all subjects; a subset was evaluated for bendamustine and rituximab PK. Population rituximab PK was assessed using nonlinear mixed-effects modeling. RESULTS: Dose-normalized plasma concentration-time bendamustine data were comparable between the arms. Systemic rituximab exposure was higher with BR-I versus BR; mean trough serum concentrations were 2- to 3-fold higher in the first three cycles and 1.2- to 1.7-fold higher subsequently. No relevant safety differences were observed. In the modeling, including treatment arm as a categorical covariate and tumor burden as a continuous time-varying covariate on overall rituximab clearance significantly improved fitting of the data. CONCLUSIONS: BR-I led to higher dose-normalized systemic rituximab exposure versus BR and more rapid steady-state achievement. The modeling data suggest that rituximab disposition is, at least in part, target mediated. Determining the clinical significance of these findings requires further assessments. TRIAL REGISTRATION: This study is registered at https://clinicaltrials.gov/ct2/show/NCT01611090 .

Population Pharmacokinetic Modeling of Canagliflozin in Healthy Volunteers and Patients with Type 2 Diabetes Mellitus.

BACKGROUND AND OBJECTIVES: Canagliflozin is an orally active, reversible, selective sodium-glucose co-transporter-2 inhibitor. A population pharmacokinetic (popPK) model of canagliflozin, including relevant covariates as sources of inter-individual variability, was developed to describe phase I, II, and III data in healthy volunteers and in patients with type 2 diabetes mellitus (T2DM). METHODS: The final analysis included 9061 pharmacokinetic (PK) samples from 1616 volunteers enrolled in nine phase I, two phase II, and three phase III studies and was performed using NONMEM(®) 7.1. Inter-individual variability was evaluated using an exponential model and the residual error model was additive in the log domain. The first-order conditional estimation method with interaction was applied and the model was parameterized in terms of rate constants. Covariate effects were explored graphically on empirical Bayes estimates of PK parameters, as shrinkage was low. Clinical relevance of statistically significant covariates was evaluated. The predictive properties of the model were illustrated by prediction-corrected visual predictive checks. RESULTS: A two-compartment PK model with lag-time and sequential zero- and first-order absorption and first-order elimination best described the observed data. Sex, age, and weight on apparent volume of distribution of the central compartment, body mass index on first-order absorption rate constant, and body mass index and over-encapsulation on lag-time, and estimated glomerular filtration rate (eGFR, by MDRD equation), dose, and genetic polymorphism (carriers of UGT1A9*3 allele) on elimination rate constant were identified as statistically significant covariates. The prediction-corrected visual predictive checks revealed acceptable predictive performance of the model. CONCLUSION: The popPK model adequately described canagliflozin PK in healthy volunteers and in patients with T2DM. Because of the small magnitude of statistically significant covariates, they were not considered clinically relevant. However, dosage adjustments are recommended for T2DM patients with renal impairment (eGFR ≥60 mL/min/1.73 m(2): 100 or 300 mg/day; eGFR of 45 to <60 mL/min/1.73 m(2): 100 mg/day).

The method of averaging applied to pharmacokinetic/pharmacodynamic indirect response models.

The computational effort required to fit the pharmacodynamic (PD) part of a pharmacokinetic/pharmacodynamic (PK/PD) model can be considerable if the differential equations describing the model are solved numerically. This burden can be greatly reduced by applying the method of averaging (MAv) in the appropriate circumstances. The MAv gives an approximate solution, which is expected to be a good approximation when the PK profile is periodic (i.e. repeats its values in regular intervals) and the rate of change of the PD response is such that it is approximately constant over a single period of the PK profile. This paper explains the basis of the MAv by means of a simple mathematical derivation. The NONMEM® implementation of the MAv using the abbreviated FORTRAN function FUNCA is described and explained. The application of the MAv is illustrated by means of an example involving changes in glycated hemoglobin (HbA1c%) following administration of canagliflozin, a selective sodium glucose co-transporter 2 inhibitor. The PK/PD model applied to these data is fitted with NONMEM® using both the MAv and the standard method using a numerical differential equation solver (NDES). Both methods give virtually identical results but the NDES method takes almost 8 h to run both the estimation and covariance steps, whilst the MAv produces the same results in less than 30 s. An outline of the NONMEM® control stream and the FORTRAN code for the FUNCA function is provided in the appendices.

Population pharmacokinetics of JNJ-37822681, a selective fast-dissociating dopamine D2-receptor antagonist, in healthy subjects and subjects with schizophrenia and dose selection based on simulated D2-receptor occupancy.

BACKGROUND AND OBJECTIVES: JNJ-37822681 is a selective, fast-dissociating dopamine D2-receptor antagonist currently in development as a candidate antipsychotic. The aim of the analyses was to develop a population pharmacokinetic model to describe the pharmacokinetics of JNJ-37822681 in healthy subjects and patients with schizophrenia and to identify covariates of interest. The model was then used to simulate D2-receptor occupancy in support of dose selection for subsequent studies. METHODS: Data were obtained from 378 subjects enrolled in three phase I and two phase II studies. Nonlinear mixed effects modeling of pooled data was conducted using NONMEM(®) to estimate population pharmacokinetic parameters and the effect of covariates on these parameters. The model was evaluated on a subset of data that was not used for model building and was subsequently used to simulate steady-state exposure for each subject in the phase IIb study. Striatal D2-receptor occupancy was predicted using simulated exposure combined with pharmacodynamic parameters from a sigmoid maximum effect model established from previous [(11)C]raclopride positron emission tomography studies. RESULTS: A two-compartment disposition model with zero-order input in a depot compartment followed by first-order absorption into and first-order elimination from the central compartment combined with a transit compartment provided the best fit to the data. Significant covariates were sex and bioavailability on apparent clearance and food intake on the absorption rate constant. Clearance was 11 % higher in females compared with males. The model passed external evaluation. The estimated pharmacokinetic parameters for the phase IIb study were similar to those observed in the phase IIa study. D2-receptor occupancy was predicted to be in the 65-80 % range at 10 mg twice daily and partially or fully reaching the 80 % threshold at doses of 20 and 30 mg twice daily. CONCLUSION: The population pharmacokinetic model of JNJ-37822681 successfully described the pharmacokinetics of JNJ-37822681 and allowed the reliable determination of individual exposure parameters in a phase IIb study. It was concluded that 5 or 7.5 mg twice daily would likely be minimal- or no-effect doses, whereas 10 mg twice daily was expected to provide an optimal balance of efficacy and tolerability.

D2-receptor occupancy measurement of JNJ-37822681, a novel fast off-rate D2-receptor antagonist, in healthy subjects using positron emission tomography: single dose versus steady state and dose selection.

RATIONALE: JNJ-37822681 is a highly selective, fast dissociating dopamine D2-receptor antagonist being developed for the treatment of schizophrenia. A single dose [¹¹C]raclopride positron emission tomography (PET) imaging study had yielded an estimated clinical dose range. Receptor occupancy at steady state was explored to test the validity of the single-dose estimates during chronic treatment. OBJECTIVES: The aims of this study are to characterize single and multiple dose pharmacokinetics and obtain striatal D2-receptor occupancies to predict doses for efficacy studies and assess the safety and tolerability of JNJ-37822681. METHODS: An open-label single- and multiple-dose study with 10 mg JNJ-37822681 (twice daily for 13 doses) was performed in 12 healthy men. Twenty [¹¹C]raclopride PET scans (up to 60 h after the last dose) from 11 subjects were used to estimate D2-receptor occupancy. A direct effect O (max) model was applied to explore the relationship between JNJ-37822681 plasma concentration and striatal D2-receptor occupancy. RESULTS: Steady state was reached after 4-5 days of twice daily dosing. JNJ-37822681 plasma concentrations of 3.17 to 63.0 ng/mL resulted in D2 occupancies of 0 % to 62 %. The concentration leading to 50 % occupancy was 18.5 ng/mL (coefficient of variation 3.9 %) after single dose and 26.0 ng/mL (8.2 %) at steady state. JNJ-37822681 was well tolerated. CONCLUSIONS: Receptor occupancy after single dose and at steady state differed for JNJ-37822681 and the robustness of the estimates at steady state will be tested in phase 2 studies. Dose predictions indicated that 10, 20, and 30 mg JNJ-37822681 twice daily could be suitable for these studies.