twitter

Saturday, 5 January 2019

Curare alkaloids from Matis Dart Poison: Comparison with d-tubocurarine in interactions with nicotinic, 5-HT3 serotonin and GABAAreceptors

Abstract

Several novel bisbenzylisoquinoline alkaloids (BBIQAs) have recently been isolated from a Matis tribe arrow poison and shown by two-electrode voltage-clamp to inhibit mouse muscle nicotinic acetylcholine receptors (nAChR). Here, using radioligand assay with Aplysia californica AChBP and radioiodinated α-bungarotoxin ([125I]-αBgt), we show that BBIQA1, BBIQA2, and d-tubocurarine (d-TC) have similar affinities to nAChR orthosteric site. However, a competition with [125I]-αBgt for binding to the Torpedo californica muscle-type nAChR revealed that BBIQAs1, 2, and 3 are less potent (IC50s = 26.3, 8.75, and 17.0 μM) than d-TC (IC50 = 0.39 μM), while with α7 nAChR in GH4C1 cells, BBIQA1 was less potent that d-TC (IC50s = 162 μM and 7.77 μM, respectively), but BBIQA2 was similar (IC50 = 5.52 μM). In inhibiting the Ca2+responses induced by acetylcholine in Neuro2a cells expressing the mouse adult α1β1εδ nAChR or human α7 nAChR, BBIQAs1 and 2 had similar potencies to d-TC (IC50s in the range 0.75–3.08 μM). Our data suggest that BBIQA1 and BBIQA2 can inhibit adult muscle α1β1εδ nAChR by both competitive and noncompetitive mechanisms. Further experiments on neuronal α3β2, α4β2, and α9α10 nAChRs, expressed in Xenopus laevis oocytes, showed that similar potencies for BBIQAs1, 2, and d-TC. With α3β2γ2 GABAAR currents were almost completely inhibited by d-TC at a high (100 μM) concentration, but BBIQAs1 and 2 were less potent (only 40–50% inhibition), whereas in competition with Alexa Fluor 546-α-cobratoxin for binding to α1β3γ2 GABAAR in Neuro2a cells, d-TC and these analogs had comparable affinities. Especially interesting effects of BBIQAs1 and 2 in comparison with d-TC were observed for 5-HT3AR: BBIQA1 and BBIQA2 were 5- and 87-fold less potent than d-TC (IC50 = 22.63 nM). Thus, our results reveal that these BBIQAs differ from d-TC in their potencies towards certain Cys-loop receptors, and we suggest that understanding the reasons behind this might be useful for future drug design.

Introduction

The nicotinic acetylcholine (nAChR), glycine (GlyR), serotonin (5-HT3R), and γ-aminobutyric acid (GABAAR) receptors belong to the Cys-loop family of ligand-gated ion channels (LGICs) [15]. All members of this family have a similar structure: five subunits form a pentamer surrounding a selective ion channel; for nAChR and 5-HT3R, the ion channel is permeable to cations, while in GABAAR and GlyR it is anion-permeable. Each subunit is composed of an N-terminal extracellular domain, four transmembrane fragments (M1-M4), and a large cytoplasmic loop between M3 and M4 fragments. In the N-terminal extracellular domain, there is a highly conserved 13 amino acid loop tethered at each end by a cysteine disulfide bridge, hence the name “Cys-loop receptors” [68].
Ligands for Cys-loop receptors are very diverse and include both synthetic and natural compounds. Typically, agonists and competitive antagonists bind to the orthosteric sites located between two adjacent subunits, whereas allosteric binders target different sites [911]. Cys-loop receptors are widely expressed in the central and peripheral nervous systems, where they mediate fast synaptic neurotransmission, and also in other systems, for example, in immune and epithelial cells [1216]. Drugs that act at Cys-loop receptors are used as muscle relaxants, chronic pain therapeutics, anti-emetics, neurodegenerative and psychiatric therapeutics, smoking cessation aids and new compounds are being actively sought [1720]. Muscle relaxants and/or neuromuscular blocking agents are widely used in anesthesia, where they act by interrupting transmission at the neuromuscular junction [21]. The development of clinical muscle relaxants began with the examination of the arrow poison produced from plant Chondodenron tomentosum by South American Indians [2122] and by isolation of d-tubocurarine (d-TC) [23]. d-TC is an antagonist of the end-plate muscle type nAChRs, which are composed by two α1, one β1, one δ, and either γ (fetal) or ɛ (adult) subunits. The orthosteric α/ε (or α/γ) site binds d-TC with 100-fold higher affinity than the orthosteric α/δ site [24]. In the middle of the last century, d-TC was a popular choice in surgery, although there were some adverse effects on blood pressure and the cardiovascular system [2225]. Today the compounds used are much safer, but still so there are some side effects and there is a strong need in understanding the key structural elements important for high selectivity of the muscle relaxants. In particular, in addition to high selectivity, the perfect relaxant should be characterized by rapid onset and short or intermediate duration of action [212627].
We have recently described an analysis of a poison which the Matis tribe in South America is using for anointing their arrows [28]. In addition to a number of well-known compounds, such as magniflorine, lindoldhamine and some others [2930], a series of novel bisbenzyltetrahydroisoquinoline alkaloids (BBIQAs) (Fig 1) were discovered. Due to their structural similarity with d-TC, it was reasonable to expect from them the inhibition of neuromuscular transmission. Indeed, it was found that both the crude venom (at a concentration of 1.0 mg/mL), and individual BBIQAs at the relatively high concentration of 20 μM, efficiently decreased the amplitudes of ACh-induced currents in mouse muscle nAChR heterologously expressed in Xlaevis oocytes [28].
thumbnail
Fig 1. Chemical structures of d-tubocurarine (d-TC) and bisbenzyltetrahydroisoquinoline alkaloids (BBIQAs).
The aim of the present work was a detailed in vitro comparison of d-TC and BBIQAs activities on a set of muscle-type and neuronal nAChRs. Previous work has shown that d-TC is able to block nAChRs competitively and noncompetitively, depending on the nAChR subtype: inhibition of neuronal α4β2 nAChR is competitive, while muscle-type and neuronal α3β4 nAChRs are blocked by a mixed mechanism involving binding to both orthosteric and allosteric sites [31]. d-TC binds with a high affinity to acetylcholine-binding protein (AChBP), a good model for studying agonist and competitive antagonist actions on nAChR [32]. Interestingly, d-TC adopts three different orientations in the AChBP ligand binding sites [33], providing a possible reason for the wide range of receptors targeted by d-TC: in addition to nAChRs, d-TC is an effective blocker of 5-HT3R (with IC50 values of 11.4–13 nM and 1.8 μM for mouse and human 5-HT3AR, respectively [3435]), and at higher concentrations (30 μM) it also inhibits GABAAR [36]. In the current research we compared the effects of d-TC and the novel BBIQAs on muscle-type and neuronal nAChRs, as well as on 5-HT3AR and GABAAR. Our data suggest that BBIQAs do not completely mimic the effects of d-TC and thus could help in the identification of structural motifs responsible for selectivity of distinct compounds, which is valuable for future drug design.

Materials and methods

Materials

The Matis poison (5.29 g) was obtained in Leticia, at the Colombian-Brazilian border, by Dr. Rainer W. Bussmann. The material was authenticated and deposited at the William L. Brown Center, Missouri Botanical Garden, in November 2011. Powdered material was maintained at room temperature (22–25°C) and protected from light until required for extraction and analysis. Tcalifornica electric organ membranes were a generous gift from Ferdinand Hucho (Free University of Berlin, Berlin, Germany) and GH4C1 cells–from EliLilly (London, UK). Mouse neuroblastoma Neuro2a cells were purchased from the Russian collection of cell cultures (Institute of Cytology, Russian Academy of Sciences, Saint Petersburg, Russia). Mature Xlaevis female frogs used in this study were obtained commercially (NASCO, Fort Atkinson WI, USA) and housed in a facility with 12:12 hours light:dark cycles, 18–20°C ambient temperature. Animals were fed twice a week and maintained according to supplier recommendations (https://www.enasco.com/page/xen_care). This study was carried out in strict accordance with the World Health Organization’s International Guiding Principles for Biomedical Research Involving Animals. The protocol was approved by the Institutional Policy on the Use of Laboratory Animals of the Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS (Protocol Number: 251/2018 26.02.18). All frog surgery was performed under benzocaine anesthesia, and all the appropriate actions were taken to minimize discomfort to frogs.

Samples purification

Preparative purification was carried out with a Gilson HPLC system (322 pump with GX 271 liquid handler) equipped with an RP column and 155 UV-Vis detector, set at 210 and 254 nm. Compounds were eluted with a H2O-MeCN gradient with 0.1% CF3COOH. UPLC-MS analysis was performed using Thermo Finnigan LCQ Deca XP Plus ion trap instrument with Thermo Accela UPLC system equipped with YMC Triart column (C-18 150 × 2 mm, 1.9 μm). Detection was achieved by UV-VIS DAD (190–600 nm) and full scan MS (ESI+, 150–2000 au). Crude samples were dissolved in mixture of water/methanol/acetic acid 88:10:2 to a final concentration of 3 mg/mL, filtered through a 45 μm nylon filter and injected into the LC system. Chromatography was carried out using a Phenomenex Luna C18(2) 5 μm 21.2x150 mm column in a linear gradient from 9 to 55% of acetonitrile in 15 Column Volumes. Fractions of interest were collected, analyzed by UPLC-MS and lyophilized.

Computer modeling with Acalifornica AChBP

Docking experiments were performed using Autodock 4.2 and analyzed with MGL Tools 1.5.6 [37] with parameters set as follows: genetic algorithm population size 150, number of evaluations 25000000 and number of runs 100. The structure of Acalifornica AChBP co-crystallized with d-TC was used as a receptor (PDB 2XYT) which contains d-TC in three distinct binding positions [33]. To address the possibility of different binding modes we performed docking of BBIQA1 and BBIQA2 to all three types of binding sites. Results were clustered and inspected using MGL Tools 1.5.6 and visualized using UCSF Chimera [38].

Radioligand assay

In competition experiments with [125I]-αBgt, the compounds were pre-incubated 3 h at room temperature with Acalifornica AChBP at final concentrations of 140 nM, Tcalifornica electric organ membranes (final concentration 1.25 nM of toxin-binding sites) or GH4C1 cells (6.5 μg of total protein with final concentration of 0.4 nM of toxin-binding sites) in 50 μL of binding buffer (20 mM Tris-HCl buffer, 1 mg/mL of bovine serum albumin, pH 8.0). [125I]-αBgt was then added to cells or membranes to a final concentration of 0.1–0.2 nM and the mixtures were incubated for 5 min. Binding was stopped by rapid filtration on GF/C filters (Whatman, UK) pre-soaked in 0.25% polyethylenimine, unbound radioactivity being removed from the filters by washing (3×3 mL) with binding buffer. Non-specific binding was determined in all cases using a 3 h pre-incubation with 30 μM α-cobratoxin (αCtx). The results were analyzed using OriginPro 2017 (OriginLab Corporation, Northampton, MA, USA) fitting to a dose-response curve with a variable Hill slope using the equation y = A1 + (A2-A1)/(1 + 10(LOGx0-x)*p): where A1 and A2 are bottom and top asymptotes, respectively; p–Hill slope; LOGx0 –log10[IC50, M].

Cell culturing and transfection

Mouse neuroblastoma Neuro2a cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Paneco, Russia) supplemented with 10% FBS (PAA Laboratories, Austria). Cells were sub-cultured 24 h before transfection and plated at a density of 10,000 cells per well (black 96-well plate, Corning, USA), followed by lipofectamine (Invitrogen, USA)-mediated transient co-transfection of human α7 nAChR-pCEP4, fluorescent calcium sensor pCase12-cyto (Evrogen, Russia) and chaperone Ric3-pCMV6-XL5 or NACHO TMEM35-pCMV6-XL5 plasmid constructs (OriGene, USA). Mouse muscle α1, β1, ε, and δ nAChR-pRBG4 plasmid constructs were expressed similarly, but without a chaperone, as well as mouse α1, β3, γ2 GABAAR-lab-pCI plasmids.

Calcium imaging

Calcium imaging was performed as described previously [39]. Briefly, cell medium was removed and cells were washed with external buffer containing (in mM) 140 NaCl, 2 CaCl2, 2.8 KCl, 4 MgCl2, 20 HEPES, 10 glucose at pH 7.4. Cells were pre-incubated with d-TC, BBIQA1, or BBIQA2 for 20 min at room temperature, or with αBgt for 5 min at room temperature. To potentiate α7 nAChR responses, PNU120596 (10 μM) was added to the pre-incubation solution. Cells were excited at 485 nm and emitted fluorescence was detected at 535±10 nm, using a multimodal microplate reader Hidex Sense (Hidex, Turku, Finland). Fluorescence was recorded every 2 s for 3 min following agonist addition. Responses were measured as peak intensity minus basal fluorescence level, and are expressed as a percentage of a maximal response to agonist. Negative controls were run in the presence of 5 μM αCtx. Data files were analyzed using HidexSence software (Hidex, Turku, Finland) and then the results were analyzed using OriginPro 2017 (OriginLab Corporation, Northampton, MA, USA) fitting to a dose-response curve with a variable Hill slope using the equation y = A1 + (A2-A1)/(1 + 10(LOGx0-x)*p): A1 and A2 are bottom and top asymptotes, respectively; p–Hill slope; LOGx0 –log10[IC50 or EC50, M].

Fluorescence assay

Fluorescence assays were performed as previously described [40]. Neuro2a cells transiently expressing α1β3γ2 GABAAR were washed with the external buffer. Cells were then pre-incubated with 50 μM d-TC, BBIQA1, or BBIQA2 for 15 min at room temperature followed by 20 min incubation with 50 nM Alexa Fluor 546 αCtx conjugate in a final volume of 100 μL. Cells were then washed 3 times with two-fold excess buffer. Non-specific fluorescence was determined using 3 μM αCtx. Pictures of 3 independent fields chosen in blind mode in each plate well were taken using an epifluorescent microscope IX71 (Olympus, Japan) equipped with a digital CCD camera. Fluorescence intensity was counted using CellX and ImageJ open-source software. Fluorescence intensity was normalized to the mean integral intensity of the field incubated in presence of 50 nM Alexa Fluor 546 αCtx conjugate. Each experimental point is an average of integral intensity independently measured for 12 independent fields.

Electrophysiology

Plasmid pT7TS constructs of human nAChR α3, α9, α10, and β2 subunits were linearized with Xba I (NEB, USA), and the plasmid pGEMHE construct with the mouse 5-HT3A subunit–with Nhe I (NEB, USA). Linearized plasmid constructs were subjected to in vitro cRNA transcription using the T7 mMessage mMachine transcription kit (AMBION, USA). Stage V-VI Xlaevisoocytes were defolliculated with 2 mg/mL collagenase Type I (Life Technologies, USA) at room temperature (21–24°C) for 2 h in Ca2+-free Barth’s solution composed of (in mM) 88 NaCl, 1.1 KCl, 2.4 NaHCO3, 0.8 MgSO4, and 15 HEPES-NaOH at pH 7.6. Oocytes were injected with 9.2 ng of cRNAs of human nAChR α3β2, α9α10 (in a ratio 1:1), or mouse 5-HT3AR, or 2–3 ng of cDNAs of rat α4β2 (in a ratio 1:1, pcDNA3.1 vector), or mouse α3, β2, and γ2 GABAAR subunits in pCI vector. Oocytes were incubated at 18ºC in regular Barth’s solution composed of (in mM) 88 NaCl, 1.1 KCl, 2.4 NaHCO3, 0.3 Ca(NO3)2, 0.4 CaCl2, 0.8 MgSO4, and 15 HEPES-NaOH at pH 7.6, supplemented with 40 μg/mL gentamicin and 100 μg/mL ampicillin for 4 days before electrophysiological recordings. Two-electrode voltage clamp recordings were made using a turbo TEC-03X amplifier (Npi electronic, Germany) and WinWCP recording software (University of Strathclyde, UK), at a holding potential of -60 mV. Oocytes were briefly washed with Barth’s solution (for human nAChR α3β2, rat nAChR α4β2, mouse 5-HT3AR, and mouse α3β2γ2 GABAAR) or Ba2+ Ringer’s solution composed of (in mM) 115 NaCl, 2.5 KCl, 1.8 BaCl2, 10 HEPES at pH 7.2 (for human α9α10 nAChR) followed by 3 applications of agonist. Washout with Barth’s or Ba2+ Ringer’s solution was done for 5 min between agonist applications. Oocytes were pre-incubated with d-TC, BBIQA1, or BBIQA2 for 5 min followed by co-application with agonist (3 s). Peak current amplitudes of agonist-evoked responses were measured before and after pre-incubation of oocytes with d-TC, BBIQA1, or BBIQA2. The ratio between these two measurements was used to assess the activity of d-TC, BBIQA1, or BBIQA2 on human nAChR α3β2, α9α10, rat α4β2, mouse 5-HT3AR, and mouse α3β2γ2 GABAAR. The results for 5-HT3AR were analyzed using OriginPro 2017 (OriginLab Corporation, Northampton, MA, USA) fitting to a dose-response curve with a variable Hill slope using the equation y = A1 + (A2-A1)/(1 + 10(LOGx0-x)*p): A1 and A2 are bottom and top asymptotes, respectively; p–Hill slope; LOGx0 –log10[IC50, M].

Data and statistical analysis

Data are presented as mean with 95% confidence interval (CI) or mean ± SEM for the indicated number of independent experiments (n). Statistical analysis (One-way ANOVA with Tukey’s HSD test) was performed using OriginPro 2017 software (OriginLab Corporation, Northampton, MA, USA). In all the tests, p < 0.05 was taken as significant.

Results

Computer modeling with Acalifornica AChBP

Molecular docking studies were performed to detect possible activity against nAChRs. It was not unexpected for BBIQAs to bind to nAChR because they have structures similar to that of d-TC. However, such intuitive suggestion should be tested with Acalifornica AChBP, co-crystallized with the d-TC, an excellent nAChR model for this work [33]. Docking results showed the potential for high affinity binding of BBIQA1 and BBIQA2 at the Acalifornica AChBP binding site (S1 Table). BBIQA1 docking poses differed substantially from the experimentally observed d-TC orientations, possibly due to greater structure flexibility (S1A Fig), while predicted BBIQA2 binding poses (S1B Fig) were similar. All docking simulations predicted BBIQA1 (S1C Fig) and BBIQA2 (S1D Fig) to bind at the classic orthosteric site under loop C of AChBP.

Samples purification

All experiments were carried out on BBIQAs1-3 after additional HPLC purification and using individual peaks according to analytical HPLC (S2S7 Figs). The mass spectra given for each of these compounds provide further confirmation of their purity and structure.

Radioligand Assay

Radioligand analysis was performed using [125I]-α-bungarotoxin (αBgt) with a specific activity of 500 Ci/mmol and either Acalifornica AChBP, membrane-bound nAChRs from the electric organ of the Tcalifornica ray, or human α7 nAChR in GH4C1 cells as previously described [41].
BBIQA1 and BBIQA2 demonstrate a little lower ability to compete with [125I]-αBgt for binding to Acalifornica AChBP in comparison with d-TC, their IC50 values were 3.80 μM, 7.63 μM, and 2.41 μM, respectively (Fig 2A and Table 1). In case of Tcalifornica nAChR (Fig 2B and Table 1) the lowest affinity was registered for BBIQA1 (IC50 26.3 μM), with slightly higher values for BBIQA3 (17 μM), and BBIQA2 (8.75 μM). However, even BBIQA2 was 22-fold less potent than d-TC (IC50 0.39 μM). Inhibition of specific [125I]-αBgt binding to α7 nAChR (Fig 2C) revealed that the affinity of BBIQA1 is low (IC50 162 μM), but that of BBIQA2 is similar to d-TC (IC50 5.52 μM and 7.77 μM, respectively).
thumbnail
Fig 2. Inhibition of [125I]-labeled α-bungarotoxin binding.
(a) to Acalifornica AChBP, (b) to nAChR from Torpedo californica electric organ membranes, and (c) to human α7 nAChR in GH4C1 cells with BBIQA1 (open blue triangles), BBIQA2 (open red squares), BBIQA3 (open green rhombuses), and d-TC (open black circles). Data are mean ± SEM of two biological replicates with duplicates for each point (the number of technical replicates is 2), i.e. n = 4. IC50 values derived from these data are shown in Table 1.
thumbnail
Table 1. Inhibitory effects of the alkaloids from Matis Dart Poison and d-TC on [125I]-αBgt binding, agonist-induced calcium responses, and/or currents mediated by Acalifornica AChBP, Tcalifornica nAChR, mouse adult α1β1εδ nAChR, human α7 nAChR, and mouse 5-HT3AR.

Calcium imaging

The next step was the analysis of BBIQAs interactions with mouse adult muscle nAChR (α1β1εδ) (Fig 3A and 3C) which differs by one subunit (ε) from the Tcalifornica nAChR (α1β1γδ). This receptor was heterologously expressed in Neuro2a cells and the inhibitory effects of BBIQAs on the ACh-induced increase in intracellular Сa2+ concentration ([Ca2+]i) were monitored based on the fluorescence of the co-expressed calcium ion sensor Case12 as described in [39]. As seen from Fig 3A and Table 1, BBIQA1 and d-TC have virtually identical IC50 values (0.75 and 0.81 μM, respectively) and BBIQA2 is ~2-fold less potent (1.75 μM). The same relationship was observed for the human α7 nAChR expressed in Neuro2a cells and analyzed in the presence of a positive allosteric modulator PNU120596 (Fig 3B). Due to extremely rapid desensitization rate of α7 nAChR, its kinetics is too fast to be measured using any conventional calcium imaging techniques [4245]. Instead, the application of its positive allosteric modulator PNU120596 significantly decreases α7 nAChR desensitization rate providing the possibility to monitor as the receptor-mediated intracellular calcium raises [4647]. PNU120596 binds to the transmembrane domain of the α7 receptor leading to an increase of the maximal agonist-elicited response, agonist potency and apparent cooperativity [4850]. In general, positive allosteric modulators have little to no effect on equilibrium binding of α7 competitive antagonists, because these modulators bind in a non-competitive manner, away from the traditional agonist binding site [51]. However, we must be cautious with calcium imaging data interpretation, as PNU120596 could affect the data on inhibition of α7 nAChR by BBIQAs due to their complex character of interaction.
thumbnail
Fig 3. Dose-response curves of inhibitory activity of d-TC (open black circles), BBIQA1 (open blue triangles), and BBIQA2 (open red squares).
(a) on the 30 μM (approx. EC50) acetylcholine-evoked intracellular calcium ion concentration ([Ca2+]i) rise in Neuro2a cells expressing mouse adult α1β1εδ nAChRs; and (b) on the 10 μM (approx. EC50) acetylcholine-evoked [Ca2+]i rise in Neuro2a cells expressing human α7 nAChRs in the presence of 10 μM PNU120596. Data are presented as mean ± SEM, n = 3. The respective IC50 values are shown in the Table 1. (c) Dose-response curves of acetylcholine (ACh)-evoked [Ca2+]i rise in the absence (grey circles, EC50 = 8.4 ± 0.6 μM) and presence of d-TC (black circles) or its analogs, BBIQA1 (blue triangles), and BBIQA2 (red squares) at different concentrations in Neuro2a cells expressing mouse adult α1β1εδ nAChRs. Data are presented as mean ± SEM, n = 3. EC50 values are shown in the Table 2.
The IC50 values for d-TC and BBIQA1 were 2.03 and 1.70 μM, respectively, while for BBIQA2 this value was 3.08 μM. Data on the inhibitory activity of BBIQAs from the calcium imaging experiments and from the above-described data on the muscle-type and α7 nAChRs obtained in radioligand experiments are shown in Table 1.
To determine possible differences between d-TC and its BBIQA1 / BBIQA2 analogs in binding to orthosteric or allosteric sites, we performed a series of calcium imaging experiments with Neuro2a cells expressing mouse adult α1β1εδ nAChR (Fig 3C). d-TC or its analogs at different concentrations (at IC50 and 3x IC50, see Table 1: Calcium Imaging for α1β1εδ nAChR) were co-applied along with acetylcholine to Neuro2a cells expressing α1β1εδ nAChR. As shown in Fig 3C, the presence of d-TC, BBIQA1, or BBIQA2 shifts the acetylcholine dose response curve to the right, with an increased acetylcholine EC50 value (Table 2), and reduces the maximal acetylcholine response (Fig 3C).
thumbnail
Table 2. Inhibitory effects of the alkaloids from Matis Dart Poison and d-TC at different concentrations (IC50 and 3x IC50) on the acetylcholine dose–response curve at mouse adult α1β1εδ nAChR expressed in Neuro2a cells.

Electrophysiology

We also analyzed the interactions of d-TC, BBIQAs1 and 2 with several subtypes of heteromeric neuronal nAChRs, namely α3β2, α4β2, and α9α10, following their expression in Xlaevis oocytes. As is seen in Fig 4A and 4B, both BBIQAs and d-TC have similar inhibitory activities against all three neuronal nAChR subtypes. We also wanted to check whether d-TC and BBIQAs differ in their activity towards other Cys-loop receptors. Fig 4C shows that a strong (90%) inhibition of the mouse α3β2γ2 GABAAR is observed with a high concentration of d-TC (100 μM), while BBIQAs1 and 2 at this concentration produce only a 40–50% inhibition.
thumbnail
Fig 4. Activity of d-TC, BBIQAs1 and 2 against heteromeric neuronal nAChRs, GABAAR, and 5-HT3AR.
(a) Representative current traces of human α3β2 nAChR, showing inhibition of nicotine (50 μM)-induced current by 10 μM d-TC, BBIQA1, or BBIQA2. (b) Bar graph for d-TC and BBIQAs (1 and 10 μM) inhibition of agonist-evoked currents mediated by human α3β2 (50 μM Nicotine), rat α4β2 (10 μM Nicotine), and human α9α10 (25 μM Acetylcholine) nAChRs. (c) Bar graph for d-TC and BBIQAs (100 μM) inhibition of agonist-evoked currents mediated by mouse α3β2γ2 GABAAR (100 μM GABA). (d) Inhibition of Alexa Fluor 546 α-cobratoxin (αCtx, 50 nM) binding to α1β3γ2 GABAAR expressed in Neuro2a cells by 50 μM d-TC, BBIQAs 1 and 2. The bar graph represents the remaining fluorescence of Alexa Fluor 546 αCtx (50 nM). In (a, b, c, d) data are presented as mean ± SEM, n = 3–6. One-way ANOVA with Tukey’s HSD test, (black asterisksp<0.05, normalized current evoked by agonist in the presence of d-TC, BBIQA1, or BBIQA2 vs normalized current induced by agonist in the absence of antagonists). (e) Dose-response curves of d-TC (open black circles), BBIQA1 (open blue triangles), or BBIQA2 (open red squares) inhibitory action on 1 μM 5-HT-evoked ion currents mediated by mouse 5-HT3AR. Data are presented as mean ± SEM, n = 3–5. IC50 values determined from these data are shown in Table 1.
As d-TC is known to potently inhibit 5-HTR [3435], we checked whether such a property is also inherent in BBIQAs (Fig 4E). A considerable difference was detected between BBIQAs1 and 2: the IC50 value for the former was 119.4 nM, while for the latter was 1975 nM, and these two compounds were 5- and 87-fold less active than d-TC, the IC50 value of which is 22.63 nM (IC50 values and 95% confidence intervals are shown in Table 1).

Fluorescence assay

For the GABAAR, we compared functional effects of the compounds on GABA-induced currents with their binding ability: Fig 4D shows that d-TC and BBIQAs 1 and 2 have similar potency in competition assays with the fluorescent derivative of αCtx, which may bind to both orthosteric and allosteric sites [40].

Discussion

In this study we compared the effects of d-TC and the structurally-related BBIQAs on the muscle-type and neuronal nAChRs, as well as on 5-HT3AR and GABAAR. The data obtained suggest that the BBIQAs do not completely mimic the effects of d-TC, and thus could provide useful tools for future studies.
Computer modeling with Acalifornica AChBP revealed that the BBIQAs may bind to Cys-loop receptors with high affinity at the orthosteric sites, as d-TC does. Indeed, radioligand binding assay with Acalifornica AChBP showed that both BBIQA1 and BBIQA2 had the affinity almost indistinguishable from that of d-TC. However, radioligand assays with Tcalifornica nAChR revealed that BBIQAs are significantly less potent than d-TC in their ability to inhibit binding of [125I]-αBgt. It is well-established that αBgt interacts with the binding sites for agonists and competitive antagonists (in other words, binds to the orthosteric sites [5253]), so our results suggest that the BBIQAs are less potent than d-TC in binding to that site. In calcium imaging experiments on the adult muscle nAChR, the IC50 values for BBIQA1 and BBIQA2 were roughly equal to those of d-TC, in being in the micromolar range and indicating that d-TC, BBIQAs 1 and 2 cannot be considered as highly potent drugs. However, the analysis of the relationship between the neuromuscular blocking dose and duration of action revealed that it is reciprocal, and the rapid onset and short duration of action are associated with drugs of low potency [265455]. Therefore, low potency of BBIQAs, as well as of d-TC, may be a positive characteristic in terms of their drug-like properties, justifying further structure-functional analysis of these compounds.
Interestingly, our calcium imaging data for muscle α1β1εδ nAChR revealed the BBIQAs potencies similar to that of d-TC, in contrast to radioligand assay data where all three BBIQAs were much weaker than d-TC, indicating that inhibition of a functional response was more complex than a simple competition. We assumed that these data could be explained by the BBIQAs binding to multiple sites on the receptor: nAChRs are large multisubunit proteins with various binding sites, orthosteric and allosteric, both in the LBD and in the transmembrane domain [5658]. Indeed, it is well known that d-TC is able to target both the orthosteric and allosteric sites [245962].
Calcium imaging is a convenient method for probing nAChR inhibition [394763]. In [47] the researchers used it with genetically encoded calcium ion sensor to study the mechanism of nicotinic receptors (α7 and α4β2 nAChRs) inhibition by well-known competitive antagonists (methyllycaconitine and dihydro-β-erythroidine, respectively). In the presence of competitive antagonist, the agonist dose-response curve shifts to the right with no reduction in the amplitude of the maximal response. Noncompetitive inhibition could as well be detected by calcium imaging with genetically encoded calcium ion sensor Case12. It was previously detected that baptide 2, a short peptide from puff adder Bitis arietans venom, inhibits α1β1εδ nAChR in noncompetitive way: the pre-incubation with baptide 2 doesn’t shift the acetylcholine dose-response curve to the right, but reduces the maximal amplitude of the responses [63]. Here we also checked the applicability of the calcium imaging with Case12 on the example of α1β1εδ nAChR inhibition by a well-known competitive antagonist, αBgt. A 5 min pre-incubation of Neuro2a cells with αBgt (100 nM) shifts the acetylcholine dose-response curve to the right, and the agonist EC50 value increases from 11.5 μM to 44.2 μM (S8 Fig). Moreover, we did not observe any reduction in the maximal amplitude of the cellular responses (S8 Fig). This indicates clearly that under conditions described in this article calcium imaging is applicable for the detection of competitive antagonism.
We compared BBIQA1, BBIQA2, and d-TC by this calcium imaging test and detected quite unusual effect of all three compounds on the acetylcholine dose-response curve for the α1β1εδ receptor: in addition to the curve shift to the right (with increasing agonist EC50 value, see Table 2), for all three compounds there was a concomitant reduction of peak acetylcholine responses (Fig 3C). This finding, along with the results of radioligand assay, supports our hypothesis that BBIQA1 and BBIQA2 inhibit adult muscle α1β1εδ nAChR by mixed mechanisms, in which both competitive and noncompetitive modes are involved. For d-TC this mixed mechanism was confirmed by the earlier data [60], where the same effect on the muscle-type nAChR was observed in the electrophysiology experiments.
In GH4C1 cells expressing α7 nAChR, d-TC competitively inhibits binding as previously reported [31]. BBIQA2 is also able to compete with [125I]-αBgt with a potency similar to d-TC, but BBIQA1 is considerably less potent. On the contrary, in calcium imaging experiments, we did not see any strong differences in the inhibitory potencies of all three compounds. Thus BBIQA1 might also interact with allosteric sites in α7 nAChR. A similar example has been reported for the muscle-type receptor [60], where the use of nondesensitizing agonist (e.g. DMPP) resulted in competitive inhibition by d-TC, instead of noncompetitive inhibition in the case of acetylcholine-evoked responses. Thus, the mechanism is dependent on the receptor activation mode [60]. However, we must be cautious, as the calcium imaging experiments were performed in the presence of positive allosteric modulator PNU120596 (see above), which decreases the receptor desensitization.
The data for α3β2, α4β2, and α9α10 nAChRs show that against them BBIQAs and d-TC have very similar inhibitory potencies. A comparable inhibitory activity of d-TC and BBIQAs against α9α10 nAChR is of considerable importance because these receptors are targets for the development of new analgesics, among which several α-conotoxins are being explored [64]. Both BBIQA1 and BBIQA2 demonstrate some decrease in their inhibiting potencies at GABAAR in comparison with d-TC in electrophysiological experiments, but there were no significant differences between them in competition with Alexa Fluor 546 αCtx, which can bind both to orthosteric and allosteric sites [40]. Thus, the situation with GABAAR might also indicate a more complicated inhibitory mechanism, in which binding to orthosteric and allosteric sites could be involved.
At 5-HT3AR, however, BBIQA1 and BBIQA2 behave quite differently, being 5- and 87-fold less active than d-TC, thus demonstrating that compared to d-TC BBIQAs are less selective for 5-HT3AR than nAChRs. This finding may be important since some α7 nAChR targeting drugs are problematic because they also act on 5-HT3R. Further investigating of key structural elements responsible for mechanisms of action of these compounds might be very useful for future drug design.

Supporting information

pone.0210182.s001.docx
Predicted affinity of the alkaloids from Matis Dart Poison and d-TC for AChBP from Acalifornica.
(DOCX)

S1 Table. Predicted affinity of the alkaloids from Matis Dart Poison and d-TC for AChBP from Acalifornica.

(DOCX)

S1 Fig. Computer docking of the BBIQA1 and BBIQA2 to Acalifornica AChBP.

Docked structures overlay of (a) BBIQA1 (purple) and d-TC (black) revealing differences in the predicted binding modes and (b) BBIQA2 (purple) and d-TC (black) revealing similarity in the predicted binding modes. Despite the difference in binding configurations, both BBIQA1 (c) and BBIQA2 (d) were docked at the classic orthosteric site under the loop C of the AChBP.
(TIF)

S2 Fig. Analytical HPLC data for BBIQA1 purification.

(TIF)

S4 Fig. Analytical HPLC data for BBIQA2 purification.

(TIF)

S6 Fig. Analytical HPLC data for BBIQA3 purification.

(TIF)

S8 Fig. Dose-response curves of acetylcholine (ACh)-evoked [Ca2+]i rise in the absence (grey circles, EC50 = 11.5 ± 1.43 μM) and presence of 100 nM αBgt (violet circles, EC50 = 44.2 ± 8.38 μM) in Neuro2a cells expressing mouse adult α1β1εδ nAChRs.

Data are presented as mean ± SEM, n = 3.
(TIF)

References

  1. 1.Bertrand D, Terry AV Jr. The wonderland of neuronal nicotinic acetylcholine receptors. Biochem Pharmacol. 2018;151:214–25. pmid:29248596
  2. 2.Burgos CF, Yevenes GE, Aguayo LG. Structure and Pharmacologic Modulation of Inhibitory Glycine Receptors. Molecular pharmacology. 2016;90(3):318–25. pmid:27401877; PubMed Central PMCID: PMC4998662.
  3. 3.Chua HC, Chebib M. GABAA Receptors and the Diversity in their Structure and Pharmacology. Adv Pharmacol. 2017;79:1–34. pmid:28528665.
  4. 4.Dani JA. Neuronal Nicotinic Acetylcholine Receptor Structure and Function and Response to Nicotine. International review of neurobiology. 2015;124:3–19. pmid:26472524; PubMed Central PMCID: PMC4795468.
  5. 5.Wu ZS, Cheng H, Jiang Y, Melcher K, Xu HE. Ion channels gated by acetylcholine and serotonin: structures, biology, and drug discovery. Acta pharmacologica Sinica. 2015;36(8):895–907. pmid:26238288; PubMed Central PMCID: PMC4564887.
  6. 6.Chua HC, Chebib M. GABAA Receptors and the Diversity in their Structure and Pharmacology. Adv Pharmacol. 2017;79:1–34. Epub 2017/05/23. pmid:28528665.
  7. 7.Hassaine G, Deluz C, Grasso L, Wyss R, Tol MB, Hovius R, et al. X-ray structure of the mouse serotonin 5-HT3 receptor. Nature. 2014;512(7514):276–81. Epub 2014/08/15. pmid:25119048.
  8. 8.Morales-Perez CL, Noviello CM, Hibbs RE. X-ray structure of the human alpha4beta2 nicotinic receptor. Nature. 2016;538(7625):411–5. Epub 2016/10/21. pmid:27698419; PubMed Central PMCID: PMCPMC5161573.
  9. 9.Nys M, Kesters D, Ulens C. Structural insights into Cys-loop receptor function and ligand recognition. Biochem Pharmacol. 2013;86(8):1042–53. Epub 2013/07/16. pmid:23850718.
  10. 10.Spurny R, Debaveye S, Farinha A, Veys K, Vos AM, Gossas T, et al. Molecular blueprint of allosteric binding sites in a homologue of the agonist-binding domain of the alpha7 nicotinic acetylcholine receptor. Proc Natl Acad Sci U S A. 2015;112(19):E2543–52. Epub 2015/04/29. pmid:25918415; PubMed Central PMCID: PMCPMC4434711.
  11. 11.Spurny R, Ramerstorfer J, Price K, Brams M, Ernst M, Nury H, et al. Pentameric ligand-gated ion channel ELIC is activated by GABA and modulated by benzodiazepines. Proc Natl Acad Sci U S A. 2012;109(44):E3028–34. Epub 2012/10/05. pmid:23035248; PubMed Central PMCID: PMCPMC3497736.
  12. 12.Braun M, Ramracheya R, Bengtsson M, Clark A, Walker JN, Johnson PR, et al. Gamma-aminobutyric acid (GABA) is an autocrine excitatory transmitter in human pancreatic beta-cells. Diabetes. 2010;59(7):1694–701. Epub 2010/04/24. pmid:20413510; PubMed Central PMCID: PMCPMC2889769.
  13. 13.Gahring LC, Myers EJ, Dunn DM, Weiss RB, Rogers SW. Nicotinic alpha 7 receptor expression and modulation of the lung epithelial response to lipopolysaccharide. PLoS One. 2017;12(4):e0175367. Epub 2017/04/07. pmid:28384302; PubMed Central PMCID: PMCPMC5383308.
  14. 14.Kawashima K, Fujii T, Moriwaki Y, Misawa H. Critical roles of acetylcholine and the muscarinic and nicotinic acetylcholine receptors in the regulation of immune function. Life Sci. 2012;91(21–22):1027–32. Epub 2012/06/05. pmid:22659391.
  15. 15.Matsuda M, Imaoka T, Vomachka AJ, Gudelsky GA, Hou Z, Mistry M, et al. Serotonin regulates mammary gland development via an autocrine-paracrine loop. Dev Cell. 2004;6(2):193–203. Epub 2004/02/13. pmid:14960274.
  16. 16.Zoli M, Pucci S, Vilella A, Gotti C. Neuronal and extraneuronal nicotinic acetylcholine receptors. Curr Neuropharmacol. 2017. Epub 2017/09/14. pmid:28901280.
  17. 17.Dineley KT, Pandya AA, Yakel JL. Nicotinic ACh receptors as therapeutic targets in CNS disorders. Trends Pharmacol Sci. 2015;36(2):96–108. Epub 2015/02/03. pmid:25639674; PubMed Central PMCID: PMCPMC4324614.
  18. 18.Foucault-Fruchard L, Antier D. Therapeutic potential of alpha7 nicotinic receptor agonists to regulate neuroinflammation in neurodegenerative diseases. Neural Regen Res. 2017;12(9):1418–21. Epub 2017/11/02. pmid:29089979; PubMed Central PMCID: PMCPMC5649454.
  19. 19.Lemoine D, Jiang R, Taly A, Chataigneau T, Specht A, Grutter T. Ligand-gated ion channels: new insights into neurological disorders and ligand recognition. Chem Rev. 2012;112(12):6285–318. Epub 2012/09/20. pmid:22988962.
  20. 20.Thompson AJ, Lummis SC. The 5-HT3 receptor as a therapeutic target. Expert Opin Ther Targets. 2007;11(4):527–40. Epub 2007/03/22. pmid:17373882; PubMed Central PMCID: PMCPMC1994432.
  21. 21.Tuba Z, Maho S, Vizi ES. Synthesis and structure-activity relationships of neuromuscular blocking agents. Current medicinal chemistry. 2002;9(16):1507–36. pmid:12171561.
  22. 22.Heier T. [Muscle relaxants]. Tidsskrift for den Norske laegeforening: tidsskrift for praktisk medicin, ny raekke. 2010;130(4):398–401. pmid:20220868.
  23. 23.Wintersteiner O, Dutcher JD. Curare Alkaloids from Chondodendron Tomentosum. Science. 1943;97(2525):467–70. pmid:17789877.
  24. 24.Sine SM. End-plate acetylcholine receptor: structure, mechanism, pharmacology, and disease. Physiol Rev. 2012;92(3):1189–234. Epub 2012/07/20. pmid:22811427; PubMed Central PMCID: PMCPMC3489064.
  25. 25.Bowman WC. Neuromuscular block. Br J Pharmacol. 2006;147 Suppl 1:S277–86. pmid:16402115; PubMed Central PMCID: PMC1760749.
  26. 26.Meakin GH. Recent advances in myorelaxant therapy. Paediatric anaesthesia. 2001;11(5):523–31. pmid:11696115.
  27. 27.Savarese JJ, Kitz RJ. Does clinical anesthesia need new neuromuscular blocking agents? Anesthesiology. 1975;42(3):236–9. pmid:123132.
  28. 28.Malca Garcia GR, Hennig L, Shelukhina IV, Kudryavtsev DS, Bussmann RW, Tsetlin VI, et al. Curare Alkaloids: Constituents of a Matis Dart Poison. Journal of natural products. 2015;78(11):2537–44. pmid:26496427.
  29. 29.Murebwayire S, Ingkaninan K, Changwijit K, Frederich M, Duez P. Triclisia sacleuxii (Pierre) Diels (Menispermaceae), a potential source of acetylcholinesterase inhibitors. The Journal of pharmacy and pharmacology. 2009;61(1):103–7. pmid:19126303.
  30. 30.Semwal DK, Semwal RB, Vermaak I, Viljoen A. From arrow poison to herbal medicine—the ethnobotanical, phytochemical and pharmacological significance of Cissampelos (Menispermaceae). Journal of ethnopharmacology. 2014;155(2):1011–28. pmid:24997389.
  31. 31.Jonsson M, Gurley D, Dabrowski M, Larsson O, Johnson EC, Eriksson LI. Distinct pharmacologic properties of neuromuscular blocking agents on human neuronal nicotinic acetylcholine receptors: a possible explanation for the train-of-four fade. Anesthesiology. 2006;105(3):521–33. pmid:16931985.
  32. 32.Brejc K, van Dijk WJ, Klaassen RV, Schuurmans M, van Der Oost J, Smit AB, et al. Crystal structure of an ACh-binding protein reveals the ligand-binding domain of nicotinic receptors. Nature. 2001;411(6835):269–76. pmid:11357122.
  33. 33.Brams M, Pandya A, Kuzmin D, van Elk R, Krijnen L, Yakel JL, et al. A structural and mutagenic blueprint for molecular recognition of strychnine and d-tubocurarine by different cys-loop receptors. PLoS biology. 2011;9(3):e1001034. pmid:21468359; PubMed Central PMCID: PMC3066128.
  34. 34.Paul M, Callahan R, Au J, Kindler CH, Yost CS. Antiemetics of the 5-hydroxytryptamine 3A antagonist class inhibit muscle nicotinic acetylcholine receptors. Anesthesia and analgesia. 2005;101(3):715–21, table of contents. pmid:16115980.
  35. 35.Zhang R, Wen X, Militante J, Hester B, Rhubottom HE, Sun H, et al. The role of loop F residues in determining differential d-tubocurarine potencies in mouse and human 5-hydroxytryptamine 3A receptors. Biochemistry. 2007;46(5):1194–204. pmid:17260949.
  36. 36.Wotring VE, Yoon KW. The inhibitory effects of nicotinic antagonists on currents elicited by GABA in rat hippocampal neurons. Neuroscience. 1995;67(2):293–300. pmid:7675170.
  37. 37.Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. Journal of computational chemistry. 2009;30(16):2785–91. pmid:19399780; PubMed Central PMCID: PMC2760638.
  38. 38.Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. Journal of computational chemistry. 2004;25(13):1605–12. pmid:15264254.
  39. 39.Shelukhina I, Spirova E, Kudryavtsev D, Ojomoko L, Werner M, Methfessel C, et al. Calcium imaging with genetically encoded sensor Case12: Facile analysis of alpha7/alpha9 nAChR mutants. PLoS One. 2017;12(8):e0181936. pmid:28797116; PubMed Central PMCID: PMC5552293.
  40. 40.Kudryavtsev DS, Shelukhina IV, Son LV, Ojomoko LO, Kryukova EV, Lyukmanova EN, et al. Neurotoxins from snake venoms and alpha-conotoxin ImI inhibit functionally active ionotropic gamma-aminobutyric acid (GABA) receptors. The Journal of biological chemistry. 2015;290(37):22747–58. pmid:26221036; PubMed Central PMCID: PMC4566246.
  41. 41.Kudryavtsev D, Makarieva T, Utkina N, Santalova E, Kryukova E, Methfessel C, et al. Marine natural products acting on the acetylcholine-binding protein and nicotinic receptors: from computer modeling to binding studies and electrophysiology. Marine drugs. 2014;12(4):1859–75. pmid:24686559; PubMed Central PMCID: PMC4012456.
  42. 42.Bouzat C, Bartos M, Corradi J, Sine SM. The interface between extracellular and transmembrane domains of homomeric Cys-loop receptors governs open-channel lifetime and rate of desensitization. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2008;28(31):7808–19. pmid:18667613; PubMed Central PMCID: PMC3844810.
  43. 43.Dickinson JA, Hanrott KE, Mok MH, Kew JN, Wonnacott S. Differential coupling of alpha7 and non-alpha7 nicotinic acetylcholine receptors to calcium-induced calcium release and voltage-operated calcium channels in PC12 cells. Journal of neurochemistry. 2007;100(4):1089–96. pmid:17181555.
  44. 44.Inserra MC, Kompella SN, Vetter I, Brust A, Daly NL, Cuny H, et al. Isolation and characterization of alpha-conotoxin LsIA with potent activity at nicotinic acetylcholine receptors. Biochem Pharmacol. 2013;86(6):791–9. pmid:23924607.
  45. 45.Roncarati R, Seredenina T, Jow B, Jow F, Papini S, Kramer A, et al. Functional properties of alpha7 nicotinic acetylcholine receptors co-expressed with RIC-3 in a stable recombinant CHO-K1 cell line. Assay and drug development technologies. 2008;6(2):181–93. pmid:18471073.
  46. 46.Bertrand D, Gopalakrishnan M. Allosteric modulation of nicotinic acetylcholine receptors. Biochem Pharmacol. 2007;74(8):1155–63. pmid:17707779.
  47. 47.Yamauchi JG, Nemecz A, Nguyen QT, Muller A, Schroeder LF, Talley TT, et al. Characterizing ligand-gated ion channel receptors with genetically encoded Ca2++ sensors. PLoS One. 2011;6(1):e16519. pmid:21305050; PubMed Central PMCID: PMC3030600.
  48. 48.Hurst RS, Hajos M, Raggenbass M, Wall TM, Higdon NR, Lawson JA, et al. A novel positive allosteric modulator of the alpha7 neuronal nicotinic acetylcholine receptor: in vitro and in vivo characterization. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2005;25(17):4396–405. pmid:15858066.
  49. 49.Szabo AK, Pesti K, Mike A, Vizi ES. Mode of action of the positive modulator PNU-120596 on alpha7 nicotinic acetylcholine receptors. Neuropharmacology. 2014;81:42–54. pmid:24486377.
  50. 50.Young GT, Zwart R, Walker AS, Sher E, Millar NS. Potentiation of alpha7 nicotinic acetylcholine receptors via an allosteric transmembrane site. Proc Natl Acad Sci U S A. 2008;105(38):14686–91. pmid:18791069; PubMed Central PMCID: PMC2535569.
  51. 51.Williams DK, Wang J, Papke RL. Positive allosteric modulators as an approach to nicotinic acetylcholine receptor-targeted therapeutics: advantages and limitations. Biochem Pharmacol. 2011;82(8):915–30. pmid:21575610; PubMed Central PMCID: PMC3162128.
  52. 52.daCosta CJ, Free CR, Sine SM. Stoichiometry for alpha-bungarotoxin block of alpha7 acetylcholine receptors. Nature communications. 2015;6:8057. pmid:26282895; PubMed Central PMCID: PMC4544739.
  53. 53.Huang S, Li SX, Bren N, Cheng K, Gomoto R, Chen L, et al. Complex between alpha-bungarotoxin and an alpha7 nicotinic receptor ligand-binding domain chimaera. The Biochemical journal. 2013;454(2):303–10. pmid:23800261; PubMed Central PMCID: PMC3920732.
  54. 54.Bowman WC, Rodger IW, Houston J, Marshall RJ, McIndewar I. Structure:action relationships among some desacetoxy analogues of pancuronium and vecuronium in the anesthetized cat. Anesthesiology. 1988;69(1):57–62. pmid:2898902.
  55. 55.Lien CA. Development and potential clinical impairment of ultra-short-acting neuromuscular blocking agents. British journal of anaesthesia. 2011;107 Suppl 1:i60–71. pmid:22156271.
  56. 56.Chatzidaki A, D'Oyley JM, Gill-Thind JK, Sheppard TD, Millar NS. The influence of allosteric modulators and transmembrane mutations on desensitisation and activation of alpha7 nicotinic acetylcholine receptors. Neuropharmacology. 2015;97:75–85. pmid:25998276; PubMed Central PMCID: PMC4548482.
  57. 57.Gill JK, Chatzidaki A, Ursu D, Sher E, Millar NS. Contrasting properties of alpha7-selective orthosteric and allosteric agonists examined on native nicotinic acetylcholine receptors. PLoS One. 2013;8(1):e55047. pmid:23383051; PubMed Central PMCID: PMC3558472.
  58. 58.Vemparala S, Domene C, Klein ML. Computational studies on the interactions of inhalational anesthetics with proteins. Accounts of chemical research. 2010;43(1):103–10. pmid:19788306.
  59. 59.Groot-Kormelink PJ, Ferrand S, Kelley N, Bill A, Freuler F, Imbert PE, et al. High Throughput Random Mutagenesis and Single Molecule Real Time Sequencing of the Muscle Nicotinic Acetylcholine Receptor. PLoS One. 2016;11(9):e0163129. Epub 2016/09/21. pmid:27649498; PubMed Central PMCID: PMCPMC5029940 BioMedical Research.
  60. 60.Jonsson Fagerlund M, Dabrowski M, Eriksson LI. Pharmacological characteristics of the inhibition of nondepolarizing neuromuscular blocking agents at human adult muscle nicotinic acetylcholine receptor. Anesthesiology. 2009;110(6):1244–52. Epub 2009/05/07. pmid:19417616.
  61. 61.Karpen JW, Hess GP. Acetylcholine receptor inhibition by d-tubocurarine involves both a competitive and a noncompetitive binding site as determined by stopped-flow measurements of receptor-controlled ion flux in membrane vesicles. Biochemistry. 1986;25(7):1786–92. Epub 1986/04/08. pmid:2423117.
  62. 62.Wang HL, Gao F, Bren N, Sine SM. Curariform antagonists bind in different orientations to the nicotinic receptor ligand binding domain. The Journal of biological chemistry. 2003;278(34):32284–91. Epub 2003/06/12. pmid:12799358.
  63. 63.Vulfius CA, Spirova EN, Serebryakova MV, Shelukhina IV, Kudryavtsev DS, Kryukova EV, et al. Peptides from puff adder Bitis arietans venom, novel inhibitors of nicotinic acetylcholine receptors. Toxicon: official journal of the International Society on Toxinology. 2016;121:70–6. pmid:27576061.
  64. 64.Hone AJ, Servent D, McIntosh JM. alpha9-containing nicotinic acetylcholine receptors and the modulation of pain. Br J Pharmacol. 2017. pmid:28662295.