- 1Department of Anesthesiology, University of Arizona, Tucson, AZ, United States
- 2Superior Institute of Biomedical Sciences, State University of Ceará, Fortaleza, Brazil
- 3Department of Physiology and Biophysics, Biomedical Sciences Institute, University of Sao Paulo, São Paulo, Brazil
Eugenol (EUG) is a bioactive monoterpenoid used as an analgesic, preservative, and flavoring agent. Our new data show EUG as a voltage-gated Na+ channel (VGSC) inhibitor, comparable but not identical to lidocaine (LID). EUG inhibits both total and only TTX-R voltage-activated Na+ currents (INa) recorded from VGSCs naturally expressed on dorsal root ganglion sensory neurons in rats. Inhibition is quick, fully reversible, and dose-dependent. Our biophysical and pharmacological analyses showed that EUG and LID inhibit VGSCs with different mechanisms. EUG inhibits VGSCs with a dose–response relationship characterized by a Hill coefficient of 2, while this parameter for the inhibition by LID is 1. Furthermore, in a different way from LID, EUG modified the voltage dependence of both the VGSC activation and inactivation processes and the recovery from fast inactivated states and the entry to slow inactivated states. In addition, we suggest that EUG, but not LID, interacts with VGSC pre-open–closed states, according to our data.
Introduction
Phytochemicals are a good source of molecules to drive the discovery of new bioactive compounds. Terpenes, molecules formed by isoprene molecule condensation, are the largest class of phytochemicals, with over 25,000 molecules identified (Harrewijn et al., 2012; Kozioł et al., 2014). Cyclic monoterpenoids, found as components of the aromatic essential oil of many plants, are cycled and oxygenated molecules based on two isoprene unit monoterpenes (Clarke and Clarke, 2008). Cyclic monoterpenoids possess a broad range of biological activities, both in vitro and in vivo. These activities include analgesic (Jorkjend and Skoglund, 1990; Ohkubo and Shibata, 1997; Park et al., 2011; Guimarães et al., 2013), anti-arrhythmogenic (Sensch et al., 2000; Magyar et al., 2004; Tomasova et al., 2015; Binu et al., 2017), antiepileptogenic (Müller et al., 2006; Jeong et al., 2015; Sucher and Carles, 2015), anticonvulsant (Dallmeier and Carlini, 1981; Dallmeier et al., 1983; de Almeida et al., 2011; Nóbrega et al., 2014; Sancheti et al., 2014), and myorelaxant effects (Beer et al., 2007; Lima et al., 2011; Olivoto et al., 2014; Peixoto-Neves et al., 2014; Peixoto-Neves et al., 2015), and they were corroborated by mechanisms of action studies strongly suggesting that the effects mentioned above, at least in part, are caused by interactions between these cyclic monoterpenoid molecules and ion channels expressed in the membranes of excitable cells. Vast literature studies show that cyclic monoterpenoids block and/or modulate voltage-gated sodium channels (VGSCs) (Haeseler et al., 2002; Park et al., 2006; Cho et al., 2008; Park et al., 2009; Moreira-Lobo et al., 2010; Gaudioso et al., 2012; Joca et al., 2012; Silva-Alves et al., 2013; Wang et al., 2015; Teixeira-Fonseca et al., 2021), voltage-gated calcium channels (Magyar et al., 2004; Lee et al., 2005; Soares et al., 2007; Chung et al., 2008; Seo et al., 2013), and voltage-gated potassium channels (Sensch et al., 2000; Magyar et al., 2002; dos Santos-Nascimento et al., 2015). In addition, somatosensory-related transient receptor potential cation channel subfamily A member 1 (TRPA1) (Karashima et al., 2007; Lee et al., 2008; Chung et al., 2014; Takaishi et al., 2014) and other channels are also reported to be modulated by cyclic monoterpenoids (Hall et al., 2004; Li et al., 2008; Yeon et al., 2011; Lee et al., 2015).
VGSCs are membrane proteins comprising an ion-conductive and voltage-sensitive alpha subunit and a varying number of regulatory beta subunits (Hille, 2001; Brackenbury and Isom, 2011; Bouza and Isom, 2018). Most of the known VGSC blockers and inhibitors are directed to the alpha subunit since it possesses the voltage-sensing mechanism and the coupled Na+-permeable conductive pore. In mammalian organisms, there are nine different VGSC alpha subunits that share more than 75% identity. Expressed in the cell membrane, each VGSC alpha subunit contains more than 2,000 amino acids and 24 membrane-spanning segments organized in 4 homologous domains radially disposed for a Na+-selective conducting pore in the center. In each domain, transmembrane segments S1–S4 form the voltage sensor that provides voltage dependence to the open probability of the central pore that is formed by the S5 and S6 segments from all domains (Pan et al., 2018). VGSCs are the most important ion channels for cell excitability since they are responsible for initiating action potentials in neurons and muscle cells (Hille, 2001; Catterall, 2012). VGSCs are crucial to maintaining physiologic cell excitability in neurons, myocytes, and endocrine cells. Natural variants of VGSCs are major causes of diseases such as epilepsy, periodic paralysis, arrhythmias, and pain disorders (Hodgkin and Huxley, 1952; Hille, 2001; George, 2005; Catterall, 2017; Jiang et al., 2022; Hernandez and Richards, 2023).
Therefore, the development of new blockers and inhibitors for VGSCs could potentially lead to the discovery of new therapies for many diseases of cell excitability (Theile and Cummins, 2011; Meng et al., 2014; Bagal et al., 2015). Currently, there is a need for VGSC inhibitors or modulators that would target the right subtype of these channels for tissue and cell specificity.
The present work aimed at speculating a putative state-dependent inhibition of VGSCs expressed by cultured dorsal root ganglion (DRG) neurons by eugenol (4-allyl-2-methoxyphenol, EUG), a bioactive cyclic monoterpenoid that is traditionally used as a dentistry material in humans in dental temporary dressings, where it provides analgesic and anti-inflammatory activities (Hume, 1939; Hume, 1984).
Our data show that EUG reversibly inhibits voltage-activated sodium currents (INa) in a dose-dependent manner without altering membrane resistance. Next, we investigated a possible state-dependent inhibition of VGSCs by EUG. To this end, we investigated the inhibition of the total INa and of tetrodotoxin-resistant (TTX-R) INa separately. TTX is one of the few partially specific modulators of VGSCs since it intensively blocks seven of the nine different VGSCs. Therefore, TTX provides a strategy to study different VGSCs when they are expressed in somatic cells. We studied total or TTX-R INa only, in response to EUG, and we compared the results with the inhibition and modulation induced by lidocaine (LID), a classical local anesthetic that inhibits VGSCs. Our data suggest that EUG may interact with the pre-open–closed states of the channels to enhance the inhibition of total INa, in addition to the inhibition of VGSCs in their resting states, which is remarkably different from the effects induced by LID.
Additionally, we looked into a possible inhibiting interaction between EUG and LID, when inhibiting VGSCs, in an attempt to speculate about the unknown EUG VGSC-inhibiting mechanisms. The inhibitory mechanisms of EUG and LID are similar but not the same, suggesting hypothetically different inhibition mechanisms on different VGSCs that are impossible to determine with studies like the one presented here. The currents we used here to study these drugs are mediated by many types of VGSCs that are expressed in DRG neurons.
We propose EUG as a new structural scaffold for the development of new target-specific VGSC modulators. The new suggestive inhibiting mechanism, based on the notion of the interaction between EUG and the pre-open–closed states of VGSCs that we have possibly identified in our studies as attributable to EUG, has not been described before for a cyclic monoterpenoid such as EUG, and it is currently a matter of further studies by our group.
Results
EUG quickly inhibits INa, and the inhibition is fully reversible
DRG neurons express a variety of VGSCs. Despite lacking a specific population of channels for new inhibitor studies, these cells provide an outstanding platform to study VGSCs natively expressed in their physiological environment. Additional benefits of studying VGSCs in DRG neurons are the practicality of cell preparation, the cell sizes ranging from 20 to 50 μm, and the absence of dendrites that may result in space clamp issues when these structures are present. We used voltage-gated K+ channel blockers, voltage gated-Ca2+ channel blockers, and online leak subtraction to isolate the voltage-activated INa for our experiments (see Methods).
We tested the fast inhibition by EUG on INa expressed in DRG neurons and compared the results with the inhibition produced by LID, a local anesthetic used in clinics that inhibits VGSCs to produce anesthesia. In a separate batch of experiments, we used 300 nM of tetrodotoxin (TTX) during the whole experiment to inhibit only the TTX-sensitive part of the total INa in DRG neurons (Roy and Narahashi, 1992; Elliott and Elliott, 1993; Ogata and Tatebayashi, 1993; Tan et al., 2014). This procedure enables the recording of only the TTX-R fraction of INa that is important in the peripheral nervous system as it relates to sensorial information such as the ones produced by noxious stimuli (Djouhri et al., 2003; Gudes et al., 2015; Bennett et al., 2019).
EUG and LID quickly inhibit the total INa and TTX-R INa in a concentration-dependent manner, and the inhibition is fully reversible after drug washout (Figure 1A, B). By using an in-house fast single-cell superfusion system, we show that EUG inhibits INa in less than 5 s upon its addition to the experiments. The inhibition may be even faster than 5 s, but we used a time series of 0.2 Hz to avoid inhibition accumulation during faster series. In addition, we used a holding potential of −110 mV in all our experiments in an effort to remove VGSCs from their inactivated states, to maximize the currents we recorded.
Figure 1. EUG inhibits the total INa and TTX-R INa in a dose-dependent and fully reversible manner. (A) Representative traces showing the inhibition of total INa and only TTX-R INa activated at +20-mV membrane potential, from a holding potential of −110 mV, by EUG and LID for comparison. (B) Fast inhibition and full recovery by EUG or LID on the peak of total INa or only TTX-R INa, followed by full recovery upon drug washout in a depolarization time series of 0.2 Hz. (C) Averaged fractional currents (symbols) were plotted for dose–response curves. Vertical bars denote SEM (n > 6 per point). Continuous lines are the best fit using Eq. 1 (see Methods for details). Filled symbols are from the total INa, and empty symbols are from TTX-R INa. The datasets plotted in red relate to EUG, and those in blue relate to LID. EUG inhibits total INa and TTX-R INa with IC50 of 2.27 ± 0.07 mM (n = 30) and 2.21 ± 0.08 mM (n = 32), respectively. For comparison, LID inhibits the total INa and TTX-R INa with IC50 values of 1.42 ± 0.17 mM (n = 24) and 0.44 ± 0.09 mM (n = 26), respectively. (D) Paired data show that 3 mM EUG does not significantly alter the membrane resistance of the cells. The mean of differences is shown, and it is not significantly different from zero (two-tailed paired t-test, n = 15; ns: not significant). (E) Typical action potential waveform recorded from a random DRG neuron upon a 1-ms 1-nA depolarizing current injection. EUG at a concentration of 1 mM promptly inhibits action potential firing, and the effect is quickly removed upon drug washout.
EUG was solubilized with ethanol and then diluted in experimental (external) solution to the desired concentration. The maximal final concentration of ethanol was 0.73% (vol/vol), and this amount was used to solubilize EUG for 10 and 30 mM working solutions only. For lower final concentrations of EUG, we used a molar solubilization ratio of 1:21 (EUG:ethanol). At 0.73% (vol/vol), ethanol reversibly blocks 10% of the total INa peak (Figure 2). For reference, the concentration of ethanol used to solubilize EUG at 2 mM was 0.2%.
Figure 2. Effect of ethanol 0.73% (vol/vol) on the total INa peaks. (A) Representative traces showing the inhibition of the current peaks (inset) by ethanol 0.73% (green traces). Recovered currents are also shown in blue. (B) Average ± SEM (n = 10) INa peak inhibition by ethanol 0.73% is shown as black lines, and individual values are shown as gray circles.
EUG inhibits the total INa with an IC50 value of 2.27 ± 0.07 mM (n = 30) and TTX-R INa with an IC50 value of 2.21 ± 0.08 mM (n = 32). Using an identical approach, we found that LID inhibited the total INa with an IC50 value of 1.42 ± 0.17 mM (n = 24) and TTX-R INa with an IC50 value of 0.44 ± 0.09 mM (n = 26) (Figure 1C). All IC50 values for the inhibition of INa and TTX-R INa by EUG and LID are shown in Table 1. It is noteworthy to mention that none of the EUG concentrations utilized in the present study affected the membrane resistance of the neurons. As an example, we show paired data with absolute membrane resistance values before and after 3 mM EUG is added to the experiment (Figure 1D). For these experiments, we used a hyperpolarizing pulse to −130 mV from the usual −110 mV holding potential. In addition, we tested EUG on action potential firing in a representative neuron. As predicted, 1 mM EUG prevents neurons from firing action potentials, a neuronal capability fully recovered a few seconds after EUG washout (Figure 1E).
Table 1. Dose–response parameters (IC50 and Hill coefficients) for the inhibition of total and TTX-R INa by EUG, LID, and combinations as indicated.
Remarkably, the inhibition profile of the total INa by EUG differs from that of LID in the slope of the dose-dependent inhibition curves, the Hill coefficient (EUG slope = 1.95 ± 0.120 vs. LID slope = 1.02 ± 0.169; extra sum-of-squares F-test p < 0.0001), suggesting that at least part of the inhibiting mechanism is different between these two drugs.
EUG inhibits total INa and TTX-R INa at all activating membrane potentials
From the holding potential of −110 mV, we applied depolarizing voltage steps from −90 mV to +50 mV in increments of +5 mV every 5 s (0.2 Hz). EUG at a concentration of 2 mM inhibits the total INa and TTX-R INa, as activated by all voltages without significantly affecting the reversal potential of the currents. Comparable results were found when LID at a concentration of 1 mM was tested as an inhibitor using identical procedures (Figure 3).
Figure 3. EUG inhibits INa at any depolarized membrane potential stimulus. (A) Typical families of the total INa under control conditions (absence of EUG) and in the presence of EUG at 2 mM. (B–E) Current-to-voltage (I-V) relationships for the peaks of the total INa, and only TTX-R INa are shown in the absence and presence of EUG at a concentration of 2 mM or LID at a concentration of 1 mM, as indicated. Average I-V relationships are shown as symbols, and vertical bars denote SEM (n > 9). Filled symbols are from the total INa, and empty symbols are from TTX-R INa, as indicated. The datasets plotted in red relate to EUG, and those in blue relate to LID.
Interestingly enough, the persistent total INa after 50-ms depolarization was remarkably inhibited by EUG but not by LID (Figures 4A–D). Nevertheless, neither EUG nor LID changed the kinetics of the inactivation process of total INa or TTX-R INa during a test pulse (Figures 4E–H), as discussed later.
Figure 4. EUG inhibits persistent total INa without affecting the time constant of INa inactivation. (A,B) Typical total INa traces showing the last 30 ms of 50-ms depolarizing pulses before (control) and after EUG or LID is added to the experiment at the indicated concentrations. Note the persistent total INa at the end of the depolarizing pulses right before the tail currents. (C,D) Typical I-V relationships of persistent total INa recorded at the end of 50-ms depolarizing pulses. EUG effectively inhibits persistent total INa (two-way ANOVA and Šídák’s multiple comparisons test; *: p < 0.05). Note that despite the consistent inhibition of the total INa peaks (Figure 3C), LID fails to inhibit persistent Na+ currents. Average persistent INa I-V relationships are shown as symbols, and vertical bars denote SEM (n > 9). (E–H) Total INa and TTX-R INa inactivation time constants kinetics are unaffected by EUG or LID. Inactivation time constants were taken by a double exponential fit (Eq. 2) of the decaying phase of the INa during a wide range of different depolarizing pulses, as indicated in the graphs. Average weighted inactivation time constants (Eq. 3) at different membrane potentials are shown as symbols, and vertical bars denote SEM (n > 9). Filled symbols are from the total INa, and empty symbols are from TTX-R INa, as indicated. The datasets plotted in red relate to EUG, and those in blue relate to LID.
EUG shifts the voltage dependence of Na+ conductance activation to more depolarized membrane potentials
We transformed INa peak values from each separate cell into Na+ conductance by using Ohm’s law (Eq. 4; see Methods) before normalizing all values by their maxima. We then plotted these cell-specific data against the activating membrane potential to produce conductance–voltage (G-V) curves in the absence and presence of 2 mM EUG or 1 mM LID for comparison. We also averaged data to highlight voltage-dependence shifts when they existed (Figures 5A, D, G, J). To the individual cell normalized Na+ conductance curves, we fitted Boltzmann’s equation (inactivation kinetics; Eq. 5; see Methods) for the individual voltage dependence of the activation process (V1/2-act) and the voltage sensitivity of the process (max slope-act). On average, EUG significantly shifted V1/2-act of the total INa but not of the TTX-R INa to more depolarized membrane potentials, and it decreased the max slope-act of both total INa and TTX-R INa (Figures 5A–F). Data from the inhibition of INa by LID show no changes in V1/2-act of both total INa and TTX-R INa and a change in the max slope-act of TTX-R INa only (Figures 5G–L). Figures 5A, D, G, J show the average curves that hide the individual variations in V1/2-act and max slope-act. Averaged fitting parameters and statistical analysis are shown in Table 2.
Figure 5. EUG changes the voltage dependence of total INa activation. (A–C) Normalized and averaged Na+ conductance–voltage (G-V) curves in the absence (control) and presence of EUG at a concentration of 2 mM. Data from individual cells were plotted as individual G-V curves and fitted with Eq. 5 (See Methods) for voltage dependence of the INa activation (V1/2-act) (B) and voltage sensitivity (Max slope-act) (C) parameters, before and after EUG is added. Averaged G-V curves are shown as symbols, with vertical bars denoting SEM (n > 9). (D–F) EUG did not shift the voltage dependence of TTX-R INa (E) but decreased its voltage sensitivity (F). (G–L) LID also did not shift the voltage dependence of either total (H) or TTX-R INa (K), but it decreased the voltage sensitivity of both currents [(I,L), respectively]. Filled symbols are from total INa, and empty symbols are from TTX-R INa, as indicated. The datasets plotted in red relate to EUG, and those in blue relate to LID. (B,C,E,F,H,I,K,L) Note that the statistically significant difference, at different levels, between conditions before and after EUG or LID is added, is indicated by asterisks (two-tailed paired t-test; ns: not significant; *p < 0.05; **p < 0.01; and ****p < 0.0001).
Table 2. Voltage activation and inactivation of total and TTX-R INa in the absence and presence of EUG or LID.
EUG shifts the voltage dependence of total INa and TTX-R INa inactivation to more negative membrane potentials
In addition to activation, one of the greatest properties of VGSCs is their fast inactivation process. This process limits the action of a VGSC to a few milliseconds, thus preventing long depolarizing periods that would impair the right functioning of cells in most cases. We tested EUG on the inactivation process of INa by studying the level of inactivation of the total INa and TTX-R INa after many depolarized (−140 mV–0 mV at every 5 mV) 100-ms conditioning pre-pulses. The inactivation level was tested by a pulse to +20 mV right after the conditioning pulse (Figure 6A). Graphically, we expressed the normalized and average levels of non-inactivated INa against the conditioning pre-pulse voltage to build steady-state inactivation curves (Figures 6B, E, H, K). Individual inactivation curves were fitted with Boltzmann’s equation (Eq. 6), similar to what was done with G-V curves for their voltage dependence (V1/2-inact) and voltage sensitivity (max slope-inact). EUG enhances the voltage-dependent inactivation of the total INa and TTX-R INa by shifting their voltage dependence to more negative voltages. The shift in the curves from the total INa was highly significant, nearly −15 mV on average. The shift in the curves from TTX-R INa was also significant and −5 mV on average. LID also shifted the inactivation curves of both the total INa and TTX-R INa to more negative potentials but not to the extent that EUG did, with −5-mV and −4-mV shifts on average, respectively. Neither EUG nor LID effectively changed the voltage sensitivity of the fast inactivation process of INa (Figure 6). Averaged fitting parameters and statistical analysis results for INa inactivation analysis are shown in Table 2.
Figure 6. EUG changes the voltage dependence of INa inactivation. (A) Typical families of total INa recorded at +20 mV after 100-ms conditioning pre-pulses at voltages ranging from −140 to 0 mV (Vmp-p, upper panel) under control conditions (absence of EUG, left panel) and in the presence of EUG at 2 mM (right panel). In both families of currents, the red trace relates to a Vmp-p value of −50 mV. (B,E,H,K) Normalized and averaged available total INa or TTX-R INa in the absence (control) and presence of EUG or LID at indicated concentrations were plotted against Vmp-p for inactivation curves. Averaged values are shown as symbols, and vertical bars denote SEM (n > 9). Data from individual cells were plotted as individual curves and fitted with Eq. 6 (see Methods) for V1/2-inact (the voltage dependence of the inactivation process) and Max slope-inact (the voltage sensitivity of the process). EUG shifts the inactivation curves of both the total INa (C) and TTX-R INa (F) to more negative voltages without altering their voltage sensitivity (D,G). (H–M) LID shifted the voltage dependence of both the total (I) and TTX-R INa (L) without changing the voltage sensitivity (J,M). All filled symbols relate to the total INa, and the empty symbols relate to TTX-R INa, as indicated. The datasets plotted in red relate to EUG, and those in blue relate to LID. (C,D,F,G,I,J,L,M) Note that the statistically significant difference, at different levels, between conditions before and after EUG or LID is added, is indicated by asterisks (two-tailed paired t-test; ns: not significant; *p < 0.05; **p < 0.01; and ****p < 0.0001).
EUG affects the recovery from the inactivation of total INa
We evaluated the speed of the recovery from the fast inactivation of the total INa and TTX-R INa in DRG neurons. To do so, we used the classic three-pulse voltage-clamp protocol, which consists of a pulse P1 to +20 mV lasting 50 ms to fully fast inactivate INa, a duration-varying pulse P2 at a holding potential of −110 mV to recover the currents from inactivation, and a pulse P3 to +20 mV once again to activate the fraction of the INa recovered during P2. Typical current recordings show small differences in the kinetics of recovery from inactivation before and after EUG is added to the experiment (Figure 7A). All INa peaks recorded during P3 were normalized by the corresponding (same sweep) INa peak during P1, averaged, and plotted against the duration of P2 (Figures 7B, F, J, N). Individual plots from the same cell before and after EUG was added to the experiment were fitted with two exponentials for a fast and slow kinetic component and a proportion amplitude component that we chose to represent as a percent of the fast component (Eq. 7). The statistical analyses of these parameters show that in the experiments with the total INa, but not with TTX-R INa, the percentage of fast component of the recovery from inactivation was significantly decreased from 80% to 74% (Figure 7C). In this same experiment, the slow recovery component of the recovery from inactivation was significantly delayed by EUG (Figures 7B–E). Because most of the recovery of the inactivation process happens at a fast pace (fast component), the delayed slow component did not strongly affect the overall recovery from inactivation. EUG did not affect any parameter of recovery from the inactivation of the TTX-R INa (Figures 7F–I). For comparison, we ran the same experiments with LID, which is known to delay INa recovery from inactivation. LID substantially and significantly delays the recovery from the inactivation of both total INa and TTX-R INa, mainly by affecting the fast component of the process. In the total INa, LID significantly decreases the percent fast component, and it delays both the fast and slow components of the process (Figures 7J–Q). In TTX-R INa, LID acts by significantly delaying the fast component of the process only. Averaged fitting parameters and statistical analysis results for INa recovery from the inactivation analysis are shown in Table 3.
Figure 7. EUG affects INa recovery from inactivation. (A) Typical families of the total INa recorded with the classic three-pulse voltage-clamp protocol with P1, P2, and P3 are shown in the upper panel. INa are activated by P1 and P3 under control conditions (absence of EUG, left panel) and in the presence of EUG at 2 mM (right panel) (see text for details). (B,F,J,N) Normalized and averaged P3/P1 current ratios at different P2 durations in the absence (control) and presence of EUG or LID at indicated concentrations are plotted for recovery curves. Data from individual cells were plotted as individual curves and fitted with Eq. 7 for respective percent of fast component and kinetics of fast and slow components (time constants τfast and τslow) (see Methods). (C, D, E) Individual cell parameters in the absence and presence of EUG (before and after, respectively) are shown as percent fast component and fast and slow recovery time constants. (F–I) EUG did not affect any recovery from the inactivation parameters of TTX-R INa. Recovery from inactivation data is shown as symbols, and vertical bars denote SEM (n > 9). (J–Q) For comparison, similar experiments were performed using LID as an INa inhibitor. As expected, LID delays recovery from inactivation in both total and TTX-R INa. (B,F,J,N) Filled symbols are from the total INa, and empty symbols are from TTX-R INa, as indicated. The datasets plotted in red relate to EUG, and those in blue relate to LID. (C–E, G–I, K–M, O–Q) Note that the statistically significant difference, at different levels, between conditions before and after EUG or LID is added, is indicated by asterisks (two-tailed paired t-test; ns: not significant; *p < 0.05; **p < 0.01; and **p < 0.001).
Table 3. Recovery from the inactivation of total and TTX-R INa in the absence and presence of EUG or LID.
EUG-inhibiting effects are not potentiated at high frequencies of depolarization
We tested the inhibition intensity of the total INa and TTX-R INa by EUG in a time series of depolarizations at 2 and 5 Hz (Figure 8). This test shows the inability of the VGSCs to recover from inactivation in a period between two depolarization events. The voltage clamp protocols used here consisted of 50-ms depolarizing pulses; therefore, they were characterized by a duty cycle of 10% when stimulation was at 2 Hz and 25% when stimulation was at 5 Hz. In other words, at 2 Hz, the INa represented the recovered currents during a 450-ms period, and at 5 Hz, the recovered currents during a 150-ms period. It is noteworthy that the first pulse in the series was taken after at least 1 min at a holding potential of −110 mV. Typically, 2 mM EUG does not have its fractional inhibition of total INa significantly increased at 2 Hz, from 0.53 ± 0.049 in the first pulse of the series to 0.55 ± 0.061 in the 20th pulse of the series. At a stimulation of 5 Hz, the fractional inhibition of the total INa was significantly increased from 0.49 ± 0.056 in the first pulse to 0.55 ± 0.062 in the 20th pulse (two-tailed paired t-tests, p = 0.20 for 2 Hz and p < 0.01 for 5 Hz). When EUG was used to inhibit TTX-R INa, the fractional inhibition was significantly increased from 0.37 ± 0.051 to 0.44 ± 0.064 at 2 Hz and from 0.37 ± 0.046 to 0.49 ± 0.061 at 5 Hz (two-tailed paired t-tests, p < 0.01 for 2 Hz and 5 Hz). For comparison, we also tested LID under the same conditions. At a stimulation of 2 Hz, LID performs similar to EUG both on total and TTX-R INa. With LID as an inhibitor, the total INa inhibition was significantly increased from 0.46 ± 0.034 to 0.54 ± 0.044 at a stimulation of 2 Hz and from 0.47 ± 0.021 to 0.81 ± 0.055 at 5 Hz (two-tailed paired t-tests, p < 0.01 for 2 Hz and 5 Hz). On TTX-R INa, inhibition by LID was significantly potentiated from 0.40 ± 0.025 to 0.55 ± 0.040 at 2 Hz and from 0.37 ± 0.027 to 0.64 ± 0.054 at 5 Hz (two-tailed paired t-tests, p < 0.01 for 2 Hz and p < 0.001 for 5 Hz). Fractional inhibition data and statistical analysis results for the INa inhibition series at 2 Hz and 5 Hz are shown in Table 4.
Figure 8. INa inhibition induced by EUG is weakly potentiated at high-frequency depolarizations. (A) Inhibition of total INa induced by EUG or LID over 20 successive depolarizations to +20 mV at 2 Hz and (B) at 5 Hz. It is noteworthy that the first pulse was taken after at least 1 min at a holding potential of −110 mV. (C, D) Inhibition of TTX-R INa induced by EUG or LID in a similar protocol as in (A,B). At 2-Hz stimulation but not at 5-Hz stimulation, EUG performs similar to LID both on total and TTX-R INa. Filled symbols are from total INa, and empty symbols are from TTX-R INa as indicated. At 5-Hz stimulation, LID performs better than EUG on both the total and TTX-R INa. The datasets plotted in red relate to EUG, and those in blue relate to LID. Note that a statistically significant difference between pulses 1 and 20 in each plot after two-tailed paired t-tests is indicated by asterisks. ns: not significant; **:p < 0.01 and ***:p < 0.001.
Table 4. Fractional inhibition of total and TTX-R INa by EUG or LID at different simulating frequencies.
EUG does not increase the rate of slow inactivation entry
VGSCs continue to undergo conformational changes after fast inactivation is completed. After long depolarizations, the channels are in the commonly called slow inactivated states (Vilin and Ruben, 2001). Because of the fast inactivation, slow inactivation cannot be observed directly as current decay like the fast inactivation. We studied slow inactivation with a three-pulse protocol consisting of a duration-varying P1 pulse to +20 mV intended to activate and inactivate VGSCs, followed by a 20-ms period at holding P2 to recover INa from fast inactivation, and finally, a P3 pulse to +20 mV once again that is intended to activate the available INa (Figure 9). We plotted averaged P3/P1 current ratios against the duration of P1 to estimate the kinetics of inactivation entry. P1 pulses up to 50 ms in duration mostly measure fast inactivation entry. Longer depolarizations drive the channels into slow inactivated states. We used two exponentials (Eq. 8) to fit the data to extract the level of inactivation when P1 = 2 ms, as well as the percent of the fast component and the time constants of the P3/P1 current ratio decay that we use to infer about the slow-inactivation entry kinetics. EUG does not significantly affect any of these two parameters from the total INa or TTX-R INa. In turn, our data using LID as an inhibitor show it significantly affects both the level of inactivation when P1 = 2 ms and the fast and the slow inactivating components of the slow inactivation process. All parameters from the analysis above are summarized in Table 5.
Figure 9. EUG does not affect the slow inactivation state entry. A modified three-pulse voltage-clamp protocol (upper panel) was used to estimate the amount of channels inactivated after a time-varying activating/inactivating P1 pulse to +20 mV. A fixed 20-ms P2 pulse at holding followed, and finally, a P3 pulse to +20 mV again serves to estimate the amount of slow inactivation that occurred during P1. (A,C) Total INa and TTX-R INa represented as the P3/P1 current fraction were plotted against P1 duration in the absence of EUG (control condition) and in the presence of EUG or LID at indicated concentrations. (B,D) Data for inhibition with LID are also shown for comparison. EUG does not change the kinetics of the slow inactivation process, and, as expected, LID greatly accelerates slow inactivation kinetics. Filled symbols are from total INa, and empty symbols are from TTX-R INa, as indicated. The data sets plotted in red relate to EUG, and those in blue relate to LID. All plots were fitted with double exponentials (continuous lines; Eq. 8), and the fitting parameters and statistical analysis are shown in Table 5.
Table 5. Amplitudes and kinetic parameters of the slow inactivation entry of total INa and TTX-R INa in the absence and presence of EUG or LID.
EUG may interact with pre-open states of VGSCs
Our data showed that EUG significantly and consistently shifts the total INa activation curves to more depolarized potentials (Figure 5). In addition, we showed that EUG shifts the steady-state inactivation curve to more negative potentials (Figure 6), although without inducing remarkable changes in the recovery from fast inactivation (Figure 7) and on the use-dependent inhibition that is a hallmark for local anesthetics like LID (Figure 8). One of our hypotheses that would explain the abovementioned data is that EUG interacts with and stabilizes pre-open–closed states of VGSCs (Sunami and Hiraoka, 1996). These pre-open–closed states are short-lived non-conducting states of the channels that are, therefore, not trivial to study (Armstrong, 2006; Lenaeus et al., 2017; Xiao et al., 2021; Gamal El-Din and Lenaeus, 2022). We designed a protocol, which, to the best of our knowledge, has not been used before exactly as is to study the possible interaction between EUG and the pre-open states of VGSCs. The protocol consists of depolarizing the membrane of the neurons from a −110-mV holding potential to −70 mV for various different durations (0–420 m) in a time series. The −70 mV was chosen based on the steady-state inactivation curves as the membrane potential that does not lead to current activation (Figure 5) and that, however, produces little inactivation (Figure 6) in the absence of drugs. Upon depolarizing the membrane potential to −70 mV, from our standard holding potential of −110 mV, we quickly shift the equilibrium of VGSCs to pre-open–closed states (Armstrong, 2006). The equilibrium quickly favors pre-open–closed states over the resting states upon stepping the voltage to −70 mV since the voltage sensors of VGSCs move within a few milliseconds at −70 mV (Lacroix et al., 2013). In addition, a much slower process takes channels from those pre-open–closed states to an inactive state of the VGSCs (Figure 10) (Armstrong, 2006). Our rationale was that if EUG interacts with such pre-open–closed states, the available current to be activated by P2 in our custom protocol would decrease in amplitude with the extended duration of P1. We plotted the normalized average inhibition intensity of the total INa recorded with P2, induced by EUG, against P1 duration (Figure 10). For comparison, we repeated the same experiment and data analysis but using LID as the total INa inhibitor. Our data show that EUG, but not LID, inhibited INa after a conditioning pulse (P1) to −70 mV shorter than 100 ms in an exponential time course with a time constant of approximately 72 ms. For comparison, LID induces similar effects in decreasing the current at pulse P2 of this protocol, with kinetics of 709 ms. The kinetics of EUG inhibition potentiation by P1 was statistically different from that of LID (Student’s t-test, p = 0.04). In fact, at −70 mV, some level of inactivation exists (Figure 10A), but we cannot rule out the notion that EUG may bind to pre-open–closed states, and this effect would be reflected in the steady-state inactivation curves as well. This matter is currently being worked out in terms of a kinetic model to further strengthen this notion in a separate paper. There is currently limited information about drug binding to pre-open–closed states of mammalian VGSCs (Gilliam et al., 1989; Chernoff, 1990; Carmeliet and Mubagwa, 1998; Nesterenko et al., 2011), and this is one of our lines of investigation moving forward.
Figure 10. EUG may interact with pre-open–closed states of VGSCs. A time-varying conditioning pulse P1 to −70 mV was applied before a INa activating pulse P2 to +20 mV was applied to fully activate available currents, as shown in the upper panel. (A) In the absence of drugs, the −70-mV conditioning period (P1) decreases the available current activated by P2 in a time-dependent manner. The control (drug-free) data in both cases are not statistically different (multiple Student’s t-tests, p > 0.05). (B) Normalized P2 current in the presence of the drug (EUG or LID) relative to control (drug-free) values was plotted against the P1 duration. The inhibiting effect of EUG during P2 was rapidly increased as P1 lasted longer. We fitted a single exponential decay to individual cell data, and the average time constant of the inhibiting effect increase was 72.80 ± 14.926 ms (n = 7). For comparison, LID was used instead as an inhibitor, and, despite the inhibition strength increasing with the duration of P1, the time course of the process is much slower, and it has a time constant of 709.10 ± 307.018 ms (n = 6). The kinetics of EUG inhibition potentiation by P1 was statistically different from that of LID (Student’s t-test, p = 0.04).
The inhibitory mechanisms of EUG and LID on VGSCs are similar but not the same
Our data suggest that EUG and LID may inhibit INa with a different mechanism. To further corroborate this notion, we built additional INa-inhibiting dose–response curves by EUG and LID but in the presence of LID and EUG, respectively. First, we established a pre-inhibition of 25% of total INa with LID at 0.5 mM. This was planned to leave the other 75% of the current amplitude for the building of an inhibiting dose–response curve using EUG. All currents inhibited in these experiments were fully recovered upon drug washout. A ∼25% pre-inhibition of the total INa with LID was associated with a statistically significant change in the IC50 value for the inhibition by EUG. In the presence of LID at 0.5 mM, EUG inhibits the total INa with an IC50 value of 1.84 ± 0.274 mM (vs. 2.27 ± 0.075 mM in the absence of LID; extra sum-of-squares F-test, p < 0.01). In addition, the Hill coefficient for the inhibition of INa by EUG was reduced to 0.98 ± 0.140 mM (n = 12) by pre-inhibiting the total INa with LID (vs. 1.95 ± 0.120 with EUG alone; extra sum-of-squares F-test, p < 0.0001; Table 1; Figures 11A, B).
Figure 11. EUG and LID effects on VGSCs are interacting. (A) Time course of the double inhibition of the total INa by LID and EUG, as indicated. (B) After an approximate 25% INa pre-inhibition by LID at 0.5 mM, a dose–response curve for the inhibition by EUG was built. The IC50 value for EUG when inhibiting the total INa was 2.27 ± 0.07 mM. After a 25% pre-blockade by LID, EUG inhibits the total INa with a significantly lower IC50 of 1.84 ± 0.275 mM (extra sum-of-squares F-test, p < 0.01). The Hill coefficient of the EUG inhibition was also significantly reduced from 1.95 ± 0.120 without LID to 0.98 ± 0.140 after pre-inhibition with LID (extra sum-of-squares F-test, p < 0.0001). (C) Time course of the double inhibition of total INa by LID and EUG as indicated. (D) Similarly, after an approximate 25% INa pre-inhibition by EUG at 1.3 mM, a dose–response curve for the inhibition by LID was built. The IC50 value for LID when inhibiting the total INa was 1.42 ± 0.169 mM. After a 25% pre-inhibition by EUG, LID inhibits the total INa with a significantly lower IC50 value of 0.73 ± 0.126 mM (extra sum-of-squares F-test, p < 0.0001), with the Hill coefficient unchanged. Fractional inhibition data plotted against inhibitor concentrations were fitted using the Hill equation (Eq. 1, see text for details). Continuous lines plotted on b and d are best fits of the Hill equation to the color-coded data.
Next, we carried out the opposite experiment with regard to inhibitors. First, we pre-inhibited around 25% of the total INa with 1.3 mM EUG and proceeded with a dose–response curve for the inhibition of the remaining 75% of the currents with LID. This time, the remaining 75% of INa was inhibited by LID with a significantly lower IC50 value of 0.73 ± 0.126 mM (n = 10) vs. 1.42 ± 0.169 with LID alone (extra sum-of-squares F-test, p < 0.0001) and with no significant changes in the Hill coefficient compared with LID alone (Hill coefficient of 1) (Table 1; Figures 11C, D).
Discussion
Here, we show, for the first time, a fast and fully reversible inhibition of VGSCs by EUG. The inhibition of INa (total and TTX-R) by EUG includes an interaction with the resting states of VGSCs with an affinity that is comparable but statistically different from the more potent inhibition produced by LID on the same currents (Figure 1; Table 1; dose–response curves compared by two-way ANOVA, p < 0.0001). Our data also suggest that EUG, different from LID, may interact with the pre-open–closed states of the VGSCs. Although speculative, this assumption is our best hypothesis to explain the hallmarks of the EUG inhibition of INa: i) INa activation curves shift to more depolarized potentials (Figure 5); ii) steady-state inactivation curves shift to more negative potentials (Figure 6); iii) light changes in the recovery from fast inactivation (Figure 7); iv) light use-dependent inhibition of INa (Figure 8); and most importantly, v) the time-dependent enhanced inhibition of INa by a −70-mV membrane potential conditioning period (Figure 10).
EUG (1-allyl-4-hydroxy-3-methoxybenzene) is a naturally occurring compound present in the essential oil of many aromatic plants like clove, wormwood, and sweet basil. EUG has been used for a long time as a medical agent, food preservative, and flavoring agent. The excellent therapeutic index of EUG has been recently addressed in a comprehensive review that points out EUG as a generally acknowledged as safe (GRAS) chemical by the World Health Organization (WHO) (Tavvabi-Kashani et al., 2024). EUG is the ingredient responsible for analgesia in zinc eugenolate chelate, a dental cement material used in dentistry (FADM, 2001). EUG has been previously studied by many groups aiming at finding its ion channel-modulating properties, and it has been studied as a modulator of VGSCs (Park et al., 2006; Cho et al., 2008; Park et al., 2009; Moreira-Lobo et al., 2010), voltage-gated calcium channels (Sensch et al., 2000; Magyar et al., 2004; Lee et al., 2005; Chung et al., 2008; Seo et al., 2013), voltage-gated potassium channels (Erdélyi, 1999; Magyar et al., 2004; Li et al., 2007; Chung et al., 2008), GABAA channels, purinergic channels (Li et al., 2008), and TRP channels (Yang et al., 2003; Chung et al., 2014; Latorre and Baldion, 2020; Takahashi et al., 2021). These previous studies suggest EUG as a pore blocker of voltage-gated ion channels since no evidence has been found that can attest to any change in the biophysical properties of the channel voltage dependence or voltage sensitivity.
DRG neurons express several biophysically distinct VGSC isoforms, including Nav1.1, Nav1.6, Nav1.7, Nav1.8, and Nav1.9 (Berta et al., 2008; Ho and O’Leary, 2011; Shiers et al., 2020). Subunits Nav1.7, Nav1.8, and Nav1.9 are preferentially expressed in DRG neurons, and they are associated with mechanisms related to neuropathic and inflammatory pain (Chevrier et al., 2004; Maruyama et al., 2004; Liu and Wood, 2011; Bennett et al., 2019). Studies using VGSCs natively expressed in mammalian cells have the advantage of the channels being in their physiologic environment, a condition that is difficult to replicate when the channels are heterologously expressed in HEK cells, for instance. Studying VGSCs in their native membranes further increases the significance of the outcome data since only in that environment might the channels be inserted in their multiprotein complexes. Recent studies focus on ion channels and other membrane proteins as protein structures that are part of the multiprotein complexes in the membrane. Within these complexes, proteins are functionally coupled (Vacher et al., 2008; Abriel et al., 2015; Baronas et al., 2015; Subramanyam and Colecraft, 2015; Heijman and Dobrev, 2018). In the present study, we took advantage of those aspects in exchange for a more precise pharmacological study that a heterologous system of expression would provide in terms of studying an INa that is mediated by a single VGSC subtype.
The affinity of INa inhibition by EUG
Previous studies have reported the inhibiting properties of EUG on VGSCs (Park et al., 2006; Cho et al., 2008; Park et al., 2009; Moreira-Lobo et al., 2010). In one of these studies, EUG was tested on the total INa and TTX-R INa, both activated from a holding membrane potential of −80 mV (Cho et al., 2008). Their recordings show a considerable rundown of currents during the experiments. In addition, the inhibition of INa induced by EUG shown in that study was not fully recovered to its initial values. Finally, that study exposed the neurons to EUG for 10 min, which, together with the drug application technique, could have caused INa rundown during the course of the experiments. We propose that these four drawbacks contributed to a possible overestimation of EUG affinity to block INa in the study conducted by Cho et al. (2008).
Our dose–response curves were built with the fast inhibition of INa that was activated from a holding potential of −110 mV chosen to maximize VGSC recovery from inactivation, therefore INa maximization (Figures 1B–C). Our data demonstrate an immediate inhibition of the total INa or TTX-R INa by EUG or LID upon inhibitor perfusion onset, with 100% recovery upon inhibitor washout. Note that the INa amplitude of the currents we studied is stable throughout the time series, and they are rundown-free (Figure 1B). Our data shown in Figure 3 also demonstrate a rundown-free recording both under control conditions and in the presence of EUG or LID as INa inhibitors: when at maximal Na+ conductance activation at test potentials of +20 mV or more, all I-V curves are linear. Our approach was intended to estimate the affinity of the fast inhibition of INa by EUG. Our experiments using LID as the positive control serve as a validation for our method.
Our data show that EUG inhibits the total INa in DRG neurons with an IC50 value of 2.27 mM and a Hill coefficient of 1.95. EUG inhibits TTX-R INa with virtually the same properties, with an IC50 value of 2.21 mM and a Hill coefficient of 1.84 (Figure 1C; Table 1). The affinity of the binding of EUG to the resting states of VGSCs (by using a holding potential of −110 mV and a test pulse to +20 mV) was statistically different from the inhibition induced by LID, with an IC50 value of 1.42 mM, and TTX-R INa with an IC50 value of 0.44 mM, which agrees with previous reports (Kistner et al., 2010). EUG affinity to inhibit VGSCs is far enough for local anesthesia. For comparison, a popular 2% LID formulation for human local anesthesia is a solution of LID at 85 mM. Interestingly, the Hill coefficient of the dose responses when EUG is used to inhibit INa, ∼2, is remarkably different from that when LID is used to inhibit the same currents (∼1). The Hill coefficient of dose–response curves is often related to binding sites, but the use of that knowledge is limited to orthosteric inhibitors such as pore blockers and not allosteric or state-dependent inhibitors (Bindslev, 2008; Prinz, 2010). Therefore, the only suggestion we can make, based on our data presented here and the Hill coefficients, is that EUG and LID may have different VGSC-inhibiting mechanisms. Additionally, in the case of a Hill coefficient of 2, for the inhibition of INa by EUG, we can rationalize that the coefficient does not arise from different EUG affinities to different VGSCs expressed in the neurons studied here. If that was the case, the slope of the dose–response curve would, instead, be shallower, thus showing an apparent Hill coefficient smaller than 1. Hence, we reinforce that since we deal with different VGSC isoforms here, with different kinetics and voltage dependencies, we simply speculate that the Hill coefficient we found for the INa inhibition by EUG, ∼2, indicates a different inhibiting mechanism compared to LID.
EUG and LID may interact to inhibit INa
Our data shown in Figure 11 suggest that EUG and LID may interact when blocking INa. The 25% INa pre-inhibition with either drug induces a significant improvement on the inhibition by the other, which blocks the remaining 75% of the currents with a lower IC50 value compared to the drug alone (Table 1). In addition, pre-blocking INa with LID reduces the Hill coefficient of the EUG inhibition curve from ∼2 with EUG alone to ∼1 in the presence of LID. We are extremely cautious when interpreting these data, especially because they come from experiments with a natively expressed variety of different VGSCs that could be inhibited by EUG or LID with different mechanisms. Therefore, we took a very conservative approach and limited the conclusion of these experiments to report that EUG and LID inhibit INa with likely different mechanisms. New experiments using a single VGSC subunit expressed heterologously are necessary to continue the investigation of this aspect of INa inhibition by EUG and LID simultaneously.
EUG and LID interact differently with distinct kinetic states of VGSCs
Our dose–response curves of inhibition of the total INa and TTX-R INa by EUG relate to the interaction of the inhibitor with the resting states of VGSCs. However, our detailed investigation revealed that EUG interacts with other states of VGSCs. Our data suggest that EUG may interact with VGSC pre-open–closed states possibly with a higher affinity than that of resting states (Figure 10). This notion is suggested by our experiments with a conditioning period at a membrane potential of −70 mV, which is known to populate pre-open–closed states of VGSCs (Armstrong, 2006), before INa is activated by a membrane potential step of +20 mV. Our data demonstrate that INa inhibition by EUG is intensified by such conditioning voltage steps, with a time constant of ∼70 ms, and we propose that this finding may be due to an interaction of EUG with the pre-open–closed states of VGSCs that are associated with a membrane potential of −70 mV. Specific experiments are still needed, but this piece of data alone suggests that such an interaction may take place with higher affinity than that with the resting states of the channels (when we activate the currents straight from a membrane potential of −110 mV). Our data using LID, an inhibitor known to bind to the inactivated states of VGSCs with higher affinity than to resting states, show the remarkable difference from that produced with EUG after time-varying −70-mV conditioning pulses. LID has its inhibiting effects on INa enhanced by the −70-mV conditioning pulse with a time constant of ∼700 ms. We propose that LID has its inhibiting potency increased in these experiments as more channels become inactivated, a slow process at −70-mV membrane potential (Armstrong, 2006). Differently, EUG would have its inhibiting potency increased faster (with a ∼70-ms time constant) as the pre-open–closed states, quickly achieved after stepping the membrane potential to −70 mV, are the states to which EUG would have an increased affinity. This notion will be a matter of further investigation in the near future using computation biology for kinetic models of this system.
We also propose that the right shift in the activation curves induced by EUG might be due to the interaction of EUG with pre-open–closed states of the VGSCs (Figure 5). Weak depolarizations would provide the channels with time to populate their pre-open–closed states, and that would be enough for EUG to further interact with the channels if the affinity to those states is higher than to resting states. LID does not show such a right shift in the voltage dependence of INa activation, in agreement with the notion that this drug would not interact with those pre-open–closed states. The remarkable effect of EUG on the steady-state inactivation curve of the total INa might also have been caused by the binding of EUG to pre-open–closed states. We propose the novel idea that a 100-ms conditioning pre-pulse, as shown in Figure 6A, is not a perfect measure of the binding of drugs to the fast inactivated states of VGSCs, even when we consider inactivation from pre-open–closed states. Other studies show that 1,000-ms or 10-s (10-fold or 100-fold longer than ours) conditioning pulses indeed show a remarkable shift in the V1/2-inact of the available currents, as induced by LID (Clarkson et al., 1988; Sheets et al., 2010). In particular, with a 10-s conditioning pulse, LID shifted the available curve by 10 mV only. For comparison, EUG shifted the inactivation curve by 16 mV with a 100-ms conditioning pulse. Our 100-ms conditioning pre-pulses may not last enough to inactivate the channels and for an interaction with LID after that. Once again, LID, known to interact with the inactivated states of VGSCs, does not shift the inactivation curves with the intensity or with the short conditioning pulse as EUG does (Figure 6). We suggest that this voltage clamp protocol with a 100-ms conditioning pre-pulse may also evidence the drug binding to the pre-open–closed states of the channels, and we think that EUG induces a consistent 16-mV shift in these curves because it might interact with the pre-open–closed states of the channels during the conditioning pre-pulses around −70 mV.
In our view, protocols to test the recovery from inactivation at a −110-mV membrane potential (Figure 7), high frequencies of depolarizing time series (Figure 8), and rate of slow inactivation entry (Figure 9) are better ways to study the interaction of inhibitors with the inactivated states of VGSCs (Clarkson et al., 1988). We sustain this notion because in all these protocols, the membrane potential is held at −110 mV (the holding potential) or at +20 mV when INa are activated and inactivated. Neither of these voltages is associated with the pre-open–closed states of VGSCs (Lenaeus et al., 2017). EUG underperformed in all three voltage clamp protocols compared with LID. These results show that EUG is not an inhibitor that binds to the inactivated states of VGSCs, fast and slow, with a higher affinity than binds to the resting states of VGSCs.
Our data show that EUG does not possess the ability to interact with pre-open–closed states of TTX-R VGSCs. Therefore, this EUG property must be related to the TTX-S VGSCs that are expressed in DRG neurons, i.e., Nav1.1, Nav1.6, and Nav1.7. Interestingly, since we observed the state-dependent effects of EUG on the total INa, i.e., on TTX-S plus TTX-R currents, we hypothesize that the observable effects on a TTX-S VGSC subunit would be even more clear.
EUG interacts with VGSCs, and this interaction modifies the voltage dependence of the channels. EUG may interact directly with the voltage sensors of the channels, as shown before for small molecules binding on extracellular residues of the homologous domain IV voltage sensor to cause a shift in the inactivation curve of the currents mediated by Nav1.3 (McCormack et al., 2013). EUD may also bind to the canonical local anesthetic binding site on S6 segments (the activation gate) that is allosterically connected to the voltage sensors of homologous domains III and IV (Sheets and Hanck, 2003; Sheets and Hanck, 2007). It is known that the activation of the voltage sensor of domain IV in VGSCs is necessary for the fast inactivation of these VGSCs (Kirsch et al., 1989; Campos et al., 2008). It seems that EUG does not interact with those regions since our data do not indicate stabilization of the inactivated states of VGSCs. Nevertheless, EUG may interact with pre-open–closed states of VGSCs as configured by the activation of voltage sensors of domains I and II or either one individually. Importantly, we cannot rule out that the voltage-dependent effect of EUG on the total INa reported here might be due to EUG interactions with an allosteric biding site elsewhere in the channel, including in the pore domain of the channels, a pathway we previously reported to be possible in Shaker K+ channels (Bassetto et al., 2021).
Another interesting and surprising piece of data in the present work is that EUG, but not LID, inhibits persistent total INa in DRG neurons (Figures 4A–D). It is important to mention that the inhibition of persistent INa could only be apparent if the inactivation kinetics of INa were accelerated. For example, we showed that neither EUG nor LID changed the kinetics of the inactivation process during a test pulse of the currents we studied (Figures 4E–H). Therefore, it is possible that EUG inhibits the specific channels that yield persistent currents more intensively than LID. This notion can lead to the inhibition of specific channels by EUG, and it will be investigated in our next studies. A previous study expressing the predominant cardiac VGSCs showed LID as an effective inhibitor of persistent voltage-activated Na+ currents (Dumaine and Kirsch, 1998). We propose that EUG might have a similar effect on the persistent INa observed in the present study.
Conclusion
Our data show that EUG is a VGSC inhibitor with an affinity comparable with LID. The inhibiting mechanism of EUG on the resting states of the VGSCs overlaps with that of LID. EUG modifies the voltage-dependent activation and inactivation processes of INa, and, different from LID, that effect might be due to the possible interaction of EUG with the pre-open–closed states of the VGSCs. Altogether, our data point out EUG as a state-dependent VGSC blocker different from LID, suggestively for those sensitive to TTX.
Materials and methods
General experimental procedures
All salts and drugs used in this work were purchased from Sigma (Sigma, St. Louis, MO). EUG was at least 98% pure, as stated by vendors. The pipette solution and the bath solutions were prepared directly from the salts and maintained at −20 °C until the day of the experiment, when they were thawed. A solution of eugenol (stock solution) was also prepared in advance to a concentration of 1 M in ethanol and stored at −20 °C. Just before the experiment, the stock solution was diluted in the bath solution to the desired final concentration of eugenol and sonicated for 5 min. The maximal final concentration of ethanol used in this study was 0.73% vol/vol, which, according to our own data (Figure 2) and data obtained by others (Wu and Kendig, 1998; Horishita and Harris, 2008), inhibits INa in approximately 10% and without any change in the voltage dependence of the activation and inactivation processes. This solution containing eugenol was applied over a single cell through a home-made pressurized perfusion system connected to a double-barreled capillary (Harvard Apparatus, Holliston, MA) to deliver the solutions with and without eugenol in the vicinity of the tested cell, with very fast interchangeable capability.
Cell preparation
DRG neurons from 1–3-day-old rats were used as a model of natively expressed VGSCs. The cell preparation procedure was adapted from the study by Kostyuk et al. (1981). In brief, the rats were euthanized by decapitation, the vertebral canal was opened, and DRGs were quickly removed. All animals were handled in compliance with the Guide for the Care and Use of Laboratory Animals by the US National Institutes of Health (NIH Publication 85–23, revised 1996; http://www.nap.edu/.readingroom/books/labrats/index.html). In addition, all experimental protocols were approved by the Comissão de Ética no Uso de Animais (CEUA), a collegiate body linked to the Congregation of the Institute of Biomedical Sciences of the University of São Paulo (ICB-USP). Ganglia were digested with trypsin 0.25% in a Ca2+, Mg2+-free Earle´s balanced salt solution (EBSS) containing (mM) 132.8 NaCl, 5.3 KCl, 1 NaH2PO4, 5.5 glucose, and 10 HEPES, pH 7.4. After digestion, the ganglia were mechanically reduced using a fire-polished Pasteur pipette in a Ca2+, Mg2+-free EBSS containing 5 U/mL DNAse (type I; Sigma) and 0.15% trypsin inhibitor (type IS; Sigma), supplemented with 10% fetal calf serum. After pelleting, the cells were resuspended in Dulbecco's modified Eagle's medium (DMEM; Sigma) supplemented with 10% fetal calf serum, 100 UI/mL penicillin, and 100 mg/mL streptomycin and seeded on glass coverslips previously treated with poly-L-lysine. The cells were kept in a water-jacket incubator at 37 °C and a 5% CO2 atmosphere until just before experiments, which were carried out during the first 7 days after cell isolation.
Electrophysiology
Voltage-activated Na+ currents passing through VGSCs were recorded using voltage clamping with the conventional whole-cell patch-clamp configuration (Hamill et al., 1981). Patch electrodes were fabricated from borosilicate glass capillaries using a model PB-7 micropipette puller (Narishige, Tokyo, Japan). Recording pipettes were pulled from borosilicate glass to achieve initial bath resistances averaging 2 MΩ and were filled with an intracellular solution containing (mM) 10 NaCl, 150 CsF, 10 TEA chloride, 1 ATP, 4.5 MgCl2, 9 EGTA, and 10 HEPES, pH 7.3. CsCl and TEA chloride were used in order to eliminate K+ currents. Cells were bathed during the recordings in an extracellular solution containing (mM) 82 choline chloride, 50 NaCl, 1.2 MgCl2, 1.8 CaCl2, 1 CoCl2, 4 KCl, 5 glucose, and 10 HEPES, pH 7.4. Na+ was partially replaced by choline (82 mM) to decrease the sodium driving force, avoid overload of the amplifier, and ensure a good voltage clamp. Co2+ in the extracellular solution was used to eliminate Ca2+ currents. After achieving a high-resistance seal, the whole-cell configuration was usually established by applying negative pressure to the pipette. Typical access resistance values were below 5 MΩ. Command voltage waveforms were generated in a computer using Clampex 10 software (Molecular Devices, Foster City, CA), and a DAC interface (model 1322; Molecular Devices, Foster City, CA) delivered the analog voltage signal to the Axopatch 200B patch-clamp amplifier (Molecular Devices, Sunnyvale, CA), which maintained the pipette voltage. The P/-4 protocol was used to remove voltage-independent leak currents and uncompensated capacitance (Bezanilla and Armstrong, 1977). The current recordings were low-pass filtered at 2 kHz by a built-in Bessel filter in the amplifier, sampled at 25 kHz in 16-bit digital levels by an ADC, and recorded on a hard disk for ulterior analysis. The cell membrane capacitance was canceled, and access resistance was routinely compensated (85% for both prediction and compensation; lag set to 10 ms). The holding potential in all experiments was set to −110 mV to avoid inactivation of the VGSC. All recordings were performed between 20°C and 23 °C. The recording chamber was continuously perfused with the bath solution to avoid unstirred layers. Recording cells exhibiting an INa peak of less than 500 pA throughout several depolarizing voltages were discarded. AP were recorded using the same cell preparation but in the current-clamp mode and using a pipette solution containing (mM) 10 NaCl, 150 KCl, 1 ATP, 4.5 MgCl2, 9 EGTA, and 10 HEPES, pH 7.3, and a bath solution, including the solution where eugenol was dissolved, containing (mM) 132 NaCl, 1.2 MgCl2, 1.8 CaCl2, 4 KCl, 5 glucose, and 10 HEPES, pH 7.4.
Data analysis and graphs
Scientific data were processed, analyzed, and plotted using Clampfit (Molecular Devices, Foster City, CA), GraphPad Prism (GraphPad Software, LLC, La Jolla, CA), Origin (OriginLab, Northampton, MA), and Microsoft Excel (Microsoft, Redmond, WA). All the graphs represent mean values from at least six independent experiments, unless otherwise noted. The vertical bars in the graphs are the standard error of the mean (SEM).
The fitting functions mentioned in the Results section are as follows:
Dose–response curves were fitted with the Hill formalism:
where the Normalized (uninhibited) current is the remaining INa after inhibition by EUG or LID (the Fractional current in the graphs), IC50 is the concentration of eugenol that inhibits 50% of the INa, and nH is the Hill coefficient.
INa inactivation kinetics were fitted by the following exponentials:
where INa is the Na+ current at a given moment t, A1-2 are the weights of the respective exponentials, τ1-2 are their time constants, and y0 is an adjusting factor for persistent currents. Inactivation time constants were calculated as follows:
Current–voltage relationships were transformed into Na+ conductance activation by voltage (G-V) curves by using Ohm’s law:
where Gm is the equivalent conductance at Ipeak, which, in turn, is the current peak value at the given voltage, Vm is the membrane potential, and Vr is the reversal potential for Na+.
Na+ conductance activation by voltage (G-V) curves was fitted by the following equation:
where Normalized Na+ conductance is the fractional conductance activated at a given membrane potential Vm, V1/2-act is the membrane potential for half-maximal Na+ conductance activation (the midpoint), and Max slope-act is the voltage sensitivity of the activation by voltage process.
INa inactivation by voltage curves (inactivation curves) was fitted by the following equation:
where Available Na+ current is the Na+ current after the conditioning pre-pulse voltage period Vmp-p, V1/2-inact is the Vmp-p that inactivates half of INa, and Max slope-inact is the voltage sensitivity of the inactivation by voltage process.
The kinetics of INa recovery from inactivated states were fitted with the following sum of exponentials:
where P3/P1 current ratio is the fractional INa recovered from inactivation during P2, %fast is the % fraction of the fast recovery from inactivation component amplitude, τ1 is the faster time constant of the INa recovery from inactivation, and τ2 is the slower time constant of the INa recovery from inactivation.
where P3/P1 current ratio is the fractional INa inactivated during P1 and not recovered from inactivated during P2 (see text for details), Plateau is the current level at P1>>10000 ms, YP1=2ms is the P3/P1 current ratio when P1 = 2 ms, %fast is the % fraction of the fast inactivation component, τ1 is the faster time constant of the slow inactivation component, and τ2 is the slower time constant of the slow inactivation component.
Statistical analysis
Data from individual cells were treated individually, including for fitting purposes. Pooled fitting parameters from different groups, e.g., EUG vs. control (its absence), were compared using a paired t-test to detect consistent changes in the parameters that relate to the drugs. Levels of significance were *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001.
Summary
Eugenol is an aromatic substance obtained from the essential oil of many plants that produces analgesia by a still uncertain mechanism. Here, we show comprehensive data indicating that eugenol inhibits voltage-gated Na+ channels with a mechanism that is different from lidocaine. We propose, based on the interpretation of our findings, that eugenol inhibits voltage-gated Na+ channels by interacting with their resting, pre-open–closed, and inactivated states.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Colegio Brasileiro de Experimentação Animal (COBEA). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
LM-J: formal analysis, and writing–review and editing. JL-C: funding acquisition, supervision, and writing–original draft. AC: conceptualization, data curation, funding acquisition, investigation, project administration, resources, supervision, and writing–original draft. JC-d-S: conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by Fundacao de Amparo a Pesquisa do Estado de Sao Paulo, Brazil–Governo do Estado de São Paulo, by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior–Capes, Brazil–Governo Federal, and by the Department of Anesthesiology at the University of Arizona. This work was supported by Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq), Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), Coordenaccao de Aperfeicoamento de Pessoal de Nivel Superior (CAPES), and Universidade de Sao Paulo (USP).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1354737/full#supplementary-material
References
Abriel, H., Rougier, J.-S., and Jalife, J. (2015). Ion channel macromolecular complexes in cardiomyocytes: roles in sudden cardiac death. Circ. Res. 116 (12), 1971–1988. doi:10.1161/CIRCRESAHA.116.305017
Armstrong, C. M. (2006). Na channel inactivation from open and closed states. Proc. Natl. Acad. Sci. U. S. A. 103 (47), 17991–17996. doi:10.1073/pnas.0607603103
Bagal, S. K., Marron, B. E., Owen, R. M., Storer, R. I., and Swain, N. A. (2015). Voltage gated sodium channels as drug discovery targets. Channels Austin Tex 9 (6), 360–366. doi:10.1080/19336950.2015.1079674
Baronas, V. A., McGuinness, B. R., Brigidi, G. S., Gomm, K. R. N., Vilin, Y. Y., Kim, R. Y., et al. (2015). Use-dependent activation of neuronal Kv1.2 channel complexes. J. Neurosci. Off. J. Soc. Neurosci. 35 (8), 3515–3524. doi:10.1523/JNEUROSCI.4518-13.2015
Bassetto, C. A., Carvalho-de-Souza, J. L., and Bezanilla, F. (2021). Molecular basis for functional connectivity between the voltage sensor and the selectivity filter gate in Shaker K+ channels. eLife 10, e63077. doi:10.7554/eLife.63077
Beer, A.-M., Lukanov, J., and Sagorchev, P. (2007). Effect of Thymol on the spontaneous contractile activity of the smooth muscles. Phytomedicine Int. J. Phytother. Phytopharm. 14 (1), 65–69. doi:10.1016/j.phymed.2006.11.010
Bennett, D. L., Clark, A. J., Huang, J., Waxman, S. G., and Dib-Hajj, S. D. (2019). The role of voltage-gated sodium channels in pain signaling. Physiol. Rev. 99 (2), 1079–1151. doi:10.1152/physrev.00052.2017
Berta, T., Poirot, O., Pertin, M., Ji, R.-R., Kellenberger, S., and Decosterd, I. (2008). Transcriptional and functional profiles of voltage-gated Na(+) channels in injured and non-injured DRG neurons in the SNI model of neuropathic pain. Mol. Cell Neurosci. 37 (2), 196–208. doi:10.1016/j.mcn.2007.09.007
Bezanilla, F., and Armstrong, C. (1977). Inactivation of the sodium channel. I. Sodium current experiments. J. Gen. Physiol. 70 (5), 549–566. doi:10.1085/jgp.70.5.549
Binu, P., Priya, N., Abhilash, S., Vineetha, R. C., and Nair, R. H. (2017). Studies on curative efficacy of monoterpene eugenol on anti-leukemic drug arsenic trioxide induced cardiotoxicity. Biomed. Pharmacother. Biomedecine Pharmacother. 91, 559–566. doi:10.1016/j.biopha.2017.04.087
Bouza, A. A., and Isom, L. L. (2018). Voltage-gated sodium channel β subunits and their related diseases. Handb. Exp. Pharmacol. 246, 423–450. doi:10.1007/164_2017_48
Brackenbury, W. J., and Isom, L. L. (2011). Na channel β subunits: overachievers of the ion channel family. Front. Pharmacol. 2, 53. doi:10.3389/fphar.2011.00053
Campos, F. V., Chanda, B., Beirão, P. S. L., and Bezanilla, F. (2008). Alpha-scorpion toxin impairs a conformational change that leads to fast inactivation of muscle sodium channels. J. Gen. Physiol. 132 (2), 251–263. doi:10.1085/jgp.200809995
Carmeliet, E., and Mubagwa, K. (1998). Antiarrhythmic drugs and cardiac ion channels: mechanisms of action. Prog. Biophys. Mol. Biol. 70 (1), 1–72. doi:10.1016/s0079-6107(98)00002-9
Catterall, W. A. (2012). Voltage-gated sodium channels at 60: structure, function and pathophysiology. J. Physiol. 590 (11), 2577–2589. doi:10.1113/jphysiol.2011.224204
Catterall, W. A. (2017). Forty years of sodium channels: structure, function, Pharmacology, and epilepsy. Neurochem. Res. 42 (9), 2495–2504. doi:10.1007/s11064-017-2314-9
Chernoff, D. M. (1990). Kinetic analysis of phasic inhibition of neuronal sodium currents by lidocaine and bupivacaine. Biophys. J. 58 (1), 53–68. doi:10.1016/S0006-3495(90)82353-5
Chevrier, P., Vijayaragavan, K., and Chahine, M. (2004). Differential modulation of Nav1.7 and Nav1.8 peripheral nerve sodium channels by the local anesthetic lidocaine. Br. J. Pharmacol. 142 (3), 576–584. doi:10.1038/sj.bjp.0705796
Cho, J. S., Kim, T. H., Lim, J.-M., and Song, J.-H. (2008). Effects of eugenol on Na+ currents in rat dorsal root ganglion neurons. Brain Res. 1243, 53–62. doi:10.1016/j.brainres.2008.09.030
Chung, G., Im, S. T., Kim, Y. H., Jung, S. J., Rhyu, M.-R., and Oh, S. B. (2014). Activation of transient receptor potential ankyrin 1 by eugenol. Neuroscience 261, 153–160. doi:10.1016/j.neuroscience.2013.12.047
Chung, G., Rhee, J. N., Jung, S. J., Kim, J. S., and Oh, S. B. (2008). Modulation of CaV2.3 calcium channel currents by eugenol. J. Dent. Res. 87 (2), 137–141. doi:10.1177/154405910808700201
Clarke, S. (2008). “Chapter 3 - families of compounds that occur in essential oils,” in Essential chemistry for aromatherapy Editor S. Clarke 2 (Edinburgh: Churchill Livingstone), 41–77.
Clarkson, C. W., Follmer, C. H., Ten Eick, R. E., Hondeghem, L. M., and Yeh, J. Z. (1988). Evidence for two components of sodium channel block by lidocaine in isolated cardiac myocytes. Circ. Res. 63 (5), 869–878. doi:10.1161/01.res.63.5.869
Dallmeier, K., and Carlini, E. A. (1981). Anesthetic, hypothermic, myorelaxant and anticonvulsant effects of synthetic eugenol derivatives and natural analogues. Pharmacology 22 (2), 113–127. doi:10.1159/000137479
Dallmeier, Z. K. R., Zelger, J. L., and Carlini, E. A. (1983). New anticonvulsants derived from 4-allyl-2-methoxyphenol (Eugenol): comparison with common antiepileptics in mice. Pharmacology 27 (1), 40–49. doi:10.1159/000137828
de Almeida, R. N., Agra, M. de F., Maior, F. N. S., and de Sousa, D. P. (2011). Essential oils and their constituents: anticonvulsant activity. Mol. Basel Switz. 16 (3), 2726–2742. doi:10.3390/molecules16032726
Djouhri, L., Fang, X., Okuse, K., Wood, J. N., Berry, C. M., and Lawson, S. N. (2003). The TTX-resistant sodium channel Nav1.8 (SNS/PN3): expression and correlation with membrane properties in rat nociceptive primary afferent neurons. J. Physiol. 550 (Pt 3), 739–752. doi:10.1113/jphysiol.2003.042127
dos Santos-Nascimento, T., Veras, K. M., Cruz, J. S., and Leal-Cardoso, J. H. (2015). Inhibitory effect of terpinen-4-ol on voltage-dependent potassium currents in rat small sensory neurons. J. Nat. Prod. 78 (2), 173–180. doi:10.1021/np4009249
Dumaine, R., and Kirsch, G. E. (1998). Mechanism of lidocaine block of late current in long Q-T mutant Na+ channels. Am. J. Physiol. 274 (2), H477–H487. doi:10.1152/ajpheart.1998.274.2.H477
Elliott, A. A., and Elliott, J. R. (1993). Characterization of TTX-sensitive and TTX-resistant sodium currents in small cells from adult rat dorsal root ganglia. J. Physiol. 463, 39–56. doi:10.1113/jphysiol.1993.sp019583
Erdélyi, L. (1999). Guaiacol and vanilloid compounds modulate the A-type potassium currents in molluscan neurons. Acta Biol. Hung 50 (1–3), 65–79. doi:10.1007/bf03543032
Gamal El-Din, T. M., and Lenaeus, M. J. (2022). Fenestropathy of voltage-gated sodium channels. Front. Pharmacol. 13, 842645. doi:10.3389/fphar.2022.842645
Gaudioso, C., Hao, J., Martin-Eauclaire, M.-F., Gabriac, M., and Delmas, P. (2012). Menthol pain relief through cumulative inactivation of voltage-gated sodium channels. Pain 153 (2), 473–484. doi:10.1016/j.pain.2011.11.014
George, A. L. (2005). Inherited disorders of voltage-gated sodium channels. J. Clin. Invest. 115 (8), 1990–1999. doi:10.1172/JCI25505
Gilliam, F. R., Starmer, C. F., and Grant, A. O. (1989). Blockade of rabbit atrial sodium channels by lidocaine. Characterization of continuous and frequency-dependent blocking. Circ. Res. 65 (3), 723–739. doi:10.1161/01.res.65.3.723
Gudes, S., Barkai, O., Caspi, Y., Katz, B., Lev, S., and Binshtok, A. M. (2015). The role of slow and persistent TTX-resistant sodium currents in acute tumor necrosis factor-α-mediated increase in nociceptors excitability. J. Neurophysiol. 113 (2), 601–619. doi:10.1152/jn.00652.2014
Guimarães, A. G., Quintans, J. S. S., and Quintans, L. J. (2013). Monoterpenes with analgesic activity--a systematic review. Phytother. Res. PTR 27 (1), 1–15. doi:10.1002/ptr.4686
Haeseler, G., Maue, D., Grosskreutz, J., Bufler, J., Nentwig, B., Piepenbrock, S., et al. (2002). Voltage-dependent block of neuronal and skeletal muscle sodium channels by thymol and menthol. Eur. J. Anaesthesiol. 19 (8), 571–579. doi:10.1017/s0265021502000923
Hall, A. C., Turcotte, C. M., Betts, B. A., Yeung, W.-Y., Agyeman, A. S., and Burk, L. A. (2004). Modulation of human GABAA and glycine receptor currents by menthol and related monoterpenoids. Eur. J. Pharmacol. 506 (1), 9–16. doi:10.1016/j.ejphar.2004.10.026
Hamill, O. P., Marty, A., Neher, E., Sakmann, B., and Sigworth, F. J. (1981). Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflüg Arch. 391 (2), 85–100. doi:10.1007/BF00656997
Harrewijn, P., Oosten, A. M. van., and Piron, P. G. (2012). Natural Terpenoids as Messengers: a multidisciplinary study of their production, biological functions and practical applications. Springer Science and Business Media.
Heijman, J., and Dobrev, D. (2018). Ion channels as part of macromolecular multiprotein complexes: clinical significance. Herzschrittmachertherapie Elektrophysiologie 29 (1), 30–35. doi:10.1007/s00399-017-0542-y
Hernandez, C. M., and Richards, J. R. (2023). Physiology, sodium channels. Treasure island: StatPearls Publishing.
Hille, B. (2001). Ion channels of excitable membranes. 3. Sunderland, Mass: Sinauer Associates is an imprint of Oxford University Press.
Ho, C., and O’Leary, M. E. (2011). Single-cell analysis of sodium channel expression in dorsal root ganglion neurons. Mol. Cell Neurosci. 46 (1), 159–166. doi:10.1016/j.mcn.2010.08.017
Hodgkin, A. L., and Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117 (4), 500–544. doi:10.1113/jphysiol.1952.sp004764
Horishita, T., and Harris, R. A. (2008). n-Alcohols inhibit voltage-gated Na+ channels expressed in Xenopus oocytes. J. Pharmacol. Exp. Ther. 326 (1), 270–277. doi:10.1124/jpet.108.138370
Hume, W. R. (1939). The pharmacologic and toxicological properties of zinc oxide-eugenol. J. Am. Dent. Assoc. 113 (5), 789–791. doi:10.14219/jada.archive.1986.0256
Hume, W. R. (1984). An analysis of the release and the diffusion through dentin of eugenol from zinc oxide-eugenol mixtures. J. Dent. Res. 63 (6), 881–884. doi:10.1177/00220345840630061301
Jeong, K. H., Lee, D.-S., and Kim, S. R. (2015). Effects of eugenol on granule cell dispersion in a mouse model of temporal lobe epilepsy. Epilepsy Res. 115, 73–76. doi:10.1016/j.eplepsyres.2015.06.001
Jiang, D., Zhang, J., and Xia, Z. (2022). Structural advances in voltage-gated sodium channels. Front. Pharmacol. 13, 908867. doi:10.3389/fphar.2022.908867
Joca, H. C., Cruz-Mendes, Y., Oliveira-Abreu, K., Maia-Joca, R. P. M., Barbosa, R., Lemos, T. L., et al. (2012). Carvacrol decreases neuronal excitability by inhibition of voltage-gated sodium channels. J. Nat. Prod. 75 (9), 1511–1517. doi:10.1021/np300050g
Jorkjend, L., and Skoglund, L. A. (1990). Effect of non-eugenol- and eugenol-containing periodontal dressings on the incidence and severity of pain after periodontal soft tissue surgery. J. Clin. Periodontol. 17 (6), 341–344. doi:10.1111/j.1600-051x.1990.tb00028.x
Karashima, Y., Damann, N., Prenen, J., Talavera, K., Segal, A., Voets, T., et al. (2007). Bimodal action of menthol on the transient receptor potential channel TRPA1. J. Neurosci. Off. J. Soc. Neurosci. 27 (37), 9874–9884. doi:10.1523/JNEUROSCI.2221-07.2007
Kirsch, G. E., Skattebøl, A., Possani, L. D., and Brown, A. M. (1989). Modification of Na channel gating by an alpha scorpion toxin from Tityus serrulatus. J. Gen. Physiol. 93 (1), 67–83. doi:10.1085/jgp.93.1.67
Kistner, K., Zimmermann, K., Ehnert, C., Reeh, P. W., and Leffler, A. (2010). The tetrodotoxin-resistant Na+ channel Na (v)1.8 reduces the potency of local anesthetics in blocking C-fiber nociceptors. Pflugers Arch. 459 (5), 751–763. doi:10.1007/s00424-010-0785-5
Kostyuk, P. G., Veselovsky, N. S., and Tsyndrenko, A. Y. (1981). Ionic currents in the somatic membrane of rat dorsal root ganglion neurons-I. Sodium currents. Neuroscience 6 (12), 2423–2430. doi:10.1016/0306-4522(81)90088-9
Kozioł, A., Stryjewska, A., Librowski, T., Sałat, K., Gaweł, M., Moniczewski, A., et al. (2014). An overview of the pharmacological properties and potential applications of natural monoterpenes. Mini-Rev Med. Chem. 14 (14), 1156–1168. doi:10.2174/1389557514666141127145820
Lacroix, J. J., Campos, F. V., Frezza, L., and Bezanilla, F. (2013). Molecular bases for the asynchronous activation of sodium and potassium channels required for nerve impulse generation. Neuron 79 (4), 651–657. doi:10.1016/j.neuron.2013.05.036
Latorre, K. L., and Baldion, P. A. (2020). Polymodal activation and desensitization of TRPV1 receptor in human odontoblasts-like cells with eugenol. Int. J. Dent. 2020, 8813979. doi:10.1155/2020/8813979
Lee, M. H., Yeon, K.-Y., Park, C.-K., Li, H.-Y., Fang, Z., Kim, M. S., et al. (2005). Eugenol inhibits calcium currents in dental afferent neurons. J. Dent. Res. 84 (9), 848–851. doi:10.1177/154405910508400913
Lee, S. H., Moon, J. Y., Jung, S. J., Kang, J. G., Choi, S. P., and Jang, J. H. (2015). Eugenol inhibits the GABAA current in trigeminal ganglion neurons. PloS One 10 (1), e0117316. doi:10.1371/journal.pone.0117316
Lee, S. P., Buber, M. T., Yang, Q., Cerne, R., Cortés, R. Y., Sprous, D. G., et al. (2008). Thymol and related alkyl phenols activate the hTRPA1 channel. Br. J. Pharmacol. 153 (8), 1739–1749. doi:10.1038/bjp.2008.85
Lenaeus, M. J., Gamal El-Din, T. M., Ing, C., Ramanadane, K., Pomès, R., Zheng, N., et al. (2017). Structures of closed and open states of a voltage-gated sodium channel. Proc. Natl. Acad. Sci. U. S. A. 114 (15), E3051–E3060. doi:10.1073/pnas.1700761114
Li, H. Y., Lee, B. K., Kim, J. S., Jung, S. J., and Oh, S. B. (2008). Eugenol inhibits ATP-induced P2X currents in trigeminal ganglion neurons. Korean J. Physiol. Pharmacol. Off. J. Korean Physiol. Soc. Korean Soc. Pharmacol. 12 (6), 315–321. doi:10.4196/kjpp.2008.12.6.315
Li, H. Y., Park, C.-K., Jung, S. J., Choi, S.-Y., Lee, S. J., Park, K., et al. (2007). Eugenol inhibits K+ currents in trigeminal ganglion neurons. J. Dent. Res. 86 (9), 898–902. doi:10.1177/154405910708600918
Lima, F. C., Peixoto-Neves, D., Gomes, M. D. M., Coelho-de-Souza, A. N., Lima, C. C., Araújo, Z. W., et al. (2011). Antispasmodic effects of eugenol on rat airway smooth muscle. Fundam. Clin. Pharmacol. 25 (6), 690–699. doi:10.1111/j.1472-8206.2010.00892.x
Liu, M., and Wood, J. N. (2011). The roles of sodium channels in nociception: implications for mechanisms of neuropathic pain. Pain Med. Malden Mass 12 (Suppl. 3), S93–S99. doi:10.1111/j.1526-4637.2011.01158.x
Magyar, J., Szentandrássy, N., Bányász, T., Fülöp, L., Varró, A., and Nánási, P. P. (2002). Effects of thymol on calcium and potassium currents in canine and human ventricular cardiomyocytes. Br. J. Pharmacol. 136 (2), 330–338. doi:10.1038/sj.bjp.0704718
Magyar, J., Szentandrássy, N., Bányász, T., Fülöp, L., Varró, A., and Nánási, P. P. (2004). Effects of terpenoid phenol derivatives on calcium current in canine and human ventricular cardiomyocytes. Eur. J. Pharmacol. 487 (1–3), 29–36. doi:10.1016/j.ejphar.2004.01.011
Maruyama, H., Yamamoto, M., Matsutomi, T., Zheng, T., Nakata, Y., Wood, J. N., et al. (2004). Electrophysiological characterization of the tetrodotoxin-resistant Na+ channel, Na(v)1.9, in mouse dorsal root ganglion neurons. Pflugers Arch. 449 (1), 76–87. doi:10.1007/s00424-004-1315-0
McCormack, K., Santos, S., Chapman, M. L., Krafte, D. S., Marron, B. E., West, C. W., et al. (2013). Voltage sensor interaction site for selective small molecule inhibitors of voltage-gated sodium channels. Proc. Natl. Acad. Sci. U. S. A. 110 (29), E2724–E2732. doi:10.1073/pnas.1220844110
Meng, E., Cai, T.-F., Zhang, H., Tang, S., Li, M.-J., Li, W.-Y., et al. (2014). Screening for voltage-gated sodium channel interacting peptides. Sci. Rep. 4, 4569. doi:10.1038/srep04569
Moreira-Lobo, D. C. A., Linhares-Siqueira, E. D., Cruz, G. M. P., Cruz, J. S., Carvalho-de-Souza, J. L., Lahlou, S., et al. (2010). Eugenol modifies the excitability of rat sciatic nerve and superior cervical ganglion neurons. Neurosci. Lett. 472 (3), 220–224. doi:10.1016/j.neulet.2010.02.009
Müller, M., Pape, H.-C., Speckmann, E.-J., and Gorji, A. (2006). Effect of eugenol on spreading depression and epileptiform discharges in rat neocortical and hippocampal tissues. Neuroscience 140 (2), 743–751. doi:10.1016/j.neuroscience.2006.02.036
Nesterenko, V. V., Zygmunt, A. C., Rajamani, S., Belardinelli, L., and Antzelevitch, C. (2011). Mechanisms of atrial-selective block of Na+ channels by ranolazine: II. Insights from a mathematical model. Am. J. Physiol. Heart Circ. Physiol. 301 (4), H1615–H1624. doi:10.1152/ajpheart.00243.2011
Nóbrega, F. F. F., Salvadori, MGSS., Masson, C. J., Mello, C. F., Nascimento, T. S., Leal-Cardoso, J. H., et al. (2014). Monoterpenoid terpinen-4-ol exhibits anticonvulsant activity in behavioural and electrophysiological studies. Oxid. Med. Cell Longev. 2014, 703848. doi:10.1155/2014/703848
Ogata, N., and Tatebayashi, H. (1993). Kinetic analysis of two types of Na+ channels in rat dorsal root ganglia. J. Physiol. 466, 9–37. doi:10.1113/jphysiol.1993.sp019706
Ohkubo, T., and Shibata, M. (1997). The selective capsaicin antagonist capsazepine abolishes the antinociceptive action of eugenol and guaiacol. J. Dent. Res. 76 (4), 848–851. doi:10.1177/00220345970760040501
Olivoto, R. R., Damiani, C. E. N., Kassouf, S. I., Lofrano-Alves, M. S., Oliveira, M. A., and Fogaça, R. T. H. (2014). Effects of eugenol on resting tension of rat atria. Braz J. Med. Biol. Res. Rev. Bras. Pesqui. Medicas E Biol. 47 (4), 328–333. doi:10.1590/1414-431x20143472
Pan, X., Li, Z., Zhou, Q., Shen, H., Wu, K., Huang, X., et al. (2018). Structure of the human voltage-gated sodium channel Nav1.4 in complex with β1. Science 362 (6412), eaau2486. doi:10.1126/science.aau2486
Park, C.-K., Kim, K., Jung, S. J., Kim, M. J., Ahn, D. K., Hong, S.-D., et al. (2009). Molecular mechanism for local anesthetic action of eugenol in the rat trigeminal system. Pain 144 (1–2), 84–94. doi:10.1016/j.pain.2009.03.016
Park, C.-K., Li, H. Y., Yeon, K.-Y., Jung, S. J., Choi, S.-Y., Lee, S. J., et al. (2006). Eugenol inhibits sodium currents in dental afferent neurons. J. Dent. Res. 85 (10), 900–904. doi:10.1177/154405910608501005
Park, S.-H., Sim, Y.-B., Lee, J.-K., Kim, S.-M., Kang, Y.-J., Jung, J.-S., et al. (2011). The analgesic effects and mechanisms of orally administered eugenol. Arch. Pharm. Res. 34 (3), 501–507. doi:10.1007/s12272-011-0320-z
Peixoto-Neves, D., Leal-Cardoso, J. H., and Jaggar, J. H. (2014). Eugenol dilates rat cerebral arteries by inhibiting smooth muscle cell voltage-dependent calcium channels. J. Cardiovasc Pharmacol. 64 (5), 401–406. doi:10.1097/FJC.0000000000000131
Peixoto-Neves, D., Wang, Q., Leal-Cardoso, J. H., Rossoni, L. V., and Jaggar, J. H. (2015). Eugenol dilates mesenteric arteries and reduces systemic BP by activating endothelial cell TRPV4 channels. Br. J. Pharmacol. 172 (14), 3484–3494. doi:10.1111/bph.13156
Prinz, H. (2010). Hill coefficients, dose-response curves and allosteric mechanisms. J. Chem. Biol. 3 (1), 37–44. doi:10.1007/s12154-009-0029-3
Roy, M. L., and Narahashi, T. (1992). Differential properties of tetrodotoxin-sensitive and tetrodotoxin-resistant sodium channels in rat dorsal root ganglion neurons. J. Neurosci. Off. J. Soc. Neurosci. 12 (6), 2104–2111. doi:10.1523/JNEUROSCI.12-06-02104.1992
Sancheti, J., Shaikh, M. F., Chaudhari, R., Somani, G., Patil, S., Jain, P., et al. (2014). Characterization of anticonvulsant and antiepileptogenic potential of thymol in various experimental models. Naunyn Schmiedeb. Arch. Pharmacol. 387 (1), 59–66. doi:10.1007/s00210-013-0917-5
Sensch, O., Vierling, W., Brandt, W., and Reiter, M. (2000). Effects of inhibition of calcium and potassium currents in Guinea-pig cardiac contraction: comparison of beta-caryophyllene oxide, eugenol, and nifedipine. Br. J. Pharmacol. 131 (6), 1089–1096. doi:10.1038/sj.bjp.0703673
Seo, H., Li, H. Y., Perez-Reyes, E., and Lee, J.-H. (2013). Effects of eugenol on T-type Ca2+ channel isoforms. J. Pharmacol. Exp. Ther. 347 (2), 310–317. doi:10.1124/jpet.113.207936
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 (1), 16–21. doi:10.1016/j.tcm.2010.03.002
Sheets, M. F., and Hanck, D. A. (2003). Molecular action of lidocaine on the voltage sensors of sodium channels. J. Gen. Physiol. 121 (2), 163–175. doi:10.1085/jgp.20028651
Sheets, M. F., and Hanck, D. A. (2007). Outward stabilization of the S4 segments in domains III and IV enhances lidocaine block of sodium channels. J. Physiol. 582 (Pt 1), 317–334. doi:10.1113/jphysiol.2007.134262
Shiers, S., Klein, R. M., and Price, T. J. (2020). Quantitative differences in neuronal subpopulations between mouse and human dorsal root ganglia demonstrated with RNAscope in situ hybridization. Pain 161 (10), 2410–2424. doi:10.1097/j.pain.0000000000001973
Silva-Alves, K. S., Ferreira-da-Silva, F. W., Peixoto-Neves, D., Viana-Cardoso, K. V., Moreira-Júnior, L., Oquendo, M. B., et al. (2013). Estragole blocks neuronal excitability by direct inhibition of Na+ channels. Braz J. Med. Biol. Res. 46 (12), 1056–1063. doi:10.1590/1414-431X20133191
Soares, P. M. G., Lima, R. F., de Freitas, P. A., Souza, E. P., Assreuy, A. M. S., and Criddle, D. N. (2007). Effects of anethole and structural analogues on the contractility of rat isolated aorta: involvement of voltage-dependent Ca2+-channels. Life Sci. 81 (13), 1085–1093. doi:10.1016/j.lfs.2007.08.027
Subramanyam, P., and Colecraft, H. M. (2015). Ion channel engineering: perspectives and strategies. J. Mol. Biol. 427 (1), 190–204. doi:10.1016/j.jmb.2014.09.001
Sucher, N. J., and Carles, M. C. (2015). A pharmacological basis of herbal medicines for epilepsy. Epilepsy Behav. EB 52 (Pt B), 308–318. doi:10.1016/j.yebeh.2015.05.012
Sunami, A., and Hiraoka, M. (1996). Blockade of cardiac Na+ channels by a charged class I antiarrhythmic agent, bisaramil: possible interaction of the drug with a pre-open closed state. Eur. J. Pharmacol. 312 (2), 245–255. doi:10.1016/0014-2999(96)00461-x
Takahashi, K., Yoshida, T., and Wakamori, M. (2021). Mode-selective inhibitory effects of eugenol on the mouse TRPV1 channel. Biochem. Biophys. Res. Commun. 556, 156–162. doi:10.1016/j.bbrc.2021.03.126
Takaishi, M., Uchida, K., Fujita, F., and Tominaga, M. (2014). Inhibitory effects of monoterpenes on human TRPA1 and the structural basis of their activity. J. Physiol. Sci. JPS 64 (1), 47–57. doi:10.1007/s12576-013-0289-0
Tan, Z.-Y., Piekarz, A. D., Priest, B. T., Knopp, K. L., Krajewski, J. L., McDermott, J. S., et al. (2014). Tetrodotoxin-resistant sodium channels in sensory neurons generate slow resurgent currents that are enhanced by inflammatory mediators. J. Neurosci. Off. J. Soc. Neurosci. 34 (21), 7190–7197. doi:10.1523/JNEUROSCI.5011-13.2014
Tavvabi-Kashani, N., Hasanpour, M., Baradaran Rahimi, V., Vahdati-Mashhadian, N., and Askari, V. R. (2024). Pharmacodynamic, pharmacokinetic, toxicity, and recent advances in Eugenol’s potential benefits against natural and chemical noxious agents: a mechanistic review. Toxicon Off. J. Int. Soc. Toxinology 238, 107607. doi:10.1016/j.toxicon.2024.107607
Teixeira-Fonseca, J. L., Santos-Miranda, A., da Silva, J. B., Marques, L. P., Joviano-Santos, J. V., Nunes, P. I. C., et al. (2021). Eugenol interacts with cardiac sodium channel and reduces heart excitability and arrhythmias. Life Sci. 282, 119761. doi:10.1016/j.lfs.2021.119761
Theile, J. W., and Cummins, T. R. (2011). Recent developments regarding voltage-gated sodium channel blockers for the treatment of inherited and acquired neuropathic pain syndromes. Front. Pharmacol. 2, 54. doi:10.3389/fphar.2011.00054
Tomasova, L., Pavlovicova, M., Malekova, L., Misak, A., Kristek, F., Grman, M., et al. (2015). Effects of AP39, a novel triphenylphosphonium derivatised anethole dithiolethione hydrogen sulfide donor, on rat haemodynamic parameters and chloride and calcium Cav3 and RyR2 channels. Nitric Oxide Biol. Chem. 46, 131–144. doi:10.1016/j.niox.2014.12.012
Vacher, H., Mohapatra, D. P., and Trimmer, J. S. (2008). Localization and targeting of voltage-dependent ion channels in mammalian central neurons. Physiol. Rev. 88 (4), 1407–1447. doi:10.1152/physrev.00002.2008
Vilin, Y. Y., and Ruben, P. C. (2001). Slow inactivation in voltage-gated sodium channels: molecular substrates and contributions to channelopathies. Cell Biochem. Biophys. 35 (2), 171–190. doi:10.1385/CBB:35:2:171
Wang, Z.-J., Tabakoff, B., Levinson, S. R., and Heinbockel, T. (2015). Inhibition of Nav1.7 channels by methyl eugenol as a mechanism underlying its antinociceptive and anesthetic actions. Acta Pharmacol. Sin. 36 (7), 791–799. doi:10.1038/aps.2015.26
Wu, J. V., and Kendig, J. J. (1998). Differential sensitivities of TTX-resistant and TTX-sensitive sodium channels to anesthetic concentrations of ethanol in rat sensory neurons. J. Neurosci. Res. 54 (4), 433–443. doi:10.1002/(SICI)1097-4547(19981115)54:4<433::AID-JNR1>3.0.CO;2-A
Xiao, J., Bondarenko, V., Wang, Y., Suma, A., Wells, M., Chen, Q., et al. (2021). Regulation and drug modulation of a voltage-gated sodium channel: pivotal role of the S4-S5 linker in activation and slow inactivation. Proc. Natl. Acad. Sci. U. S. A. 118 (28), e2102285118. doi:10.1073/pnas.2102285118
Yang, B. H., Piao, Z. G., Kim, Y.-B., Lee, C.-H., Lee, J. K., Park, K., et al. (2003). Activation of vanilloid receptor 1 (VR1) by eugenol. J. Dent. Res. 82 (10), 781–785. doi:10.1177/154405910308201004
Keywords: eugenol, lidocaine, voltage-gated sodium channels, inhibition, state-dependent interactions, dorsal root ganglion neurons, patch-clamp technique
Citation: Moreira-Junior L, Leal-Cardoso JH, Cassola AC and Carvalho-de-Souza JL (2024) Eugenol and lidocaine inhibit voltage-gated Na+ channels from dorsal root ganglion neurons with different mechanisms. Front. Pharmacol. 15:1354737. doi: 10.3389/fphar.2024.1354737
Received: 12 December 2023; Accepted: 08 April 2024;
Published: 26 June 2024.
Edited by:
Klaus Schicker, Medical University of Vienna, AustriaReviewed by:
Verena Burtscher, Washington University in St. Louis, United StatesHannes Todt, Medical University of Vienna, Austria
Copyright © 2024 Moreira-Junior, Leal-Cardoso, Cassola and Carvalho-de-Souza. 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) and the copyright owner(s) 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: Joao Luis Carvalho-de-Souza, jlvet2005@gmail.com
†Present address: Joao Luis Carvalho-de-Souza, Department of Anesthesiology, University of Arizona, Tucson, AZ, United States