Use of Acetylcholine Binding Protein in the Search for Novel α7 Nicotinic Receptor Ligands. In Silico Docking, Pharmacological Screening, and X-ray Analysis

Protein Painting Mass Spectrometry in the Discovery of Interaction Sites within the Acetylcholine Binding Protein

Nicotinic acetylcholine receptors (nAChRs) are a family of ligand-gated ion channel receptors that contribute to cognition, memory, and motor control in many organisms. The pharmacological targeting of these receptors, using small molecules or peptides, presents an important strategy for the development of drugs that can treat important human diseases, including neurodegenerative disorders. The Aplysia californica acetylcholine binding protein (Ac-AChBP) is a structural surrogate of the nAChR with high homology to the extracellular ligand binding domain of homopentameric nAChRs. In this study, we optimized protein-painting-based mass spectrometry to identify regions of interaction between the Ac-AChBP and several nAChR ligands. Using molecular dyes that adhere to the surface of a solubilized Ac-AChBP complex, we identified amino acid residues that constitute a contact site within the Ac-AChBP for α-bungarotoxin, choline, nicotine, and amyloid-β 1-42. By integrating innovation in protein painting mass spectrometry with computational structural modeling, we present a new experimental tool for analyzing protein interactions of the nAChR.

Structure and function meet at the nicotinic acetylcholine receptor-lipid interface

The nicotinic acetylcholine receptor (nAChR) is a transmembrane protein that mediates fast intercellular communication in response to the endogenous neurotransmitter acetylcholine. It is the best characterized and archetypal molecule in the superfamily of pentameric ligand-gated ion channels (pLGICs). As a typical transmembrane macromolecule, it interacts extensively with its vicinal lipid microenvironment. Experimental evidence provides a wealth of information on receptor-lipid crosstalk: the nAChR exerts influence on its immediate membrane environment and conversely, the lipid moiety modulates ligand binding, affinity state transitions and gating of ion translocation functions of the receptor protein. Recent cryogenic electron microscopy (cryo-EM) studies have unveiled the occurrence of sites for phospholipids and cholesterol on the lipid-exposed regions of neuronal and electroplax nAChRs, confirming early spectroscopic and affinity labeling studies demonstrating the close contact of lipid molecules with the receptor transmembrane segments. This new data provides structural support to the postulated "lipid sensor" ability displayed by the outer ring of M4 transmembrane domains and their modulatory role on nAChR function, as we postulated a decade ago. Borrowing from the best characterized nAChR, the electroplax (muscle-type) receptor, and exploiting new structural information on the neuronal nAChR, it is now possible to achieve an improved depiction of these sites. In combination with site-directed mutagenesis, single-channel electrophysiology, and molecular dynamics studies, the new structural information delivers a more comprehensive portrayal of these lipid-sensitive loci, providing mechanistic explanations for their ability to modulate nAChR properties and raising the possibility of targetting them in disease.

Speculation on How RIC-3 and Other Chaperones Facilitate α7 Nicotinic Receptor Folding and Assembly

The process of how multimeric transmembrane proteins fold and assemble in the endoplasmic reticulum is not well understood. The alpha7 nicotinic receptor (α7 nAChR) is a good model for multimeric protein assembly since it has at least two independent and specialized chaperones: Resistance to Inhibitors of Cholinesterase 3 (RIC-3) and Nicotinic Acetylcholine Receptor Regulator (NACHO). Recent cryo-EM and NMR data revealed structural features of α7 nAChRs. A ser-ala-pro (SAP) motif precedes a structurally important but unique "latch" helix in α7 nAChRs. A sampling of α7 sequences suggests the SAP motif is conserved from C. elegans to humans, but the latch sequence is only conserved in vertebrates. How RIC-3 and NACHO facilitate receptor subunits folding into their final pentameric configuration is not known. The artificial intelligence program AlphaFold2 recently predicted structures for NACHO and RIC-3. NACHO is highly conserved in sequence and structure across species, but RIC-3 is not. This review ponders how different intrinsically disordered RIC-3 isoforms from C. elegans to humans interact with α7 nAChR subunits despite having little sequence homology across RIC-3 species. Two models from the literature about how RIC-3 assists α7 nAChR assembly are evaluated considering recent structural information about the receptor and its chaperones.

Mode of Action of Neonicotinoid Insecticides Imidacloprid and Thiacloprid to the Cockroach Pameα7 Nicotinic Acetylcholine Receptor

The functional expression of the cockroach Pameα7 nicotinic acetylcholine receptor subunit has been previously studied, and was found to be able to form a homomeric receptor when expressed in Xenopus laevis oocytes. In this study, we found that the neonicotinoid insecticide imidacloprid is unable to activate the cockroach Pameα7 receptor, although thiacloprid induces low inward currents, suggesting that it is a partial agonist. In addition, the co-application or 5 min pretreatment with 10 µM imidacloprid increased nicotine current amplitudes, while the co-application or 5 min pretreatment with 10 µM thiacloprid decreased nicotine-evoked current amplitudes by 54% and 28%, respectively. This suggesting that these two representatives of neonicotinoid insecticides bind differently to the cockroach Pameα7 receptor. Interestingly, the docking models demonstrate that the orientation and interactions of the two insecticides in the cockroach Pameα7 nAChR binding pocket are very similar. Electrophysiological results have provided evidence to suggest that imidacloprid and thiacloprid could act as modulators of the cockroach Pameα7 receptors.

Discovery of fragments inducing conformational effects in dynamic proteins using a second-harmonic generation biosensor

Biophysical screening of compound libraries for the identification of ligands that interact with a protein is efficient, but does typically not reveal if (or how) ligands may interfere with its functional properties. For this a biochemical/functional assay is required. But for proteins whose function is dependent on a conformational change, such assays are typically complex or have low throughput. Here we have explored a high-throughput second-harmonic generation (SHG) biosensor to detect fragments that induce conformational changes upon binding to a protein in real time and identify dynamic regions. Multiwell plate format SHG assays were developed for wild-type and six engineered single-cysteine mutants of acetyl choline binding protein (AChBP), a homologue to ligand gated ion channels (LGICs). They were conjugated with second harmonic-active labels via amine or maleimide coupling. To validate the assay, it was confirmed that the conformational changes induced in AChBP by nicotinic acetyl choline receptor (nAChR) agonists and antagonists were qualitatively different. A 1056 fragment library was subsequently screened against all variants and conformational modulators of AChBP were successfully identified, with hit rates from 9-22%, depending on the AChBP variant. A subset of four hits was selected for orthogonal validation and structural analysis. A time-resolved grating-coupled interferometry-based biosensor assay confirmed the interaction to be a reversible 1-step 1 : 1 interaction, and provided estimates of affinities and interaction kinetic rate constants (K D = 0.28-63 μM, k a = 0.1-6 μM-1 s-1, k d = 1 s-1). X-ray crystallography of two of the fragments confirmed their binding at a previously described conformationally dynamic site, corresponding to the regulatory site of LGICs. These results reveal that SHG has the sensitivity to identify fragments that induce conformational changes in a protein. A selection of fragment hits with a response profile different to known LGIC regulators was characterized and confirmed to bind to dynamic regions of the protein.

A Computational Analysis of the Factors Governing the Dynamics of α7 nAChR and Its Homologs

The α7 nicotinic acetylcholine receptor is a homopentameric ion channel from the Cys-loop receptor superfamily targeted for psychiatric indications and inflammatory pain. Molecular dynamics studies of the receptor have focused on residue mobility and global conformational changes to address receptor function. However, a comparative analysis of α7 with its homologs that cannot trigger channel opening has not been made so far. To identify the residues involved in α7 activation, we ran triplicate 500-ns molecular dynamics simulations with an α7 extracellular domain homology model and two acetylcholine-binding protein homologs. We tested the effect of ligand binding and amino acid sequence on the structure and dynamics of the three proteins. We found that mobile regions identified based on root mean-square deviation and root mean-square fluctuation values are not always consistent among the individual α7 extracellular domain simulations. Comparison of the replica-average properties of the three proteins based on dynamic cross-correlation maps showed that ligand binding affects the coupling between the C-loop and the Cys-loop, vestibular loop, and β1-β2 loops. In addition, the main-immunogenic-region-like domain of α7 went through correlated motions with multiple domains of the receptor. These correlated motions were absent or diminished in α7 homologs, suggesting a unique role in α7 activation.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Di- and heptavalent nicotinic analogues to interfere with α7 nicotinic acetylcholine receptors

In the field of nicotinic acetylcholine receptors (nAChRs), recognized as important therapeutic targets, much effort has been dedicated to the development of nicotinic analogues to agonize or antagonize distinct homo- and heteropentamers nAChR subtypes, selectively. In this work we developed di- and heptavalent nicotinic derivatives based on ethylene glycol (EG) and cyclodextrin cores, respectively. The compounds showed a concentration dependent inhibition of acetylcholine-induced currents on α7 nAChR expressed by Xenopus oocytes. Interesting features were observed with the divalent nicotinic derivatives, acting as antagonists with varied inhibitory concentrations (IC₅₀) in function of the spacer arm length. The best divalent compounds showed a 16-fold lowered IC₅₀ compared to the monovalent reference (12 vs 195 µM). Docking investigations provide guidelines to rationalize these experimental findings.

Identification by virtual screening and functional characterisation of novel positive and negative allosteric modulators of the α7 nicotinic acetylcholine receptor

Several previous studies have demonstrated that the activity of neurotransmitters acting on ligand-gated ion channels such as the nicotinic acetylcholine receptor (nAChR) can be altered by compounds binding to allosteric modulatory sites. In the case of α7 nAChRs, both positive and negative allosteric modulators (PAMs and NAMs) have been identified and have attracted considerable interest. A recent study, employing revised structural models of the transmembrane domain of the α7 nAChR in closed and open conformations, has provided support for an inter-subunit transmembrane allosteric binding site (Newcombe et al 2017). In the present study, we have performed virtual screening of the DrugBank database using pharmacophore queries that were based on the predicted binding mode of PAMs to α7 nAChR structural models. A total of 81 compounds were identified in the DrugBank database, of which the 25 highest-ranked hits corresponded to one of four previously-identified therapeutic compound groups (carbonic anhydrase inhibitors, cyclin-dependent kinase inhibitors, diuretics targeting the Na⁺-K⁺-Cl^- cotransporter, and fluoroquinolone antibiotics targeting DNA gyrase). The top-ranked compound from each of these four groups (DB04763, DB08122, furosemide and pefloxacin, respectively) was tested for its effects on human α7 nAChR expressed in Xenopus oocytes using two-electrode voltage-clamp electrophysiology. These studies, conducted with wild-type, mutant and chimeric receptors, resulted in all four compounds exerting allosteric modulatory effects. While DB04763, DB08122 and pefloxacin were antagonists, furosemide potentiated ACh responses. Our findings, supported by docking studies, are consistent with these compounds acting as PAMs and NAMs of the α7 nAChR via interaction with a transmembrane site.

•

Identification of α7 nAChR positive and negative allosteric modulators.

•

Furosemide is a positive allosteric modulator of α7 nAChRs.

•

DB04763, DB08122 and pefloxacin are negative allosteric modulators of α7 nAChRs.

•

Modulation of α7 nAChRs by an allosteric transmembrane site.

Designing selective modulators for the nicotinic receptor subtypes: challenges and opportunities

Nicotinic receptors are membrane proteins involved in several physiological processes. They are considered suitable drug targets for various CNS disorders or conditions, as shown by the large number of compounds which have entered clinical trials. In recent years, nonconventional agonists have been discovered: positive allosteric modulators, allosteric agonists, site-specific agonists and silent desensitizers are compounds able to modulate the receptor interacting at sites different from the orthodox one, or to desensitize the receptor without prior opening. While these new findings can further complicate the pharmacology of these proteins and the design and optimization of ligands, they undoubtedly offer new opportunities to find drugs for the many therapeutic indications involving nicotinic receptors.

Discovery of peptide ligands through docking and virtual screening at nicotinic acetylcholine receptor homology models

Venom peptide toxins such as conotoxins play a critical role in the characterization of nicotinic acetylcholine receptor (nAChR) structure and function and have potential as nervous system therapeutics as well. However, the lack of solved structures of conotoxins bound to nAChRs and the large size of these peptides are barriers to their computational docking and design. We addressed these challenges in the context of the α4β2 nAChR, a widespread ligand-gated ion channel in the brain and a target for nicotine addiction therapy, and the 19-residue conotoxin α-GID that antagonizes it. We developed a docking algorithm, ToxDock, which used ensemble-docking and extensive conformational sampling to dock α-GID and its analogs to an α4β2 nAChR homology model. Experimental testing demonstrated that a virtual screen with ToxDock correctly identified three bioactive α-GID mutants (α-GID[A10V], α-GID[V13I], and α-GID[V13Y]) and one inactive variant (α-GID[A10Q]). Two mutants, α-GID[A10V] and α-GID[V13Y], had substantially reduced potency at the human α7 nAChR relative to α-GID, a desirable feature for α-GID analogs. The general usefulness of the docking algorithm was highlighted by redocking of peptide toxins to two ion channels and a binding protein in which the peptide toxins successfully reverted back to near-native crystallographic poses after being perturbed. Our results demonstrate that ToxDock can overcome two fundamental challenges of docking large toxin peptides to ion channel homology models, as exemplified by the α-GID:α4β2 nAChR complex, and is extendable to other toxin peptides and ion channels. ToxDock is freely available at rosie.rosettacommons.org/tox_dock.

Exploring the Relationship between Nicotinic Acetylcholine Receptor Ligand Size, Efficiency, Efficacy, and C-Loop Opening

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels mediating fundamental physiological activities in the nervous system and have become important targets for drug design. For a long time, the acetylcholine binding protein (AChBP) has been used as a surrogate to study the nAChR structure-function. Taking advantage of more than 100 AChBP crystal structures in the Protein DataBank (PDB), we explored the relationship between the size, efficiency, and efficacy of nAChR ligands and the C-loop movement. We found that the size of the ligand is correlated with the opening of the C-loop, which can be used in selecting AChBP crystal structures with appropriate C-loop opening to be used for nAChR ligand docking. Ligand size and C-loop opening are reversely correlated with the ligand efficiency rather than the binding affinity. Ligand efficiency could be accurately predicted using simple computational docking, giving a correlation coefficients (R²) up to 0.73. The efficacy of nAChR ligands might be related to ligand size, C-loop opening, and ligand efficiency. Results from this study are useful for engineering the binding affinity and efficacy of nAChR ligands.

The binding orientations of structurally-related ligands can differ; A cautionary note

Crystal structures can identify ligand-receptor interactions and assist the development of novel therapeutics, but experimental challenges sometimes necessitate the use of homologous proteins. Tropisetron is an orthosteric ligand at both 5-HT₃ and α7 nACh receptors and its binding orientation has been determined in the structural homologue AChBP (pdbid: 2WNC). Co-crystallisation with a structurally-related ligand, granisetron, reveals an almost identical orientation (pdbid; 2YME). However, there is a >1000-fold difference in the affinity of tropisetron at 5-HT₃ versus α7 nACh receptors, and α7 nACh receptors do not bind granisetron. These striking pharmacological differences prompt questions about which receptor the crystal structures most closely represent and whether the ligand orientations are correct. Here we probe the binding orientation of tropisetron and granisetron at 5-HT₃ receptors by in silico modelling and docking, radioligand binding on cysteine-substituted 5-HT₃ receptor mutants transiently expressed in HEK 293 cells, and synthetic modification of the ligands. For 15 of the 23 cysteine substitutions, the effects on tropisetron and granisetron were different. Structure-activity relationships on synthesised derivatives of both ligands were also consistent with different orientations, revealing that contrary to the crystallographic evidence from AChBP, the two ligands adopt different orientations in the 5-HT₃ receptor binding site. Our results show that even quite structurally similar molecules can adopt different orientations in the same binding site, and that caution may be needed when using homologous proteins to predict ligand binding.

•

The drugs granisetron and tropisetron are structurally similar.

•

Crystals of them bound to AChBP suggest they have similar binding orientations.

•

At 5-HT₃R, the effects of mutagenesis indicate that their orientations differ.

•

SAR on both of these drugs also supports different orientations.

Insight into the Binding Mode of Agonists of the Nicotinic Acetylcholine Receptor from Calculated Electron Densities

Insect nicotinic acetylcholine receptors (nAChRs) are among the most prominent and most economically important insecticide targets. Thus, an understanding of the modes of binding of respective agonists is important for the design of specific compounds with favorable vertebrate profiles. In the case of nAChRs, the lack of available high-resolution X-ray structures leaves theoretical considerations as the only viable option. Starting from classical homology and docking approaches, binding mode hypotheses are created for five agonists of the nAChR, covering insecticides in the main group 4 of the Insecticide Resistance Action Committee (IRAC) mode of action (MoA) classification, namely, neonicotinoids, nicotine, sulfoxaflor, and butenolides. To better understand these binding modes, the topologies of calculated electron densities of small-model systems are analyzed in the framework of the quantum theory of atoms in molecules. The theoretically obtained modes of binding are very much in line with the biology-driven IRAC MoA classification of the investigated ligands.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Structural basis for cooperative interactions of substituted 2-aminopyrimidines with the acetylcholine binding protein

The nicotinic acetylcholine receptor (nAChR) and the acetylcholine binding protein (AChBP) are pentameric oligomers in which binding sites for nicotinic agonists and competitive antagonists are found at selected subunit interfaces. The nAChR spontaneously exists in multiple conformations associated with its activation and desensitization steps, and conformations are selectively stabilized by binding of agonists and antagonists. In the nAChR, agonist binding and the associated conformational changes accompanying activation and desensitization are cooperative. AChBP, which lacks the transmembrane spanning and cytoplasmic domains, serves as a homology model of the extracellular domain of the nAChRs. We identified unique cooperative binding behavior of a number of 4,6-disubstituted 2-aminopyrimidines to Lymnaea AChBP, with different molecular variants exhibiting positive, nH > 1.0, and negative cooperativity, nH < 1.0. Therefore, for a distinctive set of ligands, the extracellular domain of a nAChR surrogate suffices to accommodate cooperative interactions. X-ray crystal structures of AChBP complexes with examples of each allowed the identification of structural features in the ligands that confer differences in cooperative behavior. Both sets of molecules bind at the agonist-antagonist site, as expected from their competition with epibatidine. An analysis of AChBP quaternary structure shows that cooperative ligand binding is associated with a blooming or flare conformation, a structural change not observed with the classical, noncooperative, nicotinic ligands. Positively and negatively cooperative ligands exhibited unique features in the detailed binding determinants and poses of the complexes.

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.

Marine natural products acting on the acetylcholine-binding protein and nicotinic receptors: from computer modeling to binding studies and electrophysiology

For a small library of natural products from marine sponges and ascidians, in silico docking to the Lymnaea stagnalis acetylcholine-binding protein (AChBP), a model for the ligand-binding domains of nicotinic acetylcholine receptors (nAChRs), was carried out and the possibility of complex formation was revealed. It was further experimentally confirmed via competition with radioiodinated α-bungarotoxin ([¹²⁵I]-αBgt) for binding to AChBP of the majority of analyzed compounds. Alkaloids pibocin, varacin and makaluvamines С and G had relatively high affinities (K(i) 0.5-1.3 μM). With the muscle-type nAChR from Torpedo californica ray and human neuronal α7 nAChR, heterologously expressed in the GH4C1 cell line, no competition with [¹²⁵I]-αBgt was detected in four compounds, while the rest showed an inhibition. Makaluvamines (K(i) ~ 1.5 μM) were the most active compounds, but only makaluvamine G and crambescidine 359 revealed a weak selectivity towards muscle-type nAChR. Rhizochalin, aglycone of rhizochalin, pibocin, makaluvamine G, monanchocidin, crambescidine 359 and aaptamine showed inhibitory activities in electrophysiology experiments on the mouse muscle and human α7 nAChRs, expressed in Xenopus laevis oocytes. Thus, our results confirm the utility of the modeling studies on AChBPs in a search for natural compounds with cholinergic activity and demonstrate the presence of the latter in the analyzed marine biological sources.

Miniaturized bioaffinity assessment coupled to mass spectrometry for guided purification of bioactives from toad and cone snail

A nano-flow high-resolution screening platform, featuring a parallel chip-based microfluidic bioassay and mass spectrometry coupled to nano-liquid chromatography, was applied to screen animal venoms for nicotinic acetylcholine receptor like (nAChR) affinity by using the acetylcholine binding protein, a mimic of the nAChR. The potential of this microfluidic platform is demonstrated by profiling the Conus textile venom proteome, consisting of over 1,000 peptides. Within one analysis (<90 min, 500 ng venom injected), ligands are detected and identified. To show applicability for non-peptides, small molecular ligands such as steroidal ligands were identified in skin secretions from two toad species (Bufo alvarius and Bufo marinus). Bioactives from the toad samples were subsequently isolated by MS-guided fractionation. The fractions analyzed by NMR and a radioligand binding assay with α7-nAChR confirmed the identity and bioactivity of several new ligands.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Insights into the binding of GABA to the insect RDL receptor from atomistic simulations: a comparison of models

The resistance to dieldrin (RDL) receptor is an insect pentameric ligand-gated ion channel (pLGIC). It is activated by the neurotransmitter γ-aminobutyric acid (GABA) binding to its extracellular domain; hence elucidating the atomistic details of this interaction is important for understanding how the RDL receptor functions. As no high resolution structures are currently available, we built homology models of the extracellular domain of the RDL receptor using different templates, including the widely used acetylcholine binding protein and two pLGICs, the Erwinia Chrysanthemi ligand-gated ion channel (ELIC) and the more recently resolved GluCl. We then docked GABA into the selected three dimensional structures, which we used as starting points for classical molecular dynamics simulations. This allowed us to analyze in detail the behavior of GABA in the binding sites, including the hydrogen bond and cation-π interaction networks it formed, the conformers it visited and the possible role of water molecules in mediating the interactions; we also estimated the binding free energies. The models were all stable and showed common features, including interactions consistent with experimental data and similar to other pLGICs; differences could be attributed to the quality of the models, which increases with increasing sequence identity, and the use of a pLGIC template. We supplemented the molecular dynamics information with metadynamics, a rare event method, by exploring the free energy landscape of GABA binding to the RDL receptor. Overall, we show that the GluCl template provided the best models. GABA forming direct salt-bridges with Arg211 and Glu204, and cation-π interactions with an aromatic cage including Tyr109, Phe206 and Tyr254, represents a favorable binding arrangement, and the interaction with Glu204 can also be mediated by a water molecule.

Comparative study on the use of docking and Bayesian categorization to predict ligand binding to nicotinic acetylcholine receptors (nAChRs) subtypes

We have carried out a comparative study between docking into homology models and Bayesian categorization, as applied to virtual screening of nicotinic ligands for binding at various nAChRs subtypes (human and rat α4β2, α7, α3β4, and α6β2β3). We found that although results vary with receptor subtype, Bayesian categorization exhibits higher accuracy and enrichment than unconstrained docking into homology models. However, docking accuracy is improved when one sets up a hydrogen-bond (HB) constraint between the cationic center of the ligand and the main-chain carbonyl group of the conserved Trp-149 or its homologue (a residue involved in cation-π interactions with the ligand basic nitrogen atom). This finding suggests that this HB is a hallmark of nicotinic ligands binding to nAChRs. Best predictions using either docking or Bayesian were obtained with the human α7 nAChR, when 100 nM was used as cutoff for biological activity. We also found that ligand-based Bayesian-derived enrichment factors and structure-based docking-derived enrichment factors highly correlate to each other. Moreover, they correlate with the mean molecular fractional polar surface area of actives ligands and the fractional hydrophobic/hydrophilic surface area of the binding site, respectively. This result is in agreement with the fact that hydrophobicity strongly contributes in promoting nicotinic ligands binding to their cognate nAChRs.

Nicotinic acetylcholine receptor agonists: a milestone for modern crop protection

The destruction of crops by invertebrate pests is a major threat against a background of a continuously rising demand in food supply for a growing world population. Therefore, efficient crop protection measures in a vast range of agricultural settings are of utmost importance to guarantee sustainable yields. The discovery of synthetic agonists selectively addressing the nicotinic acetylcholine receptors (nAChRs), located in the central nervous system of insects, for use as insecticides was a major milestone in applied crop protection research. These compounds, as a result of their high target specificity and versatility in application methods, opened a new innovative era in the control of some of the world's most devastating insect pests. These insecticides also contributed massively to extending our knowledge of the biochemistry of insect nicotinic acetylcholine receptors. The global economic success of synthetic nAChR agonists as insecticides renders the nicotinic acetylcholine receptor still one of the most attractive target sites for exploration in insecticide discovery.

Structure-activity relationships of acetylcholine derivatives with Lucilia cuprina nicotinic acetylcholine receptor α1 and α2 subunits in chicken β2 subunit hybrid receptors in comparison with chicken nicotinic acetylcholine receptor α4/β2

Insect nicotinic acetylcholine (ACh) receptors (nAChRs) are the targets of several insecticide classes. In the present study, we report the gene identification and cloning of nAChR α1 and α2 subunits (Lcα1 and Lcα2) from the sheep blowfly Lucilia cuprina. Xenopus oocytes voltage clamp experiments as hybrids with the chicken β2 nAChR (Ggβ2) subunit resulted in ACh-gated ion channels with distinct dose-response curves for Lcα1/Ggβ2 (effective concentration 50% [EC50 ] = 80 nM; nH  = 1.05), and Lcα2/Ggβ2 (EC50  = 5.37 μM, nH  = 1.46). The neonicotinoid imidacloprid was a potent agonist for the α-bungarotoxin-sensitive Lcα1/Ggβ2 (EC50 ∼ 20 nM), while the α-bungarotoxin-resistant Lcα2/Ggβ2 showed a 30-fold lower sensitivity to this insecticide (EC50  = 0.62 μM). Thirteen close derivatives of ACh were analysed in EC50 , Hill coefficient and maximum current (relative to ACh) determinations for Lcα1/Ggβ2 and Lcα2/Ggβ2 and the chicken Ggα4/Ggβ2 nAChRs, and comparisons relative to ACh allowed the definition of novel structure-activity and structure-selectivity relationships. In the case of N-ethyl-acetylcholine, the EC50 of the chicken Ggα4/Ggβ2 rose by a factor of 1000, while for both Lcα1/Ggβ2 and Lcα2/Ggβ2, potency remained unchanged. Further derivatives with insect nAChR selectivity potential were acetyl-α-methylcholine and trimethyl-(3-methoxy-3-oxopropyl)ammonium, followed by acetylhomocholine and trimethyl-(4-oxopentyl) ammonium. Our results may provide guidance for the identification or design of insect-specific nAChR agonists using structure-based or in silico methods.

Development of a solid-phase receptor-based assay for the detection of cyclic imines using a microsphere-flow cytometry system

Biologically active macrocycles containing a cyclic imine were isolated for the first time from aquaculture sites in Nova Scotia, Canada, during the 1990s. These compounds display a "fast-acting" toxicity in the traditional mouse bioassay for lipophilic marine toxins. Our work aimed at developing a receptor-based detection method for spirolides using a microsphere/flow cytometry Luminex system. For the assay, two alternatives were considered as binding proteins, the Torpedo marmorata nicotinic acetylcholine receptor (nAChR) and the Lymnaea stagnalis acetylcholine binding protein (Ls-AChBP). A receptor-based inhibition assay was developed using the immobilization of nAChR or Ls-AChBP on the surface of carboxylated microspheres and the competition of cyclic imines with biotin-α-bungarotoxin (α-BTX) for binding to these proteins. The amount of biotin-α-BTX bound to the surface of the microspheres was quantified using phycoerythrin (PE)-labeled streptavidin, and the fluorescence was analyzed in a Luminex 200 system. AChBP and nAChR bound to 13-desmethyl spirolide C efficiently; however, the cross-reactivity profile of the nAChR for spirolides and gymnodimine more closely matched the relative toxic potencies reported for these toxins. The nAChR was selected for further assay development. A simple sample preparation protocol consisting of an extraction with acetone yielded a final extract with no matrix interference on the nAChR/microsphere-based assay for mussels, scallops, and clams. This cyclic imine detection method allowed the detection of 13-desmethyl spirolide C in the range of 10-6000 μg/kg of shellfish meat, displaying a higher sensitivity and wider dynamic range than other receptor-based assays previously published. This microsphere-based assay provides a rapid, sensitive, and easily performed screening method that could be multiplexed for the simultaneous detection of several marine toxins.

An integrated catch-and-hold mechanism activates nicotinic acetylcholine receptors

In neuromuscular acetylcholine (ACh) receptor channels (AChRs), agonist molecules bind with a low affinity (LA) to two sites that can switch to high affinity (HA) and increase the probability of channel opening. We measured (by using single-channel kinetic analysis) the rate and equilibrium constants for LA binding and channel gating for several different agonists of adult-type mouse AChRs. Almost all of the variation in the equilibrium constants for LA binding was from differences in the association rate constants. These were consistently below the limit set by diffusion and were substantially different even though the agonists had similar sizes and the same charge. This suggests that binding to resting receptors is not by diffusion alone and, hence, that each binding site can undergo two conformational changes ("catch" and "hold") that connect three different structures (apo-, LA-bound, and HA-bound). Analyses of ACh-binding protein structures suggest that this binding site, too, may adopt three discrete structures having different degrees of loop C displacement ("capping"). For the agonists we tested, the logarithms of the equilibrium constants for LA binding and LA↔HA gating were correlated. Although agonist binding and channel gating have long been considered to be separate processes in the activation of ligand-gated ion channels, this correlation implies that the catch-and-hold conformational changes are energetically linked and together comprise an integrated process having a common structural basis. We propose that loop C capping mainly reflects agonist binding, with its two stages corresponding to the formation of the LA and HA complexes. The catch-and-hold reaction coordinate is discussed in terms of preopening states and thermodynamic cycles of activation.

Exploring chemical space for drug discovery using the chemical universe database

Herein we review our recent efforts in searching for bioactive ligands by enumeration and virtual screening of the unknown chemical space of small molecules. Enumeration from first principles shows that almost all small molecules (>99.9%) have never been synthesized and are still available to be prepared and tested. We discuss open access sources of molecules, the classification and representation of chemical space using molecular quantum numbers (MQN), its exhaustive enumeration in form of the chemical universe generated databases (GDB), and examples of using these databases for prospective drug discovery. MQN-searchable GDB, PubChem, and DrugBank are freely accessible at www.gdb.unibe.ch.

Size matters in activation/inhibition of ligand-gated ion channels

Cys loop, glutamate, and P2X receptors are ligand-gated ion channels (LGICs) with 5, 4, and 3 protomers, respectively. There is now growing atomic level understanding of their gating mechanisms. Although each family is unique in the architecture of the ligand-binding pocket, the pathway for motions to propagate from ligand-binding domain to transmembrane domain, and the gating motions of the transmembrane domain, there are common features among the LGICs, which are the focus of the present review. In particular, agonists and competitive antagonists apparently induce opposite motions of the binding pocket. A simple way to control the motional direction is ligand size. Agonists, usually small, induce closure of the binding pocket, leading to opening of the channel pore, whereas antagonists, usually large, induce opening of the binding pocket, thereby stabilizing the closed pore. A cross-family comparison of the gating mechanisms of the LGICs, focusing in particular on the role played by ligand size, provides new insight on channel activation/inhibition and design of pharmacological compounds.

Identification of novel α7 nicotinic receptor ligands by in silico screening against the crystal structure of a chimeric α7 receptor ligand binding domain

A hierarchical in silico screening procedure using the crystal structure of an agonist bound chimeric α7/Ls-AChBP protein was successfully applied to both proprietary and commercial databases containing drug-like molecules. An overall hit rate of 26% (pK_i ⩾5.0) was obtained, with an even better hit rate of 35% for the commercial compound collection. Structurally novel and diverse ligands were identified. Binding studies with [³H]epibatidine on chimeric α7/5-HT₃ receptors yielded submicromolar inhibition constants for identified hits. Compared to a previous screening procedure that utilized the wild type Ls-AChBP crystal structure, the current study shows that the recently obtained α7/Ls-AChBP chimeric protein crystal structure is a better template for the identification of novel α7 receptor ligands.

The effect of crebanine on memory and cognition impairment via the alpha-7 nicotinic acetylcholine receptor

AIMS

The aims of the present study were to investigate the effect of crebanine on memory and cognition impairment in mice and to elucidate the underlying molecular mechanisms.

MAIN METHODS

The memory-enhancing effects of crebanine were assessed with a water maze test using scopolamine-induced amnesic mice. The molecular mechanism was explored in silico by docking crebanine against acetylcholine binding proteins (AChBPs) and in vitro with a radioligand competition assay using (±)-[(3)H]-epibatidine. The pharmacological behavior was assessed by observing changes to the functional activity of α7-nAChRs expressed in Xenopus oocytes and by fluorescent assays on recombinant ligand gated ion channel (LGIC) receptors expressed in mammalian cells.

KEY FINDINGS

The administration of crebanine significantly improved the cognitive deficits induced by scopolamine, as measured by the water maze test. The docking results demonstrated that crebanine bound to the active binding site of the AChBP template with a good docking energy. Crebanine significantly inhibited the binding of (±)-[(3)H]-epibatidine to AChBPs with K(i) values of 179 nM and 538 nM for Ls and Ac, respectively. Further functional assays performed using two separate protocols indicated that crebanine is an antagonist of the α7-nAChR with an IC(50) of 19.1μM.

SIGNIFICANCE

The observed actions of crebanine against amnesia and its effect on α7-nAChRs will be beneficial for target-based drug design; crebanine or its scaffold can be used as the starting point to develop a drug for Alzheimer's disease. The cognition-enhancing effects of crebanine and the underlying mechanism based on α7-nAChRs are consistent with its traditional use. These findings demonstrate the potential utility of crebanine in the development of neurodegenerative therapy.

Enhanced meta-analysis of acetylcholine binding protein structures reveals conformational signatures of agonism in nicotinic receptors

The soluble acetylcholine binding protein (AChBP) is the default structural proxy for pentameric ligand-gated ion channels (LGICs). Unfortunately, it is difficult to recognize conformational signatures of LGIC agonism and antagonism within the large set of AChBP crystal structures in both apo and ligand-bound states, primarily because AChBP conformations in this set are nearly superimposable (root mean square deviation < 1.5 Å). We have undertaken a systematic, alignment-free approach to elucidate conformational differences displayed by AChBP that cleanly differentiate apo/antagonist-bound from agonist-bound states. Our approach uses statistical inference based on both crystallographic states and conformations sampled during long molecular dynamics simulations to select important inter-C(α) distances and map their collective values onto functional states. We observe that binding of (nAChR) agonists to AChBP elicits clockwise rotation of the inner β-sheet with respect to the outer β-sheet, causing tilting of the cys-loop away from the five-fold axis, in a manner quite similar to that speculated for α-subunits of the heteromeric nAChR structure (Unwin, J Mol Biol 2005;346:967), making this motion potentially important in transmission of the gating signal to the transmembrane domain of a LGIC. The method is also successful at discriminating partial from full agonists and supports the hypothesis that a particularly controversial ligand, lobeline, is in fact an LGIC antagonist.

Intersubunit bridge formation governs agonist efficacy at nicotinic acetylcholine α4β2 receptors: unique role of halogen bonding revealed

The α4β2 subtype of the nicotinic acetylcholine receptor has been pursued as a drug target for treatment of psychiatric and neurodegenerative disorders and smoking cessation aids for decades. Still, a thorough understanding of structure-function relationships of α4β2 agonists is lacking. Using binding experiments, electrophysiology and x-ray crystallography we have investigated a consecutive series of five prototypical pyridine-containing agonists derived from 1-(pyridin-3-yl)-1,4-diazepane. A correlation between binding affinities at α4β2 and the acetylcholine-binding protein from Lymnaea stagnalis (Ls-AChBP) confirms Ls-AChBP as structural surrogate for α4β2 receptors. Crystal structures of five agonists with efficacies at α4β2 from 21-76% were determined in complex with Ls-AChBP. No variation in closure of loop C is observed despite large efficacy variations. Instead, the efficacy of a compound appears tightly coupled to its ability to form a strong intersubunit bridge linking the primary and complementary binding interfaces. For the tested agonists, a specific halogen bond was observed to play a large role in establishing such strong intersubunit anchoring.

Structure-based design, synthesis and structure-activity relationships of dibenzosuberyl- and benzoate-substituted tropines as ligands for acetylcholine-binding protein

Using structure-based optimization procedures on in silico hits, dibenzosuberyl- and benzoate substituted tropines were designed as ligands for acetylcholine-binding protein (AChBP). This protein is a homolog to the ligand binding domain of the nicotinic acetylcholine receptor (nAChR). Distinct SAR is observed between two AChBP species variants and between the α7 and α4β2 nAChR subtype. The AChBP species differences are indicative of a difference in accessibility of a ligand-inducible subpocket. Hereby, we have identified a region that can be scrutinized to achieve selectivity for nicotinic receptor subtypes.