PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(3): e0173551.
Published online 2017 March 16. doi:  10.1371/journal.pone.0173551
PMCID: PMC5354290

The structure, dynamics and selectivity profile of a NaV1.7 potency-optimised huwentoxin-IV variant

Robert Blum, Editor

Abstract

Venom-derived peptides have attracted much attention as potential lead molecules for pharmaceutical development. A well-known example is Huwentoxin-IV (HwTx-IV), a peptide toxin isolated from the venom of the Chinese bird-eating spider Haplopelma schmitdi. HwTx-IV was identified as a potent blocker of a human voltage-gated sodium channel (hNaV1.7), which is a genetically validated analgesic target. The peptide was promising as it showed high potency at NaV1.7 (IC50 ~26 nM) and selectivity over the cardiac NaV subtype (NaV1.5). Mutagenesis studies aimed at optimising the potency of the peptide resulted in the development of a triple-mutant of HwTx-IV (E1G, E4G, Y33W, m3-HwTx-IV) with significantly increased potency against hNaV1.7 (IC50 = 0.4 ± 0.1 nM) without increased potency against hNaV1.5. The activity of m3-HwTx-IV against other NaV subtypes was, however, not investigated. Similarly, the structure of the mutant peptide was not characterised, limiting the interpretation of the observed increase in potency. In this study we produced isotope-labelled recombinant m3-HwTx-IV in E. coli, which enabled us to characterise the atomic-resolution structure and dynamics of the peptide by NMR spectroscopy. The results show that the structure of the peptide is not perturbed by the mutations, whilst the relaxation studies reveal that residues in the active site of the peptide undergo conformational exchange. Additionally, the NaV subtype selectivity of the recombinant peptide was characterised, revealing potent inhibition of neuronal NaV subtypes 1.1, 1.2, 1.3, 1.6 and 1.7. In parallel to the in vitro studies, we investigated NaV1.7 target engagement of the peptide in vivo using a rodent pain model, where m3-HwTx-IV dose-dependently suppressed spontaneous pain induced by the NaV1.7 activator OD1. Thus, our results provide further insight into the structure and dynamics of this class of peptides that may prove useful in guiding the development of inhibitors with improved selectivity for analgesic NaV subtypes.

Introduction

Voltage-gated ion channels are ubiquitous signalling molecules that respond to changes in the cell membrane polarisation caused by ion concentration gradients between the intra- and extracellular environment [1]. They are particularly important in the initiation and propagation of electrical signals in excitable cells. Dysfunction of these transmembrane proteins can lead to a range of severe channelopathies including epilepsy, cardiac arrhythmias and various pain disorders. Their central role in these pathologies makes them attractive pharmaceutical targets [2].

One of the most attractive drug targets to emerge in recent years in the ion-channel field is the NaV1.7 subtype of voltage-gated sodium channels (NaVs). Mutations in SCN9A, the gene encoding NaV1.7, which lead to a loss of function, result in congenital insensitivity to pain (CIP) with anosmia as the only other sensory deficit, whilst gain-of-function mutations are associated with extreme sensitivity to pain (erythromelalgia and paroxysmal extreme pain disorder) [3, 4]. Thus, inhibitors of NaV1.7 are promising analgesic leads, and indeed, local anasthetics such as lidocaine are known to non-selectively inhibit NaV channels [2]. Unfortunately, off-target effects at other homologous NaVs (1.1, 1.2, 1.4, 1.5 and 1.6) are likely to lead to severe side-effects including seizures, arrhythmias and impaired motor function [5].

Small disulfide-rich peptides isolated from the venom of animals are natural modulators of voltage-gated ion channels. The larger footprint of these peptides, compared to small molecule drugs, is thought to provide an increased opportunity for development of more selective ligands [6]. Peptides are therefore considered a goldilocks compromise between small molecules and antibodies, with the former providing limited selectivity and the latter being associated with high production costs. Thus, there are currently significant efforts in industry and academia to discover novel peptides as lead molecules for development of biologics that target ion channels [614].

Perhaps the best-studied venom peptide inhibitor of NaV1.7 is the 35-residue HwTx-IV, isolated from the venom of the Chinese bird-eating spider Haplopelma schmitdi [15]. The peptide completely inhibits NaV1.7 with an IC50 in the low nanomolar range (~26 nM) [16]. Detailed characterisation of the toxin-channel interaction revealed that the peptide binds to one of the four voltage sensor domains (VSD) of the channel [16, 17]. This is in contrast to extant small molecule drugs that bind to the central pore region, such as local anaesthetics [18]. The channel pore is more highly conserved among the different NaV subtypes compared to the VSD, making the VSD an attractive target for development of selective drugs, and small molecules that target the VSD are being developed [7, 11].

NaVs contain four homologous but non-identical VSDs that control the gating of the channel. Each VSD consists of four transmembrane helical segments (S1-S4) connected via intra- and extra-cellular loops. VSDs of domains II and IV (VSDII and VSDIV) have been shown to be attractive targets as they control channel opening and inactivation, respectively. Elegant work by Xiao et al. showed that HwTx-IV specifically binds to acidic residues in the S1-S2 and S3-S4 loops of VSDII; rat NaV1.2, NaV1.3 and human NaV1.7, all of which contain E818 in S3-S4, are sensitive to the toxin whilst NaV1.4 and NaV1.5, which lack an acidic residue in this position, are insensitive to HwTx-IV. Subsequent work by others defining the pharmacophore of the peptide through extensive mutagenesis identified that in particular residues W30 and K32 are required for NaV activity, whilst mutations in the N and C terminal parts of the peptide were found to improve binding to NaV1.7 without improvement in activity against NaV1.5 [12, 14]. This led to the development of a triple mutant of HwTx-IV (E1G, E4G, Y33W HwTx-IV, hereafter m3-HwTx-IV), which is one of the most potent blockers of NaV1.7 reported to date (IC50 ~ 0.5 nM) [14]. Despite the extensive work carried out on HwTx-IV, its complete subtype selectivity is not known. Similarly, although the structure of the wild type (WT) peptide and several variants have been determined [19, 20], the effect of the triple mutation on the pharmacophore is not known and a high quality structure would provide valuable information for rationalising the observed gain in potency.

To this end, we produced a recombinant, isotopically labelled version of m3-HwTx-IV, which allowed us to determine (1) a high-resolution 3D structure of the peptide using heteronuclear NMR and (2) the dynamics of the peptide using 15N relaxation experiments. Although the overall fold of m3-HwTx-IV remained the same as the parent toxin [14], our dynamics data show large amplitude motion in the pharmacophore region indicating that key residues undergo conformational exchange in the μs-ms time scale. In vitro and in vivo functional data indicate that m3-HwTx-IV has the same NaV selectivity profile as the WT toxin, and we show for the first time that the peptide is a potent inhibitor of NaV1.6 and NaV1.1 channels. The full selectivity profile of the channel reveals that NaV isoforms 1.1, 1.2, 1.3, 1.6 and 1.7 are sensitive to the toxin whilst isoforms 1.4, 1.5 and 1.8 are not. Finally, we show that the peptide has NaV1.7 on-target activity in vivo using a recently developed NaV1.7 target-engagement mouse model [21]. Our findings provide further insight into the structure, function and activity of this fascinating class of peptides in the pursuit of turning venoms into analgesic drugs.

Materials and methods

Recombinant peptide production

A codon-optimised synthetic gene encoding m3-HwTx-IV was synthesised by GeneArt (Thermo Fisher Scientific) and subcloned into an in-house, modified version of pLIC expression vector, between the KpnI and SacI restriction sites, in frame with a MalE signal sequence to direct secretion of the fusion-protein to the periplasmic space of E. coli, where disulfide-bond formation and folding of the molecule takes place [22]. Competent E. coli BL21 (DE3) cells were transformed with the plasmid encoding His6-MBP-m3-HwTx-IV fusion protein. Cells were grown on a LB agar plate supplemented with 100 μg/mL ampicillin (Amp). LB/Amp (100 μg/mL) medium was inoculated from a single colony or glycerol stock and the cells were grown for ~18 h (overnight) at 30°C whilst shaking at 120 rpm until reaching an OD600 of ~2. Fresh LB/Amp solution was inoculated with the pre-culture to an OD600 of ~0.05 under agitation at 120 rpm in 5 litre baffled flasks at 37°C until reaching an OD600 of ~0.4, at which time the growth temperature was reduced to 16°C. Protein production was induced at OD600 of ~0.8 by addition of isopropyl β-D-1-thiogalactopyranoside (IPTG) to a final concentration of 0.2 mM and growth was continued overnight at 16°C. The cells were harvested by centrifugation (6,000 g for 15 min at 4°C) and the cell pellets stored at –80°C for subsequent purification. In order to produce uniformly 13C/15N-labelled m3-HwTx-IV peptide, LB broth was replaced with M9 minimal medium supplemented with 15NH4Cl and 13C6-glucose as the sole nitrogen and carbon sources, respectively [23].

Peptide purification

Cells overexpressing His6-MBP-m3-HwTx-IV fusion protein were resuspended by stirring in lysis buffer (40 mM Tris-HCl, pH 8.0, and 400 mM NaCl) at 4°C for 1 h, or until reaching complete homogeneity. Cells were lysed by continuous flow cell disruption (Constant Systems Ltd., TS Series Benchtop; triple passage at constant pressure of 32 kpsi). The cell lysate was centrifuged (40,000 g, 30 min, 4°C) to remove cell debris. The supernatant was loaded onto a Ni-NTA column equilibrated with lysis buffer. The column was washed with 3–5 column volumes (CV) of lysis buffer, then non-specifically bound proteins were eluted with the same buffer containing 40 mM imidazole. The fusion protein was then eluted from the column with lysis buffer supplemented with 250 mM imidazole. The peak fractions were pooled and then desalted and concentrated using an Amicon Ultra-15 centrifugal filter unit (Merck Millipore) with a 30 kDa cut-off membrane. The protein was then digested to yield His6-MBP and m3-HwTx-IV peptide using 50 μg/ml tobacco etch virus (TEV) protease in a freshly prepared redox buffer containing 3 mM reduced glutathione (GSH) and 0.3 mM oxidised glutathione (GSSG), gently shaking overnight, at room temperature. Fusion protein cleavage was confirmed by SDS-PAGE, then the cleavage reaction mixture was acidified by adding 1% trifluoroacetic acid (TFA). Precipitated protein was removed by centrifugation followed by membrane filtration (0.45 μm). Reverse-phase (RP) HPLC purification of m3-HwTx-IV was performed using a C4 semi-preparative column (Phenomenex). The bound peptide was eluted at 3 ml/min with a linear gradient of 5–80% solvent B (90% Acetonitrile (ACN), 0.043% TFA) in solvent A (0.05% TFA in water). UV absorbance was monitored at 214 and 280 nm and fractions corresponding to the peptide of interest were pooled, lyophilised and resuspended in pure water. UV absorbance was measured (Thermo Scientific™ NanoDrop 2000) at 280 nm, and the concentration was calculated using the predicted extinction coefficient of the peptide. Identity and purity of the peptide was further characterised by liquid chromatography mass spectrometry (LC-MS), matrix-assisted laser desorption/ionisation time-of-flight (MALDI-TOF) mass spectrometry and analytical RP-HPLC with a C18 analytical column (Phenomenex) using a gradient of 1% solvent B/min at a flow rate of 1 ml/min.

Mass spectrometry

MALDI-TOF MS performed on a Bruker autoflex speed High-Performance MALDI-TOF System (Bruker, MA, USA) was employed to determine peptide mass. RP-HPLC fractions were mixed (2:1 v/v) with α-cyano-4-hydroxycinnamic acid (CHCA) matrix (5 mg/ml in 50/50 ACN/H2O) on a MALDI plate spot and dried at room temperature. The MALDI-TOF mass spectrum was acquired in positive reflector mode. Applied Biosystems Data Explorer™ was used for analysis of monoisotopic [M + H]+ ions MS data.

NMR structure determination

All NMR data were recorded using a 300 μl samples, added into a susceptibility-matched 5 mm micro tube (Shigemi Inc. Japan), prepared by dissolving 15N or 15N/13C-labeled m3-HwTx-IV to a final concentration of 400 μM in 20 mM sodium acetate solution, pH 5 with 5% (v/v) D2O. All spectra were acquired at 25 oC using a 900 MHz Bruker Avance NMR spectrometer equipped with a cryogenically cooled probe (cryoprobe). Resonance assignments were obtained using two-dimentional (2D) 1H-15N heteronuclear single quantum coherence (HSQC), 2D 1H-13C-HSQC, 3D CBCA(CO)NH, 3D HNCACB, 3D HNCO, 3D HBHA(CO)NH and four-dimensional (4D) HCC(CO)NH-TOCSY spectra. 3D and 4D spectra were acquired using non-uniform sampling [24, 25] and processed using the Rowland NMR Toolkit (http://rnmrtk.uchc.edu/rnmrtk/RNMRTK.html). Maximum entropy parameters were selected automatically as described previously [26]. 3D 15N and 13C (separate spectra for aliphatic and aromatic regions) filtered NOESY-HSQC spectra (200 ms mixing time) were acquired using linear sampling and processed using the Rowland NMR Toolkit. Data were analysed in CcpNmr Analysis v2.4 (University of Cambridge, UK) [27]. NOESY cross peaks were manually picked and supplemented with dihedral angles derived from chemical shifts using the TALOS software [28], and used for structure calculation using CYANA v3.0. [29] The NOEs were automatically assigned by CYANA. 200 structures were calculated and an ensemble of 20 structures with the lowest CYANA target function were selected.

Dynamics of the peptide

NMR spin relaxation

NMR experiments for measuring longitudinal (R1) and transverse relaxation rates (R2) and 1H-15N steady-state NOEs were acquired using a single-scan interleaved pulse sequence on a 700 MHz Bruker spectrometer equipped with a z-gradient cryoprobe. Data were acquired at 25°C using a 200 μM sample of 15N labeled m3-HwTx-IV. R1 spectra were recorded with 7 delay times (10*, 50, 300*, 500, 700, 900*, 1100 ms) as were R2 spectra (16*, 32, 64, 96*, 128, 160, 192 ms), where delays marked with an asterisk were recorded twice to confirm data reproducibility. Steady-state 1H-15N heteronuclear NOEs were acquired with a relaxation delay of 9 s. Spectra were processed using the Rowland NMR toolkit and analysed using CcpNMR [27].

hNaV1 activity using patch-clamp electrophysiology

Patch-clamp experiments used HEK293 cells expressing either of hNaV1.1 –hNaV1.5 and the β1 auxiliary subunit (SB Drug Discovery, Glasgow, UK) or Chinese Hamster Ovary (CHO) cells expressing hNaV1.6, hNaV1.7 and hNaV1.8/β3 (ChanTest Corp). Na+ currents were measured by automated whole-cell patch clamp (QPatch 16X; Biolin Scientific A/S, Ballerup, Denmark). The extracellular solution comprised (in mM) 1 CaCl2, 1 MgCl2, 5 HEPES, 3 KCl, 140 NaCl, 0.1 CdCl2 and 20 TEA-Cl at pH 7.3 and 320 mOsm, and the intracellular solution comprised (in mM) 140 CsF, 1/5 EGTA/CsOH, 10 HEPES and 10 NaCl at pH 7.3 and 320 mOsm. The elicited currents were sampled at 25 kHz and filtered at 4 kHz. Cells were maintained at a holding potential –80 mV and Na+ currents elicited by 20 ms voltage steps to 0 mV (+10 mV for NaV1.8) from a –120 mV conditioning pulse applied for 200 ms. To obtain concentration–response curves, cells were incubated for 5 min with increasing concentrations of toxin. All experimental data was analysed using QPatch Assay Software v5.0 (Biolin Scientific A/S, Ballerup, Denmark).

hNaV1 activity using FLIPR membrane potential assay

Human embryonic kidney (HEK) 293 cells heterologously expressing human NaV1.1–1.8 (SB Drug Discovery) were cultured in minimal essential medium (MEM) containing 10% v/v fetal bovine serum (FBS) and selection antibiotics as recommended by the manufacturer. Cells were grown in a humidified 5% CO2 incubator at 37 °C, grown to 70–80% confluence, and passaged every 3–4 days using TrypLE Express (Invitrogen). Cells were plated in 384-well clear-bottom black imaging plates (Corning, NY, USA) at a density of 10,000 to 15,000 cells per well. After 48 h, cells were loaded with (20 μl per well) red membrane potential dye (Molecular Devices, Sunnyvale CA, USA) diluted in physiological salt solution [PSS; composition (in mM) 140 NaCl, 11.5 glucose, 5.9 KCl, 1.4 MgCl2, 1.2 NaH2PO4, 5 NaHCO3, 1.8 CaCl2 and 10 HEPES pH 7.4] and incubated for 30 min at 37°C in a humidified 5% CO2 incubator. m3-HwTx-IV (10 nM– 30 μM) was diluted in PSS/0.1% bovine serum albumin (BSA) and pre-incubated for 5 min before stimulating NaV1.1–1.7 using veratridine (60 μM) or NaV1.8 using deltamethrin (150 μM). Changes in membrane potential were assessed using a fluorescence imaging plate reader (FLIPRTETRA, Molecular Devices) (excitation 515 to 545 nm, emission 565 to 625 nm) every second for 300 s after adding agonists.

In vivo activity assay

Animals

Adult male C57BL/6J mice aged 6–8 weeks were used for behavioural assessments. Mice were housed in groups of 2–4 per cage, under 12 h light-dark cycles with standard rodent chow and water ad libitum. Ethical approval for in vivo animal experiments was obtained from the University of Queensland animal ethics committee. Experiments involving animals were conducted in accordance with the Animal Care and Protection Regulation Qld (2012), the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 8th edition (2013) and the International Association for the Study of Pain Guidelines for the Use of Animals in Research.

OD1 model

The NaV1.7 activator OD1 (300 nM) was co-administered alone (control; n = 5) or with m3-HwTx-IV at a concentration of 100 nM (n = 4) or 300 nM (n = 3) diluted in saline/0.1% BSA by intraplantar (i.pl.) injection in a volume of 40 μL under light isoflurane (3%) anesthesia in mice as previously described [21]. Mice were then individually placed in polyvinyl boxes (10 x 10 x 10 cm) and allowed to recover from anesthesia. Immediately upon recovery, spontaneous pain behaviours (number of paw lifts, licks, shakes, and flinches) were recorded on video, and later counted by an investigator, unaware of the treatment each individual animal received, for 40 min post-injection at 5-min intervals. Data were plotted and analysed using GraphPad Prism. Statistical significance was defined *P < 0.05 and was determined using one-way ANOVA with Dunnett’s post-test compared to control.

Results

Production of recombinant m3-HwTx-IV

An efficient method was required to provide sufficient amounts of m3-HwTx-IV to perform the structural and functional studies. Therefore, we created a construct that enables expression of a His6-MBP-m3-HwTx-IV fusion protein in the periplasm of E. coli, where the disulfide-bond machinery of the cell is located [22]. The fusion protein was purified from the soluble cell fraction using Ni-NTA affinity chromatography (Fig 1C). m3-HwTx-IV was cleaved from the His6-MBP fusion tag using TEV protease, then semi-preparative RP-HPLC (Fig 1D) was employed to yield a single dominant disulfide-bond isoform with a purity of >98% as assessed by analytical RP-HPLC and MALDI-TOF MS (Fig 1E). The activity of the peptide against hNaV1.7 was initially confirmed using a FLIPR assay, then by automated patch-clamp electrophysiology as described below. The final yield of recombinant m3-HwTx-IV was ~0.6 mg per litre of bacterial culture.

Fig 1
Recombinant production of m3-HwTx-IV.

Structure of m3-HwTx-IV

The structure of m3-HwTx-IV was determined using multidimensional heteronuclear NMR spectroscopy. NMR has been the dominant approach for determining the structure of peptides smaller than 5 kDa, including disulfide-rich venom peptides [22, 30]. Employing bacterial expression enabled us to produce uniformly 15N/13C labelled m3-HwTx-IV. Two-dimensional (2D) 1H-15N HSQC (Fig 2A) and 3D HNCACB, CBCA(CO)NH, and HNCO experiments were used to assign the backbone resonances of m3-HwTx-IV. A 4D HCC(CO)NH-TOCSY experiment, which has the advantage of providing sidechain 1H-13C correlations [25], was used to provide side-chain assignments. Ultimately, structure calculations were based on interproton distances derived from analysis of 3D 15N-edited and 13C-edited HSQC NOESY spectra and backbone dihedral-angle restraints derived from analysis of 1H, 13C, and 15N chemical shifts.

Fig 2
(A) 1H-15N HSQC spectrum of m3-HwTx-IV. Peaks represent backbone amide NH correlations. Sidechain correlations are annotated with italic font and the folded arginine sidechains are marked with a dotted box. NH2 groups are connected using a dashed line. ...

All non-exchanging and slow exchanging hydrogen atoms were assigned, as well as all protonated carbon and nitrogen atoms (BMRB-ID: 30160). These assignments were supplemented with assignments of all C’ atoms from the HNCO experiment, which excludes assignments from residues preceding the prolines and the C-terminal residue (N10, I35). The backbone chemical shifts were used to derive 46 dihedral angle predictions (23 ϕ, 23 ψ) by TALOS-N [31]. A total of 1575 NOEs were manually picked from the three NOESY spectra (15N, 13C aromatic, 13C aliphatic) and 1537 (98%) automatically assigned by CYANA [29]. These were translated to 790 unique distance limits. The disulfide bond connectivities were derived from unambiguous NOE assignments across the disulfide bond [32] and introduced into the structure calculations as three upper and lower distance limits across the disulfide bond (Sγ-Cβ’, Cβ-Sγ’, Sγ-Sγ’).

The structure calculations generated an ensemble of 20 structures (see Table 1 for statistics, PDB-ID: 5T3M). The ensemble was well defined with low backbone and sidechain RMSDs (0.02 Å and 0.48 Å respectively for residues 2–34). This compares favourably with the previously determined structure of the WT peptide (PDB ID 1MB6; RMSDs of 0.34 Å and 1.01 Å for the backbone and sidechain atoms of residues 2–34 [19]). The m3-HwTx-IV structure shows that the aromatic ring of W30 is directly above the Hβ atom of T28, explaining the unusual upfield chemical shift of this atom at 2.48 ppm (compared with the random coil value of 4.16 ppm).

Table 1
NMR structure statisticsa.

Dynamics of m3-HwTx-IV

Heteronuclear spin relaxation is a powerful technique to probe the overall and internal dynamics of macromolecules [33]. In order to study the dynamics of m3-HwTx-IV, we determined the 15N longitudinal (R1) and transverse (R2) relaxation rates and 1H-15N heteronuclear NOEs. The R1 and 1H-15N heteronuclear NOE values provide information about fast time-scale motions in the range of ns to ps while the R2 values probe slow time-scale motions in the range of ms to μs (Fig 3). The experimentally determined R1 values were found to be largely constant along the polypeptide chain with an average value of 2.25 ± 0.01 s–1. Residues at the N- and C-termini (C2 and I35) had shorter R1 values indicative of faster motions. In contrast to R1, the R2 values varied considerably along the polypeptide chain, ranging from 1.79 ± 0.01 to 7.65 ± 0.02 s–1 with an average value of 3.83 ± 0.02 s–1. Residues in the S25–T28 stretch, displayed higher or lower values than the average R2, suggestive of conformational exchange. The 1H-15N heteronuclear NOE was almost constant along the polypeptide chain (0.59 ± 0.03) with the exception of the flexible C-terminal residue (I35), which has a negative NOE, indicating a high level of disorder.

Fig 3
1H-15N steady state NOE and longitudinal (R1) and transverse (R2) 15N relaxation rates obtained for m3-HwTx-IV at a field-strength of 700 MHz.

NaV1 subtype selectivity of m3-HwTx-IV

The full selectivity profile of m3-HwTx-IV at all human NaV channel subtypes (hNaV1.1 –hNaV1.8) stably expressed in HEK293 (hNaV1.1 –hNaV1.5) and CHO (hNaV1.6 –hNaV1.8) cells was measured using whole-cell patch-clamp electrophysiology; fitting of the concentration-response data yielded IC50 values of: hNaV1.1, 8.4 ± 1.8 nM; hNaV1.2, 11.9 ± 2.2 nM; hNaV1.3, 7.2 ± 1.6 nM; hNaV1.4, 369 ± 196 nM; hNaV1.5, >1 μM; hNaV1.6, 6.8 ± 1.5 nM; hNaV1.7, 3.3 ± 1.1 nM; hNaV1.8, >1 μM (Fig 4 and Table 2), indicating that the neuronal NaV isoforms 1.1, 1.2, 1.3, 1.6 and 1.7 are sensitive to the peptide whilst isoforms 1.4, 1.5 and 1.8 are not. The selectivity profile matched that of the parent (HwTx-IV) peptide, with the exception that activity against the NaV1.1, NaV1.6 and NaV1.8 had previously not been measured. The effects of m3-HwTx-IV on NaV channel fast inactivation at 0 mV are shown in S1 Fig for the toxin sensitive channels (tested at ~ IC50 concentrations, expect for NaV1.4 which was assessed at 1 μM). m3-HwTx-IV had no substantial effect on the time constant of fast inactivation (τ) of any of the channels tested (S1 Fig).

Fig 4
Potency and subtype selectivity of m3-HwTx-IV.
Table 2
Potency of m3-HwTx-IV on hNaV1.1–1.8 stably expressed in HEK239 or CHO cells using whole-cell patch-clamp electrophysiology in QPatch and heterologously expressed in HEK293 cells assessed using FLIPR membrane potential assay.

In vivo activity of m3-HwTx-IV

The on-target effect of m3-HwTx-IV at NaV1.7 was examined in vivo via behavioural assessment in a recently described NaV1.7 target engagement model [21]. m3-HwTx-IV significantly attenuated pain behaviours induced by local injection of the NaV1.7-selective activator OD1 [34] into the hind paw of the mouse model of NaV1.7-mediated pain (pain behaviour count: control (n = 5), 243 ± 22; m3-HwTx-IV (100 nM; n = 4), 120 ± 17; m3-HwTx-IV (300 nM; n = 3), 25 ± 18; P < 0.05; Fig 5), confirming on-target activity at NaV1.7.

Fig 5
Dose-dependent suppression of OD1-induced spontaneous pain by m3-HwTx-IV.

Discussion

Structure of HwTx-IV and mutants

m3-HwTx-IV was developed by Revell et al. [14] to improve potency of the parent compound (HwTx-IV) at NaV1.7. Here we determined whether the improved activity affected the structure or selectivity of HwTx-IV. The structure of WT HwTx-IV was previously solved by two independent groups (PDB IDs 1mb6 and 2m4x for the original and more recent HwTx-IV structures, respectively). Notably, when these two structures are aligned over the backbone atoms in the structured region (residues 5–33), the RMSD was found to be 1.84 Å, which is remarkably high (i.e., indicative of significant differences), considering that both structures were solved using primarily 2D homonuclear NOESY data.

To structurally characterise m3-HwTx-IV we used a recombinant expression method [22] to isotope label the peptide with 13C and 15N, which allowed application of high-resolution, heteronuclear 3D/4D NMR methods as well as 15N relaxation measurments. The structure obtained here using heteronuclear NMR aligns well with the original native peptide structure, 1mb6, with an RMSD of 0.7 Å over the backbone atoms (5–33), but it aligns poorly with the more recent 2m4x structure, with an RMSD of 1.75 Å over the same range of backbone atoms. The large difference between our structure and the more recent of the two structures of the WT peptide is apparent in 3 of the 4 loops of the peptide as well as the termini. It is unclear why the 2m4x structure deviates to such an extent from the other two structures. We note, however, that for the 1mb6 structure the authors recorded separate 2D NOESY spectra in H2O and D2O and it may be that this approach provided additional 1H NOEs from Hα atoms in the β-sheet regions which are often obscured in 2D NOESY spectra by the strong H2O signal. In our 3D filtered NOESY data these NOEs are obscured along the direct dimension but can often be assigned along the indirect dimension. Given the discrepancy, and the additional experimental data used in the 1mb6 structure, we have below only made comparisons to 1mb6.

Structural consequences of mutations

In the m3-HwTx-IV structure the position of W33 with respect to K21 is well defined, whereas in the WT structure the aromatic ring of Y33 is poorly defined and predominantly positioned facing the sidechains of K32 and F6. In our structure there are clear NOEs between the sidechain Hβ protons of K21 and the Hδ1 and Hε1 protons of W33 in the 3D 13C aromatic NOESY, 3D 13C aliphatic NOESY and the 3D 15N NOESY spectra. The chemical shifts of these atoms are unique in the 3D NOESY spectra and the NOE assignments are therefore unambiguous, strongly supporting the alternative sidechain conformation in our model (Fig 6).

Fig 6
Interaction of K21 and W33 sidechains.

The NMR relaxation studies revealed no evidence of structural or conformational flexibility in the C-terminal tail, except for the last residue, which appears less ordered. It is therefore unclear wether the W33/K21 interactions are not present in the WT peptide or whether these interactions were simply not resolved in the previous homonuclear NMR studies. We note that in the 2D NOESY data the Hβ atoms of K21 are in a crowded region, which may have obscured NOEs in the 1mb6 data. It had previously been speculated that cation-π interactions may serve to reduce the loss in energy upon binding due to the desolvation of lysine residues. Here, the conformation of these two residues would not provide such stabilisation, as the charged NH3+ group of the lysine is instead proximal to the carboxyl-terminus of the peptide.

The other residues that were mutated in m3-HwTx-IV are two glutamic acid residues in the N-terminal section of the peptide. Curiously the E4 residue appears to form a salt bridge with K7 in the WT toxin, and it might be expected that this charge pairing stabilises the fold of the peptide. However, in our structure K7 assumes a nearly identical position as in the WT peptide whilst there are no major changes to the conformation of loop 1. Thus, the improved NaV1.7 activity from mutating these acidic residues to Gly likely results from removal of unfavourable interactions with either VSDII of NaV1.7 or the surrounding membrane. Indeed, previous work has shown that loop 1 residues are important in lipid binding and it is possible that the acidic residues either interact unfavourably with anionic components of the membrane leaflet or reduce the permeability of the peptide in the lipid bilayer [35].

The C-terminus of HwTx-IV is amidated and previous structure/activity studies have shown that removal of this amidation causes a 20–50 fold reduction in activity at NaV1.7 [14, 36]. Since m3-HwTx-IV was produced recombinantly, the C-terminus is not amidated and we would expect the activity of the peptide to be less than that of the amidated peptide. Indeed we found that the recombinant peptide inhibited NaV1.7 with an IC50 of 3.3 nM, which is 8-fold higher than that previously reported for the synthetic amidated version but still more potent than WT HwTx-IV (IC50 26 nM). Thus, m3-HwTx-IV is one of the most potent recombinantly produced NaV1.7 blockers reported, and in contrast to ProTx-II, it does not require refolding [37].

Dynamics of the pharmacophore region

The NaV pharmacophore in this class of peptides has been extensively characterised and it is clear that the W30/K32 pair are the most functionally critical residues [14, 20]. In our structure these residues assume a conformation similar to that in WT peptide, and we conclude that the three mutations do not affect the overall configuration of these critical residues. Instead, our relaxation data show clear evidence of conformational mobility in loop 4 residues including W30. Remarkably, we find that the largest discrepancy between our mutant structure and the WT structure is in this loop, and when this loop is omitted in the structural alignment the RMSD falls from 0.7 Å to 0.4 Å. This suggests that this loop is flexible and that the difference between the structures may simply reflect different conformational states of loop 4 (see also Fig 2B). Given that this dynamic region coincides with the active site of the peptide we speculate that this mobility provides additional conformational states that may stabilise binding of the peptide to VSDII. These data also suggest that the absolute configuration of the pharmacophore residues in the solution structure of HwTx-IV may differ from that of the channel-bound peptide. Should this be the case, the solution structures may not be suitable design templates for development of peptide mimetics or small molecules inhibitors, unless sufficient flexibility is introduced with respect to the position of these key functional residues.

Subtype selectivity

The binding site of HwTx-IV was characterised in detail by Xiao et al. [16]. The authors found that HwTx-IV inhibits rat (r) NaV1.2 and rNaV1.3 as well as human NaV1.7, but is far less potent on rNaV1.4 and rNaV1.5. Based on these results and careful mutagenesis work, Xiao et al. were able to identify residues on NaVs that conferred selectivity [17]. They showed that the peptide binds to VSDII and that one residue on the S1-S2 loop (E753), together with four residues on the S3-S4 loop (E811, L814, D816, E819), are particularly important for peptide binding. This recognition motif “EELDE” is conserved in NaV1.6, whilst only differing in a D-to-N change (EELNE) in NaV1.1, NaV1.2, NaV1.3. In contrast, NaV1.4 and NaV1.5 contain two changes (EELNQ and EELRS respectively), and NaV1.8 contains three changes (EEVKG) with NaV1.9 having the biggest difference with four changes (EDVQR).

From the above analysis it should therefore be expected that NaV1.6 should be equally sensitive to m3-HwTx-IV, and since NaV1.2 is sensitive to the toxin, we would expect NaV1.1, which has the same recognition sequence, to also be m3-HwTx-IV sensitive. Finally, based on the sequence analysis, NaV1.8 and NaV1.9 are least likely to be sensitive to the toxin. Using a FLIPR membrane potential assay we compared subtype selectivity profiles across NaV isoforms and as predicted by the work in the Cummins group [17, 38] the previously uncharacterised hNaV1.6 and hNaV1.1 are indeed sensitive to the toxin whilst hNaV1.8 is, again as predicted, insensitive to the toxin. To further characterise the selectivity of m3-HwTx-IV we measured its activity against NaV1.1–1.8 using whole cell patch-clamp electrophysiology. We found that the peptide is very potent on NaV1.1–1.3, NaV1.6 and NaV1.7. Thus, it appears that the D-to-N mutation at position 816 is well tolerated. Comparing the activity on NaV1.4/5 and NaV1.1 we can conclude that final E-to-Q/S mutation is not tolerated and the presence of a negative charge is essential in this position. The full range of NaV subtypes sensitive to HwTx-IV only differs from tetrodotoxin (TTX) in its activity on NaV1.4. This may prove useful as differential inhibition of sodium currents by HwTx-IV and TTX can be attributed to NaV1.4.

Our results show that, despite claims of improved selectivity and potency, the subtype selectivity of m3-HwTx-IV is not significantly altered compared with native HwTx-IV. From the sequence analysis our data would support a model where the peptide binds to residues within the cavity formed between the four helices of the VSDII [39]. This model would suggest that there is a large interface between the peptide and the channel, indicating that it should in theory be possible to develop subtype selective peptides. Current peptide engineering efforts, however, rely on exhaustive mutagenesis of the peptide and have unfortunately shown limited success, with the common observation that improvements in affinity against one of the sensitive subtypes often leads to a correlated increase in affinity against other subtypes [40]. To overcome these challenges further structural details of the VSDII in its peptide bound form are required to enable rational peptide design.

In vivo activity

The in vivo activity of HwTx-IV was previously reported in rat models of formalin-induced pain and spared nerve injury, where it blocked pain behaviours at doses between 25 and 200 μg/kg [41]. This exquisite in vivo potency may at least in part be due to the strain of animals used, as morphine also was effective at doses at least 20 times lower than those usually required to achieve analgesia in the same models [4245]. Given the poor selectivity of HwTx-IV, in particular over NaV1.6, we used a recently developed NaV1.7 target engagement assay to assess the biological activity of recombinant m3-HwTx-IV following local intraplantar injection, as systemic administration of peptides with activity at NaV1.6 typically leads to adverse side-effects. Intraplantar injection of OD1 at the concentration used here evokes spontaneous pain behaviours that are reduced by ~70% in NaV1.7 knockout mice or after treatment with the selective NaV1.7 inhibitor PF-04856264, with residual pain responses likely attributable to NaV1.6 [21]. Accordingly, the observed in vivo effects of m3-HwTx-IV are entirely consistent with our pharmacological characterization of m3-HwTx-IV showing potent effects at NaV1.6 and NaV1.7, and indicate target engagement of NaV1.7 and to lesser extent NaV1.6 on peripheral nerve terminals in vivo. The in vivo effects of m3-HwTx-IV were only assessed in male mice, as sex-related differences in NaV1.7-mediated pain are not expected [46].

Conclusions

We report here an efficient method of recombinant production of a highly potent NaV1.7 inhibitor, m3-HwTx-IV. Our high-resolution solution-state NMR structure revealed that the mutations introduced into HwTx-IV do not alter the structure of the peptide. Thus, with perhaps the exception of the Y33W mutation, the gain in potency for the mutant peptide appears to be due to removal of unfavourable contacts from the acidic moieties (E1G, E4G). Furthermore, we provide the first NMR relaxation studies of this family of sodium channel inhibitors, which revealed significant conformational flexibility in the critical pharmacophore region. Detailed functional studies showed that the subtype selectivity of m3-HwTx-IV differs from TTX only at NaV1.4, which was not sensitive to m3-HwTx-IV. Finally, we found that m3-HwTx-IV provides analgesia in a mouse model of NaV1.7-mediated pain. While recombinant m3-HwTx-IV is likely to be a useful pharmacological tool for structural and pharmacological studies, its off-target activity makes it a poor clinical candidate. Further structural details of the peptide-channel complex are likely required to enable design of HwTx-IV analogues with improved selectivity profiles.

Supporting information

S1 Fig

m3-HwTx-IV has no effect on NaV channel fast inactivation when tested near the IC50 value.

Na+ currents were elicited by a 20 ms depolarising pulse to 0 mV from a 200 ms pre-pulse at -120 mV. Holding potential was -80 mV. The time constant of fast inactivation (τ) was calculated from a fit of the inactivation phase to a mono exponential equation and is presented as mean +/- SEM (n = 3 repetitions).

(DOCX)

Acknowledgments

We thank the Queensland NMR Network for access to the 900 MHz spectrometer.

Funding Statement

The authors would like to acknowledge financial support from the Australian National Health & Medical Research Council (Project Grant APP1034958 [GFK and MM], APP1080405 [MM], Program Grant APP569927 [RJL and GFK], Principal Research Fellowship to GFK and RJL) and the Australian Research Council (Future Fellowship to MM and IV). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

NMR chemical shifts are available from the bio-magnetic resonance bank (BMRB), with accession number (BMRB-ID:30160). The 3D structure is available form the PDB (PDBID: 5TM3).

References

1. Yu FH, Yarov-Yarovoy V, Gutman GA, Catterall WA (2005) Overview of molecular relationships in the voltage-gated ion channel superfamily. Pharmacol Rev 57: 387–395. doi: 10.1124/pr.57.4.13 [PubMed]
2. McGivern JG (2007) Advantages of voltage-gated ion channels as drug targets. Expert Opin Ther Tar 11: 265–271. [PubMed]
3. Cox JJ, Reimann F, Nicholas AK, Thornton G, Roberts E, Springell K, et al. (2006) An SCN9A channelopathy causes congenital inability to experience pain. Nature 444: 894–898. doi: 10.1038/nature05413 [PubMed]
4. Estacion M, Dib-Hajj SD, Benke PJ, Te Morsche RH, Eastman EM, Macala LJ, et al. (2008) NaV1.7 gain-of-function mutations as a continuum: A1632E displays physiological changes associated with erythromelalgia and paroxysmal extreme pain disorder mutations and produces symptoms of both disorders. J Neurosci 28: 11079–11088. doi: 10.1523/JNEUROSCI.3443-08.2008 [PubMed]
5. Catterall WA, Goldin AL, Waxman SG (2005) Nomenclature and structure-function relationships of voltage-gated sodium channels. Pharmacol Rev 57: 397–409. doi: 10.1124/pr.57.4.4 [PubMed]
6. 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 [PMC free article] [PubMed]
7. Ahuja S, Mukund S, Deng L, Khakh K, Chang E, Ho H, et al. (2015) Structural basis of NaV1.7 inhibition by an isoform-selective small-molecule antagonist. Science 350. [PubMed]
8. Bagal SK, Chapman ML, Marron BE, Prime R, Storer RI, Swain NA (2014) Recent progress in sodium channel modulators for pain. Bioorg Med Chem Lett 24: 3690–3699. doi: 10.1016/j.bmcl.2014.06.038 [PubMed]
9. Klint JK, Chin YKY, Mobli M (2015) Rational engineering defines a molecular switch that is essential for activity of spider-venom peptides against the analgesics target NaV1.7. Mol Pharmacol 88: 1002–1010. doi: 10.1124/mol.115.100784 [PubMed]
10. Klint JK, Smith J, Vetter I, Rupasinghe DB, Senff S, Er SY, et al. (2015) Seven novel modulators of the analgesic target NaV1.7 uncovered using a high-throughput venoms-based discovery approach. Br J Pharmacol 172: 2445–2458. doi: 10.1111/bph.13081 [PMC free article] [PubMed]
11. McCormack K, Santos S, Chapman ML, Krafte DS, Marron BE, West CW, et al. (2013) Voltage sensor interaction site for selective small molecule inhibitors of voltage-gated sodium channels. Proc Nat Acad Sci USA 110: 2724–2732. [PubMed]
12. Minassian NA, Gibbs A, Shih AY, Liu Y, Neff RA, Sutton SW, 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 [PMC free article] [PubMed]
13. Murray JK, Long J, Zou A, Ligutti J, Andrews KL, Poppe L, et al. (2016) Single residue substitutions that confer NaV subtype selectivity in the NaV1.7 inhibitory peptide GpTx-1. J Med Chem 59: 2704–2017. doi: 10.1021/acs.jmedchem.5b01947 [PubMed]
14. Revell JD, Lund PE, Linley JE, 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]
15. Peng K, Shu Q, Liu Z, Liang S (2002) Function and solution structure of huwentoxin-IV, a potent neuronal tetrodotoxin (TTX)-sensitive sodium channel antagonist from Chinese bird spider Selenocosmia huwena. J Biol Chem 277: 47564–47571. doi: 10.1074/jbc.M204063200 [PubMed]
16. Xiao Y, Bingham JP, Zhu W, Moczydlowski E, Liang S, Cummins TR (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 [PMC free article] [PubMed]
17. Xiao Y, Jackson JO, Liang SP, Cummins TR (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 [PMC free article] [PubMed]
18. Ragsdale DS, McPhee JC, Scheuer T, Catterall WA (1996) Common molecular determinants of local anesthetic, antiarrhythmic, and anticonvulsant block of voltage-gated Na+ channels. Proc Natl Acad Sci 93: 9270–9275. [PubMed]
19. Peng K, Shu Q, Liu ZH, Liang SP (2002) Function and solution structure of huwentoxin-IV, a potent neuronal tetrodotoxin (TTX)-sensitive sodium channel antagonist from Chinese bird spider Selenocosmia huwena. J Biol Chem 277: 47564–47571. doi: 10.1074/jbc.M204063200 [PubMed]
20. Minassian NA, Gibbs A, Shih AY, Liu Y, Neff RA, Sutton SW, et al. (2013) Analysis of the structural and molecular basis of voltage-sensitive sodium channel inhibition by the spider toxin, huwentoxin-IV (mu-TRTX-Hh2a). J Biol Chem 113. [PMC free article] [PubMed]
21. Deuis JR, Wingerd JS, Winter Z, Durek T, Dekan Z, Sousa SR, et al. (2016) Analgesic effects of GpTx-1, PF-04856264 and CNV1014802 in a mouse model of NaV1.7-mediated pain. Toxins (Basel) 8: 8030078. [PMC free article] [PubMed]
22. Klint JK, Senff S, Saez NJ, Seshadri R, Lau HY, Bende NS, et al. (2013) Production of recombinant disulfide-rich venom peptides for structural and functional analysis via expression in the periplasm of E. coli. PLoS One 8. [PMC free article] [PubMed]
23. Marley J, Lu M, Bracken C (2001) A method for efficient isotopic labeling of recombinant proteins. J Biomol NMR 20: 71–75. [PubMed]
24. Mobli M. (2015) Reducing seed dependent variability of non-uniformly sampled multidimensional NMR data. J Magn Reson 256: 60–69. doi: 10.1016/j.jmr.2015.04.003 [PubMed]
25. Mobli M, Stern AS, Bermel W, King GF, Hoch JC (2010) A non-uniformly sampled 4D HCC(CO)NH-TOCSY experiment processed using maximum entropy for rapid protein sidechain assignment. J Magn Reson 204: 160–164. doi: 10.1016/j.jmr.2010.02.012 [PMC free article] [PubMed]
26. Mobli M, Maciejewski MW, Gryk MR, Hoch JC (2007) An automated tool for maximum entropy reconstruction of biomolecular NMR spectra. Nat Methods 4: 467–468. doi: 10.1038/nmeth0607-467 [PubMed]
27. Vranken WF, Boucher W, Stevens TJ, Fogh RH, Pajon A, Llinas M, et al. (2005) The CCPN data model for NMR spectroscopy: development of a software pipeline. Proteins 59: 687–696. doi: 10.1002/prot.20449 [PubMed]
28. Shen Y, Delaglio F, Cornilescu G, Bax A (2009) TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts. J Biomol NMR 44: 213–223. doi: 10.1007/s10858-009-9333-z [PMC free article] [PubMed]
29. Guntert P (2004) Automated NMR structure calculation with CYANA. Methods Mol Biol 278: 353–578. doi: 10.1385/1-59259-809-9:353 [PubMed]
30. Vetter I, Davis JL, Rash LD, Anangi R, Mobli M, Alewood PF, et al. (2011) Venomics: a new paradigm for natural products-based drug discovery. Amino Acids 40: 15–28. doi: 10.1007/s00726-010-0516-4 [PubMed]
31. Shen Y, Bax A (2013) Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks. J Biomol NMR 56: 227–241. doi: 10.1007/s10858-013-9741-y [PMC free article] [PubMed]
32. Mobli M, King GF (2010) NMR methods for determining disulfide-bond connectivities. Toxicon 56: 849–54. doi: 10.1016/j.toxicon.2010.06.018 [PubMed]
33. Saez NJ, Mobli M, Bieri M, Chassagnon IR, Malde AK, Gamsjaeger R, et al. (2011) A dynamic pharmacophore drives the interaction between psalmotoxin-1 and the putative drug target acid-sensing ion channel 1a. Mol Pharmacol 80: 796–808. doi: 10.1124/mol.111.072207 [PubMed]
34. Maertens C, Cuypers E, Amininasab M, Jalali A, Vatanpour H, 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]
35. Klint JK, Chin YKY, Mobli M. (2015) Rational engineering defines a molecular switch that Is essential for activity of spider-venom peptides against the analgesics target Na(V)1.7. Mol Pharmacol 88: 1002–1010. doi: 10.1124/mol.115.100784 [PubMed]
36. Sermadiras I, Revell J, Linley JE, Sandercock A, Ravn P (2013) Recombinant expression and in vitro characterisation of active huwentoxin-IV. Plos One 8. [PMC free article] [PubMed]
37. Smith JJ, Cummins TR, Alphy S, Blumenthal KM (2007) Molecular interactions of the gating modifier toxin ProTx-II with NaV1.5: implied existence of a novel toxin binding site coupled to activation. J Biol Chem 282: 12687–12697. doi: 10.1074/jbc.M610462200 [PubMed]
38. Patel RR, Barbosa C, Brustovetsky T, Brustovetsky N, Cummins TR (2016) Aberrant epilepsy-associated mutant NaV1.6 sodium channel activity can be targeted with cannabidiol. Brain 139: 2164–2181. doi: 10.1093/brain/aww129 [PMC free article] [PubMed]
39. Lau CH, King GF, Mobli M (2016) Molecular basis of the interaction between gating modifier spider toxins and the voltage sensor of voltage-gated ion channels. Sci Rep 6. [PMC free article] [PubMed]
40. 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 [PMC free article] [PubMed]
41. 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. [PubMed]
42. Munro G (2009) Pharmacological assessment of the rat formalin test utilizing the clinically used analgesic drugs gabapentin, lamotrigine, morphine, duloxetine, tramadol and ibuprofen: Influence of low and high formalin concentrations. Eur J Pharmacol 605: 95–102. doi: 10.1016/j.ejphar.2009.01.004 [PubMed]
43. Rode F, Thomsen M, Brolos T, Jensen DG, Blackburn-Munro G, Bjerrum OJ (2007) The importance of genetic background on pain behaviours and pharmacological sensitivity in the rat spared serve injury model of peripheral neuropathic pain. Eur J Pharmacol 564: 103–111. doi: 10.1016/j.ejphar.2007.02.009 [PubMed]
44. Sevostianova N, Danysz W, Bespalov AY (2005) Analgesic effects of morphine and loperamide in the rat formalin test: Interactions with NMDA receptor antagonists. Eur J Pharmacol 525: 83–90. doi: 10.1016/j.ejphar.2005.10.010 [PubMed]
45. Sevostianova N, Zvartau E, Bespalov A, Danysz W (2003) Effects of morphine on formalin-induced nociception in rats. Eur J Pharmacol 462: 109–113. [PubMed]
46. Minett MS, Pereira V, Sikandar S, Matsuyama A, Lolignier S, Kanellopoulos AH, et al. (2015) Endogenous opioids contribute to insensitivity to pain in humans and mice lacking sodium channel NaV1.7. Nat Commun 6: 8967 doi: 10.1038/ncomms9967 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science