The NMDA receptor (NMDAR) is an ionotropic glutamate receptor formed from the tetrameric assembly of GluN1 and GluN2 subunits. Within the flexible linker between the agonist binding domain (ABD) and the M1 helix of the pore-forming transmembrane helical bundle lies a two-turn, extracellular pre-M1 helix positioned parallel to the plasma membrane and in van der Waals contact with the M3 helix thought to constitute the channel gate. The pre-M1 helix is tethered to the bilobed ABD, where agonist-induced conformational changes initiate activation. Additionally, it is a locus for de novo mutations associated with neurological disorders, is near other disease-associated de novo sites within the transmembrane domain, and is a structural determinant of subunit-selective modulators. To investigate the role of the pre-M1 helix in channel gating, we performed scanning mutagenesis across the GluN2A pre-M1 helix and recorded whole-cell macroscopic and single channel currents from HEK293 cell-attached patches. We identified two residues at which mutations perturb channel open probability, the mean open time, and the glutamate deactivation time course. We identified a subunit-specific network of aromatic amino acids located in and around the GluN2A pre-M1 helix to be important for gating. Based on these results, we are able to hypothesize about the role of the pre-M1 helix in other NMDAR subunits based on sequence and structure homology. Our results emphasize the role of the pre-M1 helix in channel gating, implicate the surrounding amino acid environment in this mechanism, and suggest unique subunit-specific contributions of pre-M1 helices to GluN1 and GluN2 gating.
NMDA receptors (NMDARs) are a subfamily of ionotropic glutamate receptors whose calcium-permeable channels mediate the slow component of fast excitatory neurotransmission (Lester et al., 1990; Jahr and Stevens, 1993; Clements, 1996). These receptors, although necessary for synaptic plasticity, learning, and memory, have been implicated in the etiology of several neurological disorders, including epilepsy, intellectual disabilities, and schizophrenia (Coyle, 2006; Liu et al., 2004; Rice and DeLorenzo, 1998; Olney et al., 1999; Lemke et al., 2014; Traynelis et al., 2010; Yuan et al., 2015; XiangWei et al., 2018). NMDARs are formed from the tetramerization of two glycine-binding GluN1 subunits and two glutamate-binding GluN2 subunits. The GluN1 subunit is encoded by a single gene with eight possible splice variants. The GluN2 subunits, however, can be one of four discretely expressed gene products (A–D), each providing the receptor with unique pharmacological, physiological, and functional properties (Monyer et al., 1992; Cull-Candy et al., 2001; Traynelis et al., 2010; Paoletti et al., 2013).
Each subunit contains four semiautonomous domains: the amino terminal domain (ATD), the agonist binding domain (ABD), the transmembrane domain (TMD), and the carboxy-terminal domain. During a process referred to as gating, agonist binding to the ABD promotes channel opening within the TMD. Crystallographic data suggest that the bilobed ABDs initiate activation by adopting a closed-cleft conformation due to atomic contacts between the agonist, D1 upper lobe, and D2 lower lobe (Armstrong and Gouaux, 2000; Sun et al., 2002; Jin et al., 2003; Inanobe et al., 2005; Furukawa et al., 2005; Vance et al., 2011; Mayer, 2011; Karakas and Furukawa, 2014; Lee et al., 2014). Investigation of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor (AMPAR) by single molecule fluorescence resonance energy transfer has also contributed toward the understanding of site-specific transitions between multiple structurally distinct states. The ATD, for example, transitions between at least three distinct conformations in the desensitized state with interdimer distances ranging from 35 to 53 Å (MacLean et al., 2019). Moreover, single molecule fluorescence resonance energy transfer has been used to reveal conformational changes within the NMDAR transmembrane region associated with channel opening, which could account for the multiple receptor states observed in single channel recordings (Dolino et al., 2017). Gating of the AMPAR, a structurally similar member of the ionotropic glutamate receptor family, involves a “kinking” of two nonadjacent, pore-forming M3 transmembrane helices away from the central axis of the pore (Twomey et al., 2017; Twomey and Sobolevsky, 2018). The AMPAR M3 helices contain a highly conserved amino acid sequence (SYTANLAAF), which can be found in all members of the ionotropic glutamate receptor family, including NMDARs. Finally, NMDARs exhibit similar asymmetry to that observed during AMPAR opening (Lee et al., 2014; Tajima et al., 2016), suggesting that this kinking motion also occurs for NMDARs.
Although the mechanism of NMDAR gating remains poorly understood, it has been hypothesized that the linker tethering the ABD to the TMD is responsible for facilitating communication between these two domains (Schorge et al., 2005; Sobolevsky et al., 2009; Talukder et al., 2010; Talukder and Wollmuth, 2011; Karakas and Furukawa, 2014; Lee et al., 2014; Gibb et al., 2018; Amin et al., 2020). Thus far, the flexible and dynamic nature of the ABD-TMD linker has prevented high resolution structure of this region, and the structure of the fully open NMDAR pore remains elusive (Tajima et al., 2016). Lower resolution analyses have revealed novel structural features within the linker that are uniquely positioned to transduce agonist binding into channel opening (Sobolevsky et al., 2009; Karakas and Furukawa, 2014; Lee et al., 2014).
Specifically, the linker between the ABD and the first transmembrane helix (M1) contains a short, two-turn pre-M1 helix that is thought to be a key contributor to the gating mechanism. The pre-M1 helix, previously described as a “cuff” helix, is hypothesized to participate in channel gating based on its orientation perpendicular to the channel pore, its position downstream of the ABD, and its van der Waals contact with the extracellular end of the M3 helical bundle (Sobolevsky et al., 2009; Karakas et al., 2015; Tajima et al., 2016). It stands to reason that agonist-induced closure of the bilobed ABD alters the position of the ABD-TMD linkers, leading to a reorientation of both the pre-M1 helices and the M3 helices of the GluN1 and GluN2 subunits (Sobolevsky et al., 2002; Furukawa and Gouaux, 2003; Furukawa et al., 2005; Inanobe et al., 2005; Vance et al., 2011; Hansen et al., 2013). In support of this hypothesis, it has been shown that cysteine substitutions within this linker region are sufficient to alter NMDAR responses, and MTS covalent modification of cysteine-substituted residues in this region altered receptor leak currents (Beck et al., 1999; Sobolevsky et al., 2007; Chang and Kuo, 2008). Therefore, if contact between the pre-M1 helix and the M3 helix is essential for stabilizing a closed channel state, then agonist-induced reorientation of the pre-M1 helix or its residue side chains may constitute a required, rate-limiting pregating step that must be traversed before the channel can open (Gibb et al., 2018).
The probable involvement of the pre-M1 helix in the gating mechanism is emphasized by human genetics. In a population of >140,000 healthy individuals (available from gnomad.broadinstitute.org), the GluN2A and GluN2B pre-M1 helices are devoid of missense variants, suggesting strong selection against variation in this region. By contrast, the pre-M1 helices of these subunits harbor multiple de novo mutations in patients with epilepsy, intellectual disabilities, or developmental delays (Swanger et al., 2016; Ogden et al., 2017; XiangWei et al., 2018; Amin et al., 2020). Within the GluN2A and GluN2B subunits, functional analysis of these disease-associated mutations revealed changes to NMDAR kinetics, further implicating the pre-M1 helix in the gating mechanism (Ogden et al., 2017; Gibb et al., 2018).
It has been proposed that the pre-M1 helix is one third of a triad also comprising the M3 pore-forming helix and the pre-M4 helix of the adjacent subunit (Chen et al., 2017; Ogden et al., 2017; Gibb et al., 2018). That is, two distinct triads exist in the tetrameric assembly: GluN1–pre-M1/GluN1-M3/GluN2–pre-M4 and GluN2–pre-M1/GluN2-M3/GluN1–pre-M4. As such, there might arise structural differences in each triad as a result of the different sequences for the GluN1 or GluN2 pre-M1 and pre-M4 helices. That is, global receptor asymmetry and disparate amino acid environments may contribute to subunit-specific pre-M1 contributions to the gating mechanism.
Functional data suggest that the GluN2A pre-M1 helix plays a larger role in gating than its GluN1 counterpart. Analysis of the GluN2A-P552R mutation showed a significant delay in receptor activation and deactivation, a finding that could not be recapitulated by the GluN1 equivalent mutation GluN1-P557R (Ogden et al., 2017). Additionally, data have shown that alanine mutations within the pre-M1 helix of the GluN2D subunit significantly alter receptor open probability (Ogden and Traynelis, 2013). Because open probability is a property of the NMDAR that differs greatly between GluN2A/2B and GluN2C/2D, we might expect that residues within the pre-M1 helix, specifically those that differ between the two groups, might account for these differences. As such, we explored the GluN2A pre-M1 helix for its role in gating to better understand both its mechanistic function and how this differs from that of the pre-M1 helix in other subunits. We conducted site-directed mutagenesis coupled with two-electrode voltage clamp, single channel, and whole-cell patch clamp recordings to directly monitor changes in receptor pharmacology, response time course, and microscopic properties. We also ran MD simulations of the receptor in the presence of agonist to investigate atomic interactions between the pre-M1 helix and surrounding residues. Our results identified GluN2A-F553 as a significant contributor to receptor gating and suggest that local intra- and inter-subunit amino acid side chain interactions may be important for subunit-specific pre-M1 contributions to this mechanism.
Materials and methods
cDNAs for recombinant rat NMDAR subunits GluN1-1a (hereafter GluN1; GenBank accession no. U08261 for GluN1 lacking exon 5 but containing exons 21 and 22) and GluN2A (GenBank accession no. D13211) were subcloned into a mammalian expression vector (pcIneo). We introduced point mutations into the GluN2A cDNA using Quikchange reactions (Agilent Technologies). For expression in Xenopuslaevis oocytes, DNA constructs were linearized by restriction enzymes and used as templates for in vitro cRNA transcription using the mMessage mMachine kit (Ambion) according to the manufacturer’s recommended protocol.
Two-electrode voltage clamp recordings from Xenopus oocytes
Whole-cell current recordings from transfected HEK293 cells
Single channel recordings from transfected HEK293 cells
Cell-attached patch recordings were made at room temperature (23°C) from HEK293 cells transiently transfected with rat GluN1, GluN2, and GFP, as described above. The extracellular recording solution was the same as in whole-cell recordings, except the pH was adjusted to 8.0. Micropipettes were made using thin-walled glass as for whole-cell recordings, coated with Sylgard silicone elastomer (Dow), and fire polished to a resistance of 4–6 MΩ. Patch pipettes were filled with recording solution supplemented with 1 mM glutamate and 50 µM glycine. The holding potential was +80 mV for cell-attached patches, which corresponded to a membrane potential of approximately −100 to −120 mV, as zero current was observed at −20 to −40 mV and the reversal potential for NMDARs is 0 mV. Currents were low pass filtered at 8 kHz (−3 dB, 8-pole Bessel) and digitized at 40 kHz. Patches were determined to contain a single channel when no double openings were observed and analysis of opening frequency suggested a low probability (P < 0.001) of a multichannel recording (Colquhoun and Hawkes, 1990). Single channel records were analyzed off-line using QUB software (Qin et al., 1996; Qin et al., 1997). The critical time (Tcrit) used to chop the idealized dwells into bursts was determined as previously described (Jackson et al., 1983; Magleby and Pallotta, 1983). Individual open and closed dwell times were determined using a dead time of 0.05 ms and fitted with the sum of multiple exponential components using the maximum likelihood method. All values reported are mean ± SEM.
Homology modeling and MD
A human GluN1/GluN2A dihetero-tetrameric homology model was constructed from a glutamate/glycine-bound, exon 5–lacking GluN1a/GluN2B crystal structure (Protein Data Bank [PDB] identifier, 4PE5; 3.96 Å resolution; Karakas and Furukawa 2014), an exon 5–containing GluN1b/GluN2B crystal structure in the nonactive conformation (PDB identifier, 5FXH; 5 Å resolution; Tajima et al., 2016), and the TMD of the GluA2 crystal structure (PDB identifier, 5L1B; 4.0 Å resolution; Yelshanskaya et al., 2016). The alignment of the target and template sequences was performed using multiple sequence comparison by log-expectation (Edgar, 2004). Five homology models were generated using Modeller v9.19 (Sali and Blundell, 1993) from which the lowest energy model was selected and subjected to quality analysis using the PDBsum generator (Laskowski et al., 2018). The selected model structure is provided in Table S1. Hydrogen bond assignments were performed, followed by subsequent optimization and visual inspection. All titratable residues were assigned their dominant protonation state at pH 7.0 (Schrödinger Release 2018–4; Protein Preparation Wizard; Epik version 3.7; Impact version 7.2; Prime version 4.5, 2017; Schrödinger, LLC). Energy minimization was performed to relieve unfavorable interactions. The minimization process first involved minimizing only hydrogen atoms, followed by a restrained minimization using the imperf tool with a convergence of the RMSD of heavy atoms to 0.3 Å. Using a separate copy of this WT model, GluN2A-Phe553 was mutated to alanine for use in a second MD simulation (Schrödinger Release 2019–1, Schrödinger Suite 2019–1 Protein Preparation Wizard; Epik, 2016; Impact, 2016; Prime, 2019; Schrödinger, LLC).
To accelerate the MD simulations, each structure was truncated to exclude the ATD (residues removed up to residue 393 for the GluN1 subunits and up to residue 424 for the GluN2A subunits; Fig. 1; the structure is provided in Table S1). For both the WT and GluN2A-F553A model, the optimized and truncated receptor was inserted into an equilibrated palmitoyl oleoyl phosphatidyl choline bilayer. Membrane placement was assigned manually using the Orientations of Proteins in Membranes database (Lomize et al., 2012) entry for the 4PE5 template as a visual guide, ensuring charged residues at the intracellular end of the TMD fell outside of the hydrophobic section of the membrane. Solvation of the protein structure was done using the simple point charge water model; the solvation box included a 10-Å buffer between the protein and periodic boundary. The overall WT (F553A values in parentheses) system was neutralized by the addition of 10 Na+ ions, with the NaCl concentration set to 150 mM. The system, consisting of 1,696 amino acids, 391 (391) lipid molecules, 47,571 (47,559) waters (simple point charge), and 132 (132) NaCl atoms, totaling 222,177 (222,121) atoms with a box size of 129 × 120 × 140 Å (2,167,200 Å3; same dimensions for both the WT and GluN2A-F553A systems) was treated using the OPLS3e force field within the Desmond module of the Schrödinger software package (Desmond Molecular Dynamics System, D. E. Shaw Research). Following the Desmond membrane protein relaxation protocol, two independent 250-ns production runs—one for each receptor (WT and GluN2A-F553A)—were performed at a constant temperature of 310°K and pressure of 1.013 bar with particle mesh Ewald electrostatics using a cutoff of 9 Å for each system. Time-step calculations were performed every 2 fs, and frames were captured every 100 ps. Analysis of the MD trajectory was performed in Desmond and visual MD (VMD; Humphrey et al., 1996). RMSD analysis of the simulations is presented in Fig. S1 (the resulting trajectories are available from Burger and Kell, 2020). A discrepancy is present between the model residue numbering and that of the sequence numbering in the original trajectory files. We have provided a key and useful Python commands at the end of Table S1 to help with renumbering. Additional files with instructions on how to renumber the residues, as well as a TCL script for use in VMD, are also available at the link provided above.
For each of the two trajectories, frames were extracted from the 250-ns simulation at 200-ps intervals. These frames were clustered using the “cluster” function within VMD, an implementation of the quality threshold algorithm (Heyer et al., 1999), with an RMSD-based distance function. PDB files for the frames closest to the center of the largest cluster for each trajectory, which were used in Fig. 6, B–D, are provided in Table S1; the residues in these files are correctly numbered and do not require renumbering. The residues considered for clustering for each simulation include only those shown in Fig. 6, B and D. Additional details and analysis of the model and MD runs are provided in Figs. S2., S3., S4., S5., and S6.
Statistical analysis was conducted using OriginPro 9 or GraphPad Prism 7. Measurements are given as mean ± SEM, with the exception of the concentrations of agonist that produces a half-maximal response (EC50), which are reported as mean with the 95% confidence intervals determined from the log (EC50). For glutamate and glycine potency, no overlap in confidence intervals was considered significant. For the deactivation time course, data were analyzed using one-way ANOVA with Tukey’s post hoc test. Significance was set at P < 0.05 compared with WT, and the number of observations was selected to give power >0.95 for a minimum detectable difference of 50% for weighted τ. When multiple parameters were compared from the same recording (e.g., Fig. 3), we corrected for familywise error.
Online supplemental material
RMSD analysis of the MD simulations is presented in Fig. S1. Fig. S2 provides distributions of distances between pairs of aromatic residues across each of the MD simulations as shown in Fig. 6. Fig. S3 and Fig. S4 show the sequence alignments used for homology model building, along with the sequence identity and similarity between the templates and the target. Ramachandran plots and structure quality metrics for our homology model and the two NMDAR template structures are provided in Fig. S5. An overlay of the model used for MD simulations with a recently published cryo-EM structure of a GluN1/GluN2A receptor is shown in Fig. S6. Table S1 provides the selected model structure and useful Python commands to help with renumbering of model residues.
Sobolevsky et al. (2009) first proposed the idea that the pre-M1 helix could play a role in gating by acting as a cuff, based on the observation that it is in contact with the gate of the channel—the M3 helical bundle that occludes the pore in the closed state of a crystal structure of homomeric GluA2 (Fig. 1). Additionally, human de novo mutations are concentrated in this region (Ogden et al., 2017), suggesting that these residues play a central role in NMDAR gating that cannot be performed by other amino acids. To assess whether the residues of the pre-M1 helix play a role in controlling kinetically distinct pregating steps that reflect conformational changes that precede receptor opening, we designed a set of experiments guided by both functional studies of clinically relevant human mutations and limited structural data. We performed scanning mutagenesis of five residues constituting two turns of the GluN2A pre-M1 helix in which we substituted each residue for alanine and recorded the resultant macroscopic current response time course and single channel unitary currents to identify the contributions of individual residues to overall receptor function.
Effects of pre-M1 mutations on agonist potency and response time course
We expressed WT and mutant GluN2A subunits with the WT GluN1 subunit in Xenopus oocytes and determined the concentration–response relationship for glutamate and glycine. Of the mutants, GluN2A-F553A produced the largest effect on glutamate potency, showing a >10-fold decrease in the EC50 from 3.5 to 0.34 µM. The other mutations tested showed minimal effects on glutamate EC50, with L550A, E551A, P552A, F553Y, and A555P changing glutamate potency by twofold or less (Fig. 2 A and Table 1). These results were mirrored by the findings for glycine potency, whereby GluN2A-F553A increased the glycine potency ∼18-fold, shifting EC50 from 1.4 to 0.077 µM. Other pre-M1 mutations shifted potency by less than twofold (Fig. 2 B and Table 1). These results, where removing the aromatic side chain of Phe553 drastically alters the receptor response to agonist, suggest that this residue in particular may dictate the pre-M1 helix effects on channel gating.
If a conformational change in any part of the pre-M1 helix in response to agonist binding comprises a rate-limiting step that precedes pore dilation, then it might be expected that mutations in this region would alter the NMDAR activation or deactivation time course following rapid application and removal of glutamate. To test this idea, mutant and WT GluN2A subunits were coexpressed with the WT GluN1 subunit in HEK293 cells and whole-cell patch clamp recordings were conducted. HEK293 cells were maintained in solution with saturating glycine, and the response time course following brief coapplication of maximally effective glutamate plus glycine was measured for each GluN2A mutant. None of the mutants tested showed a significant difference in rise time compared with the WT GluN2A receptor. However, analysis of the deactivation time course following 1-s application and rapid removal of glutamate revealed that GluN2A-F553A prolonged the deactivation time course by ninefold, from 46 to 417 ms, whereas GluN2A-L550A accelerated deactivation by 3.5-fold, from 46 to 13 ms. The other GluN2A mutations produced only modest effects on the deactivation time course (Fig. 3 and Table 2). Additionally, GluN2A-L550A and GluN2A-F553A drastically reduced peak amplitude, while the other mutations had no observable effect.
Pre-M1 mutations alter NMDAR single channel properties
To assess the effects of pre-M1 mutations on channel properties, we recorded single ion channel currents in response to maximally effective concentrations of glutamate and glycine from cell-attached patches that contained one active channel (Colquhoun and Hawkes, 1990). Alanine substitutions at several positions caused a wide range of changes in the duration of prolonged shut periods, suggesting an impact on channel gating (Fig. 4, A–F and Table 3). Single channel activity was idealized into a sequence of openings and closings. The shut time histogram was fitted using the sum of four or five exponential components, and the slowest components used to identify a Tcrit with which we could differentiate bursts from long closed states likely to reflect desensitization. This approach allowed for a more accurate interpretation of burst channel data without the potential for confounding results arising from mutation-specific effects on desensitization. When selecting bursts of channel openings based on a critical shut time (Jackson et al., 1983; Magleby and Pallotta, 1983), most of the alanine-substituted channels showed intraburst open probabilities that were similar to WT GluN2A (0.34; Table 3). By contrast, GluN2A-L550A and GluN2A-F553A, which reside on the same side of the helix, decreased open probability nearly 10-fold. Evaluation of the mean open time for each mutation revealed that GluN2A-L550A and GluN2A-F553A reduced the mean open time from 1.49 ms for the WT GluN2A receptor to 0.53 ms and 0.20 ms, respectively (Table 4). Evaluation of the mean shut times showed that both GluN2A-L550A and GluN2A-F553A increased the mean shut time, with GluN2A-L550A producing the largest shift from 29 to 900 ms (Table 5). These results suggest that multiple residues on one side of the short pre-M1 helix impact channel gating and implicate this region as critical for the gating control mechanism.
Previously, WT GluN1/GluN2A receptors have been reported to display multiple gating modes, higher open probability, and a longer mean open time when studied in the absence of extracellular calcium. (Popescu and Auerbach, 2004; Popescu et al., 2004). In this study, the single channel recordings were performed in the presence of 0.5 mM Ca2+. The mean open time and open probabilities reported here for WT GluN1/GluN2A receptors are consistent with those reported previously for the same receptors studied in the presence of extracellular calcium at pH 8.0 (Dravid et al., 2007; Chopra et al., 2015; Gibb et al., 2018).
The pre-M1 helix imparts distinguishing properties to GluN2C and GluN2D
While several residues of the pre-M1 helix are highly conserved across GluN2 subunits, the phenylalanine residue in position 553 is conserved in GluN1, GluN2A, and GluN2B, but corresponds to a tyrosine in GluN2C and GluN2D. Because the aromaticity of the phenylalanine in this position appears to be critical for receptor function, and this residue is the only difference between these subunits in this region, we questioned whether it could be responsible for their differences in open probability. To address this, we recorded single channel currents from GluN2A-F553Y, the GluN2C/D equivalent residue of GluN2A-Phe553. The GluN2A-F553Y receptors showed a slight reduction in glycine EC50 but showed no detectable differences from WT for glutamate EC50 (Fig. 5, A and B) or receptor response time course (Fig. 5 C) during brief application of a maximally effective concentration of glutamate in saturating glycine. Single channel analysis of GluN2A-F553Y revealed a 30% decrease in open probability compared with GluN2A WT (Fig. 5 D and Table 3), suggesting that the tyrosine at this pre-M1 site in GluN2C and GluN2D could contribute to the lower open probability of GluN2C- and GluN2D-containing receptors (Dravid et al., 2008; Vance et al., 2012), perhaps through potential hydrogen bond capability of tyrosine compared with phenylalanine.
MTSEA modification of GluN1-A652C has been shown to lock NMDARs in an open state. Based on the assumption that open probability approaches 1 following covalent modification, we used the reciprocal relationship between the degree of potentiation and open probability to calculate the open probability (see Materials and methods). Introduction of the GluN2A-F553Y mutation was shown to significantly reduce the open probability 55% from 0.22 to 0.12 (Table 6). To support our hypothesis that this specific residue could account for differences in open probability between GluN2A and GluN2D receptors, we measured the changes in agonist potency and open probability of GluN2D-Y578F, the reciprocal substitution to GluN2A-F553Y. Both glutamate and glycine potency were reduced for the mutant receptor with glutamate EC50 shifting from 0.24 to 0.43 µM and glycine EC50 shifting from 0.12 to 0.19 µM (Table 1). Additionally, GluN2D-Y578F increased the open probability 33%, from 0.0055 to 0.0086 (Table 6). Together, these results support our hypothesis that the disparate residue identity at this position between GluN2A/GluN2B and GluN2C/GluN2D contributes in part to the difference in open probability between these subunits.
Subunit-specific pre-M1 helix interactions
A homology model of a GluN1a/GluN2A di-heteromeric structure was used to interpret the results presented thus far. Recently, several GluN1a/GluN2A structures were resolved using x-ray and cryo-EM techniques (Jalali-Yazdi et al., 2018; Zhang et al., 2018); however, none of these structures were available at the start of this study and were therefore not included. Three template structures were used in modeling of the GluN1a/GluN2A NMDAR. Two of the templates were that of the closely related GluN1/GluN2B NMDAR—PDB entries 4PE5 and 5FXH. PDB entry 4PE5 is an x-ray structure resolved at a resolution of 3.96 Å and was selected because it has a well-resolved ABD and ATD. However, there are two limitations of this structure: (1) the linker region between ABD and the TMD is not resolved; and (2), the TMD region only has the backbone residues resolved with no density for the side chains. The second template, PDB entry 5FXH, a cryo-EM structure resolved at 5 Å, was selected because it has a resolved linker region between the ABD and TMD. This template is limited in that only the backbone residues are resolved, but it allowed us to build in the linker region. The third template used, an x-ray PDB entry 5L1B (resolution 4.0 Å), is an AMPAR structure from which the ATD and ABD were removed in the model building process. This template was selected to build in missing residues of GluN2B structures in the TMD. The TMDs of GluN1/GluN2B and AMPA show a 3.5-Å RMSD when aligned on the backbone structures and shares a 30.3% (GluN1) and 27.6% (GluN2) sequence identity within this region. The GluN1 template structures share 92.2% and 94.5% sequence identity with GluN1a target sequence used to build the GluN1a/GluN2A (Fig. S3). The GluN2B template structures share a 64.0% and 69.1% sequence identity with the GluN2A target structure used in model generation (Fig. S4).
The five homology models had molpdf values ranging from 87,071 to 92,375, from which the lowest relative energy model was selected for further analysis. This selected model showed a G-factor of −0.06 (values less than −0.5 are unusual), suggesting good overall stereochemistry of the model. The Ramachandran plot was found to be similar to that of templates PDB identifiers 5FXH and 4PE5 and showed 89.7% of the residues in the most favored regions, 8.7% in additional allowed areas, 0.8% in generously allowed regions, and 0.7% in disallowed regions (none of the latter were found in the pre-M1 helix; total of 3,196 amino acids; Fig. S5). The slightly elevated value for residues in the disallowed region highlights that caution should be taken when interpreting results from the models and subsequent MD simulations even though the value is similar to the observed frequency (0.5%) in the PDB (Gore et al., 2017). These results give us confidence in the quality of the model generated in this study. Finally, we measured the backbone RMSD of the GluN2A model lacking the ATD from a newly released cryo-EM structure (PDB identifier 6IRA; Zhang et al., 2018) which gave a value of 3.7 Å (across ∼4,700 atoms) when aligned on the entire model and had a range between 2.5 and 2.9 Å (across ∼1,500 atoms) when aligned on each chain independently (Fig. S6).
The results presented thus far suggest that the aromaticity of the amino acid residue at position 553 of GluN2A is critical for normal function of GluN2A-containing NMDARs. To explore the mechanism by which this side chain contributes to channel gating, we performed MD simulations over a 250-ns time scale to allow residues in the model to relax to a steady-state position. Analysis revealed a complex network of aromatic residues that were coordinated around GluN2A-Phe553 and included residues Phe810 of the GluN1 pre-M4 linker, Phe641, Tyr645 of the GluN2A M3 transmembrane helix, and Trp558 of the GluN2A M1 transmembrane helix (Fig. 6 B). The simulation places Phe553 near the middle of this aromatic network within a range permissible to both edge-to-face and face-to-face π–π stacking. Specifically, our simulation showed that the centroid of the aromatic ring of GluN2A-Phe553 is a mean distance of 6.6, 5.9, and 5.4 Å from the ring centroids of GluN1-Phe810, GluN2A-Tyr645, and GluN2A-Trp558, respectively. In addition, GluN2A-Phe641 is within 5.7 Å of GluN2A-Phe637 and 5.7 Å from GluN1-Phe817. Histograms showing the full range of distances observed during the 250-ns simulation are provided in Fig. S2. All these reported distances are shorter than the “strict” π–π interaction distance of 6.5 Å for large biomolecules derived from a survey of the entire PDB (Piovesan et al., 2016), with the exception of the mean distance between GluN2A-Phe553 and GluN1-Phe810, which still falls within the “relaxed” threshold of 7.0 Å. Distances between the β-carbons can be found in Fig. S2.
To explore how this aromatic network might be disrupted upon substitution of GluN2A-Phe553 with alanine, we repeated the 250-ns MD simulation with the introduction of this single mutation (Fig. 6 D). Without the presence of the phenylalanine aromatic side chain at position 553, the aromatic network is disrupted and shows increased distances between the Cβ of GluN2A-Trp558, GluN2A-Tyr645, and GluN1-Phe810 from that of the central residue, GluN2A-Ala553, compared with the WT (1.5, 2.5, and 3.2 Å, respectively). We hypothesize that this disruption interferes with communication between the pre-M1 helix and both the M3 helix and the adjacent GluN1 subunit. This finding supports the electrophysiological data in that the receptor’s ability to open or close is altered if interactions between the ligand-binding domain–tethered pre-M1 helix and the channel pore are disrupted.
Many of these aromatic residues are conserved across NMDAR subunits (Fig. 1 B), implying that the interactions between these residues is similarly conserved. However, the GluN1 pre-M1 helix appears to participate in gating in a manner distinct from that of GluN2A, as suggested by results showing that mutations in the GluN1 pre-M1 helix do not slow receptor activation, while their GluN2 counterparts do (Ogden et al., 2017; Gibb et al., 2018). We therefore investigated the amino acid environment surrounding the GluN1-Phe558 residue (Fig. 6 C). We found that GluN1-Trp563 is within 7.4, 6.2, and 5.8 Å of GluN1-Tyr647, GluN1-Phe554, and GluN1-Phe558, respectively. GluN1-Tyr647 corresponds to GluN2A-Tyr645 and GluN1-Trp563 corresponds to GluN2A-Trp558. Additionally, the GluN1 equivalent of GluN2-Phe651 and the GluN2 equivalent of GluN1-Phe810 are both replaced by isoleucine, which is incapable of participating in the aromatic network. Comparison of the set of GluN1-Phe558 and GluN2A-Phe553 aromatic interactions might explain why disease-associated mutations in GluN2A do not have the same functional consequences when expressed in GluN1 (Ogden et al., 2017; Gibb et al., 2018). Because the aromatic network that potentially dictates channel gating is less extensive in GluN1 pre-M1 helix, mutations in this region might have reduced functional consequences.
In this study, we investigated the role of the GluN2A pre-M1 helix in channel gating. We found the GluN2A residue Phe553 to be particularly critical for channel function as indicated by significant changes to both macroscopic and single channel properties. Moreover, we used MD simulations to probe the mechanism by which this residue might contribute to gating. Our analysis revealed a network of aromatic amino acids within the π-stacking range surrounding GluN2A-Phe553 that was disrupted by alanine substitution at this position.
Specifically, we conducted scanning mutagenesis of the pre-M1 helix to investigate the effects on various receptor properties. Of the mutants tested, GluN2A-F553A presented the greatest effect on glutamate and glycine EC50, showing an ∼11-fold and an ∼18-fold increase in potency, respectively. That is, GluN2A-F553A shifted the potency of both NMDAR agonists, despite only being introduced to the GluN2A subunit. In addition, GluN2A-F553A prolonged the weighted deactivation time course by approximately ninefold, reduced peak amplitude by ∼16-fold, decreased open probability by ∼10-fold, and reduced the mean open time by approximately sevenfold. GluN2A-L550A, situated on the same side of the pre-M1 helix as GluN2A-Phe553, also altered single channel properties, decreasing open probability by ∼10-fold, increasing the mean open time by approximately threefold, and increasing the mean closed time by ∼30-fold. Together, these results suggest that these two amino acids in GluN2A (Leu550 and Phe553) control channel gating in a manner that is likely dependent upon their structural contacts with residues adjacent to their side of the pre-M1 helix. Additionally, these results may also provide mechanistic context for the previously characterized disease-associated mutations GluN1-P557R, GluN2A-P552R, and GluN2B-P553L in this region (Ogden et al., 2017).
Our single channel recordings were performed in 0.5 mM Ca2+, which can have a profound effect on single channel properties. The open probability we report here for the WT GluN2A receptor (0.34) is consistent with those determined in the presence of Ca2+ in cell-attached patches (Dravid et al., 2007; Iacobucci and Popescu, 2017) and in outside-out patches (Gibb et al., 2018; Rycroft and Gibb, 2002; Bhattacharya et al., 2018). Moreover, Iacobucci and Popescu (2017) show the open probability drops by threefold from 0.62 (0 mM Ca2+) to 0.21 (0.5 mM Ca2+), suggesting that the lack of divalent ions in previously published cell-attached patch studies contributes to the higher open probability observed in the absence of extracellular Ca2+ (Popescu et al., 2004).
GluN2A-Phe553 has previously been investigated for its role in desensitization. In AMPARs, a leucine resides at this position, perhaps suggesting that the amino acid at this position contributes to the differences in desensitization between these two families of ionotropic glutamate receptors. When GluN1-Phe558—the equivalent of GluN2A-Phe553—was substituted with leucine and expressed with WT GluN2A, the receptor retained its GluN2A-like desensitization phenotype. However, when GluN2A-F553L or GluN2B-F554L was expressed with WT GluN1, the desensitization resembled that of the AMPAR (Alsaloum et al., 2016). It was proposed, therefore, that the mechanism of desensitization in NMDARs differs from that of AMPARs as a result of differences between the composition of a hydrophobic box made up of aromatic residues within pre-M1 and M3. While this finding emphasizes the importance of receptor-specific hydrophobic boxes in functional divergence between the ionotropic glutamate receptors, it does not address the functional differences between GluN1 and GluN2 subunits suggested by the finding that the homologous F553L mutation in GluN1 was not sufficient to alter desensitization.
Interestingly, when GluN2A-Phe553 was substituted with tyrosine, the GluN2C and GluN2D amino acid equivalent, we found no measurable effect on desensitization as previously reported (Krupp et al., 1998; Alsaloum et al., 2016), suggesting that the pre-M1 segment cannot account for differences in desensitization between different GluN2 subunits. Our results show that, while GluN2A-F553Y had no significant effect on macroscopic properties, this mutation significantly reduced open probability, consistent with the reduced open probability of GluN2C- and GluN2D-containing receptors. We found that GluN2D-Y578F reduced glutamate potency, reduced glycine potency, and increased open probability as determined by MTSEA relative to the WT GluN2D receptor. These findings suggest that, in addition to being critical for the gating mechanism, the identity of the residue at position 553 can impart subunit-specific single channel properties to the receptor. The hydroxyl group of the tyrosine in this position for GluN2C and GluN2D receptors may introduce a hydrogen bond that is absent from GluN2A and GluN2B receptors. One possibility for such a hydrogen bond interaction, as predicted from sidechain positions in the MD simulations, could be between the hydroxyl group of the tyrosine at this site in the pre-M1 and the tyrosine of the SYTANLAAF motif (Y645 in GluN2A; Fig. 6 B). Additionally, the GluN2A-Ala555 residue is not conserved across the NMDA subunits. In GluN2C and GluN2D, the residue in this position is a proline, as tested in our experiments. However, as our data show, the GluN2A-A555P mutation was not sufficient on its own to produce significant changes to macroscopic or single channel properties of the receptor, suggesting that pre-M1 control of single channel properties is specific to the residue at position 553.
The side of the helix that harbors GluN2A-Leu550 and GluN2A-Phe553 appears to be oriented toward the SYTANLAAF-conserved helical bundle and the pre-M4 helix of the adjacent subunit (Chen et al., 2017). That is, the residues that were of consequence during our scanning mutagenesis study are positioned toward the triad that has been proposed to control channel gating (Gibb et al., 2018). Our findings, taken with the results from Alsaloum et al. (2016), suggest that the phenylalanine of the pre-M1 helix, despite being critical for channel gating and desensitization, differs in function among GluN subunits, likely as a result of differences in local residues with which this region can interact.
We used MD to explore the hydrophobic network in the context of pre-M1 scanning mutagenesis to identify potential roles for pre-M1 that could account for differences in gating control between GluN1 and GluN2. Using a GluN2A homology model built from a GluN2B crystal and cryo-EM structures, as well as the TMD of a closed AMPA structure, we examined the amino acid environment surrounding the pre-M1 helix. Within the GluN2A subunit, Phe553 is surrounded by a network of aromatic amino acids composed of GluN2A-Phe641 and GluN2A-Tyr645 of the GluN2A M3 helix, GluN2A-Trp558 of the GluN2A M1 helix, and GluN1-Phe810 of the neighboring GluN1 pre-M4 helix. When we repeated the simulation with an alanine substituted at the position of GluN2A-Phe553, we observed a disruption of the aromatic network and an increase in the distances from the ring centroids of GluN2A-Tyr645 and GluN2A-Trp558 to the 553 position. This supports the hypothesis that the aromatic network is critical for channel function and disruption of this network can perturb channel function. Finally, a comparison of the GluN1 aromatic network with that observed for GluN2A revealed a reduction in the number of interacting residues in GluN1, suggesting that the network plays a unique role in the GluN2A subunit.
Previous findings show that some of the residues close enough to contact GluN2A-Phe553 also influence channel function and gating, which lends credibility to the idea that the aromatic amino acid network controls gating in the NMDAR. First, the disease-associated mutation GluN2A-W558S was identified in a patient with epilepsy (ClinVar; available from https://www.ncbi.nlm.nih.gov/clinvar). Within the M3 helix, GluN2A-Phe637 has been shown to influence agonist potency and channel gating (Ren et al., 2007), and GluN1-Y647S has been identified in a patient with infantile spasms (Allen et al., 2013). In the M4 helix, GluN1-F817L has been identified in a patient with intellectual disability, developmental delay, and movement disorder (Lemke et al., 2016; Zhu et al., 2015). Although these disease-associated mutations have not been functionally characterized, their potential role in these neurological disorders suggest that they likely disrupt normal NMDAR function.
In 2018, 23 full-length structures of the heteromeric GluN2A receptor were resolved (4.5–16.5 Å) using cryo-EM, none of which had resolution covering the complete transmembrane and linker regions (Zhang et al., 2018; Jalali-Yazdi et al., 2018). Some structures showed resolution within the pre-M1 helix region; however, they showed high B-factor values (358–515). When we compared the two GluN2A chains, Phe553 and Leu550 had different positions on the chains (B and D), suggesting uncertainty within the structures. We therefore opted to use homology models and not the cryo-EM GluN2A structures for calculations. Future studies however could benefit from incorporating the cryo-EM GluN2A structures in conjunction with AMPAR structures as templates to build homology models. The homology models generated were from crystal and cryo-EM structures that had resolution within the pre-M1 and linker regions. We used the well-resolved transmembrane region of the AMPAR, which has high sequence homology and structural similarity to the transmembrane regions of GluN1 and GluN2A.
Overall, the results from this study emphasize the previously suggested role of the pre-M1 helix in channel gating by demonstrating how specific residues within this region are critical for channel function. In addition, these results implicate the surrounding amino acid environment in the mechanism of channel gating by revealing a network of aromatic residues that is disrupted by an alanine substitution at the central phenylalanine residue. Finally, these results suggest unique subunit-specific contributions of pre-M1 helices in GluN1 and GluN2 to channel gating due to nonequivalent gating triads.
Merritt C. Maduke served as editor.
This work was supported by the National Institute of Neurological Disorders and Stroke, National Institutes of Health (grants R01-NS036654 and R35-NS111619-01 to S.F. Traynelis and grant F31-NS106896 to M.J. McDaniel).
S.F. Traynelis is a primary investigator on research grants from Biogen and Janssen to Emory University School of Medicine; is a member of the Scientific Advisory Board for Sage Therapeutics, the GRIN2B Foundation, and the CureGRIN Foundation; is cofounder of NeurOp, Inc.; and receives royalties for software. D.C. Liotta is a member of the Board of Directors for NeurOp, Inc. D.C. Liotta and S.F. Traynelis are coinventors of Emory-owned intellectual property that includes allosteric modulators of NMDA receptor function. P.B. Burger is a cofounder of Avicenna Biosciences, LLC.
Author contributions: M.J. McDaniel designed and performed experiments, performed analysis, and wrote the paper. K.K. Ogden designed and performed experiments, performed analysis, and wrote the paper. S.A. Kell designed experiments, performed simulations, performed analysis, and wrote the paper. P.B. Burger designed experiments, performed simulations, performed analysis, and wrote the paper. D.C. Liotta designed experiments and wrote the paper. S.F. Traynelis designed experiments, performed analysis, and wrote the paper.