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

Electrodiffusion Phenomena in Neuroscience and the Nernst–Planck–Poisson Equations

This work is aimed to give an electrochemical insight into the ionic transport phenomena in the cellular environment of organized brain tissue. The Nernst–Planck–Poisson (NPP) model is presented, and its applications in the description of electrodiffusion phenomena relevant in nanoscale neurophysiology are reviewed. These phenomena include: the signal propagation in neurons, the liquid junction potential in extracellular space, electrochemical transport in ion channels, the electrical potential distortions invisible to patch-clamp technique, and calcium transport through mitochondrial membrane. The limitations, as well as the extensions of the NPP model that allow us to overcome these limitations, are also discussed.

Simulations of Membrane-Disrupting Peptides I: Alamethicin Pore Stability and Spontaneous Insertion

An all-atom molecular dynamics simulation of the archetype barrel-stave alamethicin (alm) pore in a 1,2-dioleoyl-sn-glycero-3-phosphocholine bilayer at 313 K indicates that ∼7 μs is required for equilibration of a preformed 6-peptide pore; the pore remains stable for the duration of the remaining 7 μs of the trajectory, and the structure factors agree well with experiment. A 5 μs simulation of 10 surface-bound alm peptides shows significant peptide unfolding and some unbinding, but no insertion. Simulations at 363 and 413 K with a -0.2 V electric field yield peptide insertion in 1 μs. Insertion is initiated by the folding of residues 3-11 into an α-helix, and mediated by membrane water or by previously inserted peptides. The stability of five alm pore peptides at 413 K with a -0.2 V electric field demonstrates a significant preference for a transmembrane orientation. Hence, and in contrast to the cationic antimicrobial peptide described in the following article, alm shows a strong preference for the inserted over the surface-bound state.

Cyclodextrin-scaffolded alamethicin with remarkably efficient membrane permeabilizing properties and membrane current conductance

Bacterial resistance to classical antibiotics is a serious medical problem, which continues to grow. Small antimicrobial peptides represent a potential solution and are increasingly being developed as novel therapeutic agents. Many of these peptides owe their antibacterial activity to the formation of trans-membrane ion-channels resulting in cell lysis. However, to further develop the field of peptide antibiotics, a thorough understanding of their mechanism of action is needed. Alamethicin belongs to a class of peptides called peptaibols and represents one of these antimicrobial peptides. To examine the dynamics of assembly and to facilitate a thorough structural evaluation of the alamethicin ion-channels, we have applied click chemistry for the synthesis of templated alamethicin multimers covalently attached to cyclodextrin-scaffolds. Using oriented circular dichroism, calcein release assays, and single-channel current measurements, the α-helices of the templated multimers were demonstrated to insert into lipid bilayers forming highly efficient and remarkably stable ion-channels.

Metal-assisted channel stabilization: disposition of a single histidine on the N-terminus of alamethicin yields channels with extraordinarily long lifetimes

Alamethicin, a member of the peptaibol family of antibiotics, is a typical channel-forming peptide with a helical structure. The self-assembly of the peptide in the membranes yields voltage-dependent channels. In this study, three alamethicin analogs possessing a charged residue (His, Lys, or Glu) on their N-termini were designed with the expectation of stabilizing the transmembrane structure. A slight elongation of channel lifetime was observed for the Lys and Glu analogs. On the other hand, extensive stabilization of certain channel open states was observed for the His analog. This stabilization was predominantly observed in the presence of metal ions such as Zn(2+), suggesting that metal coordination with His facilitates the formation of a supramolecular assembly in the membranes. Channel stability was greatly diminished by acetylation of the N-terminal amino group, indicating that the N-terminal amino group also plays an important role in metal coordination.

Radio frequency rectification on membrane bound pores

Probing the interaction of biological systems with radio frequencies holds great promise for research and drug screening applications. While a common assumption is that biological systems do not operate at radio frequencies, we find that currents due to ion transport through channels and pores in cell membranes are in the pA to nA range. These values translate via the average current I = ne/tau(d) = nef to frequencies in the range of 1 MHz-1 GHz, where n is the average number of ions transported and tau(d) is the dwell time of the ions in the channel. It is thus desirable to have circuitry available which facilitates radio frequency spectroscopy of ion transport. This will yield real-time in vitro information on ion channel operation. Here we present measurements on the local interaction of a radio frequency signal with single ion channels and pores. We find radio frequency rectification and pumping on the channels and pores embedded in suspended bilipid membranes, recorded in direct current measurements. This electromagnetic modulation can be used to probe the dynamics of ion channel conformational changes.

Three-dimensional structure/hydrophobicity of latarcins specifies their mode of membrane activity

Latarcins, linear peptides from the Lachesana tarabaevi spider venom, exhibit a broad-spectrum antimicrobial activity, likely acting on the bacterial cytoplasmic membrane. We study their spatial structures and interaction with model membranes by a combination of experimental and theoretical methods to reveal the structure-activity relationship. In this work, a 26 amino acid peptide, Ltc1, was investigated. Its spatial structure in detergent micelles was determined by (1)H nuclear magnetic resonance (NMR) and refined by Monte Carlo simulations in an implicit water-octanol slab. The Ltc1 molecule was found to form a straight uninterrupted amphiphilic helix comprising 8-23 residues. A dye-leakage fluorescent assay and (31)P NMR spectroscopy established that the peptide does not induce the release of fluorescent marker nor deteriorate the bilayer structure of the membranes. The voltage-clamp technique showed that Ltc1 induces the current fluctuations through planar membranes when the sign of the applied potential coincides with the one across the bacterial inner membrane. This implies that Ltc1 acts on the membranes via a specific mechanism, which is different from the carpet mode demonstrated by another latarcin, Ltc2a, featuring a helix-hinge-helix structure with a hydrophobicity gradient along the peptide chain. In contrast, the hydrophobic surface of the Ltc1 helix is narrow-shaped and extends with no gradient along the axis. We have also disclosed a number of peptides, structurally homologous to Ltc1 and exhibiting similar membrane activity. This indicates that the hydrophobic pattern of the Ltc1 helix and related antimicrobial peptides specifies their activity mechanism. The latter assumes the formation of variable-sized lesions, which depend upon the potential across the membrane.

Using ion channel-forming peptides to quantify protein-ligand interactions

This paper proposes a method for sensing affinity interactions by triggering disruption of self-assembly of ion channel-forming peptides in planar lipid bilayers. It shows that the binding of a derivative of alamethicin carrying a covalently attached sulfonamide ligand to carbonic anhydrase II (CA II) resulted in the inhibition of ion channel conductance through the bilayer. We propose that the binding of the bulky CA II protein (MW approximately 30 kD) to the ion channel-forming peptides (MW approximately 2.5 kD) either reduced the tendency of these peptides to self-assemble into a pore or extracted them from the bilayer altogether. In both outcomes, the interactions between the protein and the ligand lead to a disruption of self-assembled pores. Addition of a competitive inhibitor, 4-carboxybenzenesulfonamide, to the solution released CA II from the alamethicin-sulfonamide conjugate and restored the current flow across the bilayer by allowing reassembly of the ion channels in the bilayer. Time-averaged recordings of the current over discrete time intervals made it possible to quantify this monovalent ligand binding interaction. This method gave a dissociation constant of approximately 2 microM for the binding of CA II to alamethicin-sulfonamide in the bilayer recording chamber: this value is consistent with a value obtained independently with CA II and a related sulfonamide derivative by isothermal titration calorimetry.

Channel-Forming Activity of Alamethicin: Effects of Covalent Tethering

In this short review article, the effects of covalent tethering of alamethicin molecules on channel-forming behavior are described. Broadly speaking, these chemical modifications have provided insight into all three aspects of channel behavior: the structure of the conducting state, the ion-selectivity and ion-permeation properties, and the voltage dependence. Each of these aspects are discussed in turn.

Structure of magainin and alamethicin in model membranes studied by x-ray reflectivity

We have investigated the structure of lipid bilayers containing varied molar ratios of different lipids and the antimicrobial peptides magainin and alamethicin. For this structural study, we have used x-ray reflectivity on highly aligned solid-supported multilamellar lipid membranes. The reflectivity curves have been analyzed by semi-kinematical reflectivity theory modeling the bilayer density profile rho(z). Model simulations of the reflectivity curves cover a large range of vertical momentum transfer q(z), and yield excellent agreement between data and theory. The structural changes observed as a function of the molar peptide/lipid concentration P/L are discussed in a comparative way.

Ion channels of N-terminally linked alamethicin dimers: enhancement of cation-selectivity by substitution of Glu for Gln at position 7

Alamethicin forms voltage-gated ion channels that have moderate cation-selectivity. The enhancement of the cation-selectivity by introducing negatively charged residues at positions 7 and 18 has been studied using the tethered homodimers of alamethicin with Q7 and E18 (di-alm-Q7E18) and its analog with E7 and Q18 (di-alm-E7Q18). In the dimeric peptides, monomer peptides are linked at the N-termini by a disulfide bond. Both the peptides formed long lasting ion channels at cis-positive voltages when added to the cis-side membrane. Their long open duration enabled us to obtain current-voltage (I-V(m)) relations and reversal potentials at the single-channel level by applying a voltage ramp during the channel opening. The reversal potentials measured in asymmetric KCl solutions indicated that ionized E7 provided strong cation-selectivity, whereas ionized E18 little influenced the charge selectivity. This was also the case for the macroscopic charge selectivity determined from the reversal potentials obtained by the macroscopic I-V(m) measurements. The results are accounted for by stronger electrostatic interactions between permeant ions and negatively charged residues at the narrowest part of the pore than at the pore mouth.

The role of extra-membranous inter-helical loops in helix-helix interactions

The effect of a short loop connecting two transmembrane alpha-helices was studied using molecular dynamics simulations. Helices F and G from bacteriorhodopsin and two corresponding polyalanine helices were embedded in octane and POPC membranes in a transmembrane configuration both with and without the inter-helical loop. The results indicate that the membrane environment and the sequence of the loop are more influential on the dynamics and structure of the motif than the presence of a loop as such, at least for the time-scales investigated. The four residues in the FG loop are stabilized by four hydrogen bonds. These hydrogen bonds are not present in the polyalanine loop, causing it to be more flexible than the FG loop. This effect was observed independently of the protein environment, stressing the importance of the sequence. The structural analysis indicates that the loop has weak stabilizing properties in all environments. The stabilization due to the presence of the loop was strongest in a simulation of the FG fragment in a membrane-mimetic octane slab. In the simulations of the helix-loop-helix motif embedded in an explicit lipid bilayer model, the lipid bilayer interface compensates to a large extent for the absence of the loop.

Conformation of peptides in lipid membranes studied by x-ray grazing incidence scattering

Although the antimicrobial, fungal peptide alamethicin has been extensively studied, the conformation of the peptide and the interaction with lipid bilayers as well as the mechanism of channel gating are still not completely clear. As opposed to studies of the crystalline state, the polypeptide structures in the environment of fluid bilayers are difficult to probe. We have investigated the conformation of alamethicin in highly aligned stacks of model lipid membranes by synchrotron-based x-ray scattering. The (wide-angle) scattering distribution has been measured by reciprocal space mappings. A pronounced scattering signal is observed in samples of high molar peptide/lipid ratio which is distinctly different from the scattering distribution of an ideal helix in the transmembrane state. Beyond simple models of ideal helices, the data is analyzed in terms of models based on atomic coordinates from the Brookhaven Protein Data Bank, as well as from published molecular dynamics simulations. The results can be explained by assuming a wide distribution of helix tilt angles with respect to the membrane normal and a partial insertion of the N-terminus into the membrane.

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.

The implementation of slab geometry for membrane-channel molecular dynamics simulations

Slab geometric boundary conditions are applied in the molecular dynamics simulation of a simple membrane-channel system. The results of the simulation were compared to those of an analogous system using normal three-dimensional periodic boundary conditions. Analysis of the dynamics and electrostatics of the system show that slab geometric periodicity eliminates the artificial bulk water orientational polarization that is present while using normal three-dimensional periodicity. Furthermore, even though the water occupancy and volume of our simple channel is the same when using either method, the electrostatic properties are considerably different when using slab geometry. In particular, the orientational polarization of water is seen to be different in the interior of the channel. This gives rise to a markedly different electric field within the channel. We discuss the implications of slab geometry for the future simulation of this type of system and for the study of channel transport properties.

N-terminal insertion of alamethicin in channel formation studied using its covalent dimer N-terminally linked by disulfide bond

Alamethicin is supposed to form helix-bundle-type channels by inserting the N terminus into bilayer lipid membranes under sufficient voltages. The N-terminal insertion has been studied with an alamethicin dimer (di-alm) N-terminally linked by a disulfide bond and by the asymmetric addition of dithiothreitol (DTT) and tetrathionate (TT) to the membrane. When di-alm was added to the cis-side membrane, it forms long-lasting channels with the lifetime tau of about 100 ms at cis-positive voltages. The lifetime was reduced to a few milliseconds by addition of DTT to the cis-side membrane, indicating that most of the channels were formed by the monomers (alm-SH) that resulted from the cleavage of the disulfide bond in di-alm. The succeeding addition of TT to the trans-side produced channels of tau=10-20 ms besides the channels of alm-SH. The results suggested that TT reacted with the N-terminal thiol group of alm-SH located at the trans-side of the membrane to alter the lifetime. The N-terminal insertion of alamethicin helices by voltage activation, therefore, was confirmed.

Understanding pH-dependent selectivity of alamethicin K18 channels by computer simulation

Alamethicin K18 is a covalently linked alamethicin dimer in which the glutamine residue at position 18 in each helix has been replaced by a lysine residue. As described in previous work, channels formed by this peptide show pH-dependent selectivity. The maximum anion selectivity of the putative octameric conducting state is obtained at pH 7 or lower. Inasmuch as no change in selectivity is seen between pH 7 and pH 3, and because protons are expected to be in equilibrium with the open state of the channel during a selectivity measurement, the channel is believed to be fully charged (i.e., all eight lysines protonated) at pH 7. In an effort to understand how such a highly charged channel structure is stable in membranes and why it is not more selective for anions, we have performed a number of computer simulations of the system. Molecular dynamics simulations of 10 ns each of the octameric bundle in a lipid bilayer environment are presented, with either zero, four, or eight lysines charged in the absence of salt, and with eight lysines charged in the presence of 0.5 M and 1 M KCl. When no salt is present and all lysines are charged, on average 1.9 Cl(-) ions are inside the channel and the channel significantly deforms. With 0.5 M KCl present, 2.9 Cl(-) ions are inside the channel. With 1 M KCl present, four Cl(-) ions are present and the channel maintains a regular structure. Poisson-Boltzmann calculations on models of the octameric channel also predict an average of 2-4 Cl(-) ions near the lysine residues as a function of ionic strength. These counterions lower the apparent charge of the channel, which may underlie the decrease in selectivity observed experimentally with increasing salt concentrations. We suggest that to increase the selectivity of Alm K18 channels, positive charges could be engineered in a narrower part of the channel.

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.

Artificial ion channels showing rectified current behavior

Voltage-dependent artificial ion channels 3 and 4 were synthesized. Two cholic acid derivatives were connected through a m-xylylene dicarbamate unit at 3-hydroxyl groups. Asymmetries were introduced by terminal hydrophilic groups, carboxylic acid and phosphoric acid for 3 and hydroxyl and carboxylic acid for 4. Under basic conditions, these headgroups in 3 and 4 are expected to be dissociate into -1/-2 (pH 8.2) and 0/-1 (pH 7.2), respectively. Single ion channel properties were examined by a planar bilayer lipid membrane method under symmetrical 500 mM KCl at pH 8.2 or 7.2. When 3 and 4 were introduced into the bilayer membrane under application of positive voltage (a positive-shift method), the current values at positive applied voltage were larger than the corresponding ones at the negative applied voltage. The current-voltage plots were fitted by curves through a zero point to show clear rectification properties. The direction of rectification could be controlled by positive- or negative-shift methods. Vectorial alignment of terminal headgroup charges by the voltage-shift incorporation is essential for giving voltage-dependent rectified ion channels.

Voltage-dependent insertion of alamethicin at phospholipid/water and octane/water interfaces

Understanding the binding and insertion of peptides in lipid bilayers is a prerequisite for understanding phenomena such as antimicrobial activity and membrane-protein folding. We describe molecular dynamics simulations of the antimicrobial peptide alamethicin in lipid/water and octane/water environments, taking into account an external electric field to mimic the membrane potential. At cis-positive potentials, alamethicin does not insert into a phospholipid bilayer in 10 ns of simulation, due to the slow dynamics of the peptide and lipids. However, in octane N-terminal insertion occurs at field strengths from 0.33 V/nm and higher, in simulations of up to 100 ns duration. Insertion of alamethicin occurs in two steps, corresponding to desolvation of the Gln7 side chain, and the backbone of Aib10 and Gly11. The proline induced helix kink angle does not change significantly during insertion. Polyalanine and alamethicin form stable helices both when inserted in octane and at the water/octane interface, where they partition in the same location. In water, both polyalanine and alamethicin partially unfold in multiple simulations. We present a detailed analysis of the insertion of alamethicin into the octane slab and the influence of the external field on the peptide structure. Our findings give new insight into the mechanism of channel formation by alamethicin and the structure and dynamics of membrane-associated helices.

Simulations of ion channels--watching ions and water move

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.

Protonation of lysine residues inverts cation/anion selectivity in a model channel

A dimeric alamethicin analog with lysine at position 18 in the sequence (alm-K18) was previously shown to form stable anion-selective channels in membranes at pH 7.0 [Starostin, A. V., R. Butan, V. Borisenko, D. A. James, H. Wenschuh, M. S. Sansom, and G. A. Woolley. 1999. Biochemistry. 38:6144-6150]. To probe the charge state of the conducting channel and how this might influence cation versus anion selectivity, we performed a series of single-channel selectivity measurements at different pH values. At pH 7.0 and below, only anion-selective channels were found with P(K(+))/P(Cl(-)) = 0. 25. From pH 8-10, a mixture of anion-selective, non-selective, and cation-selective channels was found. At pH > 11 only cation-selective channels were found with P(K(+))/P(Cl(-)) = 4. In contrast, native alamethicin-Q18 channels (with Gln in place of Lys at position 18) were cation-selective (P(K(+))/P(Cl(-)) = 4) at all pH values. Continuum electrostatics calculations were then carried out using an octameric model of the alm-K18 channel embedded in a low dielectric slab to simulate a membrane. Although the calculations can account for the apparent pK(a) of the channel, they fail to correctly predict the degree of selectivity. Although a switch from cation- to anion-selectivity as the channel becomes protonated is indicated, the degree of anion-selectivity is severely overestimated, suggesting that the continuum approach does not adequately represent some aspect of the electrostatics of permeation in these channels. Side-chain conformational changes upon protonation, conformational changes, and deprotonation caused by permeating cations and counterion binding by lysine residues upon protonation are considered as possible sources of the overestimation.

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.

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 alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.

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.

A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel

A lattice relaxation algorithm is developed to solve the Poisson-Nernst-Planck (PNP) equations for ion transport through arbitrary three-dimensional volumes. Calculations of systems characterized by simple parallel plate and cylindrical pore geometries are presented in order to calibrate the accuracy of the method. A study of ion transport through gramicidin A dimer is carried out within this PNP framework. Good agreement with experimental measurements is obtained. Strengths and weaknesses of the PNP approach are discussed.

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.

Electrostatics and the ion selectivity of ligand-gated channels

The nicotinic acetylcholine receptor (nAChR) is a cation-selective ion channel that opens in response to acetylcholine binding. The related glycine receptor (GlyR) is anion selective. The pore-lining domain of each protein may be modeled as a bundle of five parallel M2 helices. Models of the pore-lining domains of homopentameric nAChR and GlyR have been used in continuum electrostatics calculations to probe the origins of ion selectivity. Calculated pKA values suggest that "rings" of acidic or basic side chains at the mouths of the nAChR or GlyR M2 helix bundles, respectively, may not be fully ionized. In particular, for the nAChR the ring of glutamate side chains at the extracellular mouth of the pore is predicted to be largely protonated at neutral pH, whereas those glutamate side chains in the intracellular and intermediate rings (at the opposite mouth of the pore) are predicted to be fully ionized. Inclusion of the other domains of each protein represented as an irregular cylindrical tube in which the M2 bundles are embedded suggests that both the M2 helices and the extramembrane domains play significant roles in determining ion selectivity.

Ionophoric activity of the antibiotic peptaibol trichorzin PA VI: a 23Na- and 35Cl-NMR study

Trichorzin PA VI (Ac Aib1 Ser Ala Aib Iva Gln Aib Val Aib Gly10 Leu Aib Pro Leu Aib Aib Gln Pheol18) is one of the seven main peptaibols forming the natural antibiotic 18-residue peptide mixture biosynthesised by a Trichoderma harzianum strain. Trichorzins exhibit antimycoplasmic activity resulting from membrane permeability perturbations. The membrane permeabilisation process by trichorzin PA VI has been examined in egg yolk phosphatidylcholine large unilamellar vesicles (LUV) and under conditions of ionic equilibrium by 23Na- and 35Cl-NMR experiments conducted in the presence of a chemical shift reagent and a relaxation agent, respectively. In such conditions, trichorzin PA VI exchanges both cations and anions across the vesicle bilayers, indicating the absence of ion- and charge-selectivity, in contrast to antibiotic ionophores, such as monensin or nigericin; the Na+ exchange is not influenced by the ionic strength. The kinetics of the Na+ exchange have been found to be third to fourth order with respect to the peptide concentration. The permeabilisation process of liposomes has been shown to be due to the formation of aggregates of three to four helical peptide monomers arranged into a supramolecular complex including presumably lipid molecules and forming a badly-defined pore in the bilayer. The major mechanism by which ions may exchange through the bilayer involves a long-lasting opening of the pores allowing complete exchange of the internal and external media in an 'all or nothing mode'.

Modelling and simulation of ion channels: applications to the nicotinic acetylcholine receptor

Molecular dynamics simulations with experimentally derived restraints have been used to develop atomic models of M2 helix bundles forming the pore-lining domains of the nicotinic acetylcholine receptor and related ligand-gated ion channels. M2 helix bundles have been used in microscopic simulations of the dynamics and energetics of water and ions within an ion channel. Translational and rotational motion of water are restricted within the pore, and water dipoles are aligned relative to the pore axis by the surrounding helix dipoles. Potential energy profiles for translation of a Na+ ion along the pore suggest that the protein and water components of the interaction energy exert an opposing effect on the ion, resulting in a relatively flat profile which favors cation permeation. Empirical conductance calculations based on a pore radius profile suggest that the M2 helix model is consistent with a single channel conductance of ca. 50 pS. Continuum electrostatics calculations indicate that a ring of glutamate residues at the cytoplasmic mouth of the alpha 7 nicotinic receptor M2 helix bundle may not be fully ionized. A simplified model of the remainder of the channel protein when added to the M2 helix bundle plays a significant role in enhancing the ion selectivity of the channel.

Models and simulations of ion channels and related membrane proteins

The past year has seen major advances in our understanding of ion channels, resulting from molecular dynamics simulations and modelling studies. Simulations of gramicidin have revealed that proton conduction along a water wire is limited by the dynamics of water reorientation. Plausible models are now available for a number of other channels, including alamethicin, the influenza A virus M2 protein, and the pore domains of the nicotinic acetylcholine receptor and Kv channels. Molecular dynamics simulations and continuum calculations have revealed some of the subtleties of the interactions between transmembrane helices and their lipid bilayer environment.