Two models of the influenza A M2 channel domain: verification by comparison

BACKGROUND

The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. It is representative of a very large family of single-transmembrane helix proteins. The functional protein is a tetramer, with the four transmembrane helices forming a proton-permeable channel across the bilayer. Two independently derived models of the M2 channel domain are compared, in order to assess the success of applying molecular modelling approaches to simple membrane proteins.

RESULTS

The Calpha RSMD between the two models is 1.7 A. Both models are composed of a left-handed bundle of helices, with the helices tilted roughly 15 degrees relative to the (presumed) bilayer normal. The two models have similar pore radius profiles, with a pore cavity lined by the Ser31 and Gly34 residues and a pore constriction formed by the ring of His37 residues.

CONCLUSIONS

Independent studies of M2 have converged on the same structural model for the channel domain. This model is in agreement with solid state NMR data. In particular, both model and NMR data indicate that the M2 helices are tilted relative to the bilayer normal and form a left-handed bundle. Such convergence suggests that, at least for simple membrane proteins, restraints-directed modelling might yield plausible models worthy of further computational and experimental investigation.

Alamethicin helices in a bilayer and in solution: molecular dynamics simulations

Alamethicin is an alpha-helical channel-forming peptide, which inserts into lipid bilayers in a voltage-dependent, asymmetrical fashion. Nanosecond molecular dynamics simulations have been used to compare alamethicin conformation and dynamics in three different environments: 1) in water; 2) in methanol; and 3) inserted into a lipid (palmitoyl-oleoyl-phosphatidylcholine) bilayer to form a transmembrane helix. In the bilayer and in methanol, there was little change (Calpha RMSD approximately 0.2 nm over 2 ns and 1 ns) from the initial helical conformation of the peptide. In water there were substantial changes (Calpha RMSD approximately 0.4 nm over 1 ns), especially in the C-terminal segment of the peptide, which lost its alpha-helical conformation. In the bilayer and in methanol, the alamethicin molecule underwent hinge-bending motion about its central Gly-X-X-Pro sequence motif. Analysis of H-bonding interactions revealed that the polar C-terminal side chains of alamethicin provided an "anchor" to the bilayer/water interface via formation of multiple H-bonds that persisted throughout the simulation. This explains why the preferred mode of helix insertion into the bilayer is N-terminal, which is believed to underlie the asymmetry of voltage activation of alamethicin channels.

Lipid properties and the orientation of aromatic residues in OmpF, influenza M2, and alamethicin systems: molecular dynamics simulations

Molecular dynamics simulations allow a direct study of the structure and dynamics of membrane proteins and lipids. We describe the behavior of aromatic residues and lipid properties in POPE and POPC bilayer models with the Escherichia coli OmpF trimer, single alamethicin and Influenza M2 helices, 4-helix M2 bundles, and two alamethicin 6-helix channel models. The total simulation time is over 24 ns, of systems containing solvent, protein, and between 104 and 318 lipids. Various types of adjustment between lipids and proteins occur, depending on the size of the protein and the degree of hydrophobic mismatch between lipid and protein. Single helices cause little measurable effect on nearby lipids whereas the 4-helix bundles, 6-helix channel models, and OmpF cause a significant lowering of order parameters in nearby lipid chains, an increased difference between odd and even chain dihedrals in the magnitude of the trans dihedral fractions and dihedral transition rates, and in most cases a decreased gauche population and a decrease in bilayer thickness. An increased tilt of the lipid chains near the proteins can account for most of the observed decrease in order parameters. The orientation of tryptophans and tyrosines on the outside of the proteins is determined by packing at the protein exterior and non-specific hydrogen bonding with lipids and solvent. The tyrosines in the broad bands that delimit the hydrophobic exterior of OmpF show little change in orientation over one nanosecond. Their rings are oriented predominantly perpendicular to the bilayer plane, with the hydroxyl group pointing toward the lipid-water interface. Phenylalanines in OmpF, alamethicin, and Influenza M2 are more mobile and assume a variety of orientations.

Viral ion channels: molecular modeling and simulation

In a number of membrane-bound viruses, ion channels are formed by integral membrane proteins. These channel proteins include M2 from influenza A, NB from influenza B, and, possibly, Vpu from HIV-1. M2 is important in facilitating uncoating of the influenza A viral genome and is the target of amantadine, an anti-influenza drug. The biological roles of NB and Vpu are less certain. In all cases, the protein contains a single transmembrane alpha-helix close to its N-terminus. Channels can be formed by homo-oligomerization of these proteins, yielding bundles of transmembrane helices that span the membrane and surround a central ion-permeable pore. Molecular modeling may be used to integrate and interpret available experimental data concerning the structure of such transmembrane pores. This has proved successful for the M2 channel domain, where two independently derived models are in agreement with one another, and with solid-state nuclear magnetic resonance (NMR) data. Simulations based on channel models may yield insights into possible ion conduction and selectivity mechanisms.

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.

The structure and organization within the membrane of the helices composing the pore-forming domain of Bacillus thuringiensis delta-endotoxin are consistent with an "umbrella-like" structure of the pore

The aim of this study was to elucidate the mechanism of membrane insertion and the structural organization of pores formed by Bacillus thuringiensis delta-endotoxin. We determined the relative affinities for membranes of peptides corresponding to the seven helices that compose the toxin pore-forming domain, their modes of membrane interaction, their structures within membranes, and their orientations relative to the membrane normal. In addition, we used resonance energy transfer measurements of all possible combinatorial pairs of membrane-bound helices to map the network of interactions between helices in their membrane-bound state. The interaction of the helices with the bilayer membrane was also probed by a Monte Carlo simulation protocol to determine lowest-energy orientations. Our results are consistent with a situation in which helices alpha4 and alpha5 insert into the membrane as a helical hairpin in an antiparallel manner, while the other helices lie on the membrane surface like the ribs of an umbrella (the "umbrella model"). Our results also support the suggestion that alpha7 may serve as a binding sensor to initiate the structural rearrangement of the pore-forming domain.

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.

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.