4. Discussion
A broad number of medicinal plants have been used for centuries as therapeutic tools in TCM. However, neither the active compounds of many of these plants nor their mechanisms of action are yet well-understood. We have now studied the effect of Pm, an isosteroidal alkaloid considered one of the main bioactive molecules of Fb, on nAChRs.
Remarkably, Pm decreased
IACh in a dose-dependent manner, with an IC
50 in the low micromolar range (circa 3 and 1 µM for
Ip and
Iss, respectively;
Figure 1). These IC
50s are markedly lower than the values previously reported for Pm blockade of voltage-dependent potassium channels. Thus, Pm IC
50 for blocking Kv1.2 was 472 µM, it was 354 µM for Kv1.3 (142 µM if measured 150 ms after the current peak), and even much higher for Kv1.4 to Kv1.8 channels
[10]; additionally, Pm inhibited the potassium channel hERG, with an IC
50 of 44 µM, most likely by enhancing its inactivation
[11]. Pm also blocked the Nav1.7 channel (IC
50, 47 µM), demonstrating use-dependent inhibition, like the blocking mechanism of lidocaine on this channel
[10]. Pm effects on hERG channels are particularly remarkable, since these channels play a key role in myocardial repolarization and, therefore, their inhibition might cause serious cardiac arrhythmias. Despite this, humans have used
Fb as a therapeutic herb for centuries, being considered safe for consumption. Consequently, intake of
Fb should not markedly affect the activity of hERG channels, despite that Pm, one of its main bioactive compounds, blocks these channels with an IC
50 of 44 µM
[11]. Nevertheless, Pm has a very low oral bioavailability
[9] and thus its plasma concentration after
Fb intake should be fairly low. Actually, Pm content in bulbs of
Fritillaria ussuriensis and
thunbergii ranged from 0.58 to 1.2 mg/g and oral administration of powder from these
Fritillaria plants to dogs (1 g/kg) raised Pm plasma concentration to a maximum of 100–200 nM
[24]. Accordingly, a similar Pm bioavailability was reported after oral administration of
Fritillaria thunbergii extracts in rats, with peak plasma concentrations of Pm of roughly 100 nM
[25]. These Pm concentrations are several orders of magnitude lower than the IC
50s reported for sodium or potassium channels (including hERG), muscarinic receptors, or acethylcholinesterase
[8]. Interestingly, submicromolar Pm concentrations elicit a significant inhibition of nAChRs (roughly 20%;
Figure 1) and therefore this family of LGIC might be relevant targets of its actions. Of note, we have assessed Pm actions on muscle-type nAChRs, because they are the prototype member of this family of receptors, but Pm might have different affinities for related receptors, as the homomeric α7. In fact, α7 nAChRs are largely expressed in non-neuronal tissues, including macrophages, and exert powerful anti-inflammatory actions
[16]. Moreover, other nAChR subtypes, such as α4β2 and/or α9α10, may also play a role in modulating inflammatory processes and even in chronic pain
[15]. Noticeably, Pm displayed a differential affinity for different LGICs, even of the same family. Thus, whereas Pm inhibits muscle-type nAChRs with and IC
50 close to 1 µM, GABA
A receptors were almost not affected by Pm at concentrations up to 100 µM.
Pm exerted a non-competitive inhibition on muscle-type nAChRs, since
IAChs halved when co-applying Pm, at its IC
50, with different ACh concentrations (
Figure 2A,B). Several blockade mechanisms seem involved in this non-competitive inhibition of nAChRs by Pm, which is coherent with the multiple binding sites predicted by the docking simulations. First, open-channel blockade, as
IACh inhibition by Pm was voltage-dependent, the more hyperpolarized the cell, the more pronounced the blockade (
Figure 5). Actually, this is what would be expected for a positively charged molecule plugging the channel pore. In agreement with this, the docking assays predicted some Pm clusters located within the channel pore, both in the open and closed conformations (
Figure 8). Thus, the open-channel blockade of nAChRs mediated by Pm resembles that mediated by some LAs, such as lidocaine
[26] or tetracaine
[27]. Nevertheless, the kinetics of open-channel blockade of nAChR elicited by Pm was rather slow. Thus, at −60 mV, the τ of open channel blockade elicited by 1 µM Pm (close to its IC
50) was over 2 s and even by 5 µM Pm (eliciting roughly 70% of
Iss blockade) the time constant was above 1 s (
Figure 6). In contrast, at the same membrane potential, 0.7 µM tetracaine (close to its IC
50) blocked open nAChRs with a time course of roughly 300 ms
[27]. These differences in time constants between Pm and tetracaine are most likely related to their distinct molecular sizes (Pm molecular weight is over 60% greater than that of tetracaine). Second, Pm enhanced nAChR desensitization, as evidenced by: (i) acceleration of
IACh decay when co-applying ACh with Pm at concentrations of 0.5 µM, or above (
Figure 3A
2,B
2), and (ii) shortening of the
IACh aTtp, which correlated with the acceleration of
IACh decay (
Figure 3A
1,B
1). Likewise, lidocaine decreased the
IACh aTtp only at concentrations that enhanced
IACh decay
[26]. Furthermore, DMA, a lidocaine analog, both sped up
IACh decay and shortened the aTtP
[28]. By contrast, diethylamine (DEA), a lidocaine analog that mainly blocks nAChR by open channel blockade neither accelerates
IACh decay nor decreases aTtP
[21], and (iii) deceleration of
IACh deactivation, which was dependent on Pm concentration and displayed a good correlation with the rate of
IACh decay (
Figure 4A, B). The deceleration of
IACh deactivation when Pm remained in the solution strongly supports that Pm enhanced nAChR desensitization, since desensitized nAChRs display higher affinity for the agonist
[22][27][29]. Third, Pm elicited the blockade of resting (closed) nAChRs. This effect was unraveled by applying Pm before challenging the cells with ACh alone. This protocol, which allowed Pm to act only on resting (closed) nAChRs, elicited a mild nAChR blockade, mostly at negative potentials. Actually, at positive potentials,
IACh only decreased by Pm pre-application when rising its concentration to 5 µM (
Figure 7). In agreement with this,
IACh inhibition by 5 µM Pm was slightly higher when Pm was pre-applied and then co-applied with ACh than when solely co-applied with ACh (
Figure 7C).
Figure 3. Pm accelerates IACh decay and shortens the time to reach Ip. (A1,A2) Superimposed IAChs elicited by 10 µM ACh either alone (black and grey recordings) or together with different Pm concentrations (shown at right). IAChs were scaled to the same Ip amplitude to better compare the differences in time to reach Ip (aTtP; A1) and kinetics of IACh decay after Ip (A2). (B1,B2) Column graphs displaying Pm effects on aTtP (B1) and τ-values of IACh-decay (B2). (*) indicates significant differences among IAChs in presence of Pm (colored columns; same color code as in (A1,A2)) and their control values (Ctr, black column; p < 0.05, ANOVA and Bonferroni t-test). Note that post-control values (after Pm applications; grey column) were similar to control ones. Each point is the average of 4–24 cells (N = 3–12).
Figure 4. IACh decay and deactivation kinetics depend on Pm concentration. (A1,A2,A3) Representative IAChs elicited by 100 µM ACh either alone (black recording) or together with 1 (orange) or 5 (purple) µM Pm (A1). Pm superfusion remained 12 s after ACh washout (as indicated by the application bars). These recordings were normalized to either the same Ip, to better compare their IACh decay (A2), or the same Iss, to facilitate comparisons of deactivation kinetics (A3). (B1,B2) Column bar plots displaying the effect of 1 (orange) or 5 µM (purple) Pm on the IACh decay time-constant (τDesensitization; (B1)) and the deactivation kinetics (τDeactivation; (B2)), as compared to control IAChs (in the presence of ACh alone; black). (*) indicates significant differences with the control group (p < 0.05, paired t-test) and (#) indicates differences between 1 and 5 µM Pm groups (n = 9, N = 3; the same cells for all comparisons; p < 0.05, paired t-test). Notice that Pm accelerated the desensitization rate and slowed down the deactivation kinetics.
Figure 5. nAChR blockade by Pm is voltage-dependent. (A) IAChs evoked by 10 µM ACh alone (black recording) or in the presence of either 1 (orange) or 5 µM (purple) Pm when applying voltage pulses from −120 to +60 mV, as illustrated underneath. (B) Net i/v relationship of IAChs elicited by the protocol shown in (A). Black symbols are for control IAChs, whereas those evoked in the presence of Pm are drawn in either orange (+1 µM Pm) or purple (+5 µM Pm). Net IAChs were normalized as the percentage of their control IACh at −60 mV (n = 5–11; N = 2–3). (C) Plot displaying the IACh remnant after co-application of either 1 (orange) or 5 µM (purple) Pm (IACh+Pm), normalized to their control (IACh), versus the membrane potential (same cells as in (B)). Notice that 1 µM Pm, in contrast to 5 µM, did not significantly decrease IACh at +60 mV.
Figure 6. Kinetics of the voltage-dependent blockade of nAChR by Pm. (A1,B1) IAChs elicited by 10 µM ACh either alone (black recordings) or in the presence of 1 (orange; (A1)) or 5 µM (purple; (B1)) Pm, at −60 mV. A 2 s voltage jump to +40 mV was given at the IACh plateau to unplug the channel pore of the positively-charged Pm. Membrane leak-currents (recorded in the absence of ACh) have been subtracted. (A2,B2) Zooming in to the area indicated by the arrows in panels (A1,B1) (just after the voltage jump). The τ of the voltage-dependent blockade of nAChRs by Pm was determined by fitting an exponential function (green curve over the recording) to the net IACh decay. Before fitting, the smaller and slower IAChs evoked by ACh alone (black recordings of panels (A1,B1)) were subtracted from the IAChs in the presence of Pm. (C) Column-graph of τ values of the voltage-dependent blockade of nAChR by 1 and 5 µM Pm (same color code as in panels (A,B)). (*) indicates significant differences of τ values between both Pm concentrations (p < 0.05, t-test). Data are for 10 (N = 3) and 7 (N = 2) oocytes for 1 µM and 5 µM Pm, respectively.
Figure 7. Effect of Pm-application timing and holding potential on nAChR blockade. (A1–A3) IAChs elicited at −60 mV (downward deflections) and at +40 mV (upward deflections) by co-application of 10 µM ACh and 1 µM Pm (A1), solely Pm pre-application before superfusing the agonist (A2) or Pm pre-application followed by its co-application with ACh (A3). (B1–B3) As in panels (A1–A3), but in the presence of 5 µM Pm instead of 1 µM. (C) Column graph shows the percentages of Ip inhibition by Pm when applied as indicated in panels (A1–A3,B1–B3), at −60 mV (on the left) and +40 mV (on the right). (*) indicates significant differences between Ip inhibition elicited by 1 and 5 µM Pm (p < 0.05, t-test). (#) denotes significant differences, for each Pm concentration, among the percentages of Ip inhibition elicited by ACh and Pm co-application and other Pm-application protocols, at the same holding potential (p < 0.05, ANOVA and Bonferroni t-test). Each column is the average of 5–11 and 5–17 oocytes, for −60 mV and +40 mV, respectively.
Figure 8. Predicted binding sites for Pm-nAChR complexes. (A,B) Lateral view of the nAChR displaying the main Pm (labelled in cyan) clusters bound to the open (A) and closed (B) conformations. The predicted loci are numbered consecutively, beginning in the transmembrane (TMD) and later in the extracellular (ECD) domain. (C, D) Top view of the nAChR (from the synaptic cleft) displaying representative Pm clusters binding to residues located within the channel pore, TMD, and ECD in the open (C) and closed (D) conformations. The inset, in the upper right corner, displays the nAChR subunits with their color code.
Our virtual docking and MD simulations used as a template the structure of
Torpedo nAChRs in the open and closed conformations released by Unwin’s group
[30][31]. However, a new structural model of the nAChR from
Torpedo, at higher resolution and stabilized in the closed conformation by α-bungarotoxin, has been recently disclosed
[32]. The new structural model (pdb entry 6UWZ) share large similarities with the Unwin’s model for the closed conformation, though there are certain differences between them. Mostly, they differ in the upper portion of the pore, which is more constricted in the new model, and in the δ subunit arrangement
[33]. It seems that these discrepancies arise because of differences in the lipid matrix surrounding the nAChR. Actually, cholesterol interactions with the nAChR are apparently essential for stabilizing its structure and the absence of cholesterol (as in the model of Rahman et al.
[32]) leads to a more compact arrangement of TM helices (displacement of helices circa 1–3 Å;
[33]). Noticeably, the major lipid present in electroplax membranes rich in nAChRs is cholesterol
[34] and purified nAChRs from
T. marmorata and
E. electricus interact preferentially with cholesterol rather than with either phospholipid monolayers or other sterols
[35]. Moreover, nAChRs in native electroplax membranes are arranged as dimers, linked by their δ-subunits. This interaction between neighboring nAChRs might account for the differences in the δ-subunit between both structural models since dimers were reduced to monomeric receptors in the Rahman’s model. In fact, we chose for our structural studies Unwin’s models because of: (i) the structures for the open and closed conformations are available, (ii) the nAChR is present in their original membrane, and (iii) we have significant experience correlating structural and functional results using these commonly accepted models; actually, Unwin’s models have so far provided a coherent correlation with our functional results
[21][28][27][36].
The virtual docking assays predicted Pm binding to the nAChR at different sites of the TMD and ECD in the open conformation. Most Pm clusters were located at the TMD, at inter- and intra-subunit crevices, although some of them located into the channel pore (
Figure 8A,C). The binding energies estimated for these clusters were rather high (from −9.3 to −12.87 kcal/mol; see
Supplementary Table S1), pointing out that Pm has a high affinity for the nAChR. Remarkably, MD simulations of nAChRs in the open conformation indicate that Pm binding to the nAChR either at the TMD (i.e., cluster 6) or at the ECD (i.e., cluster 14) markedly decreased both the volume and the number of water molecules at the hydrophobic gate region of the channel pore (
Figure 9). Furthermore, virtual docking and MD assays pointed out that Pm also interacts with the nAChR in the resting conformation, binding to residues located at both the TMD and the ECD (
Figure 8B,D and
Figure 9E,F). Interestingly, the structural changes of the nAChR induced by Pm, as predicted by docking and MD assays, are in good agreement with the functional changes elicited by Pm on this receptor. Actually, the structural and functional results can be correlated as follows: (i) the high binding energies computed accounted for the high potency of Pm blocking nAChRs, (ii) Pm interaction with residues located within the channel pore should trigger the open-channel blockade, (iii) Pm binding to different sites at the nAChR might explain both its heterogeneity of actions on nAChRs and the effect-dependence on Pm concentration, and (iv) Pm binding to the nAChR in the closed conformation might underlie the blockade of resting nAChR. Consequently, the good correlation between structural simulations and electrophysiological results strongly suggests that Pm actually blocks nAChRs by the different aforementioned mechanisms.
Figure 9. Empty volume (left panels) and number of water molecules (right panels) within the hydrophobic gate of the nAChR pore region (between Val255 and Glu262 of the alpha subunit; see inset in the upper right corner), through the period of 40 to 100 ns of MD simulations. Top panels show the empty volume (A) and the number of water molecules (B) at the hydrophobic gate region in control nAChR (in the absence of Pm), both in the open (yellow) and the closed (blue) conformations; also, it displays the effect of ACh on the closed conformation (green). Middle panels (C,D) display the effect of Pm on these parameters when located at representative sites of the nAChR in the open conformation: within the channel pore (Pm 1), at the TMD (Pm 6), and the ECD (Pm 14). Lower panels (E,F) demonstrate the effect of Pm at representative loci of the nAChR in the closed conformation: inside the channel (Pm 1), at TMD (Pm 2), and at ECD (Pm 13).
Since
Fb has been used as therapeutic herb for thousands of years, there is a strong support for its beneficial effects and its weak (or lack) of toxicity. However, neither the identity of all its bioactive compounds nor their mechanisms of action are yet well known. Now, we report here that Pm, considered one of the main bioactive compounds from
Fritillaria, exerts a powerful inhibition of muscle-type nAChRs, which, as far as we know, is the first report demonstrating that Pm might modulate LGICs, besides acting on other targets as voltage-dependent channels or metabotropic receptors. It remains to be unraveled if Pm might modulate other nAChRs, including the homomeric α7, which is broadly expressed in immune cells and has been related to powerful anti-inflammatory actions
[16]. Furthermore, both mecamylamine, a non-competitive antagonist of α7 nAChRs, and 1-ethyl-4-(3-(bromo)phenyl)piperazine, which promotes α7 desensitization, reduce pro-inflammatory responses
[37]. We have now demonstrated that Pm inhibits muscle-type nAChRs and it can be hypothesized that it might modulate, though in a different way or with a different potency, either α7 or other nAChRs. In this regard, lidocaine exerted similar inhibitory actions on muscle-type nAChRs
[26] and on neuronal nAChRs expressed in autonomic-ganglia neurons
[38]. Alternatively, it could be that different bioactive compounds from
Fb accounted for its anti-inflammatory actions, as anthocyanin pigments, which are flavonoids present in
Fb
[39]. Noticeably, some flavonoids act as positive allosteric modulators of α7 nAChRs, although without affecting desensitization
[6], and their enhancement of α7 activity has been proposed as a therapeutic strategy for inflammatory disorders
[40].