The Nav1. 7 channel critically contributes to the excitability of sensory neurons, and gain-of-function mutations of this channel have been shown to cause inherited erythromelalgia (IEM) with neuropathic pain. In this study, we report a case of a severe phenotype of IEM caused by p. V1316A mutation in the Nav1. 7 channel. Mechanistically, we first demonstrate that the Navβ4 peptide acts as a gating modifier rather than an open channel blocker competing with the inactivating peptide to give rise to resurgent currents in the Nav1. 7 channel. Moreover, there are two distinct open and two corresponding fast inactivated states in the genesis of resurgent Na+ currents. One is responsible for the resurgent route and practically existent only in the presence of Navβ4 peptide, whereas the other is responsible for the “silent” route of recovery from inactivation. In this regard, the p. V1316A mutation makes hyperpolarization shift in the activation curve, and depolarization shift in the inactivation curve, vividly uncoupling inactivation from activation. In terms of molecular gating operation, the most important changes caused by the p. V1316A mutation are both acceleration of the transition from the inactivated states to the activated states and deceleration of the reverse transition, resulting in much larger sustained as well as resurgent Na+ currents. In summary, the genesis of the resurgent currents in the Nav1. 7 channel is ascribable to the transient existence of a distinct and novel open state promoted by the Navβ4 peptide. In addition, S4–5 linker in domain III where V1316 is located seems to play a critical role in activation–inactivation coupling, chiefly via direct modulation of the transitional kinetics between the open and the inactivated states. The sustained and resurgent Na+ currents may therefore be correlatively enhanced by specific mutations involving this linker and relevant regions, and thus marked hyperexcitability in corresponding neural tissues as well as IEM symptomatology.
Voltage-gated sodium channels, such as the Nav1. 7 channel (~1,980 amino acids; ~225 kDa), are heteromeric protein complexes composed chiefly of a pore loop-forming alpha subunit, which is a large transmembrane protein containing four similar domains, and an auxiliary beta subunit [1–3]. The Nav1. 7 channel is abundantly expressed in the neurons of trigeminal, sympathetic, and dorsal root ganglia (DRG) [2,4, 5]. Mutations in the Nav1. 7 channel may result in severe disorders involving the peripheral nervous system such as paroxysmal extreme pain disorder (PEPD, OMIM 167400), congenital insensitivity to pain (CIP, OMIM 243000), and inherited erythromelalgia (IEM, OMIM 133020) [6–13]. Recently, two sporadic cases of the p. V1316A (located in the DIII/S4–5 linker) mutation were reported [14,15], demonstrating a very severe form of IEM, which is characterized by extremely enhanced activity in relevant neural tissues [16–18]. It would be desirable to explore the molecular and biophysical basis how this mutant Nav1. 7 channels could be responsible for such exteme neural activities. Structurally, each alpha subunit of the Nav1. 7 channel is composed of four homologous domains (D1–D4) and with six transmembrane segments (S1–S6) in each domain [1–3]. The S4 segment acts as a voltage sensor while the loop between S5 and S6 segments lines at the external of the pore containing the selectivity filter [3,19]. Upon depolarization, the S4 voltage sensor moves outward to pull on the S4–S5 linker in each domain, and then on the bundle-crossing region of the S6 segment (the “S6 gate”) to open the channel gate [1–3]. IEM and PEPD are both characterized by episodes of severe pain and “gain-of-function” mutations of the Nav1. 7 channel [8,9, 12,20], although IEM more likely involves extremities than rectal or ocular areas in clinical considerations. Previous studies have shown that quite a few IEM-causing mutations are located in the S4–S5 linker in each domain, such as p. I234T, p. I848T, p. G856D, p. L858F, p. L858H, p. P1308L, and p. V1316A mutations [11]. Accordingly, these mutations have been reported to induce quite a few gating changes including depolarization shift of the activation curve, slowing of slow deactivation time constants, and increase of the response to slow depolarization of ramp currents [14,15,21–24]. Recently, the T1464I (in D3–D4 linker), M1627K (in S4–S5 linker /D4), and p. A1632E (in S4–S5 linker/D4) mutant channels showed enhanced resurgent Na+ currents [25,26], which may play a critical role in the genesis of repetitive or burst neuronal discharges [25]. Conventionally, the genesis of resurgent currents has been ascribed to direct competition between the Navβ4 peptide and the inactivating peptides for the same open state of the channel, and consequently re-opening of the Navβ4-blocked but not the inactivated channel during the repolarization phase just following a depolarization [27–30]. Why and how the resurgent currents are enhanced by these point mutations, however, have remained unexplored. We investigated the changes in molecular behavior of the p. V1316A mutant channel which causes the severe clinical IEM symptomatology. We found that the Navβ4 peptide acts as a gating modifier rather than an open channel blocker to generate resurgent Na+ currents. There is a novel open state that is responsible for the genesis of resurgent current and is significantly existent only in the presence of Navβ4 peptide. The p. V1316A mutation markedly increases both the sustained and the resurgent Na+ currents, mostly ascribable to the destabilized inactivated states by both acceleration of the transition from the inactivated to the open states and deceleration of the reverse transition. These findings not only account well for the molecular mechanism underlying the most severe clinical presentations of IEM but also strongly implicate that S4–5 linker/D3, where V1316 is located, plays a critical role in the molecular operations of recovery from fast inactivation and thus genesis of resurgent currents in Na+ channels.
An 18-yr-old Taiwanese girl suffered from severe burning pain and reddish erythema on both feet (Fig 1A) with the diagnosis of primary erythromelalgia [15]. Genetic analysis showed a missense mutation (p. V1316A) of human Nav1. 7 (SCN9A) gene. The excruciating pain made her immerse her feet into ice-cold water from time to time. The local skin infection was associated with lesions, such as blisters, and the wound healing was poor despite intensive antibiotic therapy and local surgical debridement. Laser Doppler study revealed a profound reduction in the perfusion unit (P. U.) of skin capillary in the feet at resting state (Fig 1B). The phenomenon was also found at both thermal stimulation and post-stimulation state, indicative of a profound reduction of capillary perfusion in the index patient even in the remission or “pain-free” state. In line with the perfusion study, the baseline skin temperature of the patient was markedly lower than normal control (Fig 1C). The skin temperature reached roughly the same level for both the patient and control subjects during thermal stimulation (44°C). However, the patient’s skin temperature was once again significantly lower than that of the control’s post cooling at room temperature for 10 min (Fig 1C). These findings strongly implicate markedly enhanced sympathetic activities associated with IEM. Clinically, infection and inflammation very much aggravated the neuropathic pain. Unfortunately, she developed bouts of extreme pain followed by hypotension and shock with subsequent mortality despite exhaustive resuscitation. To investigate the molecular functional changes of the p. V1316A mutant Nav1. 7 channel, we first characterized the basic key biophysical properties of the WT channel. Fig 2A shows that the resurgent Na+ currents are present only in the presence, but not in the absence of the Navβ4 peptide when there is a prepulse at +40 mV. In contrast, there are much smaller or little resurgent Na+ currents with a prepulse at 0 mV, either in the presence or in the absence of 0. 1 mM Navβ4 peptide. It has been proposed that the Navβ4 competes with the fast inactivating peptide to generate resurgent currents [27,28,31,32]. If the Navβ4 peptide is indeed an open channel blocker acting fast enough to compete with fast inactivation (and thus keeping the channel in an open but non-conducting state responsible for the genesis of resurgent currents upon subsequent repolarization), then the presence of Navβ4 peptide should accelerate the decay of the transient Na+ currents at least at the +40 mV step depolarization (prepulse), where prominent resurgent currents are generated in the following repolarization phase. The cumulative results in Fig 2B show that the decay time constants of the transient Na+ currents in the WT channel remain essentially the same in either the presence or absence of the Navβ4 peptide at both +40 and 0 mV prepulse. These findings strongly implicate that the genesis of resurgent Na+ currents in the presence of the Navβ4 peptide cannot be explained simply with the conventional model inferring the competition between the open channel blockers Navβ4 peptide and the inactivating peptides, and that the open state responsible for the resurgent currents probably is not the same as that giving rise to the transient Na+ currents in the depolarization prepulse. Fig 2C and 2D further shows that the relative sustained Na+ currents (the ratio between sustained and peak of transient Na+ currents) also remain the same whether the Navβ4 peptide is present or not. The unaltered sustained and decay kinetics of transient Na+ currents are further substantiated by the essentially superimposable average Na+ currents in the presence and absence of 0. 1 mM Navβ4 peptide (Fig 2E). These findings once more argue against the open channel blocked by Navβ4 peptide and are consistent with the transient nature of the gating states responsible for the genesis of resurgent currents (see Discussion). We therefore explored the role of the Navβ4 peptide as a gating modifier in addition to the open channel pore-blocking effect. Fig 3A and 3B shows that both the activation and inactivation curves of the WT channel are evidently shifted leftward (~12 mV) on the voltage axis in the presence of the Navβ4 peptide. The slope of the curves, however, is essentially unchanged. These results once again support the role of Navβ4 peptide as a gating modifier, which does not alter the voltage dependence of WT channel opening and inactivation but changes the non-electric free energy difference between the deactivated (closed) and the activated conformations by a fixed amount. Fig 4 compares the basic gating behaviors of the p. V1316A mutant with the WT Nav1. 7 channel. The activation and inactivation curves of p. V1316A mutant channel are negatively and positively shifted in the voltage axis, respectively. On the other hand, the slope of both gating curves is again grossly unchanged or only minimally changed by the mutation (Fig 4B). These results indicate that p. V1316A mutant channel can be more activated but less inactivated especially at a specific range of membrane potentials. Fig 4C shows the presumable range of this predicted “window” or sustained current, which could be defined by the area under both the activation and inactivation curves. Consistently, the increase in sustained currents in the p. V1316A mutant channel could be demonstrated with direct measurement of the late currents during depolarizing pulse (Fig 4D), and may therefore account for part of the origin of nerve hyperexcitability in IEM. In theory, the sustained Na+ currents may signal a decrease in the inactivation rate (slowed transition from the open to the inactivation states) or an increase in the reverse rate (accelerated transition from the inactivated to the open state) or both. Fig 4E shows that in the p. V1316A mutant channel, the decay of transient currents (the macroscopic inactivation rate) is 50% slowed by the Navβ4 peptide. On the other hand, only in the presence of the Navβ4 peptide, not in the absence of the Navβ4 peptide, the macroscopic inactivation rates are ~50% slowed by the p. V1316A mutation. In other words, the Navβ4 peptide slows rather than accelerates decay of transient currents in the p. V1316A mutant channel, and the p. V1316A mutation slows the macroscopic inactivation rates of the channel only in the presence but not in the absence of the Navβ4 peptide, although the Navβ4 peptide by itself has no apparent effect on the inactivation kinetics of the WT channel (Fig 1B). These findings further substantiate that both the Navβ4 peptide and the p. V1316A mutation are gating modifiers of the Nav1. 7 channel, and p. V1316A mutation may at least decrease the inactivation rate of the channel. Resurgent Na+ currents are only discernible in the presence of Navβ4 peptide, which is more likely a gating modifier than a fast open channel blocker (Figs 2 and 3). The voltage dependence of the kinetics and amplitude of the resurgent current were thus further assessed. Fig 5 shows that prominent resurgent currents are elicited at 0~–120 mV following a prepulse at depolarization of +40 mV but not 0 mV in both WT and p. V1316A mutant channels, again demonstrating the requirement of much more positive potentials to elicit resurgent than transient currents (see also Fig 6C below). We then investigated the possibility of accelerated transitions from the inactivated to the open states, which would also increase the resurgent in addition to the sustained Na+ currents (see Discussion). Fig 5C shows that the resurgent currents indeed are significantly larger in the p. V1316A mutant channel than in the WT channel in the presence of the Navβ4 peptide. These results suggest that DIII/S4–5 linker, where p. V1316A mutation is located, may play a critical role in the genesis of resurgent currents by modulation of the kinetics of transitions between the correlative inactivated and open states of the channel in both directions. These altered gating properties by the p. V1316A mutation could be in turn responsible for the extremely heightened neural activities and clinical manifestations of IEM. We then investigated the voltage-dependent occupation of the open state that gives rise to the resurgent Na+ currents. Depolarization prepulses with 3 and 30 ms yield similar “activation curves” of the resurgent currents (Fig 6A and 6B), although the 3 ms prepulses tend to result in higher proportion of the open state toward the more negative voltages (see below), suggesting that most of the open state responsible for the resurgent currents is well accessible with a 3 ms depolarization. The most striking features of the activation of the resurgent currents are the much shallower voltage dependence but a marked positive shift in the voltage axis if compared to the activation curves of the transient currents (Fig 6C). This is not compatible with the evident negative shift of the activation/inactivation curves by Navβ4 peptide (Fig 3B) based on the conventional model. If the resurgent particle or Navβ4 could effectively bind to the open channel pore only upon that strong depolarizing prepulses and thus gives rise to resurgent current at the following repolarization phase, then the activation/inactivation curves which are located at a more negative range in the voltage axis should not be effectively shifted (especially shifted to more negative potentials). On the other hand, the activation curve of resurgent currents in the p. V1316A mutant channel is very similar to or even slightly positively shifted than that in the WT channel. Given the significant increase of both sustained and resurgent currents by the p. V1316A mutation (Figs 4D and 5C), this finding is again difficult to envisage with the conventional model. The competition between resurgent and inactivating particles should become more favorable for the resurgent particle because of significantly decreased tendency of inactivating particle binding by the mutation (so that the resurgent and sustained current could be both increased). The relatively weakened competition of inactivating particle more likely should shift the activation curve of resurgent currents negatively (rather than unchanged or even shifted positively). The characteristics of the activation curve of the resurgent currents are thus not in accordance with the conventional pore-block model, but much more consistent with the proposed new two-open-state or gating modification model. These findings not only lend strong support to a new open state responsible for the resurgent currents, but also suggest a distinct position of the gating voltage sensors of the “resurgent” open state from the “conventional” one (see Discussion). In addition to amplitude, we also investigated the changes in the kinetics of resurgent currents caused by the p. V1316A mutation (Fig 7A and 7B). The kinetics of resurgent current decay at –60 mV remain unaltered over a wide range of prepulse depolarization (i. e. , +40~+160 mV). In contrast, the time to peak resurgent currents is significantly longer in the p. V1316A mutant than in the WT channel (Fig 7B and 7C). This could be consistent with the increased sustained and resurgent currents in the p. V1316A mutant channel, which may be chiefly and straightforwardly explained by an acceleration transition from the inactivation to open state (i. e. , “destabilized” inactivated states). Consistent with the findings in Fig 5, the “re-opening” of the inactivated channel then very likely involves conformational changes of the domain III, S4–S5 linker where V1316 is located (see Discussion). Fig 7D further shows that the kinetics of decay of the resurgent currents are very similarly accelerated by membrane hyperpolarization in both WT and p. V1316A mutant channels. Interestingly, the kinetics of the decay phase of the resurgent currents are always markedly (~10-fold) slower than that of the deactivating tails of the transient currents (Fig 7D–7F). Given the very short time to peak resurgent currents (<~1 ms in the WT channel, <~2 ms in the p. V1316A mutant channel, Fig 7C), this finding is hard to envisage with the conventional model, which has only one open state. This finding further strengthens the views that there are two distinct open states responsible for the transient and resurgent currents, respectively, and that the exit rates from the open state responsible for the resurgent currents remain unchanged by the p. V1316A mutation (see Discussion). In view of the differences in the activation curves of the resurgent currents between 3 ms and 30 ms depolarization prepulses (Fig 6B), we examined the changes in the resurgent currents with a gradually lengthened prepulse. It is intriguing that the resurgent current gets smaller with lengthening of the depolarization prepulse in the WT and p. V1316A mutant channels (Fig 8A and 8B). These findings indicate that with depolarization prepulse 3 ms in length, the WT and p. V1316A mutant channels have not yet reached a steady-state distribution. Moreover, it is easy to reach the new “resurgent” open and the corresponding inactivated states in 3 ms, but subsequent distribution of the channel protein favors the conventional inactivated state, which does not have to go through an open state to deactivate during the following repolarization (see Discussion). The voltage dependence of the kinetics and the “steady-state” relative residual resurgent currents are quite the same in the WT and p. V1316A mutant channels, although the absolute speed of decay is ~2–3-fold slower in the p. V1316A mutant than in the WT channel (Fig 8C and 8D).
We have seen that in the Nav1. 7 channel, resurgent currents are discernible only in the presence of the Navβ4 peptide (Fig 2). According to the conventional model of genesis of resurgent currents [27–30], this is ascribable to the competition between the Navβ4 peptide and the inactivating peptide for the open channel pore (Fig 9A, scheme 1). However, we have also seen that amplitude of the sustained currents and the kinetics of decay of the transient currents (at prepulses to +40 mV where resurgent currents could be effectively generated) remain unchanged with the addition of Navβ4 peptide (Fig 2B and 2C). In the p. V1316A mutant channel, the addition of Navβ4 peptide even significantly slows the decay of the transient currents (Fig 4E), although the subsequent resurgent currents are larger than that in the wild-type channel (Fig 5C). Similar findings can be seen in the Nav1. 1 channel [33], whereas acceleration of the decay phase of the transient current, to our knowledge, has never been reported. These findings are incompatible with the conventional model, which should show an increase in the rate of macroscopic “inactivation” (the rate of decay of the transient currents in the prepulse) and a decrease of sustained currents with the addition of an effective competitive blocker to the system. This would be especially so considering that the resurgent current always very quickly and effectively reaches its peak amplitude within a prepulse just 3 ms in length (Figs 6 and 8). Quite a few other findings may further argue against the hypothesis that the Navβ4 peptide competes with the inactivating peptide for the same open state to generate the resurgent currents: (1) decay of resurgent currents with lengthening of the depolarizing prepulse and a time constant from a few to few tens of millisecond, which is shorter at less depolarized prepulses (Figs 8B, 8C, 10D and 10E); (2) the very much positive shift activation curve of the resurgent current but effective negative shift of the activation/inactivation curve by the Navβ4 peptide (Figs 3 and 6); and (3) the very short time to peak resurgent currents but evidently much slower decay rates of the resurgent currents than the tail currents at the same negative voltages (Fig 7D and 7F). In addition, the major findings with the p. V1316A mutation, such as the concomitant increase of resurgent and sustained currents, the slowed decay (inactivation) phase of transient currents only in the presence but not absence of Navβ4 peptide, and the lengthened time to peak but essentially unaltered decay kinetics as well as the (quite positively located) activation curve of resurgent currents (Figs 4D, 4E, 6C, 7B–7D, 10A and 10B and S1 Fig), would be quite difficult to envisage with the conventional model. We therefore propose that there are at least two distinct open states (and corresponding inactivated states) of the Nav1. 7 channel, each responsible for the transient (O1) and resurgent (O2) currents, respectively (Figs 6,9B and S1). Significant occupancy of O2 and I2 is possible only in the presence of the Navβ4 peptide. In this regard, the Navβ4 peptide is chiefly a gating modifier, which induces new gating conformations rather than being a pore blocker that competes with the inactivating peptide (Fig 10 and S1 Fig). If I2 is chiefly ascribable to a conformational change in the pore rather than a simplistic pore-blocking phenomenon, it might be interesting to re-visit the role of conformational changes in the conduction pathway in the genesis of I1 (versus a pore-blocking phenomenon based on the ball and chain or hinged-lid models), especially considering the decisive and delicate control of Na+ channel inactivation and conformation of the conduction pathway (i. e. , different single channel conductances) by the movement or positioning of S4 in domain IV (S4/D4) [34–37]. The activation curve of resurgent currents is markedly shifted toward the positive voltages compared with the transient currents (~60 mV difference in Vh) but has a much shallower slope (slope factors, ~30–40 versus ~8, Fig 6C). It seems that only in the presence of appropriate gating modifiers, such as the Navβ4 peptide, may gating voltage sensors move further from their position responsible for O1/I1 to O2/I2, this transition is characterized by an apparent transfer of 0. 7–0. 8 equivalent gating charge and a large non-electrical work, considering the much shallower slope and positively shifted Vh of the activation curve of resurgent currents than that of the transient currents. Most intriguingly, the resurgent currents always reach it peak within 3 ms of prepulse initiation, and then decays with a slightly voltage-dependent time constant of 12–40 ms in the wild-type channel (Fig 8). This cannot be explained with the conventional model, where the binding of the competing inactivating and resurgent peptides should be both very fast (the decay time constant of the macroscopic transient currents is only ~0. 5 ms, Fig 2B) and thus very fast kinetics approaching the steady-state distribution of the (blocker-) bound states. On the other hand, this finding could be easily envisioned with a very fast transition from O1 to O2 during prepulse depolarization based on the new model (scheme 2, Fig 9B). In view of the very short time constant of decay of the macroscopic transient currents (~0. 5 ms, Fig 2B), very fast O1 to O2 is actually reasonable or even mandatory for resurgent currents genesis so that O1 to O2 transition is not negligible compared with O1 to I1 transition for a channel protein at state O1. Therefore, strong rather than weak depolarization (e. g. , +40 versus 0 mV) is necessary for sufficient acceleration of the weakly voltage-dependent O1 to O2 transition to produce discernible resurgent currents. Thus, O2 and I2 are transitional states in the “detour” route (from O1 and finally to I1) favored by strong depolarization [34]. Consistently, the “steady-state” relative resurgent currents (fo; Figs 8 and 10) are larger with stronger prepulse depolarization (e. g. , +60> +40> +20 mV). Under such circumstances, the decrease of peak resurgent currents with lengthening of the prepulse depolarization would signal redistribution of the channel back to O1/I1 conformations from O2/I2 conformations during the prepulse (Figs 8 and 10D). The relatively slow kinetics of this redistribution suggests that at depolarizing prepulses, a substantial number of channels first go from state O1 to O2, then to I2, and finally to I1, with I2 to I1 being the slowest or rate-limiting step, requiring at least tens of milliseconds to complete (Figs 9B and 10A). The unaltered sustained currents in Fig 2C by the Navβ4 peptide could be straightforwardly represented by unaltered distributions between O1 and I1. The meanings of all of the major findings could be numerically recapitulated or illustrated by computer simulation based on the two-open-state model (Figs 10 and S1), further substantiating that there are two distinct open and corresponding fast inactivated states underlying the genesis of transient and resurgent Na+ currents in the human Nav1. 7 channel. The p. V1316A mutant Nav1. 7 channel is characterized by a marked increase in both sustained and resurgent Na+ currents (Figs 4D, 5C, 10B and S1). This is the first quantitative report for a specific IEM–causing mutation which concomitantly increases both sustained and resurgent currents. Moreover, the increase in both currents is evident over a wide range of membrane voltages (Figs 4 and 5), suggesting an effective promotion of repetitive discharges by p. V1316A mutation in different firing patterns (e. g. , burst discharges with a plateau or not). Based on scheme 2 (Fig 9B), both increased sustained and resurgent currents can be chiefly and straightforwardly explained by an accelerated transition from the new inactivated to new open state (increased I2 to O2; i. e. , a “destabilized” inactivated state). Also, the increase in sustained currents may be envisioned by the shift of activation and inactivation curves in opposite directions (Figs 4B, 6C and S1), signaling uncoupling between activation and inactivation by the p. V1316A mutation. These findings would provide imperative insight into the molecular pathogenesis of IEM, and is well in line with the findings that the S4–S5 linker in domain III (S4-5 linker/D3) interacts with S6 and/or S4–S5 linker in domain IV (S6/D4 or S4-5 linker/D4) to make the most important structural element relaying movement of voltage sensors to conformational changes in the activation and inactivation gates at the internal pore mouth, and thus critically contributes to Na+ channel gating including activation/inactivation coupling [14,15,34,38,39]. In this regard, it is especially interesting to note the evident correlation between the voltage-dependent entry kinetics into a subconductance level and the movement of voltage sensor in domain 4 (S4/D4) in an inactivation-deficient Nav1. 4 channel [34], a very intriguing finding providing a strong support for a critical conformation change in the conduction pathway induced by the movement of S4/D4. A natural extension of such a molecular rationale of Na+ channel gating may well underlie the creation of a different open state responsible for the genesis of resurgent currents with 0. 7–0. 8 more charge transfer or slightly further outward movement of S4/D4 (Fig 6). Consistently, homology modeling of Nav1. 7 channel shows that the p. V1316A mutation may increase the inter-residue distances between residue p. V1316A of the S4–S5 linker/D3 and the nearby residues (e. g. , L1610, R1611, V1722, and G1723) of the S5–6 segments in domain IV (Figs 9E, 9F and S2). Here we do not mean to make any rigorous argument based on these distances from modeling. Instead, it is a demonstration of the possibility that the conformational changes associated with p. V1316A mutation may well involve the pivotal points of channel activation/inactivation, and thus be responsible for inactivation weakening and activation/inactivation uncoupling (Fig 4E). IEM and paroxysmal extreme pain disorder (PEPD) are characterized by episodic excruciating regional pain associated with hyperemic swelling (Fig 1), presumably manifestations of sensory (dorsal root ganglia) and sympathetic neuronal hyperexcitability. Electrophysiologically, it has been suggested that IEM mutations tend to make negative shift of the Nav1. 7 channel activation curve and also slow the macroscopic deactivation kinetics, whereas PEPD mutations tend to make a positive shift of the inactivation curve, slow macroscopic activation kinetics, and generate more persistent as well as resurgent Nav1. 7 currents [40–42]. In this study, we demonstrated that the p. V1316A mutation, which causes very prominent and typical clinical features of IEM, resulting in not only negative and positive shift of the activation and inactivation curves, respectively, but also evidently larger sustained and resurgent currents. IEM and PEPD thus may be more a clinical classification than a legitimate electrophysiological differentiation. An “appropriate” amount of increased sustained currents might depolarize the neuron to a level not to inactivate too many Na+ channels and thus easier for an external stimulus to fire an action potential. Once the first action potential is fired, the subsequent discharges are powered by the resurgent Na+ currents. The largest amplitude of resurgent currents with very short preceding depolarization (Figs 8,10D and 10E), would also these currents especially suitable for the genesis of densely repetitive spikes or “burst-like” discharges [41,43]. The dorsal root ganglia sensory neurons and sympathetic ganglion neurons carrying IEM-causing Nav1. 7 mutations then could be hyperexcitable, overly responding to an ordinarily innocuous stimulus with excessive bursts of discharges. The episodes of painful attacks of IEM/PEPD therefore carry some of the essential electrophysiological features of epileptic seizures. From a therapeutic point of view, these “peripheral seizures” in IEM/PEPD patients are most different from classical seizures in the general ineffectiveness of conventional anticonvulsants acting on the Na+ channel [15,42,44]. We have noted that the gating states responsible for resurgent currents (O2 and I2 in scheme 2, Figs 9B and 10A) exist only transiently. The anticonvulsants with slow binding rates onto the inactivated Na+ channel (e. g. , phenytoin, carbamazepine, and lamotrigine [45–47]) therefore may not have enough time to bind to the specific gating states to inhibit resurgent currents. It would be interesting to explore the molecular interaction between different anticonvulsants and the novel “resurgent” gating states (i. e. , O2 and I2) in more detail, and thereby search for more effective pharmacotherapy of IEM and related disorders.
The clinical investigations were approved by the Institutional Review Board (IRB) of National Taiwan University Hospital Ethics Committee (Permit Number: NTUH–REC No. 9461700723). The control subjects for the clinical investigations are six (three female and three male) age-matched individuals without any symptoms and signs in both past and present medical histories. Also, each of the control subjects must have normal results in the nerve conduction velocity study prior to recruitment. The patient and control individuals were allowed to rest quietly in an air-conditioned room (20–25°C, humidity 50%) for at least 20 min before initiating the flowmetry measurement using a laser Doppler flowmeter (PeriFlux, PF3, Perimed, Sweden). The principles underlying laser Doppler flowmetry are described elsewhere [48]. A laser Doppler output signal of more than 90% by flow in the subpapillary vessels was continuously monitored [49]. The signals at baseline, during and after thermal stimulation were recorded. Thermal stimulation was achieved by directly heating with laser probe for 1 min in order to raise the skin temperature to 44°C. Skin blood flow was measured immediately after thermal stimulation and again after a 10 min cooling rest period in the air condition room. All measures were expressed in perfusion unit (P. U.). Theoretically, the flow is determined by the product of blood flow velocity and the number of moving red cell corpuscles within the surface capillaries of the skin. The mean value of the perfusion unit in a defined period was analyzed by Perisoft software (the Perimed analysis program for PeriFlux). Each measurement was repeated six times in both the patient and seven normal age-matched controls. The intra-individual coefficient of variations was ~9. 8%. Just prior to the blood flow test, the temperature of the skin was taken using an electric thermometer (Takara, Japan) to measure the temperature at baseline, during, and after thermal stimulation. Mann-Whitney U-test was used for non-parametric comparison among the measurements of the index patient and the control group. The human Nav1. 7 cDNA (SCN9A gene) was subcloned into pTracer–EF/V5–His vector [15]. The GFP sequence, which is driven by a separate promoter in pTracer vector, allows identification of the transfected cells under fluorescent microscope (Nikon, Inc. Japan). The p. V1316A mutation was made using Quikchange site-directed mutagenesis kit (Stratagene, La Jolla, CA, US), and confirmed by the automatic DNA sequencing (Applied Biosystems, 3730xl DNA, Analyzer Foster, CA, US) [15,50–53]. Chinese hamster ovary (CHO–K1) cells were obtained from Food Industry Research and Development Institute (Hsinchu, Taiwan) [15]. CHO–K1 cells were approved by the Institutional Biosafety and Use Committee (IBUC) and conformed to the ethical information guidelines of the National Taiwan University College of Medicine. CHO–K1 cells were cultured under the standard conditions (37°C, 5% CO2) in F12–K medium (Thermo Fisher Scientific Inc. , US) with 10% fetal bovine serum (FBS, Thermo Fisher Scientific Inc. , US) and 1% Penicillin–Streptomycin (Thermo Fisher Scientific Inc. , US). The cells were transfected with either wild-type (WT) or p. V1316A mutant cDNA constructs of the human Nav1. 7 channel using Lipofectamine 3000 (Thermo Fisher Scientific Inc. , US) [15]. The cells were then maintained at 37°C, 5% CO2 incubator before electrophysiological recordings. For recordings, CHO–K1 cells were treated with 1. 0 mg/ml protease type XXIII (Sigma Chemical Co. , St Lois, MO) [54], dissociated in F12–K medium with 10% FBS, and plated on glass cover-slips at 37°C, 5% CO2 incubator for 30–40 min. The enzyme action was terminated by wash with F12–K medium. Usually the recordings were carried out within 2–3 d following DNA transfection, and the isolated CHO–K1 cells were used within 5 h of preparation. Whole-cell patch recordings from CHO–K1 cells expressing WT or p. V1316A mutant channels were conducted at room temperature using the pClamp 6. 0 software and an Axopatch 200A amplifier (Axon Instruments, Inc. Sunnyvale, CA, US). Data were filtered at ~10 kHz and digitized at ~50 kHz through Digidata 1200A interface (Axon Instruments, Inc. Sunnyvale, CA, US). The pipettes were pulled from borosilicate glass (Warner Instruments, MA, US) using a horizontal puller (Zeitz Instruments, Inc. Martinsried, Germany) and fired polished (Narishige scientific instruments, Inc, Japan) [50–53]. The pipette resistance was 1–2 MΩ when filled with intracellular solutions containing (in mM): 75 CsCl, 75 CsF, 5 HEPES, 2 CaCl2 and 2. 5 EGTA, titrate to pH 7. 4 by 1 M CsOH. Whole-cell configuration was obtained and giga-seal formed in the extracellular solution containing (in mM): 150 NaCl, 2 MgCl2,2 CaCl2,10 HEPES, titrated to pH7. 4 by 1 M NaOH. The cell was lifted from the cover-glass and moved in front of a linear array of pipes, each containing a different extracellular solution. The current sweeps in Tetrodotoxin (TTX, Torcris Cookson, Langford, UK) were used for the subtraction of leak currents and capacities to obtain the TTX-sensitive currents for subsequent analysis [55]. The Navβ4 peptide (KKLITFILKKTREK–OH, the C–terminal half of the whole sequence of the Navβ4 subunits [56,57]) remains a standard and effective tool for the investigation of resurgent currents in culture cells or even native neurons in most studies, because the entire β4 subunit is much less effective in the induction of resurgent currents for so-far unknown reasons [32,33]. The Navβ4 peptide (Genomics Bioscience and Technology Co. , Ltd, Taiwan) was dissolved in distilled water to make a 10 mM stock solution and then diluted into intracellular solution for a final concentration of 0. 1 mM [42,58]. Data were analyzed using Clmapfit 9. 0 (Axon Instruments, Inc. Sunnyvale, CA, US), and Sigmaplot 10. 0 software (Systat software, Inc. , Germany). For current-voltage relationship, cells were hold at –120 mV and stepped to a range of test voltages (–140 to +40 mV in 5 mV increments) for 100 ms. Maximal inward currents were plotted as a function of test voltage to generate the current-voltage plot. We made a regression line of the currents points between +10 and +40 mV in the current-voltage plots of each cell [22,35,36,59–64], while the maximal Na+ currents usually appear at a test voltage between 0 and –20 mV. The reversal potential of Na+ (Vrev; Na+) was determined by extrapolating the regression line to the transverse axis, and the maximal Na+ conductance (Gmax) was given by the product of ΔV (the difference between each voltage and the reversal potential) and the slope of the regression line from each individual cell. The measured currents at each voltage were divided by the maximal currents at each voltage in the same cell to give G/Gmax. The activation data were then fitted with a Boltzmann function: G/Gmax = 1/ (1+exp (Vh–V) /k) ), where Gmax is the maximal Na+ conductance, Vh is the potential at which activation is half-maximal, V is the test voltage in mV, and k is the slope factor [22,35,36,59–64]. For the steady-state fast inactivation curve, we measured the currents at a +10 mV test pulse after 100 ms prepulse at different potentials from a holding potential of –120 mV. The current measured from the test pulse was normalized to the maximal currents (I/Imax) in the series from the same cell and plotted against the prepulse voltage. The inactivation data are also fitted with a Boltzmann function I/Imax = 1/ (1+exp (V–Vh) /k) ), where Imax is the maximal Na+ currents, Vh is the potential at which inactivation is half-maximal, V is the prepulse potential in mV, and k is the slope factor. The human Nav1. 7 channel (SCN9A gene) homology modeling was built based on the X-ray crystal structure data of voltage-gated Na+ channel in Arcobacter butzleri (NavAb; PDB code: 3rvy) [19,65,66]. The human Nav1. 7 channel sequence was obtained from the UniProt database (Q15858; http: //www. uniprot. org/). The homology modeling was performed in a similar way to that described in previous studies [50–53]. Sequences of the four domains (DI–DIV) of the human Nav1. 7 channel are very different from those of the NavAb channel [19,65,66]. We first aligned critical S4–S5 linker residues in each domain (D1, V247–R266; D2, F902–M921; D3, G1305–F1334; and D4, Y1591–L1610) with corresponding S4–5 linker residues in the NavAb channel (D1: A112–L131, D2: A351–S371, D3: A593–S613, and D4: A872–S892) according to BLAST research results. We then constructed multiple alignments of S5 and S6 sequences of the four domains (DI–DIV) of the human Nav1. 7 and NavAb channels (S3 Fig). We aligned residues in S5 of each domain (S5: D1, Y356–F376; D2, D933–I943; D3, G1382–Y1402; and D4, L1610–L1630), and S6 of each domain (S6: D1, C725–V746; D2, I1169–M1190; D3, Y1430–F1441; and D4, S1715–I1736) with corresponding S5–6 residues in the NavAb channel (S5: S132–F152, and S6: F201–M222) according to BLAST research results. The crystal structure of the NavAb template sequence that showed higher score in the homologous sequence alignment was chosen as the structure conserved region (SCR) [67]. The identity of S5–6 sequences in each domain between human Nav1. 7 and NavAb channels is greater than 35%, with the root mean square (r. m. s.) of the backbone 0. 01 Å based on a model including only crystallized proteins [68]. Particularly, full sequences of the three interdomain linkers of the human Nav1. 7 channel (DI–DII linker, C740–P789; DII–DIII linker, M1190–Y1239; and DIII–DIV linker, D1447–V1519) were inserted into corresponding positions between different domains (DI–DIV) by homology modeling based on the NavAb channel structure template. The aligned sequences were then presented to Discovery Studio V3. 0 (DS V3. 0) client program to generate relative positions and the secondary structures in the vicinity of selected residues of the human Nav1. 7 channel [50–53]. There were 10 human Nav1. 7 channel homology modelings [50–53], which were subjected to the steepest descent energy minimization with a tolerance of 100 KJ/mol/nm [69], followed by a second conjugate gradient energy minimization with a tolerance of 10 KJ/mol/nm [69] using Discovery Studio V3. 0 (DS V3. 0) client program. The conformation with the lowest free energy potential estimated by the Discrete Optimized Protein Energy (DOPE; ~132,970) Score was chosen for further analysis. The scheme of sodium channel without Navβ4 peptide is based on a previous model [54]. The kinetic parameters were adjusted to fit the experimental data, including activation and inactivation curves of transient currents of WT channels in presence or absence of the Navβ4 peptide, and p. V1316A mutant channels in presence of Navβ4 peptide, peaks of resurgent currents, time to peak of resurgent currents, changes of resurgent currents with lengthening of the depolarization prepulse, and the activation curve of resurgent currents of both WT and p. V1316A mutant channels. The ionic currents are assumed to have Ohmic relationship, and the reversal potential was set as +100 mV. Euler method was used for numeric integration. The maximal time step was set to be 10−6 ms. The computations and graphic constructions were performed with QuB express suite (https: //www. qub. buffalo. edu/) [70] and MATLAB R2015 suite (The MathWorks, Inc. US). All data are presented as mean ± standard error of the mean (S. E. M.). Statistical significance is assessed using Student’s independent t test and accepted at p < 0. 05 (except for Fig 1 in which Mann-Whitney U-test were applied).
The gain-of-function mutation (p. V1316A) of the Nav1. 7 channel causes inherited erythromelalgia (IEM), a disease characterized by extremely enhanced activity in relevant neural tissues that results in neuropathic pain. We found that the p. V1316A mutation alters the basic gating properties of the channel, leading to increased sustained currents during membrane depolarization and resurgent currents during repolarization. Neurons expressing these mutant channels are more difficult to maintain in a hyperpolarized state and are thus more excitable. We demonstrate that there is very likely a distinct set of open/inactivated (O/I) states responsible for the genesis of resurgent currents. We show that the p. V1316A mutation chiefly accelerates the I to O transition in this set, but also decelerates the transitions between different sets of O/I states, to cause the channel gating and cellular excitability changes. Contrary to the conventional view, we find that the Navβ4 peptide, a key element responsible for sizable resurgent currents, does not seem to act as a pore blocker that competes with the inactivation peptide. Instead, we show that it acts as a gating modifier of the Nav1. 7 channel. Thus, the DIII/S4–5 linker, where V1316 is located, may play a critical role not only in O/I coupling but also in the couplings between different sets of O/I in the Nav1. 7 channel.
sequencing techniques
ion channel gating
medicine and health sciences
depolarization
membrane potential
electrophysiology
neuroscience
mutation
ion channels
physiological parameters
body temperature
molecular biology techniques
research and analysis methods
curve fitting
sequence analysis
mathematical functions
sequence alignment
proteins
mathematical and statistical techniques
skin temperature
biophysics
molecular biology
hyperpolarization
physics
biochemistry
point mutation
physiology
genetics
biology and life sciences
physical sciences
neurophysiology
The Biophysical Basis Underlying Gating Changes in the p.V1316A Mutant Nav1.7 Channel and the Molecular Pathogenesis of Inherited Erythromelalgia