Conserved gating hinge in ligand- and voltage-dependent K+ 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.

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.

A mutation in the S6 segment of the KvAP channel changes the secondary structure and alters ion channel activity in a lipid bilayer membrane

The peptide segment S6 is known to form the inner lining of the voltage-gated K⁺ channel KvAP (potassium channel of archaea-bacterium, Aeropyrum pernix). In our previous work, it has been demonstrated that S6 itself can form an ion channel on a bilayer lipid membrane (BLM). In the present work, the role of a specific amino acid sequence 'LIG' in determining the secondary structure of S6 has been investigated. For this purpose, 22-residue synthetic peptides named S6-Wild (S6W) and S6-Mutant (S6M) were used. Sequences of these peptides are similar except that the two amino acids isoleucine and glycine of the wild peptide interchanged in the mutant peptide. Channel forming capabilities of both the peptides were checked electro-physiologically on BLM composed of DPhPC and cholesterol. Bilayer electrophysiological experiments showed that the conductance of S6M is higher than that of S6W. Significant differences in the current versus voltage (I-V) plot, open probability, and gating characteristics were observed. Interestingly, two sub-types of channels, S6M Type 1 and Type 2, were identified in S6M differing in conductances and open probability patterns. Circular dichroism (CD) spectroscopy indicated that the secondary structures of the two peptides are different in phosphatidyl choline/asolectin liposomes and 1% SDS detergent. Reduced helicity of S6M was also noticed in membrane mimetic liposomes and 1% SDS detergent micelles. These results are interpreted in view of the difference in hydrophobicity of the two amino acids, isoleucine and glycine. It is concluded that the 'LIG' stretch regulates the structure and pore-forming ability of the S6 peptide.

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.

Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein

Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.

Inferring functional units in ion channel pores via relative entropy

Coarse-grained protein models approximate the first-principle physical potentials. Among those modeling approaches, the relative entropy framework yields promising and physically sound results, in which a mapping from the target protein structure and dynamics to a model is defined and subsequently adjusted by an entropy minimization of the model parameters. Minimization of the relative entropy is equivalent to maximization of the likelihood of reproduction of (configurational ensemble) observations by the model. In this study, we extend the relative entropy minimization procedure beyond parameter fitting by a second optimization level, which identifies the optimal mapping to a (dimension-reduced) topology. We consider anisotropic network models of a diverse set of ion channels and assess our findings by comparison to experimental results.

Calcium channel gating

Tuned calcium entry through voltage-gated calcium channels is a key requirement for many cellular functions. This is ensured by channel gates which open during membrane depolarizations and seal the pore at rest. The gating process is determined by distinct sub-processes: movement of voltage-sensing domains (charged S4 segments) as well as opening and closure of S6 gates. Neutralization of S4 charges revealed that pore opening of CaV1.2 is triggered by a "gate releasing" movement of all four S4 segments with activation of IS4 (and IIIS4) being a rate-limiting stage. Segment IS4 additionally plays a crucial role in channel inactivation. Remarkably, S4 segments carrying only a single charged residue efficiently participate in gating. However, the complete set of S4 charges is required for stabilization of the open state. Voltage clamp fluorometry, the cryo-EM structure of a mammalian calcium channel, biophysical and pharmacological studies, and mathematical simulations have all contributed to a novel interpretation of the role of voltage sensors in channel opening, closure, and inactivation. We illustrate the role of the different methodologies in gating studies and discuss the key molecular events leading CaV channels to open and to close.

Transmembrane Segment XI of the Na+/H+ Antiporter of S. pombe is a Critical Part of the Ion Translocation Pore

The Na⁺/H⁺ exchanger of the plasma membrane of S. pombe (SpNHE1) removes intracellular sodium in exchange for an extracellular proton. We examined the structure and functional role of amino acids 360–393 of putative transmembrane (TM) segment XI of SpNHE1. Structural analysis suggested that it had a helical propensity over amino acids 360–368, an extended region from 369–378 and was helical over amino acids 379–386. TM XI was sensitive to side chain alterations. Mutation of eight amino acids to alanine resulted in loss of one or both of LiCl or NaCl tolerance when re-introduced into SpNHE1 deficient S. pombe. Mutation of seven other amino acids had minor effects. Analysis of structure and functional mutations suggested that Glu³⁶¹ may be involved in cation coordination on the cytoplasmic face of the protein with a negative charge in this position being important. His³⁶⁷, Ile³⁷¹ and Gly³⁷² were important in function. Ile³⁷¹ may have important hydrophobic interactions with other residues and Gly³⁷² may be important in maintaining an extended conformation. Several residues from Val³⁷⁷ to Leu³⁸⁴ are important in function possibly involved in hydrophobic interactions with other amino acids. We suggest that TM XI forms part of the ion translocation core of this Na⁺/H⁺ exchanger.

Side chain flexibility and coupling between the S4-S5 linker and the TRP domain in thermo-sensitive TRP channels: Insights from protein modeling

The transient receptor potential (TRP) superfamily is subdivided into several subfamilies on the basis of sequence similarity, which is highly heterogeneous but shows a molecular architecture that resembles the one present in members of the Kv channel superfamily. Because of this diversity, they produce a large variety of channels with different gating and permeability properties. Elucidation of these particular features necessarily requires comparative studies based on structural and functional data. The present study aims to compilate, analyze, and determine, in a coherent way, the relationship between intrinsic side-chain flexibility and the allosteric coupling in members of the TRPV, TRPM, and TRPC families. Based on the recently determined structures of TRPV1 and TRPV2, we have generated protein models for single subunits of TRPV5, TRPM8, and TRPC5 channels. With these models, we focused our attention on the apparently crucial role of the GP dipeptide at the center of the S4-S5 linker and discussed its role in the interaction with the TRP domain, specifically with the highly-conserved Trp during this coupling. Our analysis suggests an important role of the S4-S5L flexibility in the thermosensitivity, where heat-activated channels possess rigid S4-S5 linkers, whereas cold-activated channels have flexible ones. Finally, we also present evidence of the key interaction between the conserved Trp residue of the TRP box and of several residues in the S4-S5L, importantly the central Pro. Proteins 2017; 85:630-646. © 2016 Wiley Periodicals, Inc.

Conformational Response of Influenza A M2 Transmembrane Domain to Amantadine Drug Binding at Low pH (pH 5.5)

The M2 protein from influenza A plays important roles in its viral cycle. It contains a single transmembrane helix, which oligomerizes into a homotetrameric proton channel that conducts in the low-pH environment of the host-cell endosome and Golgi apparatus, leading to virion uncoating at an early stage of infection. We studied conformational rearrangements that occur in the M2 core transmembrane domain residing on the lipid bilayer, flanked by juxtamembrane residues (M2TMD21-49 fragment), upon its interaction with amantadine drug at pH 5.5 when M2 is conductive. We also tested the role of specific mutation and lipid chain length. Electron spin resonance (ESR) spectroscopy and electron microscopy were applied to M2TMD21-49, labeled at the residue L46C with either nitroxide spin-label or Nanogold® reagent, respectively. Electron microscopy confirmed that M2TMD21-49 reconstituted into DOPC/POPS at 1:10,000 peptide-to-lipid molar ratio (P/L) either with or without amantadine, is an admixture of monomers, dimers, and tetramers, confirming our model based on a dimer intermediate in the assembly of M2TMD21-49. As reported by double electron-electron resonance (DEER), in DOPC/POPS membranes amantadine shifts oligomer equilibrium to favor tetramers, as evidenced by an increase in DEER modulation depth for P/L's ranging from 1:18,000 to 1:160. Furthermore, amantadine binding shortens the inter-spin distances (for nitroxide labels) by 5-8 Å, indicating drug induced channel closure on the C-terminal side. No such effect was observed for the thinner membrane of DLPC/DLPS, emphasizing the role of bilayer thickness. The analysis of continuous wave (cw) ESR spectra of spin-labeled L46C residue provides additional support to a more compact helix bundle in amantadine-bound M2TMD 21-49 through increased motional ordering. In contrast to wild-type M2TMD21-49, the amantadine-bound form does not exhibit noticeable conformational changes in the case of G34A mutation found in certain drug-resistant influenza strains. Thus, the inhibited M2TMD21-49 channel is a stable tetramer with a closed C-terminal exit pore. This work is aimed at contributing to the development of structure-based anti-influenza pharmaceuticals.

Allosteric coupling between proximal C-terminus and selectivity filter is facilitated by the movement of transmembrane segment 4 in TREK-2 channel

TREK-2, a member of two-pore-domain potassium channel family, regulates cellular excitability in response to diverse stimuli. However, how such stimuli control channel function remains unclear. Here, by characterizing the responses of cytosolic proximal C-terminus deletant (ΔpCt) and transmembrane segment 4 (M4)-glycine hinge mutant (G312A) to 2-Aminoethoxydiphenyl borate (2-APB), an activator of TREK-2, we show that the transduction initiated from pCt domain is allosterically coupled with the conformation of selectivity filter (SF) via the movements of M4, without depending on the original status of SF. Moreover, ΔpCt and G312A also exhibited blunted responses to extracellular alkalization, a model to induce SF conformational transition. These results suggest that the coupling between pCt domain and SF is bidirectional, and M4 movements are involved in both processes. Further mechanistic exploration reveals that the function of Phe316, a residue close to the C-terminus of M4, is associated with such communications. However, unlike TREK-2, M4-hinge of TREK-1 only controls the transmission from pCt to SF, rather than SF conformational changes triggered by pHo changes. Together, our findings uncover the unique gating properties of TREK-2, and elucidate the mechanisms for how the extracellular and intracellular stimuli harness the pore gating allosterically.

Extrapolating microdomain Ca(2+) dynamics using BK channels as a Ca(2+) sensor

Ca²⁺ ions play crucial roles in mediating physiological and pathophysiological processes, yet Ca²⁺ dynamics local to the Ca²⁺ source, either from influx via calcium permeable ion channels on plasmic membrane or release from internal Ca²⁺ stores, is difficult to delineate. Large-conductance calcium-activated K⁺ (BK-type) channels, abundantly distribute in excitable cells and often localize to the proximity of voltage-gated Ca²⁺ channels (VGCCs), spatially enabling the coupling of the intracellular Ca²⁺ signal to the channel gating to regulate membrane excitability and spike firing patterns. Here we utilized the sensitivity and dynamic range of BK to explore non-uniform Ca²⁺ local transients in the microdomain of VGCCs. Accordingly, we applied flash photolysis of caged Ca²⁺ to activate BK channels and determine their intrinsic sensitivity to Ca²⁺. We found that uncaging Ca²⁺ activated biphasic BK currents with fast and slow components (time constants being τ_f ≈ 0.2 ms and τ_s ≈ 10 ms), which can be accounted for by biphasic Ca²⁺ transients following light photolysis. We estimated the Ca²⁺-binding rate constant k_b (≈1.8 × 10⁸ M⁻¹s⁻¹) for mSlo1 and further developed a model in which BK channels act as a calcium sensor capable of quantitatively predicting local microdomain Ca²⁺ transients in the vicinity of VGCCs during action potentials.

Paxilline inhibits BK channels by an almost exclusively closed-channel block mechanism

Paxilline, a tremorogenic fungal alkaloid, potently inhibits large conductance Ca(2+)- and voltage-activated K(+) (BK)-type channels, but little is known about the mechanism underlying this inhibition. Here we show that inhibition is inversely dependent on BK channel open probability (Po), and is fully relieved by conditions that increase Po, even in the constant presence of paxilline. Manipulations that shift BK gating to more negative potentials reduce inhibition by paxilline in accordance with the increase in channel Po. Measurements of Po times the number of channels at negative potentials support the idea that paxilline increases occupancy of closed states, effectively reducing the closed-open equilibrium constant, L(0). Gating current measurements exclude an effect of paxilline on voltage sensors. Steady-state inhibition by multiple paxilline concentrations was determined for four distinct equilibration conditions, each with a distinct Po. The IC50 for paxilline shifted from around 10 nM when channels were largely closed to near 10 µM as maximal Po was approached. Model-dependent analysis suggests a mechanism of inhibition in which binding of a single paxilline molecule allosterically alters the intrinsic L(0) favoring occupancy of closed states, with affinity for the closed conformation being >500-fold greater than affinity for the open conformation. The rate of inhibition of closed channels was linear up through 2 µM paxilline, with a slope of 2 × 10(6) M(-1)s(-1). Paxilline inhibition was hindered by either the bulky cytosolic blocker, bbTBA, or by concentrations of cytosolic sucrose that hinder ion permeation. However, paxilline does not hinder MTSET modification of the inner cavity residue, A313C. We conclude that paxilline binds more tightly to the closed conformation, favoring occupancy of closed-channel conformations, and propose that it binds to a superficial position near the entrance to the central cavity, but does not hinder access of smaller molecules to this cavity.

Identification of the Conformational transition pathway in PIP2 Opening Kir Channels

The gating of Kir channels depends critically on phosphatidylinositol 4,5-bisphosphate (PIP2), but the detailed mechanism by which PIP2 regulates Kir channels remains obscure. Here, we performed a series of Targeted molecular dynamics simulations on the full-length Kir2.1 channel and, for the first time, were able to achieve the transition from the closed to the open state. Our data show that with the upward motion of the cytoplasmic domain (CTD) the structure of the C-Linker changes from a loop to a helix. The twisting of the C-linker triggers the rotation of the CTD, which induces a small downward movement of the CTD and an upward motion of the slide helix toward the membrane that pulls the inner helix gate open. At the same time, the rotation of the CTD breaks the interaction between the CD- and G-loops thus releasing the G-loop. The G-loop then bounces away from the CD-loop, which leads to the opening of the G-loop gate and the full opening of the pore. We identified a series of interaction networks, between the N-terminus, CD loop, C linker and G loop one by one, which exquisitely regulates the global conformational changes during the opening of Kir channels by PIP2.

Channel Gating Dependence on Pore Lining Helix Glycine Residues in Skeletal Muscle Ryanodine Receptor

Type 1 ryanodine receptors (RyR1s) release Ca(2+) from the sarcoplasmic reticulum to initiate skeletal muscle contraction. The role of RyR1-G4934 and -G4941 in the pore-lining helix in channel gating and ion permeation was probed by replacing them with amino acid residues of increasing side chain volume. RyR1-G4934A, -G4941A, and -G4941V mutant channels exhibited a caffeine-induced Ca(2+) release response in HEK293 cells and bound the RyR-specific ligand [(3)H]ryanodine. In single channel recordings, significant differences in the number of channel events and mean open and close times were observed between WT and RyR1-G4934A and -G4941A. RyR1-G4934A had reduced K(+) conductance and ion selectivity compared with WT. Mutations further increasing the side chain volume at these positions (G4934V and G4941I) resulted in reduced caffeine-induced Ca(2+) release in HEK293 cells, low [(3)H]ryanodine binding levels, and channels that were not regulated by Ca(2+) and did not conduct Ca(2+) in single channel measurements. Computational predictions of the thermodynamic impact of mutations on protein stability indicated that although the G4934A mutation was tolerated, the G4934V mutation decreased protein stability by introducing clashes with neighboring amino acid residues. In similar fashion, the G4941A mutation did not introduce clashes, whereas the G4941I mutation resulted in intersubunit clashes among the mutated isoleucines. Co-expression of RyR1-WT with RyR1-G4934V or -G4941I partially restored the WT phenotype, which suggested lessening of amino acid clashes in heterotetrameric channel complexes. The results indicate that both glycines are important for RyR1 channel function by providing flexibility and minimizing amino acid clashes.

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels

Prokaryotic inwardly rectifying (KirBac) potassium channels are homologous to mammalian Kir channels. Their activity is controlled by dynamical conformational changes that regulate ion flow through a central pore. Understanding the dynamical rearrangements of Kir channels during gating requires high-resolution structure information from channels crystallized in different conformations and insight into the transition steps, which are difficult to access experimentally. In this study, we use MD simulations on wild type KirBac1.1 and an activatory mutant to investigate activation gating of KirBac channels. Full atomistic MD simulations revealed that introducing glutamate in position 143 causes significant widening at the helix bundle crossing gate, enabling water flux into the cavity. Further, global rearrangements including a twisting motion as well as local rearrangements at the subunit interface in the cytoplasmic domain were observed. These structural rearrangements are similar to recently reported KirBac3.1 crystal structures in closed and open conformation, suggesting that our simulations capture major conformational changes during KirBac1.1 opening. In addition, an important role of protein-lipid interactions during gating was observed. Slide-helix and C-linker interactions with lipids were strengthened during activation gating.

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.

Transmembrane helix straightening and buckling underlies activation of mechanosensitive and thermosensitive K(2P) channels

Mechanical and thermal activation of ion channels is central to touch, thermosensation, and pain. The TRAAK/TREK K(2P) potassium channel subfamily produces background currents that alter neuronal excitability in response to pressure, temperature, signaling lipids, and anesthetics. How such diverse stimuli control channel function is unclear. Here we report structures of K(2P)4.1 (TRAAK) bearing C-type gate-activating mutations that reveal a tilting and straightening of the M4 inner transmembrane helix and a buckling of the M2 transmembrane helix. These conformational changes move M4 in a direction opposite to that in classical potassium channel activation mechanisms and open a passage lateral to the pore that faces the lipid bilayer inner leaflet. Together, our findings uncover a unique aspect of K(2P) modulation, indicate a means for how the K(2P) C-terminal cytoplasmic domain affects the C-type gate which lies ∼40Å away, and suggest how lipids and bilayer inner leaflet deformations may gate the channel.

Conservation analysis of residues in the S4-S5 linker and the terminal part of the S5-P-S6 pore modulus in Kv and HCN channels: flexible determinants for the electromechanical coupling

Protein mobility is important to achieve protein function. Intrinsic flexibility associated with motion underlies this important issue and the analysis of side chain flexibility gives insights to understand it. In this work, the S5-P-S6 pore modulus (PM) of members of Kv and HCN channels was examined by a combination of sequence alignment, residue composition analysis, and intrinsic side chain flexibility. The PM sequences were organized as a database that was used to reveal and correlate the functional diversity of each analyzed family. Specifically, we focused our attention on the crucial role of the S4-S5 linker and its well-described interaction with the S6 T during the electromechanical coupling. Our analysis suggests the presence of a Gly-hinge in the middle of the S4-S5 linkers. This apparent Gly-hinge links a flexible N-terminal segment with a rigid C-terminal one, although in Kv7 channels, the latter segment is even more flexible. Instead, HCN channels exhibit a putative Thr-hinge and is rich in aromatic residues, in consequence, their linker is more rigid. Concerning S6, we confirm the presence of the two flexible kinks previously described and we provide the complete segmental flexibility profiles for the different families. Our results are discussed in terms of the relation between residue composition, conservation, and local conformational flexibility. This provides important insights to understand and differentiate the characteristic gating properties of these channels as well as their implications in cell physiology.

Coupling between the voltage-sensing and pore domains in a voltage-gated potassium channel

Voltage-dependent potassium (Kv), sodium (Nav), and calcium channels open and close in response to changes in transmembrane (TM) potential, thus regulating cell excitability by controlling ion flow across the membrane. An outstanding question concerning voltage gating is how voltage-induced conformational changes of the channel voltage-sensing domains (VSDs) are coupled through the S4-S5 interfacial linking helices to the opening and closing of the pore domain (PD). To investigate the coupling between the VSDs and the PD, we generated a closed Kv channel configuration from Aeropyrum pernix (KvAP) using atomistic simulations with experiment-based restraints on the VSDs. Full closure of the channel required, in addition to the experimentally determined TM displacement, that the VSDs be displaced both inwardly and laterally around the PD. This twisting motion generates a tight hydrophobic interface between the S4-S5 linkers and the C-terminal ends of the pore domain S6 helices in agreement with available experimental evidence.

G protein modulation of K2P potassium channel TASK-2 : a role of basic residues in the C terminus domain

TASK-2 (K2P5.1) is a background K(+) channel opened by extra- or intracellular alkalinisation that plays a role in renal bicarbonate handling, central chemoreception and cell volume regulation. Here, we present results that suggest that TASK-2 is also modulated by Gβγ subunits of heterotrimeric G protein. TASK-2 was strongly inhibited when GTP-γ-S was used as a replacement for intracellular GTP. No inhibition was present using GDP-β-S instead. Purified Gβγ introduced intracellularly also inhibited TASK-2 independently of whether GTP or GDP-β-S was present. The effects of GTP-γ-S and Gβγ subunits were abolished by neutralisation of TASK-2 C terminus double lysine residues K257-K258 or K296-K297. Use of membrane yeast two hybrid (MYTH) experiments and immunoprecipitation assays using tagged proteins gave evidence for a physical interaction between Gβ1 and Gβ2 subunits and TASK-2, in agreement with expression of these subunits in proximal tubule cells. Co-immunoprecipitation was impeded by mutating C terminus K257-K258 (but not K296-K297) to alanines. Gating by extra- or intracellular pH was unaltered in GTP-γ-S-insensitive TASK-2-K257A-K258A mutant. Shrinking TASK-2-expressing cells in hypertonic solution decreased the current to 36 % of its initial value. The same manoeuvre had a significantly diminished effect on TASK-2-K257A-K258A- or TASK-2-K296-K297-expressing cells, or in cells containing intracellular GDP-β-S. Our data are compatible with the concept that TASK-2 channels are modulated by Gβγ subunits of heterotrimeric G protein. We propose that this modulation is a novel way in which TASK-2 can be tuned to its physiological functions.

On the origin of large flexibility of P-glycoprotein in the inward-facing state

P-glycoprotein (Pgp) is one of the most biomedically relevant transporters in the ATP binding cassette (ABC) superfamily due to its involvement in developing multidrug resistance in cancer cells. Employing molecular dynamics simulations and double electron-electron resonance spectroscopy, we have investigated the structural dynamics of membrane-bound Pgp in the inward-facing state and found that Pgp adopts an unexpectedly wide range of conformations, highlighted by the degree of separation between the two nucleotide-binding domains (NBDs). The distance between the two NBDs in the equilibrium simulations covers a range of at least 20 Å, including, both, more open and more closed NBD configurations than the crystal structure. The double electron-electron resonance measurements on spin-labeled Pgp mutants also show wide distributions covering both longer and shorter distances than those observed in the crystal structure. Based on structural and sequence analyses, we propose that the transmembrane domains of Pgp might be more flexible than other structurally known ABC exporters. The structural flexibility of Pgp demonstrated here is not only in close agreement with, but also helps rationalize, the reported high NBD fluctuations in several ABC exporters and possibly represents a fundamental difference in the transport mechanism between ABC exporters and ABC importers. In addition, during the simulations we have captured partial entrance of a lipid molecule from the bilayer into the lumen of Pgp, reaching the putative drug binding site. The location of the protruding lipid suggests a putative pathway for direct drug recruitment from the membrane.

Investigations of the contribution of a putative glycine hinge to ryanodine receptor channel gating

Ryanodine receptor channels (RyR) are key components of striated muscle excitation-contraction coupling, and alterations in their function underlie both inherited and acquired disease. A full understanding of the disease process will require a detailed knowledge of the mechanisms and structures involved in RyR function. Unfortunately, high-resolution structural data, such as exist for K⁺-selective channels, are not available for RyR. In the absence of these data, we have used modeling to identify similarities in the structural elements of K⁺ channel pore-forming regions and postulated equivalent regions of RyR. This has identified a sequence of residues in the cytosolic cavity-lining transmembrane helix of RyR (G⁴⁸⁶⁴LIIDA⁴⁸⁶⁹ in RyR2) analogous to the glycine hinge motif present in many K⁺ channels. Gating in these K⁺ channels can be disrupted by substitution of residues for the hinge glycine. We investigated the involvement of glycine 4864 in RyR2 gating by monitoring properties of recombinant human RyR2 channels in which this glycine is replaced by residues that alter gating in K⁺ channels. Our data demonstrate that introducing alanine at position 4864 produces no significant change in RyR2 function. In contrast, function is altered when glycine 4864 is replaced by either valine or proline, the former preventing channel opening and the latter modifying both ion translocation and gating. Our studies reveal novel information on the structural basis of RyR gating, identifying both similarities with, and differences from, K⁺ channels. Glycine 4864 is not absolutely required for channel gating, but some flexibility at this point in the cavity-lining transmembrane helix is necessary for normal RyR function.

Being flexible: the voltage-controllable activation gate of kv channels

Kv channels form voltage-dependent potassium selective pores in the outer cell membrane and are composed out of four α-subunits, each having six membrane-spanning α-helices (S1–S6). The α-subunits tetramerize such that the S5–S6 pore domains co-assemble into a centrally located K⁺ pore which is surrounded by four operational voltage-sensing domains (VSD) that are each formed by the S1–S4 segments. Consequently, each subunit is capable of responding to changes in membrane potential and dictates whether the pore should be conductive or not. K⁺ permeation through the pore can be sealed off by two separate gates in series: (a) at the inner S6 bundle crossing (BC gate) and (b) at the level of the selectivity filter (SF gate) located at the extracellular entrance of the pore. Within the last years a general consensus emerged that a direct communication between the S4S5-linker and the bottom part of S6 (S6_c) constitutes the coupling with the VSD thus making the BC gate the main voltage-controllable activation gate. While the BC gate listens to the VSD, the SF changes its conformation depending on the status of the BC gate. Through the eyes of an entering K⁺ ion, the operation of the BC gate apparatus can be compared with the iris-like motion of the diaphragm from a camera whereby its diameter widens. Two main gating motions have been proposed to create this BC gate widening: (1) tilting of the helix whereby the S6 converts from a straight α-helix to a tilted one or (2) swiveling of the S6_c whereby the S6 remains bent. Such motions require a flexible hinge that decouples the pre- and post-hinge segment. Roughly at the middle of the S6 there exists a highly conserved glycine residue and a tandem proline motif that seem to fulfill the role of a gating hinge which allows for tilting/swiveling/rotations of the post-hinge S6 segment. In this review we delineate our current view on the operation of the BC gate for controlling K⁺ permeation in Kv channels.

hERG K+ Channels: Structure, Function, and Clinical Significance

The human ether-a-go-go related gene (hERG) encodes the pore-forming subunit of the rapid component of the delayed rectifier K+ channel, Kv11.1, which are expressed in the heart, various brain regions, smooth muscle cells, endocrine cells, and a wide range of tumor cell lines. However, it is the role that Kv11.1 channels play in the heart that has been best characterized, for two main reasons. First, it is the gene product involved in chromosome 7-associated long QT syndrome (LQTS), an inherited disorder associated with a markedly increased risk of ventricular arrhythmias and sudden cardiac death. Second, blockade of Kv11.1, by a wide range of prescription medications, causes drug-induced QT prolongation with an increase in risk of sudden cardiac arrest. In the first part of this review, the properties of Kv11.1 channels, including biogenesis, trafficking, gating, and pharmacology are discussed, while the second part focuses on the pathophysiology of Kv11.1 channels.

On the simple random-walk models of ion-channel gate dynamics reflecting long-term memory

Several approaches to ion-channel gating modelling have been proposed. Although many models describe the dwell-time distributions correctly, they are incapable of predicting and explaining the long-term correlations between the lengths of adjacent openings and closings of a channel. In this paper we propose two simple random-walk models of the gating dynamics of voltage and Ca(2+)-activated potassium channels which qualitatively reproduce the dwell-time distributions, and describe the experimentally observed long-term memory quite well. Biological interpretation of both models is presented. In particular, the origin of the correlations is associated with fluctuations of channel mass density. The long-term memory effect, as measured by Hurst R/S analysis of experimental single-channel patch-clamp recordings, is close to the behaviour predicted by our models. The flexibility of the models enables their use as templates for other types of ion channel.

Structure of a KirBac potassium channel with an open bundle crossing indicates a mechanism of channel gating

KirBac channels are prokaryotic homologs of mammalian inwardly rectifying (Kir) potassium channels, and recent crystal structures of both Kir and KirBac channels have provided major insight into their unique structural architecture. However, all of the available structures are closed at the helix bundle crossing, and therefore the structural mechanisms that control opening of their primary activation gate remain unknown. In this study, we engineered the inner pore-lining helix (TM2) of KirBac3.1 to trap the bundle crossing in an apparently open conformation and determined the crystal structure of this mutant channel to 3.05 Å resolution. Contrary to previous speculation, this new structure suggests a mechanistic model in which rotational 'twist' of the cytoplasmic domain is coupled to opening of the bundle-crossing gate through a network of inter- and intrasubunit interactions that involve the TM2 C-linker, slide helix, G-loop and the CD loop.

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.

Uncoupling charge movement from channel opening in voltage-gated potassium channels by ruthenium complexes

The Kv2.1 channel generates a delayed-rectifier current in neurons and is responsible for modulation of neuronal spike frequency and membrane repolarization in pancreatic β-cells and cardiomyocytes. As with other tetrameric voltage-activated K(+)-channels, it has been proposed that each of the four Kv2.1 voltage-sensing domains activates independently upon depolarization, leading to a final concerted transition that causes channel opening. The mechanism by which voltage-sensor activation is coupled to the gating of the pore is still not understood. Here we show that the carbon-monoxide releasing molecule 2 (CORM-2) is an allosteric inhibitor of the Kv2.1 channel and that its inhibitory properties derive from the CORM-2 ability to largely reduce the voltage dependence of the opening transition, uncoupling voltage-sensor activation from the concerted opening transition. We additionally demonstrate that CORM-2 modulates Shaker K(+)-channels in a similar manner. Our data suggest that the mechanism of inhibition by CORM-2 may be common to voltage-activated channels and that this compound should be a useful tool for understanding the mechanisms of electromechanical coupling.

Gating of a pH-sensitive K(2P) potassium channel by an electrostatic effect of basic sensor residues on the selectivity filter

K⁺ channels share common selectivity characteristics but exhibit a wide diversity in how they are gated open. Leak K_2P K⁺ channels TASK-2, TALK-1 and TALK-2 are gated open by extracellular alkalinization. The mechanism for this alkalinization-dependent gating has been proposed to be the neutralization of the side chain of a single arginine (lysine in TALK-2) residue near the pore of TASK-2, which occurs with the unusual pK_a of 8.0. We now corroborate this hypothesis by transplanting the TASK-2 extracellular pH (pH_o) sensor in the background of a pH_o-insensitive TASK-3 channel, which leads to the restitution of pH_o-gating. Using a concatenated channel approach, we also demonstrate that for TASK-2 to open, pH_o sensors must be neutralized in each of the two subunits forming these dimeric channels with no apparent cross-talk between the sensors. These results are consistent with adaptive biasing force analysis of K⁺ permeation using a model selectivity filter in wild-type and mutated channels. The underlying free-energy profiles confirm that either a doubly or a singly charged pH_o sensor is sufficient to abolish ion flow. Atomic detail of the associated mechanism reveals that, rather than a collapse of the pore, as proposed for other K_2P channels gated at the selectivity filter, an increased height of the energetic barriers for ion translocation accounts for channel blockade at acid pH_o. Our data, therefore, strongly suggest that a cycle of protonation/deprotonation of pH_o-sensing arginine 224 side chain gates the TASK-2 channel by electrostatically tuning the conformational stability of its selectivity filter.

Physicochemical properties of pore residues predict activation gating of Ca V1.2: a correlation mutation analysis

Single point mutations in pore-forming S6 segments of calcium channels may transform a high-voltage-activated into a low-voltage-activated channel, and resulting disturbances in calcium entry may cause channelopathies (Hemara-Wahanui et al., Proc Natl Acad Sci U S A 102(21):7553–7558, 16). Here we ask the question how physicochemical properties of amino acid residues in gating-sensitive positions on S6 segments determine the threshold of channel activation of Ca_V1.2. Leucine in segment IS6 (L434) and a newly identified activation determinant in segment IIIS6 (G1193) were mutated to a variety of amino acids. The induced leftward shifts of the activation curves and decelerated current activation and deactivation suggest a destabilization of the closed and a stabilisation of the open channel state by most mutations. A selection of 17 physicochemical parameters (descriptors) was calculated for these residues and examined for correlation with the shifts of the midpoints of the activation curve (ΔV_act). ΔV_act correlated with local side-chain flexibility in position L434 (IS6), with the polar accessible surface area of the side chain in position G1193 (IIIS6) and with hydrophobicity in position I781 (IIS6). Combined descriptor analysis for positions I781 and G1193 revealed that additional amino acid properties may contribute to conformational changes during the gating process. The identified physicochemical properties in the analysed gating-sensitive positions (accessible surface area, side-chain flexibility, and hydrophobicity) predict the shifts of the activation curves of Ca_V1.2.

Glycine311, a determinant of paxilline block in BK channels: a novel bend in the BK S6 helix

The tremorogenic fungal metabolite, paxilline, is widely used as a potent and relatively specific blocker of Ca²⁺- and voltage-activated Slo1 (or BK) K⁺ channels. The pH-regulated Slo3 K⁺ channel, a Slo1 homologue, is resistant to blockade by paxilline. Taking advantage of the marked differences in paxilline sensitivity and the homology between subunits, we have examined the paxilline sensitivity of a set of chimeric Slo1/Slo3 subunits. Paxilline sensitivity is associated with elements of the S5–P loop–S6 module of the Slo1 channel. Replacement of the Slo1 S5 segment or the second half of the P loop results in modest changes in paxilline sensitivity. Replacing the Slo1 S6 segment with the Slo3 sequence abolishes paxilline sensitivity. An increase in paxilline affinity and changes in block kinetics also result from replacing the first part of the Slo1 P loop, the so-called turret, with Slo3 sequence. The Slo1 and Slo3 S6 segments differ at 10 residues. Slo1-G311S was found to markedly reduce paxilline block. In constructs with a Slo3 S6 segment, S300G restored paxilline block, but most effectively when paired with a Slo1 P loop. Other S6 residues differing between Slo1 and Slo3 had little influence on paxilline block. The involvement of Slo1 G311 in paxilline sensitivity suggests that paxilline may occupy a position within the central cavity or access its blocking position through the central cavity. To explain the differences in paxilline sensitivity between Slo1 and Slo3, we propose that the G311/S300 position in Slo1 and Slo3 underlies a structural difference between subunits in the bend of S6, which influences the occupancy by paxilline.

Separate gating mechanisms mediate the regulation of K2P potassium channel TASK-2 by intra- and extracellular pH

TASK-2 (KCNK5 or K(2P)5.1) is a background K(+) channel that is opened by extracellular alkalinization and plays a role in renal bicarbonate reabsorption and central chemoreception. Here, we demonstrate that in addition to its regulation by extracellular protons (pH(o)) TASK-2 is gated open by intracellular alkalinization. The following pieces of evidence suggest that the gating process controlled by intracellular pH (pH(i)) is independent from that under the command of pH(o). It was not possible to overcome closure by extracellular acidification by means of intracellular alkalinization. The mutant TASK-2-R224A that lacks sensitivity to pH(o) had normal pH(i)-dependent gating. Increasing extracellular K(+) concentration acid shifts pH(o) activity curve of TASK-2 yet did not affect pH(i) gating of TASK-2. pH(o) modulation of TASK-2 is voltage-dependent, whereas pH(i) gating was not altered by membrane potential. These results suggest that pH(o), which controls a selectivity filter external gate, and pH(i) act at different gating processes to open and close TASK-2 channels. We speculate that pH(i) regulates an inner gate. We demonstrate that neutralization of a lysine residue (Lys(245)) located at the C-terminal end of transmembrane domain 4 by mutation to alanine abolishes gating by pH(i). We postulate that this lysine acts as an intracellular pH sensor as its mutation to histidine acid-shifts the pH(i)-dependence curve of TASK-2 as expected from its lower pK(a). We conclude that intracellular pH, together with pH(o), is a critical determinant of TASK-2 activity and therefore of its physiological function.

Characterization and functional restoration of a potassium channel Kir6.2 pore mutation identified in congenital hyperinsulinism

The inwardly rectifying potassium channel Kir6.2 assembles with sulfonylurea receptor 1 to form the ATP-sensitive potassium (K(ATP)) channels that regulate insulin secretion in pancreatic beta-cells. Mutations in K(ATP) channels underlie insulin secretion disease. Here, we report the characterization of a heterozygous missense Kir6.2 mutation, G156R, identified in congenital hyperinsulinism. Homomeric mutant channels reconstituted in COS cells show similar surface expression as wild-type channels but fail to conduct potassium currents. The mutated glycine is in the pore-lining transmembrane helix of Kir6.2; an equivalent glycine in other potassium channels has been proposed to serve as a hinge to allow helix bending during gating. We found that mutation of an adjacent asparagine, Asn-160, to aspartate, which converts the channel from a weak to a strong inward rectifier, on the G156R background restored ion conduction in the mutant channel. Unlike N160D channels, however, G156R/N160D channels are not blocked by intracellular polyamines at positive membrane potential and exhibit wild-type-like nucleotide sensitivities, suggesting the aspartate introduced at position 160 interacts with arginine at 156 to restore ion conduction and gating. Using tandem Kir6.2 tetramers containing G156R and/or N160D in designated positions, we show that one mutant subunit in the tetramer is insufficient to abolish conductance and that G156R and N160D can interact in the same or adjacent subunits to restore conduction. We conclude that the glycine at 156 is not essential for K(ATP) channel gating and that the Kir6.2 gating defect caused by the G156R mutation could be rescued by manipulating chemical interactions between pore residues.

Examining cooperative gating phenomena in voltage-dependent potassium channels: taking the energetic approach

Allosteric regulation of protein function is often achieved by changes in protein conformation induced by changes in chemical or electrical potential. In multisubunit proteins, such conformational changes may give rise to cooperativity in ligand binding. Conformational changes between open and closed states are central to the function of voltage-activated potassium (Kv) channel proteins, homotetrameric pore-forming membrane proteins involved in generating and shaping action potentials in excitable cells. Accessible to extremely high signal-to-noise ratio in functional measurements, combined with the availability of high-resolution structural data for different conformations of the protein, the Kv channel represents an excellent allosteric model system to further understand the aspects of synergism and cooperative effects in protein function. In this chapter, we demonstrate how the use of the simple law of mass action combined with thermodynamic mutant cycle energetic coupling analysis of Kv channel gating can be used to provide valuable information regarding (1) how cooperativity in Kv channel pore opening can be assessed; (2) how one can directly discriminate whether conformational transitions during Kv channel pore opening occur in a concerted or sequential manner; and (3) how mechanistically, the coupling between distant activation gate and selectivity filter functional elements of the prototypical Shaker Kv channel protein might be achieved. In addition to providing valuable insight into the function of this important protein, the conclusions reached at using high-order thermodynamic energetic coupling analysis applied to the Kv channel allosteric model system reveal much about the function of allosteric proteins, in general.

Dynamic aspects of functional regulation of the ATP receptor channel P2X2

The P2X(2) channel is a ligand-gated channel activated by ATP. Functional features that reflect the dynamic flexibility of the channel include time-dependent pore dilatation following ATP application and direct inhibitory interaction with activated nicotinic acetylcholine receptors on the membrane. We have been studying the mechanisms by which P2X(2) channel functionality is dynamically regulated. Using a Xenopus oocyte expression system, we observed that the pore properties, including ion selectivity and rectification, depend on the open channel density on the membrane. Pore dilatation was apparent when the open channel density was high and inward rectification was modest. We also observed that P2X(2) channels show voltage dependence, despite the absence of a canonical voltage sensor. At a semi-steady state after ATP application, P2X(2) channels were activated upon membrane hyperpolarization. This voltage-dependent activation was also [ATP] dependent. With increases in [ATP], the speed of hyperpolarization-induced activation was increased and the conductance-voltage relationship was shifted towards depolarized potentials. Based on analyses of experimental data and various simulations, we propose that these phenomena can be explained by assuming a fast ATP binding step and a rate-limiting voltage-dependent gating step. Complete elucidation of these regulatory mechanisms awaits dynamic imaging of functioning P2X(2) channels.

Intersubunit coupling in the pore of BK channels

The structural basis underlying the gating of large conductance Ca(2+)-activated K(+) (BK) channels remains elusive. We found that substitution of Leu-312 in the S6 transmembrane segment of mSlo1 BK channels with hydrophilic amino acids of smaller side-chain volume favored the open state. The sensitivities of channels to calcium and voltage were modified by some mutations and completely abolished by others. Interpretation of the results in terms of an allosteric model suggests that the calcium-insensitive mutants greatly destabilize the closed relative to the open conformation and may also disrupt the allosteric coupling between Ca(2+) or voltage sensors and the gate. Some Phe-315 mutations also favor the open state, suggesting that Leu-312 and Phe-315 may interact in the closed state, forming a major energy barrier that the channel has to overcome to open. Homology modeling and molecular dynamic simulations further support that the side chain of Leu-312 can couple strongly with the aromatic ring of Phe-315 in neighboring subunits (L-F coupling) to maintain the channel closed. Additionally, single-channel recordings indicate that the calcium-insensitive mutants, whose kinetics can be approximately characterized by a two-state closed-open (C-O) model, exhibit nearly 100% open probability under physiological conditions without alterations in single-channel conductance. These findings provide a basis for understanding the structure and gating of the BK channel pore.

Structural determinants of gating in the TRPV1 channel

Transient receptor potential vanilloid 1 (TRPV1) channels mediate several types of physiological responses. Despite the importance of these channels in pain detection and inflammation, little is known about how their structural components convert different types of stimuli into channel activity. To localize the activation gate of these channels, we inserted cysteines along the S6 segment of mutant TRPV1 channels and assessed their accessibility to thiol-modifying agents. We show that access to the pore of TRPV1 is gated by S6 in response to both capsaicin binding and increases in temperature, that the pore-forming S6 segments are helical structures and that two constrictions are present in the pore: one that impedes the access of large molecules and the other that hampers the access of smaller ions and constitutes an activation gate of these channels.

The evolutionarily conserved residue A653 plays a key role in HERG channel closing

Human ether-a-go-go-related gene (HERG) encodes the rapid, outwardly rectifying K(+) current I(Kr) that is critical for repolarization of the cardiac action potential. Congenital HERG mutations or unintended pharmaceutical block of I(Kr) can lead to life-threatening arrhythmias. Here, we assess the functional role of the alanine at position 653 (HERG-A653) that is highly conserved among evolutionarily divergent K(+) channels. HERG-A653 is close to the 'glycine hinge' implicated in K(+) channel opening, and is flanked by tyrosine 652 and phenylalanine 656, which contribute to the drug binding site. We substituted an array of seven (I, C, S, G, Y, V and T) amino acids at position 653 and expressed individual variants in heterologous systems to assess changes in gating and drug binding. Substitution of A653 resulted in negative shifts of the V(1/2) of activation ranging from -23.6 (A653S) to -62.5 (A653V) compared to -11.2 mV for wild-type (WT). Deactivation was also drastically altered: channels with A653I/C substitutions exhibited delayed deactivation in response to test potentials above the activation threshold, while A653S/G/Y/V/T failed to deactivate under those conditions and required hyperpolarization and prolonged holding potentials at -130 mV. While A653S/G/T/Y variants showed decreased sensitivity to the I(Kr) inhibitor dofetilide, these changes could not be correlated with defects in channel closure. Homology modelling suggests that in the closed state, A653 forms tight contacts with several residues from the neighbouring subunit in the tetramer, playing a key role in S6 helix packing at the narrowest part of the vestibule. Our study suggests that A653 plays an important functional role in the outwardly rectifying gating behaviour of HERG, supporting channel closure at membrane potentials negative to the channel activation threshold.

Gating the pore of potassium leak channels

A key feature of potassium channel function is the ability to switch between conducting and non-conducting states by undergoing conformational changes in response to cellular or extracellular signals. Such switching is facilitated by the mechanical coupling of gating domain movements to pore opening and closing. Two-pore domain potassium channels (K(2P)) conduct leak or background potassium-selective currents that are mostly time- and voltage-independent. These channels play a significant role in setting the cell resting membrane potential and, therefore modulate cell responsiveness and excitability. Thus, K(2P) channels are key players in numerous physiological processes and were recently shown to also be involved in human pathologies. It is well established that K(2P) channel conductance, open probability and cell surface expression are significantly modulated by various physical and chemical stimuli. However, in understanding how such signals are translated into conformational changes that open or close the channels gate, there remain more open questions than answers. A growing line of evidence suggests that the outer pore area assumes a critical role in gating K(2P) channels, in a manner reminiscent of C-type inactivation of voltage-gated potassium channels. In some K(2P) channels, this gating mechanism is facilitated in response to external pH levels. Recently, it was suggested that K(2P) channels also possess a lower activation gate that is positively coupled to the outer pore gate. The purpose of this review is to present an up-to-date summary of research describing the conformational changes and gating events that take place at the K(2P) channel ion-conducting pathway during the channel regulation.

Molecular mechanisms of BK channel activation

Large conductance, Ca(2+)-activated potassium (BK) channels are widely expressed throughout the animal kingdom and play important roles in many physiological processes, such as muscle contraction, neural transmission and hearing. These physiological roles derive from the ability of BK channels to be synergistically activated by membrane voltage, intracellular Ca(2+) and other ligands. Similar to voltage-gated K(+) channels, BK channels possess a pore-gate domain (S5-S6 transmembrane segments) and a voltage-sensor domain (S1-S4). In addition, BK channels contain a large cytoplasmic C-terminal domain that serves as the primary ligand sensor. The voltage sensor and the ligand sensor allosterically control K(+) flux through the pore-gate domain in response to various stimuli, thereby linking cellular metabolism and membrane excitability. This review summarizes the current understanding of these structural domains and their mutual interactions in voltage-, Ca(2+)- and Mg(2+)-dependent activation of the channel.

Voltage- and [ATP]-dependent gating of the P2X(2) ATP receptor channel

P2X receptors are ligand-gated cation channels activated by extracellular adenosine triphosphate (ATP). Nonetheless, P2X₂ channel currents observed during the steady-state after ATP application are known to exhibit voltage dependence; there is a gradual increase in the inward current upon hyperpolarization. We used a Xenopus oocyte expression system and two-electrode voltage clamp to analyze this “activation” phase quantitatively. We characterized the conductance–voltage relationship in the presence of various [ATP], and observed that it shifted toward more depolarized potentials with increases in [ATP]. By analyzing the rate constants for the channel's transition between a closed and an open state, we showed that the gating of P2X₂ is determined in a complex way that involves both membrane voltage and ATP binding. The activation phase was similarly recorded in HEK293 cells expressing P2X₂ even by inside-out patch clamp after intensive perfusion, excluding a possibility that the gating is due to block/unblock by endogenous blocker(s) of oocytes. We investigated its structural basis by substituting a glycine residue (G344) in the second transmembrane (TM) helix, which may provide a kink that could mediate “gating.” We found that, instead of a gradual increase, the inward current through the G344A mutant increased instantaneously upon hyperpolarization, whereas a G344P mutant retained an activation phase that was slower than the wild type (WT). Using glycine-scanning mutagenesis in the background of G344A, we could recover the activation phase by introducing a glycine residue into the middle of second TM. These results demonstrate that the flexibility of G344 contributes to the voltage-dependent gating. Finally, we assumed a three-state model consisting of a fast ATP-binding step and a following gating step and estimated the rate constants for the latter in P2X₂-WT. We then executed simulation analyses using the calculated rate constants and successfully reproduced the results observed experimentally, voltage-dependent activation that is accelerated by increases in [ATP].

Molecular template for a voltage sensor in a novel K+ channel. III. Functional reconstitution of a sensorless pore module from a prokaryotic Kv channel

KvLm is a prokaryotic voltage-gated K⁺ (Kv) channel from Listeria monocytogenes. The sequence of the voltage-sensing module (transmembrane segments S1-S4) of KvLm is atypical in that it contains only three of the eight conserved charged residues known to be deterministic for voltage sensing in eukaryotic Kv's. In contrast, the pore module (PM), including the S4-S5 linker and cytoplasmic tail (linker-S5-P-S6-C-terminus) of KvLm, is highly conserved. Here, the full-length (FL)-KvLm and the KvLm-PM only proteins were expressed, purified, and reconstituted into giant liposomes. The properties of the reconstituted FL-KvLm mirror well the characteristics of the heterologously expressed channel in Escherichia coli spheroplasts: a right-shifted voltage of activation, micromolar tetrabutylammonium-blocking affinity, and a single-channel conductance comparable to that of eukaryotic Kv's. Conversely, ionic currents through the PM recapitulate both the conductance and blocking properties of the FL-KvLm, yet the KvLm-PM exhibits only rudimentary voltage dependence. Given that the KvLm-PM displays many of the conduction properties of FL-KvLm and of other eukaryotic Kv's, including strict ion selectivity, we conclude that self-assembly of the PM subunits in lipid bilayers, in the absence of the voltage-sensing module, generates a conductive oligomer akin to that of the native KvLm, and that the structural independence of voltage sensing and PMs observed in eukaryotic Kv channels was initially implemented by nature in the design of prokaryotic Kv channels. Collectively, the results indicate that this robust functional module will prove valuable as a molecular template for coupling new sensors and to elucidate PM residue–specific contributions to Kv conduction properties.

Inverse coupling in leak and voltage-activated K+ channel gates underlies distinct roles in electrical signaling

Voltage-activated (Kv) and leak (K(2P)) K(+) channels have key, yet distinct, roles in electrical signaling in the nervous system. Here we examine how differences in the operation of the activation and slow inactivation pore gates of Kv and K(2P) channels underlie their unique roles in electrical signaling. We report that (i) leak K(+) channels possess a lower activation gate, (ii) the activation gate is an important determinant controlling the conformational stability of the K(+) channel pore, (iii) the lower activation and upper slow inactivation gates of leak channels cross-talk and (iv) unlike Kv channels, where the two gates are negatively coupled, these two gates are positively coupled in K(2P) channels. Our results demonstrate how basic thermodynamic properties of the K(+) channel pore, particularly conformational stability and coupling between gates, underlie the specialized roles of Kv and K(2P) channel families in electrical signaling.

High-resolution structure of the open NaK channel

We report the crystal structure of the non-selective cation channel NaK from b. cereus at a resolution of 1.6 Å. The structure reveals the intracellular gate in an open state compared to the closed form reported previously, making NaK the only channel for which the three-dimensional structures of both conformations are known. Channel opening follows a conserved mechanism of inner helix bending utilizing a flexible glycine residue, the gating hinge, seen in MthK and most other tetrameric cation channels. Additionally, distinct inter and intra-subunit rearrangements involved in channel gating are seen and characterized for the first time along with inner helix twisting motions. Furthermore, we identify a residue deeper within the cavity of the channel pore, Phe92, which likely forms a constriction point within the open pore, restricting ion flux through the channel. Mutating this residue to Ala causes a subsequent increase in ion conduction rates as measured by ⁸⁶Rb flux assays. The structures of both the open and closed conformations of the NaK channel correlate very well with those of equivalent K⁺ channel conformations, namely MthK and KcsA, respectively.

Direct analysis of cooperativity in multisubunit allosteric proteins

Allosteric regulation of protein function is a fundamental phenomenon of major importance in many cellular processes. Such regulation is often achieved by ligand-induced conformational changes in multimeric proteins that may give rise to cooperativity in protein function. At the heart of allosteric mechanisms offered to account for such phenomenon, involving either concerted or sequential conformational transitions, lie changes in intersubunit interactions along the ligation pathway of the protein. However, structure-function analysis of such homooligomeric proteins by means of mutagenesis, although it provides valuable indirect information regarding (allosteric) mechanisms of action, it does not define the contribution of individual subunits nor interactions thereof to cooperativity in protein function, because any point mutation introduced into homooligomeric proteins will be present in all subunits. Here, we present a general strategy for the direct analysis of cooperativity in multisubunit proteins that combines measurement of the effects on protein function of all possible combinations of mutated subunits with analysis of the hierarchy of intersubunit interactions, assessed by using high-order double-mutant cycle-coupling analysis. We show that the pattern of high-order intersubunit coupling can serve as a discriminative criterion for defining concerted versus sequential conformational transitions underlying protein function. This strategy was applied to the particular case of the voltage-activated potassium channel protein (Kv) to provide compelling evidence for a concerted all-or-none activation gate opening of the Kv channel pore domain. An direct and detailed analysis of the contribution of high-order intersubunit interactions to cooperativity in the function of an allosteric protein has not previously been presented.