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.     

Computational Observation of an Ion Permeation Through a Channel Protein

The ion permeation process, driven by a membrane potential through an outer membrane protein, OmpF porin of Escherichia coli, was simulated by molecular dynamics. A Na⁺ ion, initially placed in the solvent region at the outer side of the porin channel, moved along the electric field passing through the porin channel in a 1.3 nsec simulation; the permeation rate was consistent with the experimentally estimated channel activity (10⁸10⁹/sec). In this simulation, it was indicated that the ion permeation through the porin channel proceeds by a push-out mechanism, and that Asp113 is an important residue for the channel activity.     

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.     

Permeation of water through the KcsA K+ channel

Previous studies have reported that the KcsA potassium channel has an osmotic permeability coefficient of 4.8 x 10(-12) cm3/s, giving it a significantly higher osmotic permeability coefficient than that of some membrane channels specialized in water transport. This high osmotic permeability is proposed to occur when the channel is depleted of potassium ions, the presence of which slow down the water permeation process. The atomic structure of the potassium-depleted KcsA channel and the mechanisms of water permeation have not been well characterized so far. Here, all-atom molecular dynamics simulations, in conjunction with an umbrella sampling strategy and a nonequilibrium approach to simulate pressure gradients are employed to illustrate the permeation of water in the absence of ions through the KcsA K+ channel. Equilibrium molecular dynamics simulations (95 ns combined total length) identified a possible structure of the potassium-depleted KcsA channel, and umbrella sampling calculations (160 ns combined total length) revealed that this structure is not permeable by water molecules moving along the channel axis. The simulation of a pressure gradient across the channel (30 ns combined total length) identified an alternative permeation pathway with a computed osmotic permeability of approximately (2.7 +/- 0.9) x 10(-13) cm3/s. Water fluxes along this pathway did not proceed through collective water motions or transitions to vapor state. All of the major results of this study were robust against variations in a wide set of simulation parameters (force field, water model, membrane model, and channel conformation).     

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.     

Macroscopic Na+ Currents in the “Nonconducting” Shaker Potassium Channel Mutant W434F

C-type inactivation in Shaker potassium channels inhibits K⁺ permeation. The associated structural changes appear to involve the outer region of the pore. Recently, we have shown that C-type inactivation involves a change in the selectivity of the Shaker channel, such that C-type inactivated channels show maintained voltage-sensitive activation and deactivation of Na⁺ and Li⁺ currents in K⁺-free solutions, although they show no measurable ionic currents in physiological solutions. In addition, it appears that the effective block of ion conduction produced by the mutation W434F in the pore region may be associated with permanent C-type inactivation of W434F channels. These conclusions predict that permanently C-type inactivated W434F channels would also show Na⁺ and Li⁺ currents (in K⁺-free solutions) with kinetics similar to those seen in C-type-inactivated Shaker channels. This paper confirms that prediction and demonstrates that activation and deactivation parameters for this mutant can be obtained from macroscopic ionic current measurements. We also show that the prolonged Na⁺ tail currents typical of C-type inactivated channels involve an equivalent prolongation of the return of gating charge, thus demonstrating that the kinetics of gating charge return in W434F channels can be markedly altered by changes in ionic conditions.     

Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process.     

Competing Lipid-Protein and Protein-Protein Interactions Determine Clustering and Gating Patterns in the Potassium Channel from Streptomyces lividans (KcsA)

There is increasing evidence to support the notion that membrane proteins, instead of being isolated components floating in a fluid lipid environment, can be assembled into supramolecular complexes that take part in a variety of cooperative cellular functions. The interplay between lipid-protein and protein-protein interactions is expected to be a determinant factor in the assembly and dynamics of such membrane complexes. Here we report on a role of anionic phospholipids in determining the extent of clustering of KcsA, a model potassium channel. Assembly/disassembly of channel clusters occurs, at least partly, as a consequence of competing lipid-protein and protein-protein interactions at nonannular lipid binding sites on the channel surface and brings about profound changes in the gating properties of the channel. Our results suggest that these latter effects of anionic lipids are mediated via the Trp⁶⁷–Glu⁷¹–Asp⁸⁰ inactivation triad within the channel structure and its bearing on the selectivity filter.     

Switching substrate specificity of AMT/MEP/ Rh proteins

In organisms from all kingdoms of life, ammonia and its conjugated ion ammonium are transported across membranes by proteins of the AMT/Rh family. Efficient and successful growth often depends on sufficient ammonium nutrition. The proteins mediating this transport, the so called Ammonium Transporter (AMT) or Rhesus like (Rh) proteins, share a very similar trimeric overall structure and a high sequence similarity even throughout the kingdoms. Even though structural components of the transport mechanism, like an external substrate recruitment site, an essential twin histidine pore motif, a phenylalanine gate and the hydrophobic pore are strongly conserved and have been analyzed in detail by molecular dynamic simulations and mutational studies, the substrate(s), which pass the central pores of the AMT/Rh subunits, NH₄⁺, NH₃ + H⁺, NH₄⁺ + H⁺ or NH₃, are still a matter of debate for most proteins, including the best characterized AmtB protein from Escherichia coli. The lack of a robust expression system for functional analysis has hampered proof of structural and mutational studies, although the NH₃ transport function for Rh-like proteins is rarely disputed. In plant transporters belonging to the subfamily AMT1, transport is associated with electrical currents, while some plant transporters, notably of the AMT2 type, were suggested to transport NH₃ across the membrane, without associated ionic currents. Here we summarize data in favor of each substrate for the distinct AMT/Rh classes, discuss mutants and how they differ in structure and functionality. A common mechanism with deprotonation and subsequent NH₃ transport through the central subunit pore is suggested.     

Switching substrate specificity of AMT/MEP/ Rh proteins

In organisms from all kingdoms of life, ammonia and its conjugated ion ammonium are transported across membranes by proteins of the AMT/Rh family. Efficient and successful growth often depends on sufficient ammonium nutrition. The proteins mediating this transport, the so called Ammonium Transporter (AMT) or Rhesus like (Rh) proteins, share a very similar trimeric overall structure and a high sequence similarity even throughout the kingdoms. Even though structural components of the transport mechanism, like an external substrate recruitment site, an essential twin histidine pore motif, a phenylalanine gate and the hydrophobic pore are strongly conserved and have been analyzed in detail by molecular dynamic simulations and mutational studies, the substrate(s), which pass the central pores of the AMT/Rh subunits, NH₄⁺, NH₃ + H⁺, NH₄⁺ + H⁺ or NH₃, are still a matter of debate for most proteins, including the best characterized AmtB protein from Escherichia coli. The lack of a robust expression system for functional analysis has hampered proof of structural and mutational studies, although the NH₃ transport function for Rh-like proteins is rarely disputed. In plant transporters belonging to the subfamily AMT1, transport is associated with electrical currents, while some plant transporters, notably of the AMT2 type, were suggested to transport NH₃ across the membrane, without associated ionic currents. Here we summarize data in favor of each substrate for the distinct AMT/Rh classes, discuss mutants and how they differ in structure and functionality. A common mechanism with deprotonation and subsequent NH₃ transport through the central subunit pore is suggested.     

Theoretical and simulation studies on voltage-gated sodium channels

Voltage-gated sodium (Na_v) channels are indispensable membrane elements for the generation and propagation of electric signals in excitable cells. The successes in the crystallographic studies on prokaryotic Na_v channels in recent years greatly promote the mechanistic investigation of these proteins and their eukaryotic counterparts. In this paper, we mainly review the progress in computational studies, especially the simulation studies, on these proteins in the past years.     

Subthreshold Voltage Noise Due to Channel Fluctuations in Active Neuronal Membranes

Voltage-gated ion channels in neuronal membranes fluctuate randomly between different conformational states due to thermal agitation. Fluctuations between conducting and nonconducting states give rise to noisy membrane currents and subthreshold voltage fluctuations and may contribute to variability in spike timing. Here we study subthreshold voltage fluctuations due to active voltage-gated Na⁺ and K⁺ channels as predicted by two commonly used kinetic schemes: the Mainen et al. (1995) (MJHS) kinetic scheme, which has been used to model dendritic channels in cortical neurons, and the classical Hodgkin-Huxley (1952) (HH) kinetic scheme for the squid giant axon. We compute the magnitudes, amplitude distributions, and power spectral densities of the voltage noise in isopotential membrane patches predicted by these kinetic schemes. For both schemes, noise magnitudes increase rapidly with depolarization from rest. Noise is larger for smaller patch areas but is smaller for increased model temperatures. We contrast the results from Monte Carlo simulations of the stochastic nonlinear kinetic schemes with analytical, closed-form expressions derived using passive and quasi-active linear approximations to the kinetic schemes. For all subthreshold voltage ranges, the quasi-active linearized approximation is accurate within 8% and may thus be used in large-scale simulations of realistic neuronal geometries.     

Potassium transport and intracellular potassium activities in rabbit gallbladder

Summary Intracellular K activities, (K) _c , in rabbit gallbladder were determined using conventional and ion-selective microelectrodes. (K) _c averaged 73mm and was 1.5 times that predicted for an equilibrium distribution of the ion across both apical and basolateral membranes. Thus, K must be actively transported into the cell, and the responsible mechanism is almost certainly the Na–K exchange pump in the basolateral membrane.Measurements of the bidirectional transepithelial fluxes of⁴²K indicate that K is secreted into the mucosal solution at a rate of 0.8 eq/cm² hr; this value is only 6% of the rate of transcellular Na absorption by this epithelium.Calculation of the conductance of the basolateral membrane,G ^s, reveals that it is too low to account for the maintenance of the steady-state (K) _c by a 3 Na2 K pump mechanism at the basolateral membrane if K exit across that barrier is entirely electrodiffusional.Our results together with those of others strongly suggest that a significant fraction of downhill K exit from the cell across the basolateral membrane is nonconductive and coupled to the movement of some other ion, perhaps Cl.     

Ion selectivity of epithelial Na channels

Epithelial Na channels are apparently pore-forming membrane proteins which conduct Na much better than any other biologically abundant ion. The conductance to Na can be 100 to 1000 times higher than that to K. The only other ions that can readily get through this channel are protons and Li. Small organic cations cannot pass through the channel, and water may also be impermeant. The selectivity properties of epithelial Na channels appear to be determined by at least three factors: A high field-strength anionic site, most likely a carboxyl residue of glutamic or aspartic acid residues on the channel protein, probably accounts for the high conductance through these channels of Na and Li and to the low conductance of K, Rb and Cs. A restriction in the size of the pore at its narrowest point probably accounts for the low conductance of organic cations as well as the possible exclusion of water molecules. The outer mouth of the channel appears to be negatively charged and may control access to the region of highest selectivity and may serve as a preliminary selectivity filter, attracting cations over anions. These conclusions are illustrated by the cartoon of the channel in Fig. 3. This picture is obviously both fanciful and simplified, but its general points will hopefully be testable. It leaves open a number of important questions, including: does amiloride block the channel by binding within the outer mouth? what does the inner mouth of the channel look like, and does this part of the channel contribute to selectivity? and what, if any, are the interactions between the features of the channel that impart selectivity and those that control the regulation of the channel by hormonal and other factors?     

On the cyclotron resonance mechanism for magnetic field effects on transmembrane ion conductivity

The cyclotron resonance model, recently proposed to account for physiological response to weak environmental magnetic fields, is shown to violate the laws of classical mechanics. Further, it is argued that the ubiquitous presence of dynamic friction in fluid media precludes significant magnetic effects on membrane ion transport.     

驱动蛋白结构与运动机制

驱动蛋白是一类能够利用ATP水解释放的化学能驱动其所携带的“货物”分子沿着微管（microtubule,MT）定向运动的分子马达,在细胞器运输、有丝分裂、轴突运输等方面有着重要的生理作用。随着驱动蛋白结合ADP、ATP和未结合核苷酸（APO）三种特征状态的晶体结构的解析,驱动蛋白构象变化的研究得到了进一步发展,而在力产生机制和运动模型方面仍然存在较大争议。本文以kinesin-1家族为例,分析了驱动蛋白三种特征状态结构的特点、状态结构间的构象转变,论述了驱动蛋白的力产生机制和整个迈步过程。并探讨了驱动蛋白的运动模型,同时采用分子动力学模拟比较了驱动蛋白的两种迈步方式,为深入研究驱动蛋白提供了一定的理论计算。最后,基于本课题组对复杂体系的研究,对驱动蛋白体系的控制机制提出了新的假设,并对未来的研究方向进行了展望。     

Uphill transport of water by electroosmosis

Summary Lepismatidae are able to gain water from subsaturated atmosphere above a relative humidity of 45%, surmounting a water potential difference of at least 1.1×10⁸ Pa (1,100 bar). This extraordinary task is performed by the monolayered epithelium of the posterior rectum. The particle coat of the folded apical membrane of this epithelium suggests the presence of the electrogenic, lumen-directed cation transport, which is commonly found in insects. Assuming this kind of transport and considering the anatomy of the organ, a working hypothesis for this hyposmotic water transport has been developed: The electrogenic cation transport maintains the circulation of the transported ion species across the apical membrane; the voltagedriven inward current transfers water by electroosmosis against its chemical potential from the extracellular space into the cytoplasm. Voltage and current measurements and synchronous measurements of water flow across the epithelium of the posterior rectum ofLepisma saccharina strongly corroborate this hypothesis. The transepithelial voltage is up to 200 mV (lumen positive); the short-circuit current averages 200 A per cm² of the epithelium. Both depend acutely on oxidative metabolism as does spontaneous water uptake. Exogenous transepithelial current (I) induces, independently of anoxia, a proportional change in volume flow (J _v). The induced flow has the direction of the cation flow. Its mean coupling ratio (J _v/I) is 1.5×10^–9m³/A·sec corresponding to 7 to 8 H₂O per positive unit charge. Critical evaluation of experimental data reveals that water uptake by electroosmosis may quantitatively account forin vivo performance without requiring any unusual assumption.