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

Inter- and Intra-Subunit Butanol/Isoflurane Sites of Action in the Human Glycine Receptor

Glycine receptors (GlyRs) mediate inhibitory neurotransmission and are targets for alcohols and anesthetics in brain. GlyR transmembrane (TM) domains contain critical residues for alcohol/anesthetic action: amino acid A288 in TM3 forms crosslinks with TM1 (I229) in the adjacent subunit as well as TM2 (S267) and TM4 (Y406, W407, I409, Y410) in the same subunit. We hypothesized that these residues may participate in intra-subunit and inter-subunit sites of alcohol/anesthetic action. The following double and triple mutants of GLRA1 cDNA (encoding human glycine receptor alpha 1 subunit) were injected into Xenopus laevis oocytes: I229C/A288C, I229C/A288C/C290S, A288C/Y406C, A288C/W407C, A288C/I409C, and A288C/Y410C along with the corresponding single mutants and wild-type GLRA1. Butanol (22 mM) or isoflurane (0.6 mM) potentiation of GlyR-mediated currents before and after application of the cysteine crosslinking agent HgCl₂ (10 μM) was measured using two-electrode voltage clamp electrophysiology. Crosslinking nearly abolished butanol and isoflurane potentiation in the I229C/A288C and I229C/A288C/C290S mutants but had no effect in single mutants or wild-type. Crosslinking also inhibited butanol and isoflurane potentiation in the TM3-4 mutants (A288C/Y406C, A288C/W407C, A288C/I409C, A288C/Y410C) with no effect in single mutants or wild-type. We extracted proteins from oocytes expressing I229C/288C, A288C/Y410C, or wild-type GlyRs, used mass spectrometry to verify their expression and possible inter-subunit dimerization, plus immunoblotting to investigate the biochemical features of proposed crosslinks. Wild-type GlyR subunits measured about 50 kDa; after crosslinking, the dimeric/monomeric 100:50 kDa band ratio was significantly increased in I229C/288C but not A288C/Y410C mutants or wild-type, providing support for TM1-3 inter-subunit and TM3-4 intra-subunit crosslinking. A GlyR homology model based on the GluCl template provides further evidence for a multi-site model for alcohol/anesthetic interaction with human GLRA1.

Assessment of homology templates and an anesthetic binding site within the γ-aminobutyric acid receptor

BACKGROUND

Anesthetics mediate portions of their activity via modulation of the γ-aminobutyric acid receptor (GABAaR). Although its molecular structure remains unknown, significant progress has been made toward understanding its interactions with anesthetics via molecular modeling.

METHODS

The structure of the torpedo acetylcholine receptor (nAChRα), the structures of the α4 and β2 subunits of the human nAChR, the structures of the eukaryotic glutamate-gated chloride channel (GluCl), and the prokaryotic pH-sensing channels, from Gloeobacter violaceus and Erwinia chrysanthemi, were aligned with the SAlign and 3DMA algorithms. A multiple sequence alignment from these structures and those of the GABAaR was performed with ClustalW. The Modeler and Rosetta algorithms independently created three-dimensional constructs of the GABAaR from the GluCl template. The CDocker algorithm docked a congeneric series of propofol derivatives into the binding pocket and scored calculated binding affinities for correlation with known GABAaR potentiation EC50s.

RESULTS

Multiple structure alignments of templates revealed a clear consensus of residue locations relevant to anesthetic effects except for torpedo nAChR. Within the GABAaR models generated from GluCl, the residues notable for modulating anesthetic action within transmembrane segments 1, 2, and 3 converged on the intersubunit interface between α and β subunits. Docking scores of a propofol derivative series into this binding site showed strong linear correlation with GABAaR potentiation EC50.

CONCLUSION

Consensus structural alignment based on homologous templates revealed an intersubunit anesthetic binding cavity within the transmembrane domain of the GABAaR, which showed a correlation of ligand docking scores with experimentally measured GABAaR potentiation.

Inhibition versus potentiation of ligand-gated ion channels can be altered by a single mutation that moves ligands between intra- and intersubunit sites

Pentameric ligand-gated ion channels (pLGICs) are similar in structure but either inhibited or potentiated by alcohols and anesthetics. This dual modulation has previously not been understood, but the determination of X-ray structures of prokaryotic GLIC provides an ideal model system. Here, we show that a single-site mutation at the F14' site in the GLIC transmembrane domain turns desflurane and chloroform from inhibitors to potentiators, and that this is explained by competing allosteric sites. The F14'A mutation opens an intersubunit site lined by N239 (15'), I240 (16'), and Y263. Free energy calculations confirm this site is the preferred binding location for desflurane and chloroform in GLIC F14'A. In contrast, both anesthetics prefer an intrasubunit site in wild-type GLIC. Modulation is therefore the net effect of competitive binding between the intersubunit potentiating site and an intrasubunit inhibitory site. This provides direct evidence for a dual-site model of allosteric regulation of pLGICs.

Evidence for modulatory sites at the lipid-protein interface of the human multidrug transporter P-glycoprotein

The human multidrug transporter P-glycoprotein (Pgp or ABCB1) sets up pharmacological barriers to many clinically important drugs, a therapeutic remedy for which has yet to be formulated. For the rational design of mechanism-based inhibitors (or modulators), it is necessary to map the potential sites for modulator interaction and understand their modes of communication with the other functional domains of Pgp. In this study, combining directed mutagenesis with homology modeling, we provide evidence of two modulator-specific sites at the lipid protein interface of Pgp. Targeting 21 variant positions in the COOH-terminal transmembrane (TM) regions, we find residues M948 (in TM11) and F983, M986, V988, and Q990 (all four in TM12) critically involved in substrate-site modulation by a thioxanthene-based allosteric modulator cis-(Z)-flupentixol. Interestingly, for ATP-site modulation by the same modulator, only two (M948 and Q990) of those four residues appear indispensable, together with two additional residues, T837 and I864 in TM9 and TM10, respectively, suggesting independent modes of communication linking the allosteric site with the substrate binding and ATPase domains. None of the seven residues identified prove to be critical for modulation of the substrate or ATP sites by Pgp modulators that are transported by the pump, such as cyclosporin A or verapamil, indicating their specificity for cis-(Z)-flupentixol. On the other hand, ATP-site modulation by verapamil proves to be highly sensitive to replacement at positions F716 (in TM7) and I765 (in TM8), and to a more moderate extent at I764 and L772 (both in TM8). Homology modeling based on the known crystal structures of the bacterial multidrug transporter SAV1866 and the mouse Pgp homologue maps the identified residues primarily at the lipid-protein interface of Pgp, in two spatially distinct modulator-specific clusters. The two modulatory sites demonstrate negative synergism in influencing ATP hydrolysis, consolidating their spatial distinctness. Because Pgp is known to recruit drug molecules directly from the lipid bilayer, identification of modulatory sites at the lipid-protein interface and at the same time outside the conventional central drug binding cavity is mechanistically revealing.

Structural basis for alcohol modulation of a pentameric ligand-gated ion channel

Despite its long history of use and abuse in human culture, the molecular basis for alcohol action in the brain is poorly understood. The recent determination of the atomic-scale structure of GLIC, a prokaryotic member of the pentameric ligand-gated ion channel (pLGIC) family, provides a unique opportunity to characterize the structural basis for modulation of these channels, many of which are alcohol targets in brain. We observed that GLIC recapitulates bimodal modulation by n-alcohols, similar to some eukaryotic pLGICs: methanol and ethanol weakly potentiated proton-activated currents in GLIC, whereas n-alcohols larger than ethanol inhibited them. Mapping of residues important to alcohol modulation of ionotropic receptors for glycine, γ-aminobutyric acid, and acetylcholine onto GLIC revealed their proximity to transmembrane cavities that may accommodate one or more alcohol molecules. Site-directed mutations in the pore-lining M2 helix allowed the identification of four residues that influence alcohol potentiation, with the direction of their effects reflecting α-helical structure. At one of the potentiation-enhancing residues, decreased side chain volume converted GLIC into a highly ethanol-sensitive channel, comparable to its eukaryotic relatives. Covalent labeling of M2 positions with an alcohol analog, a methanethiosulfonate reagent, further implicated residues at the extracellular end of the helix in alcohol binding. Molecular dynamics simulations elucidated the structural consequences of a potentiation-enhancing mutation and suggested a structural mechanism for alcohol potentiation via interaction with a transmembrane cavity previously termed the "linking tunnel." These results provide a unique structural model for independent potentiating and inhibitory interactions of n-alcohols with a pLGIC family member.

Modeling anesthetic binding sites within the glycine alpha one receptor based on prokaryotic ion channel templates: the problem with TM4

Ligand-gated ion channels (LGICs) significantly modulate anesthetic effects. Their exact molecular structure remains unknown. This has led to ambiguity regarding the proper amino acid alignment within their 3D structure and, in turn, the location of any anesthetic binding sites. Current controversies suggest that such a site could be located in either an intra- or intersubunit locale within the transmembrane domain of the protein. Here, we built a model of the glycine alpha one receptor (GlyRa1) based on the open-state structures of two new high-resolution ion channel templates from the prokaryote, Gloebacter violaceus (GLIC). Sequence scoring suggests reasonable homology between GlyRa1 and GLIC. Three of the residues notable for modulating anesthetic action are on transmembrane segments 1-3 (TM1-3): (ILE229, SER 267, and ALA 288). They line an intersubunit interface, in contrast to previous models. However, residues from the fourth transmembrane domain (TM4) that are known to modulate a variety of anesthetic effects are quite distant from this putative anesthetic binding site. While this model can account for a large proportion of the physicochemical data regarding such proteins, it cannot readily account for the alterations on anesthetic effects that are due to mutations within TM4.

Normal mode gating motions of a ligand-gated ion channel persist in a fully hydrated lipid bilayer model

We have previously used molecular modeling and normal-mode analyses combined with experimental data to visualize a plausible model of a transmembrane ligand-gated ion channel. We also postulated how the gating motion of the channel may be affected by the presence of various ligands, especially anesthetics. As is typical for normal-mode analyses, those studies were performed in vacuo to reduce the computational complexity of the problem. While such calculations constitute an efficient way to model the large scale structural flexibility of transmembrane proteins, they can be criticized for neglecting the effects of an explicit phospholipid bilayer or hydrated environment. Here, we show the successful calculation of normal-mode motions for our model of a glycine α-1 receptor, now suspended in a fully hydrated lipid bilayer. Despite the almost uniform atomic density, the introduction of water and lipid does not grossly distort the overall gating motion. Normal-mode analysis revealed that even a fully immersed glycine α-1 receptor continues to demonstrate an iris-like channel gating motion as a low-frequency, high-amplitude natural harmonic vibration consistent with channel gating. Furthermore, the introduction of periodic boundary conditions allows the examination of simultaneous harmonic vibrations of lipid in synchrony with the protein gating motions that are compatible with reasonable lipid bilayer perturbations. While these perturbations tend to influence the overall protein motion, this work provides continued support for the iris-like motion model that characterizes gating within the family of ligand-gated ion channels.

5-HT(3) receptors: role in disease and target of drugs

Serotonin type 3 (5-HT(3)) receptors are pentameric ion channels belonging to the superfamily of Cys-loop receptors. Receptor activation either leads to fast excitatory responses or modulation of neurotransmitter release depending on their neuronal localisation. 5-HT(3) receptors are known to be expressed in the central nervous system in regions involved in the vomiting reflex, processing of pain, the reward system, cognition and anxiety control. In the periphery they are present on a variety of neurons and immune cells. 5-HT(3) receptors are known to be involved in emesis, pain disorders, drug addiction, psychiatric and GI disorders. Progress in molecular genetics gives direction to personalised medical strategies for treating complex diseases such as psychiatric and functional GI disorders and unravelling individual drug responses in pharmacogenetic approaches. Here we discuss the molecular basis of 5-HT(3) receptor diversity at the DNA and protein level, of which our knowledge has greatly extended in the last decade. We also evaluate their role in health and disease and describe specific case-control studies addressing the involvement of polymorphisms of 5-HT3 subunit genes in complex disorders and responses to drugs. Furthermore, we focus on the actual state of the pharmacological knowledge concerning not only classical 5-HT(3) antagonists--the setrons--but also compounds of various substance classes targeting 5-HT(3) receptors such as anaesthetics, opioids, cannabinoids, steroids, antidepressants and antipsychotics as well as natural compounds derived from plants. This shall point to alternative treatment options modulating the 5-HT(3) receptor system and open new possibilities for drug development in the future.

The Molecular Mechanisms of Anesthetic Action: Updates and Cutting Edge Developments from the Field of Molecular Modeling

For over 160 years, general anesthetics have been given for the relief of pain and suffering. While many theories of anesthetic action have been purported, it has become increasingly apparent that a significant molecular focus of anesthetic action lies within the family of ligand-gated ion channels (LGIC’s). These protein channels have a transmembrane region that is composed of a pentamer of four helix bundles, symmetrically arranged around a central pore for ion passage. While initial and some current models suggest a possible cavity for binding within this four helix bundle, newer calculations postulate that the actual cavity for anesthetic binding may exist between four helix bundles. In either scenario, these cavities have a transmembrane mode of access and may be partially bordered by lipid moieties. Their physicochemical nature is amphiphilic. Anesthetic binding may alter the overall motion of a ligand-gated ion channel by a “foot-in-door” motif, resulting in the higher likelihood of and greater time spent in a specific channel state. The overall gating motion of these channels is consistent with that shown in normal mode analyses carried out both in vacuo as well as in explicitly hydrated lipid bilayer models. Molecular docking and large scale molecular dynamics calculations may now begin to show a more exact mode by which anesthetic molecules actually localize themselves and bind to specific protein sites within LGIC’s, making the design of future improvements to anesthetic ligands a more realizable possibility.

Linking of Glycine Receptor Transmembrane Segments Three and Four Allows Assignment of Intrasubunit-Facing Residues

Glycine receptors (GlyRs) are pentameric ligand-gated ion channels that mediate inhibitory neurotransmission in the brain and spinal cord and are targets of alcohols and anesthetics. The transmembrane (TM) domain of GlyR subunits is composed of four α-helical segments (TM1-4), but there are conflicting data about the orientation of TM3 and TM4 and, therefore, also the proximity of residues (e.g., A288) that are important for alcohol and anesthetic effects. In the present study, we investigated the proximity of A288 in TM3 to residues in TM4 from M404 to K411. We generated eight double mutant GlyRs (A288C/M404C, A288C/F405C, A288C/Y406C, A288C/W407C, A288C/I408C, A288C/I409C, A288C/Y410C, and A288C/K411C), as well as the corresponding single mutants, and expressed them in Xenopus laevis oocytes. To measure glycine responses, we used two-electrode voltage clamp electrophysiology. We built homology models of the GlyR using structures of the nicotinic acetylcholine receptor (nAChR) and a prokaryotic ion channel (Gloeobacter violaceus, GLIC) as templates, and asked which model best fit our experimental data. Application of the cross-linking reagent HgCl(2) in the closed state produced a leftward shift in the glycine concentration-response curves of the A288C/W407C and A288C/Y410C mutants, suggesting they are able to form cross-links. In addition, when HgCl(2) was coapplied with glycine, responses were changed in the A288C/Y406C, A288C/I409C, and A288C/Y410C double mutants, suggesting that agonist-induced rotation of TM4 allows A288C/Y406C and A288C/I409C to cross-link. These results are consistent with a model of GlyR, based on nAChR, in which A288, Y406, W407, I409, and Y410 face into a four-helical bundle.

Molecular targets and mechanisms for ethanol action in glycine receptors

Glycine receptors (GlyRs) are recognized as the primary mediators of neuronal inhibition in the spinal cord, brain stem and higher brain regions known to be sensitive to ethanol. Building evidence supports the notion that ethanol acting on GlyRs causes at least a subset of its behavioral effects and may be involved in modulating ethanol intake. For over two decades, GlyRs have been studied at the molecular level as targets for ethanol action. Despite the advances in understanding the effects of ethanol in vivo and in vitro, the precise molecular sites and mechanisms of action for ethanol in ligand-gated ion channels in general, and in GlyRs specifically, are just now starting to become understood. The present review focuses on advances in our knowledge produced by using molecular biology, pressure antagonism, electrophysiology and molecular modeling strategies over the last two decades to probe, identify and model the initial molecular sites and mechanisms of ethanol action in GlyRs. The molecular targets on the GlyR are covered on a global perspective, which includes the intracellular, transmembrane and extracellular domains. The latter has received increasing attention in recent years. Recent molecular models of the sites of ethanol action in GlyRs and their implications to our understanding of possible mechanism of ethanol action and novel targets for drug development in GlyRs are discussed.

Anesthetic sensitivity of the Gloeobacter violaceus proton-gated ion channel

A prokaryotic member of the gamma-aminobutyric acid type A receptor superfamily (GLIC) was recently cloned from the cyanobacterium Gloeobacter violaceus, its function characterized, and its 3-dimensional x-ray diffraction crystal structure determined. We report its modulation by 9 anesthetics using 2-electrode voltage clamping in Xenopus laevis oocytes. Desflurane, halothane, isoflurane, sevoflurane, and propofol inhibited currents through GLIC at and below concentrations used clinically. Hill numbers averaged 0.3, indicating negative cooperativity or multiple sites or mechanisms of action. A 2-site model fit the data for desflurane and halothane better than a 1-site model. Xenon and etomidate modulated GLIC at or above clinical concentrations, with no cooperativity. Ethanol and nitrous oxide did not modulate GLIC at surgical anesthetic concentrations. These investigations lay the groundwork for further structural and functional studies of anesthetic actions on GLIC.

Anesthetic-like modulation of receptor function by surfactants: a test of the interfacial theory of anesthesia

INTRODUCTION

Inhaled anesthetics are interfacially active, concentrating at interfaces such as the protein/water or bilayer/water interfaces. We tested the hypothesis that interfacial activity was a sufficient condition for anesthetic-like modulation of receptor function by applying surfactants to gamma-aminobutyric acid type A (GABA(A)), glycine, and N-methyl-d-aspartate (NMDA) receptors. We defined anesthetic-like modulation as an increase in currents through native channels that isoflurane and ethanol increased currents through, and a decrease in currents through channels that isoflurane and ethanol decreased currents through. We also tested the null hypothesis that there would be no difference in modulation of channel currents by surfactants in receptors with point mutations that diminished their response to isoflurane and ethanol compared to the native version of these receptors.

METHODS

The effect of seven surfactants with different head group charges (anionic, cationic, zwitterionic, and uncharged) and tail lengths (8 carbons and 12 carbons) on homomeric wild type alpha1 and mutant alpha(1) (S267I) glycine receptors, wild type alpha(1)beta(2)gamma(2s) and mutant alpha(1)(S270I)beta(2)gamma(2s) GABA(A) receptors, and wild type NR1/NR2A and mutant NR1(F639A)/NR2A NMDA receptors was studied. Receptors were expressed in Xenopus laevis oocytes and studied using two-electrode voltage clamping.

RESULTS

All seven surfactants, isoflurane, and ethanol enhanced GABA(A) receptor function. Six of seven surfactants, isoflurane, and ethanol enhanced glycine receptor function. Six of seven surfactants, isoflurane, and ethanol inhibited NMDA receptor function. For the mutant receptors, five of seven surfactants increased currents through GABA(A) receptors, whereas six of seven surfactants increased currents through glycine receptors. Six of seven surfactants decreased currents through the NMDA receptor. In contrast to isoflurane and ethanol, surfactants as a group did not diminish modulation of mutant compared to wild type receptors.

CONCLUSION

These findings identify another large class of compounds (surfactants) that modulate the function of GABA(A), glycine, and NMDA receptors in a manner that is qualitatively similar to inhaled anesthetics. We cannot reject the hypothesis that interfacial activity is a sufficient condition for anesthetic-like modulation of these receptors. Mutations that diminish the modulatory effect of isoflurane and ethanol did not diminish the modulatory effect of the surfactants.

Evidence That Xenon Does Not Produce Open Channel Blockade of the NMDA Receptor

Previous studies had not excluded the possibility that the mechanism by which Xenon (Xe) blocks N-methyl-d-aspartate (NMDA) receptors might be that of an open-channel blocker. We tested this possibility on mutant NMDA receptors carrying an alanine (A) to cysteine (C) mutation located within the SYTANLAAF-motif of the third transmembrane region (TM3). This mutation was shown to yield constitutively open ion channels after modification with a thiol-modifying reagent. We expressed such mutant channels in Neuro2A cells and recorded glutamate (50 μM)-induced currents in the whole cell recording mode. Although Xe (3.5 mM) blocked the currents through the wild-type receptor NR1-1a/NR2A and NR1-1a/NR2B by ∼40% and those through the mutant receptors NR1-1a/NR2A(A650C) or NR1-1a/NR2B(A651C) by ∼30%, it was unable to block the currents through the methane thiosulfonate etyhlammonium-modified mutant receptors. On the other hand, established open-channel blockers of the NMDA receptor such as MK-801 (1 μM) or Mg ions (Mg2+; 1 mM) were able to block these permanently open channels. These results suggest that Xe does not act as a classical open-channel blocker at the NMDA receptor.

The role of multiple hydrogen-bonding groups in specific alcohol binding sites in proteins: insights from structural studies of LUSH

It is now generally accepted that many of the physiological effects of alcohol consumption are a direct result of binding to specific sites in neuronal proteins such as ion channels or other components of neuronal signaling cascades. Binding to these targets generally occurs in water-filled pockets and leads to alterations in protein structure and dynamics. However, the precise interactions required to confer alcohol sensitivity to a particular protein remain undefined. Using information from the previously solved crystal structures of the Drosophila melanogaster protein LUSH in complexes with short-chain alcohols, we have designed and tested the effects of specific amino acid substitutions on alcohol binding. The effects of these substitutions, specifically S52A, T57S, and T57A, were examined using a combination of molecular dynamics, X-ray crystallography, fluorescence spectroscopy, and thermal unfolding. These studies reveal that the binding of ethanol is highly sensitive to small changes in the composition of the alcohol binding site. We find that T57 is the most critical residue for binding alcohols; the T57A substitution completely abolishes binding, while the T57S substitution differentially affects ethanol binding compared to longer-chain alcohols. The additional requirement for a potential hydrogen-bond acceptor at position 52 suggests that both the presence of multiple hydrogen-bonding groups and the identity of the hydrogen-bonding residues are critical for defining an ethanol binding site. These results provide new insights into the detailed chemistry of alcohol's interactions with proteins.

Molecular modeling and mutagenesis reveals a tetradentate binding site for Zn2+ in GABA(A) alphabeta receptors and provides a structural basis for the modulating effect of the gamma subunit

Gamma-aminobutyric acid type A receptors (GABA(A)-R) containing alpha1beta2gamma2 subunits are weakly inhibited by Zn2+, whereas receptors containing only the alpha1beta2 subunits are strongly inhibited. We built homology models of the ion pores of alpha1beta2 and alpha1beta2gamma2 GABA(A)-R using coordinates of the nicotinic acetylcholine receptor as a template. Threading the GABA(A)-R beta2 sequence onto this template placed the 17' histidine and the 20' glutamate residues at adjacent locations in the mouth of the pore, such that a nearly ideal tetradentate site for Zn2+ was formed from two histidine and two glutamate residues between adjacent beta subunits in the alpha1beta2 GABA(A)-R. Following optimization with CHARMM, the distance between the alpha-carbons of the adjacent histidine residues was approximately 9.2 A, close to the ideal distance for a Zn2+ binding site. Loss of inhibition by Zn2+ in alpha1beta2gamma2 GABA(A)-R can be explained by the geometry of these residues in the arrangement alpha1beta2gamma2alpha1beta2, in which the nearest C-alpha-C-alpha distance between the histidine residues is 15.5 A, too far apart for an energetically optimal Zn2+ binding site. We then mutated the gamma subunit at the 17' and/or 20' positions. Zn2+ inhibition was not restored in alpha1beta2gamma2 (I282H) receptors. A novel finding is that the modeling shows the native 20' lysine in gamma2 can compete with Zn2+ for binding to the inserted 17' histidine. Sensitivity to Zn2+ was restored in the double mutant receptor, alpha1beta2gamma2 (I282H; K285E), in which the competition with lysine was removed and a more favorable Zn2+ binding site was formed.

Cross-linking of sites involved with alcohol action between transmembrane segments 1 and 3 of the glycine receptor following activation

The glycine receptor is a member of the Cys-loop, ligand-gated ion channel family and is responsible for inhibition in the CNS. We examined the orientation of amino acids I229 in transmembrane 1 (TM1) and A288 in TM3, which are both critical for alcohol and volatile anesthetic action. We mutated these two amino acids to cysteines either singly or in double mutants and expressed the receptors in Xenopus laevis oocytes. We tested whether disulfide bonds could form between A288C in TM3 paired with M227C, Y228C, I229C, or S231C in TM1. Application of cross-linking (mercuric chloride) or oxidizing (iodine) agents had no significant effect on the glycine response of wild-type receptors or the single mutants. In contrast, the glycine response of the I229C/A288C double mutant was diminished after application of either mercuric chloride or iodine only in the presence of glycine, indicating that channel gating causes I229C and A288C to fluctuate to be within 6 A apart and form a disulfide bond. Molecular modeling was used to thread the glycine receptor sequence onto a nicotinic acetylcholine receptor template, further demonstrating that I229 and A288 are near-neighbors that can cross-link and providing evidence that these residues contribute to a single binding cavity.

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.

Normal-mode analysis of the glycine alpha1 receptor by three separate methods

Predicting collective dynamics and structural changes in biological macromolecules is pivotal toward a better understanding of many biological processes. Limitations due to large system sizes and inaccessible time scales have prompted the development of alternative techniques for the calculation of such motions. In this work, we present the results of a normal-mode analysis technique based on molecular mechanics that enables the calculation of accurate force-field based vibrations of extremely large molecules and compare it with two elastic network approximate models. When applied to the glycine alpha1 receptor, all three normal-mode analysis algorithms demonstrate an "iris-like" gating motion. Such gating motions have implications for understanding the effects of anesthetic and other ligand binding sites and for the means of transducing agonist binding into ion channel opening. Unlike the more approximate methods, molecular mechanics based analyses can also reveal approximate vibrational frequencies. Such analyses may someday allow the use of protein dynamics elucidated via normal-mode calculations as additional endpoints for future drug design.

The Common Chemical Motifs Within Anesthetic Binding Sites

BACKGROUND

It is not yet possible to obtain crystal structures of anesthetic molecules bound to proteins that are plausible neuronal targets; for example, ligand-gated ion channels. However, there are x-ray crystal structures in which anesthetics are complexed with proteins that are not directly related to anesthetic action. Much useful information about anesthetic-protein interactions can be derived from the x-ray crystal structures of halothane-cholesterol oxidase, bromoform-luciferase, halothane-albumin, and dichloroethane-dehalogenase. These structures show anesthetic-protein interactions at the atomic level.

METHODS

We obtained the known coordinate files for bromoform-luciferase, halothane- albumin, dichloroethane-dehalogenase, and halothane-cholesterol oxidase. These were then modified by adding hydrogens, edited into subsets, and underwent a series of restrained molecular mechanics optimizations. Final analysis of anesthetic polarization within the anesthetic binding site occurred via combined molecular mechanics-quantum mechanics calculations.

RESULTS

The anesthetic binding sites within these well-characterized anesthetic-protein complexes possess a set of common characteristics that we refer to as "binding motifs." The common features of these motifs are polar and nonpolar interactions within an amphiphilic binding cavity, including the presence of weak hydrogen bond interactions with amino acids and water molecules. Calculations also demonstrated the polarizing effect of the amphipathic binding sites on what are otherwise considered quite hydrophobic anesthetics. This polarization appears energetically favorable.

CONCLUSIONS

Anesthetic binding to proteins involves amphipathic interactions.

Glycine receptors: recent insights into their structural organization and functional diversity

Strychnine-sensitive glycine receptors (GlyRs) are known to mediate synaptic inhibition in spinal cord, brainstem and other regions of the CNS. During the past 5 years, considerable progress has been made in delineating structural determinants of ligand binding and channel activation in recombinant GlyRs. Furthermore, immunohistochemical and gene inactivation studies have disclosed distinct distributions and functions of differentially expressed GlyR subtypes in retina, hippocampus and the dorsal horn of the spinal cord. Accordingly, GlyRs regulate not only the excitability of motor and sensory neurones, but are also essential for the processing of photoreceptor signals, neuronal development and inflammatory pain sensitization. Hence, these receptors constitute promising targets for the development of clinically useful compounds.

Molecular modelling of the GABAA ion channel protein

The GABAA ion channel protein is central to the mechanism of action of general anaesthetics and thus to the phenomenon of human consciousness. A molecular model of the alpha1beta2gamma2 gamma-aminobutyric acid type-A (GABAA) ligand-gated ion channel protein has been constructed. The cryo-electron microscopy structure of the nicotinic acetylcholine receptor (nAChR) from Torpedo marmorata and the X-ray crystal structure of the acetylcholine binding protein (AChBP) from Lymnaea stagnalis were used as starting templates for comparative modelling. Features of the modelling approach used in the development of this GABAA model include: (1) multiple sequence alignment of members of the Cys-loop superfamily; (2) the design and implementation of a quasi-ab initio loop modelling algorithm; (3) expansion of the transmembrane domain (TMD) ion pore to model the open-state of the GABAA channel; (4) hydrophobicity analysis of the TMD to refine the structure in regions involved in general anaesthetic binding. The final model of the alpha1beta2gamma2 GABAA protein agrees with available experimental data concerning general anaesthetics.

Accessibility to residues in transmembrane segment four of the glycine receptor

Glycine receptors (GlyRs) are members of the ligand-gated ion channel superfamily. Each subunit has four transmembrane segments (TM1-TM4). Several studies suggest that amino acids in all four TMs face into a water-filled, alcohol and anesthetic binding cavity in the extracellular portion of the transmembrane domain. TM4 should contribute a "wall" to this cavity, but the residues involved are unknown. Here, we determined the ability of an alcohol analog, propyl methanethiosulfonate (propyl MTS), to covalently react with twelve GlyR TM4 positions (I401-I412) after mutating the original amino acids to cysteines. Reactivity of a cysteine with propyl MTS implies that the cysteine is exposed to water. W407C, I409C, Y410C, and K411C showed altered receptor function following reaction with propyl MTS in the presence or absence of glycine. The cysteine mutations alone eliminated the effects of ethanol for I409C, Y410C, and K411C, and reduced the effects of octanol for I409C and isoflurane for K411C. The ability of propyl MTS to reduce isoflurane and chloroform potentiation was examined in the reactive mutants. Potentiation by isoflurane was significantly reduced for I409C after reaction. These data demonstrate water-accessibility of specific TM4 positions in the GlyR and suggest involvement of these residues with alcohol and anesthetic action.