Structure-based prediction of the conductance properties of ion channels

Fusion pores with low conductance are cation selective

Many neurotransmitters are organic ions that carry a net charge, and their release from secretory vesicles is therefore an electrodiffusion process. The selectivity of early exocytotic fusion pores is investigated by combining electrodiffusion theory, measurements of amperometric foot signals from chromaffin cells with anion substitution, and molecular dynamics simulation. The results reveal that very narrow fusion pores are cation selective, but more dilated fusion pores become anion permeable. The transition occurs around a fusion pore conductance of ~300 pS. The cation selectivity of a narrow fusion pore accelerates the release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

For transmission, a fusion pore forms when vesicle and target membranes are brought together by SNARE proteins. Delacruz et al. demonstrate that selectivity of the pore accelerates release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

Molecular mechanism of fusion pore formation driven by the neuronal SNARE complex

Release of neurotransmitters from synaptic vesicles begins with a narrow fusion pore, the structure of which remains unresolved. To obtain a structural model of the fusion pore, we performed coarse-grained molecular dynamics simulations of fusion between a nanodisc and a planar bilayer bridged by four partially unzipped SNARE complexes. The simulations revealed that zipping of SNARE complexes pulls the polar C-terminal residues of the synaptobrevin 2 and syntaxin 1A transmembrane domains to form a hydrophilic core between the two distal leaflets, inducing fusion pore formation. The estimated conductances of these fusion pores are in good agreement with experimental values. Two SNARE protein mutants inhibiting fusion experimentally produced no fusion pore formation. In simulations in which the nanodisc was replaced by a 40-nm vesicle, an extended hemifusion diaphragm formed but a fusion pore did not, indicating that restricted SNARE mobility is required for rapid fusion pore formation. Accordingly, rapid fusion pore formation also occurred in the 40-nm vesicle system when SNARE mobility was restricted by external forces. Removal of the restriction is required for fusion pore expansion.

The fusion pore, 60 years after the first cartoon

Neurotransmitter release occurs in the form of quantal events by fusion of secretory vesicles with the plasma membrane, and begins with the formation of a fusion pore that has a conductance similar to that of a large ion channel or gap junction. In this review, we propose mechanisms of fusion pore formation and discuss their implications for fusion pore structure and function. Accumulating evidence indicates a direct role of soluble N-ethylmaleimide-sensitive-factor attachment receptor proteins in the opening of fusion pores. Fusion pores are likely neither protein channels nor purely lipid, but are of proteolipidic composition. Future perspectives to gain better insight into the molecular structure of fusion pores are discussed.

Protein-lipid interactions and non-lamellar lipidic structures in membrane pore formation and membrane fusion

Pore-forming proteins and peptides act on their targeted lipid bilayer membranes to increase permeability. This approach to the modulation of biological function is relevant to a great number of living processes, including; infection, parasitism, immunity, apoptosis, development and neurodegeneration. While some pore-forming proteins/peptides assemble into rings of subunits to generate discrete, well-defined pore-forming structures, an increasing number is recognised to form pores via mechanisms which co-opt membrane lipids themselves. Among these, membrane attack complex-perforin/cholesterol-dependent cytolysin (MACPF/CDC) family proteins, Bax/colicin family proteins and actinoporins are especially prominent and among the mechanisms believed to apply are the formation of non-lamellar (semi-toroidal or toroidal) lipidic structures. In this review I focus on the ways in which lipids contribute to pore formation and contrast this with the ways in which lipids are co-opted also in membrane fusion and fission events. A variety of mechanisms for pore formation that involve lipids exists, but they consistently result in stable hybrid proteolipidic structures. These structures are stabilised by mechanisms in which pore-forming proteins modify the innate capacity of lipid membranes to respond to their environment, changing shape and/or phase and binding individual lipid molecules directly. In contrast, and despite the diversity in fusion protein types, mechanisms for membrane fusion are rather similar to each other, mapping out a pathway from pairs of separated compartments to fully confluent fused membranes. Fusion proteins generate metastable structures along the way which, like long-lived proteolipidic pore-forming complexes, rely on the basic physical properties of lipid bilayers. Membrane fission involves similar intermediates, in the reverse order. I conclude by considering the possibility that at least some pore-forming and fusion proteins are evolutionarily related homologues. This article is part of a Special Issue entitled: Pore-Forming Toxins edited by Mauro Dalla Serra and Franco Gambale.

Effects of MACPF/CDC proteins on lipid membranes

Recent work on the MACPF/CDC superfamily of pore-forming proteins has focused on the structural analysis of monomers and pore-forming oligomeric complexes. We set the family of proteins in context and highlight aspects of their function which the direct and exclusive equation of oligomers with pores fails to explain. Starting with a description of the distribution of MACPF/CDC proteins across the domains of life, we proceed to show how their evolutionary relationships can be understood on the basis of their structural homology and re-evaluate models for pore formation by perforin, in particular. We furthermore highlight data showing the role of incomplete oligomeric rings (arcs) in pore formation and how this can explain small pores generated by oligomers of proteins belonging to the family. We set this in the context of cell biological and biophysical data on the proteins' function and discuss how this helps in the development of an understanding of how they act in processes such as apicomplexan parasites gliding through cells and exiting from cells.

Ion selectivity of alpha-hemolysin with beta-cyclodextrin adapter. II. Multi-ion effects studied with grand canonical Monte Carlo/Brownian dynamics simulations

In a previous study of ion selectivity of alpha-hemolysin (alphaHL) in complex with beta-cyclodextrin (betaCD) adapter, we calculated the potential of mean force (PMF) and characterized the self-diffusion coefficients of isolated K(+) and Cl(-) ions using molecular dynamics simulations (Y. Luo et al., "Ion Selectivity of alpha-Hemolysin with beta-Cyclodextrin Adapter: I. Single Ion Potential of Mean Force and Diffusion Coefficient"). In the present effort, these results pertaining to single isolated ions in the wide aqueous pore are extended to take into account multi-ion effects. The grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) algorithm is used to simulate ion currents through the wild-type alphaHL ion channel, as well as two engineered alphaHL mutants, with and without the cyclic oligosaccaride betaCD lodged in the lumen of the pore. The GCMC/BD current-voltage curves agree well with experimental results and show that betaCD increases the anion selectivity of alphaHL. Comparisons between multi-ion PMFs from GCMC/BD simulations and single-ion PMFs demonstrate that multi-ion effects and pore shape are crucial for explaining this behavior. It is concluded that the narrow betaCD adapter increases the anion selectivity of alphaHL because it reduces the pore radius locally, which decreases the ionic screening and the dielectric shielding of the strong electrostatic field induced by a nearby ring of positively charged alphaHL side chains.

Model-based prediction of the alpha-hemolysin structure in the hexameric state

The alpha-hemolysin toxin self-assembles in lipid bilayers to form water-filled pores. In recent years, alpha-hemolysin has received great attention, mainly due to its possible usage as a sensing element. We measured the ion currents through single alpha-hemolysin channels and confirmed the presence of two different subpopulations of channels with conductance levels of 465 +/- 30 pS and 280 +/- 30 pS. Different oligomerization states could be responsible for these two conductances. In fact, a heptameric structure of the channel was revealed by x-ray crystallography, whereas atomic force microscopy revealed a hexameric structure. Due to the low resolution of atomic force microscopy the atomic details of the hexameric structure are still unknown, and are here predicted by computational methods. Several possible structures of the hexameric channel were defined, and were simulated by molecular dynamics. The conductances of these channel models were computed by a numerical method based on the Poisson-Nernst-Planck electrodiffusion theory, and the values were compared to experimental data. In this way, we identified a model of the alpha-hemolysin hexameric state with conductance characteristics consistent with the experimental data. Since the oligomerization state of the channel may affect its behavior as a molecular sensor, knowing the atomic structure of the hexameric state will be useful for biotechnological applications of alpha-hemolysin.

Refolding of Escherichia coli outer membrane protein F in detergent creates LPS-free trimers and asymmetric dimers

The Escherichia coli OmpF (outer-membrane protein F; matrix porin) is a homotrimeric beta-barrel and a member of the bacterial porin superfamily. It is the best characterized porin protein, but has resisted attempts to refold it efficiently in vitro. In the present paper, we report the discovery of detergent-based folding conditions, including dodecylglucoside, which can create pure samples of trimeric OmpF. Whereas outer membrane LPS (lipopolysaccharide) is clearly required for in vivo folding, the artificially refolded and LPS-free trimer has properties identical with those of the outer-membrane-derived form. Thus LPS is not required either for in vitro folding or for structural integrity. Dimeric forms of OmpF have been observed in vivo and are proposed to be folding intermediates. In vitro, dimers occur transiently in refolding of trimeric OmpF and, in the presence of dodecylmaltoside, pure dimer can be prepared. This form has less beta-structure by CD and shows lower thermal stability than the trimer. Study of these proteins at the single-molecule level is possible because each OmpF subunit forms a distinct ion channel. Whereas each trimer contains three channels of equal conductance, each dimer always contains two distinct channel sizes. This provides clear evidence that the two otherwise identical monomers adopt different structures in the dimer and indicates that the asymmetric interaction, characteristic of C3 symmetry, is formed at the dimer stage. This asymmetric dimer may be generally relevant to the folding of oligomeric proteins with odd numbers of subunits such as aspartate transcarbamoylase.

Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map

alpha-Hemolysin of Staphylococcus aureus is a self-assembling toxin that forms a water-filled transmembrane channel upon oligomerization in a lipid membrane. Apart from being one of the best-studied toxins of bacterial origin, alpha-hemolysin is the principal component in several biotechnological applications, including systems for controlled delivery of small solutes across lipid membranes, stochastic sensors for small solutes, and an alternative to conventional technology for DNA sequencing. Through large-scale molecular dynamics simulations, we studied the permeability of the alpha-hemolysin/lipid bilayer complex for water and ions. The studied system, composed of approximately 300,000 atoms, included one copy of the protein, a patch of a DPPC lipid bilayer, and a 1 M water solution of KCl. Monitoring the fluctuations of the pore structure revealed an asymmetric, on average, cross section of the alpha-hemolysin stem. Applying external electrostatic fields produced a transmembrane ionic current; repeating simulations at several voltage biases yielded a current/voltage curve of alpha-hemolysin and a set of electrostatic potential maps. The selectivity of alpha-hemolysin to Cl(-) was found to depend on the direction and the magnitude of the applied voltage bias. The results of our simulations are in excellent quantitative agreement with available experimental data. Analyzing trajectories of all water molecule, we computed the alpha-hemolysin's osmotic permeability for water as well as its electroosmotic effect, and characterized the permeability of its seven side channels. The side channels were found to connect seven His-144 residues surrounding the stem of the protein to the bulk solution; the protonation of these residues was observed to affect the ion conductance, suggesting the seven His-144 to comprise the pH sensor that gates conductance of the alpha-hemolysin channel.

Computing numerically the access resistance of a pore

The access resistance (AR) of a channel is an important component of the conductance of ion channels, particularly in wide and short channels, where it accounts for a substantial fraction of the total resistance to the movement of ions. The AR is usually calculated by using a classical and simple expression derived by Hall from electrostatics (J.E. Hall 1975 J. Gen. Phys. 66:531-532), though other expressions, both analytical and numerical, have been proposed. Here we report some numerical results for the AR of a channel obtained by solving the Poisson-Nernst-Planck equations at the entrance of a circular pore. Agreement is found between numerical calculations and analytical results from Hall's equation for uncharged pores in neutral membranes. However, for channels embedded in charged membranes, Hall's expression overestimates the AR, which is much lower and can even be neglected in some cases. The weak dependence of AR on the pore radius for charged membranes at low salt concentration can be exploited to separate the channel and the access contributions to the measured conductance.

Molecular dynamics investigation of an oriented cyclic peptide nanotube in DMPC bilayers

Nanotubes resulting from the self-assembly of cyclic peptides formed by eight alpha-amino acids and inserted into lipid bilayers have been shown to function as synthetic, integral transmembrane channels. A nanotube consisting of eight cyclo[(L-Trp-D-Leu)(3)-L-Gln-D-Leu] subunits, organized in an antiparallel, beta-sheetlike channel embedded in a hydrated dimyristoylphosphatidylcholine bilayer was investigated in an 8-ns molecular dynamics trajectory. This large-scale statistical simulation brings to light not only the atomic-level structural features of the synthetic channel, but also its dynamical properties. Overall, the nanotube conserves its hollow tubular structure. The calculation reproduces the tilt of the channel with respect to the normal of the bilayer, in reasonable agreement with experiment. The results show a dislocation of the nanotube indicative of a possible disassembly process that may influence the channel conduction. The dynamics of the water in the hollow tubular structure has been characterized, and the conductance of the channel has been estimated. Transport properties of the peptide nanotube are discussed in comparison with other transporters.

Shape anisotropy of lipid molecules and voids

Biological polymers, viz., proteins, membranes and micelles exhibit structural discontinuities in terms of spaces unfilled by the polymeric phase, termed voids. These voids exhibit dynamics and lead to interesting properties which are experimentally demonstrable. In the specific case of phospholipid membranes, numerical simulations on a two-dimensional model system showed that voids are induced primarily due to the shape anisotropy in binary mixtures of interacting disks. The results offer a minimal description required to explain the unusually large permeation seen in liposomes made up of specific lipid mixtures (Mathai & Sitaramam, 1994). The results are of wider interest, voids being ubiquitous in biopolymers.

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.

Analysis and evaluation of channel models: simulations of alamethicin

Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, 5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20 ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.

OmpA: a pore or not a pore? Simulation and modeling studies

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).

Membrane proteins: Aquaporins--channels without ions

Recently determined structures have shed new light on the way that aquaporins act as passive, but selective, pores for the transport of small molecules--such as water or glycerol--across membranes.

Structural investigations of pneumolysin/lipid complexes

Pneumolysin, a virulence factor from the human pathogen Streptococcus pneumoniae, is a water-soluble protein which forms ring-shaped oligomeric structures upon binding to cholesterol-containing lipid membranes. It induces vesicle aggregation, membrane pore formation and withdrawal of lipid material into non-bilayer proteolipid complexes. Solid-state magic angle spinning and wideline static NMR, together with freeze-fracture electron microscopy, are used to characterize the phase changes in fully hydrated cholesterol-containing lipid membranes induced by the addition of pneumolysin. A structural model for the proteolipid complexes is proposed where a 30-50-meric pneumolysin ring lines the inside of a lipid torus. Cholesterol is found to be essential to the fusogenic action of pneumolysin.

Structural basis of pore formation by cholesterol-binding toxins

In this paper we describe reconstructions by electron cryo-microscopy of two oligomeric states of the pore-forming toxin pneumolysin. The results are interpreted by the fitting of atomic models of separated domains to the 3-dimensional electron density maps, revealing two steps in the mechanism of pore formation by the family of cholesterol-binding toxins. We briefly describe the observation of the toxin pore in model membranes and contrast the apparent mechanism of pneumolysin with that of other pore-forming toxins.

Molecular dynamics of synthetic leucine-serine ion channels in a phospholipid membrane

Molecular dynamics calculations were carried out on models of two synthetic leucine-serine ion channels: a tetrameric bundle with sequence (LSLLLSL)(3)NH(2) and a hexameric bundle with sequence (LSSLLSL)(3)NH(2). Each protein bundle is inserted in a palmitoyloleoylphosphatidylcholine bilayer membrane and solvated by simple point charge water molecules inside the pore and at both mouths. Both systems appear to be stable in the absence of an electric field during the 4 ns of molecular dynamics simulation. The water motion in the narrow pore of the four-helix bundle is highly restricted and may provide suitable conditions for proton transfer via a water wire mechanism. In the wider hexameric pore, the water diffuses much more slowly than in bulk but is still mobile. This, along with the dimensions of the pore, supports the observation that this peptide is selective for monovalent cations. Reasonable agreement of predicted conductances with experimentally determined values lends support to the validity of the simulations.

Molecular basis of the charge selectivity of nicotinic acetylcholine receptor and related ligand-gated ion channels

Nicotinic acetylcholine receptors are homo- or heteropentameric proteins belonging to the superfamily of receptor channels including the glycine and GABA-A receptors. Affinity labelling and mutagenesis experiments indicated that the M2 transmembrane segment of each subunit lines the ion channel and is coiled into an alpha-helix. Comparison of the M2 sequence of the cation-selective alpha 7 nicotinic receptor to that of the anion-selective alpha 1 glycine receptor identified amino acids involved in charge selectivity. Mutations of the alpha 7 homo-oligomeric receptor within (or near) M2, namely E237A, V251T and a proline insertion P236' were shown to convert the ionic selectivity of alpha 7 from cationic to anionic. Systematic analysis of each of these three mutations supports the notion that the conversion of ionic selectivity results from a local structural reorganization of the 234-238 loop. The 234-238 coiled loop, previously shown to lie near the narrowest portion of the channel, is thus proposed to contribute directly to the charge selectivity filter. A possible functional analogy with the voltage-gated ion channels and related receptors is discussed.