Journal of Pharmacokinetics and Pharmacodynamics最新文献

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.

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.

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.

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.

Visual predictive checks (VPC) are commonly used to evaluate pharmacometrics models. However their performance may be hampered if patients with worse outcomes drop out earlier, as often occurs in clinical trials, especially in oncology. While methods accounting for dropouts have appeared in literature, they vary in assumptions, flexibility, and performance, and the differences between them are not widely understood. This manuscript aims to elucidate which methods can be used to handle VPC with dropout and when, along with a more informative VPC approach using confidence intervals. Additionally, we propose constructing the confidence interval based on the observed data instead of the simulated data. The theoretical framework for incorporating dropout in VPCs is developed and applied to propose two approaches: full and conditional. The full approach is implemented using a parametric time-to-event model, while the conditional approach is implemented using both parametric and Cox proportional-hazard (CPH) models. The practical performances of these approaches are illustrated with an application to the tumor growth dynamics (TGD) modeling of data from two cancer clinical trials of nivolumab and docetaxel, where patients were followed until disease progression. The dataset consisted of 3504 tumor size measurements from 855 subjects, which were described by a TGD model. The dropout of subjects was described by a Weibull or CPH model. Simulated datasets were also used to further illustrate the properties of the VPC methods. The results showed that the more familiar full approach might not provide meaningful improvement for TGD model evaluation over the naive approach of not adjusting for dropout, and could be outperformed by the conditional approach using either the Weibull model or the Cox proportional hazard model. Overall, including confidence intervals in VPC should improve interpretation, the conditional approach was shown to be more generally applicable when dropout occurs, and the nonparametric approach could provide additional robustness.

Maribavir is approved for management of post-transplant cytomegalovirus (CMV) infections refractory and/or resistant to CMV therapies at a dose of 400 mg twice daily (BID). Population pharmacokinetic (PopPK) and exposure-response analyses were conducted to support the appropriateness of 400 mg BID dosing. A PopPK model was developed using non-linear mixed-effects modeling with pooled maribavir plasma concentration-time data from phase 1 and 2 studies (from 100 mg up to 1200 mg as single or repeated doses) and the phase 3 SOLSTICE study (400 mg BID). Exposure-response analyses were performed for efficacy, safety, and viral resistance based on data collected in the SOLSTICE study. Maribavir PK after oral administration was adequately described by a two-compartment model with first-order elimination, first-order absorption, and an absorption lag-time. There was no evidence that maribavir PK was affected by age, sex, race, diarrhea, vomiting, disease characteristics, or concomitant use of histamine H₂ blockers, or proton pump inhibitors. In the SOLSTICE study, higher maribavir exposure was not associated with increased probability of achieving CMV DNA viremia clearance, nor with reduced probability of treatment-emergent maribavir-resistant CMV mutations. A statistically significant association with maribavir exposure was identified for taste disturbance, fatigue, and treatment-emergent serious adverse events, while transplant type, enrollment region, CMV DNA level at baseline, and/or CMV resistance at baseline were identified as additional risk factors for these safety outcomes. In conclusion, the findings of these PopPK and exposure-response analyses provide further support for the recommended maribavir dose of 400 mg BID.

Clinical studies have found there still exists a lack of gene therapy dose-toxicity and dose-efficacy data that causes gene therapy dose selection to remain elusive. Model informed drug development (MIDD) has become a standard tool implemented throughout the discovery, development, and approval of pharmaceutical therapies, and has the potential to inform dose-toxicity and dose-efficacy relationships to support gene therapy dose selection. Despite this potential, MIDD approaches for gene therapy remain immature and require standardization to be useful for gene therapy clinical programs. With the goal to advance MIDD approaches for gene therapy, in this review we first provide an overview of gene therapy types and how they differ from a bioanalytical, formulation, route of administration, and regulatory standpoint. With this biological and regulatory background, we propose how MIDD can be advanced for AAV-based gene therapies by utilizing physiological based pharmacokinetic modeling and quantitative systems pharmacology to holistically inform AAV and target protein dynamics following dosing. We discuss how this proposed model, allowing for in-depth exploration of AAV pharmacology, could be the key the field needs to treat these unmet disease populations.

Protein therapeutics have revolutionized the treatment of a wide range of diseases. While they have distinct physicochemical characteristics that influence their absorption, distribution, metabolism, and excretion (ADME) properties, the relationship between the physicochemical properties and PK is still largely unknown. In this work we present a minimal physiologically-based pharmacokinetic (mPBPK) model that incorporates a multivariate quantitative relation between a therapeutic's physicochemical parameters and its corresponding ADME properties. The model's compound-specific input includes molecular weight, molecular size (Stoke's radius), molecular charge, binding affinity to FcRn, and specific antigen affinity. Through derived and fitted empirical relationships, the model demonstrates the effect of these compound-specific properties on antibody disposition in both plasma and peripheral tissues using observed PK data in mice and humans. The mPBPK model applies the two-pore hypothesis to predict size-based clearance and exposure of full-length antibodies (150 kDa) and antibody fragments (50-100 kDa) within a onefold error. We quantitatively relate antibody charge and PK parameters like uptake rate, non-specific binding affinity, and volume of distribution to capture the relatively faster clearance of positively charged mAb as compared to negatively charged mAb. The model predicts the terminal plasma clearance of slightly positively and negatively charged antibody in humans within a onefold error. The mPBPK model presented in this work can be used to predict the target-mediated disposition of a drug when compound-specific and target-specific properties are known. To our knowledge, a combined effect of antibody weight, size, charge, FcRn, and antigen has not been incorporated and studied in a single mPBPK model previously. By conclusively incorporating and relating a multitude of protein's physicochemical properties to observed PK, our mPBPK model aims to contribute as a platform approach in the early stages of drug development where many of these properties can be optimized to improve a molecule's PK and ultimately its efficacy.