Evolutionary trace analysis of ionotropic glutamate receptor sequences and modeling the interactions of agonists with different NMDA receptor subunits

Metazoan evolution and diversity of glutamate receptors and their auxiliary subunits

Glutamate is the major excitatory neurotransmitter in vertebrate and invertebrate nervous systems. Proteins involved in glutamatergic neurotransmission, and chiefly glutamate receptors and their auxiliary subunits, play key roles in nervous system function. Thus, understanding their evolution and uncovering their diversity is essential to comprehend how nervous systems evolved, shaping cognitive function. Comprehensive phylogenetic analysis of these proteins across metazoans have revealed that their evolution is much more complex than what can be anticipated from vertebrate genomes. This is particularly true for ionotropic glutamate receptors (iGluRs), as their current classification into 6 classes (AMPA, Kainate, Delta, NMDA1, NMDA2 and NMDA3) would be largely incomplete. New work proposes a classification of iGluRs into 4 subfamilies that encompass 10 classes. Vertebrate AMPA, Kainate and Delta receptors would belong to one of these subfamilies, named AKDF, the NMDA subunits would constitute another subfamily and non-vertebrate iGluRs would be organised into the previously unreported Epsilon and Lambda subfamilies. Similarly, the animal evolution of metabotropic glutamate receptors has resulted in the formation of four classes of these receptors, instead of the three currently recognised. Here we review our current knowledge on the animal evolution of glutamate receptors and their auxiliary subunits. This article is part of the special issue on 'Glutamate Receptors - Orphan iGluRs'.

Ligand binding domain interface: A tipping point for pharmacological agents binding with GluN1/2A subunit containing NMDA receptors

N-methyl D-aspartate (NMDA) receptors play a crucial role in normal brain function, pathogenesis of neurodegenerative and psychiatric disorders. Functional tetra-heteromeric NMDA receptor contains two obligatory GluN1 subunits and two identical or different non-GluN1 subunits that evolve from six different genes including four GluN2 (A-D) and two GluN3 (A-B) subunits. Since NMDA receptors confer varied physiological properties and spatiotemporal distributions in the brain, pharmacological agents that target NMDA receptors with specific GluN2 subunits have significant potential for therapeutic applications. In the present work, by using electrophysiology techniques, we have studied the role of ligand binding domain (LBD) interactions in determining the effect of well-characterized pharmacological agents including agonists, competitive antagonists, channel blockers and an allosteric modulator. Remarkably, point mutations at the distal end (site-II&III) of GluN1 LBD interface increased memantine potency up to sevenfold when co-expressed with wild type GluN2A receptors but exhibit no effect on Mg²⁺ activity. Conversely, mutations at the proximal end (site-I) of the LBD interface did not affect the memantine but altered Zn²⁺ and Mg²⁺ potency towards opposite directions. These results indicate that GluN1/2A LBD interface interactions play a key role in determining channel function. Further, subtle changes in LBD interaction can be readily translated to the downstream extracellular vestibule of channel pore to adopt a conformation that may affect memantine, Zn²⁺ and Mg²⁺ binding. Further studies on NMDA receptor LBD to transmembrane domain signal propagation mechanisms will help develop GluN2 subunit selective biomolecules that can be used for the treatment of neurological and psychiatric disorders.

Positive Modulatory Interactions of NMDA Receptor GluN1/2B Ligand Binding Domains Attenuate Antagonists Activity

N-methyl D-aspartate receptors (NMDAR) play crucial role in normal brain function and pathogenesis of neurodegenerative and psychiatric disorders. Functional tetra-heteromeric NMDAR contains two obligatory GluN1 subunits and two identical or different non-GluN1 subunits that include six different gene products; four GluN2 (A–D) and two GluN3 (A–B) subunits. The heterogeneity of subunit combination facilities the distinct function of NMDARs. All GluN subunits contain an extracellular N-terminal Domain (NTD) and ligand binding domain (LBD), transmembrane domain (TMD) and an intracellular C-terminal domain (CTD). Interaction between the GluN1 and co-assembling GluN2/3 subunits through the LBD has been proven crucial for defining receptor deactivation mechanisms that are unique for each combination of NMDAR. Modulating the LBD interactions has great therapeutic potential. In the present work, by amino acid point mutations and electrophysiology techniques, we have studied the role of LBD interactions in determining the effect of well-characterized pharmacological agents including agonists, competitive antagonists, and allosteric modulators. The results reveal that agonists (glycine and glutamate) potency was altered based on mutant amino acid sidechain chemistry and/or mutation site. Most antagonists inhibited mutant receptors with higher potency; interestingly, clinically used NMDAR channel blocker memantine was about three-fold more potent on mutated receptors (N521A, N521D, and K531A) than wild type receptors. These results provide novel insights on the clinical pharmacology of memantine, which is used for the treatment of mild to moderate Alzheimer's disease. In addition, these findings demonstrate the central role of LBD interactions that can be exploited to develop novel NMDAR based therapeutics.

Identification of novel allosteric modulator binding sites in NMDA receptors: A molecular modeling study

The dysfunction of N-methyl-d-Aspartate receptors (NMDARs), a subtype of glutamate receptors, is correlated with schizophrenia, stroke, and many other neuropathological disorders. However, not all NMDAR subtypes equally contribute towards these disorders. Since NMDARs composed of different GluN2 subunits (GluN2A-D) confer varied physiological properties and have different distributions in the brain, pharmacological agents that target NMDARs with specific GluN2 subunits have significant potential for therapeutic applications. In our previous research, we have identified a family of novel allosteric modulators that differentially potentiate and/or inhibit NMDARs of differing GluN2 subunit composition. To further elucidate their molecular mechanisms, in the present study, we have identified four potential binding sites for novel allosteric modulators by performing molecular modeling, docking, and in silico mutations. The molecular determinants of the modulator binding sites (MBS), analysis of particular MBS electrostatics, and the specific loss or gain of binding after mutations have revealed modulators that have strong potential affinities for specific MBS on given subunits and the role of key amino acids in either promoting or obstructing modulator binding. These findings will help design higher affinity GluN2 subunit-selective pharmaceuticals, which are currently unavailable to treat psychiatric and neurological disorders.

Cloning and characterization of the N-methyl-D-aspartate receptor subunit NR1 gene from chum salmon, Oncorhynchus keta (Walbaum, 1792)

Here, we report the information about molecular and expression characterization of NR1 gene in chum salmon for the first time. The complete NR1 subunit showed a large open-reading frame of 2844 bp in the total length of 3193 bp, and this cDNA contained a coding region encoding 948 amino acids and a stop codon. The organization of the NR1 subunit of chum salmon were similar of most other fishes, except C’ terminal. The expression of NR1 subunit was to show higher in the natal river near to the hatchery than near to the coast. We expect that the information reported herein may facilitate further investigations on the relationship between memory factors of natal rivers and homing mechanisms in Salmonidae.

Pharmacological insights obtained from structure-function studies of ionotropic glutamate receptors

Ionotropic glutamate receptors mediate the vast majority of fast excitatory synaptic transmission in the CNS. Elucidating the structure of these proteins is central to understanding their overall function and in the last few years a tremendous amount of knowledge has been gained from the crystal structures of the ligand-binding domains of the receptor protein. These efforts have enabled us to unravel the possible mechanisms of partial agonism, agonist selectivity and desensitization. This review summarizes recent data obtained from structural studies of the binding pockets of the GluR2, GluR5/6, NR1 and NR2A subunits and discusses these studies together with homology modelling and molecular dynamics simulations that have suggested possible binding modes for full and partial agonists as well as antagonists within the binding pocket of various ionotropic glutamate receptor subunits. Comparison of the ligand-binding pockets suggests that the ligand-binding mechanisms may be conserved throughout the glutamate receptor family, although agonist selectivity may be explained by a number of features inherent to the AMPA, kainate and NMDA receptor-binding pockets such as steric occlusion, cavity size and the contribution of water-bridged interactions.

Characterisation of the human NMDA receptor subunit NR3A glycine binding site

In this study, we characterise the binding site of the human N-methyl-d-aspartate (NMDA) receptor subunit NR3A. Saturation radioligand binding of the NMDA receptor agonists [(3)H]-glycine and [(3)H]-glutamate showed that only glycine binds to human NR3A (hNR3A) with high affinity (K(d)=535nM (277-793nM)). Eight amino acids, which correspond to amino acids that are critical for ligand binding to other NMDA receptor subunits, situated within the S1S2 predicted ligand binding domain of hNR3A were mutated, which resulted in complete or near complete loss of [(3)H]-glycine binding to hNR3A. The NMDA NR1 glycine site agonist d-serine and partial agonist HA-966 (3-amino-1-hydroxypyrrolid-2-one), similarly to glycine displaced [(3)H]-glycine monophasically, suggesting a single common binding site. However, neither the partial agonist d-cycloserine nor the antagonist 7-chlorokynurenic acid displaced [(3)H]-glycine. Using homology modelling, a model of the NR3A binding pocket was generated which we suggest can be used to identify candidate agonists and antagonists. Our data show that glycine is a ligand, and most probably the endogenous ligand, for native NR3A at a binding site with unique pharmacological characteristics.

A two-entropies analysis to identify functional positions in the transmembrane region of class A G protein-coupled receptors

Residues in the transmembrane region of G protein-coupled receptors (GPCRs) are important for ligand binding and activation, but the function of individual positions is poorly understood. Using a sequence alignment of class A GPCRs (grouped in subfamilies), we propose a so-called "two-entropies analysis" to determine the potential role of individual positions in the transmembrane region of class A GPCRs. In our approach, such positions appear scattered, while largely clustered according to their biological function. Our method appears superior when compared to other bioinformatics approaches, such as the evolutionary trace method, entropy-variability plot, and correlated mutation analysis, both qualitatively and quantitatively.