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.     

Molecular determinants for G protein betagamma modulation of ionotropic glycine receptors

The ligand-gated ion channel superfamily plays a critical role in neuronal excitability. The functions of glycine receptor (GlyR) and nicotinic acetylcholine receptor are modulated by G protein betagamma subunits. The molecular determinants for this functional modulation, however, are still unknown. Studying mutant receptors, we identified two basic amino acid motifs within the large intracellular loop of the GlyR alpha(1) subunit that are critical for binding and functional modulation by Gbetagamma. Mutations within these sequences demonstrated that all of the residues detected are important for Gbetagamma modulation, although both motifs are necessary for full binding. Molecular modeling predicts that these sites are alpha-helixes near transmembrane domains 3 and 4, near to the lipid bilayer and highly electropositive. Our results demonstrate for the first time the sites for G protein betagamma subunit modulation on GlyRs and provide a new framework regarding the ligand-gated ion channel superfamily regulation by intracellular signaling.     

The atypical M2 segment of the beta subunit confers picrotoxinin resistance to inhibitory glycine receptor channels.

Purified preparations of the inhibitory glycine receptor (GlyR) contain alpha and beta subunits, which share homologous primary structures and a common transmembrane topology with other members of the ligand-gated ion channel superfamily. Here, a beta subunit-specific antiserum was shown to precipitate the [3H]strychnine binding sites localized on alpha subunits from membrane extracts of both rat spinal cord and mammalian cells co-transfected with alpha and beta cDNAs. Further, inhibition of alpha homo-oligomeric GlyRs by picrotoxinin, a non-competitive blocker of ion flow, was reduced 50- to 200-fold for alpha/beta hetero-oligomeric receptors generated by cotransfection. Site-directed mutagenesis identified residues within the second predicted transmembrane segment (M2) of the beta subunit as major determinants of picrotoxinin resistance. These data implicate the M2 segment in blocker binding to and lining of the GlyR chloride channel.     

Identification of Novel Specific Allosteric Modulators of the Glycine Receptor Using Phage Display

The glycine receptor (GlyR) is a member of the Cys-loop superfamily of ligand-gated ion channels and the major mediator of inhibitory neurotransmission in the spinal cord and brainstem. Many allosteric modulators affect the functioning of members of this superfamily, with some such as benzodiazepines showing great specificity and others such as zinc, alcohols, and volatile anesthetics acting on multiple members. To date, no potent and efficacious allosteric modulator acting specifically at the GlyR has been identified, hindering both experimental characterization of the receptor and development of GlyR-related therapeutics. We used phage display to identify novel peptides that specifically modulate GlyR function. Peptide D12-116 markedly enhanced GlyR currents at low micromolar concentrations but had no effects on the closely related γ-aminobutyric acid type A receptors. This approach can readily be adapted for use with other channels that currently lack specific allosteric modulators.     

The inhibitory glycine receptor: a ligand-gated chloride channel of the central nervous system

The postsynaptic glycine receptor (GlyR) is a major inhibitory chloride channel protein in the central nervous system. The affinity-purified receptor contains polypeptides of 48 kDa, 58 kDa, and 93 kDa. The 48-kDa (alpha) and 58 kDa (beta) subunits span the postsynaptic membrane in a pentameric arrangement to form the anion channel of the receptor. The 93-kDa polypeptide is cytoplasmically localized and may have an anchoring function. Molecular cloning revealed that different structural characteristics are shared by the membrane-spanning subunits of the GlyR and those of other ligand-gated ion channel proteins. Developmental regulation of the GlyR is characterized by alterations in antagonist binding, heterogeneity of alpha subunits, and increased levels of the 93-kDa polypeptide. Glycine receptor function can be reconstituted by expression of cloned alpha subunits in heterologous cell systems. Positive charges found at the presumed mouths of the GlyR channel appear to be important determinants of ion selectivity. These data establish the anion-conducting GlyR as a homolog of other ligand-gated ion channel proteins and suggest that the diversity of these channels originates from divergent evolution of a primordial channel protein early in phylogeny.     

Microsecond Simulations Indicate that Ethanol Binds between Subunits and Could Stabilize an Open-State Model of a Glycine Receptor

Cys-loop receptors constitute a superfamily of ion channels gated by ligands such as acetylcholine, serotonin, glycine, and γ-aminobutyric acid. All of these receptors are thought to share structural characteristics, but due to high sequence variation and limited structure availability, our knowledge about allosteric binding sites is still limited. These sites are frequent targets of anesthetic and alcohol molecules, and are of high pharmacological importance. We used molecular simulations to study ethanol binding and equilibrium exchange for the homomeric α1 glycine receptor (GlyRα1), modeled on the structure of the Gloeobacter violaceus pentameric ligand-gated channel. Ethanol has a well-known potentiating effect and can be used in high concentrations. By performing two microsecond-scale simulations of GlyR with/without ethanol, we were able to observe spontaneous binding in cavities and equilibrium ligand exchange. Of interest, it appears that there are ethanol-binding sites both between and within the GlyR transmembrane subunits, with the intersubunit site having the highest occupancy and slowest exchange (∼200 ns). This model site involves several residues that were previously identified via mutations as being crucial for potentiation. Finally, ethanol appears to stabilize the GlyR model built on a presumably open form of the ligand-gated channel. This stabilization could help explain the effects of allosteric ligand binding in Cys-loop receptors.     

Homology modeling and molecular dynamics simulations of the alpha1 glycine receptor reveals different states of the channel

Homology modeling is used to build initial models of the transmembrane domain of the human alpha1 glycine receptor (GlyR) based on the most recently published refined structure of nAChR (PDB ID: 2BG9). Six preliminary GlyR models are constructed using two different approaches. In one approach, five different homopentamers are built by symmetric assembly of alpha1 GlyR subunits using only one of the five unique chains of nAChR as a template. In a second approach, each nAChR subunit serves as a template for an alpha1 GlyR subunit. All six initial GlyR constructs are then embedded into a hydrated POPC lipid bilayer and subjected to molecular dynamics simulation for at least six nanoseconds. Each model is stable throughout the simulation, and the final models fall into three distinct categories. Homopentameric GlyR bundles using a single alpha nAChR subunit as a template appear to be in an open conformation. Under an applied external potential, permeation of Cl(-) ions is observed within several ns in a channel built on an alpha chain. Model channels built on non-alpha chains have a constriction either near the intracellular mouth or more centrally located in the pore domain, both of which may be narrow enough to close the channel and whose locations correspond to putative gates observed in nicotinicoid receptors. The differences between these three general models suggest that channel closure may be effected by either rotation or tangential tilting of TM2.     

Functional reconstitution and characterization of recombinant human alpha 1-glycine receptors

By utilizing a baculoviral expression system described previously (Cascio, M., Schoppa, N. E., Grodzicki, R. L., Sigworth, F. J., and Fox, R. O. (1993) J. Biol. Chem. 268, 22135-22142), functional recombinant homomeric human alpha(1)-glycine receptors (GlyR) were overexpressed in insect cell culture, solubilized, purified, and reconstituted into lipid vesicles via gel filtration. Reconstituted GlyR channels were observed to retain native-like activity in single channel recordings of planar bilayers and in flux assays of small unilamellar vesicles, providing evidence that the recombinant homomeric receptor may be functionally reconstituted. This reconstitution is significant in that it indicates that the overexpressed homomeric receptor is an appropriate substrate for subsequent biophysical characterization aimed at the general elucidation of structure-function. Circular dichroism spectroscopy of reconstituted GlyR indicated a low alpha-helical content and a significant fraction of polyproline structure. The small fraction of observed alpha-helix is insufficient to accommodate the four helical transmembrane domains proposed in models for this receptor. By inference, other members of the homologous ligand-gated channel superfamily, which include the ionotropic gamma-aminobutyric acid, acetylcholine, and serotonin receptors, may also be erroneously modeled, and alternate models should be considered.     

Charged residues at the 2' position of human GABAC rho 1 receptors invert ion selectivity and influence open state probability

The ability of members of the nicotinicoid superfamily of ligand-gated ion channels to selectively conduct anions or cations is critical to their function within the central nervous system. Recent work has demonstrated that residues at the intracellular end of the second transmembrane domain, between the -3' and 2' positions, form the ion selectivity filter of these receptors. In this study, the proline residue at the 2' position (Pro-2') at the intracellular end of the second transmembrane domain of the gamma-aminobutyric acid type C rho 1 subunit was mutated to glutamate (rho 1P2'E) and arginine (rho 1P2'R). Dilution potential experiments indicated that the charge selectivity of the rho 1P2'E receptor channels had been inverted, with the channels now becoming predominantly cation selective, indicating the ability of negatively charged residues at this 2' position to control charge selectivity. The mutation was also seen to have significantly decreased agonist potency and intrinsic efficacy. In contrast, the rho 1P2'R receptor channels were anion-selective but were now found to be constitutively open with high holding currents (inhibited by low gamma-aminobutyric acid doses and the competitive antagonist, 1,2,5,6-tetrahydropyridine-4-yl)methylphosphinic acid alone) and increased agonist activity. Hill coefficients of both mutants were decreased, but competitive antagonist studies indicated that their binding sites were not significantly affected.     

Charge selectivity of the Cys-loop family of ligand-gated ion channels

The determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels have been studied for more than a decade. The investigations have mainly covered homomeric receptors e.g. the nicotinic acetylcholine receptor alpha7, the glycine receptor alpha1 and the serotonin receptor 5-HT(3A). Only recently, the determinants of charge selectivity of heteromeric receptors have been addressed for the GABA(A) receptor alpha2beta3gamma2. For all receptor subtypes, the selectivity determinants have been located to an intracellular linker between transmembrane domains M1 and M2. Two features of the M1-M2 linker appear to control ion selectivity. A central role for charged amino acid residues in selectivity has been almost universally observed. Furthermore, recent studies point to an important role of the size of the narrowest constriction in the pore. In the present review, these determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels will be discussed in detail.     

NMR structures of the second transmembrane domain of the human glycine receptor alpha(1) subunit: model of pore architecture and channel gating

Glycine receptors (GlyR) are the primary inhibitory receptors in the spinal cord and belong to a superfamily of ligand-gated ion channels (LGICs) that are extremely sensitive to low-affinity neurological agents such as general anesthetics and alcohols. The high-resolution pore architecture and the gating mechanism of this superfamily, however, remain unclear. The pore-lining second transmembrane (TM2) segments of the GlyR alpha(1) subunit are unique in that they form functional homopentameric channels with conductance characteristics nearly identical to those of an authentic receptor (Opella, S. J., J. Gesell, A. R. Valente, F. M. Marassi, M. Oblatt-Montal, W. Sun, A. F. Montiel, and M. Montal. 1997. Chemtracts Biochem. Mol. Biol. 10:153-174). Using NMR and circular dichroism (CD), we determined the high-resolution structures of the TM2 segment of human alpha(1) GlyR and an anesthetic-insensitive mutant (S267Y) in dodecyl phosphocholine (DPC) and sodium dodecyl sulfate (SDS) micelles. The NMR structures showed right-handed alpha-helices without kinks. A well-defined hydrophilic path, composed of side chains of G2', T6', T10', Q14', and S18', runs along the helical surfaces at an angle approximately 10-20 degrees relative to the long axis of the helices. The side-chain arrangement of the NMR-derived structures and the energy minimization of a homopentameric TM2 channel in a fully hydrated DMPC membrane using large-scale computation suggest a model of pore architecture in which simultaneous tilting movements of entire TM2 helices by a mere 10 degrees may be sufficient to account for the channel gating. The model also suggests that additional residues accessible from within the pore include L3', T7', T13', and G17'. A similar pore architecture and gating mechanism may apply to other channels in the same superfamily, including GABA(A), nACh, and 5-HT(3) receptors.     

A family of acetylcholine-gated chloride channel subunits in Caenorhabditis elegans

The genome of the nematode Caenorhabditis elegans encodes a surprisingly large and diverse superfamily of genes encoding Cys loop ligand-gated ion channels. Here we report the first cloning, expression, and pharmacological characterization of members of a family of anion-selective acetylcholine receptor subunits. Two subunits, ACC-1 and ACC-2, form homomeric channels for which acetylcholine and arecoline, but not nicotine, are efficient agonists. These channels are blocked by d-tubocurarine but not by alpha-bungarotoxin. We provide evidence that two additional subunits, ACC-3 and ACC-4, interact with ACC-1 and ACC-2. The acetylcholine-binding domain of these channels appears to have diverged substantially from the acetylcholine-binding domain of nicotinic receptors.     

Single channel analysis of conductance and rectification in cation-selective, mutant glycine receptor channels

Members of the ligand-gated ion channel superfamily mediate fast synaptic transmission in the nervous system. In this study, we investigate the molecular determinants and mechanisms of ion permeation and ion charge selectivity in this family of channels by characterizing the single channel conductance and rectification of alpha1 homomeric human glycine receptor channels (GlyRs) containing pore mutations that impart cation selectivity. The A-1'E mutant GlyR and the selectivity double mutant ([SDM], A-1'E, P-2' Delta) GlyR, had mean inward chord conductances (at -60 mV) of 7 pS and mean outward conductances of 11 and 12 pS (60 mV), respectively. This indicates that the mutations have not simply reduced anion permeability, but have replaced the previous anion conductance with a cation one. An additional mutation to neutralize the ring of positive charge at the extracellular mouth of the channel (SDM+R19'A GlyR) made the conductance-voltage relationship linear (14 pS at both 60 and -60 mV). When this external charged ring was made negative (SDM+R19'E GlyR), the inward conductance was further increased (to 22 pS) and now became sensitive to external divalent cations (being 32 pS in their absence). The effects of the mutations to the external ring of charge on conductance and rectification could be fit to a model where only the main external energy barrier height for permeation was changed. Mean outward conductances in the SDM+R19'A and SDM+R19'E GlyRs were increased when internal divalent cations were absent, consistent with the intracellular end of the pore being flanked by fixed negative charges. This supports our hypothesis that the ion charge selectivity mutations have inverted the electrostatic profile of the pore by introducing a negatively charged ring at the putative selectivity filter. These results also further confirm the role of external pore vestibule electrostatics in determining the conductance and rectification properties of the ligand-gated ion channels.     

The nicotinic acetylcholine receptor gene family of the nematode Caenorhabditis elegans: an update on nomenclature

The simple nematode, Caenorhabditis elegans, possesses the most extensive known gene family of nicotinic acetylcholine receptor (nAChR)-like subunits. Whilst all show greatest similarity with nAChR subunits of both invertebrates and vertebrates, phylogenetic analysis suggests that just over half of these (32) may represent other members of the cys-loop ligand-gated ion channel superfamily. We have introduced a novel nomenclature system for these “Orphan” subunits, designating them as lgc genes (ligand-gated ion channels of the cys-loop superfamily), which can also be applied in future to unnamed and uncharacterised members of the cys-loop ligand-gated ion channel superfamily. We present here the resulting updated version of the C. elegans nAChR gene family and related ligand-gated ion channel genes.     

Alternative splicing generates two isoforms of the alpha 2 subunit of the inhibitory glycine receptor.

The inhibitory glycine receptor (GlyR) is a ligand-gated chloride channel protein which displays developmental heterogeneity in the mammalian central nervous system. Here we describe 2 novel cDNA variants of the rat GlyR alpha 2 subunit and demonstrate that alternative splicing generates these 2 isoforms. The deduced protein sequences (alpha 2A and alpha 2B) exhibit 99% identity with the previously characterized human alpha 2 subunit. In situ hybridization revealed expression of both alpha 2A and alpha 2B mRNAs in the prenatal rat brain, suggesting that these variant proteins may have a role in synaptogenesis. Heterologous expression in Xenopus oocytes showed that the more abundantly expressed alpha 2A subunit forms strychnine-sensitive ion channels which resemble human alpha 2 subunit GlyRs in their electrophysiological properties.     

Molecular basis for zinc potentiation at strychnine-sensitive glycine receptors

The divalent cation Zn(2+) is a potent potentiator at the strychnine-sensitive glycine receptor (GlyR). This occurs at nanomolar concentrations, which are the predicted endogenous levels of extracellular neuronal Zn(2+). Using structural modeling and functional mutagenesis, we have identified the molecular basis for the elusive Zn(2+) potentiation site on GlyRs and account for the differential sensitivity of GlyR alpha(1) and GlyR alpha(2) to Zn(2+) potentiation. In addition, juxtaposed to this Zn(2+) site, which is located externally on the N-terminal domain of the alpha subunit, another residue was identified in the nearby Cys loop, a region that is critical for receptor gating in all Cys loop ligand-gated ion channels. This residue acted as a key control element in the allosteric transduction pathway for Zn(2+) potentiation, enabling either potentiation or overt inhibition of receptor activation depending upon the moiety resident at this location. Overall, we propose that Zn(2+) binds to a site on the extracellular outer face of the GlyR alpha subunit and exerts its positive allosteric effect via an interaction with the Cys loop to increase the efficacy of glycine receptor gating.     

Predicting the transmembrane secondary structure of ligand-gated ion channels

Recent mutational analyses of ligand-gated ion channels (LGICs) have demonstrated a plausible site of anesthetic action within their transmembrane domains. Although there is a consensus that the transmembrane domain is formed from four membrane-spanning segments, the secondary structure of these segments is not known. We utilized 10 state-of-the-art bioinformatics techniques to predict the transmembrane topology of the tetrameric regions within six members of the LGIC family that are relevant to anesthetic action. They are the human forms of the GABA alpha 1 receptor, the glycine alpha 1 receptor, the 5HT3 serotonin receptor, the nicotinic AChR alpha 4 and alpha 7 receptors and the Torpedo nAChR alpha 1 receptor. The algorithms utilized were HMMTOP, TMHMM, TMPred, PHDhtm, DAS, TMFinder, SOSUI, TMAP, MEMSAT and TOPPred2. The resulting predictions were superimposed on to a multiple sequence alignment of the six amino acid sequences created using the CLUSTAL W algorithm. There was a clear statistical consensus for the presence of four alpha helices in those regions experimentally thought to span the membrane. The consensus of 10 topology prediction techniques supports the hypothesis that the transmembrane subunits of the LGICs are tetrameric bundles of alpha helices.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Evolutionary Relationship of the Ligand-Gated Ion Channels and the Avermectin-Sensitive,Glutamate-Gated Chloride Channels

Two cDNAs, GluClα and GluClβ, encoding glutamate-gated chloride channel subunits that represent targets of the avermectin class of antiparasitic compounds, have recently been cloned from Caenorhabditis elegans (Cully et al., Nature, 371, 707–711, 1994). Expression studies in Xenopus oocytes showed that GluClα and GluClβ have pharmacological profiles distinct from the glutamate-gated cation channels as well as the γ-aminobutyric acid (GABA)- and glycine-gated chloride channels. Establishing the evolutionary relationship of related proteins can clarify properties and lead to predictions about their structure and function. We have cloned and determined the nucleotide sequence of the GluClα and GluClβ genes. In an attempt to understand the evolutionary relationship of these channels with the members of the ligand-gated ion channel superfamily, we have performed gene structure comparisons and phylogenetic analyses of their nucleotide and predicted amino acid sequences. Gene structure comparisons reveal the presence of several intron positions that are not found in the ligand-gated ion channel superfamily, outlining their distinct evolutionary position. Phylogenetic analyses indicate that GluClα and GluClβ form a monophyletic subbranch in the ligand-gated ion channel superfamily and are related to vertebrate glycine channels/receptors. Glutamate-gated chloride channels, with electrophysiological properties similar to GluClα and GluClβ, have been described in insects and crustaceans, suggesting that the glutamate-gated chloride channel family may be conserved in other invertebrate species. The gene structure and phylogenetic analyses in combination with the distinct pharmacological properties demonstrate that GluClα and GluClβ belong to a discrete ligand-gated ion channel family that may represent genes orthologous to the vertebrate glycine channels. Received: 30 September 1996 / Accepted: 15 November 1996     

Structural and biophysical properties of a synthetic channel-forming peptide: Designing a clinically relevant anion selective pore

The design, synthesis, modeling and in vitro testing of channel-forming peptides derived from the cys-loop superfamily of ligand-gated ion channels are part of an ongoing research focus. Over 300 different sequences have been prepared based on the M2 transmembrane segment of the spinal cord glycine receptor α-subunit. A number of these sequences are water-soluble monomers that readily insert into biological membranes where they undergo supramolecular assembly, yielding channels with a range of selectivities and conductances. Selection of a sequence for further modifications to yield an optimal lead compound came down to a few key biophysical properties: low solution concentrations that yield channel activity, greater ensemble conductance, and enhanced ion selectivity. The sequence NK(4)-M2GlyR T19R, S22W (KKKKPARVGLGITTVLTMRTQW) addressed these criteria. The structure of this peptide has been analyzed by solution NMR as a monomer in detergent micelles, simulated as five-helix bundles in a membrane environment, modified by cysteine-scanning and studied for insertion efficiency in liposomes of selected lipid compositions. Taken together, these results define the structural and key biophysical properties of this sequence in a membrane. This model provides an initial scaffold from which rational substitutions can be proposed and tested to modulate anion selectivity. This article is part of a Special Issue entitled: Protein Folding in Membranes.