Model-informed dose optimization of mycophenolic acid in pediatric kidney transplant patients.

PURPOSE: We aimed to develop and evaluate a population PK model of mycophenolic acid (MPA) in pediatric kidney transplant patients to aid MPA dose optimization. METHODS: Data were collected from pediatric kidney transplant recipients from a Dutch academic hospital (Radboudumc, the Netherlands). Pharmacokinetic model-building and model-validation analyses were performed using NONMEM. Subsequently, we externally evaluated the final model using data from another academic hospital. The final model was used to develop an optimized dosing regimen. RESULTS: Thirty pediatric patients were included of whom 266 measured MPA plasma concentrations, including 20 full pharmacokinetic (PK) curves and 24 limited sampling curves, were available. A two-compartment model with a transition compartment for Erlang-type absorption best described the data. The final population PK parameter estimates were Ktr (1.48 h-1; 95% CI, 1.15-1.84), CL/F (16.0 L h-1; 95% CI, 10.3-20.4), Vc/F (24.9 L; 95% CI, 93.0-6.71E25), Vp/F (1590 L; 95% CI, 651-2994), and Q/F (36.2 L h-1; 95% CI, 9.63-74.7). The performance of the PK model in the external population was adequate. An optimized initial dose scheme based on bodyweight was developed. With the licensed initial dose, 35% of patients were predicted to achieve the target AUC, compared to 42% using the optimized scheme. CONCLUSION: We have successfully developed a pharmacokinetic model for MPA in pediatric renal transplant patients. The optimized dosing regimen is expected to result in better target attainment early in treatment. It can be used in combination with model-informed follow-up dosing to further individualize the dose when PK samples become available.

A population pharmacokinetics model of balovaptan to support dose selection in adult and pediatric populations.

Balovaptan is a brain-penetrating vasopressin receptor 1a antagonist previously investigated for the core symptoms of autism spectrum disorder (ASD). A population pharmacokinetic (PK) model of balovaptan was developed, initially to assist clinical dosing for adult and pediatric ASD studies and subsequently for new clinical indications including malignant cerebral edema (MCE) and post-traumatic stress disorder. The final model incorporates one-compartment disposition and describes time- and dose-dependent non-linear PK through empirical drug binding and a gut extraction component with turnover. An age effect on clearance observed in children was modeled by an asymptotic function that predicts adult-equivalent exposures at 40% of the adult dose for children aged 2-4 years, 70% for 5-9 years, and at the full adult dose for ≥ 10 years. The model was adapted for intravenous (IV) balovaptan dosing and combined with in vitro and ex vivo pharmacodynamic data to simulate brain receptor occupancy as a guide for dosing in a phase II trial of MCE prophylaxis after acute ischemic stroke. A sequence of three stepped-dose daily infusions of 50, 25 and 15 mg over 30 or 60 min was predicted to achieve a target occupancy of ≥ 80% in ≥ 95% of patients over a 3-day period. This model predicts both oral and IV balovaptan exposure across a wide age range and will be a valuable tool to analyze and predict its PK in new indications and target populations, including pediatric patients.

Go beyond the limits of genetic algorithm in daily covariate selection practice.

Covariate identification is an important step in the development of a population pharmacokinetic/pharmacodynamic model. Among the different available approaches, the stepwise covariate model (SCM) is the most used. However, SCM is based on a local search strategy, in which the model-building process iteratively tests the addition or elimination of a single covariate at a time given all the others. This introduces a heuristic to limit the searching space and then the computational complexity, but, at the same time, can lead to a suboptimal solution. The application of genetic algorithms (GAs) for covariate selection has been proposed as a possible solution to overcome these limitations. However, their actual use during model building is limited by the extremely high computational costs and convergence issues, both related to the number of models being tested. In this paper, we proposed a new GA for covariate selection to address these challenges. The GA was first developed on a simulated case study where the heuristics introduced to overcome the limitations affecting currently available GA approaches resulted able to limit the selection of redundant covariates, increase replicability of results and reduce convergence times. Then, we tested the proposed GA on a real-world problem related to remifentanil. It obtained good results both in terms of selected covariates and fitness optimization, outperforming the SCM.

Population Pharmacokinetic and Concentration-QTc Analysis of Delamanid in Pediatric Participants with Multidrug-Resistant Tuberculosis.

A population pharmacokinetic analysis of delamanid and its major metabolite DM-6705 was conducted to characterize the pharmacokinetics of delamanid and DM-6705 in pediatric participants with multidrug-resistant tuberculosis (MDR-TB). Data from participants between the ages of 0.67 and 17 years, enrolled in 2 clinical trials, were utilized for the analysis. The final data set contained 634 delamanid and 706 DM-6705 valid plasma concentrations from 37 children. A transit model with three compartments best described the absorption of delamanid. Two-compartment models for each component with linear elimination were selected to characterize the dispositions of delamanid and DM-6705, respectively. The covariates included in the model were body weight on the apparent volume of distribution and apparent clearance (for both delamanid and DM-6705); formulation (dispersible versus film-coated tablet) on the mean absorption time; age, formulation, and dose on the bioavailability of delamanid; and age on the fraction of delamanid metabolized to DM-6705. Based on the simulations, doses for participants within different age/weight groups that result in delamanid exposure comparable to that in adults following the approved adult dose were calculated. By concentration-QTc (QTcB [QT corrected by Bazett's formula]) analysis, a significant positive correlation was detected with concentrations of DM-6705. However, the model-predicted upper bounds of the 90% confidence intervals of ΔQTc values were <10 ms at the simulated maximum concentration (Cmax) of DM-6705 following the administration of the maximum doses simulated. This suggests that the effect on the QT interval following the proposed dosing is unlikely to be clinically meaningful in children with MDR-TB who receive delamanid.

Population Pharmacokinetic-Pharmacodynamic Model of Oxfendazole in Healthy Adults in a Multiple Ascending Dose and Food Effect Study and Target Attainment Analysis.

Oxfendazole is a potent veterinary antiparasitic drug undergoing development for human use to treat multiple parasitic infections. Results from two recently completed phase I clinical trials conducted in healthy adults showed that the pharmacokinetics of oxfendazole is nonlinear, affected by food, and, after the administration of repeated doses, appeared to mildly affect hemoglobin concentrations. To facilitate oxfendazole dose optimization for its use in patient populations, the relationship among oxfendazole dose, pharmacokinetics, and hemoglobin concentration was quantitatively characterized using population pharmacokinetic-pharmacodynamic modeling. In fasting subjects, oxfendazole pharmacokinetics was well described by a one-compartment model with first-order absorption and elimination. The change in oxfendazole pharmacokinetics when administered following a fatty meal was captured by an absorption model with one transit compartment and increased bioavailability. The effect of oxfendazole exposure on hemoglobin concentration in healthy adults was characterized by a life span indirect response model in which oxfendazole has positive but minor inhibitory effect on red blood cell synthesis. Further simulation indicated that oxfendazole has a low risk of posing a safety concern regarding hemoglobin concentration, even at a high oxfendazole dose of 60 mg/kg of body weight once daily. The final model was further used to perform comprehensive target attainment simulations for whipworm infection and filariasis at various dose regimens and target attainment criteria. The results of our modeling work, when adopted appropriately, have the potential to greatly facilitate oxfendazole dose regimen optimization in patient populations with different types of parasitic infections.

Population pharmacokinetic model selection assisted by machine learning.

A fit-for-purpose structural and statistical model is the first major requirement in population pharmacometric model development. In this manuscript we discuss how this complex and computationally intensive task could benefit from supervised machine learning algorithms. We compared the classical pharmacometric approach with two machine learning methods, genetic algorithm and neural networks, in different scenarios based on simulated pharmacokinetic data. Genetic algorithm performance was assessed using a fitness function based on log-likelihood, whilst neural networks were trained using mean square error or binary cross-entropy loss. Machine learning provided a selection based only on statistical rules and achieved accurate selection. The minimization process of genetic algorithm was successful at allowing the algorithm to select plausible models. Neural network classification tasks achieved the most accurate results. Neural network regression tasks were less precise than neural network classification and genetic algorithm methods. The computational gain obtained by using machine learning was substantial, especially in the case of neural networks. We demonstrated that machine learning methods can greatly increase the efficiency of pharmacokinetic population model selection in case of large datasets or complex models requiring long run-times. Our results suggest that machine learning approaches can achieve a first fast selection of models which can be followed by more conventional pharmacometric approaches.

Population Pharmacokinetic Model of Oxfendazole and Metabolites in Healthy Adults following Single Ascending Doses.

Oxfendazole is a potent veterinary benzimidazole anthelmintic under transition to humans for the treatment of multiple parasitic infectious diseases. The first-in-human study evaluating the disposition of oxfendazole and its metabolites in healthy adults following single ascending oral doses from 0.5 to 60 mg/kg of body weight shows that oxfendazole pharmacokinetics is substantially nonlinear, which complicates correlating oxfendazole dose to exposure. To quantitatively capture the relation between oxfendazole dose and exposure, a population pharmacokinetic model for oxfendazole and its metabolites, oxfendazole sulfone and fenbendazole, in humans was developed using a nonlinear mixed-effect modeling approach. Our final model incorporated mechanistic characterization of dose-limited bioavailability as well as different oxfendazole metabolic processes and provided insight into the significance of presystemic metabolism in oxfendazole and metabolite disposition. Oxfendazole clinical pharmacokinetics was best described by a one-compartment model with nonlinear absorption and linear elimination. Oxfendazole apparent clearance and apparent volume of distribution were estimated to be 2.57 liters/h and 35.2 liters, respectively, at the lowest dose (0.5 mg/kg), indicating that oxfendazole is a low extraction drug with moderate distribution. The disposition of both metabolites was adequately characterized by a one-compartment model with formation rate-limited elimination. Fenbendazole formation from oxfendazole was primarily through systemic metabolism, while both presystemic and systemic metabolism were critical to the formation of oxfendazole sulfone. Our model adequately captured the concentration-time profiles of both oxfendazole and its two metabolites in healthy adults over a wide dose range. The model can be used to predict oxfendazole disposition under new dosing regimens to support dose optimization in humans.

PK/PD modeling of a clazosentan thorough QT study with hysteresis in concentration-QT and RR-QT.

Clazosentan's potential QT liability was investigated in a thorough QT study in which clazosentan was administered intravenously as a continuous infusion of 20 mg/h immediately followed by 60 mg/h. Clazosentan prolonged the placebo-corrected change-from-baseline QT interval corrected for RR with Fridericia's formula (ΔΔQTcF) with the maximum QT effect occurring 4 h after the maximum drug concentration, apparently associated with vomiting. The delayed effect precluded the standard linear modeling approach. This analysis aimed at characterizing the concentration-QT relationship in consideration of RR-QT hysteresis, concentration-ΔΔQTcF hysteresis, and the influence of vomiting. Nonlinear mixed-effects modeling was applied to characterize pharmacokinetics and pharmacodynamics, i.e., ΔΔQTcF. Simulations were used to predict ΔΔQTcF for expected therapeutic dose used in Phase 3 clinical development. Correction for RR-QT hysteresis did not influence ΔΔQTcF to a relevant extent. Pharmacokinetics of clazosentan were best described by a linear two-compartment model. The delayed QT prolongation was characterized by an indirect-response model with loglinear drug effect. Vomiting had no statistically significant influence on QT prolongation despite apparent differences between subjects vomiting and not vomiting, probably since vomiting occurred mostly after the main QT prolongation. Following a simulated 3-h infusion of 15 mg/h of clazosentan, the upper bound of the predicted 90% CI for mean ΔΔQTcF was expected to exceed the 10-ms regulatory threshold of concern with maximum effect 3.5 h after end of infusion. TRN: NCT03657446, 05 Sep 2018.

A Population Pharmacokinetic Analysis of Continuous Infusion of Cloxacillin during Staphylococcus aureus Bone and Joint Infections.

Intravenous administration of antibiotics is recommended during the early phase of methicillin-susceptible S. aureus (MSSA) bone and joint infection (BJI). We sought to compare the plasma concentrations of cloxacillin administered alternately by continuous and intermittent infusion (CI and ItI) in patients with MSSA BJI. In this prospective crossover trial, patients were randomly assigned to receive either 3 days of CI (two 75-mg/kg 12-h cloxacillin infusions per day) and then 3 days of ItI (four 37.5-mg/kg 1-h cloxacillin infusions per day) or vice versa. The drug concentration measurement was performed on day 3 of each type of administration at 1, 6, and 11 h and at 1, 2, 3, 4, and 6 h after the beginning of CI and ItI, respectively. We used the nonparametric algorithm NPAG to estimate population pharmacokinetic (PK) parameters. The final model was used to perform pharmacokinetic/pharmacodynamic (PK/PD) simulations and calculate the probabilities of target attainment (PTA) for several ItI and CI dosing regimens. We considered two PK/PD targets of time spent above the MIC for free cloxacillin concentrations (fT>MIC): 50 and 100%. Eighty-four concentrations from 11 patients were analyzed. A two-compartment model adequately described the data. ItI with q6h regimens and short 1-h infusions of 2,000 or 3,000 mg were associated with low PTA, even for the low target (50% fT>MIC) while 3-h infusions and continuous infusions (6 to 12 g/day) were associated with a PTA of >90% for an MIC up to 0.5 mg/liter. These results support the use of prolonged or continuous infusion of cloxacillin in patients with BJI.

Novel Approach for the Bioequivalence Assessment of Topical Cream Formulations: Model-Based Analysis of Tape Stripping Data Correctly Concludes BE and BIE.

PURPOSE: The purpose of this study was (a) to suggest a novel dermatopharmacokinetic (DPK) approach from which pharmacokinetic parameters relevant to the bioequivalence (BE) assessment of a topical formulation can be deduced while circumventing the need for numerous measurements and assumptions, and (b) to investigate whether this approach enables the correct conclusion of BE and bioinequivalence (BIE). METHODS: Bioequivalent and bioinequivalent formulations of acyclovir were compared versus a reference product (Zovirax®). Tape Stripping was conducted at only one dose duration during the uptake phase to generate drug content in stratum corneum versus time profiles, each time point corresponding to one stripped layer. Nonlinear mixed effect modeling (ADAPT5®) (MLEM algorithm) was used to fit the DPK data and to estimate the rate (Kin) and extent (FS) of drug absorption/input into the skin. Results were evaluated using the average BE approach. RESULTS: Estimated exposure metrics were within the usual BE limits for the bioequivalent formulation (FS: 102.4 [90%CI: 97.5-107.7]; Kin: 94.2 [90%CI: 83.7-106.0]), but outside those limits for the bioinequivalent formulation (FS: 43.4 [90%CI: 27.9-67.6]; Kin: 54.5 [90%CI: 36.6-81.1]). CONCLUSIONS: The proposed novel DPK approach was shown to be successful, robust and applicable to assess BE and BIE correctly between topical formulations.