Proof-of-concept for incorporating mechanistic insights from multi-omics analyses of polymyxin B in combination with chloramphenicol against Klebsiella pneumoniae

Carbapenemase-resistant Klebsiella pneumoniae (KP) resistant to multiple antibiotic classes necessitates optimized combination therapy. Our objective is to build a workflow leveraging omics and bacterial count data to identify antibiotic mechanisms that can be used to design and optimize combination regimens. For pharmacodynamic (PD) analysis, previously published static time-kill studies (J Antimicrob Chemother 70, 2015, 2589) were used with polymyxin B (PMB) and chloramphenicol (CHL) mono and combination therapy against three KP clinical isolates over 24 h. A mechanism-based model (MBM) was developed using time-kill data in S-ADAPT describing PMB-CHL PD activity against each isolate. Previously published results of PMB (1 mg/L continuous infusion) and CHL (C-max: 8 mg/L; bolus q6h) mono and combination regimens were evaluated using an in vitro one-compartment dynamic infection model against a KP clinical isolate (10(8) CFU/ml inoculum) over 24 h to obtain bacterial samples for multi-omics analyses. The differentially expressed genes and metabolites in these bacterial samples served as input to develop a partial least squares regression (PLSR) in R that links PD responses with the multi-omics responses via a multi-omics pathway analysis. PMB efficacy was increased when combined with CHL, and the MBM described the observed PD well for all strains. The PLSR consisted of 29 omics inputs and predicted MBM PD response (R-2 = 0.946). Our analysis found that CHL downregulated metabolites and genes pertinent to lipid A, hence limiting the emergence of PMB resistance. Our workflow linked insights from analysis of multi-omics data with MBM to identify biological mechanisms explaining observed PD activity in combination therapy.