Ion transport across transmembrane pores

To study the pore-mediated transport of ionic species across a lipid membrane, a series of molecular dynamics simulations have been performed of a dipalmitoyl-phosphatidyl-choline bilayer containing a preformed water pore in the presence of sodium and chloride ions. It is found that the stability of the transient water pores is greatly reduced in the presence of the ions. Specifically, the binding of sodium cations at the lipid/water interface increases the pore line tension, resulting in a destabilization of the pore. However, the application of mechanical stress opposes this effect. The flux of ions through these mechanically stabilized pores has been analyzed. Simulations indicate that the transport of the ions through the pores depends strongly on the size of the water channel. In the presence of small pores (radius <1.5 nm) permeation is slow, with both sodium and chloride permeating at similar rates. In the case in which the pores are larger (radius >1.5 nm), a crossover is observed to a regime where the anion flux is greatly enhanced. Based on these observations, a mechanism for the basal membrane permeability of ions is discussed.     

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

Voltage-sensing arginines in a potassium channel permeate and occlude cation-selective pores

Voltage-gated ion channels sense voltage by shuttling arginine residues located in the S4 segment across the membrane electric field. The molecular pathway for this arginine permeation is not understood, nor is the filtering mechanism that permits passage of charged arginines but excludes solution ions. We find that substituting the first S4 arginine with smaller amino acids opens a high-conductance pathway for solution cations in the Shaker K(+) channel at rest. The cationic current does not flow through the central K(+) pore and is influenced by mutation of a conserved residue in S2, suggesting that it flows through a protein pathway within the voltage-sensing domain. The current can be carried by guanidinium ions, suggesting that this is the pathway for transmembrane arginine permeation. We propose that when S4 moves it ratchets between conformations in which one arginine after another occupies and occludes to ions the narrowest part of this pathway.     

Simplified bacterial "pore" channel provides insight into the assembly, stability, and structure of sodium channels

Eukaryotic sodium channels are important membrane proteins involved in ion permeation, homeostasis, and electrical signaling. They are long, multidomain proteins that do not express well in heterologous systems, and hence, structure/function and biochemical studies on purified sodium channel proteins have been limited. Bacteria produce smaller, homologous tetrameric single domain channels specific for the conductance of sodium ions. They consist of N-terminal voltage sensor and C-terminal pore subdomains. We designed a functional pore-only channel consisting of the final two transmembrane helices, the intervening P-region, and the C-terminal extramembranous region of the sodium channel from the marine bacterium Silicibacter pomeroyi. This sodium "pore" channel forms a tetrameric, folded structure that is capable of supporting sodium flux in phospholipid vesicles. The pore-only channel is more thermally stable than its full-length counterpart, suggesting that the voltage sensor subdomain may destabilize the full-length channel. The pore subdomains can assemble, fold, and function independently from the voltage sensor and exhibit similar ligand-blocking characteristics as the intact channel. The availability of this simple pore-only construct should enable high-level expression for the testing of potential new ligands and enhance our understanding of the structural features that govern sodium selectivity and permeability.     

Membrane permeability changes at early stages of influenza hemagglutinin-mediated fusion

While biological membrane fusion is classically defined as the leak-free merger of membranes and contents, leakage is a finding in both experimental and theoretical studies. The fusion stages, if any, that allow membrane permeation are uncharted. In this study we monitored membrane ionic permeability at early stages of fusion mediated by the fusogenic protein influenza hemagglutinin (HA). HAb2 cells, expressing HA on their plasma membrane, fused with human red blood cells, cultured liver cells PLC/PRF/5, or planar phospholipid bilayer membranes. With a probability that depended upon the target membrane, an increase of the electrical conductance of the fusing membranes (leakage) by up to several nS was generally detected. This leakage was recorded at the initial stages of fusion, when fusion pores formed. This leakage usually accompanied the "flickering" stage of the early fusion pore development. As the pore widened, the leakage reduced; concomitantly, the lipid exchange between the fusing membranes accelerated. We conclude that during fusion pore formation, HA locally and temporarily increases the permeability of fusing membranes. Subsequent rearrangement in the fusion complex leads to the resealing of the leaky membranes and enlargement of the pore.     

Positive charges at the intracellular mouth of the pore regulate anion conduction in the CFTR chloride channel

Many different ion channel pores are thought to have charged amino acid residues clustered around their entrances. The so-called surface charges contributed by these residues can play important roles in attracting oppositely charged ions from the bulk solution on one side of the membrane, increasing effective local counterion concentration and favoring rapid ion movement through the channel. Here we use site-directed mutagenesis to identify arginine residues contributing important surface charges in the intracellular mouth of the cystic fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel pore. While wild-type CFTR was associated with a linear current-voltage relationship with symmetrical solutions, strong outward rectification was observed after mutagenesis of two arginine residues (R303 and R352) located near the intracellular ends of the fifth and sixth transmembrane regions. Current rectification was dependent on the charge present at these positions, consistent with an electrostatic effect. Furthermore, mutagenesis-induced rectification was more pronounced at lower Cl(-) concentrations, suggesting that these mutants had a reduced ability to concentrate Cl(-) ions near the inner pore mouth. R303 and R352 mutants exhibited reduced single channel conductance, especially at negative membrane potentials, that was dependent on the charge of the amino acid residue present at these positions. However, the very low conductance of both R303E and R352E-CFTR could be greatly increased by elevating intracellular Cl(-) concentration. Modification of an introduced cysteine residue at position 303 by charged methanethiosulfonate reagents reproduced charge-dependent effects on current rectification. Mutagenesis of arginine residues in the second and tenth transmembrane regions also altered channel permeation properties, however these effects were not consistent with changes in channel surface charges. These results suggest that positively charged arginine residues act to concentrate Cl(-) ions at the inner mouth of the CFTR pore, and that this contributes to maximization of the rate of Cl(-) ion permeation through the pore.     

Spontaneous transmembrane pore formation by short-chain synthetic peptide

Amphiphilic β-peptides, which are synthetically designed short-chain helical foldamers of β-amino acids, are established potent biomimetic alternatives of natural antimicrobial peptides. An intriguing question is how the distinct molecular architecture of these short-chain and rigid synthetic peptides translates to its potent membrane-disruption ability. Here, we address this question via a combination of all-atom and coarse-grained molecular dynamics simulations of the interaction of mixed phospholipid bilayer with an antimicrobial 10-residue globally amphiphilic helical β-peptide at a wide range of concentrations. The simulation demonstrates that multiple copies of this synthetic peptide, initially placed in aqueous solution, readily self-assemble and adsorb at membrane interface. Subsequently, beyond a threshold peptide/lipid ratio, the surface-adsorbed oligomeric aggregate moves inside the membrane and spontaneously forms stable water-filled transmembrane pores via a cooperative mechanism. The defects induced by these pores lead to the dislocation of interfacial lipid headgroups, membrane thinning, and substantial water leakage inside the hydrophobic core of the membrane. A molecular analysis reveals that despite having a short architecture, these synthetic peptides, once inside the membrane, would stretch themselves toward the distal leaflet in favor of potential contact with polar headgroups and interfacial water layer. The pore formed in coarse-grained simulation was found to be resilient upon structural refinement. Interestingly, the pore-inducing ability was found to be elusive in a non-globally amphiphilic sequence isomer of the same β-peptide, indicating strong sequence dependence. Taken together, this work puts forward key perspectives of membrane activity of minimally designed synthetic biomimetic oligomers relative to the natural antimicrobial peptides.     

Bilayer lipid composition modulates the activity of dermaseptins, polycationic antimicrobial peptides

The primary targets of defense peptides are plasma membranes, and the induced irreversible depolarization is sufficient to exert antimicrobial activity although secondary modes of action might be at work. Channels or pores underlying membrane permeabilization are usually quite large with single-channel conductances two orders of magnitude higher than those exhibited by physiological channels involved, e.g., in excitability. Accordingly, the ion specificity and selectivity are quite low. Whereas, e.g., peptaibols favor cation transport, polycationic or basic peptides tend to form anion-specific pores. With dermaseptin B2, a 33 residue long and mostly α-helical peptide isolated from the skin of the South American frog Phyllomedusa bicolor, we found that the ion specificity of its pores induced in bilayers is modulated by phospholipid-charged headgroups. This suggests mixed lipid–peptide pore lining instead of the more classical barrel–stave model. Macroscopic conductance is nearly voltage independent, and concentration dependence suggests that the pores are mainly formed by dermaseptin tetramers. The two most probable single-channel events are well resolved at 200 and 500 pS (in 150 mM NaCl) with occasional other equally spaced higher or lower levels. In contrast to previous molecular dynamics previsions, this study demonstrates that dermaseptins are able to form pores, although a related analog (B6) failed to induce any significant conductance. Finally, the model of the pore we present accounts for phospholipid headgroups intercalated between peptide helices lining the pore and for one of the most probable single-channel conductance.     

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

Molecular dynamics investigation of the ω-current in the Kv1.2 voltage sensor domains

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S.     

Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations.

To test the hypothesis that water pores in a lipid membrane mediate the proton transport, molecular dynamic simulations of a phospholipid membrane, in which the formation of a water pore is induced, are reported. The probability density of such a pore in the membrane was obtained from the free energy of formation of the pore, which was computed from the average force needed to constrain the pore in the membrane. It was found that the free energy of a single file of water molecules spanning the bilayer is 108(+/-10) kJ/mol. From unconstrained molecular dynamic simulations it was further deduced that the nature of the pore is very transient, with a mean lifetime of a few picoseconds. The orientations of water molecules within the pore were also studied, and the spontaneous translocation of a turning defect was observed. The combined data allowed a permeability coefficient for proton permeation across the membrane to be computed, assuming that a suitable orientation of the water molecules in the pore allows protons to permeate the membrane relatively fast by means of a wirelike conductance mechanism. The computed value fits the experimental data only if it is assumed that the entry of the proton into the pore is not rate limiting.     

Gating pore currents,a new pathological mechanism underlying cardiac arrhythmias associated with dilated cardiomyopathy

Voltage-gated ion channels (VGIC) are transmembrane proteins responsible for the generation of electrical signals in excitable cells. VGIC were first described in 1952 by Hodgkin and Huxley,¹ and have since been associated with various physiological functions such as propagating nerve impulses, locomotion, and cardiac excitability. VGIC include channels specialized in the selective passage of K⁺, Ca²⁺ Na⁺, or H⁺. They are composed of 2 main structures: the pore domain (PD) and the voltage sensor domain (VSD). The PD ensures the physiological flow of ions and is typically composed of 8 transmembrane segments (TM). The VSD detects voltage variations and is composed of 4 TM (S1-S4). Given their crucial physiological role, VGIC dysfunctions are associated with diverse pathologies known as ion channelopathies. These dysfunctions usually affect the membrane expression of ion channels or voltage-dependent conformational changes of the pore. However, an increasing number of ion channelopathies, including periodic paralysis, dilated cardiomyopathy (DCM) associated with cardiac arrhythmias, and peripheral nerve hyperexcitability (PNH), have been linked to the appearance of a new pathological mechanism involving the creation of an alternative permeation pathway through the normally non-conductive VSD of VGIC. This permeation pathway is called the gating pore or omega pore.     

Plants do it differently. A new basis for potassium/sodium selectivity in the pore of an ion channel

Understanding of the molecular architecture necessary for selective K(+) permeation through the pore of ion channels is based primarily on analysis of the crystal structure of the bacterial K(+) channel KcsA, and structure:function studies of cloned animal K(+) channels. Little is known about the conduction properties of a large family of plant proteins with structural similarities to cloned animal cyclic nucleotide-gated channels (CNGCs). Animal CNGCs are nonselective cation channels that do not discriminate between Na(+) and K(+) permeation. These channels all have the same triplet of amino acids in the channel pore ion selectivity filter, and this sequence is different from that of the selectivity filter found in K(+)-selective channels. Plant CNGCs have unique pore selectivity filters; unlike those found in any other family of channels. At present, the significance of the unique pore selectivity filters of plant CNGCs, with regard to discrimination between Na(+) and K(+) permeation is unresolved. Here, we present an electrophysiological analysis of several members of this protein family; identifying the first cloned plant channel (AtCNGC1) that conducts Na(+). Another member of this ion channel family (AtCNGC2) is shown to have a selectivity filter that provides a heretofore unknown molecular basis for discrimination between K(+) and Na(+) permeation. Specific amino acids within the AtCNGC2 pore selectivity filter (Asn-416, Asp-417) are demonstrated to facilitate K(+) over Na(+) conductance. The selectivity filter of AtCNGC2 represents an alternative mechanism to the well-known GYG amino acid triplet of K(+) channels that has been identified as the critical basis for K(+) over Na(+) permeation through the pore of ion channels.     

Molecular Model of a Cell Plasma Membrane With an Asymmetric Multicomponent Composition: Water Permeation and Ion Effects

We present molecular dynamics simulations of a multicomponent, asymmetric bilayer in mixed aqueous solutions of sodium and potassium chloride. Because of the geometry of the system, there are two aqueous solution regions in our simulations: one mimics the intracellular region, and one mimics the extracellular region. Ion-specific effects are evident at the membrane/aqueous solution interface. Namely, at equal concentrations of sodium and potassium, sodium ions are more strongly adsorbed to carbonyl groups of the lipid headgroups. A significant concentration excess of potassium is needed for this ion to overwhelm the sodium abundance at the membrane. Ion-membrane interactions also lead to concentration-dependent and cation-specific behavior of the electrostatic potential in the intracellular region because of the negative charge on the inner leaflet. In addition, water permeation across the membrane was observed on a timescale of ∼100 ns. This study represents a step toward the modeling of realistic biological membranes at physiological conditions in intracellular and extracellular environments.     

Design of Peptide-Membrane Interactions to Modulate Single-File Water Transport through Modified Gramicidin Channels

Water permeability through single-file channels is affected by intrinsic factors such as their size and polarity and by external determinants like their lipid environment in the membrane. Previous computational studies revealed that the obstruction of the channel by lipid headgroups can be long-lived, in the range of nanoseconds, and that pore-length-matching membrane mimetics could speed up water permeability. To test the hypothesis of lipid-channel interactions modulating channel permeability, we designed different gramicidin A derivatives with attached acyl chains. By combining extensive molecular-dynamics simulations and single-channel water permeation measurements, we show that by tuning lipid-channel interactions, these modifications reduce the presence of lipid headgroups in the pore, which leads to a clear and selective increase in their water permeability.     

Calcium influx into phospholipid vesicles caused by dynorphin neuropeptides

Dynorphins, endogeneous opioid peptides, function as ligands to the opioid kappa receptors but also induce non-opioid excitotoxic effects. Dynorphin A can increase the intra-neuronal calcium concentration through a non-opioid and non-NMDA mechanism. In this investigation, we show that big dynorphin, dynorphin A and to some extent dynorphin A (1-13), but not dynorphin B, allow calcium to enter into large unilamellar phospholipid vesicles with partly negative headgroups. The effects parallel the previously studied potency of dynorphins to translocate through biological membranes and to cause calcein leakage from large unilamellar phospholipid vesicles. There is no calcium ion influx into vesicles with zwitterionic headgroups. We have also investigated if the dynorphins can translocate through the vesicle membranes and estimated the relative strength of interaction of the peptides with the vesicles by fluorescence resonance energy transfer. The results show that dynorphins do not translocate in this membrane model system. There is a strong electrostatic contribution to the interaction of the peptides with the membrane model system.     

Calcium influx into phospholipid vesicles caused by dynorphin neuropeptides

Dynorphins, endogeneous opioid peptides, function as ligands to the opioid kappa receptors but also induce non-opioid excitotoxic effects. Dynorphin A can increase the intra-neuronal calcium concentration through a non-opioid and non-NMDA mechanism. In this investigation, we show that big dynorphin, dynorphin A and to some extent dynorphin A (1-13), but not dynorphin B, allow calcium to enter into large unilamellar phospholipid vesicles with partly negative headgroups. The effects parallel the previously studied potency of dynorphins to translocate through biological membranes and to cause calcein leakage from large unilamellar phospholipid vesicles. There is no calcium ion influx into vesicles with zwitterionic headgroups. We have also investigated if the dynorphins can translocate through the vesicle membranes and estimated the relative strength of interaction of the peptides with the vesicles by fluorescence resonance energy transfer. The results show that dynorphins do not translocate in this membrane model system. There is a strong electrostatic contribution to the interaction of the peptides with the membrane model system.     

Water and ion permeability of a claudin model: A computational study

At present, the three‐dimensional structure of the multimeric paracellular claudin pore is unknown. Using extant biophysical data concerning the size of the pore and permeation of water and cations through it, two three‐dimensional models of the pore are constructed in silico. Molecular Dynamics (MD) calculations are then performed to compute water and sodium ion permeation fluxes under the influence of applied hydrostatic pressure. Comparison to experiment is made, under the assumption that the hydrostatic pressure applied in the simulations is equivalent to osmotic pressure induced in experimental measurements of water/ion permeability. One model, in which pore‐lining charged is distributed evenly over a selectivity filter section 10–16 Å in length, is found to be generally consistent with experimental data concerning the dependence of water and ion permeation on channel pore diameter, pore length, and the sign and magnitude of pore lining charge. The molecular coupling mechanism between water and ion flow under conditions where hydrostatic pressure is applied is computationally elucidated. Proteins 2016; 84:305–315. © 2016 Wiley Periodicals, Inc.