AChBP-targeted alpha-conotoxin correlates distinct binding orientations with nAChR subtype selectivity

Basic Residues at Position 11 of α-Conotoxin LvIA Influence Subtype Selectivity between α3β2 and α3β4 Nicotinic Receptors via an Electrostatic Mechanism

Understanding the determinants of α-conotoxin (α-CTX) selectivity for different nicotinic acetylcholine receptor (nAChR) subtypes is a prerequisite for the design of tool compounds to study nAChRs. However, selectivity optimization of these small, disulfide-rich peptides is difficult not only because of an absence of α-CTX/nAChR co-structures but also because it is challenging to predict how a mutation to an α-CTX will alter its potency and selectivity. As a prototypical system to investigate selectivity, we employed the α-CTX LvIA that is 25-fold selective for the α3β2 nAChR over the related α3β4 nAChR subtype, which is a target for nicotine addiction. Using two-electrode voltage clamp electrophysiology, we identified LvIA[D11R] that is 2-fold selective for the α3β4 nAChR, reversing the subtype preference. This effect is specifically due to the change in charge and not shape of LvIA[D11R], as substitution of D11 with citrulline retains selectivity for the α3β2 nAChR. Furthermore, LvIA[D11K] shows a stronger reversal, with 4-fold selectivity for the α3β4 nAChR. Motivated by these findings, using site-directed mutagenesis, we found that β2[K79A] (I79 on β4), but not β2[K78A] (N78 on β4), largely restores the potency of basic mutants at position 11. Finally, to understand the structural basis of this effect, we used AlphaFold2 to generate models of LvIA in complex with both nAChR subtypes. Both models confirm the plausibility of an electrostatic mechanism to explain the data and also reproduce a broad range of potency and selectivity structure-activity relationships for LvIA mutants, as measured using free energy perturbation simulations. Our work highlights how electrostatic interactions can drive α-CTX selectivity and may serve as a strategy for optimizing the selectivity of LvIA and other α-CTXs.

Cysteine-Rich α-Conotoxin SII Displays Novel Interactions at the Muscle Nicotinic Acetylcholine Receptor

α-Conotoxins that target muscle nicotinic acetylcholine receptors (nAChRs) commonly fall into two structural classes, frameworks I and II containing two and three disulfide bonds, respectively. Conotoxin SII is the sole member of the cysteine-rich framework II with ill-defined interactions at the nAChRs. Following directed synthesis of α-SII, NMR analysis revealed a well-defined structure containing a 3₁₀-helix frequently employed by framework I α-conotoxins; α-SII acted at the muscle nAChR with half-maximal inhibitory concentrations (IC₅₀) of 120 nM (adult) and 370 nM (fetal) though weakly at neuronal nAChRs. Truncation of α-SII to a two disulfide bond amidated peptide with framework I disulfide connectivity led to similar activity. Surprisingly, the more constrained α-SII was less stable under mild reducing conditions and displayed a unique docking mode at the nAChR.

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.

Rigidity of loop 1 contributes to equipotency of globular and ribbon isomers of α-conotoxin AusIA

α-Conotoxins are small disulfide-rich peptides targeting nicotinic acetylcholine receptors (nAChRs) characterised by a C^IC^II-X_m-C^III-X_n-C^IV framework that invariably adopt the native globular conformations which is typically most potent. α-Conotoxins are divided into several structural subgroups based on the number of residues within the two loops braced by the disulfide bonds (m/n), with the 4/7 and 4/3 subgroups dominating. AusIA is a relatively rare α5/5-conotoxin isolated from the venom of Conus australis. Surprisingly, the ribbon isomer displayed equipotency to the wild-type globular AusIA at human α7-containing nAChR. To understand the molecular basis for equipotency, we determined the co-crystal structures of both isomers at Lymnea stagnalis acetylcholine binding protein. The additional residue in the first loop of AusIA was found to be a critical determinant of equipotency, with 11-fold and 86-fold shifts in potency in favour of globular AusIA over ribbon AusIA observed following deletion of Ala4 or Arg5, respectively. This divergence in the potency between globular AusIA and ribbon AusIA was further enhanced upon truncation of the non-conserved Val at the C-termini. Conversely, equipotency could be replicated in LsIA and TxIA [A10L] following insertion of an Ala in the first loop. These findings provide a new understanding of the role the first loop in ribbon and globular α-conotoxins can play in directing α-conotoxin nAChR pharmacology.

α-Conotoxin Bt1.8 from Conus betulinus selectively inhibits α6/α3β2β3 and α3β2 nicotinic acetylcholine receptor subtypes

α-Conotoxins are small disulfide-rich peptides found in the venom of marine cone snails and are potent antagonists of nicotinic acetylcholine receptors (nAChRs). They are valuable pharmacological tools and have potential therapeutic applications for the treatment of chronic pain or neurological diseases and disorders. In the present study, we synthesized and functionally characterized a novel α-conotoxin Bt1.8, which was cloned from Conus betulinus. Bt1.8 selectively inhibited ACh-evoked currents in Xenopus oocytes expressing rat(r) α6/α3β2β3 and rα3β2 nAChRs with an IC₅₀ of 2.1 nM and 9.4 nM, respectively, and similar potency for human (h) α6/α3β2β3 and hα3β2 nAChRs. Additionally, Bt1.8 had higher binding affinity with a slower dissociation rate for the rα6/α3β2β3 subtype compared to rα3β2. The amino acid sequence of Bt1.8 is significantly different from other reported α-conotoxins targeting the two nAChR subtypes. Further Alanine scanning analyses demonstrated that residues Ile9, Leu10, Asn11, Asn12 and Asn14 are critical for its inhibitory activity at the α6/α3β2β3 and α3β2 subtypes. Moreover, the NMR structure of Bt1.8 indicated the presence of a relatively larger hydrophobic zone than other α4/7-conotoxins which may explain its potent inhibition at α6/α3β2β3 nAChRs.

Structure-Function of Neuronal Nicotinic Acetylcholine Receptor Inhibitors Derived From Natural Toxins

Neuronal nicotinic acetylcholine receptors (nAChRs) are prototypical cation-selective, ligand-gated ion channels that mediate fast neurotransmission in the central and peripheral nervous systems. nAChRs are involved in a range of physiological and pathological functions and hence are important therapeutic targets. Their subunit homology and diverse pentameric assembly contribute to their challenging pharmacology and limit their drug development potential. Toxins produced by an extensive range of algae, plants and animals target nAChRs, with many proving pivotal in elucidating receptor pharmacology and biochemistry, as well as providing templates for structure-based drug design. The crystal structures of these toxins with diverse chemical profiles in complex with acetylcholine binding protein (AChBP), a soluble homolog of the extracellular ligand-binding domain of the nAChRs and more recently the extracellular domain of human α9 nAChRs, have been reported. These studies have shed light on the diverse molecular mechanisms of ligand-binding at neuronal nAChR subtypes and uncovered critical insights useful for rational drug design. This review provides a comprehensive overview and perspectives obtained from structure and function studies of diverse plant and animal toxins and their associated inhibitory mechanisms at neuronal nAChRs.

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.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Interactions of the α3β2 Nicotinic Acetylcholine Receptor Interfaces with α-Conotoxin LsIA and its Carboxylated C-terminus Analogue: Molecular Dynamics Simulations

Notably, α-conotoxins with carboxy-terminal (C-terminal) amidation are inhibitors of the pentameric nicotinic acetylcholine receptors (nAChRs), which are therapeutic targets for neurological diseases and disorders. The (α3)₂(β2)₃ nAChR subunit arrangement comprises a pair of α3(+)β2(−) and β2(+)α3(−) interfaces, and a β2(+)β2(−) interface. The β2(+)β2(−) interface has been suggested to have higher agonist affinity relative to the α3(+)β2(−) and β2(+)α3(−) interfaces. Nevertheless, the interactions formed by these subunit interfaces with α-conotoxins are not well understood. Therefore, in order to address this, we modelled the interactions between α-conotoxin LsIA and the α3β2 subtype. The results suggest that the C-terminal carboxylation of LsIA predominantly influenced the enhanced contacts of the conotoxin via residues P7, P14 and C17 on LsIA at the α3(+)β2(−) and β2(+)α3(−) interfaces. However, this enhancement is subtle at the β2(+)β2(−) site, which can compensate the augmented interactions by LsIA at α3(+)β2(−) and β2(+)α3(−) binding sites. Therefore, the divergent interactions at the individual binding interface may account for the minor changes in binding affinity to α3β2 subtype by C-terminal carboxylation of LsIA versus its wild type, as shown in previous experimental results. Overall, these findings may facilitate the development of new drug leads or subtype-selective probes.

ConoMode, a database for conopeptide binding modes

ConoMode is a database for complex three-dimensional (3D) structures of conopeptides binding with their target proteins. Conopeptides, a large family of peptides from the venom of marine snails of the Conus genus, have exceptionally diverse sequences, and their high specificity to block ion channels makes them crucial as drug leads and tools for physiological studies. ConoMode is a specialized archive for the collection of 3D coordinate data for the conopeptides and their binding target proteins from published literature and the Protein Data Bank. These 3D structures can be determined using experimental methods such as X-ray crystallography and electron microscopy and computational methods including docking, homology modeling and molecular dynamics simulations. The binding modes for the conopeptides determined using computational modeling must be validated based on experimental data. The 3D coordinate data from ConoMode can be searched, visualized, downloaded and uploaded. Currently, ConoMode manages 19 conopeptide sequences (from 10 Conus species), 15 protein sequences and 37 3D structures. ConoMode utilizes a modern technical framework to provide a good user experience on mobile devices with touch interaction features. Furthermore, the database is fully optimized for unstructured data and flexible data models.

Database URL: http://conomode.qnlm.ac/conomode/conomode/index

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.

Oligoarginine Peptides, a New Family of Nicotinic Acetylcholine Receptor Inhibitors

Many peptide ligands of nicotinic acetylcholine receptors (nAChRs) contain a large number of positively charged amino acid residues, a striking example being conotoxins RgIA and GeXIVA from marine mollusk venom, with an arginine content of >30%. To determine whether peptides built exclusively from arginine residues will interact with different nAChR subtypes or with their structural homologs such as the acetylcholine-binding protein and ligand-binding domain of the nAChR α9 subunit, we synthesized a series of R3, R6, R8, and R16 oligoarginines and investigated their activity by competition with radioiodinated α-bungarotoxin, two-electrode voltage-clamp electrophysiology, and calcium imaging. R6 and longer peptides inhibited muscle-type nAChRs, α7 nAChRs, and α3β2 nAChRs in the micromolar range. The most efficient inhibition of ion currents was detected for muscle nAChR by R16 (IC₅₀ = 157 nM) and for the α9α10 subtype by R8 and R16 (IC₅₀ = 44 and 120 nM, respectively). Since the R8 affinity for other tested nAChRs was 100-fold lower, R8 appears to be a selective antagonist of α9α10 nAChR. For R8, the electrophysiological and competition experiments indicated the existence of two distinct binding sites on α9α10 nAChR. Since modified oligoarginines and other cationic molecules are widely used as cell-penetrating peptides, we studied several cationic polymers and demonstrated their nAChR inhibitory activity. SIGNIFICANT STATEMENT: By using radioligand analysis, electrophysiology, and calcium imaging, we found that oligoarginine peptides are a new group of inhibitors for muscle nicotinic acetylcholine receptors (nAChRs) and some neuronal nAChRs, the most active being those with 16 and 8 Arg residues. Such compounds and other cationic polymers are cell-penetrating tools for drug delivery, and we also demonstrated the inhibition of nAChRs for several of the latter. Possible positive and negative consequences of such an action should be taken into account.

Periplasmic Expression of 4/7 α-Conotoxin TxIA Analogs in E. coli Favors Ribbon Isomer Formation - Suggestion of a Binding Mode at the α7 nAChR

Peptides derived from animal venoms provide important research tools for biochemical and pharmacological characterization of receptors, ion channels, and transporters. Some venom peptides have been developed into drugs (such as the synthetic ω-conotoxin MVIIA, ziconotide) and several are currently undergoing clinical trials for various clinical indications. Challenges in the development of peptides include their usually limited supply from natural sources, cost-intensive chemical synthesis, and potentially complicated stereoselective disulfide-bond formation in the case of disulfide-rich peptides. In particular, if extended structure–function analysis is performed or incorporation of stable isotopes for NMR studies is required, the comparatively low yields and high costs of synthesized peptides might constitute a limiting factor. Here we investigated the expression of the 4/7 α-conotoxin TxIA, a potent blocker at α3β2 and α7 nicotinic acetylcholine receptors (nAChRs), and three analogs in the form of maltose binding protein fusion proteins in Escherichia coli. Upon purification via nickel affinity chromatography and release of the toxins by protease cleavage, HPLC analysis revealed one major peak with the correct mass for all peptides. The final yield was 1–2 mg of recombinant peptide per liter of bacterial culture. Two-electrode voltage clamp analysis on oocyte-expressed nAChR subtypes demonstrated the functionality of these peptides but also revealed a 30 to 100-fold potency decrease of expressed TxIA compared to chemically synthesized TxIA. NMR spectroscopy analysis of TxIA and two of its analogs confirmed that the decreased activity was due to an alternative disulfide linkage rather than the missing C-terminal amidation, a post-translational modification that is common in α-conotoxins. All peptides preferentially formed in the ribbon conformation rather than the native globular conformation. Interestingly, in the case of the α7 nAChR, but not the α3β2 subtype, the loss of potency could be rescued by an R5D substitution. In conclusion, we demonstrate efficient expression of functional but alternatively folded ribbon TxIA variants in E. coli and provide the first structure–function analysis for a ribbon 4/7-α-conotoxin at α7 and α3β2 nAChRs. Computational analysis based on these data provide evidence for a ribbon α-conotoxin binding mode that might be exploited to design ligands with optimized selectivity.

Crystal Structure of the Monomeric Extracellular Domain of α9 Nicotinic Receptor Subunit in Complex With α-Conotoxin RgIA: Molecular Dynamics Insights Into RgIA Binding to α9α10 Nicotinic Receptors

The α9 subunit of nicotinic acetylcholine receptors (nAChRs) exists mainly in heteropentameric assemblies with α10. Accumulating data indicate the presence of three different binding sites in α9α10 nAChRs: the α9(+)/α9(−), the α9(+)/α10(−), and the α10(+)/α9(−). The major role of the principal (+) side of the extracellular domain (ECD) of α9 subunit in binding of the antagonists methyllylcaconitine and α-bungarotoxin was shown previously by the crystal structures of the monomeric α9-ECD with these molecules. Here we present the 2.26-Å resolution crystal structure of α9-ECD in complex with α-conotoxin (α-Ctx) RgIA, a potential drug for chronic pain, the first structure reported for a complex between an nAChR domain and an α-Ctx. Superposition of this structure with those of other α-Ctxs bound to the homologous pentameric acetylcholine binding proteins revealed significant similarities in the orientation of bound conotoxins, despite the monomeric state of the α9-ECD. In addition, ligand-binding studies calculated a binding affinity of RgIA to the α9-ECD at the low micromolar range. Given the high identity between α9 and α10 ECDs, particularly at their (+) sides, the presented structure was used as template for molecular dynamics simulations of the ECDs of the human α9α10 nAChR in pentameric assemblies. Our results support a favorable binding of RgIA at α9(+)/α9(−) or α10(+)/α9(−) rather than the α9(+)/α10(−) interface, in accordance with previous mutational and functional data.

Effects of C-Terminal Carboxylation on α-Conotoxin LsIA Interactions with Human α7 Nicotinic Acetylcholine Receptor: Molecular Simulation Studies

α-Conotoxins selectively bind to nicotinic acetylcholine receptors (nAChRs), which are therapeutic targets due to their important role in signaling transmission in excitable cells. A previous experimental study has demonstrated that carboxylation of the C-terminal of α-conotoxin LsIA reduces its potency to inhibit human α7 nAChR relative to naturally amidated LsIA. However, little is known about the contribution of conformational changes in the receptor and interactions, induced by C-terminal amidation/carboxylation of conotoxins, to selective binding to nAChRs, since most conotoxins and some disulfide-rich peptides from other conotoxin subfamilies possess a naturally amidated C-terminal. In this study, we employ homology modeling and molecular dynamics (MD) simulations to propose the determinants for differential interactions between amidated and carboxylated LsIAs with α7 nAChR. Our findings indicate an overall increased number of contacts favored by binding of amidated LsIA versus its carboxylated counterpart. Toxin-receptor pairwise interactions, which may play a role in enhancing the potency of the former, include ARG10-TRP77, LEU141 and CYS17-GLN79 via persistent hydrogen bonds and cation-π interactions, which are weakened in the carboxylated form due to a strong intramolecular salt-bridge formed by ARG10 and carboxylated C-terminus. The binding of amidated LsIA also induces enhanced movements in loop C and the juxtamembrane Cys-loop that are closely associated with receptor function. Additionally, the impacts of binding of LsIA on the overall structure and inter-subunit contacts were examined using inter-residue network analysis, suggesting a clockwise tilting of the α7 C and F loops upon binding to carboxylated LsIA, which is absent for amidated LsIA binding. The predicted molecular mechanism of LsIA binding to the α7 receptor may provide new insights into the important role of the C-terminal in the binding potency of conotoxins at neuronal nAChRs for pharmacological purposes.

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.

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 (K_i = 0.86 nM) on Ls-AChBP and partial agonism (pEC₅₀ = 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.

Orthosteric and/or Allosteric Binding of α-Conotoxins to Nicotinic Acetylcholine Receptors and Their Models

α-Conotoxins from Conus snails are capable of distinguishing muscle and neuronal nicotinic acetylcholine receptors (nAChRs). α-Conotoxin RgIA and αO-conotoxin GeXIVA, blocking neuronal α9α10 nAChR, are potential analgesics. Typically, α-conotoxins bind to the orthosteric sites for agonists/competitive antagonists, but αO-conotoxin GeXIVA was proposed to attach allosterically, judging by electrophysiological experiments on α9α10 nAChR. We decided to verify this conclusion by radioligand analysis in competition with α-bungarotoxin (αBgt) on the ligand-binding domain of the nAChR α9 subunit (α9 LBD), where, from the X-ray analysis, αBgt binds at the orthosteric site. A competition with αBgt was registered for GeXIVA and RgIA, IC₅₀ values being in the micromolar range. However, high nonspecific binding of conotoxins (detected with their radioiodinated derivatives) to His₆-resin attaching α9 LBD did not allow us to accurately measure IC₅₀s. However, IC₅₀s were measured for binding to Aplysia californica AChBP: the RgIA globular isomer, known to be active against α9α10 nAChR, was more efficient than the ribbon one, whereas all three GeXIVA isomers had similar potencies at low µM. Thus, radioligand analysis indicated that both conotoxins can attach to the orthosteric sites in these nAChR models, which should be taken into account in the design of analgesics on the basis of these conotoxins.

Molecular determinants of α-conotoxin potency for inhibition of human and rat α6β4 nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChRs) containing α6 and β4 subunits are expressed by dorsal root ganglion neurons and have been implicated in neuropathic pain. Rodent models are often used to evaluate the efficacy of analgesic compounds, but species differences may affect the activity of some nAChR ligands. A previous candidate α-conotoxin-based therapeutic yielded promising results in rodent models, but failed in human clinical trials, emphasizing the importance of understanding species differences in ligand activity. Here, we show that human and rat α6/α3β4 nAChRs expressed in Xenopus laevis oocytes exhibit differential sensitivity to α-conotoxins. Sequence homology comparisons of human and rat α6β4 nAChR subunits indicated that α6 residues forming the ligand-binding pocket are highly conserved between the two species, but several residues of β4 differed, including a Leu-Gln difference at position 119. X-ray crystallography of α-conotoxin PeIA complexed with the Aplysia californica acetylcholine-binding protein (AChBP) revealed that binding of PeIA orients Pro¹³ in close proximity to residue 119 of the AChBP complementary subunit. Site-directed mutagenesis studies revealed that Leu¹¹⁹ of human β4 contributes to higher sensitivity of human α6/α3β4 nAChRs to α-conotoxins, and structure-activity studies indicated that PeIA Pro¹³ is critical for high potency. Human and rat α6/α3β4 nAChRs displayed differential sensitivities to perturbations of the interaction between PeIA Pro¹³ and residue 119 of the β4 subunit. These results highlight the potential significance of species differences in α6β4 nAChR pharmacology that should be taken into consideration when evaluating the activity of candidate human therapeutics in rodent models.

Stoichiometry dependent inhibition of rat α3β4 nicotinic acetylcholine receptor by the ribbon isomer of α-conotoxin AuIB

The ribbon isomer of α-conotoxin AuIB has 10-fold greater potency than the wild-type globular isomer at inhibiting nicotinic acetylcholine receptors (nAChRs) in rat parasympathetic neurons, and unlike its globular isoform, ribbon AuIB only targets a specific stoichiometry of the α3β4 nAChR subtype. Previous electrophysiological recordings of AuIB indicated that ribbon AuIB binds to the α3(+)α3(-) interface within the nAChR extracellular domain, which is displayed by the (α3)₃(β4)₂ stoichiometry but not by (α3)₂(β4)₃. This specificity for a particular stoichiometry is remarkable and suggests that ribbon isoforms of α-conotoxins might have great potential in drug design. In this study, we investigated the binding mode and structure-activity relationships of ribbon AuIB using a combination of molecular modeling and electrophysiology recording to determine the features that underpin its selectivity. An alanine scan showed that positions 4 and 9 of ribbon AuIB are the main determinants of the interaction with (α3)₃(β4)₂ nAChR. Our computational models indicate that the first loop of ribbon AuIB binds in the "aromatic box" of the acetylcholine orthosteric binding site, similar to that of globular AuIB. In contrast, the second loop and the termini of the ribbon isomer have different orientations and interactions in the binding sites to those of the globular isomer. The structure-activity relationships reported herein should be useful to design peptides displaying a ribbon α-conotoxin scaffold for inhibition of nAChR subtypes that have hitherto been difficult to selectively target.

Characterization of Disulfide Linkages in Proteins by 193 nm Ultraviolet Photodissociation (UVPD) Mass Spectrometry

Deciphering disulfide bond patterns in proteins remains a significant challenge. In the present study, interlinked disulfide bonds connecting peptide chains are homolytically cleaved with 193 nm ultraviolet photodissociation (UVPD). Analysis of insulin showcased the ability of UVPD to cleave multiple disulfide bonds and provide sequence coverage of the peptide chains in the same MS/MS event. For proteins containing more complex disulfide bonding patterns, an approach combining partial reduction and alkylation mitigated disulfide scrambling and allowed assignment of the array of disulfide bonds. The 4 disulfide bonds of lysozyme and the 19 disulfide bonds of serotransferrin were characterized through LC/UVPD-MS analysis of nonreduced and partially reduced protein digests.

Effect of Methionine Oxidation and Substitution of α-Conotoxin TxID on α3β4 Nicotinic Acetylcholine Receptor

α-Conotoxin TxID was discovered from Conus textile by gene cloning, which has 4/6 inter-cysteine loop spacing and selectively inhibits α3β4 nicotinic acetylcholine receptor (nAChR) subtype. However, TxID is susceptible to modification due to it containing a methionine (Met) residue that easily forms methionine sulfoxide (MetO) in oxidative environment. In this study, we investigated how Met-11 and its derivatives affect the activity of TxID using a combination of electrophysiological recordings and molecular modelling. The results showed most TxID analogues had substantially decreased activities on α3β4 nAChR with more than 10-fold potency loss and 5 of them demonstrated no inhibition on α3β4 nAChR. However, one mutant, [M11I]TxID, displayed potent inhibition at α3β4 nAChR with an IC₅₀ of 69 nM, which only exhibited 3.8-fold less compared with TxID. Molecular dynamics simulations were performed to expound the decrease in the affinity for α3β4 nAChR. The results indicate replacement of Met with a hydrophobic moderate-sized Ile in TxID is an alternative strategy to reduce the impact of Met oxidation, which may help to redesign conotoxins containing methionine residue.

Neuronal Nicotinic Acetylcholine Receptor Modulators from Cone Snails

Marine cone snails are a large family of gastropods that have evolved highly potent venoms for predation and defense. The cone snail venom has exceptional molecular diversity in neuropharmacologically active compounds, targeting a range of receptors, ion channels, and transporters. These conotoxins have helped to dissect the structure and function of many of these therapeutically significant targets in the central and peripheral nervous systems, as well as unravelling the complex cellular mechanisms modulated by these receptors and ion channels. This review provides an overview of α-conotoxins targeting neuronal nicotinic acetylcholine receptors. The structure and activity of both classical and non-classical α-conotoxins are discussed, along with their contributions towards understanding nicotinic acetylcholine receptor (nAChR) structure and function.

Venom based neural modulators

Different types of neuronal nicotinic acetylcholine receptors (nAChRs) are expected to occur in vivo, most structure-activity relationship studies have been carried out for just a few neuronal subtypes. The present review enlightens current aspects of venom modulators of nAChRs. Important electronic databases such as PubMed or Google scholar were explored for the collection of latest studies in the field. Clinical and basic research has shown that cholinergic receptors play a role in several disorders of the nervous system such as chronic pain, Alzheimers disease and addiction to nicotine, alcohol and drugs. Unfortunately, the lack of selective modulators for each subtype of nAChR makes their pharmacological characterization difficult, which has slowed the development of therapeutic nAChR modulators with high selectivity and absence of off-target side-effects. Animal venoms have proven to be an excellent natural source of bioactive molecules with activity against ion channels. The present review concludes that the presence of small-molecule nAChR modulators in spider venoms support the use of venoms as a potential source of novel modulators.

α-Conotoxins to explore the molecular, physiological and pathophysiological functions of neuronal nicotinic acetylcholine receptors

The vast diversity of neuronal nicotinic acetylcholine subunits expressed in the central and peripheral nervous systems, as well as in non-neuronal tissues, constitutes a formidable challenge for researchers and clinicians to decipher the role of particular subtypes, including complex subunit associations, in physiological and pathophysiological functions. Many natural products target the nAChRs, but there is no richer source of nicotinic ligands than the venom of predatory gastropods known as cone snails. Indeed, every single species of cone snail was shown to produce at least one type of such α-conotoxins. These tiny peptides (10-25 amino acids), constrained by disulfide bridges, proved to be unvaluable tools to investigate the structure and function of nAChRs, some of them having also therapeutic potential. In this review, we provide a recent update on the pharmacology and subtype specificity of several major α-conotoxins.

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.

Modelling the interactions between animal venom peptides and membrane proteins

The active components of animal venoms are mostly peptide toxins, which typically target ion channels and receptors of both the central and peripheral nervous system, interfering with action potential conduction and/or synaptic transmission. The high degree of sequence conservation of their molecular targets makes a range of these toxins active at human receptors. The high selectivity and potency displayed by some of these toxins have prompted their use as pharmacological tools as well as drugs or drug leads. Molecular modelling has played an essential role in increasing our molecular-level understanding of the activity and specificity of animal toxins, as well as engineering them for biotechnological and pharmaceutical applications. This review focuses on the biological insights gained from computational and experimental studies of animal venom toxins interacting with membranes and ion channels. A host of recent X-ray crystallography and electron-microscopy structures of the toxin targets has contributed to a dramatic increase in the accuracy of the molecular models of toxin binding modes greatly advancing this exciting field of study. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

α-Conotoxin [S9A]TxID Potently Discriminates between α3β4 and α6/α3β4 Nicotinic Acetylcholine Receptors

α3β4 nAChRs have been implicated in various pathophysiological conditions. However, the expression profile of α3β4 nAChRs and α6/α3β4 nAChRs overlap in a variety of tissues. To distinguish between these two subtypes, we redesigned peptide 1 (α-conotoxin TxID), which inhibits α3β4 and α6/α3β4 nAChR subtypes. We systematically mutated 1 to evaluate analogue selectivity for α3β4 vs α6/α3β4 nAChRs expressed in Xenopus laevis oocytes. One analogue, peptide 7 ([S9A]TxID), had 46-fold greater potency for α3β4 versus α6/α3β4 nAChRs. Peptide 7 had IC₅₀s > 10 μM for other nAChR subtypes. Molecular dynamics simulations suggested that Ser-9 of TxID was involved in a weak hydrogen bond with β4 Lys-81 in the α6β4 binding site but not in the α3β4 binding site. When Ser-9 was substituted by an Ala, this hydrogen bond interaction was disrupted. These results provide further molecular insights into the selectivity of 7 and provide a guide for designing ligands that block α3β4 nAChRs.

Nicotinic acetylcholine receptor inhibitors derived from snake and snail venoms

The nicotinic acetylcholine receptor (nAChR) represents the prototype of ligand-gated ion channels. It is vital for neuromuscular transmission and an important regulator of neurotransmission. A variety of toxic compounds derived from diverse species target this receptor and have been of elemental importance in basic and applied research. They enabled milestone discoveries in pharmacology and biochemistry ranging from the original formulation of the receptor concept, the first isolation and structural analysis of a receptor protein (the nAChR) to the identification, localization, and differentiation of its diverse subtypes and their validation as a target for therapeutic intervention. Among the venom-derived compounds, α-neurotoxins and α-conotoxins provide the largest families and still represent indispensable pharmacological tools. Application of modified α-neurotoxins provided substantial structural and functional details of the nAChR long before high resolution structures were available. α-bungarotoxin represents not only a standard pharmacological tool and label in nAChR research but also for unrelated proteins tagged with a minimal α-bungarotoxin binding motif. A major advantage of α-conotoxins is their smaller size, as well as superior selectivity for diverse nAChR subtypes that allows their development into ligands with optimized pharmacological and chemical properties and potentially novel drugs. In the following, these two groups of nAChR antagonists will be described focusing on their respective roles in the structural and functional characterization of nAChRs and their development into research tools. In addition, we provide a comparative overview of the diverse α-conotoxin selectivities that can serve as a practical guide for both structure activity studies and subtype classification. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

The crystal structure of Ac-AChBP in complex with α-conotoxin LvIA reveals the mechanism of its selectivity towards different nAChR subtypes

The α3* nAChRs, which are considered to be promising drug targets for problems such as pain, addiction, cardiovascular function, cognitive disorders etc., are found throughout the central and peripheral nervous system. The α-conotoxin (α-CTx) LvIA has been identified as the most selective inhibitor of α3β2 nAChRs known to date, and it can distinguish the α3β2 nAChR subtype from the α6/α3β2β3 and α3β4 nAChR subtypes. However, the mechanism of its selectivity towards α3β2, α6/α3β2β3, and α3β4 nAChRs remains elusive. Here we report the co-crystal structure of LvIA in complex with Aplysia californica acetylcholine binding protein (Ac-AChBP) at a resolution of 3.4 Å. Based on the structure of this complex, together with homology modeling based on other nAChR subtypes and binding affinity assays, we conclude that Asp-11 of LvIA plays an important role in the selectivity of LvIA towards α3β2 and α3/α6β2β3 nAChRs by making a salt bridge with Lys-155 of the rat α3 subunit. Asn-9 lies within a hydrophobic pocket that is formed by Met-36, Thr-59, and Phe-119 of the rat β2 subunit in the α3β2 nAChR model, revealing the reason for its more potent selectivity towards the α3β2 nAChR subtype. These results provide molecular insights that can be used to design ligands that selectively target α3β2 nAChRs, with significant implications for the design of new therapeutic α-CTxs.

Structural mechanisms for α-conotoxin activity at the human α3β4 nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors (nAChR) are therapeutic targets for a range of human diseases. α-Conotoxins are naturally occurring peptide antagonists of nAChRs that have been used as pharmacological probes and investigated as drug leads for nAChR related disorders. However, α-conotoxin interactions have been mostly characterised at the α7 and α3β2 nAChRs, with interactions at other subtypes poorly understood. This study provides novel structural insights into the molecular basis for α-conotoxin activity at α3β4 nAChR, a therapeutic target where subtype specific antagonists have potential to treat nicotine addiction and lung cancer. A co-crystal structure of α-conotoxin LsIA with Lymnaea stagnalis acetylcholine binding protein guided the design and functional characterisations of LsIA analogues that identified the minimum pharmacophore regulating α3β4 antagonism. Interactions of the LsIA R10F with β4 K57 and the conserved -NN- α-conotoxin motif with β4 I77 and I109 conferred α3β4 activity to the otherwise inactive LsIA. Using these structural insights, we designed LsIA analogues with α3β4 activity. This new understanding of the structural basis of protein-protein interactions between α-conotoxins and α3β4 may help rationally guide the development of α3β4 selective antagonists with therapeutic potential.

A Strategy for Production of Correctly Folded Disulfide-Rich Peptides in the Periplasm of E. coli

Recombinant expression of disulfide-reticulated peptides and proteins is often challenging. We describe a method that exploits the periplasmic disulfide-bond forming machinery of Escherichia coli and combines this with a cleavable, solubility-enhancing fusion tag to obtain higher yields of correctly folded target protein than is achievable via cytoplasmic expression. The protocols provided herein cover all aspects of this approach, from vector construction and transformation to purification of the cleaved target protein and subsequent quality control.

Genetic Algorithm Managed Peptide Mutant Screening: Optimizing Peptide Ligands for Targeted Receptor Binding

This study demonstrates the utility of genetic algorithms to search exceptionally large and otherwise intractable mutant libraries for sequences with optimal binding affinities for target receptors. The Genetic Algorithm Managed Peptide Mutant Screening (GAMPMS) program was used to search an α-conotoxin (α-CTx) MII mutant library of approximately 41 billion possible peptide sequences for those exhibiting the greatest binding affinity for the α₃β₂-nicotinic acetylcholine receptor (nAChR) isoform. A series of top resulting peptide ligands with high sequence homology was obtained, with each mutant having an estimated ΔG_bind approximately double that of the potent native α-CTx MII ligand. A consensus sequence from the top GAMPMS results was subjected to more rigorous binding free energy calculations by molecular dynamics and compared to α-CTx MII and other related variants for binding with α₃β₂-nAChR. In this study, the efficiency of GAMPMS to substantially reduce the sample population size through evolutionary selection criteria to produce ligands with higher predicted binding affinity is demonstrated.

High-Affinity α-Conotoxin PnIA Analogs Designed on the Basis of the Protein Surface Topography Method

Despite some success for small molecules, elucidating structure-function relationships for biologically active peptides - the ligands for various targets in the organism - remains a great challenge and calls for the development of novel approaches. Some of us recently proposed the Protein Surface Topography (PST) approach, which benefits from a simplified representation of biomolecules' surface as projection maps, which enables the exposure of the structure-function dependencies. Here, we use PST to uncover the "activity pattern" in α-conotoxins - neuroactive peptides that effectively target nicotinic acetylcholine receptors (nAChRs). PST was applied in order to design several variants of the α-conotoxin PnIA, which were synthesized and thoroughly studied. Among the best was PnIA[R9, L10], which exhibits nanomolar affinity for the α7 nAChR, selectivity and a slow wash-out from this target. Importantly, these mutations could hardly be delineated by "standard" structure-based drug design. The proposed combination of PST with a set of experiments proved very efficient for the rational construction of new bioactive molecules.

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.

Escherichia coli Protein Expression System for Acetylcholine Binding Proteins (AChBPs)

Nicotinic acetylcholine receptors (nAChR) are ligand gated ion channels, identified as therapeutic targets for a range of human diseases. Drug design for nAChR related disorders is increasingly using structure-based approaches. Many of these structural insights for therapeutic lead development have been obtained from co-crystal structures of nAChR agonists and antagonists with the acetylcholine binding protein (AChBP). AChBP is a water soluble, structural and functional homolog of the extracellular, ligand-binding domain of nAChRs. Currently, AChBPs are recombinantly expressed in eukaryotic expression systems for structural and biophysical studies. Here, we report the establishment of an Escherichia coli (E. coli) expression system that significantly reduces the cost and time of production compared to the existing expression systems. E. coli can efficiently express unglycosylated AChBP for crystallography and makes the expression of isotopically labelled forms feasible for NMR. We used a pHUE vector containing an N-terminal His-tagged ubiquitin fusion protein to facilitate AChBP expression in the soluble fractions, and thus avoid the need to recover protein from inclusion bodies. The purified protein yield obtained from the E. coli expression system is comparable to that obtained from existing AChBP expression systems. E. coli expressed AChBP bound nAChR agonists and antagonists with affinities matching those previously reported. Thus, the E. coli expression system significantly simplifies the expression and purification of functional AChBP for structural and biophysical studies.

Design and synthesis of α-conotoxin GID analogues as selective α4β2 nicotinic acetylcholine receptor antagonists

The α4β2 nicotinic acetylcholine receptor (nAChR) is an important target for currently approved smoking cessation therapeutics. However, the development of highly selective α4β2 nAChR antagonists remains a significant challenge. α-Conotoxin GID is an antagonist of α4β2 nAChRs, though it is significantly more potent toward the α3β2 and α7 subtypes. With the goal of obtaining further insights into α-conotoxin GID/nAChR interactions that could lead to the design of GID analogues with improved affinity for α4β2 nAChRs, we built a homology model of the GID/α4β2 complex using an X-ray co-crystal structure of an α-conotoxin/acetylcholine binding protein (AChBP) complex. Several additional interactions that could potentially enhance the affinity of GID for α4β2 nAChRs were observed in our model, which led to the design and synthesis of 22 GID analogues. Seven analogues displayed inhibitory activity toward α4β2 nAChRs that was comparable to GID. Significantly, both GID[A10S] and GID[V13I] demonstrated moderately improved selectivity toward α4β2 over α3β2 when compared with GID, while GID[V18N] exhibited no measurable inhibitory activity for the α3β2 subtype, yet retained inhibitory activity for α4β2. In this regard, GID[V18N] is the most α4β2 nAChR selective α-conotoxin analogue identified to date.

From crystal structure of α-conotoxin GIC in complex with Ac-AChBP to molecular determinants of its high selectivity for α3β2 nAChR

Acetylcholine binding proteins (AChBPs) are unique spatial homologs of the ligand-binding domains of nicotinic acetylcholine receptors (nAChRs), and they reproduce some pharmacological properties of nAChRs. X-ray crystal structures of AСhBP in complex with α-conotoxins provide important insights into the interactions of α-conotoxins with distinct nAChR subtypes. Although considerable efforts have been made to understand why α-conotoxin GIC is strongly selective for α3β2 nAChR, this question has not yet been solved. Here we present the structure of α-conotoxin GIC in complex with Aplysia californica AChBP (Ac-AChBP) at a resolution of 2.1 Å. Based on this co-crystal structure complemented with molecular docking data, we suggest the key residues of GIC in determining its high affinity and selectivity for human α3β2 vs α3β4 nAChRs. These suggestions were checked by radioligand and electrophysiology experiments, which confirmed the functional role of detected contacts for GIC interactions with Ac-AChBP and α3β2 nAChR subtypes. While GIC elements responsible for its high affinity binding with Ac-AChBP and α3β2 nAChR were identified, our study also showed the limitations of computer modelling in extending the data from the X-ray structures of the AChBP complexes to all nAChR subtypes.

The role of defensive ecological interactions in the evolution of conotoxins

Venoms comprise of complex mixtures of peptides evolved for predation and defensive purposes. Remarkably, some carnivorous cone snails can inject two distinct venoms in response to predatory or defensive stimuli, providing a unique opportunity to study separately how different ecological pressures contribute to toxin diversification. Here, we report the extraordinary defensive strategy of the Rhizoconus subgenus of cone snails. The defensive venom from this worm-hunting subgenus is unusually simple, almost exclusively composed of αD-conotoxins instead of the ubiquitous αA-conotoxins found in the more complex defensive venom of mollusc- and fish-hunting cone snails. A similarly compartmentalized venom gland as those observed in the other dietary groups facilitates the deployment of this defensive venom. Transcriptomic analysis of a Conus vexillum venom gland revealed the αD-conotoxins as the major transcripts, with lower amounts of 15 known and four new conotoxin superfamilies also detected with likely roles in prey capture. Our phylogenetic and molecular evolution analysis of the αD-conotoxins from five subgenera of cone snails suggests they evolved episodically as part of a defensive strategy in the Rhizoconus subgenus. Thus, our results demonstrate an important role for defence in the evolution of conotoxins.

Recombinant Expression and Characterization of α-Conotoxin LvIA in Escherichia coli

α-Conotoxin LvIA is derived from Conus lividus, native to Hainan, and is the most selective inhibitor of α3β2 nicotinic acetylcholine receptors (nAChRs) known to date. In this study, an efficient approach for the production of recombinant α-Conotoxin LvIA is described. Tandem repeats of a LvIA gene fragment were constructed and fused with a KSI gene and a His₆ tag in a Escherichia coli (E. coli) expression vector pET-31b(+). The recombinant plasmids were transformed into E. coli and were found to express well. The KSI-(LvIA)_n-His₆ fusion protein was purified by metal affinity chromatography and then cleaved with CNBr to release recombinant LvIA (rLvIA). High yields of fusion protein ranging from 100 to 500 mg/L culture were obtained. The pharmacological profile of rLvIA was determined by two-electrode voltage-clamp electrophysiology in Xenopus laevis oocytes expressing rat nAChR subtypes. The rLvIA antagonized the α3β2 nAChR subtype selectively with a nano-molar IC₅₀. The rLvIA was analgesic in a mouse hot-plate test model of pain. Overall, this study provides an effective method to synthesize α-conotoxin LvIA in an E. coli recombinant expression system, and this approach could be useful to obtain active conopeptides in large quantity and at low cost.

Computational exploration of a protein receptor binding space with student proposed peptide ligands

Computational molecular docking is a fast and effective in silico method for the analysis of binding between a protein receptor model and a ligand. The visualization and manipulation of protein to ligand binding in three-dimensional space represents a powerful tool in the biochemistry curriculum to enhance student learning. The DockoMatic tutorial described herein provides a framework by which instructors can guide students through a drug screening exercise. Using receptor models derived from readily available protein crystal structures, docking programs have the ability to predict ligand binding properties, such as preferential binding orientations and binding affinities. The use of computational studies can significantly enhance complimentary wet chemical experimentation by providing insight into the important molecular interactions within the system of interest, as well as guide the design of new candidate ligands based on observed binding motifs and energetics. In this laboratory tutorial, the graphical user interface, DockoMatic, facilitates docking job submissions to the docking engine, AutoDock 4.2. The purpose of this exercise is to successfully dock a 17-amino acid peptide, α-conotoxin TxIA, to the acetylcholine binding protein from Aplysia californica-AChBP to determine the most stable binding configuration. Each student will then propose two specific amino acid substitutions of α-conotoxin TxIA to enhance peptide binding affinity, create the mutant in DockoMatic, and perform docking calculations to compare their results with the class. Students will also compare intermolecular forces, binding energy, and geometric orientation of their prepared analog to their initial α-conotoxin TxIA docking results.

A novel 4/6-type alpha-conotoxin ViIA selectively inhibits nAchR α3β2 subtype

Conotoxins (CTxs) are typically small peptides composed of 12-50 amino acid residues with 2-5 disulfide bridges. Most of them potently and selectively target a wide variety of ion channels and membrane receptors. They are highly valued as neuropharmacological probes and in pharmaceutical development. In this work, a novel α4/6-CTx named ViIA (RDCCSNPPCAHNNPDC-NH2) was identified from a cDNA library of the venom ducts of Conus virgo (C. virgo). ViIA was then synthesized chemically and its disulfide connectivity was identified as 'C(1)-C(3), C(2)-C(4)'. Its molecular targets were further assessed using two-electrode voltage clamping. The results indicated that ViIA selectively inhibited nicotinic acetylcholine receptor (nAChR) α3β2 subtype with an IC50 of 845.5 nM, but did not target dorsal root ganglion sodium (Na(+))-, potassium (K(+))- or calcium (Ca(2+))-ion channels. Further structure-activity relationship analysis demonstrated that Arg(1) and His(11) but not Asp(2) were the functional residues. To the best of our knowledge, ViIA is the first 4/6 α-CTx that selectively inhibits nAChR α3β2 subtype. This finding expands the knowledge of targets of α4/6-family CTxs.

Comparative Venomics Reveals the Complex Prey Capture Strategy of the Piscivorous Cone Snail Conus catus

Venomous marine cone snails produce a unique and remarkably diverse range of venom peptides (conotoxins and conopeptides) that have proven to be invaluable as pharmacological probes and leads to new therapies. Conus catus is a hook-and-line fish hunter from clade I, with ∼20 conotoxins identified, including the analgesic ω-conotoxin CVID (AM336). The current study unravels the venom composition of C. catus with tandem mass spectrometry and 454 sequencing data. From the venom gland transcriptome, 104 precursors were recovered from 11 superfamilies, with superfamily A (especially κA-) conotoxins dominating (77%) their venom. Proteomic analysis confirmed that κA-conotoxins dominated the predation-evoked milked venom of each of six C. catus analyzed and revealed remarkable intraspecific variation in both the intensity and type of conotoxins. High-throughput FLIPR assays revealed that the predation-evoked venom contained a range of conotoxins targeting the nAChR, Cav, and Nav ion channels, consistent with α- and ω-conotoxins being used for predation by C. catus. However, the κA-conotoxins did not act at these targets but induced potent and rapid immobilization followed by bursts of activity and finally paralysis when injected intramuscularly in zebrafish. Our venomics approach revealed the complexity of the envenomation strategy used by C. catus, which contains a mix of both excitatory and inhibitory venom peptides.

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.