Voltage-gated ion channels are responsible for the electrical excitability of neurons and cardiomyocytes. Thus, they are obvious targets for pharmaceuticals aimed to modulate excitability. Compounds activating voltage-gated potassium (KV) channels are expected to reduce excitability. To search for new KV-channel activators, we performed a high-throughput screen of 10,000 compounds on a specially designed Shaker KV channel. Here, we report on a large family of channel-activating compounds with a carboxyl (COOH) group as the common motif. The most potent COOH activators are lipophilic (4 < LogP <7) and are suggested to bind at the interface between the lipid bilayer and the channel’s positively charged voltage sensor. The negatively charged form of the COOH-group compounds is suggested to open the channel by electrostatically pulling the voltage sensor to an activated state. Several of the COOH-group compounds also activate the therapeutically important KV7.2/7.3 channel and can thus potentially be developed into antiseizure drugs. The COOH-group compounds identified in this study are suggested to act via the same site and mechanism of action as previously studied COOH-group compounds, such as polyunsaturated fatty acids and resin acids, but distinct from sites for several other types of potassium channel–activating compounds.

Voltage-gated potassium (KV) channels play a crucial role in cardiac and neuronal excitability (Hille, 2001). Inherited or acquired mutations in KV channel genes can cause cardiovascular, neurological, and metabolic disorders known as “channelopathies” (Ashcroft, 2000), and pharmaceuticals, such as the antiseizure drug retigabine, target a KV channel (Miceli et al., 2018).

KV channels consist of two types of domains: a pore domain surrounded by four voltage-sensor domains (VSD) (Fig. 1 A, for clarity only two VSDs are shown). The VSD contains a mobile voltage sensor S4, which responds to alterations in the voltage across the membrane; when the intracellular side becomes more positive, the positively charged voltage sensor moves from a down-state to an up-state to pull the gate open and lets potassium ions flow out from the cell (Fig. 1 A) (Delemotte et al., 2011; Henrion et al., 2012; Jensen et al., 2012; Long et al., 2007). While pharmaceuticals targeting the pore domain to block the ion-conducting pore have been known and used for decades (Hille, 1977; Ragsdale et al., 1994; Strichartz, 1976; Zhou et al., 2001), few drugs activate voltage-gated channels, and no drug targeting the VSD is available in clinical practice.

Previous research has shown that lipophilic and negatively charged small-molecule compounds, such as some resin acids and derivatives thereof, and polyunsaturated fatty acids (PUFAs) are effective activators of several KV channels and have the potential to suppress neuronal and cardiac excitability (Börjesson et al., 2008; Cui et al., 2016; Imaizumi et al., 2002; Liin et al., 2016; Ottosson et al., 2015, 2021; Salari et al., 2018). One proposed mechanism is that the negatively charged compound binds at the interface between the membrane’s lipid bilayer and the channel’s VSD to electrostatically interact with the positively charged arginines of S4 in the VSD, which drives the voltage sensor S4 from a resting (down) state to an activated (up) state, to open the ion-conducting pore (Site 1 in Fig. 1 A) (Börjesson and Elinder, 2011; Ottosson et al., 2017). We refer to these types of compounds as “lipoelectric,” because the membrane-exposed site requires lipophilicity and the electrostatic pulling requires electrostatic interactions between charges (Börjesson et al., 2008).

In contrast, other types of compounds bind to other parts of the VSD to open the channel via other mechanisms. For instance, biaryl sulfonamides (Fig. 1 B) bind from the extracellular solution to the top of the VSD (Site 2 in Fig. 1 A); they bind in a crevice between the transmembrane segments of the VSD to keep the voltage sensor S4 in an activated (up) state in both Na and K channels (Ahuja et al., 2015; Liin et al., 2018). Recently, a large family of warfarin-like tautomers was reported to bind to the intracellular side of the human KV1.5 and the Drosophila Shaker KV channel, to the linker between the VSD and pore domain (Site 3 in Fig. 1 A), and to stabilize this linker in an activated state (Silverå Ejneby et al., 2020). In addition, compounds can bind to the pore domain to activate the channel (4 in Fig. 1 A). A well-known example is the discontinued antiseizure drug retigabine (Wuttke et al, 2005).

A crucial point in the development of pharmaceuticals against specific diseases is to find drugs that selectively act on one type of channel while leaving other channels unaffected. In this study, we searched for new lipoelectric compounds. To do so, we used the Shaker KV channel supplemented with two additional gating charges at the top of the voltage sensor S4 to increase the electrostatic effect (Ottosson et al., 2014). This channel will be referred to as the 3R Shaker KV channel because three arginines (R; the native top arginine in addition to the two added ones) are important for the effect. The reason for using this supplemented channel was that we expected it to identify even weak lipoelectric channel activators, which might be activators of other KV channels too.

Here, we report on a high-throughput screen of 10,096 compounds with follow-up experiments on the site and mechanism of action. We describe a large family of carboxyl (COOH)-group compounds (Fig. 1 B) acting via the lipoelectric mechanism, distinct from several other KV-channel activating compounds.

High-throughput screen

The details of the high-throughput screen have been described elsewhere (Liin et al., 2018; Silverå Ejneby et al., 2020; and references therein). In brief, the 3R Shaker KV channel (Ottosson et al., 2014) was expressed in a CHO-K1 stable cell line and studied with a 384-well IonWorks equipment. The holding voltage was ‒80 mV, and 150 ms voltage-clamp steps were applied to +10 mV, +30 mV, and +50 mV. The cells were kept at ‒80 mV for 300 ms between each step. The current signal was sampled at 10 kHz. No external electronic filter was used. All recordings were made at room temperature (∼21°C). Modulation of the KV-mediated current was assessed by dividing the post-scan steady-state KV current (for the last 10 ms of each pulse) by the respective prescan KV current for each well and corrected for run down. For the quantitative analysis presented in this paper, we used the data from +10 mV. All data between 0.75 and 1.25 (mean ± 3 SD) were considered to have no effect. If multiple experiments were carried out, we adjusted the significance level by dividing the deviation from 1 by n0.5, where n is the number of experiments.

Compounds were searched for in the following chemical libraries (library name [company or organization]): Acids (Vitas-M Laboratory Ltd.); Consortium set (Specs)*; Drug Like (Enamine)*; Hit Finder (Maybridge)*; Ion Channel Ligand and Nuclear Receptor Ligand (Enzo Life Sciences Inc.); Known Drugs (Prestwick Chemical); LCBKI Primary Screening Set (Chemical Biology Consortium Sweden)*; Natural derivatives (TimTec); NIH Clinical Collection (National Institutes of Health); Peptidomimetics (ChemDiv Inc.)*; Tocriscreen Plus (Tocris Bioscience). Apart from LCBKI, provided by the Chemical Biology Consortium Sweden, and the NIH clinical collection, provided by the Libraries Roadmap Initiative, all compound libraries were commercially available. Whenever a selection was made (denoted by *), a chemist was instructed to search for small, hydrophobic, and partially charged compounds. The final concentration, in the test solution, was 10 µM of the compound. The tested compounds were added to the extracellular solution.

Analysis of high-throughput data—frequency diagrams

Frequency diagrams (shown in Fig. 3) were constructed by calculating the floating average frequencies of a specific family of compounds among all compounds. All data on the effects of the compounds were binned in intervals of 0.01. The number of compounds per bin is plotted in Fig. 2 A. An effect value of 1 means no effect and the grey area (0.75–1.25) means no significant effect. To calculate the frequency of a specific family of compounds among all compounds for a specific effect, we included all data points (one per compound) within a range of ± 0.05. In the effect range 0.54–1.47, the number of included data for each calculated frequency point was 50–2,394 (Fig. 2 B, red line). Outside this range, the effect window was increased to make sure that at least 50 data/compounds were used for each plotted point (Fig. 2 B black curves). The effect value (the relative current plotted on the x-axis in Fig. 3) was the average effect of all included data/compounds. The original screen included 10,096 compounds (Silverå Ejneby et al., 2020). In this previous investigation, we reported on a small family (259 compounds) of a distinct warfarin-like tautomer motif, which turned out to be a very effective KV channel activator. In the presented frequency analysis, we excluded this family, which left us with 9,837 compounds.

KV channels expressed in oocytes

The Shaker H4 channel (Kamb et al., 1988) was modified to remove the fast N-type inactivation (ShH4IR) (Hoshi et al., 1990). This channel will be referred to as the wild-type (WT) Shaker channel. To enhance the sensitivity to lipophilic and negatively charged compounds, we introduced two positively charged arginines (M356R and A359R) in the extracellular end of the voltage sensor S4; these two arginines in addition to one native arginine (R362) made the channel more sensitive to the PUFA docosahexaenoic acid and resin acids (Ottosson et al., 2014). Because all three arginines are important for the effect, we refer to this channel as the 3R Shaker K channel. The human KV7.2/7.3 channel was expressed through coinjection of KV7.2 (Gene bank accession number NM_004518) and hKV7.3 (Gene bank accession number NM_004519) in a 1:1 ratio (Liin et al., 2016). Point mutations in KV7.2 were introduced using site-directed mutagenesis (QuikChange II, Agilent) and verified by sequencing at a core facility at Linköping University, Sweden. cRNA (from linearized DNA) was prepared using the mMessage mMachine T7 kit (Invitrogen) (Larsson et al., 2020).

Preparation and injection of oocytes

All animal experiments were approved by Linköping’s local Animal Care and Use Committee. The experiments were performed in accordance with relevant guidelines and regulations. African clawed frogs (Xenopus laevis) were anesthetized with 1.4 g/liter ethyl 3-aminobenzoate methanesulfonate salt (tricaine). After an incision through the abdomen, a batch of oocytes was removed. Clusters of oocytes were separated by incubation for ∼1 h in a Ca-free O-R2 solution (in mM: 82.5 NaCl, 2 KCl, 5 HEPES, and 1 MgCl2; pH adjusted to 7.4 by NaOH) containing Liberase Blendzyme. The oocytes were then incubated at 8°C in a modified Barth’s solution (MBS; in mM: 88 NaCl, 1 KCl, 2.4 NaHCO3, 15 HEPES, 0.33 Ca(NO3)2, 0.41 CaCl2 and 0.82 MgSO4; pH adjusted to 7.6 by NaOH) supplemented with penicillin (25 U/ml), streptomycin (25 μg/ml), and sodium pyruvate (2.5 mM) for 2–24 h before injection. 50 nl of cRNA (50 pg) were injected into each oocyte using a Nanoject injector (Drummond Scientific). Injected oocytes were kept at 8°C in MBS until 1 day before electrophysiological recordings, after which they were incubated at 16°C. All chemicals were supplied by Sigma-Aldrich unless stated otherwise.

Manual electrophysiology

All manual two-electrode voltage clamp recordings were performed at room temperature (20–23°C) using a GeneClamp 500B amplifier and a Digidata 1440A digitizer with pClamp 10 software (all from Molecular Devices, Inc.). Data were sampled at 10 µs and the signal was filtered at 2 kHz by a Bessel filter. The amplifiers’ leakage subtraction was not used, but digital leakage subtraction was used during the analysis to remove any leakage current between the two most negative test pulses.

Compounds were initially dissolved at 30 mM in 99.5% EtOH and stored at −20°C. Compounds were subsequently diluted to the desired test concentration in an extracellular solution. The oocyte was placed in a bath surrounded by 1 K extracellular solution that contained (in mM): 88 NaCl, 1 KCl, 0.4 CaCl2, 0.8 MgCl2, and 15 MES (for pH 5.5 and 6.5), HEPES (for pH 7.4 and 8.0), or CHES (for pH 9.0), pH adjusted by NaOH (reaching a sodium concentration of ∼100 mM). The control solution and the compound solution were added to the bath with a peristaltic pump-driven perfusion system. Two glass microelectrodes were inserted into the oocyte using micromanipulators. The microelectrodes were pulled from borosilicate glass, filled with 3 M KCl, and had a resistance of 0.5–2 MΩ.

For the Shaker potassium channel recordings, all channels were closed when the membrane voltage was clamped to −80 mV, and this voltage was set as the holding voltage. Currents were evoked from the holding voltage of −80 mV by 100-ms long, 5-mV steps ranging from −80 up to +50 mV (WT) or +70 mV (3R). For the WT human KV7.2/7.3 channel, the holding voltage was set to −100 mV and currents were evoked by 3,000-ms long 10 mV steps ranging from −120 up to +50 mV. This was followed by a 2,000-ms-long step at −30 mV, producing a tail current from which the initial currents were measured and analyzed. For the R198Q KV7.2 mutant channel, the holding voltage was set to −120 mV and currents were evoked by 6,000-ms long 10 mV steps ranging from −140 up to +10 mV, which was followed by a tail current-producing 3,000-ms long step at −60 mV from which the initial currents were measured and analyzed.

To tune the charge of some compounds, when studying the effect on the WT KV7.2/7.3 channel we altered the pH in the extracellular solution. The effect of the compound at a certain pH was always compared with a control solution at the same pH and the compound solution. The pH alteration per se had a modest effect on the apparent voltage sensed by the channel, seen as a shift of voltage-dependent parameters along the voltage axis (Prole et al., 2003).

Analysis of manual electrophysiology data

The manual electrophysiological data from oocytes were processed and analyzed with Clampfit 11.2 (Molecular Devices, LLC.) and GraphPad Prism 10 (GraphPad Software, Inc.).

The conductance GK(V) was calculated as
(1)
where IK is the average steady-state current at the end of a 100-ms pulse, V is the absolute membrane voltage, and Vrev is the reversal voltage for K+ (set to −80 mV).
These data were fitted to a Boltzmann equation,
(2)
where A is the amplitude of the curve, V1/2 is the midpoint, s is the slope, and n is an exponent set to 4 for Shaker KV channels and to 1 for the human KV7.2/7.3 and the mutant R198Q KV7.2 channels. The compound-induced G(V) shift was determined at the 10% level of the maximum conductance in the control solution for Shaker KV channels, and the V1/2 was calculated for the human KV7.2/7.3 channel with a shared slope for the control and compound curves.
Data for concentration and pH-dependent G(V) shifts and ∆VG(V) in oocytes were fit with a dose–response equation,
(3)
where ∆Vmax is the maximal shift, c½ is half maximal effective concentration or the pKa value, and c is the concentration. The main purpose of using Eq. 3 was to quantify the effect to compare the effects of different compounds. Because we did not use concentrations above 100 µM (because of solubility problems of the most lipophilic compounds), a correct determination of both c½ and ∆Vmax is not possible. Therefore, we constrained the ∆Vmax to be equal for all compounds and used the obtained c½ as a measure of the potency of the compounds.

Calculated chemical properties

JChem for Excel version 23.5.0.480 (ChemAxon Ltd.) was used for drawing chemical structures and calculating the molecular weight, the acidic dissociation constant (pKa), and the octanol/water partitioning coefficients (LogP and LogD) for uncharged and charged compounds, respectively, as previously described (Liin et al., 2018; Silverå Ejneby et al., 2020).

Statistical analysis

Average values are expressed as mean ± SEM. Statistical analyses were done using an unpaired t test when comparing two data points against each other. P < 0.05 was considered statistically significant. Most statistical analyses were performed in GraphPad Prism 10 (GraphPad Software, Inc.). The statistical analyses of the frequency diagrams (Fig. 3) and associated Table 1 were done with skewed binomial distributions calculated in Microsoft Excel 16 (Microsoft Corporation).

A high-throughput screen identified KV channel activators

To search for KV channel activators, we screened 10,096 compounds on the 3R Shaker KV channel. To increase the possibility of finding “lipoelectric” compounds, the channel was supplemented with two extra positively charged residues in the extracellular portion of the voltage-sensor S4, referred to as the 3R Shaker KV channel (Ottosson et al., 2014). The methodology of the screen has been described before (Liin et al., 2018; Silverå Ejneby et al., 2020) and we have previously found and described two families of compounds acting as KV-channel activators: (1) one family has a distinct warfarin-like motif (compounds bind to Site 3 in Fig. 1 A [Silverå Ejneby et al., 2020]). In total, 259 compounds belonged to this family, and because data from this family was analyzed in its entirety, we excluded these compounds from the present analysis, leaving us with a total of 9,837 compounds. (2) The other family had a biaryl sulfonamide motif (compounds bind to Site 2 in Fig. 1 A [Liin et al., 2018]). Because we only analyzed data from a small fraction of these compounds (26 compounds in total), compounds from this family will be included in the initial analysis of this study.

The 9,837 analyzed compounds were divided into six classes (Fig. 1 B and Table 1): compounds containing a (1) COOH group (2,298 compounds in total), (2) amide group (2,021 compounds; excluding compounds in [3]), (3) biaryl amide groups (BAA; 733 compounds; an amide group flanked by, in most cases, aromatic rings, but in some cases also other rings), (4) sulfonamide groups (270 compounds, excluding compounds in [5]), (5) biaryl sulfonamides groups (BASA; 605 compounds, a sulfonamide group flanked by, in most cases, aromatic ring, but in some cases also other rings), and (6) none of (1)–(5) (other; 5,012 compounds). Because there is some overlap between the groups, the sum is larger than the total number of compounds. All combinations and numbers are presented in Table 1. To see if any of the groups constitute more activators or inhibitors than expected from a random distribution, we constructed frequency diagrams (Fig. 3). All compounds with an effect E ± 0.05 (x-axis) were included in the frequency analysis of effect E. For effects E < 0.54 and E > 1.47, the range was increased to make sure that data from 50 compounds were included in each data point (see Materials and methods).

The COOH compounds were clearly overrepresented among the clear activators (E ≥ 1.60, Fig. 3 A and Table 1); while constituting only 23% of all compounds, they constitute almost 60% of the best activators, suggesting that the COOH group might be a critical part of a KV-channel activator. Amides are neither activators nor inhibitors, while BAAs are clearly overrepresented among the activators but not among the inhibitors (Fig. 3 A and Table 1). Sulfonamides are neither activators nor inhibitors, while BASAs are overrepresented among weak activators and weak inhibitors (see analysis below; Fig. 3 A and Table 1). The last group, containing all other compounds, constitutes 51% of all compounds. They are clearly overrepresented among the best inhibitors but underrepresented among the activators (Fig. 3 A and Table 1).

However, while COOHs, BAAs, and BASAs are overrepresented among the activators, we cannot rule out that it is the combination of these motifs that makes them activators. Therefore, we also analyzed the subgroups. COOHs in isolation, in combination with BAAs, and in total were significantly increased among the activators (Fig. 3 B and Table 1). In contrast, BAAs were only increased in combination with COOH and in total, but not in isolation (Fig. 3 C and Table 1). This suggests that it is the COOH group per se that is critical and possibly a COOH group in combination with BAA that makes a fruitful mix. The BASA motif in any combination is neither a significant activator nor an inhibitor if similarly analyzed (Table 1), but if we perform a statistical analysis of BASA occurrence in more narrow windows, we find that BASAs in isolation is a significant weak activator and weak inhibitor (Fig. 3 D); data points outside the vertical dashed lines above the continuous horizontal line denote statistically significant overrepresentation of the motif. In a previous study, we showed that BASAs bind in a groove at the top of the VSD to keep the voltage sensor S4 in an up-state (Liin et al., 2018). The current analysis suggests that BASAs also can be weak inhibitors, possibly by binding to the same site to press down S4. This will not be analyzed further in this article.

Carboxyl compounds are clear KV channel activators

Above we described that a COOH group is the most significative functional group of KV channel activators, possibly in combination with BAA. While most of the 2,298 COOH compounds did not have any effect or were weak inhibitors, 157 of them (6.8% of all COOH compounds) significantly increased the current (see Materials and methods for significance calculation), but most of them were only weak activators. 10 µM of 38 different COOH compounds increased the current of the 3R Shaker KV channel at +10 mV by at least 60% (Effect E ≥ 1.60; Table 2). These strong activators were slightly heavier than the average COOH compound (384 g/mol (red) versus 348 g/mol (black); Fig. 4 A and Table 2), they were clearly more hydrophobic (LogP = 4.51 versus 2.36, Fig. 4 B and Table 2; LogD = 1.82 versus −0.30, Fig. 4 C and Table 2) and had a marginally higher pKa value (3.89 versus 3.54, Fig. 4 D and Table 2).

Of the 38 strong activators 24 were run in duplicates in the screen (Table 2), and the mean deviation from the average value for these 24 compounds was 20%, suggesting robust data from the screen. Even though these 38 compounds do not unexpectedly show great variability in structure, some common motifs can be found (Fig. 5). In most of the compounds (29/38), the COOH group is directly attached to an aromatic ring of six or five atoms (columns 1–4, Fig. 5). 15 of these have a characteristic motif with an amide group (in most cases in ortho position) directly connected to the aromatic ring (columns 2–3, Fig. 4). Most of these also have an aromatic or non-aromatic ring directly connected to the amide motif (the BAA compounds). In other cases, there is an adjacent ring to the COOH connected ring (column 4, Fig. 4). In the last group, the COOH group is connected to a linker or a free carbon chain (column 5, Fig. 4). These structures are not surprising; in previous studies, we have reported that other COOH compounds of the types mentioned here can activate the 3R Shaker KV channel: polyunsaturated fatty acids (Ottosson et al., 2014), resin acids and a derivative thereof (Ottosson et al., 2015, 2017), and resin acid derivatives where the COOH group is attached via a carbon chain (a stalk), optimally of two atoms (Ottosson et al., 2021; Silvera Ejneby et al., 2018). While the shapes of the activators are diverse, the COOH/BAA motif is striking. 12 of the 38 top compounds have this motif (in columns 2–3, Fig. 3). In the remainder of this article, we will focus on a selection of compounds from this top group and study the effects in more detail to decide on the site of action and possible selectivity.

COOH-group compounds electrostatically act on the voltage sensor to open Shaker KV channels

To explore the site and mechanism of action of the COOH compounds and whether there is a common site and mechanism of action, we explored the effect of nine of the compounds with two-electrode voltage clamp on 3R Shaker KV channels expressed in Xenopus oocytes. Compound #29 at 100 µM increased the current at 0 mV by a factor of 58 (Fig. 6, A and B). This current increase was mainly caused by a shift of the conductance versus voltage, G(V), a curve in the negative direction along the voltage axis, and to a much smaller degree by an increase of the maximum conductance (Fig. 6 B). The shift of the G(V) curve was −45.9 ± 2.1 mV (n = 3). Also, 10 µM of compound #29 had a large effect on the current (Fig. 6 C) and the shift was −14.4 ± 1.5 mV (n = 3). All nine explored compounds’ concentration dependently shifted the G(V) along the voltage axis (Fig. 6 D). To quantify the data in a simple way, we assumed that all compounds shifted the G(V) equally much when fully bound (see Materials and methods). The common maximum shift was −47.0 mV (dotted line in Fig. 6 D). Eight of the compounds had estimated c½ values in the range of 16–51 µM, while one (#54) had a significantly higher c½ value of 141 µM.

Because the COOH compounds were identified in the screen with the 3R Shaker KV channel, our hypothesis was that these compounds act from the lipid bilayer/S4 interface via an electrostatic mechanism to pull the voltage sensor to an activated position. To test this hypothesis, we explored the effect of the compounds on the WT Shaker KV channel, which lacks the two supplemented arginines at the top of the voltage sensor S4 of the 3R Shaker KV channel. In previous research on PUFAs and resin-acid derivatives, suggested to electrostatically act on S4, we have found a reduction by a factor of about 2–3 for WT compared with 3R in the concentration range 10–100 µM (Ottosson et al., 2014, 2015, 2017; Silverå Ejneby et al., 2021). 10 µM of compound #29 significantly opened the WT Shaker KV channel (Fig. 7 A) by shifting the G(V) curve in a negative direction along the voltage axis (Fig. 7 B; −7.7 ± 0.5 mV, n = 3), but the G(V) shift, as well as the increase in maximum conductance, is reduced (Fig. 7 B, compare with Fig. 6 B). For most compounds, the G(V) shift at 10 µM was reduced by a factor of 2–3 (dashed lines in Fig. 7 C). (Note that only eight of the original compounds were studied on the WT channel because the last compound was not available to order for these experiments.) Thus, these data suggest that all studied COOH-group compounds act via an electrostatic mechanism. This contrasts sharply with other chemical motifs, previously published from the screen: the warfarin-like tautomers (Silverå Ejneby et al., 2020) and the biaryl sulfonamides (Liin et al., 2018).

COOH-group compounds open the human KV7.2/7.3 channel

To investigate if the COOH compounds also act on a pharmaceutically more relevant channel, we explored the effects on the human KV7.2/7.3 channel, which is a validated target for antiseizure drugs (Main et al., 2000; Schenzer et al., 2005). 100 µM of compound #51 clearly increased the current at −60 mV by a factor of about nine (red traces in Fig. 8, A and B). This increase is mainly caused by a shift of the G(V) curve in a negative direction along the voltage axis (Fig. 8 C). In addition, there is a small increase in the maximum conductance (Fig. 8 C). To estimate apparent affinities to the channel, we employed the same strategy as above for the 3R Shaker KV channel. The common maximum shift was −17.2 mV (Fig. 8 D). Two of the compounds had c½ values in the range of 30–50 µM (#29 and #51), while the rest had c½ values >100 µM.

How does the effect on the KV7.2/7.3 channel compare with the effect on the WT Shaker KV channel? Four of the compounds (#13, #42, #51, and #54) had almost the same effect at 10 µM on the KV7.2/7.3 channel and the WT Shaker KV channel (Fig. 9, A, B, and G), one compound (#29) had half the effect on KV7.2/7.3 compared with the WT Shaker KV channel (Fig. 9, C and G), while two compounds (#7 and #27) had almost no effect on the KV7.2/7.3 channel despite having significant effects on the WT Shaker KV channel (Fig. 9, E, B, and G). In particular, compound #7 is interesting; it was the most potent compound in the screen of 9,837 compounds (excluding the warfarin-like tautomers) and it was the most potent COOH compound (at 10 µM) in the studies on both the WT Shaker KV channel and the 3R Shaker channel expressed in oocytes (Fig. 3 F). However, the compound had no effect on the KV7.2/7.3 channel (at 10 µM, Fig. 8 D). This suggests that it should be possible to reach selectivity between these two channels by modulating the structure of the compounds. However, based on the available data in this study, we think it is premature to speculate about which molecular elements are critical to reach selectivity.

Electrostatics is the major player in compound-induced activation of the KV7.2/7.3 channel

The data on the Shaker KV channel suggest that compound-induced channel activation is caused by an electrostatic effect. To test if this also is the reason for the effect on the WT KV7.2/7.3 channel, we performed two types of experiments: (1) an alteration of the charge of the compound and (2) an alteration of the charge of the channel.

First, we altered the charge of the compound by altering the pH in the extracellular solution. Although the COOH compounds have theoretical pKa values around 4 (Table 1) and should be almost fully charged in a water solution at neutral pH, the proposed site of action in the low dielectric environment in the lipid bilayer pushes the functional pKa value up to around 7 (Börjesson et al., 2008; Rooney et al., 1983). The G(V) shift induced by 100 µM of compound #51 changed from −19.1 ± 2.6 mV (n = 4) at pH 8.0 to +5.5 ± 1.1 mV (n = 4) at pH 5.5; the apparent pKa value was 6.8 ± 0.2, the maximum (asymptotic) negative shift, for a fully charged compound at high pH values was −20.2 ± 1.7 mV, and the minimum (asymptotic) shift and for a non-charged compound at low pH values it was +6.7 ± 2.6 mV (Fig. 10 A). This pH dependence supports the hypothesis that the effect is electrostatic. The small shift in positive direction along the voltage axis at low pH suggests that the neutral form of the compound has a slightly higher affinity for the closed, non-activated, channel than for the open, activated, channel. We failed to record reliable currents at pH 9 with 100 µM of compound #51 (the oocytes did not tolerate this treatment), but for 10 µM of compound #51, we succeeded in obtaining data for the complete pH range from 5.5 to 9. The pKa value was not changed (7.0 ± 0.7 mV) but the (total) amplitude was reduced by 50% (from 26.9 to 13.7 mV; Fig. 10 A). To control for possible compound precipitation at low pH values, we prepared a solution of 100 µM of compound #51 at pH 5.5, which was centrifuged at 16 g for 1 h. The supernatant was extracted and brought to pH 8.0 using 10 M NaOH, after which experiments were conducted according to the standard protocol. The resulting G(V) shifts were −16.6 ± 0.9 mV (n = 5), thus similar to the results from the standard experiments at pH 8.0, suggesting that compound precipitation is not the reason for the reduction of the G(V) shift at low pH values.

Second, we altered the charge of the voltage sensor by neutralizing the top charge of S4 in KV7.2 (R198Q) and expressed it as a homotetramer (Ottosson et al., 2021). The induced G(V) shift at 100 µM of compound #51 was halved compared with the WT KV7.2 homotetramer, from −12.9 ± 1.9 (n = 4) to −6.6 ± 0.6 mV (n = 3) (Fig. 10, B and C), suggesting that at least part of the G(V)-shifting effect of compound #51 was caused by electrostatic interaction with R1 in the voltage-sensing domain. There was no difference in effect between KV7.2/7.3 and KV7.2 (Fig. 10 C). It is possible that neutralization of the two most extracellularly located gating charges of S4 abolishes channel opening, but this double mutated channel did not express.

In this study, we have shown that a lipophilic compound (4 < LogP < 7) with a COOH group is a common motif among strong activators of the 3R Shaker KV channel. The mechanism of action is possibly common among the compounds, and the data are consistent with a mechanism where the lipophilic compounds are located at the interface between the lipid bilayer and the voltage sensor (Site 1 in Fig. 1 A). From this position, the compound can electrostatically activate the voltage sensor to an up-state to pull the gate open and thereby open the ion-conducting pore.

In our screen, we only identified COOH compounds activating the channel via the lipoelectric mechanism. It is possible that other channel-activating COOH compounds were missing because the screen was designed to find lipoelectric compounds. However, what speaks against this is that on the screen we in fact identified other types of compounds acting via other sites and mechanisms (Liin et al., 2018; Silverå Ejneby et al., 2020). Thus, we suggest that COOH compounds in general act via the lipoelectric site and mechanism and not via other sites and mechanisms.

While this site and mechanism of action on the Shaker KV and KV7.2/7.3 channels clearly overlap with other COOH compounds, such as PUFAs and resin acids (Ottosson et al., 2017; Silverå Ejneby et al., 2021; Yazdi et al., 2021; Yazdi et al., 2016), it is distinct from the sites and mechanisms of action of other channel-activating compounds such as biaryl sulfonamides (NaV 1.7, [Ahuja et al., 2015]; Shaker KV, [Liin et al., 2018]; site 2 in Fig. 1 A), warfarin-like tautomers (Shaker KV and KV1.5, [Silverå Ejneby et al., 2020]; site 3 in Fig. 1 A), and the antiseizure drug retigabine (KV7.2–7.5, [Main et al., 2000]; site 4 in Fig. 1 A). Thus, we propose, as a general mechanism, that COOH-containing compounds in most cases activate KV channel via the lipoelectric mechanism. But, this does not mean that all lipophilic compounds with a COOH group activate all KV channels. As is evident from the screen, most lipophilic COOH compounds do not activate the 3R Shaker KV channel (Fig. 3 B), and some channel-activating compounds were selective for specific channels (e.g., Fig. 7 E). Thus, we suggest that refined selectivity between different types of channels will be possible to develop for these lipoelectric compounds.

The data are available from the corresponding author upon reasonable request.

Christopher J. Lingle served as editor.

We thank Per-Eric Lund for running the high-throughput experiments, Gunnar Nordvall for chemistry advice and compound selection, Anders B. Eriksson (Science for Life Laboratory) for help with compound handling, Andreas Nolting for construction of the cell line, Ulrika Yngve (at the Drug Discovery and Development platform at Science for Life Laboratory) for identification of compounds in the screen, and Johan Brask, Damon Frampton, Nina Ottosson, and Hani Abdelrahman for aid, discussions, and insights during the experimental and analytical process. We thank the Chemical Biology Consortium Sweden (CBCS) at the Science for Life Laboratory (Solna, Sweden) for providing the chemical libraries tested on the screen. The NIH Clinical Collection was provided through the National Institutes of Health Molecular Libraries Roadmap Initiative. The LCBKI primary screening set was provided by CBCS. Part of this work was assisted by Karolinska High Throughput Center (KHTC), a core facility at Karolinska Institutet with affiliation to Science for Life Laboratory (https://www.scilifelab.se/facilities/khtc/).

This work was supported by grants from the Swedish Research Council (2020-01019), the Swedish Heart-Lung Foundation (20210596), the Swedish Brain Foundation (2022-0219), and the County of Östergötland.

Author contributions: O. Rönnelid: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Visualization, Writing—original draft, Writing—review and editing, F. Elinder: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing—original draft, Writing—review and editing.

Ahuja
,
S.
,
S.
Mukund
,
L.
Deng
,
K.
Khakh
,
E.
Chang
,
H.
Ho
,
S.
Shriver
,
C.
Young
,
S.
Lin
,
J.P.
Johnson
Jr.
, et al
.
2015
.
Structural basis of Nav1.7 inhibition by an isoform-selective small-molecule antagonist
.
Science
.
350
:
aac5464
.
Ashcroft
,
F.M.
2000
.
Ion Channels and Disease
.
Academic Press
,
San Diego, CA
.
Börjesson
,
S.I.
, and
F.
Elinder
.
2011
.
An electrostatic potassium channel opener targeting the final voltage sensor transition
.
J. Gen. Physiol.
137
:
563
577
.
Börjesson
,
S.I.
,
S.
Hammarström
, and
F.
Elinder
.
2008
.
Lipoelectric modification of ion channel voltage gating by polyunsaturated fatty acids
.
Biophys. J.
95
:
2242
2253
.
Cui
,
Y.M.
,
X.L.
Liu
,
W.M.
Zhang
,
H.X.
Lin
,
T.
Ohwada
,
K.
Ido
, and
K.
Sawada
.
2016
.
The synthesis and BK channel-opening activity of N-acylaminoalkyloxime derivatives of dehydroabietic acid
.
Bioorg. Med. Chem. Lett.
26
:
283
287
.
Delemotte
,
L.
,
M.
Tarek
,
M.L.
Klein
,
C.
Amaral
, and
W.
Treptow
.
2011
.
Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations
.
Proc. Natl. Acad. Sci. USA
.
108
:
6109
6114
.
Henrion
,
U.
,
J.
Renhorn
,
S.I.
Börjesson
,
E.M.
Nelson
,
C.S.
Schwaiger
,
P.
Bjelkmar
,
B.
Wallner
,
E.
Lindahl
, and
F.
Elinder
.
2012
.
Tracking a complete voltage-sensor cycle with metal-ion bridges
.
Proc. Natl. Acad. Sci. USA
.
109
:
8552
8557
.
Hille
,
B.
1977
.
Local anesthetics: Hydrophilic and hydrophobic pathways for the drug-receptor reaction
.
J. Gen. Physiol.
69
:
497
515
.
Hille
,
B.
2001
.
Ion Channels of Excitable Membranes
. Third edition.
Sinauer Associates, Inc.
,
Sunderland, MA
.
Hoshi
,
T.
,
W.N.
Zagotta
, and
R.W.
Aldrich
.
1990
.
Biophysical and molecular mechanisms of Shaker potassium channel inactivation
.
Science
.
250
:
533
538
.
Imaizumi
,
Y.
,
K.
Sakamoto
,
A.
Yamada
,
A.
Hotta
,
S.
Ohya
,
K.
Muraki
,
M.
Uchiyama
, and
T.
Ohwada
.
2002
.
Molecular basis of pimarane compounds as novel activators of large-conductance Ca(2+)-activated K(+) channel alpha-subunit
.
Mol. Pharmacol.
62
:
836
846
.
Jensen
,
M.O.
,
V.
Jogini
,
D.W.
Borhani
,
A.E.
Leffler
,
R.O.
Dror
, and
D.E.
Shaw
.
2012
.
Mechanism of voltage gating in potassium channels
.
Science
.
336
:
229
233
.
Kamb
,
A.
,
J.
Tseng-Crank
, and
M.A.
Tanouye
.
1988
.
Multiple products of the Drosophila Shaker gene may contribute to potassium channel diversity
.
Neuron
.
1
:
421
430
.
Larsson
,
J.E.
,
U.
Karlsson
,
X.
Wu
, and
S.I.
Liin
.
2020
.
Combining endocannabinoids with retigabine for enhanced M-channel effect and improved KV7 subtype selectivity
.
J. Gen. Physiol.
152
:e202012576.
Liin
,
S.I.
,
U.
Karlsson
,
B.H.
Bentzen
,
N.
Schmitt
, and
F.
Elinder
.
2016
.
Polyunsaturated fatty acids are potent openers of human M-channels expressed in Xenopus laevis oocytes
.
Acta Physiol.
218
:
28
37
.
Liin
,
S.I.
,
P.E.
Lund
,
J.E.
Larsson
,
J.
Brask
,
B.
Wallner
, and
F.
Elinder
.
2018
.
Biaryl sulfonamide motifs up- or down-regulate ion channel activity by activating voltage sensors
.
J. Gen. Physiol.
150
:
1215
1230
.
Long
,
S.B.
,
X.
Tao
,
E.B.
Campbell
, and
R.
MacKinnon
.
2007
.
Atomic structure of a voltage-dependent K+ channel in a lipid membrane-like environment
.
Nature
.
450
:
376
382
.
Main
,
M.J.
,
J.E.
Cryan
,
J.R.
Dupere
,
B.
Cox
,
J.J.
Clare
, and
S.A.
Burbidge
.
2000
.
Modulation of KCNQ2/3 potassium channels by the novel anticonvulsant retigabine
.
Mol. Pharmacol.
58
:
253
262
.
Miceli
,
F.
,
M.V.
Soldovieri
,
P.
Ambrosino
,
L.
Manocchio
,
I.
Mosca
, and
M.
Taglialatela
.
2018
.
Pharmacological targeting of neuronal Kv7.2/3 channels: A focus on chemotypes and receptor sites
.
Curr. Med. Chem.
25
:
2637
2660
.
Ottosson
,
N.E.
,
S.I.
Liin
, and
F.
Elinder
.
2014
.
Drug-induced ion channel opening tuned by the voltage sensor charge profile
.
J. Gen. Physiol.
143
:
173
182
.
Ottosson
,
N.E.
,
M.
Silverå Ejneby
,
X.
Wu
,
A.
Estrada-Mondragón
,
M.
Nilsson
,
U.
Karlsson
,
M.
Schupp
,
S.
Rognant
,
T.A.
Jepps
,
P.
Konradsson
, and
F.
Elinder
.
2021
.
Synthetic resin acid derivatives selectively open the hKV 7.2/7.3 channel and prevent epileptic seizures
.
Epilepsia
.
62
:
1744
1758
.
Ottosson
,
N.E.
,
M.
Silverå Ejneby
,
X.
Wu
,
S.
Yazdi
,
P.
Konradsson
,
E.
Lindahl
, and
F.
Elinder
.
2017
.
A drug pocket at the lipid bilayer-potassium channel interface
.
Sci. Adv.
3
:e1701099.
Ottosson
,
N.E.
,
X.
Wu
,
A.
Nolting
,
U.
Karlsson
,
P.E.
Lund
,
K.
Ruda
,
S.
Svensson
,
P.
Konradsson
, and
F.
Elinder
.
2015
.
Resin-acid derivatives as potent electrostatic openers of voltage-gated K channels and suppressors of neuronal excitability
.
Sci. Rep.
5
:
13278
.
Prole
,
D.L.
,
P.A.
Lima
, and
N.V.
Marrion
.
2003
.
Mechanisms underlying modulation of neuronal KCNQ2/KCNQ3 potassium channels by extracellular protons
.
J. Gen. Physiol.
122
:
775
793
.
Ragsdale
,
D.S.
,
J.C.
McPhee
,
T.
Scheuer
, and
W.A.
Catterall
.
1994
.
Molecular determinants of state-dependent block of Na+ channels by local anesthetics
.
Science
.
265
:
1724
1728
.
Rooney
,
E.K.
,
J.M.
East
,
O.T.
Jones
,
J.
Mcwhirter
,
A.C.
Simmonds
, and
A.G.
Lee
.
1983
.
Interaction of fatty-acids with lipid bilayers
.
Biochim. Biophys. Acta Biomembr.
728
:
159
170
.
Salari
,
S.
,
M.
Silverå Ejneby
,
J.
Brask
, and
F.
Elinder
.
2018
.
Isopimaric acid - a multi-targeting ion channel modulator reducing excitability and arrhythmicity in a spontaneously beating mouse atrial cell line
.
Acta Physiol.
222
:e12895.
Schenzer
,
A.
,
T.
Friedrich
,
M.
Pusch
,
P.
Saftig
,
T.J.
Jentsch
,
J.
Grötzinger
, and
M.
Schwake
.
2005
.
Molecular determinants of KCNQ (Kv7) K+ channel sensitivity to the anticonvulsant retigabine
.
J. Neurosci.
25
:
5051
5060
.
Silverå Ejneby
,
M.
,
A.
Gromova
,
N.E.
Ottosson
,
S.
Borg
,
A.
Estrada-Mondragón
,
S.
Yazdi
,
P.
Apostolakis
,
F.
Elinder
, and
L.
Delemotte
.
2021
.
Resin-acid derivatives bind to multiple sites on the voltage-sensor domain of the Shaker potassium channel
.
J. Gen. Physiol.
153
:e202012676.
Silverå Ejneby
,
M.
,
B.
Wallner
, and
F.
Elinder
.
2020
.
Coupling stabilizers open KV1-type potassium channels
.
Proc. Natl. Acad. Sci. USA
.
117
:
27016
27021
.
Silverå Ejneby
,
M.
,
X.
Wu
,
N.E.
Ottosson
,
E.P.
Münger
,
I.
Lundström
,
P.
Konradsson
, and
F.
Elinder
.
2018
.
Atom-by-atom tuning of the electrostatic potassium-channel modulator dehydroabietic acid
.
J. Gen. Physiol.
150
:
731
750
.
Strichartz
,
G.
1976
.
Molecular mechanisms of nerve block by local anesthetics
.
Anesthesiology
.
45
:
421
441
.
Yazdi
,
S.
,
J.
Nikesjö
,
W.
Miranda
,
V.
Corradi
,
D.P.
Tieleman
,
S.Y.
Noskov
,
H.P.
Larsson
, and
S.I.
Liin
.
2021
.
Identification of PUFA interaction sites on the cardiac potassium channel KCNQ1
.
J. Gen. Physiol.
153
:e202012850.
Wuttke
,
T.V.
,
G.
Seebohm
,
S.
Bail
,
S.
Maljevic
, and
H.
Lerche
.
2005
.
The new anticonvulsant retigabine favors voltage-dependent opening of the Kv7.2 (KCNQ2) channel by binding to its activation gate
.
Mol. Pharmacol
.
67
:
1009
1017
.
Yazdi
,
S.
,
M.
Stein
,
F.
Elinder
,
M.
Andersson
, and
E.
Lindahl
.
2016
.
The molecular basis of polyunsaturated fatty acid interactions with the shaker voltage-gated potassium channel
.
PLOS Comput. Biol.
12
:e1004704.
Zhou
,
M.
,
J.H.
Morais-Cabral
,
S.
Mann
, and
R.
MacKinnon
.
2001
.
Potassium channel receptor site for the inactivation gate and quaternary amine inhibitors
.
Nature
.
411
:
657
661
.

Author notes

Disclosures: The authors declare no competing interests exist.

This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms/). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 4.0 International license, as described at https://creativecommons.org/licenses/by-nc-sa/4.0/).