Uremia, a condition characterized by the retention of uremic toxins due to impaired renal function, may affect drug metabolism mediated by CYP3A4 enzymes. Evogliptin is a dipeptidyl peptidase-4 (DPP-4) inhibitor diabetic drug that is primarily metabolized by CYP3A4. This study aimed to construct a population pharmacokinetic (PK) and pharmacodynamic (PD) model for evogliptin in patients with varying degrees of renal disease, including end-stage renal disease on hemodialysis. A total of 688 evogliptin concentration and 598 DPP-4 activity data were available from 46 subjects. PK and PD data analyses were performed using a nonlinear mixed-effects model. The PK of evogliptin was optimally described by a two-compartment model with first-order absorption. The significant covariates in the final model included blood amylase and triglyceride on F1 (relative bioavailability). The simulation findings, together with previously reported PK data, provided evidence of a significant inhibition of the first-pass effect of evogliptin in patients with renal impairment. A direct link sigmoidal Emax model was developed to describe the relationship between evogliptin concentration and DPP-4 inhibition. The PD model predicted significant inhibition of DPP-4 at maximum effect (Emax: 88.9%) and a low EC50 value (1.08 μg/L), indicating the high potency and efficacy of evogliptin. The developed PK/PD model accurately predicted exposure and the resulting DPP-4 activity of evogliptin in renal impairment. The findings of this study suggest that renal impairment and associated biochemical changes may impact the bioavailability of CYP3A4-metabolized drugs.
Stan is a powerful probabilistic programming language designed mainly for Bayesian data analysis. Torsten is a collection of Stan functions that handles the events (e.g., dosing events) and solves the ODE systems that are frequently present in pharmacometric models. To perform a Bayesian data analysis, most models in pharmacometrics require Markov Chain Monte Carlo (MCMC) methods to sample from the posterior distribution. However, MCMC is computationally expensive and can be time-consuming, enough so that people will often forgo Bayesian methods for a more traditional approach. This paper shows how to speed up the sampling process in Stan by within-chain parallelization through both multi-threading using Stan's reduce_sum() function and multi-processing using Torsten's group ODE solver. Both methods show substantial reductions in the time necessary to sufficiently sample from the posterior distribution compared with a basic approach with no within-chain parallelization.
Breastfeeding is important in childhood development, and medications are often necessary for lactating individuals, yet information on the potential risk of infant drug exposure through human milk is limited. Establishing a lactation modeling framework can advance our understanding of this topic and potentiate clinical decision making. We expanded the modeling framework previously developed for sotalol using pregabalin as a second prototypical probe compound with similar absorption, distribution, metabolism, and elimination (ADME) properties. Adult oral models were developed in PK-Sim® and used to build a lactation model in MoBi® to simulate drug transfer into human milk. The adult model was applied to breastfeeding pediatrics (ages 1 to 23 months) and subsequently integrated with the lactation model to simulate infant drug exposure according to age, size, and breastfeeding frequency. Physiologically based pharmacokinetic (PBPK) model simulations captured the data used for verification both in adults and pediatrics. Lactation simulations captured observed milk and plasma data corresponding to doses of 150 mg administered twice daily to lactating individuals, and estimated a relative infant dose (RID) of approximately 7% of the maternal dose. The infant drug exposure simulations showed peak plasma concentrations of 0.44 μg/mL occurring within the first 2 weeks of life, followed by gradual decline with age after week four. The modeling framework performs well for this second prototypical drug and warrants expansion to other drugs for further validation. PBPK modeling and simulation approaches together with clinical lactation data could ultimately help inform infant drug exposure risk assessments to guide clinical decision making.
The ICH E9 (R1) guidance and the related estimand framework propose to clearly define and separate the clinical question of interest formulated as estimand from the estimation method. With that it becomes important to assess the validity of the estimation method and the assumptions that must be made. When going beyond the intention to treat analyses that can rely on randomization, causal inference is usually used to discuss the validity of estimation methods for the estimand of interest. In pharmacometrics, mixed-effects models are routinely used to analyze longitudinal clinical trial data; however, they are rarely discussed as a method for causal inference. Here, we evaluate nonlinear mixed-effects modeling and simulation (NLME M&S) in the context of causal inference as a standardization method for longitudinal data in the presence of confounders. Standardization is a well-known method in causal inference to correct for confounding by analyzing and combining results from subgroups of patients. We show that nonlinear mixed-effects modeling is a particular implementation of standardization that conditions on individual parameters described by the random effects of the mixed-effects model. As an example, we use a simulated clinical trial with within-subject dose titration. Being interested in the outcome of the hypothetical situation that patients adhere to the planned treatment schedule, we put assumptions in a causal diagram. From the causal diagram, conditional independence assumptions are derived either by conditioning on the individual parameters or on earlier outcomes. With both conditional independencies unbiased estimates can be obtained.
Aldehyde oxidase (AO) contributes to the clearance of many approved and investigational small molecule drugs, which are often dual substrates of AO and drug-metabolizing enzymes such as cytochrome P450s (CYPs). As such, the lack of established framework for quantitative translation of the clinical pharmacologic correlates of AO-mediated clearance represents an unmet need. This study aimed to evaluate the utility of physiologically based pharmacokinetic (PBPK) modeling in the development of AO and dual AO-CYP substrates. PBPK models were developed for capmatinib, idelalisib, lenvatinib, zaleplon, ziprasidone, and zoniporide, incorporating in vitro functional data from human liver subcellular fractions and human hepatocytes. Prediction of metabolic elimination with/without the additional empirical scaling factors (ESFs) was assessed. Clinical pharmacokinetics, human mass balance, and drug–drug interaction (DDI) studies with CYP3A4 modulators, where available, were used to refine/verify the models. Due to the lack of clinically significant AO-DDIs with known AO inhibitors, the fraction metabolized by AO (fmAO) was verified indirectly. Clearance predictions were improved by using ESFs (GMFE ≤1.4-fold versus up to fivefold with physiologically-based scaling only). Observed fmi from mass balance studies were crucial for model verification/refinement, as illustrated by capmatinib, where the fmAO (40%) was otherwise underpredicted up to fourfold. Subsequently, independent DDI studies with ketoconazole, itraconazole, rifampicin, and carbamazepine verified the fmCYP3A4, with predicted ratios of the area under the concentration–time curve (AUCR) within 1.5-fold of the observations. In conclusion, this study provides a novel PBPK-based framework for predicting AO-mediated pharmacokinetics and quantitative assessment of clinical DDI risks for dual AO-CYP substrates within a totality-of-evidence approach.
The number of quantitative systems pharmacology (QSP) submissions to the U.S. Food and Drug Administration has continued to increase over the past decade. This report summarizes the landscape of QSP submissions as of December 2023. QSP was used to inform drug development across various therapeutic areas and throughout the drug development process of small molecular drugs and biologics and has facilitated dose finding, dose ranging, and dose optimization studies. Though the majority of QSP submissions (>66%) focused on drug effectiveness, QSP was also utilized to simulate drug safety including liver toxicity, risk of cytokine release syndrome (CRS), bone density, and others. This report also includes individual contexts of use from a handful of new drug applications (NDAs) and biologics license applications where QSP modeling was used to demonstrate the utility of QSP modeling in regulatory drug development. According to the models submitted in QSP submissions, an anonymous case was utilized to illustrate how QSP informed development of a bispecific monoclonal antibody with respect to CRS risk. QSP submissions for informing pediatric drug development were summarized along with highlights of a case in inborn errors of metabolism. Furthermore, simulations of response variability with QSP were described. In summary, QSP continues to play a role in informing drug development.
The cytokine interleukin-12 (IL-12) is a potential immunotherapy because of its ability to induce a Th1 immune response. However, success in the clinic has been limited due to a phenomenon called IL-12 desensitization – the trend where repeated exposure to IL-12 leads to reduced IL-12 concentrations (pharmacokinetics) and biological effects (pharmacodynamics). Here, we investigated IL-12 pharmacokinetic desensitization via a modeling approach to (i) validate proposed mechanisms in literature and (ii) develop a mathematical model capable of predicting IL-12 pharmacokinetic desensitization. Two potential causes of IL-12 pharmacokinetic desensitization were identified: increased clearance or reduced bioavailability of IL-12 following repeated doses. Increased IL-12 clearance was previously proposed to occur due to the upregulation of IL-12 receptor on T cells that causes increased receptor-mediated clearance in the serum. However, our model with this mechanism, the accelerated-clearance model, failed to capture trends in clinical trial data. Alternatively, our novel reduced-bioavailability model assumed that upregulation of IL-12 receptor on T cells in the lymphatic system leads to IL-12 sequestration, inhibiting the transport to the blood. This model accurately fits IL-12 pharmacokinetic data from three clinical trials, supporting its biological relevance. Using this model, we analyzed the model parameter space to illustrate that IL-12 desensitization occurs over a robust range of parameter values and to identify the conditions required for desensitization. We next simulated local, continuous IL-12 delivery and identified several methods to mitigate systemic IL-12 exposure. Ultimately, our results provide quantitative validation of our proposed mechanism and allow for accurate prediction of IL-12 pharmacokinetics over repeated doses.
The use of synthetic data in pharmacology research has gained significant attention due to its potential to address privacy concerns and promote open science. In this study, we implemented and compared three synthetic data generation methods, CT-GAN, TVAE, and a simplified implementation of Avatar, for a previously published pharmacogenetic dataset of 253 patients with one measurement per patient (non-longitudinal). The aim of this study was to evaluate the performance of these methods in terms of data utility and privacy trade off. Our results showed that CT-GAN and Avatar used with k = 10 (number of patients used to create the local model of generation) had the best overall performance in terms of data utility and privacy preservation. However, the TVAE method showed a relatively lower level of performance in these aspects. In terms of Hazard ratio estimation, Avatar with k = 10 produced HR estimates closest to the original data, whereas CT-GAN slightly underestimated the HR and TVAE showed the most significant deviation from the original HR. We also investigated the effect of applying the algorithms multiple times to improve results stability in terms of HR estimation. Our findings suggested that this approach could be beneficial, especially in the case of small datasets, to achieve more reliable and robust results. In conclusion, our study provides valuable insights into the performance of CT-GAN, TVAE, and Avatar methods for synthetic data generation in pharmacogenetic research. The application to other type of data and analyses (data driven) used in pharmacology should be further investigated.
Quantitative systems toxicology (QST) models are increasingly being applied for predicting and understanding toxicity liabilities in pharmaceutical research and development. A European Federation of Pharmaceutical Industries and Associations (EFPIA)-wide survey was completed by 15 companies. The results provide insights into the current use of QST models across the industry. 73% of responding companies with more than 10,000 employees utilize QST models. The most applied QST models are for liver, cardiac electrophysiology, and bone marrow/hematology. Responders indicated particular interest in QST models for the central nervous system (CNS), kidney, lung, and skin. QST models are used to support decisions in both preclinical and clinical stages of pharmaceutical development. The survey suggests high demand for QST models and resource limitations were indicated as a common obstacle to broader use and impact. Increased investment in QST resources and training may accelerate application and impact. Case studies of QST model use in decision-making within EFPIA companies are also discussed. This article aims to (i) share industry experience and learnings from applying QST models to inform decision-making in drug discovery and development programs, and (ii) share approaches taken during QST model development and validation and compare these with recommendations for modeling best practices and frameworks proposed in the literature. Discussion of QST-specific applications in relation to these modeling frameworks is relevant in the context of the recently proposed International Council for Harmonization (ICH) M15 guideline on general principles for Model-Informed Drug Development (MIDD).