Activation of Shaker potassium channels. III. An activation gating model for wild-type and V2 mutant channels

How much does TRPV1 deviate from an ideal MWC-type protein?

Many ion channels are known to behave as an allosteric protein, coupling environmental stimuli captured by specialized sensing domains to the opening of a central pore. The classic Monod-Wyman-Changeux (MWC) model, originally proposed to describe binding of gas molecules to hemoglobin, has been widely used as a framework for analyzing ion channel gating. Here, we address the issue of how accurately the MWC model predicts activation of the capsaicin receptor TRPV1 by vanilloids. Taking advantage of a concatemeric design that makes it possible to lock TRPV1 in states with zero to four bound vanilloid molecules, we showed quantitatively that the overall gating behavior is satisfactorily predicted by the MWC model. There is, however, a small yet detectable subunit position effect: ligand binding to two kitty-corner subunits is 0.3-0.4 kcal/mol more effective in inducing opening than binding to two neighbor subunits. This difference-less than 10% of the overall energetic contribution from ligand binding-might be due to the restriction on subunit arrangement imposed by the planar membrane; if this is the case, then the position effect is not expected in hemoglobin, in which each subunit is related equivalently to all the other subunits.

Structural change of the cytoplasmic N-terminus and S1 segment of voltage-sensing phosphatase reported by Anap

BACKGROUND

Voltage-sensing phosphatase contains a structurally conserved S1-S4-based voltage-sensor domain, which undergoes a conformational transition in response to membrane potential change. Unlike that of channels, it is functional even in isolation and is therefore advantageous for studying the transition mechanism, but its nature has not yet been fully elucidated. This study aimed to address whether the cytoplasmic N-terminus and S1 exhibit structural change.

METHODS

Anap, an environment-sensitive unnatural fluorescent amino acid, was site-specifically introduced to the voltage sensor domain to probe local structural changes by using oocyte voltage clamp and photometry. Tetramethylrhodamine was also used to probe some extracellularly accessible positions. In total, 51 positions were investigated.

RESULTS

We detected robust voltage-dependent signals from widely distributed positions including N-terminus and S1. In addition, response to hyperpolarization was observed at the extracellular end of S1, reflecting the local structure flexibility of the voltage-sensor domain in the down-state. We also found that the mechanical coupling between the voltage-sensor and phosphatase domains affects the depolarization-induced optical signals but not the hyperpolarization-induced signals.

CONCLUSIONS

These results fill a gap between the previous interpretations from the structural and biophysical approaches and should provide important insights into the mechanisms of the voltage-sensor domain transition as well as its coupling with the effector.

Evaluating sequential and allosteric activation models in IKs channels with mutated voltage sensors

The ion-conducting IKs channel complex, important in cardiac repolarization and arrhythmias, comprises tetramers of KCNQ1 α-subunits along with 1-4 KCNE1 accessory subunits and calmodulin regulatory molecules. The E160R mutation in individual KCNQ1 subunits was used to prevent activation of voltage sensors and allow direct determination of transition rate data from complexes opening with a fixed number of 1, 2, or 4 activatable voltage sensors. Markov models were used to test the suitability of sequential versus allosteric models of IKs activation by comparing simulations with experimental steady-state and transient activation kinetics, voltage-sensor fluorescence from channels with two or four activatable domains, and limiting slope currents at negative potentials. Sequential Hodgkin-Huxley-type models approximately describe IKs currents but cannot explain an activation delay in channels with only one activatable subunit or the hyperpolarizing shift in the conductance-voltage relationship with more activatable voltage sensors. Incorporating two voltage sensor activation steps in sequential models and a concerted step in opening via rates derived from fluorescence measurements improves models but does not resolve fundamental differences with experimental data. Limiting slope current data that show the opening of channels at negative potentials and very low open probability are better simulated using allosteric models of activation with one transition per voltage sensor, which implies that movement of all four sensors is not required for IKs conductance. Tiered allosteric models with two activating transitions per voltage sensor can fully account for IKs current and fluorescence activation kinetics in constructs with different numbers of activatable voltage sensors.

How Much Does TRPV1 Deviate from An Ideal MWC-Type Protein

Many ion channels are known to behave as an allosteric protein, coupling environmental stimuli captured by specialized sensing domains to the opening of a central pore. The classic Monod-Wyman-Changeux (MWC) model, originally proposed to describe binding of gas molecules to hemoglobin, has been widely used for analyzing ion channel gating. Here we address the issue of how accurate the MWC model predicts activation of the capsaicin receptor TRPV1 by vanilloids. Taking advantage of a concatemeric design that makes it possible to lock TRPV1 in states with zero-to-four bound vanilloid molecules, we showed quantitatively that the overall gating behavior is satisfactorily predicted by the MWC model. There is however a small yet detectable subunit position effect: ligand binding to two kitty-corner subunits is 0.4-to-0.6 kcal/mol more effective in inducing opening than binding to two neighbor subunits. This difference, less than 10% of the overall energetic contribution from ligand binding, is not expected in hemoglobin, in which each subunit is related equivalently to all the other subunits.

Significance

The MWC model, proposed more than 50 years ago, is elegantly simple yet powerful in predicting the behavior of allosteric proteins like hemoglobin. Its prediction power for ion channel gating has been beautifully demonstrated in the studies of BK channels. Our present work aims to determine how accurate the MWC model predicts TRPV1 activation induced by vanilloids. Our findings support the notion that the evolutionary drive upon allosteric proteins applies generally to multi-subunit proteins.

Fifty years of gating currents and channel gating

We celebrate this year the 50th anniversary of the first electrophysiological recordings of the gating currents from voltage-dependent ion channels done in 1973. This retrospective tries to illustrate the context knowledge on channel gating and the impact gating-current recording had then, and how it continued to clarify concepts, elaborate new ideas, and steer the scientific debate in these 50 years. The notion of gating particles and gating currents was first put forward by Hodgkin and Huxley in 1952 as a necessary assumption for interpreting the voltage dependence of the Na and K conductances of the action potential. 20 years later, gating currents were actually recorded, and over the following decades have represented the most direct means of tracing the movement of the gating charges and gaining insights into the mechanisms of channel gating. Most work in the early years was focused on the gating currents from the Na and K channels as found in the squid giant axon. With channel cloning and expression on heterologous systems, other channels as well as voltage-dependent enzymes were investigated. Other approaches were also introduced (cysteine mutagenesis and labeling, site-directed fluorometry, cryo-EM crystallography, and molecular dynamics [MD] modeling) to provide an integrated and coherent view of voltage-dependent gating in biological macromolecules. The layout of this retrospective reflects the past 50 years of investigations on gating currents, first addressing studies done on Na and K channels and then on other voltage-gated channels and non-channel structures. The review closes with a brief overview of how the gating-charge/voltage-sensor movements are translated into pore opening and the pathologies associated with mutations targeting the structures involved with the gating currents.

A fast Markovian method for modeling channel noise in neurons

Channel noise results from rapid transitions of protein channels from closed to open state and is generally considered as the most dominant source of electrical noise causing membrane-potential fluctuations even in the absence of synaptic inputs. The simulation of a realistic channel noise remains a source of possible error. Although the Markovian method is considered as the golden standard for appropriate description of channel noise, its computation time increasing exponentially with the number of channels, it is poorly suitable to simulate realistic features. We describe here a novel algorithm at discrete time unit for simulating ion channel noise based on Markov chains (MC). Although this new algorithm refers to a Monte-Carlo process, it only needs few random numbers whatever the number of channels involved. Our fast MC (FMC) model does not exhibit the drawbacks due to approximations based on stochastic differential equations and the values of spike jitter are comparable to those obtained with the true Markovian method. In fact, we show here, that these drawbacks can be highlighted in the approximation based on stochastic differential equation methods even for a high number of channels (standard deviation of the 5th spike is about two-fold larger than that of MCF or true Markovian method for 5000 sodium channels). The FMC model appears therefore as the most accurate method to simulate channel noise with a fast execution time that does not depend on the channel number.

Phosphoinositide regulates dynamic movement of the S4 voltage sensor in the second repeat in two-pore channel 3

The two-pore channels (TPCs) are voltage-gated cation channels consisting of single polypeptides with two repeats of a canonical 6-transmembrane unit. TPCs are known to be regulated by various physiological signals such as membrane voltage and phosphoinositide (PI). The fourth helix in the second repeat (second S4) plays a major role in detecting membrane voltage, whereas the first repeat contains a PI binding site. Therefore, each of these stimuli is detected by a unique repeat to regulate the gating of the TPC central pore. How these various stimuli regulate the dynamic structural rearrangement of the TPC molecule remain unknown. Here, we found that PI binding to the first repeat in TPC3 regulates the movement of the distally located second S4 helix, showing that the PI-binding signal is not confined to the pore gate but also transmitted to the voltage sensor. Using voltage clamp fluorometry, measurement of gating charges, and Cys-accessibility analysis, we observed that PI binding significantly potentiates the voltage dependence of the movement of the second S4 helix. Notably, voltage clamp fluorometry analysis revealed that the voltage-dependent movement of the second S4 helix occurred in two phases, of which the second phase corresponds to the transfer of the gating charges. This movement was observed in the voltage range where gate-opening occurs and was potentiated by PI. In conclusion, this regulation of the second S4 helix by PI indicates a tight inter-repeat coupling within TPC3, a feature which might be conserved among TPC family members to integrate various physiological signals.

Mapping Electromechanical Coupling Pathways in Voltage-Gated Ion Channels: Challenges and the Way Forward

Inter- and intra-molecular allosteric interactions underpin regulation of activity in a variety of biological macromolecules. In the voltage-gated ion channel superfamily, the conformational state of the voltage-sensing domain regulates the activity of the pore domain via such long-range allosteric interactions. Although the overall structure of these channels is conserved, allosteric interactions between voltage-sensor and pore varies quite dramatically between the members of this superfamily. Despite the progress in identifying key residues and structural interfaces involved in mediating electromechanical coupling, our understanding of the biophysical mechanisms remains limited. Emerging new structures of voltage-gated ion channels in various conformational states will provide a better three-dimensional view of the process but to conclusively establish a mechanism, we will also need to quantitate the energetic contribution of various structural elements to this process. Using rigorous unbiased metrics, we want to compare the efficiency of electromechanical coupling between various sub-families in order to gain a comprehensive understanding. Furthermore, quantitative understanding of the process will enable us to correctly parameterize computational approaches which will ultimately enable us to predict allosteric activation mechanisms from structures. In this review, we will outline the challenges and limitations of various experimental approaches to measure electromechanical coupling and highlight the best practices in the field.

Resin-acid derivatives bind to multiple sites on the voltage-sensor domain of the Shaker potassium channel

Voltage-gated potassium (KV) channels can be opened by negatively charged resin acids and their derivatives. These resin acids have been proposed to attract the positively charged voltage-sensor helix (S4) toward the extracellular side of the membrane by binding to a pocket located between the lipid-facing extracellular ends of the transmembrane segments S3 and S4. By contrast to this proposed mechanism, neutralization of the top gating charge of the Shaker KV channel increased resin-acid-induced opening, suggesting other mechanisms and sites of action. Here, we explore the binding of two resin-acid derivatives, Wu50 and Wu161, to the activated/open state of the Shaker KV channel by a combination of in silico docking, molecular dynamics simulations, and electrophysiology of mutated channels. We identified three potential resin-acid-binding sites around S4: (1) the S3/S4 site previously suggested, in which positively charged residues introduced at the top of S4 are critical to keep the compound bound, (2) a site in the cleft between S4 and the pore domain (S4/pore site), in which a tryptophan at the top of S6 and the top gating charge of S4 keeps the compound bound, and (3) a site located on the extracellular side of the voltage-sensor domain, in a cleft formed by S1-S4 (the top-VSD site). The multiple binding sites around S4 and the anticipated helical-screw motion of the helix during activation make the effect of resin-acid derivatives on channel function intricate. The propensity of a specific resin acid to activate and open a voltage-gated channel likely depends on its exact binding dynamics and the types of interactions it can form with the protein in a state-specific manner.

Na+ and K+ channels: history and structure

In this perspective, we discuss the physiological roles of Na and K channels, emphasizing the importance of the K channel for cellular homeostasis in animal cells and of Na and K channels for cellular signaling. We consider the structural basis of Na and K channel gating in light of recent structural and electrophysiological findings.

High-Resolution Structures of K+ Channels

Potassium channels are present in every living cell and essential to setting up a stable, non-zero transmembrane electrostatic potential which manifests the off-equilibrium livelihood of the cell. They are involved in other cellular activities and regulation, such as the controlled release of hormones, the activation of T-cells for immune response, the firing of action potential in muscle cells and neurons, etc. Pharmacological reagents targeting potassium channels are important for treating various human diseases linked to dysfunction of the channels. High-resolution structures of these channels are very useful tools for delineating the detailed chemical basis underlying channel functions and for structure-based design and optimization of their pharmacological and pharmaceutical agents. Structural studies of potassium channels have revolutionized biophysical understandings of key concepts in the field - ion selectivity, conduction, channel gating, and modulation, making them multi-modality targets of pharmacological regulation. In this chapter, I will select a few high-resolution structures to illustrate key structural insights, proposed allostery behind channel functions, disagreements still open to debate, and channel-lipid interactions and co-evolution. The known structural consensus allows the inference of conserved molecular mechanisms shared among subfamilies of K⁺ channels and makes it possible to develop channel-specific pharmaceutical agents.

Gating and Regulation of KCNQ1 and KCNQ1 + KCNE1 Channel Complexes

The IKs channel complex is formed by the co-assembly of Kv7.1 (KCNQ1), a voltage-gated potassium channel, with its β-subunit, KCNE1 and the association of numerous accessory regulatory molecules such as PIP2, calmodulin, and yotiao. As a result, the IKs potassium current shows kinetic and regulatory flexibility, which not only allows IKs to fulfill physiological roles as disparate as cardiac repolarization and the maintenance of endolymph K⁺ homeostasis, but also to cause significant disease when it malfunctions. Here, we review new areas of understanding in the assembly, kinetics of activation and inactivation, voltage-sensor pore coupling, unitary events and regulation of this important ion channel complex, all of which have been given further impetus by the recent solution of cryo-EM structural representations of KCNQ1 alone and KCNQ1+KCNE3. Recently, the stoichiometric ratio of KCNE1 to KCNQ1 subunits has been confirmed to be variable up to a ratio of 4:4, rather than fixed at 2:4, and we will review the results and new methodologies that support this conclusion. Significant advances have been made in understanding differences between KCNQ1 and IKs gating using voltage clamp fluorimetry and mutational analysis to illuminate voltage sensor activation and inactivation, and the relationship between voltage sensor translation and pore domain opening. We now understand that the KCNQ1 pore can open with different permeabilities and conductance when the voltage sensor is in partially or fully activated positions, and the ability to make robust single channel recordings from IKs channels has also revealed the complicated pore subconductance architecture during these opening steps, during inactivation, and regulation by 1−4 associated KCNE1 subunits. Experiments placing mutations into individual voltage sensors to drastically change voltage dependence or prevent their movement altogether have demonstrated that the activation of KCNQ1 alone and IKs can best be explained using allosteric models of channel gating. Finally, we discuss how the intrinsic gating properties of KCNQ1 and IKs are highly modulated through the impact of intracellular signaling molecules and co-factors such as PIP2, protein kinase A, calmodulin and ATP, all of which modulate IKs current kinetics and contribute to diverse IKs channel complex function.

TRP ion channels: Proteins with conformational flexibility

Ion channels display conformational changes in response to binding of their agonists and antagonists. The study of the relationships between the structure and the function of these proteins has witnessed considerable advances in the last two decades using a combination of techniques, which include electrophysiology, optical approaches (i.e. patch clamp fluorometry, incorporation of non-canonic amino acids, etc.), molecular biology (mutations in different regions of ion channels to determine their role in function) and those that have permitted the resolution of their structures in detail (X-ray crystallography and cryo-electron microscopy). The possibility of making correlations among structural components and functional traits in ion channels has allowed for more refined conclusions on how these proteins work at the molecular level. With the cloning and description of the family of Transient Receptor Potential (TRP) channels, our understanding of several sensory-related processes has also greatly moved forward. The response of these proteins to several agonists, their regulation by signaling pathways as well as by protein-protein and lipid-protein interactions and, in some cases, their biophysical characteristics have been studied thoroughly and, recently, with the resolution of their structures, the field has experienced a new boom. This review article focuses on the conformational changes in the pores, concentrating on some members of the TRP family of ion channels (TRPV and TRPA subfamilies) that result in changes in their single-channel conductances, a phenomenon that may lead to fine-tuning the electrical response to a given agonist in a cell.

Simulation of Gating Currents of the Shaker K Channel Using a Brownian Model of the Voltage Sensor

The physical mechanism underlying the voltage-dependent gating of K channels is usually addressed theoretically using molecular dynamics simulations. However, besides being computationally very expensive, this approach is presently unable to fully predict the behavior of fundamental variables of channel gating such as the macroscopic gating current, and hence, it is presently unable to validate the model. To fill this gap, here we propose a voltage-gating model that treats the S4 segment as a Brownian particle moving through a gating channel pore and adjacent internal and external vestibules. In our model, charges on the S4 segment are screened by charged residues localized on neighboring segments of the channel protein and by ions present in the vestibules, whose dynamics are assessed using a flux conservation equation. The electrostatic voltage spatial profile is consistently assessed by applying the Poisson equation to all the charges present in the system. The treatment of the S4 segment as a Brownian particle allows description of the dynamics of a single S4 segment using the Langevin stochastic differential equation or the behavior of a population of S4 segments-useful for assessing the macroscopic gating current-using the Fokker-Planck equation. The proposed model confirms the gating charge transfer hypothesis with the movement of the S4 segment among five different stable positions where the gating charges interact in succession with the negatively charged residues on the channel protein. This behavior produces macroscopic gating currents quite similar to those experimentally found.

The dynamics of K+ channel gates as a biological transistor

Potassium channels are pore-forming membrane proteins that open and close in response to changes in a chemical or electrical potential, thereby regulating the flow of potassium ions across biological membranes. Two regions of the same channels are acting in tandem and enable ion flow through the channel pore. I refer to this coupled action as a "gate linker". To closely examine the role of the gate linker in the channel function, I mutated the amino acids in the cDNA of this region, and used from knowen mutaion, either alone or together with the amino acids of adjacent regions. I have emphasized the importance of the linker between these two gates - mutations in this region may cause conformational changes that play a fundamental role in mediating the coupling between the voltage sensor, activation gate and selectivity filter elements of Kv channels. I observe that free energy considerations show the significance of the coupling between the activation and inactivation gates. Moreover, a symmetry between the coupling and sensor spring strength leads to the destruction of ion conductivity. I present a thermodynamic framework for the possible study of multiple channel blocks. The arising physical perspective of the gating process gives rise to new research avenues of the coupling mode of potassium channels and may assist in explaining the centrality of the "gate linker" to the channel function.

I Ks ion-channel pore conductance can result from individual voltage sensor movements

The I _Ks current has an established role in cardiac action potential repolarization, and provides a repolarization reserve at times of stress. The underlying channels are formed from tetramers of KCNQ1 along with one to four KCNE1 accessory subunits, but how these components together gate the I _Ks complex to open the pore is controversial. Currently, either a concerted movement involving all four subunits of the tetramer or allosteric regulation of open probability through voltage-dependent subunit activation is thought to precede opening. Here, by using the E160R mutation in KCNQ1 or the F57W mutation in KCNE1 to prevent or impede, respectively, voltage sensors from moving into activated conformations, we demonstrate that a concerted transition of all four subunits after voltage sensor activation is not required for the opening of I _Ks channels. Tracking voltage sensor movement, via [2-(trimethylammonium)ethyl]methanethiosulfonate bromide (MTSET) modification and fluorescence recordings, shows that E160R-containing voltage sensors do not translocate upon depolarization. E160R, when expressed in all four KCNQ1 subunits, is nonconducting, but if one, two, or three voltage sensors contain the E160R mutation, whole-cell and single-channel currents are still observed in both the presence and absence of KCNE1, and average conductance is reduced proportional to the number of E160R voltage sensors. The data suggest that KCNQ1 + KCNE1 channels gate like KCNQ1 alone. A model of independent voltage sensors directly coupled to open states can simulate experimental changes in I _Ks current kinetics, including the nonlinear depolarization of the conductance-voltage (G-V) relationship, and tail current acceleration as the number of nonactivatable E160R subunits is increased.

Direct Evidence for a Similar Molecular Mechanism Underlying Shaker Kv Channel Fast Inactivation and Clustering

The fast inactivation and clustering functions of voltage-dependent potassium channels play fundamental roles in electrical signaling. Recent evidence suggests that both these distinct channel functions rely on intrinsically disordered N- and C-terminal cytoplasmic segments that function as entropic clocks to time channel inactivation or scaffold protein-mediated clustering, both relying on what can be described as a "ball and chain" binding mechanism. Although the mechanisms employed in each case are seemingly analogous, both were put forward based on bulky chain deletions and further exhibit differences in reaction order. These considerations raised the question of whether the molecular mechanisms underlying Kv channel fast inactivation and clustering are indeed analogous. By taking a "chain"-level chimeric channel approach involving long and short spliced inactivation or clustering "chain" variants of the Shaker Kv channel, we demonstrate the ability of native inactivation and clustering "chains" to substitute for each other in a length-dependent manner, as predicted by the "ball and chain" mechanism. Our results thus provide direct evidence arguing that the two completely unrelated Shaker Kv channel processes of fast inactivation and clustering indeed occur according to a similar molecular mechanism.

Cholesterol-Dependent Gating Effects on Ion Channels

Biomembranes separate a live cell from its environment and keep it in an off-equilibrium, steady state. They contain both phospholipids and nonphospholipids, depending on whether there are phosphate groups in the headgroup regions. Cholesterol (CHOL) is one type of nonphospholipids, and one of the most abundant lipid molecules in humans. Its content in plasma membranes and intracellular membranes varies and is tightly regulated. Voltage-gated ion channels are universally present in every cell and are fairly diversified in the eukaryotic domain of life. Our lipid-dependent gating hypothesis postulates that the controlled switch of the voltage-sensor domains (VSDs) in a voltage-gated potassium (Kv) channel between the "down" and the "up" state (gating) is sensitive to the ratio of phospholipids:nonphospholipids in the annular layer around the channel. High CHOL content is found to exert strong inhibitory effects on Kv channels. Such effects have been observed in in vitro membranes, cultured cells, and animal models for cholesterol metabolic defects. Thermodynamic analysis of the CHOL-dependent gating suggests that the inhibitory effects of CHOL result from collective interactions between annular CHOL molecules and the channel, which appear to be a more generic principle behind the CHOL effects on other ion channels and transporters. We will review the recent progress in the CHOL-dependent gating of voltage-gated ion channels, discuss the current technical limitations, and then expand briefly the learned principles to other ion channels that are known to be sensitive to the CHOL-channel interactions.

The tarantula toxin GxTx detains K+ channel gating charges in their resting conformation

Allosteric ligands modulate protein activity by altering the energy landscape of conformational space in ligand-protein complexes. Here we investigate how ligand binding to a K⁺ channel's voltage sensor allosterically modulates opening of its K⁺-conductive pore. The tarantula venom peptide guangxitoxin-1E (GxTx) binds to the voltage sensors of the rat voltage-gated K⁺ (Kv) channel Kv2.1 and acts as a partial inverse agonist. When bound to GxTx, Kv2.1 activates more slowly, deactivates more rapidly, and requires more positive voltage to reach the same K⁺-conductance as the unbound channel. Further, activation kinetics are more sigmoidal, indicating that multiple conformational changes coupled to opening are modulated. Single-channel current amplitudes reveal that each channel opens to full conductance when GxTx is bound. Inhibition of Kv2.1 channels by GxTx results from decreased open probability due to increased occurrence of long-lived closed states; the time constant of the final pore opening step itself is not impacted by GxTx. When intracellular potential is less than 0 mV, GxTx traps the gating charges on Kv2.1's voltage sensors in their most intracellular position. Gating charges translocate at positive voltages, however, indicating that GxTx stabilizes the most intracellular conformation of the voltage sensors (their resting conformation). Kinetic modeling suggests a modulatory mechanism: GxTx reduces the probability of voltage sensors activating, giving the pore opening step less frequent opportunities to occur. This mechanism results in K⁺-conductance activation kinetics that are voltage-dependent, even if pore opening (the rate-limiting step) has no inherent voltage dependence. We conclude that GxTx stabilizes voltage sensors in a resting conformation, and inhibits K⁺ currents by limiting opportunities for the channel pore to open, but has little, if any, direct effect on the microscopic kinetics of pore opening. The impact of GxTx on channel gating suggests that Kv2.1's pore opening step does not involve movement of its voltage sensors.

Voltage-Gated Sodium Channels in Neocortical Pyramidal Neurons Display Cole-Moore Activation Kinetics

Exploring the properties of action potentials is a crucial step toward a better understanding of the computational properties of single neurons and neural networks. The voltage-gated sodium channel is a key player in action potential generation. A comprehensive grasp of the gating mechanism of this channel can shed light on the biophysics of action potential generation. However, most models of voltage-gated sodium channels assume a concerted Hodgkin and Huxley kinetic gating scheme. However, it is not clear if Hodgkin and Huxley models are suitable for use in action potential simulations of central nervous system neurons. To resolve this, we investigated the activation kinetics of voltage-gated sodium channels. Here we performed high resolution voltage-clamp experiments from nucleated patches extracted from the soma of layer 5 (L5) cortical pyramidal neurons in rat brain slices. We show that the gating mechanism does not follow traditional Hodgkin and Huxley kinetics and that much of the channel voltage-dependence is probably due to rapid closed-closed transitions that lead to substantial onset latency reminiscent of the Cole-Moore effect observed in voltage-gated potassium conductances. Thus, the classical Hodgkin and Huxley description of sodium channel kinetics may be unsuitable for modeling the physiological role of this channel. Furthermore, our results reconcile between apparently contradicting studies sodium channel activation. Our findings may have key implications for the role of sodium channels in synaptic integration and action potential generation.

Gating interaction maps reveal a noncanonical electromechanical coupling mode in the Shaker K+ channel

Membrane potential regulates the activity of voltage-dependent ion channels via specialized voltage-sensing modules but the mechanisms involved in coupling voltage-sensor movement to pore opening remain unclear due to lack of resting state structures and robust methods to identify allosteric pathways. Here, using a newly developed interaction energy analysis, we probe the interfaces of the voltage-sensing and pore modules in the drosophila Shaker K⁺ channel. Our measurements reveal unexpectedly strong equilibrium gating interactions between contacts at the S4 and S5 helices in addition to those between S6 and S4–S5 linker. Network analysis of MD trajectories shows that the voltage-sensor and pore motions are linked by two distinct pathways- canonical one through the S4–S5 linker and a hitherto unknown pathway akin to rack and pinion coupling involving S4 and S5 helices. Our findings highlight the central role of the S5 helix in electromechanical transduction in the VGIC superfamily.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

Gating energetics of a voltage-dependent K+ channel pore domain

We used targeted molecular dynamics, informed by experimentally determined inter-atomic distances defining the pore region of open and closed states of the KvAP voltage-gated potassium channel, to generate a gating pathway of the pore domain in the absence of the voltage-sensing domains. We then performed umbrella sampling simulations along this pathway to calculate a potential of mean force that describes the free energy landscape connecting the closed and open conformations of the pore domain. The resulting energetic landscape displays three minima, corresponding to stable open, closed, and intermediate conformations with roughly similar stabilities. We found that the extent of hydration of the interior of the pore domain could influence the free energy landscape for pore opening/closing. © 2017 Wiley Periodicals, Inc.

In-Depth Study of the Interaction, Sensitivity, and Gating Modulation by PUFAs on K+ Channels; Interaction and New Targets

Voltage gated potassium channels (K_V) are membrane proteins that allow selective flow of K⁺ ions in a voltage-dependent manner. These channels play an important role in several excitable cells as neurons, cardiomyocytes, and vascular smooth muscle. Over the last 20 years, it has been shown that omega-3 polyunsaturated fatty acids (PUFAs) enhance or decrease the activity of several cardiac K_V channels. PUFAs-dependent modulation of potassium ion channels has been reported to be cardioprotective. However, the precise cellular mechanism underlying the cardiovascular benefits remained unclear in part because new PUFAs targets and signaling pathways continue being discovered. In this review, we will focus on recent data available concerning the following aspects of the K_V channel modulation by PUFAs: (i) the exact residues involved in PUFAs-K_V channels interaction; (ii) the structural PUFAs determinants important for their effects on K_V channels; (iii) the mechanism of the gating modulation of K_V channels and, finally, (iv) the PUFAs modulation of a few new targets present in smooth muscle cells (SMC), K_Ca1.1, K_2P, and K_ATP channels, involved in vascular relaxation.

Upward movement of IS4 and IIIS4 is a rate-limiting stage in Cav1.2 activation

In order to specify the role of individual S4 segments in CaV1.2 gating, charged residues of segments IS4-IVS4 were replaced by glutamine and the corresponding effects on activation/deactivation of calcium channel currents were analysed. Almost all replacements of charges in IS4 and IIIS4 decreased the slope of the Boltzmann curve of channel activation (activation curve) while charge neutralisations in IIS4 and IVS4 did not significantly affect the slope. S4 mutations caused either left or rightward shifts of the activation curve, and in wild-type channels, these S4 mutations hardly affected current kinetics.In slowly gating pore (S6) mutants (G432W, A780T, G1193T or A1503G), neutralisations in S4 segments significantly accelerated current kinetics. Likewise in wild type, charge replacements in IS4 and IIIS4 of pore mutants reduced the slope of the activation curves while substitutions of charges in IIS4 and IVS4 had less or no impact. We propose a gating model where the structurally different S4 segments leave their resting positions not simultaneously. Upward movement of segments IS4 and (to a lesser extend) IIIS4 appear to be a rate-limiting stage for releasing the pore gates. These segments carry most of the effective charge for channel activation. Our study suggests that S4 segments of CaV1.2 control the closed state in domain specific manner while stabilizing the open state in a non-specific manner.

The Cole-Moore Effect: Still Unexplained?

In the first issue, on the first page of the Biophysical Journal in 1960, Cole and Moore provided the first confirmation of the Hodgkin and Huxley formulation of the sodium and potassium conductances that underlie the action potential. In addition, working with the squid giant axon, Cole and Moore noted that strong hyperpolarization preceding a depolarizing voltage-clamp pulse delayed the rise of the potassium conductance: once started, the time course of the rise was always the same but after significant hyperpolarization there was a long lag before the rise began. This phenomenon has come to be known as the Cole-Moore effect. Their article examines and disproves the hypothesis that the lag reflects the time required to refill the membrane with potassium ions after the ions are swept out of the membrane into the axoplasm by hyperpolarization. The work by Cole and Moore indirectly supports the idea of a membrane channel for potassium conductance. However, the mechanism of the Cole-Moore effect remains a mystery even now, buried in the structure of the potassium channel, which was completely unknown at the time.

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.

Allosteric substrate switching in a voltage-sensing lipid phosphatase

Allostery provides a critical control over enzyme activity, biasing the catalytic site between inactive and active states. We find the Ciona intestinalis voltage-sensing phosphatase (Ci-VSP), which modifies phosphoinositide signaling lipids (PIPs), to have not one but two sequential active states with distinct substrate specificities, whose occupancy is allosterically controlled by sequential conformations of the voltage sensing domain (VSD). Using fast FRET reporters of PIPs to monitor enzyme activity and voltage clamp fluorometry to monitor conformational changes in the VSD, we find that Ci-VSP switches from inactive to a PIP₃-preferring active state when the VSD undergoes an initial voltage sensing motion and then into a second PIP₂-preferring active state when the VSD activates fully. This novel 2-step allosteric control over a dual specificity enzyme enables voltage to shape PIP concentrations in time, and provides a mechanism for the complex modulation of PIP-regulated ion channels, transporters, cell motility and endo/exocytosis.

Functional diversity of potassium channel voltage-sensing domains

Voltage-gated potassium channels or Kv's are membrane proteins with fundamental physiological roles. They are composed of 2 main functional protein domains, the pore domain, which regulates ion permeation, and the voltage-sensing domain, which is in charge of sensing voltage and undergoing a conformational change that is later transduced into pore opening. The voltage-sensing domain or VSD is a highly conserved structural motif found in all voltage-gated ion channels and can also exist as an independent feature, giving rise to voltage sensitive enzymes and also sustaining proton fluxes in proton-permeable channels. In spite of the structural conservation of VSDs in potassium channels, there are several differences in the details of VSD function found across variants of Kvs. These differences are mainly reflected in variations in the electrostatic energy needed to open different potassium channels. In turn, the differences in detailed VSD functioning among voltage-gated potassium channels might have physiological consequences that have not been explored and which might reflect evolutionary adaptations to the different roles played by Kv channels in cell physiology.

Alkanols inhibit voltage-gated K(+) channels via a distinct gating modifying mechanism that prevents gate opening

Alkanols are small aliphatic compounds that inhibit voltage-gated K⁺ (K_v) channels through a yet unresolved gating mechanism. K_v channels detect changes in the membrane potential with their voltage-sensing domains (VSDs) that reorient and generate a transient gating current. Both 1-Butanol (1-BuOH) and 1-Hexanol (1-HeOH) inhibited the ionic currents of the Shaker K_v channel in a concentration dependent manner with an IC₅₀ value of approximately 50 mM and 3 mM, respectively. Using the non-conducting Shaker-W434F mutant, we found that both alkanols immobilized approximately 10% of the gating charge and accelerated the deactivating gating currents simultaneously with ionic current inhibition. Thus, alkanols prevent the final VSD movement(s) that is associated with channel gate opening. Applying 1-BuOH and 1-HeOH to the Shaker-P475A mutant, in which the final gating transition is isolated from earlier VSD movements, strengthened that neither alkanol affected the early VSD movements. Drug competition experiments showed that alkanols do not share the binding site of 4-aminopyridine, a drug that exerts a similar effect at the gating current level. Thus, alkanols inhibit Shaker-type K_v channels via a unique gating modifying mechanism that stabilizes the channel in its non-conducting activated state.

A Disease Mutation Causing Episodic Ataxia Type I in the S1 Links Directly to the Voltage Sensor and the Selectivity Filter in Kv Channels

The mutation F184C in Kv1.1 leads to development of episodic ataxia type I (EA1). Although the mutation has been said to alter activation kinetics and to lower expression, we show here that the underlying molecular mechanisms may be more complex. Although F184 is positioned in the "peripheral" S1 helix, it occupies a central position in the 3D fold. We show in cut-open oocyte voltage-clamp recordings of gating and ionic currents of the Shaker Kv channel expressed in Xenopus oocytes that F184 not only interacts directly with the gating charges of the S4, but also creates a functional link to the selectivity filter of the neighboring subunit. This link leads to impaired fast and slow inactivation. The effect on fast inactivation is of an allosteric nature considering that fast inactivation is caused by a linked cytosolic ball peptide. The extensive effects of F184C provide a new mechanism underlying EA.

SIGNIFICANCE STATEMENT

Episodic ataxia (EA) is an inherited disease that leads to occasional loss of motor control in combination with variable other symptoms such as vertigo or migraine. EA type I (EA1), studied here, is caused by mutations in a voltage-gated potassium channel that contributes to the generation of electrical signals in the brain. The mechanism by which mutations in voltage-gated potassium channels lead to EA is still unknown and there is no consistent pharmacological treatment. By studying in detail one disease-causing mutation in Kv1.1, we describe a novel molecular mechanism distinct from mechanisms described previously. This mechanism contributes to the understanding of potassium channel function in general and might lead to a better understanding of how EA develops.

Positive Allosteric Modulation of Kv Channels by Sevoflurane: Insights into the Structural Basis of Inhaled Anesthetic Action

Inhalational general anesthesia results from the poorly understood interactions of haloethers with multiple protein targets, which prominently includes ion channels in the nervous system. Previously, we reported that the commonly used inhaled anesthetic sevoflurane potentiates the activity of voltage-gated K⁺ (Kv) channels, specifically, several mammalian Kv1 channels and the Drosophila K-Shaw2 channel. Also, previous work suggested that the S4-S5 linker of K-Shaw2 plays a role in the inhibition of this Kv channel by n-alcohols and inhaled anesthetics. Here, we hypothesized that the S4-S5 linker is also a determinant of the potentiation of Kv1.2 and K-Shaw2 by sevoflurane. Following functional expression of these Kv channels in Xenopus oocytes, we found that converse mutations in Kv1.2 (G329T) and K-Shaw2 (T330G) dramatically enhance and inhibit the potentiation of the corresponding conductances by sevoflurane, respectively. Additionally, Kv1.2-G329T impairs voltage-dependent gating, which suggests that Kv1.2 modulation by sevoflurane is tied to gating in a state-dependent manner. Toward creating a minimal Kv1.2 structural model displaying the putative sevoflurane binding sites, we also found that the positive modulations of Kv1.2 and Kv1.2-G329T by sevoflurane and other general anesthetics are T1-independent. In contrast, the positive sevoflurane modulation of K-Shaw2 is T1-dependent. In silico docking and molecular dynamics-based free-energy calculations suggest that sevoflurane occupies distinct sites near the S4-S5 linker, the pore domain and around the external selectivity filter. We conclude that the positive allosteric modulation of the Kv channels by sevoflurane involves separable processes and multiple sites within regions intimately involved in channel gating.

Domain and interdomain energetics underlying gating in Shaker-type Kv channels

To understand gating events with a time-base many orders-of-magnitude slower than that of atomic motion in voltage-gated ion channels such as the Shaker-type KV channels, a multiscale physical model is constructed from the experimentally well-characterized voltage-sensor (VS) domains coupled to a hydrophobic gate. The four VS domains are described by a continuum electrostatic model under voltage-clamp conditions, the control of ion flow by the gate domain is described by a vapor-lock mechanism, and the simple coupling principle is informed by known experimental results and trial-and-error. The configurational energy computed for each element is used to produce a total Hamiltonian that is a function of applied voltage, VS positions, and gate radius. We compute statistical-mechanical expectation values of macroscopic laboratory observables. This approach stands in contrast with molecular-dynamic models which are challenged by increasing scale, and kinetic models which assume a probability distribution rather than derive it from the underlying physics. This generic model predicts well the Shaker charge/voltage and conductance/voltage relations; the tight constraints underlying these results allow us to quantitatively assess the underlying physical mechanisms. The total electrical work picked up by the VS domains is an order-of-magnitude larger than the work required to actuate the gate itself, suggesting an energetic basis for the evolutionary flexibility of the voltage-gating mechanism. The cooperative slide-and-interlock behavior of the VS domains described by the VS-gate coupling relation leads to the experimentally observed bistable gating. This engineering approach should prove useful in the investigation of various elements underlying gating characteristics and degraded behavior due to mutation.

Voltage Sensing in Membranes: From Macroscopic Currents to Molecular Motions

Voltage-sensing domains (VSDs) are integral membrane protein units that sense changes in membrane electric potential, and through the resulting conformational changes, regulate a specific function. VSDs confer voltage-sensitivity to a large superfamily of membrane proteins that includes voltage-gated Na[Formula: see text], K[Formula: see text], Ca[Formula: see text] ,and H[Formula: see text] selective channels, hyperpolarization-activated cyclic nucleotide-gated channels, and voltage-sensing phosphatases. VSDs consist of four transmembrane segments (termed S1 through S4). Their most salient structural feature is the highly conserved positions for charged residues in their sequences. S4 exhibits at least three conserved triplet repeats composed of one basic residue (mostly arginine) followed by two hydrophobic residues. These S4 basic side chains participate in a state-dependent internal salt-bridge network with at least four acidic residues in S1-S3. The signature of voltage-dependent activation in electrophysiology experiments is a transient current (termed gating or sensing current) upon a change in applied membrane potential as the basic side chains in S4 move across the membrane electric field. Thus, the unique structural features of the VSD architecture allow for competing requirements: maintaining a series of stable transmembrane conformations, while allowing charge motion, as briefly reviewed here.

Voltage-dependent gating and gating charge measurements in the Kv1.2 potassium channel

Kv1.2’s gating charge is less than Shaker’s, and the specific contributions of charged S4 residues differ, suggesting that the electric field distribution in the Kv1.2 voltage-sensing domain is different than Shaker’s.

Much has been learned about the voltage sensors of ion channels since the x-ray structure of the mammalian voltage-gated potassium channel Kv1.2 was published in 2005. High resolution structural data of a Kv channel enabled the structural interpretation of numerous electrophysiological findings collected in various ion channels, most notably Shaker, and permitted the development of meticulous computational simulations of the activation mechanism. The fundamental premise for the structural interpretation of functional measurements from Shaker is that this channel and Kv1.2 have the same characteristics, such that correlation of data from both channels would be a trivial task. We tested these assumptions by measuring Kv1.2 voltage-dependent gating and charge per channel. We found that the Kv1.2 gating charge is near 10 elementary charges (e_o), ∼25% less than the well-established 13–14 e_o in Shaker. Next, we neutralized positive residues in the Kv1.2 S4 transmembrane segment to investigate the cause of the reduction of the gating charge and found that, whereas replacing R1 with glutamine decreased voltage sensitivity to ∼50% of the wild-type channel value, mutation of the subsequent arginines had a much smaller effect. These data are in marked contrast to the effects of charge neutralization in Shaker, where removal of the first four basic residues reduces the gating charge by roughly the same amount. In light of these differences, we propose that the voltage-sensing domains (VSDs) of Kv1.2 and Shaker might undergo the same physical movement, but the septum that separates the aqueous crevices in the VSD of Kv1.2 might be thicker than Shaker’s, accounting for the smaller Kv1.2 gating charge.

Diffusion approximation-based simulation of stochastic ion channels: which method to use?

To study the effects of stochastic ion channel fluctuations on neural dynamics, several numerical implementation methods have been proposed. Gillespie's method for Markov Chains (MC) simulation is highly accurate, yet it becomes computationally intensive in the regime of a high number of channels. Many recent works aim to speed simulation time using the Langevin-based Diffusion Approximation (DA). Under this common theoretical approach, each implementation differs in how it handles various numerical difficulties—such as bounding of state variables to [0,1]. Here we review and test a set of the most recently published DA implementations (Goldwyn et al., 2011; Linaro et al., 2011; Dangerfield et al., 2012; Orio and Soudry, 2012; Schmandt and Galán, 2012; Güler, 2013; Huang et al., 2013a), comparing all of them in a set of numerical simulations that assess numerical accuracy and computational efficiency on three different models: (1) the original Hodgkin and Huxley model, (2) a model with faster sodium channels, and (3) a multi-compartmental model inspired in granular cells. We conclude that for a low number of channels (usually below 1000 per simulated compartment) one should use MC—which is the fastest and most accurate method. For a high number of channels, we recommend using the method by Orio and Soudry (2012), possibly combined with the method by Schmandt and Galán (2012) for increased speed and slightly reduced accuracy. Consequently, MC modeling may be the best method for detailed multicompartment neuron models—in which a model neuron with many thousands of channels is segmented into many compartments with a few hundred channels.

The differential contribution of GluN1 and GluN2 to the gating operation of the NMDA receptor channel

The Ν-methyl-D-aspartate (NMDA) receptor channel is an obligatory heterotetramer formed by two GluN1 and two GluN2 subunits. However, the differential contribution of the two different subunits to channel operation is not clear. We found that the apparent affinity of glycine to GluN1 (K gly ∼ 0.6 μM) is much higher than NMDA or glutamate to GluN2 (K NMDA ∼ 36 μM, K glu ∼ 4.8 μM). The binding rate constant (derived from the linear regression of the apparent macroscopic binding rates) of glycine to GluN1 (∼9.8 × 10(6) M(-1) s(-1)), however, is only slightly faster than NMDA to GluN2 (∼4.1 × 10(6) M(-1) s(-1)). Accordingly, the apparent unbinding rates of glycine from activated GluN1 (time constant ∼2 s) are much slower than NMDA from activated GluN2 (time constant ∼70 ms). Moreover, the decay of NMDA currents upon wash-off of both glycine and NMDA seems to follow the course of NMDA rather than glycine unbinding. But if only glycine is washed off, the current decay is much slower, apparently following the course of glycine unbinding. The apparent binding rate of glycine to the fully deactivated channel, in the absence of NMDA, is roughly the same as that measured with co-application of both ligands, whereas the apparent binding rate of NMDA to the fully deactivated channel in the absence of glycine is markedly slower. In this regard, it is interesting that the seventh residue in the highly conserved SYTANLAAF motif (A7) in GluN1 and GluN2 are so close that they may interact with each other to control the dimension of the external pore mouth. Moreover, specific mutations involving A7 in GluN1 but not in GluN2 result in channels showing markedly enhanced affinity to both glycine and NMDA and readily activated by only NMDA, as if the channel is already partially activated. We conclude that GluN2 is most likely directly responsible for the activation gate of the NMDA channel, whereas GluN1 assumes a role of more global control, especially on the gating conformational changes in GluN2. Structurally, this intersubunit regulatory interaction seems to involve the SYTANLAAF motif, especially the A7 residue.

A self-consistent approach for determining pairwise interactions that underlie channel activation

Net free-energy measurements can be combined with mutant cycle analysis to determine interaction energies between specific amino acid pairs during channel activation.

Signaling proteins such as ion channels largely exist in two functional forms, corresponding to the active and resting states, connected by multiple intermediates. Multiparametric kinetic models based on sophisticated electrophysiological experiments have been devised to identify molecular interactions of these conformational transitions. However, this approach is arduous and is not suitable for large-scale perturbation analysis of interaction pathways. Recently, we described a model-free method to obtain the net free energy of activation in voltage- and ligand-activated ion channels. Here we extend this approach to estimate pairwise interaction energies of side chains that contribute to gating transitions. Our approach, which we call generalized interaction-energy analysis (GIA), combines median voltage estimates obtained from charge-voltage curves with mutant cycle analysis to ascertain the strengths of pairwise interactions. We show that, for a system with an arbitrary gating scheme, the nonadditive contributions of amino acid pairs to the net free energy of activation can be computed in a self-consistent manner. Numerical analyses of sequential and allosteric models of channel activation also show that this approach can measure energetic nonadditivities even when perturbations affect multiple transitions. To demonstrate the experimental application of this method, we reevaluated the interaction energies of six previously described long-range interactors in the Shaker potassium channel. Our approach offers the ability to generate detailed interaction energy maps in voltage- and ligand-activated ion channels and can be extended to any force-driven system as long as associated “displacement” can be measured.

Mutations in the S6 gate isolate a late step in the activation pathway and reduce 4-AP sensitivity in shaker K(v) channel

Kv channels detect changes in the membrane potential via their voltage-sensing domains (VSDs) that control the status of the S6 bundle crossing (BC) gate. The movement of the VSDs results in a transfer of the S4 gating charges across the cell membrane but only the last 10-20% of the total gating charge movement is associated with BC gate opening, which involves cooperative transition(s) in the subunits. Substituting the proline residue P475 in the S6 of the Shaker channel by a glycine or alanine causes a considerable shift in the voltage-dependence of the cooperative transition(s) of BC gate opening, effectively isolating the late gating charge component from the other gating charge that originates from earlier VSD movements. Interestingly, both mutations also abolished Shaker's sensitivity to 4-aminopyridine, which is a pharmacological tool to isolate the late gating charge component. The alanine substitution (that would promote a α-helical configuration compared to proline) resulted in the largest separation of both gating charge components; therefore, BC gate flexibility appears to be important for enabling the late cooperative step of channel opening.

A molecular framework for temperature-dependent gating of ion channels

Perception of heat or cold in higher organisms is mediated by specialized ion channels whose gating is exquisitely sensitive to temperature. The physicochemical underpinnings of this temperature-sensitive gating have proven difficult to parse. Here, we took a bottom-up protein design approach and rationally engineered ion channels to activate in response to thermal stimuli. By varying amino acid polarities at sites undergoing state-dependent changes in solvation, we were able to systematically confer temperature sensitivity to a canonical voltage-gated ion channel. Our results imply that the specific heat capacity change during channel gating is a major determinant of thermosensitive gating. We also show that reduction of gating charges amplifies temperature sensitivity of designer channels, which accounts for low-voltage sensitivity in all known temperature-gated ion channels. These emerging principles suggest a plausible molecular mechanism for temperature-dependent gating that reconcile how ion channels with an overall conserved transmembrane architecture may exhibit a wide range of temperature-sensing phenotypes. :

cAMP control of HCN2 channel Mg2+ block reveals loose coupling between the cyclic nucleotide-gating ring and the pore

Hyperpolarization-activated cyclic nucleotide-regulated HCN channels underlie the Na⁺-K⁺ permeable I_H pacemaker current. As with other voltage-gated members of the 6-transmembrane K_V channel superfamily, opening of HCN channels involves dilation of a helical bundle formed by the intracellular ends of S6 albeit this is promoted by inward, not outward, displacement of S4. Direct agonist binding to a ring of cyclic nucleotide-binding sites, one of which lies immediately distal to each S6 helix, imparts cAMP sensitivity to HCN channel opening. At depolarized potentials, HCN channels are further modulated by intracellular Mg²⁺ which blocks the open channel pore and blunts the inhibitory effect of outward K⁺ flux. Here, we show that cAMP binding to the gating ring enhances not only channel opening but also the kinetics of Mg²⁺ block. A combination of experimental and simulation studies demonstrates that agonist acceleration of block is mediated via acceleration of the blocking reaction itself rather than as a secondary consequence of the cAMP enhancement of channel opening. These results suggest that the activation status of the gating ring and the open state of the pore are not coupled in an obligate manner (as required by the often invoked Monod-Wyman-Changeux allosteric model) but couple more loosely (as envisioned in a modular model of protein activation). Importantly, the emergence of second messenger sensitivity of open channel rectification suggests that loose coupling may have an unexpected consequence: it may endow these erstwhile “slow” channels with an ability to exert voltage and ligand-modulated control over cellular excitability on the fastest of physiologically relevant time scales.

Modeling ion channels: past, present, and future

Ion channels are membrane-bound enzymes whose catalytic sites are ion-conducting pores that open and close (gate) in response to specific environmental stimuli. Ion channels are important contributors to cell signaling and homeostasis. Our current understanding of gating is the product of 60 plus years of voltage-clamp recording augmented by intervention in the form of environmental, chemical, and mutational perturbations. The need for good phenomenological models of gating has evolved in parallel with the sophistication of experimental technique. The goal of modeling is to develop realistic schemes that not only describe data, but also accurately reflect mechanisms of action. This review covers three areas that have contributed to the understanding of ion channels: traditional Eyring kinetic theory, molecular dynamics analysis, and statistical thermodynamics. Although the primary emphasis is on voltage-dependent channels, the methods discussed here are easily generalized to other stimuli and could be applied to any ion channel and indeed any macromolecule.

Kv2 channel regulation of action potential repolarization and firing patterns in superior cervical ganglion neurons and hippocampal CA1 pyramidal neurons

Kv2 family "delayed-rectifier" potassium channels are widely expressed in mammalian neurons. Kv2 channels activate relatively slowly and their contribution to action potential repolarization under physiological conditions has been unclear. We explored the function of Kv2 channels using a Kv2-selective blocker, Guangxitoxin-1E (GxTX-1E). Using acutely isolated neurons, mixed voltage-clamp and current-clamp experiments were done at 37°C to study the physiological kinetics of channel gating and action potentials. In both rat superior cervical ganglion (SCG) neurons and mouse hippocampal CA1 pyramidal neurons, 100 nm GxTX-1E produced near-saturating block of a component of current typically constituting ∼60-80% of the total delayed-rectifier current. GxTX-1E also reduced A-type potassium current (IA), but much more weakly. In SCG neurons, 100 nm GxTX-1E broadened spikes and voltage clamp experiments using action potential waveforms showed that Kv2 channels carry ∼55% of the total outward current during action potential repolarization despite activating relatively late in the spike. In CA1 neurons, 100 nm GxTX-1E broadened spikes evoked from -70 mV, but not -80 mV, likely reflecting a greater role of Kv2 when other potassium channels were partially inactivated at -70 mV. In both CA1 and SCG neurons, inhibition of Kv2 channels produced dramatic depolarization of interspike voltages during repetitive firing. In CA1 neurons and some SCG neurons, this was associated with increased initial firing frequency. In all neurons, inhibition of Kv2 channels depressed maintained firing because neurons entered depolarization block more readily. Therefore, Kv2 channels can either decrease or increase neuronal excitability depending on the time scale of excitation.

Functional interactions of voltage sensor charges with an S2 hydrophobic plug in hERG channels

Human ether-à-go-go–related gene (hERG, Kv11.1) potassium channels have unusually slow activation and deactivation kinetics. It has been suggested that, in fast-activating Shaker channels, a highly conserved Phe residue (F290) in the S2 segment forms a putative gating charge transfer center that interacts with S4 gating charges, i.e., R362 (R1) and K374 (K5), and catalyzes their movement across the focused electric field. F290 is conserved in hERG (F463), but the relevant residues in the hERG S4 are reversed, i.e., K525 (K1) and R537 (R5), and there is an extra positive charge adjacent to R537 (i.e., K538). We have examined whether hERG channels possess a transfer center similar to that described in Shaker and if these S4 charge differences contribute to slow gating in hERG channels. Of five hERG F463 hydrophobic substitutions tested, F463W and F463Y shifted the conductance–voltage (G-V) relationship to more depolarized potentials and dramatically slowed channel activation. With the S4 residue reversals (i.e., K525, R537) taken into account, the closed state stabilization by F463W is consistent with a role for F463 that is similar to that described for F290 in Shaker. As predicted from results with Shaker, the hERG K525R mutation destabilized the closed state. However, hERG R537K did not stabilize the open state as predicted. Instead, we found the neighboring K538 residue to be critical for open state stabilization, as K538R dramatically slowed and right-shifted the voltage dependence of activation. Finally, double mutant cycle analysis on the G-V curves of F463W/K525R and F463W/K538R double mutations suggests that F463 forms functional interactions with K525 and K538 in the S4 segment. Collectively, these data suggest a role for F463 in mediating closed–open equilibria, similar to that proposed for F290 in Shaker channels.

Do lipids show state-dependent affinity to the voltage-gated potassium channel KvAP?

As all integral membrane proteins, voltage-gated ion channels are embedded in a lipid matrix that regulates their channel behavior either by physicochemical properties or by direct binding. Because manipulation of the lipid composition in cells is difficult, we investigated the influence of different lipids on purified KvAP channels reconstituted in planar lipid bilayers of known composition. Lipids developed two distinct and independent effects on the KvAP channels; lipids interacting with the pore lowered the energy barriers for the final transitions, whereas voltage sensor-bound lipids shifted the midpoint of activation dependent on their electrostatic charge. Above all, the midpoint of activation was determined only by those lipids the channels came in contact with first after purification and can seemingly only be exchanged if the channel resides in the open state. The high affinity of the bound lipids to the binding site has implications not only on our understanding of the gating mechanism but also on the general experimental design of any lipid dependence study.

Shaker IR T449 mutants separate C- from U-type inactivation

Previous studies demonstrated that slow inactivation of the Shaker potassium channel can be made ~100-fold faster or slower by point mutations at a site in the outer pore (T449). However, the discovery that two forms of slow inactivation coexist in Shaker raises the question of which inactivation process is affected by mutation. Equivalent mutations in K(V)2.1, a channel exhibiting only U-type inactivation, have minimal effects on inactivation, suggesting that mutation of Shaker T449 acts on C-type inactivation alone, a widely held yet untested hypothesis. This study reexamines mutations at Shaker T449, confirming that T449A speeds inactivation and T449Y/V slow it. T449Y and T449V exhibit U-type inactivation that is enhanced by high extracellular potassium, in contrast to C-type inactivation in T449A which is inhibited by high potassium. Automated parameter estimation for a 12-state Markov model suggests that U-type inactivation occurs mainly from closed states upon weak depolarization, but primarily from the open state at positive voltages. The model also suggests that WT channels, which in this study exhibit mostly C-type inactivation, recover from inactivation through closed-inactivated states, producing voltage-dependent recovery. This suggests that both C-type and U-type inactivation involve both open-inactivated and closed-inactivated states.

Exponential sum-fitting of dwell-time distributions without specifying starting parameters

Fitting dwell-time distributions with sums of exponentials is widely used to characterize histograms of open- and closed-interval durations recorded from single ion channels, as well as for other physical phenomena. However, it can be difficult to identify the contributing exponential components. Here we extend previous methods of exponential sum-fitting to present a maximum-likelihood approach that consistently detects all significant exponentials without the need for user-specified starting parameters. Instead of searching for exponentials, the fitting starts with a very large number of initial exponentials with logarithmically spaced time constants, so that none are missed. Maximum-likelihood fitting then determines the areas of all the initial exponentials keeping the time constants fixed. In an iterative manner, with refitting after each step, the analysis then removes exponentials with negligible area and combines closely spaced adjacent exponentials, until only those exponentials that make significant contributions to the dwell-time distribution remain. There is no limit on the number of significant exponentials and no starting parameters need be specified. We demonstrate fully automated detection for both experimental and simulated data, as well as for classical exponential-sum-fitting problems.

Minimalist Approaches to Protein Labelling: Getting the Most Fluorescent Bang for Your Steric Buck

Fluorescence methods allow one to monitor protein conformational changes, protein–protein associations, and proteolysis in real time, at the single molecule level and in living cells. The information gained in such experiments is a function of the spectroscopic techniques used and the strategic placement of fluorophore labels within the protein structure. There is often a trade-off between size and utility for fluorophores, whereby large size can be disruptive to the protein’s fold or function, but valuable characteristics, such as visible wavelength absorption and emission or brightness, require sizable chromophores. Three major types of fluorophore readouts are commonly used: (1) Förster resonance energy transfer (FRET); (2) photoinduced electron transfer (PET); and (3) environmental sensitivity. This review focuses on those probes small enough to be incorporated into proteins during ribosomal translation, which allows the probes to be placed on the interiors of proteins as they are folded during synthesis. The most broadly useful method for doing so is site-specific unnatural amino acid (UAA) mutagenesis. We discuss the use of UAA probes in applications relying on FRET, PET, and environmental sensitivity. We also briefly review other methods of protein labelling and compare their relative merits to UAA mutagenesis. Finally, we discuss small probes that have thus far been used only in synthetic peptides, but which have unusual value and may be candidates for incorporation using UAA methods.

Structure and function of the thermoTRP channel pore

Temperature-sensitive transient receptor potential (TRP) channels are structurally similar to other tetrameric cation channels, but can be potently activated by heat. Recent studies suggest that the pore-forming region directly participates in activation gating. In this chapter, we summarize major findings from both structural and functional studies concerning the gating role of the pore region, focusing in particular on TRPV1. The emerging picture is that the peripheral S1-S4 region of TRPV1 is rigid and plays a supporting role for the pore to undergo conformational rearrangements. This places the pore region in the center of activation gating.

Omega currents in voltage-gated ion channels: what can we learn from uncovering the voltage-sensing mechanism using MD simulations?

Ion channels conduct charged species through otherwise impermeable biological membranes. Their activity supports a number of physiological processes, and genetic mutations can disrupt their function dramatically. Among these channels, voltage gated cation channels (VGCCs) are ubiquitous transmembrane proteins involved in electrical signaling. In addition to their selectivity for ions, their function requires membrane-polarization-dependent gating. Triggered by changes in the transmembrane voltage, the activation and deactivation of VGCCs proceed through a sensing mechanism that prompts motion of conserved positively charged (basic) residues within the S4 helix of a four-helix bundle, the voltage sensor domain (VSD). Decades of experimental investigations, using electrophysiology, molecular biology, pharmacology, and spectroscopy, have revealed details about the function of VGCCs. However, in 2005, the resolution of the crystal structure of the activated state of one member of the mammalian voltage gated potassium (Kv) channels family (the Kv1.2) enabled researchers to make significant progress in understanding the structure-function relationship in these proteins on a molecular level. In this Account, we review the use of a complementary technique, molecular dynamics (MD) simulations, that has offered new insights on this timely issue. Starting from the "open-activated state" crystal structure, we have carried out large-scale all atom MD simulations of the Kv1.2 channel embedded in its lipidic environment and submitted to a hyperpolarizing (negative) transmembrane potential. We then used steered MD simulations to complete the full transition to the resting-closed state. Using these procedures, we have followed the operation of the VSDs and uncovered three intermediate states between their activated and deactivated conformations. Each conformational state is characterized by its network of salt bridges and by the occupation of the gating charge transfer center by a specific S4 basic residue. Overall, the global deactivation mechanism that we have uncovered agrees with proposed kinetic models and recent experimental results that point towards the presence of several intermediate states. The understanding of these conformations has allowed us to examine how mutations of the S4 basic residues analogous to those involved in genetic diseases affect the function of VGCCs. In agreement with electrophysiology experiments, mutations perturb the VSD structure and trigger the appearance of state-dependent "leak" currents. The simulation results unveil the key elementary molecular processes involved in these so-called "omega" currents. We generalize these observations to other members of the VGCC family, indicating which type of residues may generate such currents and which conditions might cause leaks that prevent proper function of the channel. Today, the understanding of the intermediate state conformations enables researchers to confidently tackle other key questions such as the mode of action of toxins or modulation of channel function by lipids.