The cardiac sodium channel alpha subunit (RHI) is less sensitive to tetrodotoxin (TTX) and saxitoxin (STX) and more sensitive to cadmium than brain and skeletal muscle (microliter) isoforms. An RHI mutant, with Tyr substituted for Cys at position 374 (as in microliter) confers three properties of TTX-sensitive channels: (i) greater sensitivity to TTX (730-fold); (ii) lower sensitivity to cadmium (28-fold); and (iii) altered additional block by toxin upon repetitive stimulation. Thus, the primary determinant of high-affinity TTX-STX binding is a critical aromatic residue at position 374, and the interaction may take place possibly through an ionized hydrogen bond. This finding requires revision of the sodium channel pore structure that has been previously suggested by homology with the potassium channel.

The crystal structure of the pore-forming part of the KcsA bacterial K(+)-selective channel suggests a possible motif for related voltage-gated channels. We examined the hypothesis that the spacial orientation of the KcsA M1 and M2 alpha-helices also predicts the backbone location of S5 and S6 helices of the voltage-gated Na(+) channel. That channel's P region structure is expected to be different because selectivity is determined by side-chain interactions rather than by main-chain carbonyls, and its outer vestibule accommodates relatively large toxin molecules, tetrodotoxin (TTX) and saxitoxin (STX), which interact with selectivity ring residues. The Na(+) channel P loop was well-modeled by the alpha-helix-turn-beta-strand motif, which preserves the relationships for toxin interaction with the Na(+) channel found experimentally. This outer vestibule was docked into the extracellular part of the inverted teepee structure formed by the S5 and S6 helices that were spacially located by coordinates of the KcsA M1 and M2 helix main chains [Doyle et al. (1998) Science 280, 69-74], but populated with side chains of the respective S5 and S6 structures. van der Waals contacts were optimized with minimal adjustment of the S5, S6, and P loop structures, forming a densely packed pore structure. Nonregular external S5-P and P-S6 segments were not modeled here, except the P-S6 segment of domain II. The resulting selectivity region structure is consistent with Na(+) channel permeation properties, offering suggestions for the molecular processes involved in selectivity. The ability to construct a Na(+) channel pore model consistent with most of the available biophysical and mutational information suggests that the KcsA structural framework may be conserved in voltage-gated channels.

While studying the adult rat skeletal muscle Na+ channel outer vestibule, we found that certain mutations of the lysine residue in the domain III P region at amino acid position 1237 of the alpha subunit, which is essential for the Na+ selectivity of the channel, produced substantial changes in the inactivation process. When skeletal muscle alpha subunits (micro1) with K1237 mutated to either serine (K1237S) or glutamic acid (K1237E) were expressed in Xenopus oocytes and depolarized for several minutes, the channels entered a state of inactivation from which recovery was very slow, i.e., the time constants of entry into and exit from this state were in the order of approximately 100 s. We refer to this process as "ultra-slow inactivation". By contrast, wild-type channels and channels with the charge-preserving mutation K1237R largely recovered within approximately 60 s, with only 20-30% of the current showing ultra-slow recovery. Coexpression of the rat brain beta1 subunit along with the K1237E alpha subunit tended to accelerate the faster components of recovery from inactivation, as has been reported previously of native channels, but had no effect on the mutation-induced ultra-slow inactivation. This implied that ultra-slow inactivation was a distinct process different from normal inactivation. Binding to the pore of a partially blocking peptide reduced the number of channels entering the ultra-slow inactivation state, possibly by interference with a structural rearrangement of the outer vestibule. Thus, ultra-slow inactivation, favored by charge-altering mutations at site 1237 in micro1 Na+ channels, may be analogous to C-type inactivation in Shaker K+ channels.

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.

Biophysical evidence has placed the binding site for the naturally occurring marine toxins tetrodotoxin (TTX) and saxitoxin (STX) in the external mouth of the Na+ channel ion permeation pathway. We developed a molecular model of the binding pocket for TTX and STX, composed of antiparallel beta-hairpins formed from peptide segments of the four S5-S6 loops of the voltage-gated Na+ channel. For TTX the guanidinium moiety formed salt bridges with three carboxyls, while two toxin hydroxyls (C9-OH and C10-OH) interacted with a fourth carboxyl on repeats I and II. This alignment also resulted in a hydrophobic interaction with an aromatic ring of phenylalanine or tyrosine residues for the brainII and skeletal Na+ channel isoforms, but not with the cysteine found in the cardiac isoform. In comparison to TTX, there was an additional interaction site for STX through its second guanidinium group with a carboxyl on repeat IV. This model satisfactorily reproduced the effects of mutations in the S5-S6 regions and the differences in affinity by various toxin analogs. However, this model differed in important ways from previously published models for the outer vestibule and the selectivity region of the Na+ channel pore. Removal of the toxins from the pocket formed by the four beta-hairpins revealed a structure resembling a funnel that terminated in a narrowed region suitable as a candidate for the selectivity filter of the channel. This region contained two carboxyls (Asp384 and Glu942) that substituted for molecules of water from the hydrated Na+ ion. Simulation of mutations in this region that have produced Ca2+ permeation of the Na+ channel created a site with three carboxyls (Asp384, Glu942, and Glu1714) in proximity.