Dynamics of cleft closure of the GluA2 ligand-binding domain in the presence of full and partial agonists revealed by hydrogen-deuterium exchange

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

Structure, Function, and Regulation of the Kainate Receptor

Neural communication and modulation are complex processes. Ionotropic glutamate receptors (iGluRs) significantly contribute to mediating the fast-excitatory branch of neurotransmission in the mammalian brain. Kainate receptors (KARs), a subfamily of the iGluRs, act as modulators of the neuronal circuitry by playing important roles at both the post- and presynaptic sites of specific neurons. The functional tetrameric receptors are formed by two different gene families, low agonist affinity (GluK1-GluK3) and high agonist affinity (GluK4-GluK5) subunits. These receptors garnered attention in the past three decades, and since then, much work has been done to understand their localization, interactome, physiological functions, and regulation. Cloning of the receptor subunits (GluK1-GluK5) in the early 1990s led to recombinant expression of kainate receptors in heterologous systems. This facilitated understanding of the functional differences between subunit combinations, splice variants, trafficking, and drug discovery. Structural studies of individual domains and recent full-length homomeric and heteromeric kainate receptors have revealed unique functional mechanisms, which have answered several long-standing questions in the field of kainate receptor biology. In this chapter, we review the current understanding of kainate receptors and associated disorders.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

Structural Mechanisms Underlying Activity Changes in an AMPA-type Glutamate Receptor Induced by Substitutions in Its Ligand-Binding Domain

α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)-type glutamate receptors produce postsynaptic current by transmitting an agonist-induced structural change in the ligand-binding domain (LBD) to the transmembrane channel. Receptors carrying T686S/A substitutions in their LBDs produce weaker glutamate-evoked currents than wild-type (WT) receptors. However, the substitutions induce little differences in the crystal structures of their LBDs. To understand the structural mechanism underlying reduced activities of these AMPAR variants, we analyzed the structural dynamics of WT, T686S, and T686A variants of LBD using nuclear magnetic resonance. The HD exchange studies of the LBDs showed that the kinetic step where the ligand-binding cleft closes was changed by the substitutions, and the substitution-induced population shift from cleft-closed to cleft-open structures is responsible for the reduced activities of the variants. The chemical shift analyses revealed another structural equilibrium between cleft-locked and cleft-partially-open conformations. The substitution-induced population shift in this equilibrium may be related to slower desensitization observed for these variants.

Structure and mechanism of AMPA receptor - auxiliary protein complexes

Ionotropic glutamate receptors in vertebrates are composed of three major subtypes - AMPA, kainate, and NMDA receptors - and mediate the majority of fast excitatory neurotransmission at chemical synapses of the central nervous system. Among the three major families, native AMPA receptors function as complexes with a variety of auxiliary subunits, which in turn modulate receptor trafficking, gating, pharmacology, and permeation. Despite the long history of structure-mechanism studies using soluble receptor domains or intact yet isolated receptors, structures of AMPA receptor-auxiliary subunit complexes have not been available until recent breakthroughs in single-particle cryo-electron microscopy. Single particle cryo-EM studies have, in turn, provided new insights into the structure and organization of AMPA receptor - auxiliary protein complexes and into the molecular mechanisms of AMPA receptor activation and desensitization.

Activation and Desensitization Mechanism of AMPA Receptor-TARP Complex by Cryo-EM

AMPA receptors mediate fast excitatory neurotransmission in the mammalian brain and transduce the binding of presynaptically released glutamate to the opening of a transmembrane cation channel. Within the postsynaptic density, however, AMPA receptors coassemble with transmembrane AMPA receptor regulatory proteins (TARPs), yielding a receptor complex with altered gating kinetics, pharmacology, and pore properties. Here, we elucidate structures of the GluA2-TARP γ2 complex in the presence of the partial agonist kainate or the full agonist quisqualate together with a positive allosteric modulator or with quisqualate alone. We show how TARPs sculpt the ligand-binding domain gating ring, enhancing kainate potency and diminishing the ensemble of desensitized states. TARPs encircle the receptor ion channel, stabilizing M2 helices and pore loops, illustrating how TARPs alter receptor pore properties. Structural and computational analysis suggests the full agonist and modulator complex harbors an ion-permeable channel gate, providing the first view of an activated AMPA receptor.

Role of the Ion Channel Extracellular Collar in AMPA Receptor Gating

AMPA subtype ionotropic glutamate receptors mediate fast excitatory neurotransmission and are implicated in numerous neurological diseases. Ionic currents through AMPA receptor channels can be allosterically regulated via different sites on the receptor protein. We used site-directed mutagenesis and patch-clamp recordings to probe the ion channel extracellular collar, the binding region for noncompetitive allosteric inhibitors. We found position and substitution-dependent effects for introduced mutations at this region on AMPA receptor gating. The results of mutagenesis suggested that the transmembrane domains M1, M3 and M4, which contribute to the ion channel extracellular collar, undergo significant relative displacement during gating. We used molecular dynamics simulations to predict an AMPA receptor open state structure and rationalize the results of mutagenesis. We conclude that the ion channel extracellular collar plays a distinct role in gating and represents a hub for powerful allosteric modulation of AMPA receptor function that can be used for developing novel therapeutics.

An improved quantum biochemistry description of the glutamate–GluA2 receptor binding within an inhomogeneous dielectric function framework

A new methodology to define the inhomogeneous dielectric constant of protein residues, to apply to the calculation of protein–ligand properties such as the electrostatic interaction.

Dynamics of the Ligand Binding Domain Layer during AMPA Receptor Activation

Ionotropic glutamate receptors are postsynaptic tetrameric ligand-gated channels whose activity mediates fast excitatory transmission. Glutamate binding to clamshell-shaped ligand binding domains (LBDs) triggers opening of the integral ion channel, but how the four LBDs orchestrate receptor activation is unknown. Here, we present a high-resolution x-ray crystal structure displaying two tetrameric LBD arrangements fully bound to glutamate. Using a series of engineered metal ion trapping mutants, we showed that the more compact of the two assemblies corresponds to an arrangement populated during activation of full-length receptors. State-dependent cross-linking of the mutants identified zinc bridges between the canonical active LBD dimers that formed when the tetramer was either fully or partially bound by glutamate. These bridges also stabilized the resting state, consistent with the recently published full-length apo structure. Our results provide insight into the activation mechanism of glutamate receptors and the complex conformational space that the LBD layer can sample.

Conformational Selection and Submillisecond Dynamics of the Ligand-binding Domain of the N-Methyl-d-aspartate Receptor

The N-methyl-d-aspartate (NMDA) receptors are heteromeric non-selective cation channels that require the binding of glycine and glutamate for gating. Based on crystal structures, the mechanism of partial agonism at the glycine-binding site is thought to be mediated by a shift in the conformational equilibrium between an open clamshell and a closed clamshell-like structure of the bilobed ligand-binding domain (LBD). Using single-molecule Förster resonance energy transfer (smFRET) and multiparameter fluorescence detection, which allows us to study the conformational states and dynamics in the submillisecond time scale, we show that there are at least three conformational states explored by the LBD: the low FRET, medium FRET, and high FRET states. The distance of the medium and low FRET states corresponds to what has been observed in crystallography structures. We show that the high FRET state, which would represent a more closed clamshell conformation than that observed in the crystal structure, is most likely the state initiating activation, as evidenced by the fact that the fraction of the protein in this state correlates well with the extent of activation. Furthermore, full agonist bound LBDs show faster dynamic motions between the medium and high FRET states, whereas they show slower dynamics when bound to weaker agonists or to antagonists.

Tweaking Subtype Selectivity and Agonist Efficacy at (S)-2-Amino-3-(3-hydroxy-5-methyl-isoxazol-4-yl)propionic acid (AMPA) Receptors in a Small Series of BnTetAMPA Analogues

A series of analogues of the (S)-2-Amino-3-(3-hydroxy-5-methyl-isoxazol-4-yl)propionic acid (AMPA) receptor agonist BnTetAMPA (5b) were synthesized and characterized pharmacologically in radioligand binding assays at native and cloned AMPA receptors and functionally by two-electrode voltage clamp electrophysiology at the four homomeric AMPA receptors expressed in Xenopus laevis oocytes. The analogues 6 and 7 exhibit very different pharmacological profiles with binding affinity preference for the subtypes GluA1 and GluA3, respectively. X-ray crystal structures of three ligands (6, 7, and 8) in complex with the agonist binding domain (ABD) of GluA2 show that they induce full domain closure despite their low agonist efficacies. Trp767 in GluA2 ABD could be an important determinant for partial agonism of this compound series at AMPA receptors, since agonist efficacy also correlated with the location of the Trp767 side chain.

Mechanism-Based Mathematical Model for Gating of Ionotropic Glutamate Receptors

We present a mathematical model for ionotropic glutamate receptors (iGluR's) that is built on mechanistic understanding and yields a number of thermodynamic and kinetic properties of channel gating. iGluR's are ligand-gated ion channels responsible for the vast majority of fast excitatory neurotransmission in the central nervous system. The effects of agonist-induced closure of the ligand-binding domain (LBD) are transmitted to the transmembrane channel (TMC) via interdomain linkers. Our model demonstrates that, relative to full agonists, partial agonists may reduce either the degree of LBD closure or the curvature of the LBD free energy basin, leading to less stabilization of the channel open state and hence lower channel open probability. A rigorous relation is derived between the channel closed-to-open free energy difference and the tension within the linker. Finally, by treating LBD closure and TMC opening as diffusive motions, we obtain gating trajectories that resemble stochastic current traces from single-channel recordings and calculate the rate constants for transitions between the channel open and closed states. Our model can be implemented by molecular dynamics simulations to realistically depict iGluR gating and may guide functional experiments in gaining deeper insight into this essential family of channel proteins.

A quantum biochemistry investigation of willardiine partial agonism in AMPA receptors

We employ quantum biochemistry methods based on the Density Functional Theory (DFT) approach to unveil detailed binding energy features of willardiines co-crystallized with the AMPA receptor.

Reduced curvature of ligand-binding domain free-energy surface underlies partial agonism at NMDA receptors

NMDA receptors are ligand-gated ion channels that mediate excitatory synaptic transmission in the central nervous system. Partial agonists elicit submaximal channel activation, but crystal structures of the ligand-binding domains (LBDs) bound with partial and full agonists show little difference. To uncover the molecular mechanism for partial agonism, here we computed the free-energy surfaces of the GluN1 (an obligatory subunit of NMDA receptors) LBD bound with a variety of ligands. The free-energy minima are similarly positioned for full and partial agonists, but the curvatures are significantly reduced in the latter case, indicating higher probabilities for sampling conformations with a not fully closed domain cleft. The free-energy surfaces for antagonists have both shifted minima and further reduced curvatures. Reduced curvature of free-energy surface appears to explain well the partial agonism at NMDA receptors and may present a unique paradigm in producing graded responses for receptors in general.

Structural dynamics of the glycine-binding domain of the N-methyl-D-aspartate receptor

N-Methyl-D-aspartate receptors mediate the slow component of excitatory neurotransmission in the central nervous system. These receptors are obligate heteromers containing glycine- and glutamate-binding subunits. The ligands bind to a bilobed agonist-binding domain of the receptor. Previous x-ray structures of the glycine-binding domain of NMDA receptors showed no significant changes between the partial and full agonist-bound structures. Here we have used single molecule fluorescence resonance energy transfer (smFRET) to investigate the cleft closure conformational states that the glycine-binding domain of the receptor adopts in the presence of the antagonist 5,7-dichlorokynurenic acid (DCKA), the partial agonists 1-amino-1-cyclobutanecarboxylic acid (ACBC) and L-alanine, and full agonists glycine and D-serine. For these studies, we have incorporated the unnatural amino acid p-acetyl-L-phenylalanine for specific labeling of the protein with hydrazide derivatives of fluorophores. The single molecule fluorescence resonance energy transfer data show that the agonist-binding domain can adopt a wide range of cleft closure states with significant overlap in the states occupied by ligands of varying efficacy. The difference lies in the fraction of the protein in a more closed-cleft form, with full agonists having a larger fraction in the closed-cleft form, suggesting that the ability of ligands to select for these states could dictate the extent of activation.

Structure of an agonist-bound ionotropic glutamate receptor

Ionotropic glutamate receptors (iGluRs) mediate most excitatory neurotransmission in the central nervous system and function by opening their ion channel in response to binding of agonist glutamate. Here, we report a structure of a homotetrameric rat GluA2 receptor in complex with partial agonist (S)-5-nitrowillardiine. Comparison of this structure with the closed-state structure in complex with competitive antagonist ZK 200775 suggests conformational changes that occur during iGluR gating. Guided by the structures, we engineered disulfide cross-links to probe domain interactions that are important for iGluR gating events. The combination of structural information, kinetic modeling, and biochemical and electrophysiological experiments provides insight into the mechanism of iGluR gating.

Deuterium content of water increases depression susceptibility: the potential role of a serotonin-related mechanism

Environmental factors can significantly affect disease prevalence, including neuropsychiatric disorders such as depression. The ratio of deuterium to protium in water shows substantial geographical variation, which could affect disease susceptibility. Thus the link between deuterium content of water and depression was investigated, both epidemiologically, and in a mouse model of chronic mild stress. We performed a correlation analysis between deuterium content of tap water and rates of depression in regions of the USA. Next, we used a 10-day chronic stress paradigm to test whether 2-week deuterium-depleted water treatment (91 ppm) affects depressive-like behavior and hippocampal SERT. The effect of deuterium-depletion on sleep electrophysiology was also evaluated in naïve mice. There was a geographic correlation between a content of deuterium and the prevalence of depression across the USA. In the chronic stress model, depressive-like features were reduced in mice fed with deuterium-depleted water, and SERT expression was decreased in mice treated with deuterium-treated water compared with regular water. Five days of predator stress also suppressed proliferation in the dentate gyrus; this effect was attenuated in mice fed with deuterium-depleted water. Finally, in naïve mice, deuterium-depleted water treatment increased EEG indices of wakefulness, and decreased duration of REM sleep, phenomena that have been shown to result from the administration of selective serotonin reuptake inhibitors (SSRI). Our data suggest that the deuterium content of water may influence the incidence of affective disorder-related pathophysiology and major depression, which might be mediated by the serotoninergic mechanisms.

The future of drug discovery: enabling technologies for enhancing lead characterization and profiling therapeutic potential

Technology often serves as a handmaiden and catalyst of invention. The discovery of safe, effective medications depends critically upon experimental approaches capable of providing high-impact information on the biological effects of drug candidates early in the discovery pipeline. This information can enable reliable lead identification, pharmacological compound differentiation and successful translation of research output into clinically useful therapeutics. The shallow preclinical profiling of candidate compounds promulgates a minimalistic understanding of their biological effects and undermines the level of value creation necessary for finding quality leads worth moving forward within the development pipeline with efficiency and prognostic reliability sufficient to help remediate the current pharma-industry productivity drought. Three specific technologies discussed herein, in addition to experimental areas intimately associated with contemporary drug discovery, appear to hold particular promise for strengthening the preclinical valuation of drug candidates by deepening lead characterization. These are: i) hydrogen-deuterium exchange mass spectrometry for characterizing structural and ligand-interaction dynamics of disease-relevant proteins; ii) activity-based chemoproteomics for profiling the functional diversity of mammalian proteomes; and iii) nuclease-mediated precision gene editing for developing more translatable cellular and in vivo models of human diseases. When applied in an informed manner congruent with the clinical understanding of disease processes, technologies such as these that span levels of biological organization can serve as valuable enablers of drug discovery and potentially contribute to reducing the current, unacceptably high rates of compound clinical failure.

Thermodynamics and mechanism of the interaction of willardiine partial agonists with a glutamate receptor: implications for drug development

Understanding the thermodynamics of binding of a lead compound to a receptor can provide valuable information for drug design. The binding of compounds, particularly partial agonists, to subtypes of the α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor is, in some cases, driven by increases in entropy. Using a series of partial agonists based on the structure of the natural product, willardiine, we show that the charged state of the ligand determines the enthalpic contribution to binding. Willardiines have uracil rings with pKa values ranging from 5.5 to 10. The binding of the charged form is largely driven by enthalpy, while that of the uncharged form is largely driven by entropy. This is due at least in part to changes in the hydrogen bonding network within the binding site involving one water molecule. This work illustrates the importance of charge to the thermodynamics of binding of agonists and antagonists to AMPA receptors and provides clues for further drug discovery.

Structure and gating of tetrameric glutamate receptors

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that open their ion-conducting pores in response to the binding of agonist glutamate. In recent years, significant progress has been achieved in studies of iGluRs by determining numerous structures of isolated water-soluble ligand-binding and amino-terminal domains, as well as solving the first crystal structure of the full-length AMPA receptor in the closed, antagonist-bound state. These structural data combined with electrophysiological and fluorescence recordings, biochemical experiments, mutagenesis and molecular dynamics simulations have greatly improved our understanding of iGluR assembly, activation and desensitization processes. This article reviews the recent structural and functional advances in the iGluR field and summarizes them in a simplified model of full-length iGluR gating.