Literature DB >> 31032878

Histidine at position 462 determines the low quinine sensitivity of ether-à-go-go channel superfamily member Kv 12.1.

Marlen Dierich1, Willem B van Ham2,3, Anna Stary-Weinzinger3, Michael G Leitner1,4.   

Abstract

BACKGROUND AND
PURPOSE: The ether-à-go-go (Eag) Kv superfamily comprises closely related Kv 10, Kv 11, and Kv 12 subunits. Kv 11.1 (termed hERG in humans) gained much attention, as drug-induced inhibition of these channels is a frequent cause of sudden death in humans. The exclusive drug sensitivity of Kv 11.1 can be explained by central drug-binding pockets that are absent in most other channels. Currently, it is unknown whether Kv 12 channels are equipped with an analogous drug-binding pocket and whether drug-binding properties are conserved in all Eag superfamily members. EXPERIMENTAL APPROACH: We analysed sensitivity of recombinant Kv 12.1 channels to quinine, a substituted quinoline that blocks Kv 10.1 and Kv 11.1 at low micromolar concentrations. KEY
RESULTS: Quinine inhibited Kv 12.1, but its affinity was 10-fold lower than for Kv 11.1. Contrary to Kv 11.1, quinine inhibited Kv 12.1 in a largely voltage-independent manner and induced channel opening at more depolarised potentials. Low sensitivity of Kv 12.1 and characteristics of quinine-dependent inhibition were determined by histidine 462, as site-directed mutagenesis of this residue into the homologous tyrosine of Kv 11.1 conferred Kv 11.1-like quinine block to Kv 12.1(H462Y). Molecular modelling demonstrated that the low affinity of Kv 12.1 was determined by only weak interactions of residues in the central cavity with quinine. In contrast, more favourable interactions can explain the higher quinine sensitivity of Kv 12.1(H462Y) and Kv 11.1 channels. CONCLUSIONS AND IMPLICATIONS: The quinoline-binding "motif" is not conserved within the Eag superfamily, although the overall architecture of these channels is apparently similar. Our findings highlight functional and pharmacological diversity in this group of evolutionary-conserved channels.
© 2019 The Authors. British Journal of Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.

Entities:  

Mesh:

Substances:

Year:  2019        PMID: 31032878      PMCID: PMC6609544          DOI: 10.1111/bph.14693

Source DB:  PubMed          Journal:  Br J Pharmacol        ISSN: 0007-1188            Impact factor:   8.739


ether‐à‐go‐go channel ether‐à‐go‐go‐gene‐like channel human ether‐à‐go‐go‐related gene channel 1 ether‐à‐go‐go‐related gene channel long QT molecular dynamics root‐mean‐square deviation

What is already known

Drug‐induced inhibition of Kv11.1 (hERG) is a frequent cause of sudden death in humans. It is not known whether closely related Kv12 channels are also sensitive to drug‐dependent inhibition.

What this study adds

Low quinine sensitivity, affinity, and characteristics of inhibition of Kv12.1 channels are determined by H462. Drug sensitivity of Kv11.1 and Kv12.1 is different but determined by homologous amino acid positions.

What is the clinical significance

Drug‐binding pockets are not conserved in Eag superfamily members, but channel architecture is similar. Our findings facilitate the understanding of the arrhythmogenic actions of quinine.

INTRODUCTION

The ether‐à‐go‐go (Eag) superfamily of Kv channels comprises the three conserved families of Eag (Kv10), ether‐à‐go‐go‐related gene (Erg; Kv11), and ether‐à‐go‐go‐gene‐like (Elk; Kv12) channels. These channels give rise to voltage‐dependent K+ currents in many cell types (Bauer & Schwarz, 2001; Bauer & Schwarz, 2018). The best characterised member, K (termed hERG for the human isoform), engenders rapidly activating K+ current I Kr responsible for membrane repolarisation in cardiac myocytes (Sanguinetti, Jiang, Curran, & Keating, 1995). As I Kr determines duration and end of heart action potentials, loss of Kv11.1 function constitutes a frequent cause of cardiac dysfunction in humans (Curran et al., 1995; Keating & Sanguinetti, 2001; Mitcheson, Chen, Lin, Culberson, & Sanguinetti, 2000; Sanguinetti et al., 1995; Trudeau, Warmke, Ganetzky, & Robertson, 1995). Mutations in KCNH2, the gene encoding Kv11.1, cause congenital long QT (LQT) syndrome‐2 characterised by a prolonged QT interval and polymorphic ventricular arrhythmias (torsade de pointes) that may lead to recurrent syncope or sudden death (Curran et al., 1995; Sanguinetti et al., 1995). More common, however, are acquired forms of LQT syndrome through drug‐induced inhibition of Kv11.1 (Keating & Sanguinetti, 2001; Roden, 1996; Sanguinetti et al., 1995). Kv11.1 channels are extremely sensitive to a wide variety of drugs including substituted quinolines (quinidine, quinine, and chloroquine), antiarrhythmic agents (e.g., MK‐499 and dofetilide), and many other substances (e.g., terfenadine, cisapride, and vesnarinone; Furutani et al., 2011; Kamiya, Mitcheson, Yasui, Kodama, & Sanguinetti, 2001; Lees‐Miller, Duan, Teng, & Duff, 2000; Mitcheson et al., 2000; Mitcheson et al., 2005; Sanchez‐Chapula, Ferrer, Navarro‐Polanco, & Sanguinetti, 2003; Sanchez‐Chapula, Navarro‐Polanco, Culberson, Chen, & Sanguinetti, 2002). All of these substances may cause undesired prolongation of cardiac action potentials by inhibiting Kv11.1 (Sanguinetti et al., 1995). The high susceptibility of human Kv11.1 to these structurally divergent drugs is governed by amino acids facing the channel's central cavity to form hydrophobic‐binding pockets (Wang & MacKinnon, 2017). Mutagenesis studies have shown that amino acids in the inner pore helix (threonine 623), selectivity filter (serine 624 and valine 625), and tyrosine 652/phenylalanine 656 of the sixth transmembrane segment (S6) determine the susceptibility of drug‐induced inhibition in Kv11.1 (Mitcheson et al., 2000; Sanchez‐Chapula et al., 2002; Sanchez‐Chapula et al., 2003; Wang & MacKinnon, 2017). Some of these residues are involved in drug binding (threonine 623, serine 624, tyrosine 652, and phenylalanine 656; Wang & MacKinnon, 2017), whereas others (e.g., valine 625 and also asparagine 588, serine 631, and serine 620) have been shown to control drug sensitivity through contributing to channel inactivation that is necessary for Kv11.1 inhibition through many high‐affinity blockers (Ficker, Jarolimek, & Brown, 2001; Ficker, Jarolimek, Kiehn, Baumann, & Brown, 1998; Kamiya, Niwa, Mitcheson, & Sanguinetti, 2006; Perrin, Kuchel, Campbell, & Vandenberg, 2008; Wu, Gardner, & Sanguinetti, 2015). Further, phenylalanine 557 located in the S5 helix is involved in binding of some drugs (Helliwell et al., 2018; Saxena et al., 2016). Most of these residues are conserved in Kv10 channels, and accordingly, the K isoform exhibits similarly high sensitivity to some substances (Schonherr, Gessner, Lober, & Heinemann, 2002). However, K is significantly less susceptible to drug‐induced inhibition, indicating that additional motifs might contribute to drug binding in certain members of the superfamily (Chen, Seebohm, & Sanguinetti, 2002; Gessner, Zacharias, Bechstedt, Schonherr, & Heinemann, 2004; Schonherr et al., 2002). In fact, it has been shown that high‐affinity blockers of Kv11.1 are less effective on Kv10 channels, as these isoforms do not inactivate. This suggests that conformational reorientation of relevant residues associated with inactivation that are necessary for inhibition of KV11.1 is absent in Kv10 channels (Chen et al., 2002; Ficker et al., 1998; Ficker et al., 2001). In contrast to Kv10 and Kv11, only little information is available on the three members of the Kv12 (Elk) family predominantly expressed in neurons (Engeland, Neu, Ludwig, Roeper, & Pongs, 1998; Miyake, Mochizuki, Yokoi, Kohda, & Furuichi, 1999; Saganich, Machado, & Rudy, 2001; Shi et al., 1998; Trudeau, Titus, Branchaw, Ganetzky, & Robertson, 1999; Zou et al., 2003). K regulates excitability of hippocampal neurons (Zhang et al., 2010), but no physiological role has been assigned to K and K channels yet. Recent studies provided significant insight into functional characteristics of Kv12.1 (Dai & Zagotta, 2017; Dierich, Evers, Wilke, & Leitner, 2018; Dierich & Leitner, 2018; Kazmierczak et al., 2013; Li et al., 2015), but it is unknown whether Kv12 channels are equipped with a high affinity drug‐binding pocket as Kv10.1 and Kv11.1. We analysed the sensitivity of human Kv12.1 to quinine, a substituted quinoline that blocks Kv11.1 channels at low micromolar concentrations (Gessner et al., 2004; Mitcheson et al., 2000; Sanchez‐Chapula et al., 2003; Schonherr et al., 2002). Quinine also inhibited recombinant Kv12.1 channels, but their sensitivity was 10‐fold lower than that of Kv11.1. Contrary to Kv11.1, quinine‐dependent block of Kv12.1 was largely voltage‐independent and the substance induced channel opening at more depolarised membrane potentials. Mutagenesis of histidine at position 462 (H462) into the homologous tyrosine (Y652) of Kv11.1 conferred Kv11.1‐like quinine sensitivity to Kv12.1(H462Y) channels. Molecular modelling analyses demonstrated that the low quinine sensitivity of Kv12.1 was determined by only weak interactions of quinine with H462 in the central cavity, whereas more favourable interactions facilitated higher quinine sensitivity of Kv12.1(H462Y) and closely related Kv11.1. Our data showed that the drug‐binding pocket in the central cavity is not conserved in the Eag superfamily, although the general architecture of these channels is apparently similar. Our findings highlight functional and pharmacological diversity within this group of evolutionary conserved K+ channels.

METHODS

Cell culture, transient transfection, and mutagenesis

CHO dhFR− cells (ATCC Cat# CRL‐9096, RRID:CVCL_1977) were maintained as previously reported (Leitner et al., 2016). In brief, cells were maintained in MEM alpha medium (with 10% fetal calf serum and 1% pen/strep; all Invitrogen GmbH, Darmstadt, Germany) at 5% CO2 and 37°C in a humidified atmosphere. Transient transfection of cultured CHO cells was performed using jetPEI (Polyplus Transfection, Illkirch, France). The following vectors for ion channel expression were used: Kv11.1 (Erg1)–pcDNA3.1 (gene: rat Kcnh2; UniProt accession number: O08962; UniProt, RRID:SCR_002380), Kv12.1 (Elk1)–pcDNA3.1–IRES–eGFP (human KCNH8; Q96L42), and pEGFP–C1 (transfection control; Addgene, Teddington, UK). An amino acid exchange (H462Y) was introduced into Kv12.1 with the QuikChangeII XL site‐directed mutagenesis kit (Stratagene, Santa Clara, CA). Site‐directed mutagenesis was confirmed by sequencing prior to the experiments (Microsynth SEQLAB, Göttingen, Germany).

Electrophysiological recordings

Whole‐cell patch clamp recordings were performed at room temperature (22–25°C) with an HEKA EPC10 USB patch clamp amplifier controlled by PatchMaster software (HEKA, Lambrecht, Germany; Patchmaster, RRID:SCR_000034). Voltage clamp recordings were low‐pass filtered at 2 kHz and sampled at 5 kHz. The series resistance (R ) was kept below 6 MΩ, and R was compensated throughout the recordings (80–90%; Dierich & Leitner, 2018; Leitner, Halaszovich, & Oliver, 2011). All voltage protocols are indicated in the figures; dashed lines highlight zero current. Borosilicate glass patch pipettes (Sutter Instrument Company, Novato, CA) had an open‐pipette resistance of 2–3 MΩ after backfilling with intracellular solution containing (in mM) 135 KCl, 2.41 CaCl2 (100‐nM free Ca2+), 3.5 MgCl2, 5 HEPES, 5 EGTA, 2.5 Na2ATP, and 0.1 Na3GTP, pH 7.3 (with KOH), 290–295 mOsm·kg−1. The extracellular solution contained (in mM) 144 NaCl, 5.8 KCl, 1.3 CaCl2, 0.9 MgCl2, 0.7 NaH2PO4, 10 HEPES, and 5.6 d‐glucose, pH 7.4 (with NaOH), 305–310 mOsm·kg−1. Liquid junction potentials were not compensated (approximately −4 mV).

Analysis of electrophysiological recordings

Patch clamp recordings were analysed with IgorPro (Wavemetrics, Lake Oswego, OR; IGOR Pro, RRID:SCR_000325) and the PatchMaster (HEKA) software. Voltage‐dependence of activation was derived from tail current amplitudes using voltage protocols indicated: Tail currents were fitted with a two‐state Boltzmann function with I = I min + (I max − I min)∕(1 + exp((V − V )∕s)), where I is the current, V is the membrane voltage, V is the voltage at half‐maximal activation, and s describes the slope of the curve (s is presented as positive values to describe the slope of voltage‐dependent channel activation; Leitner et al., 2012). Results are shown as conductance–voltage curves, obtained by normalising to (I max − I min), obtained from fits to data of individual experiments. For dose–response relationships, currents were normalised to baseline and were fitted to a Hill equation with , where I is the (normalised) current, I and I max denote minimal and maximal currents, respectively, IC50 is the concentration at the half‐maximal effect, [S] is the drug concentration, and n H is the Hill coefficient. Time constants of activation were derived from mono‐exponential fits to activating current components at indicated potentials (τ activation; c.f. Figure 3e). For presentation, whole‐cell currents were normalised to cell capacitance (current density; pA·pF−1) or to baseline current amplitudes (I/I 0; normalised current).
Figure 3

Characteristics of the Kv12.1(H462Y) channel mutant. (a) The panel shows an alignment of Kv11.1 isoforms from humans and rats (encoded by the KCNH2 gene) with human Kv12.1 channels (KCNH8 gene). Amino acids recently identified as important for quinoline‐dependent inhibition of Kv11.1 channels are highlighted in colour (green, identical amino acids; red, amino acid exchange in Kv12.1; see text for details and references). Note that the quinine‐binding “motif” of Kv12.1 differs only at one position from Kv11.1 channels. (b–h) Biophysical characteristics of Kv12.1(H462Y; red) in comparison with Kv12.1 wild‐type channels (black). (b, c) Representative recordings of currents through Kv12.1 wild‐type (black) and Kv12.1(H462Y; red) channels activated with the shown voltage protocols. (d) Steady‐state outward currents through Kv12.1(H462Y) were significantly smaller than through Kv12.1 wild‐type channels, and (e) Kv12.1(H462Y) activated significantly faster than wild‐type channels, when activated through voltage steps between −40 and 0 mV (τ derived from monoexponential fits to activating current components in recordings shown in b, c). (f–h) Kv12.1(H462Y) activated at significantly more negative membrane potentials than wild‐type channels after conditioning potential of −60 mV. (f) Summary of voltage‐dependence of Kv12.1(H462Y) channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (g) Mean V of channel activation and (h) summarised slope factors of activation (data on Kv12.1 wild‐type channels are reproduced from Figure 1). *denotes significant differences

Data and statistical analysis

The data and statistical analysis comply with the recommendations of the British Journal of Pharmacology on experimental design and analysis in pharmacology (Curtis et al., 2018). Isolated cells under investigation were randomly assigned to different treatment groups. Data analysis for experiments presented was performed in a blinded manner. For some experiments, single recordings were normalised to baseline values individually to account for baseline variations between cells. Statistical analysis was performed using Student's two‐tailed t test/Wilcoxon–Mann–Whitney test, and when appropriate comparisons between multiple groups were performed with ANOVA followed by Dunnett's test. Post hoc tests were run only if F achieved P < 0.05 and there was no significant variance inhomogeneity. Significance was assigned at P ≤ 0.05 (*P ≤ 0.05). Data subjected to statistical analysis have n over 5 per group, and data are presented as mean ± standard error of the means (SEM). In electrophysiological experiments, n represents the number of individual cells and accordingly the number of independent experiments (no pseudo‐replication).

Molecular modelling

A Kv12.1 homology model of the pore module (residues S352–Y477) in the open/inactive conformation was built using the programme modeller 9v.17 (Webb & Sali, 2016, RRID:SCR_008395) based on the cryo‐EM structure of the hERG (human Kv11.1) channel (pdb identifier: 5VA1, 3.7 Å resolution; Wang & MacKinnon, 2017). The sequence identity between target and template is 61% (sequence similarity: 76%). Docking studies were carried out with the Genetic Optimization for Ligand Docking program, version 5.6.2 (GOLD, Jones, Willett, Glen, Leach, & Taylor, 1997). The “chemscore” fitness scoring function was used, and all residues within 12 Å of residue H462 in Kv12.1 were defined as binding site. The Chemscore.DG scoring function was used to estimate free energies of binding. Side chains of residues H462 and F466 from all four domains were set as flexible. The pKa of H462 is calculated to be 3.13 (PROPKA; Olsson, Sondergaard, Rostkowski, & Jensen, 2011), indicating a hydrogen atom on the δ nitrogen. The other two tautomers were considered in docking as well; however, similar binding modes/poses were obtained. For molecular dynamics (MD) system setups and ligand parametrisation (Vanommeslaeghe et al., 2010), the CHARMM‐GUI (Jo, Kim, Iyer, & Im, 2008; CHARMM, RRID:SCR_014892) was used. The quinine molecule was protonated at the tertiary N (pKa value of the quinolone group: 9.7, O'Neil, 2013). Protein–ligand complexes obtained from docking were embedded in a 1‐palmitoyl‐2‐oleoyl‐sn‐glycero‐3‐phosphocholine bilayer and solvated with TIP3P waters. K+ ions were placed in the selectivity filter at sites S0, S2, and S4, with water molecules at sites S1 and S3 and 0.15‐M KCl added to the simulation box. Energy minimisation, 20‐ns equilibration, and three independent 50‐ns production runs (with initial atom velocities assigned independently and randomly) were performed using GROMACS v.5.1.2 (GROMACS, RRID:SCR_014565; Abraham et al., 2016) with the charmm36 force field for protein molecules (incorporating CMAP terms, Mackerell, Feig, & Brooks, 2004) and lipid molecules and salt ions (Best et al., 2012; Klauda et al., 2010; MacKerell et al., 1998). Electrostatics were modelled using particle mesh Ewald (Darden, York, & Pedersen, 1993), and LINCS was used to constrain covalent chemical bonds to hydrogens (Hess, Bekker, Berendsen, & Fraaije, 1997). Temperature was maintained at 310 K using velocity rescaling (V‐rescale; Bussi, Donadio, & Parrinello, 2007), and semi‐isotropic pressure coupling was accomplished using the Parrinello–Rahman barostat (Parrinello & Rahman, 1981). MD trajectories were analysed using VMD v.1.9.2 (VMD, RRID:SCR_001820; Humphrey, Dalke, & Schulten, 1996) and GROMACS.

Materials

Quinine was purchased from Tocris Bioscience (Bristol, UK) and was diluted in extracellular solution to concentrations indicated in Section 3. Quinine was applied locally via a glass capillary through a custom‐made application system.

Nomenclature of targets and ligands

Key protein targets and ligands in this article are hyperlinked to corresponding entries in http://www.guidetopharmacology.org, the common portal for data from the IUPHAR/BPS Guide to PHARMACOLOGY (Harding et al., 2018), and are permanently archived in the Concise Guide to PHARMACOLOGY 2017/18 (Alexander et al., 2017).

RESULTS

Voltage‐independent but mode shift‐dependent inhibition of Kv12.1 channels through quinine at high concentrations

We started our pharmacological analyses by applying increasing concentrations of quinine to CHO cells transiently expressing human Kv12.1. When activating Kv12.1 with 600‐ms voltage steps to 0 mV from a holding potential of −60 mV, quinine inhibited Kv12.1 channels with an IC50 of 0.97 mM and a Hill coefficient of 1.0 (Figure 1a). At a concentration close to the IC50 (1 mM), this inhibition developed within seconds and was partially reversible after removal of the drug (Figure S1a,b). After having established the quinine sensitivity of Kv12.1, we set out to characterise effects of the substance on voltage‐dependent activation of the channels. Noteworthy, Kv12.1 channels exhibit a mode shift of activation (also termed pre‐pulse facilitation or voltage‐dependent potentiation; Dai & Zagotta, 2017; Dierich et al., 2018; Dierich & Leitner, 2018; Li et al., 2015), which designates stabilisation of the voltage‐sensing domain in a “relaxed” open state after prolonged depolarisation of the membrane potential (Bezanilla, Taylor, & Fernandez, 1982; Villalba‐Galea, Sandtner, Starace, & Bezanilla, 2008). Most prominently, mode shift manifests through a large shift of activation voltages to hyperpolarised potentials following membrane potential depolarisation (Dai & Zagotta, 2017; Dierich et al., 2018; Dierich & Leitner, 2018; Li et al., 2015). To induce mode shift, we applied conditioning potentials of −60 or 0 mV (200 ms) before a series of voltage steps to activate Kv12.1 (pulse potentials: −140 to +20 mV; 600 ms; Figure 1b,c; c.f. Dierich & Leitner, 2018; Dierich et al., 2018). Following the hyperpolarised conditioning potential (−60 mV), half‐maximal voltage (V ) and slope factor of Kv12.1 channel activation were −30.0 ± 2.0 mV and 14.9 ± 1.1 mV, respectively (n = 6; Figure 1d–f). When we applied the conditioning potential of 0 mV, V was −75.9 ± 1.7 mV and s was 10.6 ± 0.5 mV (n = 6; Figure 1d–f). Depolarised conditioning potentials thus induced a large (and significant) shift of the voltages of activation of Kv12.1 channels to hyperpolarised potentials by about −45 mV, representing their mode shift of activation (Dai & Zagotta, 2017; Dierich et al., 2018; Dierich & Leitner, 2018; Li et al., 2015). Application of quinine (1 mM) shifted the voltage dependence of Kv12.1 to depolarised voltages for both conditioning potentials. As this shift was more pronounced after the depolarised conditioning potential (0 mV; Figure 1d,e), the degree of mode shift was significantly attenuated to about −23 mV in the presence of the substance (Figure 1e). Application of quinine (1 mM) inhibited both steady‐state currents and tail currents through Kv12.1 channels (Figure 1g,h). Following both conditioning potentials, quinine‐induced inhibition of steady‐state Kv12.1 currents was the same for all activating voltages (Figure 1i; n = 6; e.g., for conditioning potential −60 mV: not significant (ns) between −40 mV and +20 mV, P = 0.54), that is, quinine inhibited Kv12.1 channels in a voltage‐independent manner. However, quinine‐dependent block was significantly more pronounced following the hyperpolarised conditioning potential (−60 mV) for activating voltages more positive than −30 mV (indicated by “*” in Figure 1i). These data demonstrated that to some extent, Kv12.1 channels were more sensitive to quinine block after conditioning depolarisation of the membrane, and that accordingly, quinine block of Kv12.1 channels depended on mode shift.
Figure 1

Quinine inhibits human Kv12.1 channels at high concentrations in a voltage‐independent manner. (a) Dose‐dependent inhibition of recombinant Kv12.1 channels through extracellular application of quinine. The panel shows representative recordings of a CHO cell expressing human Kv12.1 treated with increasing quinine concentrations (left) and the summarised quinine dose–response relationship (right; IC50 and Hill coefficient were calculated from a Hill fit to averaged recordings as shown in the left panel; solid line represents this fit). Scale bars, voltage protocol, quinine concentrations, and number of cells as indicated. (b–f) Quinine induced the activation of recombinant Kv12.1 channels at more positive membrane potentials. Representative recordings of Kv12.1 channels activated by voltage steps between −140 and +20 mV following 200‐ms conditioning voltage steps to (b) −60 mV or (c) 0 mV before (control) and at the end of quinine application (1 mM). (d) Summary of voltage‐dependence of Kv12.1 channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (e) Mean V of channel activation and (f) summarised s factors of activation before and at the end of quinine treatment (1 mM; values derived from fits as shown in d). Quinine (1 mM) inhibited voltage‐dependent (g) outward steady‐state and (h) inward tail currents through recombinant Kv12.1 channels (current amplitudes analysed from recordings as shown in b, c). (i) Quinine‐induced inhibition of Kv12.1 was voltage‐independent (presented as % inhibition from data shown in (g); ns, no significant difference for quinine block between activating voltages of −40 and 20 mV). Note that quinine block was significantly more pronounced, when Kv12.1 channels were activated at membrane potentials more positive than −30 mV. *Significant difference of quinine block for both conditioning potentials

Quinine inhibits human Kv12.1 channels at high concentrations in a voltage‐independent manner. (a) Dose‐dependent inhibition of recombinant Kv12.1 channels through extracellular application of quinine. The panel shows representative recordings of a CHO cell expressing human Kv12.1 treated with increasing quinine concentrations (left) and the summarised quinine dose–response relationship (right; IC50 and Hill coefficient were calculated from a Hill fit to averaged recordings as shown in the left panel; solid line represents this fit). Scale bars, voltage protocol, quinine concentrations, and number of cells as indicated. (b–f) Quinine induced the activation of recombinant Kv12.1 channels at more positive membrane potentials. Representative recordings of Kv12.1 channels activated by voltage steps between −140 and +20 mV following 200‐ms conditioning voltage steps to (b) −60 mV or (c) 0 mV before (control) and at the end of quinine application (1 mM). (d) Summary of voltage‐dependence of Kv12.1 channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (e) Mean V of channel activation and (f) summarised s factors of activation before and at the end of quinine treatment (1 mM; values derived from fits as shown in d). Quinine (1 mM) inhibited voltage‐dependent (g) outward steady‐state and (h) inward tail currents through recombinant Kv12.1 channels (current amplitudes analysed from recordings as shown in b, c). (i) Quinine‐induced inhibition of Kv12.1 was voltage‐independent (presented as % inhibition from data shown in (g); ns, no significant difference for quinine block between activating voltages of −40 and 20 mV). Note that quinine block was significantly more pronounced, when Kv12.1 channels were activated at membrane potentials more positive than −30 mV. *Significant difference of quinine block for both conditioning potentials

Quinine is a voltage‐ and mode shift‐dependent inhibitor of Kv11.1 channels

We then compared our findings on Kv12.1 to quinine‐dependent inhibition of closely related Kv11.1. In line with a previous study (Sanchez‐Chapula et al., 2003), quinine inhibited recombinant Kv11.1 channels with an IC50 of 98 μM and a Hill coefficient of 1.3 (Figure 2a). Accordingly, the quinine affinity of Kv11.1 channels was about 10 times higher than that of Kv12.1 (c.f. Figure 1a). As for Kv12.1, quinine‐dependent inhibition of Kv11.1 developed in seconds and was partially reversible (50‐μM quinine; Figure S1c,d). By applying hyperpolarised (−60 mV) and depolarised (+40 mV) conditioning voltages before the activating pulse potentials, we then analysed effects of quinine on the voltage dependence of Kv11.1 channels that also exhibit mode shift of activation (Figure 2b,c). As previously reported (Dierich et al., 2018; Tan, Perry, Ng, Vandenberg, & Hill, 2012), depolarised conditioning potentials (+40 mV) significantly shifted the voltage dependence of Kv11.1 channels by about −50 mV compared with the hyperpolarised conditioning potential (Figure 2d–f). Following the hyperpolarised conditioning potential (−60 mV), quinine (50 μM) significantly shifted the voltage dependence of Kv11.1 channels by −13 mV to more negative potentials (Figure 2d–f). As the voltage dependence of Kv11.1 was not altered following the depolarised conditioning potential (+40 mV), the degree of mode shift was attenuated in the presence of quinine (50 μM; Figure 2d,e). Quinine blocked both voltage‐dependent outward and inward currents through Kv11.1 (Figure 2g,h). After the conditioning potential of −60 mV, the degree of this inhibition was significantly increased, when the channels were activated at more positive potentials (Figure 2i), that is, quinine block was voltage dependent following the hyperpolarised conditioning potential. In contrast, following the depolarised conditioning potential, quinine block of Kv11.1 channels was completely voltage independent. Thus, for Kv11.1 channels, quinine‐induced inhibition depended on mode shift (c.f. quinidine; Furutani et al., 2011), but (in contrast to Kv12.1) the Kv11.1 channels apparently were more sensitive to quinine following the depolarised conditioning potential, when activated at negative membrane potentials.
Figure 2

Voltage‐dependent inhibition of Kv11.1 channels induced by quinine. (a) Quinine inhibited Kv11.1 channels transiently expressed in CHO cells in a concentration‐dependent manner. The left panel shows representative recordings of a CHO cell expressing Kv11.1 treated with increasing quinine concentrations, the right panel depicts the summarised quinine dose–response relationship (IC50 and Hill coefficient calculated from a Hill fit to averaged data; solid line represents Hill fit). Scale bars, voltage protocol, quinine concentrations, and number of cells recorded as indicated. (b–f) After a conditioning potential of −60 mV, quinine induced activation of recombinant Kv11.1 channels at more negative membrane potentials. Representative recordings of Kv11.1 channels activated with voltage steps (as indicated) after conditioning voltage steps to (b) −60 mV or (c) +40 mV before (control) and at the end of the application of quinine (50 μM). (d) Summary of voltage‐dependence of Kv11.1 channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (e) Mean V of channel activation and (f) summarised slope factors of activation before and at the end of quinine application (50 μM; values derived from fits as shown in d). Quinine (1 mM) inhibited voltage‐dependent (g) outward steady‐state and (h) inward tail currents through recombinant Kv11.1 channels (current amplitudes analysed from recordings as shown in b, c). (i) Quinine‐induced inhibition of Kv11.1 was voltage‐dependent after conditioning potential of −60 mV but not after the depolarised conditioning potential of +40 mV (presented as % inhibition from data shown in g). *denotes significant differences

Voltage‐dependent inhibition of Kv11.1 channels induced by quinine. (a) Quinine inhibited Kv11.1 channels transiently expressed in CHO cells in a concentration‐dependent manner. The left panel shows representative recordings of a CHO cell expressing Kv11.1 treated with increasing quinine concentrations, the right panel depicts the summarised quinine dose–response relationship (IC50 and Hill coefficient calculated from a Hill fit to averaged data; solid line represents Hill fit). Scale bars, voltage protocol, quinine concentrations, and number of cells recorded as indicated. (b–f) After a conditioning potential of −60 mV, quinine induced activation of recombinant Kv11.1 channels at more negative membrane potentials. Representative recordings of Kv11.1 channels activated with voltage steps (as indicated) after conditioning voltage steps to (b) −60 mV or (c) +40 mV before (control) and at the end of the application of quinine (50 μM). (d) Summary of voltage‐dependence of Kv11.1 channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (e) Mean V of channel activation and (f) summarised slope factors of activation before and at the end of quinine application (50 μM; values derived from fits as shown in d). Quinine (1 mM) inhibited voltage‐dependent (g) outward steady‐state and (h) inward tail currents through recombinant Kv11.1 channels (current amplitudes analysed from recordings as shown in b, c). (i) Quinine‐induced inhibition of Kv11.1 was voltage‐dependent after conditioning potential of −60 mV but not after the depolarised conditioning potential of +40 mV (presented as % inhibition from data shown in g). *denotes significant differences

Functional characterisation of Kv12.1(H462Y) channels

Sensitivity of human Kv11.1 channels to several drugs has been attributed to the amino acids T623, V625, Y652, and F656, located in the channels inner pore helix, selectivity filter, or sixth transmembrane segment (S6; Figure 3a; Lees‐Miller et al., 2000; Mitcheson et al., 2000; Sanchez‐Chapula et al., 2002; Sanchez‐Chapula et al., 2003). As Kv12.1 also carries threonine, valine, and phenylalanine at homologous positions (T433, V435, and F466 in Kv12.1), we hypothesised that the amino acid exchange at position 462, where Kv12.1 has a histidine instead of the tyrosine (Figure 3a), determined the low quinine sensitivity of Kv12.1. To elaborate this hypothesis, we mutated H462 in Kv12.1 to the corresponding tyrosine (Y652) of Kv11.1 and characterised the Kv12.1(H462Y) mutant in CHO cells. Recombinant Kv12.1(H462Y) channels produced robust voltage‐dependent and outwardly rectifying K+ currents (Figure 3b,c). However, Kv12.1(H462Y)‐mediated whole‐cell currents were significantly smaller than for the wild‐type channels (Figure 3b–d), and the mutant channels activated significantly faster than the wild‐type, when activated through voltage steps between −40 and 0 mV (Figure 3e). For conditioning voltages of −60 and 0 mV, V was −38.0 ± 1.8 mV and −73.6 ± 1.6 mV, respectively, demonstrating that Kv12.1(H462Y) channels also exhibited a mode shift of activation (Figure 3f–h; n = 9). However, in the mutant, mode shift was attenuated compared with wild‐type to about −35 mV (−45 mV in wild‐type; c.f. Figure 1), mainly because Kv12.1(H462Y) activated at significantly more negative membrane potentials after the hyperpolarised conditioning potential of −60 mV (Figure 3f–h). The voltage dependence of wild‐type and mutant channels was the same after the conditioning potential of 0 mV (Figure 3f–h). Characteristics of the Kv12.1(H462Y) channel mutant. (a) The panel shows an alignment of Kv11.1 isoforms from humans and rats (encoded by the KCNH2 gene) with human Kv12.1 channels (KCNH8 gene). Amino acids recently identified as important for quinoline‐dependent inhibition of Kv11.1 channels are highlighted in colour (green, identical amino acids; red, amino acid exchange in Kv12.1; see text for details and references). Note that the quinine‐binding “motif” of Kv12.1 differs only at one position from Kv11.1 channels. (b–h) Biophysical characteristics of Kv12.1(H462Y; red) in comparison with Kv12.1 wild‐type channels (black). (b, c) Representative recordings of currents through Kv12.1 wild‐type (black) and Kv12.1(H462Y; red) channels activated with the shown voltage protocols. (d) Steady‐state outward currents through Kv12.1(H462Y) were significantly smaller than through Kv12.1 wild‐type channels, and (e) Kv12.1(H462Y) activated significantly faster than wild‐type channels, when activated through voltage steps between −40 and 0 mV (τ derived from monoexponential fits to activating current components in recordings shown in b, c). (f–h) Kv12.1(H462Y) activated at significantly more negative membrane potentials than wild‐type channels after conditioning potential of −60 mV. (f) Summary of voltage‐dependence of Kv12.1(H462Y) channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (g) Mean V of channel activation and (h) summarised slope factors of activation (data on Kv12.1 wild‐type channels are reproduced from Figure 1). *denotes significant differences

Histidine at position 462 determines the quinine sensitivity of human Kv12.1

When applying increasing quinine concentrations to CHO cells expressing Kv12.1(H462Y), we found that quinine inhibited this channel mutant with an IC50 of 56 μM and a Hill coefficient of 1.1 (Figure 4a). Thus, the quinine affinity of Kv12.1(H462Y) was about 17‐fold higher than that of the wild‐type channels and accordingly more similar to Kv11.1 (c.f. Figures 1a and 2a). As for Kv12.1 wild‐type, application of quinine (50 μM) blocked inward and outward currents through Kv12.1(H462Y) (Figure 4b–e), but in contrast to the wild‐type channels, quinine‐mediated inhibition of Kv12.1(H462Y) was rapidly reversible (Figure S1e,f). Noteworthy, we detected a transient increase of currents in the presence of quinine, when activating Kv12.1(H462Y) at membrane voltages more positive than 0 mV (following the hyperpolarised conditioning potential; Figure 4b). This observation indicated that channel opening may be required for quinine‐dependent inhibition of Kv12.1(H462Y). In line, quinine‐induced inhibition of Kv12.1(H462Y) was more pronounced, when the channels where activated at more positive membrane voltages, that is, quinine block of Kv12.1(H462Y) was voltage dependent following positive and negative conditioning potentials (Figure 4f). For Kv12.1(H462Y) channels, this voltage dependence of quinine block presumably also caused a reduction of tail currents at depolarised activating voltages after conditioning potential of 0 mV (Figure 4e,g). Of special note, in contrast to Kv12.1 wild‐type channels, quinine block of Kv12.1(H462Y) was the same following both conditioning potentials and thus was independent on mode shift. Contrary to Kv12.1 wild‐type, quinine (50 μM) shifted the voltage dependence of Kv12.1(H462Y) by −11 mV and −10 mV to hyperpolarised potentials after conditioning voltages of −60 mV and 0 mV respectively (Figure 4g–i).
Figure 4

Voltage‐dependent inhibition of Kv12.1(H462Y) channels through quinine at low concentrations. (a) Extracellular application of quinine inhibited Kv12.1(H462Y) channels overexpressed in CHO cells. The panel shows representative recordings of a CHO cell expressing Kv12.1(H462Y) treated with increasing quinine concentrations (left) and the summarised dose–response relationship (right; IC50 and Hill coefficient were calculated from a Hill fit to averaged recordings as shown in the left panel; solid line represents Hill fit). Scale bars, voltage protocol, quinine concentrations, and number or cells recorded as indicated. (b–i) Quinine inhibited inward and outward currents through Kv12.1(H462Y) channels and induced channel activation at more negative membrane potentials. Representative recordings of recombinant Kv12.1 channels activated by voltage steps between −140 and +20 mV following 200‐ms conditioning voltage steps to (b) −60 mV or (c) 0 mV before (black) and at the end of quinine application (50 μM; red). Summarised voltage‐dependent (d) steady‐state currents and (e) inward tail currents through Kv12.1(H462Y) before (black) and at the end of quinine treatment (50 μM; current amplitudes analysed from recordings as shown in b, c). (f) Quinine‐induced inhibition of Kv12.1(H462Y) was voltage‐dependent (presented as % inhibition from data shown in d). (g–i) Summary of voltage‐dependence of Kv12.1(H462Y) channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (h) Mean V of channel activation and (i) summarised slope factors of activation before and at the end of treatment with quinine (50 μM; values derived from fits as shown in g). *denotes significant differences

Voltage‐dependent inhibition of Kv12.1(H462Y) channels through quinine at low concentrations. (a) Extracellular application of quinine inhibited Kv12.1(H462Y) channels overexpressed in CHO cells. The panel shows representative recordings of a CHO cell expressing Kv12.1(H462Y) treated with increasing quinine concentrations (left) and the summarised dose–response relationship (right; IC50 and Hill coefficient were calculated from a Hill fit to averaged recordings as shown in the left panel; solid line represents Hill fit). Scale bars, voltage protocol, quinine concentrations, and number or cells recorded as indicated. (b–i) Quinine inhibited inward and outward currents through Kv12.1(H462Y) channels and induced channel activation at more negative membrane potentials. Representative recordings of recombinant Kv12.1 channels activated by voltage steps between −140 and +20 mV following 200‐ms conditioning voltage steps to (b) −60 mV or (c) 0 mV before (black) and at the end of quinine application (50 μM; red). Summarised voltage‐dependent (d) steady‐state currents and (e) inward tail currents through Kv12.1(H462Y) before (black) and at the end of quinine treatment (50 μM; current amplitudes analysed from recordings as shown in b, c). (f) Quinine‐induced inhibition of Kv12.1(H462Y) was voltage‐dependent (presented as % inhibition from data shown in d). (g–i) Summary of voltage‐dependence of Kv12.1(H462Y) channels derived from Boltzmann fits to individual recordings as shown in (b, c); solid line represents Boltzmann fit to averaged data. (h) Mean V of channel activation and (i) summarised slope factors of activation before and at the end of treatment with quinine (50 μM; values derived from fits as shown in g). *denotes significant differences

Kv12.1(H462Y) channels exhibit Kv11.1‐like quinine sensitivity

In summary, Kv12.1(H462Y) exhibited a 17 times higher quinine affinity than Kv12.1 wild‐type, that is, quinine sensitivity of Kv12.1(H462Y) was more similar to Kv11.1 than to the wild‐type channels (Figure 5a). Application of quinine induced opening of wild‐type Kv12.1 channels at more positive potentials but at more negative values for Kv12.1(H462Y) and Kv11.1 (Figure 5b). Quinine inhibited Kv12.1 channels in a voltage‐independent manner, but quinine block was completely voltage dependent for Kv12.1(H462Y) and partially voltage dependent (after hyperpolarised conditioning potentials) for Kv11.1 channels (Figure 5c,d). In contrast to Kv12.1(H462Y), Kv12.1 wild‐type channels apparently were more sensitive to quinine after conditioning hyperpolarisation of the membrane (c.f. Figure 1i), that is, for Kv12.1 wild‐type channels, quinine‐induced inhibition depended on mode shift (just as for Kv11.1 channels). Accordingly, closely related Kv12.1 and Kv11.1 channels exhibited markedly different quinine sensitivity and characteristics of inhibition, but introduction of the H462Y amino acid exchange conferred Kv11.1‐like quinine sensitivity to Kv12.1 channels. Based on these findings, we conclude that H462 determines the low quinine sensitivity and the characteristics of quinine‐dependent inhibition for Kv12.1 channels and that the binding pocket for substituted quinolines in the central cavity is not completely conserved within the Eag superfamily of Kv channels. However, as the H462Y point mutation restored Kv11.1‐like quinine sensitivity in Kv12.1, these findings demonstrated that the quinoline‐binding pocket also exists in Kv12.1 and that thus the overall architecture of the channels is similar.
Figure 5

The inhibition induced by quinine of Kv12.1(H462Y) channels shows Kv11.1‐like characteristics. This figure summarises our findings on Kv12.1 wild‐type channels and H462Y mutant as well as on Kv11.1 wild‐type channels (data are reproduced from Figures 1, 2, 3, 4). (a) The affinity of Kv12.1(H462Y) for quinine was 17‐fold higher than that of the wild‐type and accordingly very similar to Kv11.1. (b) Quinine induced activation of Kv12.1 channels at more positive potentials following the hyperpolarised (−60 mV) and the depolarised conditioning potential (0 mV). In contrast, Kv12.1(H462Y) and Kv11.1 channels activated at more negative membrane potentials in the presence of quinine following the hyperpolarised conditioning potential. Only slight changes of voltage‐dependence of Kv12.1(H462Y) and Kv11.1 channels were detected during application of quinine following the depolarised conditioning potential. (c, d) Quinine‐mediated block was voltage‐independent for Kv12.1 but voltage‐dependent for Kv12.1(H462Y). Note that quinine inhibition of Kv11.1 channels was voltage‐dependent only following the hyperpolarised conditioning potential. Summaries of quinine‐induced inhibition (% of controls) for (c) the hyperpolarised and for (d) the depolarised conditioning membrane potential. Quinine was applied at 1 mM for Kv12.1 channels and at 50 μM for Kv12.1(H462Y) and Kv11.1. Mode shift was induced through conditioning potentials of −60 and 0 mV for Kv12.1 isoforms and at −60 and +40 mV for Kv11.1 channels

The inhibition induced by quinine of Kv12.1(H462Y) channels shows Kv11.1‐like characteristics. This figure summarises our findings on Kv12.1 wild‐type channels and H462Y mutant as well as on Kv11.1 wild‐type channels (data are reproduced from Figures 1, 2, 3, 4). (a) The affinity of Kv12.1(H462Y) for quinine was 17‐fold higher than that of the wild‐type and accordingly very similar to Kv11.1. (b) Quinine induced activation of Kv12.1 channels at more positive potentials following the hyperpolarised (−60 mV) and the depolarised conditioning potential (0 mV). In contrast, Kv12.1(H462Y) and Kv11.1 channels activated at more negative membrane potentials in the presence of quinine following the hyperpolarised conditioning potential. Only slight changes of voltage‐dependence of Kv12.1(H462Y) and Kv11.1 channels were detected during application of quinine following the depolarised conditioning potential. (c, d) Quinine‐mediated block was voltage‐independent for Kv12.1 but voltage‐dependent for Kv12.1(H462Y). Note that quinine inhibition of Kv11.1 channels was voltage‐dependent only following the hyperpolarised conditioning potential. Summaries of quinine‐induced inhibition (% of controls) for (c) the hyperpolarised and for (d) the depolarised conditioning membrane potential. Quinine was applied at 1 mM for Kv12.1 channels and at 50 μM for Kv12.1(H462Y) and Kv11.1. Mode shift was induced through conditioning potentials of −60 and 0 mV for Kv12.1 isoforms and at −60 and +40 mV for Kv11.1 channels

H462 determines the low quinine affinity of human Kv12.1 channels

To elucidate the molecular principles of the low quinine affinity of Kv12.1, we performed molecular docking and subsequent MD analyses (Figure 6). Quinine binding to the pore helix/S6 domain (amino acids S353–Y477) of Kv12.1 and Kv12.1(H462Y) channels was studied using homology models based on the cryo‐EM structure of human Kv11.1 (Wang & MacKinnon, 2017). We also utilised the available structural information to study quinine binding to Kv11.1 channels. Molecular docking suggested that quinine binds below residue H462 into the central cavity of human Kv12.1 (Figure 6a, left panel), with an estimated binding affinity (Chemscore.DG) of −20.5 kJ·mol−1. To validate the stability and refine the binding pose obtained from docking, we performed three independent 50‐ns MD simulations of the highest scored docking pose. Overall, quinine remained stable in these MD simulations, as indicated by root‐mean‐square deviation (RMSD) below 2 Å (Figure 6d, left panel). Further, π–π stacking between methoxyquinoline group and H462 from one subunit was maintained throughout 50 ns of MD simulation (Figure 6b,c, left panel). Similar to wild‐type channels, quinine docked into the central cavity of Kv12.1(H462Y) without any obvious interactions to the residues of the selectivity filter (Figure 6a–c, middle panel). In Kv12.1(H462Y), π–π interactions with the methoxyquinoline moiety were also observed; however, in the mutant channel, these interactions are formed with mutated Y462 side chains from two adjacent subunits. Additionally, we observed a hydrogen bond between Y462 and quinoline moiety of quinine in the majority of docking poses that remained stable throughout the 50‐ns MD simulations. In agreement with our experimental findings, the estimated quinine‐binding affinity of −31.1 kJ·mol−1 for Kv12.1(H462Y) was considerably higher than for wild‐type Kv12.1 channels. In line, in the simulations, an additional hydrogen bond between the methoxy group of quinidine and the residue T470 formed in the mutant channels (Figure 6c, middle panel). As shown in Figure 6d, the binding pose of quinine was highly stable in MD simulations, with an RMSD of about 1 Å. In the Kv11.1 cryo‐EM structure (Wang & MacKinnon, 2017), quinine was also predicted to bind slightly below Y652, in the central cavity, with favourable hydrophobic interactions with Y652 residues of all subunits and with F656 from only one subunit. Also, in line with our experimental findings, the estimated affinity of quinine binding to Kv11.1 wild‐type channels amounts to −33.9 kJ·mol−1 in our simulations and thus was similar to Kv12.1(H462Y) and higher than that of wild‐type Kv12.1 channels. Further, π–π interactions with Y652 and F656 residues in Kv11.1 from opposite subunits were predicted by these simulations (Figure 6b,c, right panel), again leading to a very stable binding mode with RMSD values of ~1 Å (Figure 6d, right panel).
Figure 6

Molecular interactions of quinine with Kv12.1, Kv12.1(H462Y), and Kv11.1. (a) Top‐view representations of Kv12.1 (left), Kv12.1(H462Y; middle), and Kv11.1 (right) channels with docked quinine (green–cyan spheres). (b) Top and (c) side view of the quinine‐binding site after 50 ns of molecular dynamics simulations for Kv12.1 (left), Kv12.1(H462Y; middle), and Kv11.1 (right). Directly interacting residues, as well as position 462 in Kv12.1 isoforms (652 in Kv11.1), are presented as sticks, with oxygen atoms coloured in red and nitrogen atoms coloured in blue. Yellow lines indicate π‐stacking interactions, while red lines indicate H bonds. (d) Stability of the pore domains (excluding extracellular loops), measured as the root‐mean‐square deviation (RMSD) as a function of 50‐ns simulation time

Molecular interactions of quinine with Kv12.1, Kv12.1(H462Y), and Kv11.1. (a) Top‐view representations of Kv12.1 (left), Kv12.1(H462Y; middle), and Kv11.1 (right) channels with docked quinine (green–cyan spheres). (b) Top and (c) side view of the quinine‐binding site after 50 ns of molecular dynamics simulations for Kv12.1 (left), Kv12.1(H462Y; middle), and Kv11.1 (right). Directly interacting residues, as well as position 462 in Kv12.1 isoforms (652 in Kv11.1), are presented as sticks, with oxygen atoms coloured in red and nitrogen atoms coloured in blue. Yellow lines indicate π‐stacking interactions, while red lines indicate H bonds. (d) Stability of the pore domains (excluding extracellular loops), measured as the root‐mean‐square deviation (RMSD) as a function of 50‐ns simulation time In summary, molecular modelling revealed that weak interactions between quinine and particularly H462 determined the low affinity and sensitivity of Kv12.1 wild‐type channels to quinine. The higher affinity of Kv12.1(H462Y) channels and of the closely related Kv11.1 can be explained by more favourable interactions with the drug, mainly at position Y462, and in the case of Kv11.1 additionally with F656 from one subunit. Surprisingly, despite the conservation of this second aromatic side chain between Kv11 and Kv12 channels, interactions with this second aromatic residue are completely absent in the Kv12.1 (Figure 6b, left panel vs. right panel). However, also in Kv11.1 channels, these interactions are relatively weak due to a preferred orientation of the Y652/F656 side chains towards helix S5 (see Figure S2, for orientation of aromatic side chains at positions H462/F466 in Kv12.1 wild‐type and Y652/F656 in Kv11.1 wild‐type channels in simulations).

DISCUSSION

Substituted quinolines are well‐known antimalarial agents (quinine and chloroquine) and antiarrhythmic drugs (quinidine; Bozic, Uzelac, Kezic, & Bajcetic, 2018). As considerable side effect, these substances inhibit Kv11.1 channels at low micromolar concentrations, which may result in acquired (drug‐induced) LQT syndrome, syncope, and sudden death in humans (Mitcheson et al., 2005; Mitcheson et al., 2000; Sanguinetti & Tristani‐Firouzi, 2006). Several mutagenesis studies have attributed the susceptibility of Kv11.1 to drug‐dependent block to amino acids in the pore domain of these channels, threonine 623 and valine 625 of inner pore helix and selectivity filter, respectively, as well as tyrosine 652 and phenylalanine 656 in the sixth transmembrane segment (S6; Mitcheson et al., 2000; Sanchez‐Chapula et al., 2002; Sanchez‐Chapula et al., 2003). On the one hand, valine 625 (together with other residues such as asparagine 588, serine 631, and serine 620) controls drug sensitivity of Kv11.1 through its contribution to channel inactivation that has been shown to be necessary for channel inhibition by several high‐affinity blockers (e.g., Kamiya et al., 2006; Perrin et al., 2008; Wu et al., 2015). On the other hand, some of these residues form hydrophobic pockets in combination with other amino acids at the inner surface of the small central cavity of Kv11.1 channels (Wang & MacKinnon, 2017). These pockets are endowed with high electronegative potentials and therefore favour binding of positively charged drugs (Wang & MacKinnon, 2017). The exclusive sensitivity of Kv11.1 to a wide variety of drugs with diverging structures is explained by the absence of such binding pockets in other K+ channels (Wang & MacKinnon, 2017). Closely related Kv10 and Kv12 channels share high similarity to Kv11.1, which suggests that these channels also possess analogous drug‐binding pockets, similar drug sensitivity, and characteristics of inhibition (Bauer & Schwarz, 2018). Indeed, the respective amino acids determining the susceptibility of Kv11.1 to quinoline block are conserved in Kv10 channels. However, whereas Kv10.1 displayed the same sensitivity to quinidine‐dependent inhibition as Kv11.1, quinidine sensitivity of Kv10.2 channels was 100‐fold lower (Gessner et al., 2004; Lees‐Miller et al., 2000; Sanchez‐Chapula et al., 2003; Schonherr et al., 2002). The lower susceptibility of the Kv10 isoforms to drug‐induced inhibition is explained by their lack of inactivation, that is, conformational reorientations during inactivation necessary for high affinity drug block in Kv11.1 do not occur in Kv10 channels (Chen et al., 2002; Ficker et al., 1998; Ficker et al., 2001). Thus, the lower sensitivity of Kv10 channels indicates that additional drug interaction sites outside the central binding motif might determine drug block in these channels (Gessner et al., 2004; Schonherr et al., 2002). We wondered whether Kv12.1 channels possess a quinoline‐binding pocket analogous to Kv10.1 and Kv11.1 channels despite an amino acid exchange at position 462. We found that quinine inhibited Kv12.1 channels, but their sensitivity was 10 times lower than that of Kv11.1 channels. Thus, Kv12.1 constitutes a natural variant of Eag superfamily channels with low quinoline‐binding affinity.

H462 determines quinine sensitivity of Kv12.1

The quinoline‐binding “motif” of Kv11.1 is not completely conserved in the three members of the Kv12 channel family that carry a histidine at the position homologous to tyrosine 652 in Kv11.1 (H462 in Kv12.1; c.f. Figure 3a). As in Kv11.1, this aromatic residue determines the quinoline‐binding affinity (Lees‐Miller et al., 2000; Sanchez‐Chapula et al., 2003), we hypothesised that this small sequence difference explained the 10‐fold lower quinine sensitivity of Kv12.1. Indeed, replacing this histidine in Kv12.1 with the respective tyrosine of Kv11.1 through site‐directed mutagenesis dramatically increased the quinine affinity of Kv12.1(H462Y). In fact, quinine sensitivity of Kv12.1(H462Y) was more similar to Kv11.1 than to wild‐type Kv12.1. Our findings are supported by earlier studies showing that mutating Y652 dramatically changes quinidine sensitivity of Kv11.1 channels, for example, introduction of an alanine at this position caused a threefold reduction of quinidine affinity in Kv11.1(Y652A) channels (Sanchez‐Chapula et al., 2003). Utilising molecular modelling simulations, we found that the low sensitivity of Kv12.1 channels may be explained by only weak interactions between quinine and amino acids in the central cavity (including H462). Noteworthy, in Kv12.1(H462Y) channels, more favourable interactions developed that most probably rendered this channel mutant more sensitive to quinine. In line with our experimental findings, docking analyses predicted considerably lower quinine affinities for Kv12.1 (−20.5 kJ·mol−1), compared with the Kv12.1(H462Y) mutant (−31.1 kJ·mol−1) and Kv11.1 wild‐type channels (−33.9 kJ·mol−1). Our models predicted a key difference in the orientation and interactions of the first aromatic side chain, which limits hydrophobic and π–π interactions in Kv12, while providing strong interactions in Kv11 channels. Taking together, our findings are in good agreement with a hypothesis previously published by the Sanguinetti group (Chen et al., 2002). Thus, our modelling approach recapitulated experimental findings on Kv11.1 (Sanchez‐Chapula et al., 2003) and provided a ready and straightforward explanation for the higher quinine affinity of Kv11.1 and Kv12.1(H462Y) compared with Kv12.1 wild‐type channels. Yet, as mentioned above, drug sensitivity of Kv11.1 channels is determined by channel inactivation (e.g., through V625, N588, S631, and S620; Ficker et al., 1998; Ficker et al., 2001; Kamiya et al., 2006; Perrin et al., 2008; Wu et al., 2015). As shown for Kv10 channels (Chen et al., 2002; Ficker et al., 1998; Ficker et al., 2001), lack of inactivation might contribute to low drug sensitivity of Kv12.1 channels. However, above‐mentioned residues (except N588 where Kv12.1 has an E) are conserved in Kv12.1 channels, and the single amino acid exchange H462Y sufficed to significantly lower the quinine affinity of these channels (even slightly below that of Kv11.1 channels; c.f. Figure 5a). We thus estimate that relevance of these amino acids for drug interactions is low in Kv12.1, but we cannot provide any evidence for this assumption at present. Thus, further work is needed to elucidate whether these residues also determine drug sensitivity in Kv12 family members and whether the channel mutant inactivates at all (e.g., during prolonged depolarisations). Importantly, our findings once again highlight that Y652 determines quinoline sensitivity in closely related Kv11.1 channels (Lees‐Miller et al., 2000; Macdonald, Kim, Kurata, & Fedida, 2018; Sanchez‐Chapula et al., 2002; Sanchez‐Chapula et al., 2003). Although indicated by an early report analysing binding of chloroquine to Kv11.1 in silico (Sanchez‐Chapula et al., 2002), we did not detect any cation–π interactions between quinine and Kv11.1, Kv12.1, or Kv12.1(H462Y) channels. Noteworthy, our simulations support a recent study that utilised unnatural amino acid incorporation to show that cation–π interactions at position Y652 are not relevant for drug binding in Kv11.1 channels (Macdonald et al., 2018). Based on these findings, we conclude that H462 determines the low quinine sensitivity of Kv12.1 and that possibly no amino acid motifs outside the central cavity contribute to quinoline binding in Kv12.1 channels (for Kv10, see Gessner et al., 2004; Schonherr et al., 2002). Thus, as far as quinoline sensitivity is concerned, Kv12 channels are probably more similar to other Kv families that typically carry isoleucine or valine at this position than to closely related Kv11 channels.

H462 determines characteristics of quinine‐dependent inhibition in Kv12.1

Our experiments showed that quinine‐induced inhibition of Kv12.1 channels was largely voltage independent and characterised by a shift of activation voltages to depolarised potentials. Although it has been shown that many inhibitors (including quinidine or chloroquine; Sanchez‐Chapula et al., 2003; Sanchez‐Chapula et al., 2002) preferentially block activated Kv11.1 channels, we did not find any evidence for such quinine‐dependent open channel block for Kv12.1 wild‐type channels. Transient activation of Kv12.1(H462Y) channels at depolarised membrane potentials in the presence of quinine (c.f. Figure 4b), however, might indicate such open channel block for the channel mutant, but we consider that further work including other inhibitors is needed to elucidate whether the amino acid position 462 indeed determines open channel block in Kv12.1 channels. As Kv12.1 wild‐type channels were more sensitive to quinine after conditioning hyperpolarisation of the membrane, quinine block depended on mode shift, in contrast to Kv12.1(H462Y) that were inhibited by quinine independent on conditioning potentials. Although we do not have any evidence for this assumption yet, conformational changes associated with establishment of mode shift might thus render Kv12.1 channels less sensitive to quinine. Quinine inhibited Kv11.1 more significantly, when the channels were activated at depolarised potentials, but only following the hyperpolarised conditioning potential. Thus, quinine block of Kv11.1 also depended on mode shift, but the underlying mechanisms may be different between Kv11 and Kv12 channels. In contrast to Kv12.1 wild‐type channels, quinine induced Kv11.1 channel opening at more negative membrane potentials in line with a previous report (Sanchez‐Chapula et al., 2003). Interestingly, the H462Y amino acid exchange conferred Kv11.1‐type characteristics of quinine‐dependent inhibition to Kv12.1 channels. Thus, H462 determines not only quinine sensitivity but also the characteristics of quinine‐dependent inhibition in Kv12.1 channels. Likewise, mutating the homologous tyrosine substantially changed the voltage dependence of quinoline block in Kv11.1, that is, voltage dependence was abolished in Kv11.1(Y652F), but it was completely reversed in Kv11.1(Y652A) channels (Sanchez‐Chapula et al., 2003). Yet the molecular mechanisms underlying such voltage‐dependent inhibition remain elusive for the Kv12.1 mutant, as well as for Kv11.1 channels. It was hypothesised earlier for chloroquine‐induced block of Kv11.1 that Y652 might change position during voltage‐dependent gating, thereby generating a depolarisation‐induced binding pocket with higher affinity (Sanchez‐Chapula et al., 2003). This may induce pronounced channel inhibition during channel activation at more depolarised membrane potentials. Accordingly, as proposed for Kv11.1 carrying mutations at position 652 (Sanchez‐Chapula et al., 2003), the low quinine affinity of Kv12.1 wild‐type channels might prevent generation of an analogous binding site with higher affinity, which might account for lack of voltage dependence in these channels. Thus, similar molecular mechanisms may apply for Kv11.1 and Kv12.1 during drug‐induced inhibition, but we cannot provide any evidence for such processes at present.

Conclusion and significance

We showed that Kv12.1 is a natural variant of Eag superfamily channels with low quinine sensitivity. Thus, the drug‐binding pocket in the central cavity is not completely conserved in the Eag superfamily of Kv channels, which highlights functional and pharmacological diversity within this group of evolutionary conserved ion channels. However, our work also demonstrated that the drug‐binding pocket exists in Kv12.1 and that thus the overall architecture of the channels is similar. Further work is needed to elucidate whether Kv12.1 channels are also endowed with lower sensitivity to the many more drugs that block closely related Kv11.1.

CONFLICT OF INTEREST

The authors declare no conflict of interest.

AUTHOR CONTRIBUTIONS

M.D., A.S.‐W., and M.G.L. participated in the research design; M.D., W.B.V.H., A.S.‐W., and M.G.L. conducted the experiments and performed the data analysis; M.D., A.S.‐W., and M.G.L. wrote or contributed to the writing of the manuscript; and all authors approved the final version of the manuscript.

DECLARATION OF TRANSPARENCY AND SCIENTIFIC RIGOUR

This Declaration acknowledges that this paper adheres to the principles for transparent reporting and scientific rigour of preclinical research as stated in the BJP guidelines for Design & Analysis and as recommended by funding agencies, publishers, and other organisations engaged with supporting research. Figure S1. Rapid inhibition of recombinant K 12.1, K 11.1 and K 12.1(H462Y) channels through quinine. (A,C,E) Within seconds, application of quinine inhibited currents through (A) Kv12.1, (C) Kv11.1 and (E) Kv12.1(H462Y) channels transiently overexpressed in CHO cells. The left panels show representative recordings of CHO cells expressing these channel subunits before application of quinine (control; black), at the end of quinine treatment (red) and several minutes after wash‐out of the substance (grey). The right panels show analyses of the time course of current amplitudes recorded in these experiments upon application of quinine. Data are presented as normalised to base line current amplitudes (I/Icontrol). (B,D,F) In the same cells, application of quinine induced reversible depolarisation of the membrane potential. All voltage‐protocols, scale bars, quinine concentrations and time points after wash‐out as indicated. Figure S2. Molecular dynamics of residues H462/F466 in (K 12.1 and K 12.1(H462Y)) and Y652/F656 (in K 11.1). Top view of superposed snapshots from 3 independent 50 ns molecular dynamics simulations of residues H462/F466 in Kv12.1 wt (left column) and Kv12.1(H462Y) (middle column) and of corresponding residues Y652/F656 in Kv11.1 wt channels (right column). Click here for additional data file.
  60 in total

1.  Restoration of ion channel function in deafness-causing KCNQ4 mutants by synthetic channel openers.

Authors:  Michael G Leitner; Anja Feuer; Olga Ebers; Daniela N Schreiber; Christian R Halaszovich; Dominik Oliver
Journal:  Br J Pharmacol       Date:  2012-04       Impact factor: 8.739

2.  Open channel block of HERG K(+) channels by vesnarinone.

Authors:  K Kamiya; J S Mitcheson; K Yasui; I Kodama; M C Sanguinetti
Journal:  Mol Pharmacol       Date:  2001-08       Impact factor: 4.436

Review 3.  Structural determinants for high-affinity block of hERG potassium channels.

Authors:  John Mitcheson; Matthew Perry; Phillip Stansfeld; Michael C Sanguinetti; Harry Witchel; Jules Hancox
Journal:  Novartis Found Symp       Date:  2005

4.  Differential expression of genes encoding subthreshold-operating voltage-gated K+ channels in brain.

Authors:  M J Saganich; E Machado; B Rudy
Journal:  J Neurosci       Date:  2001-07-01       Impact factor: 6.167

5.  A structural basis for drug-induced long QT syndrome.

Authors:  J S Mitcheson; J Chen; M Lin; C Culberson; M C Sanguinetti
Journal:  Proc Natl Acad Sci U S A       Date:  2000-10-24       Impact factor: 11.205

6.  A mechanism underlying compound-induced voltage shift in the current activation of hERG by antiarrhythmic agents.

Authors:  Kazuharu Furutani; Yuko Yamakawa; Atsushi Inanobe; Miki Iwata; Yuko Ohno; Yoshihisa Kurachi
Journal:  Biochem Biophys Res Commun       Date:  2011-10-12       Impact factor: 3.575

7.  Deletion of the potassium channel Kv12.2 causes hippocampal hyperexcitability and epilepsy.

Authors:  Xiaofei Zhang; Federica Bertaso; Jong W Yoo; Karsten Baumgärtel; Sinead M Clancy; Van Lee; Cynthia Cienfuegos; Carly Wilmot; Jacqueline Avis; Truc Hunyh; Catherine Daguia; Christian Schmedt; Jeffrey Noebels; Timothy Jegla
Journal:  Nat Neurosci       Date:  2010-08-01       Impact factor: 24.884

8.  Voltage-sensing domain mode shift is coupled to the activation gate by the N-terminal tail of hERG channels.

Authors:  Peter S Tan; Matthew D Perry; Chai Ann Ng; Jamie I Vandenberg; Adam P Hill
Journal:  J Gen Physiol       Date:  2012-08-13       Impact factor: 4.086

9.  Probing the molecular basis of hERG drug block with unnatural amino acids.

Authors:  Logan C Macdonald; Robin Y Kim; Harley T Kurata; David Fedida
Journal:  Sci Rep       Date:  2018-01-10       Impact factor: 4.379

10.  New potential binding determinant for hERG channel inhibitors.

Authors:  P Saxena; E-M Zangerl-Plessl; T Linder; A Windisch; A Hohaus; E Timin; S Hering; A Stary-Weinzinger
Journal:  Sci Rep       Date:  2016-04-12       Impact factor: 4.996

View more
  1 in total

1.  Histidine at position 462 determines the low quinine sensitivity of ether-à-go-go channel superfamily member Kv 12.1.

Authors:  Marlen Dierich; Willem B van Ham; Anna Stary-Weinzinger; Michael G Leitner
Journal:  Br J Pharmacol       Date:  2019-06-17       Impact factor: 8.739

  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.