CPT: Pharmacometrics & Systems Pharmacology最新文献

Zastaprazan (JP-1366) is a novel potassium-competitive acid blocker for the treatment of acid-related disorders. We aimed to establish a population pharmacokinetic (PK) model of zastaprazan, thereby characterizing the PK of zastaprazan in patients with gastroesophageal reflux disease (GERD) as well as evaluating the impact of various covariates, including CYP2C19 phenotypes, on zastaprazan PK. This population PK analysis included zastaprazan plasma concentration–time data from 92 patients with erosive GERD and 68 healthy volunteers without any gastrointestinal disorders and was performed using nonlinear mixed-effect modeling. Simulations were conducted to predict zastaprazan PK under various dosing regimens in patients with GERD. The plasma PK of zastaprazan was adequately described by a two-compartment model with Erlang-type absorption (six sequential compartments) and first-order elimination. CYP2C19 phenotypes had no significant effect on zastaprazan PK. The disease status was identified as a significant covariate on apparent clearance of zastaprazan, showing lower values in patients with GERD compared to healthy volunteers. However, the model-based simulation indicated that the impact of disease status on zastaprazan exposure was not clinically meaningful. Overall, the current population PK model successfully characterized the observed zastaprazan PK in both patients with GERD and healthy volunteers.

Integrated modeling of the pharmacokinetic (PK) and target binding, by means of a TMDD model, can provide valuable insights into the expected pharmacodynamic (PD) effects of monoclonal antibodies (mAbs). Optimal characterization of the human PK and target binding for mAbs requires data obtained after intravenous (IV) administration which can be combined with subcutaneous (SC) data to further this characterization. Integration of free and/or total target measurements in a population TMDD model will allow quantification of target engagement which is the first step in the cascade leading to efficacy. However, the assays for determination of free target concentrations are analytically challenging and are inherently biased to overpredict the true concentrations in the presence of mAb:target complexes. For that reason, the objective of the current research was to evaluate the predictive value of free target concentrations in a TMDD model developed using PK and total target observations only. Further, a secondary objective was to demonstrate that prediction of SC data is feasible, based on an existing IV model and typical values of mAb parameters reported for SC absorption. GSK3772847, a human immunoglobulin G2 sigma isotype (IgG2f) mAb that binds to the extracellular domain of the interleukin-33 receptor (IL-33R or ST2) and neutralizes IL-33-mediated ST2 signaling, was used as a model compound for mAbs in this study.

Artemisinin-based combination therapy (ACT) is the first-line recommended treatment for uncomplicated malaria. Pharmacokinetic (PK) properties in pregnant women are often based on small studies and need to be confirmed and validated in larger pregnant patient populations. This study aimed to evaluate the PK properties of amodiaquine and its active metabolite, desethylamodiaquine, and piperaquine in women in their second and third trimester of pregnancy with uncomplicated P. falciparum infections. Eligible pregnant women received either artesunate-amodiaquine (200/540 mg daily, n = 771) or dihydroartemisinin-piperaquine (40/960 mg daily, n = 755) for 3 days (NCT00852423). Population PK properties were evaluated using nonlinear mixed-effects modeling, and effect of gestational age and trimester was evaluated as covariates. 1071 amodiaquine and 1087 desethylamodiaquine plasma concentrations, and 976 piperaquine plasma concentrations, were included in the population PK analysis. Amodiaquine concentrations were described accurately with a one-compartment disposition model followed by a two-compartment disposition model of desethylamodiaquine. The relative bioavailability of amodiaquine increased with gestational age (1.25% per week). The predicted exposure to desethylamodiaquine was 2.8%–32.2% higher in pregnant women than that reported in non-pregnant women, while day 7 concentrations were comparable. Piperaquine concentrations were adequately described by a three-compartment disposition model. Neither gestational age nor trimester had significant impact on the PK of piperaquine. The predicted exposure and day 7 concentrations of piperaquine were similar to that reported in non-pregnant women. In conclusion, the exposure to desethylamodiaquine and piperaquine was similar to that in non-pregnant women. Dose adjustment is not warranted for women in their second and their trimester of pregnancy.

Omaveloxolone is a nuclear factor (erythroid-derived 2)-like 2 activator approved in the United States and the European Union for the treatment of patients with Friedreich ataxia aged ≥16 years, with a recommended dosage of 150 mg orally once daily on an empty stomach. The effect of the US Food and Drug Administration (FDA) high-fat breakfast on the pharmacokinetic profile of omaveloxolone observed in study 408-C-1703 (NCT03664453) deviated from the usual linear correlation between fed/fasted maximum plasma concentration (C_max) and area under the concentration–time curve (AUC) ratios reported for various oral drugs across 323 food effect studies. Here, physiologically based biopharmaceutics modeling (PBBM) was implemented to predict and explain the effect of the FDA high-fat breakfast on a 150-mg dose of omaveloxolone. The model was developed and validated based on dissolution and pharmacokinetic data available across dose-ranging, food effect, and drug–drug interaction clinical studies. PBBM predictions support clinical observations of the unique effect of a high-fat meal on omaveloxolone pharmacokinetic profile, in which the C_max increased by 350% with only a 15% increase in the AUC. Key parameters influencing omaveloxolone pharmacokinetics in the fasted state based on a parameter sensitivity analysis included bile salt solubilization, CYP3A4 activity, drug substance particle size distribution, and permeability. Mechanistically, in vivo omaveloxolone absorption was solubility and dissolution rate limited. However, in the fed state, higher bile salt solubilization led to more rapid dissolution with predominant absorption in the upper gastrointestinal tract, resulting in increased susceptibility to first-pass gut extraction; this accounts for the lack of correlation between C_max and AUC for omaveloxolone.

Certolizumab pegol (CZP; CIMZIA™) is the only Fc-free tumor necrosis factor inhibitor with data from a clinical study demonstrating no to minimal placental transfer. The pharmacokinetics (PK) of certolizumab pegol during pregnancy and postpartum in women with chronic inflammatory diseases were assessed using a population PK model based on data from the CHERISH study (NCT04163016), a longitudinal, prospective, open-label PK phase IB study. Model development was performed in NONMEM using a frequentist prior approach, with prior information based on a population PK model for certolizumab pegol in non-pregnant adult patients (NCT04740814). A one-compartment model with first-order absorption (K_a = 0.236 1/day) and linear elimination (CL/F = 0.416 L/day) from the central compartment (V/F = 7.86 L) best described certolizumab pegol PK in the CHERISH study. The structural model parameters were estimated with good precision (RSE < 25%). Baseline BW was included as a covariate on CL/F and V/F. Pregnancy trimester and time-varying log-transformed anti-drug antibody (ADA) titer were identified as the only significant covariates for CL/F with a comparable influence on CL/F. Individuals with higher ADA titer (75th percentile) during pregnancy exhibited CL/F up to 1.43-fold higher relative to individuals postpartum that showed median levels of ADA titer. However, the confidence interval for the combined effect of pregnancy stage and ADA titer effects on CL/F overlapped with the CL/F range of the typical individual postpartum. In addition, simulations showed a large overlap in certolizumab pegol concentrations between pregnant and non-pregnant adults. The findings of this population PK analysis support the maintenance of established certolizumab pegol dosing regimens throughout pregnancy.

Lecanemab (Leqembi®) was recently approved by health authorities in the United States, Japan, and China to treat early Alzheimer's disease (AD), including patients with mild cognitive impairment (MCI) or mild dementia due to Alzheimer's disease upon confirmation of amyloid beta pathology. Extensively and sparsely sampled PK profiles from 1619 AD subjects and 21,929 serum lecanemab observations from two phase I, one phase II, and one phase III studies were well characterized using a two-compartment model with first-order elimination. The final PK model quantified covariate effects of body weight and sex on clearance and central volume of distribution, ADA-positive status, and albumin on clearance, and of Japanese ethnicity on central and peripheral volumes of distribution. Exposure to lecanemab was comparable between two lecanemab-manufacturing processes. However, none of the identified covariates in the model had a clinically relevant impact on model-predicted lecanemab C_max or AUC at steady state following 10 mg/kg bi-weekly. Importantly, age, a well-recognized risk factor for AD, was not found to significantly affect lecanemab PK. The incidence of ARIA-E as a function of lecanemab exposure was modeled using a logit function with data pooled from 2641 subjects from the phase II and phase III studies, in which a total of 177 incidences of ARIA-E were observed. The probability of ARIA-E was significantly correlated with model-predicted C_max and predicted to be higher in subjects homozygous for APOE4. The incidence of isolated ARIA-H was not associated with lecanemab exposure and was similar between placebo and lecanemab-treated subjects.

Conventional approaches for establishing bioequivalence (BE) between test and reference formulations using non-compartmental analysis (NCA) may demonstrate low power in pharmacokinetic (PK) studies with sparse sampling. In this case, model-integrated evidence (MIE) approaches for BE assessment have been shown to increase power, but may suffer from selection bias problems if models are built on the same data used for BE assessment. This work presents model averaging methods for BE evaluation and compares the power and type I error of these methods to conventional BE approaches for simulated studies of oral and ophthalmic formulations. Two model averaging methods were examined: bootstrap model selection and weight-based model averaging with parameter uncertainty from three different sources, either from a sandwich covariance matrix, a bootstrap, or from sampling importance resampling (SIR). The proposed approaches increased power compared with conventional NCA-based BE approaches, especially for the ophthalmic formulation scenarios, and were simultaneously able to adequately control type I error. In the rich sampling scenario considered for oral formulation, the weight-based model averaging method with SIR uncertainty provided controlled type I error, that was closest to the target of 5%. In sparse-sampling designs, especially the single sample ophthalmic scenarios, the type I error was best controlled by the bootstrap model selection method.

Population pharmacokinetic (PK) models are widely used to inform drug development by pharmaceutical companies and facilitate drug evaluation by regulatory agencies. Developing a population PK model is a multi-step, challenging, and time-consuming process involving iterative manual model fitting and evaluation. A tool for fully automatic model development (AMD) of common population PK models is presented here. The AMD tool is implemented in Pharmpy, a versatile open-source library for pharmacometrics. It consists of different modules responsible for developing the different components of population PK models, including the structural model, the inter-individual variability (IIV) model, the inter-occasional variability (IOV) model, the residual unexplained variability (RUV) model, the covariate model, and the allometry model. The AMD tool was evaluated using 10 real PK datasets involving the structural, IIV, and RUV modules in three sequences. The different sequences yielded generally consistent structural models; however, there were variations in the results of the IIV and RUV models. The final models of the AMD tool showed lower Bayesian Information Criterion (BIC) values and similar visual predictive check plots compared with the available published models, indicating reasonable quality, in addition to reasonable run time. A similar conclusion was also drawn in a simulation study. The developed AMD tool serves as a promising tool for fast and fully automatic population PK model building with the potential to facilitate the use of modeling and simulation in drug development.

This study addresses the critical issue of drug-induced torsades de pointes (TdP) risk assessment, a vital aspect of new drug development due to its association with arrhythmia and sudden cardiac death. Existing methodologies, particularly those reliant on a single biomarker derived from CiPA O'Hara-Rudy (CiPAORdv1.0) ventricular cell model without the hERG dynamic as input to the individual machine learning model, have limitations in capturing the complexity inherent in the comprehensive range of factors influencing drug-induced TdP risk. This study aims to overcome these limitations by proposing a stacking ensemble machine learning approach by integrating multiple in silico biomarkers derived from the CiPAORdv1.0 with hERG dynamic characteristics. The ensemble machine learning model consisted of three artificial neural network (ANN) models as baseline model and support vector machine (SVM), logistic regression (LR), random forest (RF), and extreme gradient boosting (XGBoost) models as meta-classifier. The highest AUC score of 1.00 (0.90–1.00) for high risk, 0.97 (0.84–1.00) for intermediate risk, and 1.00 (0.87–1.00) for low risk were obtained using seven biomarkers derived from the CiPAORdv1.0 with hERG dynamic characteristics. Furthering our investigation, we explored the model's robustness by incorporating interindividual variability into the generation of in silico biomarkers from a population of human ventricular cell models. This study also enabled an analysis of TdP risk classification under high clinical exposure and therapeutic scenarios for several drugs. Additionally, from a sensitivity analysis, we revealed four important ion channels, namely, CaL, NaL, Na, and Kr channels that affect significantly the important biomarkers for TdP risk prediction.

By applying nonlinear mixed-effect (NLME) models, model-integrated evidence (MIE) approaches are able to analyze bioequivalence (BE) data with pharmacokinetic end points that have sparse sampling, which is problematic for non-compartmental analysis (NCA). However, MIE approaches may suffer from inflation of type I error due to underestimation of parameter uncertainty and to the assumption of asymptotic normality. In this study, we developed a MIE BE analysis method that is based on a pre-defined model and consists of several steps including model fitting, uncertainty assessment, simulation, and BE determination. The presented MIE approach has several improvements compared with the previously reported model-integrated methods: (1) treatment, sequence, and period effects are only added to absorption parameters (such as relative bioavailability and rate of absorption) instead of all PK parameters; (2) a simulation step is performed to generate confidence intervals of the pharmacokinetic metrics for BE assessment; and (3) in an effort to maintain type I error, two more advanced parameter uncertainty evaluation approaches are explored, a nonparametric (case resampling) bootstrap, and sampling importance resampling (SIR). To evaluate the developed method and compare the uncertainty assessment methods, simulation experiments were performed for BE studies using a two-way crossover design with different amounts of information (sparse to rich designs) and levels of variability. Based on the simulation results, the method using SIR for parameter uncertainty quantification controls type I error at the nominal level of 0.05 (i.e., the significance level set for BE evaluation) even for studies with small sample size and/or sparse sampling. As expected, our MIE approach for BE assessment exhibited higher power than the NCA-based method, especially as the data becomes sparser and/or more variable.