Permeation Models and Structure-Function Relationships in Ion Channels

Recent determination of the molecular structures of potassium andmechanosensitive channels from x-ray crystallography has led to arenewed interest in ion channels. The challenge for permeation modelsis to understand the functional properties of channels from the availablestructural information. Here we give a critical review of the three maincontenders, namely, continuum theories, Brownian dynamics and moleculardynamics. Continuum theories are shown to be invalid in a narrow channel environment because they ignore the self-energy of ions arising from theinduced charges on the dielectric boundary. Brownian and moleculardynamics are thus the only physically valid methods for studying thestructure-function relations in ion channels. Applications of thesemethods to potassium and calcium channels are presented, which illustratethe multi-ion nature of the permeation mechanism in selective biologicalchannels.     

The Persistent Question of Potassium Channel Permeation Mechanisms

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K⁺ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K⁺ ions by the channel. However, the molecular mechanism of K⁺ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K⁺ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.     

Ion exclusion mechanism in aquaporin at an atomistic level

Although the mechanism of proton exclusion in aquaporin is investigated by many researchers, the detailed molecular mechanism for ion exclusion in aquaporin is still not completely understood. In the present work, a detailed mechanism for ion exclusion in aquaporin-1 (AQP1) at an atomistic level is investigated by calculating the free energy for transport of ions in AQP1 using an atomistic molecular dynamics simulation. For this purpose, sodium and chloride ions are chosen as representatives for nonprotonic ions. The simulation shows that the free energy barrier showing its maximum is located at the NPA region for sodium ion while it is located at both the front and the rear for chloride ion and that the barrier height is 18 and 9 kcal/mol, respectively, indicating that the ions are not able to pass through aquaporin. Analysis of the pair interaction energy between the permeating ion and its environment reveals that sodium ion is excluded by the positive charge generated by two alpha-helical macro-dipoles, while chloride ion is expelled by carbonyl oxygen atoms protruding from pore-making residues before it reaches the NPA motif. It is also found that the number of water molecules hydrating the ions is reduced as the ions enter the pore, implying that the energetic cost for detaching water molecules from a permeating ion also contributes to the free energy barriers of ion transport in AQP1.     

Asymmetry and Ion Selectivity Properties of Bacterial Channel NaK Mutants Derived from Ionotropic Glutamate Receptors

Ionotropic glutamate receptors are ligand-gated cation channels that play essential roles in the excitatory synaptic transmission throughout the central nervous system. A number of open-pore structures of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic-acid (AMPA)-type glutamate receptors became available recently by cryo-electron microscopy (cryo-EM). These structures provide valuable insights into the conformation of the selectivity filter (SF), the part of the ion channel that determines the ion selectivity. Nonetheless, due to the moderate resolution of the cryo-EM structures, detailed information such as ion occupancy of monovalent and divalent cations as well as positioning of the side-chains in the SF is still missing. Here, in an attempt to obtain high-resolution information about glutamate receptor SFs, we incorporated partial SF sequences of the AMPA and kainate receptors into the bacterial tetrameric cation channel NaK, which served as a structural scaffold. We determined a series of X-ray structures of NaK-CDI, NaK-SDI and NaK-SELM mutants at 1.42–2.10 Å resolution, showing distinct ion occupation of different monovalent cations. Molecular dynamics (MD) simulations of NaK-CDI indicated the channel to be conductive to monovalent cations, which agrees well with our electrophysiology recordings. Moreover, previously unobserved structural asymmetry of the SF was revealed by the X-ray structures and MD simulations, implying its importance in ion non-selectivity of tetrameric cation channels.     

Physics of Ion Channels

We review the basic physics involved in transport of ions across membrane channels in cells. Electrochemical forces that control the diffusion of ions are discussed both from microscopic and macroscopic perspectives. A case is made for use of Brownian dynamics as the minimal phenomenological model that provides a bridge between experiments and more fundamental theoretical approaches. Application of Brownian and molecular dynamics methods to channels with known molecular structures is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mechanism of Ion Permeation in Skeletal Muscle Chloride Channels

Voltage-gated Cl⁻ channels belonging to the ClC family exhibit unique properties of ion permeation and gating. We functionally probed the conduction pathway of a recombinant human skeletal muscle Cl⁻ channel (hClC-1) expressed both in Xenopus oocytes and in a mammalian cell line by investigating block by extracellular or intracellular I⁻ and related anions. Extracellular and intracellular I⁻ exert blocking actions on hClC-1 currents that are both concentration and voltage dependent. Similar actions were observed for a variety of other halide (Br⁻) and polyatomic (SCN⁻, NO₃ ⁻, CH₃SO₃ ⁻) anions. In addition, I⁻ block is accompanied by gating alterations that differ depending on which side of the membrane the blocker is applied. External I⁻ causes a shift in the voltage-dependent probability that channels exist in three definable kinetic states (fast deactivating, slow deactivating, nondeactivating), while internal I⁻ slows deactivation. These different effects on gating properties can be used to distinguish two functional ion binding sites within the hClC-1 pore. We determined K _D values for I⁻ block in three distinct kinetic states and found that binding of I⁻ to hClC-1 is modulated by the gating state of the channel. Furthermore, estimates of electrical distance for I⁻ binding suggest that conformational changes affecting the two ion binding sites occur during gating transitions. These results have implications for understanding mechanisms of ion selectivity in hClC-1, and for defining the intimate relationship between gating and permeation in ClC channels.     

Ion and pH Sensitivity of a TMBIM Ca2+ Channel

1. Download : Download high-res image (178KB)
2. Download : Download full-size image

     

Ion Channel Selectivity through Stepwise Changes in Binding Affinity

Voltage-gated Ca²⁺ channels select Ca²⁺ over competing, more abundant ions by means of a high affinity binding site in the pore. The maximum off rate from this site is ∼1,000× slower than observed Ca²⁺ current. Various theories that explain how high Ca²⁺ current can pass through such a sticky pore all assume that flux occurs from a condition in which the pore''s affinity for Ca²⁺ transiently decreases because of ion interactions. Here, we use rate theory calculations to demonstrate a different mechanism that requires no transient changes in affinity to quantitatively reproduce observed Ca²⁺ channel behavior. The model pore has a single high affinity Ca²⁺ binding site flanked by a low affinity site on either side; ions permeate in single file without repulsive interactions. The low affinity sites provide steps of potential energy that speed the exit of a Ca²⁺ ion off the selectivity site, just as potential energy steps accelerate other chemical reactions. The steps could be provided by weak binding in the nonselective vestibules that appear to be a general feature of ion channels, by specific protein structures in a long pore, or by stepwise rehydration of a permeating ion. The previous ion-interaction models and this stepwise permeation model demonstrate two general mechanisms, which might well work together, to simultaneously generate high flux and high selectivity in single file pores.     

Recent advances in ion channel research

The field of ion channels has entered into a rapid phase of development in the last few years, partly due to the breakthroughs in determination of the crystal structures of membrane proteins and advances in computer simulations of biomolecules. These advances have finally enabled the long-dreamed goal of relating function of a channel to its underlying molecular structure. Here we present simplified accounts of the competing permeation theories and then discuss their application to the potassium, gramicidin A and calcium channels.     

High-Sensitivity Fluorometry to Resolve Ion Channel Conformational Dynamics

1. Download high-res image (100KB)
2. Download full-size image

     

DICE：电压门控离子通道电生理实验数据库

针对现有相关文献中离子通道电生理数据繁多且分散的特点,开发了一套电压门控离子通道电生理实验数据库。数据库中目前主要包括钠离子通道序列数据、调制剂分子结构和序列数据,并收集整理了文献中调制剂和通道相互作用时的电生理学数据和药理学数据。系统实现了数据的收集、录入、存储和查询,为后期进行数据挖掘奠定了基础。用户可以通过网址http：／／biodb．sgst．cn／DICE对数据库进行访问。     

The passage of an ion through a synthetic ion channel

The simulated system consisted of a fatty acid bilayer membrane dividing two electrolyte layers each containing ions, and a channel composed of linked 15-crown-5 ether rings. The Na⁺ and F^?  ions in the aqueous electrolyte layers were too large to enter the channel, but the Li⁺ ions entered and were transported. Conditions that optimised the passive, electric-field-induced transport of Li⁺ ions through the channel were investigated. It was calculated and rationalised that the higher the numerical value of the electrostatic charge on the oxygen atoms of the crown ether rings, the more easily does the channel convey the Li⁺ ions.     

环核苷酸门控离子通道门控的分子机理

环核苷酸门控离子通道(CNG)最广泛地分布于神经细胞。近年来关于 CNG 通道门控的分子机制的研究取得了很大的进步。研究表明, CNG 通道的组成及组装影响通道的特性及门控。近年来有关 CNG 突变体的研究及半胱氨酸残基亲和性的分析表明, 环核苷酸首先结合到 CNG 通道 C 端的环核苷酸结合域(CNBD)上引起 CNBD 空间构像改变, 然后 4 个亚单元发生空间构像的协调改变, CNG 通道开放。本文详细讨论了 CNG 通道的门控机制、各亚单元之间的相互作用、组装的过程及其空间构想的变化, 为 CNG 通道的进一步研究, 尤其是离子通道疾病方面提供理论指导。     

Blockade of a Retinal cGMP-gated Channel by Polyamines

The cyclic nucleotide–gated (CNG) channel in retinal rods converts the light-regulated intracellular cGMP concentration to various levels of membrane potential. Blockade of the channel by cations such as Ca²⁺ and Mg²⁺ lowers its effective conductance. Consequently, the membrane potential has very low noise, which enables rods to detect light with extremely high sensitivity. Here, we report that three polyamines (putrescine, spermidine, and spermine), which exist in both the intracellular and extracellular media, also effectively block the CNG channel from both sides of the membrane. Among them, spermine has the greatest potency. Extracellular spermine blocks the channel as a permeant blocker, whereas intracellular spermine appears to block the channel in two conformations—one permeant, and the other non- (or much less) permeant. The membrane potential in rods is typically depolarized to approximately −40 mV in the dark. At this voltage, K _1/2 of the CNG channel for extracellular spermine is 3 μM, which is 100–1,000-fold higher affinity than that of the NMDA receptor-channel for extracellular spermine. Blockade of the CNG channel by polyamines may play an important role in suppressing noise in the signal transduction system in rods.     

Structural Insights into a Plant Mechanosensitive Ion Channel MSL1

1. Download : Download high-res image (262KB)
2. Download : Download full-size image

     

Phosphorylation of the Rat Skeletal Muscle Sodium Channel by Cyclic AMP-Dependent Protein Kinase

Cyclic AMP-dependent phosphorylation of the rat brain sodium channel was reported to be restricted to five sites within an approximately 210 amino acid region of the primary sequence that is deleted in the homologous sodium channel from rat skeletal muscle. We find that, in spite of this deletion, the rat muscle sodium channel alpha-subunit is also an excellent substrate for phosphorylation by this kinase both in primary muscle cells in tissue culture and in vitro after isolation from adult muscle. Sodium channel protein purified from adult rat skeletal muscle was readily phosphorylated in vitro by the catalytic subunit of the bovine cyclic AMP-dependent protein kinase (PKa). Only the 260,000 MW alpha-subunit was labeled, with a maximum level of incorporation in vitro of approximately 0.5 mol [32P]phosphate per mole of channel protein. The beta-subunit of the channel is not phosphorylated under these conditions. In primary rat skeletal muscle cells in culture, incorporation of phosphate into the channel alpha-subunit is stimulated 1.3- to 1.5-fold by treatment of the cells with forskolin. Phosphorylation of the sodium channel isolated from these cells could also be demonstrated in vitro using PKa. This in vitro phosphorylation could be inhibited 80-90% by pretreatment of the cells in culture with forskolin, suggesting that the sites labeled in vitro by PKa were the same as those phosphorylated in the intact cells by the endogenous cyclic AMP-dependent kinase. In both the adult muscle channel and the channel from muscle cells in culture, phosphorylation by PKa was limited to serine residues.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure and ion channel activity of the human respiratory syncytial virus (hRSV) small hydrophobic protein transmembrane domain

The small hydrophobic (SH) protein from the human respiratory syncytial virus (hRSV) is a glycoprotein of approximately 64 amino acids with one putative alpha-helical transmembrane domain. Although SH protein is important for viral infectivity, its exact role during viral infection is not clear. Herein, we have studied the secondary structure, orientation, and oligomerization of the transmembrane domain of SH (SH-TM) in the presence of lipid bilayers. Only one oligomer, a pentamer, was observed in PFO-PAGE. Using polarized attenuated total reflection-Fourier transform infrared (PATR-FTIR) spectroscopy, we show that the SH-TM is alpha-helical. The rotational orientation of SH-TM was determined by site-specific infrared dichroism (SSID) at two consecutive isotopically labeled residues. This orientation is consistent with that of an evolutionary conserved pentameric model obtained from a global search protocol using 13 homologous sequences of RSV. Conductance studies of SH-TM indicate ion channel activity, which is cation selective, and inactive below the predicted pK(a) of histidine. Thus, our results provide experimental evidence that the transmembrane domain of SH protein forms pentameric alpha-helical bundles that form cation-selective ion channels in planar lipid bilayers. We provide a model for this pore, which should be useful in mutagenesis studies to elucidate its role during the virus cycle.     

Protein dynamics in a family of laboratory evolved thermophilic enzymes

Molecular dynamics simulations were employed to study how protein solution structure and dynamics are affected by adaptation to high temperature. Simulations were carried out on a para-nitrobenzyl esterase (484 residues) and two thermostable variants that were generated by laboratory evolution. Although these variants display much higher melting temperatures than wild-type (up to 18 degrees C higher) they are both >97% identical in sequence to the wild-type. In simulations at 300 K the thermostable variants remain closer to their crystal structures than wild-type. However, they also display increased fluctuations about their time-averaged structures. Additionally, both variants show a small but significant increase in radius of gyration relative to wild-type. The vibrational density of states was calculated for each of the esterases. While the density of states profiles are similar overall, both thermostable mutants show increased populations of the very lowest frequency modes (<10 cm(-1)), with the more stable mutant showing the larger increase. This indicates that the thermally stable variants experience increased concerted motions relative to wild-type. Taken together, these data suggest that adaptation for high temperature stability has resulted in a restriction of large deviations from the native state and a corresponding increase in smaller scale fluctuations about the native state. These fluctuations contribute to entropy and hence to the stability of the native state. The largest changes in localized dynamics occur in surface loops, while other regions, particularly the active site residues, remain essentially unchanged. Several mutations, most notably L313F and H322Y in variant 8G8, are in the region showing the largest increase in fluctuations, suggesting that these mutations confer more flexibility to the loops. As a validation of our simulations, the fluctuations of Trp102 were examined in detail, and compared with Trp102 phosphorescence lifetimes that were previously measured. Consistent with expectations from the theory of phosphorescence, an inverse correlation between out-of-plane fluctuations on the picosecond time scale and phosphorescence lifetime was observed.     

Sodium Ions Do Not Stabilize the Selectivity Filter of a Potassium Channel

Ion conduction is an essential function for electrical activity in all organisms. The non-selective ion channel NaK was previously shown to adopt two stable conformations of the selectivity filter. Here, we present solid-state NMR measurements of NaK demonstrating a population shift between these conformations induced by changing the ions in the sample while the overall structure of NaK is not affected. We show that two K⁺-selective mutants (NaK2K and NaK2K-Y66F) suffer a complete loss of selectivity filter stability under Na⁺ conditions, but do not collapse into a defined structure. Widespread chemical shift perturbations are seen between the Na⁺ and K⁺ states of the K⁺-selective mutants in the region of the pore helix indicating structural changes. We conclude that the stronger link between the selectivity filter and the pore helix in the K⁺-selective mutants, compared to the non-selective wild-type NaK channel, reduces the ion-dependent conformational flexibility of the selectivity filter.