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.

HOLE: a program for the analysis of the pore dimensions of ion channel structural models

A method (HOLE) that allows the analysis of the dimensions of the pore running through a structural model of an ion channel is presented. The algorithm uses a Monte Carlo simulated annealing procedure to find the best route for a sphere with variable radius to squeeze through the channel. Results can be displayed in a graphical fashion or visualized with most common molecular graphical packages. Advances include a method to analyze the anisotropy within a pore. The method can also be used to predict the conductance of channels using a simple empirically corrected ohmic model. As an example the program is applied to the cholera toxin B-subunit pentamer. The compatibility of the crystal structure and conductance data is established.

Molecular dynamics simulations of isolated transmembrane helices of potassium channels

In the middle of the S6 helix in voltage-gated potassium channels there is a highly conserved Pro-Val-Pro motif, while the equivalent M2 helix of inward rectifier potassium channels contains a conserved glycine residue in a comparable position. The structural implications of these conserved motifs are of interest given the evidence that S6 and M2 are components of the lining of their respective pores. Multiple sequence alignment and TM helix prediction methods were used to define consensus regions for S6 and M2. Ensembles of 50 structures for each helix were generated by simulated annealing and restrained molecular dynamics. Time-dependent fluctuations of S6 and M2 were investigated by long time scale molecular dynamics simulations on representative members of each ensemble carried out in vacuo in the presence and absence of a hydrophobic potential that mimics a lipid bilayer. The results are discussed in terms of the structural basis of the kink in S6 and M2 and of a putative functional role for flexible helices as "molecular swivels."

The pore domain of the nicotinic acetylcholine receptor: molecular modeling, pore dimensions, and electrostatics

The pore domain of the nicotinic acetylcholine receptor has been modeled as a bundle of five kinked M2 helices. Models were generated via molecular dynamics simulations incorporating restraints derived from 9-A resolution cryoelectron microscopy data (Unwin, 1993; 1995), and from mutagenesis data that identify channel-lining side chains. Thus, these models conform to current experimental data but will require revision as higher resolution data become available. Models of the open and closed states of a homopentameric alpha 7 pore are compared. The minimum radius of the closed-state model is less than 2 A; the minimum radius of the open-state models is approximately 6 A. It is suggested that the presence of "bound" water molecules within the pore may reduce the effective minimum radii below these values by up to approximately 3 A. Poisson-Boltzmann calculations are used to obtain a first approximation to the potential energy of a monovalent cation as it moves along the pore axis. The differences in electrostatic potential energy profiles between the open-state models of alpha 7 and of a mutant of alpha 7 are consistent with the experimentally observed change in ion selectivity from cationic to anionic. Models of the open state of the heteropentameric Torpedo nicotinic acetylcholine receptor pore domain are also described. Relatively small differences in pore radius and electrostatic potential energy profiles are seen when the Torpedo and alpha 7 models are compared.

Simulation of voltage-dependent interactions of alpha-helical peptides with lipid bilayers

Pore formation in lipid bilayers by channel-forming peptides and toxins is thought to follow voltage-dependent insertion of amphipathic alpha-helices into lipid bilayers. We have developed an approximate potential for use within the CHARMm molecular mechanics program which enables one to simulate voltage-dependent interaction of such helices with a lipid bilayer. Two classes of helical peptides which interact with lipid bilayers have been studied: (a) delta-toxin, a 26 residue channel-forming peptide from Staphylococcus aureus; and (b) synthetic peptides corresponding to the alpha 5 and alpha 7 helices of the pore-forming domain of Bacillus thuringiensis CryIIIA delta-endotoxin. Analysis of delta-toxin molecular dynamics (MD) simulations suggested that the presence of a transbilayer voltage stabilized the inserted location of delta-toxin helices, but did not cause insertion per se. A series of simulations for the alpha 5 and alpha 7 peptides revealed dynamic switching of the alpha 5 helix between a membrane-associated and a membrane-inserted state in response to a transbilayer voltage. In contrast the alpha 7 helix did not exhibit such switching but instead retained a membrane associated state. These results are in agreement with recent experimental studies of the interactions of synthetic alpha 5 and alpha 7 peptides with lipid bilayers.

Structure and orientation of the mammalian antibacterial peptide cecropin P1 within phospholipid membranes

Cecropins are positively charged antibacterial peptides that act by permeating the membrane of susceptible bacteria. To gain insight into the mechanism of membrane permeation, the secondary structure and the orientation within phospholipid membranes of the mammalian cecropin P1 (CecP) was studied using attenuated total reflectance Fourier-transform infrared (ATR-FTIR) spectroscopy and molecular dynamics simulations. The shape and frequency of the amide I and II absorption peaks of CecP within acidic PE/PG multibilayers (phosphatidylethanolamine/phosphatidylglycerol) in a 7:3 (w/w) ratio (a phospholipid composition similar to that of many bacterial membranes), indicated that the peptide is predominantly alpha-helical. Polarized ATR-FTIR spectroscopy was used to determine the orientation of the peptide relative to the bilayer normal of phospholipid multibilayers. The ATR dichroic ratio of the amide I band of CecP peptide reconstituted into oriented PE/PG phospholipid membranes indicated that the peptide is preferentially oriented nearly parallel to the surface of the lipid membranes. A similar secondary structure and orientation were found when zwitterionic phosphatidylcholine phospholipids were used. The incorporation of CecP did not significantly change the order parameters of the acyl chains of the multibilayer, further suggesting that CecP does not penetrate the hydrocarbon core of the membranes. Molecular dynamics simulations were used to gain insight into possible effects of transmembrane potential on the orientation of CecP relative to the membrane. The simulations appear to confirm that CecP adopts an orientation parallel to the membrane surface and does not insert into the bilayer in response to a cis positive transmembrane voltage difference. Taken together, the results further support a "carpet-like" mechanism, rather than the formation of transmembrane pores, as the mode of action of CecP. According to this model, formation of a layer of peptide monomers on the membrane surface destablizes the phospholipid packing of the membrane leading to its eventual disintegration.

Engineering stabilized ion channels: covalent dimers of alamethicin

The peptide alamethicin forms channels with a variety of conductance states. Selective stabilization of a particular state should simplify the task of understanding conductance in terms of channel structure. We synthesized two different covalent dimers of alamethicin in which peptides were linked at their C-terminal ends by flexible tethers. Both dimeric peptides formed channels with conductances that matched those of alamethicin channels. Particular conductance states were selectively stabilized, however, with lifetimes up to 170-fold longer than the same states observed with monomers. In addition, tethering appeared to limit the size of the structures formed so that, even at higher peptide concentrations, a single predominant conductance state was obtained. We suggest this state corresponds to a channel made from six alamethicin molecules (three dimers).

Molecular dynamics simulations of water within models of ion channels

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.

Water in channel-like cavities: structure and dynamics

Ion channels contain narrow columns of water molecules. It is of interest to compare the structure and dynamics of such intrapore water with those of the bulk solvent. Molecular dynamics simulations of modified TIP3P water molecules confined within channel-like cavities have been performed and the orientation and dynamics of the water molecules analyzed. Channels were modeled as cylindrical cavities with lengths ranging from 15 to 60 A and radii from 3 to 12 A. At the end of the molecular dynamics simulations water molecules were observed to be ordered into approximately concentric cylindrical shells. The waters of the outermost shell were oriented such that their dipoles were on average perpendicular to the normal of the wall of the cavity. Water dynamics were analyzed in terms of self-diffusion coefficients and rotational reorientation rates. For cavities of radii 3 and 6 A, water mobility was reduced relative to that of simulated bulk water. For 9- and 12-A radii confined water molecules exhibited mobilities comparable with that of the bulk solvent. If water molecules were confined within an hourglass-shaped cavity (with a central radius of 3 A increasing to 12 A at either end) a gradient of water mobility was observed along the cavity axis. Thus, water within simple models of transbilayer channels exhibits perturbations of structure and dynamics relative to bulk water. In particular the reduction of rotational reorientation rate is expected to alter the local dielectric constant within a transbilayer pore.