Clockwise domain arrangement of the sodium channel revealed by (mu)-conotoxin (GIIIA) docking orientation

Functional effects of drugs and toxins interacting with NaV1.4

NaV1.4 is a voltage-gated sodium channel subtype that is predominantly expressed in skeletal muscle cells. It is essential for producing action potentials and stimulating muscle contraction, and mutations in NaV1.4 can cause various muscle disorders. The discovery of the cryo-EM structure of NaV1.4 in complex with β1 has opened new possibilities for designing drugs and toxins that target NaV1.4. In this review, we summarize the current understanding of channelopathies, the binding sites and functions of chemicals including medicine and toxins that interact with NaV1.4. These substances could be considered novel candidate compounds or tools to develop more potent and selective drugs targeting NaV1.4. Therefore, studying NaV1.4 pharmacology is both theoretically and practically meaningful.

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

Hydropathicity-based prediction of pain-causing NaV1.7 variants

Background

Mutation-induced variations in the functional architecture of the NaV1.7 channel protein are causally related to a broad spectrum of human pain disorders. Predicting in silico the phenotype of NaV1.7 variant is of major clinical importance; it can aid in reducing costs of in vitro pathophysiological characterization of NaV1.7 variants, as well as, in the design of drug agents for counteracting pain-disease symptoms.

Results

In this work, we utilize spatial complexity of hydropathic effects toward predicting which NaV1.7 variants cause pain (and which are neutral) based on the location of corresponding mutation sites within the NaV1.7 structure. For that, we analyze topological and scaling hydropathic characteristics of the atomic environment around NaV1.7’s pore and probe their spatial correlation with mutation sites. We show that pain-related mutation sites occupy structural locations in proximity to a hydrophobic patch lining the pore while clustering at a critical hydropathic-interactions distance from the selectivity filter (SF). Taken together, these observations can differentiate pain-related NaV1.7 variants from neutral ones, i.e., NaV1.7 variants not causing pain disease, with 80.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document}% sensitivity and 93.7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document}% specificity [area under the receiver operating characteristics curve = 0.872].

Conclusions

Our findings suggest that maintaining hydrophobic NaV1.7 interior intact, as well as, a finely-tuned (dictated by hydropathic interactions) distance from the SF might be necessary molecular conditions for physiological NaV1.7 functioning. The main advantage for using the presented predictive scheme is its negligible computational cost, as well as, hydropathicity-based biophysical rationalization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04119-2.

Discovery, Pharmacological Characterisation and NMR Structure of the Novel µ-Conotoxin SxIIIC, a Potent and Irreversible NaV Channel Inhibitor

Voltage-gated sodium (Na_V) channel subtypes, including Na_V1.7, are promising targets for the treatment of neurological diseases, such as chronic pain. Cone snail-derived µ-conotoxins are small, potent Na_V channel inhibitors which represent potential drug leads. Of the 22 µ-conotoxins characterised so far, only a small number, including KIIIA and CnIIIC, have shown inhibition against human Na_V1.7. We have recently identified a novel µ-conotoxin, SxIIIC, from Conus striolatus. Here we present the isolation of native peptide, chemical synthesis, characterisation of human Na_V channel activity by whole-cell patch-clamp electrophysiology and analysis of the NMR solution structure. SxIIIC displays a unique Na_V channel selectivity profile (1.4 > 1.3 > 1.1 ≈ 1.6 ≈ 1.7 > 1.2 >> 1.5 ≈ 1.8) when compared to other µ-conotoxins and represents one of the most potent human Na_V1.7 putative pore blockers (IC₅₀ 152.2 ± 21.8 nM) to date. NMR analysis reveals the structure of SxIIIC includes the characteristic α-helix seen in other µ-conotoxins. Future investigations into structure-activity relationships of SxIIIC are expected to provide insights into residues important for Na_V channel pore blocker selectivity and subsequently important for chronic pain drug development.

NMR Structure of μ-Conotoxin GIIIC: Leucine 18 Induces Local Repacking of the N-Terminus Resulting in Reduced NaV Channel Potency

μ-Conotoxins are potent and highly specific peptide blockers of voltage-gated sodium channels. In this study, the solution structure of μ-conotoxin GIIIC was determined using 2D NMR spectroscopy and simulated annealing calculations. Despite high sequence similarity, GIIIC adopts a three-dimensional structure that differs from the previously observed conformation of μ-conotoxins GIIIA and GIIIB due to the presence of a bulky, non-polar leucine residue at position 18. The side chain of L18 is oriented towards the core of the molecule and consequently the N-terminus is re-modeled and located closer to L18. The functional characterization of GIIIC defines it as a canonical μ-conotoxin that displays substantial selectivity towards skeletal muscle sodium channels (Na_V), albeit with ~2.5-fold lower potency than GIIIA. GIIIC exhibited a lower potency of inhibition of Na_V1.4 channels, but the same Na_V selectivity profile when compared to GIIIA. These observations suggest that single amino acid differences that significantly affect the structure of the peptide do in fact alter its functional properties. Our work highlights the importance of structural factors, beyond the disulfide pattern and electrostatic interactions, in the understanding of the functional properties of bioactive peptides. The latter thus needs to be considered when designing analogues for further applications.

Predicting Structural Details of the Sodium Channel Pore Basing on Animal Toxin Studies

Eukaryotic voltage-gated sodium channels play key roles in physiology and are targets for many toxins and medically important drugs. Physiology, pharmacology, and general architecture of the channels has long been the subject of intensive research in academia and industry. In particular, animal toxins such as tetrodotoxin, saxitoxin, and conotoxins have been used as molecular probes of the channel structure. More recently, X-ray structures of potassium and prokaryotic sodium channels allowed elaborating models of the toxin-channel complexes that integrated data from biophysical, electrophysiological, and mutational studies. Atomic level cryo-EM structures of eukaryotic sodium channels, which became available in 2017, show that the selectivity filter structure and other important features of the pore domain have been correctly predicted. This validates further employments of toxins and other small molecules as sensitive probes of fine structural details of ion channels.

Reverse pharmacogenomics: carbamazepine normalizes activation and attenuates thermal hyperexcitability of sensory neurons due to Nav 1.7 mutation I234T

BACKGROUND AND PURPOSE

Pharmacotherapy for pain currently involves trial and error. A previous study on inherited erythromelalgia (a genetic model of neuropathic pain due to mutations in the sodium channel, Nav 1.7) used genomics, structural modelling and biophysical and pharmacological analyses to guide pharmacotherapy and showed that carbamazepine normalizes voltage dependence of activation of the Nav 1.7-S241T mutant channel, reducing pain in patients carrying this mutation. However, whether this approach is applicable to other Nav channel mutants is still unknown.

EXPERIMENTAL APPROACH

We used structural modelling, patch clamp and multi-electrode array (MEA) recording to assess the effects of carbamazepine on Nav 1.7-I234T mutant channels and on the firing of dorsal root ganglion (DRG) sensory neurons expressing these mutant channels.

KEY RESULTS

In a reverse engineering approach, structural modelling showed that the I234T mutation is located in atomic proximity to the carbamazepine-responsive S241T mutation and that activation of Nav 1.7-I234T mutant channels, from patients who are known to respond to carbamazepine, is partly normalized with a clinically relevant concentration (30 μM) of carbamazepine. There was significantly higher firing in intact sensory neurons expressing Nav 1.7-I234T channels, compared with neurons expressing the normal channels (Nav 1.7-WT). Pre-incubation with 30 μM carbamazepine also significantly reduced the firing of intact DRG sensory neurons expressing Nav 1.7-I234T channels. Although the expected use-dependent inhibition of Nav 1.7-WT channels by carbamazepine was confirmed, carbamazepine did not enhance use-dependent inhibition of Nav 1.7-I234T mutant channels.

CONCLUSION AND IMPLICATIONS

These results support the utility of a pharmacogenomic approach to treatment of pain in patients carrying sodium channel variants.

LINKED ARTICLES

This article is part of a themed section on Recent Advances in Targeting Ion Channels to Treat Chronic Pain. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.12/issuetoc.

Guanidinium Toxins and Their Interactions with Voltage-Gated Sodium Ion Channels

Guanidinium toxins, such as saxitoxin (STX), tetrodotoxin (TTX) and their analogs, are naturally occurring alkaloids with divergent evolutionary origins and biogeographical distribution, but which share the common chemical feature of guanidinium moieties. These guanidinium groups confer high biological activity with high affinity and ion flux blockage capacity for voltage-gated sodium channels (NaV). Members of the STX group, known collectively as paralytic shellfish toxins (PSTs), are produced among three genera of marine dinoflagellates and about a dozen genera of primarily freshwater or brackish water cyanobacteria. In contrast, toxins of the TTX group occur mainly in macrozoa, particularly among puffer fish, several species of marine invertebrates and a few terrestrial amphibians. In the case of TTX and analogs, most evidence suggests that symbiotic bacteria are the origin of the toxins, although endogenous biosynthesis independent from bacteria has not been excluded. The evolutionary origin of the biosynthetic genes for STX and analogs in dinoflagellates and cyanobacteria remains elusive. These highly potent molecules have been the subject of intensive research since the latter half of the past century; first to study the mode of action of their toxigenicity, and later as tools to characterize the role and structure of NaV channels, and finally as therapeutics. Their pharmacological activities have provided encouragement for their use as therapeutants for ion channel-related pathologies, such as pain control. The functional role in aquatic and terrestrial ecosystems for both groups of toxins is unproven, although plausible mechanisms of ion channel regulation and chemical defense are often invoked. Molecular approaches and the development of improved detection methods will yield deeper understanding of their physiological and ecological roles. This knowledge will facilitate their further biotechnological exploitation and point the way towards development of pharmaceuticals and therapeutic applications.

Nav1.7-A1632G Mutation from a Family with Inherited Erythromelalgia: Enhanced Firing of Dorsal Root Ganglia Neurons Evoked by Thermal Stimuli

Voltage-gated sodium channel Nav1.7 is a central player in human pain. Mutations in Nav1.7 produce several pain syndromes, including inherited erythromelalgia (IEM), a disorder in which gain-of-function mutations render dorsal root ganglia (DRG) neurons hyperexcitable. Although patients with IEM suffer from episodes of intense burning pain triggered by warmth, the effects of increased temperature on DRG neurons expressing mutant Nav1.7 channels have not been well documented. Here, using structural modeling, voltage-clamp, current-clamp, and multielectrode array recordings, we have studied a newly identified Nav1.7 mutation, Ala1632Gly, from a multigeneration family with IEM. Structural modeling suggests that Ala1632 is a molecular hinge and that the Ala1632Gly mutation may affect channel gating. Voltage-clamp recordings revealed that the Nav1.7-A1632G mutation hyperpolarizes activation and depolarizes fast-inactivation, both gain-of-function attributes at the channel level. Whole-cell current-clamp recordings demonstrated increased spontaneous firing, lower current threshold, and enhanced evoked firing in rat DRG neurons expressing Nav1.7-A1632G mutant channels. Multielectrode array recordings further revealed that intact rat DRG neurons expressing Nav1.7-A1632G mutant channels are more active than those expressing Nav1.7 WT channels. We also showed that physiologically relevant thermal stimuli markedly increase the mean firing frequencies and the number of active rat DRG neurons expressing Nav1.7-A1632G mutant channels, whereas the same thermal stimuli only increase these parameters slightly in rat DRG neurons expressing Nav1.7 WT channels. The response of DRG neurons expressing Nav1.7-A1632G mutant channels upon increase in temperature suggests a cellular basis for warmth-triggered pain in IEM.

SIGNIFICANCE STATEMENT

Inherited erythromelalgia (IEM), a severe pain syndrome characterized by episodes of intense burning pain triggered by warmth, is caused by mutations in sodium channel Nav1.7, which are preferentially expressed in sensory and sympathetic neurons. More than 20 gain-of-function Nav1.7 mutations have been identified from IEM patients, but the question of how warmth triggers episodes of pain in IEM has not been well addressed. Combining multielectrode array, voltage-clamp, and current-clamp recordings, we assessed a newly identified IEM mutation (Nav1.7-A1632G) from a multigeneration family. Our data demonstrate gain-of-function attributes at the channel level and differential effects of physiologically relevant thermal stimuli on the excitability of DRG neurons expressing mutant and WT Nav1.7 channels, suggesting a cellular mechanism for warmth-triggered pain episodes in IEM patients.

Network topology of NaV1.7 mutations in sodium channel-related painful disorders

BACKGROUND

Gain-of-function mutations in SCN9A gene that encodes the voltage-gated sodium channel NaV1.7 have been associated with a wide spectrum of painful syndromes in humans including inherited erythromelalgia, paroxysmal extreme pain disorder and small fibre neuropathy. These mutations change the biophysical properties of NaV1.7 channels leading to hyperexcitability of dorsal root ganglion nociceptors and pain symptoms. There is a need for better understanding of how gain-of-function mutations alter the atomic structure of Nav1.7.

RESULTS

We used homology modeling to build an atomic model of NaV1.7 and a network-based theoretical approach, which can predict interatomic interactions and connectivity arrangements, to investigate how pain-related NaV1.7 mutations may alter specific interatomic bonds and cause connectivity rearrangement, compared to benign variants and polymorphisms. For each amino acid substitution, we calculated the topological parameters betweenness centrality (B ct ), degree (D), clustering coefficient (CC ct ), closeness (C ct ), and eccentricity (E ct ), and calculated their variation (Δ value  = mutant value -WT value ). Pathogenic NaV1.7 mutations showed significantly higher variation of |ΔB ct | compared to benign variants and polymorphisms. Using the cut-off value ±0.26 calculated by receiver operating curve analysis, we found that ΔB ct correctly differentiated pathogenic NaV1.7 mutations from variants not causing biophysical abnormalities (nABN) and homologous SNPs (hSNPs) with 76% sensitivity and 83% specificity.

CONCLUSIONS

Our in-silico analyses predict that pain-related pathogenic NaV1.7 mutations may affect the network topological properties of the protein and suggest |ΔB ct | value as a potential in-silico marker.

Sodium Channels and Venom Peptide Pharmacology

Venomous animals including cone snails, spiders, scorpions, anemones, and snakes have evolved a myriad of components in their venoms that target the opening and/or closing of voltage-gated sodium channels to cause devastating effects on the neuromuscular systems of predators and prey. These venom peptides, through design and serendipity, have not only contributed significantly to our understanding of sodium channel pharmacology and structure, but they also represent some of the most phyla- and isoform-selective molecules that are useful as valuable tool compounds and drug leads. Here, we review our understanding of the basic function of mammalian voltage-gated sodium channel isoforms as well as the pharmacology of venom peptides that act at these key transmembrane proteins.

The hitchhiker's guide to the voltage-gated sodium channel galaxy

Eukaryotic voltage-gated sodium (Na_v) channels contribute to the rising phase of action potentials and served as an early muse for biophysicists laying the foundation for our current understanding of electrical signaling. Given their central role in electrical excitability, it is not surprising that (a) inherited mutations in genes encoding for Na_v channels and their accessory subunits have been linked to excitability disorders in brain, muscle, and heart; and (b) Na_v channels are targeted by various drugs and naturally occurring toxins. Although the overall architecture and behavior of these channels are likely to be similar to the more well-studied voltage-gated potassium channels, eukaryotic Na_v channels lack structural and functional symmetry, a notable difference that has implications for gating and selectivity. Activation of voltage-sensing modules of the first three domains in Na_v channels is sufficient to open the channel pore, whereas movement of the domain IV voltage sensor is correlated with inactivation. Also, structure–function studies of eukaryotic Na_v channels show that a set of amino acids in the selectivity filter, referred to as DEKA locus, is essential for Na⁺ selectivity. Structures of prokaryotic Na_v channels have also shed new light on mechanisms of drug block. These structures exhibit lateral fenestrations that are large enough to allow drugs or lipophilic molecules to gain access into the inner vestibule, suggesting that this might be the passage for drug entry into a closed channel. In this Review, we will synthesize our current understanding of Na_v channel gating mechanisms, ion selectivity and permeation, and modulation by therapeutics and toxins in light of the new structures of the prokaryotic Na_v channels that, for the time being, serve as structural models of their eukaryotic counterparts.

Venom Peptides From Cone Snails: Pharmacological Probes for Voltage-Gated Sodium Channels

The venoms of cone snails provide a rich source of neuroactive peptides (conotoxins). Several venom peptide families have been identified that are either agonists (ι- and δ-conotoxins) or antagonists (μ- and μO-conotoxins) of voltage-gated sodium channels (VGSCs). Members of these conotoxin classes have been integral in identifying and characterizing specific neurotoxin binding sites on the channel. Furthermore, given the specificity of some of these peptides for one sodium channel subtype over another, conotoxins have also proven useful in exploring differences between VGSC subtypes. This chapter summarizes the current knowledge of the structure and function based on the results of conotoxin interactions with VGSCs and correlates the peptides with the phylogeny of the Conus species from which they were derived.

An autism-associated mutation in CaV1.3 channels has opposing effects on voltage- and Ca(2+)-dependent regulation

CaV1.3 channels are a major class of L-type Ca(2+) channels which contribute to the rhythmicity of the heart and brain. In the brain, these channels are vital for excitation-transcription coupling, synaptic plasticity, and neuronal firing. Moreover, disruption of CaV1.3 function has been associated with several neurological disorders. Here, we focus on the de novo missense mutation A760G which has been linked to autism spectrum disorder (ASD). To explore the role of this mutation in ASD pathogenesis, we examined the effects of A760G on CaV1.3 channel gating and regulation. Introduction of the mutation severely diminished the Ca(2+)-dependent inactivation (CDI) of CaV1.3 channels, an important feedback system required for Ca(2+) homeostasis. This reduction in CDI was observed in two major channel splice variants, though to different extents. Using an allosteric model of channel gating, we found that the underlying mechanism of CDI reduction is likely due to enhanced channel opening within the Ca(2+)-inactivated mode. Remarkably, the A760G mutation also caused an opposite increase in voltage-dependent inactivation (VDI), resulting in a multifaceted mechanism underlying ASD. When combined, these regulatory deficits appear to increase the intracellular Ca(2+) concentration, thus potentially disrupting neuronal development and synapse formation, ultimately leading to ASD.

Structure of the voltage-gated calcium channel Cav1.1 complex

The voltage-gated calcium channel Ca(v)1.1 is engaged in the excitation-contraction coupling of skeletal muscles. The Ca(v)1.1 complex consists of the pore-forming subunit α1 and auxiliary subunits α2δ, β, and γ. We report the structure of the rabbit Ca(v)1.1 complex determined by single-particle cryo-electron microscopy. The four homologous repeats of the α1 subunit are arranged clockwise in the extracellular view. The γ subunit, whose structure resembles claudins, interacts with the voltage-sensing domain of repeat IV (VSD(IV)), whereas the cytosolic β subunit is located adjacent to VSD(II) of α1. The α2 subunit interacts with the extracellular loops of repeats I to III through its VWA and Cache1 domains. The structure reveals the architecture of a prototypical eukaryotic Ca(v) channel and provides a framework for understanding the function and disease mechanisms of Ca(v) and Na(v) channels.

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

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

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

Peptide toxins often have pharmacological applications and are powerful tools for investigating the structure-function relationships of voltage-gated sodium channels (VGSCs). Although a group of potential VGSC inhibitors have been reported from tarantula venoms, little is known about the mechanism of their interaction with VGSCs. In this study, we showed that hainantoxin-IV (β-TRTX-Hn2a, HNTX-IV in brief), a 35-residue peptide from Ornithoctonus hainana venom, preferentially inhibited rNav1.2, rNav1.3 and hNav1.7 compared with rNav1.4 and hNav1.5. hNav1.7 was the most sensitive to HNTX-IV (IC50∼21nM). In contrast to many other tarantula toxins that affect VGSCs, HNTX-IV at subsaturating concentrations did not alter activation and inactivation kinetics in the physiological range of voltages, while very large depolarization above +70mV could partially activate toxin-bound hNav1.7 channel, indicating that HNTX-IV acts as a gating modifier rather than a pore blocker. Site-directed mutagenesis indicated that the toxin bound to site 4, which was located on the extracellular S3-S4 linker of hNav1.7 domain II. Mutants E753Q, D816N and E818Q of hNav1.7 decreased toxin affinity for hNav1.7 by 2.0-, 3.3- and 130-fold, respectively. In silico docking indicated that a three-toed claw substructure formed by residues with close contacts in the interface between HNTX-IV and hNav1.7 domain II stabilized the toxin-channel complex, impeding movement of the domain II voltage sensor and inhibiting hNav1.7 activation. Our data provide structural details for structure-based drug design and a useful template for the design of highly selective inhibitors of a specific subtype of VGSCs.

Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA, PIIIA, and KIIIA in a NavAb-based model of Nav1.4

Analyses of toxin binding to a homology model of Nav1.4 indicate similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Systematic study of binding of μ-conotoxins to the sodium channel NaV1.4

Voltage-gated sodium channels (NaV) are fundamental components of the nervous system. Their dysfunction is implicated in a number of neurological disorders, such as chronic pain, making them potential targets for the treatment of such disorders. The prominence of the NaV channels in the nervous system has been exploited by venomous animals for preying purposes, which have developed toxins that can block the NaV channels, thereby disabling their function. Because of their potency, such toxins could provide drug leads for the treatment of neurological disorders associated with NaV channels. However, most toxins lack selectivity for a given target NaV channel, and improving their selectivity profile among the NaV1 isoforms is essential for their development as drug leads. Computational methods will be very useful in the solution of such design problems, provided accurate models of the protein-ligand complex can be constructed. Using docking and molecular dynamics simulations, we have recently constructed a model for the NaV1.4-μ-conotoxin-GIIIA complex and validated it with the ample mutational data available for this complex. Here, we use the validated NaV1.4 model in a systematic study of binding other μ-conotoxins (PIIIA, KIIIA and BuIIIB) to NaV1.4. The binding mode obtained for each complex is shown to be consistent with the available mutation data and binding constants. We compare the binding modes of PIIIA, KIIIA and BuIIIB to that of GIIIA and point out the similarities and differences among them. The detailed information about NaV1.4-μ-conotoxin interactions provided here will be useful in the design of new NaV channel blocking peptides.

Molecular dynamics study of binding of µ-conotoxin GIIIA to the voltage-gated sodium channel Na(v)1.4

Homology models of mammalian voltage-gated sodium (NaV) channels based on the crystal structures of the bacterial counterparts are needed to interpret the functional data on sodium channels and understand how they operate. Such models would also be invaluable in structure-based design of therapeutics for diseases involving sodium channels such as chronic pain and heart diseases. Here we construct a homology model for the pore domain of the NaV1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to NaV1.4 to validate the model. The initial poses for the NaV1.4-GIIIA complex are obtained using the HADDOCK protein docking program, which are then refined in molecular dynamics simulations. The binding mode for the final complex is shown to be in broad agreement with the available mutagenesis data. The standard binding free energy, determined from the potential of mean force calculations, is also in good agreement with the experimental value. Because the pore domains of NaV1 channels are highly homologous, the model constructed for NaV1.4 will provide an excellent template for other NaV1 channels.

Mechanism of μ-conotoxin PIIIA binding to the voltage-gated Na+ channel NaV1.4

Several subtypes of voltage-gated Na+ (NaV) channels are important targets for pain management. μ-Conotoxins isolated from venoms of cone snails are potent and specific blockers of different NaV channel isoforms. The inhibitory effect of μ-conotoxins on NaV channels has been examined extensively, but the mechanism of toxin specificity has not been understood in detail. Here the known structure of μ-conotoxin PIIIA and a model of the skeletal muscle channel NaV1.4 are used to elucidate elements that contribute to the structural basis of μ-conotoxin binding and specificity. The model of NaV1.4 is constructed based on the crystal structure of the bacterial NaV channel, NaVAb. Six different binding modes, in which the side chain of each of the basic residues carried by the toxin protrudes into the selectivity filter of NaV1.4, are examined in atomic detail using molecular dynamics simulations with explicit solvent. The dissociation constants (Kd) computed for two selected binding modes in which Lys9 or Arg14 from the toxin protrudes into the filter of the channel are within 2 fold; both values in close proximity to those determined from dose response data for the block of NaV currents. To explore the mechanism of PIIIA specificity, a double mutant of NaV1.4 mimicking NaV channels resistant to μ-conotoxins and tetrodotoxin is constructed and the binding of PIIIA to this mutant channel examined. The double mutation causes the affinity of PIIIA to reduce by two orders of magnitude.

Ala-7, His-10 and Arg-12 are crucial amino acids for activity of a synthetically engineered μ-conotoxin

Cone snail toxins or conotoxins are often small cysteine-rich peptides which have shown to be highly selective ligands for a wide range of ion channels such as voltage-gated sodium channels (Na(V)s). Na(V)s participate in a wide range of electrophysiological processes. Consequently, their malfunction has been associated with numerous diseases. The development of subtype-selective modulators of Na(V)s remains highly important in the treatment of such disorders. In order to expand our knowledge in the search for novel therapeutics to treat Na(V)-related diseases, we explored the field of peptide engineering. In the current study, the impact of well considered point mutations into a bioactive peptide that was found to be a very potent and selective inhibitor of Na(V)s (i.e. Midi R2) was examined. We designed two peptides, named Midi R2[A7G] and Midi R2[H10A, R12A] which have mutations at position 7, and both 10 and 12, respectively. Electrophysiological recordings indicated that an Ala to Gly mutation at position 7 increased IC50-values from the nanomolar range to the micromolar range. For Midi R2[H10A, R12A] at a concentration of 10 μM, activity is even reduced to 0-10% for all of the tested Na(V)-channels. Circular dichroism measurements proved that overall structural conformations did not change. These findings suggest that the minimal space between the second and the third intercysteine loop of Midi R2 is the sequence RRWARDHSR and that His at position 10 and Arg at position 12 are crucial amino acids for the potency and specificity of Midi R2. In this way, new insights into the structure-activity relationships of μ-conotoxins were found.

Probing gating mechanisms of sodium channels using pore blockers

Several classes of small molecules and peptides bind at the central pore of voltage-gated sodium channels either from the extracellular or intracellular side of the membrane and block ion conduction through the pore. Biophysical studies that shed light on the chemical nature, accessibility, and kinetics of binding of these naturally occurring and synthetic compounds reveal a wealth of information about how these channels gate. Here, we discuss insights into the structural underpinnings of gating of the channel pore and its coupling to the voltage sensors obtained from pore blockers including site 1 neurotoxins and local anesthetics.

Complex structures between the N-type calcium channel (CaV2.2) and ω-conotoxin GVIA predicted via molecular dynamics

The N-type voltage-gated Ca(2+) channel CaV2.2 is one of the important targets for pain management. ω-Conotoxins isolated from venoms of cone snails, which specifically inhibit CaV2.2, are promising scaffolds for novel analgesics. The inhibitory action of ω-conotoxins on CaV2.2 has been examined experimentally, but the modes of binding of the toxins to this and other related subfamilies of Ca(2+) channels are not understood in detail. Here molecular dynamics simulations are used to construct models of ω-conotoxin GVIA in complex with a homology model of the pore domain of CaV2.2. Three different binding modes in which the side chain of Lys2, Arg17, or Lys24 from the toxin protrudes into the selectivity filter of CaV2.2 are considered. In all the modes, the toxin forms a salt bridge with an aspartate residue of subunit II just above the EEEE ring of the selectivity filter. Using the umbrella sampling technique and potential of mean force calculations, the half-maximal inhibitory concentration (IC50) values are calculated to be 1.5 and 0.7 nM for the modes in which Lys2 and Arg17 occlude the ion conduction pathway, respectively. Both IC50 values compare favorably with the values of 0.04-1.0 nM determined experimentally. The similar IC50 values calculated for the different binding modes demonstrate that GVIA can inhibit CaV2.2 with alternative binding modes. Such a multiple-binding mode mechanism may be common for ω-conotoxins.

Molecular architecture of a sodium channel S6 helix: radial tuning of the voltage-gated sodium channel 1.7 activation gate

BACKGROUND

In-frame deletion mutation (Del-L955) in NaV1.7 sodium channel from a kindred with erythromelalgia hyperpolarizes activation.

RESULTS

Del-L955 twists the S6 helix, displacing the Phe960 activation gate. Replacement of Phe960 at the correct helical position depolarizes activation.

CONCLUSION

Radial tuning of the activation gate is critical to the activation of NaV1.7 channel.

SIGNIFICANCE

Structural modeling guided electrophysiology reveals the functional importance of radial tuning of the S6 segment. Voltage-gated sodium (NaV) channels are membrane proteins that consist of 24 transmembrane segments organized into four homologous domains and are essential for action potential generation and propagation. Although the S6 helices of NaV channels line the ion-conducting pore and participate in channel activation, their functional architecture is incompletely understood. Our recent studies show that a naturally occurring in-frame deletion mutation (Del-L955) of NaV1.7 channel, identified in individuals with a severe inherited pain syndrome (inherited erythromelalgia) causes a substantial hyperpolarizing shift of channel activation. Here we took advantage of this deletion mutation to understand the role of the S6 helix in the channel activation. Based on the recently published structure of a bacterial NaV channel (NaVAb), we modeled the WT and Del-L955 channel. Our structural model showed that Del-L955 twists the DII/S6 helix, shifting location and radial orientation of the activation gate residue (Phe(960)). Hypothesizing that these structural changes produce the shift of channel activation of Del-L955 channels, we restored a phenylalanine in wild-type orientation by mutating Ser(961) (Del-L955/S961F), correcting activation by ∼10 mV. Correction of the displaced Phe(960) (F960S) together with introduction of the rescuing activation gate residue (S961F) produced an additional ∼6-mV restoration of activation of the mutant channel. A simple point mutation in the absence of a twist (L955A) did not produce a radial shift and did not hyperpolarize activation. Our results demonstrate the functional importance of radial tuning of the sodium channel S6 helix for the channel activation.

Structural modelling and mutant cycle analysis predict pharmacoresponsiveness of a Na(V)1.7 mutant channel

Sodium channel Na(V)1.7 is critical for human pain signalling. Gain-of-function mutations produce pain syndromes including inherited erythromelalgia, which is usually resistant to pharmacotherapy, but carbamazepine normalizes activation of Na(V)1.7-V400M mutant channels from a family with carbamazepine-responsive inherited erythromelalgia. Here we show that structural modelling and thermodynamic analysis predict pharmacoresponsiveness of another mutant channel (S241T) that is located 159 amino acids distant from V400M. Structural modelling reveals that Na(v)1.7-S241T is ~2.4 Å apart from V400M in the folded channel, and thermodynamic analysis demonstrates energetic coupling of V400M and S241T during activation. Atomic proximity and energetic coupling are paralleled by pharmacological coupling, as carbamazepine (30 μM) depolarizes S214T activation, as previously reported for V400M. Pharmacoresponsiveness of S241T to carbamazepine was further evident at a cellular level, where carbamazepine normalized the hyperexcitability of dorsal root ganglion neurons expressing S241T. We suggest that this approach might identify variants that confer enhanced pharmacoresponsiveness on a variety of channels.

Computational studies of marine toxins targeting ion channels

Toxins from marine animals offer novel drug leads for treatment of diseases involving ion channels. Computational methods could be very helpful in this endeavour in several ways, e.g., (i) constructing accurate models of the channel-toxin complexes using docking and molecular dynamics (MD) simulations; (ii) determining the binding free energies of toxins from umbrella sampling MD simulations; (iii) predicting the effect of mutations from free energy MD simulations. Using these methods, one can design new analogs of toxins with improved affinity and selectivity properties. Here we present a review of the computational methods and discuss their applications to marine toxins targeting potassium and sodium channels. Detailed examples from the potassium channel toxins—ShK from sea anemone and κ-conotoxin PVIIA—are provided to demonstrate capabilities of the computational methods to give accurate descriptions of the channel-toxin complexes and the energetics of their binding. An example is also given from sodium channel toxins (µ-conotoxin GIIIA) to illustrate the differences between the toxin binding modes in potassium and sodium channels.

Mechanism and molecular basis for the sodium channel subtype specificity of µ-conopeptide CnIIIC

BACKGROUND AND PURPOSE

Voltage-gated sodium channels (Na(V) channels) are key players in the generation and propagation of action potentials, and selective blockade of these channels is a promising strategy for clinically useful suppression of electrical activity. The conotoxin µ-CnIIIC from the cone snail Conus consors exhibits myorelaxing activity in rodents through specific blockade of skeletal muscle (Na(V) 1.4) Na(V) channels.

EXPERIMENTAL APPROACH

We investigated the activity of µ-CnIIIC on human Na(V) channels and characterized its inhibitory mechanism, as well as the molecular basis, for its channel specificity.

KEY RESULTS

Similar to rat paralogs, human Na(V) 1.4 and Na(V) 1.2 were potently blocked by µ-CnIIIC, the sensitivity of Na(V) 1.7 was intermediate, and Na(V) 1.5 and Na(V) 1.8 were insensitive. Half-channel chimeras revealed that determinants for the insensitivity of Na(V) 1.8 must reside in both the first and second halves of the channel, while those for Na(V) 1.5 are restricted to domains I and II. Furthermore, domain I pore loop affected the total block and therefore harbours the major determinants for the subtype specificity. Domain II pore loop only affected the kinetics of toxin binding and dissociation. Blockade by µ-CnIIIC of Na(V) 1.4 was virtually irreversible but left a residual current of about 5%, reflecting a 'leaky' block; therefore, Na(+) ions still passed through µ-CnIIIC-occupied Na(V) 1.4 to some extent. TTX was excluded from this binding site but was trapped inside the pore by µ-CnIIIC.

CONCLUSION AND IMPLICATIONS

Of clinical significance, µ-CnIIIC is a potent and persistent blocker of human skeletal muscle Na(V) 1.4 that does not affect activity of cardiac Na(V) 1.5.

Three-dimensional localization of the α and β subunits and of the II-III loop in the skeletal muscle L-type Ca2+ channel

The L-type Ca(2+) channel (dihydropyridine receptor (DHPR) in skeletal muscle acts as the voltage sensor for excitation-contraction coupling. To better resolve the spatial organization of the DHPR subunits (α(1s) or Ca(V)1.1, α(2), β(1a), δ1, and γ), we created transgenic mice expressing a recombinant β(1a) subunit with YFP and a biotin acceptor domain attached to its N- and C- termini, respectively. DHPR complexes were purified from skeletal muscle, negatively stained, imaged by electron microscopy, and subjected to single-particle image analysis. The resulting 19.1-Å resolution, three-dimensional reconstruction shows a main body of 17 × 11 × 8 nm with five corners along its perimeter. Two protrusions emerge from either face of the main body: the larger one attributed to the α(2)-δ1 subunit that forms a flexible hook-shaped feature and a smaller protrusion on the opposite side that corresponds to the II-III loop of Ca(V)1.1 as revealed by antibody labeling. Novel features discernible in the electron density accommodate the atomic coordinates of a voltage-gated sodium channel and of the β subunit in a single docking possibility that defines the α1-β interaction. The β subunit appears more closely associated to the membrane than expected, which may better account for both its role in localizing the α(1s) subunit to the membrane and its suggested role in excitation-contraction coupling.

Design of bioactive peptides from naturally occurring μ-conotoxin structures

To date, cone snail toxins ("conotoxins") are of great interest in the pursuit of novel subtype-selective modulators of voltage-gated sodium channels (Na(v)s). Na(v)s participate in a wide range of electrophysiological processes. Consequently, their malfunctioning has been associated with numerous diseases. The development of subtype-selective modulators of Na(v)s remains highly important in the treatment of such disorders. In current research, a series of novel, synthetic, and bioactive compounds were designed based on two naturally occurring μ-conotoxins that target Na(v)s. The initial designed peptide contains solely 13 amino acids and was therefore named "Mini peptide." It was derived from the μ-conotoxins KIIIA and BuIIIC. Based on this Mini peptide, 10 analogues were subsequently developed, comprising 12-16 amino acids with two disulfide bridges. Following appropriate folding and mass verification, blocking effects on Na(v)s were investigated. The most promising compound established an IC(50) of 34.1 ± 0.01 nM (R2-Midi on Na(v)1.2). An NMR structure of one of our most promising compounds was determined. Surprisingly, this structure does not reveal an α-helix. We prove that it is possible to design small peptides based on known pharmacophores of μ-conotoxins without losing their potency and selectivity. These data can provide crucial material for further development of conotoxin-based therapeutics.

Importance of position 8 in μ-conotoxin KIIIA for voltage-gated sodium channel selectivity

μ-Conotoxin KIIIA from Conus kinoshitai is a 16-residue peptide that acts as a potent pore blocker of several voltage-gated sodium channels (Na(v)). In order to obtain more selective blockers and to investigate the role of Trp at position 8, we substituted this residue with Arg, Gln and Glu. KIIIA and analogues were tested on a range of Na(v) expressed in Xenopus laevis oocytes. The rank order of potency for KIIIA was: rNa(v)1.4 ≥ rNa(v)1.2 > mNa(v)1.6 > rNa(v)1.3, with IC(50) values of 48 ± 6 nm, 61 ± 5 nm, 183 ± 31 nm and 3.6 ± 0.3 μm, respectively, whereas no effect was seen on hNa(v)1.5 and hNa(v)1.8 at a concentration of 10 μm. Replacement of Trp8 resulted in more selective blockers with a preference for neuronal sodium channels over the skeletal sodium channel. The activity on rNa(v)1.4 was reduced about 40-, 70- and 200-fold for [W8R]KIIIA, [W8Q]KIIIA and [W8E]KIIIA, respectively. All analogues showed a completely reversible block of rNa(v)1.2, as opposed to the partial reversibility of KIIIA. At saturating concentrations, complete block of rNa(v)1.2 was never achieved. The residual current was lower than 10%, except for [W8E]KIIIA. KIIIA had no effect on the voltage dependence of activation of rNa(v)1.2, whereas all analogues caused a depolarizing shift. Overall, this study shows that Trp8 is a key residue in the pharmacophore. Replacement of Trp8 enables more selective blockers to be obtained for neuronal sodium channels. Trp is a key determinant for the reversibility of block of rNa(v)1.2.

Molecular determinants for the subtype specificity of μ-conotoxin SIIIA targeting neuronal voltage-gated sodium channels

Voltage-gated sodium channels (Na(V) channels) play a pivotal role in neuronal excitability; they are specifically targeted by μ-conotoxins from the venom of marine cone snails. These peptide toxins bind to the outer vestibule of the channel pore thereby blocking ion conduction through Na(V) channels. μ-Conotoxin SIIIA from Conus striatus was shown to be a potent inhibitor of neuronal sodium channels and to display analgesic effects in mice, albeit the molecular targets are not unambiguously known. We therefore studied recombinant Na(V) channels expressed in mammalian cells using the whole-cell patch-clamp method. Synthetic μSIIIA slowly and partially blocked rat Na(V)1.4 channels with an apparent IC(50) of 0.56 ± 0.29 μM; the block was not complete, leaving at high concentration a residual current component of about 10% with a correspondingly reduced single-channel conductance. At 10 μM, μSIIIA potently blocked rat Na(V)1.2, rat and human Na(V)1.4, and mouse Na(V)1.6 channels; human Na(V)1.7 channels were only inhibited by 58.1 ± 4.9%, whereas human Na(V)1.5 as well as rat and human Na(V)1.8 were insensitive. Employing domain chimeras between rNa(V)1.4 and hNa(V)1.5, we located the determinants for μSIIIA specificity in the first half of the channel protein with a major contribution of domain-2 and a minor contribution of domain-1. The latter was largely accounted for by the alteration in the TTX-binding site (Tyr401 in rNa(V)1.4, Cys for Na(V)1.5, and Ser for Na(V)1.8). Introduction of domain-2 pore loops of all tested channel isoforms into rNa(V)1.4 conferred the μSIIIA phenotype of the respective donor channels highlighting the importance of the domain-2 pore loop as the major determinant for μSIIIA's subtype specificity. Single-site substitutions identified residue Ala728 in rNa(V)1.4 as crucial for its high sensitivity toward μSIIIA. Likewise, Asn889 at the homologous position in hNa(V)1.7 is responsible for the channel's reduced μSIIIA sensitivity. These results will pave the way for the rational design of selective Na(V)-channel antagonists for research and medical applications.

The voltage dependence of gating currents of the neuronal CA(v)3.3 channel is determined by the gating brake in the I-II loop

Low-voltage-activated Ca(v)3 Ca(2+) channels have an activation threshold around -60 mV, which is lower than the activation threshold of other voltage-dependent calcium channels (VDCCs). The kinetics of their activation at membrane voltages just above the activation threshold is much slower than the activation kinetics of other VDCCs. It was demonstrated recently that the intracellular loop connecting repeats I and II of all three Ca(v)3 channels contains a so-called gating brake. Disruption of this brake yields channels that activate at even more hyperpolarized potentials with significantly accelerated kinetics. We have compared gating of a wild-type Ca(v)3.3 channel and a mutated ID12 channel, in which the putative gating brake at the proximal part of the I-II loop was removed. Voltage dependence of the gating current activation was shifted by 34.6 mV towards more hyperpolarized potentials in ID12 channel. ON-charge movement was significantly faster in the ID12 channel, while the kinetics of the off-charge was not altered by the mutation. We conclude that the putative gating brake in I-II loop hinders not only the opening of the conducting pore but also the activating movement of voltage-sensing S4 segments, stabilizing the channel in its closed state.

Crystal structure of a potassium ion transporter, TrkH

The TrkH/TrkG/KtrB proteins mediate K⁺ uptake in bacteria and likely evolved from simple K⁺ channels by multiple gene duplications or fusions. Here we present the crystal structure of a TrkH from Vibrio parahaemolyticus. TrkH is a homodimer, and each protomer contains an ion permeation pathway. A selectivity filter, similar in architecture to those of K⁺ channels but significantly shorter, is lined by backbone and side chain oxygen atoms. Functional studies showed that the TrkH allows permeation of K⁺ and Rb⁺ but not smaller ions such as Na⁺ or Li⁺. Immediately intracellular to the selectivity filter are an intramembrane loop and an arginine residue, both highly conserved, which constrict the permeation pathway. Substituting the arginine with an alanine significantly increases the rate of K⁺ flux. These results reveal the molecular basis of K⁺ selectivity and suggest a novel gating mechanism by this large and important family of membrane transport proteins.

Evidence for lateral mobility of voltage sensors in prokaryotic voltage-gated sodium channels

Voltage-sensor domains (VSDs) in voltage-gated ion channels are thought to regulate the probability that a channel adopts an open conformation by moving vertically in the lipid bilayer. Here we characterized the movement of the VSDs of the prokaryotic voltage-gated sodium channel, NaChBac. Substitution of residue T110, which is located on the extracellular side of the fourth transmembrane helix of the VSD, by cysteine resulted in the formation of a disulfide bond between adjacent subunits in the channel. Our results suggest that T110 residues in VSDs of adjacent subunits can come into close proximity, implying that the VSDs can move laterally in the membrane and constitute a mechanism that regulates channel activity.

Cooccupancy of the outer vestibule of voltage-gated sodium channels by micro-conotoxin KIIIA and saxitoxin or tetrodotoxin

The guanidinium alkaloids tetrodotoxin (TTX) and saxitoxin (STX) are classic ligands of voltage-gated sodium channels (VGSCs). Like TTX and STX, micro-conotoxin peptides are pore blockers but with greater VGSC subtype selectivity. micro-Conotoxin KIIIA blocks the neuronal subtype Na(V)1.2 with nanomolar affinity and we recently discovered that KIIIA and its mutant with one fewer positive charge, KIIIA[K7A], could act synergistically with TTX in a ternary peptide x TTX x Na(V) complex. In the complex, the peptide appeared to trap TTX in its normal binding site such that TTX could not readily dissociate from the channel until the peptide had done so; in turn, the presence of TTX accelerated the rate at which peptide dissociated from the channel. In the present study we examined the inhibition of Na(V)1.2, exogenously expressed in Xenopus oocytes, by STX (a divalent cation) and its sulfated congener GTX2/3 (with a net +1 charge). Each could form a ternary complex with KIIIA and Na(V)1.2, as previously found with TTX (a monovalent cation), but only when STX or GTX2/3 was added before KIIIA. The KIIIA x alkaloid x Na(V) complex was considerably less stable with STX than with either GTX2/3 or TTX. In contrast, ternary KIIIA[K7A] x alkaloid x Na(V) complexes could be formed with either order of ligand addition and were about equally stable with STX, GTX2/3, or TTX. The most parsimonious interpretation of the overall results is that the alkaloid and peptide are closely apposed in the ternary complex. The demonstration that two interacting ligands ("syntoxins") occupy adjacent sites raises the possibility of evolving a much more sophisticated neuropharmacology of VGSCs.

The outer vestibule of the Na+ channel-toxin receptor and modulator of permeation as well as gating

The outer vestibule of voltage-gated Na(+) channels is formed by extracellular loops connecting the S5 and S6 segments of all four domains ("P-loops"), which fold back into the membrane. Classically, this structure has been implicated in the control of ion permeation and in toxin blockage. However, conformational changes of the outer vestibule may also result in alterations in gating, as suggested by several P-loop mutations that gave rise to gating changes. Moreover, partial pore block by mutated toxins may reverse gating changes induced by mutations. Therefore, toxins that bind to the outer vestibule can be used to modulate channel gating.

Molecular endpoints of Ca2+/calmodulin- and voltage-dependent inactivation of Ca(v)1.3 channels

Ca²⁺/calmodulin- and voltage-dependent inactivation (CDI and VDI) comprise vital prototypes of Ca²⁺ channel modulation, rich with biological consequences. Although the events initiating CDI and VDI are known, their downstream mechanisms have eluded consensus. Competing proposals include hinged-lid occlusion of channels, selectivity filter collapse, and allosteric inhibition of the activation gate. Here, novel theory predicts that perturbations of channel activation should alter inactivation in distinctive ways, depending on which hypothesis holds true. Thus, we systematically mutate the activation gate, formed by all S6 segments within Ca_V1.3. These channels feature robust baseline CDI, and the resulting mutant library exhibits significant diversity of activation, CDI, and VDI. For CDI, a clear and previously unreported pattern emerges: activation-enhancing mutations proportionately weaken inactivation. This outcome substantiates an allosteric CDI mechanism. For VDI, the data implicate a “hinged lid–shield” mechanism, similar to a hinged-lid process, with a previously unrecognized feature. Namely, we detect a “shield” in Ca_V1.3 channels that is specialized to repel lid closure. These findings reveal long-sought downstream mechanisms of inactivation and may furnish a framework for the understanding of Ca²⁺ channelopathies involving S6 mutations.

The tetrodotoxin binding site is within the outer vestibule of the sodium channel

Tetrodotoxin and saxitoxin are small, compact asymmetrical marine toxins that block voltage-gated Na channels with high affinity and specificity. They enter the channel pore’s outer vestibule and bind to multiple residues that control permeation. Radiolabeled toxins were key contributors to channel protein purification and subsequent cloning. They also helped identify critical structural elements called P loops. Spacial organization of their mutation-identified interaction sites in molecular models has generated a molecular image of the TTX binding site in the outer vestibule and the critical permeation and selectivity features of this region. One site in the channel’s domain I P loop determines affinity differences in mammalian isoforms.

Probing the bradycardic drug binding receptor of HCN-encoded pacemaker channels

I_f (or I_h), encoded by the hyperpolarization-activated, cyclic nucleotide-gated (HCN1–4) channel gene family, contributes significantly to cardiac pacing. Bradycardic agents such as ZD7288 that target HCN channels have been developed, but the molecular configuration of their receptor is poorly defined. Here, we probed the drug receptor by systematically introducing alanine scanning substitutions into the selectivity filter (C347A, I348A, G349A, Y350A, G351A in the P-loop), outer (P355A, V356A, S357A, M358A in the P-S6 linker), and inner (M377A, F378A, V379A in S6) pore vestibules of HCN1 channels. When heterologously expressed in human embryonic kidney 293 cells for patch-clamp recordings, I348A, G349A, Y350A, G351A, P355A, and V356A did not produce measurable currents. The half-blocking concentration (IC₅₀) of wild type (WT) for ZD7288 was 25.8 ± 9.7 μM. While the IC₅₀ of M358A was identical to WT, those of C347A, S357A, F378A, and V379A markedly increased to 137.6 ± 56.4, 113.3 ± 34.1, 587.1 ± 167.5, and 1726.3 ± 673.4 μM, respectively (p < 0.05). Despite the proximity of the S6 residues studied, M377A was hypersensitive (IC₅₀ = 5.1 ± 0.7 μM; p < 0.05) implicating site specificity. To explore the energetic interactions among the S6 residues, double and triple substitutions (M377A/F378A, M377A/V379A, F378A/V379A, and M377A/F378A/V379A) were generated for thermodynamic cycle analysis. Specific interactions with coupling energies (ΔΔG) >1 kT for M377–F378 and F378–V379 but not M377–V379 were identified. Based on these new data and others, we proposed a refined drug-blocking model that may lead to improved antiarrhythmics and bioartificial pacemaker designs.

Structural modeling of calcium binding in the selectivity filter of the L-type calcium channel

Calcium channels play crucial physiological roles. In the absence of high-resolution structures of the channels, the mechanism of ion permeation is unknown. Here we used a method proposed in an accompanying paper (Cheng and Zhorov in Eur Biophys J, 2009) to predict possible chelation patterns of calcium ions in a structural model of the L-type calcium channel. We compared three models in which two or three calcium ions interact with the four selectivity filter glutamates and a conserved aspartate adjacent to the glutamate in repeat II. Monte Carlo energy minimizations yielded many complexes with calcium ions bound to at least two selectivity filter carboxylates. In these complexes calcium-carboxylate attractions are counterbalanced by calcium-calcium and carboxylate-carboxylate repulsions. Superposition of the complexes suggests a high degree of mobility of calcium ions and carboxylate groups of the glutamates. We used the predicted complexes to propose a permeation mechanism that involves single-file movement of calcium ions. The key feature of this mechanism is the presence of bridging glutamates that coordinate two calcium ions and enable their transitions between different chelating patterns involving four to six oxygen atoms from the channel protein. The conserved aspartate is proposed to coordinate a calcium ion incoming to the selectivity filter from the extracellular side. Glutamates in repeats III and IV, which are most distant from the repeat II aspartate, are proposed to coordinate the calcium ion that leaves the selectivity filter to the inner pore. Published experimental data and earlier proposed permeation models are discussed in view of our model.

Docking of mu-conotoxin GIIIA in the sodium channel outer vestibule

mu-Conotoxin GIIIA (mu-CTX) is a high-affinity ligand for the outer vestibule of selected isoforms of the voltage-gated Na(+) channel. The detailed bases for the toxin's high affinity binding and isoform selectivity are unclear. The outer vestibule is lined by four pore-forming (P) loops, each with an acidic residue near the mouth of the vestibule. mu-CTX has seven positively charged residues that may interact with these acidic P-loop residues. Using pair-wise alanine replacement of charged toxin and channel residues, in conjunction with double mutant cycle analysis, we determined coupling energies for specific interactions between each P-loop acidic residue and selected toxin residues to systematically establish quantitative restraints on the toxin orientation in the outer vestibule. Xenopus oocytes were injected with the mutant or native Na(+) channel mRNA, and currents measured by two-electrode voltage clamp. Mutant cycle analysis revealed novel, strong, toxin-channel interactions between K9/E403, K11/D1241, K11/D1532, and R19/D1532. Experimentally determined coupling energies for interacting residue pairs provided restraints for molecular dynamics simulations of mu-CTX docking. Our simulations suggest a refined orientation of the toxin in the pore, with toxin basic side-chains playing key roles in high-affinity binding. This modeling also provides a set of testable predictions for toxin-channel interactions, hitherto not described, that may contribute to high-affinity binding and channel isoform selectivity.

Structure/Function Characterization of μ-Conotoxin KIIIA, an Analgesic, Nearly Irreversible Blocker of Mammalian Neuronal Sodium Channels

Peptide neurotoxins from cone snails continue to supply compounds with therapeutic potential. Although several analgesic conotoxins have already reached human clinical trials, a continuing need exists for the discovery and development of novel non-opioid analgesics, such as subtype-selective sodium channel blockers. Micro-conotoxin KIIIA is representative of micro-conopeptides previously characterized as inhibitors of tetrodotoxin (TTX)-resistant sodium channels in amphibian dorsal root ganglion neurons. Here, we show that KIIIA has potent analgesic activity in the mouse pain model. Surprisingly, KIIIA was found to block most (>80%) of the TTX-sensitive, but only approximately 20% of the TTX-resistant, sodium current in mouse dorsal root ganglion neurons. KIIIA was tested on cloned mammalian channels expressed in Xenopus oocytes. Both Na(V)1.2 and Na(V)1.6 were strongly blocked; within experimental wash times of 40-60 min, block was reversed very little for Na(V)1.2 and only partially for Na(V)1.6. Other isoforms were blocked reversibly: Na(V)1.3 (IC50 8 microM), Na(V)1.5 (IC50 284 microM), and Na(V)1.4 (IC50 80 nM). "Alanine-walk" and related analogs were synthesized and tested against both Na(V)1.2 and Na(V)1.4; replacement of Trp-8 resulted in reversible block of Na(V)1.2, whereas replacement of Lys-7, Trp-8, or Asp-11 yielded a more profound effect on the block of Na(V)1.4 than of Na(V)1.2. Taken together, these data suggest that KIIIA is an effective tool to study structure and function of Na(V)1.2 and that further engineering of micro-conopeptides belonging to the KIIIA group may provide subtype-selective pharmacological compounds for mammalian neuronal sodium channels and potential therapeutics for the treatment of pain.

Speeding the recovery from ultraslow inactivation of voltage-gated Na+ channels by metal ion binding to the selectivity filter: a foot-on-the-door?

Slow inactivated states in voltage-gated ion channels can be modulated by binding molecules both to the outside and to the inside of the pore. Thus, external K(+) inhibits C-type inactivation in Shaker K(+) channels by a "foot-in-the-door" mechanism. Here, we explore the modulation of a very long-lived inactivated state, ultraslow inactivation (I(US)), by ligand binding to the outer vestibule in voltage-gated Na(+) channels. Blocking the outer vestibule by a mutant mu-conotoxin GIIIA substantially accelerated recovery from I(US). A similar effect was observed if Cd(2+) was bound to a cysteine engineered to the selectivity filter (K1237C). In K1237C channels, exposed to 30 microM Cd(2+), the time constant of recovery from I(US) was decreased from 145.0 +/- 10.2 s to 32.5 +/- 3.3 s (P < 0.001). Recovery from I(US) was only accelerated if Cd(2+) was added to the bath solution during recovery (V = -120 mV) from I(US), but not when the channels were selectively exposed to Cd(2+) during the development of I(US) (-20 mV). These data could be explained by a kinetic model in which Cd(2+) binds with high affinity to a slow inactivated state (I(S)), which is transiently occupied during recovery from I(US). A total of 50 microM Cd(2+) produced an approximately 8 mV hyperpolarizing shift of the steady-state inactivation curve of I(S), supporting this kinetic model. Binding of lidocaine to the internal vestibule significantly reduced the number of channels entering I(US), suggesting that I(US) is associated with a conformational change of the internal vestibule of the channel. We propose a molecular model in which slow inactivation (I(S)) occurs by a closure of the outer vestibule, whereas I(US) arises from a constriction of the internal vestibule produced by a widening of the selectivity filter region. Binding of Cd(2+) to C1237 promotes the closure of the selectivity filter region, thereby hastening recovery from I(US). Thus, Cd(2+) ions may act like a foot-on-the-door, kicking the I(S) gate to close.

Insect-selective spider toxins targeting voltage-gated sodium channels

The voltage-gated sodium (Na(v)) channel is a target for a number of drugs, insecticides and neurotoxins. These bind to at least seven identified neurotoxin binding sites and either block conductance or modulate Na(v) channel gating. A number of peptide neurotoxins from the venoms of araneomorph and mygalomorph spiders have been isolated and characterized and determined to interact with several of these sites. These all conform to an 'inhibitor cystine-knot' motif with structural, but not sequence homology, to a variety of other spider and marine snail toxins. Of these, spider toxins several show phyla-specificity and are being considered as lead compounds for the development of biopesticides. Hainantoxin-I appears to target site-1 to block Na(v) channel conductance. Magi 2 and Tx4(6-1) slow Na(v) channel inactivation via an interaction with site-3. The delta-palutoxins, and most likely mu-agatoxins and curtatoxins, target site-4. However, their action is complex with the mu-agatoxins causing a hyperpolarizing shift in the voltage-dependence of activation, an action analogous to scorpion beta-toxins, but with both delta-palutoxins and mu-agatoxins slowing Na(v) channel inactivation, a site-3-like action. In addition, several other spider neurotoxins, such as delta-atracotoxins, are known to target both insect and vertebrate Na(v) channels most likely as a result of the conserved structures within domains of voltage-gated ion channels across phyla. These toxins may provide tools to establish the molecular determinants of invertebrate selectivity. These studies are being greatly assisted by the determination of the pharmacophore of these toxins, but without precise identification of their binding site and mode of action their potential in the above areas remains underdeveloped.