Binding Modes and Functional Surface of Anti-mammalian Scorpion α-Toxins to Sodium Channels

Ts17, a Tityus serrulatus β-toxin structurally related to α-scorpion toxins

Ts17 was purified from the venom of the scorpion Tityus serrulatus, the most dangerous scorpion species in Brazil. The activity on Na_v1.1-Na_v1.7 channels was electrophysiologically characterized by patch-clamp technique. Ts17 amino acid sequence indicated high similarity to alpha-scorpion toxins; however, it presented beta-toxin activity, altering the kinetics of the Na⁺-channels. The most affected subtypes during activation (with and without prepulse) and inactivation phases were Na_v1.2 and Na_v1.5, respectively. For recovery from inactivation, the most affected voltage-gated sodium channel was Na_v1.5. Circular dichroism spectra showed that Ts17 presents mainly β-sheet and unordered structures at all analyzed pHs, and the maximum value of α-helix was found at pH 4.0 (13.3 %). Based on the results, Ts17 might be used as a template to develop a new cardiac drug. Key contribution Purification of Ts17 from Tityus serrulatus, electrophysiological characterization of Ts17 on voltage-gated sodium channel subtypes, β-toxin classification.

Methionine 58 is a key residue in the modulation of BmK scorpion toxin AGP-SYPU2 activity through in silico and in vivo studies

Protein dynamic networks play an important role in the regulation of many protein systems. Some residues that are far away from the interface between proteins and their targets have a critical role in modulating the activity of some scorpion toxins. Here, conservation analysis combined with an in vivo experiment has reveals that Met58 is a key residue of BmK scorpion toxin AGP-SYPU2 in the modulation of analgesic activity. Molecular dynamics simulations clearly reveal the conformational changes that allow the loop between the β2 and β3 sheets to be exposed on the toxin surface to interact with its targets. Our results emphasize specific roles for the residue Met58 in the NC domain and our work gives valuable information for further modification of scorpion toxins to obtain new analgesic peptides with enhanced activity. Communicated by Ramaswamy H. Sarma.

Comparison of different strategies towards the chemical synthesis of long-chain scorpion toxin AaH-II

Long-chain scorpion toxin AaH-II isolated from Androctonus australis Hector can selectively inhibit mammalian voltage-gated sodium ion channel Na_v 1.7 responsible for pain sensation. Efficient chemical synthesis of AaH-II and its derivatives is beneficial to the study of the function and mechanism of Na_v 1.7 and the development of potential peptide inhibitors. Herein, we compared three different strategies, namely, direct solid-phase peptide synthesis, hydrazide-based two-segment native chemical ligation, and hydrazide-based three-segment native chemical ligation for the synthesis of AaH-II. The hydrazide-based two-segment native chemical ligation affords the target toxin with the optimal efficiency, which provides a practically robust procedure for the preparation of tool molecules derived from AaH-II to study the biological functions and modulation of Na_v 1.7. Our work highlights the importance of selecting suitable segment condensation approach in the chemical synthesis of protein toxins.

How a Scorpion Toxin Selectively Captures a Prey Sodium Channel: The Molecular and Evolutionary Basis Uncovered

The growing resistance of insects to chemical pesticides is reducing the effectiveness of conventional methods for pest control and thus, the development of novel insecticidal agents is imperative. Scorpion toxins specific for insect voltage-gated sodium channels (Navs) have been considered as one of the most promising insecticide alternatives due to their host specificity, rapidly evoked toxicity, biodegradability, and the lack of resistance. However, they have not been developed for uses in agriculture and public health, mainly because of a limited understanding of their molecular and evolutionary basis controlling their phylogenetic selectivity. Here, we show that the traditionally defined insect-selective scorpion toxin LqhIT2 specifically captures a prey Nav through a conserved trapping apparatus comprising a three-residue-formed cavity and a structurally adjacent leucine. The former serves as a detector to recognize and bind a highly exposed channel residue conserved in insects and spiders, two major prey items for scorpions; and the latter subsequently seizes the "moving" voltage sensor via hydrophobic interactions to reduce activation energy for channel opening, demonstrating its action in an enzyme-like manner. Based on the established toxin-channel interaction model in combination with toxicity assay, we enlarged the toxic spectrum of LqhIT2 to spiders and certain other arthropods. Furthermore, we found that genetic background-dependent cavity shapes determine the species selectivity of LqhIT2-related toxins. We expect that the discovery of the trapping apparatus will improve our understanding of the evolution and design principle of Nav-targeted toxins from a diversity of arthropod predators and accelerate their uses in pest control.

Slow-, Tight-Binding Inhibition of CYP17A1 by Abiraterone Redefines Its Kinetic Selectivity and Dosing Regimen

Substantial evidence underscores the clinical efficacy of inhibiting CYP17A1-mediated androgen biosynthesis by abiraterone for treatment of prostate oncology. Previous structural analysis and in vitro assays revealed inconsistencies surrounding the nature and potency of CYP17A1 inhibition by abiraterone. Here, we establish that abiraterone is a slow-, tight-binding inhibitor of CYP17A1, with initial weak binding preceding the subsequent slow isomerization to a high-affinity CYP17A1-abiraterone complex. The in vitro inhibition constant of the final high-affinity CYP17A1-abiraterone complex ( ( K i * = 0.39 nM )yielded a binding free energy of -12.8 kcal/mol that was quantitatively consistent with the in silico prediction of -14.5 kcal/mol. Prolonged suppression of dehydroepiandrosterone (DHEA) concentrations observed in VCaP cells after abiraterone washout corroborated its protracted CYP17A1 engagement. Molecular dynamics simulations illuminated potential structural determinants underlying the rapid reversible binding characterizing the two-step induced-fit model. Given the extended residence time (42 hours) of abiraterone within the CYP17A1 active site, in silico simulations demonstrated sustained target engagement even when most abiraterone has been eliminated systemically. Subsequent pharmacokinetic-pharmacodynamic (PK-PD) modeling linking time-dependent CYP17A1 occupancy to in vitro steroidogenic dynamics predicted comparable suppression of downstream DHEA-sulfate at both 1000- and 500-mg doses of abiraterone acetate. This enabled mechanistic rationalization of a clinically reported PK-PD disconnect, in which equipotent reduction of downstream plasma DHEA-sulfate levels was achieved despite a lower systemic exposure of abiraterone. Our novel findings provide the impetus for re-evaluating the current dosing paradigm of abiraterone with the aim of preserving PD efficacy while mitigating its dose-dependent adverse effects and financial burden. SIGNIFICANCE STATEMENT: With the advent of novel molecularly targeted anticancer modalities, it is becoming increasingly evident that optimal dose selection must necessarily be predicated on mechanistic characterization of the relationships between target exposure, drug-target interactions, and pharmacodynamic endpoints. Nevertheless, efficacy has always been perceived as being exclusively synonymous with affinity-based measurements of drug-target binding. This work demonstrates how elucidating the slow-, tight-binding inhibition of CYP17A1 by abiraterone via in vitro and in silico analyses was pivotal in establishing the role of kinetic selectivity in mediating time-dependent CYP17A1 engagement and eventually downstream efficacy outcomes.

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

Computational Study of Binding of μ-Conotoxin GIIIA to Bacterial Sodium Channels NaVAb and NaVRh

Structures of several voltage-gated sodium (NaV) channels from bacteria have been determined recently, but the same feat might not be achieved for the mammalian counterparts in the near future. Thus, at present, computational studies of the mammalian NaV channels have to be performed using homology models based on the bacterial crystal structures. A successful homology model for the mammalian NaV1.4 channel was recently constructed using the extensive mutation data for binding of μ-conotoxin GIIIA to NaV1.4, which was further validated through studies of binding of other μ-conotoxins and ion permeation. Understanding the similarities and differences between the bacterial and mammalian NaV channels is an important issue, and the NaV1.4-GIIIA system provides a good opportunity for such a comparison. To this end, we study the binding of GIIIA to the bacterial channels NaVAb and NaVRh. The complex structures are obtained using docking and molecular dynamics simulations, and the dissociation of GIIIA is studied through umbrella sampling simulations. The results are compared to those obtained from the NaV1.4-GIIIA system, and the differences in the binding modes arising from the changes in the selectivity filters are highlighted.

Target-Driven Evolution of Scorpion Toxins

It is long known that peptide neurotoxins derived from a diversity of venomous animals evolve by positive selection following gene duplication, yet a force that drives their adaptive evolution remains a mystery. By using maximum-likelihood models of codon substitution, we analyzed molecular adaptation in scorpion sodium channel toxins from a specific species and found ten positively selected sites, six of which are located at the core-domain of scorpion α-toxins, a region known to interact with two adjacent loops in the voltage-sensor domain (DIV) of sodium channels, as validated by our newly constructed computational model of toxin-channel complex. Despite the lack of positive selection signals in these two loops, they accumulated extensive sequence variations by relaxed purifying selection in prey and predators of scorpions. The evolutionary variability in the toxin-bound regions of sodium channels indicates that accelerated substitutions in the multigene family of scorpion toxins is a consequence of dealing with the target diversity. This work presents an example of atypical co-evolution between animal toxins and their molecular targets, in which toxins suffered from more prominent selective pressure from the channels of their competitors. Our discovery helps explain the evolutionary rationality of gene duplication of toxins in a specific venomous species.

Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics

Peptide toxins provide valuable therapeutic leads for many diseases. As they bind to their targets with high affinity, potency is usually ensured. However, toxins also bind to off-target receptors, causing potential side effects. Thus, a major challenge in generating drugs from peptide toxins is ensuring their specificity for their intended targets. Computational methods can play an important role in solving such design problems through construction of accurate models of receptor-toxin complexes and calculation of binding free energies. Here we review the computational methods used for this purpose and their application to toxins targeting ion channels. We describe ShK and HsTX1 toxins, high-affinity blockers of the voltage-gated potassium channel Kv1.3, which could be developed as therapeutic agents for autoimmune diseases.

Bioinformatics-Aided Venomics

Venomics is a modern approach that combines transcriptomics and proteomics to explore the toxin content of venoms. This review will give an overview of computational approaches that have been created to classify and consolidate venomics data, as well as algorithms that have helped discovery and analysis of toxin nucleic acid and protein sequences, toxin three-dimensional structures and toxin functions. Bioinformatics is used to tackle specific challenges associated with the identification and annotations of toxins. Recognizing toxin transcript sequences among second generation sequencing data cannot rely only on basic sequence similarity because toxins are highly divergent. Mass spectrometry sequencing of mature toxins is challenging because toxins can display a large number of post-translational modifications. Identifying the mature toxin region in toxin precursor sequences requires the prediction of the cleavage sites of proprotein convertases, most of which are unknown or not well characterized. Tracing the evolutionary relationships between toxins should consider specific mechanisms of rapid evolution as well as interactions between predatory animals and prey. Rapidly determining the activity of toxins is the main bottleneck in venomics discovery, but some recent bioinformatics and molecular modeling approaches give hope that accurate predictions of toxin specificity could be made in the near future.

General characterization of Tityus fasciolatus scorpion venom. Molecular identification of toxins and localization of linear B-cell epitopes

This communication describes the general characteristics of the venom from the Brazilian scorpion Tityus fasciolatus, which is an endemic species found in the central Brazil (States of Goiás and Minas Gerais), being responsible for sting accidents in this area. The soluble venom obtained from this scorpion is toxic to mice being the LD50 is 2.984 mg/kg (subcutaneally). SDS-PAGE of the soluble venom resulted in 10 fractions ranged in size from 6 to 10-80 kDa. Sheep were employed for anti-T. fasciolatus venom serum production. Western blotting analysis showed that most of these venom proteins are immunogenic. T. fasciolatus anti-venom revealed consistent cross-reactivity with venom antigens from Tityus serrulatus. Using known primers for T. serrulatus toxins, we have identified three toxins sequences from T. fasciolatus venom. Linear epitopes of these toxins were localized and fifty-five overlapping pentadecapeptides covering complete amino acid sequence of the three toxins were synthesized in cellulose membrane (spot-synthesis technique). The epitopes were located on the 3D structures and some important residues for structure/function were identified.

Purification, characterization, and bioactivity of a new analgesic-antitumor peptide from Chinese scorpion Buthus martensii Karsch

Scorpion venoms are complex mixtures of dozens or even hundreds of distinct proteins, many of which have diverse bioactivities. In this study, after bioassay-driven chromatographic purification, a new dual-function peptide with analgesic and antitumor activities was isolated and designated BmK AGAP-SYPU2. The first 12 amino acid residues were sequenced with Edman degradation. The cDNA was cloned by using rapid amplification of cDNA ends from the cDNA pool from scorpion glands. The amino acid sequence of BmK AGAP-SYPU2 was then deduced, and is consistent with the molecular mass measured with MALDI-TOF-MS. A preliminary pharmacological analysis revealed the following: in the dose-effect curve plotted with the mouse-twisting test, BmK AGAP-SYPU2 showed analgesic activity with an ED50 value of 1.42 mg/kg; in the time-effect curves plotted with a hot-plate procedure, BmK AGAP-SYPU2 had similar effects to those of the painkiller morphine, except for its longer duration. BmK AGAP-SYPU2 also showed antitumor activity against Ehrlich ascites tumor and S-180 fibrosarcoma models in vivo. Sequence alignment and homology modeling showed that BmK AGAP-SYPU2 is highly conserved relative to other scorpion α-toxins. However, a few different amino acids endow it with unique molecular properties, which may be responsible for its specific bioactivities. BmK AGAP-SYPU2, a new scorpion neurotoxin with dual functions, is a potential candidate drug amenable to exploitation and modification.

Characterization of a novel BmαTX47 toxin modulating sodium channels: the crucial role of expression vectors in toxin pharmacological activity

Long-chain scorpion toxins with four disulfide bridges exhibit various pharmacological features towards the different voltage-gated sodium channel subtypes. However, the toxin production still remains a huge challenge. Here, we reported the effects of different expression vectors on the pharmacological properties of a novel toxin BmαTX47 from the scorpion Buthus martensii Karsch. The recombinant BmαTX47 was obtained using the expression vector pET-14b and pET-28a, respectively. Pharmacological experiments showed that the recombinant BmαTX47 was a new α-scorpion toxin which could inhibit the fast inactivation of rNa_v1.2, mNa_v1.4 and hNa_v1.5 channels. Importantly, the different expression vectors were found to strongly affect BmαTX47 pharmacological activities while toxins were obtained by the same expression and purification procedures. When 10 µM recombinant BmαTX47 from the pET-28a vector was applied, the values of I_5ms/I_peak for rNa_v1.2, mNa_v1.4 and hNa_v1.5 channels were 44.12% ± 3.17%, 25.40% ± 4.89% and 65.34% ± 3.86%, respectively, which were better than those values of 11.33% ± 1.46%, 15.96% ± 1.87% and 5.24% ± 2.38% for rNa_v1.2, mNa_v1.4 and hNa_v1.5 channels delayed by 10 µM recombinant BmαTX47 from the pET-14b vector. The dose-response experiments further indicated the EC₅₀ values of recombinant BmαTX47 from the pET-28a vector were 7262.9 ± 755.9 nM for rNa_v1.2 channel and 1005.8 ± 118.6 nM for hNav1.5 channel, respectively. Together, these findings highlighted the important role of expression vectors in scorpion toxin pharmacological properties, which would accelerate the understanding of the structure-function relationships of scorpion toxins and promote the potential application of toxins in the near future.

Scorpion venom components that affect ion-channels function

The number and types of venom components that affect ion-channel function are reviewed. These are the most important venom components responsible for human intoxication, deserving medical attention, often requiring the use of specific anti-venoms. Special emphasis is given to peptides that recognize Na(+)-, K(+)- and Ca(++)-channels of excitable cells. Knowledge generated by direct isolation of peptides from venom and components deduced from cloned genes, whose amino acid sequences are deposited into databanks are nowadays in the order of 1.5 thousands, out of an estimate biodiversity closed to 300,000. Here the diversity of components is briefly reviewed with mention to specific references. Structural characteristic are discussed with examples taken from published work. The principal mechanisms of action of the three different types of peptides are also reviewed. Na(+)-channel specific venom components usually are modifier of the open and closing kinetic mechanisms of the ion-channels, whereas peptides affecting K(+)-channels are normally pore blocking agents. The Ryanodine Ca(++)-channel specific peptides are known for causing sub-conducting stages of the channels conductance and some were shown to be able to internalize penetrating inside the muscle cells.

Computational methods of studying the binding of toxins from venomous animals to biological ion channels: theory and applications

The discovery of new drugs that selectively block or modulate ion channels has great potential to provide new treatments for a host of conditions. One promising avenue revolves around modifying or mimicking certain naturally occurring ion channel modulator toxins. This strategy appears to offer the prospect of designing drugs that are both potent and specific. The use of computational modeling is crucial to this endeavor, as it has the potential to provide lower cost alternatives for exploring the effects of new compounds on ion channels. In addition, computational modeling can provide structural information and theoretical understanding that is not easily derivable from experimental results. In this review, we look at the theory and computational methods that are applicable to the study of ion channel modulators. The first section provides an introduction to various theoretical concepts, including force-fields and the statistical mechanics of binding. We then look at various computational techniques available to the researcher, including molecular dynamics, brownian dynamics, and molecular docking systems. The latter section of the review explores applications of these techniques, concentrating on pore blocker and gating modifier toxins of potassium and sodium channels. After first discussing the structural features of these channels, and their modes of block, we provide an in-depth review of past computational work that has been carried out. Finally, we discuss prospects for future developments in the field.

Effect of gating modifier toxins on membrane thickness: implications for toxin effect on gramicidin and mechanosensitive channels

Various gating modifier toxins partition into membranes and interfere with the gating mechanisms of biological ion channels. For example, GsMTx4 potentiates gramicidin and several bacterial mechanosensitive channels whose gating kinetics are sensitive to mechanical properties of the membrane, whereas binding of HpTx2 shifts the voltage-activity curve of the voltage-gated potassium channel Kv4.2 to the right. The detailed process by which the toxin partitions into membranes has been difficult to probe using molecular dynamics due to the limited time scale accessible. Here we develop a protocol that allows the spontaneous assembly of a polypeptide toxin into membranes in atomistic molecular dynamics simulations of tens of nanoseconds. The protocol is applied to GsMTx4 and HpTx2. Both toxins, released in water at the start of the simulation, spontaneously bind into the lipid bilayer within 50 ns, with their hydrophobic patch penetrated into the bilayer beyond the phosphate groups of the lipids. It is found that the bilayer is about 2 Å thinner upon the binding of a GsMTx4 monomer. Such a thinning effect of GsMTx4 on membranes may explain its potentiation effect on gramicidin and mechanosensitive channels.

Binding of hanatoxin to the voltage sensor of Kv2.1

Hanatoxin 1 (HaTx1) is a polypeptide toxin isolated from spider venoms. HaTx1 inhibits the voltage-gated potassium channel kv2.1 potently with nanomolar affinities. Its receptor site has been shown to contain the S3b-S4a paddle of the voltage sensor (VS). Here, the binding of HaTx1 to the VSs of human Kv2.1 in the open and resting states are examined using a molecular docking method and molecular dynamics. Molecular docking calculations predict two distinct binding modes for the VS in the resting state. In the two binding modes, the toxin binds the S3b-S4a from S2 and S3 helices, or from S1 and S4 helices. Both modes are found to be stable when embedded in a lipid bilayer. Only the mode in which the toxin binds the S3b-S4a paddle from S2 and S3 helices is consistent with mutagenesis experiments, and considered to be correct. The toxin is then docked to the VS in the open state, and the toxin-VS interactions are found to be less favorable. Computational mutagenesis calculations performed on F278R and E281K mutant VSs show that the mutations may reduce toxin binding affinity by weakening the non-bonded interactions between the toxin and the VS. Overall, our calculations reproduce a wide range of experimental data, and suggest that HaTx1 binds to the S3b-S4a paddle of Kv2.1 from S2 and S3 helices.

The role of glycine residues at the C-terminal peptide segment in antinociceptive activity: a molecular dynamics simulation

Elucidating structural determinants in the functional regions of toxins can provide useful knowledge for designing novel analgesic peptides. Glycine residues at the C-terminal region of the neurotoxin BmK AGP-SYPU2 from the scorpion Buthus martensii Karsch (BmK) have been shown to be crucial to its analgesic activity. However, there has been no research on the structure-function relationship between the C-terminal segment of this toxin and its analgesic activity. To address this issue, we performed three MD simulations: one on the native structure and the other two on mutants of that structure. Results of these calculations suggest that the existence of glycine residues at the C-terminal segment stabilizes the protruding topology of the NC domain, which is considered an important determinant of the analgesic activity of BmK AGP-SYPU2.