Physical basis for distinct basal and mechanically gated activity of the human K+ channel TRAAK

Multiple-Basin Go̅-Martini for Investigating Conformational Transitions and Environmental Interactions of Proteins

Proteins are inherently dynamic molecules, and their conformational transitions among various states are essential for numerous biological processes, which are often modulated by their interactions with surrounding environments. Although molecular dynamics (MD) simulations are widely used to investigate these transitions, all-atom (AA) methods are often limited by short time scales and high computational costs, and coarse-grained (CG) implicit-solvent Go̅-like models are usually incapable of studying the interactions between proteins and their environments. Here, we present an approach called Multiple-basin Go̅-Martini, which combines the recent Go̅-Martini model with an exponential mixing scheme to facilitate the simulation of spontaneous protein conformational transitions in explicit environments. We demonstrate the versatility of our method through five diverse case studies: GlnBP, Arc, Hinge, SemiSWEET, and TRAAK, representing ligand-binding proteins, fold-switching proteins, de novo designed proteins, transporters, and mechanosensitive ion channels, respectively. Multiple-basin Go̅-Martini offers a new computational tool for investigating protein conformational transitions, identifying key intermediate states, and elucidating essential interactions between proteins and their environments, particularly protein-membrane interactions. In addition, this approach can efficiently generate thermodynamically meaningful data sets of protein conformational space, which may enhance deep learning-based models for predicting protein conformation distributions.

Gating mechanism of the two-pore-domain potassium channel THIK1

TWIK-related halothane-inhibited potassium channel (THIK1) maintains the resting membrane potential and regulates potassium efflux in microglia. It is a potential therapeutic target for neurodegenerative disorders, neuropathic pain and inflammation. However, the mechanism underlying its function remains unclear. Here we used cryo-electron microscopy to solve the structures of full-length human THIK1, revealing two inner gates and a C-type selectivity filter gate, distinct from other two-pore-domain potassium channels. One inner gate, formed by a short helix in the distal C terminus, introduces a unique gating mechanism involving the distal cytoplasmic domain. The other, beneath the selectivity filter, is constricted by Y273 in the M4 helix, dividing the cavity. In addition, the selectivity filter gate is modulated by polyunsaturated fatty acids. These structural insights into THIK1 gating, through the distal C-terminal helices, hydrophilic residues and selectivity filter, advance our understanding of THIK1's role in microglial homeostasis and neuropathologies.

De novo KCNK4 variant caused epilepsy with febrile seizures plus, neurodevelopmental abnormalities, and hypertrichosis

Background

Epilepsy with febrile seizures plus (EFS+) is a syndrome with a strong genetic component. Previously, variants in several genes encoding ion channels have been associated with EFS+. However, the etiology in the majority of patients remains undetermined.

Methods

Trio-based whole-exome sequencing was performed on a patient with EFS+. Previously reported KCNK4 variants were systemically reviewed to analyze the phenotypic spectrum and core phenotypes.

Results

A novel de novo KCNK4 variant (c.415G>A/p.Gly139Arg) was identified in a patient with EFS+, neurodevelopmental abnormalities, and hypertrichosis. The identified variant was absent in normal populations, indicated to alter hydrogen bonds with surrounding residues by various protein modeling, predicted to be damaging for protein function by twenty algorithms, located in residues of high conservation across species, and classified as pathogenic by the ACMG guidelines. Protein modeling analyses of the variant suggested a possible gain-of-function effect. Analysis of other eight cases with KCNK4 variants outlined the phenotypic spectrums of KCNK4, ranging from mild benign epilepsy, EFS+ with neurodevelopmental abnormalities, to syndromic neurodevelopmental disorders and revealed neurodevelopmental abnormalities and epilepsy as its core phenotypes. Integrated analysis suggested that minor allele frequency and in silico meta-predictors effectively distinguish pathogenic variants.

Conclusion

This study suggested the KCNK4 gene as a novel candidate causative gene of EFS+, which would be helpful for the genetic diagnosis and clinical management of patients.

Structure of the human K2P13.1 channel reveals a hydrophilic pore restriction and lipid cofactor site

Polyunsaturated fatty acid (PUFA) lipids modulate the neuronal and microglial leak potassium channel K2P13.1 (THIK1) and other voltage-gated ion channel (VGIC) superfamily members through poorly understood mechanisms. Here we present cryo-electron microscopy structures of human THIK1 and mutants, revealing a unique two-chamber aqueous inner cavity obstructed by a hydrophilic barrier important for gating, the flow restrictor, and a P1-M4 intersubunit interface lipid at a site, the PUFA site, corresponding to the K2P small-molecule modulator pocket. This overlap, together with functional studies, indicates that PUFA site lipids are THIK1 cofactors. Comparison with a PUFA-responsive VGIC, Kv7.1, reveals a shared modulatory role for the pore domain intersubunit interface, providing a framework for understanding PUFA action on the VGIC superfamily. Our findings reveal the distinct THIK1 architecture, highlight the importance of the P1-M4 interface for K2P control by natural and synthetic ligands and should aid in the development of THIK subfamily modulators for neuroinflammation and autism.

Insights into the structure and modulation of human TWIK-2

The T andem of pore domain in a W eak I nward R ectifying K + channel 2 (TWIK-2; KCNK6) is a member of the Two-Pore Domain K + (K2P) channel family, which is associated with pulmonary hypertension, lung injury, and inflammation. The structure and regulatory mechanisms of TWIK-2 remain largely unknown. Here, we present the cryo-electron microscopy (cryo-EM) structure of human TWIK-2 at ~3.7 Å and highlight its conserved and unique features. Using automated whole-cell patch clamp recordings, we demonstrate that gating in TWIK-2 is voltage-dependent and insensitive to changes in the extracellular pH. We identify key residues that influence TWIK-2 activity by employing structure and sequence-guided site-directed mutagenesis and provide insights into the possible mechanism of TWIK-2 regulation. Additionally, we demonstrate the application of high-throughput automated whole-cell patch clamp platforms to screen small molecule modulators of TWIK-2. Our work serves as a foundation for designing high-throughput small molecule screening campaigns to identify specific high-affinity TWIK-2 modulators, including promising new anti-inflammatory therapeutics.

Programmable Lipid Bilayer Tension-Control Apparatus for Quantitative Mechanobiology

The biological membrane is not just a platform for information processing but also a field of mechanics. The lipid bilayer that constitutes the membrane is an elastic body, generating stress upon deformation, while the membrane protein embedded therein deforms the bilayer through structural changes. Among membrane-protein interplays, various channel species act as tension-current converters for signal transduction, serving as elementary processes in mechanobiology. However, in situ studies in chaotically complex cell membranes are challenging, and characterizing the tension dependency of mechanosensitive channels remains semiquantitative owing to technical limitations. Here, we developed a programmable membrane tension-control apparatus on a lipid bilayer system. This synthetic membrane system [contact bubble bilayer (CBB)] uses pressure to drive bilayer tension changes via the Young-Laplace principle, whereas absolute bilayer tension is monitored in real-time through image analysis of the bubble geometry via the Young principle. Consequently, the mechanical nature of the system permits the implementation of closed-loop feedback control of bilayer tension (tension-clamp CBB), maintaining a constant tension for minutes and allowing stepwise tension changes within a hundred milliseconds in the tension range of 0.8 to 15 mN·m-1. We verified the system performance by examining the single-channel behavior of tension-dependent KcsA and TREK-1 potassium channels under scheduled tension time courses prescribed via visual interfaces. The result revealed steady-state activity and dynamic responses to the step tension changes, which are essential to the biophysical characterization of the channels. The apparatus explores a frontier for quantitative mechanobiology studies and promotes the development of a tension-operating experimental robot.

Effect of two activators on the gating of a K2P channel

TWIK-related potassium channel 1 (TREK1), a two-pore-domain mammalian potassium (K+) channel, regulates the resting potential across cell membranes, presenting a promising therapeutic target for neuropathy treatment. The gating of this channel converges in the conformation of the narrowest part of the pore: the selectivity filter (SF). Various hypotheses explain TREK1 gating modulation, including the dynamics of loops connecting the SF with transmembrane helices and the stability of hydrogen bond (HB) networks adjacent to the SF. Recently, two small molecules (Q6F and Q5F) were reported as activators that affect TREK1 by increasing its open probability in single-channel current measurements. Here, using molecular dynamics simulations, we investigate the effect of these ligands on the previously proposed modulation mechanisms of TREK1 gating compared to the apo channel. Our findings reveal that loop dynamics at the upper region of the SF exhibit only a weak correlation with permeation events/nonpermeation periods, whereas the HB network behind the SF appears more correlated. These nonpermeation periods arise from both distinct mechanisms: a C-type inactivation (resulting from dilation at the top of the SF), which has been described previously, and a carbonyl flipping in an SF binding site. We find that, besides the prevention of C-type inactivation in the channel, the ligands increase the probability of permeation by modulating the dynamics of the carbonyl flipping, influenced by a threonine residue at the bottom of the SF. These results offer insights for rational ligand design to optimize the gating modulation of TREK1 and related K+ channels.

The epilepsy phenotype of KCNK4-related neurodevelopmental disease

INTRODUCTION

Potassium ion channels play a crucial role in maintaining cellular electrical stability and are implicated in various epilepsies. Heterozygous pathogenic variants in KCNK4 cause a recognizable neurodevelopmental syndrome with facial dysmorphism, hypertrichosis, epilepsy, intellectual disability (ID), and gingival overgrowth (FHEIG). To date, no more than nine patients with FHEIG have been described worldwide and still little is known about epileptic phenotype in KCNK4-related disease.

METHODS

We identified a novel de novo p.(Gly139Arg) variant in KCNK4 in a patient with drug-resistant nocturnal seizures, mild ID, and dysmorphic features. In silico analyses of the variant strongly suggest a gain-of-function effect. We conducted a retrospective review of previously published cases, focusing on the epileptic features and response to various treatments.

RESULTS

To date, epilepsy has been reported in 8/10 patients with KCNK4-related disease. The mean age of seizure onset was 1.8 years, and the most common seizure type was focal to bilateral tonic-clonic (5/8). Sodium channel blockers and valproate were effective in the majority of patients, but in 3/8 the epilepsy was drug-resistant. Our patient showed improved seizure control after treatment with the carbonic anhydrase inhibitor sulthiame. Interestingly, the patient showed features of peripheral nerve hyperexcitability syndrome, a phenomenon not previously described in potassium channelopathies caused by increased K+ conductance.

CONCLUSION

Gain-of-function variants in KCNK4 cause a spectrum of epilepsies, ranging from benign isolated epilepsy to epileptic encephalopathy, with focal to bilateral tonic-clonic seizures being the most commonly observed. Importantly, a subgroup of patients present with a mild extra-neurological phenotype without characteristic facial dysmorphism or generalized hypertrichosis. This report expands the phenotypic spectrum of KNCK4-associated disease and provides new insights into the clinical heterogeneity of this rare neurodevelopmental syndrome.

Probing membrane deformation energy by KcsA potassium channel gating under varied membrane thickness and tension

This study investigated how membrane thickness and tension modify the gating of KcsA potassium channels when simultaneously varied. The KcsA channel undergoes global conformational changes upon gating: expansion of the cross-sectional area and longitudinal shortening upon opening. Thus, membranes impose differential effects on the open and closed conformations, such as hydrophobic mismatches. Here, the single-channel open probability was recorded in the contact bubble bilayer, by which variable thickness membranes under a defined tension were applied. A fully open channel in thin membranes turned to sporadic openings in thick membranes, where the channel responded moderately to tension increase. Quantitative gating analysis prompted the hypothesis that tension augmented the membrane deformation energy when hydrophobic mismatch was enhanced in thick membranes.

Structure of the human K2P13.1(THIK-1) channel reveals a novel hydrophilic pore restriction and lipid cofactor site

The halothane-inhibited K2P leak potassium channel K2P13.1 (THIK-1)1-3 is found in diverse cells1,4 including neurons1,5 and microglia6-8 where it affects surveillance6, synaptic pruning7, phagocytosis7, and inflammasome-mediated interleukin-1β release6,8,9. As with many K2Ps1,5,10-14 and other voltage-gated ion channel (VGIC) superfamily members3,15,16, polyunsaturated fatty acid (PUFA) lipids modulate K2P13.1 (THIK-1)1,5,14,17 via a poorly understood mechanism. Here, we present cryo-electronmicroscopy (cryo-EM) structures of human K2P13.1 (THIK-1) and mutants in lipid nanodiscs and detergent. These reveal that, unlike other K2Ps13,18-24, K2P13.1 (THIK-1) has a two-chamber aqueous inner cavity obstructed by a M4 transmembrane helix tyrosine (Tyr273, the flow restrictor). This hydrophilic barrier can be opened by an activatory mutation, S136P25, at natural break in the M2 transmembrane helix and by intrinsic channel dynamics. The structures also reveal a buried lipid in the P1/M4 intersubunit interface at a location, the PUFA site, that coincides with the TREK subfamily K2P modulator pocket for small molecule agonists18,26,27. This overlap, together with the effects of mutation on K2P13.1 (THIK-1) PUFA responses, indicates that the PUFA site lipids are K2P13.1 (THIK-1) cofactors. Comparison with the PUFA-responsive VGIC Kv7.1 (KCNQ1)28-31 reveals a shared role for the equivalent pore domain intersubunit interface in lipid modulation, providing a framework for dissecting the effects of PUFAs on the VGIC superfamily. Our findings reveal the unique architecture underlying K2P13.1 (THIK-1) function, highlight the importance of the P1/M4 interface in control of K2Ps by both natural and synthetic agents, and should aid development of THIK subfamily modulators for diseases such as neuroinflammation6,32 and autism6.

C-type inactivation and proton modulation mechanisms of the TASK3 channel

The TWIK-related acid-sensitive K+ channel 3 (TASK3) belongs to the two-pore domain (K2P) potassium channel family, which regulates cell excitability by mediating a constitutive "leak" potassium efflux in the nervous system. Extracellular acidification inhibits TASK3 channel, but the molecular mechanism by which channel inactivation is coupled to pH decrease remains unclear. Here, we report the cryo-electron microscopy structures of human TASK3 at neutral and acidic pH. Structural comparison revealed selectivity filter (SF) rearrangements upon acidification, characteristic of C-type inactivation, but with a unique structural basis. The extracellular mouth of the SF was prominently dilated and simultaneously blocked by a hydrophobic gate. His98 protonation shifted the conformational equilibrium between the conductive and C-type inactivated SF toward the latter by engaging a cation-π interaction with Trp78, consistent with molecular dynamics simulations and electrophysiological experiments. Our work illustrated how TASK3 is gated in response to extracellular pH change and implies how physiological stimuli might directly modulate the C-type gating of K2P channels.

Tension activation of mechanosensitive two-pore domain K+ channels TRAAK, TREK-1, and TREK-2

TRAAK, TREK-1, and TREK-2 are mechanosensitive two-pore domain K+ (K2P) channels that contribute to action potential propagation, sensory transduction, and muscle contraction. While structural and functional studies have led to models that explain their mechanosensitivity, we lack a quantitative understanding of channel activation by membrane tension. Here, we define the tension response of mechanosensitive K2Ps using patch-clamp recording and imaging. All are low-threshold mechanosensitive channels (T10%/50% 0.6-2.7 / 4.4-6.4 mN/m) with distinct response profiles. TRAAK is most sensitive, TREK-1 intermediate, and TREK-2 least sensitive. TRAAK and TREK-1 are activated broadly over a range encompassing nearly all physiologically relevant tensions. TREK-2, in contrast, activates over a narrower range like mechanosensitive channels Piezo1, MscS, and MscL. We further show that low-frequency, low-intensity focused ultrasound increases membrane tension to activate TRAAK and MscS. This work provides insight into tension gating of mechanosensitive K2Ps relevant to understanding their physiological roles and potential applications for ultrasonic neuromodulation.

High-fidelity encoding of mechanostimuli by tactile food-sensing neurons requires an ensemble of ion channels

The nematode C. elegans uses mechanosensitive neurons to detect bacteria, which are food for worms. These neurons release dopamine to suppress foraging and promote dwelling. Through a screen of genes highly expressed in dopaminergic food-sensing neurons, we identify a K2P-family potassium channel-TWK-2-that damps their activity. Strikingly, loss of TWK-2 restores mechanosensation to neurons lacking the NOMPC-like channel transient receptor potential 4 (TRP-4), which was thought to be the primary mechanoreceptor for tactile food sensing. The alternate mechanoreceptor mechanism uncovered by TWK-2 mutation requires three Deg/ENaC channel subunits: ASIC-1, DEL-3, and UNC-8. Analysis of cell-physiological responses to mechanostimuli indicates that TRP and Deg/ENaC channels work together to set the range of analog encoding of stimulus intensity and to improve signal-to-noise characteristics and temporal fidelity of food-sensing neurons. We conclude that a specialized mechanosensory modality-tactile food sensing-emerges from coordination of distinct force-sensing mechanisms housed in one type of sensory neuron.

Membrane phospholipids control gating of the mechanosensitive potassium leak channel TREK1

Tandem pore domain (K2P) potassium channels modulate resting membrane potentials and shape cellular excitability. For the mechanosensitive subfamily of K2Ps, the composition of phospholipids within the bilayer strongly influences channel activity. To examine the molecular details of K2P lipid modulation, we solved cryo-EM structures of the TREK1 K2P channel bound to either the anionic lipid phosphatidic acid (PA) or the zwitterionic lipid phosphatidylethanolamine (PE). At the extracellular face of TREK1, a PA lipid inserts its hydrocarbon tail into a pocket behind the selectivity filter, causing a structural rearrangement that recapitulates mutations and pharmacology known to activate TREK1. At the cytoplasmic face, PA and PE lipids compete to modulate the conformation of the TREK1 TM4 gating helix. Our findings demonstrate two distinct pathways by which anionic lipids enhance TREK1 activity and provide a framework for a model that integrates lipid gating with the effects of other mechanosensitive K2P modulators.

The energetics of rapid cellular mechanotransduction

Cells throughout the human body detect mechanical forces. While it is known that the rapid (millisecond) detection of mechanical forces is mediated by force-gated ion channels, a detailed quantitative understanding of cells as sensors of mechanical energy is still lacking. Here, we combine atomic force microscopy with patch-clamp electrophysiology to determine the physical limits of cells expressing the force-gated ion channels (FGICs) Piezo1, Piezo2, TREK1, and TRAAK. We find that, depending on the ion channel expressed, cells can function either as proportional or nonlinear transducers of mechanical energy and detect mechanical energies as little as ~100 fJ, with a resolution of up to ~1 fJ. These specific energetic values depend on cell size, channel density, and cytoskeletal architecture. We also make the surprising discovery that cells can transduce forces either nearly instantaneously (<1 ms) or with a substantial time delay (~10 ms). Using a chimeric experimental approach and simulations, we show how such delays can emerge from channel-intrinsic properties and the slow diffusion of tension in the membrane. Overall, our experiments reveal the capabilities and limits of cellular mechanosensing and provide insights into molecular mechanisms that different cell types may employ to specialize for their distinct physiological roles.

Pressure and ultrasound activate mechanosensitive TRAAK K + channels through increased membrane tension

TRAAK is a mechanosensitive two-pore domain K ⁺ (K2P) channel found in nodes of Ranvier within myelinated axons. It displays low leak activity at rest and is activated up to one hundred-fold by increased membrane tension. Structural and functional studies have led to physical models for channel gating and mechanosensitivity, but no quantitative analysis of channel activation by tension has been reported. Here, we use simultaneous patch-clamp recording and fluorescent imaging to determine the tension response characteristics of TRAAK. TRAAK shows high sensitivity and a broad response to tension spanning nearly the entire physiologically relevant tension range. This graded response profile distinguishes TRAAK from similarly low-threshold mechanosensitive channels Piezo1 and MscS, which activate in a step-like fashion over a narrow tension range. We further use patch imaging to show that ultrasonic activation of TRAAK and MscS is due to increased membrane tension. Together, these results provide mechanistic insight into TRAAK tension gating, a framework for exploring the role of mechanosensitive K ⁺ channels at nodes of Ranvier, and biophysical context for developing ultrasound as a mechanical stimulation technique for neuromodulation.

Structural Basis for pH-gating of the K+ channel TWIK1 at the selectivity filter

TWIK1 (K2P1.1, KCNK1) is a widely expressed pH-gated two-pore domain K+ channel (K2P) that contributes to cardiac rhythm generation and insulin release from pancreatic beta cells. TWIK1 displays unique properties among K2Ps including low basal activity and inhibition by extracellular protons through incompletely understood mechanisms. Here, we present cryo-EM structures of TWIK1 in lipid nanodiscs at high and low pH that reveal a previously undescribed gating mechanism at the K+ selectivity filter. At high pH, TWIK1 adopts an open conformation. At low pH, protonation of an extracellular histidine results in a cascade of conformational changes that close the channel by sealing the top of the selectivity filter, displacing the helical cap to block extracellular ion access pathways, and opening gaps for lipid block of the intracellular cavity. These data provide a mechanistic understanding for extracellular pH-gating of TWIK1 and illustrate how diverse mechanisms have evolved to gate the selectivity filter of K+ channels.

Transition between conformational states of the TREK-1 K2P channel promoted by interaction with PIP2

Members of the TREK family of two-pore domain (K2P) potassium channels are highly sensitive to regulation by membrane lipids, including phosphatidylinositol-4,5-bisphosphate (PIP₂). Previous studies have demonstrated that PIP₂ increases TREK1 channel activity, however, the mechanistic understanding of the conformational transitions induced by PIP₂ remain unclear. Here, we used coarse-grained molecular dynamics (CG-MD) and atomistic MD simulations to model the PIP₂ binding site on both the up and down state conformations of TREK-1. We also calculated the free energy of PIP₂ binding relative to other anionic phospholipids in both conformational states using potential of mean force (PMF) and free energy perturbation (FEP) calculations. Our results identify state-dependent binding of PIP₂ to sites involving the proximal C-terminus and we show that PIP₂ promotes a conformational transition from a down state towards an intermediate that resembles the up state. These results are consistent with functional data for PIP₂ regulation and together provide evidence for a structural mechanism of TREK-1 channel activation by phosphoinositides.