Abstract
Acid-sensing ion channel 3 (ASIC3) is a proton-gated Na+ channel with important roles in pain. ASIC3 quickly desensitizes in less than a second, limiting its capacity to sense sustained acidosis during pain. RFamide neuropeptides are modulators of ASIC3 that slow its desensitization and induce a variable sustained current. The molecular mechanism of slowed desensitization and the RFamide binding site on ASIC3 are unknown. RPRFamide, a RFamide from the venom of a cone snail, has a comparatively high affinity for ASIC3 and strongly slows its desensitization. Here we show that covalent binding of a UV-sensitive RPRFamide variant to ASIC3 prevents desensitization, suggesting that RPRFamide has to unbind from ASIC3 before it can desensitize. Moreover, we show by in silico docking to a homology model of ASIC3 that a cavity in the lower palm domain, which is also known as the nonproton ligand–sensing domain, is a potential binding site of RPRFamide. Finally, using extensive mutagenesis of residues lining the nonproton ligand–sensing domain, we confirm that this domain is essential for RPRFamide modulation of ASIC3. As comparative analysis of ASIC crystal structures in the open and in the desensitized conformation suggests that the lower palm domain contracts during desensitization, our results collectively suggest that RPRFamide, and probably also other RFamide neuropeptides, bind to the nonproton ligand–sensing domain to stabilize the open conformation of ASIC3.
Introduction
Acid-sensing ion channels (ASICs) are Na+ channels gated by extracellular protons. They assemble as trimers (Bartoi et al., 2014). Each subunit has a topology with two transmembrane domains and a large extracellular domain (ECD), which is composed of 7 α-helices and 12 β-strands that fold into a structure that has been compared with a clenched hand (Jasti et al., 2007). In this image, the two transmembrane domains build the forearm and the ECD is composed of a palm domain, a finger, a knuckle, and a thumb domain that surround the central β-ball domain. A pocket in the ECD containing several conserved acidic residues has been implicated in binding protons (Jasti et al., 2007) and gating modifier toxins (Baconguis and Gouaux, 2012; Dawson et al., 2012).
ASIC3 is an isoform that is predominantly expressed in the peripheral nervous system, where it contributes to the detection of painful acidosis (Deval and Lingueglia, 2015; Sluka and Gregory, 2015). It is expressed in primary muscle afferents and mediates muscle pain in ischemic-like conditions (Sluka et al., 2003; Molliver et al., 2005; Ross et al., 2016), for example during muscle fatigue. Muscle incision also decreases pH, and ASIC3 also mediates postincisional (postoperative) pain (Deval et al., 2011). Its expression is particularly high in afferents innervating the heart muscle, where it probably mediates the pain during myocardial ischemia (angina) (Benson et al., 1999; Sutherland et al., 2001). Moreover, in rodents it contributes to cutaneous acid pain (Deval et al., 2008), to hyperalgesia associated with neuropathic pain (Omori et al., 2008), to acute arthritic pain (Ikeuchi et al., 2009), and to itch sensation in response to acid and pruritogens (Peng et al., 2015). Thus, ASIC3 has a well established role in the detection of pain associated with sustained acidosis. In addition, ASIC3 is expressed in cell models of glioblastoma (Tian et al., 2017), where it might sense the sustained acidic microenvironment of a brain tumor.
Upon rapid application of protons, however, ASIC3 completely desensitizes, with a time constant of 0.3 seconds (Sutherland et al., 2001). During a slight but prolonged acidosis, it desensitizes without apparent opening (Gründer and Pusch, 2015), clearly limiting its capacity to sense sustained acidosis. Although it has been shown that the overlap of steady-state desensitization and activation curves of ASIC3 leads to a small window current in a tiny window of pH values that roughly correspond to the pH values reached during myocardial ischemia (Yagi et al., 2006), it is probable that different modulators interact with ASIC3 to either extend the window current (Deval et al., 2008; Marra et al., 2016) or to induce sustained currents.
One group of such modulators are RFamide neuropeptides (Vick and Askwith, 2015). The prototypical RFamide is FMRFamide, an important neurotransmitter in different invertebrates. The main effect of FMRFamide is a slowing of desensitization (prolonging opening) of ASICs. In addition, it induces a variable sustained current (Askwith et al., 2000). Mammalian RFamides, neuropeptide FF (FLFQPQRFamide) and neuropeptide AF (AGEGLNSQFWSLAAPQRFamide), have similar effects that are smaller, however (Askwith et al., 2000). Although RFamides bind to the open conformation of ASICs (Chen et al., 2006b), the on-rate to their binding site is relatively slow, and therefore binding to the closed conformation is usually necessary for an efficient modulation (Askwith et al., 2000; Chen et al., 2006b). Neither the mechanism by which RFamides slow desensitization and induce a sustained current nor their binding site on ASICs are known.
Recently, RPRFamide (RPRFa), a conorfamide that has a comparatively high affinity to ASIC3 and strongly slows its desensitization, has been isolated from the venom of the cone snail Conus textile (Reimers et al., 2017). RPRFa increases the excitability of sensory neurons expressing ASIC3, and injection of RPRFa into the gastrocnemius muscle of mice enhances muscle pain in response to acidic stimuli (Reimers et al., 2017). Thus, RPRFa is an ideal tool to study the modulation of ASIC3 by RFamides. Here, we reveal that RPRFa has to unbind from ASIC3 before the channel can desensitize, explaining the slowed desensitization. Moreover, we show that a cavity in the lower ECD of ASIC3, which is surrounded by the base of the three palm domains, the so-called nonproton ligand–sensing domain (NPLSD), is important for the modulation by RPRFa. Since the lower palm domain contracts during desensitization (Baconguis and Gouaux, 2012; Baconguis et al., 2014), our results provide a mechanistic explanation for the functional effects of RFamides on ASIC3 and strongly suggest that the NPLSD is the RFamide binding site on ASIC3.
Materials and Methods
Peptides and Chemicals
RPRFa and RPR[4-azido-l-phenylalanine (azF)]a were synthesized by ProteoGenix (Schiltigheim, France) with a purity of ≥97%. All other chemicals were purchased from Carl Roth (Karlsruhe, Germany), MilliporeSigma (St. Louis, MO) or VWR (Radnor, PA).
Site-Directed Mutagenesis and RNA Synthesis
We used rat ASIC3 cDNA contained in the expression vector pRSSP (Chen et al., 2006b). Mutants of ASIC3 were generated by site-directed mutagenesis following the Stratagene QuikChange protocol; polymerase chain reaction was performed with KAPA HiFi Polymerase (Kapa Biosystems, Wilmington, MA). The sequences of mutated cDNAs was entirely verified by sequencing (Eurofins Genomics, Ebersberg, Germany).
After linearization of plasmids with MluI (New England Biolabs, Ipswich, MA), capped cRNA was synthesized using the mMessage mMachine SP6 transcription kit (Thermo Fischer Scientific, Waltham, MA). Finally, cRNA was purified using RNA Clean & Concentrator (Zymo Research, Irvine, CA). Concentrations were determined with a Nanodrop 2000 spectrophotometer (Thermo Fischer Scientific). RNA was stored at −80°C.
Preparation and Injection of Oocytes
Animal care and experiments were conducted according to protocols approved by the state office for nature, environment, and consumer protection (LANUV) of the state North Rhine-Westphalia (NRW). Surgical removal of ovaries of female Xenopus laevis frogs was described previously (Springauf and Gründer, 2010). Oocytes were separated and the follicular membrane was removed by mechanical dissection followed by enzymatic digestion for a period of 2 hours with collagenase type 2 (Worthington Biochemical Corporation, Lakewood, NJ). We injected oocytes of stages V and VI with 3 ng capped cRNA of ASIC3 wild-type or ASIC3 mutants. For ASIC3 mutated at residues E239 or E79 and for substitutions by cysteine, 8–16 ng cRNA were injected. Injected oocytes were incubated for 24–72 hours at 19°C in oocyte Ringer solution (in millimolars: 82.5 NaCl, 2.5 KCl, 1.0 Na2HPO4, 5.0 HEPES, 1.0 MgCl2, 1.0 CaCl2, 0.5 g/l polyvinylpyrrolidone, 1000 IU/l penicillin and 10 mg/l streptomycin; pH was adjusted to pH 7.3).
Electrophysiology
For most experiments, we used the Screening Tool (npi electronic, Tamm, Germany) as measurement setup, which provides a small recording chamber combined with a computer-driven dispensing robot, allowing fast and programmable solution exchange (Baburin et al., 2006). Two-electrode voltage clamp recordings of Xenopus oocytes were conducted at ambient temperature (20–23°C). A Turbo TEC-03X amplifier (npi electronic) was used to establish voltage clamp conditions at −70 mV in all experiments. Micropipettes filled with 3M KCl and a resistance of 0.5–2 MΩ were used to gain electrical access to oocytes. The current signal was filtered at 20 Hz and digitized with a sampling rate of 500 Hz using the software Cellworks (npi electronic).
For the experiments reported in Fig. 1 and in Fig. 2C-E, a setup with a larger recording chamber and a computer-controlled valve-driven solution exchange system was used. This setup allowed shorter intervals between solution applications (as short as 2 seconds) (Madeja et al., 1995). For these experiments, currents were digitized with a sampling rate of 200 Hz.
Bath solutions for Two-electrode voltage clamp recordings recordings contained (in millimolars) 140 NaCl, 1.8 CaCl2, 1.0 MgCl2, and 10 HEPES (for pH ≥ 6.6) or 10 MES (for pH < 6.6); pH was adjusted with NaOH. If not stated otherwise, pH 7.4 was used as conditioning solution. Compounds were added to the bath solution as indicated.
The desensitized state of ASIC3 with RPRFa bound is not strongly populated. (A) Simplified kinetic model of the three main states of ASIC3 with or without RPRFa bound. τ1 and τ2 are the time constants of current decay in the absence or presence of RPRFa, respectively. The transitions that describe current decay are highlighted in bold. State DR is highlighted, as its population is unsure. (B) Left, representative current trace of an ASIC3-expressing oocyte when activated by pH 6.0 (black bars) and treated with 100 μM RPRFa (red bars). Peaks are marked by roman numerals. Right, mean I2.5 seconds/Ipeak ratios for the pH-activated currents. Mean amplitude of peak i was −17.4 ± 3.4 μA; n = 12. (C) Left, representative current trace of an ASIC3-expressing oocyte when activated by pH 6.0 (black bars) and treated with 30 μM RPRFa (red bars). Peaks are marked by roman numerals. The dotted red line represents a fit of the current decline of peak i. Current decline of peak ii followed this fit. The trace shown is representative for nine similar experiments. Right, mean I2.5 seconds/Ipeak ratios. Mean amplitude of peak i was −14.0 ± 5.1 μA; n = 9. Error bars represent S.E.M.; n.s., not significant; **P < 0.01; ***P < 0.001 (Student’s paired t test).
The photo-reactive peptide RPR[azF]a abolishes desensitization when crosslinked to ASIC3. (A) Left, representative current trace illustrating that increasing concentrations of RPR[azF]a (red bars) slow the desensitization and increase current amplitude of ASIC3 activated by pH 6.3 (black bars). Right, concentration-response curve. The response was quantified as the I5 seconds/Ipeak ratio in measurements like the one shown on the left. n = 7 for 1, 10 and 100 μM and n = 5 for 0.1, 3 and 30 μM. (B) Scheme illustrating that ASIC3-expressing oocytes were incubated in RPR[azF]a and irradiated by UV light. Irradiation time and RPR[azF]a concentration were 1 minute and 30 μM, 6 minutes and 100 μM, or 10 minutes and 30 μM, respectively. When applied for 10 minutes, half of the RPR[azF]a solution was renewed every 2 minutes. The control condition was incubated in 100 μM RPRFa and irradiated for 6 minutes. (C) Representative ASIC3-expressing oocytes, which had been photo-labeled, were activated by pH 6.0 (black bar) in the absence of any peptide. (D) Mean I25 seconds/Ipeak ratios of pH 6.0 currents of UV-irradiated ASIC3-expressing oocytes. n = 5–10. Error bars represent S.E.M.; ***P < 0.001. Mean peak currents were in the range of −7.0 ± 3.1 to −17.0 ± 2.5 μA. (E) Representative current trace of an ASIC3-expressing oocyte, which had been photo-labeled with 30 μM RPR[azF]a for 4 minutes. The first activation with pH 6.0 was in the absence of peptide and the second activation was after preapplication of 30 μM RPR[azF]a. Six similar measurements were recorded in the same week.
Photo Affinity Labeling of ASIC3 by RPR[azF]a
ASIC3-expressing oocytes were incubated in oocyte Ringer medium containing either 30–100 μM RPR[azF]a or 100 μM RPRFa and were irradiated by UV light. Irradiation was for 1–10 minutes with UV light of 302 nm in a Gel Doc XR UV-transilluminator (Bio-Rad, Hercules, CA). Electrophysiological recordings of photo-crosslinked ASIC3-expressing oocytes were conducted within 2 hours after photo-crosslinking.
Data Analysis
For analysis of current recordings, we used Cellworks Reader (version 6.2.2; npi Electronic) and IgorPro (version 5.0.3.0; WaveMetrics, Lake Oswego, OR). For kinetic analysis of the current decay phase, it was fit with a double exponential function:(1)I(t) represents the current at time point t, y0 a sustained current component, A the desensitizing current component, and a1 and a2 the two fractions of current components desensitizing with fast τ1 and slower τ2, respectively; a1 + a2 = 1. The underlying assumption was that currents of ASIC3 with no RPRFa bound decay with τ1 and currents of ASIC3 with RPRFa bound with τ2 (Reimers et al., 2017).
To determine τ after photo-crosslinking with 100 μM RPR[azF]a (Fig. 1), we used a single exponential function. To analyze desensitization of mutants of the NPLSD (Fig. 6), we first obtained τ1 for each mutant by fitting with a single exponential function (a2 = 0) current decay in the absence of RPRFa. We then fitted current decay in the presence of RPRFa with a double exponential function, with τ1 fixed to the value obtained in the absence of RPRFa, yielding mutant-specific a2 and τ2. All exponential fits were visually checked and inappropriate fits (e.g., with a negative time constant) were rejected.
Statistical analysis was performed in Excel 2016 (Microsoft, Redmond, WA) and R 3.4.2 (R Core Team, Vienna, Austria). Results are reported as the mean ± S.E.M. Student’s unpaired t test was used to determine P values, with the exception of experiments reported in Fig. 1, for which Student’s paired t test was used. Significance levels were defined as follows: *P < 0.05; **P < 0.01; ***P < 0.001. In case of multiple testing on the same data set, the Bonferroni correction was used by multiplying the P value by the number of tests. Graphs and current traces were plotted with IgorPro and arranged with Adobe Illustrator CS6. All experiments were conducted with oocytes of at least two different frogs, with the exception of photo-labeling experiments, which were conducted with oocytes of one frog.
Modeling of Rat ASIC3 and RPRFa
We created a knowledge-based homology model of rat ASIC3 (UniProtKB: O35240, 533 aa) in both its open and desensitized state, on the basis of the corresponding ion channel structures of chicken ASIC1a (cASIC1). The sequences were aligned with ClustalW, revealing a sequence identity of 44%. In the open cASIC1 structure, the channel is in complex with the snake toxin MitTx [Protein Data Bank (PDB) entry code: 4NTW; resolution of 2.07 Å] (Baconguis et al., 2014). For the desensitized receptor state, we used as a template the desensitized structure of cASIC1 (PDB entry code: 2QTS; resolution of 1.9 Å) (Jasti et al., 2007).
The full biologic assemblies were retrieved from the PDB and processed with Schrödinger software (Maestro Version 11.4.011; Small-Molecule Drug Discovery Suite 2017-4; Schrödinger, LLC, New York, NY), using default settings unless otherwise noted. All entities not belonging to the receptor trimer were removed, and missing side chains/loops were modeled. With PROPKA (propka.org) we adjusted the pH of the receptors to 6 to account for the pH experienced by an open channel. Homology models are provided in Supplemental Figs. 1 and 2 in Protein Data Bank format (pdb). The tetrapeptide RPRF was built from the sequence in an extended conformation and a C-terminal amide group was manually introduced yielding RPRFa.
Binding Site Assessment, Docking, and Scoring of RPRFa Binding Poses
Binding Site.
To characterize potential binding sites on the homology models corresponding to the open and the desensitized states, we used Schrödinger SiteMap (Halgren, 2009). Site candidates were ranked on the basis of SiteScore (see Supplemental Table 1), which considers the presence of hydrophilic and -phobic features and pocket size. The required amount of enclosure and threshold for van der Waals interactions was decreased to detect also shallower binding sites known to be relevant in protein-protein interactions. SiteScore allows ranking potential binding sites but does not allow predictions concerning the relevance of these sites for a specific ligand.
Docking Protocol and Scoring.
A peptide-protein docking approach as implemented in Schrödinger Glide was performed (Friesner et al., 2006). The centroid of the interaction points in SiteMap served to center a Glide grid, suitable for peptide docking. This grid represents a cubic box (8000 Å3) that encapsulates the lower half of the three palm domains with one cube face positioned approximately parallel to the membrane plane. RPRFa docking was carried out with a standard-precision docking approach tailored for peptide ligands, SP-peptide (Tubert-Brohman et al., 2013), which ensures sufficient conformational sampling of the peptide. The Glide docking procedure penalizes ligand-receptor interactions producing steric clashes. Moreover, side chains of the receptor were allowed to relax their conformation around the docked ligand. Therefore, none of the binding poses had steric clashes. We scored the resulting binding poses with Glide Emodel (Friesner et al., 2006) since force field terms are particularly well suited to rank individual ligand conformers, rather than different ligands. A homology model of ASIC3 in its open conformation in complex with the best ranked docking pose of RPRFa is provided in Supplemental Fig. 3 in Protein Data Bank format (PDB). A protein-ligand interaction fingerprint (PLIF) was generated with MOE (Molecular Operating Environment, 2013.08; Chemical Computing Group ULC, Montreal, Canada).
Results
RPRFa Has to Unbind before ASIC3 Desensitizes.
To describe the interaction of RPRFa with ASIC3, Reimers et al. (2017) proposed a simplified kinetic model with three sequential states: closed (C), open (O) and desensitized (D) (Fig. 1A). According to this model, RPRFa binds to all of these three states with varying dissociation constants. The slower desensitization in the presence of RPRFa is then explained by a slower transition rate (τ2) from OR to DR, when RPRFa is bound. An alternative explanation would be that RPRFa has to unbind before the channel can desensitize. The slower desensitization in the presence of RPRFa would then be explained by the transition from OR to O, which would be significantly slower than the transition from O to D (Fig. 1A), leading to an indirect slowing of desensitization. The main difference between these two explanations is the existence of a desensitized conformation that has RPRFa bound (DR state), which should be readily populated according to the first explanation, but not according to the second explanation. To decide between these two explanations, we therefore tested the population of the DR state.
First, we examined whether RPRFa can bind to the desensitized conformation D of ASIC3. We desensitized ASIC3 with pH 6.0 and then applied 100 μM RPRFa for 60 seconds at pH 6.0. Next, we washed out the peptide by two short wash steps (pH 6.0 followed by pH 7.4, 2 seconds each). When we then activated ASIC3 by application of pH 6.0, we observed no slowing of desensitization, as measured by the ratio of the current 2.5 seconds post-peak (I2.5s) to the current at peak (I2.5s/Ipeak = 1.8% ± 0.2% for peak iii and I2.5s/Ipeak = 2.7% ± 0.4% for peak ii; iii vs. ii, P = 0.26; n = 12; Fig. 1B). In contrast, when we applied RPRFa for 60 seconds at pH 7.4, followed by two wash steps of identical duration as before (pH 7.4, 2 seconds each), desensitization of ASIC3 after activation by pH 6.0 was strongly slowed down (I2.5s/Ipeak = 71% ± 13% for peak iv; iv vs. iii; P < 0.001; Fig. 1B). We conclude that the affinity for RPRFa of the desensitized state D is much smaller than of the closed state C, owing to either a smaller binding rate or a larger unbinding rate.
Next, we tested whether RPRFa unbinds during current decline (Fig. 1C). We preapplied RPRFa (30 μM, 60 seconds) and then activated ASIC3 for 17 seconds with pH 6.0 in the absence of RPRFa. During this long activation current declined to <20% of its peak amplitude. It was followed by a short 2-second interval in pH 7.4 (no RPRFa), just enough to allow recovery of ASIC3. We then activated again with pH 6.0 and evaluated current decline. In case RPRFa stayed bound to ASIC3 during the long first current decline, the slowly desensitizing current component of the second activation should be larger than the current at the end of the previous activation. In contrast, in case RPRFa unbound during current decline, the slowly desensitizing current component of the second activation should be similar to or smaller than the current at the end of the previous activation. We found that the slowly desensitizing current component of the second activation was indeed similar to the current at the end of the previous activation (Fig. 1C, left panel). To rule out that RPRFa unbinding happened during the 2-second interval at pH 7.4, we preapplied RPRFa as before (60 seconds), followed by 2-second preconditioning at pH 7.4 in the absence of RPRFa. When we then activated with pH 6.0 and evaluated I2.5s/Ipeak, we found that channels were strongly modified (I2.5s/Ipeak = 47.8% ± 2.4%; n = 7). Although I2.5s/Ipeak was larger without the 2-second interval at pH 7.4 (I2.5s/Ipeak = 60.6% ± 2.6%, peak i; n = 9; iii vs. i, P = 0.004), this result still rules out that RPRFa unbinds quickly at pH 7.4, suggesting that current decline of ASIC3 after RPRFa-modification is accompanied by a substantial unbinding of RPRFa.
So far, we found no evidence that the desensitized conformation with RPRFa bound (DR) is strongly populated. If RPRFa indeed needs to unbind to allow desensitization of ASIC3, the irreversible binding of RPRFa should prevent desensitization of ASIC3. To test this prediction, we replaced the F residue of RPRFa by the photo-reactive unnatural amino acid 4-azido-l-phenylalanine (AzF) to yield RPR[azF]a. Upon irradiation by UV light, the azide moiety of AzF is converted to the highly reactive nitrene group, which allows the covalent attachment of RPR[azF]a to ASIC3.
RPR[azF]a slowed desensitization of ASIC3 likewise to RPRFa (Fig. 2A), even with a higher apparent affinity and a slower τdes than RPRFa (EC50: 1.2 ± 0.2 μM; τdes = 12.1 ± 3 seconds; n = 7; Fig. 2A) (Reimers et al., 2017), demonstrating that the peptide containing the AzF modification still binds to and modifies ASIC3. We then irradiated ASIC3-expressing oocytes with UV light for 1–10 minutes in 30–100 μM peptide (Fig. 2B). As control, we irradiated ASIC3-expressing oocytes with UV light for 6 minutes in the presence of 100 μM unmodified RPRFa. When cells treated with RPR[azF]a were subsequently activated by pH 6.0, a sustained current was observed, although no peptide was present or preapplied (Fig. 2C). This sustained current was not observed for control oocytes that had been incubated in RPRFa, suggesting that RPR[azF]a was covalently bound to ASIC3 and modulated the channel. The sustained current increased with increasing durations of photo labeling (I25s/Ipeak ratios: control, 0.3% ± 0.2%; UV treatment 1 minute, 3.7% ± 0.4%; UV 6 minutes, 20.5% ± 1.8%; UV 10 minutes, 59.8% ± 5.0%; P < 0.001 for all conditions with UV treatment; n = 5–10; Fig. 2D). But in contrast to modulation by RPRFa, the sustained current of covalently modulated channels persisted for the whole duration of pH 6.0 application (30 seconds) without any apparent desensitization. Moreover, the transient part of the current desensitized with a τdes that is typical for unmodified ASIC3; only at prolonged UV irradiation, τdes was slowed down (τdes: control, 0.44 ± 0.02 seconds; UV treatment 1 minute, 0.62 ± 0.06 seconds, P = 0.08; UV 6 minutes, 0.82 ± 0.07 seconds, P = 0.004; UV 10 minutes, 1.43 ± 0.14 seconds, P = 0.005). These results are expected when UV irradiation leads to two populations of ASIC3, one covalently modified by RPRFa and one unmodified. The modified population would be trapped in the open conformation (O-R) and would not desensitize at all, and the other fraction would desensitize with the usual τdes. The increase in the sustained current with increasing duration of UV irradiation speaks in favor of this interpretation.
To further test the presence of two populations (one covalently modified and one unmodified) in these experiments, we irradiated ASIC3-expressing oocytes with UV light for 4 minutes in 30 μM RPR[azF]a. As before, these oocytes displayed a nondesensitizing sustained current when subsequently activated by pH 6.0. Importantly, the transient current was still sensitive to RPR[azF]a (no UV irradiation), which slowed desensitization likewise to RPRFa (Fig. 2E). Thus, this result is in agreement with two populations of ASIC3 in UV-irradiated oocytes that had been incubated in RPR[azF]a.
Collectively, these results strongly argue that RPRFa has to unbind before ASIC3 can desensitize and that the kinetics of current decline of ASIC3 in the presence of RPRFa mainly reflects the slow unbinding of RPRFa (transition OR to O, Fig. 1A).
Molecular Docking Predicts the Nonproton Ligand-Sensing Domain as a Possible Binding Site for RPRFa.
To discover possible binding sites of RPRFa on ASIC3, we modeled the open and the desensitized conformations of rat ASIC3 on the open and the desensitized conformation of chicken ASIC1, respectively (PDB ID: 4NTW and 2QTS, see Materials and Methods for further details). We then performed an unbiased binding sites detection protocol (see Materials and Methods for detail) across these two in silico models and found several potential binding sites; for each model, the five highest ranking sites according to SiteScore are reported in Supplemental Table 1. Since experimental evidence, as reported above, suggests that RPRFa unbinds before ASIC3 desensitizes, the binding site of RPRFa should be located in a region accessible only in the open state but not in the desensitized state. The largest binding site in the open-state model that is not accessible in the desensitized state model is a cavity in the lower palm domain, the NPLSD. Indeed, in cASIC1 the lower palm domain slightly shrinks on passing from the open to the desensitized state (Baconguis and Gouaux, 2012) and this is also reflected in our models. The conformational changes in the lower palm domain, when transitioning from open to the desensitized state, suggest that RPRFa binding to the NPLSD might be sterically unfavorable in the desensitized state.
RPRFa was therefore docked into the NPLSD of the open-state homology model of ASIC3 (Fig. 3). Protein-peptide docking is more challenging than small-molecule docking, since even short peptides are often highly flexible, and several binding poses with comparable probability are possible. We obtained 34 distinct binding poses (Supplemental Table 2). These are summarized by a protein ligand interaction fingerprint (PLIF) showing the interaction between RPRFa and ASIC3 (Fig. 4B; Supplemental Fig. 4). High-ranking poses form cation-π interactions with the C-terminal F of RPRFa and R376 in β-strand 10, whereas the amide capping group hydrogen bonds with Q269 of β9. Salt bridges can be observed between the C-terminal R of RPRFa and E79/E423 as well as for the N-terminal R of RPRFa and E418 (Fig. 4A). Notably, the interactions between the peptide and receptor can involve more than one monomer at the same time (Supplemental Fig. 4), possibly altering monomer interaction. Supplemental Fig. 3 provides the homology model of ASIC3 in its open conformation in complex with RPRFa in its best-ranked docking pose and Supplemental Fig. 5 provides a movie of this complex.
Cartoon representation of ASIC3 in complex with RPRFa. (A) Top view. (B) Side view with one monomer removed; highest ranked RPRFa docking pose in green licorice. The cartoon was made with PyMol (The PyMOL Molecular Graphics System, Version 2.0; Schrödinger, LLC). (C) Higher magnification view of the highest ranked docking pose with key interactions highlighted.
Binding of RPRFa to the NPLSD of ASIC3. (A) Two-dimensional protein-ligand interaction diagram of the best binding pose. (B) Protein-ligand interaction fingerprint (PLIF) as population histogram. For each amino acid of the receptor as many columns are displayed as there are distinguishable noncovalent interactions in at least one pose. The individual scores are provided in Supplemental Table 2 and the barcode representation in Supplemental Fig. 4.
The Nonproton Ligand–Sensing Domain Is Involved in Modulation by RPRFa.
The NPLSD is mainly encapsulated by four β-strands of each of the three palm domains (Jasti et al., 2007). To test experimentally the prediction by in silico docking, we mutated eight residues in the NPLSD that surrounded RPRFa docked to the homology model: L77 and E79 in β-strand 1, E423 in β12, Q269 and Q271 in β9, and R376, A378, and E380 in β10 (total of 23 mutations; Fig. 5; Supplemental Table 3). First, we activated ASIC3 wild-type and each mutant with pH 6.3 and pH 4.0 and then preapplied 3 μM RPRFa followed by another activation with pH 6.3 (Fig. 5B). To quantify modulation by RPRFa, we calculated the difference of the current remaining after 2.5 seconds with and without RPRFa (I2.5s+RFa – I2.5s-RFa) and normalized it to the peak current without RPRFa ([I2.5s+RFa – I2.5s-RFa]/IPeak-RFa). Out of a total of 23 mutants, for 18 mutants [I2.5s+RFa – I2.5s-RFa]/IPeak-RFa was significantly reduced, compared with the wild-type (Fig. 5D; Supplemental Table 3). Notably, all substitutions of R376 and Q269, two amino acids that appeared to be crucial for interaction with the Famide group of RPRFa in in silico docking (Fig. 4), significantly reduced modulation by RPRFa. Likewise, substitutions of E423, which makes extensive contacts to RPRFa in most docking poses (Fig. 4B), also significantly reduced modulation by RPRFa, with the exception of the relatively conservative exchange E423Q (Fig. 5D). As an alternative way to quantify modulation by RPRFa, we also calculated the ratio I2.5s+RFa/I2.5s-RFa. Also in this type of analysis, 17 out of the 23 mutations in the NPLSD, including all substitutions of R376 and Q269, significantly reduced modulation by RPRFa (Supplemental Fig. 6; Supplemental Table 3). Together, these results show that the NPLSD is crucial for the modulation of ASIC3 by RPRFa. To test whether the mutations specifically changed the interaction of amino acids of the NPLSD with RPRFa, we additionally generated mutants Q270E and K379E, which are also located on β-strands 9 and 10 but the side chains of which point away from the cavity to which RPRFa was docked. Both mutants indeed did not significantly affect the modulation by RPRFa as measured by the ratio [I2.5s+RFa – I2.5s-RFa]/IPeak-RFa (Fig. 5D; Supplemental Table 3). For K379E, the ratio I2.5s+RFa/I2.5s-RFa was reduced, albeit less strongly than for the other mutants (Supplemental Fig. 6; Supplemental Table 3). The results are in agreement with the interpretation that RPRFa was bound to the NPLSD.
Mutation of residues in the NPLSD reduce the modulation by RPRFa. (A) Scheme illustrating the four β strands of one subunit that surround the cavity in the lower palm domain. Approximate positions of amino acids that point toward the cavity and have been mutated are indicated. (B) Representative current trace of an ASIC3-expressing oocyte when activated by pH 6.3 (black bar) or pH 4.0 (blue bar) with or without preapplication of RPRFa (3 μM, red bar). (C) Representative pH-activated currents of NPLSD mutants in the absence of (black trace) and after preincubation with (red trace) 3μM RPRFa were overlaid. (D) Mean [I2.5s+RFa – I2.5s-RFa]/IPeak-RFa ratios for all 25 mutants of the NPLSD. Bars are shown in the color of the corresponding β-strands in (A); white bars represent residues whose side chains point away from the cavity. Error bars indicate S.E.M. nmutant = 8–19; nwt = 4–28 (each condition), nwt (total) = 102. *P < 0.05; **P < 0.01; ***P < 0.001 (Student’s t test followed by Bonferroni correction). WT, wild type.
Evidence that Mutants within the Nonproton Ligand–Sensing Domain Reduce Binding of RPRFa and Fasten Unbinding.
In principle, there are two ways by which the mutations might have reduced the modulation by RPRFa: first, by reducing binding, or second, by reducing the effect of binding on desensitization. To further differentiate between these two possibilities, for 12 mutants we also tested a higher RPRFa concentration (100 μM) (Fig. 6); for wild-type ASIC3, 100 μM RPRFa is a saturating concentration (Reimers et al., 2017). When looking at a shorter time scale (0.5 seconds after peak), RPRFa was indeed able to modulate desensitization kinetics of many mutants (Fig. 6B). To quantify RPRFa effects, we determined [I0.5s+RFa – I0.5s-RFa]/IPeak-RFa. For wild-type ASIC3 this ratio was 0.90 ± 0.07, whereas it was significantly reduced to a range of −0.10 to 0.51 for 11 mutants (L77F, L77C, E79S, E423R, Q269R, Q269C, Q271R, R376E, A378R, A378C, and E380Q; Fig. 6C; Supplemental Table 4), confirming the importance of these amino acids for RPRFa modulation of ASIC3.
NPLSD mutants are affected by high concentrations of RPRFa. (A) Representative recording of an ASIC3-expressing oocyte, illustrating the application protocol. ASIC3 was first activated by pH 6.3 (black bar) without peptide and then after preapplication of 100 μM RPRFa (red bar). (B) Representative pH-activated currents of NPLSD mutants in the absence of (black trace) and after preincubation with (red trace) 3 μM RPRFa were overlaid. Gray vertical lines indicate the time points 0, 0.5, and 2.5 seconds post-peak. (C) Bar graph displaying mean ratios of current amplitudes 0.5 seconds post-peak in the presence and in the absence of RPRFa ([I0.5s+RFa – I0.5s-RFa]/IPeak-RFa). Error bars indicate S.E.M.; nwt = 17, nmutant = 6–12. (D) Current decay after preapplication of RPRFa was fit with double exponential functions yielding a2 (top) and the ratio τ2/τ1 (bottom). Mutant Q271R was not analyzed, because it could not be fit properly. nwt = 15, nmutant = 6–10. *P < 0.05; **P < 0.01; ***P < 0.001. WT, wild type.
We then fitted current decay after RPRFa application by double exponential fits. We fixed the fast time constant τ1 to the value obtained by a single exponential fit in absence of RPRFa (see Materials and Methods) and determined the fraction a2 and the slow time constant τ2. This procedure specifically detects and characterizes components of the current decay that are introduced by RPRFa. Since for some mutants τ1 was different from wild-type, we analyzed the ratio τ2/τ1 to evaluate desensitization in the presence of RPRFa. For wild-type ASIC3, 100 μM RPRFa modified 91% ± 3% of the current (a2) and decelerated desensitization about 22.0 ± 1.4-fold (τ2/τ1; n = 15). In contrast, for L77C, E79S, R376E, A378R, and A378C, RPRFa modified a significantly smaller fraction of the current (a2 ranging from 17% to 67%; Fig. 6D, top, and Supplemental Table 4), which suggests that these mutants had less RPRFa bound. The ratio τ2/τ1 was significantly decreased 2.2–11.1-fold for all mutants with the exception of E423R (Fig. 6D, bottom and Supplemental Table 4; Q271R could not be fitted). As outlined above, the speed of current decline in the presence of RPRFa probably reflects unbinding of RPRFa (transition OR to O; Fig. 1A). A decreased ratio τ2/τ1 therefore indicates faster unbinding of RPRFa (smaller τ2), and thus provides another indication that the mutations within the NPLSD reduced affinity of ASIC3 for RPRFa. Taken together, our results reveal the NPLSD as a good candidate for the RPRFa binding site.
The Acidic Pocket of ASIC3 Is Not Involved in Modulation by RPRFa.
To further prove that mutations in the NPLSD specifically changed modulation by RPRFa, we made 12 individual amino acid substitutions within the acidic pocket of ASIC3 (Fig. 7A), the binding site of a peptide toxin on ASIC1a (Baconguis and Gouaux, 2012; Dawson et al., 2012) and a highly ranked potential binding site according to SiteScore (see above), and tested their effect on RPRFa modulation (Fig. 7B). In strong contrast to mutations within the NPLSD, none of the mutations within the acidic pocket significantly reduced the ratio [I2.5s+RFa – I2.5s-RFa]/IPeak-RFa compared with wild-type (Fig. 7C). We also calculated the ratio I2.5s+RFa/I2.5s-RFa. For the wild-type, preapplication of 3 μM RPRFa increased I2.5s 31.7 ± 3.3-fold (n = 102) relative to I2.5s when no RPRFa was preapplied. For three of the mutants, I2.5s+RFa/I2.5s-RFa was significantly decreased (R183E, E212C, and R239E; Supplemental Fig. 7; Supplemental Table 3). In contrast to an effect on RPRFa modulation, many mutations in the acidic pocket (R183E, M184Y, E212R, E212C, E231R, and R362E) significantly decreased apparent proton affinity as assessed by the ratio IpH6.3/IpH4.0 (Supplemental Table 3). Taken together these results suggest that the acidic pocket is not involved in the modulation of ASIC3 by RPRFa but might be involved in pH sensing.
Mutation of residues in the acidic pocket do not impair modulation by RPRFa. (A) Scheme of the topology of the acidic pocket of ASIC3. (B) pH 6.3-activated currents with (red trace) or without (black trace) 3 μM RPRFa were taken from recordings as the one shown in Fig. 5B and overlaid. (C) Bar graph showing mean ratios [I2.5s+RFa – I2.5s-RFa]/IPeak-RFa for the 14 mutants of the acidic pocket. The dotted line represents the ratio of the wild type as reference. Error bars indicate S.E.M. nmutant = 8–18; nwt = 6–20 (each condition), nwt (total) = 102. WT, wild type.
Discussion
Our results show: 1) that RPRFa has to unbind before ASIC3 can desensitize, 2) that RPRFa docks to the NPLSD of an ASIC3 open-state model, and 3) that mutations within the NPLSD, but not within the acidic pocket, reduce modulation of ASIC3 by RPRFa. The finding that RPRFa has to unbind before ASIC3 can desensitize suggests that the slow time constant τ2 of the current decay phase, which is introduced by RPRFa (Reimers et al., 2017) and which is commonly described as slower desensitization of ASIC3, in reality reflects unbinding of RPRFa (Fig. 1A). Thus, it appears that RPRFa uses a simple mechanism to modify gating of ASIC3: by stabilizing the open conformation.
Many mutations within the NPLSD, a cavity of the lower palm domain, reduced the modulation of ASIC3 by RPRFa. Mutation of two amino acids that do not point toward the interior of the cavity, however, did not reduce RPRFa modulation. In addition, the reduced a2 in the bi-exponential fits of mutant ASIC3 and the faster slow-time constant τ2, in our interpretation indicate reduced binding of RPRFa and faster unbinding of RPRFa, respectively. These findings suggest that RPRFa indeed binds to the NPLSD. In agreement with this idea, the lower palm has been implicated in desensitization by several previous studies (Cushman et al., 2007; Springauf and Gründer, 2011; Roy et al., 2013). In particular, comparison of the open and desensitized conformation revealed that the ECD contracts in the lower palm domain that encapsulates the negatively charged central vestibule (Baconguis et al., 2014), which forms the NPLSD of ASIC3 (Yu et al., 2010). A contraction of the lower palm domain was later confirmed by functional studies; preventing the contraction impairs desensitization (Roy et al., 2013). Moreover, interfering with the conformational changes in the lower palm domain induces sustained currents by destabilizing the desensitized conformation (Springauf and Gründer, 2011). Of particular interest, it has been shown that FRRFa limits covalent modification of a residue in the lower palm domain of ASIC1a (Frey et al., 2013). All these findings are consistent with binding of RPRFa to the NPLSD of ASIC3.
In contrast to mutations within the NPLSD, mutations within the acidic pocket did not reduce modulation of ASIC3 by RPRFa. The acidic pocket has been involved in activation of ASICs by protons (Jasti et al., 2007; Gründer and Chen, 2010; Krauson et al., 2013; Ramaswamy et al., 2013) and binds psalmotoxin 1, a gating modifier toxin that acts like an agonist of ASIC1 (Chen et al., 2005, 2006a). But the acidic pocket has not been implicated in desensitization of ASICs. Our finding that the acidic pocket apparently is not involved in RPRFa modulation is, thus, in agreement with the previous findings on the role of the acidic pocket for ASIC gating.
Collectively, our results are consistent with the following mechanistic model: RPRFa slowly (within several seconds) binds to a cavity in the lower palm domain, the NPLSD, in the closed or the open conformation and sterically stabilizes the open conformation by preventing the contraction of the central vestibule. Unbinding of RPRFa with a time constant on the order of 6 seconds (Reimers et al., 2017) would then allow delayed contraction of the central vestibule and desensitization of ASIC3. Since the time constant of desensitization of ASIC3 (0.3 seconds) (Sutherland et al., 2001) is approximately 10-fold faster than unbinding of RPRFa, RPRFa unbinding would dominate the kinetics of current decay. Although this is a plausible explanation of our results, we can, at present, not exclude other models in which RPRFa binds to a site different from the NPLSD on ASIC3 and, for example, allosterically enhances unbinding of RPRFa. The availability of RPR[AzF]a as a covalent ligand provides the opportunity to unequivocally determine the binding site of RPRFa on ASIC3 in the future, by identifying the covalently modified channel fragments, for example by tandem MS analysis.
We have previously shown that the prototypical RFamide FMRFa competes with RPRFa for modulation of ASIC3 (Reimers et al., 2017), suggesting that both RFamides share a binding site. Since all RFamides have similar effects on ASICs (Vick and Askwith, 2015), mainly slowing current decay (Askwith et al., 2000; Deval et al., 2003; Chen et al., 2005; Reimers et al., 2017), it is conceivable that all RFamide neuropeptides, which modulate ASICs, share a binding site on ASICs. Thus, the cavity in the lower palm domain is a probable common binding site for RFamide neuropeptides on ASICs. Although RFamides also modulate ASICs other than ASIC3 (Askwith et al., 2000), ASIC3 and ASIC3-containing heteromers seem to be more strongly modulated by RFamides than other ASICs (Deval et al., 2003; Chen et al., 2006b; Vick and Askwith, 2015); RPRFa, in particular, modulates only ASIC3 and different ASIC3-containing heteromers but not ASIC1a (Reimers et al., 2017). It is probable that structural differences in the central vestibule are responsible for the differential sensitivity of different ASICs to RFamides. It might be that RFamides bind with lower affinity to the NPLSD of ASIC1a, for example. Since ASIC1a desensitizes more slowly than ASIC3, an enhanced unbinding of RFamides, owing to lower affinity, could result in a situation in which the kinetics of current decay would still be dominated by desensitization rather than by RFamide unbinding, effectively reducing the functional outcome of RFamide binding to the channel. Future studies are necessary to reveal the molecular basis for the subtype-specific action of RFamides.
If RFamide neuropeptides indeed bind to the NPLSD of ASIC3 they would join several other ligands for which there is evidence of binding to the NPLSD. 2-Guanidine-4-methylquinazoline (GMQ) is a synthetic ligand that already induces sustained ASIC3 currents at pH 7.4 (Yu et al., 2010) and also binds to the NPLSD (Yu et al., 2010). Also, different natural ligands bind to the NPLSD, such as the arginine metabolite agmatine (Li et al., 2010; Yu et al., 2010), serotonin (Wang et al., 2013), and the pruritogenic peptide SLIGRLamide (Peng et al., 2015). Thus, the NPLSD appears as an important binding cavity in the ECD of ASIC3. All these ligands, including RPRFa, are positively charged, containing guanidine or amine groups, which suggests that they are attracted by the negative charges of the central vestibule (Jasti et al., 2007). Moreover, the presence of two guanidine and one amino group in RPRFa might explain its particularly high affinity for ASIC3 (Reimers et al., 2017), for example compared with FMRFa, which has only one guanidine and one amino group.
Apart from the positive charges of the ligands binding to the NPLSD, their precise structures are quite diverse, however. This situation is reminiscent of the delayed rectifier K+ channel hERG, where many structurally diverse compounds bind to a common inner cavity to inhibit the channel (Mitcheson et al., 2005). This cavity is particularly large in hERG compared with other related K+ channels, explaining the drug promiscuity of hERG. Perhaps the central vestibule of ASIC3 is also larger than that of other ASICs explaining its sensitivity to different compounds, some of which are quite bulky. It might be that the central vestibule of ASIC3 evolved to a cavity where different physiologic ligands bind to modulate the activity of ASIC3.
In summary, our results reveal how RPRFa slows current decay of ASIC3 and that the NPLSD is important for the modulation of ASIC3 by RPRFa. Although we cannot exclude other models, we propose a mechanistic model in which RPRFa binds to the NPLSD preventing the contraction of the central vestibule, which appears to be sufficient to slow desensitization.
Acknowledgments
We thank D. Wiemuth for many helpful discussions.
Authorship Contributions
Participated in research design: Schmidt, Gründer.
Conducted experiments: Reiners, Margreiter, Oslender-Bujotzek, Schmidt.
Performed data analysis: Reiners, Margreiter, Schmidt, Rossetti.
Wrote or contributed to writing of the manuscript: Rossetti, Schmidt, Gründer.
Footnotes
- Received March 12, 2018.
- Accepted July 6, 2018.
Financial support was provided by the Medical Faculty of RWTH Aachen University.
↵
This article has supplemental material available at molpharm.aspetjournals.org.
Abbreviations
- ASIC
- acid-sensing ion channel
- azF
- 4-azido-l-phenylalanine
- ECD
- extracellular domain
- GMQ
- guanidine-4-methylquinazoline
- NPLSD
- nonproton ligand–sensing domain
- PDB
- protein data bank
- RPRFa
- RPRFamide
- Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics