A Population WB-PBPK Model of Colistin and its Prodrug CMS in Pigs: Focus on the Renal Distribution and Excretion

The pharmacokinetic-nephrotoxicity relationships of CMS and CMS-E2 from the perspective of plasma and kidney drug concentrations in rats

Nephrotoxicity has seriously affected the clinical application of colistin methanesulphonate (CMS). Colistin B methanesulphonate (CMS-E2) is a novel polymyxin developed and aimed to have lower nephrotoxicity. This study aimed to investigate the relationships between pharmacokinetics (PK) and nephrotoxicity of CMS and CMS-E2 and compare the toxicity of the two drugs in rats. Rats were treated intraperitoneally with a single dose of saline, CMS [10, 20 mg/kg of colistin base activity (CBA)], and CMS-E2 (20, 40 mg/kg CBA). An LC-MS/MS method was developed to determine plasma and renal tissue concentrations of CMS/CMS-E2 and colistin/colistin B. The severity of renal injuries was examined both biochemically and histologically. The PK-toxicodynamic (TD) model was evaluated to characterize the PK of CMS/CMS-E2 and colistin/colistin B in plasma as well as its relationship with nephrotoxicity. Creatinine (CR) and blood urea nitrogen (BUN) profiles were described using an indirect link PK-TD model, with linear-effect relationship. Both the slope between colistin or colistin B concentrations in the effect compartment and CR, BUN was significantly lower for CMS-E2 compared with CMS (CR: P = 0.027, BUN: P = 0.043). The concentrations of colistin and colistin B in kidneys were correlated with CR, BUN values, and histologic examination scores. The regression coefficient of CMS-E2 between the colistin B concentrations in renal tissues and CR, BUN values were lower, as well (CR: P = 0.003, BUN: P = 0.001). The renal injuries induced by CMS and CMS-E2 lagged behind the change of plasma colistin or colistin B concentrations and correlated to those in kidneys. CMS-E2 showed significantly lower nephrotoxicity compared to CMS in vivo.

Chemical risk assessment in food animals via physiologically based pharmacokinetic modeling - Part I: Veterinary drugs on human food safety assessment

Veterinary drugs and environmental pollutants can enter food animals and remain as residues in food chains threatening human food safety and health. Performing health risk and food safety assessments to derive safety levels of these xenobiotics can protect human health. Physiologically based pharmacokinetic (PBPK) modeling is a mathematical tool to quantitatively describe chemical disposition in humans and animals informing human food safety and health risk assessments. However, few reviews focus on the application of PBPK models in food animals and discuss their relationship to human food safety and health risk assessments in the last five years (2020-2024). In this series of reviews, we introduce the methodology, recent progress and challenges of PBPK modeling in food animals. The present review is Part I of this series of reviews and it focuses on applications of PBPK models of veterinary drugs in food animals, whereas Part II is a companion review focusing on environmental chemicals. Advanced strategies of PBPK modeling in risk and food safety assessment, including population PBPK, interactive PBPK web dashboard, and generic PBPK are also summarized in Part I. Additionally, we share our perspective on the existing challenges and future direction for PBPK modeling of veterinary medicines in food animals.

Quantifying combined effects of colistin and ciprofloxacin against Escherichia coli in an in silico pharmacokinetic-pharmacodynamic model

Co-administering a low dose of colistin (CST) with ciprofloxacin (CIP) may improve the antibacterial effect against resistant Escherichia coli, offering an acceptable benefit-risk balance. This study aimed to quantify the interaction between ciprofloxacin and colistin in an in silico pharmacokinetic-pharmacodynamic model from in vitro static time-kill experiments (using strains with minimum inhibitory concentrations, MICCIP 0.023-1 mg/L and MICCST 0.5-0.75 mg/L). It was also sought to demonstrate an approach of simulating concentrations at the site of infection with population pharmacokinetic and whole-body physiologically based pharmacokinetic models to explore the clinical value of the combination when facing more resistant strains (using extrapolated strains with lower susceptibility). The combined effect in the final model was described as the sum of individual drug effects with a change in drug potency: for ciprofloxacin, concentration at half maximum killing rate (EC50) in combination was 160% of the EC50 in monodrug experiments, while for colistin, the change in EC50 was strain-dependent from 54.1% to 119%. The benefit of co-administrating a lower-than-commonly-administrated colistin dose with ciprofloxacin in terms of drug effect in comparison to either monotherapy was predicted in simulated bloodstream infections and pyelonephritis. The study illustrates the value of pharmacokinetic-pharmacodynamic modelling and simulation in streamlining rational development of antibiotic combinations.

A physiologically based pharmacokinetic model to optimize the dosage regimen and withdrawal time of cefquinome in pigs

Cefquinome is widely used to treat respiratory tract diseases of swine. While extra-label dosages of cefquinome could improve clinical efficacy, they might lead to excessively high residues in animal-derived food. In this study, a physiologically based pharmacokinetic (PBPK) model was calibrated based on the published data and a microdialysis experiment to assess the dosage efficiency and food safety. For the microdialysis experiment, in vitro/in vivo relative recovery and concentration-time curves of cefquinome in the lung interstitium were investigated. This PBPK model is available to predict the drug concentrations in the muscle, kidney, liver, plasma, and lung interstitial fluid. Concentration-time curves of 1000 virtual animals in different tissues were simulated by applying sensitivity and Monte Carlo analyses. By integrating pharmacokinetic/pharmacodynamic target parameters, cefquinome delivered at 3-5 mg/kg twice daily is advised for the effective control of respiratory tract infections of nursery pig, which the bodyweight is around 25 kg. Based on the predicted cefquinome concentrations in edible tissues, the withdrawal interval is 2 and 3 days for label and the extra-label doses, respectively. This study provides a useful tool to optimize the dosage regimen of cefquinome against respiratory tract infections and predicts the concentration of cefquinome residues in edible tissues. This information would be helpful to improve the food safety and guide rational drug usage.

Unravelling the Hepatic Elimination Mechanisms of Colistin

PURPOSE

Colistin is an antibiotic which is increasingly used as a last-resort therapy in critically-ill patients with multidrug resistant Gram-negative infections. The purpose of this study was to evaluate the mechanisms underlying colistin's pharmacokinetic (PK) behavior and to characterize its hepatic metabolism.

METHODS

In vitro incubations were performed using colistin sulfate with rat liver microsomes (RLM) and with rat and human hepatocytes (RH and HH) in suspension. The uptake of colistin in RH/HH and thefraction of unbound colistin in HH (fu,hep) was determined. In vitro to in vivo extrapolation (IVIVE) was employed to predict the hepatic clearance (CLh) of colistin.

RESULTS

Slow metabolism was detected in RH/HH, with intrinsic clearance (CLint) values of 9.34± 0.50 and 3.25 ± 0.27 mL/min/kg, respectively. Assuming the well-stirred model for hepatic drug elimination, the predicted rat CLh was 3.64± 0.22 mL/min/kg which could explain almost 70% of the reported non-renal in vivo clearance. The predicted human CLh was 91.5 ± 8.83 mL/min, which was within two-fold of the reported plasma clearance in healthy volunteers. When colistin was incubated together with the multidrug resistance-associated protein (MRP/Mrp) inhibitor benzbromarone, the intracellular accumulation of colistin in RH/HH increased significantly.

CONCLUSION

These findings indicate the major role of hepatic metabolism in the non-renal clearance of colistin, while MRP/Mrp-mediated efflux is involved in the hepatic disposition of colistin. Our data provide detailed quantitative insights into the hereto unknown mechanisms responsible for non-renal elimination of colistin.

Prediction of Tissue Exposures of Meropenem, Colistin, and Sulbactam in Pediatrics Using Physiologically Based Pharmacokinetic Modeling

BACKGROUND

The combination of polymyxins, meropenem, and sulbactam demonstrated efficacy against multi-drug-resistant bacillus Acinetobacter baumannii. These three antibiotics are commonly used against major blood, skin, lung, and heart muscle infections.

OBJECTIVE

The objective of this study was to predict drug disposition and extrapolate the efficacy in these tissues using a physiologically based pharmacokinetic modeling approach that linked drug exposures to their target pharmacodynamic indices associated with antimicrobial activities against A. baumannii.

METHODS

An adult physiologically based pharmacokinetic model was developed for meropenem, colistin, and sulbactam and scaled to pediatrics accounting for both renal and non-renal clearances. The model reliability was evaluated by comparing simulated plasma and tissue drug exposures to observed data. Target pharmacodynamic indices were used to evaluate whether pediatric and adult dosing regimens provided sufficient coverage.

RESULTS

The modeled plasma drug exposures in adults and pediatric patients were consistent with reported literature data. The mean fold errors for meropenem, colistin, and sulbactam were in the range of 0.710-1.37, 0.981-1.47, and 0.647-1.39, respectively. Simulated exposures in the blood, skin, lung, and heart were consistent with reported penetration rates. In a virtual pediatric population aged from 2 to < 18 years, the interpretive breakpoints were achieved in 85-90% of subjects for their targeted pharmacodynamic indices after administration of pediatric dosing regimens consisting of 30 mg/kg of meropenem, and 40 mg/kg of sulbactam three times daily as a 3-h or continuous infusion and 5 mg/kg/day of colistin base activity.

CONCLUSIONS

The physiologically based pharmacokinetic modeling supports pediatric dosing regimens of meropenem/colistin/sulbactam in a co-administration setting against infections in the blood, lung, skin, and heart tissues due to A. baumannii.

Acute kidney injury in cystic fibrosis patients treated with intravenous colistimethate sodium or tobramycin

OBJECTIVES

Colistimethate sodium and tobramycin are important systemic antibiotics for treatment of cystic fibrosis (CF) pulmonary exacerbations but can induce acute kidney injury (AKI). We characterize the rate of AKI in CF patients treated with systemic colistimethate sodium compared with tobramycin.

METHODS

This single-centre, retrospective cohort study included hospitalized CF patients treated with IV colistimethate sodium or tobramycin. The primary outcome was AKI defined using the RIFLE criteria. Multivariate logistic regression using a mixed model was performed to identify variables that were independently associated with AKI.

RESULTS

Overall, 156 patients representing 507 care encounters were included. The OR of AKI was not increased with IV colistimethate sodium relative to IV tobramycin after adjusting for other potential predictor variables (aOR 1.00; 95% CI 0.16-6.03). The frequency of AKI was 9.5% across all encounters, 6.9% with IV colistimethate sodium and 9.9% with IV tobramycin, with RIFLE category R (risk) being the most common stage, accounting for 4.2% of encounters with IV colistimethate sodium and 9.2% with IV tobramycin. The concomitant use of another nephrotoxin (aOR 2.51; 95% CI 1.27-4.95) or the combination of vancomycin and piperacillin/tazobactam (aOR 5.95; 95% CI 2.05-17.3) were both associated with increased odds of AKI.

CONCLUSIONS

Systemic treatment with colistimethate sodium or tobramycin in the CF patient population is associated with a similar rate of nephrotoxicity. However, clinicians should be mindful of the increased risk for AKI in patients treated with either IV colistimethate sodium or IV tobramycin when used concurrently with other nephrotoxic agents, particularly the combination of vancomycin and piperacillin/tazobactam.

Application of the Adverse Outcome Pathway Concept to In Vitro Nephrotoxicity Assessment: Kidney Injury due to Receptor-Mediated Endocytosis and Lysosomal Overload as a Case Study

Application of adverse outcome pathways (AOP) and integration of quantitative in vitro to in vivo extrapolation (QIVIVE) may support the paradigm shift in toxicity testing to move from apical endpoints in test animals to more mechanism-based in vitro assays. Here, we developed an AOP of proximal tubule injury linking a molecular initiating event (MIE) to a cascade of key events (KEs) leading to lysosomal overload and ultimately to cell death. This AOP was used as a case study to adopt the AOP concept for systemic toxicity testing and risk assessment based on in vitro data. In this AOP, nephrotoxicity is thought to result from receptor-mediated endocytosis (MIE) of the chemical stressor, disturbance of lysosomal function (KE1), and lysosomal disruption (KE2) associated with release of reactive oxygen species and cytotoxic lysosomal enzymes that induce cell death (KE3). Based on this mechanistic framework, in vitro readouts reflecting each KE were identified. Utilizing polymyxin antibiotics as chemical stressors for this AOP, the dose-response for each in vitro endpoint was recorded in proximal tubule cells from rat (NRK-52E) and human (RPTEC/TERT1) in order to (1) experimentally support the sequence of key events (KEs), to (2) establish quantitative relationships between KEs as a basis for prediction of downstream KEs based on in vitro data reflecting early KEs and to (3) derive suitable in vitro points of departure for human risk assessment. Time-resolved analysis was used to support the temporal sequence of events within this AOP. Quantitative response-response relationships between KEs established from in vitro data on polymyxin B were successfully used to predict in vitro toxicity of other polymyxin derivatives. Finally, a physiologically based kinetic (PBK) model was utilized to transform in vitro effect concentrations to a human equivalent dose for polymyxin B. The predicted in vivo effective doses were in the range of therapeutic doses known to be associated with a risk for nephrotoxicity. Taken together, these data provide proof-of-concept for the feasibility of in vitro based risk assessment through integration of mechanistic endpoints and reverse toxicokinetic modelling.

Nebulized Colistin in Ventilator-Associated Pneumonia and Tracheobronchitis: Historical Background, Pharmacokinetics and Perspectives

Clinical evidence suggests that nebulized colistimethate sodium (CMS) has benefits for treating lower respiratory tract infections caused by multidrug-resistant Gram-negative bacteria (GNB). Colistin is positively charged, while CMS is negatively charged, and both have a high molecular mass and are hydrophilic. These physico-chemical characteristics impair crossing of the alveolo-capillary membrane but enable the disruption of the bacterial wall of GNB and the aggregation of the circulating lipopolysaccharide. Intravenous CMS is rapidly cleared by glomerular filtration and tubular excretion, and 20-25% is spontaneously hydrolyzed to colistin. Urine colistin is substantially reabsorbed by tubular cells and eliminated by biliary excretion. Colistin is a concentration-dependent antibiotic with post-antibiotic and inoculum effects. As CMS conversion to colistin is slower than its renal clearance, intravenous administration can lead to low plasma and lung colistin concentrations that risk treatment failure. Following nebulization of high doses, colistin (200,000 international units/24h) lung tissue concentrations are > five times minimum inhibitory concentration (MIC) of GNB in regions with multiple foci of bronchopneumonia and in the range of MIC breakpoints in regions with confluent pneumonia. Future research should include: (1) experimental studies using lung microdialysis to assess the PK/PD in the interstitial fluid of the lung following nebulization of high doses of colistin; (2) superiority multicenter randomized controlled trials comparing nebulized and intravenous CMS in patients with pandrug-resistant GNB ventilator-associated pneumonia and ventilator-associated tracheobronchitis; (3) non-inferiority multicenter randomized controlled trials comparing nebulized CMS to intravenous new cephalosporines/ß-lactamase inhibitors in patients with extensive drug-resistant GNB ventilator-associated pneumonia and ventilator-associated tracheobronchitis.

Population pharmacokinetics of colistin and the relation to survival in critically ill patients infected with colistin susceptible and carbapenem-resistant bacteria

OBJECTIVES

The aim was to analyse the population pharmacokinetics of colistin and to explore the relationship between colistin exposure and time to death.

METHODS

Patients included in the AIDA randomized controlled trial were treated with colistin for severe infections caused by carbapenem-resistant Gram-negative bacteria. All subjects received a 9 million units (MU) loading dose, followed by a 4.5 MU twice daily maintenance dose, with dose reduction if creatinine clearance (CrCL) < 50 mL/min. Individual colistin exposures were estimated from the developed population pharmacokinetic model and an optimized two-sample per patient sampling design. Time to death was evaluated in a parametric survival analysis.

RESULTS

Out of 406 randomized patients, 349 contributed pharmacokinetic data. The median (90% range) colistin plasma concentration was 0.44 (0.14-1.59) mg/L at 15 minutes after the end of first infusion. In samples drawn 10 hr after a maintenance dose, concentrations were >2 mg/L in 94% (195/208) and 44% (38/87) of patients with CrCL ≤120 mL/min, and >120 mL/min, respectively. Colistin methanesulfonate sodium (CMS) and colistin clearances were strongly dependent on CrCL. High colistin exposure to MIC ratio was associated with increased hazard of death in the multivariate analysis (adjusted hazard ratio (95% CI): 1.07 (1.03-1.12)). Other significant predictors included SOFA score at baseline (HR 1.24 (1.19-1.30) per score increase), age and Acinetobacter or Pseudomonas as index pathogen.

DISCUSSION

The population pharmacokinetic model predicted that >90% of the patients had colistin concentrations >2 mg/L at steady state, but only 66% at 4 hr after start of treatment. High colistin exposure was associated with poor kidney function, and was not related to a prolonged survival.

A Whole-Body Physiologically Based Pharmacokinetic Model for Colistin and Colistin Methanesulfonate in Rat

Colistin is a polymyxin antibiotic used to treat patients infected with multidrug-resistant Gram-negative bacteria (MDR-GNB). The objective of this work was to develop a whole-body physiologically based pharmacokinetic (WB-PBPK) model to predict tissue distribution of colistin in rat. The distribution of a drug in a tissue is commonly characterized by its tissue-to-plasma partition coefficient, K_p . Colistin and its prodrug, colistin methanesulfonate (CMS) K_p priors, were measured experimentally from rat tissue homogenates or predicted in silico. The PK parameters of both compounds were estimated fitting in vivo their plasma concentration-time profiles from six rats receiving an i.v. bolus of CMS. The variability in the data was quantified by applying a nonlinear mixed effect (NLME) modelling approach. A WB-PBPK model was developed assuming a well-stirred and perfusion-limited distribution in tissue compartments. Prior information on tissue distribution of colistin and CMS was investigated following three scenarios: K_p was estimated using in silico K_p priors (I) or K_p was estimated using experimental K_p priors (II) or K_p was fixed to the experimental values (III). The WB-PBPK model best described colistin and CMS plasma concentration-time profiles in scenario II. Colistin-predicted concentrations in kidneys in scenario II were higher than in other tissues, which was consistent with its large experimental K_p prior. This might be explained by a high affinity of colistin for renal parenchyma and active reabsorption into the proximal tubular cells. In contrast, renal accumulation of colistin was not predicted in scenario I. Colistin and CMS clearance estimates were in agreement with published values. The developed model suggests using experimental priors over in silico K_p priors for kidneys to provide a better prediction of colistin renal distribution. Such models might serve in drug development for interspecies scaling and investigate the impact of disease state on colistin disposition.