Non-linear IV pharmacokinetics of the ATR inhibitor berzosertib (M6620) in mice.

BACKGROUND: The Ataxia Telangiectasia and Rad3-related (ATR) protein complex is an apical initiator of DNA damage response pathways. Several ATR inhibitors (ATRi) are in clinical development including berzosertib (formerly M6620, VX-970). Although clinical studies have examined plasma pharmacokinetics (PK) in humans, little is known regarding dose/exposure relationships and tissue distribution. To understand these concepts, we extensively characterized the PK of berzosertib in mouse plasma and tissues. METHODS: A highly sensitive LC-MS/MS method was utilized to quantitate berzosertib in plasma and tissues. Dose proportionality was assessed in female BALB/c mice following single IV doses (2, 6, 20 or 60 mg/kg). A more extensive PK study was conducted in tumor-bearing mice following a single IV dose of 20 mg/kg to evaluate distribution to tissues. PK parameters were calculated by non-compartmental analysis (NCA). A compartmental model was developed to describe the PK behavior of berzosertib. Plasma protein binding was determined in vitro. RESULTS: Increased doses of berzosertib were associated with less than proportional increases in early plasma concentrations and greater than proportional increase in tissue exposure, attributable to saturation of plasma protein binding. Berzosertib extensively distributed into bone marrow, tumor, thymus, and lymph nodes, however; brain and spinal cord exposure was less than plasma. CONCLUSION: The nonlinear PK of berzosertib displayed here can be attributed to saturation of plasma protein binding and occurred at concentrations close to those observed in clinical trials. Our results will help to understand preclinical pharmacodynamic and toxicity data and to inform optimal dosing and deployment of berzosertib.

An in-silico modeling approach to separate exogenous and endogenous plasma insulin appearance, with application to inhaled insulin.

The aim of this study was to develop a dynamic model-based approach to separately quantify the exogenous and endogenous contributions to total plasma insulin concentration and to apply it to assess the effects of inhaled-insulin administration on endogenous insulin secretion during a meal test. A three-step dynamic in-silico modeling approach was developed to estimate the two insulin contributions of total plasma insulin in a group of 21 healthy subjects who underwent two equivalent standardized meal tests on separate days, one of which preceded by inhalation of a Technosphere® Insulin dose (22U or 20U). In the 30-120 min test interval, the calculated endogenous insulin component showed a divergence in the time course between the test with and without inhaled insulin. Moreover, the supra-basal area-under-the-curve of endogenous insulin in the test with inhaled insulin was significantly lower than that in the test without (2.1 ± 1.7 × 104 pmol·min/L vs 4.2 ± 1.8 × 104 pmol·min/L, p < 0.01). The percentage of exogenous insulin reaching the plasma, relative to the inhaled dose, was 42 ± 21%. The proposed in-silico approach separates exogenous and endogenous insulin contributions to total plasma insulin, provides individual bioavailability estimates, and can be used to assess the effect of inhaled insulin on endogenous insulin secretion during a meal.

Long-term forecasting of a motor outcome following rehabilitation in chronic stroke via a hierarchical bayesian dynamic model.

BACKGROUND: Given the heterogeneity of stroke, it is important to determine the best course of motor therapy for each patient, i.e., to personalize rehabilitation based on predictions of long-term outcomes. Here, we propose a hierarchical Bayesian dynamic (i.e., state-space) model (HBDM) to forecast long-term changes in a motor outcome due to rehabilitation in the chronic phase post-stroke. METHODS: The model incorporates the effects of clinician-supervised training, self-training, and forgetting. In addition, to improve forecasting early in rehabilitation, when data are sparse or unavailable, we use the Bayesian hierarchical modeling technique to incorporate prior information from similar patients. We use HBDM to re-analyze the Motor Activity Log (MAL) data of participants with chronic stroke included in two clinical trials: (1) the DOSE trial, in which participants were assigned to a 0, 15, 30, or 60-h dose condition (data of 40 participants analyzed), and (2) the EXCITE trial, in which participants were assigned a 60-h dose, in either an immediate or a delayed condition (95 participants analyzed). RESULTS: For both datasets, HBDM accounts well for individual dynamics in the MAL during and outside of training: mean RMSE = 0.28 for all 40 DOSE participants (participant-level RMSE 0.26 ± 0.19-95% CI) and mean RMSE = 0.325 for all 95 EXCITE participants (participant-level RMSE 0.32 ± 0.31), which are small compared to the 0-5 range of the MAL. Bayesian leave-one-out cross-validation shows that the model has better predictive accuracy than static regression models and simpler dynamic models that do not account for the effect of supervised training, self-training, or forgetting. We then showcase model's ability to forecast the MAL of "new" participants up to 8 months ahead. The mean RMSE at 6 months post-training was 1.36 using only the baseline MAL and then decreased to 0.91, 0.79, and 0.69 (respectively) with the MAL following the 1st, 2nd, and 3rd bouts of training. In addition, hierarchical modeling improves prediction for a patient early in training. Finally, we verify that this model, despite its simplicity, can reproduce previous findings of the DOSE trial on the efficiency, efficacy, and retention of motor therapy. CONCLUSIONS: In future work, such forecasting models can be used to simulate different stages of recovery, dosages, and training schedules to optimize rehabilitation for each person. Trial registration This study contains a re-analysis of data from the DOSE clinical trial ID NCT01749358 and the EXCITE clinical trial ID NCT00057018.

Monoclonal Antibody Pharmacokinetics in Cynomolgus Monkeys Following Subcutaneous Administration: Physiologically Based Model Predictions from Physiochemical Properties.

An integrated physiologically based modeling framework is presented for predicting pharmacokinetics and bioavailability of subcutaneously administered monoclonal antibodies in cynomolgus monkeys, based on in silico structure-derived metrics characterizing antibody size, overall charge, local charge, and hydrophobicity. The model accounts for antibody-specific differences in pinocytosis, transcapillary transport, local lymphatic uptake, and pre-systemic degradation at the subcutaneous injection site and reliably predicts the pharmacokinetics of five different wild-type mAbs and their Fc variants following intravenous and subcutaneous administration. Significant associations were found between subcutaneous injection site degradation rate and the antibody's local positive charge of its complementarity-determining region (R = 0.56, p = 0.0012), antibody pinocytosis rate and its overall positive charge (R = 0.59, p = 0.00063), and antibody paracellular transport and its overall charge together with hydrophobicity (R = 0.63, p = 0.00096). Based on these results, population simulations were performed to predict the relationship between bioavailability and antibody local positive charge. In addition, model simulations were conducted to calculate the relative contribution of absorption pathways (lymphatic and blood), pre-systemic degradation pathways (interstitial and lysosomal), and the influence of injection site lymph flow on antibody bioavailability and pharmacokinetics. The proposed physiologically based modeling framework integrates fundamental mechanisms governing antibody subcutaneous absorption and disposition, with structured-based physiochemical properties, to predict antibody bioavailability and pharmacokinetics in vivo.

Analysis of Complex Absorption After Multiple Dosing: Application to the Interaction Between the P-glycoprotein Substrate Talinolol and Rifampicin.

PURPOSE: In order to clarify the effect of rifampicin on the bioavailability of the P-glycoprotein substrate talinolol, its absorption kinetics was modeled after multiple-dose oral administration of talinolol in healthy subjects. METHODS: A sum of two inverse Gaussian functions was used to calculate the time course of the input rate into the systemic circulation. RESULTS: The estimated rate of drug entry into the systemic circulation revealed two distinct peaks at 1 and 3.5 h after administration. Rifampicin did not affect bioavailability of talinolol, but did shift the second peak of the input function by 1.3 h to later times. Elimination clearance and one of the intercompartmental distribution clearances increased significantly under rifampicin treatment. CONCLUSIONS: Rifampicin changes the time course of absorption rate but not the fraction absorbed of talinolol. The model suggests the existence of two intestinal absorption windows for talinolol.

Physiologically Based Modeling to Predict Monoclonal Antibody Pharmacokinetics in Humans from in vitro Physiochemical Properties.

A model-based framework is presented to predict monoclonal antibody (mAb) pharmacokinetics (PK) in humans based on in vitro measures of antibody physiochemical properties. A physiologically based pharmacokinetic (PBPK) model is used to explore the predictive potential of 14 in vitro assays designed to measure various antibody physiochemical properties, including nonspecific cell-surface interactions, FcRn binding, thermal stability, hydrophobicity, and self-association. Based on the mean plasma PK time course data of 22 mAbs from humans reported in the literature, we found a significant positive correlation (R = 0.64, p = .0013) between the model parameter representing antibody-specific vascular to endothelial clearance and heparin relative retention time, an in vitro measure of nonspecific binding. We also found that antibody-specific differences in paracellular transport due to convection and diffusion could be partially explained by antibody heparin relative retention time (R = 0.52, p = .012). Other physiochemical properties, including antibody thermal stability, hydrophobicity, cross-interaction and self-association, in and of themselves were not predictive of model-based transport parameters. In contrast to other studies that have reported empirically derived expressions relating in vitro measures of antibody physiochemical properties directly to antibody clearance, the proposed PBPK model-based approach for predicting mAb PK incorporates fundamental mechanisms governing antibody transport and processing, informed by in vitro measures of antibody physiochemical properties, and can be expanded to include more descriptive representations of each of the antibody processing subsystems, as well as other antibody-specific information.

Open-source maximum a posteriori-bayesian dosing AdDS to current therapeutic drug monitoring: Adapting to the era of individualized therapy.

Recent updates in the therapeutic drug monitoring (TDM) guidelines for vancomycin have rekindled interest in maximum a posteriori-Bayesian (MAP-Bayesian) estimation of patient-specific pharmacokinetic parameters. To create a versatile infrastructure for MAP-Bayesian dosing of vancomycin or other drugs, a freely available, R-based software package, Advanced Dosing Solutions (AdDS), was created to facilitate clinical implementation of these improved TDM methods. The objective of this study was to utilize AdDS for pre- and post-processing of data in order to streamline the therapeutic management of vancomycin in healthy and obese veterans. Patients from a local Veteran Affairs hospital were utilized to compare the process of full re-estimation versus Bayesian updating of priors on healthy adult and obese patient populations for use with AdDS. Twenty-four healthy veterans were utilized to train (14/24) and test (10/24) the base pharmacokinetic model of vancomycin while comparing the effects of updated and fully re-estimated priors. This process was repeated with a total of 18 obese veterans for both training (11/18) and testing (7/18). Comparison of MAP objective function between the original and re-estimated models for healthy adults indicated that 78.6% of the subjects in the training and 70.0% of the subjects in the testing datasets had similar or improved predictions by the re-estimated model. For obese veterans, 81.8% of subjects in the training dataset and 85.7% of subjects in the testing dataset had similar or improved predictions. Re-estimation of model parameters provided more significant improvements in objective function compared with Bayesian updating, which may be a useful strategy in cases where sufficient samples and subjects are available. The generation of bespoke regimens based on patient-specific clearance and minimal sampling may improve patient care by addressing fundamental pharmacokinetic differences in healthy and obese veteran populations.

Characteristics of Passive Solute Transport across Primary Rat Alveolar Epithelial Cell Monolayers.

Primary rat alveolar epithelial cell monolayers (RAECM) were grown without (type I cell-like phenotype, RAECM-I) or with (type II cell-like phenotype, RAECM-II) keratinocyte growth factor to assess passive transport of 11 hydrophilic solutes. We estimated apparent permeability (Papp) in the absence/presence of calcium chelator EGTA to determine the effects of perturbing tight junctions on "equivalent" pores. Papp across RAECM-I and -II in the absence of EGTA are similar and decrease as solute size increases. We modeled Papp of the hydrophilic solutes across RAECM-I/-II as taking place via heterogeneous populations of equivalent pores comprised of small (0.41/0.32 nm radius) and large (9.88/11.56 nm radius) pores, respectively. Total equivalent pore area is dominated by small equivalent pores (99.92-99.97%). The number of small and large equivalent pores in RAECM-I was 8.55 and 1.29 times greater, respectively, than those in RAECM-II. With EGTA, the large pore radius in RAECM-I/-II increased by 1.58/4.34 times and the small equivalent pore radius increased by 1.84/1.90 times, respectively. These results indicate that passive diffusion of hydrophilic solutes across an alveolar epithelium occurs via small and large equivalent pores, reflecting interactions of transmembrane proteins expressed in intercellular tight junctions of alveolar epithelial cells.

An Analysis of Glucose Effectiveness in Subjects With or Without Type 2 Diabetes via Hierarchical Modeling.

Glucose effectiveness, defined as the ability of glucose itself to increase glucose utilization and inhibit hepatic glucose production, is an important mechanism maintaining normoglycemia. We conducted a minimal modeling analysis of glucose effectiveness at zero insulin (GEZI) using intravenous glucose tolerance test data from subjects with type 2 diabetes (T2D, n=154) and non-diabetic (ND) subjects (n=343). A hierarchical statistical analysis was performed, which provided a formal mechanism for pooling the data from all study subjects, to yield a single composite population model that quantifies the role of subject specific characteristics such as weight, height, age, sex, and glucose tolerance. Based on the resulting composite population model, GEZI was reduced from 0.021 min-1 (standard error - 0.00078 min-1) in the ND population to 0.011 min-1 (standard error - 0.00045 min-1) in T2D. The resulting model was also employed to calculate the proportion of the non-insulin-dependent net glucose uptake in each subject receiving an intravenous glucose load. Based on individual parameter estimates, the fraction of total glucose disposal independent of insulin was 72.8% ± 12.0% in the 238 ND subjects over the course of the experiment, indicating the major contribution to the whole-body glucose clearance under non-diabetic conditions. This fraction was significantly reduced to 48.8% ± 16.9% in the 30 T2D subjects, although still accounting for approximately half of the total in the T2D population based on our modeling analysis. Given the potential application of glucose effectiveness as a predictor of glucose intolerance and as a potential therapeutic target for treating diabetes, more investigations of glucose effectiveness in other disease conditions can be conducted using the hierarchical modeling framework reported herein.

A whole-body circulatory neutrophil model with application to predicting clinical neutropenia from in vitro studies.

A circulatory model of granulopoiesis and its regulation is presented that includes neutrophil trafficking in the lungs, liver, spleen, bone marrow, lymph nodes, and blood. In each organ, neutrophils undergo transendothelial migration from vascular to interstitial space, clearance due to apoptosis, and recycling via the lymphatic flow. The model includes cell cycling of progenitor cells in the bone marrow, granulocyte colony-stimulating factor (G-CSF) kinetics and its neutrophil regulatory action, as well as neutrophil margination in the blood. From previously reported studies, 111 In-labeled neutrophil kinetic data in the blood and sampled organs were used to estimate the organ trafficking parameters in the model. The model was further developed and evaluated using absolute neutrophil count (ANC), band cell, and segmented neutrophil time course data from healthy volunteers following four dose levels of pegfilgrastim (r2  = 0.77-0.99), along with ANC time course responses following filgrastim (r2  = 0.96). The baseline values of various cell types in bone marrow and blood, as well as G-CSF concentration in the blood, predicted by the model are consistent with available literature reports. After incorporating the mechanism of action of both paclitaxel and carboplatin, as determined from an in vitro bone marrow studies, the model reliably predicted the observed ANC time course following paclitaxel plus carboplatin observed in a phase I trial of 46 patients (r2  = 0.70). The circulatory neutrophil model may provide a mechanistic framework for predicting multi-organ neutrophil homeostasis and dynamics in response to therapeutic agents that target neutrophil dynamics and trafficking in different organs.

Biokinetic modeling of nanoparticle interactions with lung alveolar epithelial cells: uptake, intracellular processing, and egress.

Studies on health effects of engineered nanomaterials (ENMs) in the lung have provided information on ENM toxicity and translocation across airway and alveolar epithelial barriers. Various inhaled ENMs (e.g., gold and iridium nanoparticles) have been reported to partially cross the air-blood barrier in the lung, enter the vasculature, and distribute in several end organs, including the heart, liver, spleen, and kidney. Using an in vitro primary rat alveolar epithelial cell (AEC) monolayer model, we reported transport rates of relatively nontoxic polystyrene nanoparticles (PNPs), which appear to be taken up via nonendocytic processes into AECs. PNPs internalized into cytoplasm then trigger autophagy, followed by delivery of PNPs from autophagosomes into lysosomes, from where PNPs are exocytosed. We used the data from these experiments to perform biokinetic modeling that incorporates the processes associated with internalization and intracellular distribution of PNPs, autophagy, lysosomal exocytosis of PNPs, and several putative mechanisms of action that extend our previous understanding of AEC processing of PNPs. Results suggest that entry of PNPs into AECs, subsequent activation of autophagy by cytosolic PNPs, accumulation of PNPs in lysosomes, and lysosomal exocytosis are interwoven by proposed regulatory mechanisms.

Predicting Chemotherapy-Induced Neutropenia and Granulocyte Colony-Stimulating Factor Response Using Model-Based In Vitro to Clinical Translation.

The ability to predict the incidence of chemotherapy-induced neutropenia in early drug development can inform risk monitoring and mitigation strategies, as well as decisions on advancing compounds to clinical trials. In this report, a physiological model of granulopoiesis that incorporates the drug's mechanism of action on cell cycle proliferation of bone marrow progenitor cells was extended to include the action of the cytotoxic agents paclitaxel, carboplatin, doxorubicin, and gemcitabine. In vitro bone marrow studies were conducted with each compound, and results were used to determine the model's drug effect parameters. Population simulations were performed to predict the absolute neutrophil count (ANC) and incidence of neutropenia for each compound, which were compared to results reported in the literature. In addition, using the single agent in vitro study results, the model was able to predict ANC time course in response to paclitaxel plus carboplatin in combination, which compared favorably to the results reported in a phase 1 clinical trial of 46 patients (r2 = 0.70). Model simulations were used to compare the relative risk (RR) of neutropenia in patients with high baseline ANCs for five chemotherapeutic regimens: doxorubicin (RR = 0.59), paclitaxel plus carboplatin combination (RR = 0.079), carboplatin (RR = 0.047), paclitaxel (RR = 0.031), and gemcitabine (RR = 0.013). Finally, the model was applied to quantify the reduced incidence of neutropenia with coadministration of pegfilgrastim or filgrastim, for both paclitaxel and the combination of paclitaxel plus carboplatin. The model provides a framework for predicting clinical neutropenia using in vitro bone marrow studies of anticancer agents that may guide drug development decisions.

Population Pharmacokinetic Modeling of Vancomycin in Thai Patients With Heterogeneous and Unstable Renal Function.

BACKGROUND: Vancomycin is widely used to treat gram-positive bacterial infections. However, given significant interpatient variability in its pharmacokinetics, maintaining plasma concentrations is difficult within its characteristically narrow therapeutic window. This is especially challenging in patients with unstable renal function. Thus, the aim of this study was to develop a population pharmacokinetic model for vancomycin that is suitable for Thai patients with variable renal functions, including those with unstable renal function. METHODS: Data from 213 patients, including 564 blood samples, were retrospectively collected; approximately 70% patients exhibited unstable renal function during vancomycin treatment. The model building group was randomly assigned 108 patients and the remaining 33 patients comprised the validation group. A population pharmacokinetic model was developed that incorporated drug clearance (CL) as a function of time-varying creatine clearance (CrCL). The predictive ability of the resulting population model was evaluated using the validation data set, including its ability to forecast serum concentrations within a Bayesian feedback algorithm. RESULTS: A 2-compartment model with drug CL values that changed with time-varying CrCL adequately described vancomycin pharmacokinetics in the evaluated heterogeneous patient population with unstable renal function. Vancomycin CL was related to time-varying CrCL as follows: CL (t) = 0.11 + 0.021 × CrCL (t) (CrCL <120 mL/min. Using the population model, Bayesian estimation with at least one measured serum concentration resulted in a forecasting error of small bias (-2.4%) and adequate precision (31.5%). CONCLUSIONS: In hospitals with a high incidence of unstable renal function, incorporating time-varying CrCL with Bayesian estimation and at least one measured drug concentration, along with frequent CrCL monitoring, improves the predictive performance of therapeutic drug monitoring of vancomycin.

The Efficiency, Efficacy, and Retention of Task Practice in Chronic Stroke.

In motor skill learning, larger doses of practice lead to greater efficacy of practice, lower efficiency of practice, and better long-term retention. Whether such learning principles apply to motor practice after stroke is unclear. Here, we developed novel mixed-effects models of the change in the perceived quality of arm movements during and following task practice. The models were fitted to data from a recent randomized controlled trial of the effect of dose of task practice in chronic stroke. Analysis of the models' learning and retention rates demonstrated an increase in efficacy of practice with greater doses, a decrease in efficiency of practice with both additional dosages and additional bouts of training, and fast initial decay following practice. Two additional effects modulated retention: a positive "self-practice" effect, and a negative effect of dose. Our results further suggest that for patients with sufficient arm use post-practice, self-practice will further improve use.

Predicting monoclonal antibody pharmacokinetics following subcutaneous administration via whole-body physiologically-based modeling.

Use of the subcutaneous (SC) route for administering monoclonal antibodies (mAbs) to treat chronic conditions has been hindered because of an incomplete understanding of fundamental mechanisms controlling mAb absorption from the SC site, and due to the limited translatability of preclinical studies. In this paper, we report on the development and evaluation of a whole-body physiologically-based model to predict mAb pharmacokinetics following SC administration. The circulatory model is based on the physiological processes governing mAb transport and includes two mAb-specific parameters representing differences in pinocytosis rate and the diffusive/convective transport rates among mAbs. At the SC administration site, two additional parameters are used to represent mAb differences in lymphatic capillary uptake and in pre-systemic clearance. Model development employed clinical intravenous (IV) plasma PK data from 20 mAbs and SC plasma PK data from 12 of these mAbs, as obtained from the literature. The resulting model reliably described both the IV and SC measured plasma concentration data. In addition, a metric based on the positive charge across the mAb's complementarity determining region vicinity was found to positively correlate with the model-based estimates of the mAb-specific parameter governing organ/tissue pinocytosis transport and with estimates of the mAb's SC lymphatic capillary clearance. These two relationships were incorporated into the model and accurately predicted the SC PK profiles of three out of four separate mAbs not included in model development. The whole-body physiologically-based model reported herein, provides a platform to characterize and predict the plasma disposition of monoclonal antibodies following SC administration in humans.

A physiological model of granulopoiesis to predict clinical drug induced neutropenia from in vitro bone marrow studies: with application to a cell cycle inhibitor.

Neutropenia is one of the most common dose-limiting toxocities associated with anticancer drug therapy. The ability to predict the probability and severity of neutropenia based on in vitro studies of drugs in early drug development will aid in advancing safe and efficacious compounds to human testing. Toward this end, a physiological model of granulopoiesis and its regulation is presented that includes the bone marrow progenitor cell cycle, allowing for a mechanistic representation of the action of relevant anticancer drugs based on in vitro studies. Model development used data from previously reported tracer kinetic studies of granulocyte disposition in healthy humans to characterize the dynamics of neutrophil margination in the presence of endogenous granulocyte-colony stimulating factor (G-CSF). In addition, previously published data from healthy volunteers following pegfilgrastim and filgrastim were used to quantify the regulatory effects of support G-CSF therapies on granulopoiesis. The model was evaluated for the cell cycle inhibitor palbociclib, using an in vitro system of human bone marrow mononuclear cells to quantify the action of palbociclib on proliferating progenitor cells, including its inhibitory effect on G1 to S phase transition. The in vitro results were incorporated into the physiological model of granulopoiesis and used to predict the time course of absolute neutrophil count (ANC) and the incidence of neutropenia observed in three previously reported clinical trials of palbociclib. The model was able to predict grade 3 and 4 neutropenia due to palbociclib treatment with 86% accuracy based on in vitro data.

Analytical solution of linear multi-compartment models with non-zero initial condition and its implementation with R.

The analytical solution for multi-compartment models with a non-zero initial condition is complex because of the inter-compartmental transfer. An elegant solution and its implementation in the 'wnl' R package can be useful in solving examples of textbooks and developing software of therapeutic drug monitoring, pharmacokinetic simulation, and parameter estimation. This solution uses Laplace transformation, convolution, matrix inversion, and the fact that the general solution of an inhomogeneous ordinary differential equation is the sum of a homogenous and a particular solution, together.

Pharmacokinetics of Tedizolid in Plasma and Sputum of Adults with Cystic Fibrosis.

Over the past decade, the prevalence of infections involving methicillin-resistant Staphylococcus aureus (MRSA) in patients with cystic fibrosis (CF) has increased significantly. Tedizolid (TZD) demonstrates excellent activity against MRSA and a favorable safety profile. The pharmacokinetics of several antibiotics have been shown to be altered in CF patients. The purpose of this study was to characterize the pharmacokinetics of tedizolid in this population. Eleven patients with CF were randomized to receive tedizolid phosphate at 200 mg orally or intravenously once daily for 3 doses with a minimum 2-day washout, followed by crossover to the remaining dosage form. Plasma and expectorated sputum were collected following the third dose of each dosage form for analysis. Population pharmacokinetic analysis was performed using the maximum likelihood expectation maximization method, and the disposition of TZD was described by a two-compartment model. The sputum concentrations exceeded the unbound plasma concentrations with an estimated mean sputum-to-unbound plasma penetration ratio of 2.88 (coefficient of variation, 50.3%). The estimated population mean ± standard deviation of total clearance, central volume of distribution, and bioavailability were 9.72 ± 1.62 liters/h, 61.6 ± 6.94 liters, and 1.04 ± 0.232, respectively. The total clearance was higher in CF patients than in healthy volunteers; however, it was similar to published data for patients with complicated skin and skin structure infections (cSSSIs). This study demonstrates that the oral bioavailability of tedizolid is excellent in patients with CF and that the plasma pharmacokinetics are similar to those reported for patients with cSSSIs.

Pharmacokinetic-Pharmacodynamic Target Attainment Analyses To Determine Optimal Dosing of Ceftazidime-Avibactam for the Treatment of Acute Pulmonary Exacerbations in Patients with Cystic Fibrosis.

Acute pulmonary exacerbations (APE) involving Pseudomonas aeruginosa are associated with increased morbidity and mortality in cystic fibrosis (CF) patients. Drug resistance is a significant challenge to treatment. Ceftazidime-avibactam (CZA) demonstrates excellent in vitro activity against isolates recovered from CF patients, including drug-resistant strains. Altered pharmacokinetics (PK) of several beta-lactam antibiotics have been reported in CF patients. Therefore, this study sought to characterize the PK of CZA and perform target attainment analyses to determine the optimal treatment regimen. The PK of CZA in 12 adult CF patients administered 3 intravenous doses of 2.5 g every 8 h infused over 2 h were determined. Population modeling utilized the maximum likelihood expectation method. Monte Carlo simulations determined the probability of target attainment (PTA). An exposure target consisting of the cumulative percentage of a 24-h period that the free drug concentration exceeds the MIC under steady-state pharmacokinetic conditions (fT>MIC) was evaluated for ceftazidime (CAZ), and an exposure target consisting of the cumulative percentage of a 24-h period that the free drug concentration exceeds a 1-mg/liter threshold concentration (fT>1 mg/liter) was evaluated for avibactam (AVI). Published CAZ and CZA MIC distributions were incorporated to evaluate cumulative response probabilities. CAZ and AVI were best described by one-compartment models. The values of total body clearance (CL; CAZ CL, 7.53 ± 1.28 liters/h; AVI CL, 12.30 ± 1.96 liters/h) and volume of distribution (V; CAZ V, 18.80 ± 6.54 liters; AVI V, 25.30 ± 4.43 liters) were broadly similar to published values for healthy adults. CZA achieved a PTA (fT>MIC, 50%) of >0.9 for MICs of ≤16 mg/liter. The overall likelihood of a treatment response was 0.82 for CZA, whereas it was 0.42 for CAZ. These data demonstrate improved pharmacodynamics of CZA in comparison with those of CAZ and provide guidance on the optimal dosing of CZA for future studies. (This study has been registered at ClinicalTrials.gov under registration no. NCT02504827.).