Intersubunit coupling in the pore of BK channels

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.

A gating lever and molecular logic gate that couple voltage and calcium sensor activation to opening in BK potassium channels

BK channels uniquely integrate voltage and calcium signaling in diverse cell types through allosteric activation of their K+-conducting pore by structurally distinct V and Ca2+ sensor domains. Here, we define mechanisms and interaction pathways that link V sensors to the pore by analyzing effects on allosteric coupling of point mutations in the context of Slo1 BK channel structure. A gating lever, mediated by S4/S5 segment interaction within the transmembrane domain, rotates to engage and stabilize the open conformation of the S6 inner pore helix upon V sensor activation. In addition, an indirect pathway, mediated by the carboxyl-terminal cytosolic domain (CTD) and C-linker that connects the CTD to S6, stabilizes the closed conformation when V sensors are at rest. Unexpectedly, this mechanism, which bypasses the covalent connections of C-linker to CTD and pore, also transduces Ca2+-dependent coupling in a manner that is completely nonadditive with voltage, analogous to the function of a digital logic (OR) gate.

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.

Atomic determinants of BK channel activation by polyunsaturated fatty acids

Docosahexaenoic acid (DHA), a polyunsaturated ω-3 fatty acid enriched in oily fish, contributes to better health by affecting multiple targets. Large-conductance Ca²⁺- and voltage-gated Slo1 BK channels are directly activated by nanomolar levels of DHA. We investigated DHA-channel interaction by manipulating both the fatty acid structure and the channel composition through the site-directed incorporation of unnatural amino acids. Electrophysiological measurements show that the para-group of a Tyr residue near the ion conduction pathway has a critical role. To robustly activate the channel, ionization must occur readily by a fatty acid for a good efficacy, and a long nonpolar acyl tail with a Z double bond present at the halfway position for a high affinity. The results suggest that DHA and the channel form an ion-dipole bond to promote opening and demonstrate the channel druggability. DHA, a marine-derived nutraceutical, represents a promising lead compound for rational drug design and discovery.

Single Channel Recordings Reveal Differential β2 Subunit Modulations Between Mammalian and Drosophila BKCa(β2) Channels

Large-conductance Ca²⁺- and voltage-activated potassium (BK) channels are widely expressed in tissues. As a voltage and calcium sensor, BK channels play significant roles in regulating the action potential frequency, neurotransmitter release, and smooth muscle contraction. After associating with the auxiliary β2 subunit, mammalian BK(β2) channels (mouse or human Slo1/β2) exhibit enhanced activation and complete inactivation. However, how the β2 subunit modulates the Drosophila Slo1 channel remains elusive. In this study, by comparing the different functional effects on heterogeneous BK(β2) channel, we found that Drosophila Slo1/β2 channel exhibits “paralyzed”-like and incomplete inactivation as well as slow activation. Further, we determined three different modulations between mammalian and Drosophila BK(β2) channels: 1) dSlo1/β2 doesn’t have complete inactivation. 2) β2(K33,R34,K35) delays the dSlo1/Δ3-β2 channel activation. 3) dSlo1/β2 channel has enhanced pre-inactivation than mSlo1/β2 channel. The results in our study provide insights into the different modulations of β2 subunit between mammalian and Drosophila Slo1/β2 channels and structural basis underlie the activation and pre-inactivation of other BK(β) complexes.

Modulation of BK Channels by Small Endogenous Molecules and Pharmaceutical Channel Openers

Voltage- and Ca(2+)-activated K(+) channels of big conductance (BK channels) are abundantly found in various organs and their relevance for smooth muscle tone and neuronal signaling is well documented. Dysfunction of BK channels is implicated in an array of human diseases involving many organs including the nervous, pulmonary, cardiovascular, renal, and urinary systems. In humans a single gene (KCNMA1) encodes the pore-forming α subunit (Slo1) of BK channels, but the channel properties are variable because of alternative splicing, tissue- and subcellular-specific auxiliary subunits (β, γ), posttranslational modifications, and a multitude of endogenous signaling molecules directly affecting the channel function. Initiatives to develop drugs capable of activating BK channels (channel openers) therefore need to consider the tissue-specific variability of BK channel structure and the potential interference with endogenously produced regulatory factors. The atomic structural basis of BK channel function is only beginning to be revealed. However, building on detailed knowledge of BK channel function, including its single-channel characteristics, voltage- and Ca(2+) dependence of channel gating, and modulation by diffusible messengers, a multi-tier allosteric model of BK channel gating (Horrigan and Aldrich (HA) model) has become a valuable tool in studying modulation of the channel. Using the conceptual framework of the HA model, we here review the functional impact of endogenous modulatory factors and select small synthetic compounds that regulate BK channel activity. Furthermore, we devise experimental approaches for studying BK channel-drug interactions with the aim to classify BK-modulating substances according to their molecular mode of action.

Biophysics of BK Channel Gating

BK channels are universal regulators of cell excitability, given their exceptional unitary conductance selective for K(+), joint activation mechanism by membrane depolarization and intracellular [Ca(2+)] elevation, and broad expression pattern. In this chapter, we discuss the structural basis and operational principles of their activation, or gating, by membrane potential and calcium. We also discuss how the two activation mechanisms interact to culminate in channel opening. As members of the voltage-gated potassium channel superfamily, BK channels are discussed in the context of archetypal family members, in terms of similarities that help us understand their function, but also seminal structural and biophysical differences that confer unique functional properties.

Hydrophobic interaction between contiguous residues in the S6 transmembrane segment acts as a stimuli integration node in the BK channel

Phenylalanine 380 and leucine 377 in the BK channel S6 transmembrane helix of contiguous subunits participate in a hydrophobic interaction in both the closed and open state; this interaction is important in the allosteric coupling between the Ca²⁺ and voltage sensors and pore domain.

Large-conductance Ca²⁺- and voltage-activated K⁺ channel (BK) open probability is enhanced by depolarization, increasing Ca²⁺ concentration, or both. These stimuli activate modular voltage and Ca²⁺ sensors that are allosterically coupled to channel gating. Here, we report a point mutation of a phenylalanine (F380A) in the S6 transmembrane helix that, in the absence of internal Ca²⁺, profoundly hinders channel opening while showing only minor effects on the voltage sensor active–resting equilibrium. Interpretation of these results using an allosteric model suggests that the F380A mutation greatly increases the free energy difference between open and closed states and uncouples Ca²⁺ binding from voltage sensor activation and voltage sensor activation from channel opening. However, the presence of a bulky and more hydrophobic amino acid in the F380 position (F380W) increases the intrinsic open–closed equilibrium, weakening the coupling between both sensors with the pore domain. Based on these functional experiments and molecular dynamics simulations, we propose that F380 interacts with another S6 hydrophobic residue (L377) in contiguous subunits. This pair forms a hydrophobic ring important in determining the open–closed equilibrium and, like an integration node, participates in the communication between sensors and between the sensors and pore. Moreover, because of its effects on open probabilities, the F380A mutant can be used for detailed voltage sensor experiments in the presence of permeant cations.

Mutation of I696 and W697 in the TRP box of vanilloid receptor subtype I modulates allosteric channel activation

Residues I696 and W697 are crucial to coupling between the TRPV1 ligand- and voltage-sensing domains and the channel pore.

The transient receptor potential vanilloid receptor subtype I (TRPV1) channel acts as a polymodal sensory receptor gated by chemical and physical stimuli. Like other TRP channels, TRPV1 contains in its C terminus a short, conserved domain called the TRP box, which is necessary for channel gating. Substitution of two TRP box residues—I696 and W697—with Ala markedly affects TRPV1’s response to all activating stimuli, which indicates that these two residues play a crucial role in channel gating. We systematically replaced I696 and W697 with 18 native l-amino acids (excluding cysteine) and evaluated the effect on voltage- and capsaicin-dependent gating. Mutation of I696 decreased channel activation by either voltage or capsaicin; furthermore, gating was only observed with substitution of hydrophobic amino acids. Substitution of W697 with any of the 18 amino acids abolished gating in response to depolarization alone, shifting the threshold to unreachable voltages, but not capsaicin-mediated gating. Moreover, vanilloid-activated responses of W697X mutants showed voltage-dependent gating along with a strong voltage-independent component. Analysis of the data using an allosteric model of activation indicates that mutation of I696 and W697 primarily affects the allosteric coupling constants of the ligand and voltage sensors to the channel pore. Together, our findings substantiate the notion that inter- and/or intrasubunit interactions at the level of the TRP box are critical for efficient coupling of stimulus sensing and gate opening. Perturbation of these interactions markedly reduces the efficacy and potency of the activating stimuli. Furthermore, our results identify these interactions as potential sites for pharmacological intervention.

BK channel opening involves side-chain reorientation of multiple deep-pore residues

Three deep-pore locations, L312, A313, and A316, were identified in a scanning mutagenesis study of the BK (Ca(2+)-activated, large-conductance K(+)) channel S6 pore, where single aspartate substitutions led to constitutively open mutant channels (L312D, A313D, and A316D). To understand the mechanisms of the constitutive openness of these mutant channels, we individually mutated these three sites into the other 18 amino acids. We found that charged or polar side-chain substitutions at each of the sites resulted in constitutively open mutant BK channels, with high open probability at negative voltages, as well as a loss of voltage and Ca(2+) dependence. Given the fact that multiple pore residues in BK displayed side-chain hydrophilicity-dependent constitutive openness, we propose that BK channel opening involves structural rearrangement of the deep-pore region, where multiple residues undergo conformational changes that may increase the exposure of their side chains to the polar environment of the pore.

Transduction of voltage and Ca2+ signals by Slo1 BK channels

Large-conductance Ca2+ -and voltage-gated K+ channels are activated by an increase in intracellular Ca2+ concentration and/or depolarization. The channel activation mechanism is well described by an allosteric model encompassing the gate, voltage sensors, and Ca2+ sensors, and the model is an excellent framework to understand the influences of auxiliary β and γ subunits and regulatory factors such as Mg2+. Recent advances permit elucidation of structural correlates of the biophysical mechanism.

TRPV1 channels are functionally coupled with BK(mSlo1) channels in rat dorsal root ganglion (DRG) neurons

The transient receptor potential vanilloid receptor 1 (TRPV1) channel is a nonselective cation channel activated by a variety of exogenous and endogenous physical and chemical stimuli, such as temperature (≥42 °C), capsaicin, a pungent compound in hot chili peppers, and allyl isothiocyanate. Large-conductance calcium- and voltage-activated potassium (BK) channels regulate the electric activities and neurotransmitter releases in excitable cells, responding to changes in membrane potentials and elevation of cytosolic calcium ions (Ca²⁺). However, it is unknown whether the TRPV1 channels are coupled with the BK channels. Using patch-clamp recording combined with an infrared laser device, we found that BK channels could be activated at 0 mV by a Ca²⁺ influx through TRPV1 channels not the intracellular calcium stores in submilliseconds. The local calcium concentration around BK is estimated over 10 μM. The crosstalk could be affected by 10 mM BAPTA, whereas 5 mM EGTA was ineffectual. Fluorescence and co-immunoprecipitation experiments also showed that BK and TRPV1 were able to form a TRPV1-BK complex. Furthermore, we demonstrated that the TRPV1-BK coupling also occurs in dosal root ganglion (DRG) cells, which plays a critical physiological role in regulating the “pain” signal transduction pathway in the peripheral nervous system.

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.

Perspectives on: conformational coupling in ion channels: conformational coupling in BK potassium channels

Large conductance calcium- and voltage-dependent BK potassium channels (aka BK(Ca), MaxiK, Slo1, KCa1.1, and KCNMA1) are expressed in a wide variety of tissues throughout the body and are activated by both intracellular Ca(2+) and membrane depolarization. Owing to these properties, BK channels participate in diverse physiological processes from electrical excitability in neurons and secretory cells, and regulation of smooth muscle tone to tuning of auditory hair cells (Vergara et al., 1998; Ghatta et al., 2006). The response to voltage and Ca(2+) allows BK channels to integrate electrical and calcium signaling, which is central to their physiological role. Understanding how BK and other multimodal channels are regulated by and integrate diverse stimuli is not only physiologically important but also relevant to the topic of conformational coupling. As a voltage- and ligand-dependent channel, BK channels contain both voltage-sensor and ligand-binding domains as well as a gate to regulate the flow of K(+) through the pore. Coupling of conformational changes in one domain to another provides the basis for transducing voltage and ligand binding into channel opening and, therefore, defines, together with the functional properties of the gate and sensors, the signal transduction properties of the channel. The goal of this perspective is to provide an overview on the role and molecular basis of conformational coupling between functional domains in BK channels and outline some of the questions that remain to be answered.

DiBAC₄(3) hits a "sweet spot" for the activation of arterial large-conductance Ca²⁺-activated potassium channels independently of the β₁-subunit

The voltage-sensitive dye bis-(1,3-dibutylbarbituric acid)trimethine oxonol [DiBAC₄(3)] has been reported as a novel large-conductance Ca²⁺-activated K⁺ (BK) channel activator with selectivity for its β₁- or β₄-subunits. In arterial smooth muscle, BK channels are formed by a pore-forming α-subunit and a smooth muscle-abundant regulatory β₁-subunit. This tissue specificity has driven extensive pharmacological research aimed at regulating arterial tone. Using animals with a disruption of the gene for the β₁-subunit, we explored the effects of DiBAC₄(3) in native channels from arterial smooth muscle. We tested the hypothesis that, in native BK channels, activation by DiBAC₄(3) relies mostly on its α-subunit. We studied BK channels from wild-type and transgenic β₁-knockout mice in excised patches. BK channels from brain arteries, with or without the β₁-subunit, were similarly activated by DiBAC₄(3). In addition, we found that saturating concentrations of DiBAC₄(3) (~30 μM) promote an unprecedented persistent activation of the channel that negatively shifts its voltage dependence by as much as -300 mV. This "sweet spot" for persistent activation is independent of Ca²⁺ and/or the β₁₋₄-subunits and is fully achieved when DiBAC₄(3) is applied to the intracellular side of the channel. Arterial BK channel response to DiBAC₄(3) varies across species and/or vascular beds. DiBAC₄(3) unique effects can reveal details of BK channel gating mechanisms and help in the rational design of BK channel activators.

Effects of ketamine and its metabolites on ion currents in differentiated hippocampal H19-7 neuronal cells and in HEK293T cells transfected with α-hslo subunit

Ketamine (KT), a dissociative anesthetic, is known to induce schizophrenia-like psychosis. The percentage of KT abuse has recently grown fast despite KT being a controlled drug. The mechanism of KT actions is related to the inhibition of NMDA receptors. Whether KT produces other effects on ion currents in hippocampal neurons remains unclear. In this study, we attempted to evaluate the possible effects of KT and other related compounds on ion currents in hippocampal neuron-derived H19-7 cells. This drug exerted an inhibitory effect on Ca(2+)-activated K(+) current (IK(Ca)) in these cells with an IC(50) value of 274 μM. Pimaric acid (30 μM) or abietic acid (30 μM), known to stimulate large-conductance Ca(2+)-activated K(+) channels, reversed KT-induced inhibition of I(K)(Ca). In HEK293T cells expressing a-humans low poke, KT-induced inhibition of I(K)(Ca) still existed. Dehydronorketamine (300 μM) had little or no effect on the IK(Ca) amplitude, while norketamine (300 μM) slightly but significantly suppressed it. In inside–out configuration, KT applied to the intracellular face of the membrane did not alter single channel conductance of large-conductance Ca(2+)-activated K(+) (BKCa) channels; however, it did significantly reduce the probability of channel openings. Addition of KT was effective in depressing the peak amplitude of voltage-gated Na(+) current. Moreover, the presence of KT was noted to enhance the amplitude of membrane electroporation-induced inward currents (IMEP) in differentiated H19-7 cells. KT-stimulated IMEP was reversed by further application of LaCl(3) (100 μM), but not by NMDA (30 μM). The modulations by this compound of ion channels may contribute to the underlying mechanisms through which KT and its metabolites influence the electrical behavior of hippocampal neurons if similar findings occur in vivo.

Mechanism of the modulation of BK potassium channel complexes with different auxiliary subunit compositions by the omega-3 fatty acid DHA

Large-conductance Ca(2+)- and voltage-activated K(+) (BK) channels are well known for their functional versatility, which is bestowed in part by their rich modulatory repertoire. We recently showed that long-chain omega-3 polyunsaturated fatty acids such as docosahexaenoic acid (DHA) found in oily fish lower blood pressure by activating vascular BK channels made of Slo1+β1 subunits. Here we examined the action of DHA on BK channels with different auxiliary subunit compositions. Neuronal Slo1+β4 channels were just as well activated by DHA as vascular Slo1+β1 channels. In contrast, the stimulatory effect of DHA was much smaller in Slo1+β2, Slo1+LRRC26 (γ1), and Slo1 channels without auxiliary subunits. Mutagenesis of β1, β2, and β4 showed that the large effect of DHA in Slo1+β1 and Slo1+β4 is conferred by the presence of two residues, one in the N terminus and the other in the first transmembrane segment of the β1 and β4 subunits. Transfer of this amino acid pair from β1 or β4 to β2 introduces a large response to DHA in Slo1+β2. The presence of a pair of oppositely charged residues at the aforementioned positions in β subunits is associated with a large response to DHA. The Slo1 auxiliary subunits are expressed in a highly tissue-dependent fashion. Thus, the subunit composition-dependent stimulation by DHA demonstrates that BK channels are effectors of omega-3 fatty acids with marked tissue specificity.

Targeting BK (big potassium) channels in epilepsy

INTRODUCTION

Epilepsies are disorders of neuronal excitability characterized by spontaneous and recurrent seizures. Ion channels are critical for regulating neuronal excitability and, therefore, can contribute significantly to epilepsy pathophysiology. In particular, large conductance, Ca2+-activated K+ (BKCa) channels play an important role in seizure etiology. These channels are activated by both membrane depolarization and increased intracellular Ca2+. This unique coupling of Ca2+ signaling to membrane depolarization is important in controlling neuronal hyperexcitability, as outward K+ current through BKCa channels hyperpolarizes neurons.

AREAS COVERED

BKCa channel structure-function and the role of these channels in epilepsy pathophysiology.

EXPERT OPINION

Loss-of-function BKCa channel mutations contribute to neuronal hyperexcitability that can lead to temporal lobe epilepsy, tonic-clonic seizures and alcohol withdrawal seizures. Similarly, BKCa channel blockade can trigger seizures and status epilepticus. Paradoxically, some mutations in BKCa channel subunit can give rise to channel gain-of-function that leads to development of idiopathic epilepsy (primarily absence epilepsy). Seizures themselves also enhance BKCa channel currents associated with neuronal hyperexcitability, and blocking BKCa channels suppresses generalized tonic-clonic seizures. Thus, both loss-of-function and gain-of-function BKCa channels might serve as molecular targets for drugs to suppress certain seizure phenotypes including temporal lobe seizures and absence seizures, respectively.

Cysteine scanning and modification reveal major differences between BK channels and Kv channels in the inner pore region

BK channels are regulated by two distinct physiological signals, transmembrane potential and intracellular Ca(2+), each acting through independent modular sensor domains. However, despite a presumably central role in the coupling of sensor activation to channel gating, the pore-lining S6 transmembrane segment has not been systematically studied. Here, cysteine substitution and modification studies of the BK S6 point to substantial differences between BK and Kv channels in the structure and function of the S6-lined inner pore. Gating shifts caused by introduction of cysteines define a pattern and direction of free energy changes in BK S6 distinct from Shaker. Modification of BK S6 residues identifies pore-facing residues that occur at different linear positions along aligned BK and Kv S6 segments. Periodicity analysis suggests that one factor contributing to these differences may be a disruption of the BK S6 α-helix from the unique diglycine motif at the position of the Kv hinge glycine. State-dependent MTS accessibility reveals that, even in closed states, modification can occur. Furthermore, the inner pore of BK channels is much larger than that of K(+) channels with solved crystal structures. The results suggest caution in the use of Kv channel structures as templates for BK homology models, at least in the pore-gate domain.

Coupling and cooperativity in voltage activation of a limited-state BK channel gating in saturating Ca2+

Voltage-dependent gating mechanisms of large conductance Ca²⁺ and voltage-activated (BK) channels were investigated using two-dimensional maximum likelihood analysis of single-channel open and closed intervals. To obtain sufficient data at negative as well as positive voltages, single-channel currents were recorded at saturating Ca²⁺ from BK channels mutated to remove the RCK1 Ca²⁺ and Mg²⁺ sensors. The saturating Ca²⁺ acting on the Ca²⁺ bowl sensors of the resulting BK_B channels increased channel activity while driving the gating into a reduced number of states, simplifying the model. Five highly constrained idealized gating mechanisms based on extensions of the Monod-Wyman-Changeux model for allosteric proteins were examined. A 10-state model without coupling between the voltage sensors and the opening/closing transitions partially described the voltage dependence of Po but not the single-channel kinetics. With allowed coupling, the model gave improved descriptions of Po and approximated the single-channel kinetics; each activated voltage sensor increased the opening rate approximately an additional 23-fold while having little effect on the closing rate. Allowing cooperativity among voltage sensors further improved the description of the data: each activated voltage sensor increased the activation rate of the remaining voltage sensors approximately fourfold, with little effect on the deactivation rate. The coupling factor was decreased in models with cooperativity from ∼23 to ∼18. Whether the apparent cooperativity among voltage sensors arises from imposing highly idealized models or from actual cooperativity will require additional studies to resolve. For both cooperative and noncooperative models, allowing transitions to five additional brief (flicker) closed states further improved the description of the data. These observations show that the voltage-dependent single-channel kinetics of BK_B channels can be approximated by highly idealized allosteric models in which voltage sensor movement increases Po mainly through an increase in channel opening rates, with limited effects on closing rates.