Charge selectivity of the designed uncharged peptide ion channel Ac-(LSSLLSL)3-CONH2

Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel

Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K⁺ and Cl^- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K⁺, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl^- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.

Design of self-assembling transmembrane helical bundles to elucidate principles required for membrane protein folding and ion transport

Ion transporters and channels are able to identify and act on specific substrates among myriads of ions and molecules critical to cellular processes, such as homeostasis, cell signalling, nutrient influx and drug efflux. Recently, we designed Rocker, a minimalist model for Zn²⁺/H⁺ co-transport. The success of this effort suggests that de novo membrane protein design has now come of age so as to serve a key approach towards probing the determinants of membrane protein folding, assembly and function. Here, we review general principles that can be used to design membrane proteins, with particular reference to helical assemblies with transport function. We also provide new functional and NMR data that probe the dynamic mechanism of conduction through Rocker.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Gating of Connexin Channels by transjunctional-voltage: Conformations and models of open and closed states

Voltage is an important physiologic regulator of channels formed by the connexin gene family. Connexins are unique among ion channels in that both plasma membrane inserted hemichannels (undocked hemichannels) and intercellular channels (aggregates of which form gap junctions) have important physiological roles. The hemichannel is the fundamental unit of gap junction voltage-gating. Each hemichannel displays two distinct voltage-gating mechanisms that are primarily sensitive to a voltage gradient formed along the length of the channel pore (the transjunctional voltage) rather than sensitivity to the absolute membrane potential (V_m or V_i-o). These transjunctional voltage dependent processes have been termed V_j- or fast-gating and loop- or slow-gating. Understanding the mechanism of voltage-gating, defined as the sequence of voltage-driven transitions that connect open and closed states, first and foremost requires atomic resolution models of the end states. Although ion channels formed by connexins were among the first to be characterized structurally by electron microscopy and x-ray diffraction in the early 1980's, subsequent progress has been slow. Much of the current understanding of the structure-function relations of connexin channels is based on two crystal structures of Cx26 gap junction channels. Refinement of crystal structure by all-atom molecular dynamics and incorporation of charge changing protein modifications has resulted in an atomic model of the open state that arguably corresponds to the physiologic open state. Obtaining validated atomic models of voltage-dependent closed states is more challenging, as there are currently no methods to solve protein structure while a stable voltage gradient is applied across the length of an oriented channel. It is widely believed that the best approach to solve the atomic structure of a voltage-gated closed ion channel is to apply different but complementary experimental and computational methods and to use the resulting information to derive a consensus atomic structure that is then subjected to rigorous validation. In this paper, we summarize our efforts to obtain and validate atomic models of the open and voltage-driven closed states of undocked connexin hemichannels. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

Molecular dynamics simulations of the Cx26 hemichannel: evaluation of structural models with Brownian dynamics

The recently published crystal structure of the Cx26 gap junction channel provides a unique opportunity for elucidation of the structure of the conductive connexin pore and the molecular determinants of its ion permeation properties (conductance, current-voltage [I-V] relations, and charge selectivity). However, the crystal structure was incomplete, most notably lacking the coordinates of the N-terminal methionine residue, which resides within the pore, and also lacking two cytosolic domains. To allow computational studies for comparison with the known channel properties, we completed the structure. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) simulations of the completed and the published Cx26 hemichannel crystal structure indicate that the pore is too narrow to permit significant ion flux. The GCMC/BD simulations predict marked inward current rectification and almost perfect anion selectivity, both inconsistent with known channel properties. The completed structure was refined by all-atom molecular dynamics (MD) simulations (220 ns total) in an explicit solvent and POPC membrane system. These MD simulations produced an equilibrated structure with a larger minimal pore diameter, which decreased the height of the permeation barrier formed by the N terminus. GCMC/BD simulations of the MD-equilibrated structure yielded more appropriate single-channel conductance and less anion/cation selectivity. However, the simulations much more closely matched experimentally determined I-V relations when the charge effects of specific co- and posttranslational modifications of Cx26 previously identified by mass spectrometry were incorporated. We conclude that the average equilibrated structure obtained after MD simulations more closely represents the open Cx26 hemichannel structure than does the crystal structure, and that co- and posttranslational modifications of Cx26 hemichannels are likely to play an important physiological role by defining the conductance and ion selectivity of Cx26 channels. Furthermore, the simulations and data suggest that experimentally observed heterogeneity in Cx26 I-V relations can be accounted for by variation in co- and posttranslational modifications.

Charges dispersed over the permeation pathway determine the charge selectivity and conductance of a Cx32 chimeric hemichannel

Previous studies have shown that charge substitutions in the amino terminus of a chimeric connexin, Cx32*43E1, which forms unapposed hemichannels in Xenopus oocytes, can result in a threefold difference in unitary conductance and alter the direction and amount of open channel current rectification. Here, we determine the charge selectivity of Cx32*43E1 unapposed hemichannels containing negative and/or positive charge substitutions at the 2nd, 5th and 8th positions in the N-terminus. Unlike Cx32 intercellular channels, which are weakly anion selective, the Cx32*43E1 unapposed hemichannel is moderately cation selective. Cation selectivity is maximal when the extracellular surface of the channel is exposed to low ionic strength solutions implicating a region of negative charge in the first extracellular loop of Cx43 (Cx43E1) in influencing charge selectivity analogous to that reported. Negative charge substitutions at the 2nd, 5th and 8th positions in the intracellular N-terminus substantially increase the unitary conductance and cation selectivity of the chimeric hemichannel. Positive charge substitutions at the 5th position decrease unitary conductance and produce a non-selective channel while the presence of a positive charge at the 5th position and negative charge at the 2nd results in a channel with conductance similar to the parental channel but with greater preference for cations. We demonstrate that a cysteine substitution of the 8th residue in the N-terminus can be modified by a methanthiosulphonate reagent (MTSEA-biotin-X) indicating that this residue lines the aqueous pore at the intracellular entrance of the channel. The results indicate that charge selectivity of the Cx32*43E1 hemichannel can be determined by the combined actions of charges dispersed over the permeation pathway rather than by a defined region that acts as a charge selectivity filter.

Modulation of the conductance of a 2,2'-bipyridine-functionalized peptidic ion channel by Ni2+

An α-helical amphipathic peptide with the sequence H₂N-(LSSLLSL)₃-CONH₂ was obtained by solid phase synthesis and a 2,2′-bipyridine was coupled to its N-terminus, which allows complexation of Ni²⁺. Complexation of the 2,2′-bipyridine residues was proven by UV/Vis spectroscopy. The peptide helices were inserted into lipid bilayers (nano black lipid membranes, nano-BLMs) that suspend the pores of porous alumina substrates with a pore diameter of 60 nm by applying a potential difference. From single channel recordings, we were able to distinguish four distinct conductance states, which we attribute to an increasing number of peptide helices participating in the conducting helix bundle. Addition of Ni²⁺ in micromolar concentrations altered the conductance behaviour of the formed ion channels in nano-BLMs considerably. The first two conductance states appear much more prominent demonstrating that the complexation of bipyridine by Ni²⁺ results in a considerable confinement of the observed multiple conductance states. However, the conductance levels were independent of the presence of Ni²⁺. Moreover, from a detailed analysis of the open lifetimes of the channels, we conclude that the complexation of Ni²⁺ diminishes the frequency of channel events with larger open times.

Stochastic study of the effect of ionic strength on noncovalent interactions in protein pores

Salt plays a critical role in the physiological activities of cells. We show that ionic strength significantly affects the kinetics of noncovalent interactions in protein channels, as observed in stochastic studies of the transfer of various analytes through pores of wild-type and mutant alpha-hemolysin proteins. As the ionic strength increased, the association rate constant of electrostatic interactions was accelerated, whereas those of both hydrophobic and aromatic interactions were retarded. Dramatic decreases in the dissociation rate constants, and thus increases in the overall reaction formation constants, were observed for all noncovalent interactions studied. The results suggest that with the increase of salt concentration, the streaming potentials for all the protein pores decrease, whereas the preferential selectivities of the pores for either cations or anions drop. Furthermore, results also show that the salt effect on the rate of association of analytes to a pore differs significantly depending on the nature of the noncovalent interactions occurring in the protein channel. In addition to providing new insights into the nature of analyte-protein pore interactions, the salt-dependence of noncovalent interactions in protein pores observed provides a useful means to greatly enhance the sensitivity of the nanopore, which may find useful application in stochastic sensing.

Ca2+ selectivity of a chemically modified OmpF with reduced pore volume

We studied an E. coli OmpF mutant (LECE) containing both an EEEE-like locus, typical of Ca(2+) channels, and an accessible and reactive cysteine. After chemical modification with the cysteine-specific, negatively charged (-1e) reagents MTSES or glutathione, this LECE mutant was tested for Ca(2+) versus alkali metal selectivity. Selectivity was measured by conductance and zero-current potential. Conductance measurements showed that glutathione-modified LECE had reduced conductance at Ca(2+) mole fractions <10(-3). MTSES-modified LECE did not. Apparently, the LECE protein is (somehow) a better Ca(2+) chelator after modification with the larger glutathione. Zero-current potential measurements revealed a Ca(2+) versus monovalent cation selectivity that was highest in the presence of Li(+) and lowest in the presence of Cs(+). Our data clearly show that after the binding of Ca(2+) the LECE pore (even with the bulky glutathione present) is spacious enough to allow monovalent cations to pass. Theoretical computations based on density functional theory combined with Poisson-Nernst-Planck theory and a reduced pore model suggest a functional separation of ionic pathways in the pore, one that is specific for small and highly charged ions, and one that accepts preferentially large ions, such as Cs(+).

Chemical modification of the bacterial porin OmpF: gain of selectivity by volume reduction

OmpF is an essentially nonselective porin isolated from the outer membrane of Escherichia coli. Here we report on the manipulation of the ion selectivity of OmpF by chemical modification with MTS reagents (MTSET, MTSEA, and MTSES) and the (rather bulky) tripeptide glutathione, all cysteine specific. When recorded in a gradient of 0.1//1 M CaCl2 or 0.1//1 M NaCl, pH 7.4 solutions, measured reversal potentials of the most cation-selective modified mutants were (virtually) identical to the Nernst potential of Ca2+ or Na+. Compared to this full cation selectivity, the anion-selective modified mutants performed somewhat less but nevertheless showed high anion selectivity. We conclude that a low permanent charge in combination with a narrow pore can render the same selectivity as a highly charged but wider pore. These results favor the view that both the electrostatic potential arising form the fixed charge in the pore and the space available at the selectivity filter contribute to the charge selection (i.e., cation versus anion selectivity) of a biological ion channel.

Topological Equilibria of Ion Channel Peptides in Oriented Lipid Bilayers Revealed by 15N Solid-State NMR Spectroscopy

Ion channel peptides have been prepared by solid-phase peptide synthesis, labeled with 15N at selected sites, and reconstituted into oriented lipid bilayers. The (Leu-Ser-Ser-Leu-Leu-Ser-Leu)3-CONH2 peptide has previously been shown to exhibit well-defined and discrete ionic conductances when investigated by single-channel measurements [Lear, J. D., et al. (1988) Science 240, 1177]. Proton-decoupled 15N solid-state NMR spectroscopy indicates that (Leu-Ser-Ser-Leu-Leu-Ser-Leu)3-CONH2 preferentially aligns parallel to the membrane surface in excellent agreement with its amphipathic helical structure. However, by carefully choosing the conditions of the membrane environment, significant contributions that are indicative of transmembrane alignments become obvious in the 15N chemical shift solid-state NMR spectra. The data thereby provide experimental evidence for an equilibrium between in-plane and transmembrane-oriented helix configurations where the transmembrane and surface-oriented peptide fractions are in slow exchange. Similar topological equilibria are observed when the N-terminus of the LS21 peptide is acetylated. These observations provide experimental support for previous models, suggesting that the channels observed in single-channel conductance measurements are indeed formed by hexameric transmembrane helical bundles. In contrast, the shorter peptide (Leu-Ser-Ser-Leu-Leu-Ser-Leu)2-CONH2 is oriented parallel to the membrane surface under all conditions tested. This peptide exhibits erratic conductance changes when investigated by electrophysiological methods, probably because it is too short to span the lipid bilayer.

Engineering charge selectivity in model ion channels

Most ion channel proteins exhibit some degree of charge selectivity, that is, an ability to conduct ions of one charge more efficiently than ions of the opposite charge. The structural origins of charge selectivity remain incompletely understood despite recent advances in the determination of cation-selective and anion-selective channel protein structures. Helix bundle channels formed via self-assembly of the peptide alamethicin provide a tractable model system for exploring the structural basis of charge selectivity. We synthesized covalently-linked alamethicin dimers, with amino acid substitutions at position 18 [lysine (Lys), arginine (Arg), glutamine (Gln), 2,3-diaminopropionic acid (Dpr)] in each helix, to assess the role of this position as a charge-selectivity determinant in alamethicin channels. Of the position 18 substitutions investigated, the Lys derivative exhibited the greatest degree of anion selectivity. Arg-containing channels were slightly less anion-selective than Lys. Interestingly, Dpr channels showed cation selectivity nearly equivalent to that exhibited by the neutral Gln derivative. We suggest that this result is due to a wider pore diameter that permits a greater number of counter-ions leading to enhanced charge screening and a lower effective side-chain positive charge.

Proton conduction through the M2 protein of the influenza A virus; a quantitative, mechanistic analysis of experimental data

The M2 proton channel from influenza A virus forms proton-selective ion channels, which are the target of the drug amantadine. Here, existing experimental data are quantitatively examined for insights into mechanisms to account for the pH- and voltage-dependences of M2 proton conduction. The analysis shows that a model involving protonation equilibria of His37, including pH-dependent changes in the relative rates of diffusion on either side of the pore, is quantitatively able to account for recently reported electrophysiological data examining the pH- and voltage-dependences of Rostock and Weybridge strain M2 proton conduction.

Ion channels of alamethicin dimer N-terminally linked by disulfide bond

A covalent dimer of alamethicin Rf30 was synthesized by linking the N-termini by a disulfide bond. When the dimer peptides were added to the cis-side of a diphytanoyl PC membrane, macroscopic channel current was induced only at cis positive voltages. The single-channel recordings showed several conductance levels that were alternately stabilized. These results indicate that the dimer peptides form stable channels by N-terminal insertion like alamethicin and that most of the pores are assembled from even numbers of helices. Taking advantages of the long open duration of the dimer peptide channels, the current-voltage (I-V) relations of the single-channels were obtained by applying fast voltage ramps during the open states. The I-V relations showed rectification, such that current from the cis-side toward the trans-side is larger than that in the opposite direction. The intrinsic rectification is mainly attributed to the macro dipoles of parallel peptide helices surrounding a central pore.

Protein engineering modulates the transport properties and ion selectivity of the pores formed by staphylococcal gamma-haemolysins in lipid membranes

Staphylococcal gamma-haemolysins are bicomponent toxins in a family including other leucocidins and alpha-toxin. Two active toxins are formed combining HlgA or HlgC with HlgB. Both open pores in lipid membranes with conductance, current voltage characteristics and stability similar to alpha-toxin, but different selectivity (cation instead of anion). Structural analogies between gamma-haemolysins and alpha-toxin indicate the presence, at the pore entry, of a conserved region containing four positive charges in alpha-toxin, but either positive or negative in gamma-haemolysins. Four mutants were produced (HlgA D44K, HlgB D47K, HlgB D49K and HlgB D47K/D49K) converting those negative charges to positive in HlgA and HlgB. When all charges were positive, the pores had the same selectivity and conductance as alpha-toxin, suggesting that the cluster may form an entrance electrostatic filter. As mutated HlgC-HlgB pores were less affected, additional charges in the lumen of the pore were changed (HlgB E107Q, HlgB D121N, HlgB T136D and HlgA K108T). Removing a negative charge from the lumen made the selectivity of both HlgA-HlgB D121N and HlgC-HlgB D121N more anionic. Residue D121 of HlgB is compensated by a positive residue (HlgA K108) in the HlgA-HlgB pore, but isolated in the more cation-selective HlgC-HlgB pore. Interestingly, the pore formed by HlgA K108T-HlgB, in which the positive charge of HlgA was removed, was as cation selective as HlgC-HlgB. Meanwhile, the pore formed by HlgA K108T-HlgB D121N, in which the two charge changes compensated, retrieved the properties of wild-type HlgA-HlgB. We conclude that the conductance and selectivity of the gamma-haemolysin pores depend substantially on the presence and location of charged residues in the channel.

The first extracellular loop domain is a major determinant of charge selectivity in connexin46 channels

Intercellular channels formed of members of the gene family of connexins (Cxs) vary from being substantially cation selective to being anion selective. We took advantage of the ability of Cx46 to function as an unopposed hemichannel to examine the basis of Cx charge selectivity. Previously we showed Cx46 hemichannels to be large pores that predominantly conduct cations and inwardly rectify in symmetric salts, properties suggesting selectivity is influenced by fixed negative charges located toward the extracellular end of the pore. Here we demonstrate that high ionic strength solutions applied to the extracellular, but not the intracellular, side of Cx46 hemichannels substantially reduce the ratio of cation to anion permeability. Substitution of the first extracellular loop (E1) domain of Cx32, an anion-preferring Cx, reduces conductance, converts Cx46 from cation to anion preferring, and changes the I-V relation form inwardly to outwardly rectifying. These data suggest that fixed negative charges influencing selectivity in Cx46 are located in E1 and are substantially reduced and/or are replaced with positive charges from the Cx32 E1 sequence. Extending studies to Cx46 cell-cell channels, we show that they maintain a strong preference for cations, have a conductance nearly that expected by the series addition of hemichannels, but lack rectification in symmetric salts. These properties are consistent with preservation of the fixed charge region in E1 of hemichannels, which upon docking, become symmetrically placed near the center of the cell-cell channel pore. Furthermore, heterotypic cell-cell channels formed by pairing Cx46 with Cx32 or Cx43 rectify in symmetric salts in accordance with the differences in the charges we ascribed to E1. These data are consistent with charged residues in E1 facing the channel lumen and playing an important role in determining Cx channel conductance and selectivity.

Reversal of charge selectivity in transmembrane protein pores by using noncovalent molecular adapters

In this study, the charge selectivity of staphylococcal alpha-hemolysin (alphaHL), a bacterial pore-forming toxin, is manipulated by using cyclodextrins as noncovalent molecular adapters. Anion-selective versions of alphaHL, including the wild-type pore and various mutants, become more anion selective when beta-cyclodextrin (betaCD) is lodged within the channel lumen. By contrast, the negatively charged adapter, hepta-6-sulfato-beta-cyclodextrin (s(7)betaCD), produces cation selectivity. The cyclodextrin adapters have similar effects when placed in cation-selective mutant alphaHL pores. Most probably, hydrated Cl(-) ions partition into the central cavity of betaCD more readily than K(+) ions, whereas s(7)betaCD introduces a charged ring near the midpoint of the channel lumen and confers cation selectivity through electrostatic interactions. The molecular adapters generate permeability ratios (P(K+)/P(Cl-)) over a 200-fold range and should be useful in the de novo design of membrane channels both for basic studies of ion permeation and for applications in biotechnology.

Engineering charge selectivity in alamethicin channels

The peptide alamethicin provides a system for engineering ion channel charge selectivity. To define alamethicin charge selectivity experimentally, we measured single-channel current-voltage relationships in KCl gradients using covalently linked peptide dimers. Two factors were found to contribute to the charge selectivity of these channels: (i) the ionic strength of the surrounding solutions; and (ii) the distribution of fixed charge on the peptide. Native alamethicin channels exhibited either cation selectivity or anion selectivity depending on which end of the channel was at the low salt side of the membrane. When the glutamine residue at position 18 in the sequence was replaced with a lysine residue, an anion-selective channel was obtained regardless of which end of the channel was at the low salt side of the membrane.

Effective diffusion coefficients of K+ and Cl- ions in ion channel models

We have used molecular dynamics simulations, corresponding to a total simulation time of 11 ns, to investigate the effective short-time local diffusion coefficient of potassium and chloride ions in a series of model ion channels. These models, which include channels formed by the fungal peptide alamethicin, by a synthetic leucine-serine peptide, and by the pore-lining M2 helix bundle of the nicotinic acetylcholine receptor, have a range of different secondary structures, diameters and hydrophobicities. We find that the diffusion coefficients of both ions are appreciably reduced in the narrower channels, the extent of the reduction being similar for both the anionic and cationic species. This suggests that a difference in mobility cannot be the source of the ion selectivity exhibited by some of the channels (for example, the leucine-serine peptide). We find no evidence for a reduction in mobility of either ion in the nAChR model. These results are broadly in line with a previous similar study of Na+ ions, and may be useful in Poisson-Nernst-Planck, Eyring rate theory or Brownian dynamics calculations of channel conductance.

An anion-selective analogue of the channel-forming peptide alamethicin

The peptide alamethicin self-assembles to form helix bundle ion channels in membranes. Previous macroscopic measurements have shown that these channels are mildly cation-selective. Models indicate that a source of cation selectivity is a zone of partial negative charge toward the C-terminal end of the peptide. We synthesized an alamethicin derivative with a lysine in this zone (replacing the glutamine at position 18 in the sequence). Microscopic (single-channel) measurements demonstrate that dimeric alamethicin-lysine18 (alm-K18) forms mildly anion-selective channels under conditions where channels formed by the parent peptide are cation-selective. Long-range electrostatic interactions can explain the inversion of ion selectivity and the conductance properties of alamethicin channels.

Exploration of the structural features defining the conduction properties of a synthetic ion channel

The finite-difference Poisson-Boltzmann methodology was applied to a series of parallel, alpha-helical bundle models of the designed ion channel peptide Ac-(LSSLLSL)3-CONH2. This method is able to fully describe the current-voltage curves for this channel and quantitatively explains their cation selectivity and rectification. We examined a series of energy-minimized models representing different aggregation states, side-chain rotamers, and helical rotations, as well as an ensemble of structures from a molecular dynamics trajectory. Potential energies were computed for single, permeating K+ and Cl- ions at a series of positions along a central pathway through the models. A variable-electric-field Nernst-Planck electrodiffusion model was used, with two adjustable parameters representing the diffusion coefficients of K+ and Cl- to scale the individual ion current magnitudes. The ability of a given DelPhi potential profile to fit the experimental data depended strongly on the magnitude of the desolvation of the permeating ion. Below a pore radius of 3.8 A, the predicted profiles showed large energy barriers, and the experimental data could be fit only with unrealistically high values for the K+ and Cl- diffusion coefficients. For pore radii above 3.8 A, the desolvation energies were 2kT or less. The electrostatic calculations were sensitive to positioning of the Ser side chains, with the best fits associated with maximum exposure of the Ser side-chain hydroxyls to the pore. The backbone component was shown to be the major source of asymmetry in the DelPhi potential profiles. Only two of the energy-minimized structures were able to explain the experimental data, whereas an average of the dynamics structures gave excellent agreement with experimental results. Thus this method provides a promising approach to prediction of current-voltage curves from three-dimensional structures of ion channel proteins.

Dynamic properties of Na+ ions in models of ion channels: a molecular dynamics study

We present simulation results for the effective diffusion coefficients of a sodium ion in a series of model ion channels of different diameters and hydrophobicities, including models of alamethicin, a leucine-serine peptide, and the M2 helix bundle of the nicotinic acetylcholine receptor. The diffusion coefficient, which in the simulations has a value of 0.15(2) A2ps-1 in bulk water, is found to be reduced to as little as 0.02(1) A2ps-1 in the narrower channels, and to about 0.10(5) A2ps-1 in wider channels such as the nicotinic acetylcholine receptor. It is anticipated that this work will be useful in connection with calculations of channel conductivity using such techniques as the Poisson-Nernst-Planck equation, Eyring rate theory, or Brownian dynamics.

Synthesis and purification of hydrophobic peptides for use in biomimetic ion channels

The synthesis and subsequent purification of several hydrophobic peptides is described. These peptides include the 24-residue M3 transmembrane domain of the rat connexin 32 protein, a peptide sequence that contains only seven amino acids with hydrophilic side-chains (71% hydrophobic). Moreover, for comparison, a much smaller hydrophobic octapeptide, designed to exist with alpha-helical secondary structure, was also studied. Optimum conditions for the RP-HPLC purification of these peptides was dependent on peptide length and solubility properties.

Molecular dynamics study of the LS3 voltage-gated ion channel

Molecular dynamics calculations have been carried out on a model of the LS3 synthetic ion channel in a membrane mimetic environment. In the absence of an external electrostatic field, the LS3 channel, which consists of a bundle of six alpha-helices with sequence Ac-(LSSLLSL)3-CONH2, exhibits large structural fluctuations. However, in the presence of the field, the bundle adopts a well defined coiled-coil structure with an inner pore of water. The observed structural changes induced by the applied field are consistent with the proposed gating mechanism of the ion channel.

Molecular dynamics simulation of a synthetic ion channel

A molecular dynamics simulation has been performed on a synthetic membrane-spanning ion channel, consisting of four alpha-helical peptides, each of which is composed of the amino acids leucine (L) and serine (S), with the sequence Ac-(LSLLLSL)3-CONH2. This four-helix bundle has been shown experimentally to act as a proton-conducting channel in a membrane environment. In the present simulation, the channel was initially assembled as a parallel bundle in the octane portion of a phase-separated water/octane system, which provided a membrane-mimetic environment. An explicit reversible multiple-time-step integrator was used to generate a dynamical trajectory, a few nanoseconds in duration for this composite system on a parallel computer, under ambient conditions. After more than 1 ns, the four helices were found to adopt an associated dimer state with twofold symmetry, which evolved into a coiled-coil tetrameric structure with a left-handed twist. In the coiled-coil state, the polar serine side chains interact to form a layered structure with the core of the bundle filled with H2O. The dipoles of these H2O molecules tended to align opposite the net dipole of the peptide bundle. The calculated dipole relaxation function of the pore H2O molecules exhibits two reorientation times. One is approximately 3.2 ps, and the other is approximately 100 times longer. The diffusion coefficient of the pore H2O is about one-third of the bulk H2O value. The total dipole moment and the inertia tensor of the peptide bundle have been calculated and reveal slow (300 ps) collective oscillatory motions. Our results, which are based on a simple united atom force-field model, suggest that the function of this synthetic ion channel is likely inextricably coupled to its dynamical behavior.

Structure-function relationships in helix-bundle channels probed via total chemical synthesis of alamethicin dimers: effects of a Gln7 to Asn7 mutation

Alamethicin channels are prototypical helix bundles that may serve as tractable models for more complex protein ion channels. Solid-phase peptide synthesis of alamethicin analogues using FMOC-amino acid fluorides followed by chemical dimerization of these peptides facilitates structure-function studies of particular channel states in bilayer membranes. State 3 in particular, tentatively assigned to a hexameric helix bundle, is sufficiently long-lived that current-voltage measurements can be made during the lifetime of an individual channel opening. Molecular models of hexameric helix bundles, generated using restrained molecular dynamics with simulated annealing, indicate that a Gln7-->Asn7 (Q7-->N7) mutation will increase channel diameter locally. Experimentally, the conductance of state 3 of the N7-alm channel is found to be larger than that of the Q7-alm channel when ion flow is in the usual direction (cations entering the C-terminal end of the channel). When ion flow is in the opposite direction, no difference in the conductances of state 3 of Q7 and state 3 of N7 channels is observed. These results indicate that the effect of a change in pore diameter at position 7 is dependent on the magnitude of other barriers to permeation and that these barriers are voltage-dependent.

Intrinsic rectification of ion flux in alamethicin channels: studies with an alamethicin dimer

Covalent dimers of alamethicin form conducting structures with gating properties that permit measurement of current-voltage (I-V) relationships during the lifetime of a single channel. These I-V curves demonstrate that the alamethicin channel is a rectifier that passes current preferentially, with voltages of the same sign as that of the voltage that induced opening of the channel. The degree of rectification depends on the salt concentration; single-channel I-V relationships become almost linear in 3 M potassium chloride. These properties may be qualitatively understood by using Poisson-Nernst-Planck theory and a modeled structure of the alamethicin pore.

Ion-channels formed by hypelcins, antibiotic peptides, in planar bilayer lipid membranes

Ion-channel properties of native hypelcins (HP) A-I, A-V and B-V isolated from Hypocrea peltata and a synthetic analog, HP-A-Pheol, were studied in planar bilayer lipid membranes by a single-channel recording technique. The native and synthetic hypelcins formed ion-channels with three conductance levels for 3 mole dm(-3) KCl: < or = 0.09 nS at 225 mV (level 0, only detectable at voltages above 200 mV), approximately 0.6 nS at 150 mV (level 1, most common level) and approximately 3 nS at 150 mV (level 2). The effects of the C-terminal aminoalcohol on the channel properties were examined with HP-A-I, HP-A-V and HP-A-Pheol, whose C-termini are leucinol (Leuol), isoleucinol (Ileol) and phenylalaninol (Pheol), respectively. The substitution of Pheol for Leuol and Ileol prolonged the open channel lifetime. A comparison of HP-A-V (Gln18) and HP-B-V (Glu18) indicated that the carboxyl group at position 18 increased both the open channel lifetime and the magnitude of unitary channel conductance at each conductance level. The pores of level 1 showed poor ion-selectivity for K+ over Cl-. The selectivity order of alkali metal cations was Rb > or = Cs > or = K > Na > Li for level 1 and Cs > Rb > K > Na > Li for level 0. The unitary current-voltage characteristics showed non-linear relationships, which were simulated by a Nernst-Planck approach with a simple barrier model.

Ion selectivity predictions from a two-site permeation model for the cyclic nucleotide-gated channel of retinal rod cells

We developed a two-site, Eyring rate theory model of ionic permeation for cyclic nucleotide-gated channels (CNGCs). The parameters of the model were optimized by simultaneously fitting current-voltage (IV) data sets from excised photoreceptor patches in electrolyte solutions containing one or more of the following ions: Na+, Ca2+, Mg2+, and K+. The model accounted well for 1) the shape of the IV relations; 2) the binding affinity for Na+; 3) reversal potential values with single-sided additions of Ca2+ or Mg2+ and biionic KCl; and 4) the K1 and voltage dependence for divalent block from the cytoplasmic side of the channel. The differences between the predicted K1's for extracellular block by Ca2+ and Mg2+ and the values obtained from heterologous expression of only the alpha-subunit of the channel suggest that the beta-subunit or a cell-specific factor affects the interaction of divalent cations at the external but not the internal face of the channel. The model predicts concentration-dependent permeability ratios with single-sided addition of Ca2+ and Mg2+ and anomalous mole fraction effects under a limited set of conditions for both monovalent and divalent cations. Ca2+ and Mg2+ are predicted to carry 21% and 10%, respectively, of the total current in the retinal rod cell at -60 mV.

Permeation through an open channel: Poisson-Nernst-Planck theory of a synthetic ionic channel

The synthetic channel [acetyl-(LeuSerSerLeuLeuSerLeu)3-CONH2]6 (pore diameter approximately 8 A, length approximately 30 A) is a bundle of six alpha-helices with blocked termini. This simple channel has complex properties, which are difficult to explain, even qualitatively, by traditional theories: its single-channel currents rectify in symmetrical solutions and its selectivity (defined by reversal potential) is a sensitive function of bathing solution. These complex properties can be fit quantitatively if the channel has fixed charge at its ends, forming a kind of macrodipole, bracketing a central charged region, and the shielding of the fixed charges is described by the Poisson-Nernst-Planck (PNP) equations. PNP fits current voltage relations measured in 15 solutions with an r.m.s. error of 3.6% using four adjustable parameters: the diffusion coefficients in the channel's pore DK = 2.1 x 10(-6) and DCl = 2.6 x 10(-7) cm2/s; and the fixed charge at the ends of the channel of +/- 0.12e (with unequal densities 0.71 M = 0.021e/A on the N-side and -1.9 M = -0.058e/A on the C-side). The fixed charge in the central region is 0.31e (with density P2 = 0.47 M = 0.014e/A). In contrast to traditional theories, PNP computes the electric field in the open channel from all of the charges in the system, by a rapid and accurate numerical procedure. In essence, PNP is a theory of the shielding of fixed (i.e., permanent) charge of the channel by mobile charge and by the ionic atmosphere in and near the channel's pore. The theory fits a wide range of data because the ionic contents and potential profile in the channel change significantly with experimental conditions, as they must, if the channel simultaneously satisfies the Poisson and Nernst-Planck equations and boundary conditions. Qualitatively speaking, the theory shows that small changes in the ionic atmosphere of the channel (i.e., shielding) make big changes in the potential profile and even bigger changes in flux, because potential is a sensitive function of charge and shielding, and flux is an exponential function of potential.

A semi-microscopic Monte Carlo study of permeation energetics in a gramicidin-like channel: the origin of cation selectivity

The influence of a gramicidin-like channel former on ion free energy barriers is studied using Monte Carlo simulation. The model explicitly describes the ion, the water dipoles, and the peptide carbonyls; the remaining degrees of freedom, bulk electrolyte, non-polar lipid and peptide regions, and electronic (high frequency) permittivity, are treated in continuum terms. Contributions of the channel waters and peptide COs are studied both separately and collectively. We found that if constrained to their original orientations, the COs substantially increase the cationic permeation free energy; with or without water present, CO reorientation is crucial for ion-CO interaction to lower cation free energy barriers; the translocation free energy profiles for potassium-, rubidium-, and cesium-like cations exhibit no broad barriers; the lipid-bound peptide interacts more effectively with anions than cations; anionic translocation free energy profiles exhibit well defined maxima. Using experimental data to estimate transfer free energies of ions and water from bulk electrolyte to a non-polar dielectric (continuum lipid), we found reasonable ion permeation profiles; cations bind and permeate, whereas anions cannot enter the channel. Cation selectivity arises because, for ions of the same size and charge, anions bind hydration water more strongly.

Molecular mimicry in channel-protein structure

The approach to probing the sequence-structure relationship of ion-channel proteins using small peptides stems from the abundance of sequence information and the virtual absence of structures at atomic resolution. It is anticipated that model peptides may fold predictably into stable structures and reproduce functional properties of specific proteins. Model peptides are well suited to the application of NMR methods to determine protein structure in a membrane environment or to high-resolution X-ray diffraction analysis. It is timely to ask what we have learned through this strategy and where it may lead in our quest to understand the sequence-structure determinism.