Transcriptome profiles of macrophages upon infection by morphotypic smooth and rough variants of Mycobacterium abscessus.

Mycobacterium abscessus (Mab) infection can be deadly in patients with chronic lung diseases like cystic fibrosis (CF). In vitro and in vivo, Mab may adopt a smooth (S) or rough (R) morphotype, the latter linked to more severe disease conditions. In vitro studies revealed differences in pathogenicity and immune response to S and R morphotypes. We propose that in vivo both morphotypes exist and may transiently switch depending on the environment, having important pathogenic and immunologic consequences. This can be modeled by morphotypic S and R variants of Mab selected based on in vitro growth conditions. Here, we report the first analysis of early transcriptional events in mouse bone marrow derived macrophages (BMDMs) upon infection with media-selected interchangeable Mab-S and Mab-R morphotypes. The early transcriptional events after infection with both morphotypes showed considerable overlap of the pro-inflammatory genes that were differentially regulated compared to the uninfected macrophages. We also observed signature genes significantly differentially regulated in macrophages during infection of media-selected morphotypic Mab-S and Mab-R variants. In conclusion, media-selected Mab-S and Mab-R behave in a similar fashion to stable S and R types with respect to pathogenesis and immune response, serving as a useful model for environmentally influenced morphotype selection.

Time-dependent pharmacodynamics of amikacin on Mycobacterium abscessus growth and resistance emergence.

Mycobacterium abscessus pulmonary disease is increasing in prevalence globally, particularly for individuals with cystic fibrosis. These infections are challenging to treat due to a high rate of resistance. Amikacin is critical to treatment, but the development of toxicity, amikacin resistance, and treatment failure are significant challenges. Amikacin has been characterized previously as peak-dependent and extended-interval dosing is commonly used. In our hollow fiber infection model of M. abscessus, amikacin exhibited time-dependent rather than the expected peak-dependent pharmacodynamics. Humanized amikacin exposures with more frequent, short-interval dosing (continuous infusion or every 12 hours) yielded improved microbiological response compared to extended-interval dosing (every 24 hours or 1-3 times per week). Short-interval dosing inhibited growth with a mean (SD) maximum Δlog10 colony forming units of -4.06 (0.52), significantly more than extended-interval dosing (P = 0.0013) every 24 hours, -2.40 (0.58), or 1-3 times per week, -2.39 (0.38). Growth recovery, an indicator of resistance emergence, occurred at 6.56 (0.70) days with short-interval dosing but was significantly earlier with extended-interval dosing (P = 0.0032) every 24 hours, 3.88 (0.85) days, and 1-3 times per week, 3.27 (1.72) days. Microbiological response correlated best with the pharmacodynamic index of %T > minimum inhibitory concentration (MIC), with an EC80 for growth inhibition of ~40%T > MIC. We used a previously published population model of amikacin to determine the probability of achieving 40%T > MIC and show that current dosing strategies are far below this target, which may partially explain why treatment failure remains so high for these infections. These data support a cautious approach to infrequent amikacin dosing for the treatment of M. abscessus.IMPORTANCEPulmonary disease caused by Mycobacterium abscessus complex (MABSC) is increasing worldwide, particularly in patients with cystic fibrosis. MABSC is challenging to treat due to high levels of antibiotic resistance. Treatment requires 2-4 antibiotics over more than 12 months and has a significant risk of toxicity but still fails to eradicate infection in over 50% of patients with cystic fibrosis. Antibiotic dosing strategies have been largely informed by common bacteria such as Pseudomonas aeruginosa. The "pharmacodynamic" effects of amikacin, a backbone of MABSC treatment, were thought to be related to maximum "peak" drug concentration, leading to daily or three times weekly dosing. However, we found that amikacin MABSC kill and growth recovery, an indicator of antibiotic resistance, are dependent on how long amikacin concentrations are above the minimum inhibitory concentration, not how high the peak concentration is. Therefore, we recommend a re-evaluation of amikacin dosing to determine if increased frequency can improve efficacy.

Intravenous N-Acetylcysteine to Prevent Cisplatin-Induced Hearing Loss in Children: A Nonrandomized Controlled Phase I Trial.

PURPOSE: Cisplatin-induced hearing loss (CIHL) is common and permanent. As compared with earlier otoprotectants, we hypothesized N-acetylcysteine (NAC) offers potential for stronger otoprotection through stimulation of glutathione (GSH) production. This study tested the optimal dose, safety, and efficacy of NAC to prevent CIHL. PATIENTS AND METHODS: In this nonrandomized, controlled phase Ia/Ib trial, children and adolescents newly diagnosed with nonmetastatic, cisplatin-treated tumors received NAC intravenously 4 hours post-cisplatin. The trial performed dose-escalation across three dose levels to establish a safe dose that exceeded the targeted peak serum NAC concentration of 1.5 mmol/L (as identified from preclinical models). Patients with metastatic disease or who were otherwise ineligible were enrolled in an observation-only/control arm. To evaluate efficacy, serial age-appropriate audiology assessments were performed. Integrated biology examined genes involved in GSH metabolism and post-NAC GSH concentrations. RESULTS: Of 52 patients enrolled, 24 received NAC and 28 were in the control arm. The maximum tolerated dose was not reached; analysis of peak NAC concentration identified 450 mg/kg as the recommended phase II dose (RP2D). Infusion-related reactions were common. No severe adverse events occurred. Compared with the control arm, NAC decreased likelihood of CIHL at the end of cisplatin therapy [OR, 0.13; 95% confidence interval (CI), 0.021-0.847; P = 0.033] and recommendations for hearing intervention at end of study (OR, 0.082; 95% CI, 0.011-0.60; P = 0.014). NAC increased GSH; GSTP1 influenced risk for CIHL and NAC otoprotection. CONCLUSIONS: NAC was safe at the RP2D, with strong evidence for efficacy to prevent CIHL, warranting further development as a next-generation otoprotectant.

Optimizing Vancomycin Therapy in Critically Ill Children: A Population Pharmacokinetics Study to Inform Vancomycin Area under the Curve Estimation Using Novel Biomarkers.

Area under the curve (AUC)-directed vancomycin therapy is recommended, but Bayesian AUC estimation in critically ill children is difficult due to inadequate methods for estimating kidney function. We prospectively enrolled 50 critically ill children receiving IV vancomycin for suspected infection and divided them into model training (n = 30) and testing (n = 20) groups. We performed nonparametric population PK modeling in the training group using Pmetrics, evaluating novel urinary and plasma kidney biomarkers as covariates on vancomycin clearance. In this group, a two-compartment model best described the data. During covariate testing, cystatin C-based estimated glomerular filtration rate (eGFR) and urinary neutrophil gelatinase-associated lipocalin (NGAL; full model) improved model likelihood when included as covariates on clearance. We then used multiple-model optimization to define the optimal sampling times to estimate AUC24 for each subject in the model testing group and compared the Bayesian posterior AUC24 to AUC24 calculated using noncompartmental analysis from all measured concentrations for each subject. Our full model provided accurate and precise estimates of vancomycin AUC (bias 2.3%, imprecision 6.2%). However, AUC prediction was similar when using reduced models with only cystatin C-based eGFR (bias 1.8%, imprecision 7.0%) or creatinine-based eGFR (bias -2.4%, imprecision 6.2%) as covariates on clearance. All three model(s) facilitated accurate and precise estimation of vancomycin AUC in critically ill children.

Growth medium and nitric oxide alter Mycobacterium abscessus morphotype and virulence.

Mycobacterium abscessus complex (MABC) infections cause significant morbidity and mortality among patients with chronic lung disease, like cystic fibrosis. MABC exists in smooth (S) and rough (R) morphotypes, but triggers of morphotype switching and associated pathogenicity or antimicrobial susceptibility are poorly understood. We demonstrate that M. abscessus subspecies abscessus (Mab), massiliense (Mms), and bolletii (Mbl) cultured in Middlebrook (MB) broth exhibit S morphotype, whereas the bacteria grown in Luria Bertani (LB) broth adopt the R morphotype, characterized by low glycopeptidolipid (GPL) expression. The components of broth that mediate this selection are complex, with albumin supplementation promoting growth of S morphotype, but not sufficient for complete selection. Consistent with the findings of other groups, R forms of Mab, Mms and Mbl selected by LB broth were internalized in RAW 264.7 macrophages with higher efficiency than S. Intracellular survival of broth-selected organisms was variable and was higher for S Mab, but lower for S Mms and Mbl. It is proposed that growth in R morphotype is induced during stress conditions, such as nutrient poor environments or during inflammation. One key component of inflammation is release of nitric oxide. We demonstrated that a nitric oxide donor (DETA-NONOate) appears to induce growth in an R morphotype, as indicated by reduced GPL expression of Mab. Mab treated with DETA-NONOate also enhanced susceptibility to azithromycin at sub-MIC concentrations. In conclusion, morphotype and macrophage intracellular bacterial load of MABC subspecies can be manipulated by growing the bacteria in different culture conditions. Nitric oxide may also drive morphotype selection and enhanced azithromycin activity against Mab and macrophage killing.

Effect of Body Fat on Population Pharmacokinetics of High-Dose Methotrexate in Pediatric Patients With Acute Lymphoblastic Leukemia.

Nearly all international regimens for pediatric acute lymphoblastic leukemia (ALL) incorporate intravenous "high-dose" methotrexate (HDMTX, ≥1 g/m2 ) to penetrate the central nervous system. Dosing is routinely adjusted for body surface area (BSA), but limited data describe the pharmacokinetics of HDMTX, particularly in obese and/or large patients. To understand the impact of body size (BSA) and body fat percentage (BFP) on HDMTX pharmacokinetics, we performed a secondary analysis of 36 children and adolescents 10-21 years old treated for newly diagnosed ALL and who were enrolled in a prospective study examining body composition. All patients received 5 g/m2 of HDMTX infused over 24 hours. Plasma methotrexate concentrations were measured at 24, 42, and 48 hours. At 48 hours, ≥0.4 µmol/L was defined as "delayed elimination," necessitating prolonged supportive care. BFP was measured using dual-energy x-ray absorptiometry. A nonparametric population pharmacokinetic model was constructed with subsequent simulations to explore effects of BSA and BFP extremes. Despite standard BSA-adjusted dosing, we found significant intrapatient variability in mean MTX concentration (38%; range, 1.2%-86%). BSA and BFP were not linearly associated with increased area under the curve (AUC, P = 0.74 and P = 0.12), but both larger size (BSA) and greater obesity (BFP) were associated with an approximately 2-fold higher risk for delayed elimination at 48 hours. HDMTX AUC was not associated with toxicity. MTX pharmacokinetics vary among and even within patients despite BSA-adjusted dosing. Obesity and large size are identified as new risk factors for delayed elimination, requiring further investigation.

Intra-individual Pharmacokinetic Variability of Intravenous Busulfan in Hematopoietic Stem Cell-Transplanted Children.

BACKGROUND: Busulfan therapeutic drug monitoring (TDM) is necessary to better achieve the target exposure in children before hematopoietic stem cell transplantation (HSCT). However, TDM-based dosing may be challenging if intra-individual pharmacokinetic variability (also denoted inter-occasion variability [IOV]) occurs during therapy. OBJECTIVES: The objectives of this study were to describe and quantify busulfan IOV in children, and to investigate its potential determinants. METHODS: We performed a new analysis of published data from children who received intravenous busulfan over 4 days before HSCT. We calculated individual pharmacokinetic parameters on each day of therapy using a published population pharmacokinetic model of busulfan and analyzed their changes. Population estimation of IOV was also performed with non-linear mixed effects (NLME) modeling. Potential predictors of significant decrease in busulfan clearance (CL) were assessed by using machine learning approaches. RESULTS: IOV could be assessed in 136 children. Between day (D) 1 and D2, most patients (80%) experienced a decrease in busulfan CL, with a median change of - 7.9%. However, both large decreases (minimum, - 48.5%) and increases in CL (maximum, + 44%) were observed. Over D1-D3 of therapy, mean CL significantly decreased (- 15%), with a decrease of ≥ 20% in 22% of patients. Some patients also showed unstable CL from day to day. NLME modeling of IOV provided a coefficient of variation of 10.6% and 13.1% for volume of distribution (Vd) and CL, respectively. Some determinants of significant decreases in busulfan CL were identified, but predictive performance of the models was limited. CONCLUSIONS: Significant busulfan intra-individual variability may occur in children who receive a HSCT and is hardly predictable. The main risk is busulfan overexposure. Performing TDM repeatedly over therapy appears to be the best way to accurately estimate busulfan exposure and perform precision dosing.

Development of Population and Bayesian Models for Applied Use in Patients Receiving Cefepime.

BACKGROUND AND OBJECTIVE: Understanding pharmacokinetic disposition of cefepime, a ß-lactam antibiotic, is crucial for developing regimens to achieve optimal exposure and improved clinical outcomes. This study sought to develop and evaluate a unified population pharmacokinetic model in both pediatric and adult patients receiving cefepime treatment. METHODS: Multiple physiologically relevant models were fit to pediatric and adult subject data. To evaluate the final model performance, a withheld group of 12 pediatric patients and two separate adult populations were assessed. RESULTS: Seventy subjects with a total of 604 cefepime concentrations were included in this study. All adults (n = 34) on average weighed 82.7 kg and displayed a mean creatinine clearance of 106.7 mL/min. All pediatric subjects (n = 36) had mean weight and creatinine clearance of 16.0 kg and 195.6 mL/min, respectively. A covariate-adjusted two-compartment model described the observed concentrations well (population model R2, 87.0%; Bayesian model R2, 96.5%). In the evaluation subsets, the model performed similarly well (population R2, 84.0%; Bayesian R2, 90.2%). CONCLUSION: The identified model serves well for population dosing and as a Bayesian prior for precision dosing.

Population Pharmacokinetics of Imipenem in Critically Ill Patients: A Parametric and Nonparametric Model Converge on CKD-EPI Estimated Glomerular Filtration Rate as an Impactful Covariate.

BACKGROUND: Population pharmacokinetic (popPK) models for antibiotics are used to improve dosing strategies and individualize dosing by therapeutic drug monitoring. Little is known about the differences in results of parametric versus nonparametric popPK models and their potential consequences in clinical practice. We developed both parametric and nonparametric models of imipenem using data from critically ill patients and compared their results. METHODS: Twenty-six critically ill patients treated with intravenous imipenem/cilastatin were included in this study. Median estimated glomerular filtration rate (eGFR) measured by the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation was 116 mL/min/1.73 m2 (interquartile range 104-124) at inclusion. The usual dosing regimen was 500 mg/500 mg four times daily. On average, five imipenem levels per patient (138 levels in total) were drawn as peak, intermediate, and trough levels. Imipenem concentration-time profiles were analyzed using parametric (NONMEM 7.2) and nonparametric (Pmetrics 1.5.2) popPK software. RESULTS: For both methods, data were best described by a model with two distribution compartments and the CKD-EPI eGFR equation unadjusted for body surface area as a covariate on the elimination rate constant (Ke). The parametric population parameter estimates were Ke 0.637 h-1 (between-subject variability [BSV]: 19.0% coefficient of variation [CV]) and central distribution volume (Vc) 29.6 L (without BSV). The nonparametric values were Ke 0.681 h-1 (34.0% CV) and Vc 31.1 L (42.6% CV). CONCLUSIONS: Both models described imipenem popPK well; the parameter estimates were comparable and the included covariate was identical. However, estimated BSV was higher in the nonparametric model. This may have consequences for estimated exposure during dosing simulations and should be further investigated in simulation studies.

Maximal concentration of intravenous busulfan as a determinant of veno-occlusive disease: a pharmacokinetic-pharmacodynamic analysis in 293 hematopoietic stem cell transplanted children.

Veno-occlusive disease (VOD) is a severe adverse reaction to busulfan-containing regimens used in the preparation of children for hematopoietic stem cell transplantation (HSCT). We conducted a retrospective analysis of data to examine determinants of VOD in children who received IV busulfan for HSCT conditioning. Busulfan PK parameters as well as various indices (maximal concentration-Cmax, area under the concentration-time curve-AUC) were estimated using a validated Bayesian approach. The influence of available PK, demographic, and clinical variables on the incidence of VOD was evaluated by using logistic regression and classification and regression tree (CART) analyses. Among the 293 patients included, the mean age was 6.5 years and the mean actual body weight was 26.3 kg. The incidence of VOD was 25.6%. Busulfan Cmax as well as weight <9 kg or age <3 years were identified as independent predictors of VOD in logistic regression analysis. CART analysis identified busulfan Cmax over the entire regimen as the strongest predictor of VOD. This study suggests that busulfan-associated VOD is in part a concentration-dependent reaction. In addition, the youngest children showed the highest risk of VOD. These findings may have important implications for busulfan dosing and therapeutic drug monitoring practice in HSCT children.

The effect of medication nonadherence on progression-free survival among patients with renal cell carcinoma.

OBJECTIVE: To examine how observed medication nonadherence to 2 second-line, oral anticancer medications (axitinib and everolimus) affects progression-free survival (PFS) among patients with renal cell carcinoma. METHODS: We used an adherence-exposure-outcome model to simulate the impact of adherence on PFS. Using a pharmacokinetic/pharmacodynamic (PK/PD) population model, we simulated drug exposure measured by area under the plasma concentration-time curve (AUC) and minimum blood or trough concentration (Cmin) under 2 scenarios: 1) optimal adherence and 2) real-world adherence. Real-world adherence was measured using the medication possession ratios as calculated from health insurance claims data. A population PK/PD model was simulated on individuals drawn from the Medical Expenditure Panel Survey (MEPS), a large survey broadly representative of the US population. Finally, we used previously published PK/PD models to estimate the effect of drug exposure (i.e., Cmin and AUC) on PFS outcomes under optimal and real-world adherence scenarios. RESULTS: Average adherence measured using medication possession ratios was 76%. After applying our simulation model to 2164 individuals in MEPS, drug exposure was significantly higher among adherent patients compared with nonadherent patients for axitinib (AUC: 249.5 vs. 159.8 ng×h/mL, P<0.001) and everolimus (AUC: 185.4 vs. 118.0 µg×h/L, P<0.001). Patient nonadherence in the real world decreased the expected PFS from an optimally adherent population by 29% for axitinib (8.4 months with optimal adherence vs. 6.0 months using real-world adherence, P<0.001) and by 5% (5.5 vs. 5.2 months, P<0.001) for everolimus. CONCLUSION: Nonadherence by renal cell carcinoma patients to second-line oral therapies significantly decreased the expected PFS.

Adipocytes Sequester and Metabolize the Chemotherapeutic Daunorubicin.

Obesity is associated with poorer outcome for many cancers. Previously, we observed that adipocytes protect acute lymphoblastic leukemia (ALL) cells from the anthracycline, daunorubicin. In this study, it is determined whether adipocytes clear daunorubicin from the tumor microenvironment (TME). Intracellular daunorubicin concentrations were evaluated using fluorescence. Daunorubicin and its largely inactive metabolite, daunorubicinol, were analytically measured in media, cells, and tissues using liquid chromatography/mass spectrometry (LC/MS). Expression of daunorubicin-metabolizing enzymes, aldo-keto reductases (AKR1A1, AKR1B1, AKR1C1, AKR1C2, AKR1C3, and AKR7A2) and carbonyl reductases (CBR1, CBR3), in human adipose tissue, were queried using public databases and directly measured by quantitative PCR (qPCR) and immunoblot. Adipose tissue AKR activity was measured by colorimetric assay. Adipocytes absorbed and efficiently metabolized daunorubicin to daunorubicinol, reducing its antileukemia effect in the local microenvironment. Murine studies confirmed adipose tissue conversion of daunorubicin to daunorubicinol in vivo Adipocytes expressed high levels of AKR and CBR isoenzymes that deactivate anthracyclines. Indeed, adipocyte protein levels of AKR1C1, AKR1C2, and AKR1C3 are higher than all other human noncancerous cell types. To our knowledge, this is the first demonstration that adipocytes metabolize and inactivate a therapeutic drug. Adipocyte-mediated daunorubicin metabolism reduces active drug concentration in the TME. These results could be clinically important for adipocyte-rich cancer microenvironments such as omentum, breast, and marrow. As AKR and CBR enzymes metabolize several drugs, and can be expressed at higher levels in obese individuals, this proof-of-principle finding has important implications across many diseases.Implications: Adipocyte absorption and metabolism of chemotherapies can reduce cytotoxicity in cancer microenvironments, potentially contributing to poorer survival outcomes. Mol Cancer Res; 15(12); 1704-13. ©2017 AACR.

A Nonparametric Method to Optimize Initial Drug Dosing and Attainment of a Target Exposure Interval: Concepts and Application to Busulfan in Pediatrics.

The traditional approach for model-based initial dosing is based on the use of a single vector of typical population parameters for targeting a specific exposure. This approach is theoretically ill-suited for targeting a range of exposure. The objective of this work was to develop a general approach for optimal (OPT) targeting of a drug exposure interval. After methodological purposes, we applied our method to the busulfan case. We used a nonparametric population pharmacokinetic model of intravenous busulfan to estimate the individual pharmacokinetic parameters of 163 bone marrow-transplanted children. Then, an array of 151 doses of busulfan ranging from 0.5 to 2 mg/kg was simulated a priori in each patient. For each dose, 29 possible busulfan plasma concentration profiles, corresponding to the nonparametric prior, each associated with a probability, were obtained. The multiple-model-based, OPT dose was identified as the dose maximizing the a priori probability of achieving the busulfan target area under the concentration-time curve (AUC). Two AUC targets were considered: 900-1500 (conventional) or <1500 µM min-1. Finally, the OPT dose was individually simulated in each patient. We compared the ability of this method to achieve the target exposure interval with that of three other traditional model-based methods and one based on the non-parametric approach. When targeting the busulfan conventional AUC range, the OPT dose provided better attainment than the best of the three other methods after one dose (82.2 vs. 41.7 %, p < 0.005), two doses (79.1 vs. 65.0 %, p < 0.005), and at the end of therapy (80.4 vs. 76.7 %, p < 0.42). The approach provided a balanced distribution between under- (10.4 %) and overexposure (9.2 %), while other approaches showed higher rates of underexposure (≥19 %). When targeting an AUC <1500 µM min, the OPT dose was successful in minimizing overexposure as 0 % of children showed simulated AUC >1500 µM min-1. Our approach has been designed to optimize the targeting of an exposure interval. When applied to busulfan in children, it outperformed the traditional model-based dosing approach, with earlier and better achievement of busulfan target AUC. The approach can be applied for OPT dosing of many drugs, when the target objective is an interval.

Evaluation of Vancomycin Exposures Associated with Elevations in Novel Urinary Biomarkers of Acute Kidney Injury in Vancomycin-Treated Rats.

Vancomycin has been associated with acute kidney injury (AKI). However, the pharmacokinetic/toxicodynamic relationship for AKI is not well defined. Allometrically scaled vancomycin exposures were used to assess the relationship between vancomycin exposure and AKI. Male Sprague-Dawley rats received clinical-grade vancomycin in normal saline (NS) as intraperitoneal (i.p.) injections for 24- to 72-h durations with doses ranging 0 to 200 mg/kg of body weight divided once or twice daily. Urine was collected over the protocol's final 24 h. Renal histopathology was qualitatively scored. Urinary biomarkers (e.g., cystatin C, clusterin, kidney injury molecule 1 [KIM-1], osteopontin, lipocalin 2/neutrophil gelatinase-associated lipocalin 2) were assayed using a Luminex xMAP system. Plasma vancomycin concentrations were assayed by high-performance liquid chromatography with UV detection. A three-compartment vancomycin pharmacokinetic model was fit to the data with the Pmetrics package for R. The exposure-response in the first 24 h was evaluated using Spearman's nonparametric correlation coefficient (rs) values for the area under the concentration-time curve during the first 24 h (AUC0-24), the maximum concentration in plasma during the first 24 h (Cmax0-24 ), and the lowest (minimum) concentration in plasma after the dose closest to 24 h (Cmin0-24 ). A total of 52 rats received vancomycin (n = 42) or NS (n = 10). The strongest exposure-response correlations were observed between AUC0-24 and Cmax0-24 and urinary AKI biomarkers. Exposure-response correlations (rs values) for AUC0-24, Cmax0-24 , and Cmin0-24 were 0.37, 0.39, and 0.22, respectively, for clusterin; 0.42, 0.45, and 0.26, respectively, for KIM-1; and 0.52, 0.55, and 0.42, respectively, for osteopontin. However, no differences in histopathological scores were observed. Optimal sampling times after administration of the i.p. dose were 0.25, 0.75, 2.75, and 8 h for the once-daily dosing schemes and 0.25, 1.25, 14.5, and 17.25 h for the twice-daily dosing schemes. Our observations suggest that AUC0-24 or Cmax0-24 correlates with increases in urinary AKI biomarkers.

Accurately Achieving Target Busulfan Exposure in Children and Adolescents With Very Limited Sampling and the BestDose Software.

BACKGROUND: Busulfan dose adjustment is routinely guided by plasma concentration monitoring using 4-9 blood samples per dose adjustment, but a pharmacometric Bayesian approach could reduce this sample burden. METHODS: The authors developed a nonparametric population model with Pmetrics. They used it to simulate optimal initial busulfan dosages, and in a blinded manner, they compared dosage adjustments using the model in the BestDose software to dosage adjustments calculated by noncompartmental estimation of area under the time-concentration curve at a national reference laboratory in a cohort of patients not included in model building. RESULTS: Mean (range) age of the 53 model-building subjects was 7.8 years (0.2-19.0 years) and weight was 26.5 kg (5.6-78.0 kg), similar to nearly 120 validation subjects. There were 16.7 samples (6-26 samples) per subject to build the model. The BestDose cohort was also diverse: 10.2 years (0.25-18 years) and 46.4 kg (5.2-110.9 kg). Mean bias and imprecision of the 1-compartment model-predicted busulfan concentrations were 0.42% and 9.2%, and were similar in the validation cohorts. Initial dosages to achieve average concentrations of 600-900 ng/mL were 1.1 mg/kg (≤12 kg, 67% in the target range) and 1.0 mg/kg (>12 kg, 76% in the target range). Using all 9 concentrations after dose 1 in the Bayesian estimation of dose requirements, the mean (95% confidence interval) bias of BestDose calculations for the third dose was 0.2% (-2.4% to 2.9%, P = 0.85), compared with the standard noncompartmental method based on 9 concentrations. With 1 optimally timed concentration 15 minutes after the infusion (calculated with the authors' novel MMopt algorithm) bias was -9.2% (-16.7% to -1.5%, P = 0.02). With 2 concentrations at 15 minutes and 4 hours bias was only 1.9% (-0.3% to 4.2%, P = 0.08). CONCLUSIONS: BestDose accurately calculates busulfan intravenous dosage requirements to achieve target plasma exposures in children up to 18 years of age and 110 kg using only 2 blood samples per adjustment compared with 6-9 samples for standard noncompartmental dose calculations.

Optimal timing of oral fosfomycin administration for pre-prostate biopsy prophylaxis.

OBJECTIVES: As the optimal administration time for fosfomycin peri-procedural prophylaxis is unclear, we sought to determine optimal administration times for fosfomycin peri-procedural prophylaxis. METHODS: Plasma, peripheral zone and transition zone fosfomycin concentrations were obtained from 26 subjects undergoing transurethral resection of the prostate (TURP), following a single oral dose of 3 g of fosfomycin. Population pharmacokinetic modelling was completed with the Nonparametric Adaptive Grid (NPAG) algorithm (Pmetrics package for R), with a four-compartment model. Plasma and tissue concentrations were simulated during the first 24 h post-dose, comparing these with EUCAST susceptibility breakpoints for Escherichia coli, a common uropathogen. RESULTS: Non-compartmental-determined pharmacokinetic values in our population were similar to those reported in the package insert. Predicted plasma concentrations rapidly increased after the first hour, giving more than 90% population coverage for organisms with an MIC ≤4 mg/L over the first 12 h post-dose. Organisms with higher MICs fared much worse, with organisms at the EUCAST breakpoint being covered for <10% of the population at any time. Transitional zone prostate concentrations exceeded 4 mg/L for 90% of the population between hours 1 and 9. Peripheral zone prostate concentrations were much lower and only exceeded 4 mg/L for 70% of the population between hours 1 and 4. CONCLUSIONS: Until more precise plasma and tissue data are available, we recommend that fosfomycin prophylaxis be given 1-4 h prior to prostate biopsy. We do not recommend fosfomycin prophylaxis for subjects with known organisms with MICs >4 mg/L.

Pharmacokinetic tools for the dose adjustment of ciclosporin in haematopoietic stem cell transplant patients.

AIMS: Ciclosporin A (CsA) is used in the prophylaxis and treatment of acute and chronic graft vs. host disease after haematopoietic stem cell (HSCT) transplantation. Our objective was to build and compare three independent Bayesian estimators of CsA area under the curve (AUC) using a limited sampling strategy (LSS), to assist in dose adjustment. METHODS: The Bayesian estimators were developed using in parallel: two independent parametric modelling approaches (nonmem® and iterative two stage (ITS) Bayesian modelling) and the non-parametric adaptive grid method (Pmetrics®). Seventy-two full pharmacokinetic profiles (at pre-dose and 0.33, 0.66, 1, 2, 3, 4, 6, 8 and 12h after dosing) collected from 40 HSCT patients given CsA were used to build the pharmacokinetic models, while 15 other profiles (n = 7) were kept for validation. For each Bayesian estimator, AUCs estimated using the full profiles were compared with AUCs estimated using three samples. RESULTS: The pharmacokinetic profiles were well fitted using a two compartment model with first order elimination, combined with a gamma function for the absorption phase with ITS and Pmetrics or an Erlang distribution with nonmem. The derived Bayesian estimators based on a C0-C1 h-C4 h sampling schedule (best LSS) accurately estimated CsA AUC(0,12 h) in the validation group (n = 15; nonmem: bias (mean ± SD)/RMSE 2.05% ± 13.31%/13.02%; ITS: 4.61% ± 10.56%/11.20%; Pmetrics: 0.30% ± 10.12%/10.47%). The dose chosen confronting the three results led to a pertinent dose proposal. CONCLUSIONS: The developed Bayesian estimators were all able to predict ciclosporin AUC(0,12 h) in HSCT patients using only three blood with minimal bias and may be combined to increase the reliability of CsA dose adjustment in routine.

Use of nucleoside reverse transcriptase inhibitor-only regimens in HIV-infected children and adolescents.

OBJECTIVE: In adults, nucleoside reverse transcriptase inhibitor-only antiretroviral regimens (NOARs) with ≥3 nucleoside reverse transcriptase inhibitors are less potent than highly active antiretroviral therapy (HAART). Published pediatric experience with NOARs is limited; thus, we wished to better define the virological, immunological and toxicological effects of NOARs in children and adolescents. METHODS: We analyzed data from NOAR-treated participants in LEGACY, a multicenter observational cohort study of HIV-infected children and adolescents. NOAR-treated case-participants were matched to participants without prior NOAR who initiated HAART during the same year for comparison. RESULTS: Of 575 participants with data from time of HIV diagnosis through 2006, 67 (12%) received NOARs for at least 24 weeks; most (46%) received the fixed dose combination of zidovudine/lamivudine/abacavir. NOAR use peaked in 2001 to 2002. NOAR-treated participants were significantly older and more treatment experienced than HAART-treated participants. Virologic outcomes, including the percentage of participants with a plasma HIV RNA viral load <400 copies/mL at week 24 (47% versus 34%) and the mean 24-week change in log10 plasma HIV RNA viral load from baseline (-0.63 versus -1.02), were similar between NOAR- and HAART-treated participants, but virologic rebound was more likely in NOAR-treated participants (77% versus 54%, P = 0.02). Increase in CD4 percentage points from baseline to 24 weeks was negligible in NOAR-treated participants compared with HAART-treated participants (0.95% versus 10.1%, P < 0.001). Anemia and leukopenia were more commonly reported with NOARs than HAART. DISCUSSION: Week 24 virologic outcomes were similar between NOAR- and HAART-treated participants, but NOAR durability was poorer and their use was associated with less immunologic reconstitution. NOARs should play a limited role in pediatric and adolescent antiretroviral therapy.

Pharmacodynamics of liposomal amphotericin B and flucytosine for cryptococcal meningoencephalitis: safe and effective regimens for immunocompromised patients.

BACKGROUND: Cryptococcal meningoencephalitis is a lethal infection with relatively few therapeutic options. The optimal dosage of liposomal amphotericin B (LAmB) alone or in combination with flucytosine is not known. METHODS: A murine model of cryptococcal meningoencephalitis was used. The fungal density in the brain was determined using quantitative cultures. Pharmacokinetic-pharmacodynamic relationships were determined for LAmB and flucytosine administered alone. The effect of the combination was described using the Greco model and a mathematical model. The results were bridged to humans. RESULTS: Inoculation resulted in hematogenous dissemination and logarithmic growth within the central nervous system. There was histological evidence of multifocal infection throughout the brain. Both LAmB and flucytosine produced a dose-dependent reduction in fungal burden. The effect of the combination of agents in the brain was additive. Bridging studies suggested that a human dosage of LAmB 3 mg/kg/d resulted in a submaximal antifungal effect. Regimens of LAmB 6 mg/kg/d alone, LAmB 3 mg/kg/d plus flucytosine 50 mg/kg/d, and LAmB 3 mg/kg/d plus flucytosine 100 mg/kg/d all resulted in near-maximal antifungal activity. CONCLUSIONS: Potential regimens for further study in clinical trials include LAmB 6 mg/kg/d alone, LAmB 3 mg/kg/d plus flucytosine 50 mg/kg/d, and LAmB 3 mg/kg/d plus flucytosine 100 mg/kg/d.

Software for dosage individualization of voriconazole for immunocompromised patients.

The efficacy of voriconazole is potentially compromised by considerable pharmacokinetic variability. There are increasing insights into voriconazole concentrations that are safe and effective for treatment of invasive fungal infections. Therapeutic drug monitoring is increasingly advocated. Software to aid in the individualization of dosing would be an extremely useful clinical tool. We developed software to enable the individualization of voriconazole dosing to attain predefined serum concentration targets. The process of individualized voriconazole therapy was based on concepts of Bayesian stochastic adaptive control. Multiple-model dosage design with feedback control was used to calculate dosages that achieved desired concentration targets with maximum precision. The performance of the software program was assessed using the data from 10 recipients of an allogeneic hematopoietic stem cell transplant (HSCT) receiving intravenous (i.v.) voriconazole. The program was able to model the plasma concentrations with a high level of precision, despite the wide range of concentration trajectories and interindividual pharmacokinetic variability. The voriconazole concentrations predicted after the last dosages were largely concordant with those actually measured. Simulations provided an illustration of the way in which the software can be used to adjust dosages of patients falling outside desired concentration targets. This software appears to be an extremely useful tool to further optimize voriconazole therapy and aid in therapeutic drug monitoring. Further prospective studies are now required to define the utility of the controller in daily clinical practice.