- Open Access
In silico analysis of binding interaction of conantokins with NMDA receptors for potential therapeutic use in Alzheimer’s disease
Journal of Venomous Animals and Toxins including Tropical Diseasesvolume 23, Article number: 42 (2017)
The N-methyl-D-aspartate (NMDA) receptors are glutamate receptors that play vital roles in central nervous system development and are involved in synaptic plasticity, which is an essential process for learning and memory. The subunit N-methyl D-aspartate receptor subtype 2B (NR2B) is the chief excitatory neurotransmitter receptor in the mammalian brain. Disturbances in the neurotransmission mediated by the NMDA receptor are caused by its overexposure to glutamate neurotransmitter and can be treated by its binding to an antagonist. Among several antagonists, conantokins from cone snails are reported to bind to NMDA receptors.
This study was designed to analyze the binding mode of conantokins with NMDA receptors in both humans and rats. To study interactions, dockings were performed using AutoDock 4.2 and their results were further analyzed using various computational tools.
Detailed analyses revealed that these ligands can bind to active site residues of both receptors as reported in previous studies.
In light of the present results, we suggest that these conantokins can act as antagonists of those receptors and play an important role in understanding the importance of inhibition of NMDA receptors for treatment of Alzheimer’s disease.
The N-methyl-D-aspartate (NMDA) receptors are inotropic glutamate receptors that are gated cation channels [1, 2]. The NMDA receptors (NMDAR) play vital roles in central nervous system (CNS) development . These receptors are highly permeable to Ca2+ ions and the calcium flux is critically important for synaptic plasticity, which is an essential neurochemical process for learning and memory [2,3,4]. The receptor itself has many subunits and their variants have numerous functions in the brain. The subunit N-methyl-D-aspartate receptor subtype 2B (NR2B) is the chief excitatory neurotransmitter receptor in the mammalian brain . The glutamate neurotransmitter allows for a transmembrane ion flow through the receptor to increase the action potential of the neuron. This characteristic makes the synapsis among these neurons to be the main memory storage unit and hence associates them with learning and memory . Due to their functioning in the CNS, the potential of these receptors as drug targets for various neurodegenerative diseases has been highlighted in the literature.
NMDAR antagonists have emerged as potential lead compounds for Alzheimer’s patients . The cognitive symptoms associated with deficits in learning and menory have been attributed to disturbances in glutaminergic neurotransmission . The excessive stimulation by the glutamate neurotransmitter of neurons causes excitotoxicity and results in damage and death of neurons . Blocking the glutaminergic neurotransmission mediated by NMDA receptors can alleviate the excitotoxicity and prevent further neuron damage and death.
Several venom toxins have made their way in scientific studies and clinical trials for their therapeutic potential against various diseases. Having mostly inhibitory effects, these toxin peptides target different receptors across the body, the hyperactivity of which is associated with the pathophysiology of many diseases. Numerous toxins have been reported to target and block receptors used for the treatment of Alzheimer’s disease, such as neurotoxins that target acetylcholinesterase enzymes and certain toxins from Mamba snakes that were reported to inhibit the muscarinic acetylcholine receptors [9, 10]. A few toxins with antagonistic properties against NMDA receptor have been reported as well. The conantokins are powerful and potent blockers of the NMDA receptor, with particularly high selectivity for the NR2B subunit . These toxins have been reported to block these glutamine receptors and therefore have therapeutic potential for treating Alzheimer’s disease.
The objective of this study was to perform and look into the in silico analysis of the binding interaction of conantokins with the NMDA receptor NR2B subunit. The mode of interaction and the binding residues for both the ligand dataset and the receptor dataset were collected. Due to unavailability of the crystal structures of the NMDA receptor in humans and most of conantokins, their three-dimensional structures were predicted via computational homology modeling methods and the predicted models were validated to continue their further use.
Docking studies provided insights into binding pattern of receptors and ligands. A number of in silico studies investigated many computational approaches ranging from construction of structural models to investigation and discovery of potential drug candidates [12,13,14]. The analysis of the binding interactions of the receptor and the ligand peptides produced results that helped us demonstrate the pharmacological importance of conantokins and their potential use as NMDA receptor antagonists for treatment of Alzheimer’s disease.
Receptor dataset collection
The first step of methodology included collection of receptor proteins. As the structure of NMDA receptor in humans is not available, it had to be predicted via computational homology modeling, which allowed the construction of a three-dimensional structure of a protein based on the known structures of similar protein templates. SWISS-MODEL was used for this purpose, as it is a fully automated protein structure homology modeling server . The protein template used for structure prediction of NMDAR in humans was the structure of NMDAR in Rattus norvegicus (brown rat) . The structure of the template (pdb id: 3JPW) was retrieved from the Research Collaboratory for Structural Bioinformatics (RSCB) . After structure prediction, the Structure Analysis and Verification Server (SAVES) was used for validation by generating Ramachandran plot  and ERRAT , which gives a factor of overall quality of the predicted structure. Whereas Verify_3D  was employed to analyze the compatibility of the atomic model of the protein with its own amino acid sequence.
Ligand dataset collection
The ligand dataset comprised reported conantokins that are a class of conopeptides (17–27 amino acids) without cysteine residues that selectively influence NMDA receptors . Among them, only three-dimensional structures of conantokin G and conantokin T are available on Protein Data Bank (PDB ID: 1ONU and 1ONT, respectively). The crystal structures of other conantokins Br, L, P, R, E, Pr1, Pr2, Pr3, R1A, R1B and R1C are not available on the Protein Data Bank and therefore were subjected to three-dimensional structure prediction. Homology modeling was used initially, but due to absence of any homologues for the conantokins, the structure prediction was carried out via fold recognition (threading) on the Iterative Threading ASSEmbly Refinement (iTASSER) server , which detects structure templates from the Protein Data Bank and constructs full-length structure models by reassembling structural fragments from threading templates. The predicted structures were then validated on the SAVES metaserver using Procheck, ERRAT and Verify_3D. Multiple sequence alignment (MSA) using PRALINE  multiple sequence alignment toolbox was performed on all the conantokins to find out residues and secondary structure conservation. The residues responsible for the signal peptides, peptide precursors and active protein peptides were also highlighted for each conantokin showing considerable conservation in these peptide regions as well, as shown in Fig. 1.
The conserved glutamic acid residues in each conantokin are the positions for residue modification, where the glutamic acid is converted to 4-carboxyglutamate. Sufficient conservation of secondary structures was observed among all the conantokins.
Experimental background reported for conantokins and NMDAR
Conantokins have been the subject of interest for their potential as NMDA receptor antagonists in studying many neuropathologies. Wet lab experiments have been conducted to understand the basis for interaction between conantokins and NMDA receptor. These wet lab techniques were carried out using NMDA receptors in rats by employing electrophysiological techniques to understand the diversity in the functioning of these toxins . Techniques such as polyamine enhancement and NMR spectroscopy have allowed the understanding of how similar the inhibition of NMDAR by conantokins in humans is to previously reported inhibition of NMDAR in rats by conantokins . Using chemically synthesized conantokin variants (created by techniques of point mutation and chimeric proteins), the subtype selectivity of NMDAR has been extensively studied for most conantokins . In addition, conantokins represent a class of NMDA antagonists with an improved safety profile compared to other antagonists that show psychotomimetic, amnesic and motor-impairing actions and neurotoxicity, which limit their usefulness in humans . These wet lab studies revealed the pharmacological importance of conantokins as antagonists of NMDAR for their therapeutic applications in many neuropathologies and have led us to investigate this property computationally.
Molecular dockings are performed to predict the binding orientation between a receptor and its ligand to form a stable molecular complex . This allows a detailed insight of the three-dimensional structure formed between two biomolecules and to correlate this information to find a potential drug candidate for these receptors . Molecular docking was performed for NMDA receptors for both humans and rats. AutoDock 4.2  was used to perform automated docking runs, in order to find the binding mode of each receptor with each conantokin ligand. The number of runs for each docking was set to 50 to allow each ligand 50 different conformations with the receptors in order to let them bind freely anywhere on the receptor. The grid size was set to cover the entire receptor in order to find a potential binding site for each ligand and to analyze if each ligand indeed occupies the site on the receptor that is crucial for its functioning in Alzheimer’s and could potentially inhibit it. The docking parameters used are shown in Table 1.
Post-docking analysis of the results was carried out on AutoDock 4.2 and later visualized on Chimera . Ligplot+  was used for further validation of the protein-protein interactions, which generates two-dimensional schematics on the basis of hydrogen bonds and hydrophobic interactions. The three-dimensional structural analysis was performed on PyMol  that allowed the residues in the protein-protein interaction of the receptor-ligand complex to be highlighted. The energy of each complex that was achieved via docking allowed the receptor-ligand bindings to be confirmed and refined and used for further analysis and results discussion.
Results and discussion
Structure prediction and validation of NMDAR in human
The predicted three-dimensional structure of the NMDA receptor in humans was superimposed on Chimera with its template, a NMDA receptor in rat yielding a root mean square deviation (RSMD) of 0.181 Å. In addition, the results for structure evaluation proved the quality of the predicted structure. As indicated by the Ramachandran plot, 92% of the residues were in favored regions. Moreover, the main chain parameters such as peptide bond planarity, measure of non-bonded interactions, α carbon tetrahedral distortion, H-bond energy and overall G factor for the structure were found within favorable regions. ERRAT and Verify_3D further validated the structure by scoring it with an overall quality factor of 79.155 and by passing it with 80% of the amino acids having scored ≥0.2 respectively (Fig. 2). The validation results suggested that the predicted model of the NMDA receptor in humans was of good quality and can be used for further analysis.
Structure prediction and validation of conantokins
The predicted structure of each conantokin was also validated in order to check the quality of their three-dimensional structures. All predicted structures were passed by Ramachandran plot, ERRAT and Verify_3D, suggesting that these structures were of good quality and could be used further study.
Binding site residue information
After structure prediction and evaluation, binding site residue information for NMDA receptor in humans was gathered. Since this is a computational based analysis, it is very important to identify binding site residues correctly and verify the results. As no previous information for binding residues of NMDAR in humans was reported, the binding site was retrieved by using the binding site residues of NMDAR reported in the literature that showed residues specific for the NR2B subunit in rats. The two structures were aligned and it was observed that the binding site residues reported for NMDAR in rats are conserved in humans. The binding residues for NMDAR in both humans and rats are shown in Fig. 3. These identified residues have been taken as references to further evaluate the docking results.
Binding residue information for conantokin G only suggests Glu2, Gln6, Asn8, Arg13, Asn17 and Lys15 as predominant residues that are involved in interactions with the NMDAR subtype selectivity for NR2B subunit. Moreover, the first five amino acids of conantokins, especially Glu2 and Gla4 and a hydrophobic residue at position 12 are critical for functional activity in vitro . It has been observed that the peptide region has shown importance in interactions. These residues are also found to be conserved in other conantokins as well.
After structure prediction and binding site identification, molecular dockings were performed on NMDAR receptors with conantokin ligands. Docking results were run based on their energy values and were ranked as such with each run showing results for receptor-ligand complex with the lowest energy. Detailed analysis of each run for each complex showed that all the conantokins showed bindings with the conserved active sites of NMDA receptor in both humans and rats. All the conantokins seemed to bind with the glutamine and glutamic acid residues in the NMDA receptors in both humans and rats. Table 2 shows the detailed analysis retrieved from plotting the docking results on LigPlot+, showing the binding of NMDA receptor in humans with conantokins. These data include information about respective hydrogen bond residues of both receptors and ligands, their bond distances, bond atoms and the interacting hydrophobic residues.
Glutamic acid at position 236 and glutamine at position 110 were predominantly involved in the binding of the NR2B subunit of NMDAR in human with all the conantokins. Isoleucine at position 111, phenylalanine at position 114 and proline at position 177 were mostly found in hydrophobic interactions. The binding patterns of all the conantokin ligands individually in complex with the NR2B subunit of NMDAR in human are shown in Fig. 4.
It is clear in Fig. 4 that due to difference in sizes and structure, each conantokin occupies the binding site on NMDAR in humans in its own respective orientation to form the most stable complex. However, each ligand is shown to bind to the same binding residues as reported. The position of binding pocket of the NR2B subunit of NMDAR in humans resided by conantokins is shown in Fig. 5. It is shown clearly that all the ligands occupied the same binding pocket in structure of NMDAR as they showed binding with the same residues that have been reported.
Similar binding pattern was observed for the NR2B subunit of NMDAR in rat. Table 3 shows the binding information for NMDAR in rats, highlighting residues involved in hydrogen bonding, their atoms and their bond distances, along with the hydrophobic residues for both the receptor and the ligands.
It has been observed that likewise for human NR2B, Glu236 and Gln110 residues of NR2B subunit in rats were found interacting with the residues of every conantokin ligand. Ile111, Phe114 and Pro177 were mostly found in hydrophobic interactions. This binding pattern was conserved among all the conantokins showing their similarity in potency for the NR2B subunit. The binding pattern of all conantokins for NMDAR in rats is shown in Fig. 6.
The binding pocket of conantokins with the receptor is highlighted in Fig. 7. As in the case of NMDAR in rats, Fig. 7 clearly demonstrates that each conantokin occupies the same binding pocket as NMDAR in humans.
The overall stability of a complex is directly associated with the free energy of that complex. Lower energy values suggest the presence of a high binding affinity between the ligand and the receptor. It is of substantial importance to highlight each complex’s free energy in order to evaluate the successful complex formation. The successive energy values for the ligand-receptor complex formed by each toxin with the NMDA receptor in both humans and rats are shown in Fig. 8.
The energy values for all the complexes fall below the value of 0, as clearly observed from the graphical representation. This suggests the formation of relatively stable complexes between NMDA receptor in both rats and humans with all conantokins. NMDAR in humans showed most affinity for conantokin P, conantokin E, conantokin R and conantokin Pr3, whereas NMDAR in rats showed more affinity for Conantokin Br and Conantokin Pr3.
Comparison of reported residues for NMDAR to docking results
A comparative analysis was performed of the docking results of NMDA receptors in both humans and rats with all the conantokins. The reported residues were inspected to verify if the docking results of this study were indeed in line with the reported information regarding the binding site residues of NR2B subunit of NMDAR. It is clearly displayed in Table 4 that the binding site residues of NMDAR reported in the literature were achieved in the docking experiments for both humans and rats.
Conantokins had been reported as selective antagonists of NR2B subtype of the NMDA receptor. The NR2B subunit is directly involved in the excitotoxicity caused by the over-exposure to glutamate neurotransmitter. Therefore, inhibiting its activity has become a possible target for treatment of Alzheimer’s. We designed an in silico study and analysis of the binding interaction of conantokins with the NMDA receptor in both humans and rats. After the binding site residues for the receptors were understood, docking studies were performed and the residues achieved via docking were compared to the binding residues reported in the literature. The similarity of the results suggests the potential successful binding of the toxin ligands with both the receptors, and their subsequent function as their antagonists. Our results helped us to demonstrate the potential of these receptors in complex with the conantokins for the symptomatic treatment of Alzheimer’s patients.
Central nervous system
Multiple sequence alignment
N-methyl D-aspartate receptor subtype 2B
Research collaboratory for structural bioinformatics
Root mean square deviation
Structure analysis and verification server
Blanke ML, VanDongen AMJ. Activation mechanisms of the NMDA receptor. In: Van Dongen AM, editor. Biology of the NMDA receptor. Boca Raton (FL): CRC Press; 2009.
Li ST, Ju JG. Functional roles of synaptic and extrasynaptic NMDA receptors in physiological and pathological neuronal activities. Curr Drug Targets. 2012;13(2):207–21.
Hunt DL, Castillo PE. Synaptic plasticity of NMDA receptors: mechanisms and functional implications. Curr Opin Neurobiol. 2012;22(3):496–508.
Martin SJ, Grimwood PD, Morris RG. Synaptic plasticity and memory: an evaluation of the hypothesis. Annu Rev Neurosci. 2000;23:649–711.
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.
Yang JL, Sykora P, Wilson DM, Mattson MP, Bohr VA. The excitatory neurotransmitter glutamate stimulates DNA repair to increase neuronal resiliency. Mech Ageing Dev. 2011;132(8–9):405–11.
Farlow MR. NMDA receptor antagonists. A new therapeutic approach for Alzheimer's disease. Geriatrics. 2004;59(6):22–7.
Danysz W, Pasons CG. Alzheimer’s disease, β-amyloid, glutamate, NMDA receptors and memantine--searching for the connections. Br J Pharmacol. 2012;167(2):324–52.
Servent D, Blanchet G, Mourier G, Marquer C, Marcon E, Fruchart-Gaillard C. Muscarinic toxins. Toxicon. 2011;58(6–7):455–63.
Waqar M, Batool S. In silico analysis of binding of neurotoxic venom ligands with acetylcholinesterase for therapeutic use in treatment of Alzheimer’s disease. J Theor Biol. 2015;372:107–17.
Waqar M, Kamal MA, Batool S. In silico analysis of binding interaction of mamba toxins with M4 and M2 muscarinic acetylcholine receptors for therapeutic use in Alzheimer’s disease. CNS Neurol Disord Drug Targets. 2015;14(8):1031–40.
Amin SA, Bhattacharya P, Basak S, Gayen S, Nandy A, Saha A. Pharmacoinformatics study of Piperolactam a from Piper Betle root as new lead for non steroidal anti fertility drug development. Comput Biol Chem. 2017;67:213–24.
Amin SA, Adhikari N, Agrawal RK, Jha T, Gayen S. Possible binding mode analysis of Pyrazolo-triazole hybrids as potential anticancer agents through validated molecular docking and 3D-QSAR modeling approaches. Letters in Drug Design & Discovery. 2017;14(5):515–27. www.eurekaselect.com/145590
Amin SA, Bhargava S, Adhikari N, Gayen S, Jha T. Exploring pyrazolo [3, 4-d] pyrimidine phosphodiesterase 1 (PDE1) inhibitors: a predictive approach combining comparative validated multiple molecular modelling techniques. J Biomol Struct Dyn. 2017;13:1–19.
Schwede T, Kopp J, Guex N, Peitsch MC. SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003;31(3):3381–5.
Karakas E, Simorowski N, Furukawa H. Structure of the zinc-bound amino-terminal domain of the NMDA receptor NR2B subunit. EMBO J. 2009;28(24):3910–20.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42.
Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996;8(4):477–86.
Colovos C, Yeates TO. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993;2(9):1511–9.
Bowie JU, Lüthy R, Eisenberg D. A method to identify protein sequences that fold into a known three-dimensional structure. Science. 1991;253(5016):164–70.
Lebbe EKM, Tytgat J. In the picture: disulfide-poor conopeptides, a class of pharmacologically interesting compounds. J Venom Anim Toxins incl Trop Dis. 2016;22:30. doi: 10.1186/s40409-016-0083-6.
Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5(4):725–38.
Simossis VA, Heringa J. PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Res. 2005;33:W289–94.
Teichert RW, Jimenez EC, Twede V, Watkins M, Hollmann M, Bulaj G, et al. Novel conantokins from Conus parius venom are specific antagonists of N-methyl-D-aspartate receptors. J Biol Chem. 2007;282(51):36905–13.
Nielsen KJ, Skjaerbaek N, Dooley M, Adams DA, Mortensen M, Dodd PR, et al. Structure-activity studies of conantokins as human N-methyl-D-aspartate receptor modulators. J Med Chem. 1999;42(3):415–26.
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.
Lengauer T, Rarey M. Computational methods for biomolecular docking. Curr Opin Struct Biol. 1996;6(3):402–6.
Mondal C, Halder AK, Adhikari N, Saha A, Saha KD, Gayen S, et al. Comparative validated molecular modeling of p53-HDM2 inhibitors as antiproliferative agents. Eur J Med Chem. 2015;90:860–75.
Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexiblity. J Comput Chem. 2009;30(16):2785–91.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995;8(2):127–34.
DeLano WL. Pymol: An open-source molecular graphics tool. CCP4 Newsletter On Protein Crystallography. 2002;40:82–92.
The authors would like to thank the facilities provided by the Department of Biosciences, COMSATS Institute of Information Technology, Pakistan.
Ethics approval and consent to participate
Consent for publication
Both authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
- Synaptic plasticity
- In silico
- Alzheimer’s dieases