Engineering of a Spider Peptide via Conserved Structure-Function Traits Optimizes Sodium Channel Inhibition In Vitro and Anti-Nociception In Vivo

Introduction

Animal venoms are an exquisite source of bio-active peptides that modulate human neurophysiology. Spider venoms, in particular, are rich in inhibitory cysteine knot (ICK) peptides that often modulate voltage-gated ion channels in pain pathways (Cardoso et al., 2018; Cardoso and Lewis, 2018). These ICK peptides bind to the voltage-sensor domains of ion channels and modify their gating properties to open and close the ion-selective pore at non-conventional membrane potentials, and are hence named gating modifiers (Cardoso and Lewis, 2019). Amongst known spider gating-modifiers are potent modulators of voltage-gated sodium (Na_V) channels such as the ProTx family, Df1a, Tap1a, Hm1a and Pn3a (Priest et al., 2007; Schmalhofer et al., 2008; Cardoso et al., 2015; Cardoso et al., 2017; Deuis et al., 2017; Cardoso et al., 2021), with a few simultaneously modulating the low-voltage T-type calcium (Ca_V3) channels (Bladen et al., 2014; Cardoso et al., 2017; Cardoso et al., 2021) and the voltage-gated potassium (K_V) channels (Middleton et al., 2002).

The Na_V channel family comprises nine subtypes (Na_V1.1–Na_V1.9) composed of an α-subunit associated with one or more auxiliary subunits β1 to β4 (Cardoso and Lewis, 2018). These channels are widely expressed in the central and peripheral nervous system, with subtypes Na_V1.4 and Na_V1.5 expressed in skeletal and cardiac muscles. Alterations in Na_V channel function contribute to a range of neurological disorders including chronic pain, epilepsy, and motor neuron disease (Kubat Oktem et al., 2016; Catterall, 2017; Cardoso and Lewis, 2018; Saba et al., 2019; Cardoso, 2020). These alterations are presented as remodelling of the expression and excitability of the subtypes Na_V1.1, Na_V1.3, Na_V1.6, Na_V1.7, Na_V1.8 and Na_V1.9 (Cardoso and Lewis, 2018; Cardoso, 2020). Chronic visceral pain is a highly prevalent problem associated to irritable bowel syndrome (IBS) that affects 11% of the global population (Chey et al., 2015; Enck et al., 2016), and in which Na_V channel subtypes were shown to participate via signalling visceral pain in IBS (Feng et al., 2015; Osteen et al., 2016; Erickson et al., 2018; Grundy et al., 2018; Salvatierra et al., 2018). Spider peptides targeting Na_V channels have shown promising pre-clinical therapeutic effects in reverting some of these Na_V-related complex conditions. For example, the peptides ProTx-II, HwTx-IV and Tap1a showed analgesic effects in painful diabetic neuropathy, spared nerve injury-induced neuropathy and chronic visceral pain (Liu et al., 2014; Tanaka et al., 2015; Flinspach et al., 2017; Cardoso et al., 2021), respectively, and the peptide Hm1a reduced seizures in Dravet syndrome (Richards et al., 2018).

Studies of the structure-function relationships of Na_V-targeting spider peptides have guided the design of optimized leads with improved potency and selectivity for Na_V channel subtypes (Cardoso and Lewis, 2019). Key residues involved in channel inhibition and selectivity and conserved amongst spider peptides of the same NaSpTx family have been unravelled (Klint et al., 2012; Cardoso and Lewis, 2019). For example, NaSpTx1 peptides often bind to the voltage-sensor domain II (VSDII) of Na_V channels to engender ion channel inhibition, as observed for HwTx-IV (Xiao et al., 2008), CcoTx-1 (Shcherbatko et al., 2016), Df1a (Cardoso et al., 2017) and Hd1a (Klint et al., 2015). Studies using residue substitutions demonstrated the positions 5 or 6 in the loop 1 of these ICK structures are typically occupied by phenylalanine and are critical for Na_V inhibition (Cardoso and Lewis, 2019). In addition, studies of activity optimization in NaSpTx1 peptides introduced a decrease in negative charge at the N-terminus (Minassian et al., 2013; Revell et al., 2013; Rong et al., 2013; Shcherbatko et al., 2016; Rahnama et al., 2017; Mueller et al., 2020), with such introduction also described to increase the ability of these peptides to bind to the cell membrane (Henriques et al., 2016; Mueller et al., 2020). At the C-terminal, residue modifications that enhanced Na_V inhibition introduced positively charged and/or hydrophobic residues (Minassian et al., 2013; Revell et al., 2013; Shcherbatko et al., 2016). Although these optimized leads have shown enhanced in vitro pharmacology, the translation into improved benefits in vivo in pre-clinical models of pain remain to be achieved for these peptide leads (Cardoso et al., 2017; Agwa et al., 2020).

In this work, we investigated the structure-function relationships of the spider peptide Tap1a by using pharmacological studies with the VSDII of Na_V1.1 and by evaluating the impact of alanine residue substitutions on the activity of Tap1a at Na_V1.1 and Na_V1.7. We designed two optimized versions of Tap1a based on the enhanced structure-function properties of other NaSpTx1 peptides and characterized their activities at Na_V channels and their binding mode to VSDII. Finally, we investigated if the optimization of Na_V inhibition observed in vitro translated in vivo in the reversal of nocifensive behaviours in a murine model of Na_V-induced spontaneous pain. This work contributes to the understanding of the structure-function relationships and optimization of spider peptides targeting Na_V channels, and the development of useful therapeutic leads with promising enhanced in vivo benefits for the treatment of neurological disorders.

Results

Tap1a Binds to VDSII S3–S4 Loop of Na_V1.1

Spider peptides belonging to NaSpTx1 are known to bind to VSDII of Na_V channels and induce gating-modifying effects that inhibit channel activation (Xiao et al., 2011; Cai et al., 2015). Tap1a displayed similar modulatory mechanism by inhibiting ion currents of the K_V2.1/Na_V1.1 VSDII chimera (Figure 1A) with a calculated IC₅₀ value of approximately 185 nM (data not shown). Tap1a at 1 μM had weak effect on currents mediated by the wild type K_V2.1 channel and reduced potassium currents by approximately 18% (area under the curve) (Figure 1B). These results suggested that the potent nanomolar inhibitory effects of Tap1a observed at the K_V2.1/Na_V1.1 VSDII chimera were due to its binding to the Na_V1.1 VSDII S3-S4 loop. Non-transfected HEK293 cells had endogenous potassium currents measuring up to 0.3 nA that were not affected by Tap1a at up to 1 μM but were inhibited by 100 μM nifedipine (Figure 1C). Cells transfected with the chimera K_V2.1/Na_V1.1 VSDII produced potassium currents which were also sensitive to nifedipine at 100 μM (Supplementary Figure S1).

FIGURE 1

FIGURE 1. Binding site of Tap1a at hNa_V1.1 and pharmacological properties of Tap1a alanine mutants for hNa_V1.1 and hNa_V1.7 channels. (A–C) A rK_V2.1/hNa_V1.1 chimera containing the paddles S3–S4 from DII from Na_V1.1 was used to explore the binding site of Tap1a over Na_V1.1. Potassium currents were elicited by a pulse at +20 mV for 500 ms from –90 mV holding potential. (A) The K⁺ currents are shown before and after addition of increasing concentrations of Tap1a. (B) Tap1a had weak effect on wild type rK_V2.1 at 1 μM. (C) Tap1a had no effect on endogenous K⁺ currents from non-transfected HEK293 cells. Data are representatives from n = 3–5 independent experiments for each condition assayed, one whole cell was considered per independent experiment. (D, E) Tap1a and alanine mutants produced by recombinant expression were characterized via automated electrophysiology using QPatch 16X and HEK293 cells expressing hNa_V1.1 (D) or hNa_V1.7 (E) co-expressed with the β1 subunit. Sodium currents were elicited by depolarization to 0 mV for 20 ms from a pre-pulse at –120 mV for 200 ms and –80 mV holding potential. Serial diluted concentrations of each alanine mutant were incubated with the cells for 5 min at holding potential before channel activation. Data are described by mean ± SEM from n = 3–5 experiments. The calculated IC₅₀ values are described in Table 1 (F) pIC₅₀ values for Tap1a and alanine mutants tested against hNa_V1.1 and hNa_V1.7 and calculated from the dose responses produced in (D) and (E). The selectivity between hNa_V1.1 and hNa_V1.7 are described above the bars when applicable. (G) Primary sequence of Tap1a highlighting in red residues which alanine substitution deplete the inhibitory activity, in blue the residue which alanine substitution led to a decrease in hNa_V1.7 activity and maintained hNa_V1.1 inhibitory activity, and in green the residue which alanine substitution led to an increase in hNa_V1.1 activity and decrease in hNa_V1.7 inhibitory activity. Cysteines are highlighted in grey boxes. (H) Model of the three-dimensional structure of Tap1a wild type highlighting in red the residues which were essential for hNa_V1.1 and hNa_V1.7 activity, in blue the residue D27 in which alanine substitution led to an increase or maintained inhibitory activity, and in green the residue Q33 in which alanine substitution led to a mixed increase and decrease in inhibitory activity.

Structure-Function Properties of Tap1a Resemble NaSpTx1 Peptides

In order to determine the structure-function relationships of individual residues in Tap1a, we applied the structure-function knowledge for NaSpTx1 (Cardoso and Lewis, 2019) and designed twelve alanine mutants of Tap1a: M6A, F7A, P12A, D15A, S24A, D27A, Q28A, K31A, Y32A, Q33A, L34A and W35A, and produced these new mutants via recombinant expression (Supplementary Table S1; Supplementary Figure S2). The recombinant expression of these peptides was achieved successfully and produced a major single isomer purified for further testing (Supplementary Figure S2). The molecular weight of the purified recombinant peptides was verified and confirmed by mass spectrometry (Supplementary Table S1; Supplementary Figure S2).

Tap1a alanine mutants were evaluated via automated whole cell patch-clamp electrophysiology with activities determined for the subtypes Na_V1.1 and Na_V1.7 (Figures 1D–G and Table 1). Residues identified as critical for Tap1a inhibitory activity were M6, F7, P12, D15, S24, Q28, K31, Y32, L34 and W35 in which alanine substitutions led to weaker or complete loss of activity at up to 1 μM peptide tested. The substitution D27A induced partial loss of activity at Na_V1.7 (4.3-fold decrease) and did not affect activity at Na_V1.1 compared to Tap1a-WT. Interestingly, the substitution Q33A increased the Tap1a inhibitory activity for Na_V1.1 by 1.8-fold (not statistically significant) and had nearly no effect on the inhibition of Na_V1.7 compared to Tap1a-WT (1.2-fold decrease).

TABLE 1

TABLE 1. Pharmacological activity of Tap1a alanine mutants for hNa_V1.1 and hNa_V1.7 determined by automated electrophysiology in QPatch 16X. The values for the Tap1a wild type (WT) were obtained from our previous study (Cardoso et al., 2021). Tap1a alanine mutants were tested at up to 1 μM concentration.

The modulation of Na_V1.1 and Na_V1.7 by these alanine mutants affected the size of the Na⁺ peak current, and no alterations in the kinetics (tau) of fast inactivation were observed in the presence of these peptides under experimental conditions (data not shown). While Tap1a-WT displayed 3.7-fold preference for Na_V1.7 over Na_V1.1, the mutant D27A had 1.2-fold preference for Na_V1.1 over Na_V1.7, and Q33A had 1.6-fold preference for Na_V1.7 over Na_V1.1 (Figure 1F). The primary structure of Tap1a shows key residues for Na_V-activity revealed from alanine scanning in the loops 1, 2 and 4 and the C-terminal regions (Figure 1G), while the three-dimensional model of Tap1a showed most of these residues clustered at one face of the peptide, while the residues P12, D15 and Q28 were located adjacent to this cluster (Figure 1H).

Optimized Tap1a Displayed Increased Potency for Na_V Channels In Vitro

We designed and produced two optimized versions of the peptide Tap1a (Tap1a-WT) here named Tap1a-OPT1 and Tap1a-OPT2 (Figure 2). Our rational design focused on residue substitutions that warranted additional or new biochemical properties introduced to Tap1a-WT and were in line with the findings for the optimized peptides here named GpTx-1-OPT (Murray et al., 2015), HwTx-IV-OPT (Revell et al., 2013; Rahnama et al., 2017) and CcoTx-1-OPT (Shcherbatko et al., 2016). As a result of the rational design, Tap1a-OPT1 comprised the substitutions D1G, D2G, M6I, N13E, P19Y, K22V, R25K, D27H, Q28R, Y32W, Q33K, while Tap1a-OPT2 included the additional substitution F7A (Figures 2A–C). These newly designed peptides were successfully produced via recombinant expression and purified as described in Figures 2D–G.

FIGURE 2

FIGURE 2. Rational design of Tap1a-OPT1 and Tap1a-OPT2, and peptide production via recombinant expression. (A) Alignment of NaSpTx1 spider peptides which were pharmacologically optimized for Na_V channels inhibition (Cardoso and Lewis, 2019), and Tap1a, Tap1a-OPT1 and Tap1a-OPT2. Residues coloured in red are positions where substitutions led to Na_V-activity loss reported in previous studies and for Tap1a in this work, and residues coloured in blue and green are substitutions that enhanced potency and selectivity for Na_Vs, respectively. Residue substitutions in Tap1a-OPT1 and Tap1a-OPT2 are highlighted in blue. Cysteines are in grey boxes. (B, C) Model of the three-dimensional structures of Tap1a-OPT1 (B) and Tap1a-OPT2 (C) highlighting in blue the substitutions introduced in Tap1a-OPT1 and in pink the additional substitution F7A introduced in Tap1a-OPT2. (D–G) Production of the peptides Tap1a-OPT1 and Tap1a-OPT2. RP-HPLC chromatogram and mass spectrometry analysis results of the purified recombinants Tap1a-OPT1 (D, E) and Tap1a-OPT2 (F, G) purified in C18 column using a TFA/ACN gradient. Mass spectrometry analysis of purified peptides displayed the expected masses of 4143 Da for Tap1a-OPT1 (E) and 4066 Da for Tap1a-OPT2 (G).

The ability of Tap1a-OPT1 and Tap1a-OPT2 to inhibit Na_V channels was measured via automated whole cell patch-clamp electrophysiology and their IC₅₀ values determined for Na_V1.1 to Na_V1.7 channel subtypes (Figure 3; Table 2). Tap1a-OPT1 displayed a significant increase in potency for Na_V1.1 and Na_V1.7 (p = 0.0005 and p = 0.001, respectively), as well as for Na_V1.2 (p = 0.01), Na_V1.3 (p = 0.009) and Na_V1.6 (p = 0.0003), and Na_V1.4 (p < 0.00002) compared to Tap1a-WT, while no significant inhibition of Na_V1.5 was observed at up to 1 μM (Figures 3A–D; Table 2). The substitution F7A in Tap1a-OPT1 produced Tap1a-OPT2 that displayed significant weaker Na_V inhibition compared to Tap1a-OPT1 for Na_V1.1 (p = 0.02) and Na_V1.7 (p = 0.0002), Na_V1.2 (p = 0.0007), Na_V1.3 (p = 0.002) and Na_V1.6 (p = 0.007) (Figures 3A–D; Table 2), while Tap1a-OPT2 did not significantly inhibit Na_V1.4 and Na_V1.5 at up to 1 μM (Figure 3C; Table 2). Nonetheless, Tap1a-OPT2 maintained significant enhanced inhibition of Na_V1.1 (p = 0.0009), Na_V1.3 (p = 0.01), Na_V1.6 (p = 0.004) and Na_V1.7 (p = 0.004) compared to Tap1a-WT (Figures 3A–D; Table 2).

FIGURE 3

FIGURE 3. Na_V channel subtypes pharmacological profile of Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 determined using automated electrophysiology. (A–C) Representative concentration-responses for the inhibition of the hNa_V1.1 to hNa_V1.7 channels in the presence of increasing concentrations of Tap1a-WT (circles), Tap1a-OPT1 (triangles) and Tap1a-OPT2 (squares) measured by automated whole-cell patch clamp in QPatch 16X. Holding potential was –80 mV and Na⁺ currents were elicited by 20 ms voltage step to 0 mV from a –120 mV conditioning pulse applied for 200 ms. The IC₅₀ values calculated are described in Table 2 (D) Comparison of the pIC₅₀ values for Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 tested against hNa_V1.1 to hNa_V1.7 and calculated from the dose responses produced in (A–C). *Statistically significant compared to Tap1a-WT; ^#Statistically significant compared to Tap1a-OPT1; Individual p values are described in the text.

TABLE 2

TABLE 2. Pharmacological activity of the spider peptide Tap1a-WT and optimized analogues Tap1a-OPT1 and Tap1a-OPT2 for hNa_V channels. The values for the Tap1a-WT were obtained from our previous study (Cardoso et al., 2021). Peptides were tested at up to 1 μM concentration.

Remarkable enhancement in inhibition by Tap1a-OPT1 was observed for the subtypes Na_V1.1–Na_V1.3, Na_V1.6 and Na_V1.7 with increases of 18 to 41-fold compared to Tap1a-WT (Table 2). Tap1a-OPT2 inhibited Na_V channels similarly to Tap1a-WT and displayed 2- to 6-fold potency increase compared to Tap1a-WT. Interestingly, Tap1a-OPT2 displayed a more selective pharmacology away from the off-targets hNa_V1.4 and hNa_V1.5 compared to Tap1a-WT as observed by the remaining Na⁺ currents in the presence of these peptides tested at up to 1 µM (Figure 3C).

Enhanced Inhibitory Properties of Tap1a-OPT Occur via Decreased Off-Rate

The rates of modulation of hNa_V1.1 and hNa_V1.7 by Tap1a-WT and Tap1a-OPT1 were measured in the presence of 10x the IC₅₀ for each Na_V subtype (Figure 4). The on-rate evaluations revealed Tap1a-OPT1 produced a slight but significant increase in the on-rate for the hNa_V1.1 subtype when compared to Tap1a-WT (Figure 4A and Table 3). Tap1a-OPT1 showed on-rate of 26 s at hNa_V1.1 and was 1.9-fold faster than Tap1a-WT with on-rate of 48.4 s, and had a significant increase in K_on observed compared to Tap1a-WT. At hNa_V1.7, Tap1a-OPT1 showed on-rate of 135 s and was 1.2-fold faster than Tap1a-WT with an on-rate of 163 s but not statistically significant (Figure 4B). Both Tap1a-WT and Tap1a-OPT1 produced statistically significant decreased K_on observed and faster on-rates to inhibit hNa_V1.1 compared to hNa_V1.7 (Figures 4A,B; Table 3). Tap1a-WT was 3.3-fold faster to inhibit hNa_V1.1 compared to hNa_V1.7, and Tap1a-OPT1 was 5.3-fold faster to inhibit hNa_V1.1 compared to hNa_V1.7.

FIGURE 4

FIGURE 4. Kinetics of Na_V1.1 and Na_V1.7 inhibition by Tap1a-WT and Tap1a-OPT1. (A, B) On-rate measurements for Tap1a-WT and Tap1a-OPT1 over Na_V1.1 (A) and Na_V1.7 (B) were calculated from fitted exponentials of experiments recorded before and after addition of peptides at 10x IC₅₀ concentration. (C, D) Off-rate measurements of Tap1a-WT and Tap1a-OPT1 over Na_V1.1 and Na_V1.7 were calculated from the fitted exponentials of experiments recorded before and after addition of peptides at 10x IC₅₀ concentration followed by washes with extracellular solution. (E) Na_V currents measured at the end of 30 min for Na_V1.l and 25 min for Na_V1.7 ligand wash-out period. **p = 0.0042. Statistical significance was determined by t-student test. All data points are mean ± SEM of 4–7 independent experiments.

TABLE 3

TABLE 3. Kinetic properties of Tap1a-WT and Tap1a-OPT1 at hNa_V1.1 and hNa_V1.7 determined via patch-clamp electrophysiology. Kinetic data were determined using peptide concentrations equivalent to 10× IC₅₀, and k_obs and k_off are given as the mean ± SEM from 4–7 independent experiments.

The rates of dissociation of Tap1a-WT and Tap1a-OPT1 from hNa_V1.1 and hNa_V1.7 revealed Tap1a-OPT1 had a statistically significant decreased off-rate at both hNa_V1.1 and hNa_V1.7 subtypes compared to Tap1a-WT (Figures 4C,D; Table 3). While Tap1a-WT was washed off from hNa_V1.7 with an off-rate of 19.7 min, Tap1a-OPT1 was nearly irreversible at the same experimental conditions and K_off was not able to be calculated. Similarly, the rates of dissociation from hNa_V1.1 showed Tap1a-WT was washed off with an off-rate of 25.41 min, and Tap1a-OPT1 was nearly irreversible at the same experimental conditions. The irreversibility of hNa_V1.1 inhibition by Tap1a-WT was more pronounced as per remaining Na⁺ currents measured at 30 min wash out and significantly higher than hNa_V1.7 as described in Figure 4E.

Tap1a-OPT1 Displayed Improved Target Engagement and Therapeutic Effects In Vivo

To assess the efficacy of Tap1a-WT and Tap1a-OPT1 to inhibit Na_V-induced nocifensive responses in vivo, we utilized a model of Na_V1.7-mediated pain based on the intraplantar injection of the Na_V1.7 activator OD1 (Deuis et al., 2016). Intraplantar injection of Tap1a-WT (WT, 1 μM) or Tap1a-OPT1 (OPT, 1 μM) significantly reduced OD1-induced spontaneous pain behaviors at the time points indicated (Figure 5A), albeit Tap1a-WT had a slower onset of action and was less effective compared to Tap1a-OPT1 (Figure 5B). In line with this, only Tap1a-OPT1 significantly reduced total pain behaviors, which is consistent with its improved potency at peripheral Na_Vs observed in vitro (cumulative flinches: control 181 ± 15; Tap1a WT 140 ± 10; Tap1a OPT1 72 ± 12; Figure 5B).

FIGURE 5

FIGURE 5. Tap1a-OPT1 has enhanced in vivo target engagement and antinociceptive effects. (A) Time course of spontaneous pain behaviors induced by intraplantar injection of the Na_V1.7 activator OD1. Intraplantar injection of Tap1a wild-type (Tap1a-WT, 1 μM) or Tap1a-OPT1 (1 μM) significantly reduced OD1-induced spontaneous pain behaviors at the time points indicated. Statistical significance was determined by two-way ANOVA with Dunnett’s post-test. (B) Total OD1-induced pain behaviors over 20 min. Only Tap1a-OPT1 (1 μM) significantly reduced cumulative pain behaviors. Statistical significance was determined by one-way ANOVA with Dunnett’s post-test. Data is presented as mean ± SEM, n = 4 per mice per group. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 compared to control.

Interactions of Tap1a-WT and Tap1a-OPT With Na_V1.7 and Na_V1.4

The interactions of Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 with Na_V channels were predicted by molecular docking using HADDOCK webserver (van Zundert et al., 2016) (Figure 6). Potential interactions between Tap1a-WT and Na_V1.7 VSDII revealed strong hydrogen bonds (h-bonds) with distances ranging from (in Å) 1.5 to 2.9 (Figures 6Ai,C; Supplementary Table S2). Most h-bond interactions occurred between the loop 4 and the C-terminal of Tap1a-WT and the S2 and S3–S4 loop of Na_V1.7 VSDII. Hydrophobic interactions with distances ranging from (in Å) 3.4 to 3.9 were observed between loop 1 and the C-terminal of Tap1a-WT and the S2 and S3–S4 loop of Na_V1.7 VSDII (Figures 6Aii,C; Supplementary Table S2).

FIGURE 6

FIGURE 6. Molecular model of the Tap1a and Tap1a-OPT1 docked on the human Na_V1.7 channel DII. (Ai) The three-dimensional representation of the interactive surface of Tap1a-WT (salmon) docked over the hNa_V1.7 channel voltage sensor domain II / Na_VAb chimera (Na_V1.7 DII) (light blue) are represented in cartoon. Side chains are represented in sticks. Hydrogen bonds with distances up to 3Å are represented by dashed lines colored in yellow. (Aii) Hydrophobic interacting surfaces of Tap1a-WT colored in yellow and the hNa_V1.7 DII in orange. Zoom in at the Tap1a-WT C-terminal W35 showing a potential hydrophobic interaction with F812 in the S3 of Na_V1.7 DII. Distance between the closet carbons of their aromatic rings is represented in yellow dashed lines. Side chains are represented by sticks. (Bi) The three-dimensional representation of the interactive surface of Tap1a-OPT1 (pink) docked over the hNa_V1.7 DII (light blue) are represented in cartoon. Side chains are represented by sticks and hydrogen bonds with distances up to 3Å are represented by dashed lines colored in yellow. (Bii) Hydrophobic interacting surfaces of Tap1a-OPT1 colored in yellow and the hNa_V1.7 DII in orange. Zoom in at the Tap1a-OPT1 C-terminal W35 showing a potential hydrophobic interaction with L770 in the S2 of Na_V1.7 DII. Distance between the closet carbons of their side chains is represented in yellow dashed lines. Side chains are represented by sticks. (C) Primary structure of Tap1a-WT and Tap1a-OPT1 colored in red the residues predicted to form the strongest interactions between Tap1a-WT and hNa_V1.7 DII, and in green the residues predicted to form the strongest interactions between Tap1a-OPT1 and hNa_V1.7 DII. (D, E) Electrostatic surface of Tap1a-WT (D) and Tap1a-OPT1 (E). Their predicted interactive ionic surface is visualized in the top right representations. The electrostatic potential of the interacting surfaces of Tap1a-WT or Tap1a-OPT1 with the Na_V1.7 DII are visualized in the bottom representation of each figure. Electron density charge bars range from –5.0 (red) to +5.0 (blue).

The predicted interactions of Tap1a-OPT1 with Na_V1.7 VSDII revealed h-bonds with distances ranging from (in Å) 1.7 to 3.1 (Figures 6Bi,C; Supplementary Table S2). These occurred at loops 3 and 4 and C-terminal of Tap1a-OPT1 binding to the S3 and S3–S4 loop of the Na_V1.7 VSDII. Hydrophobic interactions with distances ranging from (in Å) 3.2 to 4.5 were also observed between loops 3 and 4 and C-terminal of Tap1a-OPT1 and S2 and S3–S4 loop of Na_V1.7 VSDII (Figures 6Bii,C; Supplementary Table S2). In an experiment that docked Tap1a-WT over Nav1.7 DII using the optimization restrains, we observed loops 3 and 4 and C-terminal forming the strongest h-bonds and hydrophobic interactions with the S3–S4 loop of the Na_V1.7 VSDII (Supplementary Table S2).

The electrostatic surface of Tap1a-WT and Tap1a-OPT1 were calculated, as well as the electrostatic interacting surface of these peptides with the Na_V1.7 VSDII as represented in Figures 6D,E. The interacting surface of Tap1a-OPT1 was favoured via residues that simultaneously form stronger hydrophobic or h-bonds interactions and additional electrostatic interactions with Na_V1.7 VSDII, such as W32 and K33. Some of the Tap1a-WT predicted interacting residues forming hydrophobic or h-bonds are negatively charged or neutral, such as D27 and Q33, respectively, and likely produce opposed forces and weaker interactions with the negatively charged exposed surface of Na_V1.7 VSDII.

In order to predict interactions with Na_V1.4 VSDII S3–S4 loop that could explain the loss of Na_V1.4 inhibition by Tap1a-OPT2, we docked Tap1a-OPT1 and Tap1a-OPT2 over Na_V1.4 (Supplementary Table S2). Our results showed F7 in Tap1a-OPT1 or A7 in Tap1a-OPT2 likely do not interact with Na_V1.4 VSDII. Instead, we observed an overall weakening of hydrophobic interactions between Tap1a-OPT2 and Na_V1.4 VSDII (Supplementary Table S2). The docking between Tap1a-WT and Na_V1.4 using the alanine scanning results as restrains also showed a potential reduction in hydrophobic interactions, while F7 was in the non-interacting surface of the Tap1a and away from the ion channel (Supplementary Table S2).

Discussion

Venoms from species of spider, scorpion, cone snail and snake are unique cocktails of bioactive compounds with potential for novel drugs and therapies. These compounds range from small molecules to peptides and proteins, from which a few are approved and available for human use such as ziconotide derived from the ω-conopeptide MVIIA isolated from Conus magus and effective in the treatment of neuropathic pain (Deer et al., 2019), and captopril derived from proline-rich oligopeptides isolated from the Brazilian snake Bothrops jararaca and effective in the treatment of hypertension and congestive heart failure (Ferreira et al., 1970). While MVIIA was not required to undergo optimization to become ziconotide, the proline-rich oligopeptides derived from snake venom were optimized via structure-active relationship investigations that led to a small molecule derived from the terminal sulfhydryl moiety of these oligopeptides (Patlak, 2004). A number of venom-derived compounds have been investigated in pre-clinical trials (Pennington et al., 2018; Bordon et al., 2020) or are in more advanced stages of clinical trials such as chlorotoxin derivates for tumour painting in intraoperative visualization (Patil et al., 2019), and the ShK-186 analogue derived from the venom peptide ShK from sea anemone for the treatment of autoimmune diseases such as psoriatic arthritis and lupus (Pennington et al., 2015). These illustrate the exquisite potential of venom peptides and great feasibility of translation into novel therapeutics for human use.

Tap1a has Conserved Structure-Function Properties as the NaSpTx1 Family

Recent advancements in the understanding of the structure-function relationships of Na_V-targeting spider peptides allows a refined rational design of novel optimized inhibitors for the development of drugs and therapies (Cardoso and Lewis, 2019; Cardoso, 2020). Poorly treated diseases in which Na_V channels are proved to underly pathological processes, e.g., chronic pain, motor neuron disease and epilepsy, would benefit enormously from these novel optimized Na_V-targeting molecules (Cardoso, 2020). The elucidation of the binding sites of spider-peptide inhibitors to Na_V channels revealed they preferably bind to the VSDII S3-S4 loop of these channels to engender inhibitory effects (Bosmans et al., 2008; Xiao et al., 2008; Shcherbatko et al., 2016; Cardoso et al., 2017). In this work, we showed Tap1a engendered inhibition through alike mechanism via interactions with the VSDII S3–S4 loop of Na_V1.1. This mechanism was investigated using a K_V2.1/Na_V1.1 VSDII channel chimera, and our results showed Tap1a at nanomolar concentration interacted with the Na_V1.1 VSDII S3–S4 loop which is in line with our previous findings with the wild-type Na_V1.1 channel α subunit associated to the β1 subunit (Cardoso et al., 2021).

Remarkable work elucidated the structure-function and/or optimized the spider peptides CcoTx-1 (Shcherbatko et al., 2016), HwTx-IV (Minassian et al., 2013; Revell et al., 2013; Rahnama et al., 2017) and GpTx-1 (Murray et al., 2015; Murray et al., 2016), all belonging to the NaSpTx1 (Klint et al., 2012). Besides, CcoTx-1 and HwTx-IV were previously shown to preferably bind to VSDII S3-S4 loop of Na_V channels (Xiao et al., 2008; Xiao et al., 2010; Shcherbatko et al., 2016). The strong resemblance of mechanisms of action and structure-function properties of NaSpTx1 peptides encouraged studies of Tap1a, also in NaSpTx1 family, through select alanine scanning and guided rational design for the optimization of Na_V-inhibitory properties.

In line with the NaSpTx1 family, Tap1a lost Na_V inhibitory activity by individual alanine substitutions of hydrophobic residues located at the centre of loop 1 and at the C-terminal, P and D in loop 2 and S and Q in loop 4. Contrary to HwTx-IV, the alanine substitution of a Q at the C-terminal of Tap1a enhanced its activity for Na_V1.1 and nearly did not alter its activity for Na_V1.7, indicating that Tap1a still maintain intrinsic properties associated to its structure-function relationships and unique potency and selectivity for Na_V subtypes. We suggested from our molecular docking studies that most of the Tap1a wild-type interactive surface is in the loop 4 and C-terminal, with fewer potential hydrophobic interactions in loop 1. In addition, the three-dimensional model of the Tap1a structure suggested the residues P12, D15 and Q28 were not located in the binding surface and this way are likely involved in the maintenance of Tap1a’s active structure. Experiments investigating the secondary and three-dimensional structures of Tap1a and derived analogues are essential to further elucidate the interactions between Tap1a and Na_V channels.

Tap1a Binding Mode to Na_V Channels is Broadly Conserved in NaSpTx Peptides

The binding mode of NaSpTx peptides to Na_V channels has been studied in detail to unravel a conserved binding mechanism here named as “Lysin and Arginine electrostatic trap”. This mode uses strong electrostatic interactions driven by basic/positive charged lysin and arginine residues generally located in the loop 4 and C-terminal of a peptide that interacts with negatively charged residues such as glutamic acid and aspartic acid in S3–S4 loop of VSDII, and in which interactions are enhanced by vicinal hydrophobic interactions both within the channel and/or the cellular membrane. This has been shown in detail for NaSpTx1 peptides, including HwTx-IV (Wisedchaisri et al., 2021), GpTx-1 (Murray et al., 2016) and CcoTx-1 (Shcherbatko et al., 2016). For example, the binding of GpTx-1 to the hNav1.7 was studied by single substitutions and molecular docking, revealing the phenylalanine in loop 1 (F5) and the lysin in the C-terminal (K31) interacted with I767 in S2 and E811 of the S3–S4 loop of VSDII, respectively (Murray et al., 2016). In the present study, Tap1a wild-type showed a phenylalanine in loop 1 (F6) interacted with A766 in S2 and a lysin at the C-terminal (K31) interacted with E811 in the S3–S4 loop.

In Tap1a-OPT1, we hypothesized that these interactions were enhanced by the replacement of aspartic acid by lysin in position 33 which provided additional electrostatic interactions which are not present in the Tap1a wild-type. Similarly, in NaSpTx2 the peptide Pn3a showed interactions between K22 and K24 in loop 4 with the D816 and E818 in the S3–S4 segment predicted by molecular docking (Mueller et al., 2020). In NaSpTx-3, a similar mechanism was proposed for the peptide ProTx-II where R22 and K26 in loop 4 and C-terminal, respectively, interacted with E811 and D816 in S3–S4 loop of VSDII of hNa_V1.7 (Xu et al., 2019). Overall, this conserved inhibition mechanism uses the insertion of positive charges often via lysin and arginine residues in the S3–S4 cleft of the VSDII to produce an electrostatic repulsion against S4 and the requirement for more positive potentials to induce the upward movement of S4 and the opening of the pore.

Although Cryo-EM and X-ray studies of Na_V channels have disclosed their three-dimensional structure in detail, these were mostly performed with channels at the open-state in which the voltage-sensor S4 is upwards and the channel pore open (Pan et al., 2018; Pan et al., 2019; Shen et al., 2019; Pan et al., 2021; Wisedchaisri et al., 2021), or did not describe key residues in the loop S3–S4 of VSDII which are essential for peptide binding (Shen et al., 2019). We performed docking experiments with Tap1a-WT and Tap1a-OPT1 on these available channel structures, and the models that produced closer interactions between VSDII and Tap1a used the structure determined for the complex hNa_V1.7 VSDII-ProTx-II (Xu et al., 2019), and the hNa_V1.4 channel structure determined in complex with the β1 auxiliary subunit (Pan et al., 2018). Our findings with these channel structures are consistent with the structure-function relationship findings for other NaSpTx peptides (Cardoso and Lewis, 2019; Shen et al., 2019; Xu et al., 2019; Mueller et al., 2020; Wisedchaisri et al., 2021).

Tap1a was Optimized via SAR-Incorporation to Enhance Pharmacology In Vitro and Bioacivity In Vivo

Optimized Na_V-targeting spider peptides display an increased affinity for the Na_V α subunit and the cell membrane, and hence improved potency. It is known the reduction in negative charge at the N-terminal, loop 1 or loop 3 of spider peptides can improve Na_V inhibition through an increase in affinity for the cell membrane (Henriques et al., 2016; Agwa et al., 2017). These occur in the NaSpTx1 and NaSpTx3 peptides (Cardoso and Lewis, 2019). In Tap1a-OPT1, the substitutions D1G and D2G at the N-terminal, D27H and Q28R at loop 4, and Q33K at the C-terminal likely increased its affinity for the cell membrane. These suggest that some residues could function by interacting with the cell membrane to initiate the binding process, and then by forming strong electrostatic interactions with the Na_V α subunit leading to the final binding fit.

The activity kinetics observed in our study implies the optimization of Na_V inhibition by Tap1a-OPT1 occurred via a significant slower off-rate, and this mechanism was conserved in both hNa_V1.1 and hNa_V1.7 subtypes. In a previous report, we demonstrated the enhancement of pharmacological Na_V-inhibition by venom-derived bivalent ligand peptides produced via a similar slower off-rate mechanism that conferred improved potency at Na_V1.4 (Peschel et al., 2020). Additional studies investigating lipids binding and ion channel interactions are essential to unravel if the slower off-rate is mediated via interactions with the Na_V α subunit and/or via interactions with the cell membrane.

An interesting observation from this study was the key role of lysine residues in the loop 4 and C-terminal to the binding to VSDII S3–S4 loop and inhibition of Na⁺ currents. Tap1a-WT displayed strong h-bonds formed by K26 and K31 and Na_V1.7 VSDII, while the substitution K31A led a loss of activity for hNa_V1.1 and hNa_V1.7. In Tap1a-OPT1, h-bonds were formed by K31 and the new K33 in the C-terminal. Similarly, a study with the optimized peptide m3-HwTx-IV and the Na_V1.7 performed via Cryo-EM revealed K27 in loop 4 and K32 at the C-terminal as the main drivers of this inhibitory mechanism (Wisedchaisri et al., 2021). Besides forming strong h-bonds, lysin residues can effectively interact with lipids in the cell membrane as demonstrated in studies with the peptide ProTx-II (Henriques et al., 2016). The substitution of lysine residues by arginine in ProTx-II led to a loss of affinity for lipids, and the substitution E17K led to a significant increase in affinity for lipids. These observations further support the hypothesis that key residues in ICK peptides can function initially as lipid binders to facilitate the initial encounter with the ion channel and later as ion channel binders to reach the final fit and nanomolar inhibition of the channel activity.

Because Tap1a-OPT2 lost at least 10-fold activity for Na_V1.4 but only 2-fold activity for Na_V1.7 compared to Tap1a-OPT1, we suggest that the phenylalanine residue in loop 1 of Tap1a-OPT1 is critical for an effective binding to Na_V1.4. A study with GpTx-1 showed a similar trend for the analogue GpTx-F5A that showed 2.5-fold loss activity for Na_V1.7 and 32.5-fold loss of activity for Na_V1.4, and structural changes in K31 and Y32 (Lawrence et al., 2019). Interestingly, the analogue Tap1a-F7A completely lost its inhibition for Na_V1.1 and Na_V1.7 at up to 1 μM. Based on these observations, we suggest the phenylalanine in loop 1 is critical for the affinity of NaSpTx1 peptides for Na_V1.4, while the affinity of NaSpTx1 peptides for Na_V1.7 is less dependent on this residue. We must also consider the likeliness of structural changes Tap1a-F7A and Tap1a-OPT2 mutants, as observed for GpTx-1 F5A, which can allosterically affect the binding of vicinal residues to the ion channel. Further experiments via nuclear magnetic resonance to determine the three-dimensional structure of Tap1a and its derived analogues are essential to test this hypothesis.

The translation of the pharmacology optimized in vitro to enhanced bioactive effects in vivo has challenged studies with spider peptides targeting Na_V channels. For example, optimization of HwTx-IV towards lipid bilayers binding and Na_V inhibition did not result in the enhancement of bioactivity in vivo compared to HwTx-IV wild type (Agwa et al., 2020). In another study, the enhancement of Na_V inhibition observed in vitro by the addition of C-terminal amide to the spider peptide Df1a did not translate in enhanced in vivo effects (Cardoso et al., 2017). These studies suggest the C-terminal amidation could be associated to an enhancement in lipids affinity only and that this may not be sufficient to produce efficient inhibitory effects in vivo such as the interactions with the ion channel. In this present study, the rational optimization of Tap1a generating Tap1a-OPT1 produced a new therapeutic lead with remarkable and significant increase in bioactive effects in vivo compared to the Tap1a wild type. In addition, we designed Tap1a-OPT2 with selectivity away from Na_V off-targets and greater potential to provide lesser side effects in more complex animal models of diseases. Such results are very encouraging and can guide the rapid optimization of other peptide Na_V inhibitors to produce leads with useful improved in vivo therapeutic effects.

Conclusions and Perspectives

Research for the development of novel drugs to tackle unmet needs in complex neurological diseases has progressed significantly in the last few years. Advancements in the understanding of disease mechanisms (Cardoso, 2020) and the structure function properties of key ion channel drug targets (Cardoso and Lewis, 2019; Mueller et al., 2020; Wisedchaisri et al., 2021) have provided key support for this progress. Although a significant number of ion channel modulators are available (Cardoso et al., 2018; Cardoso and Lewis, 2018), spider peptides have exceled in their therapeutic potential as proved by academic and industry research.

Spider peptides interact with the extracellular domains of voltage-gated ion channels and allosterically alter the channel conformation to facilitate (activator) or impede (inhibitor) the upward movement of voltage-sensor segments and the opening of the ion-selective pore. This modulatory mechanism is useful in the context of diseases with pathologies underlay by Na_V-mediated dysfunctional neuronal signalling as observed for chronic visceral pain in IBS. For example, the spider peptide Hm1a, an allosteric activator of Na_V1.1, produced chronic visceral pain-like symptoms in mice (Osteen et al., 2016), and the spider peptide Tap1a, an allosteric inhibitor of Na_V channels, reversed chronic visceral pain mediated by Na_V hyperfunction (Cardoso et al., 2021).

By applying the available structure-function relationship findings for spider peptides, we were able to rapidly produce an optimized spider peptide with enhanced in vitro pharmacology and enhanced in vivo bioactive effects. We observed this optimization produced a significant structure-function relationship improvement that likely generated a new and more efficient binding mode to the Na_V channel. The exact mechanisms through which this beneficial enhancement was achieved are still not fully clear, and we hypothesized it involved a combination of optimized interactions with the ion channel via the conserved “lysine and arginine electrostatic trap”, as well as with the lipidic cell membrane.

As perspectives from this work, further studies with chemically synthesized Tap1a analogues with C-terminal amide and non-natural amino acids to explore the chemical space available, the determination of their three-dimensional structure and binding sites on Na_V channels, along the evaluation of novel peptide leads in more relevant pre-clinical models of neurological diseases will expand this study to further support the research into novel therapeutics to treat complex neurological disorders.

Materials and Methods

Animals

For behavioral assessment we used adult male C57BL/6J mice aged 6–10 weeks. Mice were housed in groups of 4 per cage, under 12 h light-dark cycles and had standard rodent chow and water ad libitum. Ethical approval for in vivo experiments was obtained from the University of Queensland animal ethics committee. All experiments were conducted in accordance with the International Association for the Study of Pain Guidelines for the Use of Animals in Research, and the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 8th edition (2013, updated 2021).

Cell Culture

HEK293 cells were cultured in Dulbecco’s MEM supplemented with 10% FBS, 100 U.ml⁻¹ penicillin and 100 μg.ml⁻¹ streptomycin. HEK293 cells stably expressing recombinant hNa_V subtypes and the β1 auxiliary subunit (SB Drug Discovery, Glasgow, United Kingdom) were cultured in Minimal Essential medium (MEM) (Sigma-Aldrich, MO, United States) supplemented with 10% FBS, 100 U.ml⁻¹ penicillin, 100 μg.ml⁻¹ streptomycin, 2 mM L-glutamine, and variable concentrations of blasticidin, geneticin and zeocin according to manufacturer’s instructions. All cells were maintained at 37°C in a humidified 5% CO₂ incubator, and subcultured every 3–4 days in a ratio of 1:5 using 0.25% trypsin/EDTA.

Patch-Clamp Electrophysiology of Na_V1.1/K_V2.1 Chimeric Channel

The plasmid constructs containing the voltage sensor extracellular domains II of the Na_V1.1 channel inserted into the K_V2.1 channel was kindly provided by Prof Frank Bosman from the Department of Basic and Applied Medical Sciences, Ghent University, Belgium (Osteen et al., 2016). The Na_V1.1/K_V2.1 chimera was subcloned into the mammalian vector pCDNA3.1 and new constructs confirmed by DNA sequencing. Plasmid constructs were used to transfect HEK293 cells using the FuGENE reagent as previously described (Hasan et al., 2020). Briefly, HEK293 cells were transfected with 18 μg DNA using a ration of DNA:FuGENE at 1:3 and transfection media replaced with fresh media after 16 h incubation at 37°C 5% CO₂. After 48 h transfection, cells were used in the electrophysiology assays in the automated whole-cell patch clamp system QPatch 16X (Sophion Bioscience).

The extracellular solution comprised 140 NaCl, 5 KCl, 10 CaCl₂, 2 MgCl₂, 10 glucose and 10 HEPES at pH 7.4 and 320 mOsm. The intracellular solution comprised (in mM) 150 KCl, 1 MgCl₂, 4 NaCl, 0.5 EGTA and 10 HEPES at pH 7.4 and 320 mOsm. Cells were maintained at a holding potential –90 mV and K⁺ currents elicited by +20 mV pulse for 500 ms followed by –40 mV pulse for additional 500 ms. Tap1a-WT was added at increasing concentrations from 1.4 nM to 1 μM, while nifedipine was added at a single concentration of 100 μM. Electrophysiological data were analyzed using QPatch Assay Software v5.6 (Sophion Bioscience).

Patch-Clamp Electrophysiology of Na_V Channels

Na_V channels currents were recorded from HEK293 cells expressing human Na_V subtypes co-expressed with the β1 auxiliary subunit (SB Drug Discovery) using the automated whole-cell patch clamp system QPatch 16X (Sophion Bioscience A/S, Ballerup, Denmark). The extracellular solution comprised (in mM) 1 CaCl₂, 1 MgCl₂, 5 HEPES, 3 KCl, 140 NaCl, 0.1 CdCl₂, 20 TE A-Cl, at pH 7.3 and 320 mOsm. The intracellular solution comprised (in mM) 140 CsF, 1/5 EGTA/CsOH, 10 HEPES, 10 NaCl at pH 7.3 and 320 mOsm.

Cells were maintained at a holding potential –80 mV and Na⁺ currents elicited by 20 ms voltage steps to 0 mV from a –120 mV conditioning pulse applied for 200 ms. Voltage-activation relationships were obtained by measuring steady-state Na⁺ currents elicited by step depolarizations from –110 to +80 mV using 5-mV increments. The peak conductance (G_Na) was calculated from G = I/(V–V_rev), where I, V and V_rev represent the current value, membrane potential, and reversal potential, respectively. The voltage of steady-state fast inactivation was estimated using a double-pulse protocol with currents elicited by a 20 ms depolarizing potential of 0 mV from a 500 ms pre-pulse to potentials ranging from –130 to –10 mV using 10-mV increments.

To obtain concentration–response curves, cells at holding potential were incubated for 5 min with increasing concentrations of Tap1a-OPT1 and Tap1a-OPT2. The voltage-dependence of activation and inactivation were determined either in the absence or following 5-min exposure to Tap1a-OPT1 at the respective IC₅₀ concentration. For on-rate experiments, Na⁺ currents were measured at 15 s intervals over 15 min immediately following the addition of peptide at a concentration equivalent to 10 times their IC₅₀ for the Na_V subtype being analysed. For off-rate measurements, cells were incubated with peptide for 10 min at a concentration equivalent to 10 times their IC₅₀ for the Na_V subtype being analysed, and Na⁺ currents were assessed at 10-s intervals during saline washes. The k_on, k_off and K_d values were calculated using K_d = k_off/k_on (nM), where k_off = 1/τ_off (s⁻¹) and k_on=[(1/τ_on)−K_off]/(Maertens et al., 2006) (nM⁻¹S⁻¹). Data were analysed using Assay software (Sophion Biosciences) and Na⁺ currents (I_Na) plotted as I/Imax.

Alanine Scanning of Tap1a

Alanine substitutions in the primary sequence of Tap1a were made based in structure-activity relationship (SAR) studies of other spider peptides belonging to NaSpTx1 (Cardoso and Lewis, 2019). The alanine mutants were produced using the QuikChange Lightning site-directed mutagenesis kit and the Tap1a-pLICC plasmid construct as directed by the manufacturer (Agilent, CA, United States).

Briefly, primers containing the substitutions M6A, F7A, P12A, D15A, S24A, D27A, Q28A, W29A, K31A, Y32A, Q33A, L34A and W35A were designed and used in polymerase chain reactions (PCR) to generate alanine mutants. After digestion with the restriction enzyme DpnI to eliminate the Tap1a wild-type construct from the reaction, the newly synthesized mutants were transformed into E. coli strain DH5-α and plated in LB agar containing 100 μg.ml⁻¹ ampicillin, 80 μg.ml⁻¹ X-gal and 20 mM IPTG. White colonies were collected, and plasmid DNA extracted using QIAprep Spin Miniprep kit (QIAGEN, Hilden, DEU) as directed by the manufacturer. The new Tap1a-alanine mutants constructs were confirmed by DNA sequencing.

Rational Design of Optimized Tap1a

The rational design of Tap1a used the available information on the structure-activity relationship (SAR) studies of other spider peptides belonging to NaSpTx1 (Cardoso and Lewis, 2019) (Figure 3A). The following substitutions to the primary sequence of Tap1a were made to create Tap1a-OPT1: D1G, D2G, M6I, N13E, P19Y, K22V, R25K, D27D, Q28R, Y32W and Q33K. Synthetic gene encoding Tap1a-OPT1, with codons optimized for expression in E. coli, was produced and subcloned into the pLICC vector by GeneArt (Life Technologies, CA, United States). The Tap1a-OPT2 expression plasmid was generated by the replacement F7A in Tap1a-OPT1 using the QuikChange Lightning site-directed mutagenesis kit as described above. The new Tap1a-OPT2 construct was confirmed by DNA sequencing.

Production of Recombinant Peptides

The peptide mutants were produced using the QuikChange Lightning site-directed mutagenesis kit and the Tap1a-pLICC plasmid construct as directed by the manufacturer (Agilent, CA, United States) with codons optimized for expression in Escherichia coli and subcloned into the pLICC vector by GeneArt (Life Technologies, CA, United States). This expression vector contains a MalE signal sequence for periplasmic export, a His₆ tag fused to the maltose binding protein (MBP), and a tobacco etch virus (TEV) protease recognition site preceding the peptide gene. Peptide expression was performed as described previously (Cardoso et al., 2015).

Briefly, plasmids were transformed into the protease-deficient E. coli strain BL21 (DE3), then cells were cultured at 37°C in Luria-Bertani medium. Peptide expression was induced at OD₆₀₀ = 1 with 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG), then cells were harvested after 16 h at 16°C and pelleted by centrifugation for 20 min at 10,000 rpm. Cells were lysed using FastBreak lysis buffer (Promega, WI, United States) containing 0.2 mg.ml⁻¹ lysozyme and 10 U.ml⁻¹ DNAse.

The fusion protein was captured by passing the cell lysate through a NI-NTA resin (Sigma-Aldrich, Aldrich, MO, United States), then the His₆-MBP-toxin fusion protein was eluted with 500 mM imidazole in TN buffer (25 mM Tris, 150 mM NaCl, pH 8). After desalting, the purified fusion protein was cleaved with TEV protease, then the peptide was purified using RP-HPLC performed using an Ultimate 3000 LC system (Dionex, Sunnyvale, CA, United States) with a C18 column (Vydac 4.6 mm × 250 mm, 5 μM). The peptide was eluted using a gradient of 10–50% solvent B in Solvent A, using a flow rate of 0.7 ml.min⁻¹. Peaks were collected at 0.7 ml per well, then fractions lyophilized and stored at –20°C.

OD1-Induced Spontaneous Pain

The Na_V1.7 activator OD1 (300 nM) was diluted in saline/0.1% BSA and administered by intraplantar injection into the left hind paw of adult (6–10 weeks) male C57BL/6 mice under light isoflurane (4%) anesthesia as previously described (Cardoso et al., 2015; Deuis et al., 2016). Tap1a WT (1 μM) or Tap1a OPT (1 μM) were co-administered by intraplantar injection with OD1. Immediately after injection mice were placed in polyvinyl boxes (10 × 10 × 10 cm), and spontaneous pain behaviors, which included flinching, licking, and shaking of the injected hind paw, were recorded by video and later counted by a blinded observer unaware of the treatment each individual mouse received.

Molecular Docking

The molecular docking of Tap1a-WT and Tap1a-OPT1 over the hNa_V1.7 VSDII was performed using HADDOCK2.4 (van Zundert et al., 2016) and visualized and analyzed using PyMol (DeLano, 2002). The Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 molecular models were produced using as template the NMR structure of the tarantula toxin μ-TRTX-Pre1a-W7A (PDB 512P, unpublished results, 71.4% identity with Tap1a-WT), and models generated in SWISS-MODEL (Arnold et al., 2006).

For the molecular docking experiments over hNa_V1.7, the cryo-EM structure of the chimera hNa_V1.7VSDII-Na_VAb on its deactivated state and S4 downwards was used (PDB 6N4R) (Xu et al., 2019). Residues involved in the bioactivity of Tap1a-WT for Na_V1.7 channels determined using an alanine scanning (M6, F7, P12, D15, S24, D27, Q28, K31, Y32, Q33, L34 and W35) or rational optimization (G1, G2, I6, E13, Y19, V22, K25, H27, R28, W32 and K33) were used as restrains. Additional restrains were applied at F812, D815 and E817 located in the Na_V1.7 VSDII S3-S4 loop which were previously described to participate in the interactions with spider peptides (Xiao et al., 2008; Xiao et al., 2010).

For the molecular docking experiments over hNa_V1.4, the cryo-EM structure of the human Nav1.4 was used (PDB 6AGF) (Pan et al., 2018). Modified residues in the sequences of Tap1a-OPT1 (G1, G2, I6, E13, Y19, V22, K25, H27, R28, W32 and K33) and Tap1a-OPT2 (OPT1 plus A7) were used as restrains. Additional restrains were applied at N661 and Q663 located in the Na_V1.4 voltage-sensor domain II S3–S4 loop and previously described to participate in the interactions with spider peptides (Xiao et al., 2008; Xiao et al., 2010).

The HADDOCK cluster result displaying the most negative HADDOCK Z-score was selected for further analysis. H-bond interaction distances in Å were calculated automatically in Pymol, while hydrophobic interaction distances were calculated manually in Pymol using as criteria the closest distance between the carbons located in the side chains. The electrostatic surface potential of Tap1a-WT, Tap1a-OPT1 and hNav1.7 DII were calculated using the linearized Poisson-Boltzmann equation through the APBS Electrostatics plugin for PyMol (Dolinsky et al., 2004).

Data Analysis

For the in vitro electrophysiological recordings, curve fitting was performed using GraphPad Prism Version 8 (GraphPad Software, San Diego, CA, United States) using nonlinear regression with log-inhibitor versus normalized response and variable Hill slope for dose-responses and IC₅₀ determination, and exponential one-phase association and dissociation for on and off-rate analysis, respectively. Potencies expressed in IC₅₀ values were compared amongst Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 using multiple t-student test with Welch correction. The kinetics of interaction and dissociation were compared amongst Tap1a-WT, Tap1a-OPT1 and Tap1a-OPT2 using Mann-Whitney test. Data are mean ± SEM, n = number of individual cells recorded. For the in vivo experiments, statistical significance was determined by one-way ANOVA with Dunnett’s post-test and data was presented as mean ± SEM, n = 4 per mice per group.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by The University of Queensland animal ethics committee. All experiments were conducted in accordance with the International Association for the Study of Pain Guidelines for the Use of Animals in Research, and the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 8th edition (2013, updated 2021).

Author Contributions

Conceptualization: FC; Funding acquisition: FC and RL; Performed experiments and analysis: HH, SM, TZ, JD, SJ, and FC; Wrote the manuscript: FC; All authors revised and edited the manuscript.

Funding

This work was supported by the Australian National Health and Medical Research Council (NHMRC) Ideas Grant GNT1188959 (to FC and RL) and NHMRC Program Grant APP1072113 (to RL). FC was supported by The University of Queensland. JD, RL, and IV were supported by research fellowships from the Australian National Health and Medical Research Council (APP1162503, APP1119056, and APP1139961, respectively).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors kindly thank Prof Frank Bosmans from the Department of Basic and Applied Medical Sciences, Ghent University, Belgium, for providing the plasmid construct of the Na_V1.1/K_V2.1 chimera. The FP7 WeNMR (project# 261572), H2020 West-Life (project# 675858) and the EOSC-hub (project# 777536) European e-Infrastructure projects are acknowledged for the use of their web portals, which make use of the EGI infrastructure with the dedicated support of CESNET-MetaCloud, INFN-PADOVA, NCG-INGRID-PT, TW-NCHC, SURFsara and NIKHEF, and the additional support of the national GRID Initiatives of Belgium, France, Italy, Germany, Netherlands, Poland, Portugal, Spain, United Kingdom, Taiwan and the US Open Science Grid.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.742457/full#supplementary-material

References

Agwa, A. J., Lawrence, N., Deplazes, E., Cheneval, O., Chen, R. M., Craik, D. J., et al. (2017). Spider Peptide Toxin HwTx-IV Engineered to Bind to Lipid Membranes Has an Increased Inhibitory Potency at Human Voltage-Gated Sodium Channel hNa V 1.7. Biochim. Biophys. Acta (Bba) - Biomembranes 1859, 835–844. doi:10.1016/j.bbamem.2017.01.020

CrossRef Full Text | Google Scholar

Agwa, A. J., Tran, P., Mueller, A., Tran, H. N. T., Deuis, J. R., Israel, M. R., et al. (2020). Manipulation of a Spider Peptide Toxin Alters its Affinity for Lipid Bilayers and Potency and Selectivity for Voltage-Gated Sodium Channel Subtype 1.7. J. Biol. Chem. 295, 5067–5080. doi:10.1074/jbc.ra119.012281

CrossRef Full Text | Google Scholar

Arnold, K., Bordoli, L., Kopp, J., and Schwede, T. (2006). The SWISS-MODEL Workspace: a Web-Based Environment for Protein Structure Homology Modelling. Bioinformatics 22, 195–201. doi:10.1093/bioinformatics/bti770

PubMed Abstract | CrossRef Full Text | Google Scholar

Bladen, C., Hamid, J., Souza, I. A., and Zamponi, G. W. (2014). Block of T-type Calcium Channels by Protoxins I and II. Mol. Brain 7, 36. doi:10.1186/1756-6606-7-36

PubMed Abstract | CrossRef Full Text | Google Scholar

Bordon, K. d. C. F., Cologna, C. T., Fornari-Baldo, E. C., Pinheiro-Júnior, E. L., Cerni, F. A., Amorim, F. G., et al. (2020). From Animal Poisons and Venoms to Medicines: Achievements, Challenges and Perspectives in Drug Discovery. Front. Pharmacol. 11, 1132. doi:10.3389/fphar.2020.01132

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosmans, F., Martin-Eauclaire, M.-F., and Swartz, K. J. (2008). Deconstructing Voltage Sensor Function and Pharmacology in Sodium Channels. Nature 456, 202–208. doi:10.1038/nature07473

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, T., Luo, J., Meng, E., Ding, J., Liang, S., Wang, S., et al. (2015). Mapping the Interaction Site for the Tarantula Toxin Hainantoxin-IV (β-TRTX-Hn2a) in the Voltage Sensor Module of Domain II of Voltage-Gated Sodium Channels. Peptides 68, 148–156. doi:10.1016/j.peptides.2014.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C., Castro, J., Grundy, L., Schober, G., Garcia-Caraballo, S., Zhao, T., et al. (2021). A Spider-Venom Peptide with Multitarget Activity on Sodium and Calcium Channels Alleviates Chronic Visceral Pain in a Model of Irritable Bowel Syndrome. Pain 162, 569–581. doi:10.1097/j.pain.0000000000002041

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C., Dekan, Z., Rosengren, K. J., Erickson, A., Vetter, I., Deuis, J. R., et al. (2015). Identification and Characterization of ProTx-III [μ-TRTX-Tp1a], a New Voltage-Gated Sodium Channel Inhibitor from Venom of the Tarantula Thrixopelma Pruriens. Mol. Pharmacol. 88, 291–303. doi:10.1124/mol.115.098178

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C., Dekan, Z., Smith, J. J., Deuis, J. R., Vetter, I., Herzig, V., et al. (2017). Modulatory Features of the Novel Spider Toxin μ‐TRTX‐Df1a Isolated from the Venom of the Spider Davus Fasciatus. Br. J. Pharmacol. 174, 2528–2544. doi:10.1111/bph.13865

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C., Hasan, M., Zhao, T., and Lewis, R. J. (2018). Toxins in Pain. Curr. Opin. Support. Palliat. Care 12, 132–141. doi:10.1097/spc.0000000000000335

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C., and Lewis, R. J. (2018). Sodium Channels and Pain: from Toxins to Therapies. Br. J. Pharmacol. 175, 2138–2157. doi:10.1111/bph.13962

CrossRef Full Text | Google Scholar

Cardoso, F. C., and Lewis, R. J. (2019). Structure-Function and Therapeutic Potential of Spider Venom-Derived Cysteine Knot Peptides Targeting Sodium Channels. Front. Pharmacol. 10, 366. doi:10.3389/fphar.2019.00366

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F. C. (2020). Multi-targeting Sodium and Calcium Channels Using Venom Peptides for the Treatment of Complex Ion Channels-Related Diseases. Biochem. Pharmacol. 181, 114107. doi:10.1016/j.bcp.2020.114107

PubMed Abstract | CrossRef Full Text | Google Scholar

Catterall, W. A. (2017). Forty Years of Sodium Channels: Structure, Function, Pharmacology, and Epilepsy. Neurochem. Res. 42, 2495–2504. doi:10.1007/s11064-017-2314-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chey, W. D., Kurlander, J., and Eswaran, S. (2015). Irritable Bowel Syndrome. JAMA 313, 949–958. doi:10.1001/jama.2015.0954

PubMed Abstract | CrossRef Full Text | Google Scholar

Deer, T. R., Pope, J. E., Hanes, M. C., and Mcdowell, G. C. (2019). Intrathecal Therapy for Chronic Pain: a Review of Morphine and Ziconotide as Firstline Options. Pain Med. 20, 784–798. doi:10.1093/pm/pny132

PubMed Abstract | CrossRef Full Text | Google Scholar

DeLano, W. L. (2002). The PyMOL Molecular Graphics System. San Carlos, CA: DeLano Scientific.

Deuis, J. R., Dekan, Z., Wingerd, J. S., Smith, J. J., Munasinghe, N. R., Bhola, R. F., et al. (2017). Pharmacological Characterisation of the Highly NaV1.7 Selective Spider Venom Peptide Pn3a. Sci. Rep. 7, 40883. doi:10.1038/srep40883

PubMed Abstract | CrossRef Full Text | Google Scholar

Deuis, J., Wingerd, J., Winter, Z., Durek, T., Dekan, Z., Sousa, S., et al. (2016). Analgesic Effects of GpTx-1, PF-04856264 and CNV1014802 in a Mouse Model of NaV1.7-Mediated Pain. Toxins 8, 78. doi:10.3390/toxins8030078

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolinsky, T. J., Nielsen, J. E., Mccammon, J. A., and Baker, N. A. (2004). PDB2PQR: an Automated Pipeline for the Setup of Poisson-Boltzmann Electrostatics Calculations. Nucleic Acids Res. 32, W665–W667. doi:10.1093/nar/gkh381

PubMed Abstract | CrossRef Full Text | Google Scholar

Enck, P., Aziz, Q., Barbara, G., Farmer, A. D., Fukudo, S., Mayer, E. A., et al. (2016). Irritable Bowel Syndrome. Nat. Rev. Dis. Primers 2, 16014. doi:10.1038/nrdp.2016.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Erickson, A., Deiteren, A., Harrington, A. M., Garcia-Caraballo, S., Castro, J., Caldwell, A., et al. (2018). Voltage-gated Sodium Channels: (NaV)igating the Field to Determine Their Contribution to Visceral Nociception. J. Physiol. 596, 785–807. doi:10.1113/jp273461

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, B., Zhu, Y., La, J.-H., Wills, Z. P., and Gebhart, G. F. (2015). Experimental and Computational Evidence for an Essential Role of NaV1.6 in Spike Initiation at Stretch-Sensitive Colorectal Afferent Endings. J. Neurophysiol. 113, 2618–2634. doi:10.1152/jn.00717.2014

CrossRef Full Text | Google Scholar

Ferreira, S. H., Bartelt, D. C., and Greene, L. J. (1970). Isolation of Bradykinin-Potentiating Peptides from Bothrops jararaca Venom. Biochemistry 9, 2583–2593. doi:10.1021/bi00815a005

PubMed Abstract | CrossRef Full Text | Google Scholar

Flinspach, M., Xu, Q., Piekarz, A. D., Fellows, R., Hagan, R., Gibbs, A., et al. (2017). Insensitivity to Pain Induced by a Potent Selective Closed-State Nav1.7 Inhibitor. Sci. Rep. 7, 39662. doi:10.1038/srep39662

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundy, L., Erickson, A., Caldwell, A., Garcia-Caraballo, S., Rychkov, G., Harrington, A., et al. (2018). Tetrodotoxin-sensitive Voltage-Gated Sodium Channels Regulate Bladder Afferent Responses to Distension. Pain 159, 2573–2584. doi:10.1097/j.pain.0000000000001368

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasan, M., Ragnarsson, L., Cardoso, F. C., and Lewis, R. J. (2020). Transfection Methods for High-Throughput Cellular Assays of Voltage-Gated Calcium and Sodium Channels Involved in Pain. PLoS One 16 3, e0243645. doi:10.1371/journal.pone.0243645

CrossRef Full Text | Google Scholar

Henriques, S. T., Deplazes, E., Lawrence, N., Cheneval, O., Chaousis, S., Inserra, M., et al. (2016). Interaction of Tarantula Venom Peptide ProTx-II with Lipid Membranes Is a Prerequisite for its Inhibition of Human Voltage-Gated Sodium Channel NaV1.7. J. Biol. Chem. 291, 17049–17065. doi:10.1074/jbc.m116.729095

PubMed Abstract | CrossRef Full Text | Google Scholar

Klint, J. K., Senff, S., Rupasinghe, D. B., Er, S. Y., Herzig, V., Nicholson, G. M., et al. (2012). Spider-venom Peptides that Target Voltage-Gated Sodium Channels: Pharmacological Tools and Potential Therapeutic Leads. Toxicon 60, 478–491. doi:10.1016/j.toxicon.2012.04.337

PubMed Abstract | CrossRef Full Text | Google Scholar

Klint, J. K., Smith, J. J., Vetter, I., Rupasinghe, D. B., Er, S. Y., Senff, S., et al. (2015). Seven Novel Modulators of the Analgesic Target NaV1.7 Uncovered Using a High-Throughput Venom-Based Discovery Approach. Br. J. Pharmacol. 172, 2445–2458. doi:10.1111/bph.13081

CrossRef Full Text | Google Scholar

Kubat Öktem, E., Mruk, K., Chang, J., Akin, A., Kobertz, W. R., and Brown, R. H. (2016). Mutant SOD1 Protein Increases Nav1.3 Channel Excitability. J. Biol. Phys. 42, 351–370. doi:10.1007/s10867-016-9411-x

CrossRef Full Text | Google Scholar

Lawrence, N., Wu, B., Ligutti, J., Cheneval, O., Agwa, A. J., Benfield, A. H., et al. (2019). Peptide-membrane Interactions Affect the Inhibitory Potency and Selectivity of Spider Toxins ProTx-II and GpTx-1. ACS Chem. Biol. 14, 118–130. doi:10.1021/acschembio.8b00989

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Wu, Z., Tang, D., Xun, X., Liu, L., Li, X., et al. (2014). Analgesic Effects of Huwentoxin-IV on Animal Models of Inflammatory and Neuropathic Pain. Protein Pept. Lett. 21, 153–158. doi:10.2174/09298665113206660119

PubMed Abstract | CrossRef Full Text | Google Scholar

Maertens, C., Cuypers, E., Amininasab, M., Jalali, A., Vatanpour, H., and Tytgat, J. (2006). Potent Modulation of the Voltage-Gated Sodium Channel Nav1.7 by OD1, a Toxin from the Scorpion Odonthobuthus Doriae. Mol. Pharmacol. 70, 405–414. doi:10.1124/mol.106.022970

PubMed Abstract | CrossRef Full Text | Google Scholar

Middleton, R. E., Warren, V. A., Kraus, R. L., Hwang, J. C., Liu, C. J., Dai, G., et al. (2002). Two Tarantula Peptides Inhibit Activation of Multiple Sodium Channels†. Biochemistry 41, 14734–14747. doi:10.1021/bi026546a

PubMed Abstract | CrossRef Full Text | Google Scholar

Minassian, N. A., Gibbs, A., Shih, A. Y., Liu, Y., Neff, R. A., Sutton, S. W., et al. (2013). Analysis of the Structural and Molecular Basis of Voltage-Sensitive Sodium Channel Inhibition by the Spider Toxin Huwentoxin-IV (μ-TRTX-Hh2a). J. Biol. Chem. 288, 22707–22720. doi:10.1074/jbc.m113.461392

CrossRef Full Text | Google Scholar

Mueller, A., Dekan, Z., Kaas, Q., Agwa, A. J., Starobova, H., Alewood, P. F., et al. (2020). Mapping the Molecular Surface of the Analgesic NaV1.7-Selective Peptide Pn3a Reveals Residues Essential for Membrane and Channel Interactions. ACS Pharmacol. Transl. Sci. 3, 535–546. doi:10.1021/acsptsci.0c00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Murray, J. K., Ligutti, J., Liu, D., Zou, A., Poppe, L., Li, H., et al. (2015). Engineering Potent and Selective Analogues of GpTx-1, a Tarantula Venom Peptide Antagonist of the NaV1.7 Sodium Channel. J. Med. Chem. 58, 2299–2314. doi:10.1021/jm501765v

CrossRef Full Text | Google Scholar

Murray, J. K., Long, J., Zou, A., Ligutti, J., Andrews, K. L., Poppe, L., et al. (2016). Single Residue Substitutions that Confer Voltage-Gated Sodium Ion Channel Subtype Selectivity in the NaV1.7 Inhibitory Peptide GpTx-1. J. Med. Chem. 59, 2704–2717. doi:10.1021/acs.jmedchem.5b01947

PubMed Abstract | CrossRef Full Text | Google Scholar

Osteen, J. D., Herzig, V., Gilchrist, J., Emrick, J. J., Zhang, C., Wang, X., et al. (2016). Selective Spider Toxins Reveal a Role for the Nav1.1 Channel in Mechanical Pain. Nature 534, 494–499. doi:10.1038/nature17976

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, X., Li, Z., Zhou, Q., Shen, H., Wu, K., Huang, X., et al. (2018). Structure of the Human Voltage-Gated Sodium Channel Nav1.4 in Complex with β1. Science 362, eaau2486. doi:10.1126/science.aau2486

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, X., Li, Z., Huang, X., Huang, G., Gao, S., Shen, H., et al. (2019). Molecular Basis for Pore Blockade of Human Na+ Channel Nav1.2 by the μ-conotoxin KIIIA. Science 363, 1309–1313. doi:10.1126/science.aaw2999

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, X., Li, Z., Jin, X., Zhao, Y., Huang, G., Huang, X., et al. (2021). Comparative Structural Analysis of Human Na_V1.1 and Na_V1.5 Reveals Mutational Hotspots for Sodium Channelopathies. Proc. Natl. Acad. Sci. U S A. 118, e2100066118. doi:10.1073/pnas.2100066118

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, C. G., Walker, D. G., Miller, D. M., Butte, P., Morrison, B., Kittle, D. S., et al. (2019). Phase 1 Safety, Pharmacokinetics, and Fluorescence Imaging Study of Tozuleristide (BLZ-100) in Adults with Newly Diagnosed or Recurrent Gliomas. Neurosurgery 85, E641–E649. doi:10.1093/neuros/nyz125

PubMed Abstract | CrossRef Full Text | Google Scholar

Patlak, M. (2004). From viper's Venom to Drug Design: Treating Hypertension. FASEB j. 18, 421. doi:10.1096/fj.03-1398bkt

PubMed Abstract | CrossRef Full Text | Google Scholar

Pennington, M., Chang, S., Chauhan, S., Huq, R., Tajhya, R., Chhabra, S., et al. (2015). Development of Highly Selective Kv1.3-Blocking Peptides Based on the Sea Anemone Peptide ShK. Mar. Drugs 13, 529–542. doi:10.3390/md13010529

PubMed Abstract | CrossRef Full Text | Google Scholar

Pennington, M. W., Czerwinski, A., and Norton, R. S. (2018). Peptide Therapeutics from Venom: Current Status and Potential. Bioorg. Med. Chem. 26, 2738–2758. doi:10.1016/j.bmc.2017.09.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Peschel, A., Cardoso, F. C., Walker, A. A., Durek, T., Stone, M. R. L., Braga Emidio, N., et al. (2020). Two for the price of One: Heterobivalent Ligand Design Targeting Two Binding Sites on Voltage-Gated Sodium Channels Slows Ligand Dissociation and Enhances Potency. J. Med. Chem. 63, 12773–12785. doi:10.1021/acs.jmedchem.0c01107

PubMed Abstract | CrossRef Full Text | Google Scholar

Priest, B. T., Blumenthal, K. M., Smith, J. J., Warren, V. A., and Smith, M. M. (2007). ProTx-I and ProTx-II: Gating Modifiers of Voltage-Gated Sodium Channels. Toxicon 49, 194–201. doi:10.1016/j.toxicon.2006.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahnama, S., Deuis, J. R., Cardoso, F. C., Ramanujam, V., Lewis, R. J., Rash, L. D., et al. (2017). The Structure, Dynamics and Selectivity Profile of a NaV1.7 Potency-Optimised Huwentoxin-IV Variant. PLoS One 12, e0173551. doi:10.1371/journal.pone.0173551

PubMed Abstract | CrossRef Full Text | Google Scholar

Revell, J. D., Lund, P.-E., Linley, J. E., Metcalfe, J., Burmeister, N., Sridharan, S., et al. (2013). Potency Optimization of Huwentoxin-IV on hNav1.7: a Neurotoxin TTX-S Sodium-Channel Antagonist from the Venom of the Chinese Bird-Eating Spider Selenocosmia Huwena. Peptides 44, 40–46. doi:10.1016/j.peptides.2013.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, K. L., Milligan, C. J., Richardson, R. J., Jancovski, N., Grunnet, M., Jacobson, L. H., et al. (2018). Selective NaV1.1 Activation Rescues Dravet Syndrome Mice from Seizures and Premature Death. Proc. Natl. Acad. Sci. USA 115, E8077–E8085. doi:10.1073/pnas.1804764115

PubMed Abstract | CrossRef Full Text | Google Scholar

Rong, M., Duan, Z., Chen, J., Li, J., Xiao, Y., and Liang, S. (2013). Native Pyroglutamation of Huwentoxin-IV: a post-translational Modification that Increases the Trapping Ability to the Sodium Channel. PLoS One 8, e65984. doi:10.1371/journal.pone.0065984

PubMed Abstract | CrossRef Full Text | Google Scholar

Saba, L., Viscomi, M. T., Martini, A., Caioli, S., Mercuri, N. B., Guatteo, E., et al. (2019). Modified Age-dependent Expression of NaV1.6 in an ALS Model Correlates with Motor Cortex Excitability Alterations. Neurobiol. Dis. 130, 104532. doi:10.1016/j.nbd.2019.104532

PubMed Abstract | CrossRef Full Text | Google Scholar

Salvatierra, J., Castro, J., Erickson, A., Li, Q., Braz, J., Gilchrist, J., et al. (2018). NaV1.1 Inhibition Can Reduce Visceral Hypersensitivity. JCI Insight 3, e121000. doi:10.1172/jci.insight.121000

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmalhofer, W. A., Calhoun, J., Burrows, R., Bailey, T., Kohler, M. G., Weinglass, A. B., et al. (2008). ProTx-II, a Selective Inhibitor of NaV1.7 Sodium Channels, Blocks Action Potential Propagation in Nociceptors. Mol. Pharmacol. 74, 1476–1484. doi:10.1124/mol.108.047670

PubMed Abstract | CrossRef Full Text | Google Scholar

Shcherbatko, A., Rossi, A., Foletti, D., Zhu, G., Bogin, O., Galindo Casas, M., et al. (2016). Engineering Highly Potent and Selective Microproteins against Nav1.7 Sodium Channel for Treatment of Pain. J. Biol. Chem. 291, 13974–13986. doi:10.1074/jbc.m116.725978

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, H., Liu, D., Wu, K., Lei, J., and Yan, N. (2019). Structures of Human Nav1.7 Channel in Complex with Auxiliary Subunits and Animal Toxins. Science 363, 1303–1308. doi:10.1126/science.aaw2493

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, K., Sekino, S., Ikegami, M., Ikeda, H., and Kamei, J. (2015). Antihyperalgesic Effects of ProTx-II, a Nav1.7 Antagonist, and A803467, a Nav1.8 Antagonist, in Diabetic Mice. J. Exp. Pharmacol. 7, 11–16. doi:10.2147/JEP.S79973

PubMed Abstract | CrossRef Full Text | Google Scholar

van Zundert, G. C. P., Rodrigues, J. P. G. L. M., Trellet, M., Schmitz, C., Kastritis, P. L., Karaca, E., et al. (2016). The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. J. Mol. Biol. 428, 720–725. doi:10.1016/j.jmb.2015.09.014

CrossRef Full Text | Google Scholar

Wisedchaisri, G., Tonggu, L., Gamal El-Din, T. M., Mccord, E., Zheng, N., and Catterall, W. A. (2021). Structural Basis for High-Affinity Trapping of the NaV1.7 Channel in its Resting State by Tarantula Toxin. Mol. Cell 81, 38–48. doi:10.1016/j.molcel.2020.10.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Y., Bingham, J.-P., Zhu, W., Moczydlowski, E., Liang, S., and Cummins, T. R. (2008). Tarantula Huwentoxin-IV Inhibits Neuronal Sodium Channels by Binding to Receptor Site 4 and Trapping the Domain Ii Voltage Sensor in the Closed Configuration. J. Biol. Chem. 283, 27300–27313. doi:10.1074/jbc.m708447200

CrossRef Full Text | Google Scholar

Xiao, Y., Blumenthal, K., Jackson, J. O., Liang, S., and Cummins, T. R. (2010). The Tarantula Toxins ProTx-II and Huwentoxin-IV Differentially Interact with Human Nav1.7 Voltage Sensors to Inhibit Channel Activation and Inactivation. Mol. Pharmacol. 78, 1124–1134. doi:10.1124/mol.110.066332

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Y., Jackson, J. O., Liang, S., and Cummins, T. R. (2011). Common Molecular Determinants of Tarantula Huwentoxin-IV Inhibition of Na+ Channel Voltage Sensors in Domains II and IV. J. Biol. Chem. 286, 27301–27310. doi:10.1074/jbc.m111.246876

CrossRef Full Text | Google Scholar

Xu, H., Li, T., Rohou, A., Arthur, C. P., Tzakoniati, F., Wong, E., et al. (2019). Structural Basis of Nav1.7 Inhibition by a Gating-Modifier Spider Toxin. Cell 176, 702–715. doi:10.1016/j.cell.2018.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar