Desensitization is a common feature of ligand-gated ion channels, although the molecular cause varies widely between channel types. Mutations that greatly reduce or nearly abolish desensitization have been described for many ligand-gated ion channels, including glutamate, GABA, glycine, and nicotinic receptors, but not for acid-sensing ion channels (ASICs) until recently. Mutating Gln276 to a glycine (Q276G) in human ASIC1a was reported to mostly abolish desensitization at both the macroscopic and the single channel levels, potentially providing a valuable tool for subsequent studies. However, we find that in both human and chicken ASIC1, the effect of Q276G is modest. In chicken ASIC1, the equivalent Q277G slightly reduces desensitization when using pH 6.5 as a stimulus but desensitizes, essentially like wild-type, when using more acidic pH values. In addition, steady-state desensitization is intact, albeit right-shifted, and recovery from desensitization is accelerated. Molecular dynamics simulations indicate that the Gln277 side chain participates in a hydrogen bond network that might stabilize the desensitized conformation. Consistent with this, destabilizing this network with the Q277N or Q277L mutations largely mimics the Q277G phenotype. In human ASIC1a, the Q276G mutation also reduces desensitization, but not to the extent reported previously. Interestingly, the kinetic consequences of Q276G depend on the human variant used. In the common G212 variant, Q276G slows desensitization, while in the rare D212 variant desensitization accelerates. Our data reveal that while the Q/G mutation does not abolish or substantially impair desensitization as previously reported, it does point to unexpected differences between chicken and human ASICs and the need for careful scrutiny before using this mutation in future studies.
Introduction
Desensitization is a near-ubiquitous feature of ligand-gated ion channels (LGICs), which was first described >60 yr ago (Katz and Thesleff, 1957). In general, desensitization is thought to act as a protective mechanism, terminating aberrant signaling, although other roles are possible (Gielen et al., 2020; Jones and Westbrook, 1995; Papke et al., 2011). As such, the molecular basis of desensitization has been a subject of inquiry for every type of LGIC. Mutations that essentially abolish or greatly reduce desensitization have been reported for glutamate, GABA, glycine, and nicotinic receptors (Bertrand et al., 1992; Gielen and Corringer, 2018; Gielen et al., 2015; Nayeem et al., 2009; Revah et al., 1991; Stern-Bach et al., 1998). While there have been controversies surrounding the microscopic mechanisms of particular cases (Daniels et al., 2013), these mutations have been enormously helpful in driving structure–function investigations of desensitization as well as providing insight into the physiological role (Christie et al., 2010). There are mutations that somewhat attenuate desensitization of acid-sensing ion channels (ASICs) to varying degrees (Li et al., 2010; Passero et al., 2009; Rook et al., 2020b; Roy et al., 2013; Vullo et al., 2017; Yoder et al., 2018), but no mutations purported to nearly block or abolish desensitization have been described until recently (Chen et al., 2021; Wu et al., 2019).
ASICs are sodium-selective pH-activated trimeric ion channels. They are expressed widely in the central and peripheral nervous systems, as well as other tissues (Kellenberger and Schild, 2015). Given the ubiquity of the ligand, it is unsurprising that ASICs are implicated in a host of physiological processes and disease states including ischemic cell death, fear and anxiety, learning and memory, pain, muscle fatigue, migraine, bone morphogenesis, inflammation, and cancer (Sluka and Gregory, 2015; Wemmie et al., 2013; Zhu et al., 2017). In mammals, the ASIC family includes four subunits capable of forming proton-sensitive homomers: ASIC1a, ASIC1b, ASIC2a, and ASIC3. The individual subunits all have the same topology with intracellular amino- and carboxy-terminal tails of ∼20–80 amino acid residues, separated by a large extracellular domain, two transmembrane helices, and a small amino-terminal reentrant loop (Rook et al., 2020a; Yoder and Gouaux, 2020). The extracellular domain is divided into distinct thumb, finger, knuckle, palm, and β-ball domains (Fig. 1 A). ASIC activation by acidic conditions is believed to occur through protonation of distinct residues in the interface between the thumb and finger as well as a cluster of titratable side chains in the palm domain (Jasti et al., 2007; Krauson et al., 2013; Rook et al., 2020a; Vullo et al., 2017). Protonation also triggers desensitization, either with mild acidic stimuli (pH in the 7.4–6.9 range), which leads to steady-state desensitization (SSD) in the absence of channel activation, or with strong stimuli (i.e., pH 6.8–4), which also opens the channel. Desensitization depends on the isomerization or swivel of a critical linker in the palm domain, which connects the 11th and 12th β strands (Rook et al., 2020b; Yoder et al., 2018). This linker is composed of Leu414 and Asn415 (Fig. 1 B), and in the resting and open states the Leu residue points outward, away from the central axis of the channel. However, in the desensitized state, these amino acid residues essentially switch positions, with Leu414 swiveling downward and in toward the central axis. It has been suggested that in human ASIC1a, Gln276 (Gln277 in chicken ASIC1) acts as a valve to prevent linker rotation, stabilizing the desensitized state, and, furthermore, that eliminating the 276 side chain using the Q276G mutation in effect produces a “leaky” valve that enables channels to readily escape desensitization and remain open (Wu et al., 2019).
To further study this intriguing Q276G mutant and relate it to the majority of structural data, we tested the Q276G equivalent in cASIC1 (Q277G) using piezo-driven fast perfusion in excised patches. We found that when using pH 6.5 to open the channels, cASIC1 Q277G does have slightly reduced desensitization; however, when using more acidic stimuli, Q276G behaves essentially like WT with desensitization principally intact. Moreover, we found that Q277G accelerates recovery from desensitization by orders of magnitude and reduces the apparent stability of the desensitized state. Based on molecular dynamics simulations, we hypothesize that Gln277 coordinates a series of hydrogen bonds within the palm domain, thereby stabilizing the desensitized conformation. Consistent with this electrostatic mechanism, a slight mutation to Q277N or Q277L also accelerates recovery from desensitization. Finally, we find that hASIC1a Q276G exhibits robust pH-dependent desensitization, in contrast to prior work.
Materials and methods
Human embryonic kidney (HEK) 293T ASIC knockout (KO) cell creation
A guide RNA sequence (5′-GGCTAAAGCGGAACTCGTTG-PAM-3′) targeting the coding region of ASIC1 was cloned into Bbsl-linearized pSpCas9(BB)-2A-GFP vector (a kind gift from Feng Zhang, Massachusetts Institute of Technology, Boston, MA; Addgene plasmid #48138) as previously described (Ran et al., 2013). Transfected HEK cells (HEK293T; American Type Culture Collection no. CRL-3216) were verified for GFP expression and clonally expanded following serial dilution. Clonal lines were screened for on-target genome editing by Sanger sequencing of PCR products (Fwd: 5′-TTGGAGGAACCCTGGATGTGTC-3′, and Rev: 5′-TAACTCCTCTGCTGTGAGTGGC-3′). KO was confirmed by Western blotting. Briefly, 105 cell equivalents of radioimmunoprecipitation assay buffer lysate from the parental cell line, KO clone, and clone transiently transfected with human ASIC1a cDNA were resolved on an acrylamide gel and transferred to nitrocellulose membranes. Blots were blocked with bovine serum albumin and probed with an ASIC1-specific antibody (NeuroMab clone N271/44) overnight at a 1:1,000 dilution. Blots were washed with Tris-buffered saline supplemented with 0.1% Tween-20 and probed with goat anti-mouse IgG horseradish peroxidase (HRP)–conjugated secondary antibody. Blots were imaged with an Azure 300 Imaging System. After inactivation of HRP with sodium azide, the blot was probed again with Direct-Blot HRP anti-GAPDH (Biolegend) antibody as a loading control.
Cell culture, mutagenesis, and transfection
HEK293T ASIC KO cells were maintained in Dulbecco’s modified Eagle’s medium with 4.5 g/liter glucose, L-glutamine, and sodium pyruvate (Corning/Mediatech, Inc.) or Eagle’s minimum essential medium with Glutamax & Earle’s Salts (Gibco), supplemented with 10% FBS (Atlas Biologicals) and penicillin/streptomycin (Invitrogen). Cells were passaged every 2–3 d when ∼90% confluence was achieved. HEK293T KO cells were plated on tissue culture–treated 35-mm dishes, transfected 24–48 h later, and recorded from 12–48 h after transfection. Cells were transiently transfected with the indicated ASIC construct and enhanced GFP (eGFP) using an ASIC:eGFP ratio of between 2.5 and 10:1 µg of cDNA per 10 ml of medium, depending on the construct. For hASIC1a whole cell experiments (Fig. 6), a ratio of 0.25:0.25:1 μg of hASCI1a, eGFP, and pUC empty vector was used. We used commercial gene synthesis (IDT) to create human codon optimized cDNA of hASIC1a with the G212 variant that was then inserted into pcDNA 3 at the EcoRV site. Transfections were performed using polyethylenimine 25k (Polysciences, Inc.) following the manufacturer’s instructions, with media change at 1–8 h after transfection. In general, WT recordings were obtained sooner after transfection using short transfection durations (i.e., 24 h and 1–2 h incubation), while Q276G recordings were obtained later and with longer transfection durations (i.e., 48 h and 8 h incubation). Mutations were introduced using site-directed mutagenesis PCR and confirmed by sequencing (Thermo Fisher Scientific/Eurofins Genomics).
Electrophysiology
Culture dishes were visualized with phase contrast on a Nikon Ti2 microscope using a 20× objective. GFP was excited using a 455-nm or 470-nm LED (Thorlabs) and dichroic filter cube for emission detection. Outside-out patches were excised using heat-polished, thick-walled borosilicate glass pipettes of 3–15 MΩ resistance. The pipette internal solution contained (in mM) 135 CsF, 33 CsOH, 11 EGTA, 10 HEPES, 2 MgCl2, and 1 CaCl2, pH 7.4. External solutions with pH values >7 were composed of (in mM) 150 NaCl, 10 HEPES, 1 CaCl2, and 1 MgCl2, with pH values adjusted to their respective values using NaOH. For solutions with a pH value <7, HEPES was replaced with Mes. All recordings were performed at room temperature with a holding potential of −60 mV using an Axopatch 200B amplifier (Molecular Devices). Data were acquired using AxoGraph software (Axograph) at 20 kHz, filtered at 10 kHz, and digitized using a USB-6343 DAQ (National Instruments). Series resistance was routinely compensated by 90–95% where the peak amplitude exceeded 100 pA. Rapid perfusion was performed using home-built, double- or triple-barrel application pipettes (Vitrocom), manufactured according to a prior method (MacLean, 2015). Application pipettes were translated using piezo actuators driven by voltage power supplies. The command voltages were generally low-pass–filtered (50–100 Hz, eight-pole Bessel). Whole-cell recording (Fig. S1 and Fig. 6) used identical conditions except patch pipette and application pipette diameters tended to be larger.
Molecular dynamics simulations
The systems for molecular dynamics simulations were constructed using the cASIC1 structure proposed to illustrate the desensitized state (Protein Data Bank [PDB] accession no. 4NYK; resolution: 3.00 Å; Gonzales et al., 2009). Residues 42–455 were resolved in the crystal structure; of these, 23 residues had missing atoms, which were added using Modeller v9.21 (Sali and Blundell, 1993). The model was oriented for placement in a lipid bilayer by aligning the complete structure with the corresponding structure from the Orientations of Proteins in Membranes database (Lomize et al., 2012). Protonation states of specific residues were set using pdb2gmx during the system setup in GROMACS (Abraham et al., 2015). In all cases, residues Glu98, His111, Glu239, His328, Glu354, Glu374, Asp408, and Asp433 were protonated, leaving the acidic residues neutral and histidine residues with a positive charge. This was considered the “background” protonation setup, and the purpose was to maintain an overall stable protein structure. However, the importance of the presence of these individual protons was not tested in this study as they are relatively far away from our region of interest. On the given background, Glu80, Glu412, and Glu417 were protonated or deprotonated in accordance with Table 1 to test the importance of protonation of these specific residues.
All systems were prepared from scratch from the original protein structure and with the given protonation scheme applied (“background” plus a specific combination for E80, E412, and E417) before insertion into the membrane and equilibration of each system. The Charmm36m force field was applied (Huang et al., 2017). The initial POPC lipid bilayer (120 Å × 120 Å) was generated using the membrane builder of the CHARMM-GUI with PDB accession no. 4NYK inserted using the replacement method (Jo et al., 2008). The protein structures with different protonation states were then inserted into this original lipid bilayer using the InflateGro method (Kandt et al., 2007). The crystallographic water molecules and chloride ions were retained. Water molecules (TIP3P model; Jorgensen et al., 1983) were further generated to fill the box (120 Å × 120 Å × 161 Å) with solvent, and sodium and chloride ions were added to neutralize the system at a concentration of 0.15 M NaCl.
The simulations were performed using GROMACS 2019.4 (Abraham et al., 2015). All systems were minimized until convergence or to a maximum of 5,000 steps. The systems were then equilibrated in six steps totaling 375 ps, using the standard method from the CHARMM-GUI. The first three equilibration runs used a time step of 1 fs, while the last three and the production run used a time step of 2 fs. The first three equilibration runs were each 25 ps long, and the final three were each 100 ps long. The position restraints were gradually lifted during the equilibration steps as suggested in the default CHARMM-GUI protocol. Periodic boundary conditions were applied. The Verlet cutoff scheme was used throughout with a force-switch modifier starting at 10 Å and a cutoff of 12 Å. A cutoff of 12 Å was used for short-range electrostatics, and the particle mesh Ewald method was used for long-range electrostatics (Darden et al., 1993; Essmann et al., 1995). A Berendsen thermostat was used for all steps of the equilibration, and a Nosé–Hoover thermostat (Hoover, 1985; Nosé, 1984) was used in the production run to maintain the temperature at 310.15 K for all steps. Using semi-isotropic pressure coupling, the pressure was maintained at 1 bar in the last four steps of equilibration and in the production run using the Berendsen barostat (Berendsen et al., 1984) and the Parrinello–Rahman barostat (Parrinello and Rahman, 1981), respectively. The LINCS algorithm was used to constrain covalent bonds to hydrogen atoms (Hess, 2008). The production runs were 100 ns long with a total of three repeats for each system. Each repeat had different initial velocities. The system for the Q277N mutant was prepared as above, with the exception that the Gln277 side chain was manually mutated to Asn before system setup. The same background protonation scheme was used, and additionally, Glu412 and Glu417 were protonated. The system was simulated as described with three repeats of 100 ns each.
Potential hydrogen bonds between the residues Glu80, Gln277, Glu412, Leu414, and Glu417 were sampled every 10 ps. The donor atoms include Q277NH1, Q277NH2, E80HE2 (if protonated), E412HE2 (if protonated), and E417HE2 (if protonated). The acceptor atoms include Q277OE1, E80OE1, E80OE2, E412OE1, E412OE2, L414O, E417OE1, and E417OE2. The 4.3.1 Hydrogen Bond Analysis module (Smith et al., 2019) of MDAnalysis (Michaud-Agrawal et al., 2011; Gowers, R.J., et al. 2016. SciPy2016. https://doi.org/10.25080/majora-629e541a-00e) was used for the analysis, using an updated and adapted version of M. Chavent’s Jupyter Notebook available on GitHub (https://github.com/MChavent/Hbond-analysis; Del Toro et al., 2020). Default cutoffs were used for the donor-acceptor distance (3.0 Å) and the donor-hydrogen-acceptor angle (150°). The presence of each unique hydrogen bond was calculated over the trajectory and expressed as a percentage of the total trajectory; the presence of equivalent hydrogen bonds (e.g., from 417OE1 and 417OE2 in the deprotonated state) was added to give one overall percentage for the given interaction. Plots were prepared using the Matplotlib package in Python. Figures were prepared using visual molecular dynamics (Humphrey et al., 1996).
Statistics and data analysis
Unless otherwise noted, to avoid assumptions about the underlying distribution of the data, statistical testing was done using nonparametric permutation tests with at least 100,000 iterations implemented using the mlxtend package in Python (http://rasbt.github.io/mlxtend/). Statistical comparisons of recovery from desensitization were based and reported on differences in recovery time constant.
Online supplemental material
Fig. S1 shows the validation of hASIC1 KO cells using sequencing, Western blot, and patch clamp. Fig. S2 depicts the hydrogen bond stability of all possible protonation schemes for residues Glu80, Glu412, and Glu417. Fig. S3 shows a similar hydrogen bond analysis of the Glu80/Glu412H/Glu417H scheme for both WT and Q277N simulations. Video 1 and Video 2 show trajectories of Gln277 or Asn277, respectively, upon deprotonation using the E80-/E412H/E417H scheme.
Results
Our goal was to investigate the functional properties of the Q276G mutation in a cASIC1 background, to permit easy comparison with structural data and molecular dynamic simulations. HEK293T cells are an ideal system for this as they are easily cultured and transfected, and amenable to patch clamp. However, HEK293T cells express endogenous human ASIC1, which may complicate interpretation. Therefore, we removed the endogenous human ASIC1 using CRISPR. To do this, exon 2 of the human ASIC1 gene was targeted with a guide RNA cloned into a Cas9-GFP expressing vector (see Materials and methods). Single GFP-positive HEK293T cells were clonally expanded and screened using PCR followed by sequencing (Fig. S1 A). One such clonal population was selected for further characterization. As seen in Fig. S1, this cell line had negligible ASIC1 immunoreactivity compared with either WT HEK293T cells or HEK293T cells transfected with human ASIC1a. Furthermore, whole-cell patch clamp recordings from these presumptive KO cells found no significant currents in response to pH 5.0 application. All other experiments in this study used this HEK293T ASIC1 KO cell line where endogenous ASIC1 has been removed.
To investigate the kinetic consequences of Q277G in cASIC1, we excised outside-out patches from HEK293T KO cells transfected with either WT cASIC1 or cASIC1 Q277G, along with eGFP. Patches were jumped from pH 8, to populate the resting state, into pH 6.5, 6, or 5.5, to activate and desensitize the channels. We were surprised to find that desensitization is completely intact in Q277G (Fig. 1 C). Indeed, the rate of desensitization was accelerated by about twofold (pH 5 desensitization time constant: WT, 180 ± 7 ms, n = 6; Q277G, 78 ± 8 ms, n = 7; P < 1e−6; Fig. 1 C). We also noted that there was a slight elevation of the steady-state or equilibrium current with pH 5.5 stimuli (Fig. 1, C and F; percent steady-state current: WT, 0.09 ± 0.03%, n = 6; Q277G, 2.0 ± 0.6%, n = 7; P = 0.005). To better compare with past work on hASIC1a (Wu et al., 2019), we used the same pH 6.5 stimulus in our cASIC1 experiments. Interestingly, the elevated steady-state current was more prominent with less acidic stimuli, increasing to ∼10% of the peak response using pH 6.5 (Fig. 1, D–F; percent steady-state current: pH 6.0, 1.8 ± 0.4%; pH 6.5, 10 ± 2%, n = 7, P = 0.0001). Such a pH-dependent increase in steady-state current was not detectable in WT channels, although the amplitudes of these steady-state currents are exceedingly small and hence difficult to measure (Fig. 1, D and F; percent steady-state current: pH 6.0, 0.15 ± 0.09%; pH 6.5, 0.06 ± 0.04%, n = 6, P = 0.32 versus pH 5.5 steady-state).
The robust desensitization of Q277G was unexpected given prior work. However, we did observe a small yet significant increase in the current at steady-state, particularly at more alkaline stimulating values (Fig. 1, E and F). We hypothesize that this phenotype arises from a weaker pH dependence of recovery from desensitization. ASIC recovery from desensitization is strongly dependent on the pH between the conditioning and test stimuli. Relatively alkaline inter-stimuli pH values accelerate recovery, while more acidic inter-stimuli pH values slow recovery (Alijevic et al., 2020; Gwiazda et al., 2015; MacLean and Jayaraman, 2016; Rook et al., 2020b). If one extrapolates this trend, then at more acidic values (i.e., pH 5.5) recovery is very slow, and transitions from the desensitized state to the open or resting states are very unfavorable. Consequently, there is minimal steady-state current. The elevated steady-state current of Q277G suggests that Q277G recovery may be faster than WT and/or less influenced by inter-stimuli pH values. To test this, we examined Q277G recovery from desensitization using several inter-stimuli pH values in the same patch. Consistent with our hypothesis, Q277G recovery from desensitization is substantially faster than WT cASIC1 (Fig. 2). Specifically, Q277G recovery had a time constant (and slope m) of 2.03 ± 0.05 ms (m = 12 ± 1, n = 5) at pH 8.0, which is roughly 400-fold faster than WT cASIC1 (840 ± 90 ms, m = 0.96 ± 0.05, n = 5, P < 1e−5; Rook et al., 2020b). Furthermore, the recovery time constants remained very fast at pH 7.4 and 7.0 (pH 7.4 τrec = 4.7 ± 0.2 ms, m = 6.1 ± 0.4; pH 7.0 τrec = 34 ± 2 ms, 1.8 ± 0.1, n = 5). Thus, these data support the notion that the elevated steady-state current in Q277G arises from faster recovery from desensitization in general.
SSD at any given pH value reflects a balance between channels entering and exiting the desensitized state. We hypothesized that the 400-fold faster exit from desensitization is a consequence of a less stable desensitized state. Another consequence of a less stable desensitized state is that a stronger pH stimulus is needed to induce desensitization. This could be observed as a right shift in the SSD curve. To test this, we constructed both activation and SSD curves of Q277G and WT cASIC1 in excised patches (Fig. 3). We found that the pH dependence of activation of Q277G was very slightly more alkaline compared with WT (WT pH50 = 6.51 ± 0.01, n = 5; Q277G pH50 = 6.55 ± 0.01, n = 6, P = 0.006; Fig. 3). However, the SSD of Q277G was considerably right-shifted (Fig. 3, B–D). Specifically, the pH50 of SSD shifted from 7.30 ± 0.01 in WT to 6.73 ± 0.01 in Q277G (n = 6 for both, P < 1e−5). The magnitude of the right shift was sufficiently large to induce overlap with the activation curve, leading to distinct standing currents with baseline pH values such as pH 6.8 or 6.4 (Fig. 3 B). Taken together, we have found that Q277G produces only a small reduction in desensitization (or enhanced steady-state current) while dramatically accelerating recovery from desensitization and right-shifting SSD curves with minimal effect on activation curves. We hypothesize that the recovery and SSD phenotypes all result from reducing the stability of the desensitized state. It has previously been suggested that the conformation of Gln276 (human ASIC numbering) controls the stability of the desensitized state by acting as a valve or steric barrier to regulate isomerization of the β11–12 linker (Wu et al., 2019). To gain insight into the structural mechanism, we turned to molecular dynamics simulations.
Examining the proposed desensitized state structure of cASIC1 suggested that the Gln277 side chain might form a hydrogen bond to the backbone oxygen atom of Leu414 when the linker is in the desensitized conformation (Fig. 4 A). Rather than acting as a valve, Gln277 could potentially stabilize the desensitized conformation through this hydrogen bond, as we proposed from previous simulations (Rook et al., 2020b). Additionally, three acidic residues, Glu80 and Glu417 in the lower palm domain and Glu412 in the upper palm domain, are within potential hydrogen bond distance of Gln277, partly depending on the protonation states of the acidic residues. Thus, the structure suggests that Glu277 could play a role in a larger hydrogen bond network (Fig. 4 A). Because protonation states of the acidic side chains cannot be observed but must be inferred, we first tested the relative stability of the potential hydrogen bonds in the presence of different protonation states of Glu80, Glu412, and Glu417. Each residue can be protonated or not, giving rise to eight possible protonation combinations. Using the desensitized state structure (PDB accession no. 4NYK; Gonzales et al., 2009) as a starting point, we simulated each protonation scheme for three repeats of 100 ns each. To quantify the stability of potential hydrogen bond interactions between the residues of interest, we measured the fraction of time that each potential hydrogen bond was present over the course of the simulations. Potential hydrogen bond donors considered were the side chains of Glu80, Gln277, Glu412, and Glu417, while potential hydrogen bond acceptors were the same side chains along with the backbone oxygen atom of Leu414. An interaction was considered as a hydrogen bond when the donor-acceptor distance was within 3.0 Å and the donor-hydrogen-acceptor angle >150°. The overall hydrogen bond analysis (Fig. S2) illustrated that no matter the protonation states, Gln277 very rarely acted as a hydrogen bond acceptor. On the contrary, Gln277 often participated as a hydrogen bond donor in fairly stable hydrogen bonds. Looking at all 72 chains analyzed (8 setups × 3 repeats × 3 chains), Gln277 formed hydrogen bonds of varying stability to Glu80 in 75% of cases, to Glu412 in 25% of cases, to L414 in 90% of cases, and to Glu417 in 35% of cases. Therefore, we deemed the hydrogen bonds to Glu80 and to Leu414 to be most important. These two hydrogen bonds showed the highest stability in the setup in which Glu412 and Glu417 were protonated while Glu80 was deprotonated (E80-/E412H/E417H in Fig. S2). Thus, we chose this protonation setup to be the most stable for the desensitized state. Under these conditions, the side chain conformation of Gln277 was generally stable and positioned to hydrogen bond with the side chain of Glu80 and the backbone carbonyl oxygen of Leu414 (Fig. 4 B and Video 1). These interactions are noteworthy as mutations of either Glu80 or Leu414 can profoundly alter desensitization kinetics (Cushman et al., 2007; Della Vecchia et al., 2013; MacLean and Jayaraman, 2017; Rook et al., 2020b; Roy et al., 2013). In particular, motion of Leu414 is a critical regulator of ASIC desensitization, underscoring the potential significance of these contacts. Fig. 4 and Fig. S3 illustrate this analysis, showing that Q277 spends considerable time in putative hydrogen bond interactions with both Glu80 and Leu414. We hypothesized that such a network stabilizes the desensitized state with Gln277 acting as a critical hub. This role of Gln277 as an electrostatic hub is in contrast to the purely steric “valve” model of Gln277 proposed previously (Wu et al., 2019). We reasoned that a Q277N mutation may delineate between these hypotheses. If the “steric” hypothesis is true, then shortening the side chain (Q277N) should produce only a small effect on desensitization kinetics. However, if the electrostatic hub model is more accurate, then the sub-optimal bonding distances of Q277N should result in much faster recovery from desensitization. To confirm that Q277N does attenuate hydrogen bond interactions, we repeated simulations using the Q277N mutation and observed that Q277N showed a greatly reduced capacity to participate in hydrogen bonds with Glu80 and Leu414 (Fig. 4, C–F; Fig. S3; and Video 2). Therefore, we measured the recovery from desensitization of Q277N in excised patches. Consistent with the electrostatic hub hypothesis, Q277N recovers from desensitization much faster than WT at all pH values tested. Specifically, at pH 8.0, the recovery time constant (and slope m) for Q277N is 4.0 ± 0.1 ms (m = 8.7 ± 0.6, n = 7, P < 1e−5 versus WT; Fig. 5). This is slowed to 32 ± 3 ms (m = 2.4 ± 0.1) and 1,500 ± 150 ms (m = 0.79 ± 0.07) at pH 7.4 and 7 (n = 6 and 5, respectively; Fig. 5). Next, we eliminated any residual capacity of the 277 position to participate in hydrogen bonds by using the Q277L mutation, which has identical steric factors as Q277N but no capacity for electrostatic interactions with nearby side chains. Consistent with the electrostatic hub hypothesis, Q277L has comparable recovery kinetics to Q277N (pH 8: τrec = 3.5 ± 0.1 ms, m = 4.3 ± 0.3; pH 7.4: τrec = 39 ± 2 ms, m = 2.23 ± 0.07; pH 7: τrec = 2,400 ± 90 ms, m = 0.88 ± 0.03, n = 5–7; Fig. 5, B–D). While these time constants are slower than Q277G, they are orders of magnitude faster than WT, suggesting the essential feature of Gln277’s function is as a hydrogen bond hub or coordinator and not a steric valve.
These data demonstrate that in cASIC1, Q277G does not block desensitization. Rather, Q277G induces a slight increase in steady-state current that is pH-dependent. Given this, we reexamined the Q276G mutation in human ASIC1a as was previously published (Wu et al., 2019), as well as mouse ASIC1a. In both cases, the Q276G mutation gave small, barely detectable currents in excised patches, necessitating whole-cell recording. In the case of mouse ASIC1a Q276G, even whole-cell currents were too small to resolve and examine convincingly (43 ± 13 pA, n = 10). Therefore, we confined ourselves to hASIC1a Q276G. WT hASIC1a contains a Gly at position 212, as does cASIC1, but many laboratories have been inadvertently using a rare human variant containing an Asp at that position (Vaithia et al., 2019). We measured the decay kinetics of both WT G212 and the rare variant D212. As expected (Vaithia et al., 2019), the D212 variant had faster desensitization kinetics (pH 5.5 desensitization time constant: G212, 790 ± 100 ms, n = 4; D212, 330 ± 20 ms, n = 5; P < 1e−6; Fig. 6), while both had comparable percent steady-state currents (pH 5.5 percent steady-state current: G212, 0.1 ± 0.2%, n = 4; D212, 0.6 ± 0.3%, n = 5; P = 0.27; Fig. 6). Interestingly, combining the Q276G mutation with either variant produced different results. In the G212 background, we observed rather fast rundown or inhibition when using a 5-s stimulus and 20-s intervals. This stimulus and interval duration have proven adequate for WT hASIC1a in our hands. To properly measure the desensitization time course and allow for complete recovery, we progressively extended both the stimulus and interval times. Ultimately, using a 100-s pH application spaced by 120 s, we found that hASIC1a Q276G channels desensitized very slowly (Q276G: 8,800 ± 1,400 ms, n = 5; WT: 790 ± 100 ms, n = 4; P < 1e−5). Q276G in the D212 background also required long intervals for recovery. It was impractical to fully map the recovery time course, but we estimate the pH 8.0 recovery time constant to be ∼60 s for G212 and 35 s for D212. Interestingly, and in contrast to the G212, D212/Q276G showed an acceleration of decay kinetics. Fitting the D212/Q276G decays required double exponentials where ∼80% of the current was fit by a rather fast component (150 ± 10 ms for pH 5.5, 517 ± 5 ms for pH 6.5, n = 5) and the rest with a slower component (6.0 ± 1.5 s for pH 5.5, 3.0 ± 0.4 ms for pH 6.5, n = 5). However, in both backgrounds, Q276G increased the steady-state current to ∼20% with pH 6.5 but only a few percent with pH 5.5 (Fig. 6), as was observed with cASIC1 Q277G. Taken together, we demonstrate that the Q/G mutation does not abolish desensitization as previously reported. Rather, in cASIC1 this mutation elevates the steady-state current, accelerates recovery from desensitization, and reduces the stability of the desensitized state. Molecular dynamics simulations and subsequent mutagenesis suggest that these phenotypes arise by destabilizing a critical hydrogen bond network, which in the WT stabilizes the desensitized state. In hASIC1a, this mutation also does not abolish desensitization and either slows or accelerates decay kinetics, depending on the variant.
Discussion
We explored the properties of the recently described Q276G hASIC1a mutation (Wu et al., 2019; human numbering) using a combination of fast perfusion electrophysiology and molecular dynamics simulations. In contrast to prior work on hASIC1a, we find that this mutation does not abolish cASIC1 desensitization. Rather, this mutation leads to a slight elevation in steady-state current that is more pronounced with weaker pH stimuli (Fig. 1). In cASIC1, Q277G also markedly accelerates recovery from desensitization over a wide pH range (Fig. 2) and right-shifts the pH dependence of SSD without substantially altering activation (Fig. 3). All-atom simulations of the cASIC1 desensitized state indicate that this conformation is stabilized by a network of hydrogen bonds linking the lower palm residue Glu80, through Gln277, with the β11–12 linker (Fig. 5). It is worth noting that our analysis focused on interactions centered around Gln277. Thus, we cannot exclude that effects stemming from longer-range electrostatic interactions within or between subunits also affect this hydrogen bond network and thus the stability of the desensitized conformation. Consistent with the hydrogen bond analysis, compromising the hydrogen bond network by slightly shortening the Q277 side chain (Q277N) or eliminating the capacity for hydrogen bonding (Q277G and Q277L) has a profound impact on the stability of the desensitized state as measured by recovery from desensitization (Fig. 5). Finally, we found that in hASIC1a, Q276G also desensitizes but enters the desensitized state slower or faster depending on the hASIC1a variant (Fig. 6).
Comparison with previous studies
The original report that Q276G blocks desensitization used human ASIC1a in a Xenopus laevis oocyte expression system primarily using bath perfusion, pH 6.5 as a stimulus with pH 7.4 as a baseline pH (Wu et al., 2019). Oocyte experiments have an inherently slower perfusion speed. Moreover, if recordings were obtained from the faster desensitizing phenotype of hASIC1a D212 Q276G, the peak response might have been missed or at least insufficiently resolved. Thus, the experimenter would predominantly observe the larger steady-state current produced by pH 6.5, leading to the conclusion that desensitization has been blocked.
Rather than Gln277 controlling desensitization and recovery by the proposed valve mechanism (Wu et al., 2019), we provide evidence that Gln277 is central to an important hydrogen bond network linking the influential Glu80 residue in the lower palm with the critical β11–12 linker that governs desensitization. How might such a network function in the ASIC gating cycle? In our simulations, both Glu412 and Glu417 are protonated, leaving the Gln277 amide to act as a hydrogen bond donor to the deprotonated Glu80 and the backbone carbonyl of Leu414. The interaction with the carbonyl is the most commonly observed (Fig. 4, Fig. S2, and Fig. S3). We propose that in the desensitized state, Gln277 partly contributes to the stability of Leu414 by this hydrogen bond, with Gln277 being held in this advantageous position by Glu80. Upon alkalization, Glu412, Glu417, or both tend to become deprotonated, acting as alternative hydrogen bond acceptors and thereby helping to pull the amide group of Q277 away from the backbone carbonyl of Leu414, releasing Leu414. This would enable the swivel of the β11–12 linker and thus recovery from desensitization. However, it is difficult to reconcile this hypothesis with the hASIC1a Q276G data that show an apparent slowing of recovery from desensitization.
Human versus chicken data
In our hands, the Q/G mutant gives distinct effects in cASIC1 versus hASIC1a. In every case, the channels still desensitize quite strongly when stimulated with pH 5.5. However, in the G212 (or equivalent) background, cASIC1 and hASIC1a show opposing effects of the Q/G mutation on desensitization kinetics, accelerating decays in cASIC1 but slowing them in hASIC1a (Fig. 1 versus Fig. 6). This is reminiscent of the effects of psalmotoxin, which inhibits mammalian ASICs by stabilizing a desensitized state (Chen et al., 2005) yet activates cASIC1, promoting an unusual nonselective open state (Baconguis and Gouaux, 2012; Samways et al., 2009). Another recent example is the blunted effect of mambalgin in cASIC1 compared with hASIC1a, which can largely be reversed by several point mutations (Sun et al., 2020). Presently it is unclear what the source of these differences is. Human and chicken ASIC1 contain 56 amino acid differences, 31 of which are in the extracellular domain. A number of these are concentrated in the wrist region, including a 2–amino acid insertion. Given the wrist region’s involvement in gating (Li et al., 2009), it is possible that many species-specific differences arise from here. Further differences relevant for our kinetic experiments include the TRL versus SQL substitutions around amino acids 84–86 (Coric et al., 2003) as well as Ser275Ala, Val368Leu, and Ala413Val, which are all relatively proximal to Gln277 (chicken to human differences). We hypothesize that one or more of these changes subtly alter the structure of hASIC1a, potentially imparting distinct pKa values on critical palm residues and thus changing the phenotype of Q276G. As more hASIC1a structures become available in distinct functional states (Sun et al., 2020), we hope to explore the source of these differences and the conservation of mechanisms in more detail. A similar examination may uncover why the equivalent Q269G mutation in ASIC3 does appear to inhibit desensitization even with pH 5.0 (Klipp et al., 2020). Regardless of phenotypic differences, our data clearly indicate that both cASIC1 and hASIC1a Q277G mutants desensitize to a large extent. Therefore, using these mutations to explore either biophysical mechanisms of desensitization, or its physiological consequences, may be problematic.
Acknowledgments
Christopher J. Lingle served as editor.
We thank Dr. Matthieu Chavent for initial discussions concerning the hydrogen bond analysis.
Funding for this work was provided by National Institutes of Health grant T32GM068411-15 and a Joan Wright Goodman Fellowship to M.L. Rook; National Science Foundation Graduate Reseach Fellowship Program to T. Couch; an National Science and Engineering Research Council Discovery Grant (RGPIN 2019-06864) and Canada Research Chairs grant (950-232154) to M. Musgaard; and National Institutes of Health grants R00NS094761 and R35GM137951 and a National Alliance for Research on Schizophrenia and Depression Young Investigator Award (27518) to D.M. MacLean.
The authors declare no competing financial interests.
Author contributions: M.L. Rook, M. Miaro, D.L. Kneisley, T. Couch, and D.M. MacLean conducted experiments and analyzed data. M.L. Rook, M. Miaro, M. Musgaard, and D.M. MacLean interpreted results and edited the manuscript.