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

Electronic Polarization Leads to a Drier Dewetted State for Hydrophobic Gating in the Big Potassium Channel

In the hydrophobic gating mechanism proposed for some ion channels, ion permeation is not blocked by the physical dimension of the channel pore but by its dewetted state which constitutes the energetic bottleneck. A major source of uncertainty in the mechanism is that the dewetted state was not observed in experiments and only probed in simulations using nonpolarizable force fields, which do not accurately represent the properties of confined water. Here we analyze hydration of the central cavity in the pore-gate domain of the Big Potassium channel using molecular dynamics and grand canonical Monte Carlo simulations with enhanced sampling techniques. Including polarization leads to a much drier dewetted state and a higher barrier for the transition to the wet state, suggesting more effective hydrophobic gating. The simulations also identify two backbone carbonyls at the bottom of the selectivity filter as good candidates for characterizing the dewetted state using infrared spectroscopies.

Large conductance voltage-and calcium-activated K+ (BK) channel in health and disease

Large Conductance Voltage- and Calcium-activated K+ (BK) channels are transmembrane pore-forming proteins that regulate cell excitability and are also expressed in non-excitable cells. They play a role in regulating vascular tone, neuronal excitability, neurotransmitter release, and muscle contraction. Dysfunction of the BK channel can lead to arterial hypertension, hearing disorders, epilepsy, and ataxia. Here, we provide an overview of BK channel functioning and the implications of its abnormal functioning in various diseases. Understanding the function of BK channels is crucial for comprehending the mechanisms involved in regulating vital physiological processes, both in normal and pathological conditions, controlled by BK. This understanding may lead to the development of therapeutic interventions to address BK channelopathies.

Natural allelic variation modifies acute ethanol response phenotypes in wild strains of C. elegans

BACKGROUND

Genetic variation contributes to the likelihood that an individual will develop an alcohol use disorder (AUD). Traditional laboratory studies in animal models have elucidated the molecular pharmacology of ethanol, but laboratory-derived genetic manipulations rarely model the naturally occurring genetic variation observed in wild populations. Rather, these manipulations are biased toward identifying genes of central importance in the phenotypes. Because changes in such genes can confer selective disadvantages, they are not ideal candidates for carrying AUD risk alleles in humans. We sought to exploit Caenorhabditis elegans to identify allelic variation existing in the wild that modulates ethanol response behaviors.

METHODS

We tested the acute ethanol responses of four strains recently isolated from the wild (JU1511, JU1926, JU1931, and JU1941) and 41 multiparental recombinant inbred lines (mpRILs) derived from them. We assessed locomotion at 10, 30, and 50 min on low and high ethanol concentrations. We performed principal component analyses (PCA) on the different phenotypes, tested for transgressive behavior, calculated heritability, and determined the correlations between behavioral responses.

RESULTS

We observed a range of responses to ethanol across the strains. We detected a low-concentration locomotor activation effect in some of the mpRILs not seen in the laboratory wild-type strain. PCA showed different ethanol response behaviors to be independent. We observed transgressive behavior for many of the measured phenotypes and found that multiple behaviors were uncorrelated. The average broad-sense heritability for all phenotypes was 23.2%.

CONCLUSIONS

Genetic variation significantly affects multiple acute ethanol response behaviors, many of which are independent of one another. This suggests that the genetic variation captured by these strains likely affects multiple biological mechanisms through which ethanol acts. Further study of these strains may allow these distinct mechanisms to be identified.

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.

BK channels of five different subunit combinations underlie the de novo KCNMA1 G375R channelopathy

The molecular basis of a severe developmental and neurological disorder associated with a de novo G375R variant of the tetrameric BK channel is unknown. Here, we address this question by recording from single BK channels expressed to mimic a G375R mutation heterozygous with a WT allele. Five different types of functional BK channels were expressed: 3% were consistent with WT, 12% with homotetrameric mutant, and 85% with three different types of hybrid (heterotetrameric) channels assembled from both mutant and WT subunits. All channel types except WT showed a marked gain-of-function in voltage activation and a smaller decrease-of-function in single-channel conductance, with both changes in function becoming more pronounced as the number of mutant subunits per tetrameric channel increased. The net cellular response from the five different types of channels comprising the molecular phenotype was a shift of -120 mV in the voltage required to activate half of the maximal current through BK channels, giving a net gain-of-function. The WT and homotetrameric mutant channels in the molecular phenotype were consistent with genetic codominance as each displayed properties of a channel arising from only one of the two alleles. The three types of hybrid channels in the molecular phenotype were consistent with partial dominance as their properties were intermediate between those of mutant and WT channels. A model in which BK channels randomly assemble from mutant and WT subunits, with each subunit contributing increments of activation and conductance, approximated the molecular phenotype of the heterozygous G375R mutation.

Cholesterol and PIP2 Modulation of BKCa Channels

Ca2+/voltage-gated, large conductance K+ channels (BKCa) are formed by homotetrameric association of α (slo1) subunits. Their activity, however, is suited to tissue-specific physiology largely due to their association with regulatory subunits (β and γ types), chaperone proteins, localized signaling, and the channel's lipid microenvironment. PIP2 and cholesterol can modulate BKCa activity independently of downstream signaling, yet activating Ca2+i levels and regulatory subunits control ligand action. At physiological Ca2+i and voltages, cholesterol and PIP2 reduce and increase slo1 channel activity, respectively. Moreover, slo1 proteins provide sites that seem to recognize cholesterol and PIP2: seven CRAC motifs in the slo1 cytosolic tail and a string of positively charged residues (Arg329, Lys330, Lys331) immediately after S6, respectively. A model that could explain the modulation of BKCa activity by cholesterol and/or PIP2 is hypothesized. The roles of additional sites, whether in slo1 or BKCa regulatory subunits, for PIP2 and/or cholesterol to modulate BKCa function are also discussed.

KCa-Related Neurological Disorders: Phenotypic Spectrum and Therapeutic Indications

Although potassium channelopathies have been linked to a wide range of neurological conditions, the underlying pathogenic mechanism is not always clear, and a systematic summary of clinical manifestation is absent. Several neurological disorders have been associated with alterations of calcium-activated potassium channels (KCa channels), such as loss- or gain-of-function mutations, post-transcriptional modification, etc. Here, we outlined the current understanding of the molecular and cellular properties of three subtypes of KCa channels, including big conductance KCa channels (BK), small conductance KCa channels (SK), and the intermediate conductance KCa channels (IK). Next, we comprehensively reviewed the loss- or gain-of-function mutations of each KCa channel and described the corresponding mutation sites in specific diseases to broaden the phenotypic-genotypic spectrum of KCa-related neurological disorders. Moreover, we reviewed the current pharmaceutical strategies targeting KCa channels in KCa-related neurological disorders to provide new directions for drug discovery in anti-seizure medication.

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.

Evidence for Dual Activation of IK(M) and IK(Ca) Caused by QO-58 (5-(2,6-Dichloro-5-fluoropyridin-3-yl)-3-phenyl-2-(trifluoromethyl)-1H-pyrazolol[1,5-a]pyrimidin-7-one)

QO-58 (5-(2,6-dichloro-5-fluoropyridin-3-yl)-3-phenyl-2-(trifluoromethyl)-1H-pyrazolol[1,5-a]pyrimidin-7-one) has been regarded to be an activator of K_V7 channels with analgesic properties. However, whether and how the presence of this compound can result in any modifications of other types of membrane ion channels in native cells are not thoroughly investigated. In this study, we investigated its perturbations on M-type K⁺ current (I_K(M)), Ca²⁺-activated K⁺ current (I_K(Ca)), large-conductance Ca²⁺-activated K⁺ (BK_Ca) channels, and erg-mediated K⁺ current (I_K(erg)) identified from pituitary tumor (GH₃) cells. Addition of QO-58 can increase the amplitude of I_K(M) and I_K(Ca) in a concentration-dependent fashion, with effective EC₅₀ of 3.1 and 4.2 μM, respectively. This compound could shift the activation curve of I_K(M) toward a leftward direction with being void of changes in the gating charge. The strength in voltage-dependent hysteresis (V_hys) of I_K(M) evoked by upright triangular ramp pulse (V_ramp) was enhanced by adding QO-58. The probabilities of M-type K⁺ (K_M) channels that will be open increased upon the exposure to QO-58, although no modification in single-channel conductance was seen. Furthermore, GH₃-cell exposure to QO-58 effectively increased the amplitude of I_K(Ca) as well as enhanced the activity of BK_Ca channels. Under inside-out configuration, QO-58, applied at the cytosolic leaflet of the channel, activated BK_Ca-channel activity, and its increase could be attenuated by further addition of verruculogen, but not by linopirdine (10 μM). The application of QO-58 could lead to a leftward shift in the activation curve of BK_Ca channels with neither change in the gating charge nor in single-channel conductance. Moreover, cell exposure of QO-58 (10 μM) resulted in a minor suppression of I_K(erg) amplitude in response to membrane hyperpolarization. The docking results also revealed that there are possible interactions of the QO-58 molecule with the KCNQ or K_Ca1.1 channel. Overall, dual activation of I_K(M) and I_K(Ca) caused by the presence of QO-58 eventually may have high impacts on the functional activity (e.g., anti-nociceptive effect) residing in electrically excitable cells. Care must be exercised when interpreting data generated with QO-58 as it is not entirely KCNQ/K_V7 selective.

Neuronal mechanism of a BK channelopathy in absence epilepsy and dyskinesia

BK channelopathy has been increasingly implicated in diverse neurological disorders, including epilepsy and movement, cognitive, and neurodevelopmental disorders. However, precision medicine to treat BK channelopathy is lacking. We characterized a mouse model carrying a gain-of-function BK channelopathy D434G from a large family of patients with absence epilepsy and paroxysmal dyskinesia. The BK-D434G mice manifest the clinical features of absence seizures and exhibit severe locomotor defects including involuntary dyskinesia-like behavior. Pharmacological inhibition of BK channels suppresses neuronal hyperactivity and mitigates absence seizure and the locomotor defects. The BK-D434G mice thus serve as a model to understand the pathogenic mechanisms of absence epilepsy and dyskinesia. Our study also suggests that BK inhibition is a promising strategy for treating BK gain-of-function channelopathy.

A growing number of gain-of-function (GOF) BK channelopathies have been identified in patients with epilepsy and movement disorders. Nevertheless, the underlying pathophysiology and corresponding therapeutics remain obscure. Here, we utilized a knock-in mouse model carrying human BK-D434G channelopathy to investigate the neuronal mechanism of BK GOF in the pathogenesis of epilepsy and dyskinesia. The BK-D434G mice manifest the clinical features of absence epilepsy and exhibit severe motor deficits and dyskinesia-like behaviors. The cortical pyramidal neurons and cerebellar Purkinje cells from the BK-D434G mice show hyperexcitability, which likely contributes to the pathogenesis of absence seizures and paroxysmal dyskinesia. A BK channel blocker, paxilline, potently suppresses BK-D434G–induced hyperexcitability and effectively mitigates absence seizures and locomotor deficits in mice. Our study thus uncovered a neuronal mechanism of BK GOF in absence epilepsy and dyskinesia. Our findings also suggest that BK inhibition is a promising therapeutic strategy for mitigating BK GOF-induced neurological disorders.