Semiclosed Conformations of the Ligand-Binding Domains of NMDA Receptors during Stationary Gating

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.

NMDA Receptor Opening and Closing-Transitions of a Molecular Machine Revealed by Molecular Dynamics

We report the first complete description of the molecular mechanisms behind the transition of the N-methyl-d-aspartate (NMDA) receptor from the state where the transmembrane domain (TMD) and the ion channel are in the open configuration to the relaxed unliganded state where the channel is closed. Using an aggregate of nearly 1 µs of unbiased all-atom implicit membrane and solvent molecular dynamics (MD) simulations we identified distinct structural states of the NMDA receptor and revealed functionally important residues (GluN1/Glu522, GluN1/Arg695, and GluN2B/Asp786). The role of the "clamshell" motion of the ligand binding domain (LBD) lobes in the structural transition is supplemented by the observed structural similarity at the level of protein domains during the structural transition, combined with the overall large rearrangement necessary for the opening and closing of the receptor. The activated and open states of the receptor are structurally similar to the liganded crystal structure, while in the unliganded receptor the extracellular domains perform rearrangements leading to a clockwise rotation of up to 45 degrees around the longitudinal axis of the receptor, which closes the ion channel. The ligand-induced rotation of extracellular domains transferred by LBD-TMD linkers to the membrane-anchored ion channel is responsible for the opening and closing of the transmembrane ion channel, revealing the properties of NMDA receptor as a finely tuned molecular machine.

Effect of mutant NMDA receptors on the oscillations in a model of hippocampus

A computation approach to identify the effect of missense mutations on the protein function is proposed. Using molecular dynamics simulation we have analyzed the gating kinetics of mutant NMDA synaptic receptors carrying mutations in their NR2 subunits. Analysis of channel geometry and Mg ion binding allowed to estimate the receptor conductivity. As a result, it was possible to identify the effect of these mutations on the generation of theta and gamma rhythms by the hippocampal neural network. Obtained results can be adapted for the analysis and evaluation of possible cognitive impairments caused by neurological diseases or consequences of radiation and other negative factors.

Computer Simulations Predict High Structural Heterogeneity of Functional State of NMDA Receptors

N-methyl-D-aspartate receptors (NMDARs)-i.e., transmembrane proteins expressed in neurons-play a central role in the molecular mechanisms of learning and memory formation. It is unclear how the known atomic structures of NMDARs determined by x-ray crystallography and electron cryomicroscopy (18 published Protein Data Bank entries) relate to the functional states of NMDARs inferred from electrophysiological recordings (multiple closed, open, preopen, etc. states). We address this problem by using molecular dynamics simulations at atomic resolution, a method successfully applied in the past to much smaller biomolecules. Our simulations predict that several conformations of NMDARs with experimentally determined geometries, including four "nonactive" electron cryomicroscopy structures, rapidly interconvert on submicrosecond timescales and therefore may correspond to the same functional state of the receptor (specifically, one of the closed states). This conclusion is not trivial because these conformational transitions involve changes in certain interatomic distances as large as tens of Å. The simulations also predict differences in the conformational dynamics of the apo and holo (i.e., agonist and coagonist bound) forms of the receptor on the microsecond timescale. To our knowledge, five new conformations of NMDARs, with geometries joining various features from different known experimental structures, are also predicted by the model. The main limitation of this work stems from its limited sampling (30 μs of aggregate length of molecular dynamics trajectories). Though this level significantly exceeds the sampling in previous simulations of parts of NMDARs, it is still much lower than the sampling recently achieved for smaller biomolecules (up to a few milliseconds), thus precluding, in particular, the observation of transitions between different functional states of NMDARs. Despite this limitation, such computational predictions may guide further experimental studies on the structure, dynamics, and function of NMDARs, for example by suggesting optimal locations of spectroscopic probes. Overall, atomic resolution simulations provide, to our knowledge, a novel perspective on the structure and dynamics of NMDARs, complementing information obtained by experimental methods.

Glutamate and Glycine Binding to the NMDA Receptor

At central nervous system synapses, agonist binding to postsynaptic ionotropic glutamate receptors (iGluRs) results in signaling between neurons. N-Methyl-D-aspartic acid (NMDA) receptors are a unique family of iGluRs that activate in response to the concurrent binding of glutamate and glycine. Here, we investigate the process of agonist binding to the GluN2A (glutamate binding) and GluN1 (glycine binding) NMDA receptor subtypes using long-timescale unbiased molecular dynamics simulations. We find that positively charged residues on the surface of the GluN2A ligand-binding domain (LBD) assist glutamate binding via a "guided-diffusion" mechanism, similar in fashion to glutamate binding to the GluA2 LBD of AMPA receptors. Glutamate can also bind in an inverted orientation. Glycine, on the other hand, binds to the GluN1 LBD via an "unguided-diffusion" mechanism, whereby glycine finds its binding site primarily by random thermal fluctuations. Free energy calculations quantify the glutamate- and glycine-binding processes.

Kinetic models for activation and modulation of NMDA receptor subtypes

NMDA receptors are a diverse family of excitatory channels with critical roles in central synaptic transmission, development, and plasticity. Controlled expression of seven subunits and their combinatorial assembly into tetrameric receptors produces a range of molecularly distinct receptor subtypes. Despite relatively similar atomic structures, each subtype has input-output functions with unique biophysical and pharmacologic profiles. Here, we briefly summarize recent advances in understanding how gating and allosteric modulation are similar or distinct across NMDA receptor isoforms and identify open questions that will focus research in this area going forward.

Neurotransmitter Funneling Optimizes Glutamate Receptor Kinetics

Ionotropic glutamate receptors (iGluRs) mediate neurotransmission at the majority of excitatory synapses in the brain. Little is known, however, about how glutamate reaches the recessed binding pocket in iGluR ligand-binding domains (LBDs). Here we report the process of glutamate binding to a prototypical iGluR, GluA2, in atomistic detail using unbiased molecular simulations. Charged residues on the LBD surface form pathways that facilitate glutamate binding by effectively reducing a three-dimensional diffusion process to a spatially constrained, two-dimensional one. Free energy calculations identify residues that metastably bind glutamate and help guide it into the binding pocket. These simulations also reveal that glutamate can bind in an inverted conformation and also reorient while in its pocket. Electrophysiological recordings demonstrate that eliminating these transient binding sites slows activation and deactivation, consistent with slower glutamate binding and unbinding. These results suggest that binding pathways have evolved to optimize rapid responses of AMPA-type iGluRs at synapses.

Energetics of Glutamate Binding to an Ionotropic Glutamate Receptor

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that are responsible for the majority of excitatory transmission at the synaptic cleft. Mechanically speaking, agonist binding to the ligand binding domain (LBD) activates the receptor by triggering a conformational change that is transmitted to the transmembrane region, opening the ion channel pore. We use fully atomistic molecular dynamics simulations to investigate the binding process in the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor, an iGluR subtype. The string method with swarms of trajectories was applied to calculate the possible pathways glutamate traverses during ligand binding. Residues peripheral to the binding cleft are found to metastably bind the ligand prior to ligand entry into the binding pocket. Umbrella sampling simulations were performed to compute the free energy barriers along the binding pathways. The calculated free energy profiles demonstrate that metastable interactions contribute substantially to the energetics of ligand binding and form local minima in the overall free energy landscape. Protein-ligand interactions at sites outside of the orthosteric agonist-binding site may serve to lower the transition barriers of the binding process.

Simulated Dynamics of Glycans on Ligand-Binding Domain of NMDA Receptors Reveals Strong Dynamic Coupling between Glycans and Protein Core

N-Methyl-d-aspartate (NMDA) receptors, key neuronal receptors playing the central role in learning and memory, are heavily glycosylated in vivo. Astonishingly little is known about the structure, dynamics, and physiological relevance of glycans attached to them. We recently demonstrated that certain glycans on the ligand binding domain (LBD) of NMDA receptors (NMDARs) can serve as intramolecular potentiators, changing EC50 of NMDAR coagonists. In this work, we use molecular dynamics trajectories, in aggregate 86.5 μs long, of the glycosylated LBD of the GluN1 subunit of the NMDAR to investigate the behavior of glycans on NMDARs. Though all glycans in our simulations were structurally the same (Man5), the dynamics of glycans at different locations on NMDARs was surprisingly different. The slowest-time scale motions that we detected in various glycans in some cases corresponded to a flipping of parts of glycans relative to each other, while in other cases they reduced to a head-to-tail bending of a glycan. We predict that time scales of conformational changes in glycans on the GluN1 LBD of NMDARs range from nanoseconds to at least hundreds of microseconds. Some of the conformational changes in the glycans correlate with the physiologically important clamshell-like opening and closing of the GluN1 LBD domain. Thus, glycans are an integral part of NMDARs, and computational models of NMDARs should include glycans to faithfully represent the structure and the dynamics of these receptors.

Probing the Structural Dynamics of the NMDA Receptor Activation by Coarse-Grained Modeling

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. High-resolution structures have been solved for an allosterically inhibited and agonist-bound form of a functional NMDA receptor; however, other key functional states (particularly the active open-channel state) were only resolved at moderate resolutions by cryo-electron microscopy (cryo-EM). To decrypt the mechanism of the NMDA receptor activation, structural modeling is essential to provide presently missing information about structural dynamics. We performed systematic coarse-grained modeling using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on an active-state cryo-EM map. We observed extensive conformational changes that allosterically couple the extracellular regulatory and agonist-binding domains to the pore-forming trans-membrane domain (TMD), and validated these, to our knowledge, new observations against known mutational and functional studies. Our results predict two key modes of collective motions featuring shearing/twisting of the extracellular domains relative to the TMD, reveal subunit-specific flexibility profiles, and identify functional hotspot residues at key domain-domain interfaces. Finally, by examining the conformational transition pathway between the allosterically inhibited form and the active form, we predict a discrete sequence of domain motions, which propagate from the extracellular domains to the TMD. In summary, our results offer rich structural and dynamic information, which is consistent with the literature on structure-function relationships in NMDA receptors, and will guide in-depth studies on the activation dynamics of this important neurotransmitter receptor.

Allosteric Interactions between NMDA Receptor Subunits Shape the Developmental Shift in Channel Properties

During development of the central nervous system, there is a shift in the subunit composition of NMDA receptors (NMDARs) resulting in a dramatic acceleration of NMDAR-mediated synaptic currents. This shift coincides with upregulation of the GluN2A subunit and triheteromeric GluN1/2A/2B receptors with fast deactivation kinetics, whereas expression of diheteromeric GluN1/2B receptors with slower deactivation kinetics is decreased. Here, we show that allosteric interactions occur between the glutamate-binding GluN2 subunits in triheteromeric GluN1/2A/2B NMDARs. This allosterism is dominated by the GluN2A subunit and results in functional properties not predicted by those of diheteromeric GluN1/2A and GluN1/2B NMDARs. These findings suggest that GluN1/2A/2B NMDARs may maintain some signaling properties of the GluN2B subunit while having the kinetic properties of GluN1/2A NMDARs and highlight the complexity in NMDAR signaling created by diversity in subunit composition.

NMDA receptors: linking physiological output to biophysical operation

NMDA receptors are preeminent neurotransmitter-gated channels in the CNS, which respond to glutamate in a manner that integrates multiple external and internal cues. They belong to the ionotropic glutamate receptor family and fulfil unique and crucial roles in neuronal development and function. These roles depend on characteristic response kinetics, which reflect the operation of the receptors. Here, we review biologically salient features of the NMDA receptor signal and its mechanistic origins. Knowledge of distinctive NMDA receptor biophysical properties, their structural determinants and physiological roles is necessary to understand the physiological and neurotoxic actions of glutamate and to design effective therapeutics.

Gating Motions and Stationary Gating Properties of Ionotropic Glutamate Receptors: Computation Meets Electrophysiology

Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels essential to all aspects of brain function, including higher order processes such as learning and memory. For decades, electrophysiology was the primary means for characterizing the function of iGluRs and gaining mechanistic insight. Since the turn of the century, structures of isolated water-soluble domains and transmembrane-domain-containing constructs have provided the basis for formulating mechanistic hypotheses. Because these structures only represent sparse, often incomplete snapshots during iGluR activation, significant gaps in knowledge remain regarding structures, energetics, and dynamics of key substates along the functional processes. Some of these gaps have recently been filled by molecular dynamics simulations and theoretical modeling. In this Account, I describe our work in the latter arena toward characterizing iGluR gating motions and developing a formalism for calculating thermodynamic and kinetic properties of stationary gating. The structures of iGluR subunits have a highly modular architecture, in which the ligand-binding domain and the transmembrane domain are well separated and connected by flexible linkers. The ligand-binding domain in turn is composed of two subdomains. During activation, agonist binding induces the closure of the intersubdomain cleft. The cleft closure leads to the outward pulling of a linker tethered to the extracellular terminus of the major pore-lining helix of the transmembrane domain, thereby opening the channel. This activation model based on molecular dynamics simulations was validated by residue-specific information from electrophysiological data on cysteine mutants. A further critical test was made through introducing glycine insertions in the linker. Molecular dynamics simulations showed that, with lengthening by glycine insertions, the linker became less effective in pulling the pore-lining helix, leading to weaker stabilization of the channel-open state. In full agreement, single-channel recordings showed that the channel open probability decreased progressively as the linker was lengthened by glycine insertions. Crystal structures of ligand-binding domains showing different degrees of cleft closure between full and partial agonists suggested a simple mechanism for one subtype of iGluRs, but mysteries surrounded a second subtype, where the ligand-binding domains open to similar degrees when bound with either full or partial agonists. Our free energy simulations now suggest that broadening of the free energy basin for cleft closure is a plausible solution. A theoretical basis for these mechanistic hypotheses on partial agonisms was provided by a model for the free energy surface of a full receptor, where the stabilization by cleft closure is transmitted via the linker to the channel-open state. This model can be implemented by molecular dynamics simulations to predict thermodynamic and kinetics properties of stationary gating that are amenable to direct test by single-channel recordings. Close integration between computation and electrophysiology holds great promises in revealing the conformations of key substates in functional processes and the mechanisms of disease-associated mutations.

Computationally Discovered Potentiating Role of Glycans on NMDA Receptors

N-methyl-D-aspartate receptors (NMDARs) are glycoproteins in the brain central to learning and memory. The effects of glycosylation on the structure and dynamics of NMDARs are largely unknown. In this work, we use extensive molecular dynamics simulations of GluN1 and GluN2B ligand binding domains (LBDs) of NMDARs to investigate these effects. Our simulations predict that intra-domain interactions involving the glycan attached to residue GluN1-N440 stabilize closed-clamshell conformations of the GluN1 LBD. The glycan on GluN2B-N688 shows a similar, though weaker, effect. Based on these results, and assuming the transferability of the results of LBD simulations to the full receptor, we predict that glycans at GluN1-N440 might play a potentiator role in NMDARs. To validate this prediction, we perform electrophysiological analysis of full-length NMDARs with a glycosylation-preventing GluN1-N440Q mutation, and demonstrate an increase in the glycine EC50 value. Overall, our results suggest an intramolecular potentiating role of glycans on NMDA receptors.

Advancing NMDA Receptor Physiology by Integrating Multiple Approaches

NMDA receptors (NMDARs) are ion channels activated by the excitatory neurotransmitter glutamate and are essential to all aspects of brain function, including learning and memory formation. Missense mutations distributed throughout NMDAR subunits have been associated with an array of neurological disorders. Recent structural, functional, and computational studies have generated many insights into the activation process connecting glutamate binding to ion-channel opening, which is central to NMDAR physiology and pathophysiology. The field appears poised for breakthroughs, including the exciting prospect of resolving the conformations and energetics of elementary steps in the activation process, and atomic-level modeling of the effects of missense mutations on receptor function. The most promising strategy going forward is through strong integration of multiple approaches.