The gating pore current, also called omega current, consists of a cation leak through the typically nonconductive voltage-sensor domain (VSD) of voltage-gated ion channels. Although the study of gating pore currents has refined our knowledge of the structure and the function of voltage-gated ion channels, their implication in cardiac disorders has not been established. Two Nav1.5 mutations (R222Q and R225W) located in the VSD are associated with atypical clinical phenotypes involving complex arrhythmias and dilated cardiomyopathy. Using the patch-clamp technique, in silico mutagenesis, and molecular dynamic simulations, we tested the hypothesis that these two mutations may generate gating pore currents, potentially accounting for their clinical phenotypes. Our findings suggest that the gating pore current generated by the R222Q and R225W mutations could constitute the underlying pathological mechanism that links Nav1.5 VSD mutations with human cardiac arrhythmias and dilatation of cardiac chambers.
INTRODUCTION
Voltage-gated Na+ (Nav) channels are transmembrane proteins that are responsible for action potential (AP) initiation in several excitable cells. The SCN5A gene encodes Nav1.5, the main Nav isoform expressed in the heart. This large 2,016–amino acid protein plays a pivotal role in the excitation–contraction process. Nav1.5 mediates the rapid Na+ upstroke, leading to cell depolarization and contraction. It is composed of 24 transmembrane segments organized in four homologous domains (DI–DIV) (Fig. 1). The first four transmembrane segments of each domain (S1–S4) form the voltage-sensor domain (VSD), and the assembly of the transmembrane segments S5 and S6 of the four homologous domains forms the pore domain. Usually, dysfunctions of Nav1.5 channels cause well-described cardiac arrhythmias such as type 3 long QT syndrome, Brugada syndrome, or progressive cardiac conduction defect and sick sinus syndrome (Amin et al., 2010).
The Nav1.5 R222Q and R225W mutations have recently been reported to be associated with the development of an atypical phenotype combining several cardiac arrhythmias and dilatation of cardiac chambers (Bezzina et al., 2003; Cheng et al., 2010; McNair et al., 2011; Laurent et al., 2012; Mann et al., 2012; Nair et al., 2012). These highly conserved R222 and R225 residues are located on the S4 segment of DI of the channel protein (Fig. 1). Although their clinical phenotypes share major similarities, highly divergent biophysical properties have also been observed (Table S1) (Bezzina et al., 2003; Cheng et al., 2010; Laurent et al., 2012; Mann et al., 2012; Nair et al., 2012). Indeed, the R222Q mutant channel exhibits gain of function characteristics, whereas the markedly lower current density observed for the R225W mutant channel leads to loss of function. Thus, so far, no common mechanism explains the similar phenotypes of R222Q and R225W mutations.
Using electrophysiology and molecular dynamics techniques, we investigated the biophysical properties of the R222Q and R225W mutant channels in search of a potentially unifying molecular mechanism that would explain the previously reported phenotypic manifestations. We discovered that the two mutations induce a cationic leak under depolarized condition. Furthermore, after long depolarization periods, we recorded gating pore currents under hyperpolarized conditions, probably caused by a freezing of the VSD after depolarization. As this leak through the VSD is the only alteration of the channel’s properties common to both mutations, we propose that the gating pore currents induced by these mutations contribute to the atypical phenotype that associates complex arrhythmias and cardiac dilatation.
MATERIALS AND METHODS
Cell culture
The electrophysiological experiments were conducted using tsA201 cells, which are modified human embryonic kidney HEK-293 cells stably transfected with simian virus 40 large T antigen. The cells were cultured in high glucose Dulbecco’s modified Eagle’s medium supplemented with 10% FBS, 2 mM l-glutamine, 100 U/ml penicillin G, and 10 mg/ml streptomycin (Gibco) at 37°C in a 5% CO2 humidified atmosphere. They were transfected with WT or mutant Nav1.5 cDNA (5 µg) using the calcium phosphate method as described previously (Deschênes et al., 2000). The Nav1.5 α subunit was cotransfected with the β1 subunit. In brief, 5 µg of each cDNA was mixed with 500 µl of a 250-mM CaCl2 solution to which 500 µl HeBS2x (0.28 M NaCl, 0.05 M HEPES, and 1.5 M Na2HPO4) was added.
Patch-clamp recordings
Macroscopic Na+ currents (α currents) were recorded at room temperature 48–72 h after transfection using the whole-cell configuration of the patch-clamp technique as described previously (Huang et al., 2011). The liquid junction potential between the patch pipette and the bath solution was not corrected. P/4 leak subtraction was used for the α pore current recordings but was not used for the gating pore current recordings. However, for the gating pore currents, offline linear leak subtraction at hyperpolarized voltages was performed to eliminate inherent nonspecific leaks. Currents were recorded with a sampling rate of 83.33 kHz and were low-pass filtered at 5 kHz. The estimation of the unitary conductance of the gating pore currents was based on Ohm’s law, assuming that the Nav1.5 unitary conductance was 22 pS (Gellens et al., 1992) and that the maximum open probability of the gating pore occurred at 40 mV.
Solutions
The compositions of intrapipette and bath solutions used to record the whole-cell α Na+ currents and gating pore currents are described in Table 1. The pHs of all the solutions were adjusted to 7.4 using methanethiosulfonic acid and CsOH, NaOH, or KOH as required. The bath solutions contained 10 mM TEA-Cl to block endogenous potassium currents. 1 mM niflumic acid freshly prepared in 100% ethanol was added to the bath solution before recording the gating pore currents to block endogenous chloride channels. When indicated, 10 µM tetrodotoxin (TTX; Latoxan) was added to the bath solution to block α pore current.
Sequence alignment and homology modeling of the activated state
The procedure was similar to the one described in Gosselin-Badaroudine et al. (2012a). The best available template to build a model of the activated state of the DI VSD of the Nav1.4 channel was identified by a PSI-BLAST search of the Protein Data Bank (PDB) database. We selected the NavAb high resolution structure released in 2011 (Payandeh et al., 2011) because it bears the transmembrane domains with the highest max score. To build the pairwise alignment between Nav1.5 and NavAb, the first 100 human homologous sequences to NavAb in the National Center for Biotechnology Information reference protein database, which includes Nav1.5, were identified using the PSI-BLAST server. A multiple sequence alignment of the 100 sequences was constructed using the ClustalW2 online server (Larkin et al., 2007). The BLOSUM protein weight matrix was used together with standard parameters (opening gap penalty of 10, gap extension penalty of 0.20, and gap distance of 5). The pairwise alignment of NavAb and Nav1.5 was extracted and used to build a 3-D atomistic model of the mammalian channels. The alignment of the VSD transmembrane segments used for the structural modeling is presented in Fig. S1. A standard MODELLER (Eswar et al., 2006) routine was then used to build a comparative model of the DI VSD of the Nav1.5 channel. Because the S3 of NavAb was surprisingly short and because the transmembrane topology prediction software predicted that the S3 helix of the DI of the Nav1.5 channel extended to residue 212, we imposed helicity in our modeler routine.
Sequence alignment and homology modeling of alternate VSD states
The activation/deactivation mechanism of the voltage sensor of voltage-gated ion channels is thought to proceed in a stepwise manner and involves a sliding of the S4 helix relative to a rather static S1–S3 bundle (outward/inward, respectively) (Delemotte et al., 2011; Amaral et al., 2012; Gosselin-Badaroudine et al., 2012a; Henrion et al., 2012; Jensen et al., 2012; Lacroix et al., 2012; Vargas et al., 2012; Yarov-Yarovoy et al., 2012). Each state is believed to be stabilized by several transmembrane salt bridges between S4 gating charges and negative charges provided by the S1–S3 glutamic and aspartic residues and the lipid head group phosphate moieties. During deactivation, the transitions between consecutive states involve a downward ratchet-like movement of S4 in which each of the S4 charges jumps from a negatively charged binding site to the next. This arrangement of positive S4 amino acids and negative S1–S3 charges has been identified as the gating charge transfer center (GCTC) (Tao et al., 2010). Together, S1–S3 segments strongly interact with the S4 segment, thus stabilizing its position either during its transitions or in its resting and activated conformations where R1 and K4 are, respectively, located in the GCTC (Lin et al., 2011; Schwaiger et al., 2012; Gamal El-Din et al., 2013; Lacroix et al., 2014; Moreau et al., 2014b).
To produce models of alternate states of the DI VSD of the Nav1.5 channel, we used a procedure similar to the one described in Wood et al. (2012) used to produce models of two different states of the proton voltage-gated ion channel Hv1. Building an alignment in which S4 is manually shifted by three residues toward the N terminus (corresponding to one helical turn; see Fig. S1) produced a homology model of an intermediate state downward of α, that we call β, using the nomenclature proposed in previous work (Delemotte et al., 2011; Gosselin-Badaroudine et al., 2012a). Producing another alignment in which S4 is manually shifted by six residues relative to the original alignment produced another intermediate state further along deactivation, called γ.
Note that this procedure results in states that are more stable than the steered molecular dynamics procedure devised earlier by our group (Delemotte et al., 2011; Amaral et al., 2012; Gosselin-Badaroudine et al., 2012a).
Molecular dynamic simulations
The three states of the DI VSD of Nav1.5 were then inserted in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer. The systems were equilibrated under normal constant temperature and pressure conditions (298°K, 1 atm.) in a 150-mM NaCl solution. The lipid tails were melted during the first nanosecond, restraining the position of the protein, lipid head groups, water molecules, and ions to their initial position with a strong harmonic potential. To ensure correct reorganization of the lipids and solution, the positions of all the atoms of the channel were then restrained for 2 ns. The side chains were then allowed to reorganize while the backbone was kept restrained for 8 ns. Lastly, a 100-ns unrestrained molecular dynamics simulation of the entire system was conducted, enabling the system to relax. The molecular dynamics simulations were performed using the NAMD2 program (Phillips et al., 2005). Langevin dynamics were applied to keep the temperature (300°K) fixed. The equations of motion were integrated using a multiple time-step algorithm (Izaguirre et al., 1999). Short- and long-range forces were calculated every one and two time steps, respectively, with a time step of 2.0 fs. Chemical bonds between hydrogen and heavy atoms were constrained to their equilibrium value. Long-range electrostatic forces were taken into account using the particle mesh Ewald approach (Darden et al., 1993). The water molecules were described using the TIP3P model (Jorgensen et al., 1983). The simulation used the CHARMM22-CMAP force field with torsional cross-terms for the protein (MacKerell et al., 1998, 2004) and CHARMM36 for the phospholipids (Feller and MacKerrell, 2000). The simulations were performed on the CRAY XT5 KRAKEN at the National Institute for Computer Science and on Owlsnest, Temple University’s supercomputing facility.
Data analysis and statistics
The electrophysiological data were analyzed using Clampfit (pCLAMP v10.0; Molecular Devices) and custom-written MATLAB programs (The MathWorks Inc.). Data are expressed as means ± SEM. When indicated, a t test was performed using SigmaPlot software (Jandel Scientific Software). Differences were considered significant at a P < 0.05 (*), < 0.01 (**), or < 0.001 (***).
Online supplemental material
In Fig. S1, the alignment used to build homology models of the alternate α, β, and γ states of the VSD is shown. In Fig. S2, the water density profiles along the main axis of the VSD in the WT and mutant channels is shown. In Fig. S3, a proposed pathological mechanism related to the presence of a gating pore current is shown. The published biophysical characteristics and clinical phenotypes associated to the R222Q and R225W mutations are shown in Table S1. The salt-bridge patterns in the different models are summarized in Table S2. The PDB files of each model are also available for download as a ZIP file.
RESULTS
The R222Q and R225W mutations induce divergent biophysical properties
The classical biophysical characterization of the mutant channels revealed that the R222Q mutation slightly increased the peak current density, whereas the R225W mutation markedly reduces it (Table 2 and Fig. 2). However, the maximal normalized conductance (Gmax/pF) of the R222Q channel did not increase significantly, indicating that the cell surface channel density was not affected (Fig. 2). The two mutations produced opposite shifts in steady-state activation and inactivation (Fig. 2). The recovery from inactivation was not affected by the R222Q mutation, whereas it was significantly faster for the R225W mutation (Table 2 and Fig. 2). The time to peak and the current decay were both faster for the R222Q channel and slower for the R225W channel than for the WT channel (Fig. 2). Lastly, when normalized to the peak current, the ramp current was higher for both mutations (Fig. 3), although it exhibited different voltage dependence.
The Nav1.5 R222Q and R225W mutant channels conduct a gating pore current at depolarized voltages
The nature of the R to Q or W mutation and their spatial localization on the Nav1.5 protein (Fig. 1, R2, R3, respectively; S4/DI) prompted us to examine the possible existence of a gating pore current through the mutant channels. We recorded nonleak-subtracted currents induced by 80-ms voltage steps from −100 to 40 mV. Fig. 4 shows examples of raw data of WT and mutant channel gating pore currents recorded in Cs+ solutions. The corresponding offline linear leak subtraction showed that a nonlinear component appears on the mutant channels but not on the WT channel (Fig. 4 A, bottom traces). After an offline subtraction of the linear leak, the gating pore currents were normalized to the cell capacitance to obtain the gating pore current density–voltage relationships for the WT and mutant channels (Fig. 4 B). These gating pore currents were also recorded in the presence of 10 µM TTX to completely block the central pore (Fig. 4 B). The resulting gating pore currents were indistinguishable from traces without TTX, demonstrating that the current does not flow through the usual α pore of the protein. Fig. 4 C shows the normalized voltage dependence of gating pores based on the theoretical reversal potential of Cs+ in the recordings solutions. The R222Q and R225W mutations show a V1/2 of 1.2 ± 3.7 mV (n = 7) and 13.5 ± 5.7 mV (n = 7), and a slope factor of −19.0 ± 1.7 mV (n = 7) and −15.1 ± 2.2 mV (n = 7), respectively. The gating pore current density at 40 mV was larger for the R222Q mutant (6.6 ± 1.0 pA/pF, n = 8, without TTX and 6.7 ± 0.6 pA/pF, n = 7, with TTX) than for the WT channel (0.17 ± 0.09 pA/pF, n = 5) and R225W mutant (2.7 ± 0.3 pA/pF, n = 10, without TTX and 2.4 ± 0.4 pA/pF, n = 6, with TTX) channels (Fig. 4 D). The α current at −10-mV and the gating pore current at 40-mV relationship is shown Fig. 4 E. A linear correlation was found (r2 = 0.85 for R222Q channels and r2 = 0.86 for R225W channels), indicating that the amplitude of gating pore currents is correlated to the amplitude of the α pore current. The gating pore created by the mutations allowed the passage of K+, Cs+, and Na+ ions while it excludes the passage of NMDG+ (Fig. 5). Correction for the electrochemical gradient yielded the relative permeation for the ions. The gating pores of the two mutations were less permeable to K+ and Na+ than to Cs+, whereas no permeation was found with NMDG+ (Fig. 5). Indeed, the relative permeability of the mutant channels was approximately two times higher for Cs+ than for K+ (2.2 ± 0.3 for the R222Q channel and 2.0 ± 0.2 for the R225W channel) (Fig. 5 C), whereas the relative permeability of Na+ was lower than for K+ (0.59 ± 0.08 for the R222Q channel and 0.67 ± 0.05 for the R225W channel; Fig. 5 C).
It was recently proposed that arginine mutations inside the VSD would cause a loss of mobility (or a partial freezing) of the S4 after long depolarization periods (Sokolov et al., 2008; Fan et al., 2013; Groome et al., 2014). In the presence of 10 µM TTX, a 500-ms predepolarization was used to mimic the cardiac AP duration. We then applied a ramp protocol from −140 to 0 mV at 0.72 mV/ms to measure gating pore current (see protocol in inset of Fig. 6). At hyperpolarized potentials, gating pore currents were observed after a long depolarization for both R222Q and R225W mutations (−5.6 ± 1.1 pA/pF, n = 5, for R222Q and −3.3 ± 0.6, n = 3, for R225W) and were absent for the WT (−0.9 ± 0.1 pA/pF, n = 5). Such currents could not be observed without the 500-ms predepolarization (Fig. 6). These results show that after long depolarization periods, the mutated S4 segments need more time to resume their resting position and thus partially remain in a conformation, likely enabling gating pore current conduction.
Structural models of the Nav1.5 DI VSD reveal that the R222Q and R225W mutations open an alternative permeation pathway
To gain a better molecular insight into the effect of the mutations, we constructed homology models of the Nav1.5 VSD using the high resolution structure of the partially activated state of NavAb (Payandeh et al., 2011) as a template. We used three alternate alignments in the S4 region (Fig. S1) to build models of three different states (α, β, and γ, with α being the most activated).
As reported previously, the solvent-accessible volume in the four-helix VSD has an hourglass shape in all three states (Fig. 7) (Treptow and Tarek, 2006; Jogini and Roux, 2007; Krepkiy et al., 2009; Delemotte et al., 2010; Ramsey et al., 2010; Gosselin-Badaroudine et al., 2012a; Wood et al., 2012; Li et al., 2014). In each state, a single S4 charge that is located in the constricted region in the center of the VSD, R228 in α, R225 in β, and R222 in γ (Table S2). In the WT models, water density profiles (Fig. S2, dark blue and red lines) indicated that water protrudes from both sides of the membrane, defining two crevices separated by a thin hydrophobic septum. The septum is indicated by a local number density <1.0 located around the conserved aromatic tyrosine at position 168 on S2.
The introduction of the R222Q mutation in the γ state and the R225W mutation in the β state disrupted the salt bridge connecting these residues. This leads to a swelling of the four-helix bundle and the appearance of a solvated pathway through the VSD (Fig. 7). Notably, the R225W mutation opened up a wider hydrophilic funnel than the R222Q mutation (Figs. 7 and S2), potentially accounting for the higher gating pore conductance in this mutant.
DISCUSSION
We hypothesized previously that the gating pore current might be a potential pathological mechanism in the development of complex arrhythmias associated with cardiac dilatation (Gosselin-Badaroudine et al., 2012b, 2014; Moreau et al., 2014a). Several Nav1.5 mutations, including R222Q and R225W, have been reported to be associated with very similar clinical phenotypes. However, there is no consensus on the biophysical defects induced by these mutations. Also, the possibility that these two mutations might generate gating pore currents has never been investigated (Table S1). Here, we report that the R222Q and R225W mutant channels have divergent biophysical properties, which is in agreement with previously published results. No persistent current was observed for both mutant channels (not depicted). Similar mutations are usually associated with pure electrical disorders such as type 3 long QT syndrome and Brugada syndrome (Ruan et al., 2009). However, the clinical phenotypes associated with these mutations are clearly atypical and are probably not caused only by the gain or loss of channel function. Furthermore, such an opposite gain and loss of channel function mechanism has never been associated to identical clinical phenotype. The hypothesis of arrhythmia-induced dilatation has been discarded because of the young age (1-yr-old) patient carrying the R225W mutation (Bezzina and Remme, 2008). Physiologically, both mutations are expected to impact the integrity of the AP in cardiomyocytes. Indeed, both R222Q and R225W mutations have an important impact on the biophysical properties of the α pore of Nav1.5 mutated channels. The gain of function of the R222Q mutant channels would probably favor the premature initiation of AP, notably by increasing the Na+ entry through the α pore (caused by shifts in both activation and inactivation), which increases the depolarizing currents especially at hyperpolarized voltages where WT channels are expected to be closed. The premature ventricular contractions usually observed in patients carrying this mutation strengthen this hypothesis (Laurent et al., 2012; Mann et al., 2012). In contrast, the loss of function of R225W mutant channels caused by reduced membrane expression and shifts in both activation and inactivation would decrease the ability of mutated channels to initiate APs. This would be the cause of atrio-ventricular blocks shown in patients carrying the R225W mutation (Bezzina et al., 2003). Consequently, although both mutations induce divergent biophysical defects, each defect could potentially explain some pathological phenotype aspects. However, even if the impact of the defective α pore cannot be denied, the complexity and the similarity of the pathological phenotypes would probably not be only related to these divergent biophysical defects.
Our results revealed that under depolarized conditions, a new permeation pathway different from the α pore is created. Indeed, similar currents were recorded with or without TTX. These gating pores are selective for larger cations, as shown by their relative permeability (Fig. 5). The gating pore caused by the R222Q mutation activates at less positive potentials than the one caused by the R225W mutation, as shown in Fig. 4 C. These cation leak currents are essentially a gain of function, as they induce a cation leakage that is not present when an arginine occupies the 222 or 225 position. Using Ohm’s law, the unitary conductance of the gating pore was estimated at 72 fS for the R222Q mutant and 126 fS for the R225W mutant, raising the intriguing possibility that the perturbation induced by the R225W mutation is more marked. Indeed, our structural model also provided a molecular explanation for the increased conductance. The increased hydrophobic nature of the tryptophan causes this residue to bury its side chain in the S3–S4–lipid tail interface, twisting slightly the bottom of S4 (counterclockwise when looking from the bottom), whereas glutamine remains preferentially oriented toward the hydrophilic pore (Fig. 7). This leads to the opening of a wider hydrophilic pore throughout the VSD in R225W (as shown Fig. S2). As we suggested previously (Gosselin-Badaroudine et al., 2014; Moreau et al., 2014a), the resulting gating pore current might be a common biophysical defect of the R222Q and R225W mutant channels.
Both R222Q and R225W mutations impact the channel activation and inactivation process, denoting an effect on the voltage-sensing mechanism of the mutated S4 segment. Interestingly, the published data about S4 arginine neutralization indicate a wide variability of the impact of such neutralization on the α pore gating (Chahine et al., 2004; Blanchet and Chahine, 2007; Miceli et al., 2008, 2012; Sokolov et al., 2008; Berger and Isacoff, 2011; Gosselin-Badaroudine et al., 2012b). Given the crucial role of arginines in the voltage-sensing mechanism, it can be expected that their neutralization would impair the S4 movement and thus induce a depolarizing shift of the activation. It is also important to note that the α pore channel opening might be dependent on the S4 movement but also on the coupling between the VSD and the pore. The absence of effect on channel activation observed with several mutations such as R669Q/W (R3) on Nav1.4 DII (Sokolov et al., 2008) might denote their little impact on the S4-sensing mechanism (and that the VSD can be adequately coupled to the α pore of the channel). Oppositely, impacts on channel gating of R222Q and R225W mutations might indicate that either the sensing mechanism or the coupling with the α pore of mutant channels is affected. Finally, such discrepancies on the effect of S4 neutralization could also be a marker of structural differences between channel VSDs. Indeed, although VSDs are composed of similar motifs (positives charges on the S4, the presence of a GCTC in S1–S3), each VSD possesses its own biophysical characteristics and thus its own structural identity.
Gating pores have been associated with the development of periodic paralyses (Sokolov et al., 2007, 2008; Struyk and Cannon, 2007; Moreau et al., 2014a). Periodic paralyses are caused by arginine mutations in S4 of Cav1.1 or Nav1.4, the skeletal muscle homologue of Nav1.5. Interestingly, the Nav1.4 R222W mutation (R2 of DI of S4) has been associated with the development of hypokalemic periodic paralysis (HypoPP) (Matthews et al., 2009; Park and Kim, 2010).
The properties of two HypoPP mouse models also lend credence to the association between the gating pore and the development of HypoPP. These two models also present a dilatation of transverse tubules, indicating that the gating pore current may have an impact on cell morphology (Wu et al., 2011, 2012).
Previous work on Shaker-type and Nav channels revealed that gating pore current pathways open in a state-dependent manner, i.e., in the specific state in which the mutated arginine is in the GCTC and in which the disruption of the salt bridge leads to a junction of the two water crevices (Tombola et al., 2005, 2006, 2007; Sokolov et al., 2007; Delemotte et al., 2010; Gamal El-Din et al., 2010; Gosselin-Badaroudine et al., 2012a; Khalili-Araghi et al., 2012; Moreau et al., 2014b). Two different states fulfill these conditions for the R222Q and R225W mutants (β for R225W and γ for R222Q), raising the possibility that the relative population of the various states at a given transmembrane voltage is also altered by the mutations. Indeed, although the VSD is expected to be in its fully activated state (at 40 mV), gating pore currents are observed for R222Q and R225W mutant channels, indicating that the mutated residue is still located in the GCTC. This suggests that the “fully activated” state for these mutant channels might be different from the “fully activated” state of WT channels.
Similar depolarization-activated gating pores have been shown to be permeant at hyperpolarized potentials after long depolarization periods (Sokolov et al., 2008; Fan et al., 2013; Gamal El-Din et al., 2014; Groome et al., 2014). This phenomenon would be caused by S4 segment immobilization whereas, in WT condition, after depolarization, the S4 segments are expected to rapidly reach their resting conformation (Cha et al., 1999). Our results indeed show that after the depolarization period mimicking the cardiac AP duration, the mutated S4 segments remain trapped in their conductive conformations at hyperpolarized voltages (Fig. 6).
Based on physiological ionic homeostasis, the R222Q and R225W gating pores would cause an outward K+ current and an inward Na+ current during the cardiac AP (Fig. 8). Because of S4 segment immobilization, at the end of the cardiac AP, the gating pore would be temporarily permeable to Na+ ions (Fig. 8). In this case, the Na+ leak would cause a Na+ overload, as observed in HypoPP patients (Jurkat-Rott et al., 2009; Fan et al., 2013). The potential effects of this Na+ overload are schematized in Fig. S3. The Na+ leak may counterbalance Kir channels, limiting repolarizing currents (Tricarico and Camerino, 2011; Moreau et al., 2014a). This Na+ leak may also activate reverse Na+/Ca2+ exchange (Remme and Wilde, 2014), leading to a Ca2+ overload. A Ca2+ overload would modulate Ca2+ affinity of troponins, which would disturb the contractile apparatus (Liu et al., 2012). Ca2+ overload would also disrupt intercellular communication (by blocking connexins) (Bukauskas et al., 2001), and disrupt Ca2+ homeostasis through the ryanodine and IP3 receptors. However, the exact downstream consequences of such gating pores remain to be elucidated. The study of animal models could highlight the exact pathophysiological consequences of cation leaks through gating pores.
In conclusion, our results demonstrate that R2 and R3 mutations in S4 DI of Nav1.5 channels lead to the creation of a gating pore. Such gating pores have already been observed with the Nav1.5 R219H mutation and other Nav1.4 and Cav1.1 mutations associated with Hypo or NormoPP (Gosselin-Badaroudine et al., 2012a; Moreau et al., 2014a). Collectively, our data suggest that gating pores caused by mutations in Nav1.5 VSD could constitute a novel pathological mechanism at the origin of an atypical clinical phenotype combining complex cardiac arrhythmias and dilatation.
Acknowledgments
The authors wish to thank H. Poulin and O. Theriault for many insightful discussions. Thanks are also extended to Valérie Pouliot for her technical assistance and for performing the site-directed mutagenesis and the generation of constructs.
This study was supported by grants from the Heart and Stroke Foundation of Quebec, the Canadian Institutes of Health Research (MOP-111072 and MOP-130373), and the National Science Foundation (major research instrumentation grant CNS-09-58854). L. Delemotte receives funding from the European Union Seventh Framework Program (PIOF-GA-2012-329534) “Voltsens.”
The authors declare no competing financial interests.
Kenton J. Swartz served as editor.