Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations

Budding of Asymmetric Lipid Bilayers: Effects of Cholesterol, Anionic Lipid, and Electric Field

Membrane budding is vital for various cellular processes such as synaptic activity regulation, vesicle transport and release, and endocytosis/exocytosis. Although protein-mediated membrane budding has been extensively investigated, the effects of the lipid asymmetry of the two leaflets and the asymmetrically electrical environments of the cellular membrane on membrane budding remain elusive. In this work, using coarse-grained molecular dynamics simulations, we systematically investigate the impacts of lipid bilayer asymmetry and external electric fields mimicking the asymmetric membrane potential on the membrane budding. The results show that the differential stress induced by the asymmetric distribution of lipids in the two leaflets is a crucial factor for the membrane budding. The unidirectional flip of cholesterol induced by the membrane curvature and the asymmetric ion adsorption induced by the anionic lipids promote the budding process. Furthermore, the external electric field applied perpendicularly to the bilayer plane increases the transmembrane potential and produces an additional differential stress across the leaflets by imposing an asymmetric torque on the lipid headgroups in the two leaflets, facilitating the membrane budding. These findings offer insights into how the structural and the environmental asymmetry in natural cellular membranes influence membrane budding in cellular processes.

A Deep Learning Approach to Uncover Voltage-Gated Ion Channels' Intermediate States

Owing to recent advancements in cryo-electron microscopy, voltage-gated ion channels have gained a greater comprehension of their structural characteristics. However, a significant enigma remains unsolved for a large majority of these channels: their gating mechanism. This mechanism, which encompasses the conformational changes between open and closed states, is pivotal to their proper functioning. Beyond the binary states of open and closed, an ensemble of intermediate states defines the transition path in-between. Due to the lack of experimental data, one might resort to molecular dynamics simulations as an alternative to decipher these states and the transitions between them. However, the high-energy barriers and the colossal time scales involved hinder access to the latter. We present here an application of deep learning as a reliable pipeline for a comprehensive exploration of voltage-gated ion channel conformational rearrangements during gating. We showcase the pipeline performance specifically on the Kv1.2 voltage sensor domain and confront the results with existing data. We demonstrate how our physics-based deep learning approach contributes to the theoretical understanding of these channels and how it might provide further insights into the exploration of channelopathies.

Dual mechanisms contribute to enhanced voltage dependence of an electric fish potassium channel

The voltage dependence of different voltage-gated potassium channels, described by the voltage at which half of the channels are open (V1/2), varies over a range of 80 mV and is influenced by factors such as the number of positive gating charges and the identity of the hydrophobic amino acids in the channel's voltage sensor (S4). Here we explore by experimental manipulations and molecular dynamics simulation the contributions of two derived features of an electric fish potassium channel (Kv1.7a) that is among the most voltage-sensitive Shaker family potassium channels known. These are a patch of four contiguous negatively charged glutamates in the S3-S4 extracellular loop and a glutamate in the S3b helix. We find that these negative charges affect V1/2 by separate, complementary mechanisms. In the closed state, the S3-S4 linker negative patch reduces the membrane surface charge biasing the channel to enter the open state while, upon opening, the negative amino acid in the S3b helix faces the second (R2) gating charge of the voltage sensor electrostatically biasing the channel to remain in the open state. This work highlights two evolutionary novelties that illustrate the potential influence of negatively charged amino acids in extracellular loops and adjacent helices to voltage dependence.

In Silico Methods for the Discovery of Kv7.2/7.3 Channels Modulators: A Comprehensive Review

The growing interest in Kv7.2/7.3 agonists originates from the involvement of these channels in several brain hyperexcitability disorders. In particular, Kv7.2/7.3 mutants have been clearly associated with epileptic encephalopathies (DEEs) as well as with a spectrum of focal epilepsy disorders, often associated with developmental plateauing or regression. Nevertheless, there is a lack of available therapeutic options, considering that retigabine, the only molecule used in clinic as a broad-spectrum Kv7 agonist, has been withdrawn from the market in late 2016. This is why several efforts have been made both by both academia and industry in the search for suitable chemotypes acting as Kv7.2/7.3 agonists. In this context, in silico methods have played a major role, since the precise structures of different Kv7 homotetramers have been only recently disclosed. In the present review, the computational methods used for the design of Kv.7.2/7.3 small molecule agonists and the underlying medicinal chemistry are discussed in the context of their biological and structure-function properties.

Carboxyl-group compounds activate voltage-gated potassium channels via a distinct mechanism

Voltage-gated ion channels are responsible for the electrical excitability of neurons and cardiomyocytes. Thus, they are obvious targets for pharmaceuticals aimed to modulate excitability. Compounds activating voltage-gated potassium (KV) channels are expected to reduce excitability. To search for new KV-channel activators, we performed a high-throughput screen of 10,000 compounds on a specially designed Shaker KV channel. Here, we report on a large family of channel-activating compounds with a carboxyl (COOH) group as the common motif. The most potent COOH activators are lipophilic (4 < LogP <7) and are suggested to bind at the interface between the lipid bilayer and the channel's positively charged voltage sensor. The negatively charged form of the COOH-group compounds is suggested to open the channel by electrostatically pulling the voltage sensor to an activated state. Several of the COOH-group compounds also activate the therapeutically important KV7.2/7.3 channel and can thus potentially be developed into antiseizure drugs. The COOH-group compounds identified in this study are suggested to act via the same site and mechanism of action as previously studied COOH-group compounds, such as polyunsaturated fatty acids and resin acids, but distinct from sites for several other types of potassium channel-activating compounds.

Structural Insights and Influence of Terahertz Waves in Midinfrared Region on Kv1.2 Channel Selectivity Filter

Potassium ion channels are the structural basis for excitation transmission, heartbeat, and other biological processes. The selectivity filter is a critical structural component of potassium ion channels, whose structure is crucial to realizing their function. As biomolecules vibrate and rotate at frequencies in the terahertz band, potassium ion channels are sensitive to terahertz waves. Therefore, it is worthwhile to investigate how the terahertz wave influences the selectivity filter of the potassium channels. In this study, we investigate the structure of the selectivity filter of Kv1.2 potassium ion channels using molecular dynamics simulations. The effect of an electric field on the channel has been examined at four different resonant frequencies of the carbonyl group in SF: 36.75 37.06, 37.68, and 38.2 THz. As indicated by the results, 376GLY appears to be the critical residue in the selectivity filter of the Kv1.2 channel. Its dihedral angle torsion is detrimental to the channel structural stability and the transmembrane movement of potassium ions. 36.75 THz is the resonance frequency of the carbonyl group of 376GLY. Among all four frequencies explored, the applied terahertz electric field of this frequency has the most significant impact on the channel structure, negatively impacting the channel stability and reducing the ion permeability by 20.2% compared to the absence of fields. In this study, we simulate that terahertz waves in the mid-infrared frequency region can significantly alter the structure and function of potassium ion channels and that the effects of terahertz waves differ greatly based on frequency.

Dynamics of activation in the voltage-sensing domain of Ciona intestinalis phosphatase Ci-VSP

The Ciona intestinalis voltage-sensing phosphatase (Ci-VSP) is a membrane protein containing a voltage-sensing domain (VSD) that is homologous to VSDs from voltage-gated ion channels responsible for cellular excitability. Previously published crystal structures of Ci-VSD in putative resting and active conformations suggested a helical-screw voltage sensing mechanism in which the S4 helix translocates and rotates to enable exchange of salt-bridge partners, but the microscopic details of the transition between the resting and active conformations remained unknown. Here, by combining extensive molecular dynamics simulations with a recently developed computational framework based on dynamical operators, we elucidate the microscopic mechanism of the resting-active transition at physiological membrane potential. Sparse regression reveals a small set of coordinates that distinguish intermediates that are hidden from electrophysiological measurements. The intermediates arise from a noncanonical helical-screw mechanism in which translocation, rotation, and side-chain movement of the S4 helix are only loosely coupled. These results provide insights into existing experimental and computational findings on voltage sensing and suggest ways of further probing its mechanism.

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.

Protein-Mediated Electroporation in a Cardiac Voltage-Sensing Domain Due to an nsPEF Stimulus

This study takes a step in understanding the physiological implications of the nanosecond pulsed electric field (nsPEF) by integrating molecular dynamics simulations and machine learning techniques. nsPEF, a state-of-the-art technology, uses high-voltage electric field pulses with a nanosecond duration to modulate cellular activity. This investigation reveals a relatively new and underexplored phenomenon: protein-mediated electroporation. Our research focused on the voltage-sensing domain (VSD) of the NaV1.5 sodium cardiac channel in response to nsPEF stimulation. We scrutinized the VSD structures that form pores and thereby contribute to the physical chemistry that governs the defibrillation effect of nsPEF. To do so, we conducted a comprehensive analysis involving the clustering of 142 replicas simulated for 50 ns under nsPEF stimuli. We subsequently pinpointed the representative structures of each cluster and computed the free energy between them. We find that the selected VSD of NaV1.5 forms pores under nsPEF stimulation, but in a way that significant differs from the traditional VSD opening. This study not only extends our understanding of nsPEF and its interaction with protein channels but also adds a new effect to further study.

Domain- and state-specific shape of the electric field tunes voltage sensing in voltage-gated sodium channels

The ability to sense transmembrane voltage underlies most physiological roles of voltage-gated sodium (Nav) channels. Whereas the key role of their voltage-sensing domains (VSDs) in channel activation is well established, the molecular underpinnings of voltage coupling remain incompletely understood. Voltage-dependent energetics of the activation process can be described in terms of the gating charge that is defined by coupling of charged residues to the external electric field. The shape of the electric field within VSDs is therefore crucial for the activation of voltage-gated ion channels. Here, we employed molecular dynamics simulations of cardiac Nav1.5 and bacterial NavAb, together with our recently developed tool g_elpot, to gain insights into the voltage-sensing mechanisms of Nav channels via high-resolution quantification of VSD electrostatics. In contrast to earlier low-resolution studies, we found that the electric field within VSDs of Nav channels has a complex isoform- and domain-specific shape, which prominently depends on the activation state of a VSD. Different VSDs vary not only in the length of the region where the electric field is focused but also differ in their overall electrostatics, with possible implications in the diverse ion selectivity of their gating pores. Due to state-dependent field reshaping, not only translocated basic but also relatively immobile acidic residues contribute significantly to the gating charge. In the case of NavAb, we found that the transition between structurally resolved activated and resting states results in a gating charge of 8e, which is noticeably lower than experimental estimates. Based on the analysis of VSD electrostatics in the two activation states, we propose that the VSD likely adopts a deeper resting state upon hyperpolarization. In conclusion, our results provide an atomic-level description of the gating charge, demonstrate diversity in VSD electrostatics, and reveal the importance of electric-field reshaping for voltage sensing in Nav channels.

Allosteric Modulation of Membrane Proteins by Small Low-Affinity Ligands

Membrane proteins may respond to a variety of ligands under an applied external stimulus. These ligands include small low-affinity molecules that account for functional effects in the mM range. Understanding the modulation of protein function by low-affinity ligands requires characterizing their atomic-level interactions under dilution, challenging the current resolution of theoretical and experimental routines. Part of the problem derives from the fact that small low-affinity ligands may interact with multiple sites of a membrane protein in a highly degenerate manner to a degree that it is better conceived as a partition phenomenon, hard to track at the molecular interface of the protein. Looking for new developments in the field, we rely on the classic two-state Boltzmann model to devise a novel theoretical description of the allosteric modulation mechanism of membrane proteins in the presence of small low-affinity ligands and external stimuli. Free energy stability of the partition process and its energetic influence on the protein coupling with the external stimulus are quantified. The outcome is a simple formulation that allows the description of the equilibrium shifts of the protein in terms of the grand-canonical partition function of the ligand at dilute concentrations. The model's predictions of the spatial distribution and response probability shift across a variety of ligand concentrations, and thermodynamic conjugates can be directly compared to macroscopic measurements, making it especially useful to interpret experimental data at the atomic level. Illustration and discussion of the theory is shown in the context of general anesthetics and voltage-gated channels for which structural data are available.

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.

Structural mechanisms for the activation of human cardiac KCNQ1 channel by electro-mechanical coupling enhancers

The cardiac KCNQ1 potassium channel carries the important I_Ks current and controls the heart rhythm. Hundreds of mutations in KCNQ1 can cause life-threatening cardiac arrhythmia. Although KCNQ1 structures have been recently resolved, the structural basis for the dynamic electro-mechanical coupling, also known as the voltage sensor domain-pore domain (VSD-PD) coupling, remains largely unknown. In this study, utilizing two VSD-PD coupling enhancers, namely, the membrane lipid phosphatidylinositol 4,5-bisphosphate (PIP₂) and a small-molecule ML277, we determined 2.5-3.5 Å resolution cryo-electron microscopy structures of full-length human KCNQ1-calmodulin (CaM) complex in the apo closed, ML277-bound open, and ML277-PIP₂-bound open states. ML277 binds at the "elbow" pocket above the S4-S5 linker and directly induces an upward movement of the S4-S5 linker and the opening of the activation gate without affecting the C-terminal domain (CTD) of KCNQ1. PIP₂ binds at the cleft between the VSD and the PD and brings a large structural rearrangement of the CTD together with the CaM to activate the PD. These findings not only elucidate the structural basis for the dynamic VSD-PD coupling process during KCNQ1 gating but also pave the way to develop new therapeutics for anti-arrhythmia.

The 70-year search for the voltage sensing mechanism of ion channels

This retrospective on the voltage‐sensing mechanisms and gating models of ion channels begins in 1952 with the charged gating particles postulated by Hodgkin and Huxley, viewed as charges moving across the membrane and controlling its permeability to Na⁺ and K⁺ ions. Hodgkin and Huxley postulated that their movement should generate small and fast capacitive currents, which were recorded 20 years later as gating currents. In the early 1980s, several voltage‐dependent channels were cloned and found to share a common architecture: four homologous domains or subunits, each displaying six transmembrane α‐helical segments, with the fourth segment (S4) displaying four to seven positive charges invariably separated by two non‐charged residues. This immediately suggested that this segment was serving as the voltage sensor of the channel (the molecular counterpart of the charged gating particle postulated by Hodgkin and Huxley) and led to the development of the sliding helix model. Twenty years later, the X‐ray crystallographic structures of many voltage‐dependent channels allowed investigation of their gating by molecular dynamics. Further understanding of how channels gate will benefit greatly from the acquisition of high‐resolution structures of each of their relevant functional or structural states. This will allow the application of molecular dynamics and other approaches. It will also be key to investigate the energetics of channel gating, permitting an understanding of the physical and molecular determinants of gating. The use of multiscale hierarchical approaches might finally prove to be a rewarding strategy to overcome the limits of the various single approaches to the study of channel gating.

Nanosecond Pulsed Electric Field (nsPEF): Opening the Biotechnological Pandora’s Box

Nanosecond Pulsed Electric Field (nsPEF) is an electrostimulation technique first developed in 1995; nsPEF requires the delivery of a series of pulses of high electric fields in the order of nanoseconds into biological tissues or cells. They primary effects in cells is the formation of membrane nanopores and the activation of ionic channels, leading to an incremental increase in cytoplasmic Ca2+ concentration, which triggers a signaling cascade producing a variety of effects: from apoptosis up to cell differentiation and proliferation. Further, nsPEF may affect organelles, making nsPEF a unique tool to manipulate and study cells. This technique is exploited in a broad spectrum of applications, such as: sterilization in the food industry, seed germination, anti-parasitic effects, wound healing, increased immune response, activation of neurons and myocites, cell proliferation, cellular phenotype manipulation, modulation of gene expression, and as a novel cancer treatment. This review thoroughly explores both nsPEF’s history and applications, with emphasis on the cellular effects from a biophysics perspective, highlighting the role of ionic channels as a mechanistic driver of the increase in cytoplasmic Ca2+ concentration.

The pH-dependent gating of the human voltage-gated proton channel from computational simulations

Gating of the voltage-gated proton channel H_V1 is strongly controlled by pH. There is evidence that this involves the sidechains of titratable amino acids that change their protonation state with changes of the pH. Despite experimental investigations to identify the amino acids involved in pH sensing only few progress has been made, including one histidine at the cytoplasmic side of the channel that is involved in sensing cellular pH. We have used constant pH molecular dynamics simulations in symmetrical and asymmetrical pH conditions across the membrane to investigate the pH- and ΔpH-dependent gating of the human H_V1 channel. Therefore, the pK_a of every titratable amino acids has been assessed in single simulations. Our simulations captured initial conformational changes between a deactivated and an activated state of the channel induced solely by changes of the pH. The pH-dependent gating is accompanied by an outward displacement of the three S4 voltage sensing arginines that moves the second arginine past the hydrophobic gasket (HG) which separates the inner and outer pores of the channel. H_V1 activation, when outer pH increases, involves amino acids at the extracellular entrance of the channel that extend the network of interactions from the external solution down to the HG. Whereas, amino acids at the cytoplasmic entrance of the channel are involved in activation, when inner pH decreases, and in a network of interactions that extend from the cytoplasm up to the HG.

Tracking the movement of discrete gating charges in a voltage-gated potassium channel

Positively charged amino acids respond to membrane potential changes to drive voltage sensor movement in voltage-gated ion channels, but determining the displacements of voltage sensor gating charges has proven difficult. We optically tracked the movement of the two most extracellular charged residues (R1 and R2) in the Shaker potassium channel voltage sensor using a fluorescent positively charged bimane derivative (qBBr) that is strongly quenched by tryptophan. By individually mutating residues to tryptophan within the putative pathway of gating charges, we observed that the charge motion during activation is a rotation and a tilted translation that differs between R1 and R2. Tryptophan-induced quenching of qBBr also indicates that a crucial residue of the hydrophobic plug is linked to the Cole-Moore shift through its interaction with R1. Finally, we show that this approach extends to additional voltage-sensing membrane proteins using the Ciona intestinalis voltage-sensitive phosphatase (CiVSP).

Charge substitutions at the voltage-sensing module of domain III enhance actions of site-3 and site-4 toxins on an insect sodium channel

Scorpion α-toxins bind at the pharmacologically-defined site-3 on the sodium channel and inhibit channel inactivation by preventing the outward movement of the voltage sensor in domain IV (IVS4), whereas scorpion β-toxins bind at site-4 on the sodium channel and enhance channel activation by trapping the voltage sensor of domain II (IIS4) in its outward position. However, limited information is available on the role of the voltage-sensing modules (VSM, comprising S1-S4) of domains I and III in toxin actions. We have previously shown that charge reversing substitutions of the innermost positively-charged residues in IIIS4 (R4E, R5E) increase the activity of an insect-selective site-4 scorpion toxin, Lqh-dprIT3-c, on BgNav1-1a, a cockroach sodium channel. Here we show that substitutions R4E and R5E in IIIS4 also increase the activity of two site-3 toxins, LqhαIT from Leiurusquinquestriatus hebraeus and insect-selective Av3 from Anemonia viridis. Furthermore, charge reversal of either of two conserved negatively-charged residues, D1K and E2K, in IIIS2 also increase the action of the site-3 and site-4 toxins. Homology modeling suggests that S2-D1 and S2-E2 interact with S4-R4 and S4-R5 in the VSM of domain III (III-VSM), respectively, in the activated state of the channel. However, charge swapping between S2-D1 and S4-R4 had no compensatory effects on gating or toxin actions, suggesting that charged residue interactions are complex. Collectively, our results highlight the involvement of III-VSM in the actions of both site 3 and site 4 toxins, suggesting that charge reversing substitutions in III-VSM allosterically facilitate IIS4 or IVS4 voltage sensor trapping by these toxins.

Gating current noise produced by Brownian models of a voltage sensor

The activation of voltage-dependent ion channels is associated with the movement of gating charges, which give rise to gating currents. Although gating currents from a single channel are too small to be detected, analysis of the fluctuations of macroscopic gating currents from a population of channels allows a good guess of their magnitude. The analysis of experimental gating current fluctuations, when interpreted in terms of a rate model of channel activation and assuming sufficiently high bandwidth, is in accordance with the presence of a main step along the activation pathway carrying a charge of 2.3-2.4 e₀. To give a physical interpretation to these results and to relate them to the known atomic structure of the voltage sensor domain, we used a Brownian model of voltage-dependent gating based on atomic detail structure, that follows the laws of electrodynamics. The model predicts gating currents and gating current fluctuations essentially similar to those experimentally observed. The detailed study of the model output, also performed by making several simplifications aimed at understanding the basic dependencies of the gating current fluctuations, suggests that in real channels the voltage sensor moves along a sequence of intermediate states separated by relatively low (<5 kT) energy barriers. As a consequence, crossings of successive gating charges through the gating pore become very frequent, and the corresponding current shots are often seen to overlap because of the relatively high filtering. Notably, this limited bandwidth effect is at the origin of the relatively high single-step charge experimentally detected.

Exploring the Conformational Changes Induced by Nanosecond Pulsed Electric Fields on the Voltage Sensing Domain of a Ca2+ Channel

Nanosecond Pulsed Electric Field (nsPEF or Nano Pulsed Stimulation, NPS) is a technology that delivers a series of pulses of high-voltage electric fields during a short period of time, in the order of nanoseconds. The main consequence of nsPEF upon cells is the formation of nanopores, which is followed by the gating of ionic channels. Literature is conclusive in that the physiological mechanisms governing ion channel gating occur in the order of milliseconds. Hence, understanding how these channels can be activated by a nsPEF would be an important step in order to conciliate fundamental biophysical knowledge with improved nsPEF applications. To get insights on both the kinetics and thermodynamics of ion channel gating induced by nsPEF, in this work, we simulated the Voltage Sensing Domain (VSD) of a voltage-gated Ca2+ channel, inserted in phospholipidic membranes with different concentrations of cholesterol. We studied the conformational changes of the VSD under a nsPEF mimicked by the application of a continuous electric field lasting 50 ns with different intensities as an approach to reveal novel mechanisms leading to ion channel gating in such short timescales. Our results show that using a membrane with high cholesterol content, under an nsPEF of 50 ns and E→ = 0.2 V/nm, the VSD undergoes major conformational changes. As a whole, our work supports the notion that membrane composition may act as an allosteric regulator, specifically cholesterol content, which is fundamental for the response of the VSD to an external electric field. Moreover, changes on the VSD structure suggest that the gating of voltage-gated Ca2+ channels by a nsPEF may be due to major conformational changes elicited in response to the external electric field. Finally, the VSD/cholesterol-bilayer under an nsPEF of 50 ns and E→ = 0.2 V/nm elicits a pore formation across the VSD suggesting a new non-reported effect of nsPEF into cells, which can be called a “protein mediated electroporation”.

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.

Insights into Cardiac IKs (KCNQ1/KCNE1) Channels Regulation

The delayed rectifier potassium IKs channel is an important regulator of the duration of the ventricular action potential. Hundreds of mutations in the genes (KCNQ1 and KCNE1) encoding the IKs channel cause long QT syndrome (LQTS). LQTS is a heart disorder that can lead to severe cardiac arrhythmias and sudden cardiac death. A better understanding of the IKs channel (here called the KCNQ1/KCNE1 channel) properties and activities is of great importance to find the causes of LQTS and thus potentially treat LQTS. The KCNQ1/KCNE1 channel belongs to the superfamily of voltage-gated potassium channels. The KCNQ1/KCNE1 channel consists of both the pore-forming subunit KCNQ1 and the modulatory subunit KCNE1. KCNE1 regulates the function of the KCNQ1 channel in several ways. This review aims to describe the current structural and functional knowledge about the cardiac KCNQ1/KCNE1 channel. In addition, we focus on the modulation of the KCNQ1/KCNE1 channel and its potential as a target therapeutic of LQTS.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

The interplay of excitation and electroporation in nanosecond pulse stimulation

Conventional electric stimuli of micro- and millisecond duration excite or activate cells at voltages 10-100 times below the electroporation threshold. This ratio is remarkably different for nanosecond electric pulses (nsEP), which caused excitation and activation only at or above the electroporation threshold in diverse cell lines, primary cardiomyocytes, neurons, and chromaffin cells. Depolarization to the excitation threshold often results from (or is assisted by) the loss of the resting membrane potential due to ion leaks across the membrane permeabilized by nsEP. Slow membrane resealing and the build-up of electroporation damages prevent repetitive excitation by nsEP. However, peripheral nerves and muscles are exempt from this rule and withstand multiple cycles of excitation by nsEP without the loss of function or signs of electroporation. We show that the damage-free excitation by nsEP may be enabled by the membrane charging time constant sufficiently large to (1) cap the peak transmembrane voltage during nsEP below the electroporation threshold, and (2) extend the post-nsEP depolarization long enough to activate voltage-gated ion channels. The low excitatory efficacy of nsEP compared to longer pulses makes them advantageous for medical applications where the neuromuscular excitation is an unwanted side effect, such as electroporation-based cancer and tissue ablation.

Aromatic interactions with membrane modulate human BK channel activation

Large-conductance potassium (BK) channels are transmembrane (TM) proteins that can be synergistically and independently activated by membrane voltage and intracellular Ca²⁺. The only covalent connection between the cytosolic Ca²⁺ sensing domain and the TM pore and voltage sensing domains is a 15-residue ‘C-linker’. To determine the linker’s role in human BK activation, we designed a series of linker sequence scrambling mutants to suppress potential complex interplay of specific interactions with the rest of the protein. The results revealed a surprising sensitivity of BK activation to the linker sequence. Combining atomistic simulations and further mutagenesis experiments, we demonstrated that nonspecific interactions of the linker with membrane alone could directly modulate BK activation. The C-linker thus plays more direct roles in mediating allosteric coupling between BK domains than previously assumed. Our results suggest that covalent linkers could directly modulate TM protein function and should be considered an integral component of the sensing apparatus.

Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated Ion Channels

Pulsed electric fields are increasingly used in medicine to transiently increase the cell membrane permeability via electroporation to deliver therapeutic molecules into the cell. One type of event that contributes to this increase in membrane permeability is the formation of pores in the membrane lipid bilayer. However, electrophysiological measurements suggest that membrane proteins are affected as well, particularly voltage-gated ion channels (VGICs). The molecular mechanisms by which the electric field could affects these molecules remain unidentified. In this study, we used molecular dynamics simulations to unravel the molecular events that take place in different VGICs when exposing them to electric fields mimicking electroporation conditions. We show that electric fields can induce pores in the voltage-sensor domains (VSDs) of different VGICs and that these pores form more easily in some channels than in others. We demonstrate that poration is more likely in VSDs that are more hydrated and are electrostatically more favorable for the entry of ions. We further show that pores in VSDs can expand into so-called complex pores, which become stabilized by lipid headgroups. Our results suggest that such complex pores are considerably more stable than conventional lipid pores, and their formation can lead to severe unfolding of VSDs from the channel. We anticipate that such VSDs become dysfunctional and unable to respond to changes in transmembrane voltage, which is in agreement with previous electrophysiological measurements showing a decrease in the voltage-dependent transmembrane ionic currents after pulse treatment. Finally, we discuss the possibility of activation of VGICs by submicrosecond-duration pulses. Overall, our study reveals a new, to our knowledge, mechanism of electroporation through membranes containing VGICs.

Roles for Countercharge in the Voltage Sensor Domain of Ion Channels

Voltage-gated ion channels share a common structure typified by peripheral, voltage sensor domains. Their S4 segments respond to alteration in membrane potential with translocation coupled to ion permeation through a central pore domain. The mechanisms of gating in these channels have been intensely studied using pioneering methods such as measurement of charge displacement across a membrane, sequencing of genes coding for voltage-gated ion channels, and the development of all-atom molecular dynamics simulations using structural information from prokaryotic and eukaryotic channel proteins. One aspect of this work has been the description of the role of conserved negative countercharges in S1, S2, and S3 transmembrane segments to promote sequential salt-bridge formation with positively charged residues in S4 segments. These interactions facilitate S4 translocation through the lipid bilayer. In this review, we describe functional and computational work investigating the role of these countercharges in S4 translocation, voltage sensor domain hydration, and in diseases resulting from countercharge mutations.

Molecular Insights from Conformational Ensembles via Machine Learning

Biomolecular simulations are intrinsically high dimensional and generate noisy data sets of ever-increasing size. Extracting important features from the data is crucial for understanding the biophysical properties of molecular processes, but remains a big challenge. Machine learning (ML) provides powerful dimensionality reduction tools. However, such methods are often criticized as resembling black boxes with limited human-interpretable insight. We use methods from supervised and unsupervised ML to efficiently create interpretable maps of important features from molecular simulations. We benchmark the performance of several methods, including neural networks, random forests, and principal component analysis, using a toy model with properties reminiscent of macromolecular behavior. We then analyze three diverse biological processes: conformational changes within the soluble protein calmodulin, ligand binding to a G protein-coupled receptor, and activation of an ion channel voltage-sensor domain, unraveling features critical for signal transduction, ligand binding, and voltage sensing. This work demonstrates the usefulness of ML in understanding biomolecular states and demystifying complex simulations.

Two-stage electro-mechanical coupling of a KV channel in voltage-dependent activation

In voltage-gated potassium (K_V) channels, the voltage-sensing domain (VSD) undergoes sequential activation from the resting state to the intermediate state and activated state to trigger pore opening via electro-mechanical (E-M) coupling. However, the spatial and temporal details underlying E-M coupling remain elusive. Here, utilizing K_V7.1's unique two open states, we report a two-stage E-M coupling mechanism in voltage-dependent gating of K_V7.1 as triggered by VSD activations to the intermediate and then activated state. When the S4 segment transitions to the intermediate state, the hand-like C-terminus of the VSD-pore linker (S4-S5L) interacts with the pore in the same subunit. When S4 then proceeds to the fully-activated state, the elbow-like hinge between S4 and S4-S5L engages with the pore of the neighboring subunit to activate conductance. This two-stage hand-and-elbow gating mechanism elucidates distinct tissue-specific modulations, pharmacology, and disease pathogenesis of K_V7.1, and likely applies to numerous domain-swapped K_V channels.

Voltage-dependent gating in K channels: experimental results and quantitative models

Voltage-dependent K channels open and close in response to voltage changes across the cell membrane. This voltage dependence was postulated to depend on the presence of charged particles moving through the membrane in response to voltage changes. Recording of gating currents originating from the movement of these particles fully confirmed this hypothesis, and gave substantial experimental clues useful for the detailed understanding of the process. In the absence of structural information, the voltage-dependent gating was initially investigated using discrete Markov models, an approach only capable of providing a kinetic and thermodynamic comprehension of the process. The elucidation of the crystal structure of the first voltage-dependent channel brought in a dramatic change of pace in the understanding of channel gating, and in modeling the underlying processes. It was now possible to construct quantitative models using molecular dynamics, where all the interactions of each individual atom with the surroundings were taken into account, and its motion predicted by Newton's laws. Unfortunately, this modeling is computationally very demanding, and in spite of the advances in simulation procedures and computer technology, it is still limited in its predictive ability. To overcome these limitations, several groups have developed more macroscopic voltage gating models. Their approaches understandably require a number of approximations, which must however be physically well justified. One of these models, based on the description of the voltage sensor as a Brownian particle, that we have recently developed, is able to simultaneously describe the behavior of a single voltage sensor and to predict the macroscopic gating current originating from a population of sensors. The basics of this model are here described, and a typical application using the Kv1.2/2.1 chimera channel structure is also presented.

Helix breaking transition in the S4 of HCN channel is critical for hyperpolarization-dependent gating

In contrast to most voltage-gated ion channels, hyperpolarization- and cAMP gated (HCN) ion channels open on hyperpolarization. Structure-function studies show that the voltage-sensor of HCN channels are unique but the mechanisms that determine gating polarity remain poorly understood. All-atom molecular dynamics simulations (~20 μs) of HCN1 channel under hyperpolarization reveals an initial downward movement of the S4 voltage-sensor but following the transfer of last gating charge, the S4 breaks into two sub-helices with the lower sub-helix becoming parallel to the membrane. Functional studies on bipolar channels show that the gating polarity strongly correlates with helical turn propensity of the substituents at the breakpoint. Remarkably, in a proto-HCN background, the replacement of breakpoint serine with a bulky hydrophobic amino acid is sufficient to completely flip the gating polarity from inward to outward-rectifying. Our studies reveal an unexpected mechanism of inward rectification involving a linker sub-helix emerging from HCN S4 during hyperpolarization.

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.

Nuclear Magnetic Resonance Solution Structure and Functional Behavior of the Human Proton Channel

The human voltage-gated proton channel [Hv1⁽¹⁾ or VSDO⁽²⁾] plays an important role in the human innate immune system. Its structure differs considerably from those of other cation channels. It is built solely of a voltage-sensing domain and thus lacks the central pore domain, which is essential for other cation channels. Here, we determined the solution structure of an N- and C-terminally truncated human Hv1 (Δ-Hv1) in the resting state by nuclear magnetic resonance (NMR) spectroscopy. Δ-Hv1 comprises the typical voltage-sensing antiparallel four-helix bundle (S1-S4) preceded by an amphipathic helix (S0). The solution structure corresponds to an intermediate state between resting and activated forms of voltage-sensing domains. Furthermore, Zn²⁺-induced closing of proton channel Δ-Hv1 was studied with two-dimensional NMR spectroscopy, which showed that characteristic large scale dynamics of open Δ-Hv1 are absent in the closed state of the channel. Additionally, pH titration studies demonstrated that a higher H⁺ concentration is required for the protonation of side chains in the Zn²⁺-induced closed state than in the open state. These observations demonstrate both structural and dynamical changes involved in the process of voltage gating of the Hv1 channel and, in the future, may help to explain the unique properties of unidirectional conductance and the exceptional ion selectivity of the channel.

Upgraded molecular models of the human KCNQ1 potassium channel

The voltage-gated potassium channel KCNQ1 (K_V7.1) assembles with the KCNE1 accessory protein to generate the slow delayed rectifier current, I_KS, which is critical for membrane repolarization as part of the cardiac action potential. Loss-of-function (LOF) mutations in KCNQ1 are the most common cause of congenital long QT syndrome (LQTS), type 1 LQTS, an inherited genetic predisposition to cardiac arrhythmia and sudden cardiac death. A detailed structural understanding of KCNQ1 is needed to elucidate the molecular basis for KCNQ1 LOF in disease and to enable structure-guided design of new anti-arrhythmic drugs. In this work, advanced structural models of human KCNQ1 in the resting/closed and activated/open states were developed by Rosetta homology modeling guided by newly available experimentally-based templates: X. leavis KCNQ1 and various resting voltage sensor structures. Using molecular dynamics (MD) simulations, the capacity of the models to describe experimentally established channel properties including state-dependent voltage sensor gating charge interactions and pore conformations, PIP2 binding sites, and voltage sensor–pore domain interactions were validated. Rosetta energy calculations were applied to assess the utility of each model in interpreting mutation-evoked KCNQ1 dysfunction by predicting the change in protein thermodynamic stability for 50 experimentally characterized KCNQ1 variants with mutations located in the voltage-sensing domain. Energetic destabilization was successfully predicted for folding-defective KCNQ1 LOF mutants whereas wild type-like mutants exhibited no significant energetic frustrations, which supports growing evidence that mutation-induced protein destabilization is an especially common cause of KCNQ1 dysfunction. The new KCNQ1 Rosetta models provide helpful tools in the study of the structural basis for KCNQ1 function and can be used to generate hypotheses to explain KCNQ1 dysfunction.

Conduction through a narrow inward-rectifier K+ channel pore

G-protein–gated inwardly rectifying potassium channels are important mediators of inhibitory neurotransmission. Based on microsecond-scale molecular dynamics simulations, Bernsteiner et al. propose novel gating details that may enable K⁺ flux via a direct knock-on mechanism.

Inwardly rectifying potassium (Kir) channels play a key role in controlling membrane potentials in excitable and unexcitable cells, thereby regulating a plethora of physiological processes. G-protein–gated Kir channels control heart rate and neuronal excitability via small hyperpolarizing outward K⁺ currents near the resting membrane potential. Despite recent breakthroughs in x-ray crystallography and cryo-EM, the gating and conduction mechanisms of these channels are poorly understood. MD simulations have provided unprecedented details concerning the gating and conduction mechanisms of voltage-gated K⁺ and Na⁺ channels. Here, we use multi-microsecond–timescale MD simulations based on the crystal structures of GIRK2 (Kir3.2) bound to phosphatidylinositol-4,5-bisphosphate to provide detailed insights into the channel’s gating dynamics, including insights into the behavior of the G-loop gate. The simulations also elucidate the elementary steps that underlie the movement of K⁺ ions through an inward-rectifier K⁺ channel under an applied electric field. Our simulations suggest that K⁺ permeation might occur via direct knock-on, similar to the mechanism recently shown for K_v channels.

Voltage-Dependent Profile Structures of a Kv-Channel via Time-Resolved Neutron Interferometry

Available experimental techniques cannot determine high-resolution three-dimensional structures of membrane proteins under a transmembrane voltage. Hence, the mechanism by which voltage-gated cation channels couple conformational changes within the four voltage sensor domains, in response to either depolarizing or polarizing transmembrane voltages, to opening or closing of the pore domain's ion channel remains unresolved. Single-membrane specimens, composed of a phospholipid bilayer containing a vectorially oriented voltage-gated K+ channel protein at high in-plane density tethered to the surface of an inorganic multilayer substrate, were developed to allow the application of transmembrane voltages in an electrochemical cell. Time-resolved neutron reflectivity experiments, enhanced by interferometry enabled by the multilayer substrate, were employed to provide directly the low-resolution profile structures of the membrane containing the vectorially oriented voltage-gated K+ channel for the activated, open and deactivated, closed states of the channel under depolarizing and hyperpolarizing transmembrane voltages applied cyclically. The profile structures of these single membranes were dominated by the voltage-gated K+ channel protein because of the high in-plane density. Importantly, the use of neutrons allowed the determination of the voltage-dependent changes in both the profile structure of the membrane and the distribution of water within the profile structure. These two key experimental results were then compared to those predicted by three computational modeling approaches for the activated, open and deactivated, closed states of three different voltage-gated K+ channels in hydrated phospholipid bilayer membrane environments. Of the three modeling approaches investigated, only one state-of-the-art molecular dynamics simulation that directly predicted the response of a voltage-gated K+ channel within a phospholipid bilayer membrane to applied transmembrane voltages by utilizing very long trajectories was found to be in agreement with the two key experimental results provided by the time-resolved neutron interferometry experiments.

ML277 specifically enhances the fully activated open state of KCNQ1 by modulating VSD-pore coupling

Upon membrane depolarization, the KCNQ1 potassium channel opens at the intermediate (IO) and activated (AO) states of the stepwise voltage-sensing domain (VSD) activation. In the heart, KCNQ1 associates with KCNE1 subunits to form IKs channels that regulate heart rhythm. KCNE1 suppresses the IO state so that the IKs channel opens only to the AO state. Here, we tested modulations of human KCNQ1 channels by an activator ML277 in Xenopus oocytes. It exclusively changes the pore opening properties of the AO state without altering the IO state, but does not affect VSD activation. These observations support a distinctive mechanism responsible for the VSD-pore coupling at the AO state that is sensitive to ML277 modulation. ML277 provides insights and a tool to investigate the gating mechanism of KCNQ1 channels, and our study reveals a new strategy for treating long QT syndrome by specifically enhancing the AO state of native IKs currents.

Tubulin response to intense nanosecond-scale electric field in molecular dynamics simulation

Intense pulsed electric fields are known to act at the cell membrane level and are already being exploited in biomedical and biotechnological applications. However, it is not clear if electric pulses within biomedically-attainable parameters could directly influence intra-cellular components such as cytoskeletal proteins. If so, a molecular mechanism of action could be uncovered for therapeutic applications of such electric fields. To help clarify this question, we first identified that a tubulin heterodimer is a natural biological target for intense electric fields due to its exceptional electric properties and crucial roles played in cell division. Using molecular dynamics simulations, we then demonstrated that an intense - yet experimentally attainable - electric field of nanosecond duration can affect the bβ-tubulin’s C-terminus conformations and also influence local electrostatic properties at the GTPase as well as the binding sites of major tubulin drugs site. Our results suggest that intense nanosecond electric pulses could be used for physical modulation of microtubule dynamics. Since a nanosecond pulsed electric field can penetrate the tissues and cellular membranes due to its broadband spectrum, our results are also potentially significant for the development of new therapeutic protocols.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Molecular mechanism of phosphoinositides' specificity for the inwardly rectifying potassium channel Kir2.2

We investigated the binding mechanism of PI(4,5)P₂ and variants on the inwardly rectifying potassium channel, Kir2.2. Our results not only demonstrated the molecular origin for their binding specificity, but also revealed the major driving forces.

Phosphoinositides are essential signaling lipids that play a critical role in regulating ion channels, and their dysregulation often results in fatal diseases including cardiac arrhythmia and paralysis. Despite decades of intensive research, the underlying molecular mechanism of lipid agonism and specificity remains largely unknown. Here, we present a systematic study of the binding mechanism and specificity of a native agonist, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) and two of its variants, PI(3,4)P₂ and PI(3,4,5)P₃, on inwardly rectifying potassium channel Kir2.2, using molecular dynamics simulations and free energy perturbations (FEPs). Our results demonstrate that the major driving force for the PI(4,5)P₂ specificity on Kir2.2 comes from the highly organized salt-bridge network formed between the charged inositol head and phosphodiester linker of PI(4,5)P₂. The unsaturated arachidonic chain is also shown to contribute to the stable binding through hydrophobic interactions with nearby Kir2.2 hydrophobic residues. Consistent with previous experimental findings, our FEP results confirmed that non-native ligands, PI(3,4)P₂ and PI(3,4,5)P₃, show significant loss in binding affinity as a result of the substantial shift from the native binding mode and unfavorable local solvation environment. However, surprisingly, the underlying molecular pictures for the unfavorable binding of both ligands are quite distinctive: for PI(3,4)P₂, it is due to a direct destabilization in the bound state, whereas for PI(3,4,5)P₃, it is due to a relative stabilization in its free state. Our findings not only provide a theoretical basis for the ligand specificity, but also generate new insights into the allosteric modulation of ligand-gated ion channels.

Binding of the general anesthetic sevoflurane to ion channels

The direct-site hypothesis assumes general anesthetics bind ion channels to impact protein equilibrium and function, inducing anesthesia. Despite advancements in the field, a first principle all-atom demonstration of this structure-function premise is still missing. We focus on the clinically used sevoflurane interaction to anesthetic-sensitive Kv1.2 mammalian channel to resolve if sevoflurane binds protein’s well-characterized open and closed structures in a conformation-dependent manner to shift channel equilibrium. We employ an innovative approach relying on extensive docking calculations and free-energy perturbation of all potential binding sites revealed by the latter, and find sevoflurane binds open and closed structures at multiple sites under complex saturation and concentration effects. Results point to a non-trivial interplay of site and conformation-dependent modes of action involving distinct binding sites that increase channel open-probability at diluted ligand concentrations. Given the challenge in exploring more complex processes potentially impacting channel-anesthetic interaction, the result is revealing as it demonstrates the process of multiple anesthetic binding events alone may account for open-probability shifts recorded in measurements.

General anesthetics are central to modern medicine, yet their microscopic mechanism of action is still unknown. Here, we demonstrate that a clinically used anesthetic, sevoflurane, binds the mammalian voltage-gated potassium channel Kv1.2 effecting a shift in its open probability, even at low concentrations. The results, supported by recent experimental measurements, are promising as they demonstrate that the molecular process of direct binding of anesthetic to ion channels play a relevant role in anesthesia.

A New Cardiac Channelopathy: From Clinical Phenotypes to Molecular Mechanisms Associated With Nav1.5 Gating Pores

Voltage gated sodium channels (Na_V) are broadly expressed in the human body. They are responsible for the initiation of action potentials in excitable cells. They also underlie several physiological processes such as cognitive, sensitive, motor, and cardiac functions. The Na_V1.5 channel is the main Na_V expressed in the heart. A dysfunction of this channel is usually associated with the development of pure electrical disorders such as long QT syndrome, Brugada syndrome, sinus node dysfunction, atrial fibrillation, and cardiac conduction disorders. However, mutations of Na_v1.5 have recently been linked to the development of an atypical clinical entity combining complex arrhythmias and dilated cardiomyopathy. Although several Na_v1.5 mutations have been linked to dilated cardiomyopathy phenotypes, their pathogenic mechanisms remain to be elucidated. The gating pore may constitute a common biophysical defect for all Na_V1.5 mutations located in the channel's VSDs. The creation of such a gating pore may disrupt the ionic homeostasis of cardiomyocytes, affecting electrical signals, cell morphology, and cardiac myocyte function. The main objective of this article is to review the concept of gating pores and their role in structural heart diseases and to discuss potential pharmacological treatments.

Molecular Dynamics Simulation for All

The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K_v1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of −70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Determining the molecular basis of voltage sensitivity in membrane proteins

The identification of voltage-sensing elements in membrane proteins is challenging due to the diversity of voltage-sensing mechanisms. Kasimova et al. present a computational approach to predict the elements involved in voltage sensing, which they validate using voltage-gated ion channels.

Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mechanical work. They are responsible for a spectrum of physiological processes in living organisms, including electrical signaling and cell-cycle progression. Although the mechanism of voltage-sensing has been well characterized for some membrane proteins, including voltage-gated ion channels, even the location of the voltage-sensing elements remains unknown for others. Moreover, the detection of these elements by using experimental techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estimation of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis experiments will significantly reduce the number of potential voltage-sensitive elements to be tested.

Gating currents

A voltage change across a membrane protein moves charges or dipoles producing a gating current that is an electrical expression of a conformational change.

Many membrane proteins sense the voltage across the membrane where they are inserted, and their function is affected by voltage changes. The voltage sensor consists of charges or dipoles that move in response to changes in the electric field, and their movement produces an electric current that has been called gating current. In the case of voltage-gated ion channels, the kinetic and steady-state properties of the gating charges provide information of conformational changes between closed states that are not visible when observing ionic currents only. In this Journal of General Physiology Milestone, the basic principles of voltage sensing and gating currents are presented, followed by a historical description of the recording of gating currents. The results of gating current recordings are then discussed in the context of structural changes in voltage-dependent membrane proteins and how these studies have provided new insights on gating mechanisms.