Selectivity filter residues contribute unequally to pore stabilization in voltage-gated sodium channels

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Structural basis for antiarrhythmic drug interactions with the human cardiac sodium channel

Voltage-gated sodium channels play a central role in cellular excitability and are key targets for drug development. Recent breakthroughs in high-resolution cryo-electron microscopy protein structure determination, Rosetta computational protein structure modeling, and multimicrosecond molecular dynamics simulations are empowering advances in structural biology to study the atomistic details of channel−drug interactions. We used Rosetta structural computational modeling and molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with cardiac sodium channel. Our results provide crucial atomic-scale mechanistic insights into the channel–drug interactions, necessary for the rational design of novel modulators of the human cardiac sodium channel to be used for the treatment of cardiac arrhythmias.

The human voltage-gated sodium channel, hNa_V1.5, is responsible for the rapid upstroke of the cardiac action potential and is target for antiarrhythmic therapy. Despite the clinical relevance of hNa_V1.5-targeting drugs, structure-based molecular mechanisms of promising or problematic drugs have not been investigated at atomic scale to inform drug design. Here, we used Rosetta structural modeling and docking as well as molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with hNa_V1.5. These calculations revealed several key drug binding sites formed within the pore lumen that can simultaneously accommodate up to two drug molecules. Molecular dynamics simulations identified a hydrophilic access pathway through the intracellular gate and a hydrophobic access pathway through a fenestration between DIII and DIV. Our results advance the understanding of molecular mechanisms of antiarrhythmic and local anesthetic drug interactions with hNa_V1.5 and will be useful for rational design of novel therapeutics.

Large-effect mutations generate trade-off between predatory and locomotor ability during arms race coevolution with deadly prey

Adaptive evolution in response to one selective challenge may disrupt other important aspects of performance. Such evolutionary trade‐offs are predicted to arise in the process of local adaptation, but it is unclear if these phenotypic compromises result from the antagonistic effects of simple amino acid substitutions. We tested for trade‐offs associated with beneficial mutations that confer tetrodotoxin (TTX) resistance in the voltage‐gated sodium channel (Na_V1.4) in skeletal muscle of the common garter snake (Thamnophis sirtalis). Separate lineages in California and the Pacific Northwest independently evolved TTX‐resistant changes to the pore of Na_V1.4 as a result of arms race coevolution with toxic prey, newts of the genus Taricha. Snakes from the California lineage that were homozygous for an allele known to confer large increases in toxin resistance (Na_V1.4^LVNV) had significantly reduced crawl speed compared to individuals with the ancestral TTX‐sensitive channel. Heterologous expression of native snake Na_V1.4 proteins demonstrated that the same Na_V1.4^LVNV allele confers a dramatic increase in TTX resistance and a correlated decrease in overall channel excitability. Our results suggest the same mutations that accumulate during arms race coevolution and beneficially interfere with toxin‐binding also cause changes in electrophysiological function of the channel that may affect organismal performance. This trade‐off was only evident in the predator lineage where coevolution has led to the most extreme resistance phenotype, determined by four critical amino acid substitutions. If these biophysical changes also translate to a fitness cost—for example, through the inability of T. sirtalis to quickly escape predators—then pleiotropy at this single locus could contribute to observed variation in levels of TTX resistance across the mosaic landscape of coevolution.

Distinct modulation of inactivation by a residue in the pore domain of voltage-gated Na+ channels: mechanistic insights from recent crystal structures

Inactivation of voltage-gated Na⁺ channels (VGSC) is essential for the regulation of cellular excitability. The molecular rearrangement underlying inactivation is thought to involve the intracellular linker between domains III and IV serving as inactivation lid, the receptor for the lid (domain III S4-S5 linker) and the pore-lining S6 segements. To better understand the role of the domain IV S6 segment in inactivation we performed a cysteine scanning mutagenesis of this region in rNav 1.4 channels and screened the constructs for perturbations in the voltage-dependence of steady state inactivation. This screen was performed in the background of wild-type channels and in channels carrying the mutation K1237E, which profoundly alters both permeation and gating-properties. Of all tested constructs the mutation I1581C was unique in that the mutation-induced gating changes were strongly influenced by the mutational background. This suggests that I1581 is involved in specific short-range interactions during inactivation. In recently published crystal structures VGSCs the respective amino acids homologous to I1581 appear to control a bend of the S6 segment which is critical to the gating process. Furthermore, I1581 may be involved in the transmission of the movement of the DIII voltage-sensor to the domain IV S6 segment.

Modeling the human Nav1.5 sodium channel: structural and mechanistic insights of ion permeation and drug blockade

Abnormalities in the human Nav1.5 (hNav1.5) voltage-gated sodium ion channel (VGSC) are associated with a wide range of cardiac problems and diseases in humans. Current structural models of hNav1.5 are still far from complete and, consequently, their ability to study atomistic interactions of this channel is very limited. Here, we report a comprehensive atomistic model of the hNav1.5 ion channel, constructed using homology modeling technique and refined through long molecular dynamics simulations (680 ns) in the lipid membrane bilayer. Our model was comprehensively validated by using reported mutagenesis data, comparisons with previous models, and binding to a panel of known hNav1.5 blockers. The relatively long classical MD simulation was sufficient to observe a natural sodium permeation event across the channel's selectivity filters to reach the channel's central cavity, together with the identification of a unique role of the lysine residue. Electrostatic potential calculations revealed the existence of two potential binding sites for the sodium ion at the outer selectivity filters. To obtain further mechanistic insight into the permeation event from the central cavity to the intracellular region of the channel, we further employed "state-of-the-art" steered molecular dynamics (SMD) simulations. Our SMD simulations revealed two different pathways through which a sodium ion can be expelled from the channel. Further, the SMD simulations identified the key residues that are likely to control these processes. Finally, we discuss the potential binding modes of a panel of known hNav1.5 blockers to our structural model of hNav1.5. We believe that the data presented here will enhance our understanding of the structure-property relationships of the hNav1.5 ion channel and the underlying molecular mechanisms in sodium ion permeation and drug interactions. The results presented here could be useful for designing safer drugs that do not block the hNav1.5 channel.

Structure and function of voltage-gated sodium channels at atomic resolution

Voltage-gated sodium channels initiate action potentials in nerve, muscle and other excitable cells. Early physiological studies described sodium selectivity, voltage-dependent activation and fast inactivation, and developed conceptual models for sodium channel function. This review article follows the topics of my 2013 Sharpey-Schafer Prize Lecture and gives an overview of research using a combination of biochemical, molecular biological, physiological and structural biological approaches that have elucidated the structure and function of sodium channels at the atomic level. Structural models for voltage-dependent activation, sodium selectivity and conductance, drug block and both fast and slow inactivation are discussed. A perspective for the future envisions new advances in understanding the structural basis for sodium channel function and the opportunity for structure-based discovery of novel therapeutics.

Use-dependent block of the voltage-gated Na(+) channel by tetrodotoxin and saxitoxin: effect of pore mutations that change ionic selectivity

Voltage-gated Na(+) channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca(2+) permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca(2+) or Na(+) ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca(2+) permeability, suggesting that ion-toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.

A molecular switch between the outer and the inner vestibules of the voltage-gated Na+ channel

Voltage-gated ion channels are transmembrane proteins that undergo complex conformational changes during their gating transitions. Both functional and structural data from K⁺ channels suggest that extracellular and intracellular parts of the pore communicate with each other via a trajectory of interacting amino acids. No crystal structures are available for voltage-gated Na⁺ channels, but functional data suggest a similar intramolecular communication involving the inner and outer vestibules. However, the mechanism of such communication is unknown. Here, we report that amino acid Ile-1575 in the middle of transmembrane segment 6 of domain IV (DIV-S6) in the adult rat skeletal muscle isoform of the voltage-gated sodium channel (rNa_V1.4) may act as molecular switch allowing for interaction between outer and inner vestibules. Cysteine scanning mutagenesis of the internal part of DIV-S6 revealed that only mutations at site 1575 rescued the channel from a unique kinetic state (“ultra-slow inactivation,” I_US) produced by the mutation K1237E in the selectivity filter. A similar effect was seen with I1575A. Previously, we reported that conformational changes of both the internal and the external vestibule are involved in the generation of I_US. The fact that mutations at site 1575 modulate I_US produced by K1237E strongly suggests an interaction between these sites. Our data confirm a previously published molecular model in which Ile-1575 of DIV-S6 is in close proximity to Lys-1237 of the selectivity filter. Furthermore, these functional data define the position of the selectivity filter relative to the adjacent DIV-S6 segment within the ionic permeation pathway.

Biophysical costs associated with tetrodotoxin resistance in the sodium channel pore of the garter snake, Thamnophis sirtalis

Tetrodotoxin (TTX) is a potent toxin that specifically binds to voltage-gated sodium channels (NaV). TTX binding physically blocks the flow of sodium ions through NaV, thereby preventing action potential generation and propagation. TTX has different binding affinities for different NaV isoforms. These differences are imparted by amino acid substitutions in positions within, or proximal to, the TTX-binding site in the channel pore. These substitutions confer TTX-resistance to a variety of species. The garter snake Thamnophis sirtalis has evolved TTX-resistance over the course of an arms race, allowing some populations of snakes to feed on tetrodotoxic newts, including Taricha granulosa. Different populations of the garter snake have different degrees of TTX-resistance, which is closely related to the number of amino acid substitutions. We tested the biophysical properties and ion selectivity of NaV of three garter snake populations from Bear Lake, Idaho; Warrenton, Oregon; and Willow Creek, California. We observed changes in gating properties of TTX-resistant (TTXr) NaV. In addition, ion selectivity of TTXr NaV was significantly different from that of TTX-sensitive NaV. These results suggest TTX-resistance comes at a cost to performance caused by changes in the biophysical properties and ion selectivity of TTXr NaV.

Mechanisms of action of ligands of potential-dependent sodium channels

Potential-dependent sodium channels play a leading role in generating action potentials in excitable cells. Sodium channels are the site of action of a variety of modulator ligands. Despite numerous studies, the mechanisms of action of many modulators remain incompletely understood. The main reason that many important questions cannot be resolved is that there is a lack of precise data on the structures of the channels themselves. Structurally, potential-dependent sodium channels are members of the P-loop channel superfamily, which also include potassium and calcium channels and glutamate receptor channels. Crystallization of a series of potassium channels showed that it was possible to analyze the structures of different members of the superfamily using the "homologous modeling" method. The present study addresses model investigations of the actions of ligands of sodium channels, including tetrodotoxin and batrachotoxin, as well as local anesthetics. Comparison of experimental data on sodium channel ligands with x-ray analysis data allowed us to reach a new level of understanding of the mechanisms of channel modulation and to propose a series of experimentally verifiable hypotheses.

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.

Sodium channels: ionic model of slow inactivation and state-dependent drug binding

Inactivation is a fundamental property of voltage-gated ion channels. Fast inactivation of Na(+) channels involves channel block by the III-IV cytoplasmic interdomain linker. The mechanisms of nonfast types of inactivation (intermediate, slow, and ultraslow) are unclear, although the ionic environment and P-loops rearrangement appear to be involved. In this study, we employed a TTX-based P-loop domain model of a sodium channel and the MCM method to investigate a possible role of P-loop rearrangement in the nonfast inactivation. Our modeling predicts that Na(+) ions can bind between neighboring domains in the outer-carboxylates ring EEDD, forming an ordered structure with interdomain contacts that stabilize the conducting conformation of the outer pore. In this model, the permeant ions can transit between the EEDD ring and the selectivity filter ring DEKA, retaining contacts with at least two carboxylates. In the absence of Na(+), the electrostatic repulsion between the EEDD carboxylates disrupts the permeable configuration. In this Na(+)-deficient model, the region between the EEDD and DEKA rings is inaccessible for Na(+) but is accessible for TMA. Taken together, these results suggest that Na(+)-saturated models are consistent with experimental characteristics of the open channels, whereas Na(+)-deficient models are consistent with experimentally defined properties of the slow-inactivated channels. Our calculations further predict that binding of LAs to the inner pore would depend on whether Na(+) occupies the DEKA ring. In the absence of Na(+) in the DEKA ring, the cationic group of lidocaine occurs in the focus of the pore helices' macrodipoles and would prevent occupation of the ring by Na(+). Loading the DEKA ring with Na(+) results in the electrostatic repulsion with lidocaine. Thus, there are antagonistic relations between a cationic ligand bound in the inner pore and Na(+) in the DEKA ring.