Abstract
Nicotinic receptors (acetylcholine receptors, AChRs) play key roles in synaptic transmission throughout the nervous system. AChRs mediate neuromuscular transmission in nematodes, and they are targets for antiparasitic drugs. The anthelmintic agents levamisole and pyrantel, which are potent agonists of nematode muscle AChRs, are partial agonists of mammalian muscle AChRs. To further explore the structural basis of the differential activation of AChR subtypes by anthelmintics, we studied the activation of α7 AChRs using the high-conductance form of the α7-5-hydroxytryptamine-3A receptor, which is a good model for pharmacological studies involving the extracellular region of α7. Macroscopic and single-channel current recordings show that levamisole is a weak agonist of α7. It is interesting that pyrantel is a more potent agonist of α7 than acetylcholine (ACh). To identify determinants of this differential activation, we replaced residues of the complementary face of the binding site by the homologous residues in the muscle ϵ subunit and evaluated changes in activation. The mutation Q57G does not affect the activation by either ACh or levamisole. However, it increases EC50 values and decreases the maximal response to pyrantel. Kinetic analysis shows that gating of the mutant channel activated by pyrantel is profoundly impaired. The decreased sensitivity of α7-Q57G to pyrantel agrees with its weak action at muscle AChRs, indicating that when glycine occupies position 57, as in the mammalian muscle AChR, pyrantel behaves as a partial agonist. Thus, position 57 located at the complementary face of the binding site plays a key role in the selective activation of AChRs by pyrantel.
Cysteine-loop receptors are pentameric ligand-gated ion channels that play key roles in chemical synapses throughout the nervous system. In vertebrates, the cysteine-loop receptor superfamily consists of two families of cation-selective channels: the nicotinic acetylcholine (AChR) and the 5-hydroxytryptamine type 3 receptors (5HT3Rs); and two families of anionic channels: the γ-aminobutyric acid and the glycine receptors (Le Novére and Changeux, 2001; Lester et al., 2004). In invertebrates, this superfamily also includes γ-aminobutyric acid-gated cationic channels and ACh-, serotonin-, glutamate-, and histamine-gated chloride channels (Jones and Sattelle, 2003; Putrenko et al., 2005).
AChR subunits are classified as α, which contain a disulfide bridge involved in the recognition and binding of agonists, and non-α subunits, which lack this motif. To date, several α and non-α subunits have been identified. Receptors can be either heteromeric, composed of α and non-α subunits, or homomeric, composed of five identical α subunits. Each subunit is divided into an N-terminal or ligand-binding domain and a transmembrane region. The binding sites are located at the interface between two subunits (Brejc et al., 2001; Sine, 2002). Results from affinity labeling and sitedirected mutagenesis studies on AChR, which were later supported by the structural model of the acetylcholine binding protein (Brejc et al., 2001), revealed that the binding site is formed by three (called loops A-C) and four (called loops D-G) linearly separated regions of the α (principal face) and the adjacent subunit (complementary face), respectively. The principal face is always formed by the α subunit, whereas the complementary face is formed by the γ or ϵ and δ subunits in the muscle AChR and by the α subunit in the homomeric α7 AChR (Sine, 2002).
The mammalian muscle AChR has a composition of α2βγδ in embryonic and denervated muscle and α2βϵδ in adult muscle. Nematode muscle AChRs are of considerable interest because they are targets for antiparasitic drugs. The freeliving helminth Caenorhabditis elegans, which is a model for the study of parasitic nematodes, has one of the largest and most diverse AChR gene families known (Jones and Sattelle, 2003). The subunit composition of each AChR subtype is not known. Levamisole and pyrantel are two widely used anthelmintic drugs that cause hypercontraction of body wall muscles of nematode worms. The efficacy of these drugs is based on their ability to act as full agonists of AChRs in nematodes (Martin et al., 1996). In contrast, they are very low-efficacy agonists of the mammalian muscle AChR (Rayes et al., 2001, 2004). The molecular basis of this selectivity has not been fully elucidated. To further understand the structural basis of the differential activation of AChR subtypes by anthelmintic agents we here explore how these drugs activate homomeric α7 AChRs.
Studies were carried out using the high-conductance form of the α7-5HT3A receptor, which has two advantages that make it an excellent model for pharmacological studies of α7 receptors: it shows, in contrast to α7, high surface expression in mammalian cells; and it allows currents to be detected at the single-channel level (Rayes et al., 2005). Our results reveal that, as in the mammalian muscle AChR, levamisole is a very low-efficacy agonist of α7 AChRs. However, pyrantel acts at α7 as a more potent and equally efficacious agonist with respect to ACh. To identify the structural basis of the different actions of pyrantel at muscle and α7 AChRs, we performed site-directed mutagenesis at candidate residues. Because residues that form the principal face of the binding site are highly conserved between α7 and muscle AChRs (Sine, 2002), we studied residues at the complementary face of the binding site. Our results show that position 57 located at the complementary binding site face is a main determinant of the selectivity of pyrantel for AChR subtypes.
These studies provide new insights into the differential activation of AChR subtypes by anthelmintic agents and identify specific differences that can be exploited therapeutically. They also help to understand how parasites develop resistance to antiparasitic drugs.
Materials and Methods
Site-Directed Mutagenesis and Expression of α7-5HT3A Receptors. Mutant subunits were constructed using the QuikChange Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA) and were confirmed by restriction enzyme analysis and sequencing. The highconductance form (HC) of the α7-5HT3A receptor (Eisele et al., 1993) was constructed as described in Rayes et al. (2005). In brief, three arginine residues responsible for the low conductance of the serotonin type 3A receptor (5HT3A) were mutated as described before (Kelley et al., 2003). BOSC cells were transfected with α7-5HT3A HC subunit cDNAs using calcium phosphate precipitation, as described previously (Bouzat et al., 1994, 2000, 2002). Cells were used for single-channel and macroscopic current measurements 1 or 2 days after transfection.
Patch-Clamp Recordings and Kinetic Analysis. Single-channel recordings were obtained in the cell-attached configuration (Hamill et al., 1981) at a membrane potential of -70 mV and at 20°C, essentially as described before (Bouzat et al., 2004; Rayes et al., 2005). The bath and pipette solutions contained 142 mM KCl, 5.4 mM NaCl, 0.2 mM CaCl2, and 10 mM HEPES, pH 7.4. Solutions free of magnesium and with low-calcium were used to minimize channel block by divalent cations. Patch pipettes were pulled from 7052 capillary tubes (Garner Glass, Claremont, CA) and coated with Sylgard (Dow Corning, Midland, MI). ACh and anthelmintic drugs were added to the pipette solution.
Single-channel currents were recorded using an Axopatch 200 B patch-clamp amplifier (Axon Instruments, Sunnyvale, CA), digitized at 5-μs intervals with the PCI-6111E interface (National Instruments, Austin, TX), recorded to the computer hard disk using the program Acquire (Bruxton Corporation, Seattle, WA), and detected by the half amplitude threshold criterion using the program TAC 4.0.10 (Bruxton Corporation) at a final bandwidth of 10 kHz (Bouzat et al., 2000). Open-, closed-, and burst-duration histograms were plotted using a logarithmic abscissa and a square root ordinate (Sigworth and Sine, 1987) and fitted to the sum of exponential functions by maximum likelihood using the program TACFit (Bruxton Corporation). Bursts were identified as a series of closely spaced openings preceded and followed by closed intervals longer than a critical duration, which was taken as the point of intersection of the second and third closed components in the closed-time histogram. The burst duration and number of openings per burst were determined from the burst-duration histogram, in which the component corresponding to isolated openings is clearly distinguished from that corresponding to bursts.
For whole-cell recordings, the pipette solution contained 140 mM KCl, 5 mM EGTA, 5 mM MgCl2, and 10 mM HEPES, pH 7.3. Extracellular solution contained 150 mM NaCl, 5.6 mM KCl, 0.5 mM CaCl2, and 10 mM HEPES, pH 7.3. A series of applications of extracellular solution containing ACh, levamisole, or pyrantel was applied to the cell as described before (Liu and Dilger, 1991; Spitzmaul et al., 2001). Macroscopic currents were filtered at 5 kHz and digitized at 20 kHz. Data analysis was performed using the IgorPro software (WaveMetrics Inc., Lake Oswego, OR). The ensemble mean current was calculated for 5 to 10 individual current traces. Mean currents were usually fitted by a single exponential function: I(t) = I0 exp(-t/τd) + I∞, where I0 and I∞ are the peak and the steady-state current values, respectively, and τd is the decay time constant. Current records were aligned with each other at the point at which the current had increased to 50% of its maximum level. EC50 and Hill coefficient values were calculated by nonlinear regression analysis using the Hill equation: I/Imax = 1/[1 + (EC50/L)nH], where EC50 is the agonist concentration that elicits the half-maximal response, nH is the Hill coefficient, and L is the agonist concentration.
Kinetic analysis was performed as described before (Wang et al., 1997; Bouzat et al., 2000, 2002; Rayes et al., 2005). The analysis was restricted to bursts of channel openings, each reflecting the activity of a single channel. For the kinetic analysis, bursts containing three or more openings were selected. Typically, the burst terminator was 2 ms. For each recording, kinetic homogeneity was determined by selecting bursts on the basis of distributions of mean open duration, mean closed duration, and open probability (Wang et al., 1997; Bouzat et al., 2000; Rayes et al., 2005). Bursts showing mean open channel durations, mean closed channel durations, and open probabilities within 2 S.D. of the mean of the major component were selected. Typically, more than 90% of bursts were selected. The resulting open and closed intervals from single patches were analyzed according to kinetic schemes using an interval-based maximum likelihood algorithm that corrects for missed events and computes error limits of the fitted rate constants (Qin et al., 1996) (http://www.qub.buffalo.edu; QuB suite, State University of New York, Buffalo, NY). A dead time of 30 μs was imposed in all recordings. Probability density functions of open and closed durations were calculated from the fitted rate constants and instrumentation dead time and superimposed on the experimental dwell-time histogram as described by Qin et al. (1996).
Results
Activation of α7-5HT3A HC by Pyrantel and Levamisole
Macroscopic Currents Elicited by Anthelmintic Agents. Levamisole and pyrantel are full agonists of the levamisole-sensitive subtype of nematode muscle AChRs (Martin et al., 1996; Richmond and Jorgensen, 1999), but they are extremely low-efficacy agonists of mammalian muscle AChRs (Rayes et al., 2001, 2004). To understand the structural basis of the differential sensitivity among AChR subtypes to anthelmintics, we studied the activation of α7 AChRs by these agents. We used the high-conductance form of the chimera α7-5HT3A (α7-5HT3A HC), which allows detection of activity at the single-channel level (Rayes et al., 2005).
To evaluate the efficacy of levamisole and pyrantel for activating α7-5HT3A HC, we first recorded macroscopic currents elicited by 1 mM ACh, levamisole, and pyrantel (Fig. 1A). As shown in the figure, levamisole is unable to elicit current responses, similarly to what is observed in mammalian muscle AChRs. In contrast, pyrantel is capable of activating α7-5HT3A receptors. Moreover, the amplitude of the current elicited by pyrantel is similar to that generated by ACh in the same patch. By the relationship between the peak current and pyrantel concentration, an EC50 value of 44 ± 3 μM (nH = 2.6 ± 0.4) is obtained (Fig. 1B). This value is five times lower than that calculated for ACh responses (EC50 ACh = 205 ± 10 μM; nH = 2.43 ± 0.06; Rayes et al., 2005). These results suggest that, in contrast to its weak action at mammalian muscle AChRs, pyrantel is a high-efficacy and potent agonist of α7 AChRs.
Single-Channel Currents Elicited by Anthelmintic Agents. The activation of α7-5HT3A HC by anthelmintics was also evaluated at the single-channel level at a membrane potential of -70 mV (Fig. 2). Single-channel currents are observed at ACh concentrations greater than 30 μM (Rayes et al., 2005). At 500 μM ACh, channel currents appear either as isolated openings flanked by long closings or as bursts of several openings in quick succession (Fig. 2). Open-time histograms can be fitted by three components whose mean durations and relative areas are the following: τ1 = 6.20 ± 0.20 ms (0.50 ± 0.10); τ2 = 0.60 ± 0.30 ms (0.16 ± 0.01); and τ3 = 140 ± 20 μs (0.35 ± 0.08) (Fig. 2). The areas and durations of these components remain almost unaffected at ACh concentrations greater than 30 μM. At 30 μM, the brief events are significantly dominant (Rayes et al., 2005).
Single-channel currents are observed at levamisole concentrations higher than 500 μM (Fig. 2). Openings are extremely infrequent and significantly briefer than those activated by ACh. Open-time histograms can be fitted by a single component whose duration is 50-fold briefer than that of the slowest component observed at 500 μM ACh (127 ± 25 μs). Due to the extremely low number of openings, histograms had to be constructed by adding data from at least three different recordings. In contrast, single-channel currents are readily detected in the presence of pyrantel at concentrations as low as 5 μM. Open-time distributions of α7-5HT3A HC activated by 10 μM pyrantel show two main components, whose durations and relative areas are the following: 2.60 ± 0.55 ms (0.60 ± 0.04) and 240 ± 50 μs (0.40 ± 0.04) (Fig. 2). The slight changes in the open components at 10 μM pyrantel with respect to those observed with ACh are due to the fact that pyrantel also acts as an open-channel blocker (see below). As observed for ACh, at the lowest pyrantel concentration at which channels can be detected (5 μM), the briefest open component is the predominant one (180 ± 20 μs; relative area, 0.60 ± 0.07).
Short bursts of several openings in quick succession are typically observed during activation of α7-5HT3A receptors (Rayes et al., 2005). The duration of bursts recorded at 500 μM ACh (14.8 ± 3.9 ms) and at 10 μM pyrantel (10.6 ± 1.9 ms) are not statistically different (p > 0.05, n = 6).
Kinetic Analysis of α7-5HT3A Channels Activated and Blocked by Levamisole or Pyrantel
Rates for Levamisole and Pyrantel Blockade of ACh-Activated α7-5HT3A Receptors Estimated by Single-Channel Kinetic Analysis. We have shown previously that pyrantel and levamisole act as open-channel blockers of the mammalian muscle AChRs (Rayes et al., 2001, 2004) and that nicotinic agonists act as open-channel blockers of α7-5HT3A HC (Rayes et al., 2005). We therefore examined channel blockade of α7-5HT3A by pyrantel and levamisole. In this receptor, block is produced at the 5HT3A pore. Increasing pyrantel concentration from 10 to 300 μM leads to a clear flickering effect and a concentration-dependent decrease of the mean open time. At 300 μM, open-time histograms are fitted by a single component of 165 ± 30 μs (Fig. 3A). Because openings in the presence of levamisole are very infrequent, we evaluated its action as an open-channel blocker of α7-5HT3A channels activated by 500 μM ACh. As described for pyrantel, increasing levamisole concentration progressively displaces the open-time histogram of ACh-activated channels to briefer durations (Fig. 3B). The relative area of the briefest component of the closed-time histogram increases as a function of the anthelmintic concentration for both pyrantel and levamisole (Fig. 3, A and B). The mean duration of this component, which does not change with concentration, is 42 ± 4 and 30 ± 3 μs for pyrantel and levamisole, respectively.
In summary, channel activity in the presence of pyrantel or levamisole bears all the hallmarks of open-channel block: 1) concentration-dependent decrease in mean open time; 2) concentration-dependent increase in the fractional area of the block component; and 3) constant mean duration of the blocked intervals across all blocker concentrations.
To kinetically describe the block by anthelmintic agents, we performed kinetic analysis of α7-5HT3A activated by ACh and blocked by levamisole or pyrantel. We have shown recently that burst analysis of α7-5HT3A receptor gives information about its activation kinetics (Rayes et al., 2005). Open- and closed-time distributions from bursts activated by ACh are well described by the sum of two exponentials, indicating brief and long openings and closings within each burst. We have also shown that activation within bursts can be well described by a cyclic scheme containing two open and two closed states (Scheme 1).
Here, each open state is connected to two closed states and vice versa (Rayes et al., 2005). To estimate the rates of block of ACh-activated channels by levamisole and pyrantel, we fitted Scheme 2 to the closed and open intervals of bursts recorded in the presence of 500 μM ACh and 100 to 500 μM levamisole, or 1 mM ACh and 1 to 5 μM pyrantel.
Scheme 2 is an extension of Scheme 1, containing two additional blocked states (O1B and O2B) connected to each open state (O1 and O2). Given that levamisole is an extremely weak agonist of this receptor and that the infrequent openings elicited by this drug do not occur in bursts, we can ensure that activation is mediated only by ACh. In contrast, because pyrantel is a full agonist of this receptor, we used only low concentrations of this drug (1-5 μM) to ensure that activation is elicited by ACh. Moreover, the analysis is performed with kinetically homogenous bursts. Therefore, data from either anthelmintic-activated or heteroliganded receptors, which should show different kinetics from those of ACh-activated receptors (Steinbach and Chen, 1995; Akk et al., 2005), would be discarded during the selection of bursts. Given that channel block is not observed at 1 mM ACh (Rayes et al., 2005), we can ensure that block is only produced by levamisole or pyrantel.
We constrained the rate constants between closed and open states to their previously determined values for ACh (Rayes et al., 2005; Table 1). We assumed that blocking and unblocking rates are the same for the two open states (O1 and O2). We simultaneously fitted the scheme to data obtained at multiple levamisole or pyrantel concentrations (Table 1). The probability density functions resulting from the kinetic analysis adequately fitted the experimental histograms (Fig. 3C), indicating that Scheme 2 provides a good description of the anthelmintic block of ACh-activated channels. The resulting estimates for blocking rates are the following: k+b = 10 ± 0.30 and 49 ± 2.40 μM-1 s-1 for levamisole and pyrantel, respectively; k-b = 32,000 ± 800 and 15,400 ± 640 s-1 for levamisole and pyrantel, respectively. Considering the estimated rates for k+b and k-b, the Kd values (k-b/k+b) for the blocking process are 3.2 mM and 310 μM for levamisole and pyrantel, respectively. On the basis of the Kd value for pyrantel, a decrease in the peak current should be observed at 1 mM pyrantel (Fig. 1). However, only at 2 mM did we detect a 36 ± 5% reduction of the peak current. Thus, a quantitative correlation between the Kd calculated from the single channel data and the reduction of the macroscopic currents is not observed. Discrepancies between results obtained from cell-attached and whole-cell recordings have been reported previously (Akk et al., 2004, 2005). The reasons for these discrepancies are not understood. The possible explanations include the patch excision, exposure of the receptor to different monovalent inorganic cations, statistical and biological variability, and experimental artifacts (Akk et al., 2004, 2005).
Kinetic Analysis of Pyrantel-Activated α7-5HT3A Receptors. To estimate the rate constants for activation of α7-5HT3A by pyrantel, and after having assessed that Scheme 2 can well describe activation and blockade, we fitted the open and closed times of channels activated by 10 and 50 μM pyrantel to this scheme. In this analysis, we constrained the blocking and unblocking rates to the values determined above for pyrantel block of ACh-activated α7-5HT3A receptors.
The estimated values for pyrantel activation rates are detailed in Table 1. The main difference with respect to ACh activation is observed in α12. The fact that this rate is very slow when pyrantel is the agonist (∼40,000-fold slower than for ACh) indicates that long closings within a burst occur very rarely. In close correspondence, a thorough inspection of the recordings shows very few long closings within bursts, in contrast to what is observed when ACh is the agonist (Rayes et al., 2005).
As a control, we constrained the rates for the blocking process to those determined as described above and fitted Scheme 2 to data obtained in the presence of 1 mM ACh and 5 μM pyrantel. No changes in ACh-activation rates were observed with respect to those estimated in the absence of pyrantel. We can therefore ensure that the differences in the estimated rates for pyrantel- and ACh activation do not arise from fitting the data to different schemes either containing (Scheme 2) or lacking (Scheme 1) blocked states.
Residues Involved in the Activation of α7-5HT3A by Pyrantel
Because pyrantel is a weak agonist of the mammalian muscle AChR but a potent agonist of α7-5HT3A, we studied the structural basis of such a difference by mutating residues that are probably involved in agonist binding. Residues of the principal face that have been shown to be involved in agonist binding (Sine, 2002) are highly conserved between α7 and α1 subunits. Less conservation is found at the complementary face of the binding site (Fig. 4). We therefore replaced residues of this side of α7-5HT3A by the homologous residues in the ϵ subunit, which contributes to one of the binding sites in the muscle AChR (Brejc et al., 2001). We constructed the following mutant receptors: α7-Q57G, α7-N111Y, α7-Q117T, and α7-G167A and evaluated channel activity elicited by pyrantel (Fig. 5).
The minimum concentrations of both ACh (30 μM) and pyrantel (5 μM) at which channel currents can be detected and the channel properties are not affected by the α7-G167A mutation (Fig. 5). In the α7-N111Y and α7-Q117T mutants, the minimum concentration to elicit single-channel currents increases for both ACh and pyrantel with respect to the control chimera. In these mutants, channels appear at approximately 2- and 3-fold higher concentrations of ACh and pyrantel, respectively, than in the control α7-5HT3A (60 μM for ACh and 15 μM for pyrantel). The mean open time, mean closed time and mean burst duration do not change in the mutants with respect to control α7-5HT3A for both agonists (Fig. 5). It can therefore be concluded that α7-N111Y and α7-Q117T mutations do not affect gating kinetics; they only produce a slight and quantitatively similar decrease in the sensitivity for activation by both agonists. Thus, these residues do not contribute significantly to the differential activation between pyrantel and ACh among AChR subtypes.
When Gln57 is replaced by glycine, which is the equivalent residue at the ϵ subunit, the potency of pyrantel to activate α7-5HT3A decreases significantly. Single-channel openings can be observed only at concentrations greater than 50 μM (instead of 5 μM). The duration and relative areas of the slowest component of the open time distribution corresponding to the entire recording are similar to those of the control chimera at each pyrantel concentration [810 ± 30 μs (0.64 ± 0.01) and 800 ± 90 μs (0.73 ± 0.02), 530 ± 30 μs (0.66 ± 0.05) and 560 ± 10 μs (0.70 ± 0.05), and 100 ± 10 μs (1) and 120 ± 20 μs (1) for mutant and control at 50, 100, and 500 μM pyrantel, respectively].
In contrast, the duration of bursts of openings decreases 3-fold in the mutant receptor (11.63 ± 2.02 and 4.10 ± 0.41 ms for the control and mutant α7-5HT3A, respectively). Such reduction is accompanied by a reduction in the number of openings per burst (20.5 ± 9 and 6.5 ± 0.2 events per burst for the control and Q57G mutant chimera activated by 50 μM pyrantel, respectively). To discard that the briefer bursts observed in the mutant were isolated openings blocked by pyrantel, we recorded channels at +70 mV. At this membrane potential, channel block is significantly reduced. However, openings still occur in bursts (4.8 ± 0.7 openings per burst).
It is interesting to note that when the mutant Q57G chimera is activated by ACh, no differences in the minimum concentration required to detect single events and in the mean open time, mean closed time, and mean burst duration are observed with respect to those of the control chimera. At 500 μM ACh, the mean durations and relative areas of the open components are the following: τ1 = 6.50 ± 0.80 ms (0.25 ± 0.01); τ2 = 0.80 ± 0.50 ms (0.16 ± 0.05); and τ3 = 90 ± 10 μs (0.56 ± 0.02); and the mean burst duration is 14.3 ± 6.4 ms. These results show that the mutation Q57G selectively alters activation by pyrantel.
We examined activation of Q57G mutant by levamisole, but channels were detected only in 10 to 20% of the patches of transfected cells. The few channels detected were similar to those of the control chimera. Given the low activity in both the control and mutant α7-5HT3A receptors activated by levamisole, it is not possible to determine at this step whether this position also affects levamisole activation.
To determine the overall changes in pyrantel activation produced by the mutation Q57G, we measured macroscopic currents (Fig. 6). The EC50 for pyrantel, calculated by the dose-response curve, is 2-fold greater than that of the control chimera (77 ± 4 μM instead of 44 ± 3 μM). The Hill coefficient decreases from 2.60 ± 0.4 (control chimera) to 1.83 ± 0.02 in the mutant chimera (p < 0.05, n = 4), thus indicating an impairment of channel gating (Colquhoun, 1998). The peak current induced by a saturating concentration of pyrantel (500 μM) is 61 ± 4% of that elicited by 500 μM ACh, showing that the maximum response is also affected. A similar decrease in the peak current is observed at +100 mV (Fig. 6). At this positive potential, open-channel block is almost undetected (Neher and Steinbach, 1978; Rayes et al., 2004). This observation together with the fact that a reduction in peak current by channel block is not observed in the control chimera at 500 μM pyrantel (see above) indicate that the reduction in the maximum response is originated from a reduced efficacy and not from channel block.
In contrast, neither the EC50 value (220 ± 20 μM) nor the Hill coefficient (2.34 ± 0.4) changed for ACh. As in the control chimera, levamisole is unable to elicit current responses in the mutant Q57G chimera.
In conclusion, the mutation Q57G does not affect activation by ACh, but it significantly decreases the potency and efficacy of pyrantel activation. This position is differentially conserved among subunits. It is glycine in mammalian ϵ subunits, glutamic acid in γ subunits, and glutamic acid or aspartic acid in δ subunits.
If a glutamine is found at this position, as in α7 receptors, pyrantel behaves as a full agonist. If a glycine occupies position 57, as in the muscle ϵ subunit, pyrantel behaves as a partial agonist. These results indicate that position 57, located at the complementary face of the AChR binding site, is a main determinant of the different selectivity of pyrantel between muscle and α7 AChRs.
Kinetics of Activation of Q57G Mutant Receptor by Pyrantel or ACh
To analyze the activation kinetics of the Q57G mutant, we first fitted Scheme 1 to the open and closed intervals of bursts recorded at 50 and 500 μM ACh. As expected from the single-channel and macroscopic current data, no changes are observed in the estimated rates for ACh activation between the control and the mutant chimera (Table 1; Fig. 7). In contrast, the kinetic analysis for pyrantel activation of the mutant receptor reveals important differences in the estimated rate constants with respect to the control receptor (Table 1; Fig. 7). For this analysis, we fitted Scheme 2 to data obtained in the presence of 60 and 70 μM pyrantel. These concentrations, which are higher than the minimum concentration at which channels can be detected (50 μM), allow us to maintain channel blockade as low as possible. As shown in Fig. 7, open-time histograms constructed with the selected bursts show a main component of 690 ± 90 μs with a relative area of 0.95 ± 0.02. The main changes in activation rates due to the mutation Q57G are observed in β1 and α1, which are ∼4000- and ∼180-fold, respectively, lower than those of the control chimera (Table 1). The kinetic analysis reveals that in this mutant activated by pyrantel, the cyclic scheme could be simplified to a linear one, containing one open state (O2), one closed state (C1), and the blocked state. Thus, in contrast to the nonmutant receptor that is able to oscillate between two different open states, the mutant receptor opens mainly to a single state. Assuming that the burst finishes by desensitization from the longest duration open state (O2), the absence of a second open component may explain the reduction in the number of openings per burst. An increase in the dissociation rate, which cannot be measured by kinetic analysis in the chimeric receptor (Rayes et al., 2005), may also contribute to the decrease in the number of openings per burst. In summary, the mutation Q57G selectively impairs activation by pyrantel.
Activation by Anthelmintic Agents of α7-5HT3A Receptors Carrying a Mutation at Gly153, Located at the Principal Face of the Binding Site
We have shown previously that mutation of the highly conserved Gly153 of mammalian α1 subunit to a glutamic acid, which is found at equivalent positions in several nematode α subunits, increases the efficacy of levamisole and pyrantel to activate mammalian muscle AChRs (Rayes et al., 2004). We therefore determined how this position is involved in the activation of α7 by pyrantel and levamisole. As shown in Fig. 8, single channels of the mutant α7-G153E receptor activated by pyrantel and levamisole can be detected at concentrations 5-fold lower than those required to activate the control receptor. High channel activity is observed at 100 μM levamisole in the α7-G153E mutant, and open-time histograms show two components (τ1, 795 ± 120 μs; τ2, 154 ± 40 μs). Channels are observed at pyrantel concentrations as low as 1 μMinthe α7-G153E mutant. Open-time distributions at 1 μM pyrantel are fitted by three components whose durations are 5.10 ± 0.71 ms, 1.04 ± 0.21 ms, and 105 ± 41 μs. For ACh, the minimum concentration for channel detection decreases 3-fold (10 μM). Open-time histograms are similar to those of the control chimera (τ1, 6.80 ± 0.20 ms; τ2, 0.40 ± 0.20 ms; τ3, 110 ± 20 μs).
Thus, the mutation α7-G153E affects activation by both ACh and anthelmintics, albeit changes are more pronounced for anthelmintics. In contrast, Gln57 does not contribute to ACh activation, but it significantly affects activation by pyrantel. Thus, this residue plays an important role in the selectivity for pyrantel.
To further explore the contributions of Gln57 and Gly153 to pyrantel and levamisole activation, we analyzed the double mutant Q57G-G153E. In the double-mutant receptors activated by both levamisole and ACh, the minimum concentration required to detect single channels and the channel properties are the same as those of the single G153E mutant. Channels are detected at 10 μM ACh or 100 μM levamisole. For 10 μM ACh, neither the open-time distribution (τ1, 6.80 ± 0.40 ms; τ2, 0.50 ± 0.20 ms; τ3, 74 ± 11 μs) nor the mean burst duration (13.3 ± 1.5 ms) is affected. At 100 μM levamisole, open-time distributions show two components whose durations are similar to those of the G153E single mutant (τ1, 790 ± 14 μs; τ2, 90 ± 4 μs). These results demonstrate that Gln57 does not contribute to the activation by levamisole and ACh.
It is interesting that the double mutant Q57G-G153E activated by pyrantel shows an intermediate behavior between the single Q57G, which shows decreased efficacy, and G153E, which shows increased efficacy. The minimum concentration required to detect single channels is 10 μM, which is intermediate between that of G153E (1 μM) and Q57G (50 μM) mutant receptors. As described for both single mutations, no changes in the open-time components are observed with respect to the control chimera.
Taken together, our experiments demonstrate that position 153 located at the principal face is involved in the activation by ACh, levamisole, and pyrantel whereas position 57 located at the complementary face is involved only in pyrantel activation.
Discussion
Levamisole and pyrantel are anthelmintic agents that act by disrupting neuromuscular transmission in nematodes. They are full agonists of the levamisole-sensitive AChR. In contrast, they are partial agonists of mammalian muscle AChRs (Rayes et al., 2001, 2004). Position 153 of the α subunit is involved in this differential selectivity (Rayes et al., 2004). To further explore determinants of such selectivity, we investigated the actions of these agents on the neuronal α7. We took advantage of the high-conductance form of α7-5HT3A that allows the detection of single-channel events (Rayes et al., 2005). Both single-channel and macroscopic current recordings show that levamisole produces a very weak activation. It is interesting that pyrantel is a more potent agonist than ACh for this receptor. By replacing residues located at the complementary face of the binding site in α7 by the equivalent residues in the muscle ϵ subunit, we identified Gln57 as one of the residues involved in this differential sensitivity. The same residue has been shown to have a critical role in the interaction of α7 with neonicotinoid insecticides, and it may be involved in the selective actions of these compounds at insect AChRs (Shimomura et al., 2002). Thus, this position seems to have a fundamental role in drug selectivity.
Differential Actions of Pyrantel and Levamisole at α7 AChRs. The behavior of levamisole at α7-5HT3A receptors is similar to that at mammalian muscle AChRs: it acts as an extremely low-efficacy agonist. The weak activation is revealed by the low frequency of opening events at all levamisole concentrations and the absence of macroscopic currents elicited by levamisole.
This is the first report showing that pyrantel is a full agonist of α7 AChRs. The following lines of experimental evidence support this statement: 1) the current amplitudes are similar to those elicited by ACh; 2) the dose-response curve is displaced to lower agonist concentrations, indicating that it has indeed a higher potency than ACh; 3) the minimum concentration that allows detection of single-channel events is 5-fold lower for pyrantel than for ACh; and 4) single-channel kinetic analysis shows similar rates for activation of α7-5HT3A by pyrantel and ACh.
Kinetics of Activation of α7-5HT3A by Pyrantel. Channel events of α7-5HT3A activated by ACh appear as single openings flanked by long closings or as bursts of several openings in quick succession (Rayes et al., 2005). Bursts elicited by a saturating concentration of ACh contain brief and long openings and closings, and a cyclic scheme containing two open and two closed states describes activation of receptors with the optimal number of agonist sites occupied (Rayes et al., 2005). The same cyclic scheme with two additional blocked states well describes the activation by pyrantel, which shows rate constants similar to those of ACh activation. The analysis indicates that once fully occupied by pyrantel, the receptor opens efficiently and closes and reopens several times, undergoing transitions mainly through C1, O1, and O2 states.
Structural Basis for the High Selectivity of Pyrantel for α7. The differences in the potency of pyrantel among AChR subtypes may result from differences in the affinity or in their ability to open the channel once bound (efficacy). Making the distinction between effects on binding and effects on conformational change is not easy (Colquhoun, 1998). Therefore, residues involved in the selectivity of pyrantel for α7 could be located anywhere in the receptor. We here focused on residues located at the complementary face of the binding site, which show poor conservation between α7 and the muscle ϵ subunit. Mutations at residues α7-Asn111 and α7-Gln117 slightly decrease the potency for activation by ACh and pyrantel in a quantitatively similar manner, indicating that they do not govern pyrantel selectivity. The equivalent residues in the Torpedo AChR, γTyr111 and γTyr117, were reported to interact with d-tubocurarine (Chiara et al., 1999; Willcockson et al., 2002). Thus, it is possible that mutations at these residues affect only the affinity for agonists. The mutation G167A does not affect the sensitivity to either ACh or pyrantel, but it has been shown previously to be involved in the sensitivity of α7 to imidacloprid (Matsuda et al., 2000).
The mutation Q57G selectively decreases the efficacy for activation by pyrantel without affecting activation by ACh. At the macroscopic current level, such a decrease was revealed by the decrease in the maximum peak current elicited by pyrantel, an increase in the EC50 value, and a decrease in the Hill coefficient. In contrast, dose-response curves for ACh do not change in the mutant with respect to the control chimera. The reduction in the Hill slope indicates an impairment of gating (Colquhoun, 1998), and it would not be expected if the residue were affecting only binding affinity of the resting state (Spier and Lummis, 2000). At the single-channel level, no changes are observed in Q57G channels activated by ACh with respect to the control ones. In contrast, in the presence of pyrantel, a dramatic reduction of the burst duration is observed. Such reduction may be due to changes in the agonist-dissociation rate, in channel gating, and/or desensitization. Because of the absence of closed components sensitive to agonist concentration, it is not possible to evaluate the binding steps in our kinetic analysis of α7-5HT3A (Rayes et al., 2005). Therefore, we cannot determine whether the mutation α7-Q57G affects agonist dissociation. Macroscopic currents do not show significant changes in the decay rates, thus discarding important changes in desensitization. Nevertheless, the kinetic analysis reveals profound changes in channel gating. The cyclic scheme can be simplified to a linear one containing only one open (O2), one closed (C1), and a blocked state. It is interesting to note that in the control α7-5HT3A chimera, the brief openings within bursts (O1) correspond to approximately 30% of the total openings, and that in the mutant chimera, the peak current shows a decrease of approximately 40%. Given that C1 is the most likely state linking binding to channel opening steps (Rayes et al., 2005), it is therefore possible that the reduction in the peak current is mainly due to the fact that the C1-to-O1 transition has a low probability to occur (β1 decreases 4000-fold). The remaining current may correspond to the C1-to-O2 transition, which is not significantly affected in the Q57G mutant.
Although the molecular interactions that result in partial agonism are not yet understood, we here show that a glycine at position 57, which is found in the muscle AChR, converts pyrantel into a partial agonist by affecting the way in which it couples binding to channel gating. The same position has been shown to be involved in the sensitivity of α7 to neonicotinoid insecticides (Shimomura et al., 2002).
The mutation α1-G153E increases the efficacy of levamisole as an agonist of mammalian muscle AChRs (Rayes et al., 2004). A glutamic acid is found in α subunits of nematode muscle, such as UNC-38, α1A. suum, TAR1 T. colubriformis, and it could be involved in the high efficacy of levamisole at nematode AChRs. This residue is located at loop B, which contributes to the principal face of the ACh binding pocket in the α1 subunit (Brejc et al., 2001). Our results show that, as in the muscle AChR, the α7-G153E mutation increases the potency for activation by ACh and anthelmintics and that the increase is selectively higher for the latter drugs than for ACh. These results agree with those reported by Grutter et al. (2003), showing that the mutation G153K increases the potency and the affinity for different ligands in a selective manner.
The experiments with the Q57G-G153E double mutant enlighten the contribution of both residues to anthelmintic activation. The results reveal that activation by ACh or levamisole is affected by the residue at position 153 but not by that located at position 57. In contrast, both positions are similarly involved in pyrantel activation. Extrapolated to the nematode AChR, these findings suggest that mutations at specific residues, such as Gln57, may affect the sensitivity to a specific anthelmintic, such as pyrantel. This could be a possible explanation for the observation that under pyrantel treatment of pigs infested with Oesophagostomum species, the nematodes acquire resistance to this anthelmintic but not to levamisole (Bjorn et al., 1990).
Our results reveal that the maximum sensitivity to pyrantel would be achieved when the AChR contains both a glutamic acid at position 153 (principal face) and a glutamine at position 57 (complementary face). The minimum sensitivity would be achieved in the receptor containing glycine at both positions (Gly57 and Gly153), which corresponds to the mammalian muscle AChR. It remains to be determined whether there is either a muscle levamisole-sensitive AChR, which is also sensitive to pyrantel, composed of UNC-38, UNC-63, UNC-29, and LEV-1, or more than one receptor, each made from various combinations of subunits (Culetto et al., 2004). We compared the sequences of C. elegans α and non-α subunits that may form levamisole-sensitive AChRs. On the basis of our results, we speculate that Glu153 is provided by UNC-38 and Gln57 by UNC-63, and that this combination leads to receptors highly sensitive to pyrantel.
Footnotes
-
This work was supported by grants from CONICET, Universidad Nacional del Sur, FONCYT and a fellowship from John Simon Guggenheim Memorial Foundation (to C.B.).
-
ABBREVIATIONS: AChR, nicotinic acetylcholine receptor; 5-HT3A, 5-hydroxytryptamine-3A receptor; ACh, acetylcholine; HC, high conductance.
- Received May 4, 2006.
- Accepted July 5, 2006.
- The American Society for Pharmacology and Experimental Therapeutics