- Research article
- Open Access
Antagonistic action on NMDA/GluN2B mediated currents of two peptides that were conantokin-G structure-based designed
BMC Neuroscience volume 18, Article number: 44 (2017)
The GluN2B subunit of the N-methyl-d-aspartate receptor (NMDAr) modulates many physiological processes including learning, memory, and pain. Excessive increase in NMDAr/GluN2B activity has been associated with various disorders such neuropathic pain and neuronal death following hypoxia. Thus there is an interest in identifying NMDAr antagonists that interact specifically with the GluN2B subunit. Recently based on structural analysis between the GluN2B subunit and conantokin-G, a toxin that interacts selectively with the GluN2B subunit, we designed various peptides that are predicted to act as NMDAr antagonists by interacting with the GluN2B subunit. In this study we tested this prediction for two of these peptides EAR16 and EAR18.
The effects of EAR16 and EAR18 in NMDA-evoked currents were measured in cultured rat embryonic hippocampal neurons and in HEK-293 cells expressing recombinant NMDAr comprised of GluN1a–GluN2A or GluN1a–GluN2B subunits. In hippocampal neurons, EAR16 and EAR18 reduced the NMDA-evoked calcium currents in a dose-dependent and reversible manner with comparable IC50 (half maximal inhibitory concentration) values of 241 and 176 µM, respectively. At 500 µM, EAR16 blocked more strongly the NMDA-evoked currents mediated by the GluN1a–GluN2B (84%) than those mediated by the GluN1a–GluN2A (50%) subunits. At 500 µM, EAR18 blocked to a similar extent the NMDA-evoked currents mediated by the GluN1a–GluN2B (62%) and the GluN1a–GluN2A (55%) subunits.
The newly designed EAR16 and EAR18 peptides were shown to block in reversible manner NMDA-evoked currents, and EAR16 showed a stronger selectivity for GluN2B than for GluN2A.
N-methyl-d-aspartate receptors (NMDAr) are inotropic glutamate receptors (iGluRs) that contribute to multiple neuronal functions such as [1,2,3]. However, excessive NMDAr activity can lead to neuronal dysfunction and neuronal death [4, 5]. Functional NMDAr require two GluN1 and two GluN2 subunits [4, 6,7,8]. The GluN2 subunit can be replaced by a GluN3 subunit, but when present the GluN3 subunit decreases the NMDAR activity [9, 10]. There are four GluN2 subunits (GluN2A, B, C and D) . An increase in the expression of GluN1, GluN2A, and GluN2B subunits in the hippocampus contribute to neurological problems that develop following traumatic brain injury . Increase in the activity/expression of GluN2B has been postulated to contribute to neuronal damage following stroke, as well as to the development of diseases such as Parkinson, Huntington, Alzheimer and chronic neuropathic pain .
NMDAr are similar to other iGluRs in that each subunit has four domains, which include the amino-terminal domain (ATD), an agonist/ligand binding domain (ABD), a trans-membrane domain (TMD) and a carboxyl terminal domain (CTD) [13,14,15]. In GluN2 subunits, the ATD domain contains binding sites for allosteric inhibitors such as zinc  and ifenprodil . The ABD domain in GluN1 and GluN3 binds the coagonist glycine, whereas the ABD domain in GluN2 binds the agonist glutamate [11, 13,14,15]. Some peptides isolated from marine cone snails have shown specificity for GluN2 subunits of the NMDAr . One member of the conantokin family, conantokin-G is a 17-residue peptide that shows selectivity for NMDAr containing the GluN2B subunit . Although, conantokin-G shows high selectivity for GluN2B, its reversibility appears to be very slow and incomplete on neuronal cells . We designed a number of peptides with potential GluN2B–NMDAr antagonistic activity, this was done by using point mutations on the conantokin-G sequence and evaluating their structure as well as their binding capacity to the ligand binding domain of the GluN2B subunit by using various in silico evaluations. Here we tested this prediction for two of these peptides, EAR16 and EAR18, and found that both of these peptides blocked in a highly reversible manner NMDA-evoked currents in hippocampal neurons, and that EAR16 was more selective for NMDArs containing the GluN2B than for those containing the GluN2A subunit. The high reversibility of both of these peptides and the higher selectivity of EAR16 for the NMDArs containing the GluN2B increase their potential as pharmacological agents.
The peptides EAR16 (GEDDLQDNQDLIRDKSN) and EAR18 (GEDDYQDAQDLIRDKSN) were designed in the Bioinfo Grip laboratory at the National University of Colombia (patent in progress, the peptides are commercially available) to interact with the ligand binding domain (LBD) of the GluN2B subunit. The peptide design and their in silico evaluation was achieved by using point mutations on the conantokin-G sequence and evaluating their binding capacity with AutoDock Vina  and Rosetta FlexPepDock server , by using a structural model of GluN2B that was developed in our laboratory (Reyes-Guzman et al. in press), as briefly described below. As a control, we initially selected a scramble peptide with the same length (17 a.a.) like is usual in this kind of design. However, this scramble peptide displayed agonist-like activity, potentially by interacting at a site different than the LBD. Hence we selected another non-related peptide (EAR7) as control (TKRSSRAFRE), one that was designed to interact with an intracellular site, the C-terminal domain (CTD) of GluN2B subunit . EAR16 and EAR 18 show good water solubility (http://pepcalc.com) and have similar hydrophobicity levels: 19.36 and 19.65, respectively (https://www.thermofisher.com/co/en/home/life-science/protein-biology/peptides-proteins/custom-peptide-synthesis-services/peptide-analyzing-tool.html). Permeability was predicted using CellPPD (http://crdd.osdd.net/raghava/cellppd/index.html) which provided negative CPP values for these peptides (−0.06, −0.57 and −0.57 for EAR7, EAR16 and EAR18, respectively), indicating that they are non-cell permeant (CPP less than 0, indicate that the peptide is non-cell permeant). The peptide synthesis was made by a commercial supplier (Peptide 2.0 in Sterling, VA. USA). These peptides show no modifications at either terminal, their N terminal is a free amine and the C terminal is a free acid (analysis performed by Peptide 2.0. The peptide purity and other properties are shown in Table 1.
Peptide three-dimensional structure
For each peptide a three-dimensional structure model was generated using UCSF (University of California San Francisco) Chimera software (http://www.cgl.ucsf.edu/chimera/, ), and a refined three-dimensional structure was obtained in an hydrophilic environment using PEPstr (Peptide Tertiary Structure Prediction Server; http://www.imtech.res.in/raghava/pepstr/) as previously described .
Structural model of the GluN2B
Details of this model are described in a separate publication . Briefly, this model was obtained by using the amino acids 1-817 of the reported sequence of the GluN2B subunit (Ensembl ENSRNOP00000011697, Rattus norvegicus). This amino acid sequence includes the amino terminal (ATD) and the ligand binding (LBD) domains. The GluN2B structural model was predicted using I-TASSER (http://zhanglab.ccmb.med.umich.edu/, ), this software used several templates to determine the best model. The best model (based on the C-Score and TM-score) was obtained when the template used was the inotropic glutamate AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor (PDB 3KG2) . This template has been previously used to elucidate other functions of glutamate receptors . The C-score usually ranges between −5 and 2, the most positive the value the better the model topology, the GluN2B model selected has a C-core of 0.01 . The TM-score provides a value of the structural correspondence between two structures (GluN2B-template), values higher than 0.5 indicate that the model topology is likely to be correct , the selected GluN2B model has a TM-score of 0.71.
The selected GluN2B model (using I-TASSER) was further evaluated using PDBsum (http://www.ebi.ac.uk/pdbsum/, ), that provides various parameters including hydrophobicity profile, solvent accessibility, prediction of ligand binding sites. The secondary structure of the model was also evaluated using Ramachandran plots (Procheck in PDBsum), which provides the percentage of structural regions that are either “favorable”, “allowed” or “not-allowed”. A model with good stereo-chemistry should have close to 90% of the residues in structurally favorable regions . The selected model has 86.7% of the residues located in structurally “favorable” regions, and 10.9% of the residues in “allowed” regions . The model has structural topology comparable to that described for other inotropic glutamate receptors , including the ATD, LDB and the transmembrane (TDM) domains. The ATD and LBD structural configurations are composed of mostly α helices and allows the binding of modulators an agonists .
Molecular docking was performed with the obtained three-dimensional structures of the peptides and of the GluN2B subunit by using “Autodock Vina” (http://vina.scripps.edu/, ). A Grid Box was generated at the center of the LBD of the GluN2B model to which the corresponding hydrogen atoms and chargers were assigned. The Grid Box was 1 Å, in the software “completeness” was set to 10 and the number of results per Docking was set to 10, for each peptide. The peptides corresponded to the ligands. The best conformations and orientations, which displayed the lowest values for binding energy (−6.6 and −7.7 kcal/mol for EAR16 and EAR18, respectively), were visualized using Autodock Tools (v. 1.5.4.).
The Docking results were then analyzed with the Rosetta FlexPepDock web server (http://flexpepdock.furmanlab.cs.huji.ac.il/index.php, ), which refines the peptide-receptor interaction by optimizing the peptide orientation within the binding site using the Monte-Carlo method and focusing on energy minimization of the ligand with the receptor. With the starting structure FlexPepDock realized 200 independent simulations, 100 simulations were done under high resolution. The other 100 simulation included a step of pre-optimization at low resolution, followed by a method of refinement of high resolution. This resulted in 200 models for each peptide, we selected the best one for each peptide based on the score provided by the punctuation function of Rosetta (675.461 and 755.051 for EAR16 and EAR18, respectively) . The docking images were generated using UCSF Chimera (http://www.cgl.ucsf.edu/chimera/, ).Interactions between peptide-receptor were determine using Discovery Studio Visualizer v4.1 software (http://accelrys.com/products/discovery-studio/visualization.html).
Pregnant rats were maintained and used following the guidelines approved by the New York University Langone Medical Center Institutional Animal Care and Use Committee. Pregnant rats were killed using CO2 followed by cervical dislocation. Day 18 embryos (E18) were removed by caesarian section. Dissociated hippocampal cultures were prepared as previously described [29, 30]. Cells were plated on 12 mm coverslips (7.0 × 104 cells/well, 24 well plates) with neurobasal medium supplemented with B27 (Gibco Invitrogen, ThermoFisher Scientific, USA) or with 2% bovine fetal serum and used between 7 and 37 days in vitro; the data was collected from 57 cells that had on average 14.9 days in vitro (14.9 ± 10.8 mean ± SD). Coverslips were pre-treated by placing 50 µL of poly-d-Lysine (37.5 μg/mL, Sigma-Aldrich, USA) and Laminin (2.5 μg/mL, Invitrogen) in the middle of the coverslip for 5 min at room temperature (RT).
Transfection of HEK-293 cells
The plasmids pEGFP-NR1a (GluN1a) (Addgene plasmid # 17926), pEGFP-NR2B (GluN2B) (Addgene plasmid # 17925), and pEGFP-NR2A (GluN2A) (Addgene plasmid # 17924) were a gift from Stefano Vicini . HEK-293 cells (ATCC, Manassas,Virginia, USA) were grown on 12-mm coverslips in 24 well plates. HEK-293 cells (60–70% confluent) were transfected using Lipofectamine 2000 according to the product instructions (Invitrogen), with a mixture of GluN1a and GluN2A subunits or a mixture of GluN1a and GluN2B subunits. The amount of cDNA for each plasmid was 0.25 µg/well. Cells were used 2–4 days following transfection. GFP signal was used to identify positively transfected cells.
Electrophysiological recordings were performed under voltage-clamp by using the whole cell conformation [32, 33]. Intracellular solution (in mM): 110 Cs-gluconate, 20 CsCl, 10 Hepes, 10 EGTA, 4 Mg-ATP, 0.4 Na-GTP, pH 7.3 (adjusted with CsOH) . Extracellular solution (in mM): 140 NaCl, 2.5 KCl, 2.0 CaCl2, 10 Hepes, 10 d-glucose, pH 7.4. Mg2+ was omitted to prevent the voltage-dependent block of the NMDAr ion currents [35, 36]. The following blockers were added to extracellular solution: 1.0 μM tetrodotoxin (TTX; Sigma, St. Louis, MO) to block Na+ channels ; 20 μM CNQX (6-cyano-7-nitroquinoxaline-2,3-dione, Sigma) to block AMPA and Kainate receptors, 50 μM Bicuculine (Sigma) and 100 μM Picrotoxin (Sigma) to block GABA receptors, and 1 μM Estricnin to block glycine channels. Currents were measured using borosilicate pipettes (1.5 mm outer diameter; 0.86 mm inner diameter with filaments, World Precision Instruments, Inc, Sarasota, FL) with resistances of 1–3 MΩ. Experiments were conducted at room temperature.
Currents were measured using a 2-electrode voltage clamp amplifier Axopatch 200B (Axon Instruments, Foster City, CA; Molecular Devices, LLC, Sunnyvale, CA). Output from the voltage clamp amplifier was sent to a microcomputer using a data acquisition interface (Digidata 1440A; Axon Instruments). Currents were filtered at 2.5 or 5 kHz (Axopatch 200B) and sampled at 10–50 kHz. For analysis the current traces were filtered at 1 kHz. Pclamp 10.4 (Axon Instruments) was used for data acquisition and data analyses.
Cells were held at −60 mV and the NMDAr-evoked currents were obtained by applying a 5 s pulse of NMDA (50 or 100 μM NMDA + glycine 10 μM) followed by a 10 s pulse of bath solution through a puffer (flow rate: 250 μL/min) located directly on top of the cells (perfusion manifold MMF; Scientifica, East Sussex, United Kingdom). Successive stimulations with NMDA were done every 4 min, which under control conditions resulted in NMDA-evoked Ca2+ currents of similar magnitude over time (Fig. 1a). When the effects of (+)MK-801(Dizocilpine, Tocris) a general NMDAr blocker, Ro 25-6981 (Sigma) a specific GluN2B/NMDAR blocker), and peptides (EAR16 and EAR18) were evaluated, the cells were first exposed to two or three consecutive NMDA pulses (4 min interval) to assure that the cell recording was stable, then the blocker (or peptide) was applied 10 s prior and during the following NMDA pulse. In all cases the NMDA pulses were followed by a 10 s pulse of bath solution through the puffer. The chamber was continuously perfused with the bath solution (flow rate: 400 μL/min), except when the puffer was on. Peptides were added at a final concentration of 10, 100 and 500 μM. The area of the NMDA-evoked inward currents was measured over a period of 5 s, from the onset of the NMDA-evoked inward current. The areas were normalized using the NMDA-evoked inward current measured prior to the addition of the blocker or peptide (−1: maximal inward current; 0: no current). Figures were done using GraphPad Prism v6. Results are expressed as mean and standard error of mean (mean ± S.E.M).
When comparing two groups with various stimulations we used Two-way ANOVA (analysis of variance) followed by Sidak’s multiple comparison post-test. When comparing more than two groups with various stimulations we used Two-way ANOVA followed by Dunnett’s multiple comparison post-test, as recommended by GraphPadPrism 7.02. Results are expressed as mean and standard error mean (mean ± S.E.M.).
Electrophysiological evaluation of EAR16 and EAR18 peptides on hippocampal primary cultures
Primary hippocampal cultures were used between 4 and 20 days of culture, a period when the GluN2B subunit has been shown to be highly expressed [37, 38]. When NMDA stimulations were separated by a 4 min interval the magnitude of the NMDA-evoked Ca2+ currents was similar between consecutive stimulations (Fig. 1a), as reflected by the average behavior shown by the black squares in Figs. 1b, c and 2a, b. The general NMDA blockers (+)MK-801 (10 µM), blocked the NMDA-evoked currents (Fig. 1b). The magnitude of (+)MK-801-mediated block increased with exposures to the blocker and it was irreversible (Fig. 1b), as previously reported in hippocampal neurons . In order to determine the contribution of the GluN2B to the NMDA-evoked response we exposed the hippocampal cells to a GluN2B selective blocker Ro 25-6981 (1 µM), at a concentration that is 100 fold higher that producing 50% inhibition of NMDA/GluN2B-mediated currents . At this concentration Ro 25-6981 has minimal or no effects on voltage-dependent Na+ and Ca+2 channels . Block with Ro 25-6981 showed no reversibility following 4 min of washout (Fig. 1c, top). We observed that Ro 25-6981 blocked 70% of the NMDA-evoked currents (Fig. 1c, bottom); hence under our experimental conditions about 70% of the NMDA-evoked currents in hippocampal neurons are mediated by activation of GluN2B/NMDAr. Moreover, the increase in the magnitude of the block observed following consecutive exposures to Ro 25-6981 is consistent with the blockers being use-dependent.
In hippocampal neurons, the effect of EAR16 and EAR18 on the NMDA-evoked currents was accessed by applying each peptide for four consecutive NMDA-evoked currents to determine whether these peptides blocked the NMDA-evoked current and if so whether their block displayed use-dependence behavior. EAR18 inhibited NMDA-evoked currents in a dose-dependent manner, and this inhibition was reversible (Fig. 2a, c). EAR16 had a similar inhibitory effect on NMDA-evoked currents in hippocampal neurons (Fig. 2b, d). For a given dose, the magnitude of the block mediated by EAR18 or EAR16 did not increase significantly with consecutive exposures (Fig. 2a, b); hence these peptides show no clear use-dependence (compared to that observed by Ro 25-6981 Fig. 1c). Examination of the current traces, show that following the simultaneous removal of the NMDA + EAR16 (or NMDA + EAR18) there was a relatively short-lived increase in the NMDA-evoked current, suggesting that the off rate of EAR16 (and EAR18) is fast enough that would allow binding of NMDA early on during washout (Fig. 2c, d indicated with arrows). Figure 2e shows the dose–response curves. When considering the total current magnitude the estimated IC50 values were similar for both peptides being 176 and 242 µM for EAR18 and EAR16, respectively (smooth lines). No inhibitory action was observed with the control peptide at 500 µM (data not shown).
EAR16 displays higher selectivity than EAR18 for the GluN2B subunit over the GluN2A subunit
To test the selectively of these peptides for the GluN2B subunit, the blocking effect of these peptides was tested in HEK-293 cells expressing recombinant NMDAr comprised of either GluN1a–GluN2A or GluN1a–GluN2B subunits. The rationale for doing so is that these are the predominant GluN2 subunits in hippocampal neurons, moreover the amino acid sequence corresponding to the sequence used for the structural model of GluN2B (1-818) is conserved (71.9% identity) and of the LBD is highly conserved (95.5% identity) between GluN2A and GluN2B [25, 41, 42]. We selected 500 µM of each peptide to test their potency on HEK293 cells expressing recombinant NMDAr, since this peptide concentration produced a strong but not maximal block (~60%) on NMDA-evoked currents in hippocampal neurons (Fig. 2e). This should allow to determine whether or not these peptides displayed differences in potency between the GluN2B and GluN2A subunits.
Figure 3 shows that as predicted EAR16 at 500 µM blocked more strongly NMDA-evoked currents in HEK cells transfected with the GluN1a–GluN2B subunits (Fig. 3a, e) than those in cells transfected with the GluN1a–GluN2A subunits (Fig. 3b, e). In contrast, EAR18 at 500 µM blocked to a similar level the NMDA-evoked currents in HEK-cells transfected with either GluN1a–GluN2B (Fig. 3c, f) or GluN1a–GluN2A subunits (Fig. 3d, f). The IC50 of each peptide for blocking GluN1a–GluN2B and GluN1a–GluN2A was estimated by using Eq. (1), where ‘x’ is the peptide concentration in log10 (500 µM in log10 = 2.69897); ‘y’ is the % inhibition at ‘x’ (from Fig. 3c, f); Emax: is the maximal inhibitory effect which we are assuming to be 100%; ‘n’ is the slope of the fitting, which it is assumed to be “1” since there is only a single recombinant NMDAr (containing either the GluN2B or the GluN2A).
The results are shown in Fig. 3g, and indicate that EAR16 has about 100 fold higher affinity for GluN1a–GluN2B than for GluN1a–GluN2A (3.6 vs 416.1 µM). While EAR18 displayed a comparable affinity for both GluN1a–GluN2B and GluN1a–GluN2A (30.2 and 18.4 µM).
As observed in hippocampal neurons, in HEK-293 cells transfected with recombinant NMDAr the simultaneous removal of the NMDA + EAR16 (or NMDA + EAR18) also resulted in a short-lived increase in the NMDA-evoked current, suggesting that the off rate of EAR16 (and EAR18) was fast enough that would allow binding of NMDA early on during washout.
Docking of EAR16 and EAR18 with the LBD of the GluN2B model
Figure 4 shows the docking results of the interaction between the EAR16 and the GluN2B subunit. EAR16 binds to the GluN2B through 7 hydrogen bridges, 4 of them are with residues located at the LBD in GluN2B (His 486, Ser 690, Asp 732, Tyr 762) . The EAR16 residues Asp 4 and 7 contribute to the formation of these hydrogen bridges. The residue Met 739, that has been reported to contribute to the molecular interaction between Con-G and GluN2B , is predicted to form a hydrophobic interaction (alkyl) with the Leu 11 of EAR16. The majority of the interactions involve the amino terminal of EAR16, similar to what has been described for conantokin-G , indicating that the substitution of the residue γ-carboxiglutamic acid (Gla) by Asp appears to emulate the Con-G interaction with the GluN2B.
In the case of EAR18, the interactions involve 9 hydrogen bridges (Fig. 5), 3 of them are with residues located at the LBD in GluN2B (Thr 514, Ser 690, Tyr 762). The EAR18 Asp 4 residue plays an important role in establishing these 3 hydrogen bridges (Fig. 5). The Tyr 5 of EAR18 contributes to the formation of a hydrophobic interaction (Pi–sigma) with the Ile 534 of the GluN2B, and the Ala 8 of EAR18 is forming and hydrogen bridge with the Lys 485 of the GluN2B. At the LBD, the EAR18 conformation is less extended than that of the EAR16, allowing EAR18 to interact with a larger number of residues in the GluN2B.
In this study we demonstrated that the peptides EAR16 and EAR18 can inhibit NMDA-evoked currents in hippocampal neurons in a dose-dependent and highly reversible manner. We also demonstrated that these peptides inhibited NMDA-evoked currents evoked by recombinant NMDAr containing the GluN1a–GluN2B subunits. However, only EAR16 showed higher selectivity for GluN1a–GluN2B over GluN1a–GluN2A.
Conantokin-G is a 17 amino acid peptide (MW 2264.2 Da) blocks NMDA-evoked currents in hippocampal  and cortical  neurons. In expression systems, conantokin-G has been shown to be selective for NMDAr containing the GluN2B subunit . EAR16 and EAR18 were designed based on evaluations of in silico interactions by using point mutation on the conantokin-G sequence and by evaluating their binding capacity to an extracellular site of the GluN2B receptor. EAR16 and EAR18 also have 17 amino acids, between them their sequence differs in two amino acids. Both peptides have the same calculated isoelectric point of 3.76, which is higher to that reported for conantokin-G of 2.02 . EAR16 and EAR18 peptides do not have γ-carboxiglutamic acid (Gla) which is present in conantokin-G. The latter was chosen to increase the flexibility of the peptide’s structure, with the rationale that this would facilitate their interactions with the conantokin-G binding pocket in the GluN2B subunit.
The prediction for the three-dimensional structure for EAR18 and EAR16 was performed using a hydrophilic environment. Figure 6 shows the predicted structure for EAR18 and EAR16 and the one reported for conantokin-G [47, 48]. While conantokin-G adopts a helical conformation in more than 70% of its sequence [47, 48], EAR18 and EAR16 adopt a helical conformation in 50% (EAR18) or less (EAR16) of their sequence. The high level of helical conformation in conantokin-G is due to the presence of Gla, which coordinates binding of 4 calcium ions [48,49,50]. The decrease in the proportion of helical conformation will increase the molecular flexibility of EAR18 and EAR16 and this may contribute to the observed high reversibility of EAR16 and EAR18 (present study). Conantokin-G displays a very slow and incomplete reversibility in cortical neurons, such that after 3 min of washout only about 10% of the NMDA-evoked current was recuperated . In the same study, it was found that the reversibility of conantokin-G was faster and more complete in oocytes expressing GluN1a–GluN2B, such that after 4 min of washout the NMDA-evoked currents displayed about 80% recovery . We observed that EAR16 and EAR18 both were highly reversible in both hippocampal and in HEK cells expressing recombinant NMDAr, but again it also appears that their reversibility was faster from the recombinant NMDAr than from NMDAr expressed in neuronal cells. Together, the results indicate that the conformation of the ligand binding domain (LBD) in the GluN2 subunits may be affected by their neuronal environment. Regardless of such difference we found that EAR16 and EAR18 were highly reversible. This high reversibility represents an advantage for the use of these peptides as potential pharmacological agents in comparison to other available NMDAr blockers, including general NMDAr blockers ((+)MK801) and those selective for GluN2B (conantokin-G, Ro 25-6981).
The Met 739 residue in GluN2B, is not present in the GluN2A subunit and has been reported to be a molecular determinant for the selective antagonist of conantokin-G on GluN2B over the GlunN2A subunit [44, 51]. The Docking results show the formation of a hydrophobic (alkyl type) interaction between the Leu 11 of EAR16 and the Met 739-GluN2B, such interaction is not observed for any of the residues of EAR18. This interaction may contribute to the apparent higher selectivity of EAR16 over EAR18 for the GluN2B subunit over the GluN2A subunit.
Previously by testing various conantokin analogs, it was found that the presence and location of Gla appeared to be important to confer to the conantokin analogs selectivity for the GluN2B subunit over the GluN2A subunit . We found that EAR18, in which Gla was replaced by Asp, did not display selectivity between GluN2B and GluN2A, which is consistent with the idea that Gla is important to confer selectivity for GluN2B. However, we also found that EAR16, in which Gla was also replaced by Asp, showed higher selectivity for GluN2B over GluN2A, indicating that Gla is not the only residue conferring selectivity for GluN2B over GluN2A.
There are only two amino acids difference between EAR18 and EAR16, in EAR16 these are a Leucine (Leu 5) and an asparagine (Asn 8), which in EAR18 correspond to tyrosine (Tyr 5) and alanine (Ala 8). Our results are consistent with previous observation in which it has been reported that replacement of Leu 5 with Tyr in conantokin-G results in a reduction on the conantokin-G selectivity for GluN2B [18, 52].
To determine whether the amino acid substitutions in EAR16 results in selectivity for only the GluN2B, additional studies are required measuring its selectivity for GluN2B over the other two GluN2 subunits, the GluN2C and GLuN2D.
With respect to the predicted free energy of binding, one would expect that EAR18 (−7.7 kcal/mol) would have a higher affinity and potentially a higher interaction with GluN2B than EAR16 (−6.6 kcal/mol). However, their affinity is comparable, and EAR16 shows an apparent higher selectivity for GluN2B than GluN2A, than EAR18. Hence, neither the lower free energy of binding or the higher number of interactions of EAR18 over EAR 16 predict selectivity for GluN2B with respect to GluN2A.
The IC50 value for EAR16 and EAR18 is high as compared to that displayed by conantokin-G. The reported IC50 for conantokin-G is 0.48 µM in hippocampal  and cortical  neurons, and 0.1 µM in Xenopus oocytes expressing the GluN1a–GluN2B subunits . Hence our results with EAR16 and EAR18 suggest that the Gla residue contributes to the high affinity of conantokin-G.
This study has several limitations. First, we did not measure whether the peptide block was voltage-dependent; neither whether the peptide block was affected by the presence of physiological levels of Mg2+. Second, we measured the selectivity of these peptides between GluN2B and GluN2A. This was done in part because the amino acid sequence corresponding to that used for the GluN2B model (1-818 amino acids) is conserved (71.9% identity) and of the LBD (44 amino acids) is highly conserved (95.5% identity) between GluN2A and BluN2B [25, 41, 42]. However, although to a lower extent, the corresponding sequences in the GluN2B subunit also show conservation with the GluN2C (1-818 amino acids, 60.2% identity; LBD 90.9% identity) and with the GluN2D (1-818 amino acids, 57.6% identity, LBD 93.2% identity) subunits . Hence determining the potency of these peptides on recombinant NMDAr containing other GluN2 and GluN3 subunits, would provide additional information about their level of selectivity which will help define their potential use as pharmacological agents.
We found that EAR16 and EAR18 block in a reversible manner NMDA-evoked currents in hippocampal neurons. Moreover, EAR16 showed selectivity for the GluN2B over GluN2A. The high reversibility of both of these peptides and the selectivity of EAR16 for GluN2B, make them attractive as potential pharmacological agents.
half maximal inhibitory concentration
maximal inhibitory effect
inotropic glutamate receptors
amino terminal domain
agonist/ligand binding domain
ligand binding domain
carboxyl terminal domain
- MS (M + H+):
hydrogen cation mass spectrometry analysis
- MS (M + Na+):
sodium ion mass spectrometry analysis
University of California San Francisco
high-performance liquid chromatograph
Peptide Tertiary Structure Prediction Server
analysis of variance
Barria A, Malinow R. NMDA receptor subunit composition controls synaptic plasticity by regulating binding to CaMKII. Neuron. 2005;48(2):289–301.
Gardoni F, Polli F, Cattabeni F, Di Luca M. Calcium-calmodulin-dependent protein kinase II phosphorylation modulates PSD-95 binding to NMDA receptors. Eur J Neurosci. 2006;24(10):2694–704.
Fan X, Jin WY, Wang YT. The NMDA receptor complex: a multifunctional machine at the glutamatergic synapse. Front Cell Neurosci. 2014;8:160.
Paoletti P, Neyton J. NMDA receptor subunits: function and pharmacology. Curr Opin Pharmacol. 2007;7(1):39–47.
Carvajal FJ, Mattison HA, Cerpa W. Role of NMDA receptor-mediated Glutamatergic signaling in chronic and acute neuropathologies. Neural Plast. 2016;2016:2701526.
Laube B, Kuhse J, Betz H. Evidence for a tetrameric structure of recombinant NMDA receptors. J Neurosci. 1998;18(8):2954–61.
Matsuda S, Kamiya Y, Yuzaki M. Roles of the N-terminal domain on the function and quaternary structure of the ionotropic glutamate receptor. J Biol Chem. 2005;280(20):20021–9.
Sobolevsky AI, Rosconi MP, Gouaux E. X-ray structure, symmetry and mechanism of an AMPA-subtype glutamate receptor. Nature. 2009;462(7274):745–56.
Low CM, Wee KS. New insights into the not-so-new NR3 subunits of N-methyl-d-aspartate receptor: localization, structure, and function. Mol Pharmacol. 2010;78(1):1–11.
Perez-Otano I, Larsen RS, Wesseling JF. Emerging roles of GluN3-containing NMDA receptors in the CNS. Nat Rev Neurosci. 2016;17(10):623–35.
Paoletti P, Bellone C, Zhou Q. NMDA receptor subunit diversity: impact on receptor properties, synaptic plasticity and disease. Nat Rev Neurosci. 2013;14(6):383–400.
Schumann J, Alexandrovich GA, Biegon A, Yaka R. Inhibition of NR2B phosphorylation restores alterations in NMDA receptor expression and improves functional recovery following traumatic brain injury in mice. J Neurotrauma. 2008;25(8):945–57.
Dingledine R, Borges K, Bowie D, Traynelis SF. The glutamate receptor ion channels. Pharmacol Rev. 1999;51(1):7–61.
Furukawa H, Singh SK, Mancusso R, Gouaux E. Subunit arrangement and function in NMDA receptors. Nature. 2005;438(7065):185–92.
Yao Y, Mayer ML. Characterization of a soluble ligand binding domain of the NMDA receptor regulatory subunit NR3A. J Neurosci. 2006;26(17):4559–66.
Paoletti P, Ascher P, Neyton J. High-affinity zinc inhibition of NMDA NR1-NR2A receptors. J Neurosci. 1997;17(15):5711–25.
Williams K. Ifenprodil discriminates subtypes of the N-methyl-d-aspartate receptor: selectivity and mechanisms at recombinant heteromeric receptors. Mol Pharmacol. 1993;44(4):851–9.
Sheng Z, Prorok M, Castellino FJ. Specific determinants of conantokins that dictate their selectivity for the NR2B subunit of N-methyl-d-aspartate receptors. Neuroscience. 2010;170(3):703–10.
Donevan SD, McCabe RT. Conantokin G is an NR2B-selective competitive antagonist of N-methyl-d-aspartate receptors. Mol Pharmacol. 2000;58(3):614–23.
Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.
London N, Raveh B, Cohen E, Fathi G, Schueler-Furman O. Rosetta FlexPepDock web server-high resolution modeling of peptide-protein interactions. Nucleic acids research. 2011;39(Web Server issue):W249–53.
Reyes-Guzman E, Vega-Castro NA, Reyes-Montano EA. Diseño de péptidos inhibidores de interacciones de la subunidad GluN2B del receptor NMDA en isquemia. Rev Cienc Tecnol (RECyT). 2017;27:ISSN:1851–7587 (in press).
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008;9:40.
Traynelis SF, Wollmuth LP, McBain CJ, Menniti FS, Vance KM, Ogden KK, Hansen KB, Yuan H, Myers SJ, Dingledine R. Glutamate receptor ion channels: structure, regulation, and function. Pharmacol Rev. 2010;62(3):405–96.
Laskowski RA, Hutchinson EG, Michie AD, Wallace AC, Jones ML, Thornton JM. PDBsum: a Web-based database of summaries and analyses of all PDB structures. Trends Biochem Sci. 1997;22(12):488–90.
Morris AL, MacArthur MW, Hutchinson EG, Thornton JM. Stereochemical quality of protein structure coordinates. Proteins. 1992;12(4):345–64.
Raveh B, London N, Schueler-Furman O. Sub-angstrom modeling of complexes between flexible peptides and globular proteins. Proteins. 2010;78(9):2029–40.
Longart M, Liu Y, Karavanova I, Buonanno A. Neuregulin-2 is developmentally regulated and targeted to dendrites of central neurons. J Comp Neurol. 2004;472(2):156–72.
Beaudoin GM 3rd, Lee SH, Singh D, Yuan Y, Ng YG, Reichardt LF, Arikkath J. Culturing pyramidal neurons from the early postnatal mouse hippocampus and cortex. Nat Protoc. 2012;7(9):1741–54.
Luo JH, Fu ZY, Losi G, Kim BG, Prybylowski K, Vissel B, Vicini S. Functional expression of distinct NMDA channel subunits tagged with green fluorescent protein in hippocampal neurons in culture. Neuropharmacology. 2002;42(3):306–18.
Hamill OP, Marty A, Neher E, Sakmann B, Sigworth FJ. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflugers Arch. 1981;391(2):85–100.
Shelley C. Single-channel recording of glutamate receptors. Curr Protoc Pharmacol. 2015;68:11.16.1–19. doi:10.1002/0471141755.ph1116s68.
Huang L, Balsara RD, Sheng Z, Castellino FJ. Conantokins inhibit NMDAR-dependent calcium influx in developing rat hippocampal neurons in primary culture with resulting effects on CREB phosphorylation. Mol Cell Neurosci. 2010;45(2):163–72.
Mayer ML, Westbrook GL, Guthrie PB. Voltage-dependent block by Mg2+ of NMDA responses in spinal cord neurones. Nature. 1984;309(5965):261–3.
Nowak L, Bregestovski P, Ascher P, Herbet A, Prochiantz A. Magnesium gates glutamate-activated channels in mouse central neurones. Nature. 1984;307(5950):462–5.
Loftis JM, Janowsky A. The N-methyl-d-aspartate receptor subunit NR2B: localization, functional properties, regulation, and clinical implications. Pharmacol Ther. 2003;97(1):55–85.
Thomas CG, Miller AJ, Westbrook GL. Synaptic and extrasynaptic NMDA receptor NR2 subunits in cultured hippocampal neurons. J Neurophysiol. 2006;95(3):1727–34.
Tovar KR, Westbrook GL. Mobile NMDA receptors at hippocampal synapses. Neuron. 2002;34(2):255–64.
Fischer G, Mutel V, Trube G, Malherbe P, Kew JN, Mohacsi E, Heitz MP, Kemp JA. Ro 25-6981, a highly potent and selective blocker of N-methyl-d-aspartate receptors containing the NR2B subunit. Characterization in vitro. J Pharmacol Exp Ther. 1997;283(3):1285–92.
Wyllie DJ, Livesey MR, Hardingham GE. Influence of GluN2 subunit identity on NMDA receptor function. Neuropharmacology. 2013;74:4–17.
Ishii T, Moriyoshi K, Sugihara H, Sakurada K, Kadotani H, Yokoi M, Akazawa C, Shigemoto R, Mizuno N, Masu M, et al. Molecular characterization of the family of the N-methyl-d-aspartate receptor subunits. J Biol Chem. 1993;268(4):2836–43.
Karakas E, Furukawa H. Crystal structure of a heterotetrameric NMDA receptor ion channel. Science. 2014;344(6187):992–7.
Sheng Z, Liang Z, Geiger JH, Prorok M, Castellino FJ. The selectivity of conantokin-G for ion channel inhibition of NR2B subunit-containing NMDA receptors is regulated by amino acid residues in the S2 region of NR2B. Neuropharmacology. 2009;57(2):127–36.
Klein RC, Galdzicki Z, Castellino FJ. Inhibition of NMDA-induced currents by conantokin-G and conantokin-T in cultured embryonic murine hippocampal neurons. Neuropharmacology. 1999;38(12):1819–29.
Kaas Q, Yu R, Jin AH, Dutertre S, Craik DJ. ConoServer: updated content, knowledge, and discovery tools in the conopeptide database. Nucleic Acids Res. 2012;40((Database issue)):325–30.
Cnudde SE, Prorok M, Dai Q, Castellino FJ, Geiger JH. The crystal structures of the calcium-bound con-G and con-T[K7γ] dimeric peptides demonstrate a metal-dependent helix-forming motif. J Am Chem Soc. 2007;129(6):1586–93.
Skjaerbaek N, Nielsen KJ, Lewis RJ, Alewood P, Craik DJ. Determination of the solution structures of conantokin-G and conantokin-T by CD and NMR spectroscopy. J Biol Chem. 1997;272(4):2291–9.
Chandler P, Pennington M, Maccecchini ML, Nashed NT, Skolnick P. Polyamine-like actions of peptides derived from conantokin-G, an N-methyl-d-aspartate (NMDA) antagonist. J Biol Chem. 1993;268(23):17173–8.
Lin SY, Wu K, Len GW, Xu JL, Levine ES, Suen PC, Mount HT, Black IB. Brain-derived neurotrophic factor enhances association of protein tyrosine phosphatase PTP1D with the NMDA receptor subunit NR2B in the cortical postsynaptic density. Brain Res Mol Brain Res. 1999;70(1):18–25.
Teichert RW, Jimenez EC, Twede V, Watkins M, Hollmann M, Bulaj G, Olivera BM. Novel conantokins from Conus parius venom are specific antagonists of N-methyl-d-aspartate receptors. J Biol Chem. 2007;282(51):36905–13.
Klein RC, Prorok M, Galdzicki Z, Castellino FJ. The amino acid residue at sequence position 5 in the conantokin peptides partially governs subunit-selective antagonism of recombinant N-methyl-d-aspartate receptors. J Biol Chem. 2001;276(29):26860–7.
EARG prepared the hippocampal cultures, performed the electrophysiological measurements, data analysis, contributed to making the figures and in the preparation of the manuscript. NVC contributed in training EARG in cell culture and data discussions; EARM contributed in data discussions, responsible for writing the manuscript. ERP trained and directly supervised the work of EARG, performed the transfections with recombinant NMDAr subunits, responsible for writing the manuscript, contributed to making the figures. All authors read and approved the final manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We are grateful to Dr. Monica Norcini for growing the plasmids used and to Dr. Jhon Jairo Sutachan for teaching us the preparation of hippocampal cultures.
EAR16 and EAR18 are relevant for patent application. Both fund sources (indicated under funding) contributed to the design of the study, collection, analysis, and interpretation of the data, and in writing the manuscript.
Availability of data and materials
The dataset during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Pregnant rats were maintained and used following the guidelines approved by the New York University Langone Medical Center Institutional Animal Care and Use Committee (Protocol Number: 161015-01; ERP).
Support was provided by The Anesthesia Research Fund of the New York University, Department of Anesthesiology (ERP), and COLCIENCIAS by the Grant to the Project No. 110156933188 entitled “Diseño de péptidos basados en la interacción de GluN2B con DAPK1, SRC, D2R y Conantokina G y evaluación in vitro de su posible efecto neuroprotector en eventos isquémicos”, interadministrative Contract No. 476-2013 (COLCIENCIAS-Universidad Nacional de Colombia) (NVC, EARM). Both fund sources contributed to the design of the study, collection, analysis, and interpretation of the data, and in writing the manuscript.
About this article
Cite this article
Reyes-Guzman, E.A., Vega-Castro, N., Reyes-Montaño, E.A. et al. Antagonistic action on NMDA/GluN2B mediated currents of two peptides that were conantokin-G structure-based designed. BMC Neurosci 18, 44 (2017) doi:10.1186/s12868-017-0361-4
- Hippocampal neurons