Interaction between residues in the Mg2+-binding site regulates BK channel activation

BK Channelopathies and KCNMA1-Linked Disease Models

Novel KCNMA1 variants, encoding the BK K+ channel, are associated with a debilitating dyskinesia and epilepsy syndrome. Neurodevelopmental delay, cognitive disability, and brain and structural malformations are also diagnosed at lower incidence. More than half of affected individuals present with a rare negative episodic motor disorder, paroxysmal nonkinesigenic dyskinesia (PNKD3). The mechanistic relationship of PNKD3 to epilepsy and the broader spectrum of KCNMA1-associated symptomology is unknown. This review summarizes patient-associated KCNMA1 variants within the BK channel structure, functional classifications, genotype-phenotype associations, disease models, and treatment. Patient and transgenic animal data suggest delineation of gain-of-function (GOF) and loss-of-function KCNMA1 neurogenetic disease, validating two heterozygous alleles encoding GOF BK channels (D434G and N999S) as causing seizure and PNKD3. This discovery led to a variant-defined therapeutic approach for PNKD3, providing initial insight into the neurological basis. A comprehensive clinical definition of monogenic KCNMA1-linked disease and the neuronal mechanisms currently remain priorities for continued investigation.

High-resolution structures illuminate key principles underlying voltage and LRRC26 regulation of Slo1 channels

Multi-modal regulation of Slo1 channels by membrane voltage, intracellular calcium, and auxiliary subunits enables its pleiotropic physiological functions. Our understanding of how voltage impacts Slo1 conformational dynamics and the mechanisms by which auxiliary subunits, particularly of the LRRC (Leucine Rich Repeat containing) family of proteins, modulate its voltage gating remain unresolved. Here, we used single particle cryo-electron microscopy to determine structures of human Slo1 mutants which functionally stabilize the closed pore (F315A) or the activated voltage-sensor (R207A). Our structures, obtained under calcium-free conditions, reveal that a key step in voltage-sensing by Slo1 involves a rotameric flip of the voltage-sensing charges (R210 and R213) moving them by ∼6 Å across a hydrophobic gasket. Next we obtained reconstructions of a complex of human Slo1 with the human LRRC26 (γ1) subunit in absence of calcium. Together with extensive biochemical tests, we show that the extracellular domains of γ1 form a ring of interlocked dominos that stabilizes the quaternary assembly of the complex and biases Slo1:γ1 assembly towards high stoichiometric complexes. The transmembrane helix of γ1 is kinked and tightly packed against the Slo1 voltage-sensor. We hypothesize that γ1 subunits exert relatively small effects on early steps in voltage-gating but structurally stabilize non-S4 helices of Slo1 voltage-sensor which energetically facilitate conformational rearrangements that occur late in voltage stimulated transitions.

Structural mapping of patient-associated KCNMA1 gene variants

KCNMA1-linked channelopathy is a neurological disorder characterized by seizures, motor abnormalities, and neurodevelopmental disabilities. The disease mechanisms are predicted to result from alterations in KCNMA1-encoded BK K+ channel activity; however, only a subset of the patient-associated variants have been functionally studied. The localization of these variants within the tertiary structure or evaluation by pathogenicity algorithms has not been systematically assessed. In this study, 82 nonsynonymous patient-associated KCNMA1 variants were mapped within the BK channel protein. Fifty-three variants localized within cryoelectron microscopy-resolved structures, including 21 classified as either gain of function (GOF) or loss of function (LOF) in BK channel activity. Clusters of LOF variants were identified in the pore, the AC region (RCK1), and near the Ca2+ bowl (RCK2), overlapping with sites of pharmacological or endogenous modulation. However, no clustering was found for GOF variants. To further understand variants of uncertain significance (VUSs), assessments by multiple standard pathogenicity algorithms were compared, and new thresholds for sensitivity and specificity were established from confirmed GOF and LOF variants. An ensemble algorithm was constructed (KCNMA1 meta score (KMS)), consisting of a weighted summation of this trained dataset combined with a structural component derived from the Ca2+-bound and unbound BK channels. KMS assessment differed from the highest-performing individual algorithm (REVEL) at 10 VUS residues, and a subset were studied further by electrophysiology in HEK293 cells. M578T, E656A, and D965V (KMS+;REVEL-) were confirmed to alter BK channel properties in voltage-clamp recordings, and D800Y (KMS-;REVEL+) was assessed as benign under the test conditions. However, KMS failed to accurately assess K457E. These combined results reveal the distribution of potentially disease-causing KCNMA1 variants within BK channel functional domains and pathogenicity evaluation for VUSs, suggesting strategies for improving channel-level predictions in future studies by building on ensemble algorithms such as KMS.

Structural mapping of patient-associated KCNMA1 gene variants

KCNMA1-linked channelopathy is a neurological disorder characterized by seizures, motor abnormalities, and neurodevelopmental disabilities. The disease mechanisms are predicted to result from alterations in KCNMA1-encoded BK K+ channel activity; however, only a subset of the patient-associated variants have been functionally studied. The localization of these variants within the tertiary structure or evaluation by pathogenicity algorithms has not been systematically assessed. In this study, 82 nonsynonymous patient-associated KCNMA1 variants were mapped within the BK channel protein. Fifty-three variants localized within cryo-EM resolved structures, including 21 classified as either gain-of-function (GOF) or loss-of-function (LOF) in BK channel activity. Clusters of LOF variants were identified in the pore, the AC region (RCK1), and near the Ca 2+ bowl (RCK2), overlapping with sites of pharmacological or endogenous modulation. However, no clustering was found for GOF variants. To further understand variants of uncertain significance (VUS), assessments by multiple standard pathogenicity algorithms were compared, and new thresholds for sensitivity and specificity were established from confirmed GOF and LOF variants. An ensemble algorithm was constructed (KCNMA1 Meta Score), consisting of a weighted summation of this trained dataset combined with a structural component derived from the Ca 2+ bound and unbound BK channels. KMS assessment differed from the highest performing individual algorithm (REVEL) at 10 VUS residues, and a subset were studied further by electrophysiology in HEK293 cells. M578T, E656A, and D965V (KMS+;REVEL-) were confirmed to alter BK channel properties in voltage-clamp recordings, and D800Y (KMS-;REVEL+) was assessed as benign under the test conditions. However, KMS failed to accurately assess K457E. These combined results reveal the distribution of potentially disease-causing KCNMA1 variants within BK channel functional domains and pathogenicity evaluation for VUS, suggesting strategies for improving channel-level predictions in future studies by building on ensemble algorithms such as KMS.

Incorporating physics to overcome data scarcity in predictive modeling of protein function: A case study of BK channels

Machine learning has played transformative roles in numerous chemical and biophysical problems such as protein folding where large amount of data exists. Nonetheless, many important problems remain challenging for data-driven machine learning approaches due to the limitation of data scarcity. One approach to overcome data scarcity is to incorporate physical principles such as through molecular modeling and simulation. Here, we focus on the big potassium (BK) channels that play important roles in cardiovascular and neural systems. Many mutants of BK channel are associated with various neurological and cardiovascular diseases, but the molecular effects are unknown. The voltage gating properties of BK channels have been characterized for 473 site-specific mutations experimentally over the last three decades; yet, these functional data by themselves remain far too sparse to derive a predictive model of BK channel voltage gating. Using physics-based modeling, we quantify the energetic effects of all single mutations on both open and closed states of the channel. Together with dynamic properties derived from atomistic simulations, these physical descriptors allow the training of random forest models that could reproduce unseen experimentally measured shifts in gating voltage, ∆V1/2, with a RMSE ~ 32 mV and correlation coefficient of R ~ 0.7. Importantly, the model appears capable of uncovering nontrivial physical principles underlying the gating of the channel, including a central role of hydrophobic gating. The model was further evaluated using four novel mutations of L235 and V236 on the S5 helix, mutations of which are predicted to have opposing effects on V1/2 and suggest a key role of S5 in mediating voltage sensor-pore coupling. The measured ∆V1/2 agree quantitatively with prediction for all four mutations, with a high correlation of R = 0.92 and RMSE = 18 mV. Therefore, the model can capture nontrivial voltage gating properties in regions where few mutations are known. The success of predictive modeling of BK voltage gating demonstrates the potential of combining physics and statistical learning for overcoming data scarcity in nontrivial protein function prediction.

Incorporating physics to overcome data scarcity in predictive modeling of protein function: a case study of BK channels

Machine learning has played transformative roles in numerous chemical and biophysical problems such as protein folding where large amount of data exists. Nonetheless, many important problems remain challenging for data-driven machine learning approaches due to the limitation of data scarcity. One approach to overcome data scarcity is to incorporate physical principles such as through molecular modeling and simulation. Here, we focus on the big potassium (BK) channels that play important roles in cardiovascular and neural systems. Many mutants of BK channel are associated with various neurological and cardiovascular diseases, but the molecular effects are unknown. The voltage gating properties of BK channels have been characterized for 473 site-specific mutations experimentally over the last three decades; yet, these functional data by themselves remain far too sparse to derive a predictive model of BK channel voltage gating. Using physics-based modeling, we quantify the energetic effects of all single mutations on both open and closed states of the channel. Together with dynamic properties derived from atomistic simulations, these physical descriptors allow the training of random forest models that could reproduce unseen experimentally measured shifts in gating voltage, ΔV 1/2 , with a RMSE ∼ 32 mV and correlation coefficient of R ∼ 0.7. Importantly, the model appears capable of uncovering nontrivial physical principles underlying the gating of the channel, including a central role of hydrophobic gating. The model was further evaluated using four novel mutations of L235 and V236 on the S5 helix, mutations of which are predicted to have opposing effects on V 1/2 and suggest a key role of S5 in mediating voltage sensor-pore coupling. The measured ΔV 1/2 agree quantitatively with prediction for all four mutations, with a high correlation of R = 0.92 and RMSE = 18 mV. Therefore, the model can capture nontrivial voltage gating properties in regions where few mutations are known. The success of predictive modeling of BK voltage gating demonstrates the potential of combining physics and statistical learning for overcoming data scarcity in nontrivial protein function prediction.

Author Summary

Deep machine learning has brought many exciting breakthroughs in chemistry, physics and biology. These models require large amount of training data and struggle when the data is scarce. The latter is true for predictive modeling of the function of complex proteins such as ion channels, where only hundreds of mutational data may be available. Using the big potassium (BK) channel as a biologically important model system, we demonstrate that a reliable predictive model of its voltage gating property could be derived from only 473 mutational data by incorporating physics-derived features, which include dynamic properties from molecular dynamics simulations and energetic quantities from Rosetta mutation calculations. We show that the final random forest model captures key trends and hotspots in mutational effects of BK voltage gating, such as the important role of pore hydrophobicity. A particularly curious prediction is that mutations of two adjacent residues on the S5 helix would always have opposite effects on the gating voltage, which was confirmed by experimental characterization of four novel mutations. The current work demonstrates the importance and effectiveness of incorporating physics in predictive modeling of protein function with scarce data.

An allosteric modulator activates BK channels by perturbing coupling between Ca2+ binding and pore opening

BK type Ca²⁺-activated K⁺ channels activate in response to both voltage and Ca²⁺. The membrane-spanning voltage sensor domain (VSD) activation and Ca²⁺ binding to the cytosolic tail domain (CTD) open the pore across the membrane, but the mechanisms that couple VSD activation and Ca²⁺ binding to pore opening  are not clear. Here we show that a compound, BC5, identified from in silico screening, interacts with the CTD-VSD interface and specifically modulates the Ca²⁺ dependent activation mechanism. BC5 activates the channel in the absence of Ca²⁺ binding but Ca²⁺ binding inhibits BC5 effects. Thus, BC5 perturbs a pathway that couples Ca²⁺ binding to pore opening to allosterically affect both, which is further supported by atomistic simulations and mutagenesis. The results suggest that the CTD-VSD interaction makes a major contribution to the mechanism of Ca²⁺ dependent activation and is an important site for allosteric agonists to modulate BK channel activation.

BK Channel Gating Mechanisms: Progresses Toward a Better Understanding of Variants Linked Neurological Diseases

The large conductance Ca2+-activated potassium (BK) channel is activated by both membrane potential depolarization and intracellular Ca2+ with distinct mechanisms. Neural physiology is sensitive to the function of BK channels, which is shown by the discoveries of neurological disorders that are associated with BK channel mutations. This article reviews the molecular mechanisms of BK channel activation in response to voltage and Ca2+ binding, including the recent progress since the publication of the atomistic structure of the whole BK channel protein, and the neurological disorders associated with BK channel mutations. These results demonstrate the unique mechanisms of BK channel activation and that these mechanisms are important factors in linking BK channel mutations to neurological disorders.

Coupling of Ca2+ and voltage activation in BK channels through the αB helix/voltage sensor interface

Large-conductance Ca²⁺ and voltage-activated K⁺ (BK) channels control membrane excitability in many cell types. BK channels are tetrameric. Each subunit is composed of a voltage sensor domain (VSD), a central pore-gate domain, and a large cytoplasmic domain (CTD) that contains the Ca²⁺ sensors. While it is known that BK channels are activated by voltage and Ca²⁺, and that voltage and Ca²⁺ activations interact, less is known about the mechanisms involved. We explore here these mechanisms by examining the gating contribution of an interface formed between the VSDs and the αB helices located at the top of the CTDs. Proline mutations in the αB helix greatly decreased voltage activation while having negligible effects on gating currents. Analysis with the Horrigan, Cui, and Aldrich model indicated a decreased coupling between voltage sensors and pore gate. Proline mutations decreased Ca²⁺ activation for both Ca²⁺ bowl and RCK1 Ca²⁺ sites, suggesting that both high-affinity Ca²⁺ sites transduce their effect, at least in part, through the αB helix. Mg²⁺ activation also decreased. The crystal structure of the CTD with proline mutation L390P showed a flattening of the first helical turn in the αB helix compared to wild type, without other notable differences in the CTD, indicating that structural changes from the mutation were confined to the αB helix. These findings indicate that an intact αB helix/VSD interface is required for effective coupling of Ca²⁺ binding and voltage depolarization to pore opening and that shared Ca²⁺ and voltage transduction pathways involving the αB helix may be involved.

Threading the biophysics of mammalian Slo1 channels onto structures of an invertebrate Slo1 channel

Zhou et al. consider the biophysics of large-conductance Ca²⁺-activated Slo1 channels in the context of Aplysia Slo1 structures.

For those interested in the machinery of ion channel gating, the Ca²⁺ and voltage-activated BK K⁺ channel provides a compelling topic for investigation, by virtue of its dual allosteric regulation by both voltage and intracellular Ca²⁺ and because its large-single channel conductance facilitates detailed kinetic analysis. Over the years, biophysical analyses have illuminated details of the allosteric regulation of BK channels and revealed insights into the mechanism of BK gating, e.g., inner cavity size and accessibility and voltage sensor-pore coupling. Now the publication of two structures of an Aplysia californica BK channel—one liganded and one metal free—promises to reinvigorate functional studies and interpretation of biophysical results. The new structures confirm some of the previous functional inferences but also suggest new perspectives regarding cooperativity between Ca²⁺-binding sites and the relationship between voltage- and Ca²⁺-dependent gating. Here we consider the extent to which the two structures explain previous functional data on pore-domain properties, voltage-sensor motions, and divalent cation binding and activation of the channel.

Understanding the conformational motions of RCK gating rings

A timely review of the structural basis of Ca²⁺-activated K⁺ channel modulation by regulator of conduction of K⁺ (RCK) domains

Regulator of conduction of K⁺ (RCK) domains are ubiquitous regulators of channel and transporter activity in prokaryotes and eukaryotes. In humans, RCK domains form an integral component of large-conductance calcium-activated K channels (BK channels), key modulators of nerve, muscle, and endocrine cell function. In this review, we explore how the study of RCK domains in bacterial and human channels has contributed to our understanding of the structural basis of channel function. This knowledge will be critical in identifying mechanisms that underlie BK channelopathies that lead to epilepsy and other diseases, as well as regions of the channel that might be successfully targeted to treat such diseases.

Deletion of cytosolic gating ring decreases gate and voltage sensor coupling in BK channels

Both cellular depolarization and intracellular Ca²⁺ can gate open large conductance Ca²⁺-activated K⁺ channels. Zhang et al. show that the intracellular gating ring, which forms the Ca²⁺-sensing machinery of the channel, is also required for activated voltage sensors to effectively gate open the pore.

Large conductance Ca²⁺-activated K⁺ channels (BK channels) gate open in response to both membrane voltage and intracellular Ca²⁺. The channel is formed by a central pore-gate domain (PGD), which spans the membrane, plus transmembrane voltage sensors and a cytoplasmic gating ring that acts as a Ca²⁺ sensor. How these voltage and Ca²⁺ sensors influence the common activation gate, and interact with each other, is unclear. A previous study showed that a BK channel core lacking the entire cytoplasmic gating ring (Core-MT) was devoid of Ca²⁺ activation but retained voltage sensitivity (Budelli et al. 2013. Proc. Natl. Acad. Sci. USA. http://dx.doi.org/10.1073/pnas.1313433110). In this study, we measure voltage sensor activation and pore opening in this Core-MT channel over a wide range of voltages. We record gating currents and find that voltage sensor activation in this truncated channel is similar to WT but that the coupling between voltage sensor activation and gating of the pore is reduced. These results suggest that the gating ring, in addition to being the Ca²⁺ sensor, enhances the effective coupling between voltage sensors and the PGD. We also find that removal of the gating ring alters modulation of the channels by the BK channel’s β1 and β2 subunits.

Molecular Determinants of BK Channel Functional Diversity and Functioning

Large-conductance Ca2+- and voltage-activated K+ (BK) channels play many physiological roles ranging from the maintenance of smooth muscle tone to hearing and neurosecretion. BK channels are tetramers in which the pore-forming α subunit is coded by a single gene ( Slowpoke, KCNMA1). In this review, we first highlight the physiological importance of this ubiquitous channel, emphasizing the role that BK channels play in different channelopathies. We next discuss the modular nature of BK channel-forming protein, in which the different modules (the voltage sensor and the Ca2+ binding sites) communicate with the pore gates allosterically. In this regard, we review in detail the allosteric models proposed to explain channel activation and how the models are related to channel structure. Considering their extremely large conductance and unique selectivity to K+, we also offer an account of how these two apparently paradoxical characteristics can be understood consistently in unison, and what we have learned about the conduction system and the activation gates using ions, blockers, and toxins. Attention is paid here to the molecular nature of the voltage sensor and the Ca2+ binding sites that are located in a gating ring of known crystal structure and constituted by four COOH termini. Despite the fact that BK channels are coded by a single gene, diversity is obtained by means of alternative splicing and modulatory β and γ subunits. We finish this review by describing how the association of the α subunit with β or with γ subunits can change the BK channel phenotype and pharmacology.

Structural basis for gating the high-conductance Ca2+-activated K+ channel

The precise control of an ion channel gate by environmental stimuli is crucial for the fulfilment of its biological role. The gate in Slo1 K⁺ channels is regulated by two separate stimuli, intracellular Ca²⁺ concentration and membrane voltage. Slo1 is thus central to understanding the relationship between intracellular Ca²⁺ and membrane excitability. Here we present the Slo1 structure from Aplysia californica in the absence of Ca²⁺ and compare it with the Ca²⁺-bound channel. We show that Ca²⁺ binding at two unique binding sites per subunit stabilizes an expanded conformation of the Ca²⁺ sensor gating ring. These conformational changes are propagated from the gating ring to the pore through covalent linkers and through protein interfaces formed between the gating ring and the voltage sensors. The gating ring and the voltage sensors are directly connected through these interfaces, which allow membrane voltage to regulate gating of the pore by influencing the Ca²⁺ sensors.

BK channels: multiple sensors, one activation gate

Ion transport across cell membranes is essential to cell communication and signaling. Passive ion transport is mediated by ion channels, membrane proteins that create ion conducting pores across cell membrane to allow ion flux down electrochemical gradient. Under physiological conditions, majority of ion channel pores are not constitutively open. Instead, structural region(s) within these pores breaks the continuity of the aqueous ion pathway, thereby serves as activation gate(s) to control ions flow in and out. To achieve spatially and temporally regulated ion flux in cells, many ion channels have evolved sensors to detect various environmental stimuli or the metabolic states of the cell and trigger global conformational changes, thereby dynamically operate the opening and closing of their activation gate. The sensors of ion channels can be broadly categorized as chemical sensors and physical sensors to respond to chemical (such as neural transmitters, nucleotides and ions) and physical (such as voltage, mechanical force and temperature) signals, respectively. With the rapidly growing structural and functional information of different types of ion channels, it is now critical to understand how ion channel sensors dynamically control their gates at molecular and atomic level. The voltage and Ca²⁺ activated BK channels, a K⁺ channel with an electrical sensor and multiple chemical sensors, provide a unique model system for us to understand how physical and chemical energy synergistically operate its activation gate.

A point mutation in the human Slo1 channel that impairs its sensitivity to omega-3 docosahexaenoic acid

Long-chain polyunsaturated omega-3 fatty acids such as docosahexaenoic acid (DHA) at nanomolar concentrations reversibly activate human large-conductance Ca²⁺- and voltage-gated K⁺ (Slo1 BK) channels containing auxiliary β1 or β4 subunits in cell-free patches. Here we examined the action of DHA on the Slo1 channel without any auxiliary subunit and sought to elucidate the biophysical mechanism and the molecular determinants of the DHA sensitivity. Measurements of ionic currents through human Slo1 (hSlo1) channels reveal that the stimulatory effect of DHA does not require activation of the voltage or Ca²⁺ sensors. Unlike gating of the hSlo1 channel, that of the Drosophila melanogaster Slo1 (dSlo1) channel is unaltered by DHA. Our mutagenesis study based on the differential responses of human and dSlo1 channels to DHA pinpoints that Y318 near the cytoplasmic end of S6 in the hSlo1 channel is a critical determinant of the stimulatory action of DHA. The mutation Y318S in hSlo1, which replaces Y with S as found in dSlo1, greatly diminishes the channel’s response to DHA with a 22-carbon chain whether β1 or β4 is absent or present. However, the responses to α-linolenic acid, an omegea-3 fatty acid with an 18-carbon chain, and to arachidonic acid, an omega-6 fatty acid with a 20-carbon chain, remain unaffected by the mutation. Y318 in the S6 segment of hSlo1 is thus an important determinant of the electrophysiological response of the channel to DHA. Furthermore, the mutation Y318S may prove to be useful in dissecting out the complex lipid-mediated modulation of Slo1 BK channels.

The interface between membrane-spanning and cytosolic domains in Ca²+-dependent K+ channels is involved in β subunit modulation of gating

Large-conductance, voltage-, and Ca²⁺-dependent K⁺ (BK) channels are broadly expressed in various tissues to modulate neuronal activity, smooth muscle contraction, and secretion. BK channel activation depends on the interactions among the voltage sensing domain (VSD), the cytosolic domain (CTD), and the pore gate domain (PGD) of the Slo1 α-subunit, and is further regulated by accessory β subunits (β1-β4). How β subunits fine-tune BK channel activation is critical to understand the tissue-specific functions of BK channels. Multiple sites in both Slo1 and the β subunits have been identified to contribute to the interaction between Slo1 and the β subunits. However, it is unclear whether and how the interdomain interactions among the VSD, CTD, and PGD are altered by the β subunits to affect channel activation. Here we show that human β1 and β2 subunits alter interactions between bound Mg²⁺ and gating charge R213 and disrupt the disulfide bond formation at the VSD-CTD interface of mouse Slo1, indicating that the β subunits alter the VSD-CTD interface. Reciprocally, mutations in the Slo1 that alter the VSD-CTD interaction can specifically change the effects of the β1 subunit on the Ca²⁺ activation and of the β2 subunit on the voltage activation. Together, our data suggest a novel mechanism by which the β subunits modulated BK channel activation such that a β subunit may interact with the VSD or the CTD and alter the VSD-CTD interface of the Slo1, which enables the β subunit to have effects broadly on both voltage and Ca²⁺-dependent activation.