Unveiling hidden toxin diversity: Discovery of novel venom components through manual curation of highly expressed sequences annotated as "no hits" in Phoneutria nigriventer spider venom gland transcriptome.

Spider venoms have evolved over thousands of years, optimizing feeding and defense mechanisms. Venom components show pharmacological and biotechnological potential, rising interest in their study. However, the isolation of spider toxins for experimental evaluation poses significant challenges. To address this, transcriptomic analysis combined with computational tools has emerged as an appealing approach to characterizing spider venoms. However, many sequences remain unidentified after automatic annotation. In this study, we manually curated a subset of previously unannotated sequences from the Phoneutria nigriventer transcriptome and identified new putative venom components. Our manual analysis revealed 29 % of the analyzed sequences were potential venom components, 29 % hypothetical/uncharacterized proteins, and 17 % cellular function proteins. Only 25 % of the originally unannotated dataset remained without any identification. Most reclassified components were cysteine-rich peptides, including 23 novel putative toxins. We also found glycine-rich peptides (GRP), corroborating the previous description of GRPs in Phoneutria pertyi venom glands. Furthermore, to emphasize the recurrence of the lack of annotation in spider venom glands transcripts, we provide a survey of the percentage of unidentified sequences in several published spider venom transcriptomics studies. In conclusion, our study highlights the importance of manual curation in uncovering novel venom components and underscores the need for improved annotation strategies to fully exploit the medical and biotechnological potential of spider venoms.

Venom-gland transcriptomics and venom proteomics of the Tibellus oblongus spider.

The Tibellus oblongus spider is an active hunter that does not spin webs and remains highly underinvestigated in terms of the venom composition. Here, we describe venom glands transcriptome and venom proteome analysis for unveiling the polypeptide composition of Tibellus oblongus spider venom. The resulting EST database includes 1733 records, including 1263 nucleotide sequences with ORFs, of these 942 have been identified as toxin-coding. The database of peptide sequences was built based on of the transcriptomics results. It contains 217 new toxins, 212 of them were detected in the T. oblongus venom by the proteomics.

Transcriptome sequencing of wolf spider Lycosa sp. (Araneae: Lycosidae) venom glands provides insights into the evolution and diversity of disulfide-rich toxins.

Wolf spiders in the genus Lycosa are important pest predators in agroforestry ecosystems, capable of feeding on a wide range of pests through the use of complex venom which can to quickly immobilize and kill prey. Because of these characteristics the toxins in wolf spiders venom may prove to be natural sources for novel drug development and biopesticides. To better understand the toxins in Lycosa venom we sequenced the transcriptome from venom glands from an undescribed species of Lycosa and comparatively analyzed the data using known protein motifs. A series of 19 disulfide-rich peptide (DRP) toxin sequences were identified and categorized into seven groups based on the number and arrangement of cysteine residues. Notably, we identified three peptide sequences with low identity to any known toxin, which may be toxin peptides specific to this species of Lycosa. In addition, to further understand the evolutionary relationships of disulfide-rich peptide toxins in spider venom, we constructed phylogenetic trees of DRP toxins from three spiders species and found that the Lycosa sp. DRPs are comparatively diverse with previous research results. This study reveals the toxin diversity of wolf spiders (Lycosa sp.) at the transcriptomic level and provides initial insights into the evolution of DRP toxins in spiders, enriching our knowledge of toxin diversity and providing new compounds for functional studies.

Transcriptome analysis reveals the peptide toxins diversity of Macrothele palpator venom.

Spider venom is a large pharmacological repertoire of different bioactive peptide toxins. However, obtaining crude venom from some spiders is challenging. Thus, studying individual toxins through venom purification is a daunting task. In this study, we constructed the cDNA library and transcriptomic sequencing from the Macrothele palpator venom glands. Subsequently, 718 high-quality expressed sequence tags (ESTs) were identified, and grouped into three categories, including 449 toxin-like (62.53 %), 136 cellular component (18.94 %) and 133 non-matched (18.52 %) based on the gene function annotation. Additionally, 112 non-redundant toxin-like peptides were classified into 13 families (families A-M) based on their sequence homology and cysteine framework. Bioinformatics analysis revealed a high sequence similarity between families A-J and the toxins from Macrothele gigas in the NR database. In contrast, families K-M had a generally low sequence homology with known spider peptide toxins and unpredictable biological functions. Taken together, this study adds many new members to the spider toxin superfamily and provides a basis for identifying various potential biological tools in M. palpator venom.

Heterologous Expression of an Insecticidal Peptide Obtained from the Transcriptome of the Colombian Spider Phoneutria depilate.

Spider venoms are composed, among other substances, of peptide toxins whose selectivity for certain physiological targets has made them powerful tools for applications such as bioinsecticides, analgesics, antiarrhythmics, antibacterials, antifungals and antimalarials, among others. Bioinsecticides are an environmentally friendly alternative to conventional agrochemicals. In this paper, the primary structure of an insecticidal peptide was obtained from the venom gland transcriptome of the ctenid spider Phoneutria depilata (Transcript ID PhdNtxNav24). The peptide contains 53 amino acids, including 10 Cys residues that form 5 disulfide bonds. Using the amino acid sequence of such peptide, a synthetic gene was constructed de novo by overlapping PCRs and cloned into an expression vector. A recombinant peptide, named delta-ctenitoxin (rCtx-4), was obtained. It was expressed, folded, purified and validated using mass spectrometry (7994.61 Da). The insecticidal activity of rCtx-4 was demonstrated through intrathoracic injection in crickets (LD50 1.2 µg/g insect) and it was not toxic to mice. rCtx-4 is a potential bioinsecticide that could have a broad spectrum of applications in agriculture.

Transcriptome-based analyses reveal venom diversity in two araneomorph spiders, Psechrus triangulus and Hippasa lycosina.

The spiders Psechrus triangulus and Hippasa lycosina are widely distributed in Yunnan Province, China, and are important natural enemies of agricultural pests, yet studies regarding the composition of their venom are lacking. In this study, cDNA libraries were constructed from venom gland tissue of P. triangulus and H. lycosina and used for transcriptomic analysis. From the analysis, 39 and 31 toxin-like sequences were predicted for P. triangulus and H. lycosina, respectively. The predicted neurotoxin sequences were categorized according to cysteine sequence motifs, and the predicted neurotoxin sequences of P. triangulus and H. lycosina could be classified into 9 and 6 toxin families, respectively. In addition, potential acetylcholinesterase, hyaluronidase, and astaxanthin-like metalloproteinases were identified through annotation. In summary, transcriptomic techniques were invaluable in mining the gene expression information from these two spider species to explore the toxin composition of their venom and determine how they differ. Studies of this type provide essential baseline data for studying the evolution and physiological activities of spider toxins and for the potential development of medicinal compounds.

The deep-rooted origin of disulfide-rich spider venom toxins.

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

The majority of spiders rely on their venom to defend themselves, to hunt, or both. Armed with this formidable weapon, they have managed to conquer every continent besides Antarctica since they first emerged about 495 million years ago. A closer look at spider venoms hints at an intriguing evolutionary history which has been rarely examined so far. The venom of other animals, such as snakes or scorpions, is usually formed of a wide range of unrelated toxins; in contrast, spiders rely on a single class of proteins, known as disulfide-rich peptides, to create their deadly venom cocktail. This family of molecules is impressively diverse, with each peptide having a distinct structure and mode of action. Its origins, however, have remained elusive. To fill this knowledge gap, Shaikh and Sunagar scanned the sequences of all disulfide-rich peptides generated to date, bringing together a dataset that includes 60% of all modern-day spiders. The analyses allowed the identification of 78 new superfamilies of spider toxins. They also revealed that all existing peptides originate from a single molecule, which Shaikh and Sunagar named after the powerful Hindu goddess Adi Shakti. This ancestral toxin was present 375 million years ago in the last common ancestor of modern-day spiders. The work also highlighted that disulfide-rich peptides evolved under different pressures in various groups of spiders; this may be because some species primarily use their venom for hunting, and others for defence. While the 'hunters' may need to constantly acquire toxins with new roles and structures to keep their edge over their prey, those that rely on venom to protect themselves may instead benefit from relying on tried-and-tested toxins useful against a range of infrequent predators. Finally, the analyses revealed that the disulphide-rich peptides of Mygalomorphae tarantulas, which form one of the three major groups of spiders, are much more diverse than the related toxins in other spiders. The underlying reason for this difference is still unclear. Several life-saving drugs currently on the market are based on toxins first identified in the venoms of snakes, cone sails or lizards. Similar discoveries could be unlocked by better understanding the range of deadly molecules used by spiders, and how these came to be.

A secretory system for extracellular production of spider neurotoxin huwentoxin-I in Escherichia coli.

Due to their advantages in structural stability and versatility, cysteine-rich peptides, which are secreted from the venom glands of venomous animals, constitute a naturally occurring pharmaceutical arsenal. However, the correct folding of disulfide bonds is a challenging task in the prokaryotic expression system like Escherichia coli due to the reducing environment. Here, a secretory expression plasmid pSE-G1M5-SUMO-HWTX-I for the spider neurotoxin huwentoxin-I (HWTX-I) with three disulfides as a model of cysteine-rich peptides was constructed. By utilizing the signal peptide G1M5, the fusion protein 6 × His-SUMO-HWTX-I was successfully secreted into extracellular medium of BL21(DE3). After enrichment using cation-exchange chromatography and purification utilizing the Ni-NTA column, 6 × His-SUMO-HWTX-I was digested via Ulp1 kinase to release recombinant HWTX-I (rHWTX-I), which was further purified utilizing RP-HPLC. Finally, both impurities with low and high molecular weights were completely removed. The molecular mass of rHWTX-I was identified as being 3750.8 Da, which was identical to natural HWTX-I with three disulfide bridges. Furthermore, by utilizing whole-cell patch clamp, the sodium currents of hNav1.7 could be inhibited by rHWTX-I and the IC50 value was 419 nmol/L.

Molecular Diversity of Linear Peptides Revealed by Transcriptomic Analysis of the Venom Gland of the Spider Lycosa poonaensis.

Spider venom is a complex mixture of bioactive components. Previously, we identified two linear peptides in Lycosa poonaensis venom using mass spectrometric analysis and predicted the presence of more linear peptides therein. In this study, a transcriptomic analysis of the L. poonaensis venom gland was conducted to identify other undetermined linear peptides in the venom. The results identified 87 contigs encoding peptides and proteins in the venom that were similar to those in other spider venoms. The number of contigs identified as neurotoxins was the highest, and 15 contigs encoding 17 linear peptide sequences were identified. Seven peptides that were representative of each family were chemically synthesized, and their biological activities were evaluated. All peptides showed significant antibacterial activity against Gram-positive and Gram-negative bacteria, although their selectivity for bacterial species differed. All peptides also exhibited paralytic activity against crickets, but none showed hemolytic activity. The secondary structure analysis based on the circular dichroism spectroscopy showed that all these peptides adopt an amphiphilic α-helical structure. Their activities appear to depend on the net charge, the arrangement of basic and acidic residues, and the hydrophobicity of the peptides.

Multiomics Profiling of Toxins in the Venom of the Amazonian Spider Acanthoscurria juruenicola.

Acanthoscurria juruenicola is an Amazonian spider described for the first time almost a century ago. However, little is known about their venom composition. Here, we present a multiomics characterization of A. juruenicola venom by a combination of transcriptomics, proteomics, and peptidomics approaches. Transcriptomics of female venom glands resulted in 93,979 unique assembled mRNA transcript encoding proteins. A total of 92 proteins were identified in the venom by mass spectrometry, including 14 mature cysteine-rich peptides (CRPs). Quantitative analysis showed that CRPs, cysteine-rich secretory proteins, metalloproteases, carbonic anhydrases, and hyaluronidase comprise >90% of the venom proteome. Relative quantification of venom toxins was performed by DIA and DDA, revealing converging profiles of female and male specimens by both methods. Biochemical assays confirmed the presence of active hyaluronidases, phospholipases, and proteases in the venom. Moreover, the venom promoted in vivo paralytic activities in crickets, consistent with the high concentration of CRPs. Overall, we report a comprehensive analysis of the arsenal of toxins of A. juruenicola and highlight their potential biotechnological and pharmacological applications. Mass spectrometry data were deposited to the ProteomeXchange Consortium via the PRIDE repository with the dataset identifier PXD013149 and via the MassIVE repository with the dataset identifier MSV000087777.

General mechanism of spider toxin family I acting on sodium channel Nav1.7.

Various peptide toxins in animal venom inhibit voltage-gated sodium ion channel Nav1.7, including Nav-targeting spider toxin (NaSpTx) Family I. Toxins in NaSpTx Family I share a similar structure, i.e., N-terminal, loops 1-4, and C-terminal. Here, we used Mu-theraphotoxin-Ca2a (Ca2a), a peptide isolated from Cyriopagopus albostriatus, as a template to investigate the general properties of toxins in NaSpTx Family I. The toxins interacted with the cell membrane prior to binding to Nav1.7 via similar hydrophobic residues. Residues in loop 1, loop 4, and the C-terminal primarily interacted with the S3-S4 linker of domain II, especially basic amino acids binding to E818. We also identified the critical role of loop 2 in Ca2a regarding its affinity to Nav1.7. Our results provide further evidence that NaSpTx Family I toxins share similar structures and mechanisms of binding to Nav1.7.

Analysis of NGS data from Peruvian Loxosceles laeta spider venom gland reveals toxin diversity.

Accidents involving spiders from the genus Loxosceles cause medical emergencies in several countries of South America. The species Loxosceles laeta is ubiquitously present in Peru and is responsible for severe accidents in this country. To further characterize L. laeta venom components and to unveil possible variations in the Peruvian population, we provide an overview of the toxins-related transcripts present in the venom gland of Peruvian L. laeta. A dataset from a cDNA library previously sequenced by MiSeq sequencer (Illumina) was re-analyzed and the obtained data was compared with available sequences from Loxosceles toxins. Phospholipase-D represent the majority (69,28 %) of the transcripts related to venom toxins, followed by metalloproteases (20,72 %), sicaritoxins (6,03 %), serine-proteases (2,28 %), hyaluronidases (1,80 %) and Translationally Controlled Tumor Protein (TCTP) (0,56 %). New sequences of phospholipases D,sicaritoxins, hyaluronidase, TCTP and serine proteinases were described. Differences between the here-described toxin sequences and others, previously identified in venom glands from other spiders, were visualized upon sequence alignments. In addition, an in vitro hyaluronidase activity assay was also performed to complement comparisons between Peruvian and Brazilian L. laeta venom enzymatic activities, revealing a superior activity in the venom from Brazilian specimens. These new data provide a molecular basis that can help to explain the difference in toxicity among L. laeta venoms from different countries in South America.

Improving the Annotation of the Venom Gland Transcriptome of Pamphobeteus verdolaga, Prospecting Novel Bioactive Peptides.

Spider venoms constitute a trove of novel peptides with biotechnological interest. Paucity of next-generation-sequencing (NGS) data generation has led to a description of less than 1% of these peptides. Increasing evidence supports the underestimation of the assembled genes a single transcriptome assembler can predict. Here, the transcriptome of the venom gland of the spider Pamphobeteus verdolaga was re-assembled, using three free access algorithms, Trinity, SOAPdenovo-Trans, and SPAdes, to obtain a more complete annotation. Assembler's performance was evaluated by contig number, N50, read representation on the assembly, and BUSCO's terms retrieval against the arthropod dataset. Out of all the assembled sequences with all software, 39.26% were common between the three assemblers, and 27.88% were uniquely assembled by Trinity, while 27.65% were uniquely assembled by SPAdes. The non-redundant merging of all three assemblies' output permitted the annotation of 9232 sequences, which was 23% more when compared to each software and 28% more when compared to the previous P. verdolaga annotation; moreover, the description of 65 novel theraphotoxins was possible. In the generation of data for non-model organisms, as well as in the search for novel peptides with biotechnological interest, it is highly recommended to employ at least two different transcriptome assemblers.

Toxin diversity revealed by de novo transcriptome assembly for venom gland in two species of spiders (Trichonephila clavata and Sinopoda pengi).

During long-term predator-prey coevolution, spiders have generated a vast diversity of toxins. Trichonephila clavata is a web-spinning spider whose large, well-constructed webs and venomous arsenal facilitate prey capture. In contrast, Sinopoda pengi is an ambush predator with agile locomotion and strong chelicerae for hunting. In this study, transcriptomic analysis was performed to describe the predicted toxins of S. pengi and T. clavata. A total of 43 and 47 of these unigenes from S. pengi and T. clavata, respectively, were predicted to have toxin activity. Putative neurotoxins were classified to the family level according to cysteine arrangement; 4 and 6 toxin families were produced by S. pengi and T. clavata, respectively. In addition, potential metalloproteases, acetylcholinesterases, serine proteases, hyaluronidases and phospholipases were found by annotation in databases. In summary, molecular templates with potential application value for medical and biological fields were obtained by classifying and characterizing presumed venom components, which established a foundation for further study of venom.

Conformational decoupling in acid-sensing ion channels uncovers mechanism and stoichiometry of PcTx1-mediated inhibition.

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels involved in fast synaptic transmission. Pharmacological inhibition of ASIC1a reduces neurotoxicity and stroke infarct volumes, with the cysteine knot toxin psalmotoxin-1 (PcTx1) being one of the most potent and selective inhibitors. PcTx1 binds at the subunit interface in the extracellular domain (ECD), but the mechanism and conformational consequences of the interaction, as well as the number of toxin molecules required for inhibition, remain unknown. Here, we use voltage-clamp fluorometry and subunit concatenation to decipher the mechanism and stoichiometry of PcTx1 inhibition of ASIC1a. Besides the known inhibitory binding mode, we propose PcTx1 to have at least two additional binding modes that are decoupled from the pore. One of these modes induces a long-lived ECD conformation that reduces the activity of an endogenous neuropeptide. This long-lived conformational state is proton-dependent and can be destabilized by a mutation that decreases PcTx1 sensitivity. Lastly, the use of concatemeric channel constructs reveals that disruption of a single PcTx1 binding site is sufficient to destabilize the toxin-induced conformation, while functional inhibition is not impaired until two or more binding sites are mutated. Together, our work provides insight into the mechanism of PcTx1 inhibition of ASICs and uncovers a prolonged conformational change with possible pharmacological implications.

A Deep Learning Approach with Data Augmentation to Predict Novel Spider Neurotoxic Peptides.

As major components of spider venoms, neurotoxic peptides exhibit structural diversity, target specificity, and have great pharmaceutical potential. Deep learning may be an alternative to the laborious and time-consuming methods for identifying these peptides. However, the major hurdle in developing a deep learning model is the limited data on neurotoxic peptides. Here, we present a peptide data augmentation method that improves the recognition of neurotoxic peptides via a convolutional neural network model. The neurotoxic peptides were augmented with the known neurotoxic peptides from UniProt database, and the models were trained using a training set with or without the generated sequences to verify the augmented data. The model trained with the augmented dataset outperformed the one with the unaugmented dataset, achieving accuracy of 0.9953, precision of 0.9922, recall of 0.9984, and F1 score of 0.9953 in simulation dataset. From the set of all RNA transcripts of Callobius koreanus spider, we discovered neurotoxic peptides via the model, resulting in 275 putative peptides of which 252 novel sequences and only 23 sequences showing homology with the known peptides by Basic Local Alignment Search Tool. Among these 275 peptides, four were selected and shown to have neuromodulatory effects on the human neuroblastoma cell line SH-SY5Y. The augmentation method presented here may be applied to the identification of other functional peptides from biological resources with insufficient data.

Analysis of High Molecular Mass Compounds from the Spider Pamphobeteus verdolaga Venom Gland. A Transcriptomic and MS ID Approach.

Nowadays, spider venom research focuses on the neurotoxic activity of small peptides. In this study, we investigated high-molecular-mass compounds that have either enzymatic activity or housekeeping functions present in either the venom gland or venom of Pamphobeteus verdolaga. We used proteomic and transcriptomic-assisted approaches to recognize the proteins sequences related to high-molecular-mass compounds present in either venom gland or venom. We report the amino acid sequences (partial or complete) of 45 high-molecular-mass compounds detected by transcriptomics showing similarity to other proteins with either enzymatic activity (i.e., phospholipases A2, kunitz-type, hyaluronidases, and sphingomyelinase D) or housekeeping functions involved in the signaling process, glucanotransferase function, and beta-N-acetylglucosaminidase activity. MS/MS analysis showed fragments exhibiting a resemblance similarity with different sequences detected by transcriptomics corresponding to sphingomyelinase D, hyaluronidase, lycotoxins, cysteine-rich secretory proteins, and kunitz-type serine protease inhibitors, among others. Additionally, we report a probably new protein sequence corresponding to the lycotoxin family detected by transcriptomics. The phylogeny analysis suggested that P. verdolaga includes a basal protein that underwent a duplication event that gave origin to the lycotoxin proteins reported for Lycosa sp. This approach allows proposing an evolutionary relationship of high-molecular-mass proteins among P. verdolaga and other spider species.

Neurotoxins in the venom gland of Calommata signata, a burrowing spider.

Calommata signata, a burrowing spider, represents a special type of predation mode in spiders, and its utilization of toxins is different from that of web-weaving spiders and wandering spiders. The existing researches on spider toxins are mainly focused on the web-weaving and wandering spiders, but little attention on that of the burrowing spiders. Through transcriptome sequencing of C. signata venom gland and the remaining part as the counterpart tissue, 25 putative neurotoxin precursors were identified. These most neurotoxins were novel because their low similarities with the known sequences except for that of over 50% similarities in four neuropeptide toxins. The 25 neuropeptide toxins were divided into five families according to the constitution of cysteines for the possible disulfide bonds and the similarities of the deduced amino acid sequences. Besides neuropeptide toxins, other potential toxins in the venom gland were also analyzed. Unlike web-weaving spiders and wandering spiders, only a few neurotoxin genes were significantly expressed in the venom gland of C. signata. In the non-peptide toxin genes, only CsTryp_SPc-1, CsPA2-1, CsVa5-2 and four PDI genes were abundantly expressed in the venom gland. The present study provided an improved understanding on the spider toxin diversity and useful information for the exploitation of spider toxins.

Description of a serpin toxin in Loxosceles (Brown spider) venoms: Cloning, expression in baculovirus-infected insect cells and functional characterization.

Several classes of toxins are present in the venom of Brown spiders (Loxosceles genus), some of them are highly expressed and others are less expressed. In this work, we aimed to clone the sequence of a little expressed novel toxin from Loxosceles venom identified as a serine protease inhibitor (serpin), as well as to express and characterize its biochemical and biological properties. It was named LSPILT, derived from Loxoscelesserine protease inhibitor-like toxin. Multiple alignment analysis revealed high identity between LSPILT and other serpin molecules from spiders and crab. LSPILT was produced in baculovirus-infected insect cells, resulting in a 46-kDa protein fused to a His-tag. Immunological assays showed epitopes in LSPILT that resemble native venom toxins of Loxosceles spiders. The inhibitory activity of LSPILT on trypsin was found both by reverse zymography and fluorescent gelatin-degradation assay. Additionally, LSPILT inhibited the complement-dependent lysis of Trypanosoma cruzi epimastigotes, reduced thrombin-dependent clotting and suppressed B16-F10 melanoma cells migration. Results described herein prove the existence of conserved serpin-like toxins in Loxosceles venoms. The availability of a recombinant serpin enabled the determination of its biological and biochemical properties and indicates potential applications in future studies regarding the pathophysiology of the envenoming or for biotechnological purposes.

Potency-Enhancing Mutations of Gating Modifier Toxins for the Voltage-Gated Sodium Channel NaV1.7 Can Be Predicted Using Accurate Free-Energy Calculations.

Gating modifier toxins (GMTs) isolated from venomous organisms such as Protoxin-II (ProTx-II) and Huwentoxin-IV (HwTx-IV) that inhibit the voltage-gated sodium channel NaV1.7 by binding to its voltage-sensing domain II (VSDII) have been extensively investigated as non-opioid analgesics. However, reliably predicting how a mutation to a GMT will affect its potency for NaV1.7 has been challenging. Here, we hypothesize that structure-based computational methods can be used to predict such changes. We employ free-energy perturbation (FEP), a physics-based simulation method for predicting the relative binding free energy (RBFE) between molecules, and the cryo electron microscopy (cryo-EM) structures of ProTx-II and HwTx-IV bound to VSDII of NaV1.7 to re-predict the relative potencies of forty-seven point mutants of these GMTs for NaV1.7. First, FEP predicted these relative potencies with an overall root mean square error (RMSE) of 1.0 ± 0.1 kcal/mol and an R2 value of 0.66, equivalent to experimental uncertainty and an improvement over the widely used molecular-mechanics/generalized born-surface area (MM-GB/SA) RBFE method that had an RMSE of 3.9 ± 0.8 kcal/mol. Second, inclusion of an explicit membrane model was needed for the GMTs to maintain stable binding poses during the FEP simulations. Third, MM-GB/SA and FEP were used to identify fifteen non-standard tryptophan mutants at ProTx-II[W24] predicted in silico to have a at least a 1 kcal/mol gain in potency. These predicted potency gains are likely due to the displacement of high-energy waters as identified by the WaterMap algorithm for calculating the positions and thermodynamic properties of water molecules in protein binding sites. Our results expand the domain of applicability of FEP and set the stage for its prospective use in biologics drug discovery programs involving GMTs and NaV1.7.