Original Research ARTICLE
Mutations in the Voltage Sensors of Domains I and II of Nav1.5 that are Associated with Arrhythmias and Dilated Cardiomyopathy Generate Gating Pore Currents
- 1Centre de Recherche de l’Institut Universitaire en Santé Mentale de Québec, Quebec City, QC, Canada
- 2Cardiovascular Research Program, VA New York Harbor Healthcare System, Brooklyn, NY, USA
- 3Department of Medicine, Université Laval, Quebec City, QC, Canada
Voltage gated sodium channels (Nav) are transmembrane proteins responsible for action potential initiation. Mutations mainly located in the voltage sensor domain (VSD) of Nav1.5, the cardiac sodium channel, have been associated with the development of arrhythmias combined with dilated cardiomyopathy. Gating pore currents have been observed with three unrelated mutations associated with similar clinical phenotypes. However, gating pores have never been associated with mutations outside the first domain of Nav1.5. The aim of this study was to explore the possibility that gating pore currents might be caused by the Nav1.5 R225P and R814W mutations (R3, S4 in DI and DII, respectively), which are associated with rhythm disturbances and dilated cardiomyopathy. Nav1.5 WT and mutant channels were transiently expressed in tsA201 cells. The biophysical properties of the alpha pore currents and the presence of gating pore currents were investigated using the patch-clamp technique. We confirmed the previously reported gain of function of the alpha pores of the mutant channels, which mainly consisted of increased window currents mostly caused by shifts in the voltage dependence of activation. We also observed gating pore currents associated with the R225P and R814W mutations. This novel permeation pathway was open under depolarized conditions and remained temporarily open at hyperpolarized potentials after depolarization periods. Gating pore currents could represent a molecular basis for the development of uncommon electrical abnormalities and changes in cardiac morphology. We propose that this biophysical defect be routinely evaluated in the case of Nav1.5 mutations on the VSD.
Voltage-gated sodium channels (Nav) are transmembrane proteins responsible for the initiation and propagation of action potentials in many excitable cells. Cardiac action potentials are initiated by Nav1.5, the main Nav isoform in the heart. This large 2016-amino-acid protein is encoded by the SCN5A gene and is composed of 24 transmembrane segments organized in four homologous domains (Gellens et al., 1992). The assembly of the S5–S6 segments of each domain forms the alpha pore, which is responsible for the selective permeation of Na+ ions that underlies the physiological functions of the channel. The S1–S4 segments of each domain make up the voltage sensor domain (VSD), which modulates the activity of the channel through conformational changes induced by variations in the transmembrane voltage. Due to the crucial role of these channels in the generation and regulation of the electrical signals, dysfunctions result in pure arrhythmic phenotypes. For instance, Nav1.5 dysfunctions result in cardiac pathologies such as Brugada syndrome, type 3 Long QT syndrome, sick sinus syndrome, and cardiac conduction disorders with no structural heart disease (Amin et al., 2010).
Nav1.5 mutations have recently been linked to the development of atypical phenotypes associating cardiac arrhythmias and dilated cardiomyopathy (Bezzina et al., 2003; McNair et al., 2011; Laurent et al., 2012; Mann et al., 2012; Nair et al., 2012; Beckermann et al., 2014; Gosselin-Badaroudine et al., 2014; Moreau et al., 2014a, 2015). The exact pathological mechanism is not completely understood. While the clinical phenotypes are similar, some mutations result in a gain of channel function, some in a loss of channel function, and some have no effect on normal channel function (Bezzina et al., 2003; Gosselin-Badaroudine et al., 2012b, 2014; Laurent et al., 2012; Mann et al., 2012; Nair et al., 2012; Moreau et al., 2014a, 2015). Two Nav1.5 mutations located at homologous positions on the S4 segments of DI and DII (R225P and R814W; Figure 1) have been characterized (Nguyen et al., 2008; Beckermann et al., 2014). These mutations were found in patients whose phenotypes feature complex arrhythmias associated with dilated cardiomyopathy. Atrial flutter and non-sustained ventricular tachycardia were reported for the patient carrying the R814W mutation (Olson et al., 2005). The patient Nav1.5 R225P carrier experienced atrial and ventricular tachycardia, 2:1 atrio-ventricular block, ventricular, atrial ectopy and prolonged QTc (Beckermann et al., 2014). All carriers also demonstrate dilated cardiomyopathy with reduced ventricular ejection fraction (Olson et al., 2005; Beckermann et al., 2014). There is currently no consensus on the exact pathological mechanism linking morphological changes in the heart with a channel typically associated with pure arrhythmic phenotypes.
FIGURE 1. Location of the R225P and R814W mutations that target highly conserved residues in VSDs. (A) 2D structure of Nav1.5 channels showing the 24 TM segments organized into 4 domains (DI-DIV), each containing 6 TM. The S1-S4 segments of each domain form the VSD (light gray) while the assembly of the S5–S6 segments (dark gray) form the pore of the channel. The locations of the R225P and R814W mutations are indicated by red circles. (B) 3D crystal structure of the VSD of NavAb channel showing the residues that are conserved in all VSDs. These residues were determined by calculating the Shannon entropy at every position for the 6,652 VSD sequence alignments reported by (Palovcak et al., 2014). Highly conserved residues are in blue while less conserved residues are in red. The S1 segment was removed for purposes of clarity. The positively charged residues of the S4 segment (left) are highly conserved as are the residues forming the GCTCs of the S2 and S3 segments (right and middle respectively).
Interestingly, most of the reported mutations (9 of 13) are located on the positively charged S4 segment of the VSD (McNair et al., 2011; Gosselin-Badaroudine et al., 2012b, 2014; Moreau et al., 2014a, 2015). This segment drives the conformational changes of the VSD in response to changes in voltage (Moreau et al., 2014b). In the course of its movement, the positively charged amino acids of the S4 segment sequentially interact with their negatively charged counterparts on the S1–S3 segments (Delemotte et al., 2011). This arrangement of negatively charged amino acids on the S1–S3 segments is called the gating charge transfer center (GCTC; Tao et al., 2010). Mutations of the positively charged amino acids of the S4 segment disrupt interactions between the S4 segment and the GCTC, resulting in the creation of a new permeation pathway that passes directly through the usually non-conductive VSD. This new permeation pathway is called the gating pore (or omega pore; Sokolov et al., 2005, 2007; Tombola et al., 2005; Struyk and Cannon, 2007; Gosselin-Badaroudine et al., 2012a; Moreau et al., 2014a,b, 2015).
We hypothesized that gating pores might be a common pathological mechanism of several VSD mutations linked to cardiac arrhythmias and dilatation (Gosselin-Badaroudine et al., 2014; Moreau et al., 2014a, 2015). We previously showed that the R219H, R222Q, and R225W mutations in DI of the cardiac Na+ channel cause gating pore currents (Gosselin-Badaroudine et al., 2012b; Moreau et al., 2015).
We used a heterologous expression system (tsA201 cells) to study the R225P and R814W mutations (R3, S4 in DI and DII, respectively; Figure 1; Nguyen et al., 2008; Beckermann et al., 2014) to determine whether gating pores are present. The characterization of both mutations confirmed the previously observed biophysical defects. We also showed that mutations in DI and DII both result in gating pore currents.
We propose that this permeation pathway may participate in the pathological mechanism. We further propose that this pathway should be routinely evaluated in future studies on SCN5A mutations in the VSD of the channel that are associated with arrhythmias and dilated cardiomyopathy.
Materials and Methods
TsA201 cells, which are modified HEK-293 (human embryonic kidney) cells stably transfected with simian virus 40 large T antigen, were used for patch clamp experiments (Moreau et al., 2012). The cells were cultured in high-glucose DMEM supplemented with 10% fetal bovine serum (FBS), 2 mM L-glutamine, 100 U/ml of penicillin G, and 10 mg/ml of streptomycin (GIBCO- BRL Life Technologies, Canada) at 37°C in a 5% CO2 humidified atmosphere. The cells were transiently transfected with WT or mutant Nav1.5 cDNA (5 μg) using the calcium phosphate method as previously described (Deschenes et al., 2000). The Nav1.5 α-subunit was co-transfected with the human β1-subunit. Briefly, 5 μg of each cDNA was mixed with 500 μl of a 250 mM CaCl2 solution to which 500 μl of HeBS2x (0.28 M NaCl, 0.05 M HEPES, 1.5 M Na2HPO4, pH 7.05) was added. Given the size of plasmids and genes inserted, we transfected a ratio of 2 β1 subunits for each alpha subunits.
Macroscopic Na+ currents (alpha currents) were recorded at room temperature 48–72 h after transfection using the whole-cell configuration of the patch-clamp technique, as previously described (Huang et al., 2011). The liquid junction potential between the patch pipette and the bath solution were not corrected. P/4 leak subtraction was only used for the alpha pore current recordings. For the gating pore current recordings, linear leak subtraction at hyperpolarized voltages was performed off-line to eliminate the inherent non-specific leak. Currents were recorded at a sampling rate of 83.33 kHz and were low pass filtered at 5 kHz. All currents were recorded using a pipette solution composed of (in mM): 135 CsF, 5 NaCl, 10 EGTA, and 10 HEPES (pH 7.4). The bath solution was composed of (in mM): 115 NMDG, 20 NaCl, 2 CsCl, 2 Ca(OH)2, 10 HEPES, and 10 TEA-Cl (pH 7.4). TEA-Cl was used to block endogenous potassium currents. Niflumic acid (1 mM) freshly prepared in 100% ethanol was added to the bath solution to block endogenous chloride channels before recording the gating pore currents. When indicated, 10 μM TTX (Latoxan, France) was added to the bath solution to block alpha pore currents.
We used the procedure described by Moreau and coworkers to generate models for the WT and mutated VSDs for both DI and DII (Moreau et al., 2015). Briefly, the high-resolution of the pre-activated NavAb structure released in Payandeh et al. (2011) in which all VSDs were considered to be in the activated position was used as a template (accession code 3RVY). Standard MODELER routines (Eswar et al., 2006) were used to generate comparative models of the DI and DII VSDs of the Nav1.5 channels. Homology models in which the R225 and R814 residues were in the GCTC were generated using a procedure described previously (Wood et al., 2012; Moreau et al., 2015). Since the activation/deactivation mechanism involves a sliding S4 helix relative to a static S1-S3 bundle with discrete steps during which the positive residues sequentially interact with the GCTC (Delemotte et al., 2011; Gosselin-Badaroudine et al., 2012a; Henrion et al., 2012; Jensen et al., 2012; Moreau et al., 2015). The alignment of the S4 segments used for homology modeling was shifted by three amino acids toward the C-terminus in order to align R225 and R814 to the residue designated as K4 by Payandeh et al. (2011). These alignments correspond to a meta-stable state previously named β state (Delemotte et al., 2011). In this state, R225 and R814 interact with the GCTC.
The homology models were then inserted in a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and 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 while restraining the positions of the protein, the lipid head groups, the water molecules, and the ions at their initial positions. Afterward, only the protein was restrained for 2 ns. The restraints were relaxed progressively for 6 ns. Lastly, a 100-ns unrestrained MD simulation of the entire channel was conducted, enabling the system to relax. The MD simulations were carried out using the NAMD2 program. Langevin dynamics were applied to keep the temperature (300 K) fixed. The equations of motion were integrated using a multiple time-step algorithm. Short- and long-range forces were calculated every 1 and 2 time-steps, respectively. The time steps were 2.0-fs long. The simulation used the CHARMM22-CMAP force field with torsional cross-terms for the protein and CHARMM36 for the phospholipids.
Mutant VSDs were created by introducing the mutation into the equilibrated WT VSD using MODELER. They were equilibrated using the same procedure. The simulations were performed on the Calcul Québec/Compute Canada Colosse supercomputer. The operation of the supercomputer is funded by the Canada Foundation for Innovation (CFI), NanoQuébec, RMGA, and the Fonds de recherche du Québec – Nature et technologies (FRQ-NT).
Data Analysis and Statistics
The electrophysiological data were analyzed using Clampfit (pCLAMP v10.0, Molecular Devices) and custom MATLAB programs (The MathWorks, Inc.). Data are expressed as mean ± SEM (standard error of the mean). When indicated, an ANOVA with a Dunnett post hoc test or a t-test was performed. Differences were considered significant at a ∗P < 0.05, ∗∗P < 0.01, or ∗∗∗P < 0.001.
The R225P and R814W Mutations Cause a Gain of Function of the Alpha Current
As previously described (Nguyen et al., 2008; Beckermann et al., 2014), the R225P and R814W mutations both induce a gain of channel function. All experiments are presented in Figure 2 while in Table 1 all values and number of independent experiments are listed (Figure 2 and Table 1). However, the defects induced by the R814W mutation were less pronounced than that of the R225P mutation. Both mutations increased the opening probability at hyperpolarized voltages [left shift of their half voltage (V1/2) of activation, increased slope parameter]. R225P mutant channels displayed a larger slope factor for their steady-state inactivation curve while R814W mutant channels displayed a left shift in the V1/2 of inactivation. Both mutations produced a slower recovery from inactivation, an increase in time to peak, and a slow current decay. These alterations in the biophysical properties of the channel resulted in an increase in channel function through an increase in the window current as well as an increase in channel availability as evidenced by the calculated open probability and the ramp currents (Figure 2).
FIGURE 2. Biophysical characterization of the Nav1.5 mutant channels. The values for the Nav1.5 WT channel are indicated by black symbols, those for the R225P mutant channel by orange symbols, and those for the R814W mutant channel by green symbols. (A) Representative whole-cell current traces of the WT and mutant channels (see protocol in inset). (B) Current density-voltage (I–V) relationships of the WT and mutant channels. (C) Voltage-dependence of steady-state activation of the WT and mutant channels. Activation curves were generated using a standard Boltzmann distribution [G(V)/Gmax = 1/(1+exp(-(V-V1/2)/k))] and gave the V1/2 and k values listed in Table 1. (D) Steady state inactivation of the WT and mutant channels (see protocol in inset). The recorded inactivation values were fitted to a standard Boltzmann equation [I(V)/Imax = 1/(1+exp((V-V1/2)/k))+C] and gave the values listed in Table 1. (E) Recovery from fast inactivation values were obtained using a two-pulse protocol at +30 mV to obtain maximal activation (see protocol in inset). The time constants listed in Table 1 were obtained using a single-exponential function: (Afast x (1 – exp(-t/τfast)) + C). (F) The times to peak of the WT and mutant channels were used to evaluate activation kinetics. The times to peak were measured from the same current traces used to construct the I-V relationship (A). (G) The time constants of fast inactivation decay were plotted as a function of voltage for the WT and mutant channels. The time constants were obtained using a simple-exponential function: (A x (exp(-t/τ) + C). Normalized raw data shown in the inset illustrate the current decay kinetics. (H) The overlap between activation and inactivation defines the window current. The predicted window current was obtained using the following equation: (1/(1 + exp((V1/2activation – V)/kactivation)) x ((1 – C)/(1 + exp((V – V1/2inactivation)/kinactivation)) + C). (I) Ramp protocols (see protocol in inset) were imposed (0.64 mV/ms) to study the window current. Ramp current traces were normalized to the alpha peak current. (J) Histogram showing the peak window current normalized to the alpha peak current (% of peak current). ∗Difference vs. Nav1.5/WT using ANOVA (∗∗∗P < 0.001).
The R225P and R814W Mutations both Create Gating Pores that are Activated at Depolarized Potentials
Gating pores are an alternative permeation pathway created in the usually non-conductive VSDs of voltage-gated ion channels by the disruption of interactions between the S4 segment and the GCTC (Moreau et al., 2014a,b, 2015). Due to the nature and location of the R225P and R814W mutations (R3 and S4 in DI and DII, respectively), we hypothesized that they might generate gating pore currents. Unlike WT channels, gating pore currents were induced by the R225P and R814W mutations when 80-ms voltage steps ranging from -100 to +40 mV in 5 mV increments were applied. Because P/4 leak subtraction was not used to record the gating pore currents, linear leak subtraction at hyperpolarized voltages, when gating pores should be closed, was performed off-line to eliminate the non-specific leak inherent to the patch-clamp method (Figure 3, bottom panels). The gating pore currents created by the R225P mutation were larger than those created by the R814W mutation [at +40 mV: 11.8 ± 1.0 pA/pF for R225P (n = 5), 3.4 ± 0.4 pA/pF for R814W (n = 4), and 0.4 ± 0.2 pA/pF for WT (n = 4); Figure 4]. Tetrodotoxin (TTX), an alpha pore blocker, did not modify the gating pore currents induced by the R225P [at +40 mV: 11.8 ± 1.0 pA/pF (n = 5) without TTX vs. 10.9 ± 1.1 pA/pF (n = 6) with 10 μM TTX] and R814W mutant channels [at +40 mV 3.4 ± 0.4 pA/pF (n = 4) without TTX vs. 3.2 ± 0.3 pA/pF (n = 3) with 10 μM TTX; Figure 4B]. This confirmed that these cation leak currents do not flow through the alpha pore of the channel but rather through an alternative pathway created by the mutations. The gating pores investigated here were activated by membrane depolarizations. The normalized gating pore currents (I/Imax) showed that the R225P and R814W gating pores have similar voltage dependence (Figure 4C).
FIGURE 3. Gating pore currents of the Nav1.5 R225P and R814W mutant channels. Gating pore currents were recorded in absence of TTX. The cells were held at -100 mV and currents were recorded using a voltage-step protocol from -100 to +40 mV in 5 mV increments. The top panels show examples of raw traces of gating pore currents. The currents are plotted as a function of voltage in the bottom panels. Linear non-specific leaks are indicated by dotted lines.
FIGURE 4. Voltage dependence of gating pore currents. (A) Current density-voltage relationships of gating pore currents recorded for the WT, R225P, and R814W channels are shown as filled symbols in the absence of TTX (n = 5 for R225P, n = 4 for R814W, and n = 4 for WT channels) and as empty symbols in the presence of 10 μM TTX (n = 6 for R225P, n = 3 for R814W, and n = 5 for WT channels). (B) Histograms summarizing the gating pore current densities at +40 mV for the WT, R225P, and R814W channels (0.4 ± 0.2 pA/pF, 11.8 ± 1.0 pA/pF, and 3.4 ± 0.4 pA/pF, respectively). No differences were observed with or without TTX for any of the conditions. (C) Normalized gating pore currents (I/Imax) exhibiting similar voltage dependence for both the R225P and R814W mutant channels. ∗Difference vs. Nav1.5/WT (∗∗∗P < 0.001).
Mutations at the bottom of the S4 segment near the cytoplasm amplify the naturally occurring S4 segment immobilization due to long depolarization periods (Sokolov et al., 2008; Fan et al., 2013; Groome et al., 2014; Moreau et al., 2014b, 2015). Gating pore currents were thus recorded (in the presence of 10 μM TTX) using a ramp protocol ranging from -140 to 0 mV at 0.72 mV/ms after a 500-ms pre-depolarization to +40 mV and were compared to those elicited by the same ramp protocol without pre-depolarization (see protocol shown as inset in Figure 5). To minimize the impact of capacitance artifacts, gating pore currents were measured from -135 to 0 mV. The I/V curves were constructed by averaging current values every 5 mV. After linear leak subtraction, gating pore currents at hyperpolarized potentials were observed only after a 500-ms pre-depolarization for both mutants [at -135 mV: -4.6 ± 0.8 pA/pF (n = 6) vs. -0.3 ± 0.1 pA/pF (n = 5) without pre-depolarization for R225P channels and -2.6 ± 0.4 pA/pF (n = 4) vs. -0.4 ± 0.2 pA/pF (n = 4) for R814W channels]. As expected, no currents were recorded with WT channels (-0.8 ± 0.1 pA/pF, n = 5) with a 500-ms pre-depolarization.
FIGURE 5. Gating pore currents after long depolarizations. (A) Currents generated by ramp pulses (see protocols in inset) for the WT (left), R225P (middle), and R814W (right) channels. The I–V curves were constructed by averaging current values every 5 mV. Voltages were calculated using the known time course of the ramp protocol. A linear leak subtraction around -75 to -45 mV was performed to eliminate the inherent linear leak. The insets show the currents in response to the ramp protocols. Dashed lines indicate the current obtained without pre-depolarization while solid lines indicate the response after a 500-ms pre-depolarization at +40 mV. (B) Histogram summarizing the inward gating pore current density at -135 mV recorded with or without a 500-ms pre-depolarization. ∗ Indicates a statistical difference between the condition with or without pre-depolarization tested using a t-test (∗∗P < 0.01); # Indicates differences with the WT channel (with the 500-ms pre-depolarization) tested using an ANOVA. (#P < 0.05, ###P < 0.001) .
Nav1.5 mutations are involved in pure electrical disorders such as Brugada syndrome, type 3 long QT syndrome, and cardiac conduction defects with no structural heart diseases (Amin et al., 2010; Moreau et al., 2012, 2013). However, SCN5A mutations have been recently associated with cardiac dysfunctions combining arrhythmias and dilatation (McNair et al., 2011; Moreau et al., 2014a). How this gene could be at the origin of a structural heart disease is not fully understood. Gating pores are alternative permeation pathways created in the VSDs of ion channels by the disruption of interactions between the GCTC and the S4 segment (Moreau et al., 2014b). These gating pores notably appear as a result of mutations of the charged residues of the S4 segment. Such pores are known to be permeant to cations with a selectivity sequence that is more permissive for larger ions (Moreau et al., 2014a,b, 2015).
Several mutations on the S4 segment of the Nav1.5 channel have been studied due to their association with the pathological development of arrhythmias and dilated cardiomyopathy (Bezzina et al., 2003; Nguyen et al., 2008; Laurent et al., 2012; Mann et al., 2012; Nair et al., 2012; Beckermann et al., 2014; Moreau et al., 2014a, 2015). However, no consensus has emerged on the pathological process leading to this atypical clinical phenotype. In addition, some mutations with similar clinical phenotypes cause different biophysical defects (gain or loss of function; Nguyen et al., 2008; Gosselin-Badaroudine et al., 2014; Moreau et al., 2014a, 2015). We propose that gating pore currents may be a common pathological mechanism linking all mutations on the VSDs of the Nav1.5 channel (Gosselin-Badaroudine et al., 2012b, 2014; Moreau et al., 2014a, 2015). We previously showed that gating pore currents are generated by three mutations found in unrelated patients (R219H, R222Q, and R225W) that cause complex arrhythmias and dilated cardiomyopathy (Gosselin-Badaroudine et al., 2012b; Moreau et al., 2015). All these mutations are located in the S4 segment of DI of the Nav1.5 channel. The creation of gating pores due to mutations in the DII of the Nav1.5 channel has never been investigated before. Such investigations remain of critical importance, because the association between a mutation and the creation of a gating pore cannot be only based on its location. Mutations on Nav1.4 DI, DII, and DIII VSDs have been shown to cause the creation of a gating pore. Against all expectations, similar mutations on DIV VSD do not create a gating pore. This difference appears to be related to a larger hydrophobic septum in this particular VSD (Gosselin-Badaroudine et al., 2012a). Oppositely, the study of the similar DIV VSD of Cav1.1 channels has clearly shown that mutations in this domain create gating pores (Fan et al., 2013). Nevertheless, given the high sequence homology between Nav1.4 and Nav1.5, the creation of gating pores due to mutations in Nav1.5 DI, DII, or DIII VSD’s is expected. While single mutations on Nav1.5 DI, Nav1.4 DI, and Nav1.4 DII VSDs have been shown to open a gating pore (Gosselin-Badaroudine et al., 2012a,b; Moreau et al., 2015), other studies revealed that multiple mutations are required to create a gating pore in Nav1.2 DII or in the Shaker VSD (Sokolov et al., 2005; Gamal El-Din et al., 2010). All together, these results clearly demonstrate that caution is warranted to conclude that a mutation creates or not a gating pore based on its location in the VSD.
We report that gating pores are caused by two other recently characterized mutations, R225P and R814W, which are located in DI and DII, respectively, of Nav1.5 (Figure 1). We previously studied the R222Q and R225W Nav1.5 DI VSD mutations (Moreau et al., 2015). Besides the creation of gating pores due to both mutations, several biophysical defects were also recorded. This study notably reports an increase (R222Q) or a decrease (R225W) in current density and hyperpolarized (R222Q) or depolarized (R225W) shifts in the steady state activation. Thus, biophysical defects of R222Q pointed toward a gain of function while defects of R225W mutants indicated a loss of channel function (Moreau et al., 2015). Unlike the R225W mutation, the R225P and R814W mutations both cause a gain of channel function, although the defect is less pronounced with the R814W mutant. Gating pores created by these mutations are activated by membrane depolarizations. Interestingly, the biophysical defects observed for alpha pore currents should also affect gating pore currents. The left shift of steady state activation should indicate an early movement of the VSD, resulting in an early opening of the gating pore. Nevertheless, the voltage dependence of the activation is a complicated process resulting from both the movement of the VSD and the coupling between the VSD and the pore domain. The correlation between the alpha pore voltage dependence of activation and the movement of a VSD is thus not necessarily direct. Here, both steady state activation curves were left shifted (Table 1). This could be due to an earlier movement of the VSD or a tighter coupling between the VSD and the PD. In the case of an earlier movement of the VSD, this would result in an earlier opening of the gating pore. In our previous study of R222Q and R225W, the voltage dependence of R222Q gating pore was slightly left shifted when compared to the voltage dependence of the R225W gating pore (Moreau et al., 2015). This result would indicate that shifts in the voltage dependence of alpha pore currents would be due to differences in the mutants’ VSDs movement. Indeed, the voltage dependence of the alpha pore current of the R222Q mutant was importantly hyperpolarized when compared with the voltage dependence of the alpha pore current of the R225W mutant (Moreau et al., 2015).
The gating pores currents generated by the R225P and R814W mutations were recorded using Cs+ due to its higher permeability. However, based on previously reported selectivity sequences for gating pores, K+ and Na+ are highly expected to permeate this pathological pore and should affect electrical signals and ionic homeostasis (Moreau et al., 2014a,b, 2015). Like other depolarization-activated gating pores (Sokolov et al., 2008; Fan et al., 2013; Groome et al., 2014; Moreau et al., 2015), the gating pores created by the R225P and R814W mutations remained open after long periods of depolarization, probably due to the temporary partial freezing of the S4 segment (Gamal El-Din et al., 2014; Moreau et al., 2014b). Given the physiological ion concentrations on both sides of the membrane and the biophysical properties of these gating pores, the pores should cause outward K+ and Na+ leaks under depolarized conditions (during the action potential plateau phase) and inward Na+ leak at the end of each action potential (Moreau et al., 2014a,b, 2015). Gating pore currents constitute small leak currents opened as soon as the voltage sensor has moved (Moreau et al., 2014a,b, 2015). Considering their small amplitude, their contribution to the pathologic mechanism could thus be questioned. Nevertheless, it is now well-accepted that small late Nav1.5 currents can cause the LQT3 syndrome (Ruan et al., 2009; Amin et al., 2010). Although gating pore currents or late currents are of small amplitude, they are opened for long durations. Consequently, their impact would probably be related to the continuous stress imposed to cardiac myocytes. Several pathological consequences of gating pores have been described in the case of Nav1.4 mutations associated with the development of hypokalemic periodic paralysis (HypoPP). In this case, an important elevation of the resting membrane potential (VRest) is thought to be the cause of clinical manifestations (Rudel et al., 1984; Jurkat-Rott et al., 2009). Cation leaks due to gating pores usually represent 1–6% of the alpha pore current amplitude (Moreau et al., 2014a,b). Consequently, charges flowing through the gating pores are thought to not be enough to cause this elevation in VRest. This VRest elevation would most probably be related to the ionic homeostasis imbalance (intracellular Na+ and Ca2+ overload) caused by the gating pore. The ionic homeostasis imbalance would block ion channels such as Kir channels and thus cause the VRest elevation (Hughes and Swaminathan, 2008; Struyk and Cannon, 2008; Tricarico and Camerino, 2011; Moreau et al., 2014a). A similar process is thus expected to affect patient’s cardiac myocytes, resulting in a potentially pathological depolarization of their VRest. This ionic homeostasis imbalance could then result in several other pathological consequences such as the block of connexins, impair the excitation-contraction coupling or also the function and the assembly of contractile proteins (Fabiato and Fabiato, 1978; Chandra et al., 1999; Stergiopoulos et al., 1999; Bukauskas et al., 2001; Duffy et al., 2004; Feng et al., 2008; Moreau et al., 2014a). Taken together, these defects could result in the atypical clinical phenotype associating DCM and mixed arrhythmias. Furthermore, the pathological potential of such gating pore currents is strongly strengthened by the association of a similar mutation in Nav1.4 (R1135C/H, R3 – S4, DIII) associated with HypoPP phenotype (Groome et al., 2014).
The R225P and R814W mutations should disrupt interactions between the S4 segment and the GCTC. Based on previously published results (Moreau et al., 2014a,b, 2015), these two mutations should also affect the structure of the VSD, allowing the merging of water crevices as illustrated in Figure 6. Once the water crevices merge, a new permeation pathway appears through the usually non-conductive VSD structure, allowing ions to permeate. The substitution of the arginine by a tryptophan (R814W) causes a major impediment to the S4 structure and its movement, resulting in a twist of the segment due to the affinity of tryptophan for lipid tails. The substitution of the arginine by a proline (R225P) should also cause a major modification of the S4 structure. While other reported mutations mainly result in modifications to the side chain of the amino acid, the substitution by a proline affects both the side chain and the alpha carbon. Due to these modifications and based on previously published results, the creation of a continuous water wire through the VSD would thus underlie the creation of gating pores in DI and DII of Nav1.5 channels. Indeed, although extensive perturbations to the VSD structure such as the ones possibly induced by the proline and the tryptophan may be hard to resolve in the 100 ns time, our models show that the mutant VSDs are potentially permeable to ions. The water density profile was assessed by counting the average occupancy of water molecules in the VSD along the z-axis (perpendicular to the membrane) over 6 ns (Figure 6B). For both mutants the hydrophobic septa (volumes that are not accessible to water) are drastically shorter and more hydrated than the corresponding WT VSDs. In contrast with our previous study (Moreau et al., 2015), along very few angstroms, the water accessible volume in the R225P and R814W VSDs is below one water molecule but clearly above 0. These results more likely indicate that in some frames, no water molecules can be found, while in other frames at the same location a water molecule is found. Since the water density profile is built on at least 600 frames, the mean value is neither 0 nor above 1. This behavior was not observed in WT VSDs as they featured constriction sites in which the water profile value was 0 over several angstroms. Finally, consistently with both our previous study and our experimental data, these results would indicate that R225P and R814W mutated VSDs can conduct ions due to hydrophobic septum disruption resulting in water crevices junction.
FIGURE 6. Molecular basis for the creation of gating pores. (A) Models of the DI and DII of Nav1.5 were built based on previously published results (Moreau et al., 2015). The top row shows the DI models. The bottom row shows the DII models. The left panels represent WT VSDs and the right panels represent mutant VSDs (top R225P, bottom R814W). The protein is shown as orange ribbons, and the water crevice is shown in blue. Positively charged residues in the S4 segment and residues making up the GCTC are highlighted. Positively charged residues are depicted in blue, negatively charged residues in red, aromatic residues in yellow, and the proline in gray. The complete VSD is shown (left) for each panel. A higher magnification of the GCTC is also shown to clearly see the interactions between the S4 segment and the GCTC (middle) as well as the merging of the water crevices from both sides of the membrane (right) for the mutant channels. For purposes of clarity, the S1 segment has been removed. (B) Water density profiles of the Nav1.5/DI (orange) – DII (green) WT and R225P – R814W mutant channels (dash lines). The histograms were built using a 1-A grid, and the averages were calculated from at least 6 ns of the trajectories. 0 corresponds to the position of the Cα of Y168 of S2 (DI) or the position of the Cα of F760 of S2 (DII).
Further investigation is also required concerning VSD mutations located outside the S4 segment. Indeed, mutations in S1–S3 segments are known to be associated with the development of arrhythmias and dilated cardiomyopathy (McNair et al., 2004, 2011; Olson et al., 2005; Nguyen et al., 2008; Moreau et al., 2014a; Swan et al., 2014). The creation of a gating pore due to such mutations has never been investigated. Nevertheless, mutations in the VSD but outside the S4 segment might also disrupt interactions between the S4 segment and the GCTC and thus cause the creation of a gating pore.
Gating pores are now well-accepted to be also involved in the pathological development of other muscular or neuronal disorders (periodic paralysis and peripheral nerve hyper-excitability; Sokolov et al., 2007, 2008; Struyk and Cannon, 2007; Miceli et al., 2012; Moreau et al., 2014a). The exact pathological downstream consequences of these cation leaks in myocytes require further exploration. Patient specific induced pluripotent stem cells (hiPSC) or transgenic mouse models could thus constitute models of choice to further investigate the pathologic potential of these leak currents and clearly study their pathologic importance. Nevertheless, we hypothesize that such leaks unbalance ionic homeostasis, thus impairing electrical signals, cellular communication, and excitation-contraction coupling (Moreau et al., 2014a, 2015). The combination of gating pores and defects of the alpha pore could thus lead to the atypical clinical phenotype observed.
Our study identified gating pores induced by two unrelated Nav1.5 mutations. This is the first time that a mutation generating a gating pore outside the DI of Nav1.5 has been investigated. The R225P and R814W mutations both caused an atypical clinical phenotype associating cardiac arrhythmia with dilatation, bringing the number of known mutations that create gating pores related to this clinical phenotype to five. Given that the alpha pores of several mutated channels are also defective, gating pores likely contribute to the pathogenesis rather than cause it. However, the presence of these gating pores may determine whether the mutations result in a pure electrical disorder or a pathology associating electrical and morphological abnormalities. We thus propose that gating pores that are created by Nav1.5 mutations on the VSD and that are associated with combined electrical and morphological cardiac pathologies should be investigated as a biophysical defect. Finally, while this biophysical characterization is a prerequisite to the better understanding of the pathogenesis, future studies of transgenic animals are warranted for the comprehension of the link between gating pore current and the clinical phenotype combining complex arrhythmias and DCM.
AM performed all the experiments involving the biophysical characterization of the wild-type and mutant channels in mammalian cells, he was involved in the recording of gating pore currents, analyzed the data, and helped drafted the manuscript. PG-B performed dynamic simulations and helped drafted the manuscript. MB discuss the data and helped drafted the manuscript MC developed the concept, designed the project and experiments, analyzed and interpreted the results, and drafted the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was supported by grants from the Heart and Stroke Foundation of Quebec and the Canadian Institutes of Health Research (MOP-111072 and MOP-130373). PG-B received a Postgraduate Scholarship from NSERC.
Beckermann, T. M., McLeod, K., Murday, V., Potet, F., and George, A. L. Jr. (2014). Novel SCN5A mutation in amiodarone-responsive multifocal ventricular ectopy-associated cardiomyopathy. Heart Rhythm. 11, 1446–1453. doi: 10.1016/j.hrthm.2014.04.042
Bezzina, C. R., Rook, M. B., Groenewegen, W. A., Herfst, L. J., Van Der Wal, A. C., Lam, J., et al. (2003). Compound heterozygosity for mutations (W156X and R225W) in SCN5A associated with severe cardiac conduction disturbances and degenerative changes in the conduction system. Circ. Res. 92, 159–168. doi: 10.1161/01.RES.0000052672.97759.36
Bukauskas, F. F., Bukauskiene, A., Bennett, M. V., and Verselis, V. K. (2001). Gating properties of gap junction channels assembled from connexin43 and connexin43 fused with green fluorescent protein. Biophys. J. 81, 137–152. doi: 10.1016/S0006-3495(01)75687-1
Chandra, M., Montgomery, D. E., Kim, J. J., and Solaro, R. J. (1999). The N-terminal region of troponin T is essential for the maximal activation of rat cardiac myofilaments. J. Mol. Cell Cardiol. 31, 867–880. doi: 10.1006/jmcc.1999.0928
Delemotte, L., Tarek, M., Klein, M. L., Amaral, C., and Treptow, W. (2011). Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations. Proc. Natl. Acad. Sci. U.S.A. 108, 6109–6114. doi: 10.1073/pnas.1102724108
Deschenes, I., Baroudi, G., Berthet, M., Barde, I., Chalvidan, T., Denjoy, I., et al. (2000). Electrophysiological characterization of SCN5A mutations causing long QT (E1784K) and Brugada (R1512W and R1432G) syndromes. Cardiovasc. Res. 46, 55–65. doi: 10.1016/S0008-6363(00)00006-7
Duffy, H. S., Ashton, A. W., O’donnell, P., Coombs, W., Taffet, S. M., Delmar, M., et al. (2004). Regulation of connexin43 protein complexes by intracellular acidification. Circ. Res. 94, 215–222. doi: 10.1161/01.RES.0000113924.06926.11
Eswar, N., Webb, B., Marti-Renom, M. A., Madhusudhan, M. S., Eramian, D., Shen, M. Y., et al. (2006). Comparative protein structure modeling using Modeller. Curr. Protoc. Bioinformatics Chap. 5, Unit 5.6. doi: 10.1002/0471250953.bi0506s15
Fabiato, A., and Fabiato, F. (1978). Effects of pH on the myofilaments and the sarcoplasmic reticulum of skinned cells from cardiace and skeletal muscles. J. Physiol. 276, 233–255. doi: 10.1113/jphysiol.1978.sp012231
Fan, C., Lehmann-Horn, F., Weber, M. A., Bednarz, M., Groome, J. R., Jonsson, M. K., et al. (2013). Transient compartment-like syndrome and normokalaemic periodic paralysis due to a Ca(v)1.1 mutation. Brain 136, 3775–3786. doi: 10.1093/brain/awt300
Feng, H. Z., Biesiadecki, B. J., Yu, Z. B., Hossain, M. M., and Jin, J. P. (2008). Restricted N-terminal truncation of cardiac troponin T: a novel mechanism for functional adaptation to energetic crisis. J. Physiol. 586, 3537–3550. doi: 10.1113/jphysiol.2008.153577
Gamal El-Din, T. M., Heldstab, H., Lehmann, C., and Greeff, N. G. (2010). Double gaps along Shaker S4 demonstrate omega currents at three different closed states. Channels (Austin) 4, 93–100. doi: 10.4161/chan.4.2.10672
Gamal El-Din, T. M., Scheuer, T., and Catterall, W. A. (2014). Tracking S4 movement by gating pore currents in the bacterial sodium channel NaChBac. J. Gen. Physiol. 144, 147–157. doi: 10.1085/jgp.201411210
Gellens, M. E., George, A. L. Jr., Chen, L. Q., Chahine, M., Horn, R., Barchi, R. L., et al. (1992). Primary structure and functional expression of the human cardiac tetrodotoxin-insensitive voltage-dependent sodium channel. Proc. Natl. Acad. Sci. U.S.A. 89, 554–558. doi: 10.1073/pnas.89.2.554
Gosselin-Badaroudine, P., Delemotte, L., Moreau, A., Klein, M. L., and Chahine, M. (2012a). Gating pore currents and the resting state of Nav1.4 voltage sensor domains. Proc. Natl. Acad. Sci. U.S.A. 109, 19250–19255. doi: 10.1073/pnas.1217990109
Gosselin-Badaroudine, P., Keller, D. I., Huang, H., Pouliot, V., Chatelier, A., Osswald, S., et al. (2012b). A proton leak current through the cardiac sodium channel is linked to mixed arrhythmia and the dilated cardiomyopathy phenotype. PLoS ONE 7:e38331. doi: 10.1371/journal.pone.0038331
Gosselin-Badaroudine, P., Moreau, A., and Chahine, M. (2014). Nav 1.5 mutations linked to dilated cardiomyopathy phenotypes: is the gating pore current the missing link? Channels (Austin) 8, 90–94. doi: 10.4161/chan.27179
Groome, J. R., Lehmann-Horn, F., Fan, C., Wolf, M., Winston, V., Merlini, L., et al. (2014). Nav1.4 mutations cause hypokalaemic periodic paralysis by disrupting IIIS4 movement during recovery. Brain 137, 998–1008. doi: 10.1093/brain/awu015
Henrion, U., Renhorn, J., Borjesson, S. I., Nelson, E. M., Schwaiger, C. S., Bjelkmar, P., et al. (2012). Tracking a complete voltage-sensor cycle with metal-ion bridges. Proc. Natl. Acad. Sci. U.S.A. 109, 8552–8557. doi: 10.1073/pnas.1116938109
Huang, H., Priori, S. G., Napolitano, C., O’leary, M. E., and Chahine, M. (2011). Y1767C, a novel SCN5A mutation, induces a persistent Na+ current and potentiates ranolazine inhibition of Nav1.5 channels. Am. J. Physiol. Heart. Circ. Physiol. 300, H288–H299. doi: 10.1152/ajpheart.00539.2010
Hughes, B. A., and Swaminathan, A. (2008). Modulation of the Kir7.1 potassium channel by extracellular and intracellular pH. Am. J. Physiol. Cell Physiol. 294, C423–C431. doi: 10.1152/ajpcell.00393.2007
Jurkat-Rott, K., Weber, M. A., Fauler, M., Guo, X. H., Holzherr, B. D., Paczulla, A., et al. (2009). K+-dependent paradoxical membrane depolarization and Na+ overload, major and reversible contributors to weakness by ion channel leaks. Proc. Natl. Acad. Sci. U.S.A. 106, 4036–4041. doi: 10.1073/pnas.0811277106
Laurent, G., Saal, S., Amarouch, M. Y., Beziau, D. M., Marsman, R. F., Faivre, L., et al. (2012). Multifocal ectopic Purkinje-related premature contractions: a new SCN5A-related cardiac channelopathy. J. Am. Coll. Cardiol. 60, 144–156. doi: 10.1016/j.jacc.2012.02.052
Mann, S. A., Castro, M. L., Ohanian, M., Guo, G., Zodgekar, P., Sheu, A., et al. (2012). R222Q SCN5A mutation is associated with reversible ventricular ectopy and dilated cardiomyopathy. J. Am. Coll. Cardiol. 60, 1566–1573. doi: 10.1016/j.jacc.2012.05.050
McNair, W. P., Ku, L., Taylor, M. R., Fain, P. R., Dao, D., Wolfel, E., et al. (2004). SCN5A mutation associated with dilated cardiomyopathy, conduction disorder, and arrhythmia. Circulation 110, 2163–2167. doi: 10.1161/01.CIR.0000144458.58660.BB
McNair, W. P., Sinagra, G., Taylor, M. R., Di Lenarda, A., Ferguson, D. A., Salcedo, E. E., et al. (2011). SCN5A mutations associate with arrhythmic dilated cardiomyopathy and commonly localize to the voltage-sensing mechanism. J. Am. Coll. Cardiol. 57, 2160–2168. doi: 10.1016/j.jacc.2010.09.084
Miceli, F., Vargas, E., Bezanilla, F., and Taglialatela, M. (2012). Gating currents from Kv7 channels carrying neuronal hyperexcitability mutations in the voltage-sensing domain. Biophys. J. 102, 1372–1382. doi: 10.1016/j.bpj.2012.02.004
Moreau, A., Gosselin-Badaroudine, P., Delemotte, L., Klein, M. L., and Chahine, M. (2015). Gating pore currents are defects in common with two Nav1.5 mutations in patients with mixed arrhythmias and dilated cardiomyopathy. J. Gen. Physiol. 145, 93–106. doi: 10.1085/jgp.201411304
Moreau, A., Keller, D. I., Huang, H., Fressart, V., Schmied, C., Timour, Q., et al. (2012). Mexiletine differentially restores the trafficking defects caused by two brugada syndrome mutations. Front. Pharmacol. 3:62. doi: 10.3389/fphar.2012.00062
Moreau, A., Krahn, A. D., Gosselin-Badaroudine, P., Klein, G. J., Christe, G., Vincent, Y., et al. (2013). Sodium overload due to a persistent current that attenuates the arrhythmogenic potential of a novel LQT3 mutation. Front. Pharmacol. 4:126. doi: 10.3389/fphar.2013.00126
Nair, K., Pekhletski, R., Harris, L., Care, M., Morel, C., Farid, T., et al. (2012). Escape capture bigeminy: phenotypic marker of cardiac sodium channel voltage sensor mutation R222Q. Heart. Rhythm. 9, 1681–1688. doi: 10.1016/j.hrthm.2012.06.029
Nguyen, T. P., Wang, D. W., Rhodes, T. H., and George, A. L. Jr. (2008). Divergent biophysical defects caused by mutant sodium channels in dilated cardiomyopathy with arrhythmia. Circ. Res. 102, 364–371. doi: 10.1161/CIRCRESAHA.107.164673
Olson, T. M., Michels, V. V., Ballew, J. D., Reyna, S. P., Karst, M. L., Herron, K. J., et al. (2005). Sodium channel mutations and susceptibility to heart failure and atrial fibrillation. JAMA 293, 447–454. doi: 10.1001/jama.293.4.447
Rudel, R., Lehmann-Horn, F., Ricker, K., and Kuther, G. (1984). Hypokalemic periodic paralysis: in vitro investigation of muscle fiber membrane parameters. Muscle Nerve 7, 110–120. doi: 10.1002/mus.880070205
Sokolov, S., Scheuer, T., and Catterall, W. A. (2005). Ion permeation through a voltage- sensitive gating pore in brain sodium channels having voltage sensor mutations. Neuron 47, 183–189. doi: 10.1016/j.neuron.2005.06.012
Sokolov, S., Scheuer, T., and Catterall, W. A. (2008). Depolarization-activated gating pore current conducted by mutant sodium channels in potassium-sensitive normokalemic periodic paralysis. Proc. Natl. Acad. Sci. U.S.A. 105, 19980–19985. doi: 10.1073/pnas.0810562105
Stergiopoulos, K., Alvarado, J. L., Mastroianni, M., Ek-Vitorin, J. F., Taffet, S. M., and Delmar, M. (1999). Hetero-domain interactions as a mechanism for the regulation of connexin channels. Circ. Res. 84, 1144–1155. doi: 10.1161/01.RES.84.10.1144
Struyk, A. F., and Cannon, S. C. (2007). A Na+ channel mutation linked to hypokalemic periodic paralysis exposes a proton-selective gating pore. J. Gen. Physiol. 130, 11–20. doi: 10.1085/jgp.200709755
Struyk, A. F., and Cannon, S. C. (2008). Paradoxical depolarization of BA2+- treated muscle exposed to low extracellular K+: insights into resting potential abnormalities in hypokalemic paralysis. Muscle Nerve 37, 326–337. doi: 10.1002/mus.20928
Swan, H., Amarouch, M. Y., Leinonen, J., Marjamaa, A., Kucera, J. P., Laitinen-Forsblom, P. J., et al. (2014). Gain-of-function mutation of the SCN5A gene causes exercise-induced polymorphic ventricular arrhythmias. Circ. Cardiovasc. Genet. 7, 771–781. doi: 10.1161/CIRCGENETICS.114.000703
Tombola, F., Pathak, M. M., and Isacoff, E. Y. (2005). Voltage-sensing arginines in a potassium channel permeate and occlude cation-selective pores. Neuron 45, 379–388. doi: 10.1016/j.neuron.2004.12.047
Keywords: SCN5A, gating pore current, arrhythmias, dilated cardiomyopathy, voltage sensor
Citation: Moreau A, Gosselin-Badaroudine P, Boutjdir M and Chahine M (2015) Mutations in the Voltage Sensors of Domains I and II of Nav1.5 that are Associated with Arrhythmias and Dilated Cardiomyopathy Generate Gating Pore Currents. Front. Pharmacol. 6:301. doi: 10.3389/fphar.2015.00301
Received: 20 October 2015; Accepted: 09 December 2015;
Published: 24 December 2015.
Edited by:Hugues Abriel, University of Bern, Switzerland
Reviewed by:Stephan Kellenberger, University of Lausanne, Switzerland
Peter Ruben, Simon Fraser University, Canada
Copyright © 2015 Moreau, Gosselin-Badaroudine, Boutjdir and Chahine. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mohamed Chahine, firstname.lastname@example.org