Voltage Sensing in Membranes: From Macroscopic Currents to Molecular Motions

Hydrophobic residues in S1 modulate enzymatic function and voltage sensing in voltage-sensing phosphatase

The voltage-sensing domain (VSD) is a four-helix modular protein domain that converts electrical signals into conformational changes, leading to open pores and active enzymes. In most voltage-sensing proteins, the VSDs do not interact with one another, and the S1-S3 helices are considered mainly scaffolding, except in the voltage-sensing phosphatase (VSP) and the proton channel (Hv). To investigate its contribution to VSP function, we mutated four hydrophobic amino acids in S1 to alanine (F127, I131, I134, and L137), individually or in combination. Most of these mutations shifted the voltage dependence of activity to higher voltages; however, not all substrate reactions were the same. The kinetics of enzymatic activity were also altered, with some mutations significantly slowing down dephosphorylation. The voltage dependence of VSD motions was consistently shifted to lower voltages and indicated a second voltage-dependent motion. Additionally, none of the mutations broke the VSP dimer, indicating that the S1 impact could stem from intra- and/or intersubunit interactions. Lastly, when the same mutations were introduced into a genetically encoded voltage indicator, they dramatically altered the optical readings, making some of the kinetics faster and shifting the voltage dependence. These results indicate that the S1 helix in VSP plays a critical role in tuning the enzyme's conformational response to membrane potential transients and influencing the function of the VSD.

The Pentameric Ligand-Gated Ion Channel Family: A New Member of the Voltage Gated Ion Channel Superfamily?

In this report we present seven lines of bioinformatic evidence supporting the conclusion that the Pentameric Ligand-gated Ion Channel (pLIC) Family is a member of the Voltage-gated Ion Channel (VIC) Superfamily. In our approach, we used the Transporter Classification Database (TCDB) as a reference and applied a series of bioinformatic methods to search for similarities between the pLIC family and members of the VIC superfamily. These include: (1) sequence similarity, (2) compatibility of topology and hydropathy profiles, (3) shared domains, (4) conserved motifs, (5) similarity of Hidden Markov Model profiles between families, (6) common 3D structural folds, and (7) clustering analysis of all families. Furthermore, sequence and structural comparisons as well as the identification of a 3-TMS repeat unit in the VIC superfamily suggests that the sixth transmembrane segment evolved into a re-entrant loop. This evidence suggests that the voltage-sensor domain and the channel domain have a common origin. The classification of the pLIC family within the VIC superfamily sheds light onto the topological origins of this family and its evolution, which will facilitate experimental verification and further research into this superfamily by the scientific community.

S1 hydrophobic residues modulate voltage sensing phosphatase enzymatic function and voltage sensing

The voltage sensing domain (VSD) is a four-helix modular protein domain that converts electrical signals into conformational changes, leading to open pores and active enzymes. In most voltage sensing proteins, the VSDs do not interact with one another and the S1-S3 helices are considered mainly as scaffolding. The two exceptions are the voltage sensing phosphatase (VSP) and the proton channel (Hv). VSP is a voltage-regulated enzyme and Hvs are channels that only have VSDs. To investigate the S1 contribution to VSP function, we individually mutated four hydrophobic amino acids in S1 to alanine (F127, I131, I134 and L137). We also combined these mutations to generate quadruple mutation designated S1-Q. Most of these mutations shifted the voltage dependence of activity to higher voltages though interestingly, not all substrate reactions were the same. The kinetics of enzymatic activity were also altered with some mutations significantly slowing down dephosphorylation. The voltage dependence of VSD motions were consistently shifted to lower voltages and indicated a second voltage dependent motion. Co-immunoprecipitation demonstrated that none of the mutations broke the VSP dimer indicating that the S1 impact could stem from intrasubunit and/or intersubunit interactions. Lastly, when the same alanine mutations were introduced into a genetically encoded voltage indicator, they dramatically altered the optical readings, making some of the kinetics faster and shifting the voltage dependence. These results indicate that the S1 helix in VSP plays a critical role in tuning the enzymes conformational response to membrane potential transients and influencing the function of the VSD.

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.

A novel de novo KCNB1 variant altering channel characteristics in a patient with periventricular heterotopia, abnormal corpus callosum, and mild seizure outcome

KCNB1 encodes the α-subunit of Kv2.1, the main contributor to neuronal delayed rectifier potassium currents. The subunit consists of six transmembrane α helices (S1-S6), comprising the voltage-sensing domain (S1-S4) and the pore domain (S5-P-S6). Heterozygous KCNB1 pathogenic variants are associated with developmental and epileptic encephalopathy. Here we report an individual who shows the milder phenotype compared to the previously reported cases, including delayed language development, mild intellectual disability, attention deficit hyperactivity disorder, late-onset epilepsy responsive to an antiepileptic drug, elevation of serum creatine kinase, and peripheral axonal neuropathy. On the other hand, his brain MRI showed characteristic findings including periventricular heterotopia, polymicrogyria, and abnormal corpus callosum. Exome sequencing identified a novel de novo KCNB1 variant c.574G>A, p.(Ala192Thr) located in the S1 segment of the voltage-sensing domain. Functional analysis using the whole-cell patch-clamp technique in Neuro2a cells showed that the Ala192Thr mutant reduces both activation and inactivation of the channel at membrane voltages in the range of -50 to -30 mV. Our case could expand the phenotypic spectrum of patients with KCNB1 variants, and suggested that variants located in the S1 segment might be associated with a milder outcome of seizures.

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.

Cumulative hydropathic topology of a voltage-gated sodium channel at atomic resolution

Voltage-gated sodium channels (NavChs) are biological pores that control the flow of sodium ions through the cell membrane. In humans, mutations in genes encoding NavChs can disrupt physiological cellular activity thus leading to a wide spectrum of diseases. Here, we present a topological connection between the functional architecture of a NavAb bacterial channel and accumulation of atomic hydropathicity around its pore. This connection is established via a scaling analysis methodology that elucidates how intrachannel hydropathic density variations translate into hydropathic dipole field configurations along the pore. Our findings suggest the existence of a nonrandom cumulative hydropathic topology that is organized parallel to the membrane surface so that pore's stability, as well as, gating behavior are guaranteed. Given the biophysical significance of the hydropathic effect, our study seeks to provide a computational framework for studying cumulative hydropathic topological properties of NavChs and pore-forming proteins in general.

Voltage-dependent structural models of the human Hv1 proton channel from long-timescale molecular dynamics simulations

The voltage-gated Hv1 proton channel is a ubiquitous membrane protein that has roles in a variety of cellular processes, including proton extrusion, pH regulation, production of reactive oxygen species, proliferation of cancer cells, and increased brain damage during ischemic stroke. A crystal structure of an Hv1 construct in a putative closed state has been reported, and structural models for the channel open state have been proposed, but a complete characterization of the Hv1 conformational dynamics under an applied membrane potential has been elusive. We report structural models of the Hv1 voltage-sensing domain (VSD), both in a hyperpolarized state and a depolarized state resulting from voltage-dependent conformational changes during a 10-μs-timescale atomistic molecular dynamics simulation in an explicit membrane environment. In response to a depolarizing membrane potential, the S4 helix undergoes an outward displacement, leading to changes in the VSD internal salt-bridge network, resulting in a reshaping of the permeation pathway and a significant increase in hydrogen bond connectivity throughout the channel. The total gating charge displacement associated with this transition is consistent with experimental estimates. Molecular docking calculations confirm the proposed mechanism for the inhibitory action of 2-guanidinobenzimidazole (2GBI) derived from electrophysiological measurements and mutagenesis. The depolarized structural model is also consistent with the formation of a metal bridge between residues located in the core of the VSD. Taken together, our results suggest that these structural models are representative of the closed and open states of the Hv1 channel.

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.

Cancellation of nerve excitation by the reversal of nanosecond stimulus polarity and its relevance to the gating time of sodium channels

The initiation of action potentials (APs) by membrane depolarization occurs after a brief vulnerability period, during which excitation can be abolished by the reversal of the stimulus polarity. This vulnerability period is determined by the time needed for gating of voltage-gated sodium channels (VGSC). We compared nerve excitation by ultra-short uni- and bipolar stimuli to define the time frame of bipolar cancellation and of AP initiation. Propagating APs in isolated frog sciatic nerve were elicited by cathodic pulses (200 ns-300 µs), followed by an anodic (canceling) pulse of the same duration after a 0-200-µs delay. We found that the earliest and the latest boundaries for opening the critical number of VGSC needed to initiate AP are, respectively, between 11 and 20 µs and between 100 and 200 µs after the onset of depolarization. Stronger depolarization accelerated AP initiation, apparently due to faster VGSC opening, but not beyond the 11-µs limit. Bipolar cancellation was augmented by reducing pulse duration, shortening the delay between pulses, decreasing the amplitude of the cathodic pulse, and increasing the amplitude of the anodic one. Some of these characteristics contrasted the bipolar cancellation of cell membrane electroporation (Pakhomov et al. in Bioelectrochemistry 122:123-133, 2018; Gianulis et al. in Bioelectrochemistry 119:10-19, 2017), suggesting different mechanisms. The ratio of nerve excitation thresholds for a unipolar cathodic pulse and a symmetrical bipolar pulse increased as a power function as the pulse duration decreased, in remarkable agreement with the predictions of SENN model of nerve excitation (Reilly and Diamant in Health Phys 83(3):356-365, 2002).

Influence of Lipid Saturation, Hydrophobic Length and Cholesterol on Double-Arginine-Containing Helical Peptides in Bilayer Membranes

Membrane proteins are essential for many cell processes yet are more difficult to investigate than soluble proteins. Charged residues often contribute significantly to membrane protein function. Model peptides such as GWALP23 (acetyl-GGALW⁵ LAL⁸ LALALAL¹⁶ ALW¹⁹ LAGA-amide) can be used to characterize the influence of specific residues on transmembrane protein domains. We have substituted R8 and R16 in GWALP23 in place of L8 and L16, equidistant from the peptide center, and incorporated specific ² H-labeled alanine residues within the central sequence for detection by solid-state ² H NMR spectroscopy. The resulting pattern of [² H]Ala quadrupolar splitting (Δν_q ) magnitudes indicates the core helix for R^8,16 GWALP23 is significantly tilted to give a similar transmembrane orientation in thinner bilayers with either saturated C12:0 or C14:0 acyl chains (1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)) or unsaturated C16:1 Δ9 cis acyl chains. In bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC; C18:1 Δ9 cis) multiple orientations are indicated, whereas in longer, unsaturated 1,2-dieicosenoyl-sn-glycero-3-phosphocholine (DEiPC; C20:1 Δ11 cis) bilayers, the R^8,16 GWALP23 helix adopts primarily a surface orientation. The inclusion of 10-20 mol % cholesterol in DOPC bilayers drives more of the R^8,16 GWALP23 helix population to the membrane surface, thereby allowing both charged arginines access to the interfacial lipid head groups. The results suggest that hydrophobic thickness and cholesterol content are more important than lipid saturation for the arginine peptide dynamics and helix orientation in lipid membranes.

Clinical features and outcome of 6 new patients carrying de novo KCNB1 gene mutations

Objective

To describe electroclinical features and outcome of 6 patients harboring KCNB1 mutations.

Methods

Clinical, EEG, neuropsychological, and brain MRI data analysis. Targeted next-generation sequencing of a 95 epilepsy gene panel.

Results

The mean age at seizure onset was 11 months. The mean follow-up of 11.3 years documented that 4 patients following an infantile phase of frequent seizures became seizure free; the mean age at seizure offset was 4.25 years. Epilepsy phenotypes comprised West syndrome in 2 patients, infantile-onset unspecified generalized epilepsy, myoclonic and photosensitive eyelid myoclonia epilepsy resembling Jeavons syndrome, Lennox-Gastaut syndrome, and focal epilepsy with prolonged occipital or clonic seizures in each and every one. Five patients had developmental delay prior to seizure onset evolving into severe intellectual disability with absent speech and autistic traits in one and stereotypic hand movements with impulse control disorder in another. The patient with Jeavons syndrome evolved into moderate intellectual disability. Mutations were de novo, 4 missense and 2 nonsense, 5 were novel, and 1 resulted from somatic mosaicism.

Conclusions

KCNB1-related manifestations include a spectrum of infantile-onset generalized or focal seizures whose combination leads to early infantile epileptic encephalopathy including West, Lennox-Gastaut, and Jeavons syndromes. Long-term follow-up highlights that following a stormy phase, seizures subside or cease and treatment may be eased or withdrawn. Cognitive and motor functions are almost always delayed prior to seizure onset and evolve into severe, persistent impairment. Thus, KCNB1 mutations are associated with diffuse brain dysfunction combining seizures, motor, and cognitive impairment.

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.

Atomistic Modeling of Ion Conduction through the Voltage-Sensing Domain of the Shaker K+ Ion Channel

Voltage-sensing domains (VSDs) sense changes in the membrane electrostatic potential and, through conformational changes, regulate a specific function. The VSDs of wild-type voltage-dependent K⁺, Na⁺, and Ca²⁺ channels do not conduct ions, but they can become ion-permeable through pathological mutations in the VSD. Relatively little is known about the underlying mechanisms of conduction through VSDs. The most detailed studies have been performed on Shaker K⁺ channel variants in which ion conduction through the VSD is manifested in electrophysiology experiments as a voltage-dependent inward current, the so-called omega current, which appears when the VSDs are in their resting state conformation. Only monovalent cations appear to permeate the Shaker VSD via a pathway that is believed to be, at least in part, the same as that followed by the S4 basic side chains during voltage-dependent activation. We performed μs-time scale atomistic molecular dynamics simulations of a cation-conducting variant of the Shaker VSD under applied electric fields in an experimentally validated resting-state conformation, embedded in a lipid bilayer surrounded by solutions containing guanidinium chloride or potassium chloride. Our simulations provide insights into the Shaker VSD permeation pathway, the protein-ion interactions that control permeation kinetics, and the mechanism of voltage-dependent activation of voltage-gated ion channels.

Regulation of KCNQ/Kv7 family voltage-gated K+ channels by lipids

Many years of studies have established that lipids can impact membrane protein structure and function through bulk membrane effects, by direct but transient annular interactions with the bilayer-exposed surface of protein transmembrane domains, and by specific binding to protein sites. Here, we focus on how phosphatidylinositol 4,5-bisphosphate (PIP₂) and polyunsaturated fatty acids (PUFAs) impact ion channel function and how the structural details of the interactions of these lipids with ion channels are beginning to emerge. We focus on the Kv7 (KCNQ) subfamily of voltage-gated K⁺ channels, which are regulated by both PIP₂ and PUFAs and play a variety of important roles in human health and disease. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.