Journal of Pharmacokinetics and Pharmacodynamics最新文献

In this paper, we use Time Scale Calculus (TSC) to formulate and solve pharmacokinetic models exploring multiple dose dynamics. TSC is a mathematical framework that allows the modeling of dynamical systems comprising continuous and discrete processes. This characteristic makes TSC particularly suited for multi-dose pharmacokinetic problems, which inherently feature a blend of continuous processes (such as absorption, metabolization, and elimination) and discrete events (drug intake). We use this toolkit to derive analytical expressions for blood concentration trajectories under various multi-dose regimens across several flagship pharmacokinetic models. We demonstrate that this mathematical framework furnishes an alternative and simplified way to model and retrieve analytical solutions for multi-dose dynamics. For instance, it enables the study of blood concentration responses to arbitrary dose regimens and facilitates the characterization of the long-term behavior of the solutions, such as their steady state.

Congenital adrenal hyperplasia (CAH) is characterized by impaired adrenal cortisol production. Hydrocortisone (synthetic cortisol) is the drug-of-choice for cortisol replacement therapy, aiming to mimic physiological cortisol circadian rhythm. The hypothalamic-pituitary-adrenal (HPA) axis controls cortisol production through the pituitary adrenocorticotropic hormone (ACTH) and feedback mechanisms. The aim of this study was to quantify key mechanisms involved in the HPA axis activity regulation and their interaction with hydrocortisone therapy. Data from 30 healthy volunteers was leveraged: Endogenous ACTH and cortisol concentrations without any intervention as well as cortisol concentrations measured after dexamethasone suppression and single dose administration of (i) 0.5-10 mg hydrocortisone as granules, (ii) 20 mg hydrocortisone as granules and intravenous bolus. A stepwise model development workflow was used: A newly developed model for endogenous ACTH and cortisol was merged with a refined hydrocortisone pharmacokinetic model. The joint model was used to simulate ACTH and cortisol trajectories in CAH patients with varying degrees of enzyme deficiency, with or without hydrocortisone administration, and healthy individuals. Time-dependent ACTH-driven endogenous cortisol production and cortisol-mediated feedback inhibition of ACTH secretion processes were quantified and implemented in the model. Comparison of simulated ACTH and cortisol trajectories between CAH patients and healthy individuals showed the importance of administering hydrocortisone before morning ACTH secretion peak time to suppress ACTH overproduction observed in untreated CAH patients. The developed framework allowed to gain insights on the physiological mechanisms of the HPA axis regulation, its perturbations in CAH and interaction with hydrocortisone administration, paving the way towards cortisol replacement therapy optimization.

Forward addition/backward elimination (FABE) has been the standard for population pharmacokinetic model selection (PPK) since NONMEM® was introduced. We investigated five machine learning (ML) algorithms (Genetic algorithm [GA], Gaussian process [GP], random forest [RF], gradient boosted random tree [GBRT], and particle swarm optimization [PSO]) as alternatives to FABE. These algorithms were applied to PPK model selection with a focus on comparing the efficiency and robustness of each of them. All machine learning algorithms included the combination of ML algorithms with a local downhill search. The local downhill search consisted of systematically changing one or two "features" at a time (a one-bit or a two-bit local search), alternating with the ML methods. An exhaustive search (all possible combinations of model features, N = 1,572,864 models) was the gold standard for robustness, and the number of models examined leading prior to identification of the final model was the metric for efficiency.All algorithms identified the optimal model when combined with the two-bit local downhill search. GA, RF, GBRT, and GP identified the optimal model with only a one-bit local search. PSO required the two-bit local downhill search. In our analysis, GP was the most efficient algorithm as measured by the number of models examined prior to finding the optimal (495 models), and PSO exhibited the least efficiency, requiring 1710 unique models before finding the best solution. Additionally, GP was also the algorithm that needed the longest elapsed time of 2975.6 min, in comparison with GA, which only required 321.8 min.

Physiologically based pharmacokinetic (PBPK) models can be used to leverage physiological and in vitro data to predict monoclonal antibody (mAb) concentrations in serum and tissues. However, it is currently not known how consistent predictions of mAb disposition are across PBPK modelling platforms. In this work PBPK simulations of IgG, adalimumab and infliximab were compared between three platforms (Simcyp, PK-Sim, and GastroPlus). Accuracy of predicted serum and tissue concentrations was assessed using observed data collected from the literature. Physiological and mAb related input parameters were also compared and sensitivity analyses were carried out to evaluate model behavior when input values were altered. Differences in serum kinetics of IgG between platforms were minimal for a dose of 1 mg/kg, but became more noticeable at higher dosages (> 100 mg/kg) and when reference (healthy) physiological input values were altered. Predicted serum concentrations of both adalimumab and infliximab were comparable across platforms, but were noticeably higher than observed values. Tissue concentrations differed remarkably between the platforms, both for total- and interstitial fluid (ISF) concentrations. The accuracy of total tissue concentrations was within a three-fold of observed values for all tissues, except for brain tissue concentrations, which were overpredicted. Predictions of tissue ISF concentrations were less accurate and were best captured by GastroPlus. Overall, these simulations show that the different PBPK platforms generally predict similar mAb serum concentrations, but variable tissue concentrations. Caution is therefore warranted when PBPK models are used to simulate effect site tissue concentrations of mAbs without data to verify the predictions.

Pharmacokinetics and pharmacodynamics of many biologics are influenced by their complex binding to biological receptors. Biologics consist of diverse groups of molecules with different binding kinetics to its receptors including IgG with simple one-to-one drug receptor bindings, bispecific antibody (BsAb) that binds to two different receptors, and antibodies that can bind to six or more identical receptors. As the binding process is typically much faster than elimination (or internalization) and distribution processes, quasi-equilibrium (QE) binding models are commonly used to describe drug-receptor binding kinetics of biologics. However, no general QE modeling framework is available to describe complex binding kinetics for diverse classes of biologics. In this paper, we describe novel approaches of using differential algebraic equations (DAE) to solve three QE multivalent drug-receptor binding (QEMB) models. The first example describes the binding kinetics of three-body equilibria of BsAb that binds to 2 different receptors for trimer formation. The second example models an engineered IgG variant (Multabody) that can bind to 24 identical target receptors. The third example describes an IgG with modified neonatal Fc receptor (FcRn) binding affinity that competes for the same FcRn receptor as endogenous IgG. The model parameter estimates were obtained by fitting the model to all data simultaneously. The models allowed us to study potential roles of cooperative binding on bell-shaped drug exposure-response relationships of BsAb, and concentration-depended distribution of different drug-receptor complexes for Multabody. This DAE-based QEMB model platform can serve as an important tool to better understand complex binding kinetics of diverse classes of biologics.

Incorporating realistic sets of patient-associated covariates, i.e., virtual populations, in pharmacometric simulation workflows is essential to obtain realistic model predictions. Current covariate simulation strategies often omit or simplify dependency structures between covariates. Copula models are multivariate distribution functions suitable to capture dependency structures between covariates with improved performance compared to standard approaches. We aimed to develop and evaluate a copula model for generation of adult virtual populations for 12 patient-associated covariates commonly used in pharmacometric simulations, using the publicly available NHANES database, including sex, race-ethnicity, body weight, albumin, and several biochemical variables related to organ function. A multivariate (vine) copula was constructed from bivariate relationships in a stepwise fashion. Covariate distributions were well captured for the overall and subgroup populations. Based on the developed copula model, a web application was developed. The developed copula model and associated web application can be used to generate realistic adult virtual populations, ultimately to support model-based clinical trial design or dose optimization strategies.

Hereditary angioedema (HAE) due to C1-inhibitor deficiency is a rare, debilitating, genetic disorder characterized by recurrent, unpredictable, attacks of edema. The clinical symptoms of HAE arise from excess bradykinin generation due to dysregulation of the plasma kallikrein-kinin system (KKS). A quantitative systems pharmacology (QSP) model that mechanistically describes the KKS and its role in HAE pathophysiology was developed based on HAE attacks being triggered by autoactivation of factor XII (FXII) to activated FXII (FXIIa), resulting in kallikrein production from prekallikrein. A base pharmacodynamic model was constructed and parameterized from literature data and ex vivo assays measuring inhibition of kallikrein activity in plasma of HAE patients or healthy volunteers who received lanadelumab. HAE attacks were simulated using a virtual patient population, with attacks recorded when systemic bradykinin levels exceeded 20 pM. The model was validated by comparing the simulations to observations from lanadelumab and plasma-derived C1-inhibitor clinical trials. The model was then applied to analyze the impact of nonadherence to a daily oral preventive therapy; simulations showed a correlation between the number of missed doses per month and reduced drug effectiveness. The impact of reducing lanadelumab dosing frequency from 300 mg every 2 weeks (Q2W) to every 4 weeks (Q4W) was also examined and showed that while attack rates with Q4W dosing were substantially reduced, the extent of reduction was greater with Q2W dosing. Overall, the QSP model showed good agreement with clinical data and could be used for hypothesis testing and outcome predictions.

Pozelimab, a monoclonal antibody directed against C5, is the first and only treatment for adult and pediatric patients (≥ 1 year) with CD55-deficient protein-losing enteropathy (CHAPLE) disease. A target-mediated drug disposition (TMDD) population pharmacokinetic (PopPK) model was developed using pooled data from four phase 1-3 studies to characterize the pharmacokinetics (PK) of total pozelimab and total C5, and to simulate free pozelimab and free C5 to support the dose regimen in patients with CHAPLE disease. A TMDD PopPK model was developed using total pozelimab and total C5 concentration-time data from 106 participants (82 healthy volunteers; 24 patients with paroxysmal nocturnal hemoglobinuria [PNH]). This model was refined and updated to include PK data from 10 patients with CHAPLE disease from a phase 2/3 study. Stochastic simulations predicted concentration-time profiles for total pozelimab, free pozelimab, and free C5, to obtain pozelimab exposure metrics for patients with CHAPLE disease. A two-compartment TMDD model with two binding sites based on the quasi-equilibrium approximation adequately described the concentration-time profiles of total pozelimab and total C5. Body weight was identified as the most important source of pozelimab PK variability; therefore, the dose was adjusted based on body weight for the predominantly pediatric patients with CHAPLE disease. A robust TMDD PopPK model was developed to describe the PK of total pozelimab and total C5 following pozelimab administration. Reliable predictions for individual exposures of total pozelimab and free C5 were possible and supported the 10 mg/kg weight-based dose regimen in patients with CHAPLE disease.

Clinical trial endpoints are often bounded outcome scores (BOS), which are variables having restricted values within finite intervals. Common analysis approaches may treat the data as continuous, categorical, or a mixture of both. The appearance of BOS data being simultaneously continuous and categorical easily leads to confusions in pharmacometrics regarding the appropriate domain for model evaluation and the circumstances under which data likelihoods can be compared. This commentary aims to clarify these fundamental issues and facilitate appropriate pharmacometric analyses.

The generation of synthetic patient data that reflect the statistical properties of real data plays a fundamental role in today's world because of its potential to (i) be enable proprietary data access for statistical and research purposes and (ii) increase available data (e.g., in low-density regions-i.e., for patients with under-represented characteristics). Generative methods employ a family of solutions for generating synthetic data. The objective of this research is to benchmark numerous state-of-the-art deep-learning generative methods across different scenarios and clinical datasets comprising patient covariates and several pharmacokinetic/pharmacodynamic endpoints. We did this by implementing various probabilistic models aimed at generating synthetic data, such as the Multi-layer Perceptron Conditioning Generative Adversarial Neural Network (MLP cGAN), Time-series Generative Adversarial Networks (TimeGAN), and a more traditional approach like Probabilistic Autoregressive (PAR). We evaluated their performance by calculating discriminative and predictive scores. Furthermore, we conducted comparisons between the distributions of real and synthetic data using Kolmogorov-Smirnov and Chi-square statistical tests, focusing respectively on covariate and output variables of the models. Lastly, we employed pharmacometrics-related metric to enhance interpretation of our results specific to our investigated scenarios. Results indicate that multi-layer perceptron-based conditional generative adversarial networks (MLP cGAN) exhibit the best overall performance for most of the considered metrics. This work highlights the opportunities to employ synthetic data generation in the field of clinical pharmacology for augmentation and sharing of proprietary data across institutions.