- Research Article
- Open Access
A scalable metabolite supplementation strategy against antibiotic resistant pathogen Chromobacterium violaceum induced by NAD+/NADH+ imbalance
BMC Systems Biology volume 11, Article number: 51 (2017)
The leading edge of the global problem of antibiotic resistance necessitates novel therapeutic strategies. This study develops a novel systems biology driven approach for killing antibiotic resistant pathogens using benign metabolites.
Controlled laboratory evolutions established chloramphenicol and streptomycin resistant pathogens of Chromobacterium. These resistant pathogens showed higher growth rates and required higher lethal doses of antibiotic. Growth and viability testing identified malate, maleate, succinate, pyruvate and oxoadipate as resensitising agents for antibiotic therapy. Resistant genes were catalogued through whole genome sequencing. Intracellular metabolomic profiling identified violacein as a potential biomarker for resistance. The temporal variance of metabolites captured the linearized dynamics around the steady state and correlated to growth rate. A constraints-based flux balance model of the core metabolism was used to predict the metabolic basis of antibiotic susceptibility and resistance.
The model predicts electron imbalance and skewed NAD/NADH ratios as a result of antibiotics – chloramphenicol and streptomycin. The resistant pathogen rewired its metabolic networks to compensate for disruption of redox homeostasis. We foresee the utility of such scalable workflows in identifying metabolites for clinical isolates as inevitable solutions to mitigate antibiotic resistance.
A post-antibiotic apocalypse portends the utility of antibiotics being massively compromised through evolution of antibiotic resistance. The genetic basis of several antibiotic resistant populations have been delineated via whole genome sequencing to the level of identifying resistance genes [1,2,3,4]. Antibiotic resistance genes catalogued through these efforts have rarely provided individualized therapies. Penicillin resistance in Streptococcus pneumoniae is a consequence of mutations in putative iron transport systems . Stress pathways contribute to glycopeptide resistance in Staphylococcus aureus . The evolution of antibiotic resistance in pathogens is characterized by uncontrolled proliferation even in the presence of drugs. Growth and energy generation are two principal dimensions of cell function and proliferation. This duality of cell function, orchestrated by metabolic networks, is critical for survival and governs resistance. Redox homeostasis is important to effectively harness reducing power produced through the catabolism of various substrates and to utilize this power in the anabolism of cellular components such as DNA, lipids and proteins. Metabolic regulation and gene expression modulation are now recognized as major players in antibiotic resistance [6,7,8,9]. β-Lactamases are amongst the most common causes of bacterial resistance to β-Lactam antimicrobial agents . Derepression of secondary metabolism in Nocardiopsis, identified through metabolome characterization, was a consequence of acquired resistance . In E. coli it was reported that overflow metabolism and reactive oxygen species (ROS) formation are inherent cellular response to antibiotic lethality . A recent study implicates accelerating cellular respiration rates in the bactericidal mode of action . Specific metabolites have been associated with varying degrees of killing antibiotic tolerant pathogens (persisters) by stimulating proton motive force (PMF) and increasing antibiotic uptake [12, 13]. Promoting Tricarboxylic acid cycle (TCA)/Krebs cycle by glucose/alanine activation that subsequently also increase PMF stimulating uptake of antibiotic have initiated death in multi-drug resistant Edwardisiella tarda . Stoichiometric flux balance models based on evolutionary optimality predict outcomes of single environment evolution accurately and can compute operational flux states, growth and energy phenotypes of resistant/susceptible cells [14,15,16,17,18]. Investigating altered metabolism and connecting to evolved resistant genotypes may provide simple strategies to overcome drug resistance and induce susceptibility to existing antibiotics.
In this study, we have identified benign metabolites to stimulate antibiotic action and death of streptomycin and chloramphenicol resistant pathogen Chromobacterium violaceum. Primarily a zoonotic pathogen, it is opportunistic in humans and converts the essential amino acid tryptophan to violacein, a blue-violet pigment . The oxidative dimerization of two tryptophan molecules is an essential and regulatory step for the synthesis of violacein scaffold [20,21,22]. Chromobacterium violaceum is also reported sensitive to aminoglycosides, chloramphenicol, and tetracycline and resistant to ampicillin, penicillin, and first-generation cephalosporins [23, 24]. Our work highlights for the first time that differential violacein levels could act as a potential biomarker for resistance to two distinct classes of antibiotics. Here we pioneer a scalable workflow (Fig. 1) from controlled evolution of resistance to rationally identifying metabolites that induce drug susceptibility using systems biology approaches. Constraints-based flux balance modeling, being an integral part of systems biology approach, was used for the first time to predict that disruption of redox homeostasis was causal for antibiotic action. Compensatory metabolic reprogramming to overcome redox cofactor imbalance was delineated. In addition the in silico resistant growth phenotype was predicted to be a function of the rigidity of flux network.
Adaptive Laboratory Evolution (ALE)
The detailed ALE workflow is well described in Fig. 1a. Chromobacterium violaceum strain ATCC 12472T, wild type, (C. violaceum or WT) was obtained from the American Type Culture Collection Center (ATCC), USA. The wild type C. violaceum was routinely cultured on Luria-Broth (LB, Hi-Media-M575) and maintained at 30 °C with continuous aeration in a shaker incubator set at 180 revolutions per minute (rpm). C. violaceum was tested to be susceptible to low concentrations of both antibiotics (Chloramphenicol: MIC 8 μg/mL and Streptomycin: MIC 10 μg/mL in liquid culture and 60 μg/mL in agar plates). Antibiotic resistant strains of C. violaceum were evolved separately under controlled laboratory environments using the two antibiotics, chloramphenicol (chl) and streptomycin (str) at sub-lethal concentrations (10 μg/mL) on Luria Bertani agar (LBA) plates. Clonal purification by repeated sub-culturing of the colonies obtained on LBA plates with antibiotic (10 μg/mL) resulted in single colonies. These colonies were cultured in LB with antibiotic (10 μg/mL) and further cryopreserved in 50% glycerol and all genotyping and phenotyping experiments done by thawing the frozen vials and sub-culturing in LB with antibiotic (10 μg/mL) at 30 °C with continuous aeration in a shaker incubator set at 180 rpm until mentioned otherwise (Fig. 1d).
Minimum Inhibitory Concentration (MIC) determination
Antibiotic susceptibilities were determined with EzyMICTM Strips (HiMedia Laboratories, India) on Müller-Hinton agar plates using the manufacturer’s instructions (Fig. 2c and d). The MICs were further determined using broth micro-dilution method for C. violaceum (at 30 °C, 180 rpm) according to the dilution method described previously  (Fig. 2e).
Effect of varying concentration of antibiotics on growth profiles and growth rate estimation
All strains were profiled for growth by varying concentrations of antibiotics at 30 °C, 180 rpm. Exponential-phase cultures were prepared at 30 °C, 180 rpm using a shaker incubator and further used to inoculate 3 mL of LB to an initial OD600 of 0.1. Antibiotic stock solutions were added to yield desired concentrations of the antibiotic ranging from 0 to 256 μg/mL. Cultures were incubated in a shaker incubator at 30 °C, 180 rpm and bacterial cell densities were estimated hourly using a spectrophotometer. Growth profile assays for each C. violaceum strain were performed in triplicate using independent starter cultures and antibiotic stocks. Growth rate was estimated graphically from growth curves by plotting the natural log values of OD600 for each time point and determining the slope by linear regression. A minimum of four time points were used to determine the growth rate (Fig. 2a and b).
Growth profiling on different exogenous carbon and nitrogen sources
C. violaceum was cultured overnight, then diluted to a density of 0.002 (Optical density at 600 nm, OD600), mixed with an equal volume of LB that contains different metabolites at a final concentration of 2 mg/mL, except for lactate which was at a concentration of 0.27 mg/mL. C. violaceum was cultured in 96-well flat-bottom plastic microplates at 30 °C for 30 h with or without antibiotic pressure. 30 different metabolites were tested in biological triplicates. The plates were monitored for growth, biomass and violacein using iMark™ Microplate Absorbance Reader (BIO-RAD) at 550 nm and 655 nm at regular time intervals. The amount of violacein and dry cell weight were determined using 550 nm and 655 nm readings and compared to standard calibration graphs (Additional file 1: Figure S4) for quantitation purpose. Different conditions were tested wherein one set had the respective antibiotics, to which the strains were resistant, from zeroth time point (t0) and in the other set, antibiotic was added 6 h (t6) post inoculation. After 30 h the t6 set of plates were used to plate fresh LBA plates without any antibiotic to analyze viable colonies after the 30 h incubation period. Growth rates were measured for the overall 30 h duration experiment in four different phases: pre-6 h phase, post 6 h phase, overall growth rate and a maximum growth rate. Same phase wise analysis for biomass and violacein was also performed. Curve fitting, visualization and analysis of the different plots for this experiment were done using GraphPad Prism Version 6.01 (GraphPad Software, San Diego California USA, www.graphpad.com). Nonlinear curve fitting was done using Gompertz growth equation  for the growth data points (Additional files 2, 3 and 4). All the heat maps were generated using MATLAB platform (Fig. 3).
Whole genome sequencing
In order to extract genomic DNA for whole genome re-sequencing, ChlR and StrpR cultures were revived from previously cryopreserved glycerol stocks on LBA plates with respective antibiotic at 30 °C. These cultures have been previously tested for all the phenotypic traits as described in the results section. A single colony from LBA was cultured in LB broth at 30 °C, 180 rpm and mid – log phase cells were harvested for genomic DNA extraction using the DNeasy Blood and Tissue Kit (Qiagen, USA) according to the manufacturer’s instructions. The quality of the genomic DNA was assessed for RNA contamination using A260/A280 ratio and visualized on agarose gel. The DNA was also quantitated using Qubit before library preparation. Genome sequences were obtained for the two evolved populations, ChlR and StrpR, using Ion Torrent PGMTM (Life Technologies) NGS using the 314TM chip with mean read length of 180 and 188 base pairs respectively. The sequencing data analysis for ChlR and StrpR samples showed 86 and 84% of the bases read were of ≥ Q20 quality respectively. For 77.07% of ChlR and 71.3% of StrpR samples the genome base coverage was 20x and for 18.61 and 10.01% coverage was 100x. The assemblies and comparative analysis against NCBI sequence of C. violaceum (accession number NC_005085) were also performed. The genome sequence data for this publication has been deposited at Sequence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra) under accession number SRP072862.
Confirmation of mutations by Sanger Sequencing
All NGS identified sequence variations were confirmed by Sanger sequencing. Primers (Additional file 1: Table S5) were designed to amplify around 200 to 600 base pairs (bp) amplicons such that the nucleotide of interest was at a position for easy read during Sanger sequencing. Amplification and sequencing was performed by Eurofins Genomics India Pvt. Ltd.
In silico Structure - Function analysis of the mutation acquired
Ab-initio models were made using ROBETTA server  (http://robetta.bakerlab.org) except for RpsL (homology modeling was used) and the models generated were checked for various parameters for model assessment, such as Ramachandran Plots using PROCHECK. Finally, visualization and manipulation of the three dimensional models were performed using the software PyMOL (Schrödinger, LLC, 2015). In addition 3DLigandSite  was used to get a better understanding of the structure-function change post mutation/variation (Fig. 4).
Preparation of intracellular metabolite extracts
For performing metabolomics experiments, the three different populations of C. violaceum were inoculated (10% inoculum of overnight starter culture) and incubated in a shaker incubator at 30 °C, 180 rpm for 30 h. 2 mL of cell cultures were harvested at the end of 0, 6, 12, 18, 24 and 30 h by centrifugation at 12000 g at 4 °C for removal of extracellular media. The following steps were carried out on ice. The pellets were reconstituted with ice-cold ethanol for quenching as well as maximal extraction of metabolite features as used in previous study [29, 30]. For collecting the extracts, the suspension was centrifuged at 14000 g at 4 °C. The intracellular extracts were aliquoted (100 μL) and immediately stored at −80 °C till further use.
The intracellular as well as extracellular extracts were used for violacein estimation with BioPhotometer (Eppendorf) by comparing to standard calibration graphs for quantitation purpose as previously described .
LC-MS analysis of intracellular metabolites
The intracellular extracts were dried in centrivap (Labconco) at 4 °C, followed by reconstitution in 100 μL of 10% water in acetonitrile containing 4 μM atorvastatin (internal standard). 10 μL of each sample was pooled to prepare a technical quality control (QC) sample. Metabolic profiling of samples was carried out on High Performance Liquid Chromatography – Heated Electrospray Ion Source – High-resolution Mass spectrometry (HPLC-HESI-HRMS). The separation was achieved on Sequent (Zwitter ionic hydrophilic interaction Liquid Chromatography) ZIC-HILIC column (100 mm*2.1 mm*5 μm, Merck Millipore) column using HPLC consisting of Accela quaternary gradient pump, a degasser and Accela autosampler (Thermo Fisher). The column was maintained at 45 °C using column oven (PerkinElmer). The mobile phase for elution consisted of 0.1% formic acid in deionized water (Mobile phase ‘A’) and 0.08% formic acid in acetonitrile (Mobile phase ‘B’). Gradient was set with 5% of mobile phase A (0–5.0 min, 300 μL/min), 13% A (15.0, 300 μl/min), 45% A (20.0 min, 300 μL/min), 90% A (23.0 min-25.0 min, 300 μL/min), 5% A (27.0–32.0 min, 700 μL/min). Heated electrospray ionization (HESI) source was used as an interface between Liquid Chromatography (LC) and High-resolution Mass Spectrometry (HRMS) instruments. The spray voltage of the source was set at 3.7 kV with capillary temperature 300 °C, sheath gas 45 units, auxiliary gas 10 units, heater temperature 390 °C and S-lens RF at 50 units. The data was acquired in range of 70–1050 m/z at resolution of 70,000 Full Width at Half Maximum (FWHM) with Automated Gain Control (AGC) target 1e6 and injection time of 120 ms. Two technical replicates, each of 5 μL sample volumes were injected during analysis in both positive and negative ion mode. A total of 7 QC samples were run at beginning, intermittently and end of the run.
LC-MS Data analysis
For data analysis, the Qual browser module of Xcalibur (Thermo) was used for manual inspection of presence of metabolite of interest through accurate mass - extracted ion chromatogram (AM-XIC). A mass extraction window (MEW) of 20 ppm around monoisotopic m/z of possible adduct was used to generate the XIC. After establishing the retention time (Rt) and peak width of respective metabolites as well as internal standard (IS), Tandem Mass Spectrometry (MS/MS) was carried out in respective ion modes to confirm their identities. The MS/MS data was analyzed using fragment search tool in METLIN  data base (http://metlin.scripps.edu) and Mass frontier 7.0. A data processing method from processing setup module of Xcalibur software was prepared to integrate and generate area under peak data. The peak areas of metabolites were normalized to the area response of uniformly spiked internal standard, atorvastatin (ATV). This ratio is representative of the intracellular metabolite concentration (abundance) after taking into consideration the dilution factor for each sample extract. Principal component analysis (PCA) was used to get an overview of the data and to reduce the high dimension of the data set. For each metabolite, we calculated the temporal variation, using the coefficient of variation (CV):
where σ and μ are the standard deviation and mean of the measurements across the time points respectively.
AP-MALDI of intracellular metabolites
It was also estimated on Thermo Q-Exactive mass spectrometer (MS) coupled with an atmospheric pressure - matrix-assisted laser desorption/ionization (AP-MALDI) source equipped with a solid state Nd:YAG laser operating at 355 nm. A mixture of 2,5-dihydroxybenzoic acid (2,5-DHB) and 2-cyano-4-hydroxycinnamic acid (CHCA) was used as a matrix for the analysis. Samples were mixed with internal standard (2,4-diamino-6-methyl-1,3,5-triazine) before spotting on MALDI target plate in 6 replicates. The instrument was operated in full MS scan mode within m/z 50–750 at resolution 35,000 FWHM. Spray voltage at 2.5 kV, capillary temperature at 250 °C, AGC target of 1e6 and 500 μs injection time were optimized before beginning analysis. Laser fluence was optimized at 70% and Pulse Dynamic Focusing (PDF) value of 100 μs was used with automated rastering motion chosen for data acquisition.
AP-MALDI data analysis
Data analysis was performed with Thermo XCalibur Qual Browser, mMass  for qualification and in-house software MQ v 5.0 (http://www.ldi-ms.com/services/software) for quantification  with a chosen MEW of 20 ppm.
Constraints based modeling of C. violaceum central metabolism: network reconstruction
Stoichiometric network analysis based on the constraint-based modeling framework has been proven to be a valuable tool to study cellular metabolism and phenotypic capabilities of many organisms [35,36,37,38]. -The small-scale central metabolic model of C. violaceum presented here is a manually curated stoichiometric network reconstruction and model that allows probing special characteristics of this bacterium. It was done using available literature data [39,40,41,42] as well as information from databases such as KEGG, Biocyc, Metacyc. Further the violacein biosynthesis reaction list, including reaction stoichiometry, reversibility, sub-cellular localization, and gene locus/loci for each reaction comprising core metabolism were also incorporated into the model. The biomass equation for the model was also modified in order to take into account for tryptophan contribution towards biomass production. For all simulations maximization of the biomass equation was fixed as the objective function until mentioned otherwise. The detailed description of the model can be found in Additional file 5. The model was initially validated using a set of 10 substrate utilization BIOLOG GN2 plate data existing in literature [43,44,45].
Further a set of constraints that define the antibiotic susceptible WT and differentiated the evolved populations (ChlR and StrpR) were determined. The constraints used in different simulations included (i) Substrate (Glucose) uptake rates (ii) Growth yields (iii) Violacein secretion (iv) ATP maintenance costs associated with molar growth yields of each strain as discussed (Table 3, Additional file 5). The specific growth rates were calculated using 1 g biomass as the basis. The growth yields thus calculated were compared across the three strains. The goal of the simulations was to understand the flux distribution in silico and the sensitivity of growth yield to various precursors with specific reference to the cofactors NADH, NADPH and ATP. Constraints based flux balance analysis (FBA), as described in the following section, was used to simulate for growth (maximize biomass objective function) and violacein production. Constraints-based methods were used to perform a comparative analysis between the susceptible and resistant populations to understand the connections between metabolism and resistance.
Implementation of the central metabolic reconstruction for C. violaceum and constraints-based analysis was done using Constraints Based Reconstruction and Analysis (COBRA) Toolbox 2.0.2  with MATLAB v 7.11, (R2010b) and TOMLAB/CPLEX v7.7 optimizer. MATLAB codes for all referenced COBRA functions are available at the COBRA’s website (https://opencobra.github.io/).
Flux Balance Analysis (FBA) and associated sensitivity parameters
Flux-balance analysis is a method for assessing the systemic properties and cell behaviors of a metabolic genotype. In short the primal FBA problem, equation (1) describes the steady-state mass balances of the biochemical reaction network [47,48,49,50]
where c, v, vLB, and vUB are vectors of length n, and S is the m x n stoichiometric matrix.
Mathematically, the S matrix acts as a linear transformation between the vector that defines fluxes through n reactions in the biochemical network and the vector of the time derivatives of the concentrations of m metabolites involved in these reactions. The fundamentals of FBA have been widely and critically reviewed [46, 51,52,53,54].
The function optimizeCbModel(model), in COBRA toolbox was used for maximization of pathogen growth or biomass by fixing the objective function to be the biomass equation in the model.
Two sensitivity parameters – shadow prices and reduced costs  - were assessed in order to understand the effects of changing biomass, ATPM, metabolites and reactions of the different populations of C. violaceum. Shadow price corresponds to the sensitivity of the growth rate as an objective function (Z) in response to a change in the availability of a metabolite (i), and indicates how much an increment in that metabolite will increase or decrease the growth rate. Analogous to shadow price, reduced cost is the sensitivity of the objective function in response to change in fluxes of a particular reaction and its effect on the objective. In addition to the primal solution (optimal fluxes), the LP solver provides the corresponding dual solution i.e., shadow price and reduced cost for the FBA problem.
Scaled reduced costs  and logarithmic sensitivity coefficient  were also calculated as shown in equation (2) and (3) which better assess the sensitivity taking into account the substrate and growth yield. The logarithmic sensitivity coefficient (Di) represents the sensitivity of each precursor yield to its biosynthetic demand whereas scaled reduced costs (Wi) are used to assess the limiting reactions.
In equation (2) vi is the flux through a particular reaction having a wi reduced cost associated with it. Z is the objective function value, in this case being biomass. Similarly, in equation (3) dX/dM is the associated shadow price to a particular metabolite and dM is the coefficient of the metabolite in the objective function equation.
Another sensitivity parameter, gamma redox, the difference between the shadow prices for the redox couplet (NADH/NAD) and (NADPH/NADP) were calculated .
This is an index of available reducing capacity available to the cell and whether it is limiting or in excess for biomass formation. A positive shadow price suggests available reducing capacity in excess of the optimal demand for growth and where it is negative, limiting for optimal growth. A positive shadow price doesn’t necessarily indicate accumulation but also suggests metabolic rewiring into overflow metabolism.
Flux Variability Analysis (FVA)
FVA was utilized to investigate the resulting space of feasible flux distributions . FVA can be set up in COBRA toolbox using the function fluxVariability(). One can thus determine the minimum and maximum flux value that each reaction in the model can take up while satisfying all constraints on the system for a specific objective. The objectives under consideration for this study include optimal growth. These will be considered as forced or fixed fluxes. Differences in these unique forced fixed rates in resistant populations as compared to wild type indicate metabolic reprogramming. To highlight the differences between the antibiotic sensitive and resistant populations, we classified reactions in the network based on their minimum and maximum flux values and assigned categories that reflect their rigidity or flexibility. Nine categories can thus be mapped onto the flux variability map based on the magnitude and direction of the flux (Additional file 1: Table S4a) ranging from category 1 for forward direction (positive fixed) fixed flux, i.e., minimum and maximum flux values were same, a non-zero positive value to category 9 wherein the minimum and maximum flux was zero (blocked). Specific attention was paid to reaction rates that were uniquely determined (i.e., if upper and lower boundaries as computed by FVA coincide; Categories 1 and 4). Changes between such rigid fluxes (1 and 4) to more variable flux capabilities (2,3,5,6,7 and 8) would reprogram the metabolic network by either changing the direction of equilibrium or modulating magnitude and span of the reaction rates.
Dynamic Flux Balance Analysis (dFBA)
dFBA was utilized to qualitatively predict the outcomes of growth in batch culture conditions matching our experimental condition . The resistant populations ChlR and StrpR needed to be assessed for the onset of overflow metabolism and secretion patterns as observed in FVA simulations. The dFBA can be set up in COBRA toolbox using the function dynamicFBA() . It is an implicit iterative process wherein at each iteration; FBA is used to simulate for growth, nutrient uptake and by-product secretion rates using an initial concentration for nutrients, which are in turn used to calculate biomass and nutrient concentrations in the culture at the end of the step. The same values are used to calculate maximum uptake rates of nutrients for the next time step.
Glucose levels of C. violaceum cultures (WT, ChlR and StrpR) were measured using MyQubit Amplex Red Glucose Assay Kit (Thermo Fischer Scientific) according to the manufacturer’s instructions using Qubit 2.0 Fluorometer (Life Technologies, CA, USA). The assay is based on the D-glucose dependent generation of hydrogen peroxide with glucose oxidase (GO) followed by the horse radish peroxidase (HRP) catalyzed oxidation of colorless, stable Amplex Red (10-acetyl-3,7-dihydroxyphenoxazine) to fluorescent resorufin.
NADH and NAD measurements
NAD/NADH levels of C. violaceum cultures (WT, ChlR and StrpR) were measured using a commercially available kit (MAK037, Sigma Chemical, St. Louis, MO, United States) according to the manufacturer’s instructions (Fig. 6c to e). NAD total (NAD and NADH) or NADH levels are quantified in a colorimetric assay at 450 nm using iMark™ Microplate Absorbance Reader (BIO-RAD).
Adaptive Laboratory Evolution (ALE)
The schema for the workflow followed in this study (Fig. 1) involved the evolution of Chromobacterium violaceum strain ATCC 12472T (C. violaceum or WT) from a small inoculum onto Luria Bertani agar (LBA) plates chloramphenicol (chl) and streptomycin (str), both targeting protein translation at the ribosomal subunit level. To provide strong evolutionary pressure while maintaining a sizeable population, the concentration of antibiotic was chosen such that no more than 60% of growth was inhibited. Once defined this sub-lethal concentration (SLC), of the antibiotic (10 μg/mL) was not varied throughout the evolution as well as other experiments unless specified. The adaptive evolution continued for about 3 weeks and when the first positive trait appeared, these colonies were cultured multiple times on LB agar plates with the respective antibiotics followed by colony purification. Multiple clones that were evolved in parallel were colony purified. Broth cultures of the respective antibiotic resistant population were used for Minimum Inhibitory concentrations (MIC) calculations. One of the parallel lines of evolved clones resistant to antibiotics was further used in genotypic, phenotypic and metabolic profiling studies (Fig. 1d).
Evolution of antibiotic resistance and fitness
We first analyzed that growth rates and kinetic profiles of the resistant populations support evolution towards fitness on both antibiotics. The effect of varying antibiotic concentration on growth rate was studied (Fig. 2a and b). The growth rates of the chloramphenicol resistant (ChlR) population was reduced to 50% at 32 μg/ml chloramphenicol, while the streptomycin resistant population (StrpR) growth rate was lowered down to just 15% of that without antibiotic. At even 10 μg/ml of streptomycin the growth rate of the StrpR population was reduced by only 10%. The growth rate exponentially decreased with increasing concentration of the antibiotic in the wild type and evolved resistant populations (ChlR and StrpR). The resistant populations improved in growth rate and biomass yield substantially even in the presence of higher concentrations of antibiotic in contrast to the wild type (WT). Surprisingly, no fitness costs associated with the acquired resistance were observed in the absence of antibiotics.
Next, we assessed broad-spectrum antibiotic susceptibility for ChlR and StrpR populations via estimation of MIC on 11 different antibiotics (Fig. 2c). The population evolved on chloramphenicol, ChlR, had an MIC much greater than the wild type, being able to resist high titres of 256 μg/ml of chloramphenicol. The streptomycin evolved population, StrpR was able to resist twice the amount of antibiotic as the wild type reaching a titre of 120 μg/ml. StrpR populations showed low MICs for chloramphenicol indicating higher sensitivity. Only the ChlR population showed cross-resistance to Piperacillin/Tazobactam (PTZ) combination (Fig. 2d). Similar trends were observed when MIC values were estimated by broth dilution method (Fig. 2e) for chloramphenicol (>256 μg/mL) and streptomycin (120 μg/mL) and were represented as percentage viability of the cells.
Effect of exogenous metabolites and antibiotics on growth
A systematic evaluation of benign microenvironment metabolites in excess of being limiting (Additional file 1: Table S2) showed unique fitness landscapes and associated costs for the evolved and wild type populations (Fig. 3). Wild type C. violaceum (WT) does not show capacity to utilize citrate, oxalate and glyceraldehyde-3-phosphate (Additional files 2, 3 and 4). Streptomycin (bactericidal) showed a more profound effect on growth, unable to support growth on 50% of the substrates tested while chloramphenicol (bacteriostatic) affected growth on only 7/30 (23%) of the substrates (Table 1 and Fig. 3a, c). The ChlR and StrpR populations showed fitness costs associated with growth on 13 and 17 substrates respectively. StrpR populations showed lag for extended period of times on many substrates (consistent with mutations discussed in the next section). Lowered fitness is observed on glycolytic intermediates like fructose-1,6-diphosphate and other TCA cycle intermediates. Strikingly, the ChlR population recorded almost no growth and viability on organic acids maleate (C3), pyruvate (C3), succinate (C4) and 2-oxoadipate (C6) even in the presence of chloramphenicol antibiotic (Fig. 3b and d). The StrpR population exhibited similar growth patterns with the only exception of D-malate (C3) also being able to make the resistant population susceptible again to antibiotic.
Altered genotypes and in silico protein functional analysis in resistance
In order to see how genotype was shaping the growth phenotype of the resistant populations, we re-sequenced the whole genome and confirmed 14 sequence changes (Table 2, Additional file 1: Table S5) using capillary sequencing. The sequence changes were seen to affect the protein structure and function in silico as discussed (Fig. 4).
The ChlR population acquired mutations in marC  and the transcription repressor acrR of the tripartite AcrAB-TolC multidrug efflux pump [60,61,62]. Mutations in marC were silent substitutions. The repressor protein suffered premature truncation after translation of 141 amino acids and a non-synonymous substitution R60L resulting in altered domain that is involved in ligand binding for repression (Fig. 4a and Additional file 1: Figure S5c). In case of the StrpR population, a non-synonymous substitution R86S was detected in rpsL, coding for the 30S ribosomal protein S12 [63,64,65,66] along with deleterious mutations in kdpD, encoding for the sensor kinase of the two-component signal transduction system (TCS)  and pabC, encoding the pyridoxyl 5’ phosphate (PLP)–dependent enzyme 4-amino-4-deoxychorismate lyase catalyzing PABA biosynthesis [68, 69]. Ab initio models built using ROBETTA and ligand binding analysis using 3DLigandSite showed that the mutated KdpD and PabC proteins truncated after 682 and 226 amino acids respectively resulting in loss of critical ligand binding or sensory domains (see Discussion).
The metabolic phenotype of resistant populations
To determine the fundamental genotype-phenotype relationship and establish a metabolic basis for antibiotic resistance, intracellular metabolomic profiles were measured using Liquid Chromatography – High-resolution Mass Spectrometry with Tandem Mass spectrometry (LC-HRMS with MS/MS) (Fig. 5, Additional file 6). Including both ion modes, total 126 metabolites were screened. At the end of the analysis, 59 metabolites were finally selected based upon qualification criteria of mass accuracy, MS/MS confirmation, elution profile, reproducibility of response of technical QC samples and biological relevance. The details of these 59 metabolites are mentioned in Additional file 6. Relative abundance was calculated by normalizing metabolite peak area response with that of internal standard (metabolite of interest peak area/internal standard peak area). This ratio is representative of the intracellular metabolite concentration (abundance) after taking into consideration the dilution factor for each sample extract across the three different populations of C. violaceum. Metabolite relative abundance levels in the three populations span three orders of magnitude intracellularly (Additional file 1: Figure S1). Guanosine, methylmalonate, glutamine and aspartate vary slightly, but high order of magnitude differences are seen in PABA, succinate, leucine, hypoxanthine and violacein. To validate that the metabolism is different in the wild type (WT) and resistant strains (ChlR and StrpR), Principal Component Analysis (PCA) of quantitative features of metabolites of intermediary metabolism extracted from LC-HRMS data was performed (Additional file 1: Figure S2). Score plots of principal components for both biological replicates, show trends that showcase maximum separation of data with respect to different time points (PC1) and also clustering of data points based on susceptibility or resistance to antibiotic (PC3). Differential expression of features exhibited by StrpR was captured by PC3 in both the biological replicates. A significant difference after 6 h time point in scores of StrpR shows a distinct deviation in metabolic behavior in comparison to the other populations studied. This is also observed in the growth profiles of StrpR. The ChlR strain shows separation from the WT and the StrpR populations in terms of intermediary core metabolism on glucose as also identified through flux balance modeling discussed in following sections. The significant metabolites after almost 30 h of growth include deoxyviolacein, xanthine and β-cyanoalanine while the metabolites in the early hours of growth include more core metabolite candidates like fumarate, maleate, malate, succinate and pyruvate (Additional file 1: Figure S2). The metabolite abundances also show oscillatory behavior with varying amplitude, period and phase lag (Fig. 5e). Leucine, lysine and proline had a characteristic oscillatory behavior with a period of 12 h; also there is phase lag in WT compared to resistant populations. Certain metabolites such as arginine and adenosine showed negligible changes. Intermediate of violacein biosynthesis pathway, prodeoxyviolacein was seen to increase linearly only in ChlR whereas it was very low for the other two i.e., WT and StrpR. The presence of prodeoxyviolacein, a precursor only in ChlR (Fig. 5b) potentially explains the lowered violacein through limited availability of cofactor NADPH since tryptophan levels are similar. Increased recycling of nucleotides through salvage pathways was reflected in high levels of adenosine, xanthine and hypoxanthine in the resistant populations (Fig. 5c and e). 8-oxoguanine, a major oxidized base lesion formed by reactive oxygen species, was higher in the StrpR population indicating potential oxygen radical effects.
The linearized dynamics around the steady-state level of metabolites is captured by the temporal variation (Fig. 5d, See Methods). Growth limiting metabolites are known to show relatively lower temporal variation . Temporal variation (TV) identified potential growth limitation by malate, glucose, glyceraldehyde-3-phosphate and uracil across all populations. Phenylalanine and methionine are potentially growth limiting (low TV) in ChlR while pyruvate only in StrpR. Tyrosine and serine are less growth limiting (high TV) in the resistant populations than the wild type. Pyruvate and malate showed low average intracellular concentrations or had low temporal variation in the two resistant populations (Fig. 5c and d).
In silico prediction of NAD/NADH balance and redox homeostasis
To delineate a metabolic basis for the emergent resistance, the heterogeneous components of resistant genomes and elucidated metabolic physiology were integrated into a constraints-based flux balance model. A central metabolic network reconstruction represented by the genotype of C. violaceum iDB149 was developed (Fig. 6a and b; Additional file 5). The network reconstruction represents the core metabolism for the pathogen including glycolysis, pentose phosphate pathway, TCA cycle, electron transport and basic amino acid metabolism. Although, primarily a generic reconstruction of central metabolism, it includes the virulence factor metabolism, tailoring it to mirror C. violaceum metabolism. It includes detailed amino acid metabolism of specifically tryptophan, due to its direct connection with the production of violacein, a virulence factor specific to this pathogen. A biomass equation was defined as a drain on metabolites present in intermediary metabolism and macromolecules in the precursor biomass based on E. coli and Chromobacterium legacy data. This biomass composition is kept constant throughout the analysis and between strains. In order to understand the effects of changing biomass composition, a logarithmic sensitivity coefficient  that represents the sensitivity of each precursor yield to its biosynthetic demand was calculated for each of the 14 precursors of biomass (Additional file 5). Rank ordering identified higher impact of cofactors NADH, NADPH and ATP to increasing growth. Molar growth yields calculated based on experimental glucose utilization data and the maximum growth yield for Chromobacterium (Additional file 5) was used to estimate the growth-rate independent energy ATP maintenance flux  (vATPM) that represents the energy required to sustain basal cellular activities. The values interpolated/fitted from growth yield curves differed for the wild type and the resistant populations due to varying molar growth yields (Additional file 5, Table 3). The differential violacein phenotype (represented as a production/secretion rate constraint) calculated from experimental data was used to define resistant populations in silico (Additional file 1: Table S3). Based on the sensitivity analysis of violacein production and growth yields, a trade-off exists between the production of violacein and biomass production. Fixing this biosynthetic demand as a critical constraint in the model, growth rates (via growth yields) predicted for both resistant ChlR, StrpR and WT populations were consistent with experimental data. The experimental rates used for the simulations were from the exponential growth phase of the three different populations of C. violaceum. With these constraints determined by experiments, the wild type model was tested for prediction of carbon source utilization patterns using experimental legacy data [43,44,45] of Biolog™. The model being a core metabolic model was validated for utilization of carbon sources that had transporters included and growth against a small subset of substrates was accurate (Additional file 5). Further, the sensitivity of the yields to different biosynthetic demands, maintenance and changes in fluxes were probed. Shadow prices in the solution of the linear optimization problem of Flux Balance Analysis (FBA) define the sensitivity of the objective function with respect to each constraint indicating the utility of the metabolite in accelerating growth. Growth, as an objective in FBA, is defined as multiple simultaneous demands on precursors to make macromolecules related to biomass. In this context, a scaled shadow price for metabolites and scaled reduced costs for reactions that account for substrate and the growth yield are better sensitivity indicators (See Additional file 5 and Methods). The logarithmic sensitivity for cofactor NADH is the lowest in the ChlR populations followed by StrpR and differ from wild type at molar yields (Table 4). The 25% decrease in scaled shadow price values indicate that a compensation for that particular cofactor must have taken place during evolution resulting in a higher yield of that cofactor in the evolved strains for the already achieved higher growth and biomass yield. The logarithmic coefficients show that NADH and NADPH as compared to ATP may potentially play a role in increase in biomass yields through changing biomass composition. The reduced costs of each reaction indicate their significance in increasing the objective (growth). The alpha-keto glutarate dehydrogenase (AKGDH) reaction in the ChlR strain while the isocitrate lyase (ICL) reaction in the StrpR strain have scaled reduced costs associated with them (Table 4). Flux variability analysis (FVA) assesses the entire range of cellular function and the redundancy of optimal phenotypes. Applying FVA to identify reaction rates that can be uniquely determined allow us to explore the immutable or rigid metabolic state of the cell at maximal specific growth rate consistent with experimental data. Some reactions can be assigned fixed values, while the remaining calculable fluxes remain within the extreme bounds (Additional file 1: Table S4a). Uniquely computed reaction rates that are forced or fixed fluxes (coinciding upper and lower bounds) define metabolic rigidity and govern the plasticity of growth phenotype. Differences in these unique forced fixed rates in resistant and susceptible populations overall flux distribution indicate compensatory changes in metabolism due to antibiotic selection pressure. The constraints-based model identified two major features based on alternate optima predictions. Firstly, the resistant populations showed rigid flux distribution in secretion of overflow metabolites acetate and formate (Table 5). The onset of overflow metabolism and the details of secretion patterns were probed further using dynamic flux balance analysis (Additional file 1: Figure S3a - c). Both the resistant populations identified acetate as a common overflow metabolite. Dynamic FBA (dFBA) (Additional file 1: Figure S3a - c) qualitatively identified secretion of acetate and formate in that order in the ChlR population as indicated in the FVA and ethanol on lowering the oxygen uptake rates. The second major feature included reactions changing the rigid flux distribution in the wild type to a more flexible flux in the resistant populations. The reactions included alpha-keto glutarate dehydrogenase (AKGDH) and malate dehydrogenase (MDH) (Table 5). These change from a rigid flux configuration to more flexible one that potentially could modulate the direction and magnitude of flux involving NADH. The significance of alpha-keto glutarate (AKG) in the growth of the resistant populations is also supported by the high logarithmic sensitivity with respect to growth (Additional file 5). In order to probe this further, we looked at the scaled shadow prices, γredox and logarithmic sensitivity during experimental conditions (Additional file 1: Table S4b and Additional file 5) and identified NADH and NADPH to limit growth in the ChlR strain while ATP was also growth limiting in the StrpR strain. Thus, disruption of redox homeostasis through NADH/NAD ratios (Fig. 6, Table 5) and biomass precursor anabolism through NADPH/NADP ratios were identified as central to antibiotic action by FVA  and sensitivity/shadow price analysis . A recent report uses a NADH oxidase enzyme system to delineate the role of NADH imbalance and show decoupling of electron transfer via ETC. and proton pumping for ATP synthesis . Analysis of pareto fronts and trade-off between ATP and NADH and NADPH maintenance reactions (NADH/NADPH oxidases, not generally present in Chromobacterium) was performed to understand the modulation of NADH/NAD and NADPH/NADP ratios in growth. A reaction representing the NADH/NADPH oxidase (water forming) with the right balance of protons and oxygen was added to the model. This reaction essentially acts as a drain if there is excess NADH/NADPH in the system. On probing the relation to ATPase (representing ATP maintenance, ATPM) by a pareto front analysis (Additional file 1: Figure S3d - e, the differential relation for ChlR and StrpR was established. To reduce the growth of the ChlR strain to the wild type molar yield we identified NADH oxidase as a critical constraint (Pareto front analysis). For StrpR, however, both NADH and NADPH oxidases were needed. This is also suggested by γredox analysis (Additional file 1: Table S4b), which identified both NADH and NADPH as growth limiting in StrpR and only NADH as growth limiting in ChlR. At molar growth yields this suggests that excess NADH yields are indeed responsible for the excess growth associated with the resistance to chloramphenicol, while both NADH and NADPH yields play a role in StrpR. The rigidity of the flux held through AKGDH and MDH reactions was restored, when these constraints were added. Experiments confirmed an increase in NADH levels (Fig. 6e) in the ChlR population. For the StrpR population the NAD levels go up (Fig. 6e), seen in the molar yield simulations that show a 2 fold increase in flux through the NADH16 reaction that is the quinone associated conversion of NADH to NAD (Additional file 5).
The advent of genome-scale experimentation allows acquisition of heterogeneous data-types that are critical to delineating the genotype-phenotype relationship [1,2,3,4]. However, the mechanistic basis for the killing of antibiotic resistant populations only partially emerges through such data and requires integration of multiple data-types into a predictive scalable model. In this work we show that an integrative approach is able to predict the underlying mechanism related to redox and energy homeostasis operational at the level of cellular metabolism in C. violaceum. The model explains/predicts how benign metabolites in combination with antibiotics could potentially kill antibiotic resistant Chromobacterium populations by driving metabolism in a direction causing imbalance and disruption of the delicate redox or energy balance needed for the organism to survive.
Changed kinetic parameters (Table 1) of the two resistant populations in comparison to the susceptible wild type on several substrates indicated differential utilization and metabolic patterns resultant from altered genotypes and physiology via adaptive evolution of C. violaceum. Interestingly, fitness costs associated with the acquired resistance only manifested during growth on carbon sources and not in the environment used for evolution. The identification of four metabolites pyruvate, maleate, succinate and oxoadipate, potentially all electron donors, do not support growth across both the ChlR and StrpR resistant populations. The null post treatment viability count (Fig. 3e and f) make them ideal candidates for antibiotic therapy for resistant populations.
Also, these metabolites with the exception of maleate are electron donors and enter catabolism in central metabolism as glycolytic or TCA intermediates (Additional file 1: Figure S5). Each of the substrates is taken up and metabolized via specific dehydrogenases that involve the cofactor couple NAD/NADH. The results obtained highlight the importance of measuring fitness costs under multiple micro-environmental conditions. They provide a more relevant estimate of fitness in Chromobacterium and also reveal novel physiological weaknesses exploitable for drug development such as the redox homeostasis. Knowledge of such associated fitness costs in other pathogens can identify microenvironment metabolites that in combination with the antibiotic can target reduction of pathogen fitness.
The evolutionary fitness landscape of Chromobacterium can be viewed as a random adaptive (weighted) walk dictated by the drugs in the regime of strong selection and weak mutation (SSWM) [72, 73]. Such weighted walks in space of genotypes in the presence of chloramphenicol may have resulted in mutations in acrR, the transcription repressor of the tripartite AcrAB-TolC multidrug efflux pump and marC another multidrug efflux pump and may be potentially commutative. The continuous activity of the tripartite AcrAB-TolC multidrug efflux pump, which is proton dependent, could result in membrane potential changes due to efflux of small molecules like violacein. Antibiotics are known to activate the AcRAB-TolC pumps  and hence continuous de-repression may be a mechanism for evolution. The relatively few changes in genome sequence observed with chloramphenicol as selection pressure indicate a role beyond genetic causality in the resistant phenotypes.
Previously implicated mutations in rpsL (R86S; Fig. 4b) related to streptomycin resistance [63, 74] dictate an “error-restrictive”, hyper-accurate translation phenotype, accurate ribosomal function  and could explain the long lags and lowered growth rates observed for this population on specific carbon sources. Enhanced growth of the StrpR in poor carbon sources like glycerol probably occurs due to lower levels of the transcription factor, σS . Lack of growth on pyruvate, succinate and maleate indicate potential induction of rpoS . The mutated PabC protein involved in de novo folate biosynthesis via PABA could potentially lead to excess PLP known to perturb amino acid metabolism , including the observed tryptophan synthesis as observed in LC-HRMS with MS/MS profiles for StrpR intracellular temporal extracts. Streptomycin is known to cause leakage of low molecular weight antibiotic, ions and amino acids by damaging the permeability barrier creating potential secondary selection pressure . This implies low cytoplasmic K+ could shape the evolution as a secondary selection pressure specifically based on the mutated cytoplasmic domain of KdpD TCS regulator. On genotypic profiling and protein functional analysis we reach a hypothesis that the antibiotic adapted landscape of C. violaceum with a few beneficial mutations and fitness distributions potentially supports the Gillespie [72, 73] landscape model for evolution of resistance. Whether these genome changes are causal or correlational need to be investigated.
Violacein with its differential abundances (Fig. 5a) as compared to wild type in the StrpR (>50% increase) and ChlR populations (~50% reduction) could be a potential biomarker for resistance. This could be potentially due to indole known to be shared in antibiotic resistant E. coli to provide cover for more susceptible bacteria . C. violaceum accumulates tryptophan (downstream of indole) that is converted subsequently to virulence factor violacein. Potentially, developing a synthetic operon to convert tryptophan (found in all pathogens) to violacein may be useful as a visual reporter/biosensor of emergent resistance in even other infectious diseases. Reports exist that chloramphenicol lower intracellular indole and accumulation of tryptophan as observed . Further exploration is needed to connect violacein to strategic intracellular communication and resistant growth phenotypes of the pathogen.
Using the in silico central metabolic network reconstruction representing C. violaceum, iDB149, we were able to understand emergent properties of redox and energy homeostasis across susceptible and resistant populations. On shadow price and reduced costs analysis of growth limiting metabolites and reactions, AKGDH and ICL were identified in ChlR and and StrpR respectively. The alpha-keto glutarate dehydrogenase (AKGDH) reaction is known to regulate oxidative phosphorylation, lysine and tryptophan synthesis. Isocitrate lyase (ICL) reaction on the other hand, has been implicated in pathogenesis and persistence in Salmonella and resistance in Mycobacterium . The use of ICL involves shunting isocitrate through the glyoxylate shunt and bypassing part of the TCA cycle. This potential mechanistic difference between ChlR and StrpR in glucose metabolism is also evident in the intracellular metabolite profile (Additional file 1: Figure S2). AKGDH is also a known modulator for oxidative phosphorylation and dependent on the ADP/ATP ratios . The flux through the ATP synthase is higher in the StrpR strain (Additional file 5) suggesting a potential increase in respiration rates as previously described . The apparent lowering of ATP synthase flux potentially suggests a mechanism for decoupling of electron transfer from proton pumping in oxidative phosphorylation in the resistant populations. The StrpR population with mutations in rpsL may need a model of regulation and metabolic to explain the big change in growth rates and yields in presence of streptomycin. Flux variability analysis showed reactions that have flux category change in resistant strains compared to WT as a function of antibiotic involved redox cofactor balancing (NADH/NAD and NADPH/NADP ratios); increased TCA cycle/Oxidation Phosphorylation in the StrpR population and increased overflow metabolism population diverting from TCA cycle in the ChlR population. The wild type growth on pyruvate, malate and succinate also indicates the increase in NADH levels as a mechanism for survival in the presence of antibiotics (data not shown). The metabolism of these substrates after their uptake, all involve the use of the cofactor couple NAD/NADH . The pareto front analysis of the added NADH/NADPH oxidase to the ChlR and StrpR strains potentially show an increase in NADH as a mechanism of evolution to become resistant. The addition of the above mentioned compounds potentially increase the NADH levels further causing major redox imbalance. This cofactor imbalance prevents the maintenance of a rigid core flux distribution, through certain control nodes, eventually preventing minimal cellular function of the metabolic network for energy and biomass formation. The model thus predicts the emergence of NAD/NADH ratios and electron imbalance to be critical to survival and susceptibility of the antibiotic resistant phenotype and can be leveraged to make the resistant pathogens susceptible to antibiotics again.
Taken together, our data unveiled that disruption of redox homeostasis by certain benign metabolites as key to killing antibiotic resistant pathogens. To distinguish between causal and correlational factors in the evolution of antibiotic resistance more rigorous experiments and comprehensive genome-scale models may be needed. Reversing the phenotype by perturbing metabolism through modulating micro-environments could be one way of subverting the onset of the post antibiotic era. Determining the genetic and metabolic basis of ribosome targeting antibiotic’s resistance helps address the ‘fourth dimension’ of how heterogeneous networks in cells evolve simultaneously in space and time. This could lead to scalable pipelines integrating growth/metabolite/MIC profiling and constraints-based flux balance models for clinical isolates, ultimately leading to personalized treatment and individualized therapy.
Acetaldehyde dehydrogenase (acetylating)
Automated gain control
Alcohol dehydrogenase (ethanol)
Adaptive laboratory evolutions
Accurate mass - extracted ion chromatogram
Atmospheric Pressure - Matrix-assisted Laser Desorption/ionization
American type culture collection
ATP synthase (four protons for one ATP)
Colony-forming units per milliliter
Chloramphenicol resistant population of Chromobacterium violaceum
Constraints based reconstruction and analysis
Deoxyribose nucleic acid
Flux Balance Analysis
Formate transport via diffusion
Flux variability analysis
Full width at half maximum
Grams dry weight
Heated electrospray ion source
High performance liquid chromatography
High-resolution mass spectrometry
Luria Bertani Agar
Liquid chromatography–mass spectrometry
Mass by charge ratio
Maleate or Maleic acid
Mass extraction window
Minimum inhibitory concentration
Tandem mass spectrometry
Nicotinamide adenine dinucleotide (oxidised form)
Nicotinamide adenine dinucleotide (reduced form)
Nicotinamide adenine dinucleotide phosphate (oxidised form)
Nicotinamide adenine dinucleotide phosphate (reduced form)
Next generation sequencing
- OD600 :
Optical Density at 600 nanometre
Polymerase chain reaction
Pulse dynamic focussing
Pyruvate formate lyase
Phosphate reversible transport via symport
Pyruvate or Pyruvic acid
Streptomycin resistant population of Chromobacterium violaceum
Succinate or Succinic acid
Succinyl-CoA synthetase (ADP-forming)
Tricarboxylic acid Cycle
Wild Type Chromobacterium violaceum
Zwitterionic hydrophillic interaction Liquid Chromatography.
Fani F, Leprohon P, Légaré D, Ouellette M. Whole genome sequencing of penicillin-resistant Streptococcus pneumoniae reveals mutations in penicillin-binding proteins and in a putative iron permease. Genome Biol. 2011;12:R115.
Lobritz MA, Belenky P, Porter CBM, Gutierrez A, Yang JH, Schwarz EG, Dwyer DJ, Khalil AS, Collins JJ. Antibiotic efficacy is linked to bacterial cellular respiration. Proc Natl Acad Sci. 2015;112:8173–80.
Stoesser N, Batty EM, Eyre DW, Morgan M, Wyllie DH, Del Ojo Elias C, Johnson JR, Walker AS, Peto TEA, Crook DW. Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data. J Antimicrob Chemother. 2013;68:2234–44.
Toprak E, Veres A, Michel J-B, Chait R, Hartl DL, Kishony R. Evolutionary paths to antibiotic resistance under dynamically sustained drug selection. Nat Genet. 2012;44:101–5.
Renzoni A, Andrey DO, Jousselin A, Barras C, Monod A, Vaudaux P, Lew D, Kelley WL. Whole genome sequencing and complete genetic analysis reveals novel pathways to glycopeptide resistance in Staphylococcus aureus. PLoS One. 2011;6:e21577.
Dörries K, Schlueter R, Lalk M. Impact of antibiotics with various target sites on the metabolome of Staphylococcus aureus. Antimicrob Agents Chemother. 2014;58:7151–63.
Martínez JL, Rojo F. Metabolic regulation of antibiotic resistance. FEMS Microbiol Rev. 2011;35:768–89.
Suzuki S, Horinouchi T, Furusawa C. Prediction of antibiotic resistance by gene expression profiles. Nat Commun. 2014;5:1–12.
Derewacz DK, Goodwin CR, Mcnees CR, McLean JA, Bachmann BO. Antimicrobial drug resistance affects broad changes in metabolomic phenotype in addition to secondary metabolism. Proc Natl Acad Sci U S A. 2013;110:2336–41.
Livermore DM. beta-Lactamases in laboratory and clinical resistance. Clin Microbiol Rev. 1995;8:557–84.
Dwyer DJ, Belenky PA, Yang JH, MacDonald IC, Martell JD, Takahashi N, Chan CTY, Lobritz MA, Braff D, Schwarz EG, Ye JD, Pati M, Vercruysse M, Ralifo PS, Allison KR, Khalil AS, Ting AY, Walker GC, Collins JJ. Antibiotics induce redox-related physiological alterations as part of their lethality. Proc Natl Acad Sci U S A. 2014;111:E2100–9.
Allison KR, Brynildsen MP, Collins JJ. Metabolite-enabled eradication of bacterial persisters by aminoglycosides. Nature. 2011;473:216–20.
Peng B, Su Y, Li H, Han Y, Guo C, Tian Y, Peng X. Exogenous Alanine and/or Glucose plus Kanamycin Kills Antibiotic-Resistant Bacteria Article Exogenous Alanine and/or Glucose plus Kanamycin Kills Antibiotic-Resistant Bacteria. Cell Metab. 2015;21:249–61.
Charusanti P, Fong NL, Nagarajan H, Pereira AR, Li HJ, Abate EA, Su Y, Gerwick WH, Palsson BO. Exploiting adaptive laboratory evolution of Streptomyces clavuligerus for antibiotic discovery and overproduction. PLoS One. 2012;7:e33727.
McCloskey D, Palsson BØ, Feist AM. Basic and applied uses of genome-scale metabolic network reconstructions of Escherichia coli. Mol Syst Biol. 2013;9:661.
Raghunathan A, Reed J, Shin S, Palsson B, Daefler S. Constraint-based analysis of metabolic capacity of Salmonella typhimurium during host-pathogen interaction. BMC Syst Biol. 2009;3:38.
Reznik E, Mehta P, Segrè D. Flux imbalance analysis and the sensitivity of cellular growth to changes in metabolite pools. PLoS Comput Biol. 2013;9:e1003195.
Ibarra RU, Edwards JS, Palsson BO. Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth. Nature. 2002;420(November):20–3.
Durán N, Menck CF. Chromobacterium violaceum: a review of pharmacological and industrial perspectives. Crit Rev Microbiol. 2001;27:201–22.
Füller JJ, Röpke R, Krausze J, Rennhack KE, Daniel NP, Blankenfeldt W, Schulz S, Jahn D, Moser J. Biosynthesis of violacein, structure and function of l-Tryptophan oxidase VioA from Chromobacterium violaceum. J Biol Chem. 2016;291:20068–84.
Hoshino T. Violacein and related tryptophan metabolites produced by Chromobacterium violaceum: Biosynthetic mechanism and pathway for construction of violacein core. Appl Microbiol Biotechnol. 2011;91:1463–75.
Demoss RD, Evans NR. Physiological aspects of violacein biosynthesis in nonproliferating cells. J Bacteriol. 1959;78:583–8.
Farrar WE, O’dell NM. beta-Lactamase activity in Chromobacterium violaceum. J Infect Dis. 1976;134:290–3.
Fantinatti-Garboggini F, de Almeida R, Portillo VDA, Barbosa TAP, Trevilato PB, Neto CER, Coêlho RD, Silva DW, Bartoleti LA, Hanna ES, Brocchi M, Manfio GP. Drug resistance in Chromobacterium violaceum. Genet Mol Res. 2004;3:134–47.
Wiegand I, Hilpert K, Hancock REW. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat Protoc. 2008;3:163–75.
Zwietering MH, Jongenburger I, Rombouts FM, van’t Riet K. Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990;56:1875–81.
Kim DE, Chivian D, Baker D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004;32(Web Server):W526–31.
Wass MN, Kelley LA, Sternberg MJE. 3DLigandSite: predicting ligand-binding sites using similar structures. Nucleic Acids Res. 2010;38(Web Server issue):469–73.
Meyer H, Liebeke M, Lalk M. A protocol for the investigation of the intracellular Staphylococcus aureus metabolome. Anal Biochem. 2010;401:250–9.
Letisse F, Lindley ND. An intracellular metabolite quantification technique applicable to polysaccharide-producing bacteria. Biotechnol Lett. 2000;22:1673–7.
Blosser S, Gray KM. Extraction of violacein from Chromobacterium violaceum provides a new quantitative bioassay for N -acyl homoserine lactone autoinducers. J Microbiol Methods. 2000;40:47–55.
Smith CA, O’Maille G, Want EJ, Qin C, Trauger SA, Brandon TR, Custodio DE, Abagyan RSG. METLIN: a metabolite mass spectral database. Ther Drug Monit. 2005;27:747–51.
Strohalm M, Hassman M, Košata B, Kodíček M. mMass data miner: an open source alternative for mass spectrometric data analysis. Rapid Commun Mass Spectrom. 2008;22:905–8.
Bhattacharya N, Singh A, Ghanate A, Phadke G, Parmar D, Dhaware D, Basak T, Sengupta S, Panchagnula V. Matrix-assisted laser desorption/ionization mass spectrometry analysis of dimethyl arginine isomers from urine. Anal Methods. 2014;6:4602–9.
Ruckerbauer DE, Jungreuthmayer C, Zanghellini J. Predicting genetic engineering targets with Elementary Flux Mode Analysis : a review of four current methods. N Biotechnol. 2015;32:534–46.
Agren R, Liu L, Shoaie S, Vongsangnak W, Nookaew I, Nielsen J. The RAVEN toolbox and its use for generating a genome-scale metabolic model for Penicillium chrysogenum. PLoS Comput Biol. 2013;9:e1002980.
Ruppin E, Papin JA, de Figueiredo LF, Schuster S. Metabolic reconstruction, constraint-based analysis and game theory to probe genome- scale metabolic networks theory to probe genome-scale metabolic networks. Curr Opin Biotechnol. 2010;21:1–9.
Varma A, Palsson BO. Stoichiometric flux balance models quantitatively predict growth and metabolic by product secretion in wild type Escherichia coli W3110. Appl Environ Microbiol. 1994;60:3724–31.
Balibar CJ, Walsh CT. In Vitro Biosynthesis of Violacein from L -Tryptophan by the Enzymes VioA - E from Chromobacterium V iolaceum. Biochemistry. 2006;45:15444–57.
Creczynski-pasa TB, Antônio RV. Energetic metabolism of Chromobacterium violaceum. Genet Mol Res. 2004;3:162–6.
Demoss RD, Happel ME. Nutritional requirements of Chromobacterium violaceum. J Bacteriol. 1959;77:137–41.
Ryan KS, Balibar CJ, Turo KE, Walsh CT, Drennan CL. The violacein biosynthetic enzyme VioE shares a fold with lipoprotein transporter proteins. J Biol Chem. 2008;283:6467–75.
Lima-Bittencourt CI, Costa PS, Hollatz C, Raposeiras R, Santos FR, Chartone-Souza E, Nascimento AMA. Comparative biogeography of Chromobacterium from the neotropics. Antonie Van Leeuwenhoek. 2011;99:355–70.
Young C-C, Arun AB, Lai W-A, Chen W-M, Chou J-H, Shen F-T, Rekha PD, Kämpfer P. Chromobacterium aquaticum sp. nov., isolated from spring water samples. Int J Syst Evol Microbiol. 2008;58:877–80.
Martin PAW, Gundersen-Rindal D, Blackburn M, Buyer J. Chromobacterium subtsugae sp. nov., a betaproteobacterium toxic to Colorado potato beetle and other insect pests. Int J Syst Evol Microbiol. 2007;57:993–9.
Schellenberger J, Que R, Fleming RMT, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, Kang J, Hyduke DR, Palsson BØ. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.
Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.
Price ND, Reed JL, Palsson BO. Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol. 2004;2:886–97.
Fell DA, Small JR. Fat synthesis in adipose tissue. An examination of stoichiometric constraints. Biochem J. 1986;238:781–6.
Watson MR. A discrete model of bacterial metabolism. CABIOS. 1986;2:23–7.
Acevedo A, Aroca G, Conejeros R. Genome-scale NAD(H/(+)) availability patterns as a differentiating feature between Saccharomyces cerevisiae and Scheffersomyces stipitis in relation to fermentative metabolism. PLoS One. 2014;9:e87494.
Feist AM, Herrgard MJ, Thiele I, Reed JL, Palsson BO. Reconstruction of biochemical networks in microorganisms. Nat Rev Microbiol. 2009;7:129–43.
Schuster S, Pfeiffer T, Fell DA. Is maximization of molar yield in metabolic networks favoured by evolution? J Theor Biol. 2008;252:497–504.
Reed JL, Famili I, Thiele I, Palsson BO. Towards multidimensional genome annotation. Nat Rev Genet. 2006;7:130–41.
Maarleveld TR, Khandelwal RA, Olivier BG, Teusink B, Bruggeman FJ. Basic concepts and principles of stoichiometric modeling of metabolic networks. Biotechnol J. 2013;8:997–1008.
Varma A, Palsson BO. Metabolic capabilities of Escherichia coli II. Optimal growth patterns. J Theor Biol. 1993;165:503–22.
Mahadevan RÃ, Schilling CH. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003;5:264–76.
Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc. 2007;2:727–38.
Haselkorn R, Artur L, Bataus M, Batista S, Teno C. The complete genome sequence of Chromobacterium violaceum reveals remarkable and exploitable bacterial adaptability. Proc Natl Acad Sci U S A. 2003;100:11660–5.
Elkins CA, Mullis LB, Lacher DW, Jung CM. Single nucleotide polymorphism analysis of the major tripartite multidrug efflux pump of Escherichia coli: Functional conservation in disparate animal reservoirs despite exposure to antimicrobial chemotherapy. Antimicrob Agents Chemother. 2010;54:1007–15.
Li M, Gu R, Su C, Routh MD, Harris KC, Jewell S, Mcdermott G, Yu EW. Crystal structure of the transcriptional reguator AcrR from Escherichia coli. J Mol Biol. 2008;374:591–603.
Okusu H, Nikaido H. AcrAB efflux pump plays a major role in the antibiotic resistance phenotype of Escherichia coli Multiple-Antibiotic-Resistance (Mar) mutants. J Antibiot (Tokyo). 1996;178:306–8.
Carter AP, Clemons WM, Brodersen DE, Morgan-warren RJ, Wimberly BT, Ramakrishnan V. Functional insights from the structure of the 30S ribosomal subunit and its interactions with antibiotics. Nature. 2000;407:340–8.
Panecka J, Mura C, Trylska J. Interplay of the bacterial ribosomal a-site, S12 protein mutations and paromomycin binding: a molecular dynamics study. PLoS One. 2014;9:e111811.
Agarwal D, Gregory ST, O’Connor M. Error-prone and error-restrictive mutations affecting ribosomal protein S12. J Mol Biol. 2011;410:1–9.
Demirci H, Wang L, Iv FVM, Murphy EL, Carr JF, Blanchard SC, Jogl G, Dahlberg AE, Gregory ST. The central role of protein S12 in organizing the structure of the decoding site of the ribosome. RNA. 2013;19:1791–801.
Freeman ZN, Dorus S, Waterfield NR. The KdpD/KdpE two-component system: integrating K+ homeostasis and virulence. PLoS Pathog. 2013;9:e1003201.
Green JM, Merkel WK, Nichols BP. Characterization and sequence of Escherichia coli pabC, the gene encoding aminodeoxychorismate lyase, a pyridoxal phosphate-containing enzyme. J Bacteriol. 1992;174:5317–23.
Ye QZ, Liu J, Walsh CT. p-Aminobenzoate synthesis in Escherichia coli: purification and characterization of PabB as aminodeoxychorismate synthase and enzyme X as aminodeoxychorismate lyase. Proc Natl Acad Sci U S A. 1990;87:9391–5.
Fuhrer T, Sauer U. Different biochemical mechanisms ensure network-wide balancing of reducing equivalents in microbial metabolism. J Bacteriol. 2009;191:2112–21.
Titov DV, Cracan V, Goodman RP, Peng J, Grabarek Z, Mootha VK. Complementation of mitochondrial electron transport chain by manipulation of the NAD+/NADH ratio. Science (80-). 2016;352:231–5.
Gillespie JH. A simple stochastic gene substitution model. Theor Popul Biol. 1983;23:202–15.
Gillespie J. Molecular evolution over the mutational landscape. Evolution (N Y). 1984;38:1116–29.
Kenney TJ, Churchward G. Cloning and sequence analysis of the rpsL and rpsG genes of Mycobacterium smegmatis and characterization of mutations causing resistance to streptomycin. J Bacteriol. 1994;176:6153–6.
Paulander W, Maisnier-Patin S, Andersson DI. The fitness cost of streptomycin resistance depends on rpsL mutation, carbon source and RpoS (sigmaS). Genetics. 2009;183:539–46.
Wolfe AJ. The acetate switch. Microbiol Mol Biol Rev. 2005;69:12–50.
Di Salvo ML, Budisa N, Contestabile R. PLP-dependent enzymes : a powerful tool for metabolic synthesis of non-canonical amino acids. In Beilstein Bozen Symposium on Molecular Engineering and Control; 2012:27–66. http://www.beilstein-institut.de/download/365/bozen2012.epub
Iscla I, Wray R, Wei S, Posner B, Blount P. Streptomycin potency is dependent on MscL channel expression. Nat Commun. 2014;5:4891.
Lee HH, Molla MN, Cantor CR, Collins JJ. Bacterial charity work leads to population-wide resistance Henry. Nature. 2010;467:82–5.
Gibson F, Mcdougall B, Jones MJ, Teltscher H. The action of antibiotics on indole synthesis by cell suspensions of Escherichia coli. J Gen Microbiol. 1956;446–458.
Dunn MF, Ramírez-Trujillo JA, Hernández-Lucas I. Major roles of isocitrate lyase and malate synthase in bacterial and fungal pathogenesis. Microbiology. 2009;155:3166–75.
McLain AL, Szweda PA, Szweda LI. α-Ketoglutarate dehydrogenase: a mitochondrial redox sensor. Free Radic Res. 2011;45:29–36.
We would like to acknowledge Priyanka P. Sadhukhan and Sudeepa Nandi for their initial involvement in conceptualization of the metabolic core model, FVA categorization algorithm and ALE experiments.
This work was supported in part by grants from Department of Science and Technology (GAP292826, SERB grant EMR/2016/000167), CSIR India (network grant HCP003, TFYP CMET grant CSC0110 and infrastructure grant) DB is supported through a DST INSPIRE (Ph. D.) Fellowship. VP further acknowledges Department of Biotechnology for Ramalingawami Fellowship and RGYI award.
Availability of data and materials
The genome sequence data generated for the current study is available in Sequence Read Archive (SRA) repository under accession number SRP072862 (http://www.ncbi.nlm.nih.gov/sra/?term=SRP072862). All other data generated or analysed during this study are included in this published article and its supplementary information files.
The authors declare that they have no competing interests.
AR and DB conceptualized & designed the study, analyzed & interpreted the data and wrote the final manuscript. DB performed all experiments related to growth, evolution, genotyping and phenotyping; performed all constraints-based simulations & worked on data visualization. VP, DP and NB were involved in development of methods for mass spectrometric analysis and performed the analysis. AG was involved majorly in MS data analysis that included statistical tools among others. All authors read the manuscript and approve the content.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Tables S1 to S6 and Figures S1 to S5. (PDF 2833 kb)
Growth profiles for WT, ChlR and StrpR using 30 different C/N substrates when antibiotic was added from zero hour timepoint showing WT with no antibiotic in media (blue, wt_noab), ChlR with chloramphenicol added to the media (red, DD001_t0chl), StrpR with streptomycin added to the media (green, DD006_t0strep). Plots made using GraphPad Prism v6.01 and n = 3. (TIF 1004 kb)
Growth profiles for WT, ChlR and StrpR using 30 different C/N substrates for first six hours when no antibiotic was added (antibiotic was added after these 6 h, and hence represented as ‘t6’) showing WT (blue) with no antibiotic in first 6 h, ChlR (DD001_t6chl, red), StrpR (DD006_t0strep, green) with no antibiotic in media. Plots made using GraphPad Prism v6.01 and n = 3. (TIF 778 kb)
Growth profiles for WT, ChlR and StrpR using 30 different C/N substrates post addition of antibiotic at 6 h (t6) showing WT with chloramphenicol and streptomycin, (WT_t6chl,blue, and WT_t6strep, red, respectively), ChlR with chloramphenicol added to the media (DD001_t6chl,green), StrpR with with streptomycin added to the media (DD006_t0strep, violet). Plots made using GraphPad Prism v6.01 and n = 3. (TIF 1060 kb)
This file contains information about the C. violaceum core metabolic model. The file includes metabolites, reactions, and other simulations used to understand the in silico behaviour of redox state of the resistant populations ChlR and StrpR in comparison to WT in the subsequent tabs. (XLS 384 kb)
This file contains information about 59 metabolites analyzed using HPLC-HESI-HRMS. The file includes information including the retention time information, QC data peak area, normalized peak area among others. (XLS 292 kb)
About this article
Cite this article
Banerjee, D., Parmar, D., Bhattacharya, N. et al. A scalable metabolite supplementation strategy against antibiotic resistant pathogen Chromobacterium violaceum induced by NAD+/NADH+ imbalance. BMC Syst Biol 11, 51 (2017) doi:10.1186/s12918-017-0427-z
- Antibiotic resistance
- Flux balance analysis
- Flux variability analysis
- Redox homeostasis