Introduction

Epilepsy is a common, heterogeneous and often debilitating disorder, affecting 1% of the global population. Current clinical management of patients suffering from epilepsy is primarily via medication, although it is generally accepted that 30% of newly diagnosed epilepsy patients will be resistant to common anticonvulsant drugs1,2. Moreover, such resistance is correlated to seizure type—for example, it has been reported that 60% of patients suffering from focal epilepsy will develop and retain resistance to medications3. These shortcomings have motivated the search for new antiepileptic drugs, non-conventional treatments for epilepsy (such as the ketogenic diet) and the identification and ongoing validation of new drug targets4,5,6,7.

Pharmacotherapy in epilepsy has been primarily directed towards a small subset of ion channel targets8. Overall, the vast majority of drugs developed as anticonvulsants are targeted towards voltage-gated Na+ channels or GABA-mediated neurotransmission. However, several recently developed drugs have begun to expand the diversity of molecular targets for epilepsy. A good example is the ‘first-in-class’ Kv channel opener retigabine, presently approved for use as an adjunct therapy to antiepileptic drugs for partial/focal seizures9. Among all the antiepileptic drugs in use, retigabine has a unique mechanism—it is the only compound approved for human use that acts by activating/opening Kv channels. The resulting suppression of electrical excitation likely underlies its anticonvulsant and analgesic properties10,11. Numerous studies have demonstrated that this drug class is effective and well tolerated over a wide therapeutic window, with comparable effectiveness to many classic anticonvulsants in animal models7,12,13. However, despite the therapeutic potential for management of epilepsy and disorders involving membrane excitability, there remains a paucity of molecular understanding of the binding mode(s) and mechanism of action of this new drug class.

Retigabine acts primarily by opening KCNQ2–5 (Kv7.2–7.5) channels11,14,15. In the nervous system, KCNQ2–5 channels assemble as heteromers to generate the ‘M-current’, a subthreshold voltage-gated K+ conductance that is influenced by muscarinic receptor signalling (primarily via PIP2) and regulates neuronal excitability14,16. Retigabine activates these channels by shifting their voltage dependence to hyperpolarized membrane potentials, effectively causing them to open at voltages where they would normally adopt a closed, non-conducting conformation17,18. A conserved tryptophan residue in the pore-forming S5 helix of retigabine-sensitive KCNQ channels (KCNQ2 Trp236, KCNQ3 Trp265) is known to be essential for retigabine effects and is absent in the retigabine-resistant cardiac isoform KCNQ1 (refs 19, 20). The identification of retigabine as a potent Kv channel opener has led to numerous compound library screens to generate novel analogues with altered subunit specificity and effectiveness9,21,22,23,24,25. However, the molecular details that underlie retigabine–KCNQ channel interactions have remained elusive because of the inherent lack of resolution of conventional site-directed mutagenesis techniques.

We have used unnatural amino-acid mutagenesis to subtly rearrange atoms and electrons of the Trp side chain in the retigabine binding site. Our findings demonstrate that repositioning of the indole nitrogen atom of KCNQ3 Trp265 completely abolishes retigabine effects, suggesting that retigabine interaction requires the formation of a H-bond with Trp265. The importance of this H-bond interaction is further illustrated by introduction of fluorinated Trp analogues (with increased H-bonding propensity) at KCNQ3 Trp265, causing increased drug potency. Using multiple retigabine analogues, we identify an electrostatic fingerprint around a H-bond acceptor that correlates with drug potency. These experimental constraints pinpoint specific atoms and chemical forces involved in retigabine interactions, and highlight approaches that may be used to guide rational improvement of existing drugs.

Results

Retigabine modulation of KCNQ channels through an S5 Trp

Retigabine strongly affects voltage-dependent gating of KCNQ2 and KCNQ3* channels (Fig. 1a,b) by generating a substantial hyperpolarizing shift of the V1/2 of activation by 40–60 mV (KCNQ3* refers to KCNQ3[Ala315Thr]—see Methods). These effects have also been demonstrated in KCNQ4 and KCNQ5, but are absent in KCNQ1 because of the absence of an essential Trp residue (Trp236 in KCNQ2, Fig. 1a; Trp 265 in KCNQ3, Fig. 1b: Leu246 in KCNQ1)14,26. This Trp side chain lies in the pore-forming S5 transmembrane segment, near the intracellular voltage-operated gate of the channel. As KCNQ2–5 subunits generally assemble as heteromers in the central nervous system14,16, we tested the effects of retigabine in oocytes co-injected with KCNQ2 and KCNQ3 and observed large shifts of activation to more negative voltages, although not quite as large as with KCNQ2 or KCNQ3 alone (Fig. 1c). Some KCNQ channel activator molecules also cause a marked increase in peak current; however, the effects of retigabine on KCNQ3* channels are quite modest (rarely greater than a 15% increase in peak current), and throughout this study we have focused on the large gating shifts observed in these channels.

Figure 1: Multiple retigabine molecules modulate KCNQ2 and KCNQ3 channel subunits via an S5 Trp side chain.
figure 1

(a,b) Conductance–voltage relationships for (a) KCNQ2 (n=3) and KCNQ2[Trp236Phe] (n=6), and (b) KCNQ3* (n=5) and KCNQ3*[Trp265Phe] (n=3) homomeric channels along with indicated mutants (retigabine concentration of 100 μM). (c) Conductance–voltage relationships for heteromeric combinations of KCNQ2 and KCNQ3 (1:1 ratio of injected mRNA, with or without Trp→Phe mutations as indicated, n=5 for each combination), used to generate channels with reduced numbers of retigabine binding sites. (d) Summary of V1/2 shifts in saturating 100 μM retigabine for mutations of KCNQ2 Trp236 and KCNQ3 Trp265 as indicated (*P<0.05 in a paired Students t-test comparing control versus 100 μM retigabine in each experimental oocyte, n=3–6 per mutant). Only a Trp at either position is sufficient for retigabine sensitivity. (e) Exemplar currents of KCNQ3* and KCNQ3*[Trp265Phe] mutant coexpressed with CiVSP, illustrating that the Trp side chain responsible for retigabine sensitivity is not required for PIP2 sensitivity. (f) Summary data of tail current magnitude (−20 mV) after prepulses to a range of voltages, in oocytes expressing KCNQ3* (n=5) or KCNQ3*[Trp265Phe] (n=5) channels, along with CiVSP. In all panels, error bars represent s.e.m.

Heteromeric KCNQ2/3 channels were used to alter the number of available retigabine binding sites, by co-injecting KCNQ2[Trp236Phe] with wild-type (WT) KCNQ3, or KCNQ3[Trp265Phe] with WT KCNQ2. Assuming 1:1 stoichiometry of KCNQ2/3 assembly, both conditions are estimated to eliminate roughly half of the retigabine binding sites. Consistent with this, we observed intermediate effects of retigabine with fewer available binding sites, indicating that multiple molecules of the drug may be required to generate the full effect (Fig. 1c). A surprising, currently unexplained observation, was that the retigabine-mediated gating shift in KCNQ2/3 heteromers (−40 mV) was noticeably smaller than the shift in KCNQ2 or KCNQ3* homomeric channels. This likely hints that determination of the exact details of stoichiometry and cooperativity of drug binding in heteromeric channels will require alternative experimental approaches. We measured the effects of various mutations at position Trp265 (in KCNQ3*) or Trp236 (in KCNQ2) on retigabine sensitivity (Fig. 1d). In both KCNQ2 and KCNQ3* channels, only a Trp side chain present at this position was sufficient to enable a significant gating shift by retigabine, and even conservative mutations to other aromatic side chains (Phe or Tyr, Fig. 1a,b,d) cannot replicate the retigabine-mediated gating shift, consistent with previous reports19,20.

We tested whether KCNQ3 position Trp 265 plays a role in regulating channel responses to PIP2. Using the voltage-sensitive phosphatase ciona intestinalis voltage-sensitive phosphatase (CiVSP) to hydrolyse PIP2 at depolarized voltages, we pulsed oocytes through a range of prepulse voltages followed by a test pulse to −20 mV to assess channel activity after PIP2 depletion27. We observed similar channel rundown with depolarizing voltage steps in both KCNQ3* and KCNQ3*[Trp265Phe] channels, indicating that Trp265 does not significantly influence the effects of anionic phospholipids on channel function (Fig. 1e,f). Overall, these findings illustrate that KCNQ3 residue Trp265 is essential for the effects of retigabine, but does not play a significant role in regulating channel gating by voltage or PIP2.

Unnatural amino-acid mutagenesis of KCNQ3 channels

To investigate the underlying mechanism of retigabine interactions with this essential Trp side chain, we employed unnatural amino-acid mutagenesis to introduce subtly altered Trp variants (Fig. 2a). With this method, a stop codon (TAG) is placed in the ion channel gene at a site of interest, and this mRNA is co-injected into Xenopus laevis oocytes along with a synthetic amino-acylated tRNA (carrying an unnatural amino acid) that is orthogonal to Xenopus tRNA synthetic pathways28. When the ribosomal translation machinery encounters the introduced TAG stop codon, the complementary synthetic tRNA (with a CUA anticodon) enables readthrough and incorporation of the appended amino acid into full-length ion channels (Fig. 2a)29.

Figure 2: Nonsense suppression for amino-acid incorporation in KCNQ3* channels.
figure 2

(a) Schematic diagram of the nonsense suppression method, in which mRNA (with a stop codon at Trp265) and amino-acylated tRNA are co-injected into Xenopus oocytes. Incorporation of the unnatural amino acid enables readthrough of the stop codon and expression of functional channels. (b) Current magnitudes in oocytes injected with KCNQ3*[Trp265TAG] mRNA and either an unconjugated tRNA (pdCpA; n=5) or tRNA amino-acylated with Trp (n=8, *P<0.05, Student’s t-test). (c) Activation kinetics of KCNQ3* (n=5) and KCNQ3*[Trp265TAG] (n=8) channels rescued with Trp, in the presence or absence of 100 μM retigabine (*P<0.05, Student’s t-test). (df) Exemplar currents from oocytes injected with KCNQ3*[Trp265TAG] mRNA, and indicated synthetic tRNAs. (g) Conductance–voltage relationships for Trp-rescued KCNQ3*[Trp265TAG] (n=8), with retigabine response, illustrating faithful incorporation of the desired side chain at position 265. For KCNQ3* channels V1/2=−44±1 mV, k=7.5±0.5 mV; for Trp-rescued KCNQ3*[Trp265TAG] V1/2=−43±2 mV, k=7.9±0.5 mV (± indicates s.e.). Error bars throughout represent s.e.m.

Among the KCNQ channels we tested, KCNQ3* channels exhibited large currents when expressed as homomeric channels in Xenopus oocytes, with a short latency between injection and expression, and were found to be an effective model system for incorporation of unnatural amino acids. KCNQ3* is also a useful model for retigabine binding because it does not possess a second interaction site that is targeted by some other recently studied KCNQ2 activators (see Discussion)22,30, and so our results are specifically focused on the retigabine binding site formed by Trp265. Several criteria illustrate the feasibility and fidelity of the nonsense suppression approach in KCNQ3* channels. First, co-injection of KCNQ3*[Trp265TAG] with a full-length tRNA lacking a conjugated amino acid (pdCpA control) resulted in virtually no ionic currents, even several days after injection (Fig. 2b,d). However, when KCNQ3*[Trp265TAG] channels were co-injected with Trp-conjugated tRNA, robust currents with WT-like gating were observed, indicating efficient suppression of the introduced TAG stop codon and synthesis of functional channels with WT-like properties (Fig. 2b,e). Second, although there may be differences in the efficiency of incorporation of unnatural amino acids at this position, we have restricted our investigation to homomeric channels (rather than KCNQ2/KCNQ3 heteromers), and to biophysical parameters (V1/2, activation kinetics) that are normally independent of channel expression. Last, we are confident that the synthetic Trp-tRNA was incorporated in a significant fraction of channel subunits because elicited currents responded to retigabine (Fig. 2c,f,g), and only Trp appears to be sufficient for retigabine sensitivity (Fig. 1d).

Loss of retigabine effects after removal of Trp265 H-bonding

Next, we tested retigabine sensitivity of KCNQ3* channels synthesized with the ‘Ind’ amino acid at position 265. This unnatural amino-acid side chain ablates the hydrogen-bonding capability of Trp by shifting the position of the indole nitrogen atom (Fig. 3a)31, but retains the steric and chemical compositions of Trp. KCNQ3*[Trp265TAG] currents were efficiently rescued with Ind relative to the pdCpA control, yielding relatively large voltage-activated K+ currents that closely resembled KCNQ3* (Fig. 3b,c). Remarkably, however, this single atom alteration in each of the KCNQ3* subunits entirely abolishes retigabine activation of KCNQ3* channels (Fig. 3d,e). Notably, Trp265Ind channels displayed WT-like gating, indicating that this residue does not play a significant role in conformational stabilization of channel states in the absence of drug. However, given the potent impact of such a minor manipulation, we examined the effects of the Trp265Ind substitution on channel gating more closely, to rule out the possibility of significant perturbation of channel function. We observed no statistically significant difference in the V1/2 or slope of the voltage dependence of activation of KCNQ3* and KCNQ3*[Trp265Ind] (see Table 1 and Fig. 3d,e). Last, as described earlier, the more disruptive Trp265Phe mutation did not alter the inhibition of KCNQ3* currents by PIP2 depletion (mediated by CiVSP). Taken together, these observations demonstrate a highly specific effect of the Trp265Ind substitution on drug interactions (with little disruption of intrinsic voltage-dependent gating or lipid regulation), indicating that H-bond formation with Trp265 is a crucial step for retigabine action.

Figure 3: The position of the Trp 265 indole nitrogen is essential for retigabine sensitivity.
figure 3

(a) Chemical structures of Trp and Ind side chains, illustrating the subtle change in the position of the indole nitrogen atom. (b) Exemplar currents elicited from a Xenopus oocyte with Ind-rescued KCNQ3*[Trp265TAG] channels illustrating retigabine insensitivity. (c) Current magnitudes in oocytes injected with KCNQ3*[Trp265TAG] mRNA and either an unconjugated tRNA (pdCpA; n=5) or tRNA amino-acylated with Ind (n=7, *P<0.05, Student’s t-test). (d) Conductance–voltage relationships for Ind-rescued KCNQ3*[Trp265TAG], in the presence and absence of retigabine, illustrating the importance of the correct positioning of the N–H group. For KCNQ3*[Trp265Trp], V1/2=−43±2 mV, k=7.9±0.5 mV; for KCNQ3*[Trp265Ind], V1/2=−48±2 mV, k=7.3±0.6 mV (no statistical significance, ±indicates s.e.). (e) Activation kinetics for KCNQ3*[Trp265Ind] measured at −20 mV in the presence and absence of 100 μM retigabine (n=7, no statistical significance). In all panels, error bars represent s.e.m.

Table 1 Activation gating parameters for unnatural amino-acid substitutions of KCNQ3*[Trp265TAG] channels.

Tuning the strength of the Trp265–retigabine interaction

To investigate the possible involvement of other modalities of Trp interactions with the drug, we also examined the consequences of fluorination of Trp265 on retigabine effects. Trp fluorination is typically used to modify the electrostatic surface potential and test for the effects of cation–pi interactions (Fig. 4a)32. We were unable to detect functional channels carrying the F4-Trp side chain (fluorines at ring positions 4,5,6 and 7, as numbered in Fig. 4a); however, robust current rescue was observed for less fluorinated derivatives such as F3-Trp (Fig. 4b,c). Retigabine potentiation of KCNQ3* channels was retained with F3-Trp substitution at Trp265 (Fig. 4c–e), a substitution that still potently diminishes the negative electrostatic potential on the face of the side chain (Fig. 4a), indicating that a cation–pi interaction of Trp265 with retigabine or some other entity (perhaps another channel residue) is not required for retigabine effects. It is noteworthy that progressive fluorination of Trp265 caused variable effects on KCNQ3*-gating properties (Table 1), with F1-Trp causing a modest depolarizing shift of the activation V1/2, while F3-Trp significantly shifted gating in the hyperpolarizing direction. This lack of an obvious trend also suggests that an endogenous cation–pi interaction involving Trp265 does not contribute significantly to the gating process.

Figure 4: The polarity of the Trp265 indole nitrogen modulates retigabine sensitivity.
figure 4

(a) Chemical structures of Trp (with ring positions labelled) and F3-Trp, accompanied by colorimetric representations of electrostatic surface potentials. (b) Current magnitudes in oocytes injected with KCNQ3*[Trp265TAG] mRNA and either an unconjugated tRNA (pdCpA; n=5) or tRNA amino acylated with F3-Trp (n=12, *P<0.05, Student’s t-test). (c) Exemplar currents from F3-Trp-rescued KCNQ3*[Trp265TAG] channels in the presence and absence of retigabine. (d) Conductance–voltage relationships illustrating the effects of 1 μM retigabine on Trp-rescued (n=12) and F3-Trp-rescued KCNQ3*[Trp265TAG] (n=12) channels. (e) Activation kinetics (−20 mV) for F3-Trp-rescued channels (n=3), in the presence and absence of 100 μM retigabine (*P<0.05, Student’s t-test). (f) Concentration–response curves for retigabine effects on numerous fluoro-Trp analogues (n=9–12 per Trp analogue) substituted at position Trp265, illustrating enhanced retigabine potency with increased fluorination. (g) Conductance–voltage relationships with indicated retigabine concentrations on a KCNQ3 mutant (Asn220Cys, n=4) with an intrinsic hyperpolarizing shift in gating. Conductance–voltage relationships for F3-Trp substituted at Trp265TAG are shown for comparison. (h) Concentration–response curves for retigabine effects on KCNQ3* (n=5) and KCNQ3*[Asn220Cys] (n=4) channels. In all panels, error bars represent s.e.m.

The most interesting outcome of fluoro-Trp substitutions further highlights the role of H-bond formation by the indole N–H group. Specifically, fluorination of the Trp side chain significantly increases the polarity of the indole N–H, thereby strengthening the hydrogen-bonding propensity of the indole side chain33,34. This is illustrated by the electrostatic surface potentials for Trp and F3-Trp (note the more polarized electrostatic surface potential of F3-Trp, highlighted by the blue colour around the N–H bond, in the lower right corner as oriented in Fig. 4a). Consistent with this, retigabine potency was increased with fluoro-Trp analogues at Trp265. For instance, in KCNQ3*[Trp265TAG] channels rescued with Trp (mimicking KCNQ3* channels), 1 μM retigabine elicits very little shift in the voltage dependence of activation (Fig. 4d). However, in channels with F3-Trp substituted at Trp265, 1 μM retigabine elicits a much larger response (Fig. 4d). Consistent with this effect being attributable to strengthening of the formation of a hydrogen bond by Trp265, we observed a fluorination-dependent increase in retigabine potency, with progressive fluorination shifting the concentration response to lower retigabine concentrations (Fig. 4f). In addition, we observed a small but significant position-dependent effect of mono-fluorination of Trp265, in which fluorination of ring position 7 (closer to the indole N–H) increased retigabine effects slightly more than fluorination of ring position 5 (Fig. 4f). This is consistent with the fluorine substituent at position 7 having a stronger effect on the polarity of the N–H bond.

Since F3-Trp incorporation itself caused a hyperpolarizing shift of channel activation even in the absence of retigabine (Table 1 and Fig. 4d), a possible alternative explanation was that increased retigabine potency with F3-Trp at position 265 was a secondary allosteric effect of stabilization of the channel open state, rather than a direct effect on retigabine binding. Thus, as an additional control experiment, we investigated the effects of the Asn220Cys mutation on retigabine sensitivity (in the extracellular S3–S4 linker, distant from the Trp265 putative retigabine binding residue). This mutation shifts channel opening to more negative voltages, even more than F3-Trp substitution at position 265. In the Asn220Cys mutant, no sensitization to retigabine was observed (Fig. 4g). That is, 1 μM retigabine was unable to generate a marked gating shift, in contrast to the potent retigabine effects with F3-Trp substitution at Trp265 (curve fits for the F3-Trp substitution are included for comparison, Fig. 4g). In addition, the retigabine concentration response for Asn220Cys channels was very similar (albeit shallower) to the KCNQ3* background construct, demonstrating that the Asn220Cys mutation per se does not inhibit retigabine activation (Fig. 4h). These observations suggest that an indirect effect of open-state stabilization is unlikely to account for increased retigabine potency in F3-Trp-substituted channels.

Relative contributions of putative drug-binding residues

Given the dramatic effects of the Trp265Ind substitution, we sought to clarify the contributions of other residues reported to play a role in retigabine binding. Using a chimeric approach between KCNQ1 and KCNQ3, previous work identified KCNQ3 residues Thr271 (S5 helix), Leu272 (S5 helix), Leu314 (pore helix) and Leu338 (S6 helix) as important contributors to a putative retigabine binding pocket19. We generated several mutations at each of these positions in KCNQ3* and tested retigabine effects over a broad concentration range (Fig. 5a,b). We observed that all functional mutants at these positions retained a large retigabine-mediated gating shift (Fig. 5b), although reduced potency was reflected by a shift to higher drug concentrations in many cases. Moreover, many of these mutations caused intrinsic shifts in channel gating in the absence of retigabine (Fig. 5a). Previous reports have suggested that many of these mutations strongly diminish retigabine action; however, our application of higher concentrations illustrates that, while potency of the drug is weakened, large gating shifts still occur. These data indicate that Thr271, Leu272, Leu314 and L338 make measurable contributions to retigabine binding affinity. However, the Trp265Ind substitution abolishes any response to retigabine at experimentally achievable concentrations. Thus, Trp265Ind remains the most structurally subtle, yet most disruptive, mutation we have identified in terms of retigabine sensitivity, highlighting the potential importance of H-bonding with the Trp265 side chain.

Figure 5: Detailed characterization of secondary retigabine binding residues and alternative binding site orientations.
figure 5

(a) Conductance–voltage relationships were gathered for the indicated KCNQ3* mutant channels (n=4–6 per mutant), in 0, 100 or 300 μM retigabine. (b) Maximal ΔV1/2 in 300 μM retigabine measured in each mutant channel. Error bars in a,b represent s.e.m. (c) Retigabine was docked into a molecular model of the pore-forming domain of KCNQ3 (see ref. 19). Two orientations are shown with the carbamate group in either the vicinity of Leu314 (‘original’ model) or Trp265 (‘flip’ model). The two binding models are superimposed in the ‘overlay’, showing the similar space occupied by both drug orientations.

We next used a previously published molecular model of the retigabine binding site to investigate alternative binding modes of the drug19. In the published model (‘original’ in Fig. 5c) there are no obvious H-bond acceptors in the vicinity of Trp265. There are also no apparent H-bond interactions with Thr271. However, this putative binding pocket can accommodate altered orientations (for example, ‘flip’ in Fig. 5c) in which the drug occupies a similar space (see ‘overlay’ panel illustrating both binding models in Fig. 5c), but orients the carbamate moiety of retigabine in close proximity to Trp 265. A rotameric shift of Trp265 enables a close interaction between the indole N–H and a carbamate oxygen atom of retigabine. Using molecular dynamics simulations (of 12 retigabine binding sites starting in a ‘flip’ orientation, for 10 ns each), we measured the mean distance between the retigabine carbamate oxygen and indole N–H group at the end of the simulations as 4.0±0.5 Å (±indicates s.e., four binding sites were measured between 2 and 3 Å, three sites between 3 and 4 Å and five sites were >5 Å), consistent with the possible formation of a hydrogen bond in this region.

Identification of a likely H-bond donor in retigabine

We further investigated the mechanism of H-bonding with Trp265 by seeking potential H-bond acceptors in the retigabine molecule, a task facilitated by the availability of numerous retigabine analogues. We first considered the analogue ML-213 because it has a simplified chemical scaffold compared with retigabine and a reduced number of possible H-bond acceptors (Fig. 6a,b). We found that ML-213 is a more potent activator of KCNQ3* compared with retigabine (EC50=3.6±0.2 μM, n=5, versus 11.6±0.4 μM, n=5, respectively, Fig. 6c, ± indicates s.e.). We have also included the calculated electrostatic surface potentials for retigabine and ML-213, illustrating that the negative surface potential around the carbonyl oxygen is more pronounced for ML-213 (red protrusion, Fig. 6a,b).

Figure 6: ML-213 exhibits a stronger electrostatic surface potential and higher potency than retigabine for KCNQ3* activation.
figure 6

(a,b) Chemical structures and electrostatic surface potentials for retigabine and ML-213. Note the increased negative surface potential in the vicinity of the carbonyl oxygen atom in ML-213. The scale for electrostatic surface potential representation is red: −80 kcal, yellow: 0 kcal, blue: +80 kcal. (c) Concentration–response curves for retigabine and ML-213 (n=5) effects on KCNQ3* channels (n=5). Error bars represent s.e.m.

Expanding on these observations, we examined a spectrum of KCNQ channel activators with varying potency in KCNQ3* channels (Fig. 7). We have depicted all of the drugs tested and colorimetric representations of their calculated electrostatic surface potentials, along with the observed EC50 for the shift of activation V1/2 of KCNQ3*[Trp265TAG] channels (rescued with Trp, F3-Trp or Ind). For all drugs tested, effects were abolished by the Trp265Ind substitution, indicating a common mechanism via H-bond formation with Trp265. In addition, F3-Trp substitution generally resulted in higher drug potency, although this effect was not as pronounced with the less polar analogues ICA-110381 and ICA-069673. Importantly, all of the activator compounds contain a carbonyl oxygen (either in the carbamate moiety of retigabine and flupirtine, or in the amide linker in ML-213 or the ICA compounds), and it is noteworthy that the strength of the negative surface potential correlates well with the potency of the drug for KCNQ3* activation—drugs with a weaker electrostatic surface potential trend towards weaker potency (Fig. 7). The ICA compounds are very weakly potent activators of KCNQ3* channels within our typical experimental concentration range (up to 300 μM), and notably ICA-110381 caused considerably more activation in F3-Trp-substituted channels, demonstrating that effects of an otherwise weak channel activator can be strengthened by targeting the chemical properties of Trp265. Taken together, these findings illustrate the importance of the H-bonding propensity of Trp265, likely involved in the formation of a H-bond with a carbonyl oxygen present in retigabine and its analogues.

Figure 7: Effects of retigabine analogues correlate with electrostatic surface potential.
figure 7

(a) Chemical structures and electrostatic surface potentials for a series of retigabine analogues. All structures and surface potential maps have been aligned on the basis of the location of the conserved amide-ester bond—note the gradient of the intensity of the negative surface potential around the carbonyl oxygen atom (scaling is the same as in Fig. 6). (b) Summary illustrating the EC50 of each drug on Trp-rescued, F3-Trp-rescued and Ind-rescued KCNQ3*[Trp265TAG] channels (n=4–9 per data point), and the maximal efficacy (ΔV1/2) of each drug in F3-Trp and Trp-rescued channels (effects on Ind-rescued channels are minimal and thus have been omitted). Error bars represent s.e.m.

Discussion

Retigabine is currently the only drug widely approved for clinical use that acts via activation of voltage-gated potassium channels35. Detailed investigation of this drug prototype may deepen our understanding of how drugs can act as ion channel activators, and accelerate the development of potassium channel activators with increased potency or specificity. Such drugs may have untapped therapeutic potential for management of a variety of ‘hyperexcitability’ disorders in addition to epilepsy, such as painful syndromes, hypertension and cardiac arrhythmias9,36. In addition, recent reports have highlighted retigabine effectiveness in treatment of tinnitus, Parkinson’s disease and Huntington’s disease37,38,39.

At the outset of our study, chimeric studies with the retigabine-insensitive KCNQ1 channel had identified the importance of a Trp residue (KCNQ3 Trp265) in the pore-forming region, near the presumed cytoplasmic channel gate. Retigabine interactions with this site had been rationalized as generic hydrophobic interactions with this conserved Trp residue (also present in KCNQ2, 4 and 5)19,20.

Although most often considered to be hydrophobic, Trp side chains can deploy diverse chemical forces in ligand interactions40, and these require specialized chemical-scale approaches to be distinguished29. There is growing recognition of the general importance of electrostatic cation–pi interactions with all aromatic side chains29,41. In addition, both Trp and Tyr side chains can act as H-bond donors. Unnatural amino-acid mutagenesis approaches allow us to parse out atomistic modes of interaction, and also hopefully mitigate the possibility of gross structural perturbations that might arise using conventional mutagenesis. By engineering KCNQ3 channels to include subtle Trp variants at position 265, our study identifies the very atoms required for retigabine effects and, further, the chemical forces involved in forming a direct interaction with retigabine. Surprisingly, our findings indicate that the critical property of Trp265 for retigabine activation, rather than hydrophobicity, is the ability to H-bond. Repositioning of the indole nitrogen atom to remove Trp H-bonding abolished the effects of retigabine and all analogues tested (Figs 3 and 7). Moreover, a Trp analogue (F3-Trp) with altered electrostatics but intact/enhanced H-bonding ability enhances retigabine potency (Figs 4 and 7). A valuable future direction for these experiments will be to develop an independent measure of drug binding, to distinguish whether the H-bond interaction (disrupted by the Trp265Ind mutation) is essential for drug binding, or couples drug binding to a later step involved in channel potentiation. Nevertheless, these observations highlight a rarely investigated mode of drug interactions with Trp, as Trp residues in a binding pocket are often presumed to indicate a hydrophobic site. Our findings illustrate the importance of detailed investigation of the underlying chemistry of these binding sites to better understand mechanisms of drug action.

By demonstrating the effects of H-bond strength on a range of KCNQ3 activators, our findings suggest the importance of a carbonyl oxygen (usually in a carbamate or amide moiety) for the formation of a negative electrostatic surface potential to act as a H-bond acceptor. Previous screens of compound libraries have investigated important physicochemical features important for drug activity on KCNQ2–5 channels. Several studies have reported the importance of an amide bond (or a carbonyl oxygen) as an essential element of the pharmacophore, although the relationship between this functional group and the H-bonding propensity of Trp265 has not been recognized9,42,43,44. Notably, KCNQ activators with marked structural diversity appear to act through Trp265 (several examples are shown in Fig. 8), and these share the common feature of a carbonyl oxygen that we suggest acts as an essential hydrogen bond acceptor. Supporting this notion, our findings illustrate that ML-213 is an effective KCNQ3 activator, correlated with its strong surface potential relative to other drugs tested. In the context of understanding drug interactions with ion channels and other receptors, these observations illustrate the importance of investigating specific chemical forces that enable drug interactions with aromatic side chains known to contribute to binding sites in other ion channels (for example, voltage-gated Na+ channels, human ether-a-go-go related gene product (hERG))45,46,47. These insights could guide rational ‘tuning’ of the properties of certain functional groups to alter drug properties in desirable ways.

Figure 8: Diverse structures of KCNQ openers.
figure 8

Multiple structures of KCNQ channel openers are presented to highlight the overall features of an amide group flanked by various ring structures. Our findings highlight the importance of the amide carbonyl for interaction with KCNQ3 Trp 265 and likely equivalent positions in KCNQ2, 4 and 5. Drugs depicted are (a) retigabine, (b) ztz-240 (described in ref. 24, c) acrylamide (s)-1, (d) BMS-204352 and (e) an unnamed experimental drug described in ref. 43.

On the basis of the magnitude of the retigabine-induced shift of the conductance–voltage relationship, saturating effects of retigabine are estimated to stabilize channel opening by 6 kcal mol−1 (calculated as ΔG=Δ(zFV1/2)). As a comparison, formation of a H-bond between a carbonyl and amino groups (N-H: O=C) is generally quite weak, on the order of 2 kcal mol−1. A tetrameric KCNQ channel likely contains four potential retigabine binding sites, and our characterization of heteromeric channels with fewer binding sites (Fig. 1c) indicates that multiple sites may be occupied in saturating retigabine concentrations. Thus, the energetics of the full retigabine-mediated gating shift can be reasonably attributed to multiple retigabine molecules interacting with Trp265 via a H-bond. Although not as prominent as the Trp265 indole nitrogen, other residues also make additional measurable contributions to the overall energy of retigabine binding (Fig. 5 and ref. 19). It should also be noted that cooperativity between retigabine binding sites has not been investigated in detail, and so it is not known how multiple retigabine binding events may interact and influence channel function. Investigation of libraries of retigabine analogues has highlighted that variable aromatic ring substituents can alter the apparent affinity, maximal effect and subunit specificity24,43. For example, many ring substituents (such as halogenation) alter the polarity of the carbonyl oxygen and thereby alter drug interactions with Trp265. In addition, it seems likely that diverse flanking structures (potentially by interacting with the auxiliary retigabine binding residues investigated in Fig. 5) could influence the orientation of the carbonyl oxygen relative to Trp265.

A final noteworthy detail is that, although our findings delineate an atomic basis for the retigabine interaction with KCNQ3 Trp265 (and analogous residues in KCNQ2,4,5), there is significant evidence for additional interaction sites that enable potentiation of certain KCNQ subunits. For example, the ICA compounds tested (and other related compounds such as ztz-240), along with a family of diclofenac derivatives, have been reported to influence KCNQ2 channels independent of the conserved S5 Trp residue and exhibit dramatically stronger effects in KCNQ2 relative to KCNQ3 (refs 22, 24, 25, 30). In addition, unrelated KCNQ openers such as zinc pyrithione do not appear to depend on the hydrogen bond interaction that we have identified in this study48. Therefore, ongoing examination of drug binding to different KCNQ channel subtypes may help to reveal additional sites for KCNQ potentiation.

In conclusion, we have used unnatural amino-acid mutagenesis and available retigabine analogues to localize the effects of retigabine binding to a single H-bond interaction with KCNQ3 Trp265. These findings highlight an unusual mode of drug interaction with a Trp side chain, and demonstrate the importance of careful investigation of the mechanism of drug interactions to guide rational modification of therapeutic compounds.

Methods

Molecular biology and in vivo nonsense suppression

KCNQ2 (human) and KCNQ3 (human) channel cDNAs were propagated and manipulated using the pTLN vector (gifts of Dr M. Taglialatela and Dr T. Jentsch). Introduction of point mutations was accomplished using PCR to combine two fragments each of which contained the mutation in an overlapping region of the sequence (‘two-step’ PCR method). To accomplish this, we used standard PCR approaches to amplify a 5′ fragment and a 3′ fragment (using a ‘flanking primer’ paired with a mutagenic primer). These fragments were then mixed along with the ‘flanking primer’ pair for a second round of PCR to combine the overlapping fragments. The resulting fragment was then subcloned into the appropriate parent vector using EcoRI and NotI restriction enzymes, and the sequence was verified by Sanger sequencing approaches (Genewiz). Sequences for the flanking primers and the ‘forward’ mutagenic primers are provided in Supplementary Table 1 for all mutations that were generated (‘reverse’ mutagenic primers are simply the reverse complement). In all experiments involving homomeric expression of KCNQ3 channels (including in vivo nonsense suppression experiments), the Ala315Thr mutation was introduced to enable efficient trafficking of homomeric KCNQ3. KCNQ3 does not efficiently form functional channels; however, the Ala315Thr mutation enables efficient KCNQ3 functional expression without co-injection of KCNQ2 mRNA (throughout the text we refer to KCNQ3[Ala315Thr] as KCNQ3*)49. The plasmid encoding CivSP was kindly provided by Dr Y. Okamura. Complementary RNA was transcribed from the cDNA using the mMessage mMachine Kit (Ambion). Stage V–VI Xenopus laevis oocytes were prepared as previously described and were injected with cRNA alone or cRNA plus synthetic tRNA (for nonsense suppression, described below). Oocyte preparation was carried out using a protocol approved by the University of British Columbia Animal Care Committee, in accordance with the Canadian Council for Animal Care guidelines. We used female Xenopus laevis frogs 100 g or greater in size. After injection, oocytes were incubated for 12–48 h at 18 °C before recording.

Fluorinated Trp derivatives (5-F1-Trp; 7-F1-Trp; 5,7-F2-Trp (F2-Trp) and 5,6,7-F3-Trp (F3-Trp) were purchased from Asis Chem (Watertown, MA) and Sigma-Aldrich (St Louis, MO). The ‘Ind’ Trp variant was synthesized in-house as described previously31,50. Incorporation of unnatural amino acids into ion channel proteins was carried out as described extensively in prior publications29,34,51. Briefly, unnatural amino acids (aa) were protected with nitroveratryloxycarbonyl and were activated as the cyanomethyl ester, which were then coupled to the dinucleotide pdCpA (GE Healthcare/Dharmacon, Lafayette, CO). This aminoacyl dinucleotide was subsequently ligated to a modified (G73) Tetrahymena thermophile tRNA. The amino-acylated tRNA-aa was deprotected by ultraviolet irradiation immediately before oocyte injection. In a typical experiment, 80 ng of tRNA-aa and 40 ng of KCNQ3 cRNA were injected in a 50-nl volume. In control experiments, the cRNA alone or the cRNA together with a tRNA coupled to pdCpA but without an appended amino acid were injected (these parallel control experiments were conducted every experimental day and the results are described throughout the text).

Two-electrode voltage clamp recordings

Voltage-clamped potassium currents were recorded in standard Ringers solution (in mM): 116 NaCl, 2 KCl, 1 MgCl2, 0.5 CaCl2, 5 HEPES (pH 7.4) using an OC-725C voltage clamp (Warner, Hamden, CT). Glass microelectrodes were backfilled with 3 M KCl and had resistances of 0.1–1 MΩ. Data were filtered at 5 kHz and digitized at 10 kHz using a Digidata 1440A (Molecular Devices) controlled by the pClamp 10 software (Molecular Devices). Drugs were purchased from Toronto Research Chemicals (retigabine) or Tocris (ML-213, flupirtine, ICA-069673, ICA-110381), stored as 100 mM stocks in dimethylsulphoxide and diluted to working concentrations each experimental day.

Molecular simulations of retigabine binding

A published model of retigabine docked to the KCNQ3 pore domain was used as the starting point of simulations19. In this model the retigabine molecule is positioned with its potential H-bond donor carbonyl oxygen pointed away from Trp265. Therefore, the retigabine was rotated by 180° (‘flipped’ upside down). Three further retigabine molecules were docked into a similar position resulting in a model with all four binding sites occupied. This model was energy-minimized using force field AMBER03 and the constraints for the retigabine molecules were assigned with semi-empirical quantum mechanics (QM) using force field NOVA. Subsequently, a short molecular dynamics (MD) simulation was run while constraining the distance between the Trp265 N–H group and the carbonyl oxygen of retigabine to a distance of 4±1.5 Å. The resulting model contained four retigabine molecules in the ‘flip’ position, and all four Trp265 side chains in a rotated position (see Fig. 5c). Starting with this model, zero, one, two (in opposing or adjacent subunits), three or four retigabines were deleted from the respective model generating all possible combinations of occupied binding-site configurations. Each resulting model was inserted into a virtual membrane and MD simulations were run as described previously52. In brief, models were incorporated into membranes, the simulation box was filled with 0.9 mM NaCl/H2O and, subsequently, all-atom-mobile simulations were performed for 10 ns to reach stable conformations using YASARA Structure version 10.1 (AMBER03, 2 fs time steps; YASARA Biosciences GmbH, Vienna, Austria). The average structures of all models were calculated and the distance between the Trp265 N-H group and the carbonyl oxygen of retigabine was measured individually in each occupied binding site.

Data analysis

Voltage dependence of channel activation was fitted with a standard single component Boltzmann equation of the form G/Gmax=1/(1+ e ( V _ V 1 / 2 )/k where V1/2 is the voltage where channels exhibit half-maximal activation, and k is a slope factor reflecting the voltage range over which an e-fold change in Po is observed. For generation of concentration–response curves to retigabine and derivatives, ΔV1/2 (the shift in half-activation voltage) was normalized to the maximal ΔV1/2 observed, and fit with the following equation: normalized ΔV1/2=1/(1+EC50/[drug]). Statistical tests and significance are described in figure legends throughout the text.

Additional information

How to cite this article: Kim, R. Y. et al. Atomic basis for therapeutic activation of neuronal potassium channels. Nat. Commun. 6:8116 doi: 10.1038/ncomms9116 (2015).