A diagnostic tool for population models using non-compartmental analysis: The ncappc package for R

Physiologically based pharmacokinetic model of brivaracetam to predict the exposure and dose exploration in hepatic impairment and elderly populations

Brivaracetam (BRV) is a new third-generation antiseizure medication for the treatment of focal epileptic seizures. Its use has been increasing among epileptic populations in recent years, but pharmacokinetic (PK) behavior may change in hepatic impairment and the elderly populations. Due to ethical constraints, clinical trials are difficult to conduct and data are limited. This study used PK-Sim® to develop a physiologically based pharmacokinetic (PBPK) model for adults and extrapolate it to hepatic impairment and the elderly populations. The model was evaluated with clinical PK data, and dosage explorations were conducted. For the adult population with mild hepatic impairment, the dose is recommended to be adjusted to 70% of the recommended dose, and to 60% for moderate and severe hepatic impairment. For the elderly population with mild hepatic impairment under 80 years old, it is recommended that the dose be adjusted to 60% of the recommended dose and to 50% for moderate and severe conditions. The elderly population with hepatic impairment over 80 years old is adjusted to 50% of the recommended dose for all stages. Healthy elderly do not need to adjust. The BRV PBPK model was successfully developed, studying exposure in hepatic impairment and elderly populations and optimizing dosing regimens.

Physiologically based pharmacokinetic models for predicting lamotrigine exposure and dose optimization in pediatric patients receiving combination therapy with carbamazepine or valproic acid

INTRODUCTION

Lamotrigine (LTG) is an antiepileptic drug that has been used in pediatric epilepsy as a combination therapy or monotherapy after stabilization in recent years. However, there are significant drug-drug interactions (DDI) between LTG and combined drugs such as carbamazepine (CBZ) and valproic acid (VPA). It is particularly important to consider the risk of DDI in combination therapy for intractable epilepsy in pediatric patients. Therefore, it is necessary to adjust the dosage of LTG accordingly. The aim of this study was to establish and validate a pediatric physiologically based pharmacokinetic (PBPK) model for predicting LTG exposure. The model is designed to explore the potential for quantifying pharmacokinetic (PK) DDI of LTG when administered concurrently with CBZ or VPA in pediatric patients.

METHOD

Adult and pediatric PBPK models for LTG and VPA were developed using PK-Sim® software in combination with physiological information and drug-specific parameters, and a DDI model was developed in combination with the published CBZ model. The models were validated against available PK data.

RESULTS

Predictive and observational results in adults, children, and the DDI model were in good agreement. The recommended doses of LTG for preschool children (2-6 years) and school-aged children (6-12 years) in the absence of drug interactions were 1.47 and 1.2 times higher than those for adults, respectively; 3.1 and 2.6 times higher than those for adults in combination with CBZ; and 0.67 and 0.57 times lower than those for adults in combination with VPA. In addition, plasma exposures in adolescents (12-18 years) were similar to those in adults at the same doses.

CONCLUSION

We have successfully developed PBPK models and DDI models for LTG in adults and children, which provide a reference for rational drug use in the pediatric population.

Evaluation of model-integrated evidence approaches for pharmacokinetic bioequivalence studies using model averaging methods

Conventional approaches for establishing bioequivalence (BE) between test and reference formulations using non-compartmental analysis (NCA) may demonstrate low power in pharmacokinetic (PK) studies with sparse sampling. In this case, model-integrated evidence (MIE) approaches for BE assessment have been shown to increase power, but may suffer from selection bias problems if models are built on the same data used for BE assessment. This work presents model averaging methods for BE evaluation and compares the power and type I error of these methods to conventional BE approaches for simulated studies of oral and ophthalmic formulations. Two model averaging methods were examined: bootstrap model selection and weight-based model averaging with parameter uncertainty from three different sources, either from a sandwich covariance matrix, a bootstrap, or from sampling importance resampling (SIR). The proposed approaches increased power compared with conventional NCA-based BE approaches, especially for the ophthalmic formulation scenarios, and were simultaneously able to adequately control type I error. In the rich sampling scenario considered for oral formulation, the weight-based model averaging method with SIR uncertainty provided controlled type I error, that was closest to the target of 5%. In sparse-sampling designs, especially the single sample ophthalmic scenarios, the type I error was best controlled by the bootstrap model selection method.

Development and comparison of model-integrated evidence approaches for bioequivalence studies with pharmacokinetic end points

By applying nonlinear mixed-effect (NLME) models, model-integrated evidence (MIE) approaches are able to analyze bioequivalence (BE) data with pharmacokinetic end points that have sparse sampling, which is problematic for non-compartmental analysis (NCA). However, MIE approaches may suffer from inflation of type I error due to underestimation of parameter uncertainty and to the assumption of asymptotic normality. In this study, we developed a MIE BE analysis method that is based on a pre-defined model and consists of several steps including model fitting, uncertainty assessment, simulation, and BE determination. The presented MIE approach has several improvements compared with the previously reported model-integrated methods: (1) treatment, sequence, and period effects are only added to absorption parameters (such as relative bioavailability and rate of absorption) instead of all PK parameters; (2) a simulation step is performed to generate confidence intervals of the pharmacokinetic metrics for BE assessment; and (3) in an effort to maintain type I error, two more advanced parameter uncertainty evaluation approaches are explored, a nonparametric (case resampling) bootstrap, and sampling importance resampling (SIR). To evaluate the developed method and compare the uncertainty assessment methods, simulation experiments were performed for BE studies using a two-way crossover design with different amounts of information (sparse to rich designs) and levels of variability. Based on the simulation results, the method using SIR for parameter uncertainty quantification controls type I error at the nominal level of 0.05 (i.e., the significance level set for BE evaluation) even for studies with small sample size and/or sparse sampling. As expected, our MIE approach for BE assessment exhibited higher power than the NCA-based method, especially as the data becomes sparser and/or more variable.

Indirect modeling of derived outcomes: Are minor prediction discrepancies a cause for concern?

It is often a goal of model development to predict data from which a variety of outcomes can be derived, such as threshold-based categorization or change from baseline (CFB) transformations. This approach can improve power or support multiple decisions. Because these derivations are indirectly predicted from the model, they are valuable tests for misspecification when used in visual or numeric predictive checks (V/NPCs). However, derived outcome V/NPCs (especially if primary or key secondary) are often overly scrutinized and held to an uncommon standard when comparing model predictions to point estimates, even if by conventional standards both the directly and indirectly modeled data are captured well. Here, simulations of directly modeled data were used to determine where apparent issues in V/NPCs of derived outcomes are expected. Two types of datasets were simulated: (1) a simple pre-post study and (2) pharmacokinetic/pharmacodynamic data from a dose-ranging study. A psoriasis exposure-response model case study was also assessed. V/NPCs were generated on the raw data, CFB data, and placebo-corrected CFB (dCFB) data, and binned summary statistics of the observed data for each trial were graded as being strongly or weakly supportive of a predictive model (within the interquartile range or the 95% central distribution of all simulated trials, respectively). In all cases, the strength of support in direct data V/NPCs was minimally related to that in derived outcome V/NPCs. There are myriad benefits to modeling the underlying data of a derived measure, and these results support caution in discarding adequate models based on overly strict derived measure predictive checks.

Evaluation of critical parameters in the hollow-fibre system for tuberculosis: A case study of moxifloxacin

Aims

The hollow‐fibre system for tuberculosis (HFS‐TB) is a preclinical model qualified by the European Medicines Agency to underpin the anti‐TB drug development process. It can mimic in vivo pharmacokinetic (PK)–pharmacodynamic (PD) attributes of selected antimicrobials, which could feed into in silico models to inform the design of clinical trials. However, historical data and published protocols are insufficient and omit key information to allow experiments to be reproducible. Therefore, in this work, we aim to optimize and standardize various HFS‐TB operational procedures.

Methods

First, we characterized bacterial growth dynamics with different types of hollow‐fibre cartridges, Mycobacterium tuberculosis strains and media. Second, we mimicked a moxifloxacin PK profile within hollow‐fibre cartridges, in order to check drug–fibres compatibility. Lastly, we mimicked the moxifloxacin total plasma PK profile in human after once daily oral dose of 400 mg to assess PK–PD after different sampling methods, strains, cartridge size and bacterial adaptation periods before drug infusion into the system.

Results

We found that final bacterial load inside the HFS‐TB was contingent on the studied variables. Besides, we demonstrated that drug–fibres compatibility tests are critical preliminary HFS‐TB assays, which need to be properly reported. Lastly, we uncovered that the sampling method and bacterial adaptation period before drug infusion significantly impact actual experimental conclusions.

Conclusion

Our data contribute to the necessary standardization of HFS‐TB experiments, draw attention to multiple aspects of this preclinical model that should be considered when reporting novel results and warn about critical parameters in the HFS‐TB currently overlooked.

Safety, tolerability and toxicokinetics of the novel mitochondrial drug SUL-138 administered orally to rat and minipig

Noncommunicable Chronic Diseases (NCD) are a socioeconomic burden and considered one of the major health challenges for coming decades. Mitochondrial dysfunction has been implicated mechanistically in their pathophysiology. Therefore, targeting mitochondria holds great promise to improve clinical outcomes in NCDs. SUL-138, an orally bioavailable small molecule efficacious from 0.5 mg/kg, improves mitochondrial function during disease in several preclinical animal models. As preparation for a First-in-Human (FIH) trial, SUL-138 was investigated in 30-day GLP repeated dose toxicity studies in rat and minipig, selected based on their comparability with human metabolism, to determine toxicokinetics, potential toxicity and its reversibility. Rats were allocated to either vehicle, 27, 136 or 682 mg/kg SUL-138 dose groups and minipigs were allocated to either vehicle, 16, 82 or 409 mg/kg. Treatment occurred orally for 30 days followed by a recovery period of 14 days. During these studies clinical observations, toxicokinetic, clinical pathology, necropsy and histopathology evaluations were performed. There was significant systemic exposure to SUL-138 and toxicokinetics was characterized by a rapid absorption and elimination. In the rat, toxicokinetics was dose-proportional and AUC0-tlast ratios in both species indicated that SUL-138 does not accumulate in vivo. No treatment-related adverse effects were observed for dose levels up to 136 and 82 mg/kg/day in rat and minipig respectively. In conclusion, these preclinical studies demonstrate that SUL-138 is well tolerated after repeated administration in rat and minipig, with NOAELs of 136 and 82 mg/kg/day, respectively.

Pharmacokinetics and Safety of Twice-daily Ritonavir-boosted Atazanavir With Rifampicin

BACKGROUND

Critical drug-drug interactions (DDI) and hepatotoxicity complicate concurrent use of rifampicin and protease inhibitors. We investigated whether dose escalation of atazanavir/ritonavir could safely overcome the DDI with rifampicin.

METHODS

DERIVE (NCT04121195, EDCTP) was a dose-escalation trial in people with human immunodeficiency virus (HIV) on atazanavir/ritonavir-based antiretroviral therapy (ART) in Uganda. Four intensive pharmacokinetic (PK) visits were performed: PK1 300/100 mg OD (baseline); PK2 300/100 mg OD with rifampicin 600 mg; PK3 300/100 mg twice a day (BID) with rifampicin 600 mg OD; PK4 300/100 mg BID with rifampicin 1200 mg OD. Dolutegravir 50 mg BID throughout the study period ensured participants remained protected from subtherapeutic atazanavir concentrations. The data were interpreted with noncompartmental analysis. The target minimum concentration was atazanavir's protein-adjusted IC90 (PA-IC90), 0.014 mg/L.

RESULTS

We enrolled 26 participants (23 female) with median (range) age 44 (28-61) years and weight 67 (50-75) kg. Compared with PK1, atazanavir Ctau, and AUC were significantly reduced at PK2 by 96% and 85%, respectively. The escalation to BID dosing (PK3) reduced this difference in Ctau, and AUC24 to 18% lower and 8% higher, respectively. Comparable exposures were maintained with double doses of rifampicin. Lowest Ctau during PK1, PK3, and PK4 were 12.7-, 4.8-, and 8.6-fold higher than PA-IC90, respectively, whereas 65% of PK2 Ctau were below the limit of quantification (0.03 mg/L), hence likely below PA-IC90. No participant developed significant elevation of liver enzymes, reported a serious adverse event (SAE) or experienced rebound viraemia.

CONCLUSIONS

Twice daily atazanavir/ritonavir during rifampicin co-administration was well tolerated and achieved plasma concentrations above the target.

CLINICAL TRIALS REGISTRATION

NCT04121195. Registered on 09 October 2019, https://clinicaltrials.gov/ct2/show/NCT04121195.

A physiologically-based pharmacokinetic/pharmacodynamic modeling approach for drug-drug-gene interaction evaluation of S-warfarin with fluconazole

Warfarin is a widely used anticoagulant, and its S-enantiomer has higher potency compared to the R-enantiomer. S-warfarin is mainly metabolized by cytochrome P450 (CYP) 2C9, and its pharmacological target is vitamin K epoxide reductase complex subunit 1 (VKORC1). Both CYP2C9 and VKORC1 have genetic polymorphisms, leading to large variations in the pharmacokinetics (PKs) and pharmacodynamics (PDs) of warfarin in the population. This makes dosage management of warfarin difficult, especially in the case of drug-drug interactions (DDIs). This study provides a whole-body physiologically-based pharmacokinetic/PD (PBPK/PD) model of S-warfarin for predicting the effects of drug-drug-gene interactions on S-warfarin PKs and PDs. The PBPK/PD model of S-warfarin was developed in PK-Sim and MoBi. Drug-dependent parameters were obtained from the literature or optimized. Of the 34 S-warfarin plasma concentration-time profiles used, 96% predicted plasma concentrations within twofold range compared to observed data. For S-warfarin plasma concentration-time profiles with CYP2C9 genotype, 364 of 386 predicted plasma concentration values (~94%) fell within the twofold of the observed values. This model was tested in DDI predictions with fluconazole as CYP2C9 perpetrators, with all predicted DDI area under the plasma concentration-time curve to the last measurable timepoint (AUClast) ratio within twofold of the observed values. The anticoagulant effect of S-warfarin was described using an indirect response model, with all predicted international normalized ratio (INR) within twofold of the observed values. This model also incorporates a dose-adjustment method that can be used for dose adjustment and predict INR when warfarin is used in combination with CYP2C9 perpetrators.

Effect of albumin substitution on pharmacokinetics of piperacillin/tazobactam in patients with severe burn injury admitted to the ICU

BACKGROUND

Pathophysiological changes in severely burned patients alter the pharmacokinetics (PK) of anti-infective agents, potentially leading to subtherapeutic concentrations at the target site. Albumin supplementation, to support fluid resuscitation, may affect pharmacokinetic properties by binding drugs. This study aimed to investigate the PK of piperacillin/tazobactam in burn patients admitted to the ICU before and after albumin substitution as total and unbound concentrations in plasma.

PATIENTS AND METHODS

Patients admitted to the ICU and scheduled for 4.5 g piperacillin/tazobactam administration and 200 mL of 20% albumin substitution as part of clinical routine were included. Patients underwent IV microdialysis, and simultaneous arterial plasma sampling, at baseline and multiple timepoints after drug administration. PK analysis of total and unbound drug concentrations under steady-state conditions was performed before and after albumin supplementation.

RESULTS

A total of seven patients with second- to third-degree burns involving 20%-60% of the total body surface were enrolled. Mean (SD) AUC0-8 (h·mg/L) of total piperacillin/tazobactam before and after albumin substitution were 402.1 (242)/53.2 (27) and 521.8 (363)/59.7 (32), respectively. Unbound mean AUC0-8 before and after albumin supplementation were 398.9 (204)/54.5 (25) and 456.4 (439)/64.5 (82), respectively.

CONCLUSIONS

Albumin supplementation had little impact on the PK of piperacillin/tazobactam. After albumin supplementation, there was a numerical increase in mean AUC0-8 of total and unbound piperacillin/tazobactam, whereas similar Cmax values were observed. Future studies may investigate the effect of albumin supplementation on drugs with a higher plasma protein binding.

Pharmacokinetics and bioequivalence of Withania somnifera (Ashwagandha) extracts - A double blind, crossover study in healthy adults

Introduction

Withania somnifera (WS) or ashwagandha is an adaptogenic plant used extensively in traditional medicines and as a food supplement. Despite a long history of use and numerous clinical trials, the human pharmacokinetics of withanolides, the active phytochemicals in WS extracts, have not been fully evaluated. This study evaluated the oral pharmacokinetics and bioequivalence of active withanolides in human plasma after administration of a single dose of two commercial ashwagandha extracts containing equal amounts of total withanolides.

Methods

This randomized, double-blind, single-dose crossover study of 16 healthy human volunteers evaluated the acute oral bioavailability of withanolides and the bioequivalence of two WS extracts, WS-35 and WS-2.5. WS-35 was standardized to total withanolides not less than 40% comprising not less than 35% withanolide glycosides and WS-2.5 was standardized to 2.5% withanolides. The clinical dosages were normalized to 185 mg of total withanolide in each extract at the bioequivalent dosages. The pharmacokinetic parameters of withanolide A, withanoside IV, withaferin A, and total withanolides were quantified in the blood plasma using a validated LC-MS/MS method.

Results

The half-life, C-max, and mean residence time of the total withanolides were 5.18, 5.62 and 4.13 times significantly higher and had lower systemic clearance with WS-35 than with WS-2.5 extract. Considering the plasma AUC 0-inf of total withanolides per mg of each WS extract administered orally, WS-35 was 280.74 times more bioavailable than WS-2.5.

Conclusion

The results of this study highlight the importance of withanolide glycosides in improving the pharmacokinetics of WS extracts. Owing to its superior pharmacokinetic profile, WS-35, with 35% withanolide glycosides, is a promising candidate for further studies on Withania somnifera.

Clinical trial registration

CTRI/2020/10/028397 [registered on:13/10/2020] (Trial prospectively registered) http://ctri.nic.in/Clinicaltrials/pmaindet2.php?trialid=42149&EncHid=&userName=CTRI/2020/10/028397.

Revisiting the in-vitro and in-vivo considerations for in-silico modelling of complex injectable drug products

Complex injectable drug products (CIDPs) have often been developed to modulate the pharmacokinetics along with efficacy for therapeutic agents used for remediation of chronic disorders. The effective development of CIDPs has exhibited complex kinetics associated with multiphasic drug release from the prepared formulations. Consequently, predictability of pharmacokinetic modelling for such CIDPs has been difficult and there is need for advanced complex computational models for the establishment of accurate prediction models for in-vitro-in-vivo correlation (IVIVC). The computational modelling aims at supplementing the existing knowledge with mathematical equations to develop formulation strategies for generation of predictable and discriminatory IVIVC. Such an approach would help in reduction of the burden of effect of hidden factors on preclinical to clinical translations. Computational tools like physiologically based pharmacokinetics (PBPK) modelling have combined physicochemical and physiological properties along with IVIVC characteristics of clinically used formulations. Such techniques have helped in prediction and understanding of variability in pharmacodynamic parameters of potential generic products to clinically used formulations like Doxil®, Ambisome®, Abraxane® in healthy and diseased population using mathematical equations. The current review highlights the important formulation characteristics, in-vitro, preclinical in-vivo aspects which need to be considered while developing a stimulatory predictive PBPK model in establishment of an IVIVC and in-vitro-in-vivo relationship (IVIVR).

Relative bioavailability of delamanid 50 mg tablets dispersed in water in healthy adult volunteers

AIM

Delamanid is a novel drug for the treatment of drug-resistant tuberculosis, manufactured as 50-mg solid and 25-mg dispersible tablets. We evaluated the effects of dispersing the 50-mg tablet, focusing on the relative bioavailability.

METHODS

Delamanid, 50-mg tablets administered dispersed vs swallowed whole, was investigated in a phase I, four-period, crossover study. Two of three dose strengths of delamanid (25, 50 or 100 mg) were given to healthy adult participants, in both whole and dispersed forms, with a 7-day washout period. Blood samples were collected over 168 h after each dose. Delamanid and its metabolite DM-6705 were analysed with a validated liquid chromatography tandem mass spectrometry assay. The pharmacokinetics of both analytes were analysed using nonlinear mixed-effect modelling. Palatability and acceptability were determined using a standardized questionnaire.

RESULTS

Twenty-four participants completed the study. The bioavailability of dispersed tablets was estimated to be 107% of whole tablets, with a 90% confidence interval of 99.7-114%, fulfilling bioequivalence criteria. The two formulations were not significantly different regarding either bioavailability or its variability. Bioavailability increased at lower doses, by 34% (26-42%) at 50 mg and by 74% (64-86%) at 25 mg, relative to 100 mg. The majority of participants (93%) found the dispersed formulation acceptable in palatability across all delamanid doses.

CONCLUSIONS

Dispersed 50-mg delamanid tablets have similar bioavailability to tablets swallowed whole in adult volunteers. This can be an option for children and other patients who cannot swallow whole tablets, improving access to treatment.

Prothionamide Dose Optimization Using Population Pharmacokinetics for Multidrug-Resistant Tuberculosis Patients

Prothionamide, a second-line drug for multidrug-resistant tuberculosis (MDR-TB), has been in use for a few decades. However, its pharmacokinetic (PK) profile remains unclear. This study aimed to develop a population PK model for prothionamide and then apply the model to determine the optimal dosing regimen for MDR-TB patients. Multiple plasma samples were collected from 27 MDR-TB patients who had been treated with prothionamide at 2 different study hospitals. Prothionamide was administered according to the weight-band dose regimen (500 mg/day for weight <50 kg and 750 mg/day for weight >50 kg) recommended by the World Health Organization. The population PK model was developed using nonlinear mixed-effects modeling. The probability of target attainment, based on systemic exposure and MIC, was used as a response target. Fixed-dose regimens (500 or 750 mg/day) were simulated to compare the efficacies of various dosing regimens. PK profiles adequately described the two-compartment model with first-order elimination and the transit absorption compartment model with allometric scaling on clearance. All dosing regimens had effectiveness >90% for MIC values <0.4 μg/mL in 1.0-log kill target. However, a fixed dose of 750 mg/day was the only regimen that achieved the target resistance suppression of ≥90% for MIC values of <0.2 μg/mL. In conclusion, fixed-dose prothionamide (750 mg/day), regardless of weight-band, was appropriate for adult MDR-TB patients with weights of 40 to 67 kg.

Diagnostic Accuracy of Therapeutic Drug Monitoring During Tuberculosis Treatment

Tuberculosis (TB) patients co-infected with human immunodeficiency virus (HIV) are more likely to have low blood concentrations of the first-line anti-TB drugs (associated with poor outcomes). Therapeutic drug monitoring (TDM) is recommended for certain TB patient populations at increased risk for a poor outcome. Our objective was to estimate the diagnostic accuracy of a 2-hour TDM serum sample for the first-line anti-TB drugs among HIV/TB patients, and evaluate the information gained by an additional 6-hour sample. We created a virtual (n = 1,000) HIV/TB patient population and performed pharmacokinetic (PK) simulations using published population models for isoniazid, rifampin, pyrazinamide, and ethambutol. We performed receiver-operating-characteristic (ROC) analysis to compare the diagnostic performance of a single 2-hour serum sample with samples obtained at 2- and 6-hours post-dosing. The sensitivity of a single 2-hour serum concentration to identify HIV/TB patients with adequate serum exposures was lowest for rifampin (54.9%, 95% CI 50.79-59.41%) and highest for ethambutol (70.8%, 95% CI 66.06-72.61%) for C_max targets. Diagnostic accuracy of a single 2-hour serum sample for the AUC_0-24 target was highest for isoniazid (93%, 95% CI 90.9-94.1%) and lowest for pyrazinamide (66.3%, 95% CI 62.6-70.0%). In summary, the diagnostic performance of TDM for C_max and AUC_0-24 targets demonstrated variability across the first-line anti-TB drugs. The addition of a 6-hour serum sample lead to the highest statistically significant improvement (p < 0.001) and highest increase in diagnostic accuracy (area under the ROC curve) for rifampin for C_max and AUC. The other first-line drugs had modest/negligible increases in diagnostic accuracy. This article is protected by copyright. All rights reserved.

Tutorial: Statistical Analysis and Reporting of Clinical Pharmacokinetic Studies

Pharmacokinetics is the cornerstone of understanding drug absorption, distribution, metabolism, and elimination. It is also the key to describing variability in drug response caused by drug‐drug interactions (DDIs), pharmacogenetics, impaired kidney and liver function, etc. This tutorial aims to provide a guideline and step‐by‐step tutorial on essential considerations when designing clinical pharmacokinetic studies and reporting results. This includes a comprehensive guide on how to conduct the statistical analysis and a complete code for the statistical software R. As an example, we created a mock dataset simulating a clinical pharmacokinetic DDI study with 12 subjects who were administered 2 mg oral midazolam with and without an inducer of cytochrome P450 3A. We provide a step‐by‐step guide to the statistical analysis of this clinical pharmacokinetic study, including sample size/power calculation, descriptive statistics, noncompartmental analyses, and hypothesis testing. The different analyses and parameters are described in detail, and we provide a complete R code ready to use in supplementary files. Finally, we discuss important considerations when designing and reporting clinical pharmacokinetic studies. The scope of this tutorial is not limited to DDI studies, and with minor adjustments, it applies to all types of clinical pharmacokinetic studies. This work was done by early career researchers for early career researchers. We hope this tutorial may help early career researchers when getting started on their own pharmacokinetic studies. We encourage you to use this as an inspiration and starting point and continuously evolve your statistical skills.

Postprandial Responses on Serum Metabolome to Milk and Yogurt Intake in Young and Older Men

The identification and validation of biomarkers of food intake (BFIs) is a promising approach to develop more objective and complementary tools to the traditional dietary assessment methods. Concerning dairy, their evaluation in terms of intake is not simple, given the variety of existing foods, making it difficult to establish the association between specific dairy products consumption and the effects on human health, which is also dependent on the study population. Here, we aimed at identifying BFI of both milk (M) and yogurt (Y) in 14 healthy young (20–35 years) and 14 older (65–80 years). After a 3-week run-in period of dairy exclusion from the diet, the subjects acutely consumed 600 ml of M or Y. Metabolomics analyses were conducted on serum samples during the following 6 h (LC-MS and GC-MS). Several metabolites showing increased iAUC after milk or yogurt intake were considered as potential BFI, including lactose (M &gt; Y, 2-fold), galactitol (M &gt; Y, 1.5-fold), galactonate (M &gt; Y, 1.2-fold), sphingosine-1-phosphate (M &gt; Y from 2.1-fold), as well as an annotated disaccharide (Y &gt; M, 3.6-fold). Delayed serum kinetics were also observed after Y compared to M intake lysine (+22 min), phenylalanine (+45 min), tyrosine (+30min), threonine (+38 min) 3-phenyllactic acid (+30 min), lactose (+30 min), galactitol (+45min) and galactonate (+30 min). The statistical significance of certain discriminant metabolites, such as sphingosine-1-phosphate and several free fatty acids, was not maintained in the older group. This could be related to the physiological modifications induced by aging, like dysregulated lipid metabolism, including delayed appearance of dodecanoic acid (+60 min) or altered postprandial appearance of myristic acid (+70% Cmax), 3-dehydroxycarnitine (−26% Cmin), decanoylcarnitine (−51% Cmin) and dodecanoylcarnitine (−40% Cmin). In conclusion, candidate BFI of milk or yogurt could be identified based on the modified postprandial response resulting from the fermentation of milk to yogurt. Moreover, population specificities (e.g., aging) should also be considered in future studies to obtain more accurate and specific BFI.

Model-Based Equivalent Dose Optimization to Develop New Donepezil Patch Formulation

Donepezil patch was developed to replace the original oral formulation. To accurately describe the pharmacokinetics of donepezil and investigate compatible doses between two formulations, a population pharmacokinetic model for oral and transdermal patches was built based on a clinical study. Plasma donepezil levels were analyzed via liquid chromatography/tandem mass spectrometry. Non-compartmental analyses were performed to derive the initial parameters for compartmental analyses. Compartmental analysis (CA) was performed with NLME software NONMEM assisted by Perl-speaks-NONMEM, and R. Model evaluation was proceeded via visual predictive checks (VPC), goodness-of-fit (GOF) plotting, and bootstrap method. The bioequivalence test was based on a 2 × 2 crossover design, and parameters of AUC and C_max were considered. We found that a two-compartment model featuring two transit compartments accurately describes the pharmacokinetics of nine subjects administered in oral, as well as of the patch-dosed subjects. Through evaluation, the model was proven to be sufficiently accurate and suitable for further bioequivalence tests. Based on the bioequivalence test, 114 mg/101.3 cm²–146 mg/129.8 cm² of donepezil patch per week was equivalent to 10 mg PO donepezil per day. In conclusion, the pharmacokinetic model was successfully developed, and acceptable parameters were estimated. However, the size calculated by an equivalent dose of donepezil patch could be rather large. Further optimization in formulation needs to be performed to find appropriate usability in clinical situations.

Characterizing Absorption Properties of Dispersible Pretomanid Tablets Using Population Pharmacokinetic Modelling

BACKGROUND AND INTRODUCTION

The dispersible tablet formulation (DTF) of pretomanid has been developed to facilitate future use in children. This work aimed to assess the pharmacokinetics (PK) and relative bioavailability of the DTF compared to the marketed formulation (MF) and the potential influence of dose.

METHODS

Pretomanid DTF was investigated in a single-dose, randomized, four-period, cross-over study, with 7 days of washout between doses. Forty-eight healthy volunteers were enrolled and randomized into one of two panels to receive doses either in the fasted state or after a high-fat meal. Each volunteer received doses of 10, 50, and 200 mg DTF, and 200 mg MF pretomanid. Blood samples for pharmacokinetic assessment were drawn following a rich schedule up to 96 h after each single dose. The study data from the panel receiving the high-fat meal were analyzed using a nonlinear mixed-effects modeling approach, and all data were characterized with noncompartmental methods.

RESULTS

A one-compartment model with first-order elimination and absorption through a transit compartment captured the mean and variability of the observed pretomanid concentrations with acceptable precision. No significant difference in bioavailability was found between formulations. The mean absorption time for the DTF was typically 137% (86-171%) of that for the MF. The bioavailability was found to be dose dependent with a small positive and larger negative correlation under fed and fasted conditions, respectively.

CONCLUSION

Using data from a relative bioavailability study in healthy adult volunteers, a mathematical model has been developed to inform dose selection for the investigation of pretomanid in children using the new dispersible tablet formulation. Under fed conditions and at the currently marketed adult dose of 200 mg, the formulation type was found to influence the absorption rate, but not the bioavailability. The bioavailability of the DTF was slightly positively correlated with doses when administered with food.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT04309656, first posted on 16 March 2020.

Model-based efficacy and toxicity comparisons of moxifloxacin for multi-drug-resistant tuberculosis

Background

Moxifloxacin (MOX) is used as a first-choice drug to treat multi-drug-resistant tuberculosis (MDR-TB), however, evidence-based dosing optimization should be strengthened by integrative analysis. The primary goal of this study was to evaluate MOX efficacy and toxicity using integratvie model-based approaches in MDR-TB patients.

Methods

In total, 113 MDR-TB patients from five different clinical trials were analyzed for the development of a population pharmacokinetics (PK) model. A final population PK model was merged with a previously developed lung-lesion distribution and QT prolongation model. Monte Carlo simulation was used to calculate the probability target attainment (PTA) value based on concentration. An area under the concentration-time curve (AUC)-based target was identified as the minimum inhibitory concentration (MIC) of MOX isolated from MDR-TB patients.

Results

The presence of human immunodeficiency virus (HIV) increased clearance by 32.7% and decreased the AUC by 27.4%, compared with HIV-negative MDR-TB patients. A daily dose of 800 mg or a 400 mg twice daily dose of MOX is expected to be effective in MDR-TB patients with an MIC of ≤ 0.25 µg/mL, regardless of PK differences resulting from the presence of HIV. The effect of MOX in HIV-positive MDR-TB patients tended to be decreased dramatically from 0.5 µg/mL, in contrast to the findings in HIV-negative patients. A regimen of twice-daily doses of 400 mg should be considered safer than an 800 mg once-daily dosing regimen, because of the narrow fluctuation of concentrations.

Conclusions

Our results suggest that a 400 mg twice-daily dose of MOX is an optimal dosing regimen for MDR-TB patients because it provides superior efficacy and safety.

Pharmacokinetic Evaluation of Eltrombopag in ITP Pediatric Patients

Background: Eltrombopag (EPAG) is an oral thrombopoietin receptor agonist, approved for refractory primary immune thrombocytopenia (ITP) in pediatric patients. In two pediatric RCTs, EPAG led to an improvement of platelet counts and a reduction in bleeding severity. However, a significant number of pediatric patients did not achieve the primary endpoints. We performed a pharmacokinetic evaluation of EPAG in pediatric patients with refractory ITP. Methods: Outpatients aged from 1 to 17 y, affected by refractory ITP to first-line treatment, were enrolled for a pharmacokinetic assessment. The analysis of drug plasma concentration was performed by the LC-MS/MS platform. Non-compartmental and statistical subgroup analyses were carried out using the R package ncappc. Results: Among 36 patients eligible for PK analysis, the median dose of EPAG given once daily was 50 mg. The EPAG peak occurs between 2 and 4 h with a population Cmax and AUC 0-24 geo-mean of 23, 38 μg/ml, and 275, 4 µg*h/mL, respectively. The pharmacokinetic profile of EPAG did not show a dose proportionality. Female patients showed a statistically significant increase of dose-normalized exposure parameters, increasing by 110 and 123% for Cmax and AUC 0-24, respectively, when compared to male patients. Patients aged 1-5 y showed values increased by more than 100% considering both exposure parameters, compared to older children. Furthermore, patients presenting complete response (83%), showed augmented EPAG exposure parameters compared to subjects with partial or no response. Conclusion: These data highlight the need to further explore the variability of EPAG exposure and its pharmacokinetic/pharmacodynamic profile in pediatric patients also in a real-life setting.

Towards prevention of ischemia-reperfusion kidney injury: Pre-clinical evaluation of 6-chromanol derivatives and the lead compound SUL-138✰

Acute kidney injury (AKI) is a global healthcare burden attributable to high mortality and staggering costs of dialysis. The underlying causes of AKI include hypothermia and rewarming (H/R), ischemia/reperfusion (I/R), mitochondrial dysfunction and reactive oxygen species production. Inspired by the mechanisms conferring organ protection in hibernating hamster, 6-chromanol derived compounds were developed to address the need of effective prevention and treatment of AKI. Here we report on the pre-clinical screening of 6-chromanol leads that confer protection during I/R to select compounds with favorable profiles for clinical testing in AKI. A library of 6-chromanols (n = 63) was screened in silico for pharmacochemical properties and druggability. Selected compounds (n = 15) were screened for the potency to protect HEK293 cells from H/R cell death and subjected to a panel of in vitro safety assays. Based on these parameters, SUL-138 was selected as the lead compound and was found to safeguard kidney function and decrease renal injury after I/R in rats. The compound was without cardiovascular or respiratory effects in vivo. SUL-138 pharmacokinetics of control animals (mouse, rat) and those undergoing I/R (rat) was identical, showing a two-phase elimination profile with terminal half-life of about 8 h. Collectively, our phenotype-based screening approach led to the identification of 3 candidates for pre-clinical studies (5%, 3/64). SUL-138 emerged from this small-scale library of 6-chromanols as a novel prophylactic for AKI. The presented efficacy and safety data provide a basis for future development and clinical testing. SECTION ASSIGNMENTS: : Drug discovery and translational medicine, renal, metabolism SIGNIFICANCE STATEMENT: : Based on in silico druggability parameters, a 63 compound 6-chromanol library was narrowed down to 15 compounds. These compounds were subjected to phenotypical screening of cell survival following hypothermia damage and hit compounds were identified. After subsequent assessment of in vivo efficacy, toxicity, pharmacokinetics, and cardiovascular and respiratory safety, SUL-138 emerged as a lead compound that prevented kidney injury after ischemia/reperfusion and demonstrated a favorable pharmacokinetic profile unaffected by renal ischemia.

Population Modelling of Dexmedetomidine Pharmacokinetics and Haemodynamic Effects After Intravenous and Subcutaneous Administration

Background and Objective

Dexmedetomidine is a potent agonist of α₂-adrenoceptors causing dose-dependent sedation in humans. Intravenous dexmedetomidine is commonly used perioperatively, but an extravascular route of administration would be favoured in palliative care. Subcutaneous infusions provide desired therapeutic plasma concentrations with fewer unwanted effects as compared with intravenous dosing. We aimed to develop semi-mechanistic population models for predicting pharmacokinetic and pharmacodynamic profiles of dexmedetomidine after intravenous and subcutaneous dosing.

Methods

Non-linear mixed-effects modelling was performed using previously collected concentration and haemodynamic effects data from ten (eight in the intravenous phase) healthy human subjects, aged 19–27 years, receiving 1 µg/kg of intravenous or subcutaneous dexmedetomidine during a 10-min infusion.

Results

The absorption of dexmedetomidine from the subcutaneous injection site, and distribution to local subcutaneous fat tissue was modelled using a semi-physiological approach consisting of a depot and fat compartment, while a two-compartment mammillary model explained further disposition. Dexmedetomidine-induced reductions in plasma norepinephrine concentrations were accurately described by an indirect response model. For blood pressure models, the net effect was specified as hyper- and hypotensive effects of dexmedetomidine due to vasoconstriction on peripheral arteries and sympatholysis mediated via the central nervous system, respectively. A heart rate model combined the dexmedetomidine-induced sympatholytic effect, and input from the central nervous system, predicted from arterial blood pressure levels. Internal evaluation confirmed the predictive performance of the final models, as well as the accuracy of the parameter estimates with narrow confidence intervals.

Conclusions

Our final model precisely describes dexmedetomidine pharmacokinetics and accurately predicts dexmedetomidine-induced sympatholysis and other pharmacodynamic effects. After subcutaneous dosing, dexmedetomidine is taken up into subcutaneous fat tissue, but our simulations indicate that accumulation of dexmedetomidine in this compartment is insignificant.

ClinicalTrials.org

NCT02724098 and EudraCT 2015-004698-34

Electronic supplementary material

The online version of this article (10.1007/s40262-020-00900-3) contains supplementary material, which is available to authorized users.

Rifampicin Can Be Given as Flat-Dosing Instead of Weight-Band Dosing

BACKGROUND

The weight-band dosing in tuberculosis treatment regimen has been implemented in clinical practice for decades. Patients will receive different number of fixed dose combination tablets according to their weight-band. However, some analysis has shown that weight was not the best covariate to explain variability of rifampicin exposure. Furthermore, the rationale for using weight-band dosing instead of flat-dosing becomes questionable. Therefore, this study aimed to compare the average and the variability of rifampicin exposure after weight-band dosing and flat-dosing.

METHODS

Rifampicin exposure were simulated using previously published population pharmacokinetics model at dose 10-40 mg/kg for weight-band dosing and dose 600-2400 mg for flat-dosing. The median area under the curve (AUC0-24 h) after day 7 and 14 were compared as well as the variability of each dose group between weight-band and flat-dosing.

RESULTS

The difference of median AUC0-24 h of all dose groups between flat-dosing and weight-band dosing were considered low (< 20%) except for the lowest dose. At the dose of 10 mg/kg (600 mg for flat-dosing), flat-dosing resulted in higher median AUC0-24h compared to the weight-band dosing. A marginal decrease in between-patient variability was predicted for weight-band dosing compared to flat-dosing.

CONCLUSIONS

Weight-band dosing yields a small and non-clinically relevant decrease in variability of AUC0-24h.

Physiologically Based Pharmacokinetic Modeling Framework to Predict Neonatal Pharmacokinetics of Transplacentally Acquired Emtricitabine, Dolutegravir, and Raltegravir

BACKGROUND AND OBJECTIVE

Little is understood about neonatal pharmacokinetics immediately after delivery and during the first days of life following intrauterine exposure to maternal medications. Our objective was to develop and evaluate a novel, physiologically based pharmacokinetic modeling workflow for predicting perinatal and postnatal disposition of commonly used antiretroviral drugs administered prenatally to pregnant women living with human immunodeficiency virus.

METHODS

Using previously published, maternal-fetal, physiologically based pharmacokinetic models for emtricitabine, dolutegravir, and raltegravir built with PK-Sim/MoBi^®, placental drug transfer was predicted in late pregnancy. The total drug amount in fetal compartments at term delivery was estimated and subsequently integrated as initial conditions in different tissues of a whole-body, neonatal, physiologically based pharmacokinetic model to predict drug concentrations in the neonatal elimination phase after birth. Neonatal elimination processes were parameterized according to published data. Model performance was assessed by clinical data.

RESULTS

Neonatal physiologically based pharmacokinetic models generally captured the initial plasma concentrations after delivery but underestimated concentrations in the terminal phase. The mean percentage error for predicted plasma concentrations was - 71.5%, - 33.8%, and 76.7% for emtricitabine, dolutegravir, and raltegravir, respectively. A sensitivity analysis suggested that the activity of organic cation transporter 2 and uridine diphosphate glucuronosyltransferase 1A1 during the first postnatal days in term newborns is ~11% and ~30% of that in adults, respectively.

CONCLUSIONS

These findings demonstrate the general feasibility of applying physiologically based pharmacokinetic models to predict washout concentrations of transplacentally acquired drugs in newborns. These models can increase the understanding of pharmacokinetics during the first postnatal days and allow the prediction of drug exposure in this vulnerable population.

Randomized, Double-Blind, Pharmacokinetic Equivalence Trial Comparing DRL-Rituximab With MabThera in Patients With Diffuse Large B-Cell Lymphoma

PURPOSE

We sought to compare the pharmacokinetics (PKs) of DRL-rituximab (DRL_RI; potential biosimilar) and innovator rituximab MabThera (Roche, Grenzach-Wyhlen, Germany; reference medicinal product [RMP]) in patients with diffuse large B-cell lymphoma (DLBCL). Efficacy, pharmacodynamics (PDs), safety, and immunogenicity were also compared.

PATIENTS AND METHODS

We conducted a double-blind, parallel-group study in patients with untreated DLBCL who were eligible to receive cyclophosphamide, doxorubicin, vincristine, and prednisone (CHOP) therapy. Patients were randomly assigned at a one-to-one ratio to receive DRL_RI or RMP for six 21-day cycles of rituximab plus CHOP, with 18 months of follow-up after day 1, cycle 6 (C6). Primary end point was C1 PKs, measured as area under the plasma concentration–time curve from day 0 to 21 (AUC_{0-21 days}) and maximum plasma concentration (C_max). Equivalence was defined as 90% CIs for the DRL_RI/RMP geometric mean ratios (GMRs) within 80% and 125%. Secondary end points included efficacy noninferiority measured by objective response rate (ORR) at C6 and event-free survival and overall survival at 87 weeks, PK equivalence at C6 and PD equivalence (rate of B-cell depletion and repletion), safety, and immunogenicity. The trial was stopped after sufficient patients for primary end point evaluation were enrolled. Secondary end points are reported as observed.

RESULTS

A total of 151 patients were randomly assigned (DRL_RI, n = 76; RMP, n = 75). DRL_RI/RMP GMRs for AUC_{0-21 days} and C_max in C1 were 99.77 (90% CI, 87.60 to 113.63) and 96.19 (90% CI, 88.65 to 104.38), respectively. ORR at C6 for DRL_RI and RMP were 82.0% and 84.8%, respectively. Rates of B-cell depletion/repletion, immunogenicity, and adverse events were comparable in both groups.

CONCLUSION

DRL_RI and RMP had equivalent PKs, with comparable efficacy, PDs, safety, and immunogenicity.

Pharmacokinetics, Safety, and Efficacy of an Allometric Miltefosine Regimen for the Treatment of Visceral Leishmaniasis in Eastern African Children: An Open-label, Phase II Clinical Trial

BACKGROUND

Convenient, safe, and effective treatments for visceral leishmaniasis in Eastern African children are lacking. Miltefosine, the only oral treatment, failed to achieve adequate efficacy, particularly in children, in whom linear dosing (2.5 mg/kg/day for 28 days) resulted in a 59% cure rate, with lower systemic exposure than in adults.

METHODS

We conducted a Phase II trial in 30 children with visceral leishmaniasis, aged 4-12 years, to test whether 28 days of allometric miltefosine dosing safely achieves a higher systemic exposure than linear dosing.

RESULTS

Miltefosine accumulated during treatment. Median areas under the concentration time curve from days 0-210 and plasma maximum concentration values were slightly higher than those reported previously for children on linear dosing, but not dose-proportionally. Miltefosine exposure at the start of treatment was increased, with higher median plasma concentrations on day 7 (5.88 versus 2.67 μg/mL). Concentration-time curves were less variable, avoiding the low levels of exposure observed with linear dosing. The 210-day cure rate was 90% (95% confidence interval, 73-98%), similar to that previously described in adults. There were 19 treatment-related adverse events (AEs), but none caused treatment discontinuation. There were 2 serious AEs: both were unrelated to treatment and both patients were fully recovered.

CONCLUSIONS

Allometric miltefosine dosing achieved increased and less-variable exposure than linear dosing, though not reaching the expected exposure levels. The new dosing regimen safely increased the efficacy of miltefosine for Eastern African children with visceral leishmaniasis. Further development of miltefosine should adopt allometric dosing in pediatric patients.

CLINICAL TRIALS REGISTRATION

NCT02431143.

Pharmacokinetics of meloxicam after oral administration of a granule formulation to healthy horses

Background

Nonsteroidal anti‐inflammatory drugs are administered in horses for several systemic diseases. Selective cyclooxygenase‐2 inhibitors are preferred because of lower risk of adverse effects. Several meloxicam formulations have been tested in horses, but a recently marketed granule oral formulation has not been studied.

Objective

To characterize the pharmacokinetics of a novel granule meloxicam formulation in fasted and fed horses, and to compare pharmacokinetic features with oral suspension and tablets.

Animals

Seven healthy adult horses.

Methods

Meloxicam was administered at 0.6 mg/kg in fasted or fed horses. Blood samples were collected for pharmacokinetic analysis, and vital signs, hematology, and biochemistry variables were monitored for 72 hours.

Results

No adverse effects were detected. Volume of distribution and clearance after intravenous administration of meloxicam were 0.36 L/kg and 29.12 mL/h/kg, respectively, with a 12.39 hours of terminal half‐life. Protein binding was of 97%. Bioavailability was high for every oral formulation, ranging 70%‐110%, without feed effect. Because of a slower absorption, meloxicam after administration of granules had a longer half‐life (24 and 34 hours, fasted and fed, respectively) and mean residence time (31 and 47 hours), than suspension and tablets (ranging 10‐13 and 13‐15 hours, respectively). In addition, the time above therapeutic concentration was higher for the granule formulation than other formulations.

Conclusions and Clinical Importance

Granule formulation has different PK parameters compared to other oral formulations, which could enable this formulation to be used for different dosage regimens in order to reach a desired clinical effect or decrease the risk of adverse effects.

Semimechanistic Population Pharmacokinetic Model to Predict the Drug-Drug Interaction Between S-ketamine and Ticlopidine in Healthy Human Volunteers

Low‐dose oral S‐ketamine is increasingly used in chronic pain therapy, but extensive cytochrome P450 (CYP) mediated metabolism makes it prone to pharmacokinetic drug‐drug interactions (DDIs). In our study, concentration‐time data from five studies were used to develop a semimechanistic model that describes the ticlopidine‐mediated inhibition of S‐ketamine biotransformation. A mechanistic model was implemented to account for reversible and time‐dependent hepatic CYP2B6 inactivation by ticlopidine, which causes elevated S‐ketamine exposure in vivo. A pharmacokinetic model was developed with gut wall and hepatic clearances for S‐ketamine, its primary metabolite norketamine, and ticlopidine. Nonlinear mixed effects modeling approach was used (NONMEM version 7.3.0), and the final model was evaluated with visual predictive checks and the sampling‐importance‐resampling procedure. Our final model produces biologically plausible output and demonstrates that ticlopidine is a strong inhibitor of CYP2B6 mediated S‐ketamine metabolism. Simulations from our model may be used to evaluate chronic pain therapy with S‐ketamine.

A Population Pharmacokinetic-Pharmacodynamic Model of Pegfilgrastim

Neutropenia and febrile neutropenia (FN) are serious side effects of cytotoxic chemotherapy which may be alleviated with the administration of recombinant granulocyte colony-stimulating factor (GCSF) derivatives, such as pegfilgrastim (PG) which increases absolute neutrophil count (ANC). In this work, a population pharmacokinetic-pharmacodynamic (PKPD) model was developed based on data obtained from healthy volunteers receiving multiple administrations of PG. The developed model was a bidirectional PKPD model, where PG stimulated the proliferation, maturation, and margination of neutrophils and where circulating neutrophils in turn increased the elimination of PG. Simulations from the developed model show disproportionate changes in response with changes in dose. A dose increase of 10% from the 6 mg therapeutic dose taken as a reference leads to area under the curve (AUC) increases of ~50 and ~5% for PK and PD, respectively. A full random effects covariate model showed that little of the parameter variability could be explained by sex, age, body size, and race. As a consequence, little of the secondary parameter variability (C_max and AUC of PG and ANC) could be explained by these covariates.

Pharmacokinetics of metformin in patients with gastrointestinal intolerance

AIMS

To assess potential causes of metformin intolerance, including altered metformin uptake from the intestine, increased anaerobic glucose utilization and subsequent lactate production, altered serotonin uptake, and altered bile acid pool.

METHODS

For this pharmacokinetic study, we recruited 10 severely intolerant and 10 tolerant individuals, matched for age, sex and body mass index. A single 500-mg dose of metformin was administered, with blood sampling at 12 time points over 24 hours. Blood samples were analysed for metformin, lactate, serotonin and bile acid concentrations, and compared across the phenotypes.

RESULTS

The intolerant individuals were severely intolerant to 500 mg metformin. No significant difference was identified between tolerant and intolerant cohorts in metformin pharmacokinetics: median (interquartile range [IQR]) peak concentration 2.1 (1.7-2.3) mg/L and 2.0 (1.8-2.2) mg/L, respectively (P = .76); time to peak concentration 2.5 hours; median (IQR) area under the curve (AUC)_0-24 16.9 (13.9-18.6) and 13.9 (12.9-16.8) mg/L*h, respectively (P = .72). Lactate concentration peaked at 3.5 hours, with mean peak concentration of 2.4 mmol/L in both cohorts (95% CI 2.0-2.8 and 1.8-3.0 mmol/L, respectively), and similar incremental AUC_0-24 in each cohort: tolerant cohort 6.98 (95% CI 3.03-10.93) and intolerant cohort 4.47 (95% CI -3.12-12.06) mmol/L*h (P = .55). Neither serotonin nor bile acid concentrations were significantly different.

CONCLUSIONS

Despite evidence of severe intolerance in our cohort, there was no significant difference in metformin pharmacokinetics or systemic measures of lactate, serotonin or bile acids. This suggests that metformin intolerance may be attributable to local factors within the lumen or enterocyte.

Development of R packages: 'NonCompart' and 'ncar' for noncompartmental analysis (NCA)

Noncompartmental analysis (NCA) is a primary analytical approach for pharmacokinetic studies, and its parameters act as decision criteria in bioequivalent studies. Currently, NCA is usually carried out by commercial softwares such as WinNonlin®. In this article, we introduce our newly-developed two R packages, NonCompart (NonCompartmental analysis for pharmacokinetic data) and ncar (NonCompartmental Analysis for pharmacokinetic Report), which can perform NCA and produce complete NCA reports in both pdf and rtf formats. These packages are compatible with CDISC (Clinical Data Interchange Standards Consortium) standard as well. We demonstrate how the results of WinNonlin® are reproduced and how NCA reports can be obtained. With these R packages, we aimed to help researchers carry out NCA and utilize the output for early stages of drug development process. These R packages are freely available for download from the CRAN repository.

Dicloxacillin induces CYP2C19, CYP2C9 and CYP3A4 in vivo and in vitro

AIM

The aim of this study was to study potential cytochrome P450 (CYP) induction by dicloxacillin.

METHODS

We performed an open-label, randomized, two-phase, five-drug clinical pharmacokinetic cocktail crossover study in 12 healthy men with and without pretreatment with 1 g dicloxacillin three times daily for 10 days. Plasma and urine were collected over 24 h and the concentration of all five drugs and their primary metabolites was determined using a liquid chromatography coupled to triple quadrupole mass spectrometry method. Cryopreserved primary human hepatocytes were exposed to dicloxacillin for 48 h and changes in gene expression and the activity of CYP3A4, CYP2C9, CYP2B6 and CYP1A2 were investigated. The activation of nuclear receptors by dicloxacillin was assessed using luciferase assays.

RESULTS

A total of 10 days of treatment with dicloxacillin resulted in a clinically and statistically significant reduction in the area under the plasma concentration-time curve from 0 to 24 h for omeprazole (CYP2C19) {geometric mean ratio [GMR] [95% confidence interval (CI)]: 0.33 [0.24, 0.45]}, tolbutamide (CYP2C9) [GMR (95% CI): 0.73 (0.65, 0.81)] and midazolam (CYP3A4) [GMR (95% CI): 0.54 (0.41, 0.72)]. Additionally, other relevant pharmacokinetic parameters were affected, indicating the induction of CYP2C- and CYP3A4-mediated metabolism by dicloxacillin. Investigations in primary hepatocytes showed a statistically significant dose-dependent increase in CYP expression and activity by dicloxacillin, caused by activation of the pregnane X receptor.

CONCLUSIONS

Dicloxacillin is an inducer of CYP2C- and CYP3A-mediated drug metabolism, and we recommend caution when prescribing dicloxacillin to users of drugs with a narrow therapeutic window.

Suboptimal Antituberculosis Drug Concentrations and Outcomes in Small and HIV-Coinfected Children in India: Recommendations for Dose Modifications

This work aimed to evaluate the once‐daily antituberculosis treatment as recommended by the new Indian pediatric guidelines. Isoniazid, rifampin, and pyrazinamide concentration–time profiles and treatment outcome were obtained from 161 Indian children with drug‐sensitive tuberculosis undergoing thrice‐weekly dosing as per previous Indian pediatric guidelines. The exposure–response relationships were established using a population pharmacokinetic‐pharmacodynamic approach. Rifampin exposure was identified as the unique predictor of treatment outcome. Consequently, children with low body weight (4–7 kg) and/or HIV infection, who displayed the lowest rifampin exposure, were associated with the highest probability of unfavorable treatment (therapy failure, death) outcome (P_unfavorable). Model‐based simulation of optimized (P_unfavorable ≤ 5%) rifampin once‐daily doses were suggested per treatment weight band and HIV coinfection status (33% and 190% dose increase, respectively, from the new Indian guidelines). The established dose‐exposure–response relationship could be pivotal in the development of future pediatric tuberculosis treatment guidelines.

Microdialysis as a tool to determine the local tissue concentration of dicloxacillin in man

AIMS

The most common pathogen to cause postoperative infections in Denmark is Staphylococcus aureus. Despite using prophylactic antibiotics, infections are still seen. Whether the tissue concentration is above the minimal inhibitory concentration (MIC) for the pathogen is unknown. Thus, the concentration of dicloxacillin in muscle and adipose tissue was measured after intravenous administration, in healthy men.

METHODS

MIC for dicloxacillin against S. aureus was determined using the broth macrodilution method. A microdialysis (MD) catheter was placed in the subcutaneous tissue of the abdomen and in the lateral vastus muscle of the thigh of six healthy male volunteers. They were given 2 g dicloxacillin intravenously. Samples from blood and MD fluid were collected. The unbound dicloxacillin was isolated from plasma. Samples were analysed with high performance liquid chromatography (HPLC).

RESULTS

The maximum concentration was reached in muscle tissue after 0.5 h and in adipose tissue after 0.8 h. AUC_0-6h for the dicloxacillin concentration in adipose tissue was significantly lower when compared to the unbound dicloxacillin concentration in plasma. The dicloxacillin concentration was above the MIC for sensitive S. aureus for a minimum of 2.3 h and a median of 4.1 h in muscle tissue and a minimum of 1.8 h and a median of 3.2 h in adipose tissue.

CONCLUSIONS

The unbound dicloxacillin concentration in adipose and muscle tissue remained above the MIC for sensitive S. aureus, for a period sufficient for many orthopaedic procedures. Whether this is true in patients with compromised circulation remains to be investigated.

Development of a Gastric Absorptive, Immediate Responsive, Oral Protein-Loaded Versatile Polymeric Delivery System

A multifunctional platform to deliver three diverse proteins of insulin, interferon beta (INF-β) and erythropoietin (EPO), using a novel copolymeric microparticulate system of TMC-PEGDMA-MAA, was synthesised as an intelligent pH-responsive 2-fold gastric and intestinal absorptive system. Physiochemical and physicomechanical studies proved the degree of crystallinity that supported the controlled protein delivery of the microparticulate system. The copolymer was tableted before undertaking in vitro and in vivo analysis. After 2.5 h in simulated gastric fluid (SGF), insulin showed a fractional release of 3.2% in comparison to simulated intestinal fluid (SIF), in which a maximum of 83% of insulin was released. Similarly, INF-β and EPO released 3 and 9.7% in SGF and a maximum of 74 and 81.3% in SIF, respectively. In vivo studies demonstrated a significant decrease in blood glucose by 54.19% within 4 h post-dosing, and the comparator formulation provided 74.6% decrease in blood glucose within the same time period. INF-β peak bioavailable dose in serum was calculated to be 1.3% in comparison to an SC formulation having a peak concentration of 0.9%, demonstrating steady-state release for 24 h. EPO-loaded copolymeric microparticles had a 1.6% peak bioavailable concentration, in comparison to the 6.34% peak concentration after 8 h from the SC comparator formulation.

A Population Pharmacokinetic Model Incorporating Saturable Pharmacokinetics and Autoinduction for High Rifampicin Doses

Accumulating evidence suggests that increasing doses of rifampicin may shorten tuberculosis treatment. The PanACEA HIGHRIF1 trial assessed safety, pharmacokinetics, and antimycobacterial activity of rifampicin at doses up to 40 mg/kg. Eighty-three pulmonary tuberculosis patients received 10, 20, 25, 30, 35, or 40 mg/kg rifampicin daily over 2 weeks, supplemented with standard doses of isoniazid, pyrazinamide, and ethambutol in the second week. This study aimed at characterizing rifampicin pharmacokinetics observed in HIGHRIF1 using nonlinear mixed effects modeling. The final population pharmacokinetic model included an enzyme turnover model accounting for time-dependent elimination due to autoinduction, concentration-dependent clearance, and dose-dependent bioavailability. The relationship between clearance and concentration was characterized by a Michaelis-Menten relationship. The relationship between bioavailability and dose was described using an E_max relationship. The model will be key in determining exposure-response relationships for rifampicin and should be considered when designing future trials and when treating future patients with high-dose rifampicin.

Pharmacokinetics of the anti-androgenic drug flutamide in healthy stallions

Alternatives to surgical castration are necessary for controlling the sexual behaviour of stallions with breeding potential in training and competition. Flutamide is a potent selective non-steroidal androgen receptor competitive antagonist that has been used in human beings as an anti-androgenic drug. In this study, the pharmacokinetics and bioavailability of flutamide and its main active metabolite, 2-hydroflutamide, were determined in seven healthy mature stallions. Single doses of flutamide (1mg/kg intravenously, 1mg/kg orally in fasted horses, 5mg/kg orally in fasted horses and 5mg/kg orally in fed horses) were administered randomly at intervals of 2 weeks. All horses had full physical examinations and blood samples were collected for pharmacokinetics, complete blood counts and biochemistry before and after drug administration. Administration of flutamide did not result in any abnormalities on physical examination or in blood parameters. After intravenous administration of flutamide, the volume of distribution was 0.83L/kg and clearance was 1.20L/h/kg. Flutamide and its metabolite had high protein binding values (93-97%). After oral administration, flutamide was rapidly transformed to 2-hydroxyflutamide, with areas under the concentration-time curve ratios of metabolite:drug ∼7. Oral bioavailability was 6.63% after 1mg/kg flutamide in fasted horses, 6.50% after 5mg/kg flutamide in fasted horses and 6.95% after 5mg/kg in fed horses. Half lives of flutamide were close to 1h after intravenous administration and 2h after oral administration. Half lives of 2-hydroxyflutamide were 4.79-6.84h for all routes and doses. After oral administration, oral flutamide reached plasma concentrations that could be effective as an anti-androgenic agent in horses, but further studies are needed to determine whether flutamide has clinical value as an alternative to castration for controlling sexual behaviour in stallions.

Population Pharmacokinetics of Morphine in Patients With Nonalcoholic Steatohepatitis (NASH) and Healthy Adults

Altered expression and function of transporters in nonalcoholic steatohepatitis (NASH) patients may affect the pharmacokinetics (PK), efficacy, and safety of substrate drugs. A population pharmacokinetic (PopPK) analysis was performed to assess differences in morphine and morphine-3-glucuronide (M3G) disposition in NASH and healthy subjects. A total of 315 serum and 42 urine samples from 21 subjects (14 healthy; 7 NASH) were analyzed using NONMEM. Morphine and M3G PK were described by three- and one-compartment models, respectively. After accounting for the effect of total body weight on all clearance and volume of distribution parameters using an allometric scaling approach, NASH severity score (NASF; combination of fibrosis and nonalcoholic fatty liver disease activity scores) was the most significant predictor of differences in M3G exposure. The model predicted a linear decrease in the clearance of M3G with increasing NASF scores on a natural logarithmic scale. These results may provide some insight into the potential effect of NASH on the disposition of hepatic transporter substrates.

Levodopa/carbidopa microtablets in Parkinson's disease: a study of pharmacokinetics and blinded motor assessment

Background

Motor function assessments with rating scales in relation to the pharmacokinetics of levodopa may increase the understanding of how to individualize and fine-tune treatments.

Objectives

This study aimed to investigate the pharmacokinetic profiles of levodopa-carbidopa and the motor function following a single-dose microtablet administration in Parkinson’s disease.

Methods

This was a single-center, open-label, single-dose study in 19 patients experiencing motor fluctuations. Patients received 150% of their individual levodopa equivalent morning dose in levodopa-carbidopa microtablets. Blood samples were collected at pre-specified time points. Patients were video recorded and motor function was assessed with six UPDRS part III motor items, dyskinesia score, and the treatment response scale (TRS), rated by three blinded movement disorder specialists.

Results

AUC_0–4/dose and C_max/dose for levodopa was found to be higher in Parkinson’s disease patients compared with healthy subjects from a previous study, (p = 0.0008 and p = 0.026, respectively). The mean time to maximum improvement in sum of six UPDRS items score was 78 min (±59) (n = 16), and the mean time to TRS score maximum effect was 54 min (±51) (n = 15). Mean time to onset of dyskinesia was 41 min (±38) (n = 13).

Conclusions

In the PD population, following levodopa/carbidopa microtablet administration in fasting state, the C_max and AUC_0–4/dose were found to be higher compared with results from a previous study in young, healthy subjects. A large between subject variability in response and duration of effect was observed, highlighting the importance of a continuous and individual assessment of motor function in order to optimize treatment effect.

Levodopa-entacapone-carbidopa intestinal gel in Parkinson's disease: A randomized crossover study

BACKGROUND

The addition of oral entacapone to levodopa-carbidopa intestinal gel treatment leads to less conversion of levodopa to 3-O-methyldopa, thereby increasing levodopa plasma concentration. The objective of this study was to compare systemic levodopa exposure of the newly developed levodopa-entacapone-carbidopa intestinal gel after a 20% dose reduction with levodopa exposure after the usual levodopa-carbidopa intestinal gel dose in a randomized crossover trial in advanced Parkinson's disease patients.

METHODS

In this 48-hour study, 11 patients treated with levodopa-carbidopa intestinal gel were randomized to a treatment sequence. Blood samples were drawn at prespecified times, and patient motor function was assessed according to the treatment response scale.

RESULTS

Systemic exposure of levodopa did not differ significantly between treatments (ratio, 1.10 [95% confidence interval, 0.951-1.17]). Treatment response scale scores did not significantly differ between treatments (P = 0.84).

CONCLUSIONS

Levodopa-entacapone-carbidopa intestinal gel allowed a lower amount of levodopa administration and was well tolerated. Long-term studies are needed to confirm the results. © 2016 International Parkinson and Movement Disorder Society.

Confirming model-predicted pharmacokinetic interactions between bedaquiline and lopinavir/ritonavir or nevirapine in patients with HIV and drug-resistant tuberculosis

Bedaquiline and its metabolite M2 are metabolised by CYP3A4. The antiretrovirals ritonavir-boosted lopinavir (LPV/r) and nevirapine inhibit and induce CYP3A4, respectively. Here we aimed to quantify nevirapine and LPV/r drug-drug interaction effects on bedaquiline and M2 in patients co-infected with HIV and multidrug-resistant tuberculosis (MDR-TB) using population pharmacokinetic (PK) analysis and compare these with model-based predictions from single-dose studies in subjects without TB. An observational PK study was performed in three groups of MDR-TB patients during bedaquiline maintenance dosing: HIV-seronegative patients (n = 17); and HIV-infected patients using antiretroviral therapy including nevirapine (n = 17) or LPV/r (n = 14). Bedaquiline and M2 samples were collected over 48 h post-dose. A previously developed PK model of MDR-TB patients was used as prior information to inform parameter estimation using NONMEM. The model was able to describe bedaquiline and M2 concentrations well, with estimates close to their priors and earlier model-based interaction effects from single-dose studies. Nevirapine changed bedaquiline clearance to 82% (95% CI 67-99%) and M2 clearance to 119% (92-156%) of their original values, indicating no clinically significant interaction. LPV/r substantially reduced bedaquiline clearance to 25% (17-35%) and M2 clearance to 59% (44-69%) of original values. This work confirms earlier model-based predictions of nevirapine and LPV/r interaction effects on bedaquiline and M2 clearance from subjects without TB in single-dose studies, in MDR-TB/HIV co-infected patients studied here. To normalise bedaquiline exposure in patients with concomitant LPV/r therapy, an adjusted bedaquiline dosing regimen is proposed for further study.