Mechanism-based concepts of size and maturity in pharmacokinetics

Predicting Antimicrobial Activity at the Target Site: Pharmacokinetic/Pharmacodynamic Indices versus Time-Kill Approaches

Antibiotic dosing strategies are generally based on systemic drug concentrations. However, drug concentrations at the infection site drive antimicrobial effect, and efficacy predictions and dosing strategies should be based on these concentrations. We set out to review different translational pharmacokinetic-pharmacodynamic (PK/PD) approaches from a target site perspective. The most common approach involves calculating the probability of attaining animal-derived PK/PD index targets, which link PK parameters to antimicrobial susceptibility measures. This approach is time efficient but ignores some aspects of the shape of the PK profile and inter-species differences in drug clearance and distribution, and provides no information on the PD time-course. Time–kill curves, in contrast, depict bacterial response over time. In vitro dynamic time–kill setups allow for the evaluation of bacterial response to clinical PK profiles, but are not representative of the infection site environment. The translational value of in vivo time–kill experiments, conversely, is limited from a PK perspective. Computational PK/PD models, especially when developed using both in vitro and in vivo data and coupled to target site PK models, can bridge translational gaps in both PK and PD. Ultimately, clinical PK and experimental and computational tools should be combined to tailor antibiotic treatment strategies to the site of infection.

Effective and Safe Daclatasvir Drug Exposures Predicted in Children Using Adult Formulations

BACKGROUND

Sofosbuvir (SOF)/daclatasvir (DCV) is the direct-acting antiviral regimen of choice in many low- and middle-income countries for curative treatment of chronic hepatitis C virus (HCV) infection in adults, but data on the use of DCV in children are lacking. We performed a population pharmacokinetic (PK) analysis to predict DCV exposure in children treated with available adult formulations.

METHODS

DCV concentration data from HCV-infected adolescents receiving SOF/DCV [400/60 mg, once daily (OD)] who participated in a PK study in Egypt were used for model development. PK parameters were estimated using a population approach. Monte Carlo simulations were run for virtual children weighing 10 to <35 kg receiving 60 or 30 mg OD, and DCV exposures were compared with adults ranges.

RESULTS

Seventeen HCV-infected adolescents (13 males) provided 151 DCV concentrations. Median (range) age was 14 (11-18) years and weight 50 (32-63) kg. In these adolescents receiving 60 mg DCV, median (interquartile range) DCV area under the concentration time curve 0 to 24 hours, maximum concentrations, and minimum concentrations were 11,130 (8140-14,690) ng·h/mL, 1030 (790-1220) ng/mL and 130 (110-220) ng/mL, respectively, compared with 10,343 (7661-14,095) ng·h/mL, 1132 (876-1518) ng/mL and 110 (55.7-192) ng/mL predicted in children 10 to <35 kg receiving 30 mg. The proportion of children with DCV exposures above the adult range rapidly increased for children <30 kg using 60 mg OD, similarly for children 10-14 kg using 30 mg.

CONCLUSIONS

DCV 30 mg OD was predicted to achieve effective and safe exposures in children 14 to <35 kg, perhaps down to 10 kg. These results should be validated clinically. Low-cost available adult DCV formulations together with approved pediatric doses of SOF would expand global access to HCV treatment for children.

Population pharmacokinetics of oxycodone: Premature neonates to adults

BACKGROUND

Oxycodone is used in children and adults for the control of acute postoperative pain. Covariate influences such as age, size, and fat mass on oxycodone pharmacokinetic parameters over the human lifespan are poorly quantified.

METHODS

Pooled oxycodone time-concentration profiles were available from preterm neonates to adults. Data from intravenous, intramuscular, buccal, and epidural formulations were analyzed using nonlinear mixed-effects models. Normal fat mass was used to determine the influence of fat on oxycodone pharmacokinetics. Theory-based allometry was used to scale pharmacokinetic parameters to a 70 kg individual. A maturation function described the increase in clearance in neonates and infants.

RESULTS

There were 237 subjects (24 weeks postmenstrual age to 75 years; 0.44-110 kg) providing 1317 plasma concentrations. A three-compartment model with first-order elimination best described oxycodone disposition. Population parameter estimates were clearance (CL) 48.6 L.h^-1 .70 kg^-1 (CV 71%); intercompartmental clearances (Q2) 220 L.h^-1 .70 kg^-1 (CV 64%); Q3 1.45 L.h^-1 .70 kg^-1 ; volume of distribution in the central compartment (V1) 98.2 L.70 kg^-1 (CV 76%); rapidly equilibrating peripheral compartment (V2) 90.1 L. 70 kg^-1 (CV 76%); slow equilibrating peripheral compartment (V3) 28.9 L.70 kg^-1 . Total body weight was the best size descriptor for clearances and volumes. Absorption halftimes (T_ABS ) were: 1.1 minutes for intramuscular, 70 minutes for epidural, 82 minutes for nasogastric, and 159.6 minutes for buccal administration routes. The relative bioavailability after nasogastric administration was 0.673 with a lag time of 8.7 minutes.

CONCLUSIONS

Clearance matured with age; 8% of the typical adult value at 24 weeks postmenstrual age, 33% in a term neonate and reached 90% of the adult clearance value by the end of the first year of life. Allometric scaling using total body weight was the better size descriptor of oxycodone clearance than fat-free mass.

Population pharmacokinetic-pharmacodynamic modeling of acetaminophen in preterm neonates with hemodynamically significant patent ductus arteriosus

BACKGROUND

Pharmacokinetics (PK) of intravenous acetaminophen has not been assessed in preterm neonates with hemodynamically significant patent ductus arteriosus (PDA). Moreover, there is a lack of data evaluating the association between PK and pharmacodynamics (PD) of acetaminophen in hemodynamically significant PDA. Hence, we performed a population PK-PD modeling of acetaminophen in preterm neonates with hemodynamically significant PDA.

METHODS

A prospective, observational study was carried out in preterm neonates with hemodynamically significant PDA receiving intravenous acetaminophen (15 mg/kg six hourly) for maximum of nine days. The diameter of the ductus arteriosus was measured using General Electric Vivid 7® (echocardiography) and was the PD measure. The PK-PD modeling was performed using Monolix 2019R2. We performed Monte Carlo (MC) simulations to determine the probability of ductus arteriosus closure from first to the ninth day of acetaminophen treatment.

RESULTS

Fifty-five neonates were recruited. A one-compartment model with first-order elimination described well the PK of acetaminophen. Clearance (CL) and volume of distribution (Vd) for typical neonate weighing 0.98 kg was 0.0452 L/h and 1.18 L, respectively. A combination of an Imax model with effect compartment and an exponential disease progression model described well the PD of acetaminophen. The average baseline diameter of the ductus arteriosus (E0) was 2.53 mm while IC₅₀ was 0.477 µg/mL. The disease progression rate constant (Kprog) and effect compartment transfer rate constant (ke0) were 0.00425 h^-1 and 0.000103 h^-1, respectively. MC simulations of the current dosing regimen revealed a probability of 73.7% ductus arteriosus closure compared to 83.8% with 20 mg/kg six hourly dose.

CONCLUSION

The PK-PD model developed can be used for dosing acetaminophen in premature neonates with hemodynamically significant PDA. Intravenous dose of 20 mg/kg intravenously every six hours is likely to provide a better therapeutic effect than the existing dosing regimen.

Dexmedetomidine Clearance Decreases with Increasing Drug Exposure: Implications for Current Dosing Regimens and Target-controlled Infusion Models Assuming Linear Pharmacokinetics

BACKGROUND

Numerous pharmacokinetic models have been published aiming at more accurate and safer dosing of dexmedetomidine. The vast majority of the developed models underpredict the measured plasma concentrations with respect to the target concentration, especially at plasma concentrations higher than those used in the original studies. The aim of this article was to develop a dexmedetomidine pharmacokinetic model in healthy adults emphasizing linear versus nonlinear kinetics.

METHODS

The data of two previously published clinical trials with stepwise increasing dexmedetomidine target-controlled infusion were pooled to build a pharmacokinetic model using the NONMEM software package (ICON Development Solutions, USA). Data from 48 healthy subjects, included in a stratified manner, were utilized to build the model.

RESULTS

A three-compartment mamillary model with nonlinear elimination from the central compartment was superior to a model assuming linear pharmacokinetics. Covariates included in the final model were age, sex, and total body weight. Cardiac output did not explain between-subject or within-subject variability in dexmedetomidine clearance. The results of a simulation study based on the final model showed that at concentrations up to 2 ng · ml-1, the predicted dexmedetomidine plasma concentrations were similar between the currently available Hannivoort model assuming linear pharmacokinetics and the nonlinear model developed in this study. At higher simulated plasma concentrations, exposure increased nonlinearly with target concentration due to the decreasing dexmedetomidine clearance with increasing plasma concentrations. Simulations also show that currently approved dosing regimens in the intensive care unit may potentially lead to higher-than-expected dexmedetomidine plasma concentrations.

CONCLUSIONS

This study developed a nonlinear three-compartment pharmacokinetic model that accurately described dexmedetomidine plasma concentrations. Dexmedetomidine may be safely administered up to target-controlled infusion targets under 2 ng · ml-1 using the Hannivoort model, which assumed linear pharmacokinetics. Consideration should be taken during long-term administration and during an initial loading dose when following the dosing strategies of the current guidelines.

EDITOR’S PERSPECTIVE

Optimization of Initial Dose Regimen for Sirolimus in Pediatric Patients With Lymphangioma

Sirolimus is an effective oral treatment for pediatric patients with lymphangioma. The present clinical study in 15 children (0.12–16.39 years of age) examines the effects of underlying factors on sirolimus concentrations through application of a population pharmacokinetic model. Using Monte Carlo simulation, an initial dose regimen for sirolimus in pediatric patients with lymphangioma is presented. It is found that the lower the body weight, the higher the clearance rate and sirolimus clearances are 0.31–0.17 L/h/kg in pediatric patients with lymphangioma whose weights are 5–60 kg, respectively. The doses of sirolimus, 0.07, 0.06, 0.05 mg/kg/day are recommended for weights of 5–10, 10–24.5 and 24.5–60 kg in children with lymphangioma. This study is the first to establish a population pharmacokinetic model for sirolimus and to recommend initial doses in pediatric patients with lymphangioma. Large scale, prospective studies are needed in the future.

Pharmacokinetic properties of the antimalarial combination therapy artemether-lumefantrine in normal-weight, overweight and obese healthy male adults

The component drugs in the widely used antimalarial artemisinin combination therapy artemether-lumefantrine are lipophilic, with the possibility that recommended fixed doses in adults may lead to subtherapeutic concentrations and consequent treatment failure in overweight/obese individuals with malaria. The aim of this study was to investigate the pharmacokinetic properties of artemether, lumefantrine and their active metabolites dihydroartemisinin and desbutyl-lumefantrine in 16 normal-weight, overweight and obese healthy male volunteers [body mass index (BMI) categories ≤25 kg/m², >25-≤30 kg/m² and >30 kg/m², respectively; absolute range 19.3-37.2 kg/m²]. Participants received the conventional six doses of artemether-lumefantrine over 3 days, each dose comprising 80 mg artemether plus 480 mg lumefantrine administered with 6.7 g fat, and blood samples were collected at pre-specified time-points over 14 days. Plasma drug/metabolite concentrations were measured using liquid chromatography-mass spectrometry and included in multi-compartmental population pharmacokinetic models. There was a non-significant trend to a lower area under the plasma concentration-time curve with a higher body weight or BMI for dihydroartemisinin and especially artemether which was attenuated when normalized for mg/kg dose, but this relationship was not evident in the case of the more lipophilic lumefantrine and its metabolite desbutyl-lumefantrine. Simulated Day 7 plasma lumefantrine concentrations were >200 µg/L (the threshold at which Plasmodium falciparum recrudescences are minimized) in all participants. These results indicate that there is no need for artemether-lumefantrine dose modification in overweight and obese patients with malaria.

Dosing strategies of imipenem in neonates based on pharmacometric modelling and simulation

OBJECTIVES

Imipenem is a broad-spectrum antibacterial agent used in critically ill neonates after failure of first-line treatments. Few studies have described imipenem disposition in this population. The objectives of our study were: (i) to characterize imipenem population pharmacokinetics (PK) in a cohort of neonates; and (ii) to conduct model-based simulations to evaluate the performance of six different dosing regimens aiming at optimizing PK target attainment.

METHODS

A total of 173 plasma samples from 82 neonates were collected over 15 years at the Lausanne University Hospital, Switzerland. The majority of study subjects were preterm neonates with a median gestational age (GA) of 27 weeks (range: 24-41), a postnatal age (PNA) of 21 days (2-153) and a body weight (BW) of 1.16 kg (0.5-4.1). PK data were analysed using non-linear mixed-effect modelling (NONMEM).

RESULTS

A one-compartment model best characterized imipenem disposition. Population PK parameters estimates of CL and volume of distribution were 0.21 L/h and 0.73 L, with an interpatient variability (CV%) of 20.1% on CL in a representative neonate (GA 27 weeks, PNA 21 days, BW 1.16 kg, serum creatinine, SCr 46.6 μmol/L). GA and PNA exhibited the greatest impact on PK parameters, followed by SCr. These covariates explained 36% and 15% of interindividual variability in CL, respectively.Simulated regimens using a dose of 20-25 mg/kg every 6-12 h according to postnatal age led to the highest PTA (T>MIC over 100% of time).

CONCLUSIONS

Dosing adjustment according to BW, GA and PNA optimizes imipenem exposure in neonates.

Linezolid Population Pharmacokinetics in South African Adults with Drug-Resistant Tuberculosis

Linezolid is widely used for drug-resistant tuberculosis (DR-TB) but has a narrow therapeutic index. To inform dose optimization, we aimed to characterize the population pharmacokinetics of linezolid in South African participants with DR-TB and explore the effect of covariates, including HIV coinfection, on drug exposure. Data were obtained from pharmacokinetic substudies in a randomized controlled trial and an observational cohort study, both of which enrolled adults with drug-resistant pulmonary tuberculosis. Participants underwent intensive and sparse plasma sampling. We analyzed linezolid concentration data using nonlinear mixed-effects modeling and performed simulations to estimate attainment of putative efficacy and toxicity targets. A total of 124 participants provided 444 plasma samples; 116 were on the standard daily dose of 600 mg, while 19 had dose reduction to 300 mg due to adverse events. Sixty-one participants were female, 71 were HIV-positive, and their median weight was 56 kg (interquartile range [IQR], 50 to 63). In the final model, typical values for clearance and central volume were 3.57 liters/h and 40.2 liters, respectively. HIV coinfection had no significant effect on linezolid exposure. Simulations showed that 600-mg dosing achieved the efficacy target (area under the concentration-time curve for the free, unbound fraction of the drug [[Formula: see text] at a MIC level of 0.5 mg/liter) with 96% probability but had 56% probability of exceeding safety target ([Formula: see text]. The 300-mg dose did not achieve adequate efficacy exposures. Our model characterized population pharmacokinetics of linezolid in South African patients with DR-TB and supports the 600-mg daily dose with safety monitoring.

Effects of voriconazole on population pharmacokinetics and optimization of the initial dose of tacrolimus in children with chronic granulomatous disease undergoing hematopoietic stem cell transplantation

Background

This study aimed to explore the effects of voriconazole on population pharmacokinetics and optimization of the initial dose of tacrolimus in children with chronic granulomatous disease (CGD) undergoing hematopoietic stem cell transplantation (HSCT).

Methods

Thirty-four children with CGD undergoing HSCT were assessed to establish a population pharmacokinetic model (PPM) using the non-linear mixed effect. Tacrolimus concentrations were simulated by the Monte Carlo method in children weighing <25 kg at different doses.

Results

In the final model, weight and concomitant use of voriconazole were included as covariates. With the same weight, the relative value of tacrolimus clearance was 1:0.388 in children not taking voriconazole: children taking voriconazole. Compared with children not taking voriconazole, the measured tacrolimus concentrations were all higher in children taking voriconazole (P<0.01); however, these were not corrected by dose or body weight for concentration differences. Thus, we simulated the tacrolimus concentrations using different body weights (5–25 kg) and different dose regimens (0.1–0.8 mg/kg/day) for the same body weight and dose. Tacrolimus concentrations in children taking voriconazole were higher than those in children not taking voriconazole (P<0.01). Also, in children with CGD undergoing HSCT who were not taking voriconazole, the initial dose regimen of 0.5 mg/kg/day was recommended for body weights of 5–10 kg, and 0.4 mg/kg/day was recommended for body weights of 10–25 kg. In children with CGD undergoing HSCT who were taking voriconazole, an initial dose regimen of 0.3 mg/kg/day was recommended for body weights of 5–25 kg.

Conclusions

We established, for the first time, a PPM of tacrolimus in children with CGD undergoing HSCT in which voriconazole significantly increased tacrolimus concentrations. In addition, the initial dose of tacrolimus in children with CGD undergoing HSCT was recommended.

Can Published Population Pharmacokinetic Models of Busulfan Be Used for Individualized Dosing in Chinese Pediatric Patients Undergoing Hematopoietic Stem Cell Transplantation? An External Evaluation

Busulfan is a bifunctional alkylating agent that is widely used before hematopoietic stem cell transplantation (HSCT), in combination with other chemotherapeutic drugs. As of 2020, there is no population pharmacokinetic (popPK) model for busulfan in Chinese pediatric patients. A systemic external evaluation of 11 published popPK models was conducted in Chinese pediatric patients undergoing HSCT. Forty pediatric patients were enrolled in this study, with a total of 183 blood concentrations. The relative prediction error (PE%), median PE%, median absolute PE%, and percentage of PE% within ±20% and ±30% were calculated in prediction-based diagnostics. Simulation-based diagnostics were conducted through a prediction- and variability-corrected visual predictive check and the normalized prediction distribution error. The relative individual prediction error was calculated using Bayesian forecasting with 1 to 3 concentration points. The 1-compartment open linear popPK model, which was built by Su-jin Rhee et al (model H), incorporating the patient's body surface area, age, dosing day, and aspartate aminotransferase as significant covariates had preferable predictability than other popPK models. In prediction-based diagnostics, the median PE%, percentage of PE% within ±20%, and percentage of PE% within ±30% of model H were 8.48%, 45.35%, and 59.56%, respectively. The normalized prediction distribution error of model H showed that it followed the normal distribution. Based on Bayesian forecasting, model H showed good predictive performance. Thus, model H was the most appropriate model that can be used clinically for individualized dosage adjustments in Chinese pediatric HSCT patients.

Optimizing dosing and fixed-dose combinations of rifampicin, isoniazid, and pyrazinamide in pediatric patients with tuberculosis: a prospective population pharmacokinetic study

BACKGROUND

In 2010, the WHO revised dosing guidelines for treatment of childhood tuberculosis. Our aim was to investigate first-line antituberculosis drug exposures under these guidelines, explore dose optimization using the current dispersible fixed-dose combination (FDC) table of rifampicin/isoniazid/pyrazinamide; 75/50/150 mg , and suggest a new FDC with revised weight-bands.

METHODS

Children with drug-susceptible tuberculosis in Malawi and South Africa underwent pharmacokinetic sampling while receiving first-line tuberculosis drugs as single formulations according the 2010 WHO recommended doses. Nonlinear mixed-effects modelling and simulation was used to design the optimal FDC and weight-band dosing strategy for achieving the pharmacokinetic targets based on literature-derived adult AUC0-24h for rifampicin (38.7-72.9) isoniazid (11.6-26.3) and pyrazinamide (233-429 mg∙h/L).

RESULTS

180 children (42% female; 13.9% HIV-infected; median [range] age 1.9 [0.22-12] years; weight 10.7 [3.20-28.8] kg) were administered 1, 2, 3, or 4 FDC tablets (rifampicin/isoniazid/pyrazinamide 75/50/150 mg) daily for 4-8, 8-12, 12-16, and 16-25 kg weight-bands, respectively. Rifampicin exposure (for weight and age) was up to 50% lower than in adults. Increasing the tablet number resulted in adequate rifampicin but relatively high isoniazid and pyrazinamide exposures. Administering 1, 2, 3, or 4 optimized FDC tablets (rifampicin/isoniazid/pyrazinamide 120/35/130 mg) to children <6, 6-13, 13-20 and 20-25 kg, and 0.5 tablet in <3-month-olds with immature metabolism, improved exposures to all three drugs.

CONCLUSION

Current pediatric FDC doses resulted in low rifampicin exposures. Optimal dosing of all drugs cannot be achieved with the current FDCs. We propose a new FDC formulation and revised weight-bands.

Population pharmacokinetics of vancomycin in paediatric patients with febrile neutropenia and augmented renal clearance: development of new dosing recommendations

OBJECTIVES

The purpose of this study was to evaluate the influence of augmented renal clearance (ARC) on vancomycin clearance and provide dosage recommendations for paediatric patients with febrile neutropenia following HSCT.

METHODS

A population pharmacokinetic analysis was performed based on a two-compartment model structure using a non-linear mixed-effect modelling approach. Monte Carlo simulations were conducted as a target attainment analysis of AUC between 400 mg·h/L and 650 mg·h/L for MRSA at an MIC of 1 mg/L.

RESULTS

A total of 165 paediatric patients and 276 vancomycin serum concentrations were analysed in this study. Age, body weight, estimated glomerular filtration rate (eGFR) and fever (≥38.0°C) were identified as factors that significantly influenced vancomycin clearance. The median eGFR of the population was 143 mL/min/1.73 m2 and 34% of patients showed an eGFR ≥160 mL/min/1.73 m2, which may be classified as ARC. Our simulations showed that current dosing recommendations result in poor target attainment. In particular, children aged 6 months old to 6 years old with ARC require an initial vancomycin dose up to 35%-65% higher than the current dosing guidelines.

CONCLUSIONS

ARC is frequently observed in paediatric patients with post-HSCT febrile neutropenia, resulting in a significant increase in vancomycin clearance. We propose a vancomycin dosing strategy for children with febrile neutropenia following HSCT based on eGFR, age, weight and body temperature.

Model-based assessment of the safety of community interventions with primaquine in sub-Saharan Africa

BACKGROUND

Single low-dose primaquine (SLD-PQ) is recommended in combination with artemisinin-based combination therapy to reduce Plasmodium falciparum transmission in areas threatened by artemisinin resistance or aiming for malaria elimination. SLD-PQ may be beneficial in mass drug administration (MDA) campaigns to prevent malaria transmission but uptake is limited by concerns of hemolysis in glucose-6-phosphate dehydrogenase (G6PD)-deficient individuals. The aim of this study was to improve the evidence on the safety of MDA with SLD-PQ in a sub-Saharan African setting.

METHODS

A nonlinear mixed-effects model describing the pharmacokinetics and treatment-induced hemolysis of primaquine was developed using data from an adult (n = 16, G6PD deficient) and pediatric study (n = 38, G6PD normal). The relationship between primaquine pharmacokinetics and hemolysis was modeled using an established erythrocyte lifespan model. The safety of MDA with SLD-PQ was explored through Monte Carlo simulations for SLD-PQ at 0.25 or 0.4 mg/kg using baseline data from a Tanzanian setting with detailed information on hemoglobin concentrations and G6PD status.

RESULTS

The predicted reduction in hemoglobin levels following SLD-PQ was small and returned to pre-treatment levels after 25 days. G6PD deficiency (African A- variant) was associated with a 2.5-fold (95% CI 1.2-8.2) larger reduction in hemoglobin levels. In the Tanzanian setting where 43% of the population had at least mild anemia (hemoglobin < 11-13 g/dl depending on age and sex) and 2.73% had severe anemia (hemoglobin < 7-8 g/dl depending on age and sex), an additional 3.7% and 6.0% of the population were predicted to develop at least mild anemia and 0.25% and 0.41% to develop severe anemia after 0.25 and 0.4 mg/kg SLD-PQ, respectively. Children < 5 years of age and women ≥ 15 years of age were found to have a higher chance to have low pre-treatment hemoglobin.

CONCLUSIONS

This study supports the feasibility of MDA with SLD-PQ in a sub-Saharan African setting by predicting small and transient reductions in hemoglobin levels. In a setting where a substantial proportion of the population had low hemoglobin concentrations, our simulations suggest treatment with SLD-PQ would result in small increases in the prevalence of anemia which would most likely be transient.

Pharmacokinetics and Drug-Drug Interactions of Abacavir and Lamuvudine Co-administered With Antituberculosis Drugs in HIV-Positive Children Treated for Multidrug-Resistant Tuberculosis

Given the high prevalence of multidrug-resistant (MDR)-TB in high HIV burden settings, it is important to identify potential drug-drug interactions between MDR-TB treatment and widely used nucleoside reverse transcriptase inhibitors (NRTIs) in HIV-positive children. Population pharmacokinetic models were developed for lamivudine (n = 54) and abacavir (n = 50) in 54 HIV-positive children established on NRTIs; 27 with MDR-TB (combinations of high-dose isoniazid, pyrazinamide, ethambutol, ethionamide, terizidone, fluoroquinolones, and amikacin), and 27 controls without TB. Two-compartment models with first-order elimination and transit compartment absorption described both lamivudine and abacavir pharmacokinetics, respectively. Allometric scaling with body weight adjusted for the effect of body size. Clearance was predicted to reach half its mature value ∼ 2 (lamivudine) and ∼ 3 (abacavir) months after birth, with completion of maturation for both drugs at ∼ 2 years. No significant difference was found in key pharmacokinetic parameters of lamivudine and abacavir when co-administered with routine drugs used for MDR-TB in HIV-positive children.

Wide size dispersion and use of body composition and maturation improves the reliability of allometric exponent estimates

To evaluate study designs and the influence of dispersion of body size, body composition and maturation of clearance or reliable estimation of allometric exponents. Non-linear mixed effects modeling and parametric bootstrap were employed to assess how the study sample size, number of observations per subject, between subject variability (BSV) and dispersion of size distribution affected estimation bias and uncertainty of allometric exponents. The role of covariate model misspecification was investigated using a large data set ranging from neonates to adults. A decrease in study sample size, number of observations per subject, an increase in BSV and a decrease in dispersion of size distribution, increased the uncertainty of allometric exponent estimates. Studies conducted only in adults with drugs exhibiting normal (30%) BSV in clearance may need to include at least 1000 subjects to be able to distinguish between allometric exponents of 2/3 and 1. Nevertheless, studies including both children and adults can distinguish these exponents with only 100 subjects. A marked bias of 45% (95%CI 41-49%) in the estimate of the allometric exponent of clearance was obtained when maturation and body composition were ignored in infants. A wide dispersion of body size (e.g. infants, children and adults) is required to reliably estimate allometric exponents. Ignoring differences in body composition and maturation of clearance may bias the exponent for clearance. Therefore, pharmacometricians should avoid estimating allometric exponent parameters without suitable designs and covariate models. Instead, they are encouraged to rely on the well-developed theory and evidence that clearance and volume parameters in humans scale with theory-based exponents.

Population pharmacokinetics of cefotaxime in intensive care patients

PURPOSE

To characterise the pharmacokinetics and associated variability of cefotaxime in adult intensive care unit (ICU) patients and to assess the impact of patient covariates.

METHODS

This work was based on data from cefotaxime-treated patients included in the ACCIS (Antibiotic Concentrations in Critical Ill ICU Patients in Sweden) study. Clinical data from 51 patients at seven different ICUs in Sweden, given cefotaxime (1000-3000 mg given 2-6 times daily), were collected from the first day of treatment for up to three consecutive days. In total, 263 cefotaxime samples were included in the population pharmacokinetic analysis.

RESULTS

A two-compartment model with linear elimination, proportional residual error and inter-individual variability (IIV) on clearance and central volume of distribution best described the data. The typical individual was 64 years, with body weight at ICU admission of 92 kg and estimated creatinine clearance of 94 mL/min. The resulting typical value of clearance was 11.1 L/h, central volume of distribution 5.1 L, peripheral volume of distribution 18.2 L and inter-compartmental clearance 14.5 L/h. The estimated creatinine clearance proved to be a significant covariate on clearance (p < 0.001), reducing IIV from 68 to 49%.

CONCLUSION

A population pharmacokinetic model was developed to describe cefotaxime pharmacokinetics and associated variability in adult ICU patients. The estimated creatinine clearance partly explained the IIV in cefotaxime clearance. However, the remaining unexplained IIV is high and suggests a need for dose individualisation using therapeutic drug monitoring where the developed model, after evaluation of predictive performance, may provide support.

A Model-Based Framework to Inform the Dose Selection and Study Design of Emicizumab for Pediatric Patients With Hemophilia A

Emicizumab is a bispecific antibody mimicking the cofactor function of activated coagulation factor VIII to prevent bleeds in patients with hemophilia A. The dose selection for the first-in-child phase III study of emicizumab was addressed by pediatric pharmacokinetic prediction using an adult/adolescent population pharmacokinetic model developed in phase I-I/II studies. The model was modified to incorporate functions describing the age-dependent increase in body weight (BW) with or without clearance maturation to account for the differences in emicizumab pharmacokinetics between adults/adolescents and children. A minimal dose anticipated to achieve in children the same target efficacious exposure as for adults/adolescents was identified when considering BW and clearance maturation. It was the same BW-based dose as for adults/adolescents and was selected for the starting dose for the pediatric study. Whether considering clearance maturation or not in addition to BW led to uncertainty in the pediatric pharmacokinetic prediction and dose selection, which informed implementation of a dose-adapting scheme in the study design. Exposure matching to adults/adolescents was ultimately achieved in children with the starting dose, indicating that consideration of clearance maturation in addition to BW provided adequate pediatric pharmacokinetic predictions for emicizumab. This pharmacokinetic finding in conjunction with exposure-response information served as a basis for the efficacy demonstrated in children, avoiding a time-consuming process for exploring an optimal pediatric dose of emicizumab. This experience indicates that a model-based framework helped optimize the pediatric dose selection and study design, thereby streamlining the development process with extrapolation, of emicizumab for children.

Pharmacokinetics and safety of high-dose rifampicin in children with TB: the Opti-Rif trial

BACKGROUND

Rifampicin doses of 40 mg/kg in adults are safe and well tolerated, may shorten anti-TB treatment and improve outcomes, but have not been evaluated in children.

OBJECTIVES

To characterize the pharmacokinetics and safety of high rifampicin doses in children with drug-susceptible TB.

PATIENTS AND METHODS

The Opti-Rif trial enrolled dosing cohorts of 20 children aged 0-12 years, with incremental dose escalation with each subsequent cohort, until achievement of target exposures or safety concerns. Cohort 1 opened with a rifampicin dose of 15 mg/kg for 14 days, with a single higher dose (35 mg/kg) on day 15. Pharmacokinetic data from days 14 and 15 were analysed using population modelling and safety data reviewed. Incrementally increased rifampicin doses for the next cohort (days 1-14 and day 15) were simulated from the updated model, up to the dose expected to achieve the target exposure [235 mg/L·h, the geometric mean area under the concentration-time curve from 0 to 24 h (AUC0-24) among adults receiving a 35 mg/kg dose].

RESULTS

Sixty-two children were enrolled in three cohorts. The median age overall was 2.1 years (range = 0.4-11.7). Evaluated doses were ∼35 mg/kg (days 1-14) and ∼50 mg/kg (day 15) for cohort 2 and ∼60 mg/kg (days 1-14) and ∼75 mg/kg (day 15) for cohort 3. Approximately half of participants had an adverse event related to study rifampicin; none was grade 3 or higher. A 65-70 mg/kg rifampicin dose was needed in children to reach the target exposure.

CONCLUSIONS

High rifampicin doses in children achieved target exposures and the doses evaluated were safe over 2 weeks.

Maxacalcitol Pharmacokinetic-Pharmacodynamic Modeling and Simulation for Secondary Hyperparathyroidism in Patients Receiving Maintenance Hemodialysis

BACKGROUND

Maxacalcitol was approved in Taiwan in 2018 as the first active vitamin D3 injection for secondary hyperparathyroidism (SHPT) in patients on maintenance hemodialysis. However, no data from any clinical study with maxacalcitol in Taiwanese patients is available.

OBJECTIVES

This analysis aimed to evaluate the profiles of parathyroid hormone (PTH) and calcium (Ca) concentrations in Taiwanese SHPT patients on hemodialysis and maxacalcitol.

METHODS

We developed population pharmacokinetic (PK) and pharmacodynamic (PD) models using a modeling and simulation approach. The data for these analyses were obtained from two studies: a clinical pharmacology study in Japanese patients and an ethnic comparison study in healthy Japanese and -Taiwanese volunteers. We then conducted a simulation study with a PK-PD model comprising the PK and PD models developed here.

RESULTS

Serum maxacalcitol concentration profile was modeled using a two-compartment model that took into consideration the distribution of concentrations below the lower limit of quantification. An ethnic difference in clearance was included in the PK model as a covariate. A PD model that used a PTH/Ca feedback loop best described the observed data. There were no significant differences in Ca or PTH concentrations between Taiwanese and Japanese based on the simulation results from our PK-PD model, even though maxacalcitol exposure was approximately 40% higher in Taiwanese than in Japanese.

CONCLUSIONS

On the basis of these population PK and PD analyses and the clinical study conducted in Japan, there is no clinically relevant difference between Taiwanese and Japanese in terms of serum Ca or PTH levels.

Similar Piperacillin/Tazobactam Target Attainment in Obese versus Nonobese Patients despite Differences in Interstitial Tissue Fluid Pharmacokinetics

Precision dosing of piperacillin/tazobactam in obese patients is compromised by sparse information on target-site exposure. We aimed to evaluate the appropriateness of current and alternative piperacillin/tazobactam dosages in obese and nonobese patients. Based on a prospective, controlled clinical trial in 30 surgery patients (15 obese/15 nonobese; 0.5-h infusion of 4 g/0.5 g piperacillin/tazobactam), piperacillin pharmacokinetics were characterized in plasma and at target-site (interstitial fluid of subcutaneous adipose tissue) via population analysis. Thereafter, multiple 3-4-times daily piperacillin/tazobactam short-term/prolonged (recommended by EUCAST) and continuous infusions were evaluated by simulation. Adequacy of therapy was assessed by probability of pharmacokinetic/pharmacodynamic target-attainment (PTA ≥ 90%) based on time unbound piperacillin concentrations exceed the minimum inhibitory concentration (MIC) during 24 h (%fT>MIC). Lower piperacillin target-site maximum concentrations in obese versus nonobese patients were explained by the impact of lean (approximately two thirds) and fat body mass (approximately one third) on volume of distribution. Simulated steady-state concentrations were 1.43-times, 95%CI = (1.27; 1.61), higher in plasma versus target-site, supporting targets of %fT>2×MIC instead of %fT>4×MIC during continuous infusion to avoid target-site concentrations constantly below MIC. In all obesity and renally impairment/hyperfiltration stages, at MIC = 16 mg/L, adequate PTA required prolonged (thrice-daily 4 g/0.5 g over 3.0 h at %fT>MIC = 50) or continuous infusions (24 g/3 g over 24 h following loading dose at %fT>MIC = 98) of piperacillin/tazobactam.

Precision dosing of intravenous busulfan in pediatric hematopoietic stem cell transplantation: Results from a multicenter population pharmacokinetic study

Busulfan (Bu) is a common component of conditioning regimens before hematopoietic stem cell transplantation (HSCT) and is known for high interpatient pharmacokinetic (PK) variability. This study aimed to develop and externally validate a multicentric, population PK (PopPK) model for intravenous Bu in pediatric patients before HSCT to first study the influence of glutathione‐s‐transferase A1 (GSTA1) polymorphisms on Bu's PK in a large multicentric pediatric population while accounting for fludarabine (Flu) coadministration and, second, to establish an individualized, model‐based, first‐dose recommendation for intravenous Bu that can be widely used in pediatric patients. The model was built using data from 302 patients from five transplantation centers who received a Bu‐based conditioning regimen. External model validation used data from 100 patients. The relationship between body weight and Bu clearance (CL) was best described by an age‐dependent allometric scaling of a body weight model. A stepwise covariate analysis identified Day 1 of Bu conditioning, GSTA1 metabolic groups based on GSTA1 polymorphisms, and Flu coadministration as significant covariates influencing Bu CL. The final model adequately predicted Bu first‐dose CL in the external cohort, with 81% of predicted area under the curves within the therapeutic window. The final model showed minimal bias (mean prediction error, −0.5%; 95% confidence interval [CI], −3.1% to 2.0%) and acceptable precision (mean absolute prediction error percentage, 18.7%; 95% CI, 17.0%–20.5%) in Bu CL prediction for dosing. This multicentric PopPK study confirmed the influence of GSTA1 polymorphisms and Flu coadministration on Bu CL. The developed model accurately predicted Bu CL and first doses in an external cohort of pediatric patients.

Pharmacokinetics and Target Attainment of SQ109 in Plasma and Human-Like Tuberculosis Lesions in Rabbits

SQ109 is a novel well-tolerated drug candidate in clinical development for the treatment of drug-resistant tuberculosis (TB). It is the only inhibitor of the MmpL3 mycolic acid transporter in clinical development. No SQ109-resistant mutant has been directly isolated thus far in vitro, in mice, or in patients, which is tentatively attributed to its multiple targets. It is considered a potential replacement for poorly tolerated components of multidrug-resistant TB regimens. To prioritize SQ109-containing combinations with the best potential for cure and treatment shortening, one must understand its contribution against different bacterial populations in pulmonary lesions. Here, we have characterized the pharmacokinetics of SQ109 in the rabbit model of active TB and its penetration at the sites of disease—lung tissue, cellular and necrotic lesions, and caseum. A two-compartment model with first-order absorption and elimination described the plasma pharmacokinetics. At the human-equivalent dose, parameter estimates fell within the ranges published for preclinical species. Tissue concentrations were modeled using an “effect” compartment, showing high accumulation in lung and cellular lesion areas with penetration coefficients in excess of 1,000 and lower passive diffusion in caseum after 7 daily doses. These results, together with the hydrophobic nature and high nonspecific caseum binding of SQ109, suggest that multiweek dosing would be required to reach steady state in caseum and poorly vascularized compartments, similar to bedaquiline. Linking lesion pharmacokinetics to SQ109 potency in assays against replicating, nonreplicating, and intracellular M. tuberculosis showed SQ109 concentrations markedly above pharmacokinetic-pharmacodynamic targets in lung and cellular lesions throughout the dosing interval.

Population Pharmacokinetics of Meropenem and Vaborbactam Based on Data from Noninfected Subjects and Infected Patients

Meropenem-vaborbactam is a broad-spectrum carbapenem-beta-lactamase inhibitor combination approved in the United States and Europe to treat patients with complicated urinary tract infections and in Europe for other serious bacterial infections, including hospital-acquired and ventilator-associated pneumonia. Population pharmacokinetic (PK) models were developed to characterize the time course of meropenem and vaborbactam using pooled data from two phase 1 and two phase 3 studies. Multicompartment disposition model structures with linear elimination processes were fit to the data using NONMEM 7.2. Since both drugs are cleared primarily by the kidneys, estimated glomerular filtration rate (eGFR) was evaluated as part of the base structural models. For both agents, a two-compartment model with zero-order input and first-order elimination best described the pharmacokinetic PK data, and a sigmoidal Hill-type equation best described the relationship between renal clearance and eGFR. For meropenem, the following significant covariate relationships were identified: clearance (CL) decreased with increasing age, CL was systematically different in subjects with end-stage renal disease, and all PK parameters increased with increasing weight. For vaborbactam, the following significant covariate relationships were identified: CL increased with increasing height, volume of the central compartment (Vc) increased with increasing body surface area, and CL, Vc, and volume of the peripheral compartment were systematically different between phase 1 noninfected subjects and phase 3 infected patients. Visual predictive checks demonstrated minimal bias, supporting the robustness of the final models. These models were useful for generating individual PK exposures for pharmacokinetic-pharmacodynamic (PK-PD) analyses for efficacy and Monte Carlo simulations to evaluate PK-PD target attainment.

A Semi-Mechanistic Model of the Bactericidal Activity of High-Dose Isoniazid Against Multi-Drug-Resistant Tuberculosis: Results from a Randomized Clinical Trial

RATIONALE

There is accumulating evidence that higher-than-standard doses of isoniazid are effective against low-to-intermediate-level isoniazid-resistant strains of Mycobacterium tuberculosis, but the optimal dose remains unknown.

OBJECTIVE

Characterizing the association between isoniazid pharmacokinetics (standard or high-dose) and early bactericidal activity against M. tuberculosis (drug-sensitive and inhA-mutated) and N-acetyltransferase 2 status.

METHODS

ACTG A5312/INHindsight is 7-day early bactericidal activity study with isoniazid at normal dose (5 mg/kg) for patients with drug-sensitive bacteria and 5, 10, and 15 mg/kg doses for patients with inhA mutants. Participants with pulmonary TB received daily isoniazid monotherapy and collected sputum daily. Colony-forming units (CFU) on solid culture and time-to-positivity (TTP) in liquid culture were jointly analyzed using nonlinear mixed-effects modeling.

RESULTS

Fifty-nine adults were included in this analysis. Decline in sputum CFU was described by a one-compartment model, while an exponential bacterial growth model was used to interpret TTP data. The model found bacterial kill is modulated by isoniazid concentration using an effect compartment and a sigmoidal Emax relationship. The model predicted lower potency but similar maximum-kill of isoniazid against inhA-mutated isolates compared to drug-sensitive. Based on simulations from the PK/PD model, to achieve a drop in bacterial load comparable to 5mg/kg against drug-sensitive TB, 10- and 15-mg/kg doses are necessary against inhA-mutated isolates in slow and intermediate N-acetyltransferase 2 acetylators, respectively. Fast acetylators underperformed even at 15 mg/kg.

CONCLUSIONS

Dosing of isoniazid based on N-acetyltransferase 2 acetylator status may help patients attain effective exposures against inhA-mutated isolates while mitigating toxicity risks associated with higher doses. Clinical trial registration available at www.clinicaltrials.gov, ID: NCT01936831.

Impact of a vancomycin loading dose on the achievement of target vancomycin exposure in the first 24 h and on the accompanying risk of nephrotoxicity in critically ill patients

Background

The advocated pharmacokinetic/pharmacodynamic (PK/PD) target for vancomycin, AUC/MIC ≥ 400 mg·h/L, may not be reached with a conventional fixed starting dose of 1000 mg in critically ill patients, but increasing the dose may cause nephrotoxicity.

Objectives

To evaluate the effect of a weight-based loading dose of 25 mg/kg vancomycin on PK/PD target attainment in the first 24 h (AUC_0–24) in critically ill patients and to evaluate whether this increases the risk of acute kidney injury (AKI).

Patients and methods

A prospective observational before/after study was performed in ICU patients, comparing the percentage of vancomycin courses with AUC_0–24 ≥ 400 mg·h/L and the incidence of AKI, defined as worsening of the risk, injury, failure, loss of kidney function and end-stage kidney disease (RIFLE) score. The conventional dose group received 1000 mg of vancomycin as initial dose; the loading dose group received a weight-based loading dose of 25 mg/kg. A population PK model developed using non-linear mixed-effects modelling was used to estimate AUC_0–24 in all patients.

Results

One hundred and four courses from 82 patients were included. With a loading dose, the percentage of courses achieving AUC_0–24 ≥ 400 mg·h/L increased significantly from 53.8% to 88.0% (P = 0.0006). The percentage of patients with new-onset AKI was not significantly higher when receiving a 25 mg/kg loading dose (28.6% versus 37.8%; P = 0.48). However, the risk of AKI was significantly higher in patients achieving AUC_0–24 > 400 mg·h/L compared with patients achieving AUC < 400 mg·h/L (39.0% versus 14.8%; P = 0.031).

Conclusions

A weight-based loading dose of 25 mg/kg vancomycin led to significantly more patients achieving AUC_0–24 ≥ 400 mg·h/L without increased risk of AKI. However, some harm cannot be ruled out since higher exposure was associated with increased risk of AKI.

Pharmacokinetics of Intravenous Paracetamol (Acetaminophen) and Ductus Arteriosus Closure After Premature Birth

In preterm infants, a high risk of hemodynamically significant patent ductus arteriosus (PDA) exists and its persistence is associated with an increased risk of severe morbidity. Current pharmacological options include ibuprofen or indomethacin. However, treatment by indomethacin or ibuprofen of a large PDA was shown to reduce early pulmonary hemorrhage and later medical treatment but had no effect on neonatal death or morbidity. Early prophylactic treatment of ductus arteriosus by paracetamol seems to be an attractive opportunity to reduce life-threatening morbidity. However, there are currently no data regarding the pharmacokinetics (PK) and pharmacodynamics of paracetamol in preterm neonates in this potential new indication. In this study, we aimed to develop a population PK model for paracetamol and investigate the relationship between paracetamol exposure levels and time to contraction of the ductus. Data were modeled using Monolix software. A one-compartment model adequately described the paracetamol concentration-time course. A Weibull model adequately described the time to contraction of the ductus. Our results suggest that the dosage used in this study (i.e., first day 42.5 mg/kg, then 30 mg/kg/day) allows for reaching the maximum inhibition response from paracetamol regarding the time to close the ductus. However, this study pointed out a lower effect of paracetamol on extremely preterm neonates (below 27 weeks). Therefore, a dose-finding study focusing specifically on extremely preterm neonates with treatment efficacy and toxicity is strongly needed.

The effect of isoniazid intake on ethionamide pharmacokinetics and target attainment in multidrug-resistant tuberculosis patients

Ethionamide is recommended as part of regimens to treat multidrug-resistant and rifampicin-resistant tuberculosis. The study was conducted to (i) describe the distribution of ethionamide minimum inhibitory concentrations (MICs), (ii) describe the pharmacokinetics of ethionamide, and (iii) determine the probability of attaining target AUC_0-24/MIC values associated with suppression of resistant subpopulation and microbial kill. Participants received 15-20 mg/kg of ethionamide daily (in 500 or 750 mg doses), as part of a multidrug regimen. Pretreatment MICs of ethionamide for M. tuberculosis sputum isolates were determined using Sensititre MYCOTB MIC plates. Plasma concentrations of ethionamide (measured pre-dose and at 2, 4, 6, 8 and 10 hours post-dose) were available for 84 patients. A one-compartment disposition model including a liver compartment capturing hepatic extraction, best described ethionamide pharmacokinetics. Clearance and volume were allometrically scaled using fat-free mass. Isoniazid co-administration reduced ethionamide clearance by 31% resulting in a 44% increase in AUC_0-24. The median (range) MIC (n=111) was 2.5 mg/L (<0.3 to >40 mg/L). Simulations showed increased daily doses of ethionamide (1 250 mg, 1 500 mg, and 1 750 mg for patients weighing ≤45 kg, 46-70 kg, and >70 kg, respectively) resulted in the probability of attaining a fAUC_0-24/MIC ratio ≥ 42 in more than 90% of patients, only at the lowest MIC of 0.3 mg/L. The WHO recommended doses of ethionamide do not achieve target concentrations even for the lowest MIC measured in the cohort.

Population Pharmacokinetics of Ipatasertib and Its Metabolite in Cancer Patients

Ipatasertib is a selective AKT kinase inhibitor currently in development for the treatment of several solid tumors, including breast and prostate cancers. This study was undertaken to characterize pharmacokinetic profiles of ipatasertib and its metabolite M1 (G-037720) and to understand the sources of variability. Population pharmacokinetic models of ipatasertib and M1 were developed separately using data from 342 individuals with cancer from 5 phase 1 and 2 studies. The final population pharmacokinetic models for ipatasertib and M1 were 3-compartmental, with first-order elimination and sequential zero- and first-order absorption. Ipatasertib bioavailability and M1 formation increased after multiple dosing, resulting in an increase in exposure beyond that expected from accumulation alone. Covariate effects of ipatasertib include decreased oral clearance with increasing age and with coadministration of abiraterone, as well as decreased bioavailability with increasing weight. For ages 37 and 80 years, steady-state area under the curve (AUC_ss ) was predicted to be 81% and 109%, respectively, of the typical population value (64 years). For body weight of 49 and 111 kg, AUC_ss was predicted to be 132% and 78%, respectively, of the typical population value (75 kg). The small magnitude of change in ipatasertib exposure is not likely to be clinically relevant. For M1, the peripheral distribution volume and intercompartmental clearance increased with increasing weight. Coadministration of abiraterone was estimated to increase M1 exposure by 61% at steady state. Mild and moderate renal impairment, mild hepatic impairment, and race were not identified as significant covariates in the final models for ipatasertib and M1.

Constructing a representative in-silico population for paediatric simulations: Application to HIV-positive African children

AIMS

Simulations are an essential tool for investigating scenarios in pharmacokinetics-pharmacodynamics. The models used during simulation often include the effect of highly correlated covariates such as weight, height and sex, and for children also age, which complicates the construction of an in silico population. For this reason, a suitable and representative patient population is crucial for the simulations to produce meaningful results. For simulation in paediatric patients, international growth charts from the World Health Organization (WHO) and the Centers for Disease Control and Prevention (CDC) provide a reference, but these may not always be representative for specific populations, such as malnourished children with HIV or acutely unwell children.

METHODS

We present a workflow to construct a virtual paediatric patient population using WHO and CDC growth charts, suggest piecewise linear functions to adjust the median of the growth charts by sex and age, and suggest visual diagnostics to compare with the target population. We applied this workflow in a population of 1206 HIV-positive African children, consisting of 19 742 observations with weight ranging from 3.8 to 79.7 kg, height from 55.5 to 180 cm, and an age between 0.40 and 18 years.

RESULTS

Before adjustment, the WHO and CDC charts produced weights and heights higher compared to the observed data. After applying our methodology, we could simulate weight, height, sex and age combinations in good agreement with the observed data.

CONCLUSION

The methodology presented here is flexible and may be applied to other scenarios where WHO and CDC growth standards might not be appropriate. In addition we provide R scripts and a large ready-to-use paediatric population.

Analysis of the side effect of QTc interval prolongation in the bedaquiline regimen in drug resistant tuberculosis patients

OBJECTIVES

Indonesia is one of the top 20 countries with the highest prevalence of drug resistant tuberculosis (DR-TB) worldwide with a percentage of new cases of 2.4% and retreatment of 13%. Bedaquiline (BDQ) is one of the drugs that used in the individual long regimen treating DR-TB. BDQ is also combined with levofloxacin (LFX) and/or clofazimine (CFZ) that can cause QTc interval prolongation. The aim was to study the differences in the use of BDQ regimens to the lengthening of the QTc interval and to study risk factors (diabetes, hypokalemia, sex, BMI, and age) in BDQ regimen.

METHODS

This study was an observational retrospective study with a total sampling method, which was conducted at Dr. Soetomo General Hospital Surabaya. Samples from this study were patients diagnosed with DR-TB at Dr. Soetomo General Hospital Surabaya in the period of January 2015-December 2019 who used BDQ regimen and met the inclusion criteria. The ECG data were analyzed from the mean of each group (BDQ regimen and risk factors), also analyzed using statistical analysis.

RESULTS

Data obtained from total sample in this study were 73 patients. The most widely used different regimens in this study were the combination of BDQ + LFX by 36 patients (49.3%), BDQ + LFX + CFZ by 16 patients (21.9%), BDQ by 11 patients (15.1%) and BDQ + CFZ 10 patients (13.7%). Out of 73 patients, 52 patients (71.2%) experienced lengthening of the QT interval and grade 1 of QTc interval prolongation occurred in most patients and also the onset was mostly one month after using BDQ regimen. The side effects of QTc interval prolongation from groups of combination and risk factors were no difference in each month (p>0.05).

CONCLUSIONS

This study can be concluded that there were no differences in the QTc prolongation between the groups of BDQ regimen (BDQ, BDQ + LFX, BDQ + CFZ and BDQ + LFX + CFZ) and the groups of risk factors.

A Population Pharmacokinetic Model of Macitentan and Its Active Metabolite Aprocitentan in Healthy Volunteers and Patients with Pulmonary Arterial Hypertension

BACKGROUND

Macitentan and its active metabolite, aprocitentan, are non-peptide, potent, dual endothelin receptor antagonists. Macitentan is approved for the treatment of pulmonary arterial hypertension in adults, at a dose of 10 mg/day.

OBJECTIVE

The objective of this study was to develop a comprehensive population model to describe the pharmacokinetics of macitentan and aprocitentan in healthy adults and adult subjects with pulmonary arterial hypertension.

METHODS

Pharmacokinetic data of 452 subjects in nine studies, after single and repeated doses (dose range 0.2-600 mg), were pooled for a non-linear mixed-effects analysis and the assessment of covariates, i.e., body weight, age, sex, race, renal and hepatic impairment, health status (healthy volunteers vs patients with pulmonary arterial hypertension), and formulation (capsules vs tablets) on pharmacokinetic parameters.

RESULTS

The final model was an open one-compartment disposition model, with linear elimination for macitentan and linear formation and elimination for aprocitentan. A semi-mechanistic absorption model described the dose dependency and multiple peaks observed for macitentan. For a female patient with pulmonary arterial hypertension after oral administration at 10 mg, macitentan reached a maximum concentration after 9 h and, following daily dosing, reached steady state after 3 days with a twofold accumulation factor. The apparent volume of distribution was 34 L and clearance was 1.39 L/h. Aprocitentan reached maximum concentration after 51 h and steady state after 9 days, with a 12.5-fold accumulation factor. Body weight, sex, race, renal impairment, health status, and formulation were statistically significant covariates on pharmacokinetic parameters.

CONCLUSIONS

The comprehensive population pharmacokinetic model adequately described the pharmacokinetics of macitentan and aprocitentan across different dose concentrations, regimens, and formulations. Several covariates significantly influenced the pharmacokinetics of macitentan and aprocitentan, but none was considered clinically relevant.

Population pharmacokinetics of intravenous and oral ciprofloxacin in children to optimize dosing regimens

PURPOSE

This study aimed to characterize pharmacokinetics of intravenous and oral ciprofloxacin in children to optimize dosing scheme.

METHODS

Children treated with ciprofloxacin were included. Pharmacokinetics were described using non-linear mixed-effect modelling and validated with an external dataset. Monte Carlo simulations investigated dosing regimens to achieve a target AUC_0-24 h/MIC ratio ≥ 125.

RESULTS

A total of 189 children (492 concentrations) were included. A two-compartment model with first-order absorption and elimination best described the data. An allometric model was used to describe bodyweight (BW) influence, and effects of estimated glomerular filtration rate (eGFR) and age were significant on ciprofloxacin clearance.

CONCLUSION

The recommended IV dose of 10 mg/kg q8h, not exceeding 400 mg q8h, would achieve AUC_0-24 h to successfully treat bacteria with MICs ≤ 0.25 (e.g. Salmonella, Escherichia coli, Proteus, Haemophilus, Enterobacter, and Klebsiella). A dose increase to 600 mg q8h in children > 40 kg and to 15 mg/kg q8h (max 400 mg q8h, max 600 mg q8h if augmented renal clearance, i.e., eGFR > 200 mL/min/1.73 m²) in children < 40 kg would be needed for the strains with highest MIC (16% of Pseudomonas aeruginosa and 47% of Staphylococcus aureus). The oral recommended dose of 20 mg/kg q12h (not exceeding 750 mg) would cover bacteria with MICs ≤ 0.125 but may be insufficient for bacteria with higher MIC and a dose increase according bodyweight and eGFR would be needed. These doses should be prospectively confirmed, and a therapeutic drug monitoring could be used to refine them individually.

Isoniazid Population Pharmacokinetics and Dose Recommendation for Korean Patients With Tuberculosis Based on Target Attainment Analysis

The wide variability of isoniazid (INH) pharmacokinetics is mainly attributed to the trimodal N-acetyltransferase 2 (NAT2) acetylator phenotype, that is, rapid, intermediate, and slow. Consequently, a uniform INH dose in current clinical practice may lead to treatment failure and emergence of drug resistance. There is a lack of studies on specific doses of INH for different NAT2 acetylator phenotypes among tuberculosis (TB) patients. Therefore, we aimed to provide insight into the optimal dosing of INH for each NAT2 acetylator phenotype with respect to the probability of achieving a pharmacokinetic (PK)/pharmacodynamic target. PK, the NAT2 genotype, and clinical data were collected in a multicenter prospective cohort study conducted at 13 clinical centers in Korea. Population PK modeling and simulation were carried out. Data from 454 TB patients were divided into a training data set and a test data set at a ratio of 4 to 1. The PK of the training data were best described by a 2-compartment model with allometric scaling for body size effect. Importantly, NAT2 acetylator phenotypes significantly affected the apparent clearance. Our model, which provided better predictive performance compared with previously published models, was evaluated by external validation using the test set. The simulation for assessing target efficacy and toxicity indicated that the best INH dosing regimens for Korean tuberculosis patients were once-daily doses of 400, 300, and 200 mg for rapid, intermediate, and slow acetylators, respectively. In conclusion, our study provides a step forward in precision dosing for antituberculosis management.

A Semimechanistic Pharmacokinetic Model for Depot Medroxyprogesterone Acetate and Drug-Drug Interactions With Antiretroviral and Antituberculosis Treatment

Depot medroxyprogesterone acetate is an injectable hormonal contraceptive, widely used by women of childbearing potential living with HIV and/or tuberculosis. As medroxyprogesterone acetate is a cytochrome P450 (CYP3A4) substrate, drug-drug interactions (DDIs) with antiretroviral or antituberculosis treatment may lead to subtherapeutic medroxyprogesterone acetate concentrations (< 0.1 ng/mL), resulting in contraception failure, when depot medroxyprogesterone is dosed at 12-week intervals. A pooled population pharmacokinetic analysis with 744 plasma medroxyprogesterone acetate concentrations from 138 women treated with depot medroxyprogesterone and antiretroviral/antituberculosis treatment across three clinical trials was performed. Monte Carlo simulations were performed to predict the percentage of participants with subtherapeutic medroxyprogesterone acetate concentrations and to derive alternative dosing strategies. Medroxyprogesterone acetate clearance increased by 24.7% with efavirenz coadministration. Efavirenz plus antituberculosis treatment (rifampicin + isoniazid) increased clearance by 52.4%. Conversely, lopinavir/ritonavir and nelfinavir decreased clearance (28.7% and 15.8%, respectively), but lopinavir/ritonavir also accelerated medroxyprogesterone acetate's appearance into the systemic circulation, thus shortening the terminal half-life. A higher risk of subtherapeutic medroxyprogesterone acetate concentrations at Week 12 was predicted on a typical 60-kg woman on efavirenz (4.99%) and efavirenz with antituberculosis treatment (6.08%) when compared with medroxyprogesterone acetate alone (2.91%). This risk increased in women with higher body weight. Simulations show that re-dosing every 8 to 10 weeks circumvents the risk of subtherapeutic medroxyprogesterone acetate exposure associated with these DDIs. Dosing depot medroxyprogesterone every 8 to 10 weeks should eliminate the risk of subtherapeutic medroxyprogesterone acetate exposure caused by coadministered efavirenz and/or antituberculosis treatment, thus reducing the risk of contraceptive failure.

Prospective Cohort Study of Micafungin Population Pharmacokinetic Analysis in Plasma and Peritoneal Fluid in Septic Patients with Intra-abdominal Infections

The objective of this study was to describe the pharmacokinetics (PK) of micafungin in plasma and peritoneal fluid in septic patients with intra-abdominal infections. Twelve patients with secondary peritonitis in septic shock receiving 100 mg micafungin once daily were included. Total micafungin plasma and peritoneal fluid were subjected to a population pharmacokinetic analysis using Pmetrics. Monte Carlo simulations were performed considering the total area under the curve from 0 to 24 h (AUC_0-24)/MIC ratios in plasma. Micafungin concentrations in both plasma and the peritoneal exudate were best described by a three-compartmental PK model with the fat-free mass (FFM) as a covariate of clearance (CL) and the volume of the central compartment (V_c). The mean parameter estimates (standard deviations [SD]) were 1.18 (0.40) liters/h for CL and 12.85 (4.78) liters for V_c. The mean peritoneal exudate/plasma ratios (SD) of micafungin were 25% (5%) on day 1 and 40% (8%) between days 3 and 5. Dosing simulations supported the use of standard 100-mg daily dosing for Candida albicans (FFM, <60 kg), C. glabrata (FFM, <50 kg), and C. tropicalis (FFM, <30 kg) on the second day of therapy. There is a moderate penetration of micafungin into the peritoneal cavity (25 to 40%). For empirical treatment, a dose escalation of at least a loading dose of 150 mg depending on the FFM of patients and the Candida species is suggested to be effective from the first day of therapy.

Geographical Variability in Paromomycin Pharmacokinetics Does Not Explain Efficacy Differences between Eastern African and Indian Visceral Leishmaniasis Patients

Introduction

Intramuscular paromomycin monotherapy to treat visceral leishmaniasis (VL) has been shown to be effective for Indian patients, while a similar regimen resulted in lower efficacy in Eastern Africa, which could be related to differences in paromomycin pharmacokinetics.

Methods

Pharmacokinetic data were available from two randomized controlled trials in VL patients from Eastern Africa and India. African patients received intramuscular paromomycin monotherapy (20 mg/kg for 21 days) or combination therapy (15 mg/kg for 17 days) with sodium stibogluconate. Indian patients received paromomycin monotherapy (15 mg/kg for 21 days). A population pharmacokinetic model was developed for paromomycin in Eastern African and Indian VL patients.

Results

Seventy-four African patients (388 observations) and 528 Indian patients (1321 observations) were included in this pharmacokinetic analysis. A one-compartment model with first-order kinetics of absorption and elimination best described paromomycin in plasma. Bioavailability (relative standard error) was 1.17 (5.18%) times higher in Kenyan and Sudanese patients, and 2.46 (24.5%) times higher in Ethiopian patients, compared with Indian patients. Ethiopian patients had an approximately fourfold slower absorption rate constant of 0.446 h^–1 (18.2%). Area under the plasma concentration-time curve for 24 h at steady-state (AUC_τ,SS) for 15 mg/kg/day (median [interquartile range]) was higher in Kenya and Sudan (172.7 µg·h/mL [145.9–214.3]) and Ethiopia (230.1 µg·h/mL [146.3–591.2]) compared with India (97.26 µg·h/mL [80.83–123.4]).

Conclusion

The developed model provides detailed insight into the pharmacokinetic differences among Eastern African countries and India, however the resulting differences in paromomycin exposure do not seem to explain the geographical differences in paromomycin efficacy in the treatment of VL patients.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-01036-8.

Model-Informed Optimization of a Pediatric Clinical Pharmacokinetic Trial of a New Spironolactone Liquid Formulation

Quantitative pharmacology brings important advantages in the design and conduct of pediatric clinical trials. Herein, we demonstrate the application of a model-based approach to select doses and pharmacokinetic sampling scenarios for the clinical evaluation of a novel oral suspension of spironolactone in pediatric patients with edema. A population pharmacokinetic model was developed and qualified for spironolactone and its metabolite, canrenone, using data from adults and bridged to pediatrics (2 to <17 years old) using allometric scaling. The model was then used via simulation to explore different dosing and sampling scenarios. Doses of 0.5 and 1.5 mg/kg led to target exposures (i.e., similar to 25 and 100 mg of the reference product in adults) in all the reference pediatric ages (i.e., 2, 6, 12 and 17 years). Additionally, two different sampling scenarios were delineated to accommodate patients into sparse sampling schemes informative to characterize drug pharmacokinetics while minimizing phlebotomy and burden to participating children.

Translational precision medicine: an industry perspective

In the era of precision medicine, digital technologies and artificial intelligence, drug discovery and development face unprecedented opportunities for product and business model innovation, fundamentally changing the traditional approach of how drugs are discovered, developed and marketed. Critical to this transformation is the adoption of new technologies in the drug development process, catalyzing the transition from serendipity-driven to data-driven medicine. This paradigm shift comes with a need for both translation and precision, leading to a modern Translational Precision Medicine approach to drug discovery and development. Key components of Translational Precision Medicine are multi-omics profiling, digital biomarkers, model-based data integration, artificial intelligence, biomarker-guided trial designs and patient-centric companion diagnostics. In this review, we summarize and critically discuss the potential and challenges of Translational Precision Medicine from a cross-industry perspective.

Test-dose pharmacokinetics guided melphalan dose adjustment in reduced intensity conditioning allogeneic transplant for non-malignant disorders

AIMS

We studied melphalan pharmacokinetics (PK) and feasibility of melphalan full-dose adjustment based on test-dose PK in children and young adults with non-malignant disorders (NMD) undergoing allogeneic hematopoietic cell transplantation (HCT) using reduced intensity conditioning (RIC) containing alemtuzumab, fludarabine and melphalan.

METHODS

Patients received test-dose melphalan (10% of planned full-dose) prior to conditioning. Blood samples for PK were obtained around test and full-dose melphalan (140 mg/m² or 4.7 mg/kg in patients <10 kg). Melphalan concentration was measured by liquid chromatography electrospray ionization tandem mass-spectrometry assay and data were analysed using a population-PK model and Bayesian estimation. Test and full-dose melphalan clearance estimates were evaluated by pairwise Wilcoxon test and Bland-Altman plot.

RESULTS

Twenty-four patients undergoing 25 transplants were included in the final analysis. Patients received standard full-dose melphalan in 17 transplants, with median area under the concentration-time curve (AUC) of 5.5 mg*h/L (range, 3.0-9.5 mg*h/L). Patients received test-dose melphalan in 23 transplants with a test-dose PK predicted full-dose AUC range of 2.9-16.8 mg*h/L. In seven transplants where patients had baseline organ impairment, test-dose PK predicted higher exposure for standard full-dose (median AUC 13.8 mg*h/L). Melphalan full-dose was adjusted in these patients, with achievement of desired target AUC (3.6-5.4 mg*h/L) and no excess toxicity. Mean ratio of test-dose clearance to full-dose clearance was 1.03. Twenty of 22 patients (91%) were within the 95% confidence intervals of the clearance ratio.

CONCLUSION

Melphalan test-dose PK reliably predicts full-dose PK and allows for accurate adjustment of full-dose melphalan in RIC-HCT for NMD. This approach can avoid excess toxicity from increased systemic exposure, especially in patients with organ impairment.

Prediction of levobupivacaine concentrations in neonates and infants following neuraxial rescue blocks

AIM

Pharmacokinetic simulation was used to characterize levobupivacaine disposition after regional anesthetic rescue for failed spinal anesthesia in neonates and infants.

METHODS

Population pharmacokinetics of levobupivacaine were estimated after spinal blockade in a cohort of neonates and infants (n = 25, postnatal age 5-18 weeks, gestation 21-41 weeks, weight 2.4-6 kg). Total levobupivacaine concentrations were assayed 3-4 times in the first hour after spinal levobupivacaine 1 mg kg^-1 administration. Parameters were estimated using nonlinear mixed-effects models and supported by priors. Covariates included postnatal age and total body weight. Parameter estimates were used to simulate total levobupivacaine concentrations after a primary spinal levobupivacaine 1 mg kg^-1 with rescue caudal levobupivacaine 1.5-2.5 mg kg^-1 .

RESULTS

A one-compartment model with a mature clearance 21.5 L h^-1  70 kg^-1 (CV 47.3%) and central volume 189 L 70 kg^-1 (CV 37%) adequately described time-concentration profiles. Clearance maturation was described using a maturation half-time of 11.5 weeks postnatal age. The absorption half-time for spinal levobupivacaine was 2.6 min (CV 56.8%). The upper (97.5% prediction) for peak concentrations after rescue caudal levobupivacaine were 1.5 mg kg^-1 , 2 mg kg^-1 , and 2.5 mg kg^-1 was 2.05 mg L^-1 , 2.5 mg L^-1 , and 2.9 mg L^-1 respectively.

CONCLUSION

Total bupivacaine concentrations greater than 2.5 mg L^-1 are associated with neurotoxicity in adults. Predicted concentrations after either a repeat spinal or a caudal rescue dose of levobupivacaine 1.5 mg kg^-1 1 h after spinal levobupivacaine administration are below the neurotoxic concentration threshold.

A Semi-Mechanistic Population Pharmacokinetic/Pharmacodynamic Model of Bortezomib in Pediatric Patients with Relapsed/Refractory Acute Lymphoblastic Leukemia

INTRODUCTION

The pharmacokinetics (PK) of the 20S proteasome inhibitor bortezomib are characterized by a large volume of distribution and a rapid decline in plasma concentrations within the first hour after administration. An increase in exposure was observed in the second week of treatment, which has previously been explained by extensive binding of bortezomib to proteasome in erythrocytes and peripheral tissues. We characterized the nonlinear population PK and pharmacodynamics (PD) of bortezomib in children with acute lymphoblastic leukemia.

METHODS

Overall, 323 samples from 28 patients were available from a pediatric clinical study investigating bortezomib at an intravenous dose of 1.3 mg/m² twice weekly (Dutch Trial Registry number 1881/ITCC021). A semi-physiological PK model for bortezomib was first developed; the PK were linked to the decrease in 20S proteasome activity in the final PK/PD model.

RESULTS

The plasma PK data were adequately described using a two-compartment model with linear elimination. Increased concentrations were observed in week 2 compared with week 1, which was described using a Langmuir binding model. The decrease in 20S proteasome activity was best described by a direct effect model with a sigmoidal maximal inhibitory effect, representing the relationship between plasma concentrations and effect. The maximal inhibitory effect was 0.696 pmol AMC/s/mg protein (95% confidence interval 0.664-0.728) after administration.

CONCLUSION

The semi-physiological model adequately described the nonlinear PK and PD of bortezomib in plasma. This model can be used to further optimize dosing of bortezomib.

In vivo pharmacokinetic evaluation of carprofen delivery from intra-articular nanoparticles in rabbits: A population modelling approach

Osteoarthritis is treated with COX or fosfolipase A2 inhibitors such as carprofen, a propionic acid NSAID. The enhancement of its action over the articular cartilage is mandatory to facilitate its therapeutic application. Drug uptake into the cartilage requires high synovial fluid concentrations, anticipating its rapid distribution towards bloodstream. Thus, intraarticular administration improves local targeting of the drug, lining with the site of action. A pharmacokinetic study in rabbits has been performed to evaluate carprofen nanoparticles after intraarticular administration. Pharmacokinetic analysis of plasma profiles through a modelling approach, has demonstrated the rapid distribution of drug outside of synovial chamber but mainly remaining in plasma. The data modelling has demonstrated the existence of two release-absorption processes when the nanoparticles are administered in the synovial space. Additionally, results are predictive of the PK profile of some other species such as cat, dogs or humans.

Mycobacterium tuberculosis precursor rRNA as a measure of treatment-shortening activity of drugs and regimens

There is urgent need for new drug regimens that more rapidly cure tuberculosis (TB). Existing TB drugs and regimens vary in treatment-shortening activity, but the molecular basis of these differences is unclear, and no existing assay directly quantifies the ability of a drug or regimen to shorten treatment. Here, we show that drugs historically classified as sterilizing and non-sterilizing have distinct impacts on a fundamental aspect of Mycobacterium tuberculosis physiology: ribosomal RNA (rRNA) synthesis. In culture, in mice, and in human studies, measurement of precursor rRNA reveals that sterilizing drugs and highly effective drug regimens profoundly suppress M. tuberculosis rRNA synthesis, whereas non-sterilizing drugs and weaker regimens do not. The rRNA synthesis ratio provides a readout of drug effect that is orthogonal to traditional measures of bacterial burden. We propose that this metric of drug activity may accelerate the development of shorter TB regimens.

Initial dose recommendation for sirolimus in paediatric kaposiform haemangioendothelioma patients based on population pharmacokinetics and pharmacogenomics

Objective

Sirolimus has been used to treat paediatric kaposiform haemangioendothelioma patients. However, there is considerable pharmacokinetic variability among individuals, and it is difficult to develop an initial dosing regimen. The goal of the present study is to recommend an initial sirolimus dose in paediatric kaposiform haemangioendothelioma patients based on population pharmacokinetics and pharmacogenomics.

Methods

This was a retrospective clinical study. A population pharmacokinetics model was established and population characteristics, laboratory test results, drug combinations, and pharmacogenomics were considered as potential covariates. The Monte Carlo method was used to simulate the optimal initial dosage.

Results

The final covariates that affect sirolimus clearance include weight and the CYP3A5 genotype. The initial dosage of sirolimus for individuals with CYP3A5*3/*3 was 0.20 mg/kg split into two doses for 5 to 60 kg body weight. For individuals with CYP3A5*1, the initial dose was 0.23 mg/kg split into two doses for 5 to 30 kg body weight and 0.20 mg/kg split into two doses for 30 to 60 kg body weight.

Conclusion

The recommendation for the initial sirolimus dose in paediatric kaposiform haemangioendothelioma patients was based on population pharmacokinetics and pharmacogenomics. This study may provide practical value for sirolimus clinical use in paediatric kaposiform haemangioendothelioma patients.

Informing Pharmacokinetic Models With Physiological Data: Oral Population Modeling of L-Serine in Humans

To determine how to set optimal oral L-serine (serine) dose levels for a clinical trial, existing literature was surveyed. Data sufficient to set the dose was inadequate, and so an (n = 10) phase I-A calibration trial was performed, administering serine with and without other oral agents. We analyzed the trial and the literature data using pharmacokinetic (PK) modeling and statistical analysis. The therapeutic goal is to modulate specific serine-related metabolic pathways in the liver using the lowest possible dose which gives the desired effect since the upper bound was expected to be limited by toxicity. A standard PK approach, in which a common model structure was selected using a fit to data, yielded a model with a single central compartment corresponding to plasma, clearance from that compartment, and an endogenous source of serine. To improve conditioning, a parametric structure was changed to estimate ratios (bioavailability over volume, for example). Model fit quality was improved and the uncertainty in estimated parameters was reduced. Because of the particular interest in the fate of serine, the model was used to estimate whether serine is consumed in the gut, absorbed by the liver, or entered the blood in either a free state, or in a protein- or tissue-bound state that is not measured by our assay. The PK model structure was set up to represent relevant physiology, and this quantitative systems biology approach allowed a broader set of physiological data to be used to narrow parameter and prediction confidence intervals, and to better understand the biological meaning of the data. The model results allowed us to determine the optimal human dose for future trials, including a trial design component including IV and tracer studies. A key contribution is that we were able to use human physiological data from the literature to inform the PK model and to set reasonable bounds on parameters, and to improve model conditioning. Leveraging literature data produced a more predictive, useful model.

Extrapolation and dosing recommendations for raxibacumab in children from birth to age <18 years

AIMS

The US Food and Drug Administration's Animal Rule allows for the approval of drugs when human efficacy studies are not ethical. While the therapeutic doses of raxibacumab, a monoclonal antibody for the prophylaxis and treatment of inhalational anthrax, have been based on pharmacokinetic data from adult subjects, its disposition in children has not been investigated in clinical trials. Here we evaluate the effect of demographic covariates and maturation processes on the pharmacokinetics of raxibacumab and explore opportunities for the optimisation of paediatric doses.

METHODS

A population pharmacokinetic model was used as basis for the extrapolation of raxibacumab disposition from adults to children. Different extrapolation scenarios, including weight-banded dosing regimens, were considered to assess the effect of growth and maturation on the pharmacokinetic parameters of interest. Area under the concentration-time curve, maximum plasma concentration and the time of serum raxibacumab concentrations greater than or equimolar to the highest serum protective antigen concentrations observed for at least 28 days in any monkey challenged with Bacillus anthracis that died were derived and compared with the currently approved US doses.

RESULTS

Based on practical considerations, a weight-banded dosing regimen consisting of 4 dose levels (75 mg/kg for individuals ≤1.5 kg, 55 mg/kg for individuals <10 kg, 45 mg/kg for individuals <50 kg, 40 mg/kg for all individuals >50 kg) was required to optimise target exposure across the paediatric population.

CONCLUSIONS

Age-related maturation processes may affect raxibacumab clearance in very young patients. The proposed dosing regimens take into account effects of body weight and maturation processes on the elimination of raxibacumab.

An Integrated Paediatric Population PK/PD Analysis of dDAVP: How do PK Differences Translate to Clinical Outcomes?

INTRODUCTION

The bioequivalence of two formulations of desmopressin (dDAVP), a vasopressin analogue prescribed for nocturnal enuresis treatment in children, has been previously confirmed in adults but not in children. In this study, we aimed to study the pharmacokinetics (PK) and pharmacodynamics (PD) of these two formulations, in both fasted and fed children, including patients younger than 6 years of age.

METHODS

Previously published data from one PK study and one PK/PD study in children aged between 6 and 16 years were combined with a new PK/PD study in children aged between 6 months and 8 years, and analysed using population PK/PD modelling. Simulations were performed to further explore the relative bioavailability of both formulations and evaluate current dosing strategies.

RESULTS

The complex absorption behaviour of the lyophilizate was modelled using a double input, linked to a one-compartmental model with linear elimination and an indirect response model linking dDAVP concentration to produced urine volume and osmolality. The final model described the observed data well and elucidated the complexity of bioequivalence and therapeutic equivalence of the two formulations. Simulations showed that current dosing regimens using a fixed dose of lyophilizate 120 μg is not adequate for children, assuming children to be in the fed state when taking dDAVP. A new age- and weight-based dosing regimen was suggested and was shown to lead to improved, better tailored effects.

CONCLUSIONS

Bioequivalence and therapeutic equivalence data of two formulations of the same drug in adults cannot be readily extrapolated to children. This study shows the importance of well-designed paediatric clinical trials and how they can be analysed using mixed-effects modelling to make clinically relevant inferences. A follow-up clinical trial testing the proposed dDAVP dosing regimen should be performed.

CLINICAL TRIAL REGISTRATION

This trial has been registered at www.clinicaltrials.gov (identifier NCT02584231; EudraCT 2014-005200-13).

Normal fat mass cannot be reliably estimated in typical pharmacokinetic studies

PURPOSE

An influential covariate for pharmacokinetics is (body) size. Recently, the method of estimation of normal fat mass (NFM) has been advocated. Here, the relative contribution of fat mass, estimated as a fraction fat (Ffat), is used to explain differences in pharmacokinetic parameters. This concept is more and more applied. However, it remains unclear whether NFM can be reliably estimated in these typical studies.

METHODS

We performed an evaluation of the reliability of NFM estimation in a typical study size (n = 30), otherwise best-case scenario, by means of a pharmacokinetic simulation study. Several values of Ffat were investigated.

RESULTS

In a typical pharmacokinetic study, high imprecision was observed for NFM parameter estimates over a range of scenarios. For example, in a scenario where the true value of Ffat on clearance was 0.5, we found a 95% confidence interval of - 0.1 to 2.1, demonstrating a low precision. The implications for practice are that one could conclude that fat-free mass best describes the relationship of the pharmacokinetics with body size, while the true relationship was between fat-free mass and total body weight. Consequently, this could lead to incorrect extrapolation of pharmacokinetics to extreme body sizes.

CONCLUSION

In typical pharmacokinetic studies, NFM should be used with caution because the Ffat estimates have low precision. The estimation of Ffat should always be preceded by careful study design evaluation before planning a study, to ensure that the design and sample size is sufficient to apply this potentially useful methodology.