Physical origin of selectivity in ionic channels of biological membranes

Proton Coulomb Blockade Effect Involving Covalent Oxygen-Hydrogen Bond Switching

Instead of the canonical Grotthuss mechanism, we show that a knock-on proton transport process is preferred between organic functional groups (e.g., -COOH and -OH) and adjacent water molecules in biological proton channel and synthetic nanopores through comprehensive quantum and classical molecular dynamics simulations. The knock-on process is accomplished by the switching of covalent O─H bonds of the functional group under externally applied electric fields. The proton transport through the synthetic nanopore exhibits nonlinear current-voltage characteristics, suggesting an unprecedented proton Coulomb blockade effect. These findings not only enhance the understanding of proton transport in nanoconfined systems but also pave the way for the design of a variety of proton-based nanofluidic devices.

Highly Selective Lithium Transport through Crown Ether Pillared Angstrom Channels

Biological ion channels use the synergistic effects of various strategies to realize highly selective ion sieving. For example, potassium channels use functional groups and angstrom-sized pores to discriminate rival ions and enrich target ions. Inspired by this, we constructed a layered crystal pillared by crown ether that incorporates these strategies to realize high Li+ selectivity. The pillared channels and crown ether have an angstrom-scale size. The crown ether specifically allows the low-barrier transport of Li+ . The channels attract and enrich Li+ ions by up to orders of magnitude. As a result, our material sieves Li+ out of various common ions such as Na+ , K+ , Ca2+ , Mg2+ and Al3+ . Moreover, by spontaneously enriching Li+ ions, it realizes an effective Li+ /Na+ selectivity of 1422 in artificial seawater where the Li+ concentration is merely 25 μM. We expect this work to spark technologies for the extraction of lithium and other dilute metal ions.

SARS-CoV spike proteins can compete for electrolytes in physiological fluids according to structure-based quantum-chemical calculations

The trimeric spike (S) glycoprotein is the trojan horse and the stronghold of the severe acute respiratory syndrome coronaviruses. Although several structures of the S-protein have been solved, a complete understanding of all its functions is still lacking. Our multi-approach study, based on the combination of structural experimental data and quantum-chemical DFT calculations, led to identify a sequestration site for sodium, potassium and chloride ions within the central cavity of both the SARS-CoV-1 and SARS-CoV-2 spike proteins. The same region was found as strictly conserved, even among the sequences of the bat-respective coronaviruses. Due to the prominent role of the main three electrolytes at many levels, and their possible implication in the molecular mechanisms of COVID-19 disease, our study can take the lead in important discoveries related to the SARS-CoV-2 biology, as well as in the design of novel effective therapeutic strategies.

Adaptable and Multifunctional Ion-Conducting Aquaporins

Aquaporins function as water and neutral solute channels, signaling hubs, disease virulence factors, and metabolon components. We consider plant aquaporins that transport ions compared to some animal counterparts. These are candidates for important, as yet unidentified, cation and anion channels in plasma, tonoplast, and symbiotic membranes. For those individual isoforms that transport ions, water, and gases, the permeability spans 12 orders of magnitude. This requires tight regulation of selectivity via protein interactions and posttranslational modifications. A phosphorylation-dependent switch between ion and water permeation in AtPIP2;1 might be explained by coupling between the gates of the four monomer water channels and the central pore of the tetramer. We consider the potential for coupling between ion and water fluxes that could form the basis of an electroosmotic transducer. A grand challenge in understanding the roles of ion transporting aquaporins is their multifunctional modes that are dependent on location, stress, time, and development.

Ionic Coulomb blockade and the determinants of selectivity in the NaChBac bacterial sodium channel

Mutation-induced transformations of conductivity and selectivity in NaChBac bacterial channels are studied experimentally and interpreted within the framework of ionic Coulomb blockade (ICB), while also taking account of resonant quantised dehydration (QD) and site protonation. Site-directed mutagenesis and whole-cell patch-clamp experiments are used to investigate how the fixed charge Q_f at the selectivity filter (SF) affects both valence selectivity and same-charge selectivity. The new ICB/QD model predicts that increasing ∣Q_f∣ should lead to a shift in selectivity sequences toward larger ion sizes, in agreement with the present experiments and with earlier work. Comparison of the model with experimental data leads to the introduction of an effective charge Q_f^∗ at the SF, which was found to differ between Aspartate and Glutamate charged rings, and also to depend on position within the SF. It is suggested that protonation of the residues within the restricted space of the SF is important in significantly reducing the effective charge of the EEEE ring. Values of Q_f^∗ derived from experiments on divalent blockade agree well with expectations based on the ICB/QD model and have led to the first demonstration of ICB oscillations in Ca²⁺ conduction as a function of the fixed charge. Preliminary studies of the dependence of Ca²⁺ conduction on pH are qualitatively consistent with the predictions of the model.

Design principles of ion selective nanostructured membranes for the extraction of lithium ions

It is predicted that the continuously increasing demand for the energy-critical element of lithium will soon exceed its availability, rendering it a geopolitically significant resource. The present work critically reviews recent reports on Li⁺ selective membranes. Particular emphasis has been placed on the basic principles of the materials' design for the development of membranes with nanochannels and nanopores with Li⁺ selectivity. Fundamental and practical challenges, as well as prospects for the targeted design of Li⁺ ion-selective membranes are also presented, with the goal of inspiring future critical research efforts in this scientifically and strategically important field.

The permeation mechanism of organic cations through a CNG mimic channel

Several channels, ranging from TRP receptors to Gap junctions, allow the exchange of small organic solute across cell membrane. However, very little is known about the molecular mechanism of their permeation. Cyclic Nucleotide Gated (CNG) channels, despite their homology with K+ channels and in contrast with them, allow the passage of larger methylated and ethylated ammonium ions like dimethylammonium (DMA) and ethylammonium (EA). We combined electrophysiology and molecular dynamics simulations to examine how DMA interacts with the pore and permeates through it. Due to the presence of hydrophobic groups, DMA enters easily in the channel and, unlike the alkali cations, does not need to cross any barrier. We also show that while the crystal structure is consistent with the presence of a single DMA ion at full occupancy, the channel is able to conduct a sizable current of DMA ions only when two ions are present inside the channel. Moreover, the second DMA ion dramatically changes the free energy landscape, destabilizing the crystallographic binding site and lowering by almost 25 kJ/mol the binding affinity between DMA and the channel. Based on the results of the simulation the experimental electron density maps can be re-interpreted with the presence of a second ion at lower occupancy. In this mechanism the flexibility of the channel plays a key role, extending the classical multi-ion permeation paradigm in which conductance is enhanced by the plain interaction between the ions.

Structural Heterogeneity of CNGA1 Channels Revealed by Electrophysiology and Single-Molecule Force Spectroscopy

The determination at atomic resolution of the three-dimensional molecular structure of membrane proteins such as receptors and several ion channels has been a major breakthrough in structural biology. The molecular structure of several members of the superfamily of voltage-gated ionic channels such as K⁺ and Na⁺ is now available. However, despite several attempts, the molecular structure at atomic resolution of the full cyclic nucleotide-gated (CNG) ion channel, although a member of the same superfamily of voltage-gated ion channels, has not been obtained yet, neither by X-ray crystallography nor by electron cryomicroscopy (cryo-EM). It is possible that CNG channels have a high structural heterogeneity, making difficult crystallization and single-particle analysis. To address this issue, we have combined single-molecule force spectroscopy (SMFS) and electrophysiological experiments to characterize the structural heterogeneity of CNGA1 channels expressed in Xenopus laevis oocytes. The unfolding of the cytoplasmic domain had force peaks, occurring with a probability from 0.2 to 0.96. Force peaks during the unfolding of the transmembrane domain had a probability close to 1, but the distribution of the increase in contour length between two successive force peaks had multiple maxima differing by tens of nanometers. Concomitant electrophysiological experiments showed that the rundown in mutant channels S399C is highly variable and that the effect of thiol reagents when specific residues were mutated was consistent with a dynamic structural heterogeneity. These results show that CNGA1 channels have a wide spectrum of native conformations that are difficult to detect with X-ray crystallography and cryo-EM.

Heterogeneous sub-continuum ionic transport in statistically isolated graphene nanopores

Graphene and other two-dimensional materials offer a new class of ultrathin membranes that can have atomically defined nanopores with diameters approaching those of hydrated ions. These nanopores have the smallest possible pore volumes of any ion channel, which, due to ionic dehydration and electrokinetic effects, places them in a novel transport regime and allows membranes to be created that combine selective ionic transport with ultimate permeance and could lead to separations and sensing applications. However, experimental characterization and understanding of sub-continuum ionic transport in nanopores below 2 nm is limited. Here we show that isolated sub-2 nm pores in graphene exhibit, in contrast to larger pores, diverse transport behaviours consistent with ion transport over a free-energy barrier arising from ion dehydration and electrostatic interactions. Current-voltage measurements reveal that the conductance of graphene nanopores spans three orders of magnitude and that they display distinct linear, voltage-activated or rectified current-voltage characteristics and different cation-selectivity profiles. In rare cases, rapid, voltage-dependent stochastic switching is observed, consistent with the presence of a dissociable group in the pore vicinity. A modified Nernst-Planck model incorporating ion hydration and electrostatic effects quantitatively matches the observed behaviours.

A structural, functional, and computational analysis suggests pore flexibility as the base for the poor selectivity of CNG channels

Cyclic nucleotide-gated (CNG) ion channels, despite a significant homology with the highly selective K(+) channels, do not discriminate among monovalent alkali cations and are permeable also to several organic cations. We combined electrophysiology, molecular dynamics (MD) simulations, and X-ray crystallography to demonstrate that the pore of CNG channels is highly flexible. When a CNG mimic is crystallized in the presence of a variety of monovalent cations, including Na(+), Cs(+), and dimethylammonium (DMA(+)), the side chain of Glu66 in the selectivity filter shows multiple conformations and the diameter of the pore changes significantly. MD simulations indicate that Glu66 and the prolines in the outer vestibule undergo large fluctuations, which are modulated by the ionic species and the voltage. This flexibility underlies the coupling between gating and permeation and the poor ionic selectivity of CNG channels.

Energetics of discrete selectivity bands and mutation-induced transitions in the calcium-sodium ion channels family

We use Brownian dynamics (BD) simulations to study the ionic conduction and valence selectivity of a generic electrostatic model of a biological ion channel as functions of the fixed charge Q(f) at its selectivity filter. We are thus able to reconcile the discrete calcium conduction bands recently revealed in our BD simulations, M0 (Q(f)=1e), M1 (3e), M2 (5e), with a set of sodium conduction bands L0 (0.5e), L1 (1.5e), thereby obtaining a completed pattern of conduction and selectivity bands vs Q(f) for the sodium-calcium channels family. An increase of Q(f) leads to an increase of calcium selectivity: L0 (sodium-selective, nonblocking channel) → M0 (nonselective channel) → L1 (sodium-selective channel with divalent block) → M1 (calcium-selective channel exhibiting the anomalous mole fraction effect). We create a consistent identification scheme where the L0 band is putatively identified with the eukaryotic sodium channel The scheme created is able to account for the experimentally observed mutation-induced transformations between nonselective channels, sodium-selective channels, and calcium-selective channels, which we interpret as transitions between different rows of the identification table. By considering the potential energy changes during permeation, we show explicitly that the multi-ion conduction bands of calcium and sodium channels arise as the result of resonant barrierless conduction. The pattern of periodic conduction bands is explained on the basis of sequential neutralization taking account of self-energy, as Q(f)(z,i)=ze(1/2+i), where i is the order of the band and z is the valence of the ion. Our results confirm the crucial influence of electrostatic interactions on conduction and on the Ca(2+)/Na(+) valence selectivity of calcium and sodium ion channels. The model and results could be also applicable to biomimetic nanopores with charged walls.

Competition among Ca2+, Mg2+, and Na+ for model ion channel selectivity filters: determinants of ion selectivity

Because voltage-gated ion channels play critical biological roles, understanding how they can discriminate the native metal ion from rival cations in the milieu is of great interest. Although Ca(2+), Mg(2+), and Na(+) are present in comparable concentrations outside the cell, the factors governing the competition among these cations for the selectivity filter of voltage-gated Ca(2+) ion channel remain unclear. Using density functional theory combined with continuum dielectric methods, we evaluate the effect of (1) the number, chemical type, and charge of the ligands lining the pore, (2) the pore's rigidity, size, symmetry, and solvent accessibility, and (3) the Ca(2+) hydration number outside the selectivity filter on the competition among Ca(2+), Mg(2+), and Na(+) in model selectivity filters. The calculations show how the outcome of this competition depends on the interplay between electronic and solvation effects. Selectivity for monovalent Na(+) over divalent Ca(2+)/Mg(2+) is achieved when solvation effects outweigh electrostatic effects; thus filters comprising a few weak charge-donating groups such as Ser/Thr side chains, where electrostatic effects are relatively weak and are easily overcome by solvation effects, are Na(+)-selective. In contrast, selectivity for divalent Ca(2+)/Mg(2+) over monovalent Na(+) is achieved when metal-ligand electrostatic effects outweigh solvation effects. The key differences in selectivity between Mg(2+) and Ca(2+) lie in the pore size, oligomericity, and solvent accessibility. The results, which are consistent with available experimental data, reveal how the structure and composition of the ion channel selectivity pore had adapted to the specific physicochemical properties of the native metal ion to enhance the competitiveness of the native metal toward rival cations.

Gating of cyclic nucleotide-gated channels is voltage dependent

Cyclic nucleotide-gated channels belong to the family of voltage-gated ion channels, but pore opening requires the presence of intracellular cyclic nucleotides. In the presence of a saturating agonist, cyclic nucleotide-gated channel gating is voltage independent and it is not known why cyclic nucleotide-gated channels are voltage-insensitive despite harbouring the S4-type voltage sensor. Here we report that, in the presence of Li(+), Na(+) and K(+), the gating of wild-type cyclic nucleotide-gated A1 and native cyclic nucleotide-gated channels is voltage independent, whereas their gating is highly voltage-dependent in the presence of Rb(+), Cs(+) and organic cations. Mutagenesis experiments show that voltage sensing occurs through a voltage sensor composed of charged/polar residues in the pore and of the S4-type voltage sensor. During evolution, cyclic nucleotide-gated channels lose their voltage-sensing ability when Na(+) or K(+) permeate so that the vertebrate photoreceptor cyclic nucleotide-gated channels are open at negative voltages, a necessary condition for phototransduction.

Gating in CNGA1 channels

The aminoacid sequences of CNG and K(+) channels share a significant sequence identity, and it has been suggested that these channels have a common ancestral 3D architecture. However, K(+) and CNG channels have profoundly different physiological properties: indeed, K(+) channels have a high ionic selectivity, their gating strongly depends on membrane voltage and when opened by a steady depolarizing voltage several K(+) channels inactivate, whereas CNG channels have a low ion selectivity, their gating is poorly voltage dependent, and they do not desensitize in the presence of a steady concentration of cyclic nucleotides that cause their opening. The purpose of the present review is to summarize and recapitulate functional and structural differences between K(+) and CNG channels with the aim to understand the gating mechanisms of CNG channels.

Determinants of K+ vs Na+ selectivity in potassium channels

Ion channels, specialized pore-forming proteins, are an indispensable component of the nervous system and play a crucial role in regulating cardiac, skeletal, and smooth muscle contraction. Potassium ion channels, controlling the action potential of a number of excitable cells, are characterized by a remarkable ability to select K(+) over Na(+). Although the molecular basis for this striking ion selectivity has been a subject of extensive investigations using both experimental and theoretical methods, the following outstanding questions remain: (a) To what extent is the number of water molecules bound to the permeating ion (i.e., the hydration number) important for the K(+)/Na(+) competition? (b) Are the chemical type and number of coordinating groups lining the pore critical for the selectivity process? (c) Apart from providing cation-ligating groups, do the channel walls play any other role in the selectivity process? This work reveals that the pore's selectivity for K(+) over Na(+) increases with (i) increasing hydration number of K(+) relative to that of Na(+), (ii) increasing number of K(+)-coordinating dipoles, (iii) increasing number of Na(+)-coordinating dipoles, and (iv) decreasing magnitude of the coordinating dipoles provided by the pore. Thus, a high K(+)/Na(+) selectivity in K(+) channels could be achieved from a combination of several favorable factors involving the native ion, the metal-coordinating ligands, and the protein matrix, viz., (a) an octahydrated permeating K(+), (b) a pore lined with 8 carbonyl ligands, and (c) finely tuned physicomechanical properties of the channel walls providing a low dielectric medium favoring a high hydration number for the permeating K(+) and enough stiffness to force the competing Na(+) to adopt an unfavorable 8-fold coordination. This implies that optimal K(+)/Na(+) selectivity in K(+) channels generally does not arise from solely structural or energetic consideration. The factors affecting ion selectivity revealed herein help to rationalize why valinomycin and the KcsA ion channels are highly K(+)-selective, whereas the NaK channel is nonselective. The calculations predict that other pores containing a different number/chemical type of coordinating groups from those observed in potassium channels could also select K(+) over Na(+).

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

Effects of exogenous electromagnetic fields on a simplified ion channel model

In this paper, we calculate the effect of an exogenous perturbation (an electromagnetic field [EMF] oscillating in the range of microwave frequencies in the range of 1 GHz) on the flux of two ion species through a cylindrical ion channel, implementing a continuous model, the Poisson-Smoluchowski system of equations, to study the dynamics of charged particle density inside the channel. The method was validated through comparison with Brownian dynamics simulations, supposed to be more accurate but computationally more demanding, obtaining a very good agreement. No EMF effects were observed for low field intensities below the level for thermal effects, as the highly viscous regime and the simplicity of the channel do not exhibit resonance phenomena. For high intensities of the external field (>10(5) V/m), we observed slightly different behavior of ion concentration oscillations and ion currents as a function of EMF orientation with respect to the channel axis.

Hydration of potassiated amino acids in the gas phase

The thermochemistry of stepwise hydration of several potassiated amino acids was studied by measuring the gas-phase equilibria, AAK(+)(H(2)O)(n-1) + H(2)O = AAK(+)(H(2)O)(n) (AA = Gly, AL, Val, Met, Pro, and Phe), using a high-pressure mass spectrometer. The AAK(+) ions were obtained by electrospray and the equilibrium constants K(n-1,n) were measured in a pulsed reaction chamber at 10 mbar bath gas, N(2), containing a known partial pressure of water vapor. Determination of the equilibrium constants at different temperatures was used to obtain the DeltaH(n)(o), DeltaS(n)(o), and DeltaG(n)(o) values. The results indicate that the water binding energy in AAK(+)(H(2)O) decreases as the K(+) affinity to AA increases. This trend in binding energies is explained in terms of changes in the side-chain substituent, which delocalize the positive charge from K(+) to AA in AAK(+) complexes, varying the AAK(+)-H(2)O electrostatic interaction.

How the asymmetry of internal potential influences the shape of I-V characteristic of nanochannels

Ion transport in biological and synthetic nanochannels is characterized by such phenomena as ion current fluctuations, rectification, and pumping. Recently, it has been shown that the nanofabricated synthetic pores could be considered as analogous to biological channels with respect to their transport characteristics [P. Yu. Apel et al., Nucl. Instrum. Methods Phys. Res. B 184, 337 (2001); Z. Siwy et al., Europhys. Lett. 60, 349 (2002)]. The ion current rectification is analyzed. Ion transport through cylindrical nanopores is described by the Smoluchowski equation. The model is considering the symmetric nanopore with asymmetric charge distribution. In this model, the current rectification in asymmetrically charged nanochannels shows a diodelike shape of I-V characteristic. It is shown that this feature may be induced by the coupling between the degree of asymmetry and the depth of internal electric potential well. The role of concentration gradient is discussed.

Selectivity principle of the ligand escape process from a two-gate tunnel in myoglobin: molecular dynamics simulation

We proposed a selectivity principle for the ligand escape process from two fluctuating bottlenecks in a cavity with a multigate inside a myoglobin pocket. Our previous analytical theory proposed a fluctuating bottleneck model for a Brownian particle passing through two gates on a cavity surface of an enzyme protein and has determined the escape rate in terms of the time-dependent gate function and the competition effect. It illustrated that with two (or more than two) gates on a cavity surface the gate modulation, which is controlled by protein fluctuation, dominates the ligand escape pathway. We have performed a molecular dynamics simulation to investigate the selectivity principle of the ligand escape process from two-gate tunnel in myoglobin. The simulation results confirm our theoretical conjecture. It indicates that the escape process is actually entropy driven, and the ligand escape pathway is chosen via the gate modulation. This suggests an interesting intrinsic property, that is, the oxymyoglobin tertiary structure is favorable to the departure of the ligand from one direction rather than through a biased random walk.

Alkali Cation Binding and Permeation in the Rat Organic Cation Transporter rOCT2

Organic cation transporters of the OCT family mediate downhill transport of organic cations, compatible with carrier, pore, or gate-lumen-gate mechanisms. We studied rat OCT2 expressed in Xenopus oocytes by the two-electrode voltage-clamp technique, including membrane capacitance (C(m)) monitoring. Choline, a transported cationic substrate, elicited the expected inward currents but also elicited decreases of C(m). Similar C(m) decreases were caused by the non-transported inhibitors tetrabutylammonium (a cation) and corticosterone (uncharged). Effects on C(m) were voltage-dependent, with a maximum at -140 mV. These findings suggest that the empty rOCT2 protein can undergo an electrogenic conformation change, with one conformation highly favored at physiological voltage. Moreover, alkali cations elicited considerable inward currents and inhibited uptake of [(14)C]tetraethylammonium with a sequence Cs(+) > Rb(+) > K(+) > Na(+) approximately Li(+). Cs(+) affected current and capacitance with similar affinity (K(0.5) approximately 50 mm). Tetraethylammonium inhibited Cs(+) currents in a concentration-dependent manner. Conversely, Cs(+) inhibited tetraethylammonium uptake by a competitive mechanism. Activation energy of the currents estimated from measurements between 12 degrees C and 32 degrees C was approximately 81 kJ/mol for Cs(+) and 39 kJ/mol for tetramethylammonium, compatible with permeation of Cs(+) through rOCT2 along the same path as organic substrates and by a mechanism different from simple electrodiffusion. Rationalization of Cs(+) selectivity in terms of a pore pointed to a pore diameter of approximately 4 A. Intriguingly, that value matches the known selectivity of rOCT2 for organic compounds. Our data show that selective permeability of rOCT2 is not determined by ligand affinity but might rather be understood in terms of the ion channel concept of a distinct "selectivity filter."

A homology model of the pore region of HCN channels

HCN channels are activated by membrane hyperpolarization and regulated by cyclic nucleotides, such as cyclic adenosine-mono-phosphate (cAMP). Here we present structural models of the pore region of these channels obtained by using homology modeling and validated against spatial constraints derived from electrophysiological experiments. For the construction of the models we make two major assumptions, justified by electrophysiological observations: i), in the closed state, the topology of the inner pore of HCN channels is similar to that of K(+) channels. In particular, the orientation of the S5 and S6 helices of HCN channels is very similar to that of the corresponding helices of the K(+) KcsA and K(+) KirBac1.1 channels. Thus, we use as templates the x-ray structure of these K(+) channels. ii), In the open state, the S6 helix is bent further than it is in the closed state, as suggested (but not proven) by experimental data. For this reason, the template of the open conformation is the x-ray structure of the MthK channel. The structural models of the closed state turn out to be consistent with all the available electrophysiological data. The model of the open state turned out to be consistent with all the available electrophysiological data in the filter region, including additional experimental data performed in this work. However, it required the introduction of an appropriate, experimentally derived constraint for the S6 helix. Our modeling provides a structural framework for understanding several functional properties of HCN channels: i), the cysteine ring at the inner mouth of the pore may act as a sensor of the intracellular oxidizing/reducing conditions; ii), the bending amplitude of the S6 helix upon gating appears to be significantly smaller than that found in MthK channels; iii), the reduced ionic selectivity of HCN channels, relative to that of K(+) channels, may be caused, at least in part, by the larger flexibility of the inner pore of HCN channels.

Structural basis of gating of CNG channels

Cyclic nucleotide-gated (CNG) ion channels, underlying sensory transduction in vertebrate photoreceptors and olfactory sensory neurons, require cyclic nucleotides to open. Here, we present structural models of the tetrameric CNG channel pore from bovine rod in both open and closed states, as obtained by combining homology modeling-based techniques, experimentally derived spatial constraints and structural patterns present in the PDB database. Gating is initiated by an anticlockwise rotation of the N-terminal region of the C-linker, which is then, transmitted through the S6 transmembrane helices to the P-helix, and in turn from this to the pore lumen, which opens up from 2 to 5A thus allowing for ion permeation. The approach, here presented, is expected to provide a general methodology for model ion channels and their gating when structural templates are available and an extensive electrophysiological analysis has been performed.

Modulation of plant ion channels by oxidizing and reducing agents

Ion channels are proteins forming hydrophilic pathways through the membranes of all living organisms. They play important roles in the electrogenic transport of ions and metabolites. Because of biophysical properties such as high selectivity for the permeant ion, high turnover rate, and modulation by physico-chemical parameters (e.g., membrane potential, calcium concentration), they are involved in several physiological processes in plant cells (e.g., maintenance of the turgor pressure, stomatal movements, and nutrient absorption by the roots). As plants cannot move, plant metabolism must be flexible and dynamic, to cope with environmental changes, to compete with other living species and to prevent pathogen invasion. An example of this flexibility and dynamic behavior is represented by their handling of the so-called reactive oxygen species, inevitable by-products of aerobic metabolism. Plants cope with these species on one side avoiding their toxic effects, on the other utilizing them as signalling molecules and as a means of defence against pathogens. In this review, we present the state-of-the-art of the modulation of plant ion channels by oxidizing and reducing agents.

Mechanism of anionic conduction across ClC

ClC chloride channels are voltage-gated transmembrane proteins that have been associated with a wide range of regulatory roles in vertebrates. To accomplish their function, they allow small inorganic anions to efficiently pass through, while blocking the passage of all other particles. Understanding the conduction mechanism of ClC has been the subject of many experimental investigations, but until now, the detailed dynamic mechanism was not known despite the availability of crystallographic structures. We investigate Cl(-) conduction by means of an all-atom molecular dynamics simulation of the ClC channel in a membrane environment. Based on our simulation results, we propose a king-of-the-hill mechanism for permeation, in which a lone ion bound to the center of the ClC pore is pushed out by a second ion that enters the pore and takes its place. Although the energy required to extract the single central ion from the pore is enormous, by resorting to this two-ion process, the largest free energy barrier for conduction is reduced to 4 kcal/mol. At the narrowest part of the pore, residues Tyr-445 and Ser-107 stabilize the central ion. There, the bound ion blocks the pore, disrupting the formation of a continuous water file that could leak protons, possibly preventing the passage of uncharged solutes.

Kinetic models of ion transport through a nanopore

Kinetic equations for the stationary state distribution function of ions moving through narrow pores are solved for a number of 1D models of single ion transport. Ions move through pores of length L, under the action of a constant external field and of a concentration gradient. The interaction of single ions with the confining pore surface and with water molecules inside the pore are modeled by a Fokker-Planck term in the kinetic equation, or by uncorrelated collisions with thermalizing centers distributed along the pore. The temporary binding of ions to polar residues lining the pore is modeled by stopping traps or energy barriers. Analytic expressions for the stationary ion current through the pore are derived for several versions of the model, as functions of key physical parameters. In all cases, saturation of the current at high fields is predicted. Such simple models, for which results are analytic, may prove useful in the study of the current/voltage relations of ion channels through membranes.

Pore topology of the hyperpolarization-activated cyclic nucleotide-gated channel from sea urchin sperm

The current flow through hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, referred to as I(h), plays a major role in several fundamental biological processes. The sequence of the presumed pore region of HCN channels is reminiscent of that of most known K(+)-selective channels. In the present work, the pore topology of an HCN channel from sea urchin sperm, called SpHCN, was investigated by means of the substituted-cysteine accessibility method (SCAM). The I(h) current in the wild-type (w.t.) SpHCN channel was irreversibly blocked by intracellular Cd(2+). This blockage was not observed in mutant C428S. Extracellular Cd(2+) did not cause any inhibition of the I(h) current in the w.t. SpHCN channel, but blocked the current in mutant channels K433C and F434C. Large extracellular anions blocked the current both in the w.t. and K433Q mutant channel. These results suggest that 1) cysteine in position 428 faces the intracellular medium; 2) lysine and phenylalanine in position 433 and 434, respectively, face the extracellular side of the membrane; and 3) lysine 433 does not mediate the anion blockade. Additionally, our study confirms that the K(+) channel signature sequence GYG also forms the inner pore in HCN channels.

The chloroplast protein import channel Toc75: pore properties and interaction with transit peptides

The channel properties of Toc75 (the protein import pore of the outer chloroplastic membrane) were further characterized by electrophysiological measurements in planar lipid bilayers. After improvement of the Toc75 reconstitution procedure the voltage dependence of the channel open probability resembled those observed for other beta-barrel pores. Studies concerning the pore size of the reconstituted Toc75 indicate the presence of a narrow restriction zone corresponding to the selectivity filter and a wider pore vestibule with diameters of approximately 14 A and 26 A, respectively. Interactions between Toc75 and different peptides (a genuine chloroplastic transit peptide, a synthetic peptide resembling a transit peptide, and a mitochondrial presequence) show that Toc75 itself is able to differentiate between these peptides and the recognition is based on both conformational and electrostatic interactions.

K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations

Microscopic molecular dynamics free energy perturbation calculations of the K(+)/Na(+) selectivity in the KcsA potassium channel, based on its experimental three-dimensional structure, are reported. The relative binding free energies for K(+) and Na(+) in the most relevant ion occupancy states of the four-site selectivity filter are calculated. The previously proposed mechanism for ion permeation through the KcsA channel is predicted, in agreement with available experimental data, to have a significant selectivity for K(+) over Na(+). The calculations also show that the individual 'binding site' selectivities are generally not additive and the doubly loaded states of the filter thus display cooperative effects. The only site that is not K(+) selective is that which is located at the entrance to the internal water cavity, suggesting the possibility that internal Na(+) could block outward currents.

The interaction of Na(+) and K(+) in the pore of cyclic nucleotide-gated channels

The permeability ratio between K(+) and Na(+) ions in cyclic nucleotide-gated channels is close to 1, and the single channel conductance has almost the same value in the presence of K(+) or Na(+). Therefore, K(+) and Na(+) ions are thought to permeate with identical properties. In the alpha-subunit from bovine rods there is a loop of three prolines at positions 365 to 367. When proline 365 is mutated to a threonine, a cysteine, or an alanine, mutant channels exhibit a complex interaction between K(+) and Na(+) ions. Indeed K(+), Rb(+) and Cs(+) ions do not carry any significant macroscopic current through mutant channels P365T, P365C and P365A and block the current carried by Na(+) ions. Moreover in mutant P365T the presence of K(+) in the intracellular (or extracellular) medium caused the appearance of a large transient inward (or outward) current carried by Na(+) when the voltage command was quickly stepped to large negative (or positive) membrane voltages. This transient current is caused by a transient potentiation, i.e., an increase of the open probability. The permeation of organic cations through these mutant channels is almost identical to that through the wild type (w.t.) channel. Also in the w.t. channel a similar but smaller transient current is observed, associated to a slowing down of the channel gating evident when intracellular Na(+) is replaced with K(+). As a consequence, a rather simple mechanism can explain the complex behavior here described: when a K(+) ion is occupying the pore there is a profound blockage of the channel and a potentiation of gating immediately after the K(+) ion is driven out. Potentiation occurs because K(+) ions slow down the rate constant K(off) controlling channel closure. These results indicate that K(+) and Na(+) ions do not permeate through CNG channels in the same way and that K(+) ions influence the channel gating.

Electronic excitation as a mechanism of the ion selectivity filter

A correlation between the character of pharmacological activity and the energies of electronic transitions in some biologically active molecules, affecting the nervous system, has been found. In order to explain the correlation, a new principle of the membrane ion selectivity filter has been suggested. The principle is based on the recombination process of a metal cation, passing through the filter, with an electron, when the energy quantum (equal to the metal ionization energy) is emitted. The amino acid residue group, performing the function of the channel filter, absorbs this quantum, transits itself into an electronically excited state, changes its conformation and lets, as a result, the cation pass. The process is possible only in that case when the amino acid residue group has a transition of the same energy, therefore not all of metals can pass through the filter. From the viewpoint of this conception, an active molecule acts because of its transition into an electronically excited state of the same energy and interfering, thereby, with the natural processes.

Molecular determinants of anion selectivity in the cystic fibrosis transmembrane conductance regulator chloride channel pore

Ionic selectivity in many cation channels is achieved over a short region of the pore known as the selectivity filter, the molecular determinants of which have been identified in Ca(2+), Na(+), and K(+) channels. However, a filter controlling selectivity among different anions has not previously been identified in any Cl(-) channel. In fact, because Cl(-) channels are only weakly selective among small anions, and because their selectivity has proved so resistant to site-directed mutagenesis, the very existence of a discrete anion selectivity filter has been called into question. Here we show that mutation of a putative pore-lining phenylalanine residue, F337, in the sixth membrane-spanning region of the cystic fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel, dramatically alters the relative permeabilities of different anions in the channel. Specifically, mutations that reduce the size of the amino acid side chain present at this position virtually abolish the relationship between anion permeability and hydration energy, a relationship that characterizes the anion selectivity not only of wild-type CFTR, but of most classes of Cl(-) channels. These results suggest that the pore of CFTR may indeed contain a specialized region, analogous to the selectivity filter of cation channels, at which discrimination between different permeant anions takes place. Because F337 is adjacent to another amino acid residue, T338, which also affects anion selectivity in CFTR, we suggest that selectivity is predominantly determined over a physically discrete region of the pore located near these important residues.

General anesthetic binding to gramicidin A: the structural requirements

There is a distinct possibility that general anesthetics exert their action on the postsynaptic receptor channels. The structural requirements for anesthetic binding in transmembrane channels, however, are largely unknown. High-resolution (1)H nuclear magnetic resonance and direct photoaffinity labeling were used in this study to characterize the volatile anesthetic binding sites in gramicidin A (gA) incorporated into sodium dodecyl sulfate (SDS) micelles and into dimyristoylphosphatidylcholine (DMPC) bilayers, respectively. To confirm that the structural arrangement of the peptide side chains can affect anesthetic binding, gA in nonchannel forms in methanol was also analyzed. The addition of volatile anesthetic halothane to gA in SDS with a channel conformation caused a concentration-dependent change in resonant frequencies of the indole amide protons of W9, W11, W13, and W15, with the most profound changes in W9. These frequency changes were observed only for gA carefully prepared to ensure a channel conformation and were absent for gA in methanol. For gA in DMPC bilayers, direct [(14)C]halothane photolabeling and microsequencing demonstrated dominant labeling of W9, less labeling of W11 and W13, and no significant labeling of W15. In methanol, gA showed much less labeling of any residues. Inspection of the 3-D structure of gA suggests that the spatial arrangements of the tryptophan residues in the channel form of gA, combined with the amphiphilic regions of lipid, create a favorable anesthetic binding motif.

Exploring the open pore of the potassium channel from Streptomyces lividans

The tetrameric potassium channel from Streptomyces lividans (KcsA) embedded in planar bilayers exhibits the following electrophysiological characteristics: (i) K+ ions can cross the pore in a highly hydrated state (nH2O > or = 6), (ii) the selectivity for K+ exceeds that for Na+ ions by 11 times, and both Ca2+ and Mg2+ are permeant, (iii) the internal side is blocked by Ba2+ ions in a voltage-dependent manner, (iv) intrinsic rectification is due to gating, depending on the direction of the electric field, (v) the internal side is pH-sensitive, and (vi) the open pore has a diameter of approximately 5.8 A. In conclusion, our results show that ion conduction and selectivity of KcsA cannot easily be reconciled with the properties deduced from the rigid crystal structure [Doyle et al., Science 280 (1998) 69-77], which must be concluded to have the pore trapped in its closed state.

Molecular determinants of electrical rectification of single channel conductance in gap junctions formed by connexins 26 and 32

The fully open state of heterotypic gap junction channels formed by pairing cells expressing connexin 32 (Cx32) with those expressing connexin 26 (Cx26) rectifies in a way that cannot be predicted from the current-voltage (I-V) relation of either homotypic channel. Using a molecular genetic analysis, we demonstrate that charged amino acids positioned in the amino terminus (M1 and D2) and first extracellular loop (E42) are major determinants of the current-voltage relation of the fully open state of homotypic and heterotypic channels formed by Cx26 and Cx32. The observed I-V relations of wild-type and mutant channels were closely approximated by those obtained with the electrodiffusive model of Chen and Eisenberg (Chen, D., and R. Eisenberg. 1993. Biophys. J. 64:1405-1421), which solves the Poisson-Nernst-Plank equations in one dimension using charge distribution models inferred from the molecular analyses. The rectification of the Cx32/Cx26 heterotypic channel results from the asymmetry in the number and position of charged residues. The model required the incorporation of a partial charge located near the channel surface to approximate the linear I-V relation observed for the Cx32*Cx26E1 homotypic channel. The best candidate amino acid providing this partial charge is the conserved tryptophan residue (W3). Incorporation of the partial charge of residue W3 and the negative charge of the Cx32E41 residue into the charge profile used in the Poisson-Nernst-Plank model of homotypic Cx32 and heterotypic Cx26/Cx32 channels resulted in I-V relations that closely resembled the observed I-V relations of these channels. We further demonstrate that some channel substates rectify. We suggest that the conformational changes associated with transjunctional voltage (V(j))-dependent gating to these substates involves a narrowing of the cytoplasmic entry of the channel that increases the electrostatic effect of charges in the amino terminus. The rectification that is observed in the Cx32/Cx26 heterotypic channel is similar although less steep than that reported for some rectifying electrical synapses. We propose that a similar electrostatic mechanism, which results in rectification through the open and substates of heterotypic channels, is sufficient to explain the properties of steeply rectifying electrical synapses.