- 1Instituto de Investigação e Inovação em Saúde (i3S), Universidade do Porto, Porto, Portugal
- 2Instituto de Biologia Molecular e Celular (IBMC), Universidade do Porto, Porto, Portugal
- 3Programa Doutoral em Biologia Molecular e Celular (MCbiology), Instituto de Ciências Biomédicas Abel Salazar (ICBAS), Universidade do Porto, Porto, Portugal
The bacterial K+ homeostasis machinery is widely conserved across bacterial species, and different from that in animals. Dysfunction in components of the machinery has an impact on intracellular turgor, membrane potential, adaptation to changes in both extracellular pH and osmolarity, and in virulence. Using a fluorescence-based liposome flux assay, we have performed a high-throughput screen to identify novel inhibitors of the KtrAB ion channel complex from Bacillus subtilis, a component of the K+ homeostasis machinery that is also present in many bacterial pathogens. The screen identified 41 compounds that inhibited K+ flux and that clustered into eight chemical groups. Many of the identified inhibitors were found to target KtrAB with an in vitro potency in the low μM range. We investigated the mechanisms of inhibition and found that most molecules affected either the membrane component of the channel, KtrB alone or the full KtrAB complex without a preference for the functional conformation of the channel, thus broadening their inhibitory action. A urea derivative molecule that inhibited the membrane component of KtrAB affected cell viability in conditions in which KtrAB activity is essential. With this proof-of-concept study, we demonstrate that targeting components of the K+ homeostasis machinery has the potential as a new antibacterial strategy and that the fluorescence-based flux assay is a robust tool for screening chemical libraries.
Introduction
Potassium is the most abundant inorganic cation in cells from all organisms. In bacteria, K+ homeostasis is related to intracellular turgor and membrane potential, and is involved in the mechanisms of adaptation to environmental stress such as changes in extracellular pH or osmolarity (reviewed by Epstein, 2003). For example, during adaptation to hyper-osmotic shock the concentration of K+ rises from ~300 to 600–700 mM (Richey et al., 1987; Whatmore et al., 1990).
The K+ homeostasis machinery in bacteria includes only a few types or families of K+ importers and exporters. An important family of bacterial K+ importers is the superfamily of potassium transporters (SKT), which are fast ion importers present in all kingdoms of life except animals. This superfamily includes the subgroup of prokaryotic TrK/Ktr channels, which show low selectivity for K+ (Cao et al., 2013; Levin and Zhou, 2014; Diskowski et al., 2015). SKT members evolved from a simpler bacterial KcsA-like channel, with which they share the core structural motif of transmembrane helix-pore-transmembrane helix (TM-PLoop-TM), with the pore segment including a selectivity filter (Durell and Guy, 1999; Durell et al., 1999; Corratge-Faillie et al., 2010). The ion conducting membrane part of prokaryotic TrK/Ktr channels, TrKH/KtrB, is formed by a dimer of four TM-PLoop-TM repeats assembled around a central pore. Crystal and cryo-EM structures of Trk from Vibrio parahaemolyticus and Ktr complexes from Vibrio alginolyticus and Bacillus subtilis show a dimeric membrane protein assembled with a cytoplasmic tetrameric or octameric gating ring composed by the regulators of conductance of K+ (RCK) domains TrkA and KtrA, respectively (Cao et al., 2013; Vieira-Pires et al., 2013; Diskowski et al., 2017). TrkA and KtrA are nucleotide-dependent RCK domains regulated by ATP and ADP. Binding of the nucleotides induces conformational changes in the gating ring that propagate to the transmembrane subunit and lead to an active (ATP-bound) or inactive (ADP-bound) channel (reviewed by Schrecker et al., 2019). In some bacteria and archaea, cyclic-di-AMP also binds to these RCK domains and is proposed to inhibit K+ flux (Corrigan et al., 2013; Bai et al., 2014; Kim et al., 2015; Commichau et al., 2018, 2019; Rocha et al., 2019). In addition it has been demonstrated that cyclic-di-AMP binds to the ydaO riboswitch in the ktrAB operon of B. subtilis inhibiting transcription and therefore regulating protein levels in the cell (Nelson et al., 2013; Gundlach et al., 2017).
Members of the TrK/Ktr subfamily of SKT have been implicated in resistance to antibiotics. Screening of a library of insertion mutations in Mycobacterium smegmatis revealed mutations in TrkA that resulted in changes to the proton motive force (PMF) and reduced bacterial susceptibility to antibiotic challenge (Castaneda-Garcia et al., 2011). Escherichia coli strains adapted to various aminoglycoside antibiotics revealed mutations in TrkH. It was suggested that resistance to the aminoglycosides resulted from an altered function of TrkH that affected the PMF (Lazar et al., 2013; Oz et al., 2014). Another TrkH point mutation was identified in a highly streptomycin-resistant Salmonella enterica strain (Wistrand-Yuen et al., 2018). In addition, a point mutation in TrkH was identified in a Streptococcus pneumonia pneumococcal isolate resistant to non-aminoglycoside antibiotics (Pan et al., 2018), and KtrA was found to confer resistance to aminoglycosides and fitness advantage during infection in Staphylococcus aureus (Gries et al., 2013).
Strikingly, TrK/Ktr proteins have also been implicated in virulence of pathogenic bacteria. In a Campilobacter jejuni transposon insertion screening, KtrA and KtrB were identified as virulence factors for mouse intestinal colonization. The role of KtrA in virulence was confirmed by observing that its knockout mutant was impaired in colonization, and reintroduction of the gene restores colonization capacity (Gao et al., 2017). TrK transposon mutants arose in two different screenings for virulence factors on the pathogenic bacteria Aeromonas hydrophila and Pectobacterium wasabie, and for both their role in virulence was confirmed (Valente and Xavier, 2016; Pang et al., 2017). In addition, inactivation of TrK in Mycobacterium tuberculosis leads to attenuation of host colonization (MacGilvary et al., 2019), and deletion of TrkA in S. enterica leads to attenuation of cell invasion and animal infection, which further increases with additional deletion of the other K+ importers Kup and Kdp (Liu et al., 2013). The function of K+ importers in virulence probably results from their role in adaptation mechanisms that are crucial for the ability of the pathogen to survive in different environments in the host.
In this study, we explore the viability of screening for novel inhibitors of the SKT proteins by using KtrAB from B. subtilis as a model channel. This is an extensively studied protein, with its structural and functional properties well characterized, both in vitro and in vivo (Holtmann et al., 2003; Albright et al., 2006; Vieira-Pires et al., 2013; Szollosi et al., 2016; Gundlach et al., 2017; Teixeira-Duarte et al., 2019). We have adapted an established fluorescence-based liposome flux assay (Su et al., 2016; Teixeira-Duarte et al., 2019) to monitor K+ efflux through KtrAB into a 96-well plate format and executed a proof-of-concept screening campaign in which we screened a diverse library with 10,000 compounds. In addition, we have established a secondary lead follow-up screening strategy in which we have identified a urea-based derivative that has a specific effect on channel activity.
Materials and Methods
Cloning, Protein Expression, Purification, and Incorporation Into Liposomes
Bacillus subtilis KtrAB, KtrAB with R417K mutation in KtrB, and KtrB proteins were expressed, purified, and incorporated into liposomes as previously described (Teixeira-Duarte et al., 2019), with the following modifications: EDTA was absent from KtrA purification, KtrB expression was induced for 3 h, KCl and EDTA were absent from the KtrB lysis buffer, and then liposomes were subjected to ultracentrifugation at 267,000 × g for 25 min at 4°C. The liposomes were resuspended in half the initial volume for concentration before flash freezing in liquid nitrogen.
Bacillus cereus NaK DNA was a kind gift of Y. Jiang. The N-terminally truncated NaKΔN19 was cloned into pQE60 by Gibson Assembly (Gibson et al., 2009) with appropriate primers (TAAGCTTAATTAGCTGAGCTTGGACTCC and GTTAATTTCTCCTCTTTAATGAATTCTGTGTGAAATTGTTATC to amplify pQE60, and AATTCATTAAAGAGGAGAAATTAACATGTGGAAAGATAAAGAATTTCAAGTATTATTTGTATTAACAATTTTGAC and GATCTATCAACAGGAGTCCAAGCTCAGCTAATTAAGCTTATTAGTGATGGTGATGGTGATGAGATCTG to amplify NaK∆N19), and the F92A mutation was introduced by site-directed mutagenesis with primers TTATTGGGATTGGACTAGTGGCTGGATTTATTCATAAGTTAGCAG and CTGCTAACTTATGAATAAATCCAGCCACTAGTCCAATCCC. NaKΔN19/F92A protein was expressed and purified as described (Alam and Jiang, 2009). Incorporation of B. cereus NaKΔN19/F92A protein into liposomes was performed as for B. subtilis KtrAB.
K+ Flux Fluorescence-Based 96-Well Assay and Compound Screen
A screening campaign was conducted using a Chembridge DiverSet library (Chembridge, San Diego, CA, United States) consisting of 10,000 compounds chosen for their structure diversity. Compounds were stored at −20°C in 1 mM in DMSO in a 384-well format.
The K+ flux fluorescence-based screening assay was adapted from Su et al. (2016) and Teixeira-Duarte et al. (2019). It is based on the quenching of the pH gradient-sensitive dye 9-amino-6-chloro-2-methoxyacridine (ACMA) due to a proton influx that is coupled to the K+ efflux through KtrAB. Proteoliposomes in HEPES/NMG Buffer pH8.0 [10 mM HEPES, 7 mM NMG buffer (pH 8),150 mM KCl, 0.1 mM ATP, and 0.5 mM MgCl2] were 100-fold diluted in the same buffer containing 0.56 μM ACMA and 50 mM choline chloride instead of KCl. In this assay format, we used the lowest ACMA concentration to minimize potential compound interference. The mixture was added to 96-well black, flat-bottomed microplates (Greiner Bio-One, Monroe, NC, United States), 100 μl per well, and the plate was centrifuged at 1,000 rpm for 1 min in a bench top centrifuge. A Janus automated Workstation (PerkinElmer, Waltham, MA, United States) equipped with a pin tool (V&P Scientific, San Diego, CA, United States) was used to transfer 0.1 μl of compound from the 384-well library plates to the 96-well plates already containing the proteoliposome mixture such that the assay contained a final concentration of 1 μM compound and 0.1% DMSO. Wells in column 12 of each assay plate contained no compound and were used as negative controls representing no inhibition. The plates were transferred to a Synergy MX plate reader (BioTek, Winooski, VT, United States) and the fluorescence of ACMA was recorded over time at 490/20 nm with excitation 410/20 nm, for 2 min 40 s, using the maximum number of measurements per data point; in this way, we maximized read accuracy without excessively prolonging total assay time, which could compromise liposome stability. The assay was started with the addition of the protonophore carbonyl cyanide m-chlorophenyl hydrazine (CCCP) to a final concentration of 1 μM using a Multidrop Combi dispenser (Thermo Scientific, Waltham, MA, United States), and the fluorescence was measured every 12 s for 5 min, a time point at which quenching steady-state has been reached in the absence of inhibition (Teixeira-Duarte et al., 2019).
To measure K+ flux mediated by the KtrABR417K-ADP complex ADP replaced ATP in the liposome preparation and in the assay buffer. To measure K+ flux from KtrB or NaKΔN19/F92A, nucleotides were absent from the assay buffer, ACMA concentration was 1.12 μM, and the KtrB assay was run for 15 min.
K+ Transport Fluorescence-Based Cuvette Assay
KtrB liposomes were prepared as described above except that they were not concentrated by ultracentrifugation. They were incubated with KtrA at a ratio 1:0.6 KtrB:KtrA (w:w) for at least 1 h at room temperature, diluted to a final volume of 2 ml in buffer containing ACMA, as above, plus 2 mM MgCl2 and 150 mM choline chloride. Measurements were done in polymethacrylate cuvettes (Merck, Germany) with a PTFE magnetic stirrer (ϕ 9 × 6 mm). Fluorescence at 480 nm, upon excitation at 410 nm, was recorded on a FluoroMax-4 (Horiba, Edison, NJ, United States) spectrofluorometer in 2 s intervals, with constant stirring. The assay was started by the addition of CCCP as above, and 300 nM valinomycin (Val) was added after 500 s to estimate the fraction of active liposomes.
Screening Analysis, Statistics, and Retesting
Gen5 software (BioTek) was used to analyze all screening results. Two different types of analyses were performed: initial velocity and final fluorescence intensity. Z’ values were calculated within each assay plate and were found to be 0.1–0.7 (mean 0.4) for the initial velocities, and 0.7–0.95 (mean 0.9) for the final fluorescence intensities, indicating acceptable and excellent quality screens (Zhang et al., 1999), respectively. Fluorescence curves that gave changes of 2× SD from the initial velocity mean or 3× SD from the final fluorescence intensity mean were individually inspected. Clear false positives, generally due to compound fluorescence, were discarded, and the remaining compounds were retested at 1 μM in triplicate using samples from the original screening library plates and 12 negative controls wells per plate. Analysis for confirmation of these screening hits was done as for the initial screening, using the mean of the negative controls.
Dose-Response
Compounds with confirmed activity were ordered from Chembridge, Hit2Lead and dose-response assays were performed as described above, in the presence of 50 (in plate format) or 150 mM (in cuvette) choline chloride, in triplicate, with compound concentrations ranging from 8 nM to 128 μM (depending on the compound), and at least 8 negative controls per plate. To test if compounds destabilized liposomes leading to false positive effects, empty liposomes were also assayed in the presence of high concentrations of compounds, by adding the K+ ionophore Val to a final concentration of 1 μM at the end of the assay. If, in the presence of a compound concentration, the decrease in fluorescence after Val addition was smaller than in the absence of compound, it was concluded that at this concentration, the compound destabilized liposomes, dissipating the K+ gradient, and it was removed from the dose-response analysis of proteoliposomes.
For data collected in plate format, dose-response curves were generated using final fluorescence intensities, by calculating the percentage of inhibition relative to the final fluorescence intensity before addition of CCCP.
For data collected in cuvette format, fluorescence curves were normalized using two different procedures: (1) using as maximum the initial fluorescence (before CCCP addition) and as minimum the fluorescence after Val addition and (2) using as maximum the initial fluorescence (before CCCP addition) and as minimum the fluorescence read immediately before Val addition. Initial velocities were determined from normalized curves as the slope of the straight line that fitted the first (at least five) time points in the curve, after addition of CCCP. Dose-response curves were generated using initial velocities, either by plotting the percentage of inhibition as a fractional change of maximum initial velocity, obtained in the absence of compound or by plotting the velocity value in function of compound concentration.
All dose-response curves were fitted with a Hill equation,
in which A is the minimum inhibition or maximum velocity, B is the maximum inhibition or minimum velocity, x is the compound concentration, KI is the constant of 50% inhibition, and nH is the Hill coefficient.
KtrA Thermal Shift Assay
Forty microliters of KtrA protein at a final concentration of 2.4 μM in HEPES/NMG buffer pH8.0 containing 5 mM DTT, and Sypro Orange Dye (Sigma-Aldrich, St. Louis, MO, United States) at a final concentration of 2.5×, were mixed with 10 μl of compound at a final concentration of 10 μM or 0.1% DMSO, as control. Assay was run in a 96-well white microplate (Bio-Rad, Hercules, CA, United States), in triplicate. After a centrifugation in a bench top centrifuge at 1,500 rpm for 1 min, plates were heated from 25 to 75°C in 0.5°C increments in an IQ5 Real Time Detection System (Bio-Rad) with a pause of 30 s at each temperature. Unfolding of KtrA was monitored using excitation and emission filters 545/30 nm and 585/20 nm, respectively.
Growth Curves of B. subtilis
All B. subtilis cultures were performed in modified Spizizen’s minimal medium, containing 2 g/L (NH4)2SO4, 1 g/L sodium citrate, 0.2 g/L MgSO4·7H2O, 5 g/L glucose, 50 mg/L tryptophan, 50 mg/L phenylalanine, 125 mg/L MgCl2·6H2O, 7.3 mg/L CaCl2·2H2O, 13.5 mg/L FeCl2·6H2O, 1 mg/L MnCl2·4H2O, 1.7 mg/L ZnCl2, 0.43 mg/L CuCl2·2H2O, 0.3 mg/L CoCl2, and 0.6 mg/L NaMoO4·2H2O. Concentrations of phosphate salts (K2HPO4, KH2PO4, Na2HPO4, and NaH2PO4) were adjusted for the appropriate K+ concentration, keeping the total concentration of the monovalent cations (Na+ and K+) as 150 mM.
Bacillus subtilis JH642 (pheA1 trpC2 citS642; BGSCID 1A96) results from an 18-kb deletion on B. subtilis 168 (trpC2; BGSCID 1A1) genome (Srivatsan et al., 2008). B. subtilis 168 ∆ (ktrB::kan) strain, in which the ktrB gene is disrupted by a kanamycin resistance cassette, was constructed as follows. Genomic DNA of B. subtilis 168 ∆ (ktrB::kan) (obtained from the Bacillus Genetic Stock Center) was isolated with Gene JET Genomic DNA Purification kit (Thermo Scientific). The kanamycin resistance cassette flanked by KtrB sequences was amplified with primers GGCGGCAGCATTTGTAAAGAAC and CGCCTCTATTAAAGACAGTGCCGTTC. B. subtilis 168 was then transformed with the amplified DNA fragment according to Yasbin et al. (1975). To avoid suppressor mutations due to low K+ concentrations, positive transformants were selected in the presence of 10 g/L KCl, in addition to 7.5 μg/ml kanamycin.
All strains were grown in the presence of 1 mM K+ for 6 h at 37°C before being inoculated into 96-well sterile flat-bottomed microplates (Corning, Durham, NC, United States), at an initial OD600 0.05, in a total volume of 100 μl supplemented with the appropriate concentration of compounds. Plates were covered with Breathe-Easy sealing tapes (Diversified Biotech, Dedham, MA, United States), which allow gas exchange, placed in a Synergy 2 or MX plate reader for 15 h at 37°C with medium agitation, and the OD600 was recorded every 30 min.
Results
Initial Screening and Retest Confirmation
To identify regulators of K+ flux through the KtrAB channel, we screened 10,000 compounds at a final compound concentration of 1 μM from the Chembridge DiverSet library, using a modified fluorescence-based liposome flux assay (Su et al., 2016; Teixeira-Duarte et al., 2019) in a 96-well format (Figures 1A,B). In this assay, the ATP-bound KtrAB complex is reconstituted into liposomes in the presence of a high concentration of potassium (150 mM KCl). This reconstitution method allows both orientations of the complex, and both orientations are functional in the assay. The proteoliposomes are then diluted in a K+-free buffer containing ATP, creating a K+ gradient that drives K+ efflux through KtrAB. K+ efflux is counter balanced by H+ influx facilitated by the proton ionophore CCCP, creating a pH gradient that results in quenching of ACMA fluorescence. We followed quenching of ACMA over time after the addition of CCCP, and determined two parameters, the initial velocity and the fluorescence intensity at a final time-point, both of which reflect KtrAB activity.
Figure 1. Screening assay. (A) Scheme of reconstituted KtrAB proteoliposomes in which quenching of 9-amino-6-chloro-2-methoxyacridine (ACMA) fluorescence is monitored in response to fluxes of K+ and H+ ions driven by KtrAB and CCCP. KtrAB is randomly oriented in liposomes and it is likely that more than one membrane protein complex are present in each liposome. (B) Example of a curve of ACMA fluorescence obtained from a 96-well plate format. Carbonyl cyanide m-chlorophenyl hydrazine (CCCP) and Val addition are represented by arrows. See text for a detailed explanation of the assay. (C) Scheme of the high throughput screening methodology followed in this study.
To evaluate the non-specific impact of compounds on liposome integrity, we added the K+ ionophore Val at the end of the assay, which results in the total dissipation of the K+ gradient and rapid fluorescence quenching to a baseline level if liposomes are intact. In addition, quenching curves can be normalized to the total fluorescence quenching after Val addition. This ionophore was only included at the single dose retesting and dose-response steps of the screen.
In the initial screen, compounds were defined as “hits” if either the initial velocity diverged 2× SDs from the mean initial velocity in the absence of compound or the final fluorescence intensity diverged 3× SDs from the mean final fluorescence intensity in the absence of compound. This relatively wide cutoff was deliberately used so that we would not “miss” any potential hit compounds. Using these criteria, 171 compounds were identified resulting in an overall hit rate of 1.7% (Figure 1C; Supplementary Figure S1). The initial hits were retested in triplicate from liquid samples at 1 μM and 90 compounds retest-confirmed, resulting in a 52% retest confirmation rate.
Dose-Response Confirmation
Seventy of the 90 compounds that retest-confirmed were arranged into eight structural groups defined by heterocyclic ring backbone structure. Forty-two of these compounds were selected based on their structure and percentage of inhibition, and re-ordered as powders from Chembridge Hit2Lead for dose-response confirmation (Figure 1C). From those, 41 compounds showed a dose-response based on final fluorescence intensities and in 40 cases the dose-response could be fit with a Hill equation, deriving a constant of inhibition (KI) and Hill coefficient (nH) (Supplementary Table S1). In each case, we excluded from the dose-response curves the compound concentrations that had a deleterious effect on liposome integrity, observed when measuring the response of empty liposomes (without KtrAB) to Val in the presence of compound (Supplementary Figure S2; Supplementary Table S2). Representative dose-responses are shown in Figure 2A, for compounds #21 and #39 and Supplementary Figure S3 shows all dose responses. We further validated our plate-based assay format by repeating the dose-response of these compounds in a more rigorous single cuvette format. The cuvette data were analyzed using two approaches, by determining initial velocities from Val -normalized curves (Figure 2B) and by determining initial velocities from curves normalized to the final time-point on the segment after addition of CCCP only (Supplementary Figure S4). The KI values determined varied less than 2 fold in magnitude and were consistent in both assay formats and analysis method used.
Figure 2. Dose-response of compounds #21 and #39. (A, left) – Representative experiments of K+ flux assays performed in 96-well plate with the indicated concentrations of compound #21 (top) or #39 (bottom). (A, right) – Percentage of inhibition of KtrAB activity determined using the final fluorescence from curves (on the left) as function of compound concentration, for compound #21 (top) or #39 (bottom). Mean ± standard deviation (SD) of triplicates are shown fitted with Hill equation and indicated nH and KI ± standard error (SE) parameters. See Supplementary Figure S3 for dose-response curves of all compounds. (B, left) – Representative experiments of K+ flux assays performed in cuvette with the indicated concentrations of compound #21 (top) or #39 (bottom). Curves were normalized using as maximum the initial fluorescence (before CCCP addition) and as minimum the fluorescence read after valinomycin (Val) addition. (B, right) – Percentage of inhibition of KtrAB activity determined using initial velocities determined from curves (on the left) as function of compound concentration, for compound #21 (top) or #39 (bottom). Mean ± SD of triplicates are shown fitted with Hill equation and indicated nH and KI ± SE parameters. Arrows indicate addition of CCCP (first) and Val (second).
Investigating the mechanism of compound inhibition
Many compounds showed a dose-response with KI ~ 1 μM and were distributed throughout the different structural groups identified (Supplementary Table S1). Therefore, we decided to characterize in more detail nine representative compounds, outlined in Table 1 (compounds #21, #6, #31, #25, #3, #36, #39, #1, and #41), focusing on the following effects: binding site; state-dependent binding; binding specificity, and in vivo growth phenotype.
In order to evaluate whether the compounds affect KtrAB K+ flux through binding to the cytoplasmic KtrA protein, we measured the thermal stability of purified KtrA protein in the presence or absence of 10 μM of each compound. None of them significantly changed the melting temperature of KtrA (Tm = 36.0 ±0.5°C) suggesting that they do not bind to KtrA (data not shown). Due to incompatibility of the thermal shift assay with detergents we could not use this format to evaluate compound effect on the KtrAB complex.
To investigate if the compounds affect KtrAB K+ flux directly through the membrane protein KtrB, we performed the K+ flux assay in liposomes containing the membrane subunit KtrB only. In the absence of the regulatory octameric KtrA ring, K+ flux is reduced 10-fold but a small K+ efflux can still be observed (Teixeira-Duarte et al., 2019). We, therefore, tested the functional effect of one compound from each structural group (Table 1) on KtrB alone (Figure 3). We observed inhibition of K+ efflux with the majority of compounds tested, but quantification of the effect was difficult due to the low activity of KtrB (Figure 3). Interestingly, compound #36, a urea derivative, showed a significantly higher potency of inhibition on KtrB K+ efflux (KI = 0.5 μM) compared to that seen on the KtrAB complex (KI = 8 μM).
Figure 3. Inhibition of KtrB by representative compounds. Percentage of inhibition of KtrB or KtrAB activity determined using the final fluorescence from flux curves vs. compound concentration, for the indicated compounds. Mean ± SD of triplicates are shown. KtrAB curves were built with the parameters shown in Supplementary Table S1, and KtrB data were fitted with Hill equations with the following parameters: cpd#39, KI = 3.8 μM, nH = 0.92; cpd#36, KI = 0.5 μM, nH = 2.3.
We have previously shown that K+ flux through KtrAB is conformation or state-dependent (Vieira-Pires et al., 2013; Szollosi et al., 2016; Teixeira-Duarte et al., 2019). Up to this point, we have been assaying compound inhibition of K+ flux through the KtrAB complex in the active or open state, in the presence of ATP. In order to understand if compound binding and inhibitory effect is limited to one particular conformation, we tested the capacity of the compounds also to inhibit the inactive conformation, in the presence of ADP. This was achieved by measuring K+ flux using a previously characterized mutant of KtrAB, KtrABR417K, in the presence of ADP. The KtrABR417K-ADP mutant is thought to still acquire the inactive conformation of wild type KtrAB-ADP while exhibiting a higher K+ flux, making it more amenable to assay (Szollosi et al., 2016). As shown in Figures 4A,B, all tested compounds showed inhibition with both the active and inactive forms of KtrABR417K leading us to conclude that binding of the inhibitors is not state-dependent.
Figure 4. Inhibition of K+ flux in the ATP (open) and ADP (closed) forms of KtrABR417K by representative compounds. (A) – Representative curves of K+ flux in the presence or absence of 4 μM of cpd #1, in proteoliposomes prepared in the presence of ATP or ADP, as indicated. Arrows indicate addition of CCCP (first) and Val (second). (B) – Percentage of inhibition of activity by each compound determined using final fluorescence relative to activity determined in the absence of compound. Compound concentrations used: 4 μM of #1, 6 μM of #36, and 8 μM of the remaining. Mean ± SD of triplicates are shown.
To evaluate the specificity of compounds for the KtrAB channel, we chose to look at the inhibitory effect of the same subset of nine representative compounds on the NaK cation channel from B. cereus. This protein does not belong to the SKT superfamily; however, it does have a similar membrane structural architecture. We used the NaKΔN19/F92A mutant form of the channel in our flux liposome assay, as it has been reported to show higher activity than the wild type form (Alam and Jiang, 2009). As shown in Figure 5, the majority of compounds had lower efficacy against the NaKΔN19/F92A channel (with at least 10 × higher KI), and compound #39 showed no inhibition within the tested concentration range. This data suggests that these compounds are specific in their action on the KtrAB channel. Notably, compound #36, which affected KtrAB activity through the membrane protein, KtrB, also showed a similar inhibitory effect on the NaK channel.
Figure 5. Comparison of inhibition of KtrAB and NaK by representative compounds. Percentage of inhibition of NaK or KtrAB activity was determined using the final fluorescence of flux curves in function of compound concentration, for the indicated compounds. Mean ± SD of triplicates are shown. KtrAB curves were built with the parameters shown in Supplementary Table S1, and NaK data were fit with a Hill equation with the following parameters: #1, KI = 14 μM, nH = 1.4; #6, KI = 16 μM, nH = 1.7.
Impact on Cell Viability
To investigate the action of these identified compounds on cell viability through disruption of K+ homeostasis, we utilized the fact that KtrAB is required for B. subtilis growth at low K+ concentrations (Holtmann et al., 2003). We expected that compounds that inhibited B. subtilis KtrAB activity in liposomes would also inhibit B. subtilis cell growth in low K+ (1 mM) but not in high K+ (30 mM). Initially, we tested the B. subtilis strain JH642 as it is devoid of the KimA K+ importer (Gundlach et al., 2017) and as such is dependent on KtrAB for growth in low K+ conditions. This strain was grown in low and high K+ media in the presence of different concentrations of the compounds that showed an in vitro inhibitory dose-response with KI ≤ 10 μM (34 compounds). Beyond a certain concentration threshold, we observed complete inhibition of growth with all compounds that may also result from off target effects. We focused on intermediate concentrations where we could look for a K+ specific growth phenotype. Three compounds markedly showed an inhibition of growth at 1 mM K+ compared with 30 mM K+: #37, belonging to the benzimidazol group, #18 from the piperidine group, and #36 from the urea group (Figure 6A; Supplementary Figures S5, S6). The inhibitory effect is modest, with high concentrations of compound (100–500 μM) required.
Figure 6. Impact of selected compounds on Bacillus subtilis viability. (A) – Growth curves of B. subtilis JH642 strain at 1 or 30 mM K+ for the nine compounds characterized in vitro (Table 1) plus other compounds that showed a differential effect, in the absence (control) or presence of the compound (cpd) at the following concentrations: 500 μM of compounds #37 and #18, 300 μM of #36, 300 μM of #6, 100 μM of #31, #25,#3, #39, and #41, 50 μM of #1, and 10 μM of #21. Examples shown correspond to compound concentrations where the clearest effect was observed without total cell arrest, which invariably occurred at higher compound concentrations (Supplementary Figure S5). Supplementary Figure S6 shows growth curves for the remaining 23 compounds tested. (B) – Growth curves of B. subtilis 168 strain for compounds #37 and #18 under the same conditions as in A whereas compound #36 was 500 μM. (C) – Growth curves of B. subtilis 168 and 168 Δ(ktrB::kan) strains in the presence of compound #36 at the indicated concentrations. In all panels, mean ± SD of triplicates or quadruplicates are shown, except for compound #37 in (A), which shows the results of two independent experiments.
To further test the efficacy of compounds #37, 18 and 36 we also evaluated their impact on the parental B. subtilis strain 168, which, in addition to the KtrAB channel contains the KimA K+ importer. Compound #36 maintained the same phenotype, showing an extended delay in cell growth in 1 mM K+ that was partially compensated in 30 mM K+ (Figure 6B), while the effect disappeared for compounds #37 and #18. Moreover, the amplitude of the K+ differential effect observed for compound #36 was concentration dependent, increasing with compound concentration (Figures 6B, left panel of C). The differential effect was abolished in a mutant strain that lacks KtrAB, strain 168 ∆ (ktrB::kan), confirming that compound #36 is targeting KtrAB (Figure 6C, right panel).
Discussion
We have successfully executed a fluorescence-based liposome K+ flux assay to identify novel molecules that inhibit the KtrAB channel from B. subtilis. First, this liposome assay is very robust; it is easily adapted to a plate format for screening campaigns, and it allows screening with reduced automation coupled to a relatively simple low time-resolution plate reader. The liposome flux assay can also be used for the follow-up stages, estimation of compound potency and confirmation of liposome integrity. Second, we have identified a number of chemical groups that inhibit KtrAB activity in vitro with moderate potency, in the μM range.
Among the compounds identified, #36 (a urea derivative) was unique as it showed K+-dependent cell growth inhibition in the mutant strain JH642 (where KtrAB is essential for growth in low K+ concentrations) and in strain 168, which contains other K+ importers besides KtrAB. Functional characterization of compound #36 showed that KtrAB inhibition most likely results from binding to the membrane protein (KtrB). It also revealed in vitro inhibition of the NaK channel (a cation channel from B. cereus) with similar potency. Although KtrAB and NaK have very different amino acid sequences, they are both non-selective monovalent cation channels with ion pores that adopt the TM-Ploop-TM architecture. We speculate that compound #36 may bind to and recognize a structural feature present in both channels. This is not uncommon, as a number of studies have shown that molecules can bind to a membrane-buried crevice positioned between subunits or repeats of the ion pore of TM-Ploop-TM channels with very different sequences (Finol-Urdaneta et al., 2019; Zhu et al., 2020). Importantly, knockout of KtrB from strain 168 abolished the K+-dependent inhibition by compound #36, demonstrating specific targeting of KtrAB in the cell.
The K+-dependent phenotype associated with compound #36 may appear subtle. However, this is a striking result as it reveals that inhibition of a specific component of the bacterial K+ homeostasis machinery has an impact on bacterial growth, even when there is apparent functional redundancy due to the presence of other K+ transporters. Importantly, the K+-dependent phenotype observed with inhibited cells strongly resembles that of mutant strains with deficient K+ transport (Holtmann et al., 2003), which display perturbed growth under environmental challenges (for example, low K+ or osmotic stress). Bacteria meet this type of challenges at different stages of the process of infection (Fang et al., 2016; Janakiraman and Lesser, 2017), supporting our purpose of establishing a screening strategy for compounds that target components of the K+ homeostasis machinery. In addition, the high inhibitory concentrations required (100–500 μM) indicate that compound #36 displays modest potency. This feature may be resolved through chemical expansion, using a targeted medicinal chemistry approach. In fact, urea derivatives are present in a number of pharmacological drugs, are available commercially, and have shown antimicrobial activity (Patil et al., 2019a,b).
In conclusion, this proof-of-concept study opens the door for expanded screening campaigns, covering a larger chemical space, targeting conserved components of the K+ homeostasis machinery. With few exceptions, the current antibiotics have the same few cellular targets (cell wall, protein biosynthesis, and DNA biosynthesis). Being excluded from animal cells and essential to S. aureus and other bacteria with a marked role in virulence, antibacterial resistance and adaptation to environmental changes, SKT members are promising alternative targets for antibiotic development.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
JM-C, CH, and AF designed the experiments. AF, AP, and CT-D performed the experiments. AF and CH analyzed the data. JM-C, CH, AF, AP, and CT-D wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by FEDER funds through COMPETE 2020-POCI, Portugal 2020, and FCT – Fundação para a Ciência e a Tecnologia/Ministério da Ciência, Tecnologia e Ensino Superior: POCI-01-0145-FEDER-029863 (PTDC/BIA-BQM/29863/2017), and by “Fundação Luso-Americana para o Desenvolvimento” FLAD Life Science 2020 awarded to JM-C. We acknowledge FCT fellowship SFRH/BPD/105672/2015 and contract DL 57/2016/CP1355/CT0026 awarded to AF, fellowship SFRH/BPD/107785/2015 to AP, and fellowship SFRH/BD/123761/2016 to CT-D.
Conflict of Interest
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.603700/full#supplementary-material
Acknowledgments
We acknowledge the support of the i3S scientific platforms, BioSciences Screening (member of PPBI) – Portuguese Platform of Bioimaging (PPBI-POCI-01-0145-FEDER-022122); Biointerfaces and Nanotechnology; and Cell Culture and Genotyping.
References
Alam, A., and Jiang, Y. (2009). High-resolution structure of the open NaK channel. Nat. Struct. Mol. Biol. 16, 30–34. doi: 10.1038/nsmb.1531
Albright, R. A., Ibar, J. L., Kim, C. U., Gruner, S. M., and Morais-Cabral, J. H. (2006). The RCK domain of the KtrAB K+ transporter: multiple conformations of an octameric ring. Cell 126, 1147–1159. doi: 10.1016/j.cell.2006.08.028
Bai, Y., Yang, J., Zarrella, T. M., Zhang, Y., Metzger, D. W., and Bai, G. (2014). Cyclic di-AMP impairs potassium uptake mediated by a cyclic di-AMP binding protein in Streptococcus pneumoniae. J. Bacteriol. 196, 614–623. doi: 10.1128/JB.01041-13
Cao, Y., Pan, Y., Huang, H., Jin, X., Levin, E. J., Kloss, B., et al. (2013). Gating of the TrkH ion channel by its associated RCK protein TrkA. Nature 496, 317–322. doi: 10.1038/nature12056
Castaneda-Garcia, A., Do, T. T., and Blazquez, J. (2011). The K+ uptake regulator TrkA controls membrane potential, pH homeostasis and multidrug susceptibility in Mycobacterium smegmatis. J. Antimicrob. Chemother. 66, 1489–1498. doi: 10.1093/jac/dkr165
Commichau, F. M., Gibhardt, J., Halbedel, S., Gundlach, J., and Stulke, J. (2018). A delicate connection: c-di-AMP affects cell integrity by controlling osmolyte transport. Trends Microbiol. 26, 175–185. doi: 10.1016/j.tim.2017.09.003
Commichau, F. M., Heidemann, J. L., Ficner, R., and Stulke, J. (2019). Making and breaking of an essential poison: the Cyclases and Phosphodiesterases that produce and degrade the essential second messenger cyclic di-AMP in bacteria. J. Bacteriol. 201, e00462–e00418. doi: 10.1128/JB.00462-18
Corratge-Faillie, C., Jabnoune, M., Zimmermann, S., Very, A. A., Fizames, C., and Sentenac, H. (2010). Potassium and sodium transport in non-animal cells: the Trk/Ktr/HKT transporter family. Cell. Mol. Life Sci. 67, 2511–2532. doi: 10.1007/s00018-010-0317-7
Corrigan, R. M., Campeotto, I., Jeganathan, T., Roelofs, K. G., Lee, V. T., and Grundling, A. (2013). Systematic identification of conserved bacterial c-di-AMP receptor proteins. Proc. Natl. Acad. Sci. U. S. A. 110, 9084–9089. doi: 10.1073/pnas.1300595110
Diskowski, M., Mehdipour, A. R., Wunnicke, D., Mills, D. J., Mikusevic, V., Barland, N., et al. (2017). Helical jackknives control the gates of the double-pore K(+) uptake system KtrAB. eLife 6:e24303. doi: 10.7554/eLife.24303
Diskowski, M., Mikusevic, V., Stock, C., and Hanelt, I. (2015). Functional diversity of the superfamily of K(+) transporters to meet various requirements. Biol. Chem. 396, 1003–1014. doi: 10.1515/hsz-2015-0123
Durell, S. R., and Guy, H. R. (1999). Structural models of the KtrB, TrkH, and Trk1,2 symporters based on the structure of the KcsA K(+) channel. Biophys. J. 77, 789–807. doi: 10.1016/S0006-3495(99)76932-8
Durell, S. R., Hao, Y., Nakamura, T., Bakker, E. P., and Guy, H. R. (1999). Evolutionary relationship between K(+) channels and symporters. Biophys. J. 77, 775–788. doi: 10.1016/S0006-3495(99)76931-6
Epstein, W. (2003). The roles and regulation of potassium in bacteria. Prog. Nucleic Acid Res. Mol. Biol. 75, 293–320. doi: 10.1016/s0079-6603(03)75008-9
Fang, F. C., Frawley, E. R., Tapscott, T., and Vazquez-Torres, A. (2016). Bacterial stress responses during host infection. Cell Host Microbe 20, 133–143. doi: 10.1016/j.chom.2016.07.009
Finol-Urdaneta, R. K., McArthur, J. R., Goldschen-Ohm, M. P., Gaudet, R., Tikhonov, D. B., Zhorov, B. S., et al. (2019). Batrachotoxin acts as a stent to hold open homotetrameric prokaryotic voltage-gated sodium channels. J. Gen. Physiol. 151, 186–199. doi: 10.1085/jgp.201812278
Gao, B., Vorwerk, H., Huber, C., Lara-Tejero, M., Mohr, J., Goodman, A. L., et al. (2017). Metabolic and fitness determinants for in vitro growth and intestinal colonization of the bacterial pathogen Campylobacter jejuni. PLoS Biol. 15:e2001390. doi: 10.1371/journal.pbio.2001390
Gibson, D. G., Young, L., Chuang, R. Y., Venter, J. C., Hutchison, C. A. 3rd, and Smith, H. O. (2009). Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 6, 343–345. doi: 10.1038/nmeth.1318
Gries, C. M., Bose, J. L., Nuxoll, A. S., Fey, P. D., and Bayles, K. W. (2013). The Ktr potassium transport system in Staphylococcus aureus and its role in cell physiology, antimicrobial resistance and pathogenesis. Mol. Microbiol. 89, 760–773. doi: 10.1111/mmi.12312
Gundlach, J., Herzberg, C., Kaever, V., Gunka, K., Hoffmann, T., Weiss, M., et al. (2017). Control of potassium homeostasis is an essential function of the second messenger cyclic di-AMP in Bacillus subtilis. Sci. Signal. 10:eaal3011. doi: 10.1126/scisignal.aal3011
Holtmann, G., Bakker, E. P., Uozumi, N., and Bremer, E. (2003). KtrAB and KtrCD: two K+ uptake systems in Bacillus subtilis and their role in adaptation to hypertonicity. J. Bacteriol. 185, 1289–1298. doi: 10.1128/JB.185.4.1289-1298.2003
Janakiraman, A., and Lesser, C. F. (2017). How to manage stress: Lessons from an intracellular pathogen. Virulence 8, 359–361. doi: 10.1080/21505594.2016.1256538
Kim, H., Youn, S. J., Kim, S. O., Ko, J., Lee, J. O., and Choi, B. S. (2015). Structural studies of potassium transport protein KtrA regulator of conductance of K+ (RCK) C domain in complex with cyclic diadenosine monophosphate (c-di-AMP). J. Biol. Chem. 290, 16393–16402. doi: 10.1074/jbc.M115.641340
Lazar, V., Pal Singh, G., Spohn, R., Nagy, I., Horvath, B., Hrtyan, M., et al. (2013). Bacterial evolution of antibiotic hypersensitivity. Mol. Syst. Biol. 9:700. doi: 10.1038/msb.2013.57
Levin, E. J., and Zhou, M. (2014). Recent progress on the structure and function of the TrkH/KtrB ion channel. Curr. Opin. Struct. Biol. 27, 95–101. doi: 10.1016/j.sbi.2014.06.004
Liu, Y., Ho, K. K., Su, J., Gong, H., Chang, A. C., and Lu, S. (2013). Potassium transport of Salmonella is important for type III secretion and pathogenesis. Microbiology 159, 1705–1719. doi: 10.1099/mic.0.068700-0
MacGilvary, N. J., Kevorkian, Y. L., and Tan, S. (2019). Potassium response and homeostasis in Mycobacterium tuberculosis modulates environmental adaptation and is important for host colonization. PLoS Pathog. 15:e1007591. doi: 10.1371/journal.ppat.1007591
Nelson, J. W., Sudarsan, N., Furukawa, K., Weinberg, Z., Wang, J. X., and Breaker, R. R. (2013). Riboswitches in eubacteria sense the second messenger c-di-AMP. Nat. Chem. Biol. 9, 834–839. doi: 10.1038/nchembio.1363
Oz, T., Guvenek, A., Yildiz, S., Karaboga, E., Tamer, Y. T., Mumcuyan, N., et al. (2014). Strength of selection pressure is an important parameter contributing to the complexity of antibiotic resistance evolution. Mol. Biol. Evol. 31, 2387–2401. doi: 10.1093/molbev/msu191
Pan, F., Zhang, H., Dong, X., Ye, W., He, P., Zhang, S., et al. (2018). Comparative genomic analysis of multidrug-resistant Streptococcus pneumoniae isolates. Infect. Drug Resist. 11, 659–670. doi: 10.2147/IDR.S147858
Pang, M., Xie, X., Dong, Y., Du, H., Wang, N., Lu, C., et al. (2017). Identification of novel virulence-related genes in Aeromonas hydrophila by screening transposon mutants in a Tetrahymena infection model. Vet. Microbiol. 199, 36–46. doi: 10.1016/j.vetmic.2016.12.021
Patil, M., Noonikara-Poyil, A., Joshi, S. D., Patil, S. A., Patil, S. A., and Bugarin, A. (2019a). New Urea Derivatives as Potential Antimicrobial Agents: Synthesis, Biological Evaluation, and Molecular Docking Studies. Antibiotics 8:178. doi: 10.3390/antibiotics8040178
Patil, M., Poyil, A. N., Joshi, S. D., Patil, S. A., Patil, S. A., and Bugarin, A. (2019b). Synthesis, molecular docking studies, and antimicrobial evaluation of new structurally diverse ureas. Bioorg. Chem. 87, 302–311. doi: 10.1016/j.bioorg.2019.03.024
Richey, B., Cayley, D. S., Mossing, M. C., Kolka, C., Anderson, C. F., Farrar, T. C., et al. (1987). Variability of the intracellular ionic environment of Escherichia coli. Differences between in vitro and in vivo effects of ion concentrations on protein-DNA interactions and gene expression. J. Biol. Chem. 262, 7157–7164. doi: 10.1016/S0021-9258(18)48218-X
Rocha, R., Teixeira-Duarte, C. M., Jorge, J. M. P., and Morais-Cabral, J. H. (2019). Characterization of the molecular properties of KtrC, a second RCK domain that regulates a Ktr channel in Bacillus subtilis. J. Struct. Biol. 205, 34–43. doi: 10.1016/j.jsb.2019.02.002
Schrecker, M., Wunnicke, D., and Hanelt, I. (2019). How RCK domains regulate gating of K+ channels. Biol. Chem. 400, 1303–1322. doi: 10.1515/hsz-2019-0153
Srivatsan, A., Han, Y., Peng, J., Tehranchi, A. K., Gibbs, R., Wang, J. D., et al. (2008). High-precision, whole-genome sequencing of laboratory strains facilitates genetic studies. PLoS Genet. 4:e1000139. doi: 10.1371/journal.pgen.1000139
Su, Z., Brown, E. C., Wang, W., and MacKinnon, R. (2016). Novel cell-free high-throughput screening method for pharmacological tools targeting K+ channels. Proc. Natl. Acad. Sci. U. S. A. 113, 5748–5753. doi: 10.1073/pnas.1602815113
Szollosi, A., Vieira-Pires, R. S., Teixeira-Duarte, C. M., Rocha, R., and Morais-Cabral, J. H. (2016). Dissecting the Molecular Mechanism of Nucleotide-Dependent Activation of the KtrAB K+ Transporter. PLoS Biol. 14:e1002356. doi: 10.1371/journal.pbio.1002356
Teixeira-Duarte, C. M., Fonseca, F., and Morais-Cabral, J. H. (2019). Activation of a nucleotide-dependent RCK domain requires binding of a cation cofactor to a conserved site. elife 8:e50661. doi: 10.7554/eLife.50661
Valente, R. S., and Xavier, K. B. (2016). The Trk potassium transporter is required for RsmB-mediated activation of virulence in the Phytopathogen Pectobacterium wasabiae. J. Bacteriol. 198, 248–255. doi: 10.1128/JB.00569-15
Vieira-Pires, R. S., Szollosi, A., and Morais-Cabral, J. H. (2013). The structure of the KtrAB potassium transporter. Nature 496, 323–328. doi: 10.1038/nature12055
Whatmore, A. M., Chudek, J. A., and Reed, R. H. (1990). The effects of osmotic upshock on the intracellular solute pools of Bacillus subtilis. J. Gen. Microbiol. 136, 2527–2535. doi: 10.1099/00221287-136-12-2527
Wistrand-Yuen, E., Knopp, M., Hjort, K., Koskiniemi, S., Berg, O. G., and Andersson, D. I. (2018). Evolution of high-level resistance during low-level antibiotic exposure. Nat. Commun. 9:1599. doi: 10.1038/s41467-018-04059-1
Yasbin, R. E., Wilson, G. A., and Young, F. E. (1975). Transformation and transfection in lysogenic strains of Bacillus subtilis: evidence for selective induction of prophage in competent cells. J. Bacteriol. 121, 296–304. doi: 10.1128/JB.121.1.296-304.1975
Zhang, J. H., Chung, T. D., and Oldenburg, K. R. (1999). A simple statistical parameter for use in evaluation and validation of high throughput screening assays. J. Biomol. Screen. 4, 67–73. doi: 10.1177/108705719900400206
Keywords: KtrAB, high-throughput screen, antibacterial target, K+ homeostasis, Bacillus subtilis, superfamily of potassium transporters, membrane protein
Citation: Fernandes AS, Pombinho A, Teixeira-Duarte CM, Morais-Cabral JH and Harley CA (2021) Fluorometric Liposome Screen for Inhibitors of a Physiologically Important Bacterial Ion Channel. Front. Microbiol. 12:603700. doi: 10.3389/fmicb.2021.603700
Edited by:
Natalia V. Kirienko, Rice University, United StatesReviewed by:
Peter Graumann, University of Marburg, GermanyPiotr Koprowski, Nencki Institute of Experimental Biology (PAS), Poland
Copyright © 2021 Fernandes, Pombinho, Teixeira-Duarte, Morais-Cabral and Harley. 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(s) 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.
*Correspondence: Carol A. Harley, aGFybGV5Y2FAaWJtYy51cC5wdA==