Structural insights into phenylethanolamines high-affinity binding site in NR2B from binding and molecular modeling studies
Molecular Brain volume 1, Article number: 16 (2008)
Phenylethanolamines selectively bind to NR2B subunit-containing N-methyl-D-aspartate-subtype of ionotropic glutamate receptors and negatively modulate receptor activity. To investigate the structural and functional properties of the ifenprodil binding domain on the NR2B protein, we have purified a soluble recombinant rat NR2B protein fragment comprising the first ~400 amino acid amino-terminal domain (ATD2B) expressed in E. coli. Spectral measurements on refolded ATD2B protein demonstrated specific binding to ifenprodil. We have used site-directed mutagenesis, circular dichroism spectroscopy and molecular modeling to obtain structural information on the interactions between critical amino acid residues and ifenprodil of our soluble refolded ATD2B proteins. Ligand-induced changes in protein structure were inferred from changes in the circular dichroism spectrum, and the concentration dependence of these changes was used to determine binding constants for ifenprodil and its analogues.
Ligand binding of ifenprodil, RO25,6981 and haloperidol on soluble recombinant ATD2B determined from circular dichroism spectroscopy yielded low-to-high micromolar equilibrium constants which concurred with functional IC50 measurement determined in heterologously expressed NR1/NR2B receptors in Xenopus oocytes. Amino acid residue substitutions of Asp101, Ile150 and Phe176 with alanine residue within the ATD2B protein altered the recombinant protein dissociation constants for ifenprodil, mirroring the pattern of their functional phenotypes. Molecular modeling of ATD2B as a clam-shell-like structure places these critical residues near a putative ligand binding site.
We report for the first time biochemical measurements show that the functional measurements actually reflect binding to the ATD of NR2B subunit. Insights gained from this study help advance the theory that ifenprodil is a ligand for the ATD of NR2B subunit.
N-methyl-D-aspartate (NMDA) receptors are postsynaptic receptors for L-glutamate, the major excitatory neurotransmitter in the central nervous system. They are ligand-gated ion channels that are permeable to K+, Na+ and Ca2+. They play a pivotal role in neuronal development, synaptic plasticity, and learning and memory [1–5]. However, the overactivation of these receptors has been implicated in excitotoxic neuronal cell death observed in neurological disorders such as stroke, epilepsy [6, 7] and head trauma .
NMDA receptors are heteromeric structures composed of an obligatory NR1 subunit and one or more of the four NR2 subunits (NR2A-D) and/or NR3 subunits (NR3A-B) [9–12]. Each NMDA receptor subunit shares a characteristic modular architecture with three transmembrane segments (M1, M3 and M4), a re-entrant loop (initially known as M2) which forms the pore-lining region, an intracellular C-terminal domain, and large extracellular domains [5, 10, 12]. The large extracellular domains of the mature NR2B subunit include amino acids 27–557 (Accession number U11419) situated pre-M1 and amino acids 649–817 between M3 and M4. The glutamate binding pocket is made up of two discontinuous stretches of amino acids encoded by the S1 and S2 regions of cDNA. Each segment constitutes ~150 amino acids preceding M1 and following M3. The remaining non-agonist binding segment preceding M1 is also known as the amino-terminal domain (ATD) (Fig. 1A).
Expression of isolated domains has been applied to ligand-gated glutamate receptors, voltage-gated potassium ion channels and many other plasma membrane bound ion channels and receptors with reasonable success. For NMDA receptors, the crystal structures of the ligand binding core (S1S2) for glycine (NR1) and glutamate (NR2A) have been successfully solved recently [13, 14]. These atomic structures provide, for the first time, a bridge between data obtained from biophysical, molecular and biochemical studies of NMDA receptors. There are now over 50 crystal structures of various glutamate receptor binding cores bound to different ligands [15–19].
In contrast, there is no crystal structure for the putative allosteric modulators' binding sites located on the amino terminal domain proximal to the S1 region of ionotropic glutamate receptors, including NMDA receptors. The apparent high affinity Zn2+ binding site is situated in ATD of NR2A [20–22]; the putative phenylethanolamines binding site(s) in ATD of NR2B [23–25]; the polyamine-sensitive region in ATD of NR1 [26, 27]; and possibly a Zn2+ binding site in NR2B . Indeed, Rosahl and colleagues  provided ex vivo support for the modulatory role of ATD of NR2B by ifenprodil and ifenprodil-like NMDA receptor antagonists recently. By replacing the mouse NR2B gene function by "knocking-in" a chimeric human NR2A/B cDNA containing ATD of NR2A fused to the endogenous NR2B locus, the high affinity component from the inhibition of NR2B containing receptors was absent in culture hippocampal neurons.
Like the ATD of NR2A subunit of NMDA receptor which binds the negative allosteric modulator Zn2+ [20–22, 29], NR2B ATD binds negative modulator phenylethanolamines and analogues (ifenprodil, RO25,6981, haloperidol) (Fig. 1B) [23, 24, 30–34]. The NMDA receptor subunit ATD is a large domain made up of the first ~400 amino acids (e.g. residues 27–403) of the mature polypeptide. The ATD has weak homology and shows similarities in predicted secondary structure with bacterial periplasmic leucine/isoleucine/valine-binding protein (LIVBP) [27, 35]. NMDA receptor ATD also shares homology and predicted secondary structure with the G-protein-coupled metabotropic glutamate receptor 1 (mGluR1), for which the crystal structure of the glutamate binding domain is known . There is no crystallographic structural information on the ATD of NMDA receptor or any other ionotropic glutamate receptors at atomic resolution. This has slowed advances in our understanding about how the negative modulators inhibit the opening of the NMDA receptor channel, which is triggered by glutamate binding to the agonist binding core . However, several low level electron microscopy studies suggest an orientation of ATD with respect to the glutamate binding core in AMPA receptors [38–40].
To begin to address this gap in our understanding, we have constructed and expressed the modular amino-terminal domain (proximal 367 amino acids) of NR2B subunit (ATD2B) as a 6xHis-fusion protein (6 × His-ATD2B) in Escherichia coli expression system . In this study, we demonstrate, for the first time, that soluble 6 × His-ATD2B recombinant proteins can bind various modulators with dissociation constants (KDs) close to and correlated tightly with the functional IC50 measurements. Mutations at the critical amino acids (Asp101, Ile150, Phe176) within ATD of NR2B subunit significantly decrease the binding affinity of ifenprodil to 6 × His-ATD2B. Three-dimensional homology model of the ATD of NR2B revealed that these critical amino acids are strategically positioned around the putative ligand binding cleft. Taken together, our multi-disciplinary approach supports the concept of a modular organization of ATD within the NR2B subunit as well as its inherent ability to bind negative allosteric modulating ligands.
Buffers -9, -12 and -13 refolded ATD of NR2B
Using a fractional factorial folding screen [41, 42], we had successfully refolded 6 × His-ATD2B recombinant protein in three buffers differing in pH, refolding temperature, divalent ions, polar and nonpolar additives, chaotrope, GSH:GSSH ratio, detergent, PEG, and presence/absence of a ligand. These buffers were described in our previous work [see Additional file 1], and are referred to here by the same nomenclature as that study (Buffers -9, -12 and -13 ). We had shown previously Buffer-13 refolded 6 × His-ATD2B existed in a soluble and monodispersed species which bound ifenprodil in a concentration-dependent manner (KD = 128 nM) .
In this study we examine the structural and functional properties of refolded recombinant 6 × His-ATD2B proteins. Light scattering detection analysis allows biomolecular sizing and characterization of proteins to be measured directly. Dynamic light scattering analyses showed Buffers -9 and -12 refolded 6 × His-ATD2B proteins existed in monodispersed species (data not shown). Their relative percentages of α-helix content were estimated to be 26.7 ± 0.4% (n = 5) and 27.7 ± 0.7% (n = 6), respectively. These values were congruent with the 26.8% estimated for Buffer-13 refolded protein reported previously .
Circular dichroism (CD) spectroscopy is a useful technique for studying protein-protein interactions in solution . Circular dichroism spectroscopy measures differences in the absorption of left-handed polarized light versus right-handed polarized light. Units given are typically ellipticity (theta, θ), which is related to the difference in absorbance of left and right circularly polarized light. This difference arises from the interaction between the light and the chiral nature of the protein structure. The chromophore for wavelengths in the "far-UV" spectral region (190–250 nm) is the peptide bond, and the signal arises when it is located in a regular, folded environment. Secondary structural elements such as alpha-helix, beta-sheet, and random coil structures each give rise to a characteristic shape and magnitude of CD spectrum. The CD spectrum of a protein in the "near-UV" spectral region (250–350 nm) can be sensitive to the overall tertiary structure of the protein, including signals for phenylalanine (250–270 nm), tyrosine (270–290 nm), and tryptophan (280–300 nm), as well as weak signals from disulfide bonds throughout the near-UV spectrum. Thus, evaluation of changes in the CD spectrum can be used to monitor changes in protein conformations following ligand binding, and the concentration dependence of ligand-induced changes in the CD spectrum can also be used to estimate binding affinities [44–46].
CD spectra of the Buffers -9 (Fig. 2A) and -12 refolded (Fig. 2B) 6 × His-ATD2B proteins were obtained in the presence of 5 μM ifenprodil and 0.5 μM RO25,6981. Ifenprodil and RO25,6981 in the concentration ranges 0.01–5 μM and 0.002–0.5 μM, respectively, diminished the θ over the wavelengths in the region 210–250 nm. We interpret the concentration-dependent shifts of magnitude at wavelength 220.0 nm ellipticity to suggest that ifenprodil binds to refolded soluble 6 × His-ATD2B proteins with KD values of 60 ± 18 nM for Buffer-9 (n = 5, Fig. 3A) and 61 ± 34 nM for Buffer-12 (n = 5, Fig. 3B). The concentration-dependent changes in the CD spectrum induced by RO25,6981 (218.9 nm) upon interacting with the 6 × His-ATD2B proteins suggest that the ligand binds with KD values of 7 ± 3 nM (Buffer-9), 4 ± 1 nM (Buffer-12) (n = 5, Figs. 3C, D) and 4 ± 2 nM (Buffer-13). Figure 4A summarizes the binding properties of 6 × His-ATD2B to ifenprodil and RO25,6981 between three different protein refolding buffers. In addition, in the presence of saturating concentration of 5 μM ifenprodil, RO25,6891 could displace ifenprodil and bind to 6 × His-ATD2B with comparable affinities in the absence of ifenprodil: Buffer -9 (4 ± 1 nM), Buffer -12 (2 ± 1 nM) and Buffer -13 (3 ± 1 nM).
Ligand Inhibition of Heterologously Expressed NR1/NR2B Receptors in Xenopus Oocytes
Ifenprodil and RO25,6981 are NR2B-selective antagonists. Although both molecules share similar chemical structures that are likely to act on the same putative binding site residing in the ATD of NR2B, RO25,6981 is more potent than ifenprodil in inhibiting heterologously expressed and native NR2B-containing NMDA receptors (see  for review). We heterologously expressed NR1 with NR2B in Xenopus oocytes to determine the NMDA receptor sensitivity to both ifenprodil and RO25,6981 so as to establish the ratio of IC50(ifenprodil) to IC50(RO25,6981). This ratio will determine which buffer-refolded 6 × His-ATD2B (ie. Buffer-9 vs Buffer-12 vs Buffer-13) will be used for detail binding analyses. The negative allosteric antagonists inhibited NR1/NR2B receptor current responses with IC50 values 105 nM and 59 nM at pH7.3, respectively (Fig. 4B; Table 1). The ratio of IC50(ifenprodil) to IC50(RO25,6981) is almost 2.
6 × His-ATD2B Binds Antipsychotic Drug Haloperidol
6 × His-ATD2B refolded in Buffer-9 was able to bind a number of ligands with a wide range of affinities. Haloperidol, a dopaminergic antagonist with antipsychotic properties, is also known to block NR2B-containing NMDA receptors at micromolar concentrations [47–49]. Analysis of changes in CD spectra suggested that haloperidol bound to Buffer-9 refolded 6 × His-ATD2B in a dose-dependent manner. The maximum magnitude shift in ellipticity over 0.1–30 μM of haloperidol at 219.6 nm yielded a KD of 3094 nM (Table 1).
If the KD values estimated for the various ligands tested accurately represent dissociation equilibrium constants, then the biochemically determined KD values should correlate well with the IC50 values determined in functional assays. To address this, we performed linear regression analysis of three ligands with wide apart affinities (60–440 fold, Table 1; also see  and ) on our experimentally determined IC50 and KD. The ability of our Buffer-9 refolded 6 × His-ATD2B protein to bind three ligands of widely varying affinities (ifenprodil, RO25,6981 and haloperidol) yielded a tight correlation (R = 0.96; Fig. 4C). This strong linear correlation provided additional evidence that our soluble recombinant protein assumed a conformation that could bind its specific ligands.
Computing the ratios of KD(ifenprodil) to KD(RO25,6981) for each refolding buffers showed approximately 9-fold for Buffer-9, 15-fold for Buffer-12 and 24-fold for Buffer-13 refolded 6 × His-ATD2B proteins (Fig. 4A; Table 1). Buffer-9 refolded 6 × His-ATD2B yielded a ratio closest to the relative potency (IC50 ratio of 2) of ifenprodil to RO25,6981 (Fig. 4B and Table 1). Thus, subsequent structural and functional characterizations of the modular ATD of NR2B were performed on Buffer-9 refolded 6 × His-ATD2B protein.
Impaired Binding Interactions between Recombinant ATD2B Mutants and Ifenprodil
We next sought to more directly evaluate the ifenprodil binding pocket of Buffer-9 refolded 6 × His-ATD2B. A charged (Asp101), an aliphatic (Ile150) and an aromatic (Phe176) amino acid residues, previously reported to be critical molecular determinants for ifenprodil inhibition of NR2B-containing NMDA receptors (> 60-fold increase in IC50 values compared to wild-type NR1/NR2B receptors; ), were mutated to alanine (designated as ATD2B(D101A), ATD2B(I150A), ATD2B(F176A)) using site-directed mutagenesis (see Methods). We subjected these mutant 6 × His-ATD2B proteins to increasing concentrations of ifenprodil to evaluate the concentration dependence of binding. In the presence of increasing amounts of ifenprodil, each mutant protein D101A and I150A yielded higher KD values of 173 nM and 139 nM, respectively (Table 2). By substituting phenylalanine with alanine at amino acid 176, the concentration-effect curve yielded even higher KD value of 331 nM (Fig. 5A solid circle). When compared to the wild-type (WT) ATD2B protein, these mutations caused a 2.9- (D101A), 2.3- (I150A) and 5.5- (F176A) fold increase in KD values, which were statistically significant (P < 0.05) (Fig. 5C and Table 2).
We hypothesized that ifenprodil binding would destabilize further if more than one of the several key amino acid residues, whose side chain may make Van der Waals or hydrogen bonding contact with the ligand, were removed. We selected the charged Asp101 and hydrophobic Phe176 residues to disrupt jointly in our recombinant ATD2B protein. D101A was inserted into ATD2B(F176A) plasmid using overlap extension site-directed mutagenesis (see Methods). Indeed, this double mutant recombinant protein ATD2B(D101A/F176A) displayed a much weaker binding to ifenprodil (KD 627 nM compared to 60 nM WT; P < 0.05; Fig. 5A, C; Table 2).
Molecular Modeling of the Amino-Terminal Domain of the NR2B Subunit
To better understand the structural determinants of ligand binding to ATD of NR2B, we used Prime Structure Prediction (Schrödinger) to build a homology model of ATD2B based on the crystal structure of mGluR1 agonist binding domain complexed with glutamate (1 EWK; ). The homology model of the ATD of NR2B in the apo (absence of ligand) configuration is depicted in Fig. 6A. It is built using the sequence alignment as shown in Additional file 2D. The α-helices and β-sheets were labeled according to the mGluR1 structure . Three (Asp101, Ile150 and Phe176) out of four amino acids (Phe182) reported by Perin-Dureau and colleagues  that played a critical role in the inhibition of NR2B-containing NMDA receptors (IC50(mutant)/IC50(WT) > 60) are highlighted in Fig. 6A. Examinations of these three "critical" molecular determinants (D101A, I150A and F176A) reveal these residues to be spatially distributed near a central cleft despite their locations on different regions of the bi-lobed structure (Asp101 on lobe I while Phe176 on lobe II) as well as on the hinge region (Ile150) of the NR2B ATD (see also Fig. 8A in Perin-Dureau et al. ). Figures 6B and 6C show these three critical residues to be located within the putative central cleft and in contact with the top-scoring pose  of ifenprodil. In fact, ifenprodil is predicted to coordinate with Asp101 and Phe176 via hydrogen and hydrophobic interactions. In addition, the docked pose has significant Van der Waals contact with Val42, Thr103 and Lys234, residues that all had a moderate effect in previous studies  (data not shown).
Finally we searched for an amino acid residue located outside the cleft of the clam-shell conformation of NR2B ATD. A charged residue, Asp113 located on lobe I, was mutated and probed for its ligand binding phenotype (see Fig. 6A). The D113A mutation showed nearly superimposable dose-dependent ifenprodil titration curve with a KD similar to the wild-type soluble 6 × His-ATD2B (P > 0.4; Figs. 5B, C and Table 2).
In the present work, we have used an autonomously folded protein fragment of the amino terminal domain of NR2B subunit to elucidate several critical residues that influence the binding of putative modulators on the mammalian NR2B subunit-containing NMDA receptors. The two most important conclusions to emerge from this study are that (1) the isolated ATD2B protein is sufficient to interact with various modulators with affinities that correlate well with the physiological and biochemical studies, and (2) the critical molecular determinants for the binding of these modulators reside within the proximal ~400 amino acids preceding S1 domain of NR2B subunit. These data provide direct confirmation of ligand binding to ATD as well as suggest that mutations previously described influence the binding of ifenprodil, a NR2B-selective modulatory ligand.
Comparison to previous studies
The KD values determined with CD spectroscopy of ligands for ATD2B correlated with functionally determined IC50 values (Fig. 4C), and were consistent with KD values previously reported for [3H]ifenprodil binding to rat brain cortex (KD = 36 nM)  and both NR1/NR2B transfected HEK293 and mouse L(tk-) cell line (KD = 34–79 nM) [52, 53]. Furthermore, our ligand displacement experiment using RO25,6981 (~9-fold higher KD affinity determined in this study) over ifenprodil on ATD2B protein yielded Ki values (2–4 nM for Buffers -9, -12 and -13, n = 5 per buffer; details not shown) comparable to KD values that were determined for RO25,6981 alone. Thus, our findings on refolded ATD2B agree well with previous determinations of ligand affinities determined in intact receptors.
Previous description of alanine mutagenesis scan in the NR2B ATD on key amino acid residues known to coordinate ifenprodil and RO25,6981 binding [23, 24] corroborated our structural findings. Perin-Dureau and colleagues  showed that the substitutions of Asp101, Phe176 and Phe182 produced the largest effects on ifenprodil potency for inhibition of NR1/NR2B currents, and hypothesized that these residues were involved in ifenprodil binding. They further proposed that Ile150 participated in the clam-shell closure mechanism rather than direct contact with the ligand. D101A and F176A mutations inserted individually into the full length NR2B subunit of the heterologously expressed NR1/NR2B receptors yielded > 60-fold increase in IC50 values . Our experimentally determined KD values for these two same mutations inserted into our ATD2B refolded protein were 2.3–5.5 fold higher than those determined from wild-type ATD2B. Our results are consistent with the ideas presented by Perin-Dureau et al. , and suggest that these residues can directly impact ifenprodil binding. However, the single residue mutant ATD2B proteins did not generate the same degree of shift in KD as observed for the IC50 value determined from two electrode voltage clamp current measurements. The IC50 values determined from the voltage clamp experiments describe the concentration-dependence of ifenprodil reduction in cation flux through the aqueous channel pore. This functional measure thus takes into account the effects of mutations on allosteric interactions between domains, long range intra-protein interactions, in addition to effects on ligand binding site. It seems likely that some of the differences in our KD and functional IC50 values are attributed to the study of ATD2B in isolation.
Our KD values, determined via dose-dependent ligand-binding assay using CD, are 'snap-shots' physical constants between the ligands and the isolated ATD protein. Interestingly, D101A and F176A mutations resulted in complete loss of [3H]RO25,6981 binding to membranes of HEK293 cells transiently co-transfected with NR1 and NR2B cDNAs . Since both ifenprodil and RO25,6981 only differ chemically by a methylene between the piperidine ring and the phenyl ring (see Fig. 1B), one plausible scenario would be that other surrounding charged (eg. Asp104) and hydrophobic (eg. Phe182) residues could possibly compensate partially for the substituted Asp101 and Phe176, respectively . Alternatively, both ligands may dock themselves in an opposite manner into the binding pocket within ATD2B. Notwithstanding these discrepancies, our experimental data on double-mutant ATD2B(D101A/F176A) recombinant protein, as predicted, yielded larger shift in its KD. The degree of this reduction in affinity for the double-mutant protein in association with ifenprodil is more than an additive effect when compared to single amino acid substitutions of D101A or F176A.
Molecular features of ifenprodil binding and inhibition
Our homology model shows a clam-shell configuration of NR2B ATD with a central cleft. The carboxyl of Asp101 and the hydrophobic phenyl ring of Phe176, from lobe I and lobe II respectively, are located within the microenvironment of the cleft (Fig. 6B). The carboxyl group of Asp101 can potentially interact with the basic nitrogen on the piperidine ring, or either of the hydroxyl groups of ifenprodil via electrostatic interactions or hydrogen bonding. The phenyl rings of Phe176 and ifenprodil can establish hydrophobic interaction (Fig. 6C). The predicted distance between these two residues from our model is 11.2 Å (γ-carbon to γ-carbon). In the homology model, Ile150 is located at the end of one of three flexible loops that link Lobe I and Lobe II of the ATD. Together with other critical residues such as Phe176 and Phe182, its alkyl side chain helps form a hydrophobic backing to the pocket where ifenprodil and RO25,6981 are predicted to dock (Fig. 6C). As such, it is plausible to hypothesize that Ile150 and Phe182 potentially provide contact to ligands through nonspecific lipophilic interactions . Ile150 and Phe182 may also provide a hydrophobic environment that helps organize other nearby residues such as Phe231 or Val258, or they may provide a limit on the length of the ligand (the phenol end) as McCauley and colleagues observed of large non-polar hydrophobic side chains larger than methyl at the phenol end of ifenprodil-like compounds, which have poor inhibitory properties on NR1/NR2B expressing cell lines . In contrast to elongated binding poses with the ligand stretched across the cleft as seen in Marinelli et al. , the formation of a deep pocket allows the ligand to bury some of its hydrophobic regions as well as bringing it into proximity with other residues shown to be critical for binding.
We have provided evidence for the binding of ifenprodil and its analogues directly to soluble recombinant ATD protein of NR2B subunit. We have also provided data supporting the direct involvement of intra-cleft residues predicted by a homology model of ATD2B in ifenprodil binding and developed a proposed binding mode that more clearly demonstrates the role of these residues. These show that the soluble folded ATD2B protein provides a new tool for structure-activity studies on existing and new subunit selective modulators that bind to glutamate receptor amino terminal domains .
Rat cDNAs for NR1-1a (hereafter NR1) and NR2 subunits (GenBank numbers NR1 U11418 and U08261; NR2B U11419) were provided by Dr. S. Heinemann (Salk Institute, La Jolla CA). The vector pET30b(+)EG was a gift from Dr E. Gouaux (Vollum Institute, Portland OR). Ifenprodil, RO25,6981 and haloperidol were obtained from Sigma. Oligonucleotides were purchased from Metabion (USA) and Proligo (Singapore).
Residues that are targeted by mutagenesis in ATD sequence of NR2B subunit are shown in Additional file 2A. These amino acids were mutated in pET30b(+)EG-ATD2B plasmid  by overlap extension site-directed mutagenesis (D101A, I150A, F176A) [56, 57] and QuikChange methods (D113A) . The mutagenesis primers used are listed in Additional file 3. For overlap extension mutagenesis, two addition primers (2B-PRM5 and 2B-PRM6)  were coupled to the above mutant primers to first generate two amplicons with overlapping ends. The authenticity of the full-length cloned inserts in all recombinant constructs was confirmed with dideoxy DNA sequencing in both directions.
Wild-type and Mutant 6 × His-ATD2B Protein Expression, Solubilization, Purification and Refolding
The wild-type and mutant 6 × His-ATD2B recombinant proteins were expressed, lyzed, solubilized, affinity-purified and refolded as described previously  [Additional file 2B, C]). All proteins were finally dialyzed against Buffer 2 [see Additional file 1]. The protein concentrations of purified 6 × His-ATD2B were determined by Bradford protein assay and absorbance at an extinction coefficient (ε280 = 60,430 M-1cm-1) as reported previously .
Dynamic Light Scattering (DLS), Circular Dichroism (CD) Spectroscopy and Ligand-Binding Assay
DLS measurements were made with DynaPro (Protein Solutions) and CD spectra were measured at 25°C on a Jasco J-715 Spectropolarimeter equipped with a 0.1 cm path-length cuvette (Jasco Corporation, Tokyo, Japan) as described previously . For the CD assays, the refolded protein was concentrated to 1 mg/ml using Vivaspin 2 ml concentrator (Vivascience). Ligand (haloperidol, ifenprodil or RO25,6981) was added in small aliquots from stock solution to 350 μl of the protein solution. The final change of assay volume was < 3% for all ligands investigated. The CD spectra of the wild-type and mutant 6 × His-ATD2B proteins in the absence and presence of a ligand were recorded from 190 to 250 nm. The CD spectrum of Buffer 2 was subtracted from all the CD spectra obtained. To ensure that the ligand molecules themselves did not cause changes in the molecular ellipticities, control experiments were performed before the binding assays whereby ligand was added into Buffer 2 and scanned from 190 to 250 nm. CD spectra of all ligands tested in this study did not give rise to noticeable absorbance (data not shown). All experiments were performed using at least two separate batches of induced and refolded 6 × His-ATD2B proteins (both wild-type and mutants), each batch in duplicates or triplicates. The dose-dependent binding curves of ligands to 6 × His-ATD2B were fitted by the following equation:
where Y is the shift in ellipticity, [ligand] is the concentration of ligand investigated and K D is the dissociation equilibrium constant.
Electrophoresis and Mass Spectrometry
The wild-type and all mutant recombinant 6 × His-ATD2B proteins were separated, detected by specific antibodies recognizing hexahistidine tag (Clontech; 1:10,000) and N-terminal amino acids 27–76 of the NR2B subunit (Santa Cruz Biotechnology, Inc; 1:500). Detection of the primary antibody was performed with anti-rabbit IgG horseradish peroxidase (HRP) (1 mg/ml) followed by ECLPlus Western blotting detection reagents (Amersham Pharmacia, USA). Molecular weights of proteins were confirmed by MALDI-TOF MS analysis (Voyager DE-STR, Applied Biosystems). Trypsin-digested peptides of 6 × His-ATD2B were analyzed on 4700 Proteomic Analyzer (Applied Biosystems) and Data Explorer software .
Voltage-Clamp Recordings from Xenopus Oocytes
cRNA for rat NR1-1a (hereafter NR1) and NR2B were synthesized in vitro and injected (5–10 ng) into stage V-VI Xenopus laevis oocytes, isolated as previously described [22, 26, 58]. Two-electrode voltage-clamp current recordings were made 2–7 days post injection. The recording solution contained (in mM) 90 NaCl, 3 KCl, 10 HEPES, 0.5 BaCl2, 0.01 EDTA, 0.05 glycine (23°C); pH was adjusted to 7.3 with NaOH. EDTA (10 μM) was added to chelate contaminant extracellular divalent ions, including trace amounts of Zn2+. Solution exchange was computer controlled through an 8-modular valve positioner (Digital MVP Valve, Hamilton Company, Reno, NV). Voltage and current electrodes were filled with 0.3 and 3.0 M KCl, respectively, and current responses recorded at a holding potential of -40 mV. Data acquisition and voltage control were accomplished with a two-electrode voltage-clamp amplifier (OC-725, Warner Instruments, Hamden, CT). Data were pooled among oocytes and mean composite concentration-effect curve was fitted by the following equation:
where IC50 is the concentration of antagonist that produces a 50% inhibition, [antagonist] is the concentration of antagonist investigated, minimum a residual current response fixed to 0 and n is the Hill slope.
Homology Modeling and Docking
The protein sequence of the rat NR2B ATD (GenBank accession number U11419) was aligned against that of mGluR1 ATD (1 EWK; (36)) using a previously published alignment  [Additional file 2D]. Both the open and closed forms of mGluR1 were used as a structural basis for homology models, but due to recent work the closed model was used as the basis for docking . The homology model was then built from this modified alignment using Prime (Schrödinger), and then subjected to multiple iterative rounds of energy minimization and side chain optimization to relieve any steric strain present in the model. Ifenprodil was constructed in Maestro (Schrödinger) and docked into the homology model using Glide (Schrödinger) and induced fit docking, with the box centered on residues Asp101, Asp104, Ile150, Phe176, Phe182, and Thr233. Docking was constrained by required contact with a previously identified critical residue (Asp101) . These results were then rescored using the MM-GBSA protocol, which incorporates solvation and protein relaxation.
The shifts in ellipticities at 218.9 nm (RO25,6981), 219.6 nm (haloperidol) and 220 nm (ifenprodil), where the signal change on ligand binding was maximal, were calculated. The experiment data were analyzed by nonlinear curve fitting (ORIGIN7, USA). All pooled data were expressed as mean ± SEM. The numbers of CD assay for wild-type and mutants as well as the number of oocytes recorded on heterologously expressed NR1/NR2B receptors are shown in parentheses.
dynamic light scattering
Cull-Candy S, Brickley S, Farrant M: NMDA receptor subunits: diversity, development and disease. Curr Opin Neurobiol. 2001, 11: 327-335. 10.1016/S0959-4388(00)00215-4.
Hardingham GE, Bading H: The Yin and Yang of NMDA receptor signalling. Trends Neurosci. 2003, 26: 81-89. 10.1016/S0166-2236(02)00040-1.
Massey PV, Johnson BE, Moult PR, Auberson YP, Brown MW, Molnar E, Collingridge GL, Bashir ZI: Differential roles of NR2A and NR2B-containing NMDA receptors in cortical long-term potentiation and long-term depression. J Neurosci. 2004, 24: 7821-7828. 10.1523/JNEUROSCI.1697-04.2004.
Kew JN, Kemp JA: Ionotropic and metabotropic glutamate receptor structure and pharmacology. Psychopharmacol. 2005, 179: 4-29. 10.1007/s00213-005-2200-z.
Paoletti P, Neyton J: NMDA receptor subunits: function and pharmacology. Curr Opin Pharmacol. 2007, 7: 9-47. 10.1016/j.coph.2006.08.011.
Lynch DR, Guttmann RP: Excitotoxicity: Perspectives based on N-methyl-D-aspartate receptor subtypes. Pharmacol Exp Ther. 2002, 300: 717-723. 10.1124/jpet.300.3.717.
Loftis JM, Janowsky A: The N-methyl-D-aspartate receptor subunit NR2B: localization, functional properties, regulation, and clinical implications. Pharmacol Ther. 2003, 97: 55-85. 10.1016/S0163-7258(02)00302-9.
Arundine M, Tymianski M: Molecular mechanisms of glutamate-dependent neurodegeneration in ischemia and traumatic brain injury. Cell Mol Life Sci. 2004, 61: 657-68. 10.1007/s00018-003-3319-x.
McBain CJ, Mayer ML: N-Methyl-D-aspartate acid receptor structure and function. Physiol Rev. 1994, 74: 723-760.
Dingledine R, Borges K, Bowie D, Traynelis SF: The glutamate receptor ion channels. Pharmacol Rev. 1999, 51: 7-61.
Kemp JA, McKernan RM: NMDA receptor pathways as drug targets. Nat Neurosci. 2002, 5: 1039-1042. 10.1038/nn936.
Erreger K, Chen PE, Wyllie DJ, Traynelis SF: Glutamate receptor gating. Crit Rev Neurobiol. 2004, 16: 187-224. 10.1615/CritRevNeurobiol.v16.i3.10.
Furukawa H, Gouaux E: Mechanisms of activation, inhibition and specificity: crystal structures of the NMDA receptor NR1 ligand-binding core. EMBO J. 2003, 22: 2873-2885. 10.1093/emboj/cdg303.
Furukawa H, Singh SK, Mancusso R, Gouaux E: Subunit arrangement and function in NMDA receptors. Nature. 2005, 438: 185-192. 10.1038/nature04089.
McFeeters RL, Oswald RE: Emerging structural explanations of ionotropic glutamate receptor function. FASEB J. 2004, 18: 428-438. 10.1096/fj.03-0873rev.
Wollmuth LP, Sobolevsky AI: Structure and gating of the glutamate receptor ion channel. Trends Neurosci. 2004, 27: 321-328. 10.1016/j.tins.2004.04.005.
Kubo M, Ito E: Structural dynamics of an ionotropic glutamate receptor. Proteins. 2004, 56: 411-419. 10.1002/prot.20154.
Mayer ML: Glutamate receptor ion channel. Curr Opin Neurobiol. 2005, 15: 282-288. 10.1016/j.conb.2005.05.004.
Chen PE, Wyllie DJ: Pharmacological insights obtained from structure-function studies of ionotropic glutamate receptors. Br J Pharmacol. 2006, 147: 839-853. 10.1038/sj.bjp.0706689.
Choi YB, Lipton SA: Identification and mechanism of action of two histidine residues underlying high-affinity Zn2+ inhibition of the NMDA receptor. Neuron. 1999, 23: 171-180. 10.1016/S0896-6273(00)80763-1.
Fayyazuddin A, Villarroel A, Le Goff A, Lerma J, Neyton J: Four residues of the extracellular N-terminal domain of the NR2A subunit control high-affinity Zn2+ binding to NMDA receptors. Neuron. 2000, 25: 683-694. 10.1016/S0896-6273(00)81070-3.
Low C-M, Zheng F, Lyuboslavsky P, Traynelis SF: Molecular determinants of coordinated proton and zinc inhibition of N-methyl-D-aspartate NR1/NR2A receptors. Proc Natl Acad Sci USA. 2000, 97: 11062-11067. 10.1073/pnas.180307497.
Perin-Dureau F, Rachline J, Neyton J, Paoletti P: Mapping the binding site of the neuroprotectant ifenrpodil on NMDA receptors. J Neurosci. 2002, 22: 5955-5965.
Malherbe P, Mutel V, Broger C, Perin-Dureau F, Kemp JA, Neyton J, Paoletti P, Kew JNC: Identificaiton of critical residues in the amino terminal domain of the human NR2B subunit involved in the RO25-6981 binding pocket. J Pharmacol Exp Ther. 2003, 307: 897-905. 10.1124/jpet.103.056291.
Rosahl TW, Wingrove PB, Hunt V, Fradley RL, Lawrence JM, Heavens RP, Treacey P, Usala M, Macaulay A, Bonnert TP, Whiting PJ, Wafford KA: A genetically modified mouse model probing the selective action of ifenprodil at the N-methyl-D-aspartate type 2B receptor. Mol Cell Neurosci. 2006, 33: 47-56. 10.1016/j.mcn.2006.06.006.
Traynelis SF, Burgess MF, Zheng F, Lyuboslavsky P, Powers JL: Control of voltage-independent zinc inhibition of NMDA receptors by the NR1 subunit. J Neurosci. 1998, 18: 6163-6175.
Masuko T, Kashiwagi K, Kuno T, Nguyen ND, Pahk AJ, Fukuchi J-I, Igarashi K, Williams K: A regulatory domain (R1–R2) in the amino terminus of the N-methyl-D-aspartate receptor: effects of spermine, protons, and ifenprodil, and structural similarity to baverial leucine/isoleucine/valine binding protein. Mol Pharmacol. 1999, 55: 957-969.
Rachline J, Perin-Dureau F, Le Goff A, Neyton J, Paoletti P: The micromolar zinc-binding domain on the NMDA receptor subunit NR2B. J Neurosci. 2005, 25: 308-317. 10.1523/JNEUROSCI.3967-04.2005.
Paoletti P, Perin-Dureau F, Fayyazuddin A, Le Goff A, Callebaut I, Neyton J: Molecular organization of a zinc binding N-terminal modulatory domain in a NMDA receptor subunit. Neuron. 2000, 28: 911-925. 10.1016/S0896-6273(00)00163-X.
Williams K: Ifenprodil discriminates subtypes of N-methyl-D-aspartate receptor: selectivity and mechanisms at recombinant heteromeric receptors. Mol Pharmacol. 1993, 44: 851-859.
Williams K: Ifenprodil, a novel NMDA receptor antagonist: site and mechanism of action. Curr Drug Targets. 2001, 2: 285-298. 10.2174/1389450013348489.
Gallagher MJ, Huang H, Pritchett DB, Lynch DR: Interactions between ifenprodil and the NR2B subunit of the N-methyl-D-aspartate receptor. J Biol Chem. 1996, 271: 9603-9611. 10.1074/jbc.271.16.9603.
Brimecombe JC, Gallagher MJ, Lynch DR, Aizenman E: An NR2B point mutation affective haloperidol and CP101,606 sensitivity of single recombinant N-methyl-D-aspartate receptors. J Pharmacol Exp Ther. 1998, 286: 627-634.
Lynch DR, Shim SS, Seifert KM, Kurapathi S, Mutel V, Gallagher MJ, Guttmann RP: Pharmacological characterization of interactions of RO25-6981 with the NR2B (epsilon 2) subunit. Eur J Pharmacol. 2001, 416: 185-195. 10.1016/S0014-2999(01)00868-8.
O'Hara PJ, Sheppard PO, Thøgersen H, Venezia D, Haldeman BA, McGrane V, Houamed KM, Thomsen C, Gilbert TL, Mulvihill ER: The ligand-binding domain in metabotropic glutamate receptors is related to bacterial periplasmic binding protein. Neuron. 1993, 11: 41-52. 10.1016/0896-6273(93)90269-W.
Kunishima N, Shimade Y, Tsuji Y, Sato T, Yamamoto M, Kumasaka T, Nakanishi S, Jingami H, Morikawa K: Structural basis of glutamate recognition by a dimeric metabotropic glutamate receptor. Nature. 2000, 407: 971-977. 10.1038/35039564.
Mayer ML: Glutamate receptors at atomic resolution. Nature. 2006, 440: 456-462. 10.1038/nature04709.
Safferling M, Tichelaar W, Kümmerle G, Jourppila A, Kuusinen A, Keinänen K, Madden DR: First image of a glutamate receptor ion channel: oligomeric state and molecular dimensions of GluRB homomers. Biochemistry. 2001, 40: 13948-13953. 10.1021/bi011143g.
Tichelaar W, Safferling M, Keinänen K, Stark H, Madden DR: The three-dimensional structure of an ionotropic glutamate receptor reveals a dimer-of-dimers assembly. J Mol Biol. 2004, 344: 435-442. 10.1016/j.jmb.2004.09.048.
Nakagawa T, Cheng Y, Ramm E, Sheng M, Walz T: Structure and different conformational states of native AMPA receptor complexes. Nature. 2005, 433: 545-549. 10.1038/nature03328.
Wong E, Ng F-M, Yu C-Y, Lim P, Lim L-H, Traynelis SF, Low C-M: Expression and characterization of soluble NR2B subunit N-terminus of NMDA receptor. Protein Sci. 2005, 14: 2275-2283. 10.1110/ps.051509905.
Chen GQ, Gouaux E: Overexpression of a glutamate receptor (GluR2) ligand binding domain in Escherichia coli: application of a novel protein folding screen. Proc Natl Acad Sci USA. 1997, 94: 13431-13436. 10.1073/pnas.94.25.13431.
Greenfield NJ: Circular dichroism analysis for protein-protein interactions. Methods Mol Biol. 2004, 261: 55-78.
Heyn MP, Weischet WO: Circular dichroism and fluorescence studies on the binding of ligands to the α subunit of tryptophan synthase. Biochemistry. 1975, 14: 2962-2968. 10.1021/bi00684a026.
Hebenstreit D, Ferreira F: Structural changes in calcium-binding allergens: use of circular dichroism to study binding characteristics. Allergy. 2005, 60: 1208-1211. 10.1111/j.1398-9995.2005.00907.x.
Greenfield NJ: Using circular dichroism collected as a function of temperature to determine the thermodynamics of protein unfolding and binding interactions. Nature Protocols. 2006, 1: 2527-2535. 10.1038/nprot.2006.204.
Lynch DR, Gallagher MJ: Inhibition of N-methyl-D-aspartate receptors by haloperidol: developmental and pharmacological characterization in native and recombinant receptors. J Pharmacol Exp Ther. 1996, 279: 154-161.
Ilyin VI, Whittemore ER, Guastella J, Weber E, Woodward RM: Subtype-selective inhibition of N-methyl-D-aspartate receptors by haloperidol. Mol Pharmacol. 1996, 50: 1541-1550.
Whittemore ER, Ilyin VI, Woodward RM: Antagonism of N-methyl-D-aspartate receptors by σ site ligands: potency, subtype-selectivity and mechanisms of inhibition. J Pharmacol Exp Ther. 1997, 282: 326-338.
Warren GL, Andrews CW, Capelli A-M, Clarke B, LaLonde J, Lambert MH, Lindvall M, Nevins N, Semus SF, Senger S, Tedesco G, Wall ID, Woolven JM, Peishoff CE, Head MS: Comparison of shape-matching and docking as virtual screening tools. J Med Chem. 2006, 49: 5912-5931. 10.1021/jm050362n.
Schoemaker H, Allen J, Langer SZ: Binding of [3H]ifenprodil, a novel NMDA antagonist, to a polyamine-sensitive site in the rat cerebral cortex. Eur J Pharmacol. 1990, 176: 249-250. 10.1016/0014-2999(90)90539-I.
Kiss L, Cheng G, Bednar B, Bednar RA, Bennett PB, Kane SA, McIntyre CJ, McCauley JA, Koblan KS: In vitro characterization of novel NR2B selective NMDA receptor antagonists. Neurochem Int. 2005, 46: 453-464. 10.1016/j.neuint.2004.12.008.
McCauley JA, Theberge CR, Romano JJ, Billings SB, Anderson KD, Claremon DA, Freidinger RM, Bednar RA, Mosser SD, Gaul SL, Connolly TM, Condra CL, Xia M, Cunningham ME, Bednar B, Stump GL, Lynch JJ, MaCaulay A, Wafford KA, Koblan KS, Liverton NJ: NR2B-selective N-methyl-D-aspartate antagonists: synthesis and evaluation of 5-substituted benzimidazoles. J Med Chem. 2004, 47: 2089-2096. 10.1021/jm030483s.
Marinelli L, Cosconati S, Steinbrecher T, Limongelli V, Bertamino A, Novellino E, Case DA: Homology modelling of NR2B modulatory domain of NMDA receptor and analysis of ifenprodil binding. ChemMedChem. 2007, 2: 1498-1510. 10.1002/cmdc.200700091.
Gerber AM, Vallano ML: Structural properties of the NMDA receptor and the design of neuroprotective therapies. Mini Rev Med Chem. 2006, 6: 805-815. 10.2174/138955706777698651.
Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR: Site-directed mutagenesis by overlap extension using the polymerase chain reaction. Gene. 1989, 77: 51-59. 10.1016/0378-1119(89)90358-2.
Horton RM, Hunt HD, Ho SN, Pullen JK, Pease LR: Engineering hybrid genes without the use of restriction enzymes: gene splicing by overlap extension. Gene. 1989, 77: 61-68. 10.1016/0378-1119(89)90359-4.
Low C-M, Lyuboslavsky P, French A, Le P, Wyatte K, Thiel WH, Marchan EM, Igarashi K, Kashiwagi K, Gernert K, Williams K, Traynelis SF, Zheng F: Molecular determinants of proton-sensitive N-methyl-D-aspartate receptor gating. Mol Pharmacol. 2003, 63: 1212-1222. 10.1124/mol.63.6.1212.
We thank Professor Stephen F. Heinemann (Salk Institute) for sharing NR1 and NR2B cDNAs (GenBank accession no. U11418 and U11419). We thank Professor Xiaodong Cheng and Professor Xing Zhang (Emory University, GA) for valuable advice throughout this study, including hosting C-M.L in their laboratory for 2 months. We also thank Professor Eric Gouaux (Vollum Institute) for the 6 × His-tag protein expression vector. We thank Peiqi Lim for her excellent technical assistance and Yoke-Ping Cheong for help in figure preparation. F-M.N is the recipient of a graduate research scholarship from the National University of Singapore and partially supported by the BMRC grant. The authors thank the NUH-NUS Medical Publications Support Unit, Singapore, for assistance in the preparation of this manuscript. This work was supported by NIH-NINDS, NARSAD, Michael J Fox Foundation to SFT, Biomedical Research Council grant (R184000119305), Academic Research Funds (R184000114112) and National Medical Research Council grant (0581/2001) to C-M.L.
The authors declare that they have no competing interests.
FMN carried out the recombinant protein biochemical studies. MTG carried out and JPS participated in the homology modeling and docking studies. SFT participated in the functional oocyte studies. CML conceived of the study, participated in its design and coordination, and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Refolding buffer compositions and dialysis temperatures of 6 × His-ATD2B recombinant protein.(TIFF 449 KB)
Additional file 2: Construction and expression of ATD of NR2B subunit. A. Sschematic representation of the ATD2B construct (Arg27 to Arg393) with 6xHis-tag and linker. Two stop codons are inserted (**). The four mutated residues that were investigated in this study are labeled. B. Coomassie blue-stained 11% SDS-PAGE gel showing insoluble protein fraction induced with 0.5mM IPTG (lane 1), soluble protein fraction (lane2) and purified protein (47kDa) after elution with 300mM imidazole (lane 3). C. 6xHis-ATD2B probed with anti-6xHis (left) and anti-NR2B (right) primary antibodies. D. Rat NR2B and mGluR1 ATD sequence alignment (CLUSTAL W 1.83) and secondary structure assignment based on mGluR1 . Helices are shown as red cylinders, strands as blue arrows, with the same labeling as in Fig. 2 of . Intraprotomer disulphide bridges identified in mGluR1 are indicated by thin solid black lines connecting cysteines. Red amino acid residues indicate selected critical molecular determinants to ifenprodil binding within ATD of NR2B which are mutated in this study (also see  for more key residues); blue amino acid residue indicates non-critical residue to ifenprodil binding. Residues that are identical between mGluR1 and NR2B are shaded black while conserved residues are shaded grey. (TIFF 1 MB)
Additional file 3: Oligonucleotide primers used in mutagenesis experiments. TA represents the annealing temperature. The nucleotides coding for the mutant amino acid are underlined. (TIFF 244 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ng, FM., Geballe, M.T., Snyder, J.P. et al. Structural insights into phenylethanolamines high-affinity binding site in NR2B from binding and molecular modeling studies. Mol Brain 1, 16 (2008). https://doi.org/10.1186/1756-6606-1-16