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

Peptide-based non-viral gene delivery: A comprehensive review of the advances and challenges

Gene therapy is the most effective treatment option for diseases, but its effectiveness is affected by the choice and design of gene carriers. The genes themselves have to pass through multiple barriers in order to enter the cell and therefore require additional vectors to carry them inside the cell. In gene therapy, peptides have unique properties and potential as gene carriers, which can effectively deliver genes into specific cells or tissues, protect genes from degradation, improve gene transfection efficiency, and enhance gene targeting and biological responsiveness. This paper reviews the research progress of peptides and their derivatives in the field of gene delivery recently, describes the obstacles encountered by foreign materials to enter the interior of the cell, and introduces the following classes of functional peptides that can carry materials into the interior of the cell, and assist in transmembrane translocation of carriers, thus breaking through endosomal traps to enable successful entry of genetic materials into the nucleus of the cell. The paper also discusses the combined application of peptide vectors with other vectors to enhance its transfection ability, explores current challenges encountered by peptide vectors, and looks forward to future developments in the field.

Targeting Dysregulated Ion Channels in Liver Tumors with Venom Peptides

The regulation of cellular processes by ion channels has become central to the study of cancer mechanisms. Designing molecules that can modify ion channels specific to tumor cells is a promising area of targeted drug delivery and therapy. Despite their potential in drug discovery, venom peptides-a group of natural products-have largely remained understudied and under-characterized. In general, venom peptides display high specificity and selectivity for their target ion channels. Therefore, they may represent an effective strategy for selectively targeting the dysregulation of ion channels in tumor cells. This review examines existing venom peptide therapies for different cancer types and focuses on the application of snail venom peptides in hepatocellular carcinoma (HCC), the most common form of primary liver cancer worldwide. We provide insights into the mode of action of venom peptides that have been shown to target tumors. We also explore the benefit of using new computational methods like de novo protein structure prediction to screen venom peptides and identify potential druggable candidates. Finally, we summarize the role of cell culture, animal, and organoid models in developing effective therapies against HCC and highlight the need for creating models that represent the most disproportionately affected ethnicities in HCC.

Virtual Screening of Peptide Libraries: The Search for Peptide-Based Therapeutics Using Computational Tools

Over the last few decades, we have witnessed growing interest from both academic and industrial laboratories in peptides as possible therapeutics. Bioactive peptides have a high potential to treat various diseases with specificity and biological safety. Compared to small molecules, peptides represent better candidates as inhibitors (or general modulators) of key protein-protein interactions. In fact, undruggable proteins containing large and smooth surfaces can be more easily targeted with the conformational plasticity of peptides. The discovery of bioactive peptides, working against disease-relevant protein targets, generally requires the high-throughput screening of large libraries, and in silico approaches are highly exploited for their low-cost incidence and efficiency. The present review reports on the potential challenges linked to the employment of peptides as therapeutics and describes computational approaches, mainly structure-based virtual screening (SBVS), to support the identification of novel peptides for therapeutic implementations. Cutting-edge SBVS strategies are reviewed along with examples of applications focused on diverse classes of bioactive peptides (i.e., anticancer, antimicrobial/antiviral peptides, peptides blocking amyloid fiber formation).

Computational Design of α-Conotoxins to Target Specific Nicotinic Acetylcholine Receptor Subtypes

Nicotinic acetylcholine receptors (nAChRs) are drug targets for neurological diseases and disorders, but selective targeting of the large number of nAChR subtypes is challenging. Marine cone snail α-conotoxins are potent blockers of nAChRs and some have been engineered to achieve subtype selectivity. This engineering effort would benefit from rapid computational methods able to predict mutational energies, but current approaches typically require high-resolution experimental structures, which are not widely available for α-conotoxin complexes. Herein, five mutational energy prediction methods were benchmarked using crystallographic and mutational data on two acetylcholine binding protein/α-conotoxin systems. Molecular models were developed for six nAChR subtypes in complex with five α-conotoxins that were studied through 150 substitutions. The best method was a combination of FoldX and molecular dynamics simulations, resulting in a predictive Matthews Correlation Coefficient (MCC) of 0.68 (85 % accuracy). Novel α-conotoxin mutants designed using this method were successfully validated by experimental assay with improved pharmaceutical properties. This work paves the way for the rapid design of subtype-specific nAChR ligands and potentially accelerated drug development.

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.

A practical guide to machine-learning scoring for structure-based virtual screening

Structure-based virtual screening (SBVS) via docking has been used to discover active molecules for a range of therapeutic targets. Chemical and protein data sets that contain integrated bioactivity information have increased both in number and in size. Artificial intelligence and, more concretely, its machine-learning (ML) branch, including deep learning, have effectively exploited these data sets to build scoring functions (SFs) for SBVS against targets with an atomic-resolution 3D model (e.g., generated by X-ray crystallography or predicted by AlphaFold2). Often outperforming their generic and non-ML counterparts, target-specific ML-based SFs represent the state of the art for SBVS. Here, we present a comprehensive and user-friendly protocol to build and rigorously evaluate these new SFs for SBVS. This protocol is organized into four sections: (i) using a public benchmark of a given target to evaluate an existing generic SF; (ii) preparing experimental data for a target from public repositories; (iii) partitioning data into a training set and a test set for subsequent target-specific ML modeling; and (iv) generating and evaluating target-specific ML SFs by using the prepared training-test partitions. All necessary code and input/output data related to three example targets (acetylcholinesterase, HMG-CoA reductase, and peroxisome proliferator-activated receptor-α) are available at https://github.com/vktrannguyen/MLSF-protocol , can be run by using a single computer within 1 week and make use of easily accessible software/programs (e.g., Smina, CNN-Score, RF-Score-VS and DeepCoy) and web resources. Our aim is to provide practical guidance on how to augment training data to enhance SBVS performance, how to identify the most suitable supervised learning algorithm for a data set, and how to build an SF with the highest likelihood of discovering target-active molecules within a given compound library.

High-Throughput Prediction and Design of Novel Conopeptides for Biomedical Research and Development

Cone snail venoms have been considered a valuable treasure for international scientists and businessmen, mainly due to their pharmacological applications in development of marine drugs for treatment of various human diseases. To date, around 800 Conus species are recorded, and each of them produces over 1,000 venom peptides (termed as conopeptides or conotoxins). This reflects the high diversity and complexity of cone snails, although most of their venoms are still uncharacterized. Advanced multiomics (such as genomics, transcriptomics, and proteomics) approaches have been recently developed to mine diverse Conus venom samples, with the main aim to predict and identify potentially interesting conopeptides in an efficient way. Some bioinformatics techniques have been applied to predict and design novel conopeptide sequences, related targets, and their binding modes. This review provides an overview of current knowledge on the high diversity of conopeptides and multiomics advances in high-throughput prediction of novel conopeptide sequences, as well as molecular modeling and design of potential drugs based on the predicted or validated interactions between these toxins and their molecular targets.

High throughput virtual screening (HTVS) of peptide library: Technological advancement in ligand discovery

High-throughput virtual screening (HTVS) is a leading biopharmaceutical technology that employs computational algorithms to uncover biologically active compounds from large-scale collections of chemical compound libraries. In addition, this method often leverages the precedence of screening focused libraries for assessing their binding affinities and improving physicochemical properties. Usually, developing a drug sometimes takes ages, and lessons are learnt from FDA-approved drugs. This screening strategy saves resources and time compared to laboratory testing in certain stages of drug discovery. Yet in-silico investigations remain challenging in some cases of drug discovery. For the last few decades, peptide-based drug discoveries have received remarkable momentum for several advantages over small molecules. Therefore, developing a high-fidelity HTVS platform for chemically versatile peptide libraries is highly desired. This review summarises the modern and frequently appreciated HTVS strategies for peptide libraries from 2011 to 2021. In addition, we focus on the software used for preparing peptide libraries, their screening techniques and shortcomings. An index of various HTVS methods reported here should assist researchers in identifying tools that could be beneficial for their peptide library screening projects.

PD1 CD44 antiviral peptide as an inhibitor of the protein-protein interaction in dengue virus invasion

Dengue virus (DENV) infection is mediated by the interaction between the virus envelope protein and cellular receptors of the host cells. In this study, we designed peptides to inhibit protein-protein interaction between dengue virus and CD44 receptor, which is one of the receptors used by DENV for entry. In silico model complexes were designed between domain III of the viral envelope protein of dengue virus 2 and the domain of human CD44 receptor using ClusPro 2.0, (https://cluspro.bu.edu/login.php), and inhibition peptides were designed with Rosetta Online-Server(http://rosie.rosettacommons.org/peptiderive). We identified one linear antiviral peptide of 18 amino acids derived from the human CD44 receptor, PD1 CD44. It did not show hemolysis or toxicity in HepG2 or BHK cell lines, nor did it stimulate the release of IL-1β, IL-6, TNF-α, and IFN-γ, below 100 µM. It had an IC50 of 13.8 µM and maximum effective dose of 54.9 µM evaluated in BHK cells. The decrease in plaque-forming units/mL for DENV1, DENV2, DENV3, and DENV4 was 99.60%, 99.40%, 97.80%, and 70.50%, respectively, and similar results were obtained by RT-qPCR. Non-structural protein 1 release was decreased in pre- and co-treatment but not in post-treatment. Competition assays between the DN59 peptide, envelope protein, and the fragment of domain III "MDKLQLKGMSYSMCTGKF" of the viral envelope of DENV2 and PD1 CD44 showed that our peptide lost its antiviral activity. We demonstrated that our peptide decreased endosome formation, and we propose that it binds to the envelope protein of DENV, inhibiting viral invasion/fusion.

An Overview of Antiviral Peptides and Rational Biodesign Considerations

Viral diseases have contributed significantly to worldwide morbidity and mortality throughout history. Despite the existence of therapeutic treatments for many viral infections, antiviral resistance and the threat posed by novel viruses highlight the need for an increased number of effective therapeutics. In addition to small molecule drugs and biologics, antimicrobial peptides (AMPs) represent an emerging class of potential antiviral therapeutics. While AMPs have traditionally been regarded in the context of their antibacterial activities, many AMPs are now known to be antiviral. These antiviral peptides (AVPs) have been shown to target and perturb viral membrane envelopes and inhibit various stages of the viral life cycle, from preattachment inhibition through viral release from infected host cells. Rational design of AMPs has also proven effective in identifying highly active and specific peptides and can aid in the discovery of lead peptides with high therapeutic selectivity. In this review, we highlight AVPs with strong antiviral activity largely curated from a publicly available AMP database. We then compile the sequences present in our AVP database to generate structural predictions of generic AVP motifs. Finally, we cover the rational design approaches available for AVPs taking into account approaches currently used for the rational design of AMPs.

sVmKTx, a transcriptome analysis-based synthetic peptide analogue of Vm24, inhibits Kv1.3 channels of human T cells with improved selectivity

Kv1.3 K⁺ channels play a central role in the regulation of T cell activation and Ca²⁺ signaling under physiological and pathophysiological conditions. Peptide toxins targeting Kv1.3 have a significant therapeutic potential in the treatment of autoimmune diseases; thus, the discovery of new toxins is highly motivated. Based on the transcriptome analysis of the venom gland of V. mexicanus smithi a novel synthetic peptide, sVmKTx was generated, containing 36 amino acid residues. sVmKTx shows high sequence similarity to Vm24, a previously characterized peptide from the same species, but contains a Glu at position 32 as opposed to Lys32 in Vm24. Vm24 inhibits Kv1.3 with high affinity (K_d = 2.9 pM). However, it has limited selectivity (~1,500-fold) for Kv1.3 over hKv1.2, hKCa3.1, and mKv1.1. sVmKTx displays reduced Kv1.3 affinity (K_d = 770 pM) but increased selectivity for Kv1.3 over hKv1.2 (~9,000-fold) as compared to Vm24, other channels tested in the panel (hKCa3.1, hKv1.1, hKv1.4, hKv1.5, rKv2.1, hKv11.1, hKCa1.1, hNav1.5) were practically insensitive to the toxin at 2.5 μM. Molecular dynamics simulations showed that introduction of a Glu instead of Lys at position 32 led to a decreased structural fluctuation of the N-terminal segment of sVmKTx, which may explain its increased selectivity for Kv1.3. sVmKTx at 100 nM concentration decreased the expression level of the Ca²⁺ -dependent T cell activation marker, CD40 ligand. The high affinity block of Kv1.3 and increased selectivity over the natural peptide makes sVmKTx a potential candidate for Kv1.3 blockade-mediated treatment of autoimmune diseases.

Marine Origin Ligands of Nicotinic Receptors: Low Molecular Compounds, Peptides and Proteins for Fundamental Research and Practical Applications

The purpose of our review is to briefly show what different compounds of marine origin, from low molecular weight ones to peptides and proteins, offer for understanding the structure and mechanism of action of nicotinic acetylcholine receptors (nAChRs) and for finding novel drugs to combat the diseases where nAChRs may be involved. The importance of the mentioned classes of ligands has changed with time; a protein from the marine snake venom was the first excellent tool to characterize the muscle-type nAChRs from the electric ray, while at present, muscle and α7 receptors are labeled with the radioactive or fluorescent derivatives prepared from α-bungarotoxin isolated from the many-banded krait. The most sophisticated instruments to distinguish muscle from neuronal nAChRs, and especially distinct subtypes within the latter, are α-conotoxins. Such information is crucial for fundamental studies on the nAChR revealing the properties of their orthosteric and allosteric binding sites and mechanisms of the channel opening and closure. Similar data are provided by low-molecular weight compounds of marine origin, but here the main purpose is drug design. In our review we tried to show what has been obtained in the last decade when the listed classes of compounds were used in the nAChR research, applying computer modeling, synthetic analogues and receptor mutants, X-ray and electron-microscopy analyses of complexes with the nAChRs, and their models which are acetylcholine-binding proteins and heterologously-expressed ligand-binding domains.

Conformational selection of vasopressin upon V1a receptor binding

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V_1aR, V_1bR and V₂R) are of high interest in a wide array of drug discovery programs. V_1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V_1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V_1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V_1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V_1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V_1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Interactions of Globular and Ribbon [γ4E]GID with α4β2 Neuronal Nicotinic Acetylcholine Receptor

The α4β2 nAChR is implicated in a range of diseases and disorders including nicotine addiction, epilepsy and Parkinson's and Alzheimer's diseases. Designing α4β2 nAChR selective inhibitors could help define the role of the α4β2 nAChR in such disease states. In this study, we aimed to modify globular and ribbon α-conotoxin GID to selectively target the α4β2 nAChR through competitive inhibition of the α4(+)β2(-) or α4(+)α4(-) interfaces. The binding modes of the globular α-conotoxin [γ4E]GID with rat α3β2, α4β2 and α7 nAChRs were deduced using computational methods and were validated using published experimental data. The binding mode of globular [γ4E]GID at α4β2 nAChR can explain the experimental mutagenesis data, suggesting that it could be used to design GID variants. The predicted mutational energy results showed that globular [γ4E]GID is optimal for binding to α4β2 nAChR and its activity could not likely be further improved through amino-acid substitutions. The binding mode of ribbon GID with the (α4)₃(β2)₂ nAChR was deduced using the information from the cryo-electron structure of (α4)₃(β2)₂ nAChR and the binding mode of ribbon AuIB. The program FoldX predicted the mutational energies of ribbon [γ4E]GID at the α4(+)α4(-) interface, and several ribbon[γ4E]GID mutants were suggested to have desirable properties to inhibit (α4)₃(β2)₂ nAChR.

The Allosteric Activation of α7 nAChR by α-Conotoxin MrIC Is Modified by Mutations at the Vestibular Site

α-conotoxins are 13–19 amino acid toxin peptides that bind various nicotinic acetylcholine receptor (nAChR) subtypes. α-conotoxin Mr1.7c (MrIC) is a 17 amino acid peptide that targets α7 nAChR. Although MrIC has no activating effect on α7 nAChR when applied by itself, it evokes a large response when co-applied with the type II positive allosteric modulator PNU-120596, which potentiates the α7 nAChR response by recovering it from a desensitized state. A lack of standalone activity, despite activation upon co-application with a positive allosteric modulator, was previously observed for molecules that bind to an extracellular domain allosteric activation (AA) site at the vestibule of the receptor. We hypothesized that MrIC may activate α7 nAChR allosterically through this site. We ran voltage-clamp electrophysiology experiments and in silico peptide docking calculations in order to gather evidence in support of α7 nAChR activation by MrIC through the AA site. The experiments with the wild-type α7 nAChR supported an allosteric mode of action, which was confirmed by the significantly increased MrIC + PNU-120596 responses of three α7 nAChR AA site mutants that were designed in silico to improve MrIC binding. Overall, our results shed light on the allosteric activation of α7 nAChR by MrIC and suggest the involvement of the AA site.

Potency- and Selectivity-Enhancing Mutations of Conotoxins for Nicotinic Acetylcholine Receptors Can Be Predicted Using Accurate Free-Energy Calculations

Nicotinic acetylcholine receptor (nAChR) subtypes are key drug targets, but it is challenging to pharmacologically differentiate between them because of their highly similar sequence identities. Furthermore, α-conotoxins (α-CTXs) are naturally selective and competitive antagonists for nAChRs and hold great potential for treating nAChR disorders. Identifying selectivity-enhancing mutations is the chief aim of most α-CTX mutagenesis studies, although doing so with traditional docking methods is difficult due to the lack of α-CTX/nAChR crystal structures. Here, we use homology modeling to predict the structures of α-CTXs bound to two nearly identical nAChR subtypes, α3β2 and α3β4, and use free-energy perturbation (FEP) to re-predict the relative potency and selectivity of α-CTX mutants at these subtypes. First, we use three available crystal structures of the nAChR homologue, acetylcholine-binding protein (AChBP), and re-predict the relative affinities of twenty point mutations made to the α-CTXs LvIA, LsIA, and GIC, with an overall root mean square error (RMSE) of 1.08 ± 0.15 kcal/mol and an R2 of 0.62, equivalent to experimental uncertainty. We then use AChBP as a template for α3β2 and α3β4 nAChR homology models bound to the α-CTX LvIA and re-predict the potencies of eleven point mutations at both subtypes, with an overall RMSE of 0.85 ± 0.08 kcal/mol and an R2 of 0.49. This is significantly better than the widely used molecular mechanics-generalized born/surface area (MM-GB/SA) method, which gives an RMSE of 1.96 ± 0.24 kcal/mol and an R2 of 0.06 on the same test set. Next, we demonstrate that FEP accurately classifies α3β2 nAChR selective LvIA mutants while MM-GB/SA does not. Finally, we use FEP to perform an exhaustive amino acid mutational scan of LvIA and predict fifty-two mutations of LvIA to have greater than 100X selectivity for the α3β2 nAChR. Our results demonstrate the FEP is well-suited to accurately predict potency- and selectivity-enhancing mutations of α-CTXs for nAChRs and to identify alternative strategies for developing selective α-CTXs.

Globular and ribbon isomers of Conus geographus α-conotoxins antagonize human nicotinic acetylcholine receptors

The short disulfide-rich α-conotoxins derived from the venom of Conus snails comprise a conserved C^IC^II(m)C^III(n)C^IV cysteine framework (m and n, number of amino acids) and the majority antagonize nicotinic acetylcholine receptors (nAChRs). Depending on disulfide connectivity, α-conotoxins can exist as either globular (C^I-C^III, C^II-C^IV), ribbon (C^I-C^IV, C^II-C^III) or bead (C^I-C^II, C^III-C^IV) isomers. In the present study, C. geographus α-conotoxins GI, GIB, G1.5 and G1.9 were chemically synthesized as globular and ribbon isomers and their activity investigated at human nAChRs expressed in Xenopus oocytes using the two-electrode voltage clamp recording technique. Both the globular and ribbon isomers of the 3/5 (m/n) α-conotoxins GI and GIB selectively inhibit heterologous human muscle-type α1β1δε nAChRs, whereas G1.5, a 4/7 α-conotoxin, selectively antagonizes neuronal (non-muscle) nAChR subtypes particularly human α3β2, α7 and α9α10 nAChRs. In contrast, globular and ribbon isomers of G1.9, a novel C-terminal elongated 4/8 α-conotoxin exhibited no activity at the human nAChR subtypes studied. This study reinforces earlier observations that 3/5 α-conotoxins selectively target the muscle nAChR subtypes, although interestingly, GIB is also active at α7 and α9 α10 nAChRs. The 4/7 α-conotoxins target human neuronal nAChR subtypes whereas the pharmacology of the 4/8 α-conotoxin remains unknown.

Virtual screening of peptides with high affinity for SARS-CoV-2 main protease

The current pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused more than 2,000,000 deaths worldwide. Currently, vaccine development and drug repurposing have been the main strategies to find a COVID-19 treatment. However, the development of new drugs could be the solution if the main strategies fail. Here, a virtual screening of pentapeptides was applied in order to identify peptides with high affinity to SARS-CoV-2 main protease (M^pro). Over 70,000 peptides were screened employing a genetic algorithm that uses a docking score as the fitness function. The algorithm was coupled with a RESTful API to persist data and avoid redundancy. The docking exhaustiveness was adapted to the number of peptides in each virtual screening step, where the higher the number of peptides, the lower the docking exhaustiveness. Two potential peptides were selected (HHYWH and HYWWT), which have higher affinity to M^pro than to human proteases. Albeit preliminary, the data presented here provide some basis for the rational design of peptide-based drugs to treat COVID-19.

The new COST Action European Venom Network (EUVEN)-synergy and future perspectives of modern venomics

Venom research is a highly multidisciplinary field that involves multiple subfields of biology, informatics, pharmacology, medicine, and other areas. These different research facets are often technologically challenging and pursued by different teams lacking connection with each other. This lack of coordination hampers the full development of venom investigation and applications. The COST Action CA19144–European Venom Network was recently launched to promote synergistic interactions among different stakeholders and foster venom research at the European level.

Novel Three-Finger Neurotoxins from Naja melanoleuca Cobra Venom Interact with GABAA and Nicotinic Acetylcholine Receptors

Cobra venoms contain three-finger toxins (TFT) including α-neurotoxins efficiently binding nicotinic acetylcholine receptors (nAChRs). As shown recently, several TFTs block GABA_A receptors (GABA_ARs) with different efficacy, an important role of the TFTs central loop in binding to these receptors being demonstrated. We supposed that the positive charge (Arg36) in this loop of α-cobratoxin may explain its high affinity to GABA_AR and here studied α-neurotoxins from African cobra N. melanoleuca venom for their ability to interact with GABAARs and nAChRs. Three α-neurotoxins, close homologues of the known N. melanoleuca long neurotoxins 1 and 2, were isolated and sequenced. Their analysis on Torpedocalifornica and α7 nAChRs, as well as on acetylcholine binding proteins and on several subtypes of GABA_ARs, showed that all toxins interacted with the GABA_AR much weaker than with the nAChR: one neurotoxin was almost as active as α-cobratoxin, while others manifested lower activity. The earlier hypothesis about the essential role of Arg36 as the determinant of high affinity to GABA_AR was not confirmed, but the results obtained suggest that the toxin loop III may contribute to the efficient interaction of some long-chain neurotoxins with GABA_AR. One of isolated toxins manifested different affinity to two binding sites on Torpedo nAChR.

System Biology-Guided Chemical Proteomics to Discover Protein Targets of Monoethylhexyl Phthalate in Regulating Cell Cycle

Chemical proteomics methods have been used as effective tools to identify novel protein targets for small molecules. These methods have great potential to be applied as environmental toxicants to figure out their mode of action. However, these assays usually generate dozens of possible targets, making it challenging to validate the most important one. In this study, we have integrated the cellular thermal shift assay (CETSA), quantitative proteomics, metabolomics, computer-assisted docking, and target validation methods to uncover the protein targets of monoethylhexyl phthalate (MEHP). Using the mass spectrometry implementation of CETSA (MS-CETSA), we have identified 74 possible protein targets of MEHP. The Gene Ontology (GO) enrichment integration was further conducted for the target proteins, the cellular dysregulated proteins, and the metabolites, showing that cell cycle dysregulation could be one primary change due to the MEHP-induced toxicity. Flow cytometry analysis confirmed that hepatocytes were arrested at the G1 stage due to the treatment with MEHP. Subsequently, the potential protein targets were ranked by their binding energy calculated from the computer-assisted docking with MEHP. In summary, we have demonstrated the development of interactomics workflow to simplify the redundant information from multiomics data and identified novel cell cycle regulatory protein targets (CPEB4, ANAPC5, and SPOUT1) for MEHP.

Mass spectrometric identification and denovo sequencing of novel conotoxins from vermivorous cone snail (Conus inscriptus), and preliminary screening of its venom for biological activities in vitro and in vivo

Venom of Conus inscriptus, a vermivorous cone snail found abundantly in the southern coastal waters was studied to yield conotoxins through proteomic analysis. A total of 37 conotoxins (4 with single disulfide bonds, 20 with two disulfide bonds and 11 three disulfide-bonded peptides) were identified using mass spectrometric analysis. Among them, amino acid sequences of 11 novel conopeptides with one, two and three disulfides belonging to different classes were derived through manual de novo sequencing. Based on the established primary sequence, they were pharmacologically classified into α conotoxins, µ conotoxins and contryphans. Except In1696 all other conopeptides have undergone C-terminal amidation. The natural venom exhibited 50% lethality at 304.82 µg/mL against zebrafish embryo and 130.31 µg/mL against brine shrimp nauplii. The anticonvulsant study of natural venom effectively reduced the locomotor activity against pentylenetetrazole (PTZ) treated zebrafish. This concludes that the venom peptides from Conus inscriptus exhibit potential anticonvulsant function, which leads to the discovery of lead molecules against seizures.

Water Thermodynamics of Peptide Toxin Binding Sites on Ion Channels

Peptide toxins isolated from venomous creatures, long prized as research tools due to their innate potency for ion channels, are emerging as drugs as well. However, it remains challenging to understand why peptide toxins bind with high potency to ion channels, to identify residues that are key for activity, and to improve their affinities via mutagenesis. We use WaterMap, a molecular dynamics simulation-based method, to gain computational insight into these three questions by calculating the locations and thermodynamic properties of water molecules in the peptide toxin binding sites of five ion channels. These include an acid-sensing ion channel, voltage-gated potassium channel, sodium channel in activated and deactivated states, transient-receptor potential channel, and a nicotinic receptor whose structures were recently determined by crystallography and cryo-electron microscopy (cryo-EM). All channels had water sites in the peptide toxin binding site, and an average of 75% of these sites were stable (low-energy), and 25% were unstable (medium or high energy). For the sodium channel, more unstable water sites were present in the deactivated state structure than the activated. Additionally, for each channel, unstable water sites coincided with the positions of peptide toxin residues that previous mutagenesis experiments had shown were important for activity. Finally, for the sodium channel in the deactivated state, unstable water sites were present in the peptide toxin binding pocket but did not overlap with the peptide toxin, suggesting that future experimental efforts could focus on targeting these sites to optimize potency.

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.

Rational Design of Functional Peptide-Gold Hybrid Nanomaterials for Molecular Interactions

Gold nanoparticles (AuNPs) have been extensively used for decades in biosensing-related development due to outstanding optical properties. Peptides, as newly realized functional biomolecules, are promising candidates of replacing antibodies, receptors, and substrates for specific molecular interactions. Both peptides and AuNPs are robust and easily synthesized at relatively low cost. Hence, peptide-AuNP-based bio-nano-technological approaches have drawn increasing interest, especially in the field of molecular targeting, cell imaging, drug delivery, and therapy. Many excellent works in these areas have been reported: demonstrating novel ideas, exploring new targets, and facilitating advanced diagnostic and therapeutic technologies. Importantly, some of them also have been employed to address real practical problems, especially in remote and less privileged areas. This contribution focuses on the application of peptide-gold hybrid nanomaterials for various molecular interactions, especially in biosensing/diagnostics and cell targeting/imaging, as well as for the development of highly active antimicrobial/antifouling coating strategies. Rationally designed peptide-gold nanomaterials with functional properties are discussed along with future challenges and opportunities.

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.

Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination

Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.

APETx-Like Peptides from the Sea Anemone Heteractis crispa, Diverse in Their Effect on ASIC1a and ASIC3 Ion Channels

Currently, five peptide modulators of acid-sensing ion channels (ASICs) attributed to structural class 1b of sea anemone toxins have been described. The APETx2 toxin is the first and most potent ASIC3 inhibitor, so its homologs from sea anemones are known as the APETx-like peptides. We have discovered that two APETx-like peptides from the sea anemone Heteractis crispa, Hcr 1b-3 and Hcr 1b-4, demonstrate different effects on rASIC1a and rASIC3 currents. While Hcr 1b-3 inhibits both investigated ASIC subtypes with IC₅₀ 4.95 ± 0.19 μM for rASIC1a and 17 ± 5.8 μM for rASIC3, Hcr 1b-4 has been found to be the first potentiator of ASIC3, simultaneously inhibiting rASIC1a at similar concentrations: EC₅₀ 1.53 ± 0.07 μM and IC₅₀ 1.25 ± 0.04 μM. The closest homologs, APETx2, Hcr 1b-1, and Hcr 1b-2, previously demonstrated the ability to inhibit hASIC3 with IC₅₀ 63 nM, 5.5, and 15.9 μM, respectively, while Hcr 1b-2 also inhibited rASIC1a with IC₅₀ 4.8 ± 0.3 μM. Computer modeling allowed us to describe the peculiarities of Hcr 1b-2 and Hcr 1b-4 interfaces with the rASIC1a channel and the stabilization of the expanded acidic pocket resulting from peptides binding which traps the rASIC1a channel in the closed state.

New Insights into the Type II Toxins from the Sea Anemone Heteractis crispa

Toxins modulating NaV channels are the most abundant and studied peptide components of sea anemone venom. Three type-II toxins, δ-SHTX-Hcr1f (= RpII), RTX-III, and RTX-VI, were isolated from the sea anemone Heteractis crispa. RTX-VI has been found to be an unusual analog of RTX-III. The electrophysiological effects of Heteractis toxins on nine NaV subtypes were investigated for the first time. Heteractis toxins mainly affect the inactivation of the mammalian NaV channels expressed in the central nervous system (NaV1.1-NaV1.3, NaV1.6) as well as insect and arachnid channels (BgNaV1, VdNaV1). The absence of Arg13 in the RTX-VI structure does not prevent toxin binding with the channel but it has changed its pharmacological profile and potency. According to computer modeling data, the δ-SHTX-Hcr1f binds within the extracellular region of the rNaV1.2 voltage-sensing domain IV and pore-forming domain I through a network of strong interactions, and an additional fixation of the toxin at the channel binding site is carried out through the phospholipid environment. Our data suggest that Heteractis toxins could be used as molecular tools for NaV channel studies or insecticides rather than as pharmacological agents.

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.

Ribbon α-Conotoxin KTM Exhibits Potent Inhibition of Nicotinic Acetylcholine Receptors

KTM is a 16 amino acid peptide with the sequence WCCSYPGCYWSSSKWC. Here, we present the nuclear magnetic resonance (NMR) structure and bioactivity of this rationally designed α-conotoxin (α-CTx) that demonstrates potent inhibition of rat α3β2-nicotinic acetylcholine receptors (rα3β2-nAChRs). Two bioassays were used to test the efficacy of KTM. First, a qualitative PC12 cell-based assay confirmed that KTM acts as a nAChR antagonist. Second, bioactivity evaluation by two-electrode voltage clamp electrophysiology was used to measure the inhibition of rα3β2-nAChRs by KTM (IC50 = 0.19 ± 0.02 nM), and inhibition of the same nAChR isoform by α-CTx MII (IC50 = 0.35 ± 0.8 nM). The three-dimensional structure of KTM was determined by NMR spectroscopy, and the final set of 20 structures derived from 32 distance restraints, four dihedral angle constraints, and two disulfide bond constraints overlapped with a mean global backbone root-mean-square deviation (RMSD) of 1.7 ± 0.5 Å. The structure of KTM did not adopt the disulfide fold of α-CTx MII for which it was designed, but instead adopted a flexible ribbon backbone and disulfide connectivity of C2-C16 and C3-C8 with an estimated 12.5% α-helical content. In contrast, α-CTx MII, which has a native fold of C2-C8 and C3-C16, has an estimated 38.1% α-helical secondary structure. KTM is the first reported instance of a Framework I (CC-C-C) α-CTx with ribbon connectivity to display sub-nanomolar inhibitory potency of rα3β2-nAChR subtypes.

Qualitative Assay to Detect Dopamine Release by Ligand Action on Nicotinic Acetylcholine Receptors

A pheochromocytoma of the rat adrenal medulla derived (a.k.a. PC12) cell-based assay for dopamine measurement by luminescence detection was customized for the qualitative evaluation of agonists and antagonists of nicotinic acetylcholine receptors (nAChRs). The assay mechanism begins with ligand binding to transmembrane nAChRs, altering ion flow into the cell and inducing dopamine release from the cell. Following release, dopamine is oxidized by monoamine oxidase generating hydrogen peroxide that catalyzes a chemiluminescence reaction involving luminol and horseradish peroxidase, thus producing a detectable response. Results are presented for the action of nAChR agonists (acetylcholine, nicotine, and cytisine), and antagonists (α-conotoxins (α-CTxs) MII, ImI, LvIA, and PeIA) that demonstrate a luminescence response correlating to the increase or decrease of dopamine release. A survey of cell growth and treatment conditions, including nerve growth factor, nicotine, ethanol, and temperature, led to optimal assay requirements to achieve maximal signal intensity and consistent response to ligand treatment. It was determined that PC12 cells treated with a combination of nerve growth factor and nicotine, and incubated at 37 °C, provided favorable results for a reduction in luminescence signal upon treatment of cells with α-CTxs. The PC12 assay is intended for use as a fast, efficient, and economic qualitative method to assess the bioactivity of molecules that act on nAChRs, in which testing of ligand-nAChR binding hypotheses and computational predictions can be validated. As a screening method for nAChR bioactivity, lead compounds can be assessed for their likelihood of exhibiting desired bioactivity prior to being subjected to more complex quantitative methods, such as electrophysiology or live animal studies.

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.

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

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.

CD97 Promotes Tumor Aggressiveness Through the Traditional G Protein-Coupled Receptor-Mediated Signaling in Hepatocellular Carcinoma

Cluster of differentiation 97 (CD97) is a member of the epidermal growth factor seven-transmembrane family belonging to the class B G protein-coupled receptors (GPCRs). The protein affects tumor aggressiveness through its cellular ligand CD55 stimulation and exhibits adhesive properties. Studies have demonstrated the involvement of CD97 in dedifferentiation, migration, invasiveness, and metastasis of tumors. However, little information is currently available on the specific role of CD97 in hepatocellular carcinoma (HCC). Here, we have shown that CD97 up-regulation in HCCs is positively correlated with tumor metastasis. Functionally, CD97 promoted cell migration and invasion in vitro. In an in vivo mouse model, overexpression of CD97 in HCC cells led to accelerated lung metastasis. Mechanistically, CD97 cooperated with the altered regulator, GPCR kinase 6 (GRK6), to mediate GPCR desensitization and internalization. Down-regulation of GRK6 suppressed CD97 internalization and promoted CD97 expression. Integrated regulatory interactions between CD97 and GRK6 stimulated downstream matrix metalloproteinase 2/9 secretion and, consequently, HCC metastasis. Conclusion: Our collective findings support the utility of CD97 as an effective potential prognosticator and therapeutic target for HCC.

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.

Protein-peptide docking: opportunities and challenges

Peptides have recently attracted much attention as promising drug candidates. Rational design of peptide-derived therapeutics usually requires structural characterization of the underlying protein-peptide interaction. Given that experimental characterization can be difficult, reliable computational tools are needed. In recent years, a variety of approaches have been developed for 'protein-peptide docking', that is, predicting the structure of the protein-peptide complex, starting from the protein structure and the peptide sequence, including variable degrees of information about the peptide binding site and/or conformation. In this review, we provide an overview of protein-peptide docking methods and outline their capabilities, limitations, and applications in structure-based drug design. Key challenges are also briefly discussed, such as modeling of large-scale conformational changes upon binding, scoring of predicted models, and optimal inclusion of varied types of experimental data and theoretical predictions into an integrative modeling process.

A tale of ligands big and small: an update on how pentameric ligand-gated ion channels interact with agonists and proteins

Pentameric ligand-gated ion channels (pLGICs, also known as Cys-loop receptors) are a large family of ion channels expressed in all Bilateria and in several groups of bacteria and archaea. They are activated by small-molecule neurotransmitters to mediate fast transmission at many central and peripheral nervous system synapses and are the target of several drugs and insecticides. Here we review recent advances in the field, focussing on new insights on the structure of the agonist-binding site and on newly discovered protein-protein interactions involving pLGICs.

Are Fireworms Venomous? Evidence for the Convergent Evolution of Toxin Homologs in Three Species of Fireworms (Annelida, Amphinomidae)

Amphinomids, more commonly known as fireworms, are a basal lineage of marine annelids characterized by the presence of defensive dorsal calcareous chaetae, which break off upon contact. It has long been hypothesized that amphinomids are venomous and use the chaetae to inject a toxic substance. However, studies investigating fireworm venom from a morphological or molecular perspective are scarce and no venom gland has been identified to date, nor any toxin characterized at the molecular level. To investigate this question, we analyzed the transcriptomes of three species of fireworms-Eurythoe complanata, Hermodice carunculata, and Paramphinome jeffreysii-following a venomics approach to identify putative venom compounds. Our venomics pipeline involved de novo transcriptome assembly, open reading frame, and signal sequence prediction, followed by three different homology search strategies: BLAST, HMMER sequence, and HMMER domain. Following this pipeline, we identified 34 clusters of orthologous genes, representing 13 known toxin classes that have been repeatedly recruited into animal venoms. Specifically, the three species share a similar toxin profile with C-type lectins, peptidases, metalloproteinases, spider toxins, and CAP proteins found among the most highly expressed toxin homologs. Despite their great diversity, the putative toxins identified are predominantly involved in three major biological processes: hemostasis, inflammatory response, and allergic reactions, all of which are commonly disrupted after fireworm stings. Although the putative fireworm toxins identified here need to be further validated, our results strongly suggest that fireworms are venomous animals that use a complex mixture of toxins for defense against predators.