Structure and dynamics of the pore-lining helix of the nicotinic receptor: MD simulations in water, lipid bilayers, and transbilayer bundles

Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.     

Modeling of the pore domain of the GLUR1 channel: homology with K+ channel and binding of channel blockers

Molecular models of the M2 segments of the GluR1 channel have been elaborated using a molecular mechanics approach. The models are based on the homology between pore-lining segments of AMPA receptor channels and the KcsA K+ channel and on cyclic H bonds at the Q/R site of the AMPA receptor channel. The N-terminal region of an M2 segment of the channel is assumed, like that of the K+ channel, to adopt a helical conformation. Due to a deletion, the C-terminal end of the M2 segment of the AMPA receptor is more stretched than that of the K+ channel. As a result, only a single oxygen ring may be exposed to the AMPA receptor channel pore. Data on the block of AMPA receptor channels by dicationic adamantane derivatives have been used to select the most relevant model. The model with the oxygen of a Gly residue (position +2 from the Q/R site) exposed to the pore best fits the experimental data. This model also fits experimental data for another class of AMPA receptor antagonists, the polyamine amides. According to the model, the side-chains of the C-terminal residues are involved in intra-receptor interactions that stabilize the structure of the channel rather than in interactions with ions in the pore.     

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.     

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 structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane alpha-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12 degrees relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 A at its narrowest, to 8.6 A at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

The predicted TM10 transmembrane sequence of the cardiac Ca2+ release channel (ryanodine receptor) is crucial for channel activation and gating

The predicted TM10 transmembrane sequence, (4844)IIFDITFFFFVIVILLAIIQGLII(4867), has been proposed to be the pore inner helix of the ryanodine receptor (RyR) and to play a crucial role in channel activation and gating, as with the inner helix of bacterial potassium channels. However, experimental evidence for the involvement of the TM10 sequence in RyR channel activation and gating is lacking. In the present study, we have systematically investigated the effects of mutations of each residue within the 24-amino acid TM10 sequence of the mouse cardiac ryanodine receptor (RyR2) on channel activation by caffeine and Ca(2+). Intracellular Ca(2+) release measurements in human embryonic kidney 293 cells expressing the RyR2 wild type and TM10 mutants revealed that several mutations in the TM10 sequence either abolished caffeine response or markedly reduced the sensitivity of the RyR2 channel to activation by caffeine. By assessing the Ca(2+) dependence of [(3)H]ryanodine binding to RyR2 wild type and TM10 mutants we also found that mutations in the TM10 sequence altered the sensitivity of the channel to activation by Ca(2+) and enhanced the basal activity of [(3)H]ryanodine binding. Furthermore, single I4862A mutant channels exhibited considerable channel openings and altered gating at very low concentrations of Ca(2+). Our data indicate that the TM10 sequence constitutes an essential determinant for channel activation and gating, in keeping with the proposed role of TM10 as an inner helix of RyR. Our results also shed insight into the orientation of the TM10 helix within the RyR channel pore.     

Functional and Modeling Studies of the Transmembrane Region of the TRPM8 Channel

Members of the transient receptor potential (TRP) ion channel family act as polymodal cellular sensors, which aid in regulating Ca²⁺ homeostasis. Within the TRP family, TRPM8 is the cold receptor that forms a nonselective homotetrameric cation channel. In the absence of TRPM8 crystal structure, little is known about the relationship between structure and function. Inferences of TRPM8 structure have come from mutagenesis experiments coupled to electrophysiology, mainly regarding the fourth transmembrane helix (S4), which constitutes a moderate voltage-sensing domain, and about cold sensor and phosphatidylinositol 4,5-bisphosphate binding sites, which are both located in the C-terminus of TRPM8. In this study, we use a combination of molecular modeling and experimental techniques to examine the structure of the TRPM8 transmembrane and pore helix region including the conducting conformation of the selectivity filter. The model is consistent with a large amount of functional data and was further tested by mutagenesis. We present structural insight into the role of residues involved in intra- and intersubunit interactions and their link with the channel activity, sensitivity to icilin, menthol and cold, and impact on channel oligomerization.     

The structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane α-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12° relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 Å at its narrowest, to 8.6 Å at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

The pore helix dipole has a minor role in inward rectifier channel function

Ion channels lower the energetic barrier for ion passage across cell membranes and enable the generation of bioelectricity. Electrostatic interactions between permeant ions and channel pore helix dipoles have been proposed as a general mechanism for facilitating ion passage. Here, using genetic selections to probe interactions of an exemplar potassium channel blocker, barium, with the inward rectifier Kir2.1, we identify mutants bearing positively charged residues in the potassium channel signature sequence at the pore helix C terminus. We show that these channels are functional, selective, resistant to barium block, and have minimally altered conductance properties. Both the experimental data and model calculations indicate that barium resistance originates from electrostatics. We demonstrate that potassium channel function is remarkably unperturbed when positive charges occur near the permeant ions at a location that should counteract pore helix electrostatic effects. Thus, contrary to accepted models, the pore helix dipole seems to be a minor factor in potassium channel permeation.     

Theoretical studies of the M2 transmembrane segment of the glycine receptor: models of the open pore structure and current-voltage characteristics

The pentameric glycine receptor (GlyR), a member of the nicotinicoid superfamily of ligand-gated ion channels, is an inhibitory Cl(-) channel that is gated by glycine. Using recently published NMR data of the second transmembrane segment (M2) of the human alpha1 GlyR, structural models of pentameric assemblies embedded in a lipid bilayer were constructed using a combination of experimentally determined constraints coupled with all-atom energy minimization. Based on this structure of the pentameric M2 "pore", Brownian dynamics simulations of ion permeation through this putative conducting open state of the channel were carried out. Simulated I-V curves were in good agreement with published experimental current-voltage curves and the anion/cation permeability ratio, suggesting that our open-state model may be representative of the conducting channel of the full-length receptor. These studies also predicted regions of chloride occupancy and suggested residues critical to anion permeation. Calculations of the conductance of the cation-selective mutant A251E channel are also consistent with experimental data. In addition, both rotation and untilting of the pore helices of our model were found to be broadly consistent with closing of the channel, albeit at distinct regions that may reflect alternate gates of the receptor.     

A high resolution structure of the putative hinge region in M2 channel-lining segments of the nicotinic acetylcholine receptor

Transmembrane peptide helices play key roles in signal transduction across cell membranes, yet little is known about their high-resolution structure or the role membrane composition plays in their association, structure, dynamics and ultimately their performance. Using magic angle spinning (MAS) homonuclear dipolar recoupling experiments, the backbone structure at positions L10, L11, and A12 of the M2 ion channel peptide was determined in two lipid systems. Their measurements are in agreement with M2 forming transmembrane helices, but the torsion angles vary considerably from common α−helical values. These measurements show remarkable agreement with a previous computational model of M2 peptides forming a pore domain in which their helices are kinked near the central leucine, L11 [R. Sankararamakrishnan, C. Adcock, M.S.P. Sansom, The pore domain of the nicotinic acetylcholine receptor: Molecular modeling, pore dimensions, and electrostatics, Biophys. J. 71 (1996) 1659-1671]. The generation of high resolution data for transmembrane helices is of critical importance in refining structures for membrane protein and developing models of helix packing interactions.     

A high resolution structure of the putative hinge region in M2 channel-lining segments of the nicotinic acetylcholine receptor

Transmembrane peptide helices play key roles in signal transduction across cell membranes, yet little is known about their high-resolution structure or the role membrane composition plays in their association, structure, dynamics and ultimately their performance. Using magic angle spinning (MAS) homonuclear dipolar recoupling experiments, the backbone structure at positions L10, L11, and A12 of the M2 ion channel peptide was determined in two lipid systems. Their measurements are in agreement with M2 forming transmembrane helices, but the torsion angles vary considerably from common alpha-helical values. These measurements show remarkable agreement with a previous computational model of M2 peptides forming a pore domain in which their helices are kinked near the central leucine, L11 [R. Sankararamakrishnan, C. Adcock, M.S.P. Sansom, The pore domain of the nicotinic acetylcholine receptor: Molecular modeling, pore dimensions, and electrostatics, Biophys. J. 71 (1996) 1659-1671]. The generation of high resolution data for transmembrane helices is of critical importance in refining structures for membrane protein and developing models of helix packing interactions.     

Proton inhibition of unitary currents of vanilloid receptors

Protons, which are released during inflammation and injury, regulate many receptors and ion channels involved in pain transduction, including capsaicin channels (transient receptor potential vanilloid receptors 1). Whereas extracellular acidification both sensitizes and directly activates the channel, it also causes concomitant reduction of the unitary current amplitudes. Here, we investigate the mechanisms and molecular basis of this inhibitory effect of protons on channel conductance. Single-channel recordings showed that the unitary current amplitudes decreased with extracellular pH in a dose-dependent manner, consistent with a model in which protons bind to a site within the channel with an apparent pKa of ∼6. The inhibition was voltage dependent, ∼65% at −60 mV and 37% at +60 mV when pH was reduced from 7.4 to 5.5. The unitary current amplitudes reached saturation at [K⁺] ≥ 1 M, and notably the maximum amplitudes did not converge with different pHs, inconsistent with a blockade model based on surface charge screening or competitive inhibition of permeating ions. Mutagenesis experiments uncovered two acidic residues critical for proton inhibition, one located at the pore entrance and the other on the pore helix. Based on homology to the KcsA structure, the two acidic residues, along with another basic residue also on the pore helix, could form a triad interacting with each other through extensive hydrogen bonds and electrostatic contacts, suggesting that protons may mediate the interactions between the selectivity filter and pore helix, thereby altering the local structure in the filter region and consequently the conductance of the channel.     

Molecular characterization of the inositol 1,4,5-trisphosphate receptor pore-forming segment

Specific residues in the putative pore helix, selectivity filter, and S6 transmembrane helix of the inositol 1,4,5-trisphosphate receptor were mutated in order to examine their effects on channel function. Mutation of 5 of 8 highly conserved residues in the pore helix/selectivity filter region inactivated the channel (C2533A, G2541A, G2545A, G2546A, and G2547A). Of the remaining three mutants, C2527A and R2543A were partially active and G2549A behaved like wild type receptor. Mutation of a putative glycine hinge residue in the S6 helix (G2586A) or a putative gating residue at the cytosolic end of S6 helix (F2592A) had minimal effects on function, although channel function was inactivated by G2586P and F2592D mutations. The mutagenesis data are interpreted in the context of a structural homology model of the inositol 1,4,5-trisphosphate receptor.     

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.     

Acetylcholine receptor pore permeability studied by molecular dynamics simulation

A dynamic model of the closed-state pore of an acetylcholine receptor (five M2 α-helices stabilized with a (CH₂)₁₀₅ ring) is used to examine the migration of uncharged and charged probe particles equivalent to a hexahydrated sodium ion (van der Waals diameter 7.27 Å) propelled by varied external force along the channel axis. Ion movement through the pore is hindered by steric constraints and electrostatic interactions. The van der Waals gate is formed by helix residues 13′ (A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264), whereas the negatively charged residues in the upper part of the channel are important for ion selectivity.     

Structure modeling of the acetylcholine receptor channel and related ligand gated channels

In this model-building study a model for the pore of the acetylcholine receptor channel is proposed. The pore is formed by five -helices of the M2 segment where three rings of hydrophilic side chains point into the channel lumen. This model is in agreement with most experimental data like photolabeling, drug affinity studies, single channel conductivity measurements and cryo electron microscopy known about this channel.This study predicts a strong coupling of the motion of the ions in the channel to that of the charged and highly hydrophilic amino acid side chains at the channel wall. Due to the negative net-charge in the pore more than a single cation may occupy the pore region. The resulting strong local electric fields make the commonly used constant field approximation obsolete for this type of ion channel.     

Decrypting the sequence of structural events during the gating transition of pentameric ligand-gated ion channels based on an interpolated elastic network model

Despite many experimental and computational studies of the gating transition of pentameric ligand-gated ion channels (pLGICs), the structural basis of how ligand binding couples to channel gating remains unknown. By using a newly developed interpolated elastic network model (iENM), we have attempted to compute a likely transition pathway from the closed- to the open-channel conformation of pLGICs as captured by the crystal structures of two prokaryotic pLGICs. The iENM pathway predicts a sequence of structural events that begins at the ligand-binding loops and is followed by the displacements of two key loops (loop 2 and loop 7) at the interface between the extracellular and transmembrane domain, the tilting/bending of the pore-lining M2 helix, and subsequent movements of M4, M3 and M1 helices in the transmembrane domain. The predicted order of structural events is in broad agreement with the Φ-value analysis of α subunit of nicotinic acetylcholine receptor mutants, which supports a conserved core mechanism for ligand-gated channel opening in pLGICs. Further perturbation analysis has supported the critical role of certain intra-subunit and inter-subunit interactions in dictating the above sequence of events.     

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.     

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.