# Multiparametric Profiling for Identification of Chemosensitizers against Gram-Negative Bacteria

Vincent Lôme, Jean-Michel Brunel, Jean-Marie Pagès, Jean-Michel Michel Bolla

► To cite this version:

Vincent Lôme, Jean-Michel Brunel, Jean-Marie Pagès, Jean-Michel Michel Bolla. Multiparametric Profiling for Identification of Chemosensitizers against Gram-Negative Bacteria. Frontiers in Microbiology, 2018, 9, 10.3389/fmicb.2018.00204. [https://doi.org/10.3389/fmicb.2018.00204](https://doi.org/10.3389/fmicb.2018.00204). [hal-01824504](https://HAL Id: hal-01824504)

Submitted on 27 Jun 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Multiparametric Profiling for Identification of Chemosensitizers against Gram-Negative Bacteria

Vincent Lôme¹, Jean-Michel Brunel², Jean-Marie Pagès¹ and Jean-Michel Bolla¹*

¹ UMR MD1, Aix-Marseille University, IRBA, TMCD2, Facultés de Médecine et de Pharmacie, Marseille, France, ² Centre de Recherche en Cancérologie de Marseille (CRCM), CNRS, UMR7258, Institut Paoli Calmettes, Aix-Marseille Université, UM 105, Inserm, U1068, Faculté de Pharmacie, Marseille, France

Antibiotic resistance is now a worldwide therapeutic problem. Since the beginning of anti-infectious treatment bacteria have rapidly shown an incredible ability to develop and transfer resistance mechanisms. In the last decades, the design variation of pioneer bioactive molecules has strongly improved their activity and the pharmaceutical companies partly won the race against the clock. Since the 1980s, the new classes of antibiotics that emerged were mainly directed to Gram-positive bacteria. Thus, we are now facing to multidrug-resistant Gram-negative bacteria, with no therapeutic options to deal with them. These bacteria are mainly resistant because of their double membrane that conjointly impairs antibiotic accumulation and extrudes these molecules when entered. The main challenge is to allow antibiotics to cross the impermeable envelope and reach their targets. One promising solution would be to associate, in a combination therapy, a usual antibiotic with a non-antibiotic chemosensitizer. Nevertheless, for effective drug discovery, there is a prominent lack of tools required to understand the rules of permeation and accumulation into Gram-negative bacteria. By the use of a multidrug-resistant enterobacteria, we introduce a high-content screening procedure for chemosensitizers discovery by quantitative assessment of drug accumulation, alteration of barriers, and deduction of their activity profile. We assembled and analyzed a control chemicals library to perform the proof of concept. The analysis was based on real-time monitoring of the efflux alteration and measure of the influx increase in the presence of studied compounds in an automatized bio-assay. Then, synergistic activity of compounds with an antibiotic was studied and kinetic data reduction was performed which led to the calculation of a score for each barrier to be altered.

**Keywords:** antibiotic resistance, combination therapy, automated platform, whole-cell screening, hit-to-lead

**OPEN ACCESS**

**Edited by:**
Henrietta Venter,
University of South Australia, Australia

**Reviewed by:**
Jon William Weeks,
United States Food and Drug
Administration, United States
Hiroshi Nikaido,
University of California, Berkeley,
United States

*Correspondence:
Jean-Michel Bolla
jean-michel.bolla@univ-amu.fr

**Specialty section:**
This article was submitted to *Antimicrobials, Resistance and Chemotherapy*, a section of the journal *Frontiers in Microbiology*

Received: 16 November 2017

Accepted: 29 January 2018

Published: 19 February 2018

**Citation:**
Lôme V, Brunel J-M, Pagès J-M and Bolla J-M (2018) Multiparametric Profiling for Identification of Chemosensitizers against Gram-Negative Bacteria. *Front. Microbiol.* 9:204. doi: 10.3389/fmicb.2018.00204

## INTRODUCTION

Historically, the activity improvement of antibiotics has been widely performed by design variation of pioneer bioactive molecules (Bush et al., 2011; Silver, 2011). Since the 1980s, new classes of antibiotics have emerged but mainly are active against Gram-positive bacteria (Lomovskaya et al., 2006). Recent target-based high-throughput screening programs along with *in silico* studies have led to identification of hits with high potentiality. Although this strategy appears attractive, the major drawback of target-based assays is that they fail to consider the membrane translocation
barriers, comprising of the bacterial permeation and the efflux pump issues (Payne et al., 2006; Stavenger and Winterhalter, 2014; Tommasi et al., 2015; Zgurskaya et al., 2015). In this context, one of the greatest challenges for the design of new scaffolds of interest against Gram-negative bacteria is to promote their intrabacterial accumulation (Nikaido, 1994; Lomovskaya et al., 2006; Spellberg and Shlaes, 2014; Stavenger and Winterhalter, 2014; Zgurskaya et al., 2015). Hence, combination therapy of a usual antibiotic with a non-antibiotic chemosensitizer seems to be one of the most promising solutions (Kristiansen et al., 2007; Mazumdar et al., 2009), being able to increase antibiotics accumulation through non-specific synergy mechanisms such as permeability, enhancement, and efflux impairment (Lewis, 2013). Nevertheless, in effective drug discovery, there is a prominent lack of tools required to understand the rules of permeation (Masi et al., 2017) and accumulation into Gram-negative bacteria (Harvey et al., 2015; Schneider et al., 2017).

We introduce a high-content screening method for chemosensitizers discovery by quantitative assessment of drug accumulation barriers alteration and deduction of their activity profile. We adapted a series of whole-cell-based assays (Figure 1B-D) to the multidrug-resistant *Enterobacter aerogenes* EA289 clinical isolate since this strain presents a decreased outer membrane permeability as it does not express major porins, and an enhanced efflux transport toward antibiotics (Mallea et al., 1998; Pradel and Pagès, 2002). We assembled and analyzed a control chemicals library to perform the proof of concept.

A chemosensitizer compound can promote drug accumulation by altering at least one major cellular barrier (Figure 1A). Hence, the first stage of analysis was based on two ways: real-time monitoring of the RND efflux alteration using the 1,2'-diNA dye and measurement of increased influx using the nitrocefin probe in the presence of studied compounds.

In this context, kinetic data reduction was performed taking into account the biological variability among replicates, the non-specific absorbance, or fluorescence interactions between the compound and the probe, which led to the calculation of a score for each barrier to be altered.

## MATERIALS AND METHODS

### Strains Used in This Study

Two bacterial strains were used in this study. The *E. aerogenes* EA289 strain is a Kan⁸ derivative of the MDR clinical isolate Ea27 (Mallea et al., 1998); the *acrB* mutant EA289Δ*acrB* strain was constructed from the EA289 strain (Pradel and Pagès, 2002). Strains were maintained at -80°C in 15% (v/v) glycerol for cryoprotection. Bacteria were routinely grown in Cation-adjusted Mueller Hinton Broth (CAMHB) at 37°C.

### Glucose-Triggered 1,2'-Efflux Assay

A single colony of *E. aerogenes*, EA289 or EA289Δ*acrB*, from an overnight plate, was grown in CAMHB shaking at 250 rpm at 37°C until the OD₆₀₀nm = 0.5 was reached. Bacteria were collected by centrifugation and re-suspended

at OD₆₀₀nm = 0.25 in potassium phosphate buffer, K₂HPO₄ 20 mM, MgCl₂ 1 mM, pH 7.0 (PPB), supplemented with the proton conductor carbonyl cyanide m-chlorophenyl-hydrazone (CCCP) 5 µM, that allows the inactivation of active efflux by de-energization of the cytoplasmic membrane (Nagano and Nikaido, 2009). Bacteria were then incubated overnight with 1,2'-dinaphthylamine (1,2'-diNA) 32 µM at 37°C. After overnight incubation, the cells were washed in PPB and aliquoted in microplates (100 µl/well) in wells pre-loaded with compounds 100 µM or controls. Membrane incorporated 1,2'-diNA was followed by monitoring the fluorescence every 30 s at 37°C ($λ_{ex}$ = 370 nm; $λ_{em}$ = 420 nm). Glucose 50 mM was added at 300 s to initiate bacterial energization. Plates were read on a Tecan Infinite M200 Pro multimode plate reader. Assays were performed in Greiner Bio-One 96-well plates, ref 675076 (half area, black with solid bottom). 1,2'-diNA ref D2988 was purchased from T.C.I (Zwijndrecht, Belgium) (Brunel et al., 2013).

The pipeting was performed by a Tecan Freedom EVO automated liquid handler equipped with eight fixed PTFE-coated tips, with a robotic arm to transfer plates to the Tecan Infinite M200 Pro multimode plate reader. Automation protocols were coded on Tecan EVOware and are available upon request (Tecan, Lyon, France).

### Glucose-Triggered 1,2'-diNA Efflux Assay, Kinetic Data Reduction (Table 1)

Kinetic profiles observed in raw data are outlined in Supplementary Figure 1. Typically, a-outlined curves were

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Chemosensitizer promoting drug accumulation (A) and real-time assays used in this study (B-D). Chemosensitizer possible modes of action (A), RND efflux inhibition screening (B), OM permeabilization screening (C), and transmembrane potential disruption assay (D).</figcaption></figure>
<table><caption>**TABLE 1** | Glucose-triggered 1,2'-diNA efflux assay, kinetic data reduction.</caption><thead><tr><th>Kinetic transformed data</th><th>Comments</th></tr></thead><tbody><tr><td>basisRFU<sub>A</sub>(Δ, NC)</td><td>Pre-energization fluorescence intensity of EA289ΔacrAB in the presence of DMSO 1%</td></tr><tr><td>basisRFU<sub>A</sub>(Δ, x)</td><td>Pre-energization fluorescence intensity of EA289ΔacrAB in the presence of a tested compound</td></tr><tr><td>areaRFU<sub>S<sub>B</sub></sub>(WT, NC)</td><td>Area under the curve after dye efflux in EA289 in the presence of DMSO 1%</td></tr><tr><td>finalRFU<sub>S<sub>B</sub></sub>(WT, NC)</td><td>As: areaRFU<sub>S<sub>B</sub></sub>(WT, NC) ÷ 300. Mean fluorescence intensity after dye efflux in EA289 in the presence of DMSO 1%</td></tr><tr><td>areaRFU<sub>S<sub>B</sub></sub>(WT, x)</td><td>Area under the curve after dye efflux in EA289 in the presence of a tested compound</td></tr><tr><td>finalRFU<sub>S<sub>B</sub></sub>(WT, x)</td><td>As: areaRFU<sub>S<sub>B</sub></sub>(WT, x) ÷ 300. Mean fluorescence intensity after dye efflux in EA289 in the presence of a tested compound</td></tr></tbody></table>

non-specific fluorescence controls obtained for each compound incubated with efflux deficient mutant EA289ΔacrAB (Δ), b were positive results for 1,2'-diNA efflux inhibition curves incubated with the wild-type (WT) strain EA289, c were control and negative results incubated with WT EA289. Quenched kinetic signals not to be analyzed were reported in **Supplementary Figure 2**, where a-outlined curves were fluorescence controls obtained incubated with EA289ΔacrAB, b were the same tested compounds as a-curves incubated with the WT EA289, and c were stronger quencher compounds incubated with EA289 or with EA289ΔacrAB.

Kinetic transformed data basisRFU, areaRFU.s were calculated systematically, from raw kinetic data and served as the basis for score calculation and quencher compounds detection. basisRFU_A was calculated as the mean value of the 10 first data points (A period) and represents the pre-energization fluorescence intensity [AU]; areaRFU.sB was calculated as the area under the curve over data point 19 through data point 29 (B period) and represents the area under the curve of remaining 1,2'-diNA after transport [AU.t] (**Supplementary Figure 3**).

As observed in **Supplementary Figure 1a**, and consistent with the initial hypothesis, 1,2'-diNA efflux kinetics in the presence of EA289ΔacrAB were not sensitive to the effect of potential efflux inhibitors. However, a compound leading to the translation of EA289ΔacrAB efflux kinetic was considered a fluorescence quencher. The fluorescence quenching calculation is based on the comparison of the basisRFU_A(Δ, NC) and basisRFU_A(Δ, x) values for EA289ΔacrAB, in the absence and presence of tested compound, respectively.

If the value basisRFU_A(Δ, x) < basisRFU_A(Δ, NC) × 0.7 the fluorescence signal was considered quenched too much and a proportionality factor correction could not effectively report kinetics. Such data were not further analyzed. It is important to note that a compound having such properties on fluorescence signal is not indicative of its potential inhibitor character of efflux systems. The characterization of this activity

must be pursued by other experimental methods. Otherwise, i.e., if the value basisRFU_A(Δ, x) > basisRFU_A(Δ, NC) × 0.7 we considered that the fluorescence signal may be corrected by proportionality factor and efflux inhibition score (%EIS) may be calculated.

The %EIS calculation comprised two steps:

(1) The finalRFU.sB(WT, x) values obtained for tested compounds in the presence of the WT strain were adjusted using the mean basisRFU_A(Δ, NC) value calculated for the eight negative control (NC) replicates of the mutant strain EA289ΔacrAB, and the basisRFU_A(Δ, x) values obtained for tested compounds in the presence of the mutant strain.

$$ \text{finalRFU.sB,adjusted}(\text{WT}, x) = \frac{\text{finalRFU.sB}(\text{WT}, x) \times \text{basisRFU_A}(\Delta, \text{NC})}{\text{basisRFU_A}(\Delta, x)}. $$

This step both participates in robustness of the results by correcting moderate shift in basis fluorescence observed in **Supplementary Figure 1a**, and detects strong quenchers.

(2) The mean areaRFU.sB(WT, NC) value calculated for the eight NC replicates of the WT strain was subtracted to the adjusted areaRFU.sB,adjusted (WT, x) values obtained for tested compounds in the presence of the WT strain. This step ensures results consistency between assays by correcting the biological variability among clones. The obtained values were normalized using the maximum efflux of the WT strain which was calculated by subtracting the mean basisRFU_A(WT, x) value calculated for the eight NC replicates multiplied by the duration of the A period (300 s) to the mean areaRFU.sB(WT, NC) value calculated for the eight NC replicates of the WT strain. Finally, the %EIS was obtained as:

$$ \%EIS(x) = \frac{\text{finalRFU.sB,adjusted}(\text{WT}, x) - \text{finalRFU.sB}(\text{WT}, \text{NC})}{\text{basisRFU_A}(\text{WT}, \text{NC}) - \text{finalRFU.sB}(\text{WT}, \text{NC})} \times 100. $$

The Z'-factor (Zhang et al., 1999) was calculated based on the average of four screens, using the EA289ΔacrAB strain incubated with 1% DMSO as a positive control (PC) and the WT strain incubated with 1% DMSO as NC, in the calculation.

Kinetic data transformation, score calculation, and false positive detection were programmed and performed automatically by Tecan Magellan software. Both data from Magellan and pipetting instructions from EVOware were sent to a custom LIMS (Modul-Bio, Marseille, France), which matched results to compound names.

## Outer Membrane Permeability Assay

A single colony of *E. aerogenes*, EA289, from an overnight plate, was grown in CAMHB shaking at 250 rpm at 37°C until the OD<sub>600nm</sub> = 0.5 was reached. Bacteria were collected by centrifugation and were re-suspended at OD<sub>600nm</sub> = 0.25 in PPB, supplemented with CCCP 5 µM. Bacteria were mixed with nitrocefin 50 µg ml<sup>-1</sup> before addition of compounds 100 µM. Nitrocefin hydrolysis was followed by monitoring the absorbance
<table><caption>**TABLE 2** | Outer membrane permeability assay, kinetic data reduction.</caption><thead><tr><td>Kinetic transformed data</td><td>Comments</td></tr></thead><tbody><tr><td>basisOD<sub>A</sub>(x)</td><td>OD value at origin in the presence of a tested compound</td></tr><tr><td>maxSlopeOD/hr<sub>B</sub>(x)</td><td>Initial slope in the presence of a tested compound</td></tr><tr><td>areaOD.s<sub>C</sub>(x)</td><td>Area under the curve in the presence of a tested compound</td></tr><tr><td>areaOD.s<sub>C</sub>(NC)</td><td>Area under the curve in the presence of DMSO 1%</td></tr><tr><td>maxOD<sub>D</sub>(PC)</td><td>Maximum OD for the PC</td></tr></tbody></table>

($\lambda_{abs}$ = 490 nm). Cell suspension was added at 100 µl per well and the absorbance read every 30 s at 37°C. Assays were performed in Greiner Bio-One 96-well plates, ref 675101 (half area, clear with flat bottom). Nitrocefin SR0112C was purchased from Oxoid (Basingstoke, United Kingdom) (Matsumoto et al., 2011).

The pipetting was performed by a Tecan Freedom EVO automated liquid handler equipped with eight fixed PTFE-coated tips, with a robotic arm to transfer plates to the Tecan Infinite M200 Pro multimode plate reader. Automation protocols were coded on Tecan EVOware and are available upon request.

## Outer Membrane Permeability Assay, Kinetic Data Reduction (Table 2)

Kinetic profiles observed within raw data are outlined in **Supplementary Figure 4**. Typically, a-outlined curves were positive results, b were NC and true negatives, and Y-values at origin read at position c corresponded to the absorbance of the bacterial suspension. Biased kinetic signals were reported in **Supplementary Figure 5**, where a-outlined curves were false positives, b were associated with compounds competitively hydrolyzed with the nitrocefin β-lactam (compounds belonging to β-lactam antibiotics or β-lactamase inhibitors, families). Y-values at origin read at position c were expected to be 0.25 the OD<sub>600nm</sub> for bacterial suspension, whereas d and e corresponded to higher OD<sub>600nm</sub> due to non-specific compound absorbance at the same wavelength as the peak absorbance of hydrolyzed nitrocefin. False positive detection can be anticipated by naked eye observation of the plate during pipetting, concerned wells appearing red.

Kinetic transformed data basisOD, maxSlopeOD/hr, areaOD.s, and maxOD were calculated systematically, from raw kinetic data and served as the basis for score calculation and false positive detection. basisOD<sub>A</sub> was calculated as the mean value of the five first data points (A period) and represents the Y-axis value at origin [OD], maxSlopeOD/hr<sub>B</sub> was calculated as the maximum slope using three consecutive data points within the 10 first data points (B period) and represents the mean slope at origin [OD/t], areaOD.s<sub>C</sub> was calculated as the area under the curve over the 80 first data points (C period) [OD.t], maxOD<sub>D</sub>(PC) was calculated as the maximum OD value using five consecutive data points over data point 71 through data point 81 (D period) in the presence of sodium dodecyl sulfate (SDS) 0.5% as the PC condition (Supplementary Figure 6).

Possible false positives showed kinetic curve starts that were parallel to the NC one (Supplementary Figure 5a). Data can be corrected by translation of the curve to the NC one, and data may be saved. The basisOD<sub>A</sub>(x) value can be significantly greater than what is observed with the NC both in the case of a true positive or in the case of false positives (Supplementary Figures 4a, and 5a, respectively). The distinction between true and false positives can be done through their initial slope maxSlopeOD/hr<sub>B</sub>(x): false positives showed an initial slope that was similar to the NC one, whereas true positives showed a significantly greater one. According to repeated observations the threshold value for maxSlopeOD/hr<sub>B</sub>(x) was arbitrarily set to 1 OD hr<sup>-1</sup>.

If maxSlopeOD/hr<sub>B</sub>(x) < 1 OD hr<sup>-1</sup>, the initial slope was considered to be close to the NC one. The translation may be then applied within outer membrane permeabilization score (%OPS) calculation. Otherwise, i.e., if maxSlopeOD/hr<sub>B</sub>(x) > 1 OD hr<sup>-1</sup>, the initial slope was considered significant, due to rapid hydrolytic reaction start, before reading has taken place. It is a specific absorbance of hydrolyzed nitrocefin, and translation must therefore not be applied.

The %OPS calculation comprised two steps:

(1) All curves were translated so their Y-axis value at origin equals 0, by subtracting the basisOD<sub>A</sub>(x) value multiplied by the duration of the C period (2400 s). This step both corrects false positives observed in Supplementary Figure 5a, and participates in robustness of the result by overcoming any differences in optical path length between well (optical meniscus may not be identical for all wells, even after the agitation of the plate). The differences in endpoint value may be infinitesimal but they are exacerbated by the area under the curve calculation and can alter the further ranking of low activity compounds. This calculation was performed systematically except when a rapid acting true positive was detected. Therefore, such curves did not benefit the subtraction of basisOD<sub>A</sub>(x) × 2400 adjustment nonetheless optical path variability did not interfere with rapid acting true positives ranking.

$$ \text{areaOD.s}_{\text{C,adjusted}}(x) = \text{areaOD.s}_{\text{C}}(x) - (\text{basisODA}(x) \times 2400). $$

(2) The mean areaOD.s<sub>C,adjusted</sub>(NC) value calculated for the eight NC replicates was subtracted to the mean areaOD.s<sub>C,adjusted</sub>(x) value calculated for the two replicates of each tested compound. This step ensures results consistency between assays by correcting the biological variability among clones.

Within the same calculation step, obtained values were normalized, using the PC condition corresponding to bacterial cells lysis in the presence of SDS 0.5%. The mean maxOD<sub>D</sub>(PC) value was calculated for the eight PC replicates, and the %OPS was calculated as:

$$ \%OPS(x) = \frac{\text{areaOD.s}_{\text{C,adjusted}}(x) - \text{areaOD.s}_{\text{C,adjusted}}(\text{NC})}{(\text{maxOD}_{\text{D}}(\text{PC}) \times 2400) - \text{areaOD.s}_{\text{C}}(\text{NC})} \times 100. $$
The Z'-factor (Zhang et al., 1999) was calculated based on the average of two screens, using chlorhexidine 100 µM as PC and 1% DMSO as NC, in the calculation.

Kinetic data transformation, score calculation, and false positive detection were programmed and performed automatically by Tecan Magellan software. Both data from Magellan and pipetting instructions from EVOware were sent to a custom LIMS (Modul-Bio, Marseille, France), which matched results to compound names.

## Transmembrane Potential Disruption

A single colony of *E. aerogenes*, EA289, from an overnight plate, was grown in CAMHB shaking at 250 rpm at 37°C until the OD<sub>600nm</sub> = 0.5 was reached. Bacteria were collected by centrifugation and were re-suspended at OD<sub>600nm</sub> = 0.25 in HEPES 5 mM, EDTA 10 mM, pH 7.0, and then washed in HEPES 5 mM, pH 7.0 with 3-3'-dipropylthiadicarbocyanine iodide (diSC<sub>3</sub>(5)) 8 µM. The membrane potential-sensitive cyanine dye diSC<sub>3</sub>(5) distributes between cells and the medium depending on the cytoplasmic membrane potential gradient. Released diSC<sub>3</sub>(5) was quantified by measuring the fluorescence (λ<sub>ex</sub> = 622 nm; λ<sub>em</sub> = 690 nm) 300 s after the addition of compounds 100 µM. A control experiment was performed for every tested compound where the cells were treated with chlorhexidine 300 µM as lysis solution, to normalize the results. Cell suspension was added at 100 µl/well and the fluorescence read every 30 s at 37°C. Plates were read on a Tecan Infinite M200 Pro multimode plate reader. Assays were performed in Greiner Bio-One 96-well plates, ref 675076 (half area, black with solid bottom). 3,3'-dipropylthiadicarbocyanine iodide (diSC<sub>3</sub>(5)) ref 84923 was purchased from Anaspec Inc. (Fremont, CA, United States).

The pipetting was performed by a Tecan Freedom EVO automated liquid handler equipped with eight fixed PTFE-coated tips, with a robotic arm to transfer plates to the Tecan Infinite M200 Pro multimode plate reader. Automation protocols were coded on Tecan EVOware and are available upon request.

Results are summarized in **Supplementary Table 2**. Typically, the average fluorescence intensity after the addition of compounds was between 1000 and 3000 A.U. for NC (incubated with DMSO 1%) and negative results, and between 4000 and 10,000 A.U. for positive results. The PC condition included compounds with the lysis solution and showed fluorescence intensities between 45,000 and 60,000 A.U.

## Monodose Chemosensitization Assay (Table 3)

Previously selected hits from the control chemicals library were tested at 10 µM for 48 h growth inhibition of *E. aerogenes* EA289. Compounds were tested in both conditions: in combination with doxycycline to assess combinatorial growth inhibition and alone to assess direct growth inhibition. The screening protocol was performed according to Ejim et al.'s (2011) antibacterial combination screening, with several adjustments (Ejim et al., 2011). Screening was carried out in 200 µl in 96-well F-bottom clear plates, in duplicate, using CAMHB with 0.1% DMSO and a library compound concentration of 10 µM. When used, the

concentration of doxycycline was 5 µg ml<sup>-1</sup>, corresponding to half minimal inhibitory concentration (MIC) value obtained under the same conditions of the monodose chemosensitization assay. The preliminary MIC assay was performed in duplicate using five shifted concentration gradients of doxycycline, from 100 to 60 µg ml<sup>-1</sup>. The 48 h MIC value obtained for *E. aerogenes* EA289 was 10 µg ml<sup>-1</sup>. The inoculum of approximately 5 × 10<sup>5</sup> CFU ml<sup>-1</sup> was prepared from an overnight plate. A single colony of *E. aerogenes*, EA289, was grown in CAMHB shaking at 250 rpm at 37°C until the OD<sub>600nm</sub> = 0.5 was reached. Background controls (eight wells per plate) contained only media and DMSO (these were also the sterility controls). Growth controls, also eight wells per plate, contained media, DMSO, and inoculum. Plates were incubated at 37°C for 48 h. The pipetting was performed by a Tecan Freedom EVO automated liquid handler equipped with eight fixed PTFE-coated tips, with a robotic arm to transfer plates to the Tecan Infinite M200 Pro multimode plate reader. Plates were read at 600 nm after 5 min shaking. Automation protocols were coded on Tecan EVOware and are available upon request.

For each test well, the percentage growth and the percentage growth inhibition were calculated as:

$$ \%Growth(x) = \frac{OD(x) - OD(background)}{OD(NC) - OD(background)} \times 100 $$

$$ \%Growth \text{ inhibition} = 100 - \%Growth(x). $$

The data set consisted in eight replicates in the absence of antibiotic and eight replicates in the presence of doxycycline. Results were represented as a XY-chart showing the mean value of eight replicates for combinatorial growth inhibition percentage, in the presence of doxycycline (+DOX) as Y-value, and the mean value of eight replicates for direct growth inhibition percentage in the absence of antibiotic (-DOX) as X-value (Figure 4).

Finally, the monodose chemosensitizer score (%MCS) was calculated as:

$$ \%MCS(x) = \%Growth \text{ inhibition}(x)_{+DOX} - \%Growth \text{ inhibition}(x)_{-DOX}. $$

Compounds showing %MCS > 40% were considered lead chemosensitizers. Results are shown in Figure 4.

The Z'-factor (Zhang et al., 1999) was calculated based on the average of eight screens, using NV845 10 µM as PC and tetracycline 10 µM as NC, in the calculation.

Hits showing 100% growth inhibition in the absence of antibiotic may be retested at one-tenth or one-hundredth the

<table><caption>TABLE 3 | Monodose Chemosensitization assay, kinetic data reduction.</caption><thead><tr><td>Endpoint measurement</td><td>Comments</td></tr></thead><tbody><tr><td>OD(x)</td><td>OD value for the growth of EA289 in the presence of a tested compound</td></tr><tr><td>OD(background)</td><td>OD value for sterile CAMHB in the presence of 0.1% DMSO</td></tr><tr><td>OD(NC)</td><td>OD value for the growth of EA289 in the presence of 0.1% DMSO</td></tr></tbody></table>
concentration, i.e., 1 or 0.1 µM, as it indicates such hit was
tested above its MIC value. Results are shown in **Supplementary
Figure 7**.

Score calculation was programmed and performed automatically by Tecan Magellan software. Both data from Magellan and pipetting instructions from EVOware were sent to a custom LIMS (Modul-Bio, Marseille, France), which matched results to compound names.

RESULTS

In order to demonstrate the efficiency of this multiparametric
profiling method as a tool for chemosensitizer discovery, an 80
control compounds library was assembled and tested, including:
reference permeabilizers, antibiotics, together with natural or
synthetic molecules (Supplementary Table 1). The first stage of
analysis consisted in identifying hit compounds against EA289
along with understanding their activity profile (Figure 2). These
assays were performed by using a multimode plate reader linked
to an automated liquid handler in order to minimize exposure
time variation of the tested compounds to bacteria, human error,
and the pipetting precision issues. We determined that a 100 µM
dose provided the most discriminatory dynamic range to screen
compounds on real-time assays. Growth and buffer conditions
were also homogenized among the assays to obtain consistent
data for profiling purposes.

**Efflux Inhibition Screening**

The effect on RND efflux was assessed using the 1,2'-
diNA fluorescent dye, a substrate of the AcrAB-TolC pump
(Figure 1B) which is fluorescent when loaded into the bacteria
(Bohnert et al., 2011). The dye efflux is then triggered by glucose
addition as energy source causing the fluorescence decrease. The
plate layout included each compound to be tested both against
EA289 and EA289Δ*acrAB*, which was included as a control since
it is unable to efflux 1,2'-diNA (Supplementary Figures 1a, 3).
To prevent biological variability among replicates, and to adjust
the overall fluorescence change due to the tested compound,
basis fluorescence, i.e., before energization, was calibrated for
each compound using the mean fluorescence obtained with
EA289Δ*acrAB*. The efflux inhibition score (%EIS) was calculated
(Figure 3A) by subtracting the final fluorescence, i.e., at the
end of the measurement, obtained in the absence of compound
from the value obtained in the presence of the tested compound.
%EIS accounts for the proportion of retained 1,2'-diNA over
the efflux capacity of EA289 without compound. Practically, a
total retention of the dye within the membrane as observed

with the *acrAB* mutant corresponded to a 100% EIS whereas
a 0% EIS was related to fully functional efflux pumps. The
data analysis included systematic signal adjustments and score
calculation resulting in a robust Z'-factor (Zhang et al., 1999):
0.71. Thus, we determined that a score exceeding 10% was
related to a significant real-time effect, leading to compound
selection as a hit. From this step of screening, among the tested
library, the following compounds showed a significant %EIS:
triclosan, CCCP, PAβN, polymyxins, polyamines, permeabilizers,
quinolines, fluoroquinolones, imipenem, and meropenem.

**Membrane Permeabilization Screening**

Membrane permeability was assessed using the measurement of
nitrocefin hydrolysis by periplasmic β-lactamases (Figure 1C). In
this case, the outer membrane permeabilization score (%OPS)
was calculated (Figure 3B) by subtracting the basal hydrolysis
area under the curve due to passive transport of nitrocefin from
the area obtained in the presence of the tested compound, in
order to prevent biological variability between replicates. This
analysis led to robust data as it takes into account both the
initial slope, the maximum plateau, and the mean absorbance,
of each curve. Finally, a bacterial lysis presenting a maximum
β-lactamase hydrolysis rate was associated to a 100% OPS and
a 0% OPS was related to a basal hydrolysis rate (background).
The data analysis included systematic signal adjustments and
score calculations resulting in a robust Z'-factor (Zhang et al.,
1999): 0.78. Thus, we determined that a score exceeding 10%
was related to a significant real-time effect, leading to compound
selection as a hit. From this step of screening, among the
tested library, the following compounds showed a significant
%OPS: triclosan, PAβN, polymyxins, polyamines, benzalkonium
chloride, chlorhexidine, CTAB, SDS, and Triton X-100. However,
since nitrocefin is a good substrate for AcrAB-TolC efflux
system (Nagano and Nikaido, 2009), we cannot exclude the
possibility that the apparent permeability enhancement was
caused by the inhibition of efflux, for example in the case of
PAβN.

**Transmembrane Potential Assay**

As proton-motive force drives RND efflux pumps, disruption
of transmembrane potential could lead to efflux impairment.
In this context, it was critical to understand if an observed
%EIS resulted in the inner membrane potential disruption
or by specifically targeting efflux systems. This parameter
was assessed with membrane potential sensitive diSC₃(5)
probe that becomes fluorescent once released into the
external medium following the disruption of the membrane
potential (Katayama et al., 2007; **Figure 1D**). Despite our

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Chemosensitizer discovery workflow.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Hit identification and profiling. The chemical library was tested at 100 µM on EA289. Compounds showing a score exceeding 10% on one tested activity were considered hit compounds. Hits showing significant %OPS attributed the permeabilizer profile. The anti-efflux profile was associated with hits showing %EIS-only. Data reduction for RND efflux inhibition screening (A) and for OM permeabilization screening (B). *, Efflux inhibition data resulted in 30% or more fluorescence quenching, data were not considered. **, %OPS was a negative value, due to β-lactamase inhibitors, or β-lactam competition with nitrocefin hydrolysis. NV845, {3-[(3-amino-propyl)-methyl-amino]-propyl}-(2E,6E)-3,7,11-trimethyl-dodeca-2,6,10-triynyl-amine; NV731, {3-[bis-(3-amino-propyl)-amino]-propyl}-(3,7-dimethyl-octa-2,6-dienyl)-amine; BG1023, N1-(7-chloroquinolin-4-yl)-N2,N2-diisopropylethane-1,2-diamine; BG1189, 3-(3-dimethylamino-propyl)-6-nitro-3H-quinazolin-4-one.</figcaption></figure>

efforts, the robustness of this assay was strongly impaired by repeated quenching effect. Under this consideration, this experiment was applied as a yes/no type of assay to provide supplementary information about the mechanism of action of the considered hit (Supplementary Table 1). The following hits showed a significant transmembrane potential disruption effect: triclosan, polymyxins, squalamine, permeabilizers, and quinolines.

## Monodose Chemosensitization Screening

The second stage of the approach consisted in investigating if the observed barrier alteration by a hit compound led to effective chemosensitization. For this purpose, selected hits were tested for growth inhibition of EA289, both alone, and in combination with an intrabacterial antibiotic. Doxycycline was selected as a probe to assess chemosensitization in EA289 since it has been previously demonstrated that, besides

target-based resistance mechanisms, acquired resistance to tetracyclines in Enterobacteriaceae was due to efflux pumps overproduction and to reduce outer membrane permeability (Thaker et al., 2009; Bohnert et al., 2010). We hypothesized that an effective chemosensitizer significantly increases growth inhibition when combined with this antibiotic.

We observed that a 10 µM dose provided the most discriminatory dynamic range on growth inhibition assays. Doxycycline was combined at a 5 µg ml⁻¹ concentration corresponding to one half the MIC against EA289. Under these conditions, the monodose chemosensitization score (%MCS) was calculated by comparing growth inhibition by the compound alone and when combined with the antibiotic. In this context, Z'-factor calculation for %MCS led to a 0.58 value, implying a more stringent selection to lead chemosensitizers, which required a %MCS exceeding 40%. Finally, identified lead chemosensitisers were NV731, NV845, PAßN, and PMB nonapeptide.
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Monodose chemosensitization assay. Anti-efflux-profiled hits (blue dots), permeabilizer-profiled hits (orange dots), were tested at 10 µM for 48 h growth inhibition of EA289. Compounds were tested in both conditions: in combination with doxycycline to assess combinatorial growth inhibition and alone to assess direct growth inhibition. Compounds showing 40% or more inhibition when combined with the intrabacterial targeted antibiotic were considered lead chemosensitizers.</figcaption></figure>

DISCUSSION

The approach presented herein, resulted in highly reproducible data, detected and efficiently profiled 24 hits (Supplementary Figure 8). Taking into account our data, we were able to discriminate four activity profiles:

* Non-specific permeabilizers showing both a significant outer membrane permeabilization activity and an efflux inhibition one, that was related to transmembrane potential disruption; including: polymyxines (Vaara, 1992), squalamine (Di Pasquale et al., 2010), and permeabilizers. We could expect compounds within this profile present low cytotoxicity as they were not specifically acting, and their use might be relevant as desinfectants. However, polymyxin and squalamine present a low cytotoxicity against eukaryotic cells, and their human therapeutic use is well established (Pogue et al., 2017).

* Specific dual-effect hits showed both a significant outer membrane permeabilization activity and an efflux inhibition one, that was independent from transmembrane potential disruption; including: PAβN, NV845, and NV731. As data will be accumulated, it would be of particular interest to discern which moiety or which chemical property participates in which activity. As new prospective strategies are currently developed to associate two molecules with compatible activities, one can consider the independent

development of two associated properties, before linkage of both partners to build hybrid active molecules (Gorityala et al., 2016).

* Specific efflux pump inhibition hits showed a significant efflux inhibition activity alone, including: fluoroquinolones, imipenem, and meropenem antibiotics. Although we hypothesized that antibiotics inhibited efflux by competition with the 1,2'-diNA dye, inhibitors by direct interaction may aggregate within this profile. The possible direct interactions may include binding with TolC causing to block the exit duct (Gilardi et al., 2017), covalent binding with AcrB causing to alter the conformational change of the three AcrB protomers therefore blocking the pump, direct interaction with AcrA causing the disruption of the AcrAB-TolC complex (Daury et al., 2016).

* Finally, non-specific efflux inhibition hits, that was related to transmembrane potential disruption, include triclosan and quinolines. Due to quenching effect, we were not able to reproduce the well-known proton gradient uncoupling effect of CCCP, but we expected it to aggregate within this profile (Venter et al., 2015).

A hit-compound selected after stage 1 was either able to permeabilize outer-membrane (by letting the β-lactam nitrocefin enter the periplasm), or inhibit AcrAB-mediated efflux. Hits showing a significant %EIS with no significant %OPS were grouped as anti-efflux whereas all other hits were grouped as permeabilizers. As both barriers have previously been demonstrated to be resistance mechanisms against many classes of antibiotics (Masi et al., 2017), one antibiotic was used to measure for the chemosensitizer activity. However, we cannot exclude this approach may select for doxycycline-specific chemosensitizers along with regular chemosensitizers.

In order to link barrier to accumulation alteration and effective chemosensitization, hits were further tested in a monodose chemosensitization assay. Although most of the lead chemosensitizers were found to be permeabilizer-profiled (Figure 4), we aimed to quantify the contribution of each barrier to be altered. From the data set composed of MCS, OPS, and EIS, a multiple linear regression analysis was performed. We hypothesized that %MCS was dependant on %OPS or %EIS and performed a multiple linear regression analysis which resulted in the following equation %MCS = 0.3 × %OPS + 0.2 × %EIS linking chemosensitization with the barriers to be impaired (Supplementary Table 3). As there was a very high chemical diversity, a wide variety of activities (efflux inhibitors, permeabilizers, dual effects, etc.), and too few data points, the statistical model remains to be improved, and the formula will be refined with the accumulation of further hits. However, the multiple linear regression model as presented has the purpose: to help define the best strategy for chemosensitization. At this point, the model tends to show the best strategy would be to develop a dual-effect chemosensitizer. As a part of broader screening campaigns this model will be refined and quantitatively respond more reliably
to the relative importance of the two barriers in resistance. In the long run, this model opens two perspectives. On the one hand, it will help quantifying the contribution of each barrier in the overall resistance to antibiotics. On the other hand, this equation may be a drug-design tool that will describe the objectives to build the best antibiotic adjuvant; along with the results of QSAR (structure vs. %OPS and structure vs. %EIS) pharmacophores will be identified.

In summary, this method originally complements current bacteriological assays as it provides new tools for quantitative structure-activity relationships in the development of new chemosensitizers against Gram-negative pathogens (Lewis, 2013; Schneider et al., 2017). Thus, this method could be applied in the early search for new molecules, prior to usual studies such as checkerboard assays including several bacterial species and combined with several intrabacterial antibiotics. Besides the screening purpose, we expect the accumulated data through the scoring tools developed therein will help elucidate the rules for Gram-negative bacteria cell entry and drug accumulation (Payne et al., 2006; Silver, 2008).

## AUTHOR CONTRIBUTIONS

All authors conceived the approach. VL developed and performed the high-throughput assays and developed the data analysis methods. VL, J-MBr, and J-MBo wrote the main manuscript text. VL prepared the figures. All authors reviewed the manuscript.

## FUNDING

This research was supported by Aix-Marseille Université, Service de Santé des Armées. The research leading to the results discussed here was conducted as part of the translocation consortium (http://www.translocation.eu) and has received support from the Innovative Medicines Initiative joint Undertaking under Grant Agreement No. 115525, resources which are composed of financial contribution from the European Union's seventh framework program (FP/2007-2013) and EFPIA companies in kind contributions. VL was supported by a French Ministry of Education and Research Fellowship.

## REFERENCES

* Bohnert, J. A., Karamian, B., and Nikaido, H. (2010). Optimized Nile Red efflux assay of AcrAB-TolC multidrug efflux system shows competition between substrates. *Antimicrob. Agents Chemother.* 54, 3770-3775. doi: 10.1128/AAC.00620-10

* Bohnert, J. A., Schuster, S., Szymaniak-Vits, M., and Kern, W. V. (2011). Determination of real-time efflux phenotypes in *Escherichia coli* AcrB binding pocket phenylalanine mutants using a 1,2'-dinaphthylamine efflux assay. *PLoS One* 6:e21196. doi: 10.1371/journal.pone.0021196

* Brunel, J. M., Lieutaud, A., Lome, V., Pagès, J.-M., and Bolla, J.-M. (2013). Polyamino geranic derivatives as new chemosensitizers to combat antibiotic

## ACKNOWLEDGMENTS

The authors thank S. Alibert for providing chemical samples of BG1023 and BG1189.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00204/full#supplementary-material

* **FIGURE S1** | Example raw data from an experiment plate, obtained for the RND Efflux inhibition screening. Compounds incubated with Ea289ΔacrAB (Δ) (a), compounds incubated with Ea289 (b,c).

* **FIGURE S2** | Example raw data from an experiment plate, obtained for the RND Efflux inhibition screening. Compounds incubated with Ea289ΔacrAB (Δ) (a), compounds incubated with Ea289 (b).

* **FIGURE S3** | RND Efflux inhibition kinetic data reduction. Data reduction for diNaphthylamine real-time efflux. Kinetics were selected from example raw data: NV845 as test and 1% DMSO as negative control conditions.

* **FIGURE S4** | Example raw data from one experiment plate, obtained for the outer membrane permeability screening.

* **FIGURE S5** | Example raw data from one experiment plate, obtained for the outer membrane permeability screening.

* **FIGURE S6** | Outer membrane permeabilization kinetic data reduction. Data reduction for outer membrane permeabilization screening. Kinetics were selected from example raw data: benzalkonium chloride as test and 1% DMSO as negative control conditions.

* **FIGURE S7** | Monodose Chemosensitization retest. Above-MIC compounds were re-tested at 1 & 0.1μM.

* **FIGURE S8** | Validation of control chemicals mode of action. Real-time assays results for selected hits were evaluated in a binary.

* **TABLE S1** | Components of the control chemicals library.

* **TABLE S2** | Transmembrane potential disruption assay. Showed transmembrane potential disruption (+), did not show transmembrane potential disruption (-), the measurement could not be assessed due to strong fluorescence quenching (q).

* **TABLE S3** | Relationship between %OPS, %EIS and %MCS by linear regression analysis. The relationship between chemosensitization and barriers to be altered was modeled in a multiple linear regression analysis. %MCS was hypothesized to be dependent on %OPS or %EIS and the analysis was performed on the previously selected hits except thioridazine and chlorpromazine as EIS could not be calculated. The full data set comprised %EIS, %OPS and %MCS for the 22 selected compounds. The MCS was calculated after retest at lower concentration as indicated previously for colistin, polymyxin B, meropenem and triclosan.

* resistant Gram-negative bacteria. *Bioorg. Med. Chem.* 21, 1174-1179. doi: 10.1016/j.bmc.2012.12.030

* Bush, K., Courvalin, P., Dantas, G., Davies, J., Eisenstein, B., Huovinen, P., et al. (2011). Tackling antibiotic resistance. *Nat. Rev. Microbiol.* 9, 894-896. doi: 10.1038/nrmicro2693

* Daury, L., Orange, F., Taveau, J. C., Verchère, A., Monlezun, L., Gounou, C., et al. (2016). Tripartite assembly of RND multidrug efflux pumps. *Nat. Commun.* 12:10731. doi: 10.1038/ncomms10731

* Di Pasquale, E., Salmi-Smail, C., Brunel, J.-M., Sanchez, P., Fantini, J., and Maresca, M. (2010). Biophysical studies of the interaction of squalamine and other cationic amphiphilic molecules with bacterial and eukaryotic membranes: importance of the distribution coefficient in membrane
selectivity. *Chem. Phys. Lipids* 163, 131–140. doi: 10.1016/j.chemphyslip.2009.10.006

Ejim, L., Farha, M. A., Falconer, S. B., Wildenhain, J., Coombes, B. K., Tyers, M., et al. (2011). Combinations of antibiotics and nonantibiotic drugs enhance antimicrobial efficacy. *Nat. Chem. Biol.* 7, 348–350. doi: 10.1038/nchembio.559

Gilardi, A., Bhamidimarri, S. P., Brönstrup, M., Bilitewski, U., Marreddy, R. K. R., Pos, K. M., et al. (2017). Biophysical characterization of *E. coli* TolC interaction with the known blocker hexaamminecobalt. *Biochim. Biophys. Acta* 1861, 2702–2709. doi: 10.1016/j.bbagen.2017.07.014

Gorityala, B. K., Guchhait, G., Goswami, S., Fernando, D. M., Kumar, A., Zhanel, G. G., et al. (2016). Hybrid antibiotic overcomes resistance in *P. aeruginosa* by enhancing outer membrane penetration and reducing efflux. *J. Med. Chem.* 59, 8441–8455. doi: 10.1021/acs.jmedchem.6b00867

Harvey, A. L., Edrada-Ebel, R., and Quinn, R. J. (2015). The re-emergence of natural products for drug discovery in the genomics era. *Nat. Rev. Drug Discov.* 14, 111–129. doi: 10.1038/nrd4510

Katayama, H., Kusaka, Y., Yokota, H., and Akao, T. (2007). Parasporin-1, a novel cytotoxic protein from *Bacillus thuringiensis*, induces Ca²⁺ influx and a sustained elevation of the cytoplasmic Ca²⁺ concentration in toxin-sensitive cells. *J. Biol. Chem.* 282, 7742–7752. doi: 10.1074/jbc.M611382200

Kristiansen, J. E., Hendricks, O., Delvin, T., Butterworth, T. S., Aagaard, L., Christensen, J. B., et al. (2007). Reversal of resistance in microorganisms by help of non-antibiotics. *J. Antimicrob. Chemother.* 59, 1271–1279. doi: 10.1093/jac/dkm071

Lewis, K. (2013). Platforms for antibiotic discovery. *Nat. Rev. Drug Discov.* 12, 371–387. doi: 10.1038/nrd3975

Lomovskaya, O., Zgurskaya, H. I., Totrov, M., and Watkins, W. J. (2006). Waltzing transporters and “the dance macabre” between humans and bacteria. *Nat. Rev. Drug Discov.* 6, 56–65. doi: 10.1038/nrd2200

Mallea, M., Chevalier, J., Bornet, C., Eyraud, A., Davin-Regli, A., Bollet, C., et al. (1998). Porin alteration and active efflux: two in vivo drug resistance strategies used by *Enterobacter aerogenes*. *Microbiology* 144(Pt 11), 3003–3009. doi: 10.1099/00221287-144-11-3003

Masi, M., Refregiers, M., Pos, K. M., and Pagès, J.-M. (2017). Mechanisms of envelope permeability and antibiotic influx and efflux in Gram-negative bacteria. *Nat. Microbiol.* 2:17001. doi: 10.1038/nmicrobiol.2017.1

Matsumoto, Y., Hayama, K., Sakakihara, S., Nishino, K., Noji, H., Iino, R., et al. (2011). Evaluation of multidrug efflux pump inhibitors by a new method using microfluidic channels. *PLoS One* 6:e18547. doi: 10.1371/journal.pone.0018547

Mazumdar, K., Dastidar, S. G., Park, J. H., and Dutta, N. K. (2009). The anti-inflammatory non-antibiotic helper compound diclofenac: an antibacterial drug target. *Eur. J. Clin. Microbiol. Infect. Dis.* 28, 881–891. doi: 10.1007/s10096-009-0739-z

Nagano, K., and Nikaido, H. (2009). Kinetic behavior of the major multidrug efflux pump AcrB of *Escherichia coli*. *Proc. Natl. Acad. Sci. U.S.A.* 106, 5854–5858. doi: 10.1073/pnas.0901695106

Nikaido, H. (1994). Prevention of drug access to bacterial targets: permeability barriers and active efflux. *Science* 264, 382–388. doi: 10.1126/science.8153625

Payne, D. J., Gwynn, M. N., Holmes, D. J., and Pompliano, D. L. (2006). Drugs for bad bugs: confronting the challenges of antibacterial discovery. *Nat. Rev. Drug Discov.* 6, 29–40. doi: 10.1038/nrd2201

Pogue, J. M., Ortwine, J. K., and Kaye, K. S. (2017). Clinical considerations for optimal use of the polymyxins: a focus on agent selection and dosing. *Clin. Microbiol. Infect.* 23, 229–233. doi: 10.1016/j.cmi.2017.02.023

Pradel, E., and Pagès, J. M. (2002). The AcrAB-TolC efflux pump contributes to multidrug resistance in the nosocomial pathogen *Enterobacter aerogenes*. *Antimicrob. Agents Chemother.* 46, 2640–2643. doi: 10.1128/AAC.46.8.2640-2643.2002

Schneider, E. K., Reyes-Ortega, F., Velkov, T., and Li, J. (2017). Antibiotic-non-antibiotic combinations for combating extremely drug-resistant Gram-negative 'superbugs'. *Essays Biochem.* 2017, 115–125. doi: 10.1042/EBC2016 0058

Silver, L. L. (2008). Are natural products still the best source for antibacterial discovery? The bacterial entry factor. *Expert Opin. Drug Discov.* 3, 487–500. doi: 10.1517/17460441.3.5.487

Silver, L. L. (2011). Challenges of antibacterial discovery. *Clin. Microbiol. Rev.* 24, 71–109. doi: 10.1128/CMR.00030-10

Spellberg, B., and Shlaes, D. (2014). Prioritized current unmet needs for antibacterial therapies. *Clin. Pharmacol. Ther.* 96, 151–153. doi: 10.1038/clpt.2014.106

Stavenger, R. A., and Winterhalter, M. (2014). TRANSLOCATION project: how to get good drugs into bad bugs. *Sci. Transl. Med.* 6:228ed7. doi: 10.1126/ scitranslmed.3008605

Thaker, M., Spanogiannopoulos, P., and Wright, G. D. (2009). The tetracycline resistome. *Cell. Mol. Life Sci.* 67, 419–431. doi: 10.1007/s00018-009-0172-6

Tommasi, R., Brown, D. G., Walkup, G. K., Manchester, J. I., and Miller, A. A. (2015). ESKAPEing the labyrinth of antibacterial discovery. *Nat. Rev. Drug Discov.* 14, 529–542. doi: 10.1038/nrd4572

Vaara, M. (1992). Agents that increase the permeability of the outer membrane. *Microbiol. Rev.* 56, 395–411.

Venter, H., Mowla, R., Ohene-Agyei, T., and Ma, S. (2015). RND-type drug efflux pumps from Gram-negative bacteria: molecular mechanism and inhibition. *Front. Microbiol.* 28:377. doi: 10.3389/fmicb.2015.00377

Zgurskaya, H. I., López, C. A., and Gnanakaran, S. (2015). Permeability barrier of Gram-negative cell envelopes and approaches to bypass it. *ACS Infect. Dis.* 1, 512–522. doi: 10.1021/acsinfecdis.5b00097

Zhang, J., Chung, T., and Oldenburg, K. A. (1999). Simple statistical parameter for use in evaluation and validation of high throughput screening assays. *J. Biomol. Screen.* 4, 67–73. doi: 10.1177/108705719900400206

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Lôme, Brunel, Pagès and Bolla. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.