The s2–s3 loop of kv7.4 channels is essential for calmodulin regulation of channel activation

Introduction

Kv7 (KCNQ) voltage-gated potassium channels (Kv7. 1–7. 5; Gutman et al., 2003 ) produce Iks current (Kv7. 1; Barhanin et al., 1996 ; Sanguinetti et al., 1996 ; Wang et al., 1996 ; Jentsch, 2000 ) and M current (Kv7. 2–7. 5; Wang et al., 1998 ; Delmas and Brown, 2005 ; Hernandez et al., 2008 ), playing critical roles in controlling cellular excitability in the brain, heart, and ear. Kv7 channels open at subthreshold membrane potentials and function as a brake on membrane excitation ( Delmas and Brown, 2005 ; Brown and Passmore, 2009 ). Kv7. 4 channels are expressed in the inner ear and the central auditory pathway in the brainstem ( Kubisch et al., 1999 ; Beisel et al., 2000 , 2005 ; Kharkovets et al., 2000 ). Mutations in Kv7. 4 result in an autosomal-dominant, non-syndromic, progressive, high-frequency hearing loss, denoted as DFNA2 ( Kubisch et al., 1999 ).

Calmodulin (CaM) appears to be an essential auxiliary subunit of Kv7 channels and strongly modulates Kv7 function. The Kv7 A and B helices, which are located in the cytoplasmic C-terminal domain and close to the pore, are identified as CaM binding sites ( Wen and Levitan, 2002 ; Yus-Najera et al., 2002 ; Sachyani et al., 2014 ; Strulovich et al., 2016 ; Sun and MacKinnon, 2017 ; Chang et al., 2018 ). Disrupted CaM interactions with Kv7. 1 ( Ghosh et al., 2006 ; Shamgar et al., 2006 ), Kv7. 2 ( Etxeberria et al., 2008 ; Alaimo et al., 2009 , 2018 ), heteromeric Kv7. 2/Kv7. 3 ( Liu and Devaux, 2014 ), or Kv7. 4 ( Chang et al., 2018 ) have been shown to be tightly associated with severely impaired Kv7 surface expression, which illustrates the vital role of CaM in Kv7 assembly and trafficking. Numerous studies about the effects of calcium on the interaction between CaM and Kv7 revealed that both Ca ²⁺/CaM and Apo/CaM forms bind to the channels ( Wen and Levitan, 2002 ; Yus-Najera et al., 2002 ; Ghosh et al., 2006 ; Shamgar et al., 2006 ; Bal et al., 2008 ), whereas details of the structure of the Kv7: CaM complex are calcium dependent. Apo/CaM embraces an antiparallel pair of the Kv7 A and B helices, with the Apo/C-lobe interacting with the A helix and the Apo/N-lobe interacting with the B helix ( Bernardo-Seisdedos et al., 2018 ; Chang et al., 2018 ). Calcium binding allows the Ca ²⁺/N-lobe to stay anchored on the Kv7 B helix in a similar configuration to the apo/N-lobe ( Sachyani et al., 2014 ; Strulovich et al., 2016 ; Sun and MacKinnon, 2017 ; Bernardo-Seisdedos et al., 2018 ; Chang et al., 2018 ), while there are striking differences between the Ca ²⁺/C-lobe and Apo/C-lobe conformations. Structural studies of the interaction between CaM and the Kv7. 4 A and B helices showed that calcium binding makes the C-lobe lose interaction with the A helix and bind to the B helix weakly through a much smaller interaction surface ( Xu et al., 2013 ; Chang et al., 2018 ). CaM acts as a calcium sensor and regulates Kv7 currents in a calcium-dependent manner. Ca ²⁺/CaM facilitates Kv7. 1 function ( Sachyani et al., 2014 ; Tobelaim et al., 2017 ; Chang et al., 2018 ) but inhabits Kv7. 2–7. 5 ( Gamper and Shapiro, 2003 ; Gamper et al., 2005 ; Sihn et al., 2016 ; Gomis-Perez et al., 2017 ; Chang et al., 2018 ). Nevertheless, the exact mechanism by which CaM regulates Kv7 channels remains unclear.

The cryo-EM structure of the Kv7. 1/CaM complex reports a newly discovered interaction between the third EF hand (EF3) of CaM and the Kv7. 1 S2–S3 loop ( Sun and MacKinnon, 2017 ), which is unique to Kv7 channels but absent in the other voltage-gated potassium channels ( Supplementary Figure S1 ), raising a possibility that CaM contacts the voltage sensor of Kv7 channels through the S2–S3 loop. In our study, using electrophysiological recordings and co-immunoprecipitation assays, we demonstrate that CaM regulates Kv7. 4 activation through the S2–S3 loop and propose a new model accounting for CaM-mediated regulation of Kv7 channels.

Materials and Methods

Plasmid Constructs

The plasmids (EX-U0210-M98-5) expressing human Kv7. 4 (NP_004691. 2) were purchased from GeneCopoeia. Each amino acid of the Kv7. 4 S2–S3 loop was substituted by an alanine or a valine to generate Kv7. 4 mutants. DNA encoding CaM (NP_008819. 1) was a gift from Jiahuai Han (Xiamen University, China). CaM mutants with impaired calcium affinities had an alanine substitution of the first residue, an aspartate, and in the calcium-binding loop of each EF hand (CaM ₃ , D94A; CaM ₄ , D130A; CaM ₁₂ , D21A/D57A; CaM ₃₄ , D94A/D130A; and CaM ₁₂₃₄ , D21A/D57A/D94A/D130A; Geiser et al., 1991 ; Xia et al., 1998 ). Mutations G97A, N98A, G99A, and Y100A were introduced into wild-type (WT) CaM (referred to as CaM G97A, CaM N98A, CaM G99A, and CaM Y100A) or CaM ₃ (referred to as CaM ₃ G97A, CaM ₃ N98A, CaM ₃ G99A, and CaM ₃ Y100A). All mutants were generated by homologous recombination and verified by automated sequencing. DNA segments encoding WT or mutant Kv7. 4 were subcloned into pIRES2-EGFP plasmids, and segments encoding WT or mutant CaM were subcloned into pIRES2-mCherry plasmids.

There are two alternatively spliced isoforms of human Kv7. 4, called Kv7. 4 isoform a and Kv7. 4 isoform b (NP_751895. 1). Most studies including ours used Kv7. 4 isoform a, which we and others called Kv7. 4. Our study also performed experiments on Kv7. 4 isoform b, which we called Kv7. 4b to distinguish from isoform a. DNA segments encoding Kv7. 4b were generated by deleting residues 377–430 of Kv7. 4 isoform a by homologous recombination and subcloned into pIRES2-EGFP plasmids. The sequence was verified by automated sequencing.

For surface labeling, the Kv7. 4 subunit was tagged with a modified Myc epitope that was flanked with the extracellular D1–D2 loop of ClC-5 chloride channels to increase accessibility in the extracellular loop that connects the transmembrane domains S1 and S2 as previously described ( Kim et al., 2011 ; referred to as Myc-Kv7. 4). The amino acid sequence of the Kv7. 4 S1–S2 loop was changed to STIQEHQELANENSEHEQKLISEEDLVTFEERDKCPEWNC. The Myc-Kv7. 4 segments were subcloned into pcDNA3. 1 plasmids without fluorescent tags. DNA encoding CD8α was kindly provided by Lan Bao (the Chinese Academy of Sciences, China), and an HA epitope was inserted between the signal peptide and the mature protein of CD8α (referred to as HA-CD8α). The HA-CD8α segments were subcloned into pcDNA3. 1 plasmids as well.

For measurements of Kv7. 4 and CaM total expressions and co-immunoprecipitation experiments, the Kv7. 4 subunit was tagged with a Flag epitope (DYKDDDDK) at the N-terminus (referred to as Flag-Kv7. 4), and CaM was tagged with an HA epitope (YPYDVPDYA) at the C-terminus (referred to as CaM-HA).

Electrophysiology

Chinese hamster ovary (CHO) cells were obtained from the Cell Bank of the Chinese Academy of Sciences, and cultured in DMEM/F-12 (Gibco) with 10% fetal bovine serum (FBS) at 37°C with 5% CO ₂ . Cells in 35 mm diameter wells were transfected using Lipofectamine 3000 (Invitrogen) with 1 μg of DNA encoding Kv7. 4 channels and 1 μg of either DNA encoding CaM or empty vectors per dish.

Twenty-four hours following transfection, the media were replaced with bath solution. The bath solution contained 145 mM NaCl, 4 mM KCl, 1. 8 mM CaCl ₂ , 0. 5 mM MgCl ₂ , 10 mM HEPES, and 5 mM D-glucose (pH 7. 4 adjusted with NaOH). Electrophysiology experiments were performed at room temperature by whole-cell voltage clamp recordings using polished pipettes filled with pipette solution. The pipettes were pulled from borosilicate glass capillaries (BF150-86-10, Sutter Instrument, Novato) using a Flaming/Brown micropipette puller (P-97, Sutter Instrument, Novato) and polished to obtain 3–5 MΩ resistance. The normal pipette solution contained 140 mM KCl, 1 mM MgCl ₂ , 10 mM HEPES, 10 mM EGTA, 3 mM CaCl ₂ , and 4 mM K ₂ ATP (pH 7. 2 adjusted with KOH).

Pipette solutions containing various free calcium concentrations were used to change intracellular calcium levels. Whole-cell recordings were performed on transfected cells using pipette solutions containing 140 mM KCl, 1 mM MgCl ₂ , 10 mM HEPES, 4 mM K ₂ ATP, and 10 mM EGTA plus 0, 3, or 6 mM CaCl ₂ . The free calcium concentrations of the pipette solutions were 0, 77. 6, and 275. 9 nM, respectively, as calculated by the Webmax software online. ¹

BAPTA-AM (Abcam) was used to decrease intracellular calcium levels. BAPTA-AM was prepared as a stock solution in 20 mM in DMSO. One microliter of the BAPTA-AM stock solution was added to 1 ml of DMEM/F-12 media without FBS to obtain a final concentration of 20 μM BAPTA-AM. Twenty-four hours following transfection, cells were pretreated with the media containing 20 μM BAPTA-AM for 6 h at 37°C with 5% CO ₂ . Whole-cell recordings were performed on transfected cells pretreated with BAPTA-AM using the pipette solution containing 140 mM KCl, 1 mM MgCl ₂ , 10 mM HEPES, 4 mM K ₂ ATP, and 10 mM EGTA (pH 7. 2 adjusted with KOH).

All recordings were carried out using pClamp10. 5 software [(Axon Instruments, United States), an Axopatch 200B amplifier (Axon Instruments, United States), and an Axon Digidata 1550 digitizer (Axon Instruments, United States)]. Data were filtered at 1 kHz and digitized at 10 kHz. Kv7. 4 channel current traces were recorded from a holding potential of −80 mV to depolarizing step potentials ranging from −100 to 40 mV with 20 mV increments for 1, 000 ms, followed by −50 mV pulses for 500 ms to generate tail currents. Currents were measured after at least 80% series resistance compensation. Leakage currents were digitally subtracted after P/N leak subtraction. Kv7. 4 current amplitudes at varying test potentials were measured at steady-state levels, and the currents were then divided by the cell capacitance (pF) to generate a current density-voltage relationship. Tail currents were measured immediately after pulsing to −50 mV. To obtain voltage-dependent activation curves, the normalized tail current vs. voltage was plotted and fitted with the following Boltzmann equation: I/Imax = 1/{1 + exp[(V _1/2 − Vm)/k]}, where V _1/2 is the half-activation potential, Vm is the membrane potential, and k is the slope factor. Parameters of Kv7. 4 activation kinetics were determined by fitting the activation phase of the traces with a two-component exponential function: I(t) = A ₁ exp (−t/ τ ₁ ) + A ₂ exp (−t/τ ₂ ) + C, where I is the recorded current, A is the current amplitude, τ is the time constant, and C is the amplitude at which the activation starts. A ₁ , τ ₁ and A ₂ , τ ₂ represent the parameters for the fast and slow activation components, respectively. The deactivation time constants were measured by fitting the deactivating phase to a single exponential at −50 mV, as described above.

Immunofluorescence and Quantification of Surface Expression

Chinese hamster ovary cells in 35 mm diameter wells were transfected with 1 μg of DNA encoding Myc-tagged Kv7. 4 channels, 0. 3 μg of DNA encoding HA-tagged CD8α, and 1 μg of either DNA encoding CaM or empty vectors per dish. Forty-eight hours after co-transfection with Myc-Kv7. 4, HA-CD8α, and CaM, live CHO cells were placed in phosphate-buffered saline (PBS) with 3% bovine serum albumin (BSA) for 1 h to block unspecific binding. The cells were then incubated for 30 min with primary antibodies, mouse monoclonal anti-Myc (30601ES20, Yeasen), and rabbit monoclonal anti-HA (ab236632, Abcam). After three washes with PBS, the cells were incubated for 30 min with the secondary antibodies, Alexa-594-conjugated donkey anti-mouse IgG (H + L; 34112ES60, Yeasen), and Alexa-488-conjugated donkey anti-rabbit IgG (H + L; 34206ES60, Yeasen). All procedures were performed at room temperature.

Confocal images were acquired on a laser scanning confocal microscope (FV1200, Olympus, Japan) using a 60× oil immersion objective, in multitracking mode to minimize channel crosstalk. The surface expression of Kv7. 4 channels was quantified by the ratio of Myc fluorescence intensity to HA fluorescence intensity of all the cells captured in random subsets, and fluorescence intensities were measured using ImageJ (National Institutes of Health, United States). HA-CD8α was used here as a control for its stable and exclusive expression on the cell surface ( Zhang et al., 2008 ).

Co-immunoprecipitation and Western Blot

Chinese hamster ovary cells in 100 mm diameter wells were transfected with 9 μg of DNA encoding Flag-tagged Kv7. 4 channels and 9 μg of DNA encoding HA-tagged CaM per dish. Forty-eight hours after transfection, cells were lysed in lysis buffer containing 50 mM Tris-HCl (pH 8), 150 mM NaCl, 0. 5% Nonidet P-40, 200 μM Na ₃ VO ₄ , and protease inhibitor cocktail (B14001, Bimake), rotating at 4°C for 30 min. Cell lysates were centrifuged at 14, 000 g for 10 min at 4°C. Collected supernatants were incubated with rabbit monoclonal anti-Flag antibodies (ab205606, Abcam) in the presence of either 2 mM CaCl ₂ or 2 mM EGTA overnight at 4°C. Protein A + G agarose beads (P2012, Beyotime) were added, and the mixture was incubated for additional 3 h at 4°C. The beads that pulled down the immune complexes were collected by centrifugation at 1, 000 g for 2 min, and then washed three times with lysis buffer. SDS-PAGE loading buffer was added and boiled for 5 min to elute proteins. Immunoprecipitated proteins were separated by SDS-PAGE and blotted onto nitrocellulose membranes (10600001, GE Healthcare). The blots were blocked in TBST (10 mM Tris-HCl, 150 mM NaCl, and 0. 1% Tween 20, pH 7. 5) with 5% nonfat milk at room temperature for 1 h, and then incubated with mouse anti-Flag antibodies (30502ES20, Yeasen) and mouse anti-HA antibodies (30701ES20, Yeasen) overnight. After three washes with TBST, the blots were incubated with horseradish peroxidase (HRP)-labeled goat anti-mouse IgG (H + L; 33201ES60, Yeasen) for 1 h at room temperature. Proteins were visualized using Immun-Star HRP Substrate (1705041, Bio-Rad) and ChemiDoc Imaging System (Bio-Rad, United States). Chemiluminescent signals were collected by Image Lab software (Bio-Rad, United States) and analyzed by ImageJ (National Institutes of Health, United States).

Statistics

The statistical significances were determined by unpaired Student’s t -test for comparisons between two groups and by one-way ANOVA for comparisons of more groups. All statistical data in the figures were presented as mean ± SEM.

Results

CaM Regulates Kv7. 4 Activation in a Calcium-Dependent Manner

Previous studies reported CaM conferred calcium sensitivity to Kv7. 4 channels and affected Kv7. 4 activity ( Sihn et al., 2016 ; Chang et al., 2018 ). To test the calcium-dependent regulation mediated by CaM, we first measured currents of CHO cells transfected with Kv7. 4 alone or together with WT CaM or the mutant CaM ₁₂₃₄ that lost the calcium binding ability of all four EF hands. CHO cells have endogenous CaM ( Xu et al., 2007 ), and overexpression of endogenous WT CaM or CaM ₁₂₃₄ could knock down the effects of endogenous CaM. Shown in Figure 1A are whole-cell current traces of Kv7. 4 channels with only endogenous CaM but no overexpressed exogenous CaM and Kv7. 4 channels with overexpression of exogenous WT CaM or CaM ₁₂₃₄ . Compared with Kv7. 4 channels with only endogenous CaM of CHO cells, overexpression of exogenous WT CaM and CaM ₁₂₃₄ significantly increased current amplitudes, and Kv7. 4 co-expressed with CaM ₁₂₃₄ generated the largest current amplitudes ( Figures 1B , C ; current densities at 40 mV: Kv7. 4 = 76. 7 ± 6. 6 pA/pF; Kv7. 4: CaM = 178. 2 ± 17. 9 pA/pF; and Kv7. 4: CaM ₁₂₃₄ = 377. 6 ± 35. 9 pA/pF). Kv7. 4 with overexpression of exogenous WT CaM showed activation properties similar to Kv7. 4 channels with only endogenous CaM in terms of activation kinetics ( Figure 1D ; τ values at 40 mV: Kv7. 4 = 72. 4 ± 3. 7 ms; Kv7. 4: CaM = 72. 0 ± 8. 6 ms) and the voltage dependence of activation ( Figures 1E , F ; V _1/2 values: Kv7. 4 = −19. 9 ± 0. 8 mV; Kv7. 4: CaM = −18. 8 ± 1. 6 mV), while overexpression of exogenous CaM ₁₂₃₄ sharply accelerated the activation rate ( Figure 1D ; τ value at 40 mV: Kv7. 4: CaM ₁₂₃₄ = 20. 8 ± 3. 1 ms) and shifted the half-activation voltage leftward markedly ( Figures 1E , F ; V _1/2 value: Kv7. 4: CaM ₁₂₃₄ = −47. 9 ± 1. 0 mV). Overexpression of exogenous WT CaM or CaM ₁₂₃₄ did not change deactivation kinetics ( Figure 1G ).

FIGURE 1

Calmodulin (CaM) regulates Kv7. 4 activation and trafficking.(A)Representative current traces of Kv7. 4 channels without overexpression of exogenous CaM and Kv7. 4 channels with overexpression of exogenous CaM or CaM ₁₂₃₄ using the indicated protocol.(B)Corresponding current density-voltage relations of Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(C)Average current densities at 40 mV for Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(D)Average time constants at 40 mV for the fast activation component of Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(E)Voltage-dependent activation curves for Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(F)Individual V _1/2 values for Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(G)Average time constants at 40 mV for the deactivation of Kv7. 4 ( N = 20), Kv7. 4: CaM ( N = 20), and Kv7. 4: CaM ₁₂₃₄ ( N = 12).(H)Immunofluorescence analysis of surface expression of Kv7. 4 channels without overexpression of exogenous CaM and Kv7. 4 channels with overexpression of exogenous CaM or CaM ₁₂₃₄ . Live transfected cells were double labeled with anti-HA (green) and anti-Myc (red) antibodies to detect surface expression of HA-CD8α and Myc-Kv7. 4, respectively. The scale bars indicate 10 μm.(I)Quantification of surface expression of K7. 4 channels without overexpression of exogenous CaM ( n = 30) and Kv7. 4 channels with overexpression of exogenous CaM ( n = 30) or CaM ₁₂₃₄ ( n = 30), in which the ratio (Myc intensity/HA intensity) of each subset was normalized to the average ratio of Kv7. 4 without overexpression of exogenous CaM. n indicates the number of captured subsets obtained from five replications of experiments. Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^** p < 0. 01;^*** p < 0. 001. NS indicates not significant.

We also measured the surface expression of Kv7. 4 without overexpression of exogenous CaM and Kv7. 4 with overexpression of exogenous WT CaM or CaM ₁₂₃₄ to further interpret the different current amplitudes. To quantify surface expression, we developed a new method based on immunofluorescence assays. CHO cells were co-transfected with Myc-Kv7. 4, HA-CD8α and empty vectors, WT CaM or CaM ₁₂₃₄ . The Kv7. 4 subunit was tagged with a Myc epitope in the extracellular loop that connects the transmembrane domains S1 and S2, which has been used to detect Kv7. 4 surface expression in previous studies ( Kim et al., 2011 ). Our electrophysiological data showed that the extracellular Myc tag had no effects on Kv7. 4 currents or CaM regulation of Kv7. 4 ( Supplementary Figure S2 ). An extracellular HA tag was inserted between the signal peptide and the mature protein of CD8α. Previous studies have shown that CD8α could be stably and exclusively present on the cell surface ( Zhang et al., 2008 ). Hence, HA-CD8α was used as an indicator to observe the expression of exogenous protein and as a ruler to quantify the membrane area of the transfected cells. The fluorescence intensity of Myc-Kv7. 4 was divided by the fluorescence intensity of the co-expressed HA-CD8α, and this ratio could represent an average surface expression of Kv7. 4 channels. Statistically, overexpression of exogenous WT CaM or CaM ₁₂₃₄ significantly increased Kv7. 4 surface expression about 3. 5-fold ( Figures 1H , I ). Therefore, overexpression of exogenous WT CaM increased Kv7. 4 current amplitudes by enhancing Kv7. 4 surface expression without changing activation properties, while overexpression of exogenous CaM ₁₂₃₄ increased Kv7. 4 current amplitudes by enhancing Kv7. 4 surface expression as well as facilitating channel activation. There was no significant difference in surface expression between Kv7. 4 co-expressed with WT CaM and Kv7. 4 co-expressed with CaM ₁₂₃₄ ( Figures 1H , I ), which emphasized that only facilitated activation of Kv7. 4 co-expressed with CaM ₁₂₃₄ accounted for its greater current amplitudes than Kv7. 4 co-expressed with WT CaM. Thus, we focused on the effects of WT or mutant CaM on Kv7. 4 activation in the subsequent study.

We next sought to verify the role of calcium in CaM regulation of Kv7. 4 activation to explain the differences in activation properties between Kv7. 4 co-expressed with WT CaM and Kv7. 4 co-expressed with CaM ₁₂₃₄ . We prepared different pipette solutions containing 10 mM EGTA plus 0, 3, or 6 mM CaCl ₂ , and the free calcium concentrations are calculated by the MaxChelator software to be about 0, 80, and 270 nM, respectively. We performed whole-cell recordings using the different pipette solutions to change intracellular calcium levels of the recorded cells. Increasing free calcium concentrations in the pipette solutions exerted inhibitory effects on Kv7. 4: CaM activation, decreasing the activation rate ( Figure 2A ; τ values at 40 mV: Kv7. 4: CaM with 10 mM EGTA = 52. 3 ± 1. 2 ms; Kv7. 4: CaM with 10 mM EGTA plus 3 mM CaCl ₂ = 76. 0 ± 2. 0 ms; Kv7. 4: CaM with 10 mM EGTA plus 6 mM CaCl ₂ = 102. 4 ± 9. 2 ms), and shifting the voltage-dependent activation rightwards ( Figure 2B ; V _1/2 values at 40 mV: Kv7. 4: CaM with 10 mM EGTA = −27. 8 ± 1. 4 mV; Kv7. 4: CaM with 10 mM EGTA plus 3 mM CaCl ₂ = −17. 3 ± 1. 6 mV; Kv7. 4: CaM with 10 mM EGTA plus 6 mM CaCl ₂ = −10. 1 ± 1. 1 mV). However, both activation kinetics ( Figure 2C ; τ values at 40 mV: Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA = 20. 5 ± 2. 5 ms; Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA plus 3 mM CaCl ₂ = 17. 7 ± 2. 6 ms; Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA plus 6 mM CaCl ₂ = 18. 4 ± 2. 2 ms) and voltage-dependent activation ( Figure 2B ; V _1/2 values at 40 mV: Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA = −48. 0 ± 1. 4 mV; Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA plus 3 mM CaCl ₂ = −47. 2 ± 1. 8 mV; Kv7. 4: CaM ₁₂₃₄ with 10 mM EGTA plus 6 mM CaCl ₂ = −45. 7 ± 1. 1 mV) of Kv7. 4: CaM ₁₂₃₄ remained unchanged when we used different pipette solutions. Therefore, we further confirmed that functional CaM mediated calcium-dependent inhibition to Kv7. 4 channel activation and that the calcium-insensitive CaM ₁₂₃₄ facilitated Kv7. 4 activation essentially by abolishing the calcium-dependent inhibition.

FIGURE 2

CaM mediates calcium-dependent inhibition to Kv7. 4 activation.(A)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4: CaM pretreated with BAPTA-AM ( N = 8; B) Voltage-dependent activation curves for Kv7. 4: CaM with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4: CaM pretreated with BAPTA-AM ( N = 8).(C)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ₁₂₃₄ with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4: CaM pretreated with BAPTA-AM ( N = 9)(D)Voltage-dependent activation curves for Kv7. 4: CaM ₁₂₃₄ with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4: CaM pretreated with BAPTA-AM ( N = 9). Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^* p < 0. 05;^** p < 0. 01. NS indicates not significant.

Compared with the pipette solutions with CaCl ₂ added, the pipette solution containing only 10 mM EGTA facilitated Kv7. 4: CaM activation. Nevertheless, there are still differences between Kv7. 4: CaM with the pipette solution containing only 10 mM EGTA and Kv7. 4: CaM ₁₂₃₄ in both the activation rate and the voltage dependence of activation, possibly because binding to Kv7. 4 made CaM obtain high affinities for calcium in vivo and the calcium chelating ability of 10 mM EGTA was not enough to make all CaM in Apo/CaM forms, which CaM ₁₂₃₄ mimicked, in such a limited time of the whole-cell recording. We pretreated transfected cells with 20 μM BAPTA-AM, a cell permeable calcium chelator, for 6 h and performed whole-cell recordings using the pipette solution containing only 10 mM EGTA. Pretreating with BAPTA-AM further facilitated Kv7. 4: CaM activation but had no effect on Kv7. 4: CaM ₁₂₃₄ activation in terms of the activation rate ( Figures 2A , C ; τ values at 40 mV: Kv7. 4: CaM pretreated with 20 μM BAPTA-AM = 28. 4 ± 2. 8 ms; Kv7. 4: CaM ₁₂₃₄ pretreated with 20 μM BAPTA-AM = 19. 9 ± 2. 5 ms) and the voltage dependence of activation ( Figures 2B , D ; V _1/2 values: Kv7. 4: CaM pretreated with 20 μM BAPTA-AM = −37. 6 ± 1. 0 mV; Kv7. 4: CaM2 ₂₃₄ pretreated with 20 μM BAPTA-AM = −49. 8 ± 1. 7 mV).

Taken together, our results proved that CaM regulates Kv7. 4 channels in a calcium-dependent manner, in which Apo/CaM facilitated Kv7. 4 activation and Ca ²⁺/CaM inhibited Kv7. 4 activation.

The CaM EF3 Is Vital to the Regulation of Kv7. 4 Activation

Earlier studies have reported calcium-dependent structural properties of the Kv7. 4: CaM complex, in which calcium binding to the CaM N-lobe or C-lobe led to different interactions between CaM and Kv7. 4 ( Xu et al., 2013 ; Chang et al., 2018 ). To test if the calcium-dependent structural properties had effects on the calcium-dependent regulation of Kv7. 4 activation mediated by CaM, we sought to explore the role of each lobe or EF hand of CaM in regulating Kv7. 4 channels. Hence, we co-expressed Kv7. 4 with the CaM mutants that lost the calcium binding affinity of the N-lobe (CaM ₁₂ ) or the C-lobe (CaM ₃₄ ), and the effects of CaM ₁₂ and CaM ₃₄ were compared with the effects of WT CaM and CaM ₁₂₃₄ . There was no significant difference among the expression levels of WT CaM and the CaM mutants ( Supplementary Figure S4 ). Shown in Figure 3A are whole-cell current traces of Kv7. 4 channels co-expressed with WT CaM and the CaM mutants. CaM ₃₄ mirrored the effects of CaM ₁₂₃₄ on the current amplitudes ( Figures 3B , C ; current densities at 40 mV: Kv7. 4: CaM ₁₂₃₄ = 375. 7 ± 40. 8 pA/pF; Kv7. 4: CaM ₃₄ = 408. 9 ± 28. 0 pA/pF), the activation rate ( Figure 3D ; τ values at 40 mV: Kv7. 4: CaM ₁₂₃₄ = 20. 6 ± 2. 0 ms; Kv7. 4: CaM ₃₄ = 20. 4 ± 2. 0 ms), and the voltage dependence ( Figures 3E , F ; V _1/2 values: Kv7. 4: CaM ₃₄ = −46. 3 ± 1. 0 mV; CaM ₃₄ = −47. 1 ± 1. 2 mV). Compared with CaM ₁₂₃₄ and CaM ₃₄ , CaM ₁₂ caused much smaller changes to the indicated Kv7. 4 channel properties, only generating a small increase in the current amplitudes ( Figures 3B , C ; current densities at 40 mV: Kv7. 4: CaM ₁₂ = 293. 5 ± 22. 2 pA/pF), a mildly faster activation rate ( Figure 3D ; τ values at 40 mV: Kv7. 4: CaM ₁₂ = 42. 9 ± 4. 7 ms), and a slightly leftward shift in the voltage-dependent activation curve ( Figures 3E , F ; V _1/2 values: Kv7. 4: CaM ₁₂ = −36. 9 ± 1. 4 mV). These results suggested that the CaM C-lobe played a key role in the calcium-dependent regulation of Kv7. 4 activation, although calcium binding to the N-lobe also conferred a slight inhibition to Kv7. 4.

FIGURE 3

The CaM EF3 controls the regulation of Kv7. 4 activation.(A)Representative current traces of cells co-expressing Kv7. 4 channels with CaM or CaM mutants.(B)Corresponding current density-voltage relations of Kv7. 4 channels with CaM ( N = 10) or CaM mutants: CaM ₁₂₃₄ ( N = 11), CaM ₁₂ ( N = 14), CaM ₃₄ ( N = 17), CaM ₃ ( N = 15), and CaM ₄ ( N = 14).(C)Average current densities at 40 mV for Kv7. 4 channels with CaM ( N = 10) or CaM mutants: CaM ₁₂₃₄ ( N = 11), CaM ₁₂ ( N = 14), CaM ₃₄ ( N = 17), CaM ₃ ( N = 15), and CaM ₄ ( N = 14).(D)Average time constants at 40 mV for the fast activation component of Kv7. 4 with CaM ( N = 10) or CaM mutants: CaM ₁₂₃₄ ( N = 11), CaM ₁₂ ( N = 14), CaM ₃₄ ( N = 17), CaM ₃ ( N = 15), and CaM ₄ ( N = 14).(E)Voltage-dependent activation curves for Kv7. 4 with CaM ( N = 10) or CaM mutants: CaM ₁₂₃₄ ( N = 11), CaM ₁₂ ( N = 14), CaM ₃₄ ( N = 17), CaM ₃ ( N = 15), and CaM ₄ ( N = 14).(F)Individual V _1/2 values for Kv7. 4 with CaM ( N = 10) or CaM mutants: CaM ₁₂₃₄ ( N = 11), CaM ₁₂ ( N = 14), CaM ₃₄ ( N = 17), CaM ₃ ( N = 15), and CaM ₄ ( N = 14).(G)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ₃ with pipette solutions containing only 10 mM EGTA ( N = 8), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 8) and Kv7. 4: CaM ₃ pretreated with BAPTA-AM ( N = 9).(H)Voltage-dependent activation curves for Kv7. 4: CaM ₃ with pipette solutions containing only 10 mM EGTA ( N = 8), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 8) and Kv7. 4: CaM ₃ pretreated with BAPTA-AM ( N = 9). Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^** p < 0. 01;^*** p < 0. 001. NS indicates not significant.

Next, we used the CaM mutants with disabled calcium binding ability of the EF3 (CaM ₃ ) or the EF4 (CaM ₄ ) in the C-lobe. CaM ₃ phenocopied the effects of CaM ₁₂₃₄ and CaM ₃₄ on current amplitudes ( Figures 3B , C ; current densities at 40 mV: Kv7. 4: CaM ₃ = 385. 8 ± 30. 3 pA/pF), activation kinetics ( Figure 3D ; τ values at 40 mV: Kv7. 4: CaM ₃ = 18. 5 ± 1. 0 ms), and voltage-dependent activation ( Figures 3E , F ; V _1/2 values: Kv7. 4: CaM ₃ = −48. 0 ± 0. 9 mV). However, CaM ₄ , which caused effects almost equivalent to WT CaM, made no difference to the current amplitudes ( Figures 3B , C ; current densities at 40 mV: Kv7. 4: CaM ₄ = 216. 0 ± 24. 2 pA/pF), the activation rate ( Figure 3D ; τ values at 40 mV: Kv7. 4: CaM ₄ = 75. 2 ± 5. 4 ms), or the voltage dependence of activation ( Figures 3E , F ; V _1/2 values: Kv7. 4: CaM ₄ = −22. 2 ± 1. 2 mV). Like Kv7. 4: CaM ₁₂₃₄ , Kv7. 4: CaM ₃ showed insensitivity to the changing intracellular calcium levels, and both varying free calcium concentrations of the pipette solutions and application of BAPTA-AM have no effect on the activation rate ( Figure 3G ; τ values at 40 mV: Kv7. 4: CaM ₃ pretreated with 20 μM BAPTA-AM = 19. 6 ± 2. 1 ms; Kv7. 4: CaM ₃ with 10 mM EGTA = 19. 2 ± 1. 6 ms; Kv7. 4: CaM ₃ with 10 mM EGTA plus 3 mM CaCl ₂ = 17. 9 ± 1. 4 ms; Kv7. 4: CaM ₃ with 10 mM EGTA plus 6 mM CaCl ₂ = 18. 1 ± 2. 2 ms) or the voltage-dependent activation ( Figure 3H ; V _1/2 values: Kv7. 4: CaM ₃ pretreated with 20 μM BAPTA-AM = −45. 3 ± 1. 5 mV; Kv7. 4: CaM ₃ with 10 mM EGTA = −47. 3 ± 1. 7 mV; Kv7. 4: CaM ₃ with 10 mM EGTA plus 3 mM CaCl ₂ = −48. 5 ± 1. 2 mV; Kv7. 4: CaM ₃ with 10 mM EGTA plus 6 mM CaCl ₂ = −47. 1 ± 1. 3 mV) of Kv7. 4: CaM ₃ . These results emphasized a critical role of the CaM EF3 in modulating Kv7. 4 activation, which supported the possibility that the CaM EF3 interacted with the Kv7. 4 S2–S3 loop to contact the voltage sensor ( Sun and MacKinnon, 2017 ).

The Kv7. 4 S2–S3 Loop Is Required for CaM Regulation of Kv7. 4 Activation

We sought to interpret the importance of the CaM EF3 in regulating Kv7. 4 activation based on the potential interaction between the CaM EF3 and the Kv7 S2–S3 loop ( Sun and MacKinnon, 2017 ). Each residue of the nine-amino-acid S2–S3 loop was changed to an alanine or a valine. Because the current traces from CHO cells expressing C158A or G162A were quite noisy ( Supplementary Figure S6A ), we generated the mutants C158V and G162V for further investigation. When co-expressed with WT CaM, five of the nine Kv7. 4 mutants displayed facilitated channel activation ( Figure 4A ) with a faster activation rate ( Figure 4D ; τ values at 40 mV: Kv7. 4: CaM = 72. 7 ± 4. 8 ms, Kv7. 4 C156A: CaM = 31. 3 ± 2. 9 ms, Kv7. 4 C157A: CaM = 38. 2 ± 4. 6 ms, Kv7. 4 C158V: CaM = 33. 3 ± 1. 7 ms, Kv7. 4 R159A: CaM = 28. 3 ± 3. 0 ms, and Kv7. 4 R161A: CaM = 37. 4 ± 2. 1 ms) and a leftward shift in the voltage-dependent activation ( Figure 4F ; V _1/2 values: Kv7. 4: CaM = −18. 2 ± 1. 1 mV, Kv7. 4 C156A: CaM = −33. 9 ± 1. 2 mV, Kv7. 4 C157A: CaM = −29. 6 ± 1. 1 mV, Kv7. 4 C158V: CaM = −30. 6 ± 1. 2 mV, Kv7. 4 R159A: CaM = −35. 7 ± 1. 2 mV, and Kv7. 4 R161A: CaM = −33. 0 ± 1. 1 mV). However, these five mutants co-expressed CaM ₃ showed activation properties similar to WT Kv7. 4 co-expressed with CaM ₃ in terms of both the activation kinetics ( Figure 4E ; τ values at 40 mV: Kv7. 4: CaM ₃ = 18. 8 ± 1. 0 ms, Kv7. 4 C156A: Kv7. 4: CaM ₃ = 20. 2 ± 1. 4 ms, Kv7. 4 C157A: CaM ₃ = 19. 2 ± 1. 4 ms, Kv7. 4 C158V: CaM ₃ = 19. 3 ± 1. 8 ms, Kv7. 4 R159A: CaM ₃ = 21. 0 ± 1. 7 ms, and Kv7. 4 R161A: CaM ₃ = 19. 3 ± 2. 0 ms) and the voltage dependence of activation ( Figure 4G ; V _1/2 values: Kv7. 4: CaM ₃ = −47. 9 ± 0. 8 mV, Kv7. 4 C156A: CaM ₃ = −49. 1 ± 1. 5 mV, Kv7. 4 C157A: CaM ₃ = −48. 1 ± 1. 4 mV, Kv7. 4 C158V: CaM ₃ = −50. 1 ± 1. 3 mV, Kv7. 4 R159A: CaM ₃ = −51. 2 ± 1. 5 mV, and Kv7. 4 R161A: CaM ₃ = −49. 6 ± 0. 9 mV), indicating that these five mutations had no effect on activation of Kv7. 4 co-expressed with CaM ₃ . These results suggested that mutations within the S2–S3 loop of Kv7. 4 channels decreased the calcium-dependent inhibition mediated by functional CaM ( Supplementary Figure S7 ) and proved that the S2–S3 loop of Kv7. 4 is required for CaM regulation of Kv7. 4 activation. We could not probe the role of the residues Tyr160 and Gly162 in CaM regulation, because the mutants Y160A and G162V did not yield functional currents ( Supplementary Figures S6B , S8A, B ). Moreover, the mutations W163A and Q164A did not affect the currents of Kv7. 4 co-expressed with WT CaM ( Supplementary Figures S8C , E ) or CaM ₃ ( Supplementary Figures S8D , F ), implying that the two residues were not involved in CaM regulation.

FIGURE 4

Mutations in the Kv7. 4 S2-S3 loop affect CaM regulation of channel activation.(A)A schematic depicting a possible structural model of Ca ²⁺/CaM interacting with the Kv7. 4 subunit, in which both the Ca ²⁺/N-lobe and Ca ²⁺/C-lobe interact with the B helix of Kv7. 4 channels, based on structural studies of Ca ²⁺/CaM with the Kv7. 4 A and B helices ( Xu et al., 2013 ; Chang et al., 2018 ). A magnified view of the potential interface between the Kv7. 4 S2–S3 loop and the Ca ²⁺/CaM EF3 is shown in the dashed box. The nine residues of the Kv7. 4 S2–S3 loop studied in this work are labeled with circles. Red circles indicate residues required for CaM regulation of Kv7. 4 activation (according to our electrophysiological).(B)Representative current traces of cells co-expressing Kv7. 4 WT or Kv7. 4 mutants with CaM.(C)Representative current traces of cells co-expressing Kv7. 4 WT or Kv7. 4 mutants with CaM ₃ .(D)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ( N = 12), Kv7. 4 C156A: CaM ( N = 13), Kv7. 4 C157A: CaM ( N = 8), Kv7. 4 C158V: CaM ( N = 10), Kv7. 4 R159A: CaM ( N = 9), and Kv7. 4 R161A: CaM ( N = 13).(E)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ₃ ( N = 11), Kv7. 4 C156A: CaM ₃ ( N = 11), Kv7. 4 C157A: CaM ₃ ( N = 10), Kv7. 4 C158V: CaM ₃ ( N = 9), Kv7. 4 R159A: CaM ₃ ( N = 8), and Kv7. 4 R161A: CaM ₃ ( N = 9).(F)Voltage-dependent activation curves for Kv7. 4: CaM ( N = 12), Kv7. 4 C156A: CaM ( N = 13), Kv7. 4 C157A: CaM ( N = 8), Kv7. 4 C158V: CaM ( N = 10), Kv7. 4 R159A: CaM ( N = 9), and Kv7. 4 R161A: CaM ( N = 13).(G)Voltage-dependent activation curves for Kv7. 4: CaM ₃ ( N = 11), Kv7. 4 C156A: CaM ₃ ( N = 11), Kv7. 4 C157A: CaM ₃ ( N = 10), Kv7. 4 C158V: CaM ₃ ( N = 9), Kv7. 4 R159A: CaM ₃ ( N = 8), and Kv7. 4 R161A: CaM ₃ ( N = 9). Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^*** p < 0. 001. NS indicates not significant.

The Kv7. 4 Double Mutation C156A/R159A Completely Abolishes CaM Regulation of Kv7. 4 Activation

Although the five Kv7. 4 mutations (C156A, C157A, C158V, R159A, and R161A) facilitated activation of Kv7. 4 channels co-expressed with WT CaM, the activation time constants and the half-activation voltages of the five single mutants were still different from the parameters of Kv7. 4 channels co-expressed with CaM ₃ , suggesting that the inhibitory effects of Ca ²⁺/CaM was decreased but not abolished completely. For further investigation on the necessity of the S2–S3 loop, we sought to enhance the facilitating effects by generating double mutants (C156A/R159A, C156A/R161A, C158V/R159A, and R159A/R158A). Shown in Figure 5A are whole-cell current traces of the double mutants co-expressed with WT CaM or CaM ₃ . When co-expressed with WT CaM, all the double mutants displayed further facilitated channel activation in terms of activation kinetics ( Figure 5B ; τ values at 40 mV: Kv7. 4: CaM = 71. 9 ± 5. 0 ms, Kv7. 4 C156A/R159A: CaM = 19. 1 ± 1. 7 ms, Kv7. 4 C156A/R161A: CaM = 31. 7 ± 3. 1 ms, Kv7. 4 C158V/R159A: CaM = 23. 9 ± 2. 0 ms, and Kv7. 4 R159A/R161A: CaM = 25. 5 ± 1. 9 ms) and voltage-dependent activation ( Figures 5C , D ; V _1/2 values: Kv7. 4: CaM = −21. 1 ± 1. 2 mV, Kv7. 4 C156A/R159A = −47. 3 ± 1. 2 mV, Kv7. 4 C156A/R161A: CaM = −43. 3 ± 0. 7 mV, Kv7. 4 C158V/R159A: CaM = −41. 3 ± 0. 8 mV, and Kv7. 4 C158V/R159A: CaM = −41. 7 ± 1. 1 mV). As expected, all the double mutants co-expressed with CaM ₃ generated activation rates ( Figure 4B ; τ values at 40 mV: Kv7. 4: CaM ₃ = 19. 2 ± 1. 4 ms, Kv7. 4 C156A/R159A: CaM ₃ = 20. 5 ± 1. 6 ms, Kv7. 4 C156A/R161A: CaM ₃ = 20. 8 ± 1. 7 ms, Kv7. 4 C158V/R159A: CaM ₃ = 18. 6 ± 1. 1 ms, and Kv7. 4 R159A/R161A: CaM ₃ = 19. 9 ± 1. 6 ms) and half-activation voltages ( Figures 4C , D ; V _1/2 values: Kv7. 4: CaM ₃ = −48. 6 ± 0. 9 mV, Kv7. 4 C156A/R159A: CaM ₃ = −49. 9 ± 0. 8 mV, Kv7. 4 C156A/R161A: CaM ₃ = −50. 0 ± 1. 8 mV, Kv7. 4 C158V/R159A: CaM ₃ = −49. 0 ± 0. 9 mV, and Kv7. 4 R159A/R161A: CaM ₃ = −48. 8 ± 1. 8 mV; Figures 4C , D ) similar to WT Kv7. 4 channels with CaM ₃ . Importantly, according to our statistical results, there was no significant difference among Kv7. 4 C156A/R159A: CaM, Kv7. 4 C156A/R159A: CaM ₃ , and Kv7. 4: CaM ₃ in both activation kinetics ( Figure 5B ) and voltage-dependent activation ( Figure 5C ). This indicated that the double mutation C156A/R159A generated facilitating effects equivalent to CaM ₃ and suggested that CaM regulation of Kv7. 4 activation was abolished completely by this mutation. Additionally, both varying free calcium concentrations of the pipette solutions and application of BAPTA-AM have no effect on activation kinetics ( Figure 5E ; τ values at 40 mV: Kv7. 4 C156A/R159A: CaM pretreated with 20 μM BAPTA-AM = 21. 6 ± 2. 1 ms; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA = 20. 2 ± 1. 9 ms; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA plus 3 mM CaCl2 = 18. 5 ± 3. 0 ms; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA plus 6 mM CaCl2 = 19. 0 ± 2. 3 ms) or voltage-dependent activation ( Figure 5F ; V1/2 values: Kv7. 4 C156A/R159A: CaM pretreated with 20 μM BAPTA-AM = −48. 6 ± 1. 7 mV; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA = −47. 2 ± 1. 3 mV; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA plus 3 mM CaCl2 = −50. 9 ± 1. 4 mV; Kv7. 4 C156A/R159A: CaM with 10 mM EGTA plus 6 mM CaCl2 = −48. 6 ± 1. 1 mV) of Kv7. 4 C156A/R159A: CaM. Taken together, these results of Kv7. 4 C156A/R159A further emphasized the essential role of the Kv7. 4 S2–S3 loop in CaM-mediated calcium-dependent regulation.

FIGURE 5

The Kv7. 4 double mutation C156A/R159A completely abolishes CaM regulation of channel activation.(A)Representative current traces of cells co-expressing Kv7. 4 double mutants with CaM (top) or CaM ₃ (bottom).(B)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ( N = 17), Kv7. 4: CaM ₃ ( N = 13), Kv7. 4 C156A/R159A: CaM ( N = 9), Kv7. 4 C156A/R159A: CaM ₃ ( N = 9), Kv7. 4 C156A/R161A: CaM ( N = 17), Kv7. 4 C156A/R161A: CaM ₃ ( N = 10), Kv7. 4 C158V/R159A: CaM ( N = 10), Kv7. 4 C158V/R159A: CaM ₃ ( N = 11), Kv7. 4 R159A/R161A: CaM ( N = 11), and Kv7. 4 R159A/R161A: CaM ₃ ( N = 10).(C)Individual V _1/2 values for Kv7. 4: CaM ( N = 17), Kv7. 4: CaM ₃ ( N = 13), Kv7. 4 C156A/R159A: CaM ( N = 9), Kv7. 4 C156A/R159A: CaM ₃ ( N = 9), Kv7. 4 C156A/R161A: CaM ( N = 17), Kv7. 4 C156A/R161A: CaM ₃ ( N = 10), Kv7. 4 C158V/R159A: CaM ( N = 10), Kv7. 4 C158V/R159A: CaM ₃ ( N = 11), Kv7. 4 R159A/R161A: CaM ( N = 11), and Kv7. 4 R159A/R161A: CaM ₃ ( N = 10).(D)Voltage-dependent activation curves for Kv7. 4: CaM ( N = 17), Kv7. 4: CaM ₃ ( N = 13), Kv7. 4 C156A/R159A: CaM ( N = 9), Kv7. 4 C156A/R159A: CaM ₃ ( N = 9), Kv7. 4 C156A/R161A: CaM ( N = 17), Kv7. 4 C156A/R161A: CaM ₃ ( N = 10), Kv7. 4 C158V/R159A: CaM ( N = 10), Kv7. 4 C158V/R159A: CaM ₃ ( N = 11), Kv7. 4 R159A/R161A: CaM ( N = 11), and Kv7. 4 R159A/R161A: CaM ₃ ( N = 10).(E)Average time constants at 40 mV for the fast activation component of Kv7. 4 C156A/R159A: CaM with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4 C156A/R159A: CaM pretreated with BAPTA-AM ( N = 8).(F)Voltage-dependent activation curves for Kv7. 4 C156A/R159A: CaM with pipette solutions containing only 10 mM EGTA ( N = 7), 10 mM EGTA plus 3 mM CaCl ₂ ( N = 7), or 10 mM EGTA plus 6 mM CaCl ₂ ( N = 7) and Kv7. 4 C156A/R159A: CaM pretreated with BAPTA-AM ( N = 8). N indicates the number of cells. Asterisks indicate significance:^* p < 0. 05;^** p < 0. 01;^*** p < 0. 001. NS indicates not significant.

We used co-immunoprecipitation experiments to probe the effects of the double mutation C156/R159A, which completely abolished CaM regulation of Kv7. 4 activation, on Kv7. 4 forming complex with CaM, a requirement for their interaction. A Flag tag and an HA tag were introduced to the N-terminus of Kv7. 4 and the C-terminus of CaM, respectively. Electrophysiological recordings revealed no noticeable distinction between tagged and untagged combinations of Kv7. 4 with CaM in channel currents ( Supplementary Figure S3A ) or in the effects of CaM regulation ( Supplementary Figures S3B , C ). Immunoprecipitation accomplished using anti-Flag antibodies to precipitate Flag-Kv7. 4 or Flag-Kv7. 4 C156A/R159A in the presence of either 2 mM CaCl2 or 2 mM EGTA. The immunoprecipitates were probed with anti-Flag antibodies to confirm precipitation of Flag-Kv7. 4 or Flag-Kv7. 4 C156A/R159A and were probed with anti-HA antibodies to detect CaM-HA or CaM ₃ -HA pulled down with the channels. In the presence of 2 mM CaCl ₂ , WT CaM showed a detectable but much weaker association with the Kv7. 4 C156A/R159A mutant than with WT Kv7. 4, while the signals of CaM ₃ pulled down with WT Kv7. 4 and Kv7. 4 C156A/R159A are similar ( Figure 6 ). These results implied that the Kv7. 4 S2-S3 loop interacted with Ca ²⁺/CaM but had no association with Ca ²⁺/CaM ₃ that had an empty EF3 and calcium-bound EF hands 1, 2, and 4. Of note, the signal of CaM ₃ pulled down with Kv7. 4 channels was stronger than WT CaM in the presence of 2 mM CaCl ₂ , matching the observation by Kosenko and Hoshi (2013), who reported that longer exposure to calcium decreased CaM binding to Kv7. 2 channels. This difference might result from the different interactions of the Kv7. 4 C-terminal domain with Ca ²⁺/CaM and with CaM ₃ , in which calcium binding disturbed the interaction between the CaM C-lobe and the Kv7. 4 A helix and the Ca ²⁺/C-lobe bound to the B helix weakly through a much smaller interaction surface ( Xu et al., 2013 ; Chang et al., 2018 ). Therefore, the stronger binding of CaM ₃ than Ca ²⁺/CaM is not contradictory to our hypothesis that Ca ²⁺/CaM interacted with the Kv7. 4 S2–S3 loop, whereas Ca ²⁺/CaM ₃ lost this interaction. In the presence of 2 mM EGTA through the whole immunoprecipitation overnight, both WT CaM and CaM ₃ would exist in Apo/CaM forms, and there were no significant differences among the signals of WT CaM pulled down with WT Kv7. 4, WT CaM pulled down with Kv7. 4 C156A/R159A, CaM ₃ pulled down with WT Kv7. 4, and CaM ₃ pulled down with Kv7. 4 C156A/R159A ( Figure 6 ). Additionally, the binding of Apo/CaM is stronger than Ca ²⁺/CaM and similar to Ca ²⁺/CaM ₃ , suggesting an important role of the EF3 in CaM binding. Taken together, these results of the co-immunoprecipitation experiments coincided with our electrophysiological data, supporting an potential correlation between the structure and function of the Kv7. 4: CaM complex.

FIGURE 6

The Kv7. 4 double mutation C156A/R159A decreased Ca ²⁺/CaM binding.(A)Co-immunoprecipitation analysis of CaM or CaM ₃ binding to Kv7. 4 or Kv7. 4 C156A/R159A. Chinese hamster ovary (CHO) cells were co-transfected with Flag-Kv7. 4 and CaM-HA, Flag Kv7. 4 C156A/R159A and CaM-HA, Flag-Kv7. 4 and CaM ₃ -HA, or Flag-Kv7. 4 C156A/R159A: CaM ₃ -HA. Immunoprecipitation was accomplished using anti-Flag antibodies in the presence of either 2 mM CaCl ₂ or 2 mM EGTA. The immunoprecipitates were probed with anti-Flag antibodies to confirm precipitation of Flag-Kv7. 4 or Flag-Kv7. 4 C156A/R159A and were probed with anti-HA antibodies to detect CaM-HA or CaM ₃ -HA pulled down with the channels. The experiments were replicated four times and showed similar results.(B)Quantification of CaM or CaM ₃ binding to Kv7. 4 or Kv7. 4 C156A/R159A. The binding efficiency is represented as the ratio of HA intensity to Flag intensity, and all ratios are normalized to the ratio representing CaM binding to Kv7. 4 in the presence of CaCl ₂ . Asterisks indicate significance:^* p < 0. 05;^** p < 0. 01. NS indicates not significant.

The Residues N98 and N99 in the CaM EF3 Are Required for the Regulation of Kv7. 4 Activation

Sun and MacKinnon (2017) have previously shown that four residues in the CaM EF3 (Gly97, Asn98, Gly99, and Tyr100) are involved in the newly discovered interface between CaM and the S2–S3 loop of Kv7. 1 channels ( Supplementary Figure S9 ). We generated CaM mutants in which Gly97, Asn98, Gly99, or Tyr100 were mutated to an alanine based on WT CaM (CaM G97A, CaM N98A, CaM G99A, and CaM Y100A) or CaM ₃ (CaM ₃ G97A, CaM ₃ N98A, CaM ₃ G99A, and CaM ₃ Y100A). Kv7. 4 channels co-expressed with all these CaM mutants yielded functional currents ( Figures 7A , B ). Compared with WT CaM, CaM N98A and CaM G99A accelerated the activation rate ( Figure 7C ; τ values at 40 mV: Kv7. 4: CaM = 73. 1 ± 4. 8 ms, Kv7. 4: CaM N98A = 29. 1 ± 2. 5 ms, and Kv7. 4: CaM G99A = 35. 2 ± 2. 3 ms) and shifted the voltage-dependent activation curve leftwards ( Figure 7E ; V _1/2 values: Kv7. 4: CaM = −18. 7 ± 1. 0 mV, Kv7. 4: CaM N98A = −32. 8 ± 0. 9 ms, and Kv7. 4: CaM G99A = −31. 9 ± 1. 8 ms). Kv7. 4 co-expressed with CaM G97A or CaM Y100A generated activation properties identical to Kv7. 4 co-expressed with WT CaM in terms of both the activation rate ( Figure 7C ; τ values at 40 mV: Kv7. 4: CaM G97A = 68. 6 ± 3. 1 ms and Kv7. 4: CaM Y100A = 70. 0 ± 7. 0 ms) and the voltage dependence of activation ( Figure 7E ; V _1/2 values: Kv7. 4: CaM N98A = −20. 4 ± 1. 0 ms and Kv7. 4: CaM G99A = −21. 6 ± 0. 8 ms), suggesting that the two residues of CaM were not involved in the regulation of Kv7. 4 activation. All mutants based on CaM ₃ (CaM ₃ G97A, CaM ₃ N98A, CaM ₃ G99A, and CaM ₃ Y100A) phenocopied the effects of CaM ₃ on the activation rate ( Figure 7D ; τ values at 40 mV: Kv7. 4: CaM ₃ = 19. 2 ± 1. 1 ms, Kv7. 4: CaM ₃ G97A = 19. 5 ± 1. 4 ms, Kv7. 4: CaM ₃ N98A = 21. 3 ± 2. 0 ms, Kv7. 4: CaM ₃ G99A = 20. 7 ± 1. 3 ms, and Kv7. 4: CaM ₃ Y100A = 22. 1 ± 1. 1 ms) and the voltage dependence of activation ( Figure 7F ; V _1/2 values: Kv7. 4: CaM ₃ = −48. 4 ± 1. 0 mV, Kv7. 4: CaM ₃ G97A = −46. 3 ± 1. 2 mV, Kv7. 4: CaM ₃ N98A = −45. 4 ± 1. 3 mV, Kv7. 4: CaM ₃ G99A = −48. 1 ± 1. 1 mV, and Kv7. 4: CaM ₃ Y100A = −48. 0 ± 0. 9 mV). Taken together, the effects of the two mutations N98A and G99A in the CaM EF3 were similar to the effects of the five mutations (C156A, C157A, C158V, R159A, and R161A) in the Kv7. 4 S2–S3 loop, which potentially interacted with the CaM EF3, providing further supporting evidence for the model in which CaM modulated Kv7. 4 activation through the S2-S3 loop.

FIGURE 7

CaM mutations N98A and G99A affect the regulation of Kv7. 4 activation.(A)Representative current traces of cells co-expressing Kv7. 4 with CaM-based mutants.(B)Representative current traces of cells co-expressing Kv7. 4 with CaM ₃ -based mutants.(C)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ( N = 13), Kv7. 4: CaM G97A ( N = 15), Kv7. 4: CaM N98A ( N = 12), Kv7. 4: CaM G99A ( N = 10), and Kv7. 4: CaM Y100A ( N = 17).(D)Average time constants at 40 mV for the fast activation component of Kv7. 4: CaM ₃ ( N = 11), Kv7. 4: CaM ₃ G97A ( N = 11), Kv7. 4: CaM ₃ N98A ( N = 11), Kv7. 4: CaM ₃ G99A ( N = 14), and Kv7. 4: CaM ₃ Y100A ( N = 9).(E)Voltage-dependent activation curves for Kv7. 4: CaM ( N = 13), Kv7. 4: CaM G97A ( N = 15), Kv7. 4: CaM N98A ( N = 12), Kv7. 4: CaM G99A ( N = 10), and Kv7. 4: CaM Y100A ( N = 17).(F)Voltage-dependent activation curves for Kv7. 4: CaM ₃ ( N = 11), Kv7. 4: CaM ₃ G97A ( N = 11), Kv7. 4: CaM ₃ N98A ( N = 11), Kv7. 4: CaM ₃ G99A ( N = 14), and Kv7. 4: CaM ₃ Y100A ( N = 9). Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^*** p < 0. 001. NS indicates not significant.

CaM Regulates Activation of Kv7. 4 Isoform b

There are two alternatively spliced isoforms of human Kv7. 4, called Kv7. 4 isoform a and Kv7. 4 isoform b. Most studies including ours, used Kv7. 4 isoform a, which we and others called Kv7. 4. To further the understanding of CaM regulation of both the two isoforms, we also performed experiments on Kv7. 4 isoform b, which we called Kv7. 4b to distinguish from isoform a. Whole-cell recordings were performed on CHO cells transfected with Kv7. 4b alone or together with WT CaM, CaM ₁₂₃₄ , or CaM ₃ . Kv7. 4b with only endogenous CaM and Kv7. 4b with overexpression of exogenous WT CaM showed the same activation rate ( τ values at 40 mV: Kv7. 4b = 52. 5 ± 5. 0 ms; Kv7. 4b: CaM = 48. 0 ± 4. 2 ms) and voltage dependence of activation (V _1/2 values: Kv7. 4b = −22. 2 ± 2. 0 mV; Kv7. 4b: CaM = −24. 9 ± 1. 0 mV). Overexpression of CaM ₁₂₃₄ significantly accelerated the activation rate ( τ value at 40 mV: Kv7. 4b: CaM ₁₂₃₄ = 23. 4 ± 2. 7 ms) and shifted the voltage dependent activation leftwards markedly (V _1/2 value: Kv7. 4b: CaM ₁₂₃₄ = −48. 8 ± 1. 1 mV). CaM ₃ mirrored the facilitating effects of CaM ₁₂₃₄ on both activation kinetics ( τ value at 40 mV: Kv7. 4b: CaM ₃ = 21. 8 ± 2. 2 ms) and voltage-dependence activation (V _1/2 value: Kv7. 4b: CaM ₃ = −51. 6 ± 1. 1 mV). To probe effects of the Kv7. 4b S2–S3 loop on CaM regulation, we also performed whole-cell recordings on CHO cells transfected with Kv7. 4b C156A/R159A alone or together with WT CaM, CaM ₁₂₃₄ , or CaM ₃ . According to the statistical results, there was no significant difference among the four groups in both the activation rate ( τ values at 40 mV: Kv7. 4b C156A/R156A = 23. 5 ± 3. 4 ms; Kv7. 4b C156A/R159A: CaM = 22. 5 ± 2. 0 ms; Kv7. 4b C156A/R159A: CaM ₁₂₃₄ = 23. 6 ± 1. 8 ms; Kv7. 4b C156A/R159A: CaM ₃ = 21. 6 ± 1. 9 ms) and the voltage dependence of activation (V _1/2 values: Kv7. 4b C156A/R156A = −53. 2 ± 2. 0 mV; Kv7. 4b C156A/R159A: CaM = −53. 5 ± 0. 9 mV; Kv7. 4b C156A/R159A: CaM ₁₂₃₄ = −52. 4 ± 1. 3 mV; Kv7. 4b C156A/R159A: CaM ₃ = −53. 3 ± 0. 9 mV), suggesting that the double mutation C156A/R159A in the S2–S3 loop also abolished CaM regulation of Kv7. 4b activation. Taken together, our work revealed a unified mechanism for CaM regulation of both isoforms of Kv7. 4 channels ( Figure 8 ).

FIGURE 8

CaM regulates activation of Kv7. 4 isoform b.(A)Representative current traces of Kv7. 4b without overexpression of exogenous CaM and Kv7. 4b with overexpression of exogenous CaM, CaM ₁₂₃₄ , or CaM ₃ .(B)Representative current traces of Kv7. 4b C156A/R159A without overexpression of exogenous CaM and Kv7. 4b C156A/R159A with overexpression of exogenous CaM, CaM ₁₂₃₄ , or CaM ₃ .(C)Corresponding current density-voltage relations of Kv7. 4b ( N = 7), Kv7. 4b: CaM ( N = 10), Kv7. 4b: CaM ₁₂₃₄ ( N = 9), and Kv7. 4b: CaM ₃ ( N = 10).(D)Corresponding current density-voltage relations of Kv7. 4b C156A/R159A ( N = 7), Kv7. 4b C156A/R159A: CaM ( N = 9), Kv7. 4b C156A/R159A: CaM ₁₂₃₄ ( N = 8), and Kv7. 4b C156A/R159A: CaM ₃ ( N = 8).(E)Average time constants at 40 mV for the fast activation component of Kv7. 4b ( N = 7), Kv7. 4b: CaM ( N = 10), Kv7. 4b: CaM ₁₂₃₄ ( N = 9), Kv7. 4b: CaM ₃ ( N = 10), Kv7. 4b C156A/R159A ( N = 7), Kv7. 4b C156A/R159A: CaM ( N = 9), Kv7. 4b C156A/R159A: CaM ₁₂₃₄ ( N = 8), and Kv7. 4b C156A/R159A: CaM ₃ ( N = 8).(F)Voltage-dependent activation curves for Kv7. 4b ( N = 7), Kv7. 4b: CaM ( N = 10), Kv7. 4b: CaM ₁₂₃₄ ( N = 9), and Kv7. 4b: CaM ₃ ( N = 10).(G)Voltage-dependent activation curves for Kv7. 4b C156A/R159A ( N = 7), Kv7. 4b C156A/R159A: CaM ( N = 9), Kv7. 4b C156A/R159A: CaM ₁₂₃₄ ( N = 8), and Kv7. 4b C156A/R159A: CaM ₃ ( N = 8). Error bars show SEM. N indicates the number of cells. Asterisks indicate significance:^*** p < 0. 001. NS indicates not significant.

Discussion

With the properties of opening at subthreshold membrane potentials and lacking inactivation, Kv7 channels hold a pivotal position in controlling excitability in the heart and the nervous system ( Barhanin et al., 1996 ; Sanguinetti et al., 1996 ; Wang et al., 1996 , 1998 ; Jentsch, 2000 ; Delmas and Brown, 2005 ; Hernandez et al., 2008 ; Brown and Passmore, 2009 ). CaM conferred effective regulation to Kv7 channels. Though the importance of CaM regulation of Kv7 channels has been recognized, the regulation mechanism remains elusive. Earlier studies mostly focused on the role of the C-terminal domain of Kv7 channels in CaM binding and CaM regulation. Mutations in the C-terminal domain impairing CaM binding affect channel assembly, trafficking, and gating ( Ghosh et al., 2006 ; Shamgar et al., 2006 ; Etxeberria et al., 2008 ; Alaimo et al., 2009 , 2018 ; Liu and Devaux, 2014 ; Tobelaim et al., 2017 ; Chang et al., 2018 ). Here, our data provided new insight into CaM regulation of Kv7 channels by proving that the S2–S3 loop of Kv7. 4 channels is essential for CaM-mediated calcium-dependent regulation of Kv7. 4 activation.

Previous studies have shown that CaM conferred intracellular calcium signals to Kv7 channels ( Gamper and Shapiro, 2003 ; Gamper et al., 2005 ; Sihn et al., 2016 ; Gomis-Perez et al., 2017 ; Tobelaim et al., 2017 ; Chang et al., 2018 ). To probe the calcium-dependent regulation of Kv7. 4 channels mediated by CaM, we performed electrophysiology experiments on CHO cells transfected with Kv7. 4 anole or together with WT CaM or CaM mutants with impaired calcium binding abilities. CHO cells have endogenous CaM, which is important for many cellular physiological activities, and thus cannot be knocked out. Overexpression of exogenous CaM or CaM mutants could knockdown the effects of endogenous CaM, which is a common method used by most of the previous studies to probe calcium-dependent CaM regulation. Compared with Kv7. 4 channels with only endogenous CaM of CHO cells, overexpression of exogenous WT CaM increased current amplitudes over 2-fold but did not change activation properties, and overexpression of exogenous CaM ₁₂₃₄ increased current amplitudes about 5-fold and facilitated channel activation as well. CaM ₁₂₃₄ sharply accelerated the activation rate and produced a large leftward shift in the voltage-dependent activation, similar to the effects of CaM ₁₂₃₄ on Kv7. 2 ( Gomis-Perez et al., 2017 ), Kv7. 3 ( Gomis-Perez et al., 2017 ), and Kv7. 5 ( Chang et al., 2018 ) and implying a unified mechanism in which Apo/CaM facilitates Kv7. 2–7. 5 activation, or in other words, in which Ca ²⁺/CaM inhibits Kv7. 2–7. 5 activation.

We also measured surface expression of Kv7. 4 with only endogenous CaM and Kv7. 4 with overexpression of exogenous WT CaM and CaM ₁₂₃₄ to further explain the different current amplitudes. Kv7. 4 co-expressed with WT CaM and Kv7. 4 co-expressed with CaM ₃ showed similar surface expression and about 3. 5-fold the surface expression of Kv7. 4 with only endogenous CaM, suggesting that CaM enhances Kv7. 4 surface expression in a calcium-independent manner. We also measured Kv7. 4 total expression of the three groups, and the total expression of Kv7. 4 with overexpression of exogenous WT CaM or CaM ₁₂₃₄ was about twice the total expression of Kv7. 4 channels with only endogenous CaM ( Supplementary Figures S4A , B ). Taken together, increased surface expression accounted for the greater current amplitudes of Kv7. 4 with overexpression of exogenous WT CaM than Kv7. 4 with only endogenous CaM, and Kv7. 4 with overexpression of exogenous CaM ₁₂₃₄ generated further greater current amplitudes due to both facilitated activation and increased surface expression. As the facilitating effect of the overexpression of exogenous WT CaM or CaM ₁₂₃₄ on Kv7. 4 surface expression is more significant than the facilitating effect on Kv7. 4 total expression, we deem that the increase of Kv7. 4 surface expression is likely due to not only the increase of Kv7. 4 total expression but also enhanced channel trafficking.

Of note, Kv7. 4 with overexpression of exogenous WT CaM showed activation properties equivalent to Kv7. 4 with only endogenous CaM, possibly because exogenous WT CaM is essentially the same as endogenous CaM except the much higher expression level. CaM was reported to be an auxiliary subunit of Kv7 channels ( Wen and Levitan, 2002 ). Disrupted CaM interactions with Kv7 channels have been shown to be tightly associated with severely impaired Kv7 surface expression ( Ghosh et al., 2006 ; Shamgar et al., 2006 ; Etxeberria et al., 2008 ; Alaimo et al., 2009 , 2018 ; Liu and Devaux, 2014 ; Chang et al., 2018 ), illustrating the vital role of CaM in Kv7 assembly and trafficking, which is also supported by our study. Additionally, the stoichiometry for CaM binding to the Kv7 subunit was 1: 1 ( Sun and MacKinnon, 2017 ), and Bernardo-Seisdedos et al. (2018) reported that the Kv7. 2 A and B helices could not be purified in the absence of CaM and that the complex could not be dissociated at any calcium concentration in vitro . Therefore, we deem that there would be CaM strongly and stably binding to the functional plasma membrane Kv7. 4 channels, which has been properly assembled and transported, at a stoichiometry of 1: 1, although more Kv7. 4 channels was expressed on the cell surface due to the overexpression of exogenous WT CaM. However, the already bound CaM controls the regulating activation of the plasma membrane Kv7. 4 channels, so the different expression levels of endogenous CaM and overexpressed exogenous WT CaM did not change Kv7. 4 activation.

We also studied the effects CaM ₁₂ , CaM ₃₄ , CaM ₃ , and CaM ₄ on Kv7. 4 activation to probe the role of each lobe or EF hand in regulating Kv7. 4 channels. CaM ₃₄ and CaM ₃ caused facilitating effects equivalent to CaM ₁₂₃₄ on Kv7. 4 activation kinetics and the voltage dependence of activation, supporting a key role of the EF3 in CaM regulation of Kv7. 4 activation. CaM ₁₂ did not mirror the effects of WT CaM of Kv7. 4 activation, possibly because calcium binding could enhance the binding of the N-lobe to the B helix of Kv7. 4 ( Chang et al., 2018 ), and such configuration changes might have slight effects on channel function. However, the changes of Kv7. 4 activation caused by CaM ₁₂ were much slighter than the changes caused by CaM ₁₂₃₄ , CaM ₃₄ , and CaM ₃ . Our data largely match the results by Chang et al. (2018), but there is a difference in the effects of CaM ₄ on Kv7. 4 activation. We found that CaM ₄ did not change Kv7. 4 activation kinetics or the voltage dependence of activation, whereas Chang et al. (2018) showed that CaM ₄ facilitated Kv7. 4 activation though the facilitating effects is much weaker than CaM ₃ , CaM ₃₄ , and CaM ₁₂₃₄ . This difference might result from the different experimental approaches, in which Chang et al. (2018) used perforated-patch configuration to test cells under natural resting calcium levels but we used whole-cell configuration to substitute the cell cytosolic composition by the pipette solution. Moreover, activation of Kv7. 4 co-expressed with WT CaM was inhibited by increasing intracellular calcium levels, while CaM ₁₂₃₄ and CaM ₃ rendered Kv7. 4 insensitive to the changing calcium concentrations, further emphasizing the role of calcium in CaM regulation of Kv7. 4 channels and implying that CaM ₁₂₃₄ and CaM ₃ facilitate Kv7. 4 activation essentially by abolishing the calcium-dependent inhibition.

The free calcium concentration in the regular pipette solution is 80 nM, which is within the basal cytosolic calcium concentrations (10–100 nM). Free CaM exists in its Apo-forms in the basal cytosolic calcium concentrations [K _d (C-lobe) = 3. 4 μM, K _d (N-lobe) = 14 μM; Linse et al., 1991 ]. However, CaM binding to some target peptides or proteins can dramatically increase calcium affinities ( Olwin et al., 1984 ; Olwin and Storm, 1985 ; Peersen et al., 1997 ; Black et al., 2005 ; Evans and Shea, 2009 ; Findeisen et al., 2013 ). For instance, the calcium affinities of CaM binding to the IQ motif of Cav channels, kinase skMLCK or kinase MKII are quite high with Kd values much lower than 100 nM ( Peersen et al., 1997 ; Black et al., 2005 ). In this study, the calcium chelators, EGTA and BAPTA-AM exerted effects on Kv7. 4 activation in the same facilitating direction as CaM ₁₂₃₄ , CaM ₁₂ , CaM ₃₄ , and CaM ₃ , suggesting that there are some tonic occupations by calcium in the EF hands under the basal intracellular calcium level. Previous studies have also reported changed Kv7 currents by decreasing intracellular calcium concentrations from the basal level ( Ghosh et al., 2006 ; Shamgar et al., 2006 ; Chang et al., 2018 ). These observations are possibly due to increased calcium affinities of CaM when binding to Kv7 channels.

Structural studies have reported that the CaM EF3 interacted with the Kv7. 1 S2–S3 loop, which is conserved among Kv7 isoforms ( Sun and MacKinnon, 2017 ). To explore if the important functional role of the CaM EF3 is associated with the potential interaction between the CaM EF3 and the Kv7. 4 S2–S3, we mutated amino acids either in the Kv7. 4 S2–S3 loop or in the CaM EF3 to a nonpolar amino acid. Five mutations (C156A, C157A, C158V, R159A, and R161A) in the Kv7. 4 S2–S3 loop and two mutations (N98A and G99A) in the CaM EF3 led to decreased inhibitory effects of Ca ²⁺/CaM, facilitating Kv7. 4: CaM activation in terms of both the activation rate and the voltage dependence of activation. The double mutation C156A/R159A within the S2–S3 loop completely abolished the inhibitory effects of Ca ²⁺/CaM on Kv7. 4 activation and rendered Kv7. 4 insensitive to changing intracellular calcium levels, which further emphasizes the importance of the Kv7. 4 S2–S3 loop in CaM regulation. We also performed co-immunoprecipitation experiments to probe the effect of C156A/R159A on CaM binding. Compared with WT Kv7. 4, the double mutant C156A/R159A showed weaker association with Ca ²⁺/CaM but unchanged association with CaM ₃ and Apo/CaM, suggesting that the Kv7. 4 S2–S3 loop interacts with the Ca ²⁺/EF3 but has no association with the Apo/EF3. The co-immunoprecipitation results coincide with our electrophysiological results. Mutations in the Kv7. 4 S2–S3 loop decreased Ca ²⁺/CaM binding as well as Ca ²⁺/CaM inhibition on channel activation, while these mutations had no effect on CaM ₃ binding as well as activation of Kv7. 4 co-expressed with CaM ₃ , suggesting that Ca ²⁺/CaM inhibited Kv7. 4 activation through an interaction with the Kv7. 4 S2–S3 loop.

We also probed CaM regulation of Kv7. 4 isoform b, activation of which was reported by Sihn et al. (2016) to be not regulated by CaM. Our results are strikingly different from Sihn et al. (2016) and show that CaM regulates Kv7. 4b activation in a calcium-dependent manner similarly to the regulation of Kv7. 4 isoform a. CaM ₁₂₃₄ facilitated Kv7. 4b activation in terms of both activation kinetics and the voltage dependence of activation, and CaM ₃ phenocopied the facilitating effects of CaM ₁₂₃₄ . The double mutation C156A/R159A abolished CaM regulation of Kv7. 4b activation. These results emphasized the crucial roles of the EF3 and the S2–S3 loop in CaM regulating activation of Kv7. 4b. There are cases where the earlier studies failed to find CaM regulation of Kv7. 1 and Kv7. 3 ( Gamper et al., 2005 ), both of which were later shown to be regulated by CaM ( Sachyani et al., 2014 ; Gomis-Perez et al., 2017 ; Tobelaim et al., 2017 ; Chang et al., 2018 ). Our work proved a unified rather than isoform-dependent mechanism for CaM regulation of Kv7. 4 channels.

To the best of our knowledge, the mechanism for CaM-mediated calcium-dependent regulation of voltage-gated channels has remained puzzling since the CaM binding sites were initially found. Heretofore, the most-studied structural presentation accounting for the mechanism for CaM regulation of channel activation comes from SK4 channels ( Lee and MacKinnon, 2018 ). It was proposed that Apo/CaM binds to SK4 channels through the interaction of Apo/C-lobe with the SK4 A and B helices, and calcium loading allows the Ca ²⁺/N-lobe to remain anchored on the A and B helices and enables the Ca ²⁺/N-lobe to contact the S4–S5 linker of SK4 channels and further affect pore opening. In Kv7 channels, it has been widely assumed that the Ca ²⁺-dependent CaM regulation is also driven by conformational changes derived from calcium loading, but detailed interpretation remained unresolved. Apo/CaM binds to Kv7 channels in a clamped form, in which the C-lobe interacts with the A helix of Kv7 channels and the N-lobe interacts with the B helix ( Bernardo-Seisdedos et al., 2018 ; Chang et al., 2018 ). Calcium binding allow the Ca ²⁺/N-lobe to stay anchored on the Kv7 B helix in a similar configuration to the apo/N-lobe ( Sachyani et al., 2014 ; Strulovich et al., 2016 ; Sun and MacKinnon, 2017 ; Bernardo-Seisdedos et al., 2018 ; Chang et al., 2018 ), while there are striking differences between the Ca ²⁺/C-lobe and Apo/C-lobe conformations. Structural studies of the interaction between CaM and the Kv7. 4 A and B helices showed that calcium binding makes the C-lobe lose interaction with the A helix and bind to the B helix weakly through a much smaller interaction surface ( Xu et al., 2013 ; Chang et al., 2018 ). Xu et al. (2013) proposed that the interaction of apo/CaM with the Kv7. 4 A helix may affect the neighboring S6 in some way to favor channel opening and, in contrast, when CaM switches to Ca ²⁺/CaM and binds to the helix B, loss of the interaction with the A helix causes inhibition of the currents. Our data demonstrate that the third EF hand of CaM and the S2–S3 loop of Kv7. 4 channels play crucial roles in CaM-mediated calcium-dependent regulation of Kv7. 4 activation, and lead us to propose a new regulation model ( Figure 9 ). Apo/CaM binds only the C-terminal domain of Kv7. 4 channels through the interaction of the Ca ²⁺/C-lobe with the A helix and the interaction of the Ca ²⁺/N-lobe with the B helix, and has no association with the voltage sensor domain of Kv7. 4 channels. Thus, the voltage sensor could easily move in response to depolarized membrane potentials, which is coupled to pore opening. Upon calcium loading, the CaM C-lobe releases its interaction with the Kv7. 4 A helix and binds to the B helix weakly with a much smaller interaction surface ( Xu et al., 2013 ; Chang et al., 2018 ), making it possible for the Ca ²⁺/EF3 to contact the Kv7. 4 S2–S3 loop. This interaction between the Ca ²⁺/EF3 and the Kv7. 4 S2–S3 loop makes it possible for Ca ²⁺/CaM to “ pull” the voltage sensor of Kv7. 4 channels, and thus offers a force opposing the movement of the voltage sensor in response to depolarized potentials, affecting voltage-dependent pore opening. We could not deny the model proposed by Xu et al. (2013), although our model could better explain the shift in the voltage dependence of activation caused by CaM regulation. Probably, the inhibitory effect of Ca ²⁺/CaM on Kv7. 4 activation might be a combined effect of the “ pulled” voltage sensor and the released helix A. In our model, the interaction of the S2–S3 loop with the Ca ²⁺/EF3 rather than the Apo/EF3 is more likely to be driven by calcium-dependent changes of the constitutive interactions between CaM and the Kv7. 4 C-terminal domain. However, the residues in the Kv7. 4 A and B helices interacting with Apo/CaM and the residues in the Kv7. 4 B helix interacting with Ca ²⁺/CaM are different from the residues in Kv7. 1 to a certain degree, so there might be some differences between Kv7. 4 and Kv7. 1 in the calcium-dependent interaction with the CaM EF3, accounting for the different direction of CaM regulation.

FIGURE 9

A model for CaM-mediated calcium-dependent regulation of Kv7. 4 activation. The gray box represents the voltage sensor domain containing the Kv7. 4 transmembrane helixes S1–S4, with positively charged amino acids (shown as plus signs) in the S4 helix. The S2–S3 loop is shown as a red curve. The yellow box represents the pore domain. The brown rectangles represent the A and B helices. The CaM N-lobe is colored in green, and the C-lobe is colored in purple. The gray circles represent calcium bound to CaM. Apo/CaM binds to Kv7. 4 channels with the Apo/N-lobe interacting with the B helix and the Apo/C-lobe interacting with the A helix, and Apo/CaM has no association with the Kv7. 4 S2–S3 loop. Upon calcium loading, the Ca ²⁺/N-lobe remains anchored on the B helix, while the Ca ²⁺/C-lobe releases interaction with the A helix and binds to the B helix through a much smaller interaction surface, making it possible for the Ca ²⁺/EF3 to contact the Kv7. 4 S2–S3 loop and affect motivation of the voltage sensor and opening of the coupled pore.

Limited previous studies explored the role of the Kv7 S2–S3 loop in channel function. Gamper et al. (2006) found that increasing H ₂ O ₂ concentrations enhanced Kv7. 4 currents and facilitated Kv7. 4 activation and that currents of the Kv7. 4 triple mutant C156A/C157A/C158A remained unchanged after extracellular application of H ₂ O ₂ , so they concluded that the triple cysteines in the S2–S3 loop was critical for the oxidative modification. However, their data showed that extracellular application of H ₂ O ₂ still facilitated the activation rate of Kv7. 4 C156A/C157A/C158A, which seemed contradictory to the unchanged currents. Additionally, before extracellular application of H ₂ O ₂ and with only endogenous H ₂ O ₂ in cells, Kv7. 4 C156A/C157A/C158A showed a faster activation rate and a leftward shift in the half-activation voltage compared with WT Kv7. 4, which contrast the facilitating effects of H ₂ O ₂ and suggested that the S2–S3 loop might be involved in other regulations that affect channel function as well. Ooi et al. (2013) reported that extracellular application of SNAP that forms nitric oxide (NO) inhibited Kv7. 4 currents with an inhibition of about 30%, and that the NO inhibition of the Kv7. 4 triple mutant C156A/C157A/C158A was much smaller and almost abolished. They also tested S-nitrosylation of WT Kv7. 4 and Kv7. 4 C156A/C157A/C158A treated with SNAP and found that Kv7. 4 C156A/C157A/C158A showed decreased S-nitrosylation, so they concluded that NO possibly inhibited Kv7. 4 currents through S-nitrosylation of the triple cysteines in the S2–S3 loop. Both NO and Ca ²⁺/CaM inhibit Kv7. 4 channels, and the S2–S3 loop is required for both inhibition, implying a potential correlation between the two regulations. For example, possibly S-nitrosylation of the triple cysteines could favor the interaction between Ca ²⁺/CaM and the S2–S3 loop of Kv7. 4 channels. Further studies are required to explore the role of the S2-S3 loop in coordinating the process of H ₂ O ₂ regulation, NO regulation, and Ca ²⁺/CaM regulation.

Overall, our results highlight the importance of the Kv7. 4 S2–S3 linker in CaM regulation of channel activation. This is the first study to indicate another site in Kv7. 4 channels required for CaM regulation in addition to the C-terminal domain. Furthermore, given that the S2–S3 loop sequence is highly conserved among Kv7. 2–7. 5, this loop might play a similar role in CaM regulation of all these isoforms. More studies about the S2–S3 loop of the other Kv7 isoforms including Kv7. 1 may further our understanding of a unified mechanism for CaM regulation of Kv7 channels.

Data Availability Statement

The original contributions presented in the study are included in the article/ Supplementary Material , further inquiries can be directed to the corresponding author.

Author Contributions

WZ designed and implemented the experiments, collected and analyzed data, and wrote the paper. ZY supervised the project and provided guidance. Both the authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Key R&D Program of China Project (2017YFA0103900), the National Natural Science Foundation of China (31571083 and 31970931), the Program for Professor of Special Appointment (Eastern Scholar of Shanghai, TP2014008), the Shanghai Municipal Science and Technology Major Project (No. 2017SHZDZX01 and No. 2018SHZDZX01), and the Shanghai Rising-Star Program (14QA1400800).

Conflict of Interest

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

Acknowledgments

We are grateful to Jiahuai Han at Xiamen University and Lan Bao at the Chinese Academy of Sciences for providing the plasmids. We also appreciate Yanai Mei, Wenyong Fan, Honglan Zheng, and Zhaoyang Li for useful discussion.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www. frontiersin. org/articles/10. 3389/fphys. 2020. 604134/full#supplementary-material

Footnotes

1. https://somapp. ucdmc. ucdavis. edu/pharmacology/bers/maxchelator/webmaxc/webmaxcS. htm

References

Alaimo, A., Etxeberria, A., Gómez-Posada, J. C., Gomis-Perez, C., Fernández-Orth, J., Malo, C., et al. (2018). Lack of correlation between surface expression and currents in epileptogenic AB-calmodulin binding domain Kv72 potassium channel mutants. Channels 12, 299–310. doi: 10. 1080/19336950. 2018. 1511512

Alaimo, A., Gómez-Posada, J. C., Aivar, P., Etxeberría, A., Rodriguez-Alfaro, J. A., Areso, P., et al. (2009). Calmodulin activation limits the rate of KCNQ2 K+ channel exit from the endoplasmic reticulum. J. Biol. Chem. 284, 20668–20675. doi: 10. 1074/jbc. M109. 019539

Bal, M., Zaika, O., Martin, P., and Shapiro, M. S. (2008). Calmodulin binding to M-type K+ channels assayed by TIRF/FRET in living cells. J. Physiol. 586, 2307–2320. doi: 10. 1113/jphysiol. 2008. 152777

Barhanin, J., Lesage, F., Guillemare, E., Fink, M., Lazdunski, M., and Romey, G. (1996). K(V)LQT1 and lsK (minK) proteins associate to form the I(Ks) cardiac potassium current. Nature 384, 78–80. doi: 10. 1038/384078a0

Beisel, K. W., Nelson, N. C., Delimont, D. C., and Fritzsch, B. (2000). Longitudinal gradients of KCNQ4 expression in spiral ganglion and cochlear hair cells correlate with progressive hearing loss in DFNA2. Brain Res. Mol. Brain Res. 82, 137–149. doi: 10. 1016/S0169-328X(00)00204-7

Beisel, K. W., Rocha-Sanchez, S. M., Morris, K. A., Nie, L., Feng, F., Kachar, B., et al. (2005). Differential expression of KCNQ4 in inner hair cells and sensory neurons is the basis of progressive high-frequency hearing loss. J. Neurosci. 25, 9285–9293. doi: 10. 1523/JNEUROSCI. 2110-05. 2005

Bernardo-Seisdedos, G., Nunez, E., Gomis-Perez, C., Malo, C., Villarroel, A., and Millet, O. (2018). Structural basis and energy landscape for the Ca ²⁺ gating and calmodulation of the Kv7. 2 K(+) channel. Proc. Natl. Acad. Sci. U. S. A. 115, 2395–2400. doi: 10. 1073/pnas. 1800235115

Black, D. J., Halling, D. B., Mandich, D. V., Pedersen, S. E., Altschuld, R. A., and Hamilton, S. L. (2005). Calmodulin interactions with IQ peptides from voltage-dependent calcium channels. Am. J. Phys. Cell Physiol. 288, C669–C676. doi: 10. 1152/ajpcell. 00191. 2004

Brown, D. A., and Passmore, G. M. (2009). Neural KCNQ (Kv7) channels. Br. J. Pharmacol. 156, 1185–1195. doi: 10. 1111/j. 1476-5381. 2009. 00111. x

Chang, A., Abderemane-Ali, F., Hura, G. L., Rossen, N. D., Gate, R. E., and Minor, D. L. Jr. (2018). A Calmodulin C-lobe Ca ²⁺/-dependent switch governs Kv7 channel function. Neuron 97, 836. e836–852. e836. doi: 10. 1016/j. neuron. 2018. 01. 035

Delmas, P., and Brown, D. A. (2005). Pathways modulating neural KCNQ/M (Kv7) potassium channels. Nat. Rev. Neurosci. 6, 850–862. doi: 10. 1038/nrn1785

Etxeberria, A., Aivar, P., Rodriguez-Alfaro, J. A., Alaimo, A., Villacé, P., Gómez-Posada, J. C., et al. (2008). Calmodulin regulates the trafficking of KCNQ2 potassium channels. FASEB J. 22, 1135–1143. doi: 10. 1096/fj. 07-9712com

Evans, T. I., and Shea, M. A. (2009). Energetics of calmodulin domain interactions with the calmodulin binding domain of CaMKII. Proteins 76, 47–61. doi: 10. 1002/prot. 22317

Findeisen, F., Rumpf, C. H., and Minor, D. L. Jr. (2013). Apo states of calmodulin and CaBP1 control CaV1 voltage-gated calcium channel function through direct competition for the IQ domain. J. Mol. Biol. 425, 3217–3234. doi: 10. 1016/j. jmb. 2013. 06. 024

Gamper, N., Li, Y., and Shapiro, M. S. (2005). Structural requirements for differential sensitivity of KCNQ K+ channels to modulation by Ca2+/calmodulin. Mol. Biol. Cell 16, 3538–3551. doi: 10. 1091/mbc. e04-09-0849

Gamper, N., and Shapiro, M. S. (2003). Calmodulin mediates Ca2+−dependent modulation of M-type K+ channels. J. Gen. Physiol. 122, 17–31. doi: 10. 1085/jgp. 200208783

Gamper, N., Zaika, O., Li, Y., Martin, P., Hernandez, C. C., Perez, M. R., et al. (2006). Oxidative modification of M-type K(+) channels as a mechanism of cytoprotective neuronal silencing. EMBO J. 25, 4996–5004. doi: 10. 1038/sj. emboj. 7601374

Geiser, J. R., van Tuinen, D., Brockerhoff, S. E., Neff, M. M., and Davis, T. N. (1991). Can calmodulin function without binding calcium? Cell 65, 949–959. doi: 10. 1016/0092-8674(91)90547-C

Ghosh, S., Nunziato, D. A., and Pitt, G. S. (2006). KCNQ1 assembly and function is blocked by long-QT syndrome mutations that disrupt interaction with calmodulin. Circ. Res. 98, 1048–1054. doi: 10. 1161/01. RES. 0000218863. 44140. f2

Gomis-Perez, C., Soldovieri, M. V., Malo, C., Ambrosino, P., Taglialatela, M., Areso, P., et al. (2017). Differential regulation of PI(4, 5)P ₂ sensitivity of Kv7. 2 and Kv7. 3 channels by calmodulin. Front. Mol. Neurosci. 10: 117. doi: 10. 3389/fnmol. 2017. 00117

Gutman, G. A., Chandy, K. G., Adelman, J. P., Aiyar, J., Bayliss, D. A., Clapham, D. E., et al. (2003). International union of pharmacology. XLI. Compendium of voltage-gated ion channels: potassium channels. Pharmacol. Rev. 55, 583–586. doi: 10. 1124/pr. 55. 4. 9

Hernandez, C. C., Zaika, O., Tolstykh, G. P., and Shapiro, M. S. (2008). Regulation of neural KCNQ channels: signalling pathways, structural motifs and functional implications. J. Physiol. 586, 1811–1821. doi: 10. 1113/jphysiol. 2007. 148304

Jentsch, T. J. (2000). Neuronal KCNQ potassium channels: physiology and role in disease. Nat. Rev. Neurosci. 1, 21–30. doi: 10. 1038/35036198

Kharkovets, T., Hardelin, J. P., Safieddine, S., Schweizer, M., El-Amraoui, A., Petit, C., et al. (2000). KCNQ4, a K+ channel mutated in a form of dominant deafness, is expressed in the inner ear and the central auditory pathway. Proc. Natl. Acad. Sci. U. S. A. 97, 4333–4338. doi: 10. 1073/pnas. 97. 8. 4333

Kim, H. J., Lv, P., Sihn, C. R., and Yamoah, E. N. (2011). Cellular and molecular mechanisms of autosomal dominant form of progressive hearing loss, DFNA2. J. Biol. Chem. 286, 1517–1527. doi: 10. 1074/jbc. M110. 179010

Kosenko, A., and Hoshi, N. (2013). A change in configuration of the calmodulin-KCNQ channel complex underlies Ca2+−dependent modulation of KCNQ channel activity. PLoS One 8: e82290. doi: 10. 1371/journal. pone. 0082290

Kubisch, C., Schroeder, B. C., Friedrich, T., Lütjohann, B., El-Amraoui, A., Marlin, S., et al. (1999). KCNQ4, a novel potassium channel expressed in sensory outer hair cells, is mutated in dominant deafness. Cell 96, 437–446. doi: 10. 1016/s0092-8674(00)80556-5

Lee, C. H., and MacKinnon, R. (2018). Activation mechanism of a human SK-calmodulin channel complex elucidated by cryo-EM structures. Science 360, 508–513. doi: 10. 1126/science. aas9466

Linse, S., Helmersson, A., and Forsén, S. (1991). Calcium binding to calmodulin and its globular domains. J. Biol. Chem. 266, 8050–8054.

Liu, W., and Devaux, J. J. (2014). Calmodulin orchestrates the heteromeric assembly and the trafficking of KCNQ2/3 (Kv7. 2/3) channels in neurons. Mol. Cell. Neurosci. 58, 40–52. doi: 10. 1016/j. mcn. 2013. 12. 005

Olwin, B. B., Edelman, A. M., Krebs, E. G., and Storm, D. R. (1984). Quantitation of energy coupling between Ca2+, calmodulin, skeletal muscle myosin light chain kinase, and kinase substrates. J. Biol. Chem. 259, 10949–10955.

Olwin, B. B., and Storm, D. R. (1985). Calcium binding to complexes of calmodulin and calmodulin binding proteins. Biochemistry 24, 8081–8086. doi: 10. 1021/bi00348a037

Ooi, L., Gigout, S., Pettinger, L., and Gamper, N. (2013). Triple cysteine module within M-type K+ channels mediates reciprocal channel modulation by nitric oxide and reactive oxygen species. J. Neurosci. 33, 6041–6046. doi: 10. 1523/JNEUROSCI. 4275-12. 2013

Peersen, O. B., Madsen, T. S., and Falke, J. J. (1997). Intermolecular tuning of calmodulin by target peptides and proteins: differential effects on Ca2+ binding and implications for kinase activation. Protein Sci. 6, 794–807. doi: 10. 1002/pro. 5560060406

Sachyani, D., Dvir, M., Strulovich, R., Tria, G., Tobelaim, W., Peretz, A., et al. (2014). Structural basis of a Kv7. 1 potassium channel gating module: studies of the intracellular c-terminal domain in complex with calmodulin. Structure 22, 1582–1594. doi: 10. 1016/j. str. 2014. 07. 016

Sanguinetti, M. C., Curran, M. E., Zou, A., Shen, J., Spector, P. S., Atkinson, D. L., et al. (1996). Coassembly of K(V)LQT1 and minK (IsK) proteins to form cardiac I(Ks) potassium channel. Nature 384, 80–83. doi: 10. 1038/384080a0

Shamgar, L., Ma, L., Schmitt, N., Haitin, Y., Peretz, A., Wiener, R., et al. (2006). Calmodulin is essential for cardiac IKS channel gating and assembly: impaired function in long-QT mutations. Circ. Res. 98, 1055–1063. doi: 10. 1161/01. res. 0000218979. 40770. 69

Sihn, C. R., Kim, H. J., Woltz, R. L., Yarov-Yarovoy, V., Yang, P. C., Xu, J., et al. (2016). Mechanisms of calmodulin regulation of different isoforms of Kv7. 4 K+ channels. J. Biol. Chem. 291, 2499–2509. doi: 10. 1074/jbc. M115. 668236

Strulovich, R., Tobelaim, W. S., Attali, B., and Hirsch, J. A. (2016). Structural insights into the M-channel proximal C-terminus/calmodulin complex. Biochemistry 55, 5353–5365. doi: 10. 1021/acs. biochem. 6b00477

Sun, J., and MacKinnon, R. (2017). Cryo-EM structure of a KCNQ1/CaM complex reveals insights into congenital long QT syndrome. Cell 169, 1042. e1049–1050. e1049. doi: 10. 1016/j. cell. 2017. 05. 019

Tobelaim, W. S., Dvir, M., Lebel, G., Cui, M., Buki, T., Peretz, A., et al. (2017). Competition of calcified calmodulin N lobe and PIP2 to an LQT mutation site in Kv71 channel. Proc. Natl. Acad. Sci. U. S. A. 114, E869–e878. doi: 10. 1073/pnas. 1612622114

Wang, Q., Curran, M. E., Splawski, I., Burn, T. C., Millholland, J. M., VanRaay, T. J., et al. (1996). Positional cloning of a novel potassium channel gene: KVLQT1 mutations cause cardiac arrhythmias. Nat. Genet. 12, 17–23. doi: 10. 1038/ng0196-17

Wang, H. S., Pan, Z., Shi, W., Brown, B. S., Wymore, R. S., Cohen, I. S., et al. (1998). KCNQ2 and KCNQ3 potassium channel subunits: molecular correlates of the M-channel. Science 282, 1890–1893. doi: 10. 1126/science. 282. 5395. 1890

Wen, H., and Levitan, I. B. (2002). Calmodulin is an auxiliary subunit of KCNQ2/3 potassium channels. J. Neurosci. 22, 7991–8001. doi: 10. 1523/JNEUROSCI. 22-18-07991. 2002

Xia, X. M., Fakler, B., Rivard, A., Wayman, G., Johnson-Pais, T., Keen, J. E., et al. (1998). Mechanism of calcium gating in small-conductance calcium-activated potassium channels. Nature 395, 503–507. doi: 10. 1038/26758

Xu, Q., Chang, A., Tolia, A., and Minor, D. L. Jr. (2013). Structure of a Ca ²⁺/CaM: Kv7. 4 (KCNQ4) B-helix complex provides insight into M current modulation. J. Mol. Biol. 425, 378–394. doi: 10. 1016/j. jmb. 2012. 11. 023

Xu, T., Nie, L., Zhang, Y., Mo, J., Feng, W., Wei, D., et al. (2007). Roles of alternative splicing in the functional properties of inner ear-specific KCNQ4 channels. J. Biol. Chem. 282, 23899–23909. doi: 10. 1074/jbc. M702108200

Yus-Najera, E., Santana-Castro, I., and Villarroel, A. (2002). The identification and characterization of a noncontinuous calmodulin-binding site in noninactivating voltage-dependent KCNQ potassium channels. J. Biol. Chem. 277, 28545–28553. doi: 10. 1074/jbc. M204130200

Zhang, Z. N., Li, Q., Liu, C., Wang, H. B., Wang, Q., and Bao, L. (2008). The voltage-gated Na+ channel Nav1. 8 contains an ER-retention/retrieval signal antagonized by the beta3 subunit. J. Cell Sci. 121, 3243–3252. doi: 10. 1242/jcs. 026856