The pore of voltage-gated potassium ion channels is strained when closed

Personalized structural biology reveals the molecular mechanisms underlying heterogeneous epileptic phenotypes caused by de novo KCNC2 variants

Whole-exome sequencing (WES) in the clinic has identified several rare monogenic developmental and epileptic encephalopathies (DEE) caused by ion channel variants. However, WES often fails to provide actionable insight for rare diseases, such as DEEs, due to the challenges of interpreting variants of unknown significance (VUS). Here, we describe a "personalized structural biology" (PSB) approach that leverages recent innovations in the analysis of protein 3D structures to address this challenge. We illustrate this approach in an Undiagnosed Diseases Network (UDN) individual with DEE symptoms and a de novo VUS in KCNC2 (p.V469L), the Kv3.2 voltage-gated potassium channel. A nearby KCNC2 variant (p.V471L) was recently suggested to cause DEE-like phenotypes. Computational structural modeling suggests that both affect protein function. However, despite their proximity, the p.V469L variant is likely to sterically block the channel pore, while the p.V471L variant is likely to stabilize the open state. Biochemical and electrophysiological analyses demonstrate heterogeneous loss-of-function and gain-of-function effects, as well as differential response to 4-aminopyridine treatment. Molecular dynamics simulations illustrate that the pore of the p.V469L variant is more constricted, increasing the energetic barrier for K+ permeation, whereas the p.V471L variant stabilizes the open conformation. Our results implicate variants in KCNC2 as causative for DEE and guide the interpretation of a UDN individual. They further delineate the molecular basis for the heterogeneous clinical phenotypes resulting from two proximal pathogenic variants. This demonstrates how the PSB approach can provide an analytical framework for individualized hypothesis-driven interpretation of protein-coding VUS.

Weak Cation Selectivity in HCN Channels Results From K+-Mediated Release of Na+ From Selectivity Filter Binding Sites

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels generate the pacemaker current which plays an important role in the timing of various biological processes like the heart beat. We used umbrella sampling to explore the potential of mean force for the conduction of potassium and sodium through the open HCN4 pore. Our data explain distinct functional features like low unitary conductance and weak selectivity as a result of high energetic barriers inside the selectivity filter of this channel. They exceed the 3-5 kJ/mol threshold which is presumed as maximal barrier for diffusion-limited conductance. Furthermore, simulations provide a thermodynamic explanation for the weak cation selectivity of HCN channels that contain only two ion binding sites in the selectivity filter (SF). We find that sodium ions bind more strongly to the SF than potassium and are easier released by binding of potassium than of another sodium. Hence ion transport and selectivity in HCN channels is not determined by the same mechanism as in potassium-selective channels; it rather relies on sodium as a weak blocker that can only be released by potassium.

Characterizing fenestration size in sodium channel subtypes and their accessibility to inhibitors

Voltage-gated sodium channels (Nav) underlie the electrical activity of nerve and muscle cells. Humans have nine different subtypes of these channels, which are the target of small-molecule inhibitors commonly used to treat a range of conditions. Structural studies have identified four lateral fenestrations within the Nav pore module that have been shown to influence Nav pore blocker access during resting-state inhibition. However, the structural differences among the nine subtypes are still unclear. In particular, the dimensions of the four individual fenestrations across the Nav subtypes and their differential accessibility to pore blockers is yet to be characterized. To address this, we applied classical molecular dynamics simulations to study the recently published structures of Nav1.1, Nav1.2, Nav1.4, Nav1.5, and Nav1.7. Although there is significant variability in the bottleneck sizes of the Nav fenestrations, the subtypes follow a common pattern, with wider DI-II and DIII-IV fenestrations, a more restricted DII-III fenestration, and the most restricted DI-IV fenestration. We further identify the key bottleneck residues in each fenestration and show that the motions of aromatic residue sidechains govern the bottleneck radii. Well-tempered metadynamics simulations of Nav1.4 and Nav1.5 in the presence of the pore blocker lidocaine also support the DI-II fenestration being the most likely access route for drugs. Our computational results provide a foundation for future in vitro experiments examining the route of drug access to sodium channels. Understanding the fenestrations and their accessibility to drugs is critical for future analyses of diseases mutations across different sodium channel subtypes, with the potential to inform pharmacological development of resting-state inhibitors and subtype-selective drug design.

Ion Selectivity and Permeation Mechanism in a Cyclodextrin-Based Channel

Synthetic ion channels are a promising technology in the medical and materials sciences because of their ability to conduct ions. Channels based on cyclodextrin, a cyclic oligomer of glucose, are of particular interest because of their nontoxicity and biocompatibility. Using molecular dynamics-based free energy calculations, this study identifies cyclodextrin channel types that are best suited to serve as synthetic ion channels. Free energy profiles show that the connectivity in the channel determines whether the channel is cation-selective or anion-selective. Furthermore, the energy barrier for ion transport is governed by the number of glucose molecules making up the cyclodextrin units of the channel. A detailed mechanism is proposed for ion transport through these channels. Findings from this study will aid in designing cyclodextrin-based channels that could be either cation-selective or anion-selective, by modifying the linkages of the channel or the number of glucose molecules in the cyclodextrin rings.

Structural and molecular insight into the pH-induced low-permeability of the voltage-gated potassium channel Kv1.2 through dewetting of the water cavity

Understanding the gating mechanism of ion channel proteins is key to understanding the regulation of cell signaling through these channels. Channel opening and closing are regulated by diverse environmental factors that include temperature, electrical voltage across the channel, and proton concentration. Low permeability in voltage-gated potassium ion channels (Kv) is intimately correlated with the prolonged action potential duration observed in many acidosis diseases. The Kv channels consist of voltage-sensing domains (S1-S4 helices) and central pore domains (S5-S6 helices) that include a selectivity filter and water-filled cavity. The voltage-sensing domain is responsible for the voltage-gating of Kv channels. While the low permeability of Kv channels to potassium ion is highly correlated with the cellular proton concentration, it is unclear how an intracellular acidic condition drives their closure, which may indicate an additional pH-dependent gating mechanism of the Kv family. Here, we show that two residues E327 and H418 in the proximity of the water cavity of Kv1.2 play crucial roles as a pH switch. In addition, we present a structural and molecular concept of the pH-dependent gating of Kv1.2 in atomic detail, showing that the protonation of E327 and H418 disrupts the electrostatic balance around the S6 helices, which leads to a straightening transition in the shape of their axes and causes dewetting of the water-filled cavity and closure of the channel. Our work offers a conceptual advancement to the regulation of the pH-dependent gating of various voltage-gated ion channels and their related biological functions.

Atomistic modeling and molecular dynamics analysis of human CaV1.2 channel using external electric field and ion pulling simulations

BACKGROUND

Human Ca_V1.2 (hCav1.2), a calcium selective voltage-gated channel, plays important roles in normal cardiac and neuronal functions. Calcium influx and gating mechanisms leading to the activation of hCa_V1.2 are critical for its functionalities. Lack of an experimentally resolved structure of hCa_V1.2 remains a significant impediment in molecular-level understanding of this channel. This work focuses on building atomistic hCa_V1.2 model and studying calcium influx using computational approaches.

METHODS

We employed homology modeling and molecular dynamics (MD) to build the structure of hCa_V1.2. Subsequently, we employed steered molecular dynamics (SMD) to understand calcium ion permeation in hCa_V1.2.

RESULTS

We report a comprehensive three-dimensional model of a closed state hCa_V1.2 refined under physiological membrane-bound conditions using MD simulations. Our SMD simulations on the model revealed four important barriers for ion permeation: this includes three calcium binding sites formed by the EEEE- and TTTT- rings within the selectivity filter region and a large barrier rendered by the hydrophobic internal gate. Our results also revealed that the first hydration shell of calcium remained intact throughout the simulations, thus playing an important role in ion permeation in hCa_V1.2.

CONCLUSIONS

Our results have provided some important mechanistic insights into the structure, dynamics and ion permeation in hCa_V1.2. The significant barriers for ion permeation formed by the four phenylalanine residues at the internal gate region suggest that this site is important for channel activation.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

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.

Computational Study of the Loss-of-Function Mutations in the Kv1.5 Channel Associated with Atrial Fibrillation

Atrial fibrillation (AF) is a heart disease caused by defective ion channels in the atria, which affect the action potential (AP) duration and disturb normal heart rhythm. Rapid firing of APs in neighboring atrial cells is a common mechanism of AF, and therefore, therapeutic approaches have focused on extending the AP duration by inhibiting the K⁺ channels involved in repolarization. Of these, Kv1.5 that carries the I _Kur current is a promising target because it is expressed mainly in atria and not in ventricles. In genetic studies of AF patients, both loss-of-function and gain-of-function mutations in Kv1.5 have been identified, indicating that either decreased or increased I _Kur currents could trigger AF. Blocking of already downregulated Kv1.5 channels could cause AF to become chronic. Thus, a molecular-level understanding of how the loss-of-function mutations in Kv1.5 affect I _Kur would be useful for developing new therapeutics. Here, we perform molecular dynamics simulations to study the effect of three loss-of-function mutations in the pore domain of Kv1.5 on ion permeation. Comparison of the pore structures and ion free energies in the wild-type and mutant Kv1.5 channels indicates that conformational changes in the selectivity filter could hinder ion permeation in the mutant channels.

The role of entropic potential in voltage activation and K+ transport through Kv 1.2 channels

We analyze the entropic effects of inner pore geometry changes of Kv 1.2 channel during membrane depolarization and their implications for the rate of transmembrane transport of potassium ions. We base this on the idea that spatial confinements within the channel pore give rise to entropic barriers which can both effectively affect the stability of open macroconformation and influence channel's ability to conduct the potassium ions through the membrane. First, we calculate the differences in entropy between voltage-activated and resting states of the channel. As a template, we take a set of structures of channel pore in an open state at different membrane potentials generated in our previous research. The obtained results indicate that tendency to occupy open states at membrane depolarization is entropy facilitated. Second, we describe the differences in rates of K⁺ transport through the channel pore at different voltages based on the results of appropriate random walk simulations in entropic and electric potentials. The simulated single channel currents (I) suggest that the geometry changes during membrane depolarization are an important factor contributing to the observed flow of potassium ions through the channel. Nevertheless, the charge distribution within the channel pore (especially at the extracellular entrance) seems most prominent for the observed I/I_max relation at a qualitative level at analyzed voltages.

Effects of protein-protein interactions and ligand binding on the ion permeation in KCNQ1 potassium channel

The voltage-gated KCNQ1 potassium ion channel interacts with the type I transmembrane protein minK (KCNE1) to generate the slow delayed rectifier (I_Ks) current in the heart. Mutations in these transmembrane proteins have been linked with several heart-related issues, including long QT syndromes (LQTS), congenital atrial fibrillation, and short QT syndrome. Off-target interactions of several drugs with that of KCNQ1/KCNE1 ion channel complex have been known to cause fatal cardiac irregularities. Thus, KCNQ1/KCNE1 remains an important avenue for drug-design and discovery research. In this work, we present the structural and mechanistic details of potassium ion permeation through an open KCNQ1 structural model using the combined molecular dynamics and steered molecular dynamics simulations. We discuss the processes and key residues involved in the permeation of a potassium ion through the KCNQ1 ion channel, and how the ion permeation is affected by (i) the KCNQ1-KCNE1 interactions and (ii) the binding of chromanol 293B ligand and its derivatives into the complex. The results reveal that interactions between KCNQ1 with KCNE1 causes a pore constriction in the former, which in-turn forms small energetic barriers in the ion-permeation pathway. These findings correlate with the previous experimental reports that interactions of KCNE1 dramatically slows the activation of KCNQ1. Upon ligand-binding onto the complex, the energy-barriers along ion permeation path are more pronounced, as expected, therefore, requiring higher force in our steered-MD simulations. Nevertheless, pulling the ion when a weak blocker is bound to the channel does not necessitate high force in SMD. This indicates that our SMD simulations have been able to discern between strong and week blockers and reveal their influence on potassium ion permeation. The findings presented here will have some implications in understanding the potential off-target interactions of the drugs with the KCNQ1/KCNE1 channel that lead to cardiotoxic effects.

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.

Effect of the Southeast Asian Ovalocytosis Deletion on the Conformational Dynamics of Signal-Anchor Transmembrane Segment 1 of Red Cell Anion Exchanger 1 (AE1, Band 3, or SLC4A1)

The first transmembrane (TM1) helix in the red cell anion exchanger (AE1, Band 3, or SLC4A1) acts as an internal signal anchor that binds the signal recognition particle and directs the nascent polypeptide chain to the endoplasmic reticulum (ER) membrane where it moves from the translocon laterally into the lipid bilayer. The sequence N-terminal to TM1 forms an amphipathic helix that lies at the membrane interface and is connected to TM1 by a bend at Pro403. Southeast Asian ovalocytosis (SAO) is a red cell abnormality caused by a nine-amino acid deletion (Ala400-Ala408) at the N-terminus of TM1. Here we demonstrate, by extensive (∼4.5 μs) molecular dynamics simulations of TM1 in a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine membrane, that the isolated TM1 peptide is highly dynamic and samples the structure of TM1 seen in the crystal structure of the membrane domain of AE1. The SAO deletion not only removes the proline-induced bend but also causes a "pulling in" of the part of the amphipathic helix into the hydrophobic phase of the bilayer, as well as the C-terminal of the peptide. The dynamics of the SAO peptide very infrequently resembles the structure of TM1 in AE1, demonstrating the disruptive effect the SAO deletion has on AE1 folding. These results provide a precise molecular view of the disposition and dynamics of wild-type and SAO TM1 in a lipid bilayer, an important early biosynthetic intermediate in the insertion of AE1 into the ER membrane, and extend earlier results of cell-free translation experiments.

Contributions of the membrane dipole potential to the function of voltage-gated cation channels and modulation by small molecule potentiators

The membrane dipole potential (Ψ_d) constitutes one of three electrical potentials generated by cell membranes. Ψ_d arises from the unfavorable parallel alignment of phospholipid and water dipoles, and varies in magnitude both longitudinally and laterally across the bilayer according to membrane composition and phospholipid packing density. In this work, we propose that dynamic counter-balancing between Ψ_d and the transmembrane potential (ΔΨ_m) governs the conformational state transitions of voltage-gated ion channels. Ψ_d consists of 1) static outer, and dynamic inner leaflet components (Ψ_d(extra) and Ψ_d(intra), respectively); and 2) a transmembrane component (ΔΨ_{d(inner-outer)}), ariing from differences in intra- and extracellular leaflet composition. Ψ_d(intra), which transitions between high and low energy states (Ψ_{d(intra, high)} and Ψ_{d(intra, low)}) as a function of channel conformation, is transduced by the pore domain. ΔΨ_{d(inner-outer)} is transduced by the voltage-sensing (VS) domain in summation with ΔΨ_m. Potentiation of voltage-gated ion channels is of interest for the treatment of cardiac, neuronal, and other disorders arising from inherited/acquired ion channel dysfunction. Potentiators are widely believed to alter the rates and voltage-dependencies of channel gating transitions by binding to pockets in the membrane-facing and other regions of ion channel targets. Here, we propose that potentiators alter Ψ_d(intra) and/or Ψ_d(extra), thereby increasing or decreasing the energy barriers governing channel gating transitions. We used quantum mechanical and molecular dynamics (MD) simulations to predict the overall Ψ_d-modulating effects of a series of published positive hERG potentiators partitioned into model DOPC bilayers. Our findings suggest a strong correlation between the magnitude of Ψ_d-lowering and positive hERG potentiation across the series.

TRPV1: A Target for Rational Drug Design

Transient Receptor Potential Vanilloid 1 (TRPV1) is a non-selective, Ca²⁺ permeable cation channel activated by noxious heat, and chemical ligands, such as capsaicin and resiniferatoxin (RTX). Many compounds have been developed that either activate or inhibit TRPV1, but none of them are in routine clinical practice. This review will discuss the rationale for antagonists and agonists of TRPV1 for pain relief and other conditions, and strategies to develop new, better drugs to target this ion channel, using the newly available high-resolution structures.

Helix perturbations in membrane proteins assist in inter-helical interactions and optimal helix positioning in the bilayer

Transmembrane (TM) helices in integral membrane proteins are primarily α-helical in structure. Here we analyze 1134 TM helices in 90 high resolution membrane proteins and find that apart from the widely prevalent α-helices, TM regions also contain stretches of 3₁₀ (3 to 8 residues) and π-helices (5 to 19 residues) with distinct sequence signatures. The various helix perturbations in TM regions comprise of helices with kinked geometry, as well as those with an interspersed 3₁₀/π-helical fragment and show high occurrence in a few membrane proteins. Proline is frequently present at sites of these perturbations, but it is neither a necessary nor a sufficient requirement. Helix perturbations are also conserved within a family of membrane proteins despite low sequence identity in the perturbed region. Furthermore, a perturbation influences the geometry of the TM helix, mediates inter-helical interactions within and across protein chains and avoids hydrophobic mismatch of the helix termini with the bilayer. An analysis of π-helices in the TM regions of the heme copper oxidase superfamily shows that interspersed π-helices can vary in length from 6 to 19 amino acids or be entirely absent, depending upon the protein function. The results presented here would be helpful for prediction of 3₁₀ and π-helices in TM regions and can assist the computational design of membrane proteins.

On Application Of Langevin Dynamics In Logarithmic Potential To Model Ion Channel Gate Activity

We model the activity of an ion channel gate by Langevin dynamics in a logarithmic potential. This approach enables one to describe the power-law dwell-time distributions of the considered system, and the long-term correlations between the durations of the subsequent channel states, or fractal scaling of statistical characteristics of the gate's movement with time. Activity of an ion channel gate is described as an overdamped motion of the reaction coordinate in a confining logarithmic potential, which ensures great flexibility of the model. Depending on the chosen parameters, it allows one to reproduce many types of gate dynamics within the family of non-Markovian, anomalous conformational diffusion processes. In this study we apply the constructed model to largeconductance voltage and Ca2+-activated potassium channels (BKCa). The interpretation of model assumptions and parameters is provided in terms of this biological system. Our results show good agreement with the experimental data.

Grafting voltage and pharmacological sensitivity in potassium channels

A classical voltage-gated ion channel consists of four voltage-sensing domains (VSDs). However, the roles of each VSD in the channels remain elusive. We developed a GVTDT (Graft VSD To Dimeric TASK3 channels that lack endogenous VSDs) strategy to produce voltage-gated channels with a reduced number of VSDs. TASK3 channels exhibit a high host tolerance to VSDs of various voltage-gated ion channels without interfering with the intrinsic properties of the TASK3 selectivity filter. The constructed channels, exemplified by the channels grafted with one or two VSDs from Kv7.1 channels, exhibit classical voltage sensitivity, including voltage-dependent opening and closing. Furthermore, the grafted Kv7.1 VSD transfers the potentiation activity of benzbromarone, an activator that acts on the VSDs of the donor channels, to the constructed channels. Our study indicates that one VSD is sufficient to voltage-dependently gate the pore and provides new insight into the roles of VSDs.

Kainate receptor pore-forming and auxiliary subunits regulate channel block by a novel mechanism

KEY POINTS

Kainate receptor heteromerization and auxiliary subunits, Neto1 and Neto2, attenuate polyamine ion-channel block by facilitating blocker permeation. Relief of polyamine block in GluK2/GluK5 heteromers results from a key proline residue that produces architectural changes in the channel pore α-helical region. Auxiliary subunits exert an additive effect to heteromerization, and thus relief of polyamine block is due to a different mechanism. Our findings have broad implications for work on polyamine block of other cation-selective ion channels.

ABSTRACT

Channel block and permeation by cytoplasmic polyamines is a common feature of many cation-selective ion channels. Although the channel block mechanism has been studied extensively, polyamine permeation has been considered less significant as it occurs at extreme positive membrane potentials. Here, we show that kainate receptor (KAR) heteromerization and association with auxiliary proteins, Neto1 and Neto2, attenuate polyamine block by enhancing blocker permeation. Consequently, polyamine permeation and unblock occur at more negative and physiologically relevant membrane potentials. In GluK2/GluK5 heteromers, enhanced permeation is due to a single proline residue in GluK5 that alters the dynamics of the α-helical region of the selectivity filter. The effect of auxiliary proteins is additive, and therefore the structural basis of polyamine permeation and unblock is through a different mechanism. As native receptors are thought to assemble as heteromers in complex with auxiliary proteins, our data identify an unappreciated impact of polyamine permeation in shaping the signalling properties of neuronal KARs and point to a structural mechanism that may be shared amongst other cation-selective ion channels.

Modulation of cell function by electric field: a high-resolution analysis

Regulation of cell function by a non-thermal, physiological-level electromagnetic field has potential for vascular tissue healing therapies and advancing hybrid bioelectronic technology. We have recently demonstrated that a physiological electric field (EF) applied wirelessly can regulate intracellular signalling and cell function in a frequency-dependent manner. However, the mechanism for such regulation is not well understood. Here, we present a systematic numerical study of a cell-field interaction following cell exposure to the external EF. We use a realistic experimental environment that also recapitulates the absence of a direct electric contact between the field-sourcing electrodes and the cells or the culture medium. We identify characteristic regimes and present their classification with respect to frequency, location, and the electrical properties of the model components. The results show a striking difference in the frequency dependence of EF penetration and cell response between cells suspended in an electrolyte and cells attached to a substrate. The EF structure in the cell is strongly inhomogeneous and is sensitive to the physical properties of the cell and its environment. These findings provide insight into the mechanisms for frequency-dependent cell responses to EF that regulate cell function, which may have important implications for EF-based therapies and biotechnology development.

Calcium ions open a selectivity filter gate during activation of the MthK potassium channel

Ion channel opening and closing are fundamental to cellular signalling and homeostasis. Gates that control K⁺ channel activity were found both at an intracellular pore constriction and within the selectivity filter near the extracellular side but the specific location of the gate that opens Ca²⁺-activated K⁺ channels has remained elusive. Using the Methanobacterium thermoautotrophicum homologue (MthK) and a stopped-flow fluorometric assay for fast channel activation, we show that intracellular quaternary ammonium blockers bind to closed MthK channels. Since the blockers are known to bind inside a central channel cavity, past the intracellular entryway, the gate must be within the selectivity filter. Furthermore, the blockers access the closed channel slower than the open channel, suggesting that the intracellular entryway narrows upon pore closure, without preventing access of either the blockers or the smaller K⁺. Thus, Ca²⁺-dependent gating in MthK occurs at the selectivity filter with coupled movement of the intracellular helices.

Ion channels open and close to allow the regulated passage of ions through the membrane. Here the authors use selective ion channel blockers to analyse this regulation in a potassium channel and show that the gate is in the selectivity filter, past the entrance to the channel.

Mechanism-Based Mathematical Model for Gating of Ionotropic Glutamate Receptors

We present a mathematical model for ionotropic glutamate receptors (iGluR's) that is built on mechanistic understanding and yields a number of thermodynamic and kinetic properties of channel gating. iGluR's are ligand-gated ion channels responsible for the vast majority of fast excitatory neurotransmission in the central nervous system. The effects of agonist-induced closure of the ligand-binding domain (LBD) are transmitted to the transmembrane channel (TMC) via interdomain linkers. Our model demonstrates that, relative to full agonists, partial agonists may reduce either the degree of LBD closure or the curvature of the LBD free energy basin, leading to less stabilization of the channel open state and hence lower channel open probability. A rigorous relation is derived between the channel closed-to-open free energy difference and the tension within the linker. Finally, by treating LBD closure and TMC opening as diffusive motions, we obtain gating trajectories that resemble stochastic current traces from single-channel recordings and calculate the rate constants for transitions between the channel open and closed states. Our model can be implemented by molecular dynamics simulations to realistically depict iGluR gating and may guide functional experiments in gaining deeper insight into this essential family of channel proteins.

Voltage-dependent gating of KCNH potassium channels lacking a covalent link between voltage-sensing and pore domains

Voltage-gated channels open paths for ion permeation upon changes in membrane potential, but how voltage changes are coupled to gating is not entirely understood. Two modules can be recognized in voltage-gated potassium channels, one responsible for voltage sensing (transmembrane segments S1 to S4), the other for permeation (S5 and S6). It is generally assumed that the conversion of a conformational change in the voltage sensor into channel gating occurs through the intracellular S4–S5 linker that provides physical continuity between the two regions. Using the pathophysiologically relevant KCNH family, we show that truncated proteins interrupted at, or lacking the S4–S5 linker produce voltage-gated channels in a heterologous model that recapitulate both the voltage-sensing and permeation properties of the complete protein. These observations indicate that voltage sensing by the S4 segment is transduced to the channel gate in the absence of physical continuity between the modules.

The pore of voltage-gated ion channels opens in response to membrane depolarization sensed by a separate voltage-sensing domain. Here, Lörinczi et al. show that, contrary to assumptions, no physical linker is required to transmit changes from the voltage-sensing to the permeation domain of KCNH channels.

Structural dynamics and conformational equilibria of SERCA regulatory proteins in membranes by solid-state NMR restrained simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.

Gating topology of the proton-coupled oligopeptide symporters

Proton-coupled oligopeptide transporters belong to the major facilitator superfamily (MFS) of membrane transporters. Recent crystal structures suggest the MFS fold facilitates transport through rearrangement of their two six-helix bundles around a central ligand binding site; how this is achieved, however, is poorly understood. Using modeling, molecular dynamics, crystallography, functional assays, and site-directed spin labeling combined with double electron-electron resonance (DEER) spectroscopy, we present a detailed study of the transport dynamics of two bacterial oligopeptide transporters, PepTSo and PepTSt. Our results identify several salt bridges that stabilize outward-facing conformations and we show that, for all the current structures of MFS transporters, the first two helices of each of the four inverted-topology repeat units form half of either the periplasmic or cytoplasmic gate and that these function cooperatively in a scissor-like motion to control access to the peptide binding site during transport.

Reduced curvature of ligand-binding domain free-energy surface underlies partial agonism at NMDA receptors

NMDA receptors are ligand-gated ion channels that mediate excitatory synaptic transmission in the central nervous system. Partial agonists elicit submaximal channel activation, but crystal structures of the ligand-binding domains (LBDs) bound with partial and full agonists show little difference. To uncover the molecular mechanism for partial agonism, here we computed the free-energy surfaces of the GluN1 (an obligatory subunit of NMDA receptors) LBD bound with a variety of ligands. The free-energy minima are similarly positioned for full and partial agonists, but the curvatures are significantly reduced in the latter case, indicating higher probabilities for sampling conformations with a not fully closed domain cleft. The free-energy surfaces for antagonists have both shifted minima and further reduced curvatures. Reduced curvature of free-energy surface appears to explain well the partial agonism at NMDA receptors and may present a unique paradigm in producing graded responses for receptors in general.

Structure based computational assessment of channel properties of assembled ORF-8a from SARS-CoV

ORF 8a is a short 39 amino acid bitopic membrane protein encoded by severe acute respiratory syndrome causing corona virus (SARS‐CoV). It has been identified to increase permeability of the lipid membrane for cations. Permeability is suggested to occur due to the assembly of helical bundles. Computational models of a pentameric assembly of 8a peptides are generated using the first 22 amino acids, which include the transmembrane domain. Low energy structures reveal a hydrophilic pore mantled by residues Thr‐8, and −18, Ser‐11, Cys‐13, and Arg‐22. Potential of mean force (PMF) profiles for mono (Na⁺, K⁺, Cl⁻) and divalent (Ca²⁺) ions along the pore are calculated. The data support experimental findings of a weak cation selectivity of the channel. Calculations on 8a are compared to data derived for a pentameric bundle consisting of the M2 helices of the bacterial pentameric ligand gated ion channel GLIC (3EHZ). PMF curves of both, bundles 8a and M2, show sigmoidal shaped profiles. In comparison to the data for the M2‐GLIC model, data of the 8a bundle show lower amplitude of the PMF values between maximum and minimum and less discrimination amongst ions. Proteins 2015; 83:300–308. © 2014 Wiley Periodicals, Inc.

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

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

KCNC3(R420H), a K(+) channel mutation causative in spinocerebellar ataxia 13 displays aberrant intracellular trafficking

Spinocerebellar ataxia 13 (SCA13) is an autosomal dominant disease resulting from mutations in KCNC3 (Kv3.3), a voltage-gated potassium channel. The KCNC3(R420H) mutation was first identified as causative for SCA13 in a four-generation Filipino kindred with over 20 affected individuals. Electrophysiological analyses in oocytes previously showed that this mutation did not lead to a functional channel and displayed a dominant negative phenotype. In an effort to identify the molecular basis of this allelic form of SCA13, we first determined that human KCNC3(WT) and KCNC3(R420H) display disparate post-translational modifications, and the mutant protein has reduced complex glycan adducts. Immunohistochemical analyses demonstrated that KCNC3(R420H) was not properly trafficking to the plasma membrane and surface biotinylation demonstrated that KCNC3(R420H) exhibited only 24% as much surface expression as KCNC3(WT). KCNC3(R420H) trafficked through the ER but was retained in the Golgi. KCNC3(R420H) expression results in altered Golgi and cellular morphology. Electron microscopy of KCNC3(R420H) localization further supports retention in the Golgi. These results are specific to the KCNC3(R420H) allele and provide new insight into the molecular basis of disease manifestation in SCA13.

General rules for the arrangements and gating motions of pore-lining helices in homomeric ion channels

The pore-lining helix (PLH) bundles are central to the function of all ion channels, as their conformational rearrangements dictate channel gating. Here we explore all plausible oligomeric arrangements of the PLH bundles within homomeric ion channels by building models using generic restraints. In particular, the distance between two neighbouring PLHs was bounded both below and above in order to avoid steric clash and allow proper packing. The resulting models provide a theoretical representation of the accessible space for oligomeric arrangements. While the represented space is confined, it encompasses nearly all the ion channel PLH bundles for which the structures are currently known. For a multitude of channels, gating models suggested by paths within the confined accessible space are in qualitative agreement with those established in previous structural and computational studies.

Rearrangements of the pore-lining helix (PLH) bundles of ion channels are central to their gating mechanisms. Here, Dai et al. use a modelling approach to define the general rules that govern the arrangements and gating motions of the PLHs in homomeric ion channels.

P-loop conformation governed crizotinib resistance in G2032R-mutated ROS1 tyrosine kinase: clues from free energy landscape

Tyrosine kinases are regarded as excellent targets for chemical drug therapy of carcinomas. However, under strong purifying selection, drug resistance usually occurs in the cancer cells within a short term. Many cases of drug resistance have been found to be associated with secondary mutations in drug target, which lead to the attenuated drug-target interactions. For example, recently, an acquired secondary mutation, G2032R, has been detected in the drug target, ROS1 tyrosine kinase, from a crizotinib-resistant patient, who responded poorly to crizotinib within a very short therapeutic term. It was supposed that the mutation was located at the solvent front and might hinder the drug binding. However, a different fact could be uncovered by the simulations reported in this study. Here, free energy surfaces were characterized by the drug-target distance and the phosphate-binding loop (P-loop) conformational change of the crizotinib-ROS1 complex through advanced molecular dynamics techniques, and it was revealed that the more rigid P-loop region in the G2032R-mutated ROS1 was primarily responsible for the crizotinib resistance, which on one hand, impaired the binding of crizotinib directly, and on the other hand, shortened the residence time induced by the flattened free energy surface. Therefore, both of the binding affinity and the drug residence time should be emphasized in rational drug design to overcome the kinase resistance.

Cancers can eventually confer drug resistance to the continued medication. In most cases, mutations occurred in a drug target can attenuate the binding affinity of the drugs. Here, we studied the drug resistance mechanisms of the mutations G2032R in the ROS1 tyrosine kinase in fusion-type NSCLC. It is well known that the phosphate-binding loop (P-loop) plays a vital role in the binding of competitive inhibitors in tyrosine kinases, and numerous mutations have been found occurred around the P-loop, which may affect the binding/unbinding process of a drug. Free energy surfaces were constructed to characterize the impact of the mutation to the binding/unbinding process of a well-known NSCLC drug, crizotinib. Two advanced free energy calculation methods, namely funnel based well-tempered metadynamics and umbrella sampling based absolute binding free energy calculation achieved consistent results with the experimental data, suggesting that the rigid P-loop of the mutated target was mainly responsible for the crizotinib resistance to ROS1 tyrosine kinase.

Flexible gates generate occluded intermediates in the transport cycle of LacY

The major facilitator superfamily (MFS) transporter lactose permease (LacY) alternates between cytoplasmic and periplasmic open conformations to co-transport a sugar molecule together with a proton across the plasma membrane. Indirect experimental evidence suggested the existence of an occluded transition intermediate of LacY, which would prevent leaking of the proton gradient. As no experimental structure is known, the conformational transition is not fully understood in atomic detail. We simulated transition events from a cytoplasmic open conformation to a periplasmic open conformation with the dynamic importance sampling molecular dynamics method and observed occluded intermediates. Analysis of water permeation pathways and the electrostatic free-energy landscape of a solvated proton indicated that the occluded state contains a solvated central cavity inaccessible from either side of the membrane. We propose a pair of geometric order parameters that capture the state of the pathway through the MFS transporters as shown by a survey of available crystal structures and models. We present a model for the occluded state of apo-LacY, which is similar to the occluded crystal structures of the MFS transporters EmrD, PepT_So, NarU, PiPT and XylE. Our simulations are consistent with experimental double electron spin–spin distance measurements that have been interpreted to show occluded conformations. During the simulations, a salt bridge that has been postulated to be involved in driving the conformational transition formed. Our results argue against a simple rigid-body domain motion as implied by a strict “rocker-switch mechanism” and instead hint at an intricate coupling between two flexible gates.

•

The transport mechanism of LacY is hypothesized to involve an intermediate “occluded” state.

•

Such a state is observed in computer simulations of the conformational transitions.

•

Simulation data are validated with experimental double electron–electron spin resonance measurements.

•

The structural gating elements of LacY are identified.

•

Occluded LacY is similar to known occluded structures of homologous proteins.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Detailed Examination of a Single Conduction Event in a Potassium Channel

Although extensively studied, it has proved difficult to describe in detail how potassium ion channels conduct cations and water. We present a computational study that, by using stratified umbrella sampling, examines nearly an entire conduction event of the Kv1.2/2.1 paddle chimera and thereby identifies the expected stable configurations of ions and waters in the selectivity filter of the channel. We describe in detail the motions of the ions and waters during a conduction event, focusing on how waters and ions enter the filter, the rotation of water molecules inside the filter, and how potassium ions are coordinated as they move from a water to a protein environment. Finally, we analyze the small conformational changes undergone by the protein, showing that the stable configurations are most similar to the experimental crystal structure.