A mixed helix-beta-sheet model of the transmembrane region of the nicotinic acetylcholine receptor

Roles for loop 2 residues of alpha1 glycine receptors in agonist activation

The present study tested the hypothesis that several residues in Loop 2 of alpha1 glycine receptors (GlyRs) play important roles in mediating the transduction of agonist activation to channel gating. This was accomplished by investigating the effect of cysteine point mutations at positions 50-60 on glycine responses in alpha1GlyRs using two-electrode voltage clamp of Xenopus oocytes. Cysteine substitutions produced position-specific changes in glycine sensitivity that were consistent with a beta-turn structure of Loop 2, with odd-numbered residues in the beta-turn interacting with other agonist-activation elements at the interface between extracellular and transmembrane domains. We also tested the hypothesis that the charge at position 53 is important for agonist activation by measuring the glycine response of wild type (WT) and E53C GlyRs exposed to methanethiosulfonate reagents. As earlier, E53C GlyRs have a significantly higher EC(50) than WT GlyRs. Exposing E53C GlyRs to the negatively charged 2-sulfonatoethyl methanethiosulfonate, but not neutral 2-hydroxyethyl methanethiosulfonate, positively charged 2-aminoethyl methanethiosulfonate, or 2-trimethylammonioethyl methanethiosulfonate, decreased the glycine EC(50) to resemble WT GlyR responses. Exposure to these reagents did not significantly alter the glycine EC(50) for WT GlyRs. The latter findings suggest that the negative charge at position 53 is important for activation of GlyRs through its interaction with positive charge(s) in other neighboring agonist activation elements. Collectively, the findings provide the basis for a refined molecular model of alpha1GlyRs based on the recent x-ray structure of a prokaryotic pentameric ligand-gated ion channel and offer insight into the structure-function relationships in GlyRs and possibly other ligand-gated ion channels.

Conformational dynamics of the nicotinic acetylcholine receptor channel: a 35-ns molecular dynamics simulation study

The nicotinic acetylcholine receptor (AChR) is the paradigm of ligand-gated ion channels, integral membrane proteins that mediate fast intercellular communication in response to neurotransmitters. A 35-ns molecular dynamics simulation has been performed to explore the conformational dynamics of the entire membrane-spanning region, including the ion channel pore of the AChR. In the simulation, the 20 transmembrane (TM) segments that comprise the whole TM domain of the receptor were inserted into a large dipalmitoylphosphatidylcholine (DPPC) bilayer. The dynamic behavior of individual TM segments and their corresponding AChR subunit helix bundles was examined in order to assess the contribution of each to the conformational transitions of the whole channel. Asymmetrical and asynchronous motions of the M1-M3 TM segments of each subunit were revealed. In addition, the outermost ring of five M4 TM helices was found to convey the effects exerted by the lipid molecules to the central channel domain. Remarkably, a closed-to-open conformational shift was found to occur in one of the channel ring positions in the time scale of the present simulations, the possible physiological significance of which is discussed.

Homology model of the GABAA receptor examined using Brownian dynamics

We have developed a homology model of the GABA(A) receptor, using the subunit combination of alpha1beta2gamma2, the most prevalent type in the mammalian brain. The model is produced in two parts: the membrane-embedded channel domain and the extracellular N-terminal domain. The pentameric transmembrane domain model is built by modeling each subunit by homology with the equivalent subunit of the heteropentameric acetylcholine receptor transmembrane domain. This segment is then joined with the extracellular domain built by homology with the acetylcholine binding protein. The all-atom model forms a wide extracellular vestibule that is connected to an oval chamber near the external surface of the membrane. A narrow, cylindrical transmembrane channel links the outer segment of the pore to a shallow intracellular vestibule. The physiological properties of the model so constructed are examined using electrostatic calculations and Brownian dynamics simulations. A deep energy well of approximately 80 kT accommodates three Cl(-) ions in the narrow transmembrane channel and seven Cl(-) ions in the external vestibule. Inward permeation takes place when one of the ions queued in the external vestibule enters the narrow segment and ejects the innermost ion. The model, when incorporated into Brownian dynamics, reproduces key experimental features, such as the single-channel current-voltage-concentration profiles. Finally, we simulate the gamma2 K289M epilepsy inducing mutation and examine Cl(-) ion permeation through the mutant receptor.

Homology Modeling of a Human Glycine Alpha 1 Receptor Reveals a Plausible Anesthetic Binding Site

The superfamily of ligand-gated ion channels (LGICs) has been implicated in anesthetic and alcohol responses. Mutations within glycine and GABA receptors have demonstrated that possible sites of anesthetic action exist within the transmembrane subunits of these receptors. The exact molecular arrangement of this transmembrane region remains at intermediate resolution with current experimental techniques. Homology modeling methods were therefore combined with experimental data to produce a more exact model of this region. A consensus from multiple bioinformatics techniques predicted the topology within the transmembrane domain of a glycine alpha one receptor (GlyRa1) to be alpha helical. This fold information was combined with sequence information using the SeqFold algorithm to search for modeling templates. Independently, the FoldMiner algorithm was used to search for templates that had structural folds similar to published coordinates of the homologous nAChR (1OED). Both SeqFold and Foldminer identified the same modeling template. The GlyRa1 sequence was aligned with this template using multiple scoring criteria. Refinement of the alignment closed gaps to produce agreement with labeling studies carried out on the homologous receptors of the superfamily. Structural assignment and refinement was achieved using Modeler. The final structure demonstrated a cavity within the core of a four-helix bundle. Residues known to be involved in modulating anesthetic potency converge on and line this cavity. This suggests that the binding sites for volatile anesthetics in the LGICs are the cavities formed within the core of transmembrane four-helix bundles.

Structural basis for lipid modulation of nicotinic acetylcholine receptor function

The nicotinic acetylcholine receptor (AChR) is the archetype molecule in the superfamily of ligand-gated ion channels (LGIC). Members of this superfamily mediate fast intercellular communication in response to endogenous neurotransmitters. This review is focused on the structural and functional crosstalk between the AChR and lipids in the membrane microenvironment, and the modulation exerted by the latter on ligand binding and ion translocation. Experimental approaches using Laurdan extrinsic fluorescence and Förster-type resonance energy transfer (FRET) that led to the characterization of the polarity and molecular dynamics of the liquid-ordered phase AChR-vicinal lipids and the bulk membrane lipids, and the asymmetry of the AChR-rich membrane are reviewed first. The topological relationship between protein and lipid moieties and the changes in physical properties induced by exogenous lipids are discussed next. This background information lays the basis for understanding the occurrence of lipid sites in the AChR transmembrane region, and the selectivity of the protein-lipid interactions. Changes in FRET efficiency induced by fatty acids, phospholipid and cholesterol (Chol), led to the identification of discrete sites for these lipids on the AChR protein, and electron-spin resonance (ESR) spectroscopy has recently facilitated determination of the stoichiometry and selectivity for the AChR of the shell lipid. The influence of lipids on AChR function is discussed next. Combined single-channel and site-directed mutagenesis data fostered the recognition of lipid-sensitive residues in the transmembrane region, dissecting their contribution to ligand binding and channel gating, opening and closing. Experimental evidence supports the notion that the interface between the protein moiety and the adjacent lipid shell is the locus of a variety of pharmacologically relevant processes, including the action of steroids and other lipids.

Molecular modelling of the interactions of carbamazepine and a nicotinic receptor involved in the autosomal dominant nocturnal frontal lobe epilepsy

1. The normal and a mutant (S248F) human neuronal alpha4beta2 nicotinic receptors, and their interaction with the channel blocker carbamazepine (CBZ) have been modelled. The mutant, responsible for the autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE), has an enhanced sensitivity to and a slower recovery from desensitization, a lower conductance, short open times, reduced calcium permeability, and is 3 fold more sensitive to CBZ, a drug used in the treatment of partial epilepsies. 2. Mutant channel properties are explained by the physicochemical properties of the two Phe248 side chains, including size and cation-pi interaction, and their dynamic behaviour. A defective mechanism of dehydration might be responsible for the reduced calcium influx. 3. Phe248 residues are the main component of CBZ binding sites in the mutant, while this is not true for Ser248 in the normal receptor. 4. A higher number of blocking binding sites and a predicted higher affinity found for CBZ in the mutant account for its differential sensitivity to CBZ. 5. Aromatic-aromatic interactions between CBZ and the two Phe248 account for the difference in affinity, which is at least 12 times higher for the mutant, depending on the method used for calculating K(i). 6. Normal vs mutant differences in K(i), enhanced by the higher number of blocking binding sites in the mutant, seem excessive compared to the differential sensitivities to CBZ experimentally found. The negative cooperativity suggested by a predicted overlapping of blocking and non-blocking binding sites gives an explanation, as overlapping is higher in the mutant. 7. For both types of receptors we found that the carbamyl group of the best blocking conformers of CBZ forms hydrogen bonds with serine residues, which may explain the fundamental role of that moiety for this molecule to act as antiepileptic drug.

Structural analysis of synthetic peptide fragments from EmrE, a multidrug resistance protein, in a membrane-mimetic environment

EmrE, a multidrug resistance protein from Escherichia coli, renders the bacterium resistant to a variety of cytotoxic drugs by active translocation out of the cell. The 110-residue sequence of EmrE limits the number of structural possibilities that can be envisioned for this membrane protein. Four helix bundle models have been considered [Yerushalmi, H., Lebendiker, M., and Schuldiner, S. (1996) J. Biol. Chem. 271, 31044-31048]. The validity of EmrE structural models has been probed experimentally by investigations on overlapping peptides (ranging in length from 19 to 27 residues), derived from the sequence of EmrE. The choice of peptides was made to provide sequences of two complete, predicted transmembrane helices (peptides H1 and H3) and two helix-loop-helix motifs (peptides A and B). Peptide (B) also corresponds to a putative hairpin in a speculative beta-barrel model, with the "Pro-Thr-Gly" segment forming a turn. Structure determination in SDS micelles using NMR indicates peptide H1 to be predominantly helical, with helix boundaries in the micellar environment corroborating predicted helical limits. Peptide A adopts a helix-loop-helix structure in SDS micelles, and peptide B was also largely helical in micellar environments. An analogue peptide, C, in which the central "Pro-Thr-Gly" was replaced by "(D)Pro-Gly" displays local turn conformation at the (D)Pro-Gly segment, but neither a continuous helical stretch nor beta-hairpin formation was observed. This study implies that the constraints of membrane and micellar environments largely direct the structure of transmembrane peptides and proteins and study of judiciously selected peptide fragments can prove useful in the structural elucidation of membrane proteins.

Stability of membrane proteins: relevance for the selection of appropriate methods for high-resolution structure determinations

High stability is a prominent characteristic of integral membrane proteins of known atomic structure. But rather than being an intrinsic property, it may be due to a selection exerted by biochemical procedures prior to structure determination, since solubilization results in the transient exposure of membrane proteins to solution conditions. This may cause structural perturbations that interfere with 3D crystallization and hence with X-ray analysis. This problem also affects the preparation of samples for electron crystallography and NMR studies and may account for the fact that high-resolution structures of representatives of whole groups, such as transport proteins and signal transducers, have not been elucidated so far by any method. A knowledge of the proportion of labile proteins among membrane proteins, and of the kinetics of their denaturation, is therefore necessary. Establishing stability profiles, developing methods to maintain lateral pressure, or preventing contact with water (or both) should prove significant in establishing the structures of conformationally flexible proteins.

Structure of the pore-forming transmembrane domain of a ligand-gated ion channel

The structure of the pore-forming transmembrane domain of the nicotinic acetylcholine receptor from Torpedo has been investigated by infrared spectroscopy. Treatment of affinity-purified receptor with either Pronase or proteinase K digests the extramembranous domains (roughly 75% of the protein mass), leaving hydrophobic membrane-imbedded peptides 3-6 kDa in size that are resistant to peptide (1)H/(2)H exchange. Infrared spectra of the transmembrane domain preparations exhibit relatively sharp and symmetric amide I and amide II band contours centered near 1655 and 1545 cm(-)1, respectively, in both (1)H(2)O and (2)H(2)O. The amide I band is very similar to the amide I bands observed in the spectra of alpha-helical proteins, such as myoglobin and bacteriorhodopsin, that lack beta structure and exhibit much less beta-sheet character than is observed in proteins with as little as 20% beta sheet. Curve-fitting estimates 75-80% alpha-helical character, with the remaining peptides likely adopting extended and/or turn structures at the bilayer surface. Infrared dichroism spectra are consistent with transmembrane alpha-helices oriented perpendicular to the bilayer surface. The evidence strongly suggests that the transmembrane domain of the nicotinic receptor, the most intensively studied ligand-gated ion channel, is composed of five bundles of four transmembrane alpha-helices.

Topography of nicotinic acetylcholine receptor membrane-embedded domains

The topography of nicotinic acetylcholine receptor (AChR) membrane-embedded domains and the relative affinity of lipids for these protein regions were studied using fluorescence methods. Intact Torpedo californica AChR protein and transmembrane peptides were derivatized with N-(1-pyrenyl)maleimide (PM), purified, and reconstituted into asolectin liposomes. Fluorescence mapped to proteolytic fragments consistent with PM labeling of cysteine residues in alphaM1, alphaM4, gammaM1, and gammaM4. The topography of the pyrene-labeled Cys residues with respect to the membrane and the apparent affinity for representative lipids were determined by differential fluorescence quenching with spin-labeled derivatives of fatty acids, phosphatidylcholine, and the steroids cholestane and androstane. Different spin label lipid analogs exhibit different selectivity for the whole AChR protein and its transmembrane domains. In all cases labeled residues were found to lie in a shallow position. For M4 segments, this is compatible with a linear alpha-helical structure, but not so for M1, for which "classical" models locate Cys residues at the center of the hydrophobic stretch. The transmembrane topography of M1 can be rationalized on the basis of the presence of a substantial amount of non-helical structure, and/or of kinks attributable to the occurrence of the evolutionarily conserved proline residues. The latter is a striking feature of M1 in the AChR and all members of the rapid ligand-gated ion channel superfamily.

Functional effects of periodic tryptophan substitutions in the alpha M4 transmembrane domain of the Torpedo californica nicotinic acetylcholine receptor

Previous amino acid substitutions at the M4 domain of the Torpedo californica and mouse acetylcholine receptor suggested that the location of the substitution relative to the membrane-lipid interface and perhaps to the ion pore can be critical to the channel gating mechanism [Lasalde, J. A., Tamamizu, S., Butler, D. H., Vibat, C. R. T., Hung, B., and McNamee, M. G. (1996) Biochemistry 35, 14139-14148; Ortiz-Miranda, S. I., Lasalde, J. A., Pappone, P. A., and McNamee, M. G. (1997) J. Membr. Biol. 158, 17-30; Tamamizu, S., Lee, Y. H., Hung, B., McNamee, M. G., and Lasalde-Dominicci, J. A. (1999) J. Membr. Biol. 170, 157-164]. In this study, we introduce tryptophan substitutions at 12 positions (C412W, M415W, L416W, I417W, C418W, I419W, I420W, G421W, T422W, V423W, S424W, and V425W) along this postulated lipid-exposed segment M4 so that we can examine functional consequences on channel gating. The expression levels of mutants C412W, G421W, S424W, and V425W were almost the same as that of the wild type, whereas other mutants (M415W, L416W, C418W, I419W, I420W, T422W, and V423W) had relatively lower expression levels compared to that of the wild type as measured by iodinated alpha-bungarotoxin binding ([(125)I]-alpha-BgTx). Two positions (L416W and I419W) had less than 20% of the wild type expression level. I417W gave no detectable [(125)I]BgTx binding on the surface of oocyte, suggesting that this position might be involved in the AChR assembly, oligomerization, or transport to the cell membrane. The alphaV425W mutant exhibited a significant increase in the open channel probability with a moderate increase in the macroscopic response at higher ACh concentrations very likely due to channel block. The periodicity for the alteration of receptor assembly and ion channel function seems to favor a potential alpha-helical structure. Mutants that have lower levels of expression are clustered on one side of the postulated alpha-helical structure. Mutations that display normal expression and functional activity have been shown previously to face the membrane lipids by independent labeling studies. The functional analysis of these mutations will be presented and discussed in terms of possible structural models.

Chloride channels of glycine and GABA receptors with blockers: Monte Carlo minimization and structure-activity relationships

GABA and glycine receptors (GlyRs) are pentameric ligand-gated ion channels that respond to the inhibitory neurotransmitters by opening a chloride-selective central pore lined with five M2 segments homologous to those of alpha(1) GlyR/ ARVG(2')LGIT(6')TVLTMTTQSSGSR. The activity of cyanotriphenylborate (CTB) and picrotoxinin (PTX), the best-studied blockers of the Cl(-) pores, depends essentially on the subunit composition of the receptors, in particular, on residues in positions 2' and 6' that form the pore-facing rings R(2') and R(6'). Thus, CTB blocks alpha(1) and alpha(1)/beta, but not alpha(2) GlyRs (Rundström, N., V. Schmieden, H. Betz, J. Bormann, and D. Langosch. 1994. Proc. Natl. Acad. Sci. U.S.A. 91:8950-8954). PTX blocks homomeric receptors (alpha(1) GlyR and rat rho(1) GABAR), but weakly antagonizes heteromeric receptors (alpha(1)/beta GlyR and rho(1)/rho(2) GABAR) (Pribilla, I., T. Takagi, D. Langosch, J. Bormann, and H. Betz. 1992. EMBO J. 11:4305-4311; Zhang D., Z. H. Pan, X. Zhang, A. D. Brideau, and S. A. Lipton. 1995. Proc. Natl. Acad. Sci. U.S.A. 92:11756-11760). Using as a template the kinked-helices model of the nicotinic acetylcholine receptor in the open state (Tikhonov, D. B., and B. S. Zhorov. 1998. Biophys. J. 74:242-255), we have built homology models of GlyRs and GABARs and calculated Monte Carlo-minimized energy profiles for the blockers pulled through the pore. The profiles have shallow minima at the wide extracellular half of the pore, a barrier at ring R(6'), and a deep minimum between rings R(6') and R(2') where the blockers interact with five M2s simultaneously. The star-like CTB swings necessarily on its way through ring R(6') and its activity inversely correlates with the barrier at R(6'): Thr(6')s and Ala(2')s in alpha(2) GlyR confine the swinging by increasing the barrier, while Gly(2')s in alpha(1) GlyR and Phe(6')s in beta GlyR shrink the barrier. PTX has an egg-like shape with an isopropenyl group at the elongated end and the rounded end trimmed by ether and carbonyl oxygens. In the optimal binding mode to alpha(1) GlyR and rho(1) GABAR, the rounded end of PTX accepts several H-bonds from Thr(6')s, while the elongated end enters ring R(2'). The lack of H-bond donors on the side chains of Phe(6')s (beta GlyR) and Met(6')s (rho(2) GABAR) deteriorates the binding. The hydrophilic elongated end of picrotin does not fit the hydrophobic ring of Pro(2')s/Ala(2')s in GABARs, but fit a more hydrophilic ring with Gly(2')s in GlyRs. This analysis provides explanations for structure-activity relationships of noncompetitive agonists and predicts a narrow pore of LGICs in agreement with experimental data on the permeation of organic cations.

Effective diffusion coefficients of K+ and Cl- ions in ion channel models

We have used molecular dynamics simulations, corresponding to a total simulation time of 11 ns, to investigate the effective short-time local diffusion coefficient of potassium and chloride ions in a series of model ion channels. These models, which include channels formed by the fungal peptide alamethicin, by a synthetic leucine-serine peptide, and by the pore-lining M2 helix bundle of the nicotinic acetylcholine receptor, have a range of different secondary structures, diameters and hydrophobicities. We find that the diffusion coefficients of both ions are appreciably reduced in the narrower channels, the extent of the reduction being similar for both the anionic and cationic species. This suggests that a difference in mobility cannot be the source of the ion selectivity exhibited by some of the channels (for example, the leucine-serine peptide). We find no evidence for a reduction in mobility of either ion in the nAChR model. These results are broadly in line with a previous similar study of Na+ ions, and may be useful in Poisson-Nernst-Planck, Eyring rate theory or Brownian dynamics calculations of channel conductance.

Improved secondary structure predictions for a nicotinic receptor subunit: incorporation of solvent accessibility and experimental data into a two-dimensional representation

Abstract A refined prediction of the nicotinic acetylcholine receptor (nAChR) subunits' secondary structure was computed with third-generation algorithms. The four selected programs, PHD, Predator, DSC, and NNSSP, based on different prediction approaches, were applied to each sequence of an alignment of nAChR and 5-HT3 receptor subunits, as well as a larger alignment with related subunit sequences from glycine and GABA receptors. A consensus prediction was computed for the nAChR subunits through a "winner takes all" method. By integrating the probabilities obtained with PHD, DSC, and NNSSP, this prediction was filtered in order to eliminate the singletons and to more precisely establish the structure limits (only 4% of the residues were modified). The final consensus secondary structure includes nine alpha-helices (24.2% of the residues, with an average length of 13.9 residues) and 17 beta-strands (22.5% of the residues, with an average length of 6.6 residues). The large extracellular domain is predicted to be mainly composed of beta-strands, with only two helices at the amino-terminal end. The transmembrane segments are predicted to be in a mixed alpha/beta topology (with a predominance of alpha-helices), with no known equivalent in the current protein database. The cytoplasmic domain is predicted to consist of two well-conserved amphipathic helices joined together by an unfolded stretch of variable length and sequence. In general, the segments predicted to occur in a periodic structure correspond to the more conserved regions, as defined by an analysis of sequence conservation per position performed on 152 superfamily members. The solvent accessibility of each residue was predicted from the multiple alignments with PHDacc. Each segment with more than three exposed residues was assumed to be external to the core protein. Overall, these data constitute an envelope of structural constraints. In a subsequent step, experimental data relative to the extracellular portion of the complete receptor were incorporated into the model. This led to a proposed two-dimensional representation of the secondary structure in which the peptide chain of the extracellular domain winds alternatively between the two interfaces of the subunit. Although this representation is not a tertiary structure and does not lead to predictions of specific beta-beta interaction, it should provide a basic framework for further mutagenesis investigations and for fold recognition (threading) searches.

Acetylcholine receptor M3 domain: stereochemical and volume contributions to channel gating

By defining the functional defect in a congenital myasthenic syndrome (CMS), we show that the third transmembrane domain (M3) of the muscle acetylcholine receptor governs the speed and efficiency of gating of its channel. The clinical phenotype of this CMS results from the mutation V285I in M3 of the alpha subunit, which attenuates endplate currents, accelerates their decay and causes abnormally brief acetylcholine-induced single-channel currents. Kinetic analysis of engineered alpha V285I receptors demonstrated a predominant effect on channel gating, with abnormally slow opening and rapid closing rates. Analysis of site-directed mutations revealed stereochemical and volume-dependent contributions of alpha V285 to channel gating. Thus, we demonstrate a functional role for the M3 domain as a key component of the nicotinic acetylcholine receptor channel-gating mechanism.

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.

Secondary structure of the exchange-resistant core from the nicotinic acetylcholine receptor probed directly by infrared spectroscopy and hydrogen/deuterium exchange

The spectral changes that occur in infrared spectra recorded as a function of time after exposure of the nicotinic acetylcholine receptor (nAChR) to 2H2O buffer were examined in order to investigate the secondary structure of the transmembrane domain. The resolution-enhanced amide I band in spectra recorded during the first 12 h after exposure to 2H2O exhibits subtle downshifts in frequency of alpha-helical and beta-sheet vibrations. A strong intensity of the unexchanged alpha-helical vibration near 1655 cm-1 after 3 days exposure to 2H2O suggests that a large proportion of the remaining 25% of unexchanged peptide hydrogens adopts an alpha-helical conformation. Further exposure of the nAChR to 2H2O under conditions of both increasing pH and membrane "fluidity" led to additional exchange of peptide hydrogens for deuterium. The greatest degree of peptide 1H/2H exchange (95%) under nondenaturing conditions was found for the nAChR reconstituted into the highly fluid egg phosphatidylcholine membranes lacking cholesterol and anionic lipids at pH 9.0. This enhanced exchange was accompanied by a decrease in intensity near 1655 cm-1 due to the downshift in frequency of peptides in the alpha-helical conformation, whereas no clear evidence was found for the further exchange of beta-sheet. Some unexchanged alpha-helical peptide hydrogens were still observed. As the exchange-resistant peptides likely include those found within the hydrophobic environment of the lipid bilayer, these data strongly support an alpha-helical secondary structure of the transmembrane domain.

Binding sites for exogenous and endogenous non-competitive inhibitors of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (AChR) is the paradigm of the neurotransmitter-gated ion channel superfamily. The pharmacological behavior of the AChR can be described as three basic processes that progress sequentially. First, the neurotransmitter acetylcholine (ACh) binds the receptor. Next, the intrinsically coupled ion channel opens upon ACh binding with subsequent ion flux activity. Finally, the AChR becomes desensitized, a process where the ion channel becomes closed in the prolonged presence of ACh. The existing equilibrium among these physiologically relevant processes can be perturbed by the pharmacological action of different drugs. In particular, non-competitive inhibitors (NCIs) inhibit the ion flux and enhance the desensitization rate of the AChR. The action of NCIs was studied using several drugs of exogenous origin. These include compounds such as chlorpromazine (CPZ), triphenylmethylphosphonium (TPMP+), the local anesthetics QX-222 and meproadifen, trifluoromethyl-iodophenyldiazirine (TID), phencyclidine (PCP), histrionicotoxin (HTX), quinacrine, and ethidium. In order to understand the mechanism by which NCIs exert their pharmacological properties several laboratories have studied the structural characteristics of their binding sites, including their respective locations on the receptor. One of the main objectives of this review is to discuss all available experimental evidence regarding the specific localization of the binding sites for exogenous NCIs. For example, it is known that the so-called luminal NCIs bind to a series of ring-forming amino acids in the ion channel. Particularly CPZ, TPMP+, QX-222, cembranoids, and PCP bind to the serine, the threonine, and the leucine ring, whereas TID and meproadifen bind to the valine and extracellular rings, respectively. On the other hand, quinacrine and ethidium, termed non-luminal NCIs, bind to sites outside the channel lumen. Specifically, quinacrine binds to a non-annular lipid domain located approximately 7 A from the lipid-water interface and ethidium binds to the vestibule of the AChR in a site located approximately 46 A away from the membrane surface and equidistant from both ACh binding sites. The non-annular lipid domain has been suggested to be located at the intermolecular interfaces of the five AChR subunits and/or at the interstices of the four (M1-M4) transmembrane domains. One of the most important concepts in neurochemistry is that receptor proteins can be modulated by endogenous substances other than their specific agonists. Among membrane-embedded receptors, the AChR is one of the best examples of this behavior. In this regard, the AChR is non-competitively modulated by diverse molecules such as lipids (fatty acids and steroids), the neuropeptide substance P, and the neurotransmitter 5-hydroxytryptamine (5-HT). It is important to take into account that the above mentioned modulation is produced through a direct binding of these endogenous molecules to the AChR. Since this is a physiologically relevant issue, it is useful to elucidate the structural components of the binding site for each endogenous NCI. In this regard, another important aim of this work is to review all available information related to the specific localization of the binding sites for endogenous NCIs. For example, it is known that both neurotransmitters substance P and 5-HT bind to the lumen of the ion channel. Particularly, the locus for substance P is found in the deltaM2 domain, whereas the binding site for 5-HT and related compounds is putatively located on both the serine and the threonine ring. Instead, fatty acid and steroid molecules bind to non-luminal sites. More specifically, fatty acids may bind to the belt surrounding the intramembranous perimeter of the AChR, namely the annular lipid domain, and/or to the high-affinity quinacrine site which is located at a non-annular lipid domain. Additionally, steroids may bind to a site located on the extracellular hydrophi

Mutations at lipid-exposed residues of the acetylcholine receptor affect its gating kinetics

The firmest candidate among the transmembrane portions of the nicotinic acetylcholine receptor (AChR) to be in contact with the lipid bilayer is the fourth segment, M4. To explore the contribution of alphaM4 amino acid residues of mouse AChR to channel gating, we combined site-directed mutagenesis with single-channel recordings. Two residues in alphaM4, Cys418 and Thr422, were found to significantly affect gating kinetics when replaced by alanine. AChRs containing alphaC418A and alphaT422A subunits form channels characterized by a 3- and 5-fold reduction in the mean open time, respectively, suggesting an increase in the closing rate due to the mutations. The calculated changes in the energy barrier for the channel closing process show unequal and coupled contributions of both positions to channel gating. Single-channel recordings of hybrid wild-type alpha/alphaT422A AChR show that the closing rate depends on the number of alpha subunits mutated. Each substitution of threonine to alanine changes the energy barrier of the closing process by approximately 0.5 kcal/mol. Recordings of channels activated by high agonist concentration suggest that these mutations also impair channel opening. Both Cys418 and Thr422 have been postulated to be in contact with the lipid milieu and are highly conserved among species and subunits. Our results support the involvement of lipid-exposed residues in alphaM4 in AChR channel gating mechanism.

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

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

Models and simulations of ion channels and related membrane proteins

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

Topological disposition of Cys 222 in the alpha-subunit of nicotinic acetylcholine receptor analyzed by fluorescence-quenching and electron paramagnetic resonance measurements

The structure of the nicotinic acetylcholine receptor (AChR) has been studied using a combination of fluorescence quenching and electron paramagnetic resonance (EPR) collision gradient methods. The AChR from Torpedo californica was labeled with a fluorescent probe, N-(1-pyrenyl)maleimide, specific for sulfhydryls in a hydrophobic environment, under conditions of selective labeling of Cys222 in the alpha-subunit. alphaCys222 is located in the postulated M1 transmembrane domain and predicted to be at the center of an alpha-helical secondary structure. The spatial disposition of the acetylcholine receptor-bound pyrene with respect to the membrane bilayer was assessed by fluorescence quenching measurements. Quenching of pyrene fluorescence by spin-labeled fatty acids with the doxyl group at positions C-5 and C-12 revealed that the former was more effective, suggesting that the fluorophore is located closer to the membrane-water interface than to the hydrophobic interior. Power saturation EPR spectroscopy was also used to examine the effect of molecular oxygen and water-soluble paramagnetic reagents on the saturation behavior of a nitroxide spin label, which was specifically attached to the same alphaCys222 residue. Using the gradients of these paramagnetic reagents through the membrane-solution interface, the distance for the nitroxide derivative from the membrane-solution interface was measured to be approximately 7 A from the headgroup region of the phospholipid bilayer, in agreement with fluorescence quenching results. These results suggest that the M1 transmembrane domain of the AChR probably forms an irregular structure, a beta-strand, or an alpha-helical structure that may span the membrane in a way different from a linear alpha-helix.

Where does cholesterol act during activation of the nicotinic acetylcholine receptor?

Why agonist-induced activation of the nicotinic acetylcholine receptor (nAcChoR) fails completely in the absence of cholesterol is unknown. Affinity-purified nAcChoRs from Torpedo reconstituted into 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine/1, 2-dioleoyl-sn-glycero-3-phosphate/steroid bilayers at mole ratios of 58:12:30 were used to distinguish between three regions of the membrane where cholesterol might act: the lipid bilayer, the lipid-protein interface, or sites within the protein itself. In the bilayer, the role of fluidity has been ruled out and certain neutral lipids can substitute for cholesterol [C. Sunshine, M.G. McNamee, Biochim. Biophys. Acta 1191 (1994) 59-64]; therefore, we first tested the hypothesis that flip-flop of cholesterol across the membrane is important; a plausible mechanism might be the relief of mechanical bending strain induced by a conformation change that expands the two leaflets of the bilayer asymmetrically. Cholesterol analogs prevented from flipping by charged groups attached to the 3-position's hydroxyl supported channel opening, contrary to this hypothesis. The second hypothesis is that interstitial cholesterol binding sites exist deep within the nAcChoR that must be occupied for channel opening to occur. When cholesterol hemisuccinate was covalently 'tethered' to the glycerol backbone of phosphatidylcholine, channel opening was still supported. Thus, if there are functionally important cholesterol sites, they must be very close to the lipid-protein interface and might be termed periannular.

Kinked-helices model of the nicotinic acetylcholine receptor ion channel and its complexes with blockers: simulation by the Monte Carlo minimization method

A model of the nicotinic acetylcholine receptor ion channel was elaborated based on the data from electron microscopy, affinity labeling, cysteine scanning, mutagenesis studies, and channel blockade. A restrained Monte Carlo minimization method was used for the calculations. Five identical M2 segments (the sequence EKMTLSISVL10LALTVFLLVI20V) were arranged in five-helix bundles with various geometrical profiles of the pore. For each bundle, energy profiles for chlorpromazine, QX-222, pentamethonium, and other blocking drugs pulled through the pore were calculated. An optimal model obtained allows all of the blockers free access to the pore, but retards them at the rings of residues known to contribute to the corresponding binding sites. In this model, M2 helices are necessarily kinked. They come into contact with each other at the cytoplasmic end but diverge at the synaptic end, where N-termini of M1 segments may contribute to the pore. The kinks disengage alpha-helical H-bonds between Ala12 and Ser8. The uncoupled lone electron pairs of Ser8 carbonyl oxygens protrude into the pore, forming a hydrophilic ring that may be important for the permeation of cations. A split network of H-bonds provides a flexibility to the chains Val9-Ala12, the numerous conformations of which form only two or three intrasegment H-bonds. The cross-ectional dimensions of the interface between the flexible chains vary essentially at the level of Leu11. We suggest that conformational transitions in the chains Val9-Ala12 are responsible for the channel gating, whereas rotations of more stable alpha-helical parts of M2 segments may be necessary to transfer the channel in the desensitized state.

Topology of ligand binding sites on the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (AChR) presents two very well differentiated domains for ligand binding that account for different cholinergic properties. In the hydrophilic extracellular region of both alpha subunits there exist the binding sites for agonists such as the neurotransmitter acetylcholine (ACh) and for competitive antagonists such as d-tubocurarine. Agonists trigger the channel opening upon binding while competitive antagonists compete for the former ones and inhibit its pharmacological action. Identification of all residues involved in recognition and binding of agonist and competitive antagonists is a primary objective in order to understand which structural components are related to the physiological function of the AChR. The picture for the localisation of the agonist/competitive antagonist binding sites is now clearer in the light of newer and better experimental evidence. These sites are mainly located on both alpha subunits in a pocket approximately 30-35 A above the surface membrane. Since both alpha subunits are sequentially identical, the observed high and low affinity for agonists on the receptor is conditioned by the interaction of the alpha subunit with the delta or the gamma chain, respectively. This relationship is opposite for curare-related drugs. This molecular interaction takes place probably at the interface formed by the different subunits. The principal component for the agonist/competitive antagonist binding sites involves several aromatic residues, in addition to the cysteine pair at 192-193, in three loops-forming binding domains (loops A-C). Other residues such as the negatively changed aspartates and glutamates (loop D), Thr or Tyr (loop E), and Trp (loop F) from non-alpha subunits were also found to form the complementary component of the agonist/competitive antagonist binding sites. Neurotoxins such as alpha-, kappa-bungarotoxin and several alpha-conotoxins seem to partially overlap with the agonist/competitive antagonist binding sites at multiple point of contacts. The alpha subunits also carry the binding site for certain acetylcholinesterase inhibitors such as eserine and for the neurotransmitter 5-hydroxytryptamine which activate the receptor without interacting with the classical agonist binding sites. The link between specific subunits by means of the binding of ACh molecules might play a pivotal role in the relative shift among receptor subunits. This conformational change would allow for the opening of the intrinsic receptor cation channel transducting the external chemical signal elicited by the agonist into membrane depolarisation. The ion flux activity can be inhibited by non-competitive inhibitors (NCIs). For this kind of drugs, a population of low-affinity binding sites has been found at the lipid-protein interface of the AChR. In addition, several high-affinity binding sites have been found to be located at different rings on the M2 transmembrane domain, namely luminal binding sites. In this regard, the serine ring is the locus for exogenous NCIs such as chlorpromazine, triphenylmethylphosphonium, the local anaesthetic QX-222, phencyclidine, and trifluoromethyliodophenyldiazirine. Trifluoromethyliodophenyldiazirine also binds to the valine ring, which is the postulated site for cembranoids. Additionally, the local anaesthetic meproadifen binding site seems to be located at the outer or extracellular ring. Interestingly, the M2 domain is also the locus for endogenous NCIs such as the neuropeptide substance P and the neurotransmitter 5-hydroxytryptamine. In contrast with this fact, experimental evidence supports the hypothesis for the existence of other NCI high-affinity binding sites located not at the channel lumen but at non-luminal binding domains. (ABSTRACT TRUNCATED)

Assigning functions to residues in the acetylcholine receptor channel region (review)

This review is concerned with the functional domains of the nicotinic acetylcholine receptor (AChR) involved in ion permeation. These comprise the ion pore and its gate. The latter allows the channel to be almost exclusively closed in the absence of agonist and favours ion flux in its presence. Early photoaffinity labelling experiments using open-channel blockers and site-directed mutagenesis studies identified M2 of each AChR subunit as the transmembrane domain lining the walls of the ion pore. Several biochemical, electrophysiological, and mutagenesis studies as well as molecular modelling and in vitro studies of ion channel formation with synthetic peptides corroborate these findings. Point mutations combined with electrophysiological techniques have contributed to dissecting the AChR channel region assigning functions to individual amino acid residues, thus revealing structural and functional stratification of the M2 channel domain. Specific residues have been found to be structural determinants of conductance, ion selectivity, gating, and desensitization. The three-dimensional structure of the AChR protein at 9A resolution suggests a possible arrangement of the M2 alpha-helices in the open and closed states, respectively. In spite of the current wealth of knowledge on the AChR ion channel stemming from the combination of experimental approaches discussed in this review, the mechanistic structure by which the interaction with the agonist favours the opening of the cationic channel remains unknown.

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 high-affinity quinacrine binding site is located at a non-annular lipid domain of the nicotinic acetylcholine receptor

This work deals with the localization of the high-affinity non-competitive quinacrine binding site on the muscle-type nicotinic acetylcholine receptor (AChR). Specifically, quantitative steady-state fluorescence spectroscopy is used to determine whether quinacrine binds to a site located at either the annular or the non-annular lipid domain. For this purpose, we measure the ability of spin-labelled phosphatidylcholine (SL-PC) to quench AChR-bound quinacrine, AChR-bound ethidium and membrane-partitioned 7-(9-anthroyloxy)stearate (7-AS) fluorescence. Additionally, we compare the accessibility of SL-PC which is considered to bind only to the annular lipid domain of the AChR with the accessibility of two non-annular domain-sensing lipids such as 5-doxylstearate (5-SAL) and spin-labelled androstane (ASL). Initial experiments using 7-AS established the experimental conditions for maximum SL-PC membrane partitioning. The non-specific quenching elicited by increasing turbidity of the sample after addition of SL-PC is corrected by means of parallel experiments with unlabelled egg yolk phosphatidylcholine. After correction, the SL-PC quenching experiments show the following order in quenching efficiency: 7-AS > quinacrine >> ethidium. The relative intrinsic sensitivity of quinacrine to TEMPO paramagnetic quenching in acetonitrile is considered to be approximately two times higher than that for 7-AS. Thus, SL-PC was found to be more accessible (about 5-fold) to the membrane-partitioned 7-AS than to the quinacrine locus. In addition, SL-PC was virtually not accessible to the high-affinity non-luminal binding site for ethidium. The relative capacity of SL-PC, 5-SAL, and ASL to quench AChR-bound quinacrine fluorescence indicated that the spin-labelled lipid accessibility to the quinacrine binding site follows the order: 5-SAL > ASL >> SL-PC. Examination of the effect of high concentrations of 5-SAL, of its unlabelled parent stearate, and of SL-PC on ethidium and quinacrine binding showed that: (a) both fatty acids displace quinacrine, but not ethidium, from its high-affinity binding site, however (b) 5-SAL was found to be more effective than stearate to displace quinacrine from its locus, whereas (c) SL-PC competes neither for the ethidium locus nor for the quinacrine binding site. The results suggest that the high-affinity binding site for quinacrine is located at a non-annular lipid domain of the AChR. This particular area has been considered to be located at the intramolecular interfaces of the five AChR subunits and/or at the interstices of the transmembrane domains.

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.