Development of CHARMM compatible force field parameters and molecular dynamics simulations for the pesticide flupyradifurone.

Flupyradifurone (FLU) is a novel butenolide insecticide with partial agonist activity for insect nicotinic acetylcholine receptors. Its safety for non-target organisms has been questioned in the literature, despite initial claims of its harmlessness. Detailed understanding of its toxicity and related molecular mechanisms remain under discussion. Thus, in this work, an optimized set of CHARMM compatible parameters for FLU is presented. CHARMM General Force Field program was used as a starting point while the non-bonded and bonded parameters were adjusted and optimized to reproduce MP2/6-31G(d) accuracy level results. For the validity assessment of these parameters, infrared spectrum, water-octanol partition coefficient, and normal modes were computed and compared to experimental values found in the literature. Several MD simulations of FLU in water and FLU in complex with an acetylcholine-binding protein were performed to estimate the ability of the optimized parameters to correctly describe its torsional space and reproduce observed crystallographic trends respectively.

Efficient Expression in Leishmania tarentolae (LEXSY) of the Receptor-Binding Domain of the SARS-CoV-2 S-Protein and the Acetylcholine-Binding Protein from Lymnaea stagnalis.

Leishmania tarentolae (LEXSY) system is an inexpensive and effective expression approach for various research and medical purposes. The stated advantages of this system are the possibility of obtaining the soluble product in the cytoplasm, a high probability of correct protein folding with a full range of post-translational modifications (including uniform glycosylation), and the possibility of expressing multi-subunit proteins. In this paper, a LEXSY expression system has been employed for obtaining the receptor binding domain (RBD) of the spike-protein of the SARS-CoV-2 virus and the homopentameric acetylcholine-binding protein (AChBP) from Lymnaea stagnalis. RBD is actively used to obtain antibodies against the virus and in various scientific studies on the molecular mechanisms of the interaction of the virus with host cell targets. AChBP represents an excellent structural model of the ligand-binding extracellular domain of all subtypes of nicotinic acetylcholine receptors (nAChRs). Both products were obtained in a soluble glycosylated form, and their structural and functional characteristics were compared with those previously described.

Acetylcholine-Binding Protein Affinity Profiling of Neurotoxins in Snake Venoms with Parallel Toxin Identification.

Snakebite is considered a concerning issue and a neglected tropical disease. Three-finger toxins (3FTxs) in snake venoms primarily cause neurotoxic effects since they have high affinity for nicotinic acetylcholine receptors (nAChRs). Their small molecular size makes 3FTxs weakly immunogenic and therefore not appropriately targeted by current antivenoms. This study aims at presenting and applying an analytical method for investigating the therapeutic potential of the acetylcholine-binding protein (AChBP), an efficient nAChR mimic that can capture 3FTxs, for alternative treatment of elapid snakebites. In this analytical methodology, snake venom toxins were separated and characterised using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS) and high-throughput venomics. By subsequent nanofractionation analytics, binding profiling of toxins to the AChBP was achieved with a post-column plate reader-based fluorescence-enhancement ligand displacement bioassay. The integrated method was established and applied to profiling venoms of six elapid snakes (Naja mossambica, Ophiophagus hannah, Dendroaspis polylepis, Naja kaouthia, Naja haje and Bungarus multicinctus). The methodology demonstrated that the AChBP is able to effectively bind long-chain 3FTxs with relatively high affinity, but has low or no binding affinity towards short-chain 3FTxs, and as such provides an efficient analytical platform to investigate binding affinity of 3FTxs to the AChBP and mutants thereof and to rapidly identify bound toxins.

Residues Responsible for the Selectivity of α-Conotoxins for Ac-AChBP or nAChRs.

Nicotinic acetylcholine receptors (nAChRs) are targets for developing new drugs to treat severe pain, nicotine addiction, Alzheimer disease, epilepsy, etc. α-Conotoxins are biologically and chemically diverse. With 12-19 residues and two disulfides, they can be specifically selected for different nAChRs. Acetylcholine-binding proteins from Aplysia californica (Ac-AChBP) are homologous to the ligand-binding domains of nAChRs and pharmacologically similar. X-ray structures of the α-conotoxin in complex with Ac-AChBP in addition to computer modeling have helped to determine the binding site of the important residues of α-conotoxin and its affinity for nAChR subtypes. Here, we present the various α-conotoxin residues that are selective for Ac-AChBP or nAChRs by comparing the structures of α-conotoxins in complex with Ac-AChBP and by modeling α-conotoxins in complex with nAChRs. The knowledge of these binding sites will assist in the discovery and design of more potent and selective α-conotoxins as drug leads.

Molecular recognition of thiaclopride by Aplysia californica AChBP: new insights from a computational investigation.

The binding of thiaclopride (THI), a neonicotinoid insecticide, with Aplysia californica acetylcholine binding protein (Ac-AChBP), the surrogate of the extracellular domain of insects nicotinic acetylcholine receptors, has been studied with a QM/QM' hybrid methodology using the ONIOM approach (M06-2X/6-311G(d):PM6). The contributions of Ac-AChBP key residues for THI binding are accurately quantified from a structural and energetic point of view. The importance of water mediated hydrogen-bond (H-bond) interactions involving two water molecules and Tyr55 and Ser189 residues in the vicinity of the THI nitrile group, is specially highlighted. A larger stabilization energy is obtained with the THI-Ac-AChBP complex compared to imidacloprid (IMI), the forerunner of neonicotinoid insecticides. Pairwise interaction energy calculations rationalize this result with, in particular, a significantly more important contribution of the pivotal aromatic residues Trp147 and Tyr188 with THI through CH···π/CH···O and π-π stacking interactions, respectively. These trends are confirmed through a complementary non-covalent interaction (NCI) analysis of selected THI-Ac-AChBP amino acid pairs.

Insights into a highly conserved network of hydrogen bonds in the agonist binding site of nicotinic acetylcholine receptors: a structural and theoretical study.

Structural and theoretical studies on the geometrical features of a hydrogen-bond network occurring in the binding site of nicotinic acetylcholine receptors (nAChRs) and composed of interconnected WxPD (Trp-x-Pro-Asp) and SWyz (Ser-Trp-yz) sequences from loops A and B, respectively, have been carried out. Multiple sequence alignments using as template the sequence of the apoform of Aplysia californica acetylcholine binding protein (Ac-AChBP) show the strict conservation of serine and tryptophan residues of the loop B SWyz sequence. Considering a sample of 19 high resolution AChBP structures, the strong conformational preferences of the key tryptophan residue has been pointing out, whatever the form, free or bounded, of AChBP. The geometry of the motif hydrogen-bond network has been characterized through the analyses of seven distances. The robustness of the various hydrogen-bond interactions is pointed out, the one involving the aspartate carboxylate group and the serine residue being the shortest of the network. The role of a cooperative effect involving a NH(His145)…OH (Ser142) hydrogen bond is highlighted. Density functional theory calculations on several simplified models based on the motif hydrogen-bond network allow probing the importance of the various hydrogen-bond interactions. The removal of the Ser142 hydroxyl group induces strong structural rearrangements, in agreement with the structural observations. Molecular electrostatic potential calculations on model systems highlight the importance of a cooperative effect in the whole hydrogen-bond network. More precisely, the key role of the Ser142 hydroxyl group, involved in several hydrogen bonds, is underlined.

Discovery and biological characterization of a novel mesoionic insecticide fenmezoditiaz.

BACKGROUND: Many piercing-sucking insects have developed resistance or cross-resistance to many insecticides targeting insect neural nicotinic acetylcholine receptor (nAChR). Here we are aiming to present the discovery of a novel mesoionic insecticide, fenmezoditiaz, by BASF through structure-based drug design (SBDD) approaches. It has recently been added to the Insecticide Resistance Action Committee mode of classification (IRAC 4E). It is being developed for plant protection against piercing-sucking pests, especially rice hopper complex. RESULTS: The soluble acetylcholine binding protein (AChBP) from the sea slug Aplysia californica was modified using site-directed mutagenesis and based on putative aphid nAChR subunit sequences to create soluble insect-like AChBPs. Among them, insect-like ß1 AChBP and native aphid membrane preparation showed the highest correlated biochemical affinity toward structurally diverse ligands. This mutant AChBP was used to understand how insect nAChRs structurally interact with mesoionics, which was then utilized to design novel mesoionics including fenmezoditiaz. It is an excellent systemic insecticide with diverse application methods and has a broad insecticidal spectrum, especially against piercing/sucking insects. It lacks cross-resistance for neonicotinoid resistant plant hoppers. Field-collected brown plant hopper populations from Asian countries showed high susceptibility. CONCLUSIONS: Fenmezoditiaz is a systemic insecticide with a broad spectrum, lack of cross-resistance and it could be an additional tool for integrated pest management and insecticide resistance management, especially for the rice hopper complex. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

QM-polarized ligand docking accurately predicts the trend in binding affinity of a series of arylmethylene quinuclidine-like derivatives at the α4ß2 and α3ß4 nicotinic acetylcholine receptors (nAChRs).

Compounds containing a quinuclidine scaffold are promising drug candidates for pharmacological management of the central nervous system (CNS) pathologies implicating nAChRs. We have carried out binding affinity and in-silico docking studies of arylmethylene quinuclidine-like derivatives at the α4ß2 receptor using in-vitro receptor binding assay and comparative modeling, respectively. We found that introducing a hydrogen-bond acceptor into the 3-benzylidene quinuclidine derivative resulted in a 266-fold increase in binding affinity and confers agonism properties. By contrast, addition of a phenyl group to 3-benzylidene quinuclidine derivative only results in an 18-fold increase in binding affinity, without conferring agonism. We also found that docking into the orthosteric binding site of the α4ß2 nAChR is consistent with the fact that the basic nitrogen atom donates a hydrogen-bond to the carbonyl group of the highly conserved Trp-149, as initially observed by Dougherty and co-workers.(1) The experimentally-observed trend in binding affinity at both α4ß2 and α3ß4 nAChRs was accurately and independently confirmed by quantum mechanics (QM)-polarized docking. The reduction in binding affinity to the α3ß4 subtype primarily results from a dampening of both coulombic and cation-π interactions.

Computational studies of novel carbonyl-containing diazabicyclic ligands interacting with α4ß2 nicotinic acetylcholine receptor (nAChR) reveal alternative binding modes.

We have carried out computational studies on interactions of diazabicyclic amide analogs with α4ß2 nAChR using homology modeling, docking and pharmacophore elucidation techniques. We have found alternative ligand binding modes in most cases. All these diverse poses exhibit the quintessential hydrogen-bonding interaction between the ligand basic nitrogen and the backbone carbonyl oxygen atom of the highly conserved Trp-149. This hydrogen bond was always found to be shorter than the one contracted by the ligand carbonyl group and a second hydrogen-bond made by the cationic center with Tyr-93 of the principal face of the protein. In most of the poses observed, cation-π interactions involved three aromatic residues located in the principal face of the protein: Trp-149, Tyr-190 and Tyr-197. The latter amino acid residue appears to often donate a hydrogen-bond to the ligand carbonyl oxygen atom. We also describe two rings of alternative receptor-based hydrogen-bond donor features equidistantly separated from the carbonyl oxygen of the highly conserved Trp-149 approximately by 5 and 8Å, respectively. These findings could be exploited to design diverse and selective novel chemical libraries for the treatment of diseases and conditions where the α4ß2 nAChR is disrupted, such as Alzheimer disease, Parkinson's disease and l-dopa-induced dyskinesia (LID).

Unravelling the allosteric binding mode of αD-VxXXB at nicotinic acetylcholine receptors.

αD-conotoxins are 11 kDa homodimers that potently inhibit nicotinic acetylcholine receptors (nAChRs) through a non-competitive (allosteric) mechanism. In this study, we describe the allosteric binding mode of the granulin-like C-terminal (CTD) of VxXXB bound to Lymnea stagnalis acetylcholine binding protein (Ls-AChBP), a soluble homologue of the extracellular ligand-binding domain of nAChRs. This co-crystal complex revealed a novel allosteric binding site for nAChR antagonists outside the C-loop that caps the orthosteric site defined by the nAChR agonist nicotine and the antagonist epibatidine. Mutational and docking studies on Ls-AChBP supported a two-site binding mode for full-length VxXXB, with the first CTD binding site located outside the C-loop as seen in the co-crystal complex, with a second CTD binding site located near the N-terminal end of the adjacent subunit of AChBP. These results provide new structural insight into a novel allosteric mechanism of nAChR inhibition and define the cooperative binding mode of the N-terminal domain linked granulin core domains of αD-conotoxins.

Unique Pharmacological Properties of α-Conotoxin OmIA at α7 nAChRs.

OmIA, isolated from Conus omaria venom, is a potent antagonist at α7 nAChRs. We determined the co-crystal structure of OmIA with Lymnae stagnalis acetylcholine binding protein (Ls-AChBP) that identified His5, Val10 and Asn11 as key determinants for the high potency of OmIA at α7 nAChRs. Remarkably, despite a competitive binding mode observed in the co-crystal structure, OmIA and analogues displayed functional insurmountable antagonism at α7 and α3ß4 nAChRs, except OmIA analogues having long side chain at position 10 ([V10Q]OmIA and [V10L]OmIA), which were partial insurmountable antagonist at α7 nAChRs in the presence of type II positive allosteric modulators (PAMs). A "two-state, two-step" model was used to explain these observations, with [V10Q]OmIA and [V10L]OmIA co-existing in a fast reversible/surmountable as well as a tight binding/insurmountable state. OmIA and analogues also showed biphasic-inhibition at α7 nAChRs in the presence of PNU120596, with a preference for the high-affinity binding site following prolonged exposure. The molecular basis of binding and complex pharmacological profile of OmIA at α7 nAChRs presented in here expands on the potential of α-conotoxins to probe the pharmacological properties of nAChRs and may help guide the development novel α7 modulators.

Nicotinic cholinergic system and COVID-19: In silico evaluation of nicotinic acetylcholine receptor agonists as potential therapeutic interventions.

SARS-CoV-2 infection was announced as a pandemic in March 2020. Since then, several scientists have focused on the low prevalence of smokers among hospitalized COVID-19 patients. These findings led to our hypothesis that the Nicotinic Cholinergic System (NCS) plays a crucial role in the manifestation of COVID-19 and its severe symptoms. Molecular modeling revealed that the SARS-CoV-2 Spike glycoprotein might bind to nicotinic acetylcholine receptors (nAChRs) through a cryptic epitope homologous to snake toxins, substrates well documented and known for their affinity to the nAChRs. This binding model could provide logical explanations for the acute inflammatory disorder in patients with COVID-19, which may be linked to severe dysregulation of NCS. In this study, we present a series of complexes with cholinergic agonists that can potentially prevent SARS-CoV-2 Spike glycoprotein from binding to nAChRs, avoiding dysregulation of the NCS and moderating the symptoms and clinical manifestations of COVID-19. If our hypothesis is verified by in vitro and in vivo studies, repurposing agents currently approved for smoking cessation and neurological conditions could provide the scientific community with a therapeutic option in severe COVID-19.

An Investigation of Three-Finger Toxin-nAChR Interactions through Rosetta Protein Docking.

Three-finger toxins (3FTX) are a group of peptides that affect multiple receptor types. One group of proteins affected by 3FTX are nicotinic acetylcholine receptors (nAChR). Structural information on how neurotoxins interact with nAChR is limited and is confined to a small group of neurotoxins. Therefore, in silico methods are valuable in understanding the interactions between 3FTX and different nAChR subtypes, but there are no established protocols to model 3FTX-nAChR interactions. We followed a homology modeling and protein docking protocol to address this issue and tested its success on three different systems. First, neurotoxin peptides co-crystallized with acetylcholine binding protein (AChBP) were re-docked to assess whether Rosetta protein-protein docking can reproduce the native poses. Second, experimental data on peptide binding to AChBP was used to test whether the docking protocol can qualitatively distinguish AChBP-binders from non-binders. Finally, we docked eight peptides with known α7 and muscle-type nAChR binding properties to test whether the protocol can explain the differential activities of the peptides at the two receptor subtypes. Overall, the docking protocol predicted the qualitative and some specific aspects of 3FTX binding to nAChR with reasonable success and shed light on unknown aspects of 3FTX binding to different receptor subtypes.

Mutagenesis of α-Conotoxins for Enhancing Activity and Selectivity for Nicotinic Acetylcholine Receptors.

Nicotinic acetylcholine receptors (nAChRs) are found throughout the mammalian body and have been studied extensively because of their implication in a myriad of diseases. α-Conotoxins (α-CTxs) are peptide neurotoxins found in the venom of marine snails of genus Conus. α-CTxs are potent and selective antagonists for a variety of nAChR isoforms. Over the past 40 years, α-CTxs have proven to be valuable molecular probes capable of differentiating between closely related nAChR subtypes and have contributed greatly to understanding the physiological role of nAChRs in the mammalian nervous system. Here, we review the amino acid composition and structure of several α-CTxs that selectively target nAChR isoforms and explore strategies and outcomes for introducing mutations in native α-CTxs to direct selectivity and enhance binding affinity for specific nAChRs. This review will focus on structure-activity relationship studies involving native α-CTxs that have been rationally mutated and molecular interactions that underlie binding between ligand and nAChR isoform.

Orthosteric and benzodiazepine cavities of the α1ß2γ2 GABAA receptor: insights from experimentally validated in silico methods.

γ-aminobutyric acid-type A (GABAA) receptors mediate fast synaptic inhibition in the central nervous system of mammals. They are modulated via several sites by numerous compounds, which include GABA, benzodiazepines, ethanol, neurosteroids and anaesthetics among others. Due to their potential as targets of novel drugs, a detailed knowledge of their structure-function relationships is needed. Here, we present the model of the α1ß2γ2 subtype GABAA receptor in the APO state and in complex with selected ligands, including agonists, antagonists and allosteric modulators. The model is based on the crystallographic structure of the human ß3 homopentamer GABAA receptor. The complexes were refined using atomistic molecular dynamics simulations. This allowed a broad description of the binding modes and the detection of important interactions in agreement with experimental information. From the best of our knowledge, this is the only model of the α1ß2γ2 GABAA receptor that represents altogether the desensitized state of the channel and comprehensively describes the interactions of ligands of the orthosteric and benzodiazepines binding sites in agreement with the available experimental data. Furthermore, it is able to explain small differences regarding the binding of a variety of chemically divergent ligands. Finally, this new model may pave the way for the design of focused experimental studies that will allow a deeper description of the receptor.

Molecular interactions of type I and type II positive allosteric modulators with the human α7 nicotinic acetylcholine receptor: an in silico study.

The binding site locations and structural components for type I and type II positive allosteric modulators (PAMs) of the α7 nicotinic acetylcholine receptor (nAChR) have not been fully characterized yet. In this regard, homology models of the human α7 nAChR and hα7/m5-HT3A chimera, built using the crystal structure of the serotonin type 3A receptor (5-ΗΤ3ΑR), were used for molecular docking and molecular dynamics simulations to study the molecular interactions of selected type I (5-hydroxyindol, NS-1738, and LY-2087101) and type II (PNU-120596, PAM-2, and TBS-516) PAMs. The docking results indicate: (1) a site located in the extracellular domain (ECD) for type I PAMs such as NS-1738 and LY-2087101, but not for 5-HI; (2) an overlapping site in the ECD-transmembrane domain (TMD) junction for all studied PAMs. Additional docking results on the hα7/m5-HT3A chimera supported experimental results indicating that the ECD site might be relevant for type I PAM activity; and (3) two TMD sites, an intrasubunit site that recognizes type II PAMs, and an intersubunit pocket with high specificity for 5-HI (type I PAM). The in silico α7TSLMF mutant results support the view that M1-Ser223 and M3-Ile281 are key residues for the interaction of PAM-2 and PNU-120596 with the intrasubunit cavity. Our in silico results are in agreement with experimental data showing that the intrasubunit cavity is relevant for the activity of type II PAMs, and suggest that the ECD-TMD junction and intersubunit sites could be significant for the activity of type I PAMs.

Discovery, cocrystallization and biological evaluation of novel piperidine derivatives as high affinity Ls-AChBP ligands possessing α7 nAChR activities.

A series of novel pyridine-substituted piperidine derivatives were discovered as low nanomolar Ls-AChBP ligands with α7 nAChR partial agonism or antagonism activities. A high-resolution antagonist-bound Ls-AChBP complex was successfully resolved with a classic Loop C opening phenomenon and unique sulfur-π interactions which deviated from our previous docking mode to a large extent. With the knowledge of the co-complex, 27 novel piperidine derivatives were designed and synthesized. The structure-activity relationships (SARs) of the aromatic and pyridine regions were well established and binding modes were illustrated with the help of molecular docking which indicated that interactions with Trp 143 and the "water bridge" are essential for the high binding affinities. Halogen bonding as well as the space around 5'- or 6'- position of the pyridine ring was also proposed to influence the binding conformation of the compounds. Notably, two enantiomers of compound 2 showed opposite functions toward α7 nAChR and compound (S)-2 showed sub-nanomolar affinity (Ki = 0.86 nM) on Ls-AChBP and partial agonism (pEC50 = 4.69 ±â€¯0.11,Emax = 36.1%) on α7 nAChR with reasonable pharmacokinetics (PK) properties and fine ability of blood-brain-barrier (BBB) penetration. This study provided promising hits to develop candidates targeting nAChR-related CNS diseases.

Neurotoxicity fingerprinting of venoms using on-line microfluidic AChBP profiling.

Venoms from snakes are rich sources of highly active proteins with potent affinity towards a variety of enzymes and receptors. Of the many distinct toxicities caused by envenomation, neurotoxicity plays an important role in the paralysis of prey by snakes as well as by venomous sea snails and insects. In order to improve the analytical discovery component of venom toxicity profiling, this paper describes the implementation of microfluidic high-resolution screening (HRS) to obtain neurotoxicity fingerprints from venoms that facilitates identification of the neurotoxic components of envenomation. To demonstrate this workflow, 47 snake venoms were profiled using the acetylcholine binding protein (AChBP) to mimic the target of neurotoxic proteins, in particular nicotinic acetylcholine receptors (nAChRs). In the microfluidic HRS system, nanoliquid chromatographic (nanoLC) separations were on-line connected to both AChBP profiling and parallel mass spectrometry (MS). For virtually all neurotoxic elapid snake venoms tested, we obtained bioactivity fingerprints showing major and minor bioactive zones containing masses consistent with three-finger toxins (3FTxs), whereas, viperid and colubrid venoms showed little or no detectable bioactivity. Our findings demonstrate that venom interactions with AChBP correlate with the severity of neurotoxicity observed following human envenoming by different snake species. We further, as proof of principle, characterized bioactive venom peptides from a viperid (Daboia russelli) and an elapid (Aspidelaps scutatus scutatus) snake by nanoLC-MS/MS, revealing that different toxin classes interact with the AChBP, and that this binding correlates with the inhibition of α7-nAChR in calcium-flux cell-based assays. The on-line post-column binding assay and subsequent toxin characterization methodologies described here provide a new in vitro analytic platform for rapidly investigating neurotoxic snake venom proteins.

Expression of Cys-loop receptor subunits and acetylcholine binding protein in the mechanosensory neurons, glial cells, and muscle tissue of the spider Cupiennius salei.

The central and peripheral nervous system transcriptomes of the spider Cupiennius salei have 15 Cys-loop receptor subunits and an acetylcholine-binding protein (AChBP). Twelve subunits are predicted to form anion channels gated by γ-aminobutyric acid (GABA), glutamate, histamine, or changes in pH, and three are putative ACh-gated cation channels. Spiders have a variety of mechanosensilla and proprioceptive organs that are innervated by efferents in their peripherally located parts, and efferents also innervate muscle fibers. We investigated Cys-loop gene expression in muscle tissue by qPCR and localized this expression in mechanosensilla via in situ hybridization. The cuticular mechanosensory neurons had only CsGABArdl and CspHCl2 subunits, whereas the muscle tissue expressed a wider variety of subunits, especially CsGABAgrd, CsGABAA ß, CsGluCl1 and CspHCl, but very low levels of the CsGABArdl or CsnACh subunits. An nACh non-α subunit was expressed in a group of unidentified cells in the hypodermis and at low level in the muscle tissue, but the physiological function of this subunit is unknown. The CsnAChα subunit was not expressed in sensory neurons and was expressed at extremely low level in the muscle tissue. None of the probes gave signals in proprioceptive joint receptors, suggesting that efferent innervation to this sense organ employs other receptor types. CsAChBP and a glia-specific homeodomain CsREPO were both expressed in glial cells that surround sensory neurons and also in muscle tissue, probably around the nerve endings of the neuromuscular junction. These locations have large numbers of synapses, suggesting that AChBP may have a function in modulating synaptic transmission. J. Comp. Neurol. 525:1139-1154, 2017. © 2016 Wiley Periodicals, Inc.

The binding orientation of epibatidine at α7 nACh receptors.

Epibatidine is an alkaloid toxin that binds with high affinity to nicotinic and muscarinic acetylcholine receptors, and has been extensively used as a research tool. To examine binding interactions at the nicotinic receptor, it has been co-crystallised with the structural homologue acetylcholine binding protein (AChBP; PDB ID 2BYQ), and with an AChBP chimaera (3SQ6) that shares 64% sequence identity with the α7 nACh receptor. However, the binding orientations revealed by AChBP co-crystal structures may not precisely represent their receptor homologues and experimental evidence is needed to verify the ligand poses. Here we identify potential binding site interactions between epibatidine and AChBP residues, and substitute equivalent positions in the α7 nACh receptor. The effects of these are probed by [3H]epibatidine binding following the expression α7 nACh receptor cysteine mutants in HEK 293 cells. Of the sixteen mutants created, the affinity of epibatidine was unaffected by the substitutions Q55C, L106C, L116C, T146C, D160C and S162C, reduced by C186A and C187A, increased by Q114C and S144C, and abolished by W53C, Y91C, N104C, W145C, Y184C and Y191C. These results are consistent with the predicted orientations in AChBP and suggest that epibatidine is likely to occupy a similar location at α7 nACh receptors. We speculate that steric constraints placed upon the C-5 position of the pyridine ring in 3SQ6 may account for the relatively poor affinities of epibatidine derivatives that are substituted at this position.