Exploring ligand recognition and ion flow in comparative models of the human GABA type A receptor

GABA receptor associated protein changes the electrostatic environment around the GABA type A receptor

We have performed fully atomistic molecular dynamics simulations of the intracellular domain of a model of the GABA_A receptor with and without the GABA receptor associated protein (GABARAP) bound. We have also calculated the electrostatic potential due to the receptor, in the absence and presence of GABARAP. We find that GABARAP binding changes the electrostatic properties around the GABA_A receptor and could lead to increased conductivity of chloride ions through the receptor. We also find that ion motions that would result in conducting currents are observed nearly twice as often when GABARAP binds. These results are consistent with data from electrophysiological experiments.

Functional movements of the GABA type A receptor

We have performed a Monte Carlo simulation of the GABA type A receptor. We have analysed the configurations and developed a correlation tensor method to predict receptor gating.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.

GABAA receptor family: overview on structural characterization

The pentameric γ-aminobutyric acid type A receptors are ion channels activated by ligands, which intervene in the rapid inhibitory transmission in the mammalian CNS. Due to their rich pharmacology and therapeutic potential, it is essential to understand their structure and function thoroughly. This deep characterization was hampered by the lack of experimental structural information for many years. Thus, computational techniques have been extensively combined with experimental data, in order to undertake the study of γ-aminobutyric acid type A receptors and their interaction with drugs. Here, we review the exciting journey made to assess the structures of these receptors and outline major outcomes. Finally, we discuss the brand new structure of the α1β2γ2 subtype and the amazing advances it brings to the field.

Prediction of GABARAP interaction with the GABA type A receptor

We have performed docking simulations on GABARAP interacting with the GABA type A receptor using SwarmDock. We have also used a novel method to study hydration sites on the surface of these two proteins; this method identifies regions around proteins where desolvation is relatively easy, and these are possible locations where proteins can bind each other. There is a high degree of consistency between the predictions of these two methods. Moreover, we have also identified binding sites on GABARAP for other proteins, and listed possible binding sites for as yet unknown proteins on both GABARAP and the GABA type A receptor intracellular domain.

Different Benzodiazepines Bind with Distinct Binding Modes to GABAA Receptors

Benzodiazepines are clinically relevant drugs that bind to GABAA neurotransmitter receptors at the α+/γ2- interfaces and thereby enhance GABA-induced chloride ion flux leading to neuronal hyperpolarization. However, the structural basis of benzodiazepine interactions with their high-affinity site at GABAA receptors is controversially debated in the literature, and in silico studies led to discrepant binding mode hypotheses. In this study, computational docking of diazepam into α+/γ2- homology models suggested that a chiral methyl group, which is known to promote preferred binding to α5-containing GABAA receptors (position 3 of the seven-membered diazepine ring), could possibly provide experimental evidence that supports or contradicts the proposed binding modes. Thus, we investigated three pairs of R and S isomers of structurally different chemotypes, namely, diazepam, imidazobenzodiazepine, and triazolam derivatives. We used radioligand displacement studies as well as two-electrode voltage clamp electrophysiology in α1β3γ2-, α2β3γ2-, α3β3γ2-, and α5β3γ2-containing GABAA receptors to determine the ligand binding and functional activity of the three chemotypes. Interestingly, both imidazobenzodiazepine isomers displayed comparable binding affinities, while for the other two chemotypes, a discrepancy in binding affinities of the different isomers was observed. Specifically, the R isomers displayed a loss of binding, whereas the S isomers remained active. These findings are in accordance with the results of our in silico studies suggesting the usage of a different binding mode of imidazobenzodiazepines compared to those of the other two tested chemotypes. Hence, we conclude that different chemically related benzodiazepine ligands interact via distinct binding modes rather than by using a common binding mode.

An Electrostatic Funnel in the GABA-Binding Pathway

The γ-aminobutyric acid type A receptor (GABA_A-R) is a major inhibitory neuroreceptor that is activated by the binding of GABA. The structure of the GABA_A-R is well characterized, and many of the binding site residues have been identified. However, most of these residues are obscured behind the C-loop that acts as a cover to the binding site. Thus, the mechanism by which the GABA molecule recognizes the binding site, and the pathway it takes to enter the binding site are both unclear. Through the completion and detailed analysis of 100 short, unbiased, independent molecular dynamics simulations, we have investigated this phenomenon of GABA entering the binding site. In each system, GABA was placed quasi-randomly near the binding site of a GABA_A-R homology model, and atomistic simulations were carried out to observe the behavior of the GABA molecules. GABA fully entered the binding site in 19 of the 100 simulations. The pathway taken by these molecules was consistent and non-random; the GABA molecules approach the binding site from below, before passing up behind the C-loop and into the binding site. This binding pathway is driven by long-range electrostatic interactions, whereby the electrostatic field acts as a ‘funnel’ that sweeps the GABA molecules towards the binding site, at which point more specific atomic interactions take over. These findings define a nuanced mechanism whereby the GABA_A-R uses the general zwitterionic features of the GABA molecule to identify a potential ligand some 2 nm away from the binding site.

Neurotransmitters convey signals from one neuron to the next and are indispensable to the functioning of the nervous system. These small molecules bind to receptors to exert their action. One of the most important neurotransmitters is γ-aminobutyric acid (GABA), which binds to its type A receptor to exert an inhibitory influence on the neuron. Many drugs, both medicinal and nefarious, bind to these neuroreceptors and alter the balance of neuronal signals in the brain. There is a fine balance between these drugs eliciting the desired effect, and causing unwanted and sometimes irreversible alterations in neural behavior. To study this critical binding event, we are using computational simulations to observe precisely how the GABA molecule binds to its type A receptor (GABA_A-receptor). One hundred individual simulations were carried out where GABA was placed near the binding site and then allowed to freely bind to the GABA_A-receptor. Binding occurred in 19 of these simulations. Statistical analysis of these binding simulations reveals the consistent pathway taken by GABA molecules to enter the binding site. This improved understanding of the binding event enables development of safer medicinal neuroactive drugs and countermeasures for effects of neuronal chemical trauma.

Role of the ρ1 GABA(C) receptor N-terminus in assembly, trafficking and function

The GABAC receptor and closely related GABAA receptor are members of the pentameric ligand-gated ion channels (pLGICs) superfamily and mediate inhibitory fast synaptic transmission in the nervous system. Each pLGIC subunit comprises an N-terminal extracellular agonist-binding domain followed by a channel domain and a variable intracellular domain. Available structural information shows that the core of the agonist-binding domain is a β sandwich of ten β-strands, which form the agonist-binding pocket at the subunit interface. This β-sandwich is preceded by an N-terminal α-helix in eukaryotic structures but not in prokaryotic structures. The N-terminal α-helix has been shown to be functionally essential in α7 nicotinic acetylcholine receptors. Sequence analysis of GABAC and GABAA receptors predicts an α-helix in a similar position but preceded by 8 to 46 additional residues, of unknown function, which we term the N-terminal extension. To test the functional role of both the N-terminal extension and the putative N-terminal α-helix in the ρ1 GABAC receptor, we created a series of deletions from the N-terminus. The N-terminal extension was not functionally essential, but its removal did reduce both cell surface expression and cooperativity of agonist-gated channel function. Further deletion of the putative N-terminal α-helix abolished receptor function by preventing cell-surface expression. Our results further demonstrate the essential role of the N-terminal α-helix in the assembly and trafficking of eukaryotic pLGICs. They also provide evidence that the N-terminal extension, although not essential, contributes to receptor assembly, trafficking and conformational changes associated with ligand gating.

Relative positioning of classical benzodiazepines to the γ2-subunit of GABAA receptors

GABAA receptors are the major inhibitory neurotransmitter receptors in the brain. Benzodiazepine exert their action via a high affinity-binding site at the α/γ subunit interface on some of these receptors. Diazepam has sedative, hypnotic, anxiolytic, muscle relaxant, and anticonvulsant effects. It acts by potentiating the current evoked by the agonist GABA. Understanding specific interaction of benzodiazepines in the binding pocket of different GABAA receptor isoforms might help to separate these divergent effects. As a first step, we characterized the interaction between diazepam and the major GABAA receptor isoform α1β2γ2. We mutated several amino acid residues on the γ2-subunit assumed to be located near or in the benzodiazepine binding pocket individually to cysteine and studied the interaction with three ligands that are modified with a cysteine-reactive isothiocyanate group (-NCS). When the reactive NCS group is in apposition to the cysteine residue this leads to a covalent reaction. In this way, three amino acid residues, γ2Tyr58, γ2Asn60, and γ2Val190 were located relative to classical benzodiazepines in their binding pocket on GABAA receptors.

Probing the orthosteric binding site of GABAA receptors with heterocyclic GABA carboxylic acid bioisosteres

The ionotropic GABAA receptors (GABAARs) are widely distributed in the central nervous system where they play essential roles in numerous physiological and pathological processes. A high degree of structural heterogeneity of the GABAAR has been revealed and extensive effort has been made to develop selective and potent GABAAR agonists. This review investigates the use of heterocyclic carboxylic acid bioisosteres within the GABAAR area. Several heterocycles including 3-hydroxyisoxazole, 3-hydroxyisoxazoline, 3-hydroxyisothiazole, and the 1- and 3-hydroxypyrazole rings have been employed in order to map the orthosteric binding site. The physicochemical properties of the heterocyclic moieties making them suitable for bioisosteric replacement of the carboxylic acid in the molecule of GABA are discussed. A variety of synthetic strategies for synthesis of the heterocyclic scaffolds are available. Likewise, methods for introduction of substituents into specific positions of the heterocyclic scaffolds facilitate the investigation of different regions in the orthosteric binding pocket in close vicinity of the core scaffolds of muscimol/GABA. The development of structural models, from pharmacophore models to receptor homology models, has provided more insight into the molecular basis for binding. Similar binding modes are proposed for the heterocyclic GABA analogues covered in this review by use of ligand-receptor docking into the receptor homology model and the presented structure-activity relationships. A network of interactions between the analogues and the binding pocket is leaving no room for substituents and underline the limited space in the GABAAR orthosteric binding site when in the agonist conformation.

Approaching the 5-HT₃ receptor heterogeneity by computational studies of the transmembrane and intracellular domains

5-hydroxytryptamine type-3 receptor (5-HT₃), an important target of many neuroactive drugs, is a cation selective transmembrane pentamer whose functional stoichiometries and subunit arrangements are still debated, due to the extreme complexity of the system. The three dimensional structure of the 5-HT₃R subunits has not been solved so far. Moreover, most of the available structural and functional data is related to the extracellular ligand-binding domain, whereas the transmembrane and the intracellular receptor domains are far less characterised, although they are crucial for receptor function. Here, for the first time, 3D homology models of the transmembrane and the intracellular receptor domains of all the known human 5-HT₃ subunits have been built and assembled into homopentameric (5-HT(3A)R, 5-HT(3B)R, 5-HT(3C)R, 5-HT(3D)R and 5-HT(3E)R) and heteropentameric receptors (5-HT(3AB), 5-HT(3AC), 5-HT(3AD) and 5-HT(3AE)), on the basis of the known three-dimensional structures of the nicotinic-acetylcholine receptor and of the ligand gated channel from Erwinia chrysanthemi. The comparative analyses of sequences, modelled structures, and computed electrostatic properties of the single subunits and of the assembled pentamers shed new light both on the stoichiometric composition and on the physicochemical requirements of the functional receptors. In particular, it emerges that a favourable environment for the crossing of the pore at the transmembrane and intracellular C terminus domain levels by Ca²⁺ ions is granted by the maximum presence of two B subunits in the 5-HT₃ pentamer.

The Cys-loop pentameric ligand-gated ion channel receptors: 50 years on

This year, 2011, the Department of Pharmacology at the University of Alberta celebrated its 50th anniversary. This timeframe covers nearly the entire history of Cys-loop pentameric ligand-gated ion channel (pLGIC) research. In this review we consider how major technological advancements affected our current understanding of pLGICs, and highlight the contributions made by members of our department. The individual at the center of our story is Susan Dunn; her passing earlier this year has robbed the Department of Pharmacology and the research community of a most insightful colleague. Her dissection of ligand interactions with the nAChR, together with their interpretation, was the hallmark of her extensive collaborations with Michael Raftery. Here, we highlight some electrophysiological studies from her laboratory over the last few years, using the technique that she introduced to the department in Edmonton, the 2-electrode voltage-clamp of Xenopus oocytes. Finally, we discuss some single-channel studies of the anionic GlyR and GABA(A)R that prefaced the introduction of this technique to her laboratory.

Diazepam-bound GABAA receptor models identify new benzodiazepine binding-site ligands

Benzodiazepines exert their anxiolytic, anticonvulsant, muscle-relaxant and sedative-hypnotic properties by allosterically enhancing the action of GABA at GABA(A) receptors via their benzodiazepine-binding site. Although these drugs have been used clinically since 1960, the molecular basis of this interaction is still not known. By using multiple homology models and an unbiased docking protocol, we identified a binding hypothesis for the diazepam-bound structure of the benzodiazepine site, which was confirmed by experimental evidence. Moreover, two independent virtual screening approaches based on this structure identified known benzodiazepine-site ligands from different structural classes and predicted potential new ligands for this site. Receptor-binding assays and electrophysiological studies on recombinant receptors confirmed these predictions and thus identified new chemotypes for the benzodiazepine-binding site. Our results support the validity of the diazepam-bound structure of the benzodiazepine-binding pocket, demonstrate its suitability for drug discovery and pave the way for structure-based drug design.

Insights into the structure and pharmacology of GABA(A) receptors

GABA is the major inhibitory neurotransmitter in the adult mammalian CNS. The ionotropic GABA type A receptors (GABA(A)Rs) belong to the Cys-loop family of receptors. Each member of the family is a large pentameric protein in which each subunit traverses the cell membrane four times. Within this family, the GABA type A receptors are particularly important for their rich pharmacology as they are targets for a range of therapeutically important drugs, including the benzodiazepines, barbiturates, neuroactive steroids and anesthetics. This review discusses new insights into receptor properties that allow us to begin to relate the structure of an individual receptor to its functional and pharmacological properties.

New insights into the GABA(A) receptor structure and orthosteric ligand binding: receptor modeling guided by experimental data

GABA(A) receptors (GABA(A)Rs) are ligand gated chloride ion channels that mediate overall inhibitory signaling in the CNS. A detailed understanding of their structure is important to gain insights in, e.g., ligand binding and functional properties of this pharmaceutically important target. Homology modeling is a necessary tool in this regard because experimentally determined structures are lacking. Here we present an exhaustive approach for creating a high quality model of the α(1)β(2)γ(2) subtype of the GABA(A)R ligand binding domain, and we demonstrate its usefulness in understanding details of orthosteric ligand binding. The model was constructed by using multiple templates and by incorporation of knowledge from biochemical/pharmacological experiments. It was validated on the basis of objective energy functions, its ability to account for available residue specific information, and its stability in molecular dynamics (MD) compared with that of the two homologous crystal structures. We then combined the model with extensive structure-activity relationships available from two homologous series of orthosteric GABA(A)R antagonists to create a detailed hypothesis for their binding modes. Excellent agreement with key experimental data was found, including the ability of the model to accommodate and explain a previously developed pharmacophore model. A coupling to agonist binding was thereby established and discussed in relation to activation mechanisms. Our results highlight the importance of critical evaluation and optimization of each step in the homology modeling process. The approach taken here can greatly aid in increasing the understanding of GABA(A)Rs and related receptors where structural insight is limited and reliable models are difficult to obtain.

New insights into the molecular mechanisms of general anaesthetics

This paper provides new insights of how general anaesthetic research should be carried out in the future by an analysis of what we know, what we do not know and what we would like to know. I describe previous hypotheses on the mechanism of action of general anaesthetics (GAs) involving membranes and protein receptors. I provide the reasons why the GABA type A receptor, the NMDA receptor and the glycine receptor are strong candidates for the sites of action of GAs. I follow with a review on attempts to provide a mechanism of action, and how future research should be conducted with the help of physical and chemical methods.

Novel structural determinants of single channel conductance and ion selectivity in 5-hydroxytryptamine type 3 and nicotinic acetylcholine receptors

Nicotinic acetylcholine (nACh) and 5-hydroxytryptamine type 3 (5-HT(3)) receptors are cation-selective ion channels of the pentameric ligand-gated ion channel (pLGIC) superfamily. Multiple lines of evidence adduced over the last 30 years indicate that the lining of the channel of such receptors is formed by the alpha-helical second transmembrane (TM2) domain and flanking sequences contributed by each of the five subunits present within the receptor complex. Specific amino acid residues within, and adjacent to, the TM2 domain influence single channel conductance, ion selectivity, and other aspects of receptor function that include gating and desensitization. However, more recent work has revealed important structural determinants of single channel conductance and ion selectivity that are not associated with the TM2 domain. Direct experimental evidence indicates that the intracellular domain of eukaryotic pLGICs, in particular a region of the loop linking TM3 and TM4 termed the membrane-associated (MA) stretch, exerts a strong influence upon ion channel biophysics. Moreover, recent computational approaches, complemented by experimentation, implicate the extracellular domain as an additional important determinant of ion conduction. This brief review describes how our knowledge of ion conduction and selectivity in cation-selective pLGICs has evolved beyond TM2.

Synthesis and anticonvulsant activity of a new 6-alkoxy-[1,2,4]triazolo[4,3-b]pyridazine

A series of 6-alkoxy-[1,2,4]triazolo[4,3-b]pyridazine derivatives were synthesized. In initial screening and quantitative evaluation, compound 2r was among the most active agents, exhibiting in the same time the lowest toxicity. In the anti-maximal electroshock test, it showed median effective dose (ED50) of 17.3 mg/kg and median toxicity dose (TD50) of 380.3 mg/kg, and the protective index (PI) of 22.0, which is much better than PI of the reference drugs. In a subsequent test, compound 2r had median hypnotic dose (HD50) of 746.6 mg/kg, thus demonstrating much better margin of safety compared to reference drugs. Compound 2r also showed oral activity against MES-induced seizures and lower oral neurotoxicity. For explanation of the putative mechanism of action, compound 2r was tested in chemical induced models.

Relative positioning of diazepam in the benzodiazepine-binding-pocket of GABA receptors

GABA(A) receptors are the major inhibitory neurotransmitter receptors in the brain. Some of them are targets of benzodiazepines that are widely used in clinical practice for their sedative/hypnotic, anxiolytic, muscle relaxant and anticonvulsant effects. In order to rationally separate these different drug actions, we need to understand the interaction of such compounds with the benzodiazepine-binding pocket. With this aim, we mutated residues located in the benzodiazepine-binding site individually to cysteine. These mutated receptors were combined with benzodiazepine site ligands carrying a cysteine reactive group in a defined position. Proximal apposition of reaction partners will lead to a covalent reaction. We describe here such proximity-accelerated chemical coupling reactions of alpha(1)S205C and alpha(1)T206C with a diazepam derivative modified at the C-3 position with a reactive isothiocyanate group (-NCS). We also provide new data that identify alpha(1)H101C and alpha(1)N102C as exclusive sites of the reaction of a diazepam derivative where the -Cl atom is replaced by a -NCS group. Based on these observations we propose a relative positioning of diazepam within the benzodiazepine-binding site of alpha(1)beta(2)gamma(2) receptors.

Isomerization of the proline in the M2-M3 linker is not required for activation of the human 5-HT3A receptor

Each subunit of the cation-selective members of the Cys-loop family of ligand-gated ion channels contains a conserved proline residue in the extracellular loop between the second and third transmembrane domains. In the mouse homomeric 5-hydroxytryptamine type 3A (5-HT(3)A) receptor, the effects of substitution of this proline by unnatural amino acids led to the suggestion that trans-cis isomerization of the protein backbone at this position is integral to agonist-induced channel opening [Nature (2005) vol. 438, pp. 248-252]. We explored the generality of this conclusion using natural amino acid mutagenesis of the homologous human 5-HT(3)A receptor. The conserved proline (P303) was substituted by either a histidine or tryprophan and the mutant receptors were expressed in Xenopus oocytes. These mutations did not significantly affect the magnitude of agonist-mediated currents, compromise channel gating by 5-HT or inhibition of 5-HT-induced currents by either picrotoxin or d-tubocurarine. The mutations did, however, result in altered dependence on extracellular Ca(2+) concentration and a 10-fold increase in the rate of receptor desensitization. These results demonstrate an important role for P303 in 5-HT(3)A receptor function but indicate that trans-cis isomerization at this proline is unlikely to be a general mechanism underlying the gating process.

Molecular modeling and simulation studies of ion channel structures, dynamics and mechanisms

Ion channels are integral membrane proteins that enable selected ions to flow passively across membranes. Channel proteins have been the focus of computational approaches to relate their three-dimensional (3D) structure to their physiological function. We describe a number of computational tools to model ion channels. Homology modeling may be used to construct structural models of channels based on available X-ray structures. Electrostatics calculations enable an approximate evaluation of the energy profile of an ion passing through a channel. Molecular dynamics simulations and free-energy calculations provide information on the thermodynamics and kinetics of channel function.

The subunit arrangement and assembly of ionotropic receptors

Ionotropic receptors mediate rapid communication between neurons. These receptors are oligomers and are usually assembled from multiple subunit types. Receptors built from different subunit combinations have distinct functional properties, such as single-channel conductances, rates of desensitization and sensitivities to activators and inactivators; they can also have different intracellular locations. Methods are now available for determining not only the subunit stoichiometry but also the subunit arrangement within ionotropic receptors. This information will inform experiments designed to understand the molecular basis of receptor assembly and function. It will also permit the modelling of potential ligand-binding sites at the interfaces between the subunits and should lead to a more rational approach to drug development.

Atomic force microscopy reveals the stoichiometry and subunit arrangement of the alpha4beta3delta GABA(A) receptor

The GABA(A) receptor is a chloride-selective ligand-gated ion channel of the Cys-loop superfamily. The receptor consists of five subunits arranged pseudosymmetrically around a central pore. The predominant form of the receptor in the brain contains alpha(1)-, beta(2)-, and gamma(2)-subunits in the arrangement alphabetaalphagammabeta, counter-clockwise around the pore. GABA(A) receptors containing delta-instead of gamma-subunits, although a minor component of the total receptor population, have interesting properties, such as an extrasynaptic location, high sensitivity to GABA, and potential association with conditions such as epilepsy. They are therefore attractive targets for drug development. Here we addressed the subunit arrangement within the alpha(4)beta(3)delta form of the receptor. Different epitope tags were engineered onto the three subunits, and complexes between receptors and anti-epitope antibodies were imaged by atomic force microscopy. Determination of the numbers of receptors doubly decorated by each of the three antibodies revealed a subunit stoichiometry of 2alpha:2beta:1delta. The distributions of angles between pairs of antibodies against the alpha- and beta-subunits both had peaks at around 144 degrees , indicating that these pairs of subunits were nonadjacent. Decoration of the receptor with ligands that bind to the extracellular domain (i.e., the lectin concanavalin A and an antibody that recognizes the beta-subunit N-terminal sequence) showed that the receptor preferentially binds to the mica extracellular face down. Given this orientation, the geometry of complexes of receptors with both an antibody against the delta-subunit and Fab fragments against the alpha-subunits indicates a predominant subunit arrangement of alphabetaalphadeltabeta, counter-clockwise around the pore when viewed from the extracellular space.