Article navigation

Gain-of-function (GoF) missense variants in the two-pore domain (K2P) K+ channel TASK-1 (KCNK3) result in developmental delay with sleep apnea (DDSA), a neurodevelopmental channelopathy, while loss-of-function (LoF) variants cause pulmonary arterial hypertension. However, for the related TASK-3 channel (KCNK9), both LoF and GoF variants underlie a distinct neurodevelopmental disorder, KCNK9 imprinting syndrome (KIS). The relationship between genotype and phenotype in these disorders is further complicated because TASK-1 and TASK-3 can co-assemble into heteromeric channels with distinct functional properties. Here, we report additional patients with missense variants in KCNK3 and KCNK9 and investigate the effect of four novel genetic variants on the functional properties of homomeric and heteromeric TASK channels. Interestingly, two of these new pathogenic GoF variants (R131H and L122V) are found in both TASK-1 and TASK-3 and have equivalent functional effects on heteromeric TASK-1/TASK-3 channel activity, yet result in different clinical phenotypes. We have also determined a cryo-EM structure for the pathogenic L122V mutant TASK-3 channel, which suggests that subtle changes in gating and permeation within the inner cavity are responsible for its activatory effect. Overall, these results highlight the dominant role that homomeric TASK channels likely play in defining their associated channelopathies as well as the complexity of interpreting K+ channel dysfunction in pathophysiology.

Genetic variants that affect the function of K2P K+ channels are now increasingly associated with a range of diseases, including diabetes, nephropathy, migraine, sleep apnea, hypertension, and other cardiac disorders, as well as a range of neurodevelopmental disorders (Lee et al., 2021). However, despite the pressing need, there are few effective treatments for these complex conditions, and their underlying mechanisms are often not well understood.

K2P channels represent a structurally distinct subfamily of K+ channels. In addition to their distinct extracellular Cap domain, each subunit contains four transmembrane helices (M1–M4) and two pore domains (P1 and P2). Two subunits therefore assemble as a “dimer of dimers” to create a pseudotetrameric channel that is similar to other classical K+ channel pores (Enyedi and Czirjak, 2010; Niemeyer et al., 2016). There are 15 different human K2P (KCNK) channels, and members of this subfamily of K+ channels are often referred to as simple “leak” conductances involved in establishing the resting membrane potential, but it is now clear they exhibit complex patterns of dynamic regulation by diverse stimuli, including many G-protein–coupled receptor (GPCR)-associated regulatory pathways (Enyedi and Czirjak, 2010; Niemeyer et al., 2016). Understanding K2P channel dysfunction in the disease state is therefore critical for a clearer insight into these rare channelopathies and to exploit the therapeutic potential of K2P channels in other more common disease states (Bagal et al., 2013; Mathie et al., 2021).

KCNK3 and KCNK9 encode the TASK-1 and TASK-3 K2P channels, respectively. Although separate genes with distinct tissue-specific patterns of expression, they share almost 80% amino acid identity within the main pore-forming transmembrane regions (Fig. S1). Consequently, they exhibit significant structural and functional similarity (Rodstrom et al., 2020; Lin et al., 2024; Hall et al., 2025), and in those cell types where they are coexpressed, they can co-assemble into heteromeric TASK-1/TASK-3 channels, which possess distinct functional properties (Berg et al., 2004; Kim et al., 2009; Bohnen et al., 2017; Rinné et al., 2024).

+ Expand view − Collapse view
Figure S1
Figure S1. Refer to the image caption for details. The sequences are labeled as TASK-1 and TASK-3, with corresponding amino acid positions indicated. Various regions such as transmembrane domains (M1, M2, M3, M4), extracellular cap domains (EC1, EC2), pore domains (P1, P2), selectivity filters (SF1, SF2), and the C-terminal domain (CTD) are highlighted. Asterisks mark the positions of variants, with red asterisks indicating existing variants and blue asterisks indicating new variants described in the study. The alignment shows a high degree of homology between the two sequences, particularly within the pore-forming transmembrane regions.

Homology between TASK-1 and TASK-3. An amino acid alignment between the sequences of TASK-1 (KCNK3) and TASK-3 (KCNK9) reveals the high degree of homology between the two channel subunits. Overall, the sequence identity is 54% with almost 80% identity within the pore-forming transmembrane regions. The extent of the homology means that within the TM/Pore region the equivalent amino acid numbering is identical and does not deviate except within the CTD. The various transmembrane regions (M1–M4) are indicated by blue bars above the sequence along with the Extracellular Cap domain helices and the selectivity filters. The position of the variants is indicated by asterisks (red = existing variant, blue = new variant described in this study). This clearly demonstrates that M2 is a hotspot for pathogenic GoF variants in TASK-1.

Figure S1.
Amino acid sequence alignment between TASK-1 and TASK-3. The sequences are labeled as TASK-1 and TASK-3, with corresponding amino acid positions indicated. Various regions such as transmembrane domains (M1, M2, M3, M4), extracellular cap domains (EC1, EC2), pore domains (P1, P2), selectivity filters (SF1, SF2), and the C-terminal domain (CTD) are highlighted. Asterisks mark the positions of variants, with red asterisks indicating existing variants and blue asterisks indicating new variants described in the study. The alignment shows a high degree of homology between the two sequences, particularly within the pore-forming transmembrane regions.

Homology between TASK-1 and TASK-3. An amino acid alignment between the sequences of TASK-1 (KCNK3) and TASK-3 (KCNK9) reveals the high degree of homology between the two channel subunits. Overall, the sequence identity is 54% with almost 80% identity within the pore-forming transmembrane regions. The extent of the homology means that within the TM/Pore region the equivalent amino acid numbering is identical and does not deviate except within the CTD. The various transmembrane regions (M1–M4) are indicated by blue bars above the sequence along with the Extracellular Cap domain helices and the selectivity filters. The position of the variants is indicated by asterisks (red = existing variant, blue = new variant described in this study). This clearly demonstrates that M2 is a hotspot for pathogenic GoF variants in TASK-1.

Close modal

TASK-1 (KCNK3) is broadly expressed in many cells and tissues involved in the control of breathing and ventilation, including the carotid body, lung, heart, and pulmonary artery smooth muscle, as well as many ventilatory-associated chemosensitive regions of the brain, hypoglossal motor neurons, and spinal cord motor neurons (Trapp et al., 2008; Bayliss et al., 2015; Lin et al., 2025). By contrast, TASK-3 (KCNK9) expression shows a higher degree of specialization within the central nervous system (CNS). However, their patterns of expression do overlap in some cell types; for example, both subunits are coexpressed in carotid body glomus cells (Buckler et al., 2000; Buckler, 2015), and heteromeric TASK-1/TASK-3 channels in these cells are proposed to form the major oxygen-sensitive background K+ channel, which plays a critical role in the regulation of respiration. Both TASK-1 and TASK-3 are also a target for volatile anesthetics (Conway and Cotten, 2012; Luethy et al., 2017), and given their key role in ventilation and respiratory control, inhibitors that target TASK-1 have been developed for the treatment of sleep apnea with improvements observed in a subset of patients, where topical application of the inhibitor in the upper airway was achieved via nasal spray (Osman et al., 2024).

Given their physiological role and expression profile, it was perhaps not surprising that inherited monoallelic loss-of-function (LoF) variants in KCNK3 have been shown to be associated with a form of pulmonary arterial hypertension (Ma et al., 2013; Olschewski et al., 2017; Jouen-Tachoire et al., 2021). Furthermore, we have also recently shown that de novo heterozygous missense gain-of-function (GoF) variants underlie a novel neurodevelopmental disorder (developmental delay with sleep apnea [DDSA]), where all affected individuals also have sleep apnea (Kaplanis et al., 2020; Sörmann et al., 2022). Structural analysis reveals the majority of these mutations cluster in the M2 and M4 transmembrane helices close to an intersection between which forms the “X-gate,” a structural motif analogous to the classical bundle-crossing gate in other K+ channels (Rodstrom et al., 2020; Sörmann et al., 2022). All of these variants affect channel gating and produce TASK-1 channels with a higher intrinsic single-channel open probability (Po) combined with a markedly reduced sensitivity to inhibition by Gq-protein–coupled receptor (GqPCR) pathways (Sörmann et al., 2022). We have recently shown that this defect in GPCR inhibition is directly coupled to their higher intrinsic Po and the state-dependent inhibition of TASK channels by PLC-generated diacylglycerol (Jouen-Tachoire et al., 2026).

Missense variants in TASK-3 (KCNK9) are associated with a different neurodevelopmental disorder, “KCNK9 imprinting syndrome” (KIS), previously known as Birk-Barel syndrome (Barel et al., 2008; Graham et al., 2016). Typically, KIS-associated mutations produce a LoF in TASK-3, but more recently, GoF mutations have also been reported in individuals that present similar phenotypes (Cousin et al., 2022). Sleep apnea has also been reported in some KIS patients (Cousin et al., 2022) but is not a universal feature, unlike its defining phenotype in DDSA. It remains unclear why both LoF and GoF variants in KCNK9 result in the same phenotype, i.e., KIS. This complex etiology may be due to genomic imprinting of KCNK9, which has not been reported for KCNK3 or any other K2P channels, as well as the mosaicism associated with the de novo origin of these variants (Luedi et al., 2007; Cooper et al., 2020; Cousin et al., 2022). Another complicating factor may be their ability to form heteromeric channels as the relative contribution of defects in homomeric vs. heteromeric channels to these different disorders also remains unclear.

In this study, we identify six new affected individuals with GoF mutations in TASK-1, which display a DDSA phenotype, and these include three novel missense variants in KCNK3 (F125S, Q126E, and R131H). We also identify a novel missense variant in KCNK9 (L122V) associated with KIS and solve a cryo-EM structure of this mutant channel. Overall, the structural and functional characterization of these different mutations provide important insights into the relative contributions of both homomeric and heteromeric TASK channels in the development of TASK-related disorders.

Patient identification, genetic and clinical investigation, and ethical approval

Using the GeneMatcher platform (Kaplanis et al., 2020), Genomics England data derived from the 100,000 genomes project, and extensive international data sharing, we identified seven families. A uniform clinical pro forma was distributed to collect clinical details. Parents or legal guardians of all affected individuals provided informed consent for the publication of clinical information and genetic data, in accordance with the Declaration of Helsinki. The study was approved by the relevant local Ethics Committees, including University College London (for Integrated Research Application System Project ID: 310045), while the original projects (Kaplanis et al., 2020; Sörmann et al., 2022) have UK Research Ethics Committee approval (10/H0305/83, granted by the Cambridge South REC, and GEN/284/12, granted by the Republic of Ireland REC). Trio exome or genome sequencing, as well as patient-only sequencing where applicable, was performed on DNA extracted from blood-derived leukocytes. Data analysis and variant filtration were conducted using standard pipelines, and Sanger sequencing was used for segregation analysis when indicated.

Electrophysiology

Whole-cell currents were measured using two-electrode voltage clamp (TEVC) recordings of Xenopus laevis oocytes microinjected with the relevant mRNA. This was done as previously described in accordance with local ethical and regulatory requirements (Sörmann et al., 2022). Briefly, the WT human TASK-1 gene (KCNK3) and human TASK-3 gene (KCNK9) were subcloned into pFAW. Mutations were introduced by site-directed mutagenesis and confirmed by sequencing. mRNA was transcribed using the T7 mScript Standard mRNA Production System. Each oocyte was then injected with 2 ng of RNA and incubated for 20–24 h at 17.5°C. Recordings were performed in ND96 buffer at pH 7.4 (96 mM NaCl, 2 mM KCl, 2 mM MgCl2, 1.8 mM CaCl2, and 5 mM HEPES). Currents were recorded using a 400 ms voltage step protocol from a holding potential of −80 mV delivered in 10 mV increments between −120 mV and +50 mV and 800 ms ramp protocols from −120 to +50 mV. All recorded traces were analyzed using Clampfit, and graphs were plotted using Origin 2021. Single-channel currents in cell-attached patches were recorded under quasi-symmetrical potassium concentrations with 140 mM KCl solution containing 10 mM HEPES (pH 7.4 with KOH) in the pipette. Data were filtered at 2 or 5 kHz and recorded at a 200 kHz sampling rate with program Clampex on an Axopatch 200B amplifier. Data analysis was performed using Clampfit. Channel open probability (Po) was determined from single-channel recordings with 1–4-min duration, which exhibited only one main open level. Numbers of individual recordings (n) are quoted for each value along with the test used for any statistical comparison (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; * = P < 0.05; ** = P < 0.01; *** = P < 0.001.

Cryo-EM sample preparation and collection

TASK-3 L122V protein was expressed and purified using the method previously described for WT TASK-3 (Hall et al., 2025). Samples of TASK-3 L122V (4.0 mg/ml) were then applied onto glow-discharged AuFlat 300 1.2/1.3 gold grids (Molecular Dimensions). The grids were blotted for 4–5 s at 100% humidity and 4°C and then vitrified in liquid ethane using a Vitrobot Mark IV (Thermo Fisher Scientific). Data were collected on a Titan Krios G3 (FEI) operating at 300 kV, equipped with a BioQuantum imaging filter (Gatan) and a K3 direct detection camera (Gatan). Images were acquired at a magnification of 105,000×, corresponding to a physical pixel size of 0.832 Å with a total dose rate of 43.9 e-2. A total of 13,035 initial movies were collected.

Image processing

All processing was performed in cryoSPARC (Punjani et al., 2017) unless stated otherwise. Movies were imported and underwent motion correction and patch contrast transfer function estimation. A subset of 378 micrographs was selected for picking with a particle diameter of 125 Å. The initially extracted particles underwent two rounds of 2D classification and ab initio reconstruction (n = 2), which yielded a volume used to generate templates for particle picking.

The TASK-3 L122V templates were used for particle repicking from the full dataset with a diameter of 125 Å, resulting in the extraction of 8,381,007 particles with box size of 256 px, Fourier cropped to 128 px. These particles underwent four rounds of 2D classification during which junk particles were removed, leaving 1,718,239 particles. The particles were re-extracted with a full 256 px box size, prior to two more rounds of 2D classification. Two rounds of five ab initio classes and heterogenous refinement in C1 and C2 symmetry were completed, leading to a particle stack of 332,112 particles. Homogenous refinement, nonuniform (NU) refinement, and reference-based motion correction in C2 produced a 2.90 Å map. Polished particles were processed with cryosieve (Zhu et al., 2023), and the particle stacks were reimported into cryoSPARC. To identify the highest resolution stack, each stack was subjected to one round of 2D classification, ab initio reconstruction in C1, followed by homogeneous and NU refinements in C2. The final particle set (211,335 particles) produced a map with a resolution of 2.83 Å, based on the FSC = 0.143 criterion.

Model building and refinement

WT TASK-3 (Protein Data Bank [PDB] ID: 9G9V) was used as the initial model to build TASK-3 L122V. The structure was manually mutated and adjusted within the map using Coot (Emsley et al., 2010). The model was further refined iteratively, first through manual adjustments in Coot and then using phenix.real_space_refine (Adams et al., 2010). The refined model was subsequently processed with ISOLDE (Croll, 2018) through ChimeraX (Meng et al., 2023) to enhance their accuracy and underwent a final round of refinement in Phenix, with the ISOLDE-optimized models serving as references. Structural figures were generated using PyMOL (Schrödinger, LLC).

Online supplemental material

Fig. S1 shows an annotated alignment of the sequences of TASK-1 and TASK-3. Fig. S2 shows a summary of the functional properties of the TASK-1 T121M variant. Fig. S3 illustrates the image processing workflow and associated data for the structure of TASK-3 L122V. Table S1 summarizes the phenotype and anonymized clinical details of the new DDSA and KIS patients reported in this study. Table S2 summarizes the cryo-EM data collection, refinement, and validation statistics for the TASK-2 L122V structure reported in this study.

Three novel DDSA variants in KCNK3

Our previous study identified nine individuals with DDSA (Sörmann et al., 2022), each heterozygous for one of six de novo missense variants in KCNK3. Since that original study in 2022, we have now identified another three DDSA patients with similar de novo GoF mutations in KCNK3, including two with the recurrent N133S variant and one with the L241F variant. In addition, we have identified three further affected individuals with novel heterozygous missense variants in KCNK3 who all match the reported DDSA neurodevelopmental phenotype (Table S1). Interestingly, the three new variants (F125S, Q126E, and R131H) all cluster in the M2 transmembrane helix adjacent to several previously characterized DDSA variants (Fig. 1, a and b). The location of these novel variants alongside many others within M2 indicates that this particular region of the transmembrane helix clearly represents a major hotspot for pathogenic mutations and is consistent with the proposed role for the TM helices and the X-gate in channel gating (Rodstrom et al., 2020; Sörmann et al., 2022; Hall et al., 2025). It is also consistent with the reported effects of other mutations within this region on K2P channel activity (Ben Soussia et al., 2019).

Figure 1.
A multi-panel image depicts the functional effects of new DDSA-associated variants on TASK-1 channels. Panel a shows a schematic diagram of the TASK-1 channel highlighting transmembrane domains and mutation sites. Panel b shows a three-dimensional structural model of TASK-1 with mutation positions and channel gates. Panel c shows representative current traces comparing wild-type TASK-1 and mutant channel activities. Panel d shows current-voltage plots comparing wild-type and mutant TASK-1 channel currents. Panel e shows a bar plot comparing normalized currents of homomeric and heterozygous TASK-1 channels. Panel f shows representative single-channel recordings comparing wild-type and mutant TASK-1 channel activities. Panel g shows a bar plot comparing the open probability of wild-type and mutant TASK-1 channels.

Functional effects of new DDSA-associated variants. (a) Transmembrane topology of TASK-1 shown with the position of the new and previously characterized DDSA variants labelled in blue and red, respectively. The position of the X-gate on M4 is also shown. (b) TASK-1 crystal structure (PDB: 6RV2) with the two subunits shown in green and grey, respectively. The new and previously characterized DDSA mutations are shown on one subunit as blue and red spheres, respectively. (c) Representative TEVC recordings of WT TASK-1 and DDSA mutant currents in response to voltage steps from −120 to +50 mV in 20 mV steps from a holding potential of −80 mV. (d) Current–voltage plot of WT TASK-1 (n = 8), F125S (n = 16), Q126E (n = 19), and R131H (n = 23); data are presented as mean ± SD. (e) Currents for homomeric DDSA mutants, and heterozygous channels formed from 1:1 coexpression of WT TASK-1 and DDSA mutants normalized to WT current at +50 mV: WT (n = 43), F125S (n = 26, P = 1.4 × 10−22), F125S-WT (n = 16, P = 4.4 × 10−6), Q126E (n = 22, P = 1.1 × 10−14), Q126E-WT (n = 12, P = 4.7 × 10−6), R131H (n = 17, P = 1.1 × 10−55), and R131H-WT (n = 23, P = 5.8 × 10−8); data are presented as mean ± SD All mutant currents differ from WT (one-way ANOVA with Dunnett’s multiple comparisons test). (f) Representative single-channel recordings of new DDSA mutants at a holding potential of −160 mV. (g) Single-channel Po for each mutant. WT (n = 8), F125S (n = 8, P = 5.4 × 10−7), Q126E (n = 6, P = 0.0037), and R131H (n = 6, P = 3.2 × 10−7). All mutant Po values differ from that for WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; ** = P < 0.01; *** = P < 0.001.

Figure 1.
A multi-panel image depicts the functional effects of new DDSA-associated variants on TASK-1 channels. Panel a shows a schematic diagram of the TASK-1 channel highlighting transmembrane domains and mutation sites. Panel b shows a three-dimensional structural model of TASK-1 with mutation positions and channel gates. Panel c shows representative current traces comparing wild-type TASK-1 and mutant channel activities. Panel d shows current-voltage plots comparing wild-type and mutant TASK-1 channel currents. Panel e shows a bar plot comparing normalized currents of homomeric and heterozygous TASK-1 channels. Panel f shows representative single-channel recordings comparing wild-type and mutant TASK-1 channel activities. Panel g shows a bar plot comparing the open probability of wild-type and mutant TASK-1 channels.

Functional effects of new DDSA-associated variants. (a) Transmembrane topology of TASK-1 shown with the position of the new and previously characterized DDSA variants labelled in blue and red, respectively. The position of the X-gate on M4 is also shown. (b) TASK-1 crystal structure (PDB: 6RV2) with the two subunits shown in green and grey, respectively. The new and previously characterized DDSA mutations are shown on one subunit as blue and red spheres, respectively. (c) Representative TEVC recordings of WT TASK-1 and DDSA mutant currents in response to voltage steps from −120 to +50 mV in 20 mV steps from a holding potential of −80 mV. (d) Current–voltage plot of WT TASK-1 (n = 8), F125S (n = 16), Q126E (n = 19), and R131H (n = 23); data are presented as mean ± SD. (e) Currents for homomeric DDSA mutants, and heterozygous channels formed from 1:1 coexpression of WT TASK-1 and DDSA mutants normalized to WT current at +50 mV: WT (n = 43), F125S (n = 26, P = 1.4 × 10−22), F125S-WT (n = 16, P = 4.4 × 10−6), Q126E (n = 22, P = 1.1 × 10−14), Q126E-WT (n = 12, P = 4.7 × 10−6), R131H (n = 17, P = 1.1 × 10−55), and R131H-WT (n = 23, P = 5.8 × 10−8); data are presented as mean ± SD All mutant currents differ from WT (one-way ANOVA with Dunnett’s multiple comparisons test). (f) Representative single-channel recordings of new DDSA mutants at a holding potential of −160 mV. (g) Single-channel Po for each mutant. WT (n = 8), F125S (n = 8, P = 5.4 × 10−7), Q126E (n = 6, P = 0.0037), and R131H (n = 6, P = 3.2 × 10−7). All mutant Po values differ from that for WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; ** = P < 0.01; *** = P < 0.001.

Close modal

Novel DDSA variants also produce GoF channels with impaired GPCR sensitivity

All existing DDSA variants are associated with a GqPCR-insensitive GoF effect on TASK-1 channel activity. Therefore, to confirm these new DDSA variants also fit this profile, we next examined their functional activity. Whole-cell K+ currents of WT and mutant TASK-1 channels were measured by TEVC after expression in Xenopus oocytes. These results show that all three variants produced markedly larger currents than WT TASK-1 (Fig. 1, c and d).

K2P channels assemble as dimers, and each of the new DDSA patients are also heterozygous for these KCNK3 variants. To better replicate this, we therefore also measured average whole-cell currents from oocytes co-injected with equal amounts of WT and mutant TASK-1 RNA. Similar to the results obtained with the previous DDSA variants, we found that all three novel variants also increased these “heterozygous” currents by ∼50% (Fig. 1 e).

To confirm that these increases are due to an increase in the intrinsic activity of the channel rather than just an increase in trafficking to the cell surface, we also measured single-channel currents of homomeric WT and mutant TASK-1 channels and found that all three variants resulted in a marked (∼10-fold) increase in channel open probability (Fig. 1, f and g). A modest (∼20%) increase in single-channel conductance was also observed for the F125S and Q126E variants, but this was not large enough to account for their increased whole-cell currents (>250%).

WT TASK-1 channels are sensitive to inhibition by GqPCR mediated pathways, and the known DDSA variants are characterized by a dramatic reduction in this sensitivity that exacerbates their intrinsic GoF effect by uncoupling them from their regulatory pathways (Sörmann et al., 2022). We therefore also examined the GqPCR sensitivity of these new variants and found that whereas ATP-induced stimulation of a coexpressed P2Y2 receptor results in ∼50% inhibition of WT TASK-1 currents, all three variants exhibited a markedly reduced response with only 15–20% inhibition observed in identical conditions. This effect was seen in both homomeric and heterozygous TASK-1 mutant channels (Fig. 2). Thus, although the inheritance of the Q126E variant is unknown, its clear GoF effect as a heterozygous variant suggests it may also be of de novo origin, but this would need to be confirmed.

Figure 2.
Graphs depict current measurements and relative current comparisons. Panel a shows four line graphs depicting current measurements over time with 15 micromolar ATP applied extracellularly. The x-axis represents time in seconds, and the y-axis represents current in microamperes. The graphs compare wild-type (WT/WT) TASK-1 channels and heterozygous channels from coexpressed WT and mutant subunits coexpressed with the P2Y2 receptor. Panel b shows a bar graph comparing residual current for homozygous mutants and mutants coexpressed 1:1 with WT TASK-1 after the addition of 15 micromolar ATP. The x-axis lists different conditions (WT, F125S, F125S/WT, Q126E, Q126E/WT, R131H, R131H/WT), and the y-axis represents relative current. The bar graph includes error bars representing standard deviation and significance markers indicating statistical differences.

GqPCR-mediated inhibition is reduced by DDSA-associated mutations. (a) Representative currents at +50 mV of homomeric WT TASK-1 channels (WT/WT) and heterozygous channels from coexpressed WT and mutant subunits coexpressed with the P2Y2 receptor. Currents were recorded over time with 15 µM ATP applied extracellularly where indicated. This concentration produces ∼50% inhibition of WT TASK-1. (b) Residual current for homozygous mutants and mutants coexpressed 1:1 with WT TASK-1 after addition of 15 μM ATP. WT TASK-1 (n = 17), F125S (n = 28, P = 2.7 × 10−15), F125S/WT (n = 48, P = 4.4 × 10−16), Q126E (n = 15, P = 8.6 × 10−13), Q126E/WT (n = 36, P = 1.4 × 10−18), R131H (n = 22, P = 3.1 × 10−15), and R131H/WT (n = 45, P = 2.1 × 10−13); data are presented as mean ± SD. All mutant currents differ from WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; *** = P < 0.001.

Figure 2.
Graphs depict current measurements and relative current comparisons. Panel a shows four line graphs depicting current measurements over time with 15 micromolar ATP applied extracellularly. The x-axis represents time in seconds, and the y-axis represents current in microamperes. The graphs compare wild-type (WT/WT) TASK-1 channels and heterozygous channels from coexpressed WT and mutant subunits coexpressed with the P2Y2 receptor. Panel b shows a bar graph comparing residual current for homozygous mutants and mutants coexpressed 1:1 with WT TASK-1 after the addition of 15 micromolar ATP. The x-axis lists different conditions (WT, F125S, F125S/WT, Q126E, Q126E/WT, R131H, R131H/WT), and the y-axis represents relative current. The bar graph includes error bars representing standard deviation and significance markers indicating statistical differences.

GqPCR-mediated inhibition is reduced by DDSA-associated mutations. (a) Representative currents at +50 mV of homomeric WT TASK-1 channels (WT/WT) and heterozygous channels from coexpressed WT and mutant subunits coexpressed with the P2Y2 receptor. Currents were recorded over time with 15 µM ATP applied extracellularly where indicated. This concentration produces ∼50% inhibition of WT TASK-1. (b) Residual current for homozygous mutants and mutants coexpressed 1:1 with WT TASK-1 after addition of 15 μM ATP. WT TASK-1 (n = 17), F125S (n = 28, P = 2.7 × 10−15), F125S/WT (n = 48, P = 4.4 × 10−16), Q126E (n = 15, P = 8.6 × 10−13), Q126E/WT (n = 36, P = 1.4 × 10−18), R131H (n = 22, P = 3.1 × 10−15), and R131H/WT (n = 45, P = 2.1 × 10−13); data are presented as mean ± SD. All mutant currents differ from WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; *** = P < 0.001.

Close modal

Together, these results confirm that these three novel KCNK3 variants also fit the same functional profile as other DDSA variants by producing channels with an increased activity that is markedly less sensitive to inhibition via GqPCR-coupled signalling pathways. These functional defects will therefore have profound effects on the electrical activity of cells in which KCNK3 is expressed, especially during development.

An additional inherited variant in KCNK3 was also identified in 3-year-old child with behavioral problems, later diagnosed as an autism spectrum-like disorder and speech difficulties, but without sleep apnea. Interestingly, this variant (T121M) is also located in the M2 helix adjacent to the pathogenic L122V variant associated with DDSA. However, functional analysis revealed the T121M variant to behave similarly to WT TASK-1, producing only slightly elevated levels of whole-cell currents, especially in the heterozygous state, but with no increased single-channel open probability and only a very small reduction in GqPCR sensitivity (Fig. S2). Thus, although the M2 helix clearly represents a hotspot for DDSA mutations, this variant is unlikely to underlie the phenotype observed in this patient and so also highlights a structural boundary for this specific hotspot of GoF variants.

+ Expand view − Collapse view
Figure S2
Figure S2. Refer to the image caption for details. Panel a shows a membrane topology schematic highlighting the T121M mutation in M2 and the X-gate region. Panel b shows a structural model locating the T121M mutation below the selectivity filter near the X-gate. Panel c shows representative current traces comparing WT and T121M, with larger currents observed in T121M. Panel d shows a current–voltage graph comparing WT and T121M, demonstrating increased currents across positive voltages. Panel e shows bar graphs of normalized current for WT, homomeric T121M, and heterozygous-like channels. Panel f shows a bar graph comparing initial GPCR current ratios, with increased values in T121M. Panel g shows current-versus-time plots after 15 micromolar ATP, comparing responses of WT and T121M channels. Panel h shows representative single-channel recordings of WT and T121M, displaying similar open probability values of 0.02.

Functional effects of the T121M mutation. (a) 2D topological model of a TASK-1 subunit with the position of the T121M mutation labelled in red. The X-gate is labelled in dark grey. (b) TASK-1 crystal structure (PDB: 6RV2). Subunits are shown in teal and orange, respectively. The T121M mutation is shown as red spheres. (c) Representative TEVC recordings of WT TASK-1 and T121M mutant currents in response to voltage steps from −120 to +50 mV in 20 mV steps from a holding potential of −80 mV. (d) Current–voltage plot of WT TASK-1 (n = 12) and T121M (n = 16); data are presented as mean ± SEM. (e) Currents for homomeric and heterozygous T121M channels formed from 1:1 coexpression of WT TASK-1 and T121M normalized to WT current at +50 mV: WT (n = 27), T121M (n = 28, P = 1.5 × 10−10), and T121M-WT (n = 48, P = 5.3 × 10−5). Data are presented as mean ± SD, one-way ANOVA with Dunnett’s multiple comparisons test. (f) Residual current for homomeric T121M after addition of 15 μM ATP. WT TASK-1 (n = 15), and T121M (n = 9, P = 4.3 × 10−5); data are presented as mean ± SD, one-way ANOVA with Dunnett’s multiple comparisons test. ns, not significant; ** = P < 0.01; *** = P < 0.001. (g) Representative currents at +50 mV of WT TASK-1 and T121M channels, over time while adding 15 µM ATP. (h) Single-channel Po values at −100 mV, WT (n = 8) and T121M (n = 8).

Figure S2.
A multi-panel image depicts the functional effects of the T121M mutation. Panel a shows a membrane topology schematic highlighting the T121M mutation in M2 and the X-gate region. Panel b shows a structural model locating the T121M mutation below the selectivity filter near the X-gate. Panel c shows representative current traces comparing WT and T121M, with larger currents observed in T121M. Panel d shows a current–voltage graph comparing WT and T121M, demonstrating increased currents across positive voltages. Panel e shows bar graphs of normalized current for WT, homomeric T121M, and heterozygous-like channels. Panel f shows a bar graph comparing initial GPCR current ratios, with increased values in T121M. Panel g shows current-versus-time plots after 15 micromolar ATP, comparing responses of WT and T121M channels. Panel h shows representative single-channel recordings of WT and T121M, displaying similar open probability values of 0.02.

Functional effects of the T121M mutation. (a) 2D topological model of a TASK-1 subunit with the position of the T121M mutation labelled in red. The X-gate is labelled in dark grey. (b) TASK-1 crystal structure (PDB: 6RV2). Subunits are shown in teal and orange, respectively. The T121M mutation is shown as red spheres. (c) Representative TEVC recordings of WT TASK-1 and T121M mutant currents in response to voltage steps from −120 to +50 mV in 20 mV steps from a holding potential of −80 mV. (d) Current–voltage plot of WT TASK-1 (n = 12) and T121M (n = 16); data are presented as mean ± SEM. (e) Currents for homomeric and heterozygous T121M channels formed from 1:1 coexpression of WT TASK-1 and T121M normalized to WT current at +50 mV: WT (n = 27), T121M (n = 28, P = 1.5 × 10−10), and T121M-WT (n = 48, P = 5.3 × 10−5). Data are presented as mean ± SD, one-way ANOVA with Dunnett’s multiple comparisons test. (f) Residual current for homomeric T121M after addition of 15 μM ATP. WT TASK-1 (n = 15), and T121M (n = 9, P = 4.3 × 10−5); data are presented as mean ± SD, one-way ANOVA with Dunnett’s multiple comparisons test. ns, not significant; ** = P < 0.01; *** = P < 0.001. (g) Representative currents at +50 mV of WT TASK-1 and T121M channels, over time while adding 15 µM ATP. (h) Single-channel Po values at −100 mV, WT (n = 8) and T121M (n = 8).

Close modal

Equivalent disease-causing mutations in both TASK-1 and TASK-3

Due to their high degree of sequence homology (Fig. S1), TASK-1 and TASK-3 share many structural and functional properties, and we have previously shown that GoF variants found in DDSA produce identical functional effects when equivalent mutations are introduced at the same sites in TASK-3 channels (Sörmann et al., 2022).

Interestingly, the R131H DDSA variant reported above has also been identified as a disease-causing GoF variant in TASK-3, where it is associated with KIS (Cousin et al., 2022). This therefore suggests that a pathogenic GoF variant in TASK-1 also has the potential to cause disease if the equivalent mutation occurs in TASK-3 and vice-versa. Furthermore, we also identified an affected individual who fits the classical KIS developmental phenotype with a heterozygous L122V variant in KCNK9, demonstrating that this variant not only results in DDSA in TASK-1 (KCNK3) (Sörmann et al., 2022) but also the KIS neurodevelopmental phenotype if mutated in TASK-3 (KCNK9) (Table S1).

Mutations at this particular site (TM2.6) in the M2 helix have been shown to cause an activatory GoF effect in almost every known K2P channel (Ben Soussia et al., 2019), and in TASK-1 produce channels with increased channel open probability and reduced GqPCR-sensitivity (Sörmann et al., 2022). We therefore examined its effect in TASK-3 and found similar activatory effects that result from an increased single channel open probability (L122V TASK-3 Po = 0.22 ± 0.06 (n = 8) compared with 0.058 ± 0.022 (n = 8) for WT TASK-3; Fig. 3). The fold activation produced by this mutation in homomeric TASK-3 is lower (∼4-fold) compared with its ∼10-fold effect on TASK-1 (Sörmann et al., 2022), but this is likely explained by the higher intrinsic Po of WT TASK-3 compared with WT TASK-1 (Jouen-Tachoire et al., 2026).

Figure 3.
A multi-panel image depicts the effects of TASK channel mutations on current and GPCR inhibition. Panel a shows line traces of single-channel currents at 160 millivolts for TASK-3 wild-type and TASK-3 L122V. The x-axis represents time in milliseconds, and the y-axis represents current in picoamperes. Panel b consists of bar plots displaying normalized currents for homomeric TASK-3 wild-type, homomeric TASK-3 L122V, and heteromeric TASK-3 wild-type plus TASK-3 L122V channels at positive 50 millivolts. The x-axis labels the different channel configurations, and the y-axis shows normalized current values. Panel c presents bar plots of residual currents for homomeric TASK-3 wild-type, homomeric TASK-3 L122V, and heteromeric TASK-3 wild-type plus TASK-3 L122V channels after G protein-coupled receptor-mediated inhibition with 15 micromolar ATP. The x-axis labels the different channel configurations, and the y-axis shows the ratio of residual current to initial current. Each bar plot includes data points and error bars representing the mean and standard deviation. Significant differences are indicated with asterisks.

Position-independent effect of a single L122V mutation in heteromeric TASK channels. (a) Recordings of single channel currents at –160 mV showing the marked increase in Po for TASK-3 L122V channels compared with WT. (b) Currents for homomeric and heteromeric WT and L122V TASK channels normalized to WT current at +50 mV. WT TASK-1 (n = 8), TASK-1 L122V (n = 11, P = 3.2 × 10−6), WT TASK-3 (n = 11), and TASK-3 L122V (n = 10, P = 2.0 × 10−7). For heteromeric TASK-1/TASK-3 channels: top row = TASK-1, bottom = TASK-3. WT-WT (n = 12), WT-L122V (n = 11, P = 0.0063), L122V-WT (n = 12, P = 0.0019), and L122V-L122V (n = 10, P = 0.00013). (c) Residual currents for homomeric and heteromeric TASK channels with different combinations of L122V mutations as indicated after GPCR-mediated inhibition. Residual currents measured after addition of 15 μM ATP. Labels same as shown in panel b. Homomeric: WT TASK-1 (n = 8), TASK-1 L122V (n = 9, P = 2.4 × 10−7), WT TASK-3 (n = 15), and TASK-3 L122V (n = 5, P = 6.5 × 10−5); heteromeric TASK-1/TASK-3: WT-WT (n = 17), WT-L122V (n = 23, P = 9.7 × 10−14), L122V-WT (n = 19, P = 6.7 × 10-13), and L122V-L122V (n = 20, P = 6.7 × 10−16). Data are presented as mean ± SD All mutant currents differ from their respective WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; ** = P < 0.01; *** = P < 0.001. One-way ANOVA with Tukey’s post hoc test used to show the position-independent effect of L122V in heteromeric TASK channels (all ns compared with each other P ≥ 0.2 in each case).

Figure 3.
A multi-panel image depicts the effects of TASK channel mutations on current and GPCR inhibition. Panel a shows line traces of single-channel currents at 160 millivolts for TASK-3 wild-type and TASK-3 L122V. The x-axis represents time in milliseconds, and the y-axis represents current in picoamperes. Panel b consists of bar plots displaying normalized currents for homomeric TASK-3 wild-type, homomeric TASK-3 L122V, and heteromeric TASK-3 wild-type plus TASK-3 L122V channels at positive 50 millivolts. The x-axis labels the different channel configurations, and the y-axis shows normalized current values. Panel c presents bar plots of residual currents for homomeric TASK-3 wild-type, homomeric TASK-3 L122V, and heteromeric TASK-3 wild-type plus TASK-3 L122V channels after G protein-coupled receptor-mediated inhibition with 15 micromolar ATP. The x-axis labels the different channel configurations, and the y-axis shows the ratio of residual current to initial current. Each bar plot includes data points and error bars representing the mean and standard deviation. Significant differences are indicated with asterisks.

Position-independent effect of a single L122V mutation in heteromeric TASK channels. (a) Recordings of single channel currents at –160 mV showing the marked increase in Po for TASK-3 L122V channels compared with WT. (b) Currents for homomeric and heteromeric WT and L122V TASK channels normalized to WT current at +50 mV. WT TASK-1 (n = 8), TASK-1 L122V (n = 11, P = 3.2 × 10−6), WT TASK-3 (n = 11), and TASK-3 L122V (n = 10, P = 2.0 × 10−7). For heteromeric TASK-1/TASK-3 channels: top row = TASK-1, bottom = TASK-3. WT-WT (n = 12), WT-L122V (n = 11, P = 0.0063), L122V-WT (n = 12, P = 0.0019), and L122V-L122V (n = 10, P = 0.00013). (c) Residual currents for homomeric and heteromeric TASK channels with different combinations of L122V mutations as indicated after GPCR-mediated inhibition. Residual currents measured after addition of 15 μM ATP. Labels same as shown in panel b. Homomeric: WT TASK-1 (n = 8), TASK-1 L122V (n = 9, P = 2.4 × 10−7), WT TASK-3 (n = 15), and TASK-3 L122V (n = 5, P = 6.5 × 10−5); heteromeric TASK-1/TASK-3: WT-WT (n = 17), WT-L122V (n = 23, P = 9.7 × 10−14), L122V-WT (n = 19, P = 6.7 × 10-13), and L122V-L122V (n = 20, P = 6.7 × 10−16). Data are presented as mean ± SD All mutant currents differ from their respective WT (one-way ANOVA with Dunnett’s multiple comparisons test). ns, not significant; ** = P < 0.01; *** = P < 0.001. One-way ANOVA with Tukey’s post hoc test used to show the position-independent effect of L122V in heteromeric TASK channels (all ns compared with each other P ≥ 0.2 in each case).

Close modal

The activatory effect of this L122V variant in TASK-3 suggests it may also have an equivalent effect on heteromeric TASK-1/TASK-3 channels, similar to the effect of the L122V variant in TASK-1 (Sörmann et al., 2022). We therefore examined the functional effects of the L122V mutation in heteromeric TASK-1/TASK-3 channels and found that whole-cell currents were increased ∼50% (Fig. 3 b) along with a markedly impaired sensitivity to GqPCR-mediated inhibition (Fig. 3 c) irrespective of which TASK subunit is mutated.

Structural consequences of the TASK-3 L122V variant

The common activatory effect of TM2.6 mutations on K2P channel activity suggests a conserved gating mechanism in this region (Ben Soussia et al., 2019). In a previous study, we demonstrated that hydrophilic substitutions at this position reduced a hydrophobic barrier deep within the inner cavity of TWIK-1 (KCNK1) (Aryal et al., 2014), but it remains unclear whether this is the case in all other K2P channels.

We have recently determined a structure of the WT TASK-3 channel by single-particle cryo-EM (Hall et al., 2025) and so used an identical method to determine the structure of this L122V variant to 2.8 Å resolution (Fig.4, Fig. S3, and Table S2). Comparison of this structure with WT TASK-3 reveals an identical overall fold, i.e., the channel is in the closed X-gate conformation with no obvious differences except within the inner cavity at the site of the mutation. In WT TASK-3, the closed X-gate forms a small inner cavity just below the filter with the L122 side chain lining this cavity (Fig. 4). In the L122V mutant structure, the smaller nature of the side chain increases the radius of the pore at this point from 3 to 4 Å and thus also increases the overall volume of the cavity immediately below the filter. This change is therefore perfectly situated to affect the interaction between K+ and water immediately prior to entry of K+ ions into the filter from within the inner cavity, as well as any hydrophobic gate which exists in this narrow cavity, but it may also have allosteric effects on the filter gate itself.

Figure 4.
A multi-panel image shows structural and pore radius analysis of the TASK-3 L122V mutant channel. Panel a shows a three-dimensional structural model of the TASK-3 L122V potassium channel. Panel b provides an expanded view of the inner cavity at the site of the L122V mutation, shown as a surface representation. The mutant L122V is in purple and the wild-type TASK-3 is in green. The side chains at position 122 are shown as sticks, highlighting the expanded inner cavity in the mutant channel. Panel c presents a pore radius plot comparing wild-type TASK-3 and TASK-3 L122V. The x-axis represents the pore radius in angstroms, and the y-axis represents the distance along the pore in angstroms. The expansion of the inner cavity below the filter at position 122 is indicated.

Cryo-EM structure of the TASK-3 L122V variant reveals expansion of inner cavity below filter. (a) Cartoon representation of the overall fold of the TASK-3 L122V mutant channel as determined by cryo-EM. K+ ions are shown as gold spheres and cholesterol hemisuccinate molecules as gold sticks. (b) Expanded view of the inner cavity at the site of the L122V mutation shown as a surface representation. The mutant L122V is shown in purple and WT TASK-3 in green. The side chains at position 122 are shown as sticks. Note the expanded inner cavity in the mutant channel. (c) Pore radius plot of WT TASK-3 and TASK-3 L122V. The expansion of the inner cavity below the filter at position 122 is indicated.

Figure 4.
A multi-panel image shows structural and pore radius analysis of the TASK-3 L122V mutant channel. Panel a shows a three-dimensional structural model of the TASK-3 L122V potassium channel. Panel b provides an expanded view of the inner cavity at the site of the L122V mutation, shown as a surface representation. The mutant L122V is in purple and the wild-type TASK-3 is in green. The side chains at position 122 are shown as sticks, highlighting the expanded inner cavity in the mutant channel. Panel c presents a pore radius plot comparing wild-type TASK-3 and TASK-3 L122V. The x-axis represents the pore radius in angstroms, and the y-axis represents the distance along the pore in angstroms. The expansion of the inner cavity below the filter at position 122 is indicated.

Cryo-EM structure of the TASK-3 L122V variant reveals expansion of inner cavity below filter. (a) Cartoon representation of the overall fold of the TASK-3 L122V mutant channel as determined by cryo-EM. K+ ions are shown as gold spheres and cholesterol hemisuccinate molecules as gold sticks. (b) Expanded view of the inner cavity at the site of the L122V mutation shown as a surface representation. The mutant L122V is shown in purple and WT TASK-3 in green. The side chains at position 122 are shown as sticks. Note the expanded inner cavity in the mutant channel. (c) Pore radius plot of WT TASK-3 and TASK-3 L122V. The expansion of the inner cavity below the filter at position 122 is indicated.

Close modal
+ Expand view − Collapse view
Figure S3
Figure S3. Refer to the image caption for details. Panel a: The workflow starts with 13,035 movies undergoing motion correction and Patch CTF. A blob picker selects 378 movies, yielding 148,071 particles. These particles undergo two rounds of 2D classification with 35 and 25 classes, respectively. Ab-initio reconstruction follows, producing templates. A template picker then selects 8,381,007 particles, which are re-extracted at 256 pixels and undergo further 2D classification. Heterogeneous refinement is performed on 1,724,459 particles, resulting in ab-initio models. Homogenous refinement, NU-refinement, and reference-based motion correction are applied to 331,817 particles, achieving a resolution of 2.90 angstroms. CryoSieve is then used on 211,355 particles, resulting in a final resolution of 2.83 angstroms. Panel b: A molecular structure of TASK-3 L122V with a color gradient indicating local resolution, ranging from 2.50 to 3.00 angstroms. Panel c: A line graph showing the Gold-standard Fourier Shell Correlation (FSC) curve used for global-resolution estimates. The x-axis represents resolution in angstroms, and the y-axis represents the FSC value. Different lines represent various masking conditions: No Mask, Loose, Tight, and Corrected, with corresponding resolutions of 3.4, 3.1, 2.84, and 2.84 angstroms, respectively.

Determination of cryo-EM Structure of the TASK-3 L122V mutation. (a) Image processing workflow for TASK-3 L122V. (b) Local resolution of reconstructed map as determined within cryoSPARC. (c) Gold-standard FSC curve used for global-resolution estimates within cryoSPARC.

Figure S3.
A multi-panel image of cryo-EM structure determination process. Panel a: The workflow starts with 13,035 movies undergoing motion correction and Patch CTF. A blob picker selects 378 movies, yielding 148,071 particles. These particles undergo two rounds of 2D classification with 35 and 25 classes, respectively. Ab-initio reconstruction follows, producing templates. A template picker then selects 8,381,007 particles, which are re-extracted at 256 pixels and undergo further 2D classification. Heterogeneous refinement is performed on 1,724,459 particles, resulting in ab-initio models. Homogenous refinement, NU-refinement, and reference-based motion correction are applied to 331,817 particles, achieving a resolution of 2.90 angstroms. CryoSieve is then used on 211,355 particles, resulting in a final resolution of 2.83 angstroms. Panel b: A molecular structure of TASK-3 L122V with a color gradient indicating local resolution, ranging from 2.50 to 3.00 angstroms. Panel c: A line graph showing the Gold-standard Fourier Shell Correlation (FSC) curve used for global-resolution estimates. The x-axis represents resolution in angstroms, and the y-axis represents the FSC value. Different lines represent various masking conditions: No Mask, Loose, Tight, and Corrected, with corresponding resolutions of 3.4, 3.1, 2.84, and 2.84 angstroms, respectively.

Determination of cryo-EM Structure of the TASK-3 L122V mutation. (a) Image processing workflow for TASK-3 L122V. (b) Local resolution of reconstructed map as determined within cryoSPARC. (c) Gold-standard FSC curve used for global-resolution estimates within cryoSPARC.

Close modal

In this study, we have identified six further DDSA patients, with GoF variants in KCNK3, including three with novel GoF missense variants (F125S, Q126E, and R131H). This therefore markedly expands the global cohort of individuals with DDSA and reinforces the role that defective TASK-1 channels play in this unique neurodevelopmental disorder. We also report a novel TASK-3 variant (L122V) in an individual with KIS and have functionally and structurally characterized this mutant channel. These variants all share similar GoF effects on channel activity with a markedly reduced sensitivity to inhibition via GqPCR-coupled pathways. These GoF effects and their regulatory defects were also dominant in heterozygous channels containing one mutant subunit, as well as in heteromeric TASK-1/TASK-3 channels. Furthermore, we show that two of these pathogenic variants (R131H and L122V) are found in identical conserved positions in both KCNK3 (TASK-1) and KCNK9 (TASK-3) channels, and that their effects on heteromeric TASK channel activity appear the same, irrespective of which channel subunit contains the mutation.

Other examples also exist of equivalent pathogenic mutations in conserved domains of related channels that result in similar functional defects, yet different phenotypes, particularly within the superfamily of tetrameric cation channels. For example, the mutation of equivalent residues within the highly conserved DIV S3–S4 linker of the skeletal and cardiac muscle voltage-gated Na+ channels (SCN4A and SCN5A, respectively) produce similar effects on channel inactivation in these channels, yet result in different channelopathies (Mantegazza et al., 2021). Likewise, equivalent mutations within the positive charges of the conserved voltage sensor S4-helix in voltage-gated K+ channels also result in different channelopathies. In both cases, the unique phenotypes are due to the different spatiotemporal expression of these various channels, and so the distinct expression patterns for TASK-1 and TASK-3 will also heavily influence their phenotype. However, for the equivalent TASK channel mutations described here (e.g., L122V and R131H), the overlapping expression patterns of TASK-1 and TASK-3 channels and their ability to form physiologically relevant heteromeric channels is also important because any differences in the activity of the heteromeric channels that result from these mutations could be highly significant.

However, our results show that these equivalent mutations in either TASK-1 or TASK-3 have very similar effects on heteromeric channel activity. Given the very high degree of structural and functional conservation between TASK-1 and TASK-3 (Fig. S1), this is perhaps not surprising, but the markedly different disease phenotypes that result from these mutations therefore suggest it is defective homomeric TASK channels and their distinct spatiotemporal patterns of expression that may dominate the unique features of each disorder, especially within the developing CNS (Aller and Wisden, 2008). Such GoF mutations will not only dampen neuronal excitability, but it is well known that inappropriate spatiotemporal K+ channel expression can adversely affect the development of many cell types, especially neurons (D'Adamo et al., 2020).

Intriguingly, not all DDSA-associated variants (e.g., L239P) (Sörmann et al., 2022) result in a major increase in whole-cell currents even though they all increase channel Po. Indeed, the increase in currents for the three novel variants reported here, especially in the heterozygous state (Fig. 1 e), is relatively modest. Furthermore, in KIS some of the variants result in reduced TASK-3 channel activity (Cousin et al., 2022). However, in every case, these pathogenic mutants all share a marked reduction in their sensitivity to GqPCR-mediated inhibition. Therefore, in addition to any change in whole-cell currents, this loss of GqPCR sensitivity may represent a significant factor in defining their pathogenic phenotype. In particular, this may be especially important in the respiratory phenotype seen in all DDSA (and some KIS) individuals because GqPCR-mediated inhibition of TASK channels contributes to depolarization of cells in the carotid body as well as in pacemaker neurons in the brainstem important for sleep/wakefulness transition and respiratory rhythmogenesis. Abnormal respiratory control would therefore be expected to result from such defective channel regulation (Buckler et al., 2000; Trapp et al., 2008; Bayliss et al., 2015; Buckler, 2015).

The genomic imprinting reported for KCNK9 and the mosaicism that arises from their de novo origin may also explain why only a subset of KIS patients are reported to have sleep apnea, compared with its defining presence in all known cases of DDSA (Luedi et al., 2007; Cooper et al., 2020; Cousin et al., 2022). However, there is also precedent in other genetic conditions involving defective K+ channels where, paradoxically, both LoF and GoF variants lead to neuronal hyperexcitability (Niday and Tzingounis, 2018) and/or neurodevelopmental disorders (D'Adamo et al., 2020). Consequently, much work remains to be done in order to understand this complexity.

Finally, our cryo-EM structure of the L122V variant at the TM2.6 site in TASK-3 also highlights how even subtle changes in structure can have dramatic functional consequences. Although the hydrophobicity of this side chain does not change, the slight increase in volume of the inner cavity just below the filter may be enough to affect a critical step in permeation, i.e., the dehydration of the K+ ion prior to its entry into the filter and to affect any hydrophobic gating mechanism in this region (Aryal et al., 2015). It is also consistent with the activatory effect of hydrophilic substitutions at this site (Aryal et al., 2014; Ben Soussia et al., 2019). Indeed, recent detailed simulations of K+ permeation in open state conformations of other K2P channels have shown that K+ permeation through the filter is tightly coupled to gating (Aldakul et al., 2025). This particular mutation therefore highlights the complexity of understanding the effect of missense variants on channel activity and their relationship to clinical phenotype, even when detailed functional assays and structural information are available.

In conclusion, we have characterized the functional and structural properties of a range of novel GoF missense disease-causing variants in TASK channels to expand the molecular and genetic basis of TASK channelopathies. These variants also provide insight into the relative contributions of homomeric versus heteromeric TASK channels in their respective disorders. Overall, our characterization of the structural and functional effects of these pathogenic variants has important implications for understanding genotype/phenotype correlations and also for future drug development studies involving both these rare channelopathies as well as other more common disease states involving TASK channels.

The data underlying Figs. 1, 2 and 3 are available in the published article and its online supplemental material. The atomic coordinates for the TASK-3 L122V model shown in Fig. 4 have been deposited in the PDB with accession number 28IZ. Coulomb potential maps have been deposited in the EMDB database, under the following accession code: EMDB-56538. Patient data included in Table S1 from the National Genomic Research Library (NGRL) used in this research are available within the secure Genomics England Research Environment. Access to NGRL data is restricted to adhere to consent requirements and protect participant privacy. Visit: https://www.genomicsengland.co.uk/research for more information on patient data access.

Christopher J. Lingle served as editor.

We are also grateful for the important support from patients and families, our UK and international collaborators, brainbank, and biobanks, and grateful for other essential funding from the MSA Trust, UK Dementia Research Institute, National Institute for Health Research University College London Hospitals Biomedical Research Centre (NIHR-BRC), Michael J. Fox Foundation, Fidelity Foundation, Rosetrees Trust, EAN, ERDERA: European Union’s Horizon Europe research and innovation programme, Dolby Family Fund, Alzheimer’s Research UK, Mission MSA, Parkinson’s disease UK, Parkinson’s Foundation, Muscular Dystrophy UK, Ataxia UK, CureDRPLA, ALS Association, National Ataxia Foundation, Target ALS Foundation, Medical Research Foundation, and the National Brain Appeal. JF and CFW were supported by the National Institute for Health and Care Research Exeter Biomedical Research Centre. We also gratefully acknowledge the participants of the National Genomic Research Library (NGRL). Secure access to the NGRL under project ID RR349 was provided by Genomics England in partnership with NHS England. The NGRL contains participants’ health data collected by the NHS as part of their care, along with samples and data from their participation in research, for which fully informed consent has been obtained. This includes genomic and clinical data provided through the NHS Genomic Medicine Service. The views expressed are those of the authors and not necessarily those of NHS England, the NIHR, or the Department of Health and Social Care. We also thank the Central Oxford Structural Molecular Imaging Centre (COSMIC) cryo-EM facility for data collection, assistance, and technical support.

This work was directly supported by grants from the Biotechnology and Biological Sciences Research Council and Medical Research Council to S.J. Tucker (BB/T002018/1, BB/S008608/1, and MR/W017741/1). It was also supported by the Wellcome Trust as part of the OXION Initiative (WT084655MA and 102161/B/13/Z). K.M. Crowther was funded by the UKRI-BBSRC Interdisciplinary Bioscience Doctoral Training Partnership (BB/T008784/1). Open Access funding provided by University of Oxford.

Author contributions: Kate M. Crowther: conceptualization, data curation, formal analysis, investigation, methodology, resources, validation, visualization, and writing—original draft, review, and editing. Thibault R.H. Jouen-Tachoire: data curation, formal analysis, investigation, methodology, visualization, and writing—original draft, review, and editing. Peter Proks: formal analysis, investigation, and visualization. Peter Rory Hall: Data curation, investigation, visualization, and writing—review and editing. Emma L. Veale: resources, validation, and writing—review and editing. Janina Sörmann: writing—review and editing. Karin E.J. Rödström: formal analysis and writing—review and editing. Thomas Müller: conceptualization, resources, and writing—review and editing. Saskia B. Wortmann: data curation, formal analysis, investigation, and writing—review and editing. Nina Barisic: investigation and writing—review and editing. Natalie Hauser: resources. Vincenzo Damiano Salpietro: investigation. RaeLynn Forsyth: conceptualization, data curation, resources, and writing—original draft, review, and editing. Linford Williams: data curation and resources. Nora Derrabi: resources. Carlos A. Bacino: resources and writing—review and editing. Jill A. Rosenfeld: investigation, resources, and writing—review and editing. Henry Houlden: conceptualization, data curation, funding acquisition, investigation, project administration, resources, and writing—review and editing. Simon Newstead: supervision. Caroline F. Wright: supervision and writing—review and editing. James Fasham: investigation, and writing—review and editing. Alistair A. Mathie: conceptualization, data curation, validation, and writing—review and editing. Reza Maroofian: data curation, formal analysis, investigation, methodology, resources, validation, and writing—review and editing. Stephen J. Tucker: conceptualization, formal analysis, funding acquisition, project administration, resources, supervision, validation, visualization, and writing—original draft, review, and editing.

Adams
,
P.D.
,
P.V.
Afonine
,
G.
Bunkóczi
,
V.B.
Chen
,
I.W.
Davis
,
N.
Echols
,
J.J.
Headd
,
L.W.
Hung
,
G.J.
Kapral
,
R.W.
Grosse-Kunstleve
, et al
.
2010
.
PHENIX: A comprehensive Python-based system for macromolecular structure solution
.
Acta Crystallogr. D Biol. Crystallogr.
66
:
213
221
.
Aldakul
,
Y.K.
,
M.
Schewe
,
C.
Coll-Diez
,
S.
Hwang
,
T.
Baukrowitz
, and
H.
Sun
.
2025
.
Atomistic mechanism of noncanonical voltage gating in K2P channels
.
Sci. Adv.
11
:eadx1680.
Aller
,
M.I.
, and
W.
Wisden
.
2008
.
Changes in expression of some two-pore domain potassium channel genes (KCNK) in selected brain regions of developing mice
.
Neuroscience
.
151
:
1154
1172
.
Aryal
,
P.
,
F.
Abd-Wahab
,
G.
Bucci
,
M.S.P.
Sansom
, and
S.J.
Tucker
.
2014
.
A hydrophobic barrier deep within the inner pore of the TWIK-1 K2P potassium channel
.
Nat. Commun.
5
:
4377
.
Aryal
,
P.
,
M.S.
Sansom
, and
S.J.
Tucker
.
2015
.
Hydrophobic gating in ion channels
.
J. Mol. Biol.
427
:
121
130
.
Bagal
,
S.K.
,
A.D.
Brown
,
P.J.
Cox
,
K.
Omoto
,
R.M.
Owen
,
D.C.
Pryde
,
B.
Sidders
,
S.E.
Skerratt
,
E.B.
Stevens
,
R.I.
Storer
, and
N.A.
Swain
.
2013
.
Ion channels as therapeutic targets: A drug discovery perspective
.
J. Med. Chem.
56
:
593
624
.
Barel
,
O.
,
S.A.
Shalev
,
R.
Ofir
,
A.
Cohen
,
J.
Zlotogora
,
Z.
Shorer
,
G.
Mazor
,
G.
Finer
,
S.
Khateeb
,
N.
Zilberberg
, and
O.S.
Birk
.
2008
.
Maternally inherited Birk Barel mental retardation dysmorphism syndrome caused by a mutation in the genomically imprinted potassium channel KCNK9
.
Am. J. Hum. Genet.
83
:
193
199
.
Bayliss
,
D.A.
,
J.
Barhanin
,
C.
Gestreau
, and
P.G.
Guyenet
.
2015
.
The role of pH-sensitive TASK channels in central respiratory chemoreception
.
Pflugers Arch.
467
:
917
929
.
Ben Soussia
,
I.
,
S.
El Mouridi
,
D.
Kang
,
A.
Leclercq-Blondel
,
L.
Khoubza
,
P.
Tardy
,
N.
Zariohi
,
M.
Gendrel
,
F.
Lesage
,
E.J.
Kim
, et al
.
2019
.
Mutation of a single residue promotes gating of vertebrate and invertebrate two-pore domain potassium channels
.
Nat. Commun.
10
:
787
.
Berg
,
A.P.
,
E.M.
Talley
,
J.P.
Manger
, and
D.A.
Bayliss
.
2004
.
Motoneurons express heteromeric TWIK-related acid-sensitive K+ (TASK) channels containing TASK-1 (KCNK3) and TASK-3 (KCNK9) subunits
.
J. Neurosci.
24
:
6693
6702
.
Bohnen
,
M.S.
,
D.
Roman-Campos
,
C.
Terrenoire
,
J.
Jnani
,
K.J.
Sampson
,
W.K.
Chung
, and
R.S.
Kass
.
2017
.
The impact of heterozygous KCNK3 mutations associated with pulmonary arterial hypertension on channel function and pharmacological recovery
.
J. Am. Heart Assoc.
6
:e006465.
Buckler
,
K.J.
2015
.
TASK channels in arterial chemoreceptors and their role in oxygen and acid sensing
.
Pflugers Arch.
467
:
1013
1025
.
Buckler
,
K.J.
,
B.A.
Williams
, and
E.
Honore
.
2000
.
An oxygen-acid- and anaesthetic-sensitive TASK-like background potassium channel in rat arterial chemoreceptor cells
.
J. Physiol.
525
:
135
142
.
Conway
,
K.E.
, and
J.F.
Cotten
.
2012
.
Covalent modification of a volatile anesthetic regulatory site activates TASK-3 (KCNK9) tandem-pore potassium channels
.
Mol. Pharmacol.
81
:
393
400
.
Cooper
,
A.
,
T.
Butto
,
N.
Hammer
,
S.
Jagannath
,
D.L.
Fend-Guella
,
J.
Akhtar
,
K.
Radyushkin
,
F.
Lesage
,
J.
Winter
,
S.
Strand
, et al
.
2020
.
Inhibition of histone deacetylation rescues phenotype in a mouse model of Birk-Barel intellectual disability syndrome
.
Nat. Commun.
11
:
480
.
Cousin
,
M.A.
,
E.L.
Veale
,
N.R.
Dsouza
,
S.
Tripathi
,
R.G.
Holden
,
M.
Arelin
,
G.
Beek
,
M.R.
Bekheirnia
,
J.
Beygo
,
V.
Bhambhani
, et al
.
2022
.
Gain and loss of TASK3 channel function and its regulation by novel variation cause KCNK9 imprinting syndrome
.
Genome Med.
14
:
62
.
Croll
,
T.I.
2018
.
Isolde: A physically realistic environment for model building into low-resolution electron-density maps
.
Acta Crystallogr. D Struct. Biol.
74
:
519
530
.
D'Adamo
,
M.C.
,
A.
Liantonio
,
E.
Conte
,
M.
Pessia
, and
P.
Imbrici
.
2020
.
Ion channels involvement in neurodevelopmental disorders
.
Neuroscience
.
440
:
337
359
.
Emsley
,
P.
,
B.
Lohkamp
, and
W.G.
Scott
.
2010
.
Features and development of Coot
.
Acta Crystallogr. D Biol. Crystallogr.
66
:
486
501
.
Enyedi
,
P.
, and
G.
Czirjak
.
2010
.
Molecular background of leak K+ currents: Two-pore domain potassium channels
.
Physiol. Rev.
90
:
559
605
.
Graham
,
J.M.
, Jr.
,
N.
Zadeh
,
M.
Kelley
,
E.S.
Tan
,
W.
Liew
,
V.
Tan
,
M.A.
Deardorff
,
G.N.
Wilson
,
L.
Sagi-Dain
, and
S.A.
Shalev
.
2016
.
KCNK9 imprinting syndrome-further delineation of a possible treatable disorder
.
Am. J. Med. Genet. A.
170
:
2632
2637
.
Hall
,
P.R.
,
T.
Jouen-Tachoire
,
M.
Schewe
,
P.
Proks
,
T.
Baukrowitz
,
E.P.
Carpenter
,
S.
Newstead
,
K.E.J.
Rödström
, and
S.J.
Tucker
.
2025
.
Structures of TASK-1 and TASK-3 K2P channels provide insight into their gating and dysfunction in disease
.
Structure
.
33
:
115
122 e114
.
Jouen-Tachoire
,
T.R.H.
,
P.
Proks
,
D.
Seiferth
,
K.
Crowther
,
P.C.
Biggin
,
T.
Baukrowitz
,
M.
Schewe
, and
S.J.
Tucker
.
2026
.
Structural determinants for GPCR-mediated inhibition of TASK K2P channels by diacylglycerol and its dysfunction in disease
.
EMBO J.
45
:
2400
2412
.
Jouen-Tachoire
,
T.R.H.
,
S.J.
Tucker
, and
P.
Tammaro
.
2021
.
Ion channels as convergence points in the pathology of pulmonary arterial hypertension
.
Biochem. Soc. Trans.
49
:
1855
1865
.
Kaplanis
,
J.
,
K.E.
Samocha
,
L.
Wiel
,
Z.
Zhang
,
K.J.
Arvai
,
R.Y.
Eberhardt
,
G.
Gallone
,
S.H.
Lelieveld
,
H.C.
Martin
,
J.F.
McRae
, et al
.
2020
.
Evidence for 28 genetic disorders discovered by combining healthcare and research data
.
Nature
.
586
:
757
762
.
Kim
,
D.
,
E.J.
Cavanaugh
,
I.
Kim
, and
J.L.
Carroll
.
2009
.
Heteromeric TASK-1/TASK-3 is the major oxygen-sensitive background K+ channel in rat carotid body glomus cells
.
J. Physiol.
587
:
2963
2975
.
Lee
,
L.M.
,
T.
Müntefering
,
T.
Budde
,
S.G.
Meuth
, and
T.
Ruck
.
2021
.
Pathophysiological role of K2P channels in human diseases
.
Cell Physiol. Biochem.
55
:
65
86
.
Lin
,
A.Y.
,
C.D.
Turnbull
, and
J.J.
Pandit
.
2025
.
The “TASK” of breathing: Anesthetic relevance of background two-pore domain potassium channels as therapeutic targets for respiratory control
.
Anesth. Analg.
140
:
1414
1425
.
Lin
,
H.
,
J.
Li
,
Q.
Zhang
,
H.
Yang
, and
S.
Chen
.
2024
.
C-type inactivation and proton modulation mechanisms of the TASK3 channel
.
Proc. Natl. Acad. Sci. USA
.
121
:e2320345121.
Luedi
,
P.P.
,
F.S.
Dietrich
,
J.R.
Weidman
,
J.M.
Bosko
,
R.L.
Jirtle
, and
A.J.
Hartemink
.
2007
.
Computational and experimental identification of novel human imprinted genes
.
Genome Res.
17
:
1723
1730
.
Luethy
,
A.
,
J.D.
Boghosian
,
R.
Srikantha
, and
J.F.
Cotten
.
2017
.
Halogenated ether, alcohol, and alkane anesthetics activate TASK-3 tandem pore potassium channels likely through a common mechanism
.
Mol. Pharmacol.
91
:
620
629
.
Ma
,
L.
,
D.
Roman-Campos
,
E.D.
Austin
,
M.
Eyries
,
K.S.
Sampson
,
F.
Soubrier
,
M.
Germain
,
D.A.
Tregouet
,
A.
Borczuk
,
E.B.
Rosenzweig
, et al
.
2013
.
A novel channelopathy in pulmonary arterial hypertension
.
N. Engl. J. Med.
369
:
351
361
.
Mantegazza
,
M.
,
S.
Cestele
, and
W.A.
Catterall
.
2021
.
Sodium channelopathies of skeletal muscle and brain
.
Physiol. Rev.
101
:
1633
1689
.
Mathie
,
A.
,
E.L.
Veale
,
K.P.
Cunningham
,
R.G.
Holden
, and
P.D.
Wright
.
2021
.
Two-pore domain potassium channels as drug targets: Anesthesia and beyond
.
Annu. Rev. Pharmacol. Toxicol.
61
:
401
420
.
Meng
,
E.C.
,
T.D.
Goddard
,
E.F.
Pettersen
,
G.S.
Couch
,
Z.J.
Pearson
,
J.H.
Morris
, and
T.E.
Ferrin
.
2023
.
UCSF ChimeraX: Tools for structure building and analysis
.
Protein Sci.
32
:e4792.
Niday
,
Z.
, and
A.V.
Tzingounis
.
2018
.
Potassium Channel gain of function in epilepsy: An unresolved paradox
.
Neuroscientist
.
24
:
368
380
.
Niemeyer
,
M.I.
,
L.P.
Cid
,
W.
González
, and
F.V.
Sepulveda
.
2016
.
Gating, regulation, and structure in K2P K+ channels: In varietate Concordia?
Mol. Pharmacol.
90
:
309
317
.
Olschewski
,
A.
,
E.L.
Veale
,
B.M.
Nagy
,
C.
Nagaraj
,
G.
Kwapiszewska
,
F.
Antigny
,
M.
Lambert
,
M.
Humbert
,
G.
Czirják
,
P.
Enyedi
, et al
.
2017
.
TASK-1 (KCNK3) channels in the lung: From cell biology to clinical implications
.
Eur. Respir. J.
50
:
1700754
.
Osman
,
A.M.
,
B.
Toson
,
G.R.
Naik
,
S.
Mukherjee
,
M.
Delbeck
,
M.
Hahn
,
T.
Muller
,
G.
Weimann
, and
D.J.
Eckert
.
2024
.
A novel TASK channel antagonist nasal spray reduces sleep apnea severity in physiological responders: A randomized, blinded, trial
.
Am. J. Physiol. Heart Circ. Physiol.
326
:
H715
H723
.
Punjani
,
A.
,
J.L.
Rubinstein
,
D.J.
Fleet
, and
M.A.
Brubaker
.
2017
.
cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination
.
Nat. Methods
.
14
:
290
296
.
Rinné
,
S.
,
F.
Schick
,
K.
Vowinkel
,
S.
Schutte
,
C.
Krasel
,
S.
Kauferstein
,
M.K.
Schafer
,
A.K.
Kiper
,
T.
Muller
, and
N.
Decher
.
2024
.
Potassium channel TASK-5 forms functional heterodimers with TASK-1 and TASK-3 to break its silence
.
Nat. Commun.
15
:
7548
.
Rodstrom
,
K.E.J.
,
A.K.
Kiper
,
W.
Zhang
,
S.
Rinne
,
A.C.W.
Pike
,
M.
Goldstein
,
L.J.
Conrad
,
M.
Delbeck
,
M.G.
Hahn
,
H.
Meier
, et al
.
2020
.
A lower X-gate in TASK channels traps inhibitors within the vestibule
.
Nature
.
582
:
443
447
.
Sörmann
,
J.
,
M.
Schewe
,
P.
Proks
,
T.
Jouen-Tachoire
,
S.
Rao
,
E.B.
Riel
,
K.E.
Agre
,
A.
Begtrup
,
J.
Dean
,
M.
Descartes
, et al
.
2022
.
Gain-of-function mutations in KCNK3 cause a developmental disorder with sleep apnea
.
Nat. Genet.
54
:
1534
1543
.
Trapp
,
S.
,
M.I.
Aller
,
W.
Wisden
, and
A.V.
Gourine
.
2008
.
A role for TASK-1 (KCNK3) channels in the chemosensory control of breathing
.
J. Neurosci.
28
:
8844
8850
.
Zhu
,
J.
,
Q.
Zhang
,
H.
Zhang
,
Z.
Shi
,
M.
Hu
, and
C.
Bao
.
2023
.
A minority of final stacks yields superior amplitude in single-particle cryo-EM
.
Nat. Commun.
14
:
7822
.

This work is part of a special issue on Emerging Research on Ion Channels in Health and disease.

Author notes

*

Kate M. Crowther and Thibault R.H. Jouen-Tachoire contributed equally to this paper.

Disclosures: T. Müller reported being an employee of Bayer AG Pharma, Wuppertal, Germany, which in the past pursued clinical trials with TASK inhibitors to treat sleep apnea. All TASK activities have been stopped, and all related patents been dropped before the initiation of this work. J.A. Rosenfeld reported that the Department of Molecular and Human Genetics at Baylor College of Medicine receives revenue from clinical genetic testing completed at Baylor Genetics Laboratories. No other disclosures were reported.

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

or Create an Account

Close Modal
Close Modal