Abstract
Mexiletine is a class Ib antiarrhythmic drug and is also used clinically to reduce or prevent myotonia. In addition, mexiletine has neuroprotective effects in models of brain ischemia. We compared state-dependent affinities of mexiletine for Nav1.2, Nav1.4, and Nav1.5 and examined the effects of mutations of transmembrane segment S6 residues on mexiletine block of Nav1.5. Three channel isoforms had similar affinities of mexiletine for the rested state, and Nav1.4 and Nav1.5 had similar affinities for the open and inactivated states, while Nav1.2 had lower affinity for these states than Nav1.4 and Nav1.5. Mutational studies revealed that the largest affinity change was observed for an Ala substitution of Phe in domain IV S6. In our homology modeling based on the bacterial Na+ channel, mexiletine changed its location and orientation in the pore depending on the state of the channel, irrespective of the channel isoform. Mexiletine occurred in the upper part in the pore in the open state and lower in the closed state. High-affinity binding of mexiletine in the open states of Nav1.4 and Nav1.5 was caused by a π–π interaction with Phe, whereas mexiletine was located away from Phe in the open state of Nav1.2. These results provide crucial information on the mechanism of isoform differences in state-dependent block by local anesthetics and related drugs. Mexiletine at upper locations in the open state may effectively cause an electrostatic mechanism of block.
Introduction
Mexiletine, a class Ib antiarrhythmic drug, is considered the first line for treating myotonic syndromes (Jackson et al., 1994; Hudson et al., 1995; Camerino et al., 2007). In addition, mexiletine is clinically useful to relieve neuropathic pain (Chabal et al., 1992) and protects brain neurons from injury in the models of ischemia (Stys and Lesiuk, 1996).
Mexiletine is an analog of the local anesthetic lidocaine and exerts its therapeutic action through blockade of voltage-gated Na+ channels. Although mexiletine blocks various channel isoforms, there have been few reports on isoform differences in affinities of mexiletine. Affinity is known to vary as a function of the sodium channel gating state. It has been previously reported that the affinity of mexiletine for the inactivated channel is 2-fold higher in Nav1.5 than in Nav1.2 (Weiser et al., 1999) and Nav1.5 is more sensitive to use-dependent block by mexiletine than Nav1.7 (Wang et al., 2015). On the other hand, mexiletine caused more blocking for Nav1.5 than Nav1.2 or Nav1.4 by increasing the proportion of inactivated channels in Nav1.5 (Kawagoe et al., 2002). Although the isoform differences in local anesthetic affinity seem to be less than in toxin affinity (Fozzard and Hanck, 1996), understanding the molecular mechanisms of isoform differences in affinities may contribute to rational design of isoform-selective and safer Na+ channel blockers.
Elucidating the mechanisms of state-dependent block by local anesthetics is obviously important. Phe in domain IV S6 has long been known as a critical determinant of local anesthetic block (Ragsdale et al., 1994, 1996). Other residues in S6 segments of domains I, III, and IV also affect the local anesthetic binding (Wang et al., 1998, 2000; Nau et al., 1999; Yarov-Yarovoy et al., 2001, 2002). Although mutations of these residues are well known to have different effects on the block of the closed, open, and inactivated channels (Ragsdale et al., 1994, 1996; Wang et al., 1998, 2000; Nau et al., 1999; Yarov-Yarovoy et al., 2001, 2002), the molecular mechanisms of these differences have not been clarified. Several structural models of Na+ channels have analyzed the interactions of local anesthetics with the pore residues in the closed (Tikhonov et al., 2006; Bruhova et al., 2008; Tikhonov and Zhorov, 2012) and open states (Lipkind and Fozzard, 2005; Tikhonov et al., 2006; Tikhonov and Zhorov, 2007, 2017). Lidocaine molecules adopted different locations and orientations in the closed and open states of the Nav1.5 model (Tikhonov et al., 2006), but most of their models are for either the closed or open state, and for either channel isoform of Nav1.4 or Nav1.5.
To our knowledge, there is no comparison of local anesthetics docking into two or more different channel isoforms. Previously, homology models were based on the bacterial K+ channel templates instead of the Na+ channel templates (Lipkind and Fozzard, 2005; Tikhonov et al., 2006; Tikhonov and Zhorov, 2007; Bruhova et al., 2008). Nevertheless, sequences of mammalian Na+ channels are more similar to the bacterial Na+ channels than to the bacterial K+ channels, and the Na+ channel alignment is much less ambiguous than that of K+ channels (Tikhonov and Zhorov, 2017). Actually, a novel distinct binding mode of lidocaine was found in the Na+ channel model based on the Na+ channel template (Tikhonov and Zhorov, 2017), but this mode was not revealed in the models based on the K+ channel templates (Tikhonov et al., 2006; Tikhonov and Zhorov, 2007).
In the present study, we addressed isoform differences in state-dependent block of mexiletine. First, we measured state-dependent affinities of mexiletine for Nav1.2, Nav1.4, and Nav1.5. The affinities of mexiletine for the rested state were not different among these isoforms. Nav1.4 and Nav1.5 had similar affinities for the open and inactivated states, while Nav1.2 had a lower affinity for these states than Nav1.4 and Nav1.5. We next visualized the interactions of mexiletine with the pore residues using bacterial Na+ channel templates. Mexiletine was located at a different place in the closed and open states irrespective of the channel isoform, and mexiletine was found closer to the extracellular space in the open state. High-affinity binding of mexiletine in the open states of Nav1.4 and Nav1.5 was caused by a π–π interaction with Phe, which was consistent with our mutational studies of Nav1.5. In contrast, mexiletine was not close to Phe in the open state of Nav1.2.
Materials and Methods
Site-Directed Mutagenesis and Transfection Procedure.
The human heart Na+ channel α-subunit (Nav1.5) clone in the pcDNA3.1+ plasmid (Invitrogen, Carlsbad, CA) was a generous gift from Dr. N. Makita (Nagasaki University). All mutations (N406A, L1462A, L1462F, L1462E, F1760A) of Nav1.5 were introduced using the QuikChange Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA), according to the manufacturer’s protocols and as described in our previous report (Sunami et al., 2004). Mutations were confirmed by DNA sequencing of the full-length plasmid. Human embryonic kidney cells (HEK293 cells) were transiently transfected with the plasmid encoding wild-type (WT) or mutant Nav1.5 channels using SuperFect (Qiagen, Hilden, Germany) in combination with a plasmid encoding CD8 (pIRES-CD8) to identify transfected cells with Dynabeads M-450 CD8 (Dynal, Oslo, Norway). Electrophysiologic measurements were performed 24–72 hours after transfection.
WT channels of the rat brain IIA (Nav1.2) and the rat skeletal muscle Na+ channel α-subunit (Nav1.4) were stably expressed in HEK293 cells. The pRc/CMV plasmid (Invitrogen) encoding Nav1.2 and the pCI-neo plasmid (Promega, Madison, WI) encoding Nav1.4 were generous gifts from Dr. H. A. Fozzard (University of Chicago). Selection of stablely transfected cells was made in Dulbecco’s modified Eagle’s medium containing 600 µg/ml G418 (geneticin). The Na+ channel β-subunits were not used in all experiments.
Electrophysiology and Data Analyses.
Na+ currents were recorded with the whole-cell configuration of patch clamp using an Axopatch 200B amplifier (Axon Instruments, Union City, CA). Currents were filtered at 5 kHz through a four-pole Bessel filter and digitized at 100 kHz at 16-bit resolution with a Digidata 1322A acquisition system using pCLAMP 9 software (Axon Instruments). The pipettes had resistances between 1.0 and 1.5 MΩ. Recordings were made at room temperature (20–22°C) in a bathing solution containing (millimolar): 140 NaCl, 2.5 KCl, 2 CaCl2, 1 MgCl2, 10 HEPES (pH 7.4 with NaOH). The pipette solution was (millimolar): 100 CsF, 40 NaCl, 10 HEPES (pH 7.2 with CsOH).
To determine the activation parameters, the current-voltage (I-V) relationship was fitted to a transform of a Boltzmann distribution: I = (V − Vrev)Gmax/{1 + exp[(V1/2 − V)/s]}, where I is the peak Na+ current during the test pulse of voltage V. The parameters estimated by the fitting were V1/2 (the voltage for half-activation), s (slope factor), Vrev (the reversal potential), and Gmax (the maximum peak conductance). After induction of steady-state inactivation by 1-second depolarizing prepulses from a holding potential of −130 mV, the Na+ currents were measured during test pulses to −20 or 0 mV applied every 4 seconds. Availability was described by the following Boltzmann equation: I/Imax = 1/{1 + exp[(V − V1/2)/s]}, where I is the peak current, Imax is the maximum peak current, V is the prepulse voltage, V1/2 is the voltage for half-inactivation, and s is the slope factor. The recovery from inactivation induced by 10-millisecond depolarization to −10 mV was monitored, and the time course of recovery was fitted by a single exponential. The time to decay from 90% to 10% of the peak current at 0 mV was also measured.
For mexiletine block experiments, a train of 30 pulses of 5- or 200-millisecond duration was applied at 300-millisecond interpulse intervals in the absence or presence of mexiletine. To compare the sensitivity of three channel isoforms to mexiletine block, we used holding potentials of −120 mV for Nav1.2 and Nav1.4, and −130 mV for Nav1.5, and test potentials of 0 mV for Nav1.2 and Nav1.4, and −20 mV for Nav1.5. These voltages were chosen based on the gating data (Table 1), and the voltages minimized the gating effects on mexiletine block. In Nav1.5 mutant channels, tonic (or rested) and use-dependent block by mexiletine were produced by a train of 30 pulses of 5- or 200-millisecond duration applied to 0 (N406A, L1462E), −5 (L1462A), or −20 mV (L1462F, F1760A) from a holding potential of −130 (N406A, L1462A, F1760A) or −140 mV (L1462E, L1462F). The dissociation constants of mexiletine for the rested (KR), open (KO), and inactivated (KI) states of the channels were determined by the concentration–response curves for tonic block, and for total block (tonic block + use-dependent block) caused by 5- and 200-millisecond depolarizing pulses, respectively. Mexiletine was obtained from Sigma-Aldrich (St. Louis, MO).
Gating parameters of Nav1.2, Nav1.4, Nav1.5, and Nav1.5 mutants
The voltage dependence of activation and inactivation, and recovery from fast inactivation were measured by the method described in Materials and Methods. Data represent the mean ± S.D.
Data were analyzed using pCLAMP 9 (Axon Instruments) and SigmaPlot (SPSS, Chicago, IL). Pooled data are presented as the mean ± S.D. Statistical comparisons were made by using one- or two-way analysis of variance.
Molecular Modeling and Docking.
As with our previous study (Sunami et al., 2004), we built homology models of the closed and open pores of Nav1.2, Nav1.4, and Nav1.5 based on the crystal structures of bacterial voltage-gated Na+ channels, NavAb (for closed channel; PDB accession code 3RVY; Payandeh et al., 2011) and NavMs (for open channel; PDB accession code 5HVX; Sula et al., 2017) using MODELER version 9.17 (Šali and Blundell, 1993). Sequence alignments of the pore regions of NavAb, NavMs, Nav1.2, Nav1.4, and Nav1.5 are given in Fig. 1. The geometry optimization and partial charge calculation of mexiletine was performed using the PM6 method of MOPAC2009 (http://OpenMOPAC.net). Docking mexiletine in these models was performed using the lamarckian genetic algorithm of AutoDock4.2.6 and AutoDockTools (Morris et al., 2009). Side-chain flexibility was allowed during docking for the following S6 residues: S401, V405, I408, N927, L931, F934, S1458, T1461, L1462, F1465, F1760, V1764, and Y1767 in Nav1.5, and the equivalent residues in Nav1.2 and Nav1.4. The docking grid was composed of 28 × 28 × 31 points with a point spacing of 1 Å. Autodock was run using the default parameter set except for “Number of GA Runs” set to 350. The docking results were visualized using Biovia Discovery Studio Visualizer (Dassault Systemes Biovia, 2016).
Sequence alignment of the pore regions of NavAb, NavMs, Nav1.2, Nav1.4, and Nav1.5. The proposed identities of the P1 and P2 helices, and the inner helix or the transmembrane segment S6 are indicated above the sequences. The number next to each inner helix or S6 sequence indicates the amino acid number of the end of each inner helix or S6. The amino acid residues indicated in bold are conserved among Nav1.2, Nav1.4, and Nav1.5. The residues in red (N406, L1462, F1760) were mutated in the present study, and these are conserved among these three channels.
Results
Isoform Differences in Mexiletine Block.
We first compared the sensitivity of Nav1.2, Nav1.4, and Nav1.5 to tonic block (first-pulse block or rested block) by mexiletine, which was measured at the holding potential of −120 mV for Nav1.2 and Nav1.4, and −130 mV for Nav1.5. In Nav1.5, tonic block by 40 µM mexiletine was 9.5% ± 5.8% (n = 5) and was not significantly different from the tonic block of Nav1.2 (4.5% ± 9.6%, n = 5) or Nav1.4 (8.0% ± 2.5%, n = 5) (Fig. 2, A and B).
Tonic and use-dependent block by mexiletine of Nav1.2, Nav1.4, and Nav1.5. (A) Typical current traces before (Con) and after application of 40 µM mexiletine (Mex) recorded from Nav1.2 (left), Nav1.4 (middle), and Nav1.5 (right). To produce tonic block (TB) and use-dependent block (UDB) by mexiletine, a train of 30 pulses of 200-millisecond duration was applied at 300-millisecond interpulse intervals in the absence and presence of mexiletine. Current traces are during the first pulse for the control, and the 1st and 30th pulses for mexiletine. See Materials and Methods for the voltages of the pulse protocol. (B) Tonic block of Nav1.2, Nav1.4, and Nav1.5 by 40 µM mexiletine. (C) Use-dependent block of Nav1.2, Nav1.4, and Nav1.5 by 40 µM mexiletine. Use-dependent block by mexiletine was produced by 30 pulses with (left) 5- or (right) 200-millisecond pulse duration (PD) at 300-millisecond interpulse intervals. In B and C, data represent the mean ± S.D. from five cells. *P < 0.05; ***P < 0.001 compared with Nav1.5; +++P < 0.001 compared with Nav1.4, two-way analysis of variance followed by Bonferroni test.
We next compared use-dependent block by mexiletine in three channel isoforms. In the control, a train of 30 pulses of 5- or 200-millisecond pulse duration produced no discernible block in three isoforms. After application of 40 µM mexiletine, the three isoforms exhibited use-dependent block with 5-millisecond pulses (Nav1.2, 8.1% ± 1.3%, n = 5; Nav1.4, 13.5% ± 2.7%, n = 5; Nav1.5, 15.2% ± 4.5%, n = 5), and use-dependent block became bigger with 200-millisecond pulses in all isoforms (Nav1.2, 20.2% ± 2.2%, n = 5; Nav1.4, 37.2% ± 2.2%, n = 5; Nav1.5, 38.9% ± 5.4%, n = 5). The amount of activated and inactivated channel block by mexiletine was consistent with the previous reports (Ono et al., 1994; Carmeliet and Mubagwa, 1998; Sasaki et al., 2004). When we compared the amounts of use-dependent block with 5- or 200-millisecond pulses among these three isoforms, Nav1.4 and Nav1.5 showed similar sensitivity to use-dependent block, and Nav1.2 showed significantly less sensitivity than Nav1.4 or Nav1.5 (P < 0.05 or P < 0.001) (Fig. 2C).
When we compared the dissociation constants of mexiletine for three states of Nav1.2, Nav1.4, and Nav1.5, the KR values were not much different among these three isoforms (Nav1.2, KR = 323.2 ± 126.7 µM, n = 8; Nav1.4, KR = 315.1 ± 87.1 µM, n = 8; Nav1.5, KR = 270.3 ± 70.1 µM, n = 6) (Fig. 3). On the other hand, Nav1.4 and Nav1.5 showed similar affinities for the open and inactivated states, and they had a 2- to 3-fold higher affinity for these states than Nav1.2 (Nav1.2, KO = 190.5 ± 71.3 µM, n = 8, KI = 163.6 ± 52.6 µM, n = 13; Nav1.4, KO = 95.2 ± 36.4 µM, n = 5, KI = 50.9 ± 4.5 µM, n = 5; Nav1.5, KO = 92.4 ± 25.9 µM, n = 5, KI = 48.3 ± 14.9 µM, n = 6). In the three channel isoforms, the open state affinity was 1.7- to 3.3-fold higher and the inactivated state affinity was 2.0- to 6.2-fold higher than the rested state.
Dissociation constants of mexiletine for three states of Nav1.2, Nav1.4, and Nav1.5. The dissociation constants of mexiletine for the rested (KR), open (KO), and inactivated (KI) states were determined by the concentration–response curves for tonic block and for total block (tonic block + use-dependent block) caused by 5- and 200-millisecond depolarizing pulses, respectively. Data represent the mean ± S.D. from five to eight cells except Nav1.2 KI (13 cells). See Materials and Methods and the inset shown in Fig. 2A for more details of definition of block and dissociation constants.
Mutational Effects of S6 on Mexiletine Block.
Because mutations of S6 residues, Asn in domain I (N418 in Nav1.2, N434 in Nav1.4), Leu in domain III (L1465 in Nav1.2, L1280 in Nav1.4), and Phe in domain IV (F1764 in Nav1.2, F1579 in Nav1.4) affected local anesthetic binding (Ragsdale et al., 1994, 1996; Wang et al., 1998, 2000; Nau et al., 1999; Yarov-Yarovoy et al., 2001, 2002) and because most of these observations were from Nav1.2 and Nav1.4, we examined the effects of mutations of these equivalent residues (N406, L1462, F1760) on mexiletine block in Nav1.5. Substitution of Asn-406 with Ala (N406A) significantly increased tonic block (26.7% ± 3.7%, n = 7, P < 0.001) and decreased use-dependent block by 40 µM mexiletine (12.5% ± 3.7%, n = 7, P < 0.001) compared with WT (tonic block, 9.5% ± 5.8%, n = 5; use-dependent block, 38.9% ± 5.4%, n = 5) (Fig. 4). Replacement of Phe-1760 with Ala (F1760A) reduced tonic block (0.0% ± 4.5%, n = 5) and significantly decreased use-dependent block (3.2% ± 1.1%, n = 5, P < 0.001 vs. WT). Substitution of Leu-1462 with a negatively charged residue, Glu (L1462E) significantly decreased use-dependent block (11.4% ± 4.0%, n = 7, P < 0.001 vs. WT), and the Phe mutant of Leu-1462 (L1462F) significantly increased use-dependent block (49.1% ± 5.6%, n = 7, P < 0.01 vs. WT).
Effects of Nav1.5 mutants on mexiletine block. (A) Typical current traces before (Con) and after application of 40 µM mexiletine (Mex) recorded from Nav1.5 WT, N406A, L1462E, and F1760A. Current traces in each panel are during the first pulse in a train of 30 pulses for the control, and the 1st and 30th pulses for mexiletine. (B) Summary of tonic block of WT, and mutants of N406, L1462, and F1760 by 40 µM mexiletine. (C) Summary of use-dependent block of WT, and mutants of N406, L1462, and F1760 by 40 µM mexiletine. In A through C, the tonic and use-dependent blocks were produced by a train of 30 pulses with 200-millisecond duration applied at 300-millisecond interpulse intervals in the absence or presence of mexiletine. See Materials and Methods for the voltages of the pulse protocol. In B and C, data represent the mean ± S.D. from five to seven cells. ***P < 0.001; **P < 0.01 compared with WT, one-way analysis of variance followed by Dunnett’s test.
In Fig. 5, the mutational effects of N406, L1462, and F1760 on the dissociation constants of mexiletine for three states are summarized. The largest affinity changes were observed in F1760A. F1760A decreased the affinity of mexiletine for the inactivated state 9-fold (KI = 429.5 ± 147.9 µM, n = 6). This was consistent with the previous reports showing that an important part of the local anesthetic binding site includes F1764 in Nav1.2 (equivalent to F1760 in Nav1.5) (Ragsdale et al., 1994, 1996). In addition, F1760A decreased affinities for the rested and open states by 2-fold and 5-fold, respectively (KR = 515.3 ± 183.6 µM, n = 5; KO = 480.5 ± 138.4 µM, n = 6).
Effects of Nav1.5 mutants on dissociation constants of mexiletine for three states. Data represent the mean ± S.D. from five to eight cells. See Fig. 3 for definition of dissociation constants.
Changes in affinities for the open and inactivated states were also observed in L1462 mutations. L1462E decreased affinity for the open state 1.5-fold (KO = 138.6 ± 22.8 µM, n = 5) and inactivated state 2.3-fold (KI = 113.1 ± 25.0 µM, n = 5). L1462F increased affinities for the open and inactivated states by 1.8-fold and 2.3-fold, respectively (KO = 51.3 ± 13.0 µM, n = 7; KI = 21.3 ± 9.3 µM, n = 7), but L1462F had no effect on affinity for the rested state (KR = 270.8 ± 77.4 µM, n = 6). By contrast, N406A increased affinity for the rested state by 2.7-fold (KR = 101.3 ± 9.6 µM, n = 5) and decreased affinity for the inactivated state by 1.7-fold (KI = 80.4 ± 9.4 µM, n = 5).
The Pore Models of Closed and Open Nav1.5 Channel.
We built homology models of the inner pore of the closed and open Nav1.5 channel based on the crystal structures of bacterial voltage-gated Na+ channels, NavAb (for closed channel) (Payandeh et al., 2011) and NavMs (for open channel) (Sula et al., 2017). Because the transmembrane pore regions of NavAb and NavMs have high sequence identities (McCusker et al., 2012), structure differences seen in the models are likely to be state dependent rather than a consequence of sequence dissimilarities (Bagnéris et al., 2015). Our pore model was composed of P1 and P2 helices, nonhelix linkers including the selectivity filter residues between the P1 and P2 helices, and the S6 segments (Fig. 6). In our pore model of Nav1.5, F1760 in domain IV S6 and L1462 in domain III S6 faced the inner pore in both the closed and open states of the channel, which was consistent with the previous experimental evidence using a cysteine mutant of F1579 in Nav1.4 (equivalent to F1760 in Nav1.5) and the charged methanethiosulfonate reagents (Sunami et al., 2004), and many studies of local anesthetic block (Ragsdale et al., 1994, 1996; Wang et al., 2000; Yarov-Yarovoy et al., 2001; Nau et al., 2003). On the other hand, N406 in domain I S6 did not face the inner pore irrespective of the channel states.
The pore models of closed and open Nav1.5 channels based on the structures of bacterial Na+ channels, NavAb (for closed channel), and NavMs (for open channel). (A and C) Closed Nav1.5 channels. (B and D) Open Nav1.5 channels. (A and B) Top views of channels. (C and D) Side views of channels. Each pore model was composed of P1 and P2 helices, nonhelix linkers including selectivity filter residues between P1 and P2 helices, and S6 segments. The linkers between P2 helices and S6 helices are also shown, but their location and arrangement are not necessarily based on the homology sequences of NavAb and NavMs. The closed and open channels are colored light pink and silver, respectively; the helices and the linkers are shown as ribbons and rods, respectively. Residues N406 (red) in domain I S6 (DIS6), L1462 (green) in DIIIS6, and F1760 (blue) in DIVS6 are represented as sticks.
Docking Simulations of Mexiletine in the Inner Pore of Nav1.5.
We docked mexiletine in our inner pore model of the closed and open Nav1.5 channels. The lowest energy poses of mexiletine docked inside the inner pore are shown in the closed and open channels (Fig. 7). In the closed channel, mexiletine was in the middle part of the inner pore, and the drug was not close to L1462 and F1760 (Fig. 7A). In the open channel, mexiletine moved closer to the extracellular space in the inner pore and was nearer to L1462 and F1760 compared with the closed channel (Fig. 7B). In particular, mexiletine interacted with F1760 through a π–π interaction with a distance of 3.65 Å between the aromatic ring centers of the drug and F1760. These docking results were consistent with those of the present mutational studies showing that F1760 was more involved in mexiletine binding in the open and inactivated channels than in the rested channel.
Docking simulations of mexiletine inside the inner pore of the (A) closed and (B) open Nav1.5 channels. Side views of the lowest energy binding orientation of mexiletine (yellow stick) docked inside the inner pore are shown in the closed (A) and open (B) channels. L1462 and F1760 are represented as green and blue sticks, respectively. The left-hand or right-hand inset in each panel is a magnified view showing the relative location of mexiletine to L1462 and F1760, where the distances between the aromatic ring centers of F1760 and mexiletine, and between the aromatic ring center of F1760 and the tertiary amino group of mexiletine are shown.
Docking Simulations of Mexiletine in the Inner Pore of Nav1.2 and Nav1.4.
In the open and inactivated channels, Nav1.4 had similar affinity for mexiletine compared with Nav1.5, and Nav1.2 had lower affinity than Nav1.4 and Nav1.5. In the rested channels, affinities for mexiletine were not different among these isoforms. To compare the relative location of mexiletine in these isoforms, we built the models of the inner pore of Nav1.2 and Nav1.4 in a similar manner to that of Nav1.5, and then docked mexiletine inside the pore of the closed and open channels.
In the closed Nav1.2 channel, mexiletine was located relatively close to L1465 and away from F1764 (Fig. 8A). In the open Nav1.2 channel, mexiletine moved toward the extracellular space in the pore but remained away from F1764 (Fig. 8B). On the other hand, mexiletine in Nav1.4 behaved in a similar way to that in Nav1.5. That is, mexiletine was located away from L1280 and F1579 in the closed Nav1.4 channel, but mexiletine moved closer to F1760. The drug interacted with F1760 through a π–π interaction with a distance of 4.30 Å between the aromatic ring centers of the drug and F1760 in the open Nav1.4 channel (Fig. 8, C and D).
Docking simulations of mexiletine inside the inner pore of Nav1.2 and Nav1.4. Side views of the lowest energy binding orientation of mexiletine (yellow stick) docked inside the inner pore are shown in the (A and C) closed and (B and D) open channels for Nav1.2 (A and B) and Nav1.4 (C and D). Colors and magnified views of the relative location of mexiletine to L1462 and F1760 are the same as described in the legend to Fig. 7.
These results suggest that the interaction between mexiletine and the pore residues depends on the channel isoforms and the states of the channel, and that high-affinity binding of mexiletine in the open states of Nav1.4 and Nav1.5 is caused by a π–π interaction with Phe.
Docking Simulations of Mexiletine in the Inner Pore of the Open Nav1.5 Mutant Channels.
In the present mutational studies of Nav1.5, F1760A and L1462E decreased the affinities of mexiletine for the open and inactivated states. In contrast, L1462F increased drug affinities for the open and inactivated states. To examine the relative location of mexiletine in these mutant channels, we docked mexiletine inside the inner pore of the open mutant channels. In the open F1760A channel, mexiletine was located far away from A1760, and specific interactions between mexiletine and A1760 were not observed (Fig. 9A), which resulted in decreased affinities of mexiletine for the channel. In the open L1462E channel, mexiletine moved away from F1760 and interacted with E1462 through a salt bridge or hydrogen bonds (Fig. 9B). A salt bridge was formed between the positively charged tertiary amino group of mexiletine and the negatively charged carboxyl group of E1462 (Fig. 9B), and hydrogen bonds were formed between two amino group hydrogens of mexiletine and the hydroxyl oxygen of E1462 (data not shown).
Detailed view of mexiletine binding in the open Nav1.5 mutant channels. Side views of the lowest energy binding orientation of mexiletine (yellow stick) docked inside the inner pore are shown in the open channels of F1760A (A), L1462E (B), and L1462F (C). Domain II S6 is not depicted for clarity in each panel. Residues at 1760 (A1760, F1760) are shown as blue sticks, and residues at 1462 (L1462, E1462, F1462) are shown as green sticks. The distances between F1760 and substituted residues (A1760, E1462, F1462) and the aromatic ring center and the tertiary amino group of mexiletine are depicted to examine whether these pore residues make specific interactions with the drug including π–π interaction, cation–π interaction, salt bridge, and hydrogen bond.
Despite of the specific interaction of mexiletine with E1462, the observation that L1462E decreased the mexiletine block suggests that the interaction between mexiletine and E1462 was not important in the drug block; rather, less interaction of mexiletine with F1760 appeared to decrease the mexiletine block. On the other hand, in the open L1462F channel, mexiletine was sandwiched between F1760 and a newly substituted residue, Phe at 1462 (Fig. 9C). Mexiletine interacted with F1760 through a T-shaped π–π interaction with a distance of 3.40 Å and with F1462 through a π–π interaction with a distance of 4.01 Å. The increased affinities of mexiletine caused by L1462F appeared to be due to π–π interactions between the aromatic rings of mexiletine and two phenylalanine residues. These results further support the idea that a π–π interaction of mexiletine with aromatic residues in the pore is crucial for high-affinity binding of mexiletine to the open channel of Nav1.5.
Discussion
The class Ib antiarrhythmic drug mexiletine is used clinically to reduce or prevent myotonia (Jackson et al., 1994; Hudson et al., 1995) and has neuroprotective effects in models of brain ischemia (Stys and Lesiuk, 1996). To use the drug effectively for heart, muscles, or neurons, the isoform-selectivity of the channel block must be taken into account. Therefore, we identified and compared state-dependent affinities of mexiletine for three channel isoforms under the same experimental conditions.
The rested state affinities were not different among these three isoforms. Nav1.4 and Nav1.5 exhibited similar affinities for the open and inactivated states, which likely contributes to the therapeutic effects of mexiletine on myotonia as well as on cardiac arrhythmia. Nav1.2 had a 2- to 3-fold lower affinity for these states than did Nav1.4 and Nav1.5, which was consistent with the previous report showing that affinity of mexiletine for the inactivated channel was 2-fold higher in Nav1.5 than in Nav1.2 (Weiser et al., 1999).
The rested and inactivated state affinities of Nav1.4 obtained here were consistent with the previous study (Wang et al., 2004). The rested state affinity for Nav1.5 was also in good agreement with the previous report (Weiser et al., 1999). In contrast, there were big differences in open state affinities between our data and the previous data (Wang et al., 1997, 2004). They directly measured the open state affinity from the mutant channels of Nav1.4 or Nav1.5 showing the late openings (IC50 of about 3 µM). Because the WT channels of the three isoforms have no or little persistent late currents under normal conditions, we estimated the open state affinity from the total block caused by a train of 30 pulses of 5-millisecond duration applied at 300-millisecond interpulse intervals. Consequently, the open state affinities in Nav1.4 and Nav1.5 were approximately 30- to 45-fold lower than those in the mutant channels (Wang et al., 1997, 2004). More pulses with higher frequency of stimulation may explain the experimental differences between their reports and our experiments. Nonetheless, our estimation of the open state affinity seems unlikely to alter the observation that Nav1.2 has lower affinity for the open state than Nav1.4 and Nav1.5.
In our mutational studies of Nav1.5, N406A increased the affinity of mexiletine for the rested state and decreased affinity for the inactivated state, which was similar to the previous studies of etidocaine and bupivacaine exposed to the equivalent mutations of Nav1.2 and Nav1.4 (Nau et al., 1999; Yarov-Yarovoy et al., 2002). L1462E and L1462F decreased and increased the affinities of mexiletine for the open and inactivated states, respectively. These findings were similar to the previous studies of bupivacaine using the equivalent mutations of Nav1.4 (Nau et al., 2003). In our mutational studies, the biggest impact on mexiletine block was observed in F1760A. F1760A decreased affinities for the rested, open, and inactivated states 2-fold, 5-fold, and 9-fold, respectively, which resembled the results of the previous study of mexiletine in rat Nav1.5 (Weiser et al., 1999).
In our model of Nav1.5, L1462 and F1760 faced the inner pore, and N406 faced away from the pore, irrespective of the channel states. Consistent with our model, the KcsA-based model of the closed Nav1.5 channel indicated that Asp does not face the pore (Bruhova et al., 2008). Although N406A changed affinities of mexiletine for the rested and inactivated states, this affinity change seems to be mediated by allosteric influences.
In our docking, mexiletine was located differently depending on the state. In the closed states of Nav1.4 and Nav1.5, mexiletine was not close to or was located away from Phe. Once Nav1.4 and Nav1.5 channels opened, mexiletine approached Phe and formed π–π interactions. Consistent with this, F1760A decreased the affinity of mexiletine more markedly for the open state than for the rested state, suggesting that F1760 is more involved in mexiletine block in the open state than in the closed state. In addition, it is noteworthy that mexiletine is a potent open channel blocker of Nav1.4 and Nav1.5 (Wang et al., 1997, 2004).
By contrast, in the Nav1.2 model mexiletine was located away from Phe even in the open state. The simplest interpretation of these results is that mexiletine has no or little interaction with Phe in the open state of Nav1.2. To our knowledge, there is no evidence that mexiletine directly interacts with the open state of Nav1.2. In addition, our estimation of the open state affinity of mexiletine showed that Nav1.2 has a 2-fold lower affinity for the open state than Nav1.4 and Nav1.5. These results indicated that Phe is not a key residue for mexiletine block in the open state of Nav1.2. Consistent with this, Phe was also less important to mexiletine block in the rested and inactivated states of Nav1.2 than to binding of the drug to Nav1.5 (Weiser et al., 1999). Instead, another key aromatic residue, Tyr, contributed more importantly to mexiletine block in the inactivated state of Nav1.2.
The previous studies using unnatural fluorinated phenylalanine derivatives demonstrated that Nav1.4 and Nav1.5 channels form cation–π interactions with lidocaine and mexiletine at Phe and produce use-dependent block (Ahern et al., 2008; Pless et al., 2011). A cation–π interaction likely occurs in the inactivated state because class Ib antiarrhythmic drugs including lidocaine and mexiletine, having higher affinities for the inactivated state, produce use-dependent block through a cation–π interaction, but class Ia and Ic drugs including quinidine and flecainide prefer the open channel and rely less on this interaction (Pless et al., 2011). Because mexiletine seems to change its position and orientation depending on the channel state, it is possible that mexiletine binds to Phe through a π–π interaction in the open state and through a cation–π interaction in the inactivated state. Structural models of the inactivated state are not yet available to confirm this hypothesis because prokaryotic Na+ channels have no fast inactivation (Pavlov et al., 2005) and the electrophysiological states of electric eel Na+ channels with fast inactivation (of which structures were recently reported) are yet to be defined (Yan et al., 2017).
Our models agree with the following experimental data. 1) Mutations of L1462 and F1760 affect mexiletine block, and L1462 and F1760 face the pore in the model. 2) Nav1.4 and Nav1.5 have similar affinities for the closed and open states, and the relative locations of mexiletine are similar in the models of closed and open states of Nav1.4 and Nav1.5. 3) Nav1.2 has lower affinity for the open state than Nav1.4 and Nav1.5, and the relative location of mexiletine in the model of open Nav1.2 is different from those of open Nav1.4 and Nav1.5. 4) Decrease in affinity mediated by F1760A is more marked in the open state than in the closed state, and mexiletine interacts with F1760 through a π–π interaction in the model of open Nav1.5 and is not close to F1760 in the model of closed Nav1.5. 5) F1760A decreases affinity for the open state 5-fold (−1.0 kcal/mol), and mexiletine is located far away from A1760 in the model of the open channel containing F1760A.
In several previous models, local anesthetics adopted a horizontal orientation with respect to the ion permeation pathway in the closed channel and a vertical orientation in the open channel (Lipkind and Fozzard, 2005; Tikhonov et al., 2006; Tikhonov and Zhorov, 2007, 2012; Bruhova et al., 2008). On the other hand, several local anesthetics adopt horizontal, vertical, and intermediate (angular) orientations in the closed channel (Bruhova et al., 2008), and horizontal and two distinct vertical orientations in the open channel (Tikhonov and Zhorov, 2017). In our docking, mexiletine appeared not to adopt either horizontal or vertical orientations in any state. Differences in ligand orientations among the models might be, in part, because of different channel templates used, different channel states, and different channel isoforms.
In our model, mexiletine occurred closer to the extracellular space in the open channel regardless of the channel isoform. This may facilitate the electrostatic mechanism of block (Lipkind and Fozzard, 2005; Tikhonov et al., 2006; McNulty et al., 2007) and explain a greater block in the open channel. Specific binding of the aromatic ring of mexiletine to Phe likely helps locate the ammonium group upward in Nav1.4 and Nav1.5. In addition, our docking and experimental results suggest that the amino acid residues lining the pore, which are common to Nav1.4 and Nav1.5 but not common to Nav1.2 (Fig. 1), may help mexiletine form a π–π interaction with Phe in the open states. Consequently, mexiletine may work better in Nav1.4 or Nav1.5 than Nav1.2. Further experimental studies are needed to test this hypothesis.
In conclusion, Nav1.4 and Nav1.5 had similar affinities and interactions for mexiletine. This may explain the therapeutic effects of mexiletine on myotonia as well as on cardiac arrhythmia. Visualization of isoform differences in state-dependent block may contribute to rational design of isoform-selective and safer drugs.
Acknowledgments
We are indebted to Dr. H. A. Fozzard for his support of this work and to Dr. S. C. Dudley for reading and commenting on the manuscript. Dr. N. Makita generously provided the Nav1.5 clone. The technical assistance of N. Ishimaru is also appreciated.
Authorship Contributions
Participated in research design: Nakagawa, Munakata, Sunami.
Conducted experiments: Nakagawa, Munakata, Sunami.
Performed data analysis: Nakagawa, Munakata, Sunami.
Wrote or contributed to the writing of the manuscript: Nakagawa, Munakata, Sunami.
Footnotes
- Received August 18, 2018.
- Accepted December 19, 2018.
This work was in part supported by the grant from International University of Health and Welfare (to A.S.).
Abbreviations
- G418
- geneticin
- HEK293
- human embryonic kidney 293 cells
- Nav
- voltage-gated sodium channel
- WT
- wild type
- Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics