|Home | About | Journals | Submit | Contact Us | Français|
In a voltage-dependent sodium channel, the activation of voltage-sensors upon depolarization leads to the opening of pore gates. To elucidate the principles underlying this conformational coupling, we investigated a putative gating interface in domain III of the sodium channel using voltage-clamp fluorimetry and tryptophan-scanning mutagenesis. Most mutations have similar energetic effects on voltage-sensor activation and pore opening. However, several mutations stabilized the activated voltage-sensor while concurrently destabilizing the open pore. Mapping of these mutations onto a homology model of the sodium channel showed that most localize to hinge regions of the gating interface. Our analysis shows that these residues are involved in energetic coupling of the voltage-sensor and the pore when both are in the resting or activated conformations, supporting the notion that electromechanical coupling in a voltage-dependent ion channel involves the movement of rigid segments connected by elastic hinges.
Voltage-gated sodium channels respond to membrane depolarization by undergoing a series of conformational changes that culminate in channel opening and subsequent inactivation. Ionic fluxes through these channels underlie electrical signaling in many physiological processes particularly in the nervous system and the cardiovascular system. The voltage-dependent sodium channel shares a common architecture with other members of the voltage-gated ion channel family and comprises a central ion-conducting pore surrounded by four voltage-sensing domains1–3. The “outward” movement of charged residues in the voltage-sensors in response to a depolarization triggers a conformational change that is transmitted to the pore gates4–10. The structural mechanisms underlying conformational coupling of the voltage-sensor and pore gates is one of the central questions in ion channel biology11.
The principal subunit of the voltage-gated sodium channel is a large 260 kDa polypeptide and has four homologous domains (DI to DIV) which together form a functional channel12 (Figure 1a). Each domain of the sodium channel is equivalent to a single subunit of the voltage-dependent potassium channel. The S4 transmembrane segments in the sodium channel are highly conserved across the voltage-gated ion channel family and along with the S1–S3 segments, constitute the voltage-sensing domain of the channel 13 (Figure 1b). The S5 and S6 helices of the four domains together form the central pore of the channel. The high-resolution structures of the eukaryotic potassium channels, presumably representing the open conformation of the channel, reveal that the voltage-sensing domains make limited physical contacts with the pore 14,15. Studies in the Shaker potassium channel indicate that at least two putative pathways mediate the transmission of “information” from the voltage-sensor to the central pore 16–20. Scanning mutagenesis and domain swap experiments combined with detailed electrophysiological studies have highlighted the importance of the S4, the S4–S5 linker and the intracellular ends of the S5 and S6 segments 17,20,21. This pathway likely couples the voltage-sensor and the pore within the same subunit or domain. Tryptophan scanning along with statistical coupling analysis has uncovered another interface constituted in part by the S1 segment, which may be involved in coupling the voltage-sensing domain of one subunit to the pore of the neighboring subunit 18,19. These studies are beginning to unravel the molecular mechanisms underlying the propagation of structural change from the voltage-sensor to the pore of a voltage-dependent ion channel.
In comparison to the voltage-gated potassium channels, little is known about the molecular determinants of electromechanical coupling in the sodium channels. The complexity of the sodium channel structure and an absence of functional symmetry make interpretation of structure-function data extremely challenging. Recently, Swartz and colleagues 22 have described an elegant methodology to dissect the functional behavior of parts of the sodium channel by transplanting them into a symmetric potassium channel template. Here, we have used a novel approach that combines site-specific fluorescence recordings with conductance measurements to examine the molecular basis of coupling between the voltage-sensor and the pore of the sodium channel. We have systematically mutated residues in the S4–S5 linker, parts of the S5 and S6 of domain III (Figure 1) to study the effect of these mutations on coupling interactions between voltage-sensor and pore of the sodium channel. Site-specific fluorescence recordings determined the effect of these mutations on the thermodynamic parameters of the domain III voltage-sensor whereas the effect on the pore conformation was ascertained by conductance measurements. Upon mapping a small set of critical residues on to a structural model of the sodium channel, we find that a majority of them are either in close proximity to or in the hinge regions of key helices such as the S4, S4–S5 linker, S5 and S6 segments. Two of these, arginine (R1135) at the beginning of S4–S5 linker and a phenylalanine (F1298) at the tail end of the S6 segment, are likely to form an interacting pair. Our analysis shows that these residues are involved in coupling interactions both in the closed and open conformation of the channel.
In this study, the role of the residues in a putative interaction surface likely to be involved in voltage sensor-pore coupling was assessed by tryptophan scanning mutagenesis approach. A mutation to tryptophan, which has the largest side chain volume amongst naturally occurring amino acids, has the potential to disrupt protein-protein interactions radically by introducing defects in the interfaces stabilizing the protein core 18,23–25. We introduced tryptophans in the rat skeletal muscle sodium channel (Nav 1.4), which is one of the better-characterized members of the sodium channel family. Fluorophores attached to cysteine residues proximal to the S4 segments of this channel were previously shown to track the movement of gating charges 26. Here, the conformational state of the voltage-sensor of domain III was monitored by a fluorophore attached to substituted cysteine at the L1115 position (Figure 1a) using a voltage-clamp fluorimetry setup. Thus, all the mutants were made in L1115C background, which is hereafter referred to as the wild type (WT). The conformation of the pore gates was monitored by conductance measurements in a two-electrode voltage-clamp setup. Mutations of residues in the S6-domain III are known to significantly alter the stability of the sodium channel pore 27. Furthermore, the mutations of charged residues in the voltage-sensor of domain III reduce the slope of conductance-voltage relationships suggesting that the movement of domain III gating charges is coupled to opening of the pore gates 28. Taken together, these studies suggest that any direct or indirect alteration in the conformational state of the pore helices of domain III is likely to be manifested in the conductance properties of the sodium channel.
The tryptophan substitutions were introduced in the S4–S5 linker and parts of S5 and S6 towards the intracellular end (Figure 1b). These substitutions encompassed the gating interface identified by Lu and colleagues 17 in the Shaker potassium channel. Amongst the 54 tryptophan mutants, functional channel activity was observed in all but 10 mutants (Table I). Tryptophans at these 10 positions presumably disrupt protein packing which causes either a complete loss of channel function or low plasma membrane expression. In all such instances, substitution to alanine instead of tryptophan generated a functional channel and allowed us to probe the role of side chains at those sites. Figure 2a shows typical voltage-dependent ionic current and fluorescence traces of the WT and few of the mutants, representing a broad spectrum of functionally altered channels. The conductance-voltage (G-V) relationships were obtained from peak conductance values whereas the fluorescence-voltage (F-V) relationships represent the steady state fluorescence intensities (Figure 2b).
Structure-function studies on the brain sodium channel show that the steady state activation curves of the inactivation deficient mutants are not significantly different compared to the wild type channels29. This is mainly because the activation rates of the wild type channels are about 4–5 fold faster than the inactivation rates26. The ratios of activation to inactivation time constants of the mutants were similar to those of the WT sodium channels (Supplementary Figure 4). This suggests that inactivation is unlikely to have a significant effect on our estimates of activation parameters derived from peak conductance values.
To quantitatively evaluate the effects of the perturbation, we fitted each of the G-V and the F-V curves to a single Boltzmann function and obtained the values of two parameters, the slope (z) and the voltage eliciting half-maximal response (V1/2) (Table I). The difference in V1/2 of the mutant relative to the WT was expressed as ΔV1/2 (ΔV1/2 = V1/2 Mutant − V1/2 WT) and these differences calculated from the F-V and G-V curves were summarized in Figure 2c. The plot reveals that more than half of the tryptophan mutants cause a statistically significant shift in the voltage range describing the fluorescence response compared to the WT. These data highlight the importance of the residues in this linker region in voltage-sensor activation and/or pore opening.
For the most part, the tryptophan substitutions cause much a larger effect on the F-V rather than the G-V parameters. This trend can be partly accounted by the fact that the fluorescence specifically monitors the movement of domain III voltage-sensor whereas conductance requires contribution from the pore segments of all four domains. Based on their functional effects, the mutants can be broadly categorized into two groups: a) those that shift the ΔV1/2 of both the G-V and F-V curves in the same direction (referred hereafter as class I mutations) and, b) those that shift ΔV1/2 of F-V and G-V curves in the opposite direction (class II mutations). Class I mutations have similar effects on the energetics of voltage-sensor activation and pore opening whereas the Class II mutations have opposite effects on those processes (Supplementary figure 1). For instance, tryptophan substitution at the L1133 position generates a class II mutant, causing a 12 mV rightward shift in the V1/2 of the G-V curve and a 25 mV leftward shift in the V1/2 of the F-V curve (Figure 4). This mutation stabilizes the voltage-sensor in the activated state and concurrently stabilizes the closed conformation of the pore. In the following sections, we will analyze these mutations and their implication in detail.
An implicit assumption in this sort of analysis is that the mutations do not induce a novel conformational change but only affect the thermodynamic parameters that define the coupled system. The time course (at saturating potentials) and the direction of fluorescence change (after accounting for shifts in V1/2) in the mutants were similar to the WT (Supplementary Figure 5). Furthermore, the amplitudes of fluorescence signals of the mutants are broadly correlated with channel expression as estimated from the peak outward currents obtained during fluorescence recordings. To quantitate the effect of the mutations on fluorescence amplitudes, it is necessary to estimate the number of channels per cell and the fraction of labeled channels. Although, gating current measurements have been utilized to normalize for channel density, we cannot rule out the possibility that these mutations may also modify the total gating charge per channel by altering the structure. Recent studies on the Shaker potassium channel show that the total gating currents is highly sensitive to the position of the charge in a focused electric field30 and it is conceivable that a mutation could alter the total charge per channel by reshaping the electric field31 or by hindering charge movement. Given these considerations, we cannot completely rule out the possibility that either the fluorescence amplitudes or gating charge, in some instances, was modified by these mutations. Nevertheless, our overall fluorescence data is consistent with the notion that the probe is likely to track the same structural rearrangements in the WT and mutant channels.
The voltage range of activation (V1/2) of both F-V and G-V curves were shifted to more negative voltages in fourteen mutants suggesting that the energy required to activate the voltage-sensor and pore opening is reduced (Figure 2c and Table I). In order to evaluate this data in a structural context, we generated a homology model of the sodium channel structure (Figure 3a, red; and Supplementary figure 2a and c) based on the crystal structure of Kv1.2/2.1 chimera (Figure 3a, blue; and Supplementary figure 2a) (PDB code: 2R9R; see Materials and Methods) 15. The model of Nav1.4 has two major differences: first, the sodium channel has a shorter S3–S4 linker compared to the chimeric potassium channel, and second, the break between the S4 segment and the S4–S5 linker occurs immediately following the fourth charge. In the eukaryotic potassium channel, the break is after one helical turn. Mapping those positions onto the structural model of the sodium channel (Figure 3b) shows that many of these residues (8 of the 14 residues) are clustered in the intracellular region of the S5 and S6 segments (Figure 3c). Residues in S5 (I1150, P1151, M1154, V1156 and L1157) face the clustered residues in S6 (L1282, F1283 and I1284) and away from the voltage-sensor of domain III. The side chains of these residues project towards the S4 and S4–S5 linker of the neighboring domain (Domain II) and, thus, the mutations at these positions may perturb inter-domain interactions leading to altered activation profiles. Residues in this region have been shown to be important in activation gating of the Shaker potassium channel 20. The remaining six residues are positioned throughout the S4–S5 linker and the inner portion of S6 (Figure 3d). In contrast to the fourteen residues that cause leftward shifts in the F-V and G-V curves, mutations at two positions cause rightward shifts in both curves. N1144 is at the C-terminal end of the S4–S5 linker while K1297 is at the end of the S6 helix (Figure 3d). These two mutations in all likelihood stabilize the closed conformation of the channel. While there are a few such exceptions, we assume that for most part large aromatic substitutions such as the introduction of tryptophans disrupt protein packing. Thus, our finding that a large fraction of the mutants destabilize the closed state of the channel may reflect increased steric clashes in the closed state of the channel.
Of the 54 sites probed in this study, tryptophan mutations in seven sites corresponding to L1133, S1134, R1135, A1149, F1278, N1281 and V1286 showed left shift in the F-V and a right shift in the G-V relationships. Thus, the mutations at these positions cause voltage-sensors to activate more easily while concurrently making the opening of the pore more difficult. Figure 4a represents the F-Vs and G-Vs of four of the mutants along with those of the WT. Right shifts of G-V curves were in the range of 2.6 to 12.7 mV and left shifts of F-V relationships were in the range of 10.9 to 66.9 mV (Table I). A stretch of three residues (L1133W, S1134 and R1135) occurs at the hinge linking the S4 segment with the S4–S5 linker (Figure 5a and c). Remarkably, three other residues are also either in the hinge or in close proximity. A1149 is at the hinge between the S4–S5 linker and S5 segment (Figure 5b) and V1286 is in close proximity to this bend but in the S6 helix (Figure 5c). F1278 is near the putative gating hinge on the S6 helix (Figure 5b). Localization of these mutations to the hinges in the structural model implies that these regions may play a unique role in the gating process.
Mutations at three positions (Q1294, K1296 and F1298) in the lower half of the S6 helix caused a rightward shift in the F-V and a left shift in the G-V relationships bringing their voltage-range closer to each other (Figure 4b). Of these three, the F-V curve of Q1294W showed the largest depolarizing shift. A simplistic interpretation is that the tryptophan mutations at these positions increase the strength of coupling between the voltage-sensor and the pore domain. As a result, the functional properties of these two domains resemble each other. Conversely, we can infer that the class II mutations, which stabilize the activated voltage-sensor and destabilize the open pore, decouple the voltage-sensor and pore movements. We will revisit these interpretations in the discussion section.
The crystal structure of both Kv1.2 and the Kv1.2/2.1 chimera shows that the S4–S5 linker is juxtaposed to the S6 helix of the same subunit 14,15. The S4–S5 linker is almost parallel to the membrane plane and positioned over the S6 inner helix forming an interaction surface with one face of the S6 helix. Studies in the HERG potassium channel 32 have shown that positively charged residues in the S4–S5 linker are likely to interact with the S6 segment of the potassium channel. Sequence alignment of different members of the voltage-gated channel family shows that the tail end of the S6 segment is rich in negatively charged residues or aromatic groups. In our structural model of the sodium channel, two phenylanines (F1298 and F1136) were observed to be in proximity to positive charges (R1135 and K1296) on S4–S5 linker in the tertiary structure (Figure 6c) raising the possibility that these residues form a cation-π type of interaction 33. Tryptophan substitutions at both these positions showed a trend consistent with increased coupling. Tryptophans have a higher electron cloud density compared to phenylalanines which would contribute to a stronger cation-π interaction 34. In the model, the side chains of R1135 and F1298 are pointing towards each other and their α-carbons are within 8.6 Å of each other, close enough to interact. To test this possibility, we first mutated these residues to alanine. Interestingly, the opposite shifts in F-V and G-V curves observed in F1298W were reversed in F1298A (Figure 6a and b), suggesting that F1298, but not K1296 (Figure 6b), may interact with R1135 in the S4–S5 linker. To further probe this putative interaction between R1135 and F1298, we switched the positions of aromatic and cationic side chains to create R1135W-F1298R double mutant. The F-V and G-V curves of this double mutant were similar to the R1135W single mutant and not to the WT suggesting that these two groups are unable to form an electrostatic interaction pair. Note that cation-π interaction is geometry dependent-cation interacts in en face orientation as opposed to edge-on orientation with the aromatic side chain 33. It is possible that switching the positions of the interacting side-chains does not create an optimal orientation for such an interaction. Moreover, R1135 is in close proximity to other class II (L1133 and S1134) mutations and is in the putative hinge connecting the S4 to the S4–S5 linker. Substitution of arginine by a bulky tryptophan in this crucial putative hinge may introduce steric perturbations that prevent movements necessary for voltage-sensor-pore coupling. Thus, to fully examine the interaction between R1135 and F1298, it will be necessary to use a less disruptive substitutions such as those involving the use of non-natural amino acids 35.
By combining tryptophan-scanning mutagenesis with voltage-clamp fluorimetry, we have identified a number of key residues that are likely to be involved in conformational coupling of the domain III voltage-sensor to the pore of the Nav 1.4 channel. Class I mutants modify the energetics of voltage-sensor activation and opening of the pore in the same direction. To understand how a single amino acid substitution can cause such effects on the voltage-sensor and pore, we first consider two simple canonical models of a cooperative system where the two domains interact either when both are in their resting state or in their activated state (Scheme IA and Scheme IB in Figure 7a). The effect of changing each of the three thermodynamic parameters on the activation probabilities of the voltage-sensor and pore was tested by using numerical simulations (details in Supplementary material XX). The plots in Figures 7b–e show that the probabilities of voltage-sensor activation and pore opening for Schemes IA and IB increase when either K1or K2 increases. When the interaction term (θ) is increased, the voltage-sensor and pore opening probabilities for Scheme IA also increase whereas those probabilities for Scheme IB were reduced. To evaluate the generality of these trends, we analyzed the derivatives of the probability expressions (Supplementary Table I and II). A positive value of the derivative with respect to a parameter indicates that the probability value will increase as the parameter value is increased. The derivatives of voltage-sensor activation and pore opening also show that both these probabilities move in the same direction when any one of the thermodynamic terms is modified.
The above analysis shows that if a mutation shifts both the probability curves in the same direction, we cannot attribute the effect of a mutation on any one of thermodynamic terms in Scheme IA or IB. Nonetheless, we can make some inferences regarding the role of these amino acids based on their position in the structural model of the sodium channel. A majority of the class I mutations relatively stabilize the open state of the channel and a large fraction of these are in the S5 and S6 segments which constitutes the pore domain. In KcsA structure36, which serves as a model of the closed pore, the four S6 helices are tightly packed at the bundle crossing. Studies in the Shaker and the HERG channel show that the mutations in these regions disrupt the tight packing and destabilize the closed pore5,27,37–40. These findings suggest that the introduction of tryptophan in the S5 and S6 region mostly destabilizes the packing in the closed state of the sodium channel.
In contrast to class I mutations, the class II mutations cause opposite effects on the activation of the voltage-sensor and pore opening in the sodium channel. These “Janus-faced” mutants are reminiscent of the ILT mutant in the Shaker potassium channel 16,41. The voltage-sensor in ILT mutant is stabilized in the activated state whereas the open pore is highly destabilized. Our derivative analysis (Supplementary Table I and II) unequivocally shows that the opposing effects on the voltage-sensor and pore cannot be rationalized by changing a single parameter with Schemes IA & IB. We consider an alternate model of a cooperative system (Scheme II, Figure 7a), where the voltage-sensor and the pore interact with each other when both of them are in the doubly resting and doubly activated state. θ is the coupling term between the pore and voltage-sensor. We define this type of interaction as “conserved interactions” as opposed to “unconserved interactions” of the type depicted in Scheme IA & IB. Conserved interactions result in no net gain or loss of interaction energy going from the doubly resting to doubly activated state and achieve conformational coupling by destabilizing the intermediate states.
The probability expressions of activation for the voltage-sensor and the pore in scheme II is;
The WT F-V and G-V data were fitted to these probability expressions to obtain the initial estimates of the three thermodynamic parameters for numerical simulations. The plots in Figures 7f and g show the effect of changing each of these parameters on the probabilities of voltage-sensor activation and pore opening. Note that the probability of voltage-sensor activation decreases with increasing values of θ whereas that of pore opening increases. Analysis of derivatives of the above probability expressions also shows that a change in interaction term but not the intrinsic equilibrium constants has opposite effects on the probability terms under conditions when K2 <1 and K1>1 (Supplementary Table II). Thus in Scheme II, a loss of coupling would lead to a right shift in the G-V curve and a left-shift in the F-V curve. The assumption, K2 < 1 implies that the pore domain is intrinsically more stable in the closed or resting conformation. This is consistent with scanning mutagenesis studies on the pore segments 39,40. Our second assumption, K1> 1 implies that the voltage-sensor is intrinsically more stable in its activated conformation (in absence of membrane potential). The crystal structure of the isolated voltage-sensor domain of the prokaryotic voltage-dependent potassium channel, KvAP superposes extremely well (rms deviation 2.08 Å) with the voltage-sensor domain of KV1.2/2.1 15. The channel in the Kv 1.2/2.1 structure is in the open or activated conformation, thus indicating that the isolated voltage-sensor domain of KvAP is likely to be in the activated state.
Thus, the scheme II provides a plausible framework to rationalize the effects of class II mutations on the voltage-sensor movement and pore domain with a change in a single parameter. The three schemes that we consider here encompass all the possible ways two domains can be positively coupled to each other. If we assume that a single mutation alters only one thermodynamic term, then derivative analysis demonstrates that the class II mutations must modify interactions between the voltage-sensor and pore, both in the resting and activated states. This suggests that class II mutations play a unique role in coupling the voltage-sensor to the pore domain.
Interestingly, many of our class II mutants map close to the sites in the Shaker potassium channel that cause opposite shifts in the voltage-sensor activation and pore opening (Supplementary Figure 6). For instance, A1149 maps close to the hinge connecting the S4–S5 linker to S5 segment on the sodium channel. E395 and L398 occupy the equivalent position in the Shaker potassium channel. These residues are also facing the hinge connecting the S4 with S4–S5 linker in the neighboring subunit and may be involved in intersubunit coupling. Class II mutants F1278 and N1281 in S6 helix of domain III and map in vicinity of sites L472 in the Shaker potassium channels. S479 in the Shaker channel is in the same region as V1286 in the domain III of the sodium channel. These findings suggest that many of the structural mechanisms involved in electromechanical coupling are likely to be conserved between the members of voltage-gated ion channel superfamily.
Intriguingly, the majority of mutants with class II phenotype in the Shaker channel occur near the hinge regions (supplementary Figure XX). Although there is an uncertainty in the sodium channel due to a lack of experimentally derived structure, a surprisingly large number of the class II mutations were also localized to the hinge regions connecting the various helices such as the S4, S4–S5 linker, S5 segments in the sodium channel homology model. Calculations of angle bending potentials in hinge regions of lysozyme 42,43 and a number of other proteins 44–47 have shown that these potential functions are parabolic with respect to angular displacement. Thus, bending at the hinges generates a restoring force, which acts on the neighboring segments and thereby propagates conformational change. The conserved interactions, by definition, are elastic interactions and hinge-bending potentials will be emblematic of such interactions. Introduction of bulky tryptophan groups at or near the hinge regions may prevent hinge-bending motions and may account for the large number of class II hits at those regions. Based on these findings and previous studies on the Shaker channel, we speculate that the hinge regions play an important role in conformational coupling between the voltage-sensor and pore in voltage-dependent ion channels. The importance of hinge residues in electromechanical coupling suggests that the conformational propagation underlying voltage-dependent gating is likely to involve rigid body motions of helices rather than a wavelike spread such as those envisioned for the ligand gated ion channels 48–50.
Tryptophan and alanine mutations in this study were all generated in a mutant (L1115C) rat skeletal muscle sodium channel (rNav 1.4) background. The region encompassing domain III was subcloned into the pBluescript (KS+) cloning vector by digesting the L1115C mutant sodium channel in the PBSTA vector with the enzymes Aat II and BspEI. A unique silent restriction site (Bam HI) was inserted into this region by overlap PCR method. Tryptophan or alanine mutations were generated by using QUIKCHANGE mutagenesis kit (Stratagene, CA) or by overlap PCR method. Mutations were confirmed by sequencing the entire cloning cassette. The mutated domain III region was reinserted back in to the original vector. The mutant clones were screened by using the unique restriction enzyme. For cRNA preparation, the mutant plasmids were linearized by Not I digestion. cRNA was generated by in vitro transcription using T7 RNA polymerase kit (Ambion, Inc., CA).
Ionic conductance measurements were performed with two electrode voltage clamp (OC-725C Oocyte Clamp, Warner, CT) at room temperature (22 ± 3°C). The oocytes were clamped at a holding potential of −80 mV for at least 5 min prior to recording. All current recordings were obtained within 12–36 hours of cRNA injection. The currents were recorded in an external solution containing 105 mM Na-Mes. In many instances, tetrodotoxin (Sigma Chemical Co.) was added to a final concentration of 10–50 nM in the external solution to obtain better voltage control of inward sodium currents.
Fluorescence signals were typically recorded 3–5 days after injection. To minimize series resistance errors, the external and internal solution contained 105 mM NMG-Mes 51. The channel was labeled by incubating the oocytes with 10 μM TMRM (10mM stock in DMSO) in a depolarizing solution for 30 minutes on ice (XX). After washing the excess unbound label, the oocytes were stored in dark in 105 mM NMG containing external solution. The fluorescence recordings were obtained in two setups both combining a cut-open oocyte voltage-clamp (Dagan Corp. MN) with an upright epifluorescence microscope (Olympus; BX50WI) using a LumPlanFL 40X water-immersion objective (0.8 N.A.). In one setup, the light source was a 150 W tungsten lamp (Newport Corp., CT) supplied by low noise power supply (Kepco, NY) and in the other, a 150 W super-quiet xenon arc lamp (Hamamatsu Inc., Japan) was used. The photocurrents generated by PIN photodiodes in response to emitted light were amplified using either an Axopatch 1B patch clamp amplifier (Axon Instruments, CA) or a Photomax-200 amplifier (Dagan Corp., MN).
Both ionic conductance and fluorescence signals were sampled at 250 kHz with a Digidata 1322A and 1440 interface, respectively (MDS Analytical Technologies, CA). For ionic current measurements, linear and membrane capacity currents were subtracted online using a P/4 procedure with a subtraction holding potential at −120 mV. The ionic current signals were low-pass filtered at 1 kHz. The interval between the each potential trace was 70 ms. Fluorescence traces were recorded without subtraction and low pass filtered at 10 KHz. Each fluorescence trace represents an average of 10 recordings with 1 s interpulse interval. Electrophysiology and fluorescence data were acquired using PClamp 10 acquisition software (MDS Analytical Technologies, CA), and analysis was performed using Clampfit (MDS Analytical Technologies, CA) and Excel (Microsoft, WA). The figures were prepared in Origin (MicroCal, MA), Illustrator (Microsoft, WA) or PYMOL (DeLano Scientific LLC). Statistical analysis was performed using one-way analysis of variance (ANOVA) followed by the post-hoc Dunnett’s test using GraphPad Prism (San Diego, CA).
Conductance-voltage and fluorescence-voltage relationships were fitted to a single Boltzmann function:
where z is the valence and V1/2 is the half maximal voltage. K is the Boltzmann constant, T is the temperature and e is the electronic charge.
The homology model of the sodium channel domains was developed using the crystal structure of the KV1.2/2.1 chimera (PDB ID code 2R9RB). The sequence of the putative sodium channel domains were first aligned to that of the KV1.2/2.1 chimera using ClustalW 52, operated on default parameters (Supplementary figure 3). A separate structure-sequence alignment was also performed using MODELLER 9v4 53. These alignments were combined to generate the final alignment. The homology models of the four sodium channel domains were developed using MODELLER 9v4. For each domain, 100 candidate models were generated. The final model for each domain was selected on the basis of the DOPE (discrete optimized potential energy) scores and the stereochemical quality of the candidate models, as assessed from the Ramachandran plots. The models of the four domains were organized symmetrically around the axis of K+ conduction by using appropriate transformation of coordinates to generate the final tetrameric structure. Energy minimization of this tetrameric structure was performed using the ZMM program 54. No steric clashes were observed within the domains or at the interface between the domains. In the above modeling exercise, the loop regions between the S5 helices and the pore helices were not modeled. This loop which is longest in DI (and relatively shorter in DIII) forms a gap aligned block in the multiple sequence alignment. Moreover, since these loops were not structurally important for the purpose of our study, their modeling was not considered necessary.
This project was supported by funds from the National Institutes of Health (GM084140-01), AHA Scientist Development Award (0535214N), and the Shaw Scientist award to B. C. We thank Dr. Marcel Goldschen for help with preparing the supplementary figure 6 and Kristina Schuldt for excellent technical assistance. We also thank the members of Chanda laboratory for their comments and discussions.
AUTHOR CONTRIBUTIONSY.M. and M.A-.M. contributed equally to this work. Y.M. and M.A-.M. performed all the experiments and analyzed the data. S.C. created the structural model. S.C. and B.C generated the coupling models. Y.M, M.A-.M and B.C wrote the paper.