Pore dimensions and the role of occupancy in unitary conductance of Shaker K channels

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels

The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.

Sweetening K-channels: what sugar taught us about permeation and gating

Because they enable for the modification of both viscosity and osmolarity, sugars have been used as a biophysical probe of voltage-gated K-channels for a while. Viscosity variations made it possible to measure the pore sizes in large and small conductance K-channels using techniques similar to those used in the 1980s to study the gramicidin A channel. These analyses led to the finding that the size of the internal mouth appears to be the primary cause of the conductance differences between Shaker-like channels and large conductance BK-channels. As an osmotic agent, adding sugar unilaterally causes streaming potentials that indicate H₂O/K⁺ cotransport across the BK-channel pore. Osmotic experiments on Shaker K-channels suggest that the pore gate operation and the slow inactivation displace comparable amounts of water. Functionally isolated voltage sensors allow estimation of individual osmotic work for each voltage sensing charge during voltage-activation, reporting dramatic internal and external remodeling of the Voltage Sensing Domain´s solvent exposed surfaces. Remarkably, each charge of the VSD appears to take a unique trajectory. Thus, manipulation of viscosity and osmolarity, together with 3D structures, brings in solid grounds to harmonize function and structure in membrane proteins such as K-channels and, in a wider scope, other structurally dynamic proteins.

Modulation of the IKS channel by PIP2 requires two binding sites per monomer

The phosphatidyl-inositol-4,5-bisphosphate (PIP₂) lipid has been shown to be crucial for the coupling between the voltage sensor and the pore of the potassium voltage-gated K_V7 channel family, especially the K_V7.1 channel. Expressed in the myocardium membrane, K_V7.1 forms a complex with KCNE1 auxiliary subunits to generate the I_KS current. Here we present molecular models of the transmembrane region of this complex in its three known states, namely the Resting/Closed (RC), the Intermediate/Closed (IC), and the Activated/Open (AO), robustness of which is assessed by agreement with a range of biophysical data. Molecular Dynamics (MD) simulations of these models embedded in a lipid bilayer including phosphatidyl-inositol-4,5-bisphosphate (PIP₂) lipids show that in presence of KCNE1, two PIP₂ lipids are necessary to stabilize each state. The simulations also show that KCNE1 interacts with both PIP₂ binding sites, forming a tourniquet around the pore and preventing its opening. The present investigation provides therefore key molecular elements that govern the role of PIP₂ in KCNE1 modulation of I_KS channels, possibly a common mechanism by which auxiliary KCNE subunits might modulate a variety of other ion channels.

Effects of ionic strength on gating and permeation of TREK-2 K2P channels

In addition to the classical voltage-dependent behavior mediated by the voltage-sensing-domains (VSD) of ion channels, a growing number of voltage-dependent gating behaviors are being described in channels that lack canonical VSDs. A common thread in their mechanism of action is the contribution of the permeating ion to this voltage sensing process. The polymodal K2P K⁺ channel, TREK2 responds to membrane voltage through a gating process mediated by the interaction of K⁺ with its selectivity filter. Recently, we found that this action can be modulated by small molecule agonists (e.g. BL1249) which appear to have an electrostatic influence on K⁺ binding within the inner cavity and produce an increase in the single-channel conductance of TREK-2 channels. Here, we directly probed this K⁺-dependent gating process by recording both macroscopic and single-channel currents of TREK-2 in the presence of high concentrations of internal K⁺. Surprisingly we found TREK-2 is inhibited by high internal K⁺ concentrations and that this is mediated by the concomitant increase in ionic-strength. However, we were still able to determine that the increase in single channel conductance in the presence of BL1249 was blunted in high ionic-strength, whilst its activatory effect (on channel open probability) persisted. These effects are consistent with an electrostatic mechanism of action of negatively charged activators such as BL1249 on permeation, but also suggest that their influence on channel gating is complex.

Physics of Selective Conduction and Point Mutation in Biological Ion Channels

We introduce a statistical and linear response theory of selective conduction in biological ion channels with multiple binding sites and possible point mutation. We derive an effective grand-canonical ensemble and generalized Einstein relations for the selectivity filter, assuming strongly coordinated ionic motion, and allowing for ionic Coulomb blockade. The theory agrees well with data from the KcsA K^{+} channel and a mutant. We show that the Eisenman relations for thermodynamic selectivity follow from the condition for fast conduction and find that maximum conduction requires the binding sites to be nearly identical.

Atomistic basis of opening and conduction in mammalian inward rectifier potassium (Kir2.2) channels

This paper presents the crystal structure of a forced open inward rectifier Kir2.2 channel. Molecular dynamics reveals the details of ion permeation through the open channel.

Potassium ion conduction through open potassium channels is essential to control of membrane potentials in all cells. To elucidate the open conformation and hence the mechanism of K⁺ ion conduction in the classic inward rectifier Kir2.2, we introduced a negative charge (G178D) at the crossing point of the inner helix bundle, the location of ligand-dependent gating. This “forced open” mutation generated channels that were active even in the complete absence of phosphatidylinositol-4,5-bisphosphate (PIP₂), an otherwise essential ligand for Kir channel opening. Crystal structures were obtained at a resolution of 3.6 Å without PIP₂ bound, or 2.8 Å in complex with PIP₂. The latter revealed a slight widening at the helix bundle crossing (HBC) through backbone movement. MD simulations showed that subsequent spontaneous wetting of the pore through the HBC gate region allowed K⁺ ion movement across the HBC and conduction through the channel. Further simulations reveal atomistic details of the opening process and highlight the role of pore-lining acidic residues in K⁺ conduction through Kir2 channels.

TOK channels use the two gates in classical K+ channels to achieve outward rectification

TOKs are outwardly rectifying K⁺ channels in fungi with two pore-loops and eight transmembrane spans. Here, we describe the TOKs from four pathogens that cause the majority of life-threatening fungal infections in humans. These TOKs pass large currents only in the outward direction like the canonical isolate from Saccharomyces cerevisiae (ScTOK), and distinct from other K⁺ channels. ScTOK, AfTOK1 (Aspergillus fumigatus), and H99TOK (Cryptococcus neoformans grubii) are K⁺ -selective and pass current above the K⁺ reversal potential. CaTOK (Candida albicans) and CnTOK (Cryptococcus neoformans neoformans) pass both K⁺ and Na⁺ and conduct above a reversal potential reflecting the mixed permeability of their selectivity filter. Mutations in CaTOK and ScTOK at sites homologous to those that open the internal gates in classical K⁺ channels are shown to produce inward TOK currents. A favored model for outward rectification is proposed whereby the reversal potential determines ion occupancy, and thus, conductivity, of the selectivity filter gate that is coupled to an imperfectly restrictive internal gate, permitting the filter to sample ion concentrations on both sides of the membrane.

The Selectivity Filter Is Involved in the U-Type Inactivation Process of Kv2.1 and Kv3.1 Channels

Voltage-gated potassium (Kv) channels display several types of inactivation processes, including N-, C-, and U-types. C-type inactivation is attributed to a nonconductive conformation of the selectivity filter (SF). It has been proposed that the activation gate and the channel's SF are allosterically coupled because the conformational changes of the former affect the structure of the latter and vice versa. The second threonine of the SF signature sequence (e.g., TTVGYG) has been proven to be essential for this allosteric coupling. To further study the role of the SF in U-type inactivation, we substituted the second threonine of the TTVGYG sequence by an alanine in the hKv2.1 and hKv3.1 channels, which are known to display U-type inactivation. Both hKv2.1-T377A and hKv3.1-T400A yielded channels that were resistant to inactivation, and as a result, they displayed noninactivating currents upon channel opening; i.e., hKv2.1-T377A and hKv3.1-T400A remained fully conductive upon prolonged moderate depolarizations, whereas in wild-type hKv2.1 and hKv3.1, the current amplitude typically reduces because of U-type inactivation. Interestingly, increasing the extracellular K⁺ concentration increased the macroscopic current amplitude of both hKv2.1-T377A and hKv3.1-T400A, which is similar to the response of the homologous T to A mutation in Shaker and hKv1.5 channels that display C-type inactivation. Our data support an important role for the second threonine of the SF signature sequence in the U-type inactivation gating of hKv2.1 and hKv3.1.

Trans-toxin ion-sensitivity of charybdotoxin-blocked potassium-channels reveals unbinding transitional states

In silico and in vitro studies have made progress in understanding protein–protein complex formation; however, the molecular mechanisms for their dissociation are unclear. Protein–protein complexes, lasting from microseconds to years, often involve induced-fit, challenging computational or kinetic analysis. Charybdotoxin (CTX), a peptide from the Leiurus scorpion venom, blocks voltage-gated K⁺-channels in a unique example of binding/unbinding simplicity. CTX plugs the external mouth of K⁺-channels pore, stopping K⁺-ion conduction, without inducing conformational changes. Conflicting with a tight binding, we show that external permeant ions enhance CTX-dissociation, implying a path connecting the pore, in the toxin-bound channel, with the external solution. This sensitivity is explained if CTX wobbles between several bound conformations, producing transient events that restore the electrical and ionic trans-pore gradients. Wobbling may originate from a network of contacts in the interaction interface that are in dynamic stochastic equilibria. These partially-bound intermediates could lead to distinct, and potentially manipulable, dissociation pathways.

Gating-induced large aqueous volumetric remodeling and aspartate tolerance in the voltage sensor domain of Shaker K+ channels

Neurons encode electrical signals with critically tuned voltage-gated ion channels and enzymes. Dedicated voltage sensor domains (VSDs) in these membrane proteins activate coordinately with an unresolved structural change. Such change conveys the transmembrane translocation of four positively charged arginine side chains, the voltage-sensing residues (VSRs; R1-R4). Countercharges and lipid phosphohead groups likely stabilize these VSRs within the low-dielectric core of the protein. However, the role of hydration, a sign-independent charge stabilizer, remains unclear. We replaced all VSRs and their neighboring residues with negatively charged aspartates in a voltage-gated potassium channel. The ensuing mild functional effects indicate that hydration is also important in VSR stabilization. The voltage dependency of the VSR aspartate variants approached the expected arithmetic summation of charges at VSR positions, as if negative and positive side chains faced similar pathways. In contrast, aspartates introduced between R2 and R3 did not affect voltage dependence as if the side chains moved outside the electric field or together with it, undergoing a large displacement and volumetric remodeling. Accordingly, VSR performed osmotic work at both internal and external aqueous interfaces. Individual VSR contributions to volumetric works approached arithmetical additivity but were largely dissimilar. While R1 and R4 displaced small volumes, R2 and R3 volumetric works were massive and vectorially opposed, favoring large aqueous remodeling during VSD activation. These diverse volumetric works are, at least for R2 and R3, not compatible with VSR translocation across a unique stationary charge transfer center. Instead, VSRs may follow separated pathways across a fluctuating low-dielectric septum.

Molecular Determinants of BK Channel Functional Diversity and Functioning

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

IBiSA_Tools: A Computational Toolkit for Ion-Binding State Analysis in Molecular Dynamics Trajectories of Ion Channels

Ion conduction mechanisms of ion channels are a long-standing conundrum. Although the molecular dynamics (MD) method has been extensively used to simulate ion conduction dynamics at the atomic level, analysis and interpretation of MD results are not straightforward due to complexity of the dynamics. In our previous reports, we proposed an analytical method called ion-binding state analysis to scrutinize and summarize ion conduction mechanisms by taking advantage of a variety of analytical protocols, e.g., the complex network analysis, sequence alignment, and hierarchical clustering. This approach effectively revealed the ion conduction mechanisms and their dependence on the conditions, i.e., ion concentration and membrane voltage. Here, we present an easy-to-use computational toolkit for ion-binding state analysis, called IBiSA_tools. This toolkit consists of a C++ program and a series of Python and R scripts. From the trajectory file of MD simulations and a structure file, users can generate several images and statistics of ion conduction processes. A complex network named ion-binding state graph is generated in a standard graph format (graph modeling language; GML), which can be visualized by standard network analyzers such as Cytoscape. As a tutorial, a trajectory of a 50 ns MD simulation of the Kv1.2 channel is also distributed with the toolkit. Users can trace the entire process of ion-binding state analysis step by step. The novel method for analysis of ion conduction mechanisms of ion channels can be easily used by means of IBiSA_tools. This software is distributed under an open source license at the following URL: http://www.ritsumei.ac.jp/~ktkshr/ibisa_tools/

Pore size matters for potassium channel conductance

Ion channels are membrane proteins that mediate efficient ion transport across the hydrophobic core of cell membranes, an unlikely process in their absence. K⁺ channels discriminate K⁺ over cations with similar radii with extraordinary selectivity and display a wide diversity of ion transport rates, covering differences of two orders of magnitude in unitary conductance. The pore domains of large- and small-conductance K⁺ channels share a general architectural design comprising a conserved narrow selectivity filter, which forms intimate interactions with permeant ions, flanked by two wider vestibules toward the internal and external openings. In large-conductance K⁺ channels, the inner vestibule is wide, whereas in small-conductance channels it is narrow. Here we raise the idea that the physical dimensions of the hydrophobic internal vestibule limit ion transport in K⁺ channels, accounting for their diversity in unitary conductance.

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

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

Biophysics of BK Channel Gating

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

Ion Concentration- and Voltage-Dependent Push and Pull Mechanisms of Potassium Channel Ion Conduction

The mechanism of ion conduction by potassium channels is one of the central issues in physiology. In particular, it is still unclear how the ion concentration and the membrane voltage drive ion conduction. We have investigated the dynamics of the ion conduction processes in the Kv1.2 pore domain, by molecular dynamics (MD) simulations with several different voltages and ion concentrations. By focusing on the detailed ion movements through the pore including selectivity filter (SF) and cavity, we found two major conduction mechanisms, called the III-IV-III and III-II-III mechanisms, and the balance between the ion concentration and the voltage determines the mechanism preference. In the III-IV-III mechanism, the outermost ion in the pore is pushed out by a new ion coming from the intracellular fluid, and four-ion states were transiently observed. In the III-II-III mechanism, the outermost ion is pulled out first, without pushing by incoming ions. Increases in the ion concentration and voltage accelerated ion conductions, but their mechanisms were different. The increase in the ion concentrations facilitated the III-IV-III conductions, while the higher voltages increased the III-II-III conductions, indicating that the pore domain of potassium channels permeates ions by using two different driving forces: a push by intracellular ions and a pull by voltage.

Effective pore size and radius of capture for K(+) ions in K-channels

Reconciling protein functional data with crystal structure is arduous because rare conformations or crystallization artifacts occur. Here we present a tool to validate the dimensions of open pore structures of potassium-selective ion channels. We used freely available algorithms to calculate the molecular contour of the pore to determine the effective internal pore radius (r_E) in several K-channel crystal structures. r_E was operationally defined as the radius of the biggest sphere able to enter the pore from the cytosolic side. We obtained consistent r_E estimates for MthK and Kv1.2/2.1 structures, with r_E = 5.3–5.9 Å and r_E = 4.5–5.2 Å, respectively. We compared these structural estimates with functional assessments of the internal mouth radii of capture (r_C) for two electrophysiological counterparts, the large conductance calcium activated K-channel (r_C = 2.2 Å) and the Shaker Kv-channel (r_C = 0.8 Å), for MthK and Kv1.2/2.1 structures, respectively. Calculating the difference between r_E and r_C, produced consistent size radii of 3.1–3.7 Å and 3.6–4.4 Å for hydrated K⁺ ions. These hydrated K⁺ estimates harmonize with others obtained with diverse experimental and theoretical methods. Thus, these findings validate MthK and the Kv1.2/2.1 structures as templates for open BK and Kv-channels, respectively.