Journal of Pharmacokinetics and Pharmacodynamics最新文献

The development of optimized dosing regimens plays a crucial role in oncology drug development. This study focused on the population pharmacokinetic modelling and simulation of docetaxel, comparing the pharmacokinetic exposure of oral docetaxel plus encequidar (oDox + E) with the standard of care intravenous (IV) docetaxel regimen. The aim was to evaluate the feasibility of oDox + E as a potential alternative to IV docetaxel. The article demonstrates an approach which aligns with the FDA's Project Optimus which aims to improve oncology drug development through model informed drug development (MIDD). The key question answered by this study was whether a feasible regimen of oDox + E existed. The purpose of this question was to provide an early GO / NO-GO decision point to guide drug development and improve development efficiency.

Methods:  A stepwise approach was employed to develop a population pharmacokinetic model for total and unbound docetaxel plasma concentrations after IV docetaxel and oDox + E administration. Simulations were performed from the final model to assess the probability of target attainment (PTA) for different oDox + E dose regimens (including multiple dose regimens) in relation to IV docetaxel using AUC over effective concentration (AUCOEC) metric across a range of effective concentrations (EC). A Go / No-Go framework was defined-the first part of the framework assessed whether a feasible oDox + E regimen existed (i.e., a PTA ≥ 80%), and the second part defined the conditions to proceed with a Go decision.

Results:  The overall population pharmacokinetic model consisted of a 3-compartment model with linear elimination, constant bioavailability, constant binding mechanics, and a combined error model. Simulations revealed that single dose oDox + E regimens did not achieve a PTA greater than 80%. However, two- and three-dose regimens at 600 mg achieved PTAs exceeding 80% for certain EC levels.

Conclusion:  The study demonstrates the benefits of MIDD using oDox + E as a motivating example. A population pharmacokinetic model was developed for the total and unbound concentration in plasma of docetaxel after administration of IV docetaxel and oDox + E. The model was used to simulate oDox + E dose regimens which were compared to the current standard of care IV docetaxel regimen. A GO / NO-GO framework was applied to determine whether oDox + E should progress to the next phase of drug development and whether any conditions should apply. A two or three-dose regimen of oDox + E at 600 mg was able to achieve non-inferior pharmacokinetic exposure to current standard of care IV docetaxel in simulations. A Conditional GO decision was made based on this result and further quantification of the "effective concentration" would improve the ability to optimise the dose regimen.

This is a correspondence on "Evaluation of ChatGPT and Gemini large language models for pharmacometrics with NONMEM". Additional concern on using ChatGPT and Gemini is provided.

Dose personalization improves patient outcomes for many drugs with a narrow therapeutic index and high inter-individuality variability, including busulfan. Non-compartmental analysis (NCA) and model-based methods like maximum a posteriori Bayesian (MAP) approaches are two methods routinely used for dose optimization. These approaches vary in how they estimate patient-specific pharmacokinetic parameters to inform a dose and the impact of these differences is not well-understood. Using busulfan as an example application and area under the concentration-time curve (AUC) as a target exposure metric, these estimation methods were compared using retrospective patient data (N = 246) and simulated precision dosing treatment courses. NCA was performed with or without peak extension, and MAP Bayesian estimation was performed using either the one-compartment Shukla model or the two-compartment McCune model. All methods showed good agreement on real-world data (correlation coefficients of 0.945-0.998) as assessed by Bland-Altman plots, although agreement between NCA and MAP methods was higher during the first dosing interval (0.982-0.994) compared to subsequent dosing intervals (0.918-0.938). In dose adjustment simulations, both NCA and MAP estimated high target attainment (> 98%) although true simulated target attainment was lower for NCA (63-66%) versus MAP (91-93%). The largest differences in AUC estimation were due to different assumptions for the shape of the concentration curve during the infusion phase, followed by how the methods considered time-dependent clearance and concentration-time points collected in earlier intervals. In conclusion, although AUC estimates between the two methods showed good correlation, in a simulated study, MAP lead to higher target attainment. When changing from one method to another, or changing infusion duration and other factors, optimum estimated exposure targets may require adjusting to maintain a consistent exposure.

Aprocitentan is a novel, potent, dual endothelin receptor antagonist that recently demonstrated efficacy in the treatment of difficult-to-treat (resistant) hypertension. The aim of this study was to develop a population pharmacokinetic (PK) model describing aprocitentan plasma concentration over time, to investigate relationships between subject-specific factors (covariates) and model parameters, and to quantify the influence of the identified covariates on the exposure to aprocitentan via model-based simulations, enabling judgment about the clinical relevance of the covariates.PK data from 902 subjects in ten Phase 1, one Phase 2, and one Phase 3 study were pooled to develop a joint population PK model. The concentration-time course of aprocitentan was described by a two-compartment model with absorption lag time, first-order absorption and elimination, and reduced relative bioavailability following very high doses of 300 and 600 mg.The population PK model described the observed data well. Volume and clearance parameters were associated with body weight. Renal function as reflected by estimated glomerular filtration rate (eGFR), hepatic impairment, and sex were identified as relevant covariates on clearance.The subject-specific characteristics of body weight, eGFR, hepatic impairment, and sex were shown to influence exposure parameters area under the concentration-time curve and maximum concentration in steady state to a limited extent, i.e., not more than 25% different from a reference subject, and therefore do not warrant dose adjustments.

The use of β-lactam (BL) and β-lactamase inhibitor (BLI) combinations, such as piperacillin-tazobactam (PIP-TAZ) is an effective strategy to combat infections by extended-spectrum β-lactamase-producing bacteria. However, in Gram-negative bacteria, resistance (both mutational and adaptive) to BL-BLI combination can still develop through multiple mechanisms. These mechanisms may include increased β-lactamase activity, reduced drug influx, and increased drug efflux. Understanding the relative contribution of these mechanisms during resistance development helps identify the most impactful mechanism to target in designing a treatment to counter BL-BLI resistance. This study used semi-mechanistic mathematical modeling in combination with antibiotic sensitivity assays to assess the potential impact of different resistance mechanisms during the development of PIP-TAZ resistance in a Klebsiella pneumoniae isolate expressing CTX-M-15 and SHV-1 β-lactamases. The mathematical models were used to evaluate the potential impact of several cellular changes as a sole mediator of PIP-TAZ resistance. Our semi-mechanistic model identified 2 out of the 13 inspected mechanisms as key resistance mechanisms that may independently support the observed magnitude of PIP-TAZ resistance, namely porin loss and efflux pump up-regulation. Simulation using the resulting models also suggested the possible adjustment of PIP-TAZ dose outside its commonly used 8:1 dosing ratio. The current study demonstrated how theory-based mechanistic models informed by experimental data can be used to support hypothesis generation regarding potential resistance mechanisms, which may guide subsequent experimental studies.

Currently, model-informed precision dosing uses one population pharmacokinetic model that best fits the target population. We aimed to develop a subgroup identification-based model selection approach to improve the predictive performance of individualized dosing, using vancomycin in neonates/infants as a test case. Data from neonates/infants with at least one vancomycin concentration was randomly divided into training and test dataset. Population predictions from published vancomycin population pharmacokinetic models were calculated. The single best-performing model based on various performance metrics, including median absolute percentage error (APE) and percentage of predictions within 20% (P20) or 60% (P60) of measurement, were determined. Clustering based on median APEs or clinical and demographic characteristics and model selection by genetic algorithm was used to group neonates/infants according to their best-performing model. Subsequently, classification trees to predict the best-performing model using clinical and demographic characteristics were developed. A total of 208 vancomycin treatment episodes in training and 88 in test dataset was included. Of 30 identified models from the literature, the single best-performing model for training dataset had P20 26.2-42.6% in test dataset. The best-performing clustering approach based on median APEs or clinical and demographic characteristics and model selection by genetic algorithm had P20 44.1-45.5% in test dataset, whereas P60 was comparable. Our proof-of-concept study shows that the prediction of the best-performing model for each patient according to the proposed model selection approaches has the potential to improve the predictive performance of model-informed precision dosing compared with the single best-performing model approach.

Medication nonadherence is one of the largest problems in healthcare today, particularly for patients undergoing long-term pharmacotherapy. To combat nonadherence, it is often recommended to prescribe so-called "forgiving" drugs, which maintain their effect despite lapses in patient adherence. Nevertheless, drug forgiveness is difficult to quantify and compare between different drugs. In this paper, we construct and analyze a stochastic pharmacokinetic/pharmacodynamic (PK/PD) model to quantify and understand drug forgiveness. The model parameterizes a medication merely by an effective rate of onset of effect when the medication is taken (on-rate) and an effective rate of loss of effect when a dose is missed (off-rate). Patient dosing is modeled by a stochastic process that allows for correlations in missed doses. We analyze this "on/off" model and derive explicit formulas that show how treatment efficacy depends on drug parameters and patient adherence. As a case study, we compare the effects of nonadherence on the efficacy of various antihypertensive medications. Our analysis shows how different drugs can have identical efficacies under perfect adherence, but vastly different efficacies for adherence patterns typical of actual patients. We further demonstrate that complex PK/PD models can indeed be parameterized in terms of effective on-rates and off-rates. Finally, we have created an online app to allow pharmacometricians to explore the implications of our model and analysis.

Two-pore physiologically based pharmacokinetic (PBPK) modeling has demonstrated its potential in describing the pharmacokinetics (PK) of different-size proteins. However, all existing two-pore models lack either diverse proteins for validation or interspecies extrapolation. To fill the gap, here we have developed and optimized a translational two-pore PBPK model that can characterize plasma and tissue disposition of different-size proteins in mice, rats, monkeys, and humans. Datasets used for model development include more than 15 types of proteins: IgG (150 kDa), F(ab)2 (100 kDa), minibody (80 kDa), Fc-containing proteins (205, 200, 110, 105, 92, 84, 81, 65, or 60 kDa), albumin conjugate (85.7 kDa), albumin (67 kDa), Fab (50 kDa), diabody (50 kDa), scFv (27 kDa), dAb2 (23.5 kDa), proteins with an albumin-binding domain (26, 23.5, 22, 16, 14, or 13 kDa), nanobody (13 kDa), and other proteins (110, 65, or 60 kDa). The PBPK model incorporates: (i) molecular weight (MW)-dependent extravasation through large and small pores via diffusion and filtration, (ii) MW-dependent renal filtration, (iii) endosomal FcRn-mediated protection from catabolism for IgG and albumin-related modalities, and (iv) competition for FcRn binding from endogenous IgG and albumin. The finalized model can well characterize PK of most of these proteins, with area under the curve predicted within two-fold error. The model also provides insights into contribution of renal filtration and lysosomal degradation towards total elimination of proteins, and contribution of paracellular convection/diffusion and transcytosis towards extravasation. The PBPK model presented here represents a cross-modality, cross-species platform that can be used for development of novel biologics.

The paper extended the TMDD model to drugs with more than two (N > 2) identical binding sites (N-to-one TMDD). The quasi-steady-state (N-to-one QSS), quasi-equilibrium (N-to-one QE), irreversible binding (N-to-one IB), and Michaelis–Menten (N-to-one MM) approximations of the model were derived. To illustrate properties of new equations and approximations, N = 4 case was investigated numerically. Using simulations, the N-to-one QSS approximation was compared with the full N-to-one TMDD model. As expected, and similarly to the standard TMDD for monoclonal antibodies (mAb), N-to-one QSS predictions were nearly identical to N-to-one TMDD predictions, except for times of fast changes following initiation of dosing, when equilibrium has not yet been reached. Predictions for mAbs with soluble targets (slow elimination of the complex) were simulated from the full 4-to-one TMDD model and were fitted to the 4-to-one TMDD model and to its QSS approximation. It was demonstrated that the 4-to-one QSS model provided nearly identical description of not only the observed (simulated) total drug and total target concentrations, but also unobserved concentrations of the free drug, free target, and drug-target complexes. For mAb with a membrane-bound target, the 4-to-one MM approximation adequately described the data. The 4-to-one QSS approximation converged 8 times faster than the full 4-to-one TMDD.

Polymyxin B (PMB) is considered a last-line treatment for multidrug-resistant (MDR) gram-negative bacterial infections. Model-informed precision dosing with population pharmacokinetics (PopPK) models could help to individualize PMB dosing regimens and improve therapy. However, the external prediction ability of the established PopPK models has not been fully elaborated. This study aimed to systemically evaluate eleven PMB PopPK models from ten published literature based on a new independent population, which was divided into four different populations, patients with liver dysfunction, kidney dysfunction, liver and kidney dysfunction, and normal liver and kidney function. The whole data set consisted of 146 patients with 391 PMB concentrations. The prediction- and simulation-based diagnostics and Bayesian forecasting were conducted to evaluate model predictability. In the overall evaluation process, none of the models exhibited satisfactory predictive ability in both prediction- and simulation-based diagnostic simultaneously. However, the evaluation of the models in the subgroup of patients with normal liver and kidney function revealed improved predictive performance compared to those with liver and/or kidney dysfunction. Bayesian forecasting demonstrated enhanced predictability with the incorporation of two to three prior observations. The external evaluation highlighted a lack of consistency between the prediction results of published models and the external validation dataset. Nonetheless, Bayesian forecasting holds promise in improving the predictive performance of the models, and feedback from therapeutic drug monitoring is crucial in optimizing individual dosing regimens.