Molecular insights into the local anesthetic receptor within voltage-gated sodium channels using hydroxylated analogs of mexiletine

We previously showed that the β-adrenoceptor modulators, clenbuterol and propranolol, directly blocked voltage-gated sodium channels, whereas salbutamol and nadolol did not (Desaphy et al., 2003), suggesting the presence of two hydroxyl groups on the aromatic moiety of the drugs as a molecular requisite for impeding sodium channel block. To verify such an hypothesis, we synthesized five new mexiletine analogs by adding one or two hydroxyl groups to the aryloxy moiety of the sodium channel blocker and tested these compounds on hNav1.4 channels expressed in HEK293 cells. Concentration–response relationships were constructed using 25-ms-long depolarizing pulses at −30 mV applied from an holding potential of −120 mV at 0.1 Hz (tonic block) and 10 Hz (use-dependent block) stimulation frequencies. The half-maximum inhibitory concentrations (IC₅₀) were linearly correlated to drug lipophilicity: the less lipophilic the drug, minor was the block. The same compounds were also tested on F1586C and Y1593C hNav1.4 channel mutants, to gain further information on the molecular interactions of mexiletine with its receptor within the sodium channel pore. In particular, replacement of Phe1586 and Tyr1593 by non-aromatic cysteine residues may help in the understanding of the role of π–π or π–cation interactions in mexiletine binding. Alteration of tonic block suggests that the aryloxy moiety of mexiletine may interact either directly or indirectly with Phe1586 in the closed sodium channel to produce low-affinity binding block, and that this interaction depends on the electrostatic potential of the drug aromatic tail. Alteration of use-dependent block suggests that addition of hydroxyl groups to the aryloxy moiety may modify high-affinity binding of the drug amine terminal to Phe1586 through cooperativity between the two pharmacophores, this effect being mainly related to drug lipophilicity. Mutation of Tyr1593 further impaired such cooperativity. In conclusion, these results confirm our former hypothesis by showing that the presence of hydroxyl groups to the aryloxy moiety of mexiletine greatly reduced sodium channel block, and provide molecular insights into the intimate interaction of local anesthetics with their receptor.

Introduction

Mexiletine is a class Ib antiarrhythmic drug, also considered as the first choice drug for treating myotonic syndromes (Conte Camerino et al., 2007). Moreover, recent clinical trials demonstrate its therapeutic value to relieve neuropathic pain (Challapalli et al., 2005). Mexiletine exerts its pharmacological action through blockade of voltage-gated sodium channels, thus reducing cell excitability. Preferential binding to inactivated channels and use-dependent block are the basis of the selective action of mexiletine on pathologic hyperactive tissues. The structure of mexiletine is related to that of local anesthetic (LA) drugs, presenting a protonable amine group connected to a hydrophobic aromatic ring through an intermediate ether link (Table 1). We have shown that little modification of these three components can substantially affect sodium channel blockade in vitro and antimyotonic effects in vivo (Desaphy et al., 1999, 2001; De Luca et al., 2000, 2003, 2004; De Bellis et al., 2006). The use of mexiletine analogs has allowed to get new important information on the molecular dynamic interaction between the drug and its receptor within the sodium channel pore (De Luca et al., 2000, 2003). Besides, sodium channel mutagenesis experiments have defined the amino acids most probably involved in the interaction with LAs (Ragsdale et al., 1994, 1996; Wright et al., 1998; Nau et al., 1999; Wang et al., 2000; Sunami et al., 2001; Yarov-Yarovoy et al., 2001, 2002; O’Leary and Chahine, 2002; McNulty et al., 2007; Ahern et al., 2008).

TABLE 1

Table 1. Chemical structures and physicochemical properties.

In a previous study, we showed that drugs known as modulators of β-adrenergic receptors are able to block sodium channels in a manner reminiscent of LA drugs (Desaphy et al., 2003). The β-agonist clenbuterol and the β-blocker propranolol showed efficacies comparable to mexiletine in producing use-dependent inhibition of sodium currents in rat skeletal muscle fibers as well as in tsA201 cells transiently expressing the human muscle sodium channel isoform, hNav1.4. In contrast, the β-agonist salbutamol and the β-blocker nadolol had no effect on sodium currents at 1 mM concentration. Examination of chemical structures and physicochemical properties of these drugs revealed the presence of the two pharmacophores important for sodium channel blocking activity, the amine group and the aromatic ring, which respectively confer a high pK_a and high lipophilicity. However, the two inactive compounds were characterized by the presence of two hydroxyl groups on the aromatic moiety, which we proposed to be determinant for impeding sodium channel blockade.

To verify this hypothesis, we have synthesized new analogs of mexiletine by introducing one or two hydroxyl groups on the aryloxy moiety of the drug and tested them on sodium currents in mammalian cells permanently transfected with hNav1.4. Moreover site-directed mutagenesis was performed to introduce non-aromatic cysteine residues in place of the Phe1586 and Tyr1593 residues of hNav1.4, which are thought to be part of the binding site for LAs within the domain IV segment 6 of sodium channels (Ragsdale et al., 1994).

The results demonstrate that the presence of hydroxyl groups on the aromatic moiety drastically reduces the ability of compounds to block sodium channels and provide further insight in the molecular interaction between LA drugs and their binding site within the sodium channel pore.

Materials and Methods

Site-Directed Mutagenesis and Expression of Recombinant Sodium Channels

Full-length cDNA encoding the main α-subunit of wild-type (WT) skeletal muscle isoform of voltage-gated sodium channel (hNav1.4) was subcloned in the mammalian expression vector pRc–CMV. The two hNav1.4 mutants, F1586C and Y1593C, were engineered by standard two-step PCR-based site-directed mutagenesis. All PCR reactions were performed using Pfu DNA polymerase (Stratagene, La Jolla, CA, USA) for high-fidelity amplification (Desaphy et al., 2010). The complete coding region of channel mutant cDNAs was sequenced to exclude any polymerase errors. Permanent sodium channel expression was achieved by transfecting 1 mg ml⁻¹ of plasmid in human embryonic kidney HEK293 cells using the calcium phosphate precipitation method and clone selection using geneticin (Gibco-Invitrogen, Milan, Italy).

Sodium Current Measurement in HEK293 Cells

Whole-cell sodium currents (I_Na) were recorded at room temperature (20–22°C) using an Axopatch 1D amplifier (Axon Instruments, Union City, CA, USA). Voltage clamp protocols and data acquisition were performed with pCLAMP 9.2 software (Axon Instruments) through a 12-bit A–D/D–A interface (Digidata 1340, Axon Instruments). Pipettes made with Corning 7052 glass (Garner Glass, Claremont, CA, USA) had resistance that ranged from 1 to 3 MΩ. Currents were low-pass filtered at 2 kHz (–3 dB) by the four-pole Bessel filter of the amplifier and digitized at 10–20 kHz. After the patch membrane had been ruptured, a 25-ms-long test pulse to −30 mV from a holding potential (HP) of −120 mV was applied to the cell at a low frequency until stabilization of I_Na amplitude and kinetics was achieved (typically 5 min). Only those data obtained from cells exhibiting series resistance errors <5 mV were considered for analysis. Little (<5%) or no rundown was observed within the experiments. The concentration–response relationships were produced by obtaining the peak current amplitude measured in the presence of drug (I_DRUG), normalized to the peak current amplitude measured in the same cell before drug application (I_CTRL), as a function of drug concentration [DRUG]. Each data point is the mean ± SEM from at least three cells. The relationships were fitted with a first-order binding function:

where IC₅₀ is the half-maximum inhibitory concentration, and n_H is the slope factor. The variable I_max was introduced only for dose–response relationships of compounds effects on Y1593C channel mutant at 10 Hz, otherwise I_max was fixed to 1.

Mutant Cycle Analysis

Mutant cycle analysis was performed to evaluate the cross influence between amino acid mutations in the hNav1.4 channel and hydroxyl substitutions in mexiletine, in order to characterize the molecular interactions. The coupling energy ΔG is thus calculated as RT·lnΩ, where R is the perfect gas constant (8.314 J K mol⁻¹), T is the temperature expressed in Kelvin (295.15 K = 22°C), and Ω is the coupling constant (Hidalgo and MacKinnon, 1995).

Drugs and Solutions

Patch clamp pipette solution contained in mM: 120 CsF, 10 CsCl, 10 NaCl, 5 EGTA, and 5 HEPES, and the pH was set to 7.2 with CsOH. Bath solution for patch clamp recordings contained (in mM): 150 NaCl, 4 KCl, 2 CaCl₂, 1 MgCl₂, 5 HEPES, and 5 glucose. The pH was set to 7.4 with NaOH. Mexiletine hydrochloride [Mex, 1-(2,6-dimethylphenoxy)propan-2-amine] was purchased from Sigma (Milan, Italy). The hydroxylated analogs of mexiletine {metahydroxymexiletine, mHM, 3-(2-aminopropoxy)-2,4-dimethylphenol; parahydroxymexiletine, pHM, 4-(2-aminopropoxy)-3,5-dimethylphenol; hydroxymethylmexiletine, HMM, [2-(2-aminopropoxy)-3-methylphenyl]methanol; bis(hydroxymethyl)mexiletine, bHMM, [2-(2-aminopropoxy)-1,3-phenylene]dimethanol; parahydroxyhydroxymethyl mexiletine, pHHMM, 4-(2-aminopropoxy)-3-(hydroxymethyl)-5-methylphenol} were synthesized in our laboratories (Catalano et al., 2004, 2010; Cavalluzzi et al., 2007). All drugs were dissolved directly in external patch solution at the desired final concentration. The patched cell was continuously exposed to a stream of control or drug-supplemented bath solution flowing out from a plastic capillary.

Experimental Determination of Drug Physicochemical Properties

The pK_a and Log P values of mexiletine and its analogs were determined with a potentiometric method using Sirius GlpK_a (Sirius Analytical Instrument Ltd., Forest Row, East Sussex, UK) as described previously (Catalano et al., 2010). Because the compounds showed a low water solubility, methanol was used as co-solvent (methanol–water 10–30% w/w) for the determination of pK_a.

The electrostatic potential of each analog was determined as the molecular charge distribution of the corresponding methyl ether moiety [2,6-dimethylphenyl methyl ether, for Mex; 3-methoxy-2,4-dimethylphenol, for mHM; 4-methoxy-3,5-dimethylphenol, for pHM; (2-methoxy-3-methylphenyl)methanol, for HMM; (2-methoxybenzene-1,3-diyl)dimethanol, for bHMM; 3-(hydroxymethyl)-4-methoxy-5-methylphenol, for pHHMM] was investigated. Molecular models were constructed by fragments offered by SPARTAN PRO (Wavefunction, Inc., 18401 Von Karman Avenue, Suite 370, Irvine, CA 92715, USA) inner fragment library and assuming the suggested default starting geometries. The generated geometries were optimized by the molecular mechanics MMFF routine proposed by the software and then submitted to a systematic conformer distribution analysis. Conformers were classified according to the ab initio gas phase energy content calculated at the HF/3-21G* level. All conformers falling within a window of 5 kcal/mol above global minimum were retained. After removal of conformers differing for dihedral values lower than 10°, the retained conformers were submitted to HF/3-21G* geometry optimization. The most negative electrostatic potential (i.e., the highest energy of interaction with a positive charge used as a probe) value for the most stable conformer of each analog was used for structure–activity relationship considerations.

Results

Effects of Mexiletine and Analogs on Wild-type HNav1.4 Channels

To compare the effects of mexiletine and its hydroxylated analogs, block of sodium channels was evaluated by measuring the drug-induced reduction of I_Na elicited at −30 mV for 25 ms every 10 (0.1 Hz) or 0.1 (10 Hz) s from the HP of −120 mV. Applying this protocol in the absence of drug, there was no significant change in current amplitude for hNav1.4 channels (not shown). Typical time course of peak I_Na amplitude during application of 100 μM mexiletine is illustrated in Figure 1A. At the stimulation of 0.1 Hz, application of mexiletine reduced gradually I_Na down to a plateau. This block will be called tonic block hereafter. It is due mainly to the block of closed channels with low affinity, but depends also in part on the block of closed-state inactivated channels with high affinity (Desaphy et al., 2001, 2004). Then increasing the stimulation to 10 Hz further induced the reduction of I_Na to a lower level, determining the so-called use-dependent block, which depends on the high-affinity binding to inactivated channels. These effects were fully reversible. Examples of current traces recorded at steady state before (control) and in the presence of 100 μM mexiletine at 0.1 Hz then 10 Hz are illustrated in Figure 1B. The reduction of I_Na induced by 100 μM mexiletine was 32 and 74% at 0.1 and 10 Hz, respectively.

FIGURE 1

Figure 1. Effects of 100 μM mexiletine on hNav1.4 currents. (A) Time course of sodium current amplitude measured in a representative HEK293 cell permanently transfected with hNav1.4 sodium channel subtype before, during, and after of application of 100 μM mexiletine. The sodium current was elicited by depolarizing the cell membrane to −30 mV for 25 ms every 10 (0.1 Hz) or 0.1 s (10 Hz), from a holding potential of −120 mV. (B) Typical sodium current traces are shown, which were obtained from the average of three records obtained at steady state in absence of drug using 0.1 Hz stimulation frequency (CTRL), and in presence of drug at 0.1 Hz (red) and 10 Hz (blue).

The same protocol was used to test all the exploratory compounds at various concentrations. Figure 2 shows representative current traces recorded at 0.1 and 10 Hz stimulation during application of 1 mM of each compound. Effects of mexiletine and mHM were quite similar, blocking I_Na by 77% at 0.1 Hz and 95% at 10 Hz. All the other compounds were less potent than mexiletine, following the rank order of potency: pHM > HMM > bHMM > pHHMM. The less potent pHHMM at 1 mM had no effect at 0.1 Hz and very little effect at 10 Hz. All the effects were fully reversible (not shown). The peak I_Na amplitude measured in presence of drug at 0.1 or 10 Hz was normalized with respect to control peak I_Na and reported as a function of drug concentration to draw concentration–response curves (Figure 3). All the curves were satisfactorily fitted with a first-order binding function, given the half-maximum inhibitory concentrations (IC₅₀) reported in Table 2. Addition of one hydroxyl group to mexiletine had different effects on efficacy depending on the position of substituent on the aryloxy moiety, ranging from merely no effect (mHM) to an eightfold increase of IC₅₀ value (HMM) at 0.1 Hz. The addition of two hydroxyl groups induced a dramatic reduction of efficacy by more than 45 times for bHMM and 122 times for pHHMM. The effects of chemical maneuvers were more pronounced at 10 Hz stimulation, with an increase of IC₅₀ values ranging from 1.7-fold for mHM to 567-fold for pHHMM, compared to mexiletine. Thus the 0.1- to 10-Hz IC₅₀ ratio decreased from 9.8 for mexiletine to 6.3 for mHM and about 2.5 for the other compounds.

FIGURE 2

Figure 2. Effects of 1 mM mexiletine or hydroxylated analogs on hNav1.4 currents. Sodium currents were recorded as described in Figure 1.

FIGURE 3

Figure 3. Concentration–response relationships for mexiletine or hydroxylated analogs on hNav1.4 channels. The curves were constructed at 0.1 and 10 Hz frequency stimulation, and fitted with a first-order binding function (Eq. 1 in Materials and Methods). Each point is the mean ± SEM from at least three cells. The calculated IC₅₀ values ± SE of the fit are reported in Table 2. The slope factors ranged between 0.71 and 1.22.

TABLE 2

Table 2. Half-maximum inhibitory concentration values (IC₅₀).

The most prominent effect of hydroxylation on the physicochemical properties of compounds is a reduction of the lipophilicity at the experimental pH, expressed as the Log D determined experimentally (Table 1). The relationships of IC₅₀ values as a function of Log P or Log D, plotted on semilogarithmic axes, were linearly correlated: The higher the lipophilicity, the more efficient was the compound (Figure 4). Notwithstanding, if the lipophilicity was the unique factor influencing compound efficacy, mHM should have appeared less efficient than it was.

FIGURE 4

Figure 4. Relationship between IC₅₀ values for drug action on hNav1.4 channels and drug lipophilicity expressed either as Log P or Log D. The IC₅₀ values were determined at the holding potential of −120 mV at 0.1 Hz (A) or 10 Hz (B) stimulation frequency, as described in Figure 3. The Log P and Log D values were measured as described in Table 1. The points were linearly correlated with r² > 0.9.

Effects of Mexiletine and Derivatives on F1586C HNav1.4 Mutant

Mexiletine is thought to bind the putative high-affinity receptor for LAs within the pore of sodium channels. Indeed, the non-conservative mutation of phenylalanine to cysteine in position 1586 (hNav1.4) within the 6th segment of domain IV drastically reduced mexiletine block of inactivated channels (Desaphy et al., 2009). The hydroxylated compounds, except mHM, were tested on the same sodium channel mutant F1586C. Dose–response relationships of compound effects on F1586C I_Na are shown in Figure 5, and the IC₅₀ values are reported in Table 2. On FC channels, the compound pHM appeared as much potent as mexiletine, whereas HMM was 2.2 times less efficient, bHMM 15–20 times less, and pHHMM 20–50 times less. Compared to WT channels, the mutation reduced sodium channel block for all the compounds, especially at 10 Hz stimulation. Indeed use-dependence was almost abolished for mexiletine, pHM, HMM, and bHMM, the 0.1- to 10-Hz IC₅₀ ratio being reduced to 1.2–1.6. In contrast, use-dependent block of F1586C mutant by pHHMM was maintained.

FIGURE 5

Figure 5. Concentration–response relationships for hydroxylated mexiletine analogs on F1586C hNav1.4 channel mutants. The curves were constructed at 0.1 and 10 Hz frequency stimulation, and fitted with a first-order binding function (Eq. 1 in Materials and Methods). Each point is the mean ± SEM from at least three cells. The calculated IC₅₀ values ± SE of the fit are reported in Table 2. The slope factors ranged between 0.74 and 1.11.

Effects of Mexiletine and Derivatives on Y1593C HNav1.4 Mutant

Another amino acid potentially involved in the binding of LA-like drugs is the tyrosine in position 1593 of hNav1.4 channel (Desaphy et al., 2010). This residue was proposed to interact with the hydrophobic moiety of sodium channel blockers through a π–π interaction (Ragsdale et al., 1994), although this hypothesis has been recently challenged (Ahern et al., 2008). We tested our compounds, except mHM, on the Y1593C mutant. Representative I_Na current traces recorded before and after application of 1 mM mexiletine are shown in Figure 6. In contrast to WT and F1586C channels, the Y1593C-mutant I_Na developed significant use-dependent reduction at 10 Hz in absence of drug, ranging from 15 to 30% of control current depending on the cell. The dose–response relationships are shown in Figure 6. To take into account the use-dependent reduction of Y1593C-mutant I_Na in absence of drugs, an additional variable, namely I_max, was introduced in the fitting equation. The IC₅₀ and I_max values are reported in Table 2. Regarding the IC₅₀ values at 0.1 Hz, the Y1593C mutation was less disturbing than F1586C. Compared to WT, the IC₅₀ of mexiletine on Y1593C I_Na at 0.1 Hz was twofold higher, whereas no significant change was observed for the other compounds. In contrast, the Y1593C mutation reduced dramatically use-dependent block at 10 Hz, increasing IC₅₀ sevenfold for mexiletine, and 1.5- to 2-fold for the other compounds, except for pHHMM. Compared to WT, the 0.1- to 10-Hz IC₅₀ ratio decreased to 2.8 for mexiletine and 1–1.5 for the other compounds. No significant effect of Y1593C mutation was observed for pHHMM.

FIGURE 6

Figure 6. Concentration–response relationships for mexiletine and hydroxylated analogs on Y1593C hNav1.4 channel mutants. (A) Representative sodium current traces recorded in a HEK293 cell permanently transfected with Y1593C hNav1.4 mutant, by depolarizing the cell to −30 mV from an holding potential of −120 mV at 0.1 and 10 Hz stimulation frequency, before and during application of 1 mM mexiletine. Note the use-dependent current reduction in absence of drug. (B–F) The concentration–response relationships were constructed at 0.1 and 10 Hz frequency stimulation, and fitted with a first-order binding function (Eq. 1 in Materials and Methods). Each point is the mean ± SEM from at least three cells. The calculated IC₅₀ and I_max values ± SE of the fit are reported in Table 2. The slope factors ranged between 0.69 and 1.76.

Mutant Cycle Analysis

Being modified solely on the aryloxy moiety, the derivatives of mexiletine can be considered as simple mexiletine mutants. Thus the combined evaluation of effects of mexiletine mutants on hNav1.4 channel mutants allows the calculation of mutant cycles (Hidalgo and MacKinnon, 1995; Lipkind and Fozzard, 2005). By this way, the cross influence of one mutation on the effect of the other can be quantified using a thermodynamic cycle in order to characterize the intimate molecular interactions. Such a mutant cycle analysis is illustrated for pHM in Figure 7, and coupling energies are reported in Table 3. Since the tonic block is determined mainly by interaction of drugs with the closed channel, mutant cycle calculation using the IC₅₀ values for tonic block likely reveals information about interaction of drugs with the low-affinity binding site (Figure 7A). In these conditions, the coupling energy (ΔG) for interaction with F1586 is 1 kcal/mol, whereas ΔG is 0.5 kcal/mol for interaction with Y1593, thereby suggesting that mexiletine is closer to F1586 than to Y1593 in the closed channel. Since OH substitutions reduce ΔG for F1586 but leave unchanged ΔG for Y1593, we may hypothesize that the former amino acid may interact directly or indirectly with the aryloxy mexiletine moiety whereas the latter may not be involved in low-affinity binding. On the other hand, mutant cycle calculation using the IC₅₀ values for use-dependent block (10 Hz stimulation) may allow to get information about the interaction of drugs with the high-affinity binding site (Figure 7B). In this case, the coupling energy between the drug and Y1593 is 0.9 kcal/mol, whereas ΔG is ∼1.8 kcal/mol for F1586, which suggests that mexiletine may interact with both amino acids, and that interaction with F1586 is more critical to high-affinity binding.

FIGURE 7

Figure 7. Mutant cycle analysis of mexiletine and pHM binding to the wild-type hNav1.4 channel and its mutants Y1593C and F1586C. Mutant cycle analysis for tonic block considers IC₅₀ values calculated at 0.1 Hz stimulation frequency (A), while IC₅₀ calculated at 10 Hz stimulation frequency allowed mutant cycle analysis for use-dependent block (B). The Xi and Yi values represent the changes in affinity with each mutation, calculated as the ratio of IC₅₀ values. The coupling constant Ω is calculated as the X1 to X2 ratio and the coupling energy ΔG is calculated as RT·lnΩ, where R is the perfect gas constant (8.314 J K mol⁻¹) and T the temperature expressed in Kelvin (295.15 K = 22°C). The calculated ΔG shown in red was converted to kilocalorie per mole.

TABLE 3

Table 3. Coupling energy (ΔG, kcal/mol) calculated from mutant cycles, as shown in Figure 7.

Discussion

In this study, we used chemical maneuvers on mexiletine and site-directed mutagenesis on hNav1.4 voltage-gated sodium channel to get further insight in the intimate drug-channel interaction at the LA receptor. The binding affinity of LA drugs depends on the state of the channel, being higher for the open/inactivated with respect to the closed channel (Hille, 1977; Fozzard et al., 2011). Please note here that tonic block measured at the stimulation of 0.1 Hz results from the combination of low-affinity drug binding to closed channels and, to a lesser extent, high-affinity binding to closed-state inactivated channels, while use-dependent block measured at 10 Hz stimulation depends mainly on the high-affinity binding to inactivated channels (Desaphy et al., 2001, 2004).

Lipophilicity of Drugs is a Key Determinant of Sodium Channel Blockade

The results first confirm our former hypothesis (Desaphy et al., 2003), that the presence of two hydroxyl groups on the aryloxy moiety of LA-like drugs impedes high-affinity sodium channel blockade. We also show that the presence of a single hydroxyl group is less deleterious, depending however on the position of the OH group. The good correlation between the log D of compounds and the IC₅₀ values, either at 0.1 or 10 Hz, suggests that compound lipophilicity is a major parameter. Because the putative LA receptor lies within the channel pore, lipophilicity might be determinant allowing the compound to reach/escape its receptor through the closed or inactivated channel. For instance, we previously showed that the drug pilsicainide, which has a negative log D, blocks channels inactivated from the closed state with very low efficiency compared to mexiletine, although affinities for inactivated channels of both drugs are comparable (Desaphy et al., 2010). Thus altering hydrophobicity of compounds can modify drug access to or escape from the binding site. A reduced hydrophobicity is expected to reduce drug access to closed channels, which will result in a reduced apparent affinity for the closed and closed-inactivated states. Such a mechanism likely contributes significantly to the effects observed in the present study. Nevertheless, such an explanation is not completely sufficient because inhibition of sodium current at 10 Hz stimulation is greatly reduced as well as use-dependence expressed as 0.1- to-10 Hz IC₅₀ ratio, arguing against a trapping mechanism of the more hydrophilic compounds within the closed/inactivated channel (O’Leary et al., 2003). This result indicates that lipophilicity is also a critical determinant of the intimate interaction between the compound and its receptor. We thus tried to get more information by using contemporaneously mexiletine analogs and sodium channel mutants.

Effects of F1586C Mutation on Sodium Channel Blockade by Mexiletine and Hydroxylated Derivatives

The phenylalanine in position 1586 is thought to be a key component of LA drug binding to the LA receptor (Ragsdale et al., 1994; Sheets et al., 2010). The aromaticity of this residue is required for use-dependent inhibition, while its hydrophobicity may be relevant for resting block (Li et al., 1999). By incorporating unnatural fluorinated phenylalanine analogs in place of the natural phenylalanine, it was recently demonstrated that disruption of the electrostatic potential at this position abolishes use-dependent inhibition without affecting resting block by lidocaine and mexiletine (Ahern et al., 2008; Pless et al., 2011). Since the fluorination did not affect blockade by the neutral benzocaine, the authors concluded that a π–cation interaction between phenylalanine side chain and the charged amine of drugs is a key player in high-affinity binding, at least for class Ib antiarrhythmics such as mexiletine (Pless et al., 2011). Accordingly, we observed that the mutation of F1586 by a non-aromatic cysteine nearly zeroed use-dependent inhibition by mexiletine and its derivatives, except pHHMM. In addition, we found that low-affinity block by mexiletine was also lightly affected by F1586C, because the IC₅₀ value calculated at 0.1 Hz stimulation was 50% greater than mexiletine affinity for WT channels in the closed state (∼800 μM, Desaphy et al., 2001).

Interestingly, the inhibition of F1586C by pHM and mexiletine was comparable, indicating that the presence of one hydroxyl group in pHM mainly reduced high-affinity, use-dependent block of WT channels (as in mHM). Considering the π–cation interaction with F1586 as the main determinant of high-affinity binding, this suggests that alteration of the aryloxy moiety by introduction of the OH in mHM and pHM is able to weaken the critical interaction of the charged amine with F1586, through a mechanism resembling cooperativity between two pharmacophores. Our results suggest that this may be related to the reduction of drug lipophilicity (see Log D value in Table 1). Nevertheless, it has been recently shown that addition of a methyl group in para position of lidocaine is also able to weaken the π–cation interaction with F1586 (see supplemental data in Pless et al., 2011), suggesting that changes other than lipophilicity may be also involved. As hypothesized from previous experiments using aryloxy-substituted mexiletine analogs (Franchini et al., 2003), a possibility is that the presence of substituents on the aromatic ring may slightly modify the orientation of the molecule, thereby weakening the π–cation interaction. It is also interesting that while the overall hydrophobicity does not change between mHM and pHM, the IC₅₀ for use-dependent block of WT channels decreases 10-fold with the OH in para position, thereby confirming that not only lipophilicity is involved. Thus, in function of its position, the substituent may modify more or less the orientation of the entire molecule within the binding site. It is possible to appreciate from Figure 8 how OH addition in meta have limited effects on aryloxy hindrance and on electrostatic distribution within the aryloxy (compared to mexiletine), whereas the OH in para concentrates more electronegativity close to the molecule extremity and induces a major steric hindrance. We may speculate that, in case of π–π stacking between the aryloxy and Y1593 side chain (see below), the OH in para position as in pHM may interact very negatively with the OH in Y1593, whereas the OH in meta position may be less unfavorable. The former negative interaction may result from steric hindrance or local electrostatic repulsion between the two OH.

FIGURE 8

Figure 8. Calculation of drug electrostatic potential. Mexiletine and its hydroxylated analogs are displayed as chemical structures presenting in bold the moiety studied by ab initio calculations (upper panels). For each compound, both the corresponding HF/3-21G* minimized three-dimensional structure and the electrostatic potential map is given in equatorial and lower panels, respectively. For each moiety, the scale bar indicates the electrostatic potential surface value range in kilocalorie per mole, the negative value being the one used in structure–activity evaluations. The HMM, bHMM, and pHHMM moieties (lower half of the figure) presented the most negative electrostatic potential values (around −65 kcal/mol), which most likely may allow less favorable interactions with electron-rich aromatic rings (e.g., F1586 or Y1593 rings), compared to the ones expected for Mex, mHM, and pHM moieties with a less negative electrostatic potential (around −50 kcal/mol). The Log D values experimentally measured at pH 7.4 are indicated.

When the strong interaction between the charged amine group and F1586 is disrupted by cysteine substitution, pHM has nearly the same activity than mexiletine even if it displays a reduced Log D. This latter observation suggests that an additional OH, as in mHM and pHM, or substitution of F1586 have no additive effect on low-affinity binding. Such an observation would be in accord with the proposed horizontal orientation of LA drugs in the closed channel, where the aromatic moiety of drugs may be exposed toward F1586 (Bruhova et al., 2008; Hanck et al., 2009; Sheets et al., 2010). Thus the alteration of one or the other may determine similar effect on low-affinity binding.

Interestingly, we noted that mHM and pHM show an electrostatic potential similar to that of mexiletine (Figure 8). In addition we calculated the value of energy coupling for tonic block calculated from mexiletine/pHM and F1586C mutants cycle, which resulted close to 1 kcal/mol (Table 3), being in the range expected for a through-space electrostatic coupling (Hidalgo and MacKinnon, 1995). Altogether, these results suggest that the aromatic moiety of the drug may interact with F1586 in the closed channel through π–π stacking thereby accounting for low-affinity binding. It should be noted however that, in contrast with this hypothesis, the variation of F1586 electrostatic potential by fluorination was shown to have no effect on low-affinity binding of lidocaine (Ahern et al., 2008). Thus an alternative hypothesis to be considered may be an indirect effect of F1586C mutation, through an allosteric mechanism, on the electrostatic coupling between the aryloxy drug moiety and another aromatic residue in the closed pore.

Regarding HMM and bHMM, a similar impairment of positive cooperativity between the aromatic and protonable drug extremities may account for the reduction of high-affinity blockade of WT channels, and use-dependent block was also nearly zeroed by the F1586C mutation. Moreover, tonic block of WT channels by HMM and bHMM was greatly reduced compared to that exerted by mexiletine. Thus the IC₅₀ values for tonic block by HMM and bHMM were more than 2- and 10-fold higher than the specific affinity of mexiletine for closed channels, respectively (∼800 μM, Desaphy et al., 2001), indicating that these compounds significantly affect the low-affinity binding. Interestingly, we observed that the electrostatic potentials of the aryloxy moiety of HMM and bHMM were 10-kcal/mol more negative than that of mexiletine, mHM, and pHM. We may thus hypothesized that alteration of electrostatic coupling may contribute to weaken the interaction of the aromatic tail of HMM/bHMM with the low-affinity receptor. Such a hypothesis is supported by the low energy coupling value (<1 kcal/mol) for tonic block by HMM and bHMM. Most probably, other parameters, including hydrophobicity, should be however involved to explain the large difference in tonic block exerted by HMM and bHMM. In addition, we observed that both HMM and bHMM exert a lower tonic block of F1586C channel with respect to mexiletine and pHM, indicating that substituents on HMM and bHMM have additional effects to those resulting from the suppression of aromatic side chain at position 1586.

The compound pHHMM shows an electrostatic potential similar to HMM/bHMM, and tonic block was similarly reduced twofold by F1586C. Quite surprisingly, the mutation had no significant effect on pHHMM use-dependence, indicating that the drug does not bind to F1586 at all. This latter observation suggests that F1586 is no more involved in the higher-affinity binding of pHHMM.

Effects of Y1593C Mutation on Sodium Channel Blockade by Mexiletine and Hydroxylated Derivatives

The tyrosine in position 1593 was originally individuated as an important component of the intimate drug-receptor interaction, possibility through hydrophobic interaction with the aromatic moiety of LA drugs (Ragsdale et al., 1994). This hypothesis has been somewhat challenged because Y1593 mutations can modify sodium channel gating, which rends difficult the interpretation of pharmacological results. In addition to various effects on steady-state activation and inactivation, a slowing of recovery from inactivation has been reported (O’Reilly et al., 2000; Xiao et al., 2001), which may explain the use-dependent reduction of Y1593C currents we observed in the absence of drug. Recently, the contribution of Y1593 to π–cation or π–π interactions with lidocaine has been ruled out by experiments using fluorinated amino acids (Ahern et al., 2008).

As for pilsicainide in a previous study (Desaphy et al., 2010), we observed that Y1593C reduced use-dependent sodium current blockade by mexiletine but to a lesser extent than F1586C.

The mutation may impede a direct interaction of Y1593 with the drug, or reduce the lipophilicity in proximity of the binding site. Thus considering Y1593 being involved either directly or indirectly on mexiletine binding, its mutation to C1593 may cooperatively reduce the strong interaction between the charged amine and F1586. The Y1593C mutation had a similar effect on mexiletine derivatives (except pHHMM), reducing use-dependence without altering tonic block. This observation indicates that Y1593 is not involved in the low-affinity binding, still in accord with the proposed position of LA drugs in the closed channel (Sheets et al., 2010). Tonic block and use-dependence of pHHMM was not altered by Y1593C, most probably because hydrophobicity is no more critical for a drug having a very negative Log D.

Synopsis

We made an attempt to summarize the main findings in a graphical sketch (Figure 9). Regarding the low-affinity binding to closed channels (Figure 9A), our results suggest that the aryloxy moiety of drugs may interact either directly or indirectly with F1586 and suggest that electrostatic coupling is critical in drug binding. Regarding the high-affinity binding to inactivated channels (Figure 9B), addition of OH groups on the aryloxy moiety cooperatively reduced the interaction of drug charged amine group with F1586. The amplitude of this effect depends mainly on lipophilicity and also on the position of the OH substituent.

FIGURE 9

Figure 9. Synopsis of possible interactions of mexiletine and hydroxylated analogs with the low (A) and high (B) binding site within the wild-type hNav1.4 channel and its mutants F1586C and Y1593C. (A) Regarding the low-affinity binding to closed channels, our results are in accord with an horizontal orientation of the drug, with the aryloxy moiety interacting either directly or indirectly with F1586 through electrostatic coupling. The addition of OH groups as in mHM and pHM does not modify aryloxy electrostatic potential thereby leaving this interaction unaltered. The F1586C mutation weakens interaction with Mex, mHM, and pHM to the same extent (arrow 1). Introduction of OH as in HMM and bHMM weakens the interaction due to modification of aryloxy electrostatic potential and possibly other parameters (arrow 2). The F1586C mutation produces additional weakening of the interaction with HMM and bHMM (arrow 3). Finally the Y1593C mutation has no effect on low-affinity binding (not illustrated). (B) Regarding the high-affinity binding to inactivated channels, addition of OH groups to the aryloxy moiety cooperatively reduced the interaction of charged amine group with F1586 (arrow 1). The amplitude of this effect depends on the OH position and lipophilicity. The F1586C mutation disrupts the π–cation interaction, so that OH substitutions have no more influence on binding (arrow 2 and 3). In contrast, the Y1593C mutation impairs only partially the π–cation interaction of mexiletine through cooperativity (arrow 4). The mutation also further impairs the binding of hydroxylated analogs (arrow 5).

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.

Acknowledgments

The financial support from Telethon, Italy (non-for-profit, grants GGP04140 and GGP10101) and Association Française contre les Myopathies (non-for-profit, grant #15020) to Diana Conte Camerino is gratefully acknowledged.

References

Ahern, C. A., Eastwood, A. L., Dougherty, D. A., and Horn, R. (2008). Electrostatic contributions of aromatic residues in the local anesthetic receptor of voltage-gated sodium channels. Circ. Res. 102, 86–94.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bruhova, I., Tikhonov, D. B., and Zhorov, B. S. (2008). Access and binding of local anesthetics in the closed sodium channel. Mol. Pharmacol. 74, 1033–1045.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Catalano, A., Carocci, A., Cavalluzzi, M. M., Di Mola, A., Lentini, G., Lovece, A., Dipalma, A., Costanza, T., Desaphy, J.-F., Conte Camerino, D., and Franchini, C. (2010). Hydroxylated analogs of mexiletine as tools for structural-requirements investigation of the sodium channel blocking activity. Arch. Pharm. (Weinheim) 343, 325–332.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Catalano, A., Carocci, A., Fracchiolla, G., Franchini, C., Lentini, G., Tortorella, V., De Luca, A., De Bellis, M., Desaphy, J.-F., and Conte Camerino, D. (2004). Stereospecific synthesis of “para-hydroxymexiletine” and sodium channel blocking activity evaluation. Chirality 16, 72–78.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cavalluzzi, M. M., Catalano, A., Bruno, C., Lovece, A., Carocci, A., Corbo, F., Franchini, C., Lentini, G., and Tortorella, V. (2007). Synthesis of (R)-, (S)-, and (RS)-hydroxymethylmexiletine, one of the major metabolites of mexiletine. Tetrahedron Asymmetry 18, 2409–2417.

CrossRef Full Text

Challapalli, V., Tremont-Lukats, I. W., McNicol, E. D., Lau, J., and Carr, D. B. (2005). Systemic administration of local anesthetic agents to relieve neuropathic pain. Cochrane Database Syst. Rev. 4, CD003345.

Pubmed Abstract | Pubmed Full Text

Conte Camerino, D., Tricarico, D., and Desaphy, J.-F. (2007). Ion channel pharmacology. Neurotherapeutics 4, 184–198.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Bellis, M., De Luca, A., Rana, F., Cavalluzzi, M. M., Catalano, A., Lentini, G., Franchini, C., Tortorella, V., and Conte Camerino, D. (2006). Evaluation of the pharmacological activity of the major mexiletine metabolites on skeletal muscle sodium currents. Br. J. Pharmacol. 149, 300–310.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Luca, A., Natuzzi, F., Desaphy, J.-F., Loni, G., Lentini, G., Franchini, C., Tortorella, V., and Conte Camerino, D. (2000). Molecular determinants of mexiletine structure for potent and use-dependent block of skeletal muscle sodium channels. Mol. Pharmacol. 57, 268–277.

Pubmed Abstract | Pubmed Full Text

De Luca, A., Pierno, S., Liantonio, A., Desaphy, J.-F., Natuzzi, F., Didonna, M. P., Ferrannini, E., Jockusch, H., Franchini, C., Lentini, G., Corbo, F., Tortorella, V., and Conte Camerino, D. (2004). New potent mexiletine and tocainide analogues evaluated in vivo and in vitro as antimyotonic agents on myotonic ADR mouse. Neuromuscul. Disord. 14, 405–416.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Luca, A., Talon, S., De Bellis, M., Desaphy, J.-F., Franchini, C., Lentini, G., Catalano, A., Corbo, F., Tortorella, V., and Conte Camerino, D. (2003). Inhibition of skeletal muscle sodium currents by mexiletine analogues: specific hydrophobic interactions rather than lipophilia per se account for drug therapeutic profile. Naunyn Schmiedebergs Arch. Pharmacol. 367, 318–327.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Desaphy, J.-F., Conte Camerino, D., Franchini, C., Lentini, G., Tortorella, V., and De Luca, A. (1999). Increased hindrance on the chiral carbon atom of mexiletine enhances the block of rat skeletal muscle Na+ channels in a model of myotonia induced by ATX. Br. J. Pharmacol. 128, 1165–1174.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Desaphy, J.-F., De Luca, A., Didonna, M. P., George, A. L. Jr., and Conte Camerino, D. (2004). Different flecainide sensitivity of hNav1.4 channels and myotonic mutants explained by state-dependent block. J. Physiol. 554, 321–334.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Desaphy, J.-F., De Luca, A., Tortorella, P., De Vito, D., George, A. L. Jr., and Conte Camerino, D. (2001). Gating of myotonic Na channel mutants defines the response to mexiletine and a potent derivative. Neurology 57, 1849–1857.

Pubmed Abstract | Pubmed Full Text

Desaphy, J.-F., Dipalma, A., Costanza, T., Bruno, C., Lentini, G., Franchini, C., George, A. L. Jr., and Conte Camerino, D. (2010). Molecular determinants of state-dependent block of voltage-gated sodium channels by pilsicainide. Br. J. Pharmacol. 160, 1521–1533.

Pubmed Abstract | Pubmed Full Text

Desaphy, J.-F., Dipalma, A., De Bellis, M., Costanza, T., Gaudioso, C., Delmas, P., George, A. L. Jr., and Conte Camerino, D. (2009). Involvement of voltage-gated sodium channels blockade in the analgesic effects of orphenadrine. Pain 142, 225–235.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Desaphy, J.-F., Pierno, S., De Luca, A., Didonna, P., and Conte Camerino, D. (2003). Different ability of clenbuterol and salbutamol to block sodium channels predicts their therapeutic use in muscle excitability disorders. Mol. Pharmacol. 63, 659–670.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fozzard, H. A., Sheets, M. F., and Hanck, D. A. (2011). The sodium channel as a target of local anesthetic drugs. Front. Pharmacol. 2:68.

CrossRef Full Text

Franchini, C., Carocci, A., Catalano, A., Cavalluzzi, M. M., Corbo, F., Lentini, G., Scilimati, A., Tortorella, P., Conte Camerino, D., and De Luca, A. (2003). Optically active mexiletine analogues as stereoselective blockers of voltage-gated Na+ channels. J. Med. Chem. 46, 5238–5248.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hanck, D. A., Nikitina, E., McNulty, M. M., Fozzard, H. A., Lipkind, G. M., and Sheets, M. F. (2009). Using lidocaine and benzocaine to link sodium channel molecular conformations to state-dependent antiarrhythmic drug affinity. Circ. Res. 105, 492–499.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hidalgo, P., and MacKinnon, R. (1995). Revealing the architecture of a K⁺ channel pore through mutant cycles with a peptide inhibitor. Science 268, 307–310.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hille, B. (1977). Local anesthetics: hydrophilic and hydrophobic pathways for the drug-receptor reaction. J. Gen. Physiol. 69, 497–515.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, H. L., Galue, A., Meadows, L., and Ragsdale, D. S. (1999). A molecular basis for the different local anesthetic affinities of resting versus open and inactivated states of the sodium channel. Mol. Pharmacol. 55, 134–141.

Pubmed Abstract | Pubmed Full Text

Lipkind, G. M., and Fozzard, H. A. (2005). Molecular modelling of local anesthetic drug binding by voltage-gated sodium channels. Mol. Pharmacol. 68, 1611–1622.

Pubmed Abstract | Pubmed Full Text

McNulty, M. M., Edgerton, G. B., Shah, R. D., Hanck, D. A., Fozzard, H. A., and Lipkind, G. M. (2007). Charge at the lidocaine binding site residues Phe-1759 affects permeation in human cardiac voltage-gated sodium channels. J. Physiol. 581, 741–755.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Nau, C., Wang, S.-H., Strichartz, G. R., and Wang, G. K. (1999). Point mutations at N434 in D1-S6 of m1 Na+ channels modulate binding affinity and stereoselectivity of local anesthetic enantiomers. Mol. Pharmacol. 56, 404–413.

Pubmed Abstract | Pubmed Full Text

O’Leary, M. E., and Chahine, M. (2002). Cocaine binds to a common site on open and inactivated human heart (Nav1.5) sodium channels. J. Physiol. 541, 701–716.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

O’Leary, M. E., Digregorio, M., and Chahine, M. (2003). Closing and inactivation potentiate the cocaethylene inhibition of cardiac sodium channels by distinct mechanisms. Mol. Pharmacol. 64, 1575–1585.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

O’Reilly, J. P., Wang, S.-Y., and Wang, G. K. (2000). A point mutation in domain 4-segment 6 of the skeletal muscle sodium channel produces an atypical inactivation state. Biophys. J. 78, 773–784.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pless, S. A., Galpin, J. D., Frankel, A., and Ahern, C. A. (2011). Molecular basis for class Ib anti-arrhythmic inhibition of cardiac sodium channels. Nat. Commun. 2, 351.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ragsdale, D. S., McPhee, J. C., Scheuer, T., and Catterall, W. A. (1994). Molecular determinants of state-dependent block of Na⁺ channels by local anesthetics. Science 265, 1724–1728.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ragsdale, D. S., McPhee, J. C., Scheuer, T., and Catterall, W. A. (1996). Common molecular determinants of local anesthetic, antiarrhythmic, and anticonvulsant block of voltage-gated sodium channels. Proc. Natl. Acad. Sci. U.S.A. 93, 9270–9275.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sheets, M. F., Fozzard, H. A., Lipkind, G. M., and Hanck, D. A. (2010). Sodium channel molecular conformations and antiarrhythmic drug affinity. Trends Cardiovasc. Med. 20, 16–21.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sunami, A., Glaase, J. W., and Fozzard, H. A. (2001). Structural and gating changes of the sodium channel induced by mutation of a residue in the upper third of IVS6, creating an external access path for local anesthetics. Mol. Pharmacol. 59, 684–691.

Pubmed Abstract | Pubmed Full Text

Wang, S. Y., Nau, C., and Wang, G. K. (2000). Residues in Na+ channel D3-S6 segment modulate both batrachotoxin and local anesthetic affinities. Biophys. J. 79, 1379–1387.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wright, S. N., Wang, S. Y., and Wang, G. K. (1998). Lysine point mutations in Na+ channel D4-S6 reduce inactivated channel block by local anesthetics. Mol. Pharmacol. 54, 733–739.

Pubmed Abstract | Pubmed Full Text

Xiao, Y.-F., Ke, Q., Wang, S.-Y., Yang, Y., Wang, G. K., Morgan, J. P., and Leaf, A. (2001). Point mutations in α-subunit of human cardiac Na+ channels alter Na+ current kinetics. Biochem. Biophys. Res. Commun. 281, 45–52.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yarov-Yarovoy, V., Brown, J., Sharp, E. M., Clare, J. J., Scheuer, T., and Catterall, W. A. (2001). Molecular determinants of voltage-dependent gating and binding of pore-blocking drugs in transmembrane segment IIIS6 of the Na+ channel α subunit. J. Biol. Chem. 276, 20–27.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yarov-Yarovoy, V., McPhee, J. C., Idsvoog, D., Pate, C., Scheuer, T., and Catterall, W. A. (2002). Role of amino acid residues in transmembrane segments IS6 and IIS6 of the Na+ channel α subunit in voltage-dependent gating and drug block. J. Biol. Chem. 277, 35393–35401.

CrossRef Full Text