Conformational Changes and Slow Dynamics through Microsecond Polarized Atomistic Molecular Simulation of an Integral Kv1.2 Ion Channel

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations

A Threshold Equation for Action Potential Initiation

In central neurons, the threshold for spike initiation can depend on the stimulus and varies between cells and between recording sites in a given cell, but it is unclear what mechanisms underlie this variability. Properties of ionic channels are likely to play a role in threshold modulation. We examined in models the influence of Na channel activation, inactivation, slow voltage-gated channels and synaptic conductances on spike threshold. We propose a threshold equation which quantifies the contribution of all these mechanisms. It provides an instantaneous time-varying value of the threshold, which applies to neurons with fluctuating inputs. We deduce a differential equation for the threshold, similar to the equations of gating variables in the Hodgkin-Huxley formalism, which describes how the spike threshold varies with the membrane potential, depending on channel properties. We find that spike threshold depends logarithmically on Na channel density, and that Na channel inactivation and K channels can dynamically modulate it in an adaptive way: the threshold increases with membrane potential and after every action potential. Our equation was validated with simulations of a previously published multicompartemental model of spike initiation. Finally, we observed that threshold variability in models depends crucially on the shape of the Na activation function near spike initiation (about −55 mV), while its parameters are adjusted near half-activation voltage (about −30 mV), which might explain why many models exhibit little threshold variability, contrary to experimental observations. We conclude that ionic channels can account for large variations in spike threshold

Lipid-dependent gating of a voltage-gated potassium channel

Recent studies hypothesized that phospholipids stabilize two voltage-sensing arginine residues of certain voltage-gated potassium channels in activated conformations. It remains unclear how lipids directly affect these channels. Here, by examining the conformations of the KvAP in different lipids, we showed that without voltage change, the voltage-sensor domains switched from the activated to the resting state when their surrounding lipids were changed from phospholipids to nonphospholipids. Such lipid-determined conformational change was coupled to the ion-conducting pore, suggesting that parallel to voltage gating, the channel is gated by its annular lipids. Our measurements recognized that the energetic cost of lipid-dependent gating approaches that of voltage gating, but kinetically it appears much slower. Our data support that a channel and its surrounding lipids together constitute a functional unit, and natural nonphospholipids such as cholesterol should exert strong effects on voltage-gated channels. Our first observation of lipid-dependent gating may have general implications to other membrane proteins

Stochastic Ion Channel Gating in Dendritic Neurons: Morphology Dependence and Probabilistic Synaptic Activation of Dendritic Spikes

Neuronal activity is mediated through changes in the probability of stochastic transitions between open and closed states of ion channels. While differences in morphology define neuronal cell types and may underlie neurological disorders, very little is known about influences of stochastic ion channel gating in neurons with complex morphology. We introduce and validate new computational tools that enable efficient generation and simulation of models containing stochastic ion channels distributed across dendritic and axonal membranes. Comparison of five morphologically distinct neuronal cell types reveals that when all simulated neurons contain identical densities of stochastic ion channels, the amplitude of stochastic membrane potential fluctuations differs between cell types and depends on sub-cellular location. For typical neurons, the amplitude of membrane potential fluctuations depends on channel kinetics as well as open probability. Using a detailed model of a hippocampal CA1 pyramidal neuron, we show that when intrinsic ion channels gate stochastically, the probability of initiation of dendritic or somatic spikes by dendritic synaptic input varies continuously between zero and one, whereas when ion channels gate deterministically, the probability is either zero or one. At physiological firing rates, stochastic gating of dendritic ion channels almost completely accounts for probabilistic somatic and dendritic spikes generated by the fully stochastic model. These results suggest that the consequences of stochastic ion channel gating differ globally between neuronal cell-types and locally between neuronal compartments. Whereas dendritic neurons are often assumed to behave deterministically, our simulations suggest that a direct consequence of stochastic gating of intrinsic ion channels is that spike output may instead be a probabilistic function of patterns of synaptic input to dendrites

Tailored ß-Cyclodextrin Blocks the Translocation Pores of Binary Exotoxins from C. Botulinum and C. Perfringens and Protects Cells from Intoxication

International audienceBackgroundClostridium botulinum C2 toxin and Clostridium perfringens iota toxin are binary exotoxins, which ADP-ribosylate actin in the cytosol of mammalian cells and thereby destroy the cytoskeleton. C2 and iota toxin consists of two individual proteins, an enzymatic active (A-) component and a separate receptor binding and translocation (B-) component. The latter forms a complex with the A-component on the surface of target cells and after receptor-mediated endocytosis, it mediates the translocation of the A-component from acidified endosomal vesicles into the cytosol. To this end, the B-components form heptameric pores in endosomal membranes, which serve as translocation channels for the A-components.Here we demonstrate that a 7-fold symmetrical positively charged ß-cyclodextrin derivative, per-6-S-(3-aminomethyl)benzylthio-ß-cyclodextrin, protects cultured cells from intoxication with C2 and iota toxins in a concentration-dependent manner starting at low micromolar concentrations. We discovered that the compound inhibited the pH-dependent membrane translocation of the A-components of both toxins in intact cells. Consistently, the compound strongly blocked transmembrane channels formed by the B-components of C2 and iota toxin in planar lipid bilayers in vitro. With C2 toxin, we consecutively ruled out all other possible inhibitory mechanisms showing that the compound did not interfere with the binding of the toxin to the cells or with the enzyme activity of the A-component.Conclusions/SignificanceThe described ß-cyclodextrin derivative was previously identified as one of the most potent inhibitors of the binary lethal toxin of Bacillus anthracis both in vitro and in vivo, implying that it might represent a broad-spectrum inhibitor of binary pore-forming exotoxins from pathogenic bacteria

FKBP12 Activates the Cardiac Ryanodine Receptor Ca2+-Release Channel and Is Antagonised by FKBP12.6

Changes in FKBP12.6 binding to cardiac ryanodine receptors (RyR2) are implicated in mediating disturbances in Ca2+-homeostasis in heart failure but there is controversy over the functional effects of FKBP12.6 on RyR2 channel gating. We have therefore investigated the effects of FKBP12.6 and another structurally similar molecule, FKBP12, which is far more abundant in heart, on the gating of single sheep RyR2 channels incorporated into planar phospholipid bilayers and on spontaneous waves of Ca2+-induced Ca2+-release in rat isolated permeabilised cardiac cells. We demonstrate that FKBP12 is a high affinity activator of RyR2, sensitising the channel to cytosolic Ca2+, whereas FKBP12.6 has very low efficacy, but can antagonise the effects of FKBP12. Mathematical modelling of the data shows the importance of the relative concentrations of FKBP12 and FKBP12.6 in determining RyR2 activity. Consistent with the single-channel results, physiological concentrations of FKBP12 (3 µM) increased Ca2+-wave frequency and decreased the SR Ca2+-content in cardiac cells. FKBP12.6, itself, had no effect on wave frequency but antagonised the effects of FKBP12

Selection of Inhibitor-Resistant Viral Potassium Channels Identifies a Selectivity Filter Site that Affects Barium and Amantadine Block

BACKGROUND:Understanding the interactions between ion channels and blockers remains an important goal that has implications for delineating the basic mechanisms of ion channel function and for the discovery and development of ion channel directed drugs. METHODOLOGY/PRINCIPAL FINDINGS:We used genetic selection methods to probe the interaction of two ion channel blockers, barium and amantadine, with the miniature viral potassium channel Kcv. Selection for Kcv mutants that were resistant to either blocker identified a mutant bearing multiple changes that was resistant to both. Implementation of a PCR shuffling and backcrossing procedure uncovered that the blocker resistance could be attributed to a single change, T63S, at a position that is likely to form the binding site for the inner ion in the selectivity filter (site 4). A combination of electrophysiological and biochemical assays revealed a distinct difference in the ability of the mutant channel to interact with the blockers. Studies of the analogous mutation in the mammalian inward rectifier Kir2.1 show that the T-->S mutation affects barium block as well as the stability of the conductive state. Comparison of the effects of similar barium resistant mutations in Kcv and Kir2.1 shows that neighboring amino acids in the Kcv selectivity filter affect blocker binding. CONCLUSIONS/SIGNIFICANCE:The data support the idea that permeant ions have an integral role in stabilizing potassium channel structure, suggest that both barium and amantadine act at a similar site, and demonstrate how genetic selections can be used to map blocker binding sites and reveal mechanistic features

Millisecond dynamics of an unlabeled amino acid transporter

Excitatory amino acid transporters (EAATs) are important in many physiological processes and crucial for the removal of excitatory amino acids from the synaptic cleft. Here, we develop and apply high-speed atomic force microscopy line-scanning (HS-AFM-LS) combined with automated state assignment and transition analysis for the determination of transport dynamics of unlabeled membrane-reconstituted GltPh, a prokaryotic EAAT homologue, with millisecond temporal resolution. We find that GltPh transporters can operate much faster than previously reported, with state dwell-times in the 50 ms range, and report the kinetics of an intermediate transport state with height between the outward- and inward-facing states. Transport domains stochastically probe transmembrane motion, and reversible unsuccessful excursions to the intermediate state occur. The presented approach and analysis methodology are generally applicable to study transporter kinetics at system-relevant temporal resolution