Ion channels: a first view of K+ channels in atomic glory

Crystal structures have been solved for the transbilayer pore domain of a bacterial K+ channel and the tetramerisation domain of voltage-gated K+ channel. These provide our first real structural insights into possible mechanisms of ion selectivity and permeation for K+ channels.

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.

Functional characterization of a mutated chicken alpha7 nicotinic acetylcholine receptor subunit with a leucine residue inserted in transmembrane domain 2

1. Site-directed mutagenesis was used to create an altered form of the chicken alpha7 nicotinic acetylcholine (ACh) receptor subunit (alpha7x61) in which a leucine residue was inserted between residues Leu9' and Ser10' in transmembrane domain 2. The properties of alpha7x61 receptors are distinct from those of the wild-type receptor. 2. Oocytes expressing wild-type alpha7 receptors responded to 10 microM nicotine with rapid inward currents that desensitized with a time-constant of 710+/-409 ms (mean+/-s.e.mean, n=5). However in alpha7x61 receptors 10 microM nicotine resulted in slower onset inward currents that desensitized with a time-constant of 5684+/-3403 ms (mean+/-s.e.mean, n = 4). No significant difference in the apparent affinity of nicotine or acetylcholine between mutant and wild-type receptors was observed. Dihydro-beta-erythroidine (DHbetaE) acted as an antagonist on both receptors. 3. Molecular modelling of the alpha7x61 receptor channel pore formed by a bundle of M2 alpha-helices suggested that three of the channel lining residues would be altered by the leucine insertion i.e.; Ser10 would be replaced by the leucine insertion, Val13' and Phe14' would be replaced, by Thr and Val, respectively. 4 When present in the LEV-1 nicotinic ACh receptor subunit from Caenorhabditis elegans the same alteration conferred resistance to levamisole anthelmintic drug. Levamisole blocked responses to nicotine of wild-type and alpha7x61 receptors. However, block was more dependent on membrane potential for the alpha7x61 receptors. 5. We conclude that the leucine insertion in transmembrane domain 2 has the unusual effect of slowing desensitization without altering apparent agonist affinity.

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.

Protein-water-ion interactions in a model of the pore domain of a potassium channel: a simulation study

A model of the selectivity filter of a voltage-gated K+ (Kv) channel formed by an eight-stranded beta-barrel is compared with physiological properties of the channel. Continuum electrostatic calculations suggest that only two of the eight Asp sidechains at the extracellular mouth of the pore will ionise. A ring of four Tyr sidechains forms the narrowest region of the pore. Molecular dynamic simulations of the potential energy of a K+ ion as translated along the model pore indicate that the two ionised Asp sidechains and the hydroxyl groups of the Tyr sidechains stabilise the partially desolvated ion as it passes through the narrowest region.

Ion channel stability and hydrogen bonding. Molecular modelling of channels formed by synthetic alamethicin analogues

Several analogues of the channel-forming peptaibol alamethicin have been demonstrated to exhibit faster switching between channel substates than does unmodified alamethicin. Molecular modelling studies are used to explore the possible molecular basis of these differences. Models of channels formed by alamethicin analogues were generated by restrained molecular dynamics in vacuo and refined by short molecular dynamics simulations with water molecules within and at either mouth of the channel. A decrease in backbone solvation was found to correlate with a decrease in open channel stability between alamethicin and an analogue in which all alpha-amino-isobutyric acid residues of alamethicin were replaced by leucine. A decrease in the extent of hydrogen-bonding at residue 7 correlates with lower open channel stabilities of analogues in which the glutamine at position 7 was replaced by smaller polar sidechains. These two observations indicate the importance of alamethicin/water H-bonds in stabilizing the open channel.

Channels formed by the transmembrane helix of phospholamban: a simulation study

Phospholamban is a small membrane protein which can form cation selective ion channels in lipid bilayers. Each subunit contains a single, largely hydrophobic transmembrane helix. The helices are thought to assemble as a pentameric and approximately parallel bundle surrounding a central pore. A model of this assembly (PDB code IPSL) has been used as the starting point for molecular dynamics (MD) simulations of a system consisting of the pentameric helix bundle, plus 217 water molecules located within and at either mouth of the pore. Interhelix distance restraints were employed to maintain the integrity of the helix bundle during a 500 ps MD simulation. Water molecules within the pore exhibited reduced diffusional and rotational mobility. Interactions between the alpha-helix dipoles and the water dipoles, the latter aligned anti-parallel to the former, contribute to the stability of the system. Analysis of the potential energy of interaction of a K+ ion as it was moved through the pore suggested that unfavourable interactions of the cation with the aligned helix dipoles at the N-terminal mouth were overcome by favourable ion-water interactions. Comparable analysis for a Cl ion revealed that the ion-(pore + water) interactions were unfavourable along the whole of the pore, increasingly so from the N- to the C-terminal mouth. Overall, the interaction energy profiles were consistent with a pore selective for cations over anions. Pore radius profiles were used to predict a channel conductance of 50 to 70 ps in 0.2 M KCl, which compares well with an experimental value of 100 ps.

Secondary structure of an isolated P-region from the voltage-gated sodium channel: a molecular modelling/dynamics study

Conformational studies of synthetic peptides corresponding to the pore-forming regions of voltage-gated sodium channels show a high tendency for beta-sheet conformation when interacting with lipid vesicles, as revealed by circular dichroism and infrared spectroscopy. These observations have guided our choice of possible molecular models for the P-region peptide of domain II of voltage-gated sodium channels: three alternative beta-hairpins, with differing turn assignments, or an alpha-helical hairpin. After generation of models by distance geometry-based methods, molecular dynamics (MD) simulations were run. in the absence of explicit solvent molecules but employing three different dielectric constants, to explore possible conformational preferences. The simulations in the different dielectric environments suggest that a 4-residue turn with the sequence LCGE yields more stable beta-hairpins. The MD results suggest that the SS1 part of the peptide may be more stable as an alpha-helix, whereas the SS2 part tends to adopt a beta-conformation.

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.

The dielectric properties of water within model transbilayer pores

Ion channels contain extended columns of water molecules within their transbilayer pores. The dynamic properties of such intrapore water have been shown to differ from those of water in its bulk state. In previous molecular dynamics simulations of two classes of model pore (parallel bundles of Ala20 alpha-helices and antiparallel barrels of Ala10 beta-strands), a substantially reduced translational and rotational mobility of waters was observed within the pore relative to bulk water. Molecular dynamics simulations in the presence of a transpore electrostatic field (i.e., a voltage drop along the pore axis) have been used to estimate the resultant polarization (due to reorientation) of the intrapore water, and hence to determine the local dielectric behavior within the pore. It is shown that the local dielectric constant of water within a pore is reduced for models formed by parallel alpha-helix bundles, but not by those formed by beta-barrels. This result is discussed in the context of electrostatics calculations of ion permeation through channels, and the effect of the local dielectric of water within a helix bundle pore is illustrated with a simple Poisson-Boltzmann calculation.

Molecular dynamics study of water and Na+ ions in models of the pore region of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is an integral membrane protein that forms ligand-gated and cation-selective channels. The central pore is lined by a bundle of five approximately parallel M2 helices, one from each subunit. Candidate model structures of the solvated pore region of a homopentameric (alpha7)5 nAChR channel in the open state, and in two possible forms of the closed state, have been studied using molecular dynamics simulations with restraining potentials. It is found that the mobility of the water is substantially lower within the pore than in bulk, and the water molecules become aligned with the M2 helix dipoles. Hydrogen-bonding patterns in the pore, especially around pore-lining charged and hydrophilic residues, and around exposed regions of the helix backbone, have been determined. Initial studies of systems containing both water and sodium ions together within the pore region have also been conducted. A sodium ion has been introduced into the solvated models at various points along the pore axis and its energy profile evaluated. It is found that the ion causes only a local perturbation of the water structure. The results of these calculations have been used to examine the effectiveness of the central ring of leucines as a component of a gate in the closed-channel model.

The pore-lining region of shaker voltage-gated potassium channels: comparison of beta-barrel and alpha-helix bundle models

Although there is a large body of site-directed mutagenesis data that identify the pore-lining sequence of the voltage-gated potassium channel, the structure of this region remains unknown. We have interpreted the available biochemical data as a set of topological and orientational restraints and employed these restraints to produce molecular models of the potassium channel pore region, H5. The H5 sequence has been modeled either as a tetramer of membrane-spanning beta-hairpins, thus producing an eight-stranded beta-barrel, or as a tetramer of incompletely membrane-spanning alpha-helical hairpins, thus producing an eight-staved alpha-helix bundle. In total, restraints-directed modeling has produced 40 different configurations of the beta-barrel model, each configuration comprising an ensemble of 20 structures, and 24 different configurations of the alpha-helix bundle model, each comprising an ensemble of 24 structures. Thus, over 1300 model structures for H5 have been generated. Configurations have been ranked on the basis of their predicted pore properties and on the extent of their agreement with the biochemical data. This ranking is employed to identify particular configurations of H5 that may be explored further as models of the pore-lining region of the voltage-gated potassium channel pore.

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.

The influenza A virus M2 channel: a molecular modeling and simulation study

The M2 protein of influenza virus forms ion channels activated by low pH which are proton permeable and play a key role in the life cycle of the virus. M2 is a 97-residue integral membrane protein containing a single transmembrane (TM) helix. M2 is present as disulfide-linked homotetramers. The TM domain of M2 has been modeled as a bundle of four parallel M2 helices. The helix bundle forms a left-handed supercoil surrounding a central pore. Residue H37 has been implicated in the mechanism of low-pH activation of the channel. Models generated with H37 in a fully deprotonated state exhibit a pore occluded by a ring of H37 side chains oriented toward the lumen of the pore. Models with H37 in a fully protonated state no longer exhibit such occlusion of the pore, as the H37 side chains adopt a more interfacial location. Extended molecular dynamics simulations with water molecules within and at the mouths of the pores support this distinction between the H37-deprotonated and H37-protonated models. These simulations suggest that only in the H37-protonated model is there a continuous column of water extending the entire length of the central pore. A mechanism for activation of M2 by low pH is presented in which the H37-deprotonated model corresponds to the "closed" form of the channel, while the H37-protonated model corresponds to the "open" form. A switch from the closed to the open form of the channel occurs if H37 is protonated midway through a simulation. The open channel is suggested to contain a wire of H-bonded water molecules which enables proton permeability.

Alamethicin channels - modelling via restrained molecular dynamics simulations

Alamethicin channels have been modelled as approximately parallel bundles of transbilayer helices containing between N = 4 and 8 helices per bundle. Initial models were generated by in vacuo restrained molecular dynamics (MD) simulations, and were refined by 60 ps MD simulations with water molecules present within and at the mouths of the central pore. The helix bundles were stabilized by networks of H-bonds between intra-pore water molecules and Gln-7 side-chains. Channel conductances were predicted on the basis of pore radius profiles, and suggested that the N = 4 bundle formed an occluded pore, whereas pores with N > or = 5 helices per bundle were open. Continuum electrostatics calculations suggested that the N = 6 pore is cation-selective, whereas pores with N > or = 7 helices per bundle were predicted to be somewhat less ion-selective.

Ion channels formed by HIV-1 Vpu: a modelling and simulation study

Vpu is an oligomeric integral membrane protein encoded by HIV-1 which forms ion channels, each subunit of which contains a single transmembrane helix. Models of Vpu channels formed by bundles of N = 4, 5 or 6 transmembrane helices have been developed by restrained molecular dynamics and refined by 100 ps simulations with water molecules within the pore. Pore radius profiles and conductance predictions suggest that the N = 5 model corresponds to the predominant channel conductance level of the channel. Potential energy profiles for translation of Na+ or Cl- ions along the Vpu N = 5 pore are consistent with the weak cation selectivity of Vpu channels.

A novel method for structure-based prediction of ion channel conductance properties

A rapid and easy-to-use method of predicting the conductance of an ion channel from its three-dimensional structure is presented. The method combines the pore dimensions of the channel as measured in the HOLE program with an Ohmic model of conductance. An empirically based correction factor is then applied. The method yielded good results for six experimental channel structures (none of which were included in the training set) with predictions accurate to within an average factor of 1.62 to the true values. The predictive r2 was equal to 0.90, which is indicative of a good predictive ability. The procedure is used to validate model structures of alamethicin and phospholamban. Two genuine predictions for the conductance of channels with known structure but without reported conductances are given. A modification of the procedure that calculates the expected results for the effect of the addition of nonelectrolyte polymers on conductance is set out. Results for a cholera toxin B-subunit crystal structure agree well with the measured values. The difficulty in interpreting such studies is discussed, with the conclusion that measurements on channels of known structure are required.

Ferrocenoyl derivatives of alamethicin: redox-sensitive ion channels

The synthesis and single-channel characterization of two redox-active C-terminal derivatives of alamethicin are herein described. The reduced [Fe(II)] forms of ferrocenoyl-alamethicin (Fc-ALM) and 1'-carboxyferrocenoyl-alamethicin (cFc-ALM) are shown to form voltage-dependent ion channels at cis positive potentials in planar lipid bilayers (PLB) with conductance properties similar to those of alamethicin. In situ oxidation of Fc-ALM [to Fe(III)] in the PLB apparatus causes a time-dependent elimination of channel openings, which can be restored by an increase in the transbilayer potential. In contrast, oxidation of cFc-ALM leads to the formation of shorter-lived channels. Pretreatment of the ferrocenoyl peptides with oxidizing agent alters their single-channel properties in a qualitatively similar manner, establishing that the changes in channel properties in the presence of oxidizing agents are due specifically to ferrocenoyl oxidation. We suggest that the redox sensitivity of these ferrocene-containing ion channels may be governed by a combination of the following factors: (1) changes in hydrophobicity; (2) alteration of peptide molecular dipole; and (3) alterations in tendencies toward self-association. However, oxidation induced changes in peptide conformation cannot be ruled out. Our results provide evidence that it is possible to engineer channel-forming peptides that respond to specific changes in the chemical environment.

Simulation studies of alamethicin-bilayer interactions

Alamethicin is an alpha-helical peptide that forms voltage-activated ion channels. Experimental data suggest that channel formation occurs via voltage-dependent insertion of alamethicin helices into lipid bilayers, followed by self-assembly of inserted helices to form a parallel helix bundle. Changes in the kink angle of the alamethicin helix about its central proline residue have also been suggested to play a role in channel gating. Alamethicin helices generated by simulated annealing and restrained molecular dynamics adopt a kink angle similar to that in the x-ray crystal structure, even if such simulations start with an idealized unkinked helix. This suggests that the kinked helix represents a stable conformation of the molecule. Molecular dynamics simulations in the presence of a simple bilayer model and a transbilayer voltage difference are used to explore possible mechanisms of helix insertion. The bilayer is represented by a hydrophobicity potential. An alamethicin helix inserts spontaneously in the absence of a transbilayer voltage. Application of a cis positive voltage decreases the time to insertion. The helix kink angle fluctuates during the simulations. Insertion of the helix is associated with a decrease in the mean kink angle, thus helping the alamethicin molecule to span the bilayer. The simulation results are discussed in terms of models of alamethicin channel gating.