Tegoprazan is a novel potassium-competitive acid blocker (P-CAB) developed by CJ Healthcare (Korea) for the treatment of gastroesophageal reflux disease and helicobacter pylori infections. Tegoprazan is mainly metabolized by cytochrome P450 (CYP) 3A4. Considering the therapeutic indications, tegoprazan is likely to be administered in combination with various drugs. Therefore, the investigation of drug–drug interactions (DDI) between tegoprazan and CYP3A4 perpetrators is imperative. In the present study, we first aimed to develop a physiologically based pharmacokinetic (PK) model for tegoprazan and its major metabolite, M1, using PK-Sim^®. This model was applied to predict the DDI between tegoprazan and CYP3A4 perpetrators. Clarithromycin, a potent inhibitor of CYP3A4, and rifampicin, a strong inducer of CYP3A4, were selected as case studies. Our results show that clarithromycin significantly increased the exposure of tegoprazan. The area under the concentration–time curve (AUC) and C_max of tegoprazan in the steady state increased up to 4.54- and 2.05-fold, respectively, when tegoprazan (50 mg, twice daily) was coadministered with clarithromycin (500 mg, three times daily). Rifampicin significantly reduced the exposure of tegoprazan. The AUC and C_max of tegoprazan were reduced by 5.71- and 3.51-fold when tegoprazan was coadministered with rifampicin (600 mg, once daily). Due to the high DDI potential, the comedication of tegoprazan with CYP3A4 perpetrators should be controlled. The dosage adjustment for each individual is suggested.

Although alterations in the methionine metabolism cycle (MMC) have been associated with vascular complications of diabetes, there have not been consistent results about the levels of methionine and homocysteine in type 2 diabetes mellitus (T2DM). The aim of the current study was to predict changes in plasma methionine and homocysteine concentrations after simulated consumption of methionine-rich foods, following the development of a mathematical model for MMC in Zucker Diabetic Fatty (ZDF) rats, as a representative T2DM animal model. The model building and simulation were performed using NONMEM® (ver. 7.3.0) assisted by Perl-Speaks-NONMEM (PsN, ver. 4.3.0). Model parameters were derived using first-order conditional estimation method with interactions permitted among the parameters (FOCE-INTER). NCA was conducted using Phoenix (ver. 6.4.0). For all tests, we considered a P-value < 0.05 to reflect statistical significance. Our model featured seven compartments that considered all parts of the cycle by applying non-linear mixed effects model. Conversion of S-adenosyl-L-homocysteine (SAH) to homocysteine increased and the metabolism of homocysteine was reduced under diabetic conditions, and consequently homocysteine accumulated in the elimination phase. Using our model, we performed simulations to compare the changes in plasma methionine and homocysteine concentrations between ZDF and normal rats, by multiple administrations of the methionine-rich diet of 1 mmol/kg, daily for 60 days. The levels of methionine and homocysteine were elevated approximately two- and three-fold, respectively, in ZDF rats, while there were no changes observed in the normal control rats. These results can be interpreted to mean that both methionine and homocysteine will accumulate in patients with T2DM, who regularly consume high-methionine foods.

Exploratory preclinical, as well as clinical trials, may involve a small number of patients, making it difficult to calculate and analyze the pharmacokinetic (PK) parameters, especially if the PK parameters show very high inter-individual variability (IIV). In this study, the performance of a classical first-order conditional estimation with interaction (FOCE-I) and expectation maximization (EM)-based Markov chain Monte Carlo Bayesian (BAYES) estimation methods were compared for estimating the population parameters and its distribution from data sets having a low number of subjects. In this study, 100 data sets were simulated with eight sampling points for each subject and with six different levels of IIV (5%, 10%, 20%, 30%, 50%, and 80%) in their PK parameter distribution. A stochastic simulation and estimation (SSE) study was performed to simultaneously simulate data sets and estimate the parameters using four different methods: FOCE-I only, BAYES(C) (FOCE-I and BAYES composite method), BAYES(F) (BAYES with all true initial parameters and fixed ω 2 ), and BAYES only. Relative root mean squared error (rRMSE) and relative estimation error (REE) were used to analyze the differences between true and estimated values. A case study was performed with a clinical data of theophylline available in NONMEM distribution media. NONMEM software assisted by Pirana, PsN, and Xpose was used to estimate population PK parameters, and R program was used to analyze and plot the results. The rRMSE and REE values of all parameter (fixed effect and random effect) estimates showed that all four methods performed equally at the lower IIV levels, while the FOCE-I method performed better than other EM-based methods at higher IIV levels (greater than 30%). In general, estimates of random-effect parameters showed significant bias and imprecision, irrespective of the estimation method used and the level of IIV. Similar performance of the estimation methods was observed with theophylline dataset. The classical FOCE-I method appeared to estimate the PK parameters more reliably than the BAYES method when using a simple model and data containing only a few subjects. EM-based estimation methods can be considered for adapting to the specific needs of a modeling project at later steps of modeling.

Sample sizes for single-period clinical trials, including pharmacokinetic studies, are statistically determined by within-subject variability (WSV). However, it is difficult to determine WSV without replicate-designed clinical trial data, and statisticians typically estimate optimal sample sizes using total variability, not WSV. We have developed an efficient population-based method to predict WSV accurately with single-period clinical trial data and demonstrate method performance with eperisone. We simulated 1000 virtual pharmacokinetic clinical trial datasets based on single-period and dense sampling studies, with various study sizes and levels of WSV and interindividual variabilities (IIVs). The estimated residual variability (RV) resulting from population pharmacokinetic methods were compared with WSV values. In addition, 3 × 3 bioequivalence results of eperisone were used to evaluate method performance with a real clinical dataset. With WSV of 40% or less, regardless of IIV magnitude, RV was well approximated by WSV for sample sizes greater than 18 subjects. RV was underestimated at WSV of 50% or greater, even with datasets having low IIV and numerous subjects. Using the eperisone dataset, RV was 44% to 48%, close to the true value of 50%. In conclusion, the estimated RV accurately predicted WSV in single-period studies, validating this method for sample size estimation in clinical trials.

Potassium‐competitive acid blockers (PCABs) are emerging alternatives to proton pump inhibitors for the treatment of acid‐related diseases. However, due to the complex, nonlinear interaction between drug exposure, food intake, and physiological rhythms, optimizing dosing strategies remains challenging. A multi‐leveled population analysis was conducted using published pharmacokinetic and pharmacodynamic data on four representative PCABs: tegoprazan, YH4808, fexuprazan, and vonoprazan. A semi‐mechanistic population PK/PD model was developed to account for food effects, circadian pH rhythms, and pH‐dependent drug absorption. A multi‐level nonlinear mixed‐effects modeling framework was implemented to capture both inter‐drug and inter‐study variability. The model successfully described the time course of plasma concentration and intragastric pH for all four PCABs under various conditions. The model identified differences in pharmacokinetics and pharmacodynamic potency between drugs (with the relative in vitro potency ranked as vonoprazan > fexuprazan > YH4808 > tegoprazan), and simulations demonstrated that both pre‐ and post‐meal administration enhanced pH control in early time period via potentially distinct mechanisms: the pre‐meal effect may arise from temporally separated contributions of food‐ and drug‐induced pH elevation, whereas the post‐meal effect is likely driven by temporally overlapping, additive actions, particularly under low‐dose or non‐steady‐state conditions. Predicted pH profiles and holding times above pH 4 closely matched reported clinical outcomes. The study demonstrates the application of a mechanistic, multi‐level population approach for cross‐drug PK/PD evaluation of PCABs. The findings support drug‐specific dose optimization and highlight the clinical relevance of food–drug interactions. The modeling approach provides a model platform for pharmacotherapy or model‐informed drug development (MIDD).