# Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the Gloeobacter violaceus pentameric ligand-gated ion channel

Ákos Nemecz, Haidai Hu, Zaineb Fourati, Catherine van Renterghem, Marc Delarue, Pierre-Jean Corringer

► To cite this version:

Ákos Nemecz, Haidai Hu, Zaineb Fourati, Catherine van Renterghem, Marc Delarue, et al.. Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the Gloeobacter violaceus pentameric ligand-gated ion channel. PLoS Biology, 2017, 15 (12), pp.e2004470. 10.1371/journal.pbio.2004470. pasteur-01721168

HAL Id: pasteur-01721168

https://pasteur.hal.science/pasteur-01721168

Submitted on 1 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the *Gloeobacter violaceus* pentameric ligand-gated ion channel

Ákos Nemecz¹, Haidai Hu², Zaineb Fourati², Catherine Van Renterghem¹, Marc Delarue², Pierre-Jean Corringer¹*

1 Unité Récepteurs-Canaux, Unité Mixte de Recherche 3571 du Centre National de la Recherche Scientifique, Institut Pasteur, Paris, France, 2 Unité de Dynamique Structurale des Macromolécules, Unité Mixte de Recherche 3528 du Centre National de la Recherche Scientifique, Institut Pasteur, Paris, France

* pierre-jean.corringer@pasteur.fr

<figure><img src="image_2.png" /></figure>

## Abstract

The *Gloeobacter violaceus* ligand-gated ion channel (GLIC) has been extensively studied by X-ray crystallography and other biophysical techniques. This provided key insights into the general gating mechanism of pentameric ligand-gated ion channel (pLGIC) signal transduction. However, the GLIC is activated by lowering the pH and the location of its putative proton activation site(s) still remain(s) unknown. To this end, every Asp, Glu, and His residue was mutated individually or in combination and investigated by electrophysiology. In addition to the mutational analysis, key mutations were structurally resolved to address whether particular residues contribute to proton sensing, or alternatively to GLIC-gating, independently of the side chain protonation. The data show that multiple residues located below the orthosteric site, notably E26, D32, E35, and D122 in the lower part of the extracel-lular domain (ECD), along with E222, H235, E243, and H277 in the transmembrane domain (TMD), alter GLIC activation. D122 and H235 were found to also alter GLIC expression. E35 is identified as a key proton-sensing residue, whereby neutralization of its side chain carbox-ylate stabilizes the active state. Thus, proton activation occurs allosterically to the orthos-teric site, at the level of multiple loci with a key contribution of the coupling interface between the ECD and TMD.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Nemecz Á, Hu H, Fourati Z, Van Renterghem C, Delarue M, Corringer P-J (2017) Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the *Gloeobacter violaceus* pentameric ligand-gated ion channel. PLoS Biol 15(12): e2004470. https://doi.org/10.1371/journal.pbio.2004470

**Academic Editor:** Raimund Dutzler, University of Zurich, Switzerland

**Received:** October 9, 2017

**Accepted:** December 4, 2017

**Published:** December 27, 2017

**Copyright:** © 2017 Nemecz et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files. Crystallographic structures have been deposited with the RCSB protein data bank. The PDB ID's are: 6F0I, 6F0J, 6F0U, 6F0M, 6F0V, 6F13, 6F15, 6F0N, 6F0R, 6F11, 6F10, 6F0Z, 6F12, 6F16 for the mutations: E26A, E26Q, E35A, E35Q, E67A, E75A, H127Q, E82A, E82Q, D86A, D88A, D88N, E181A, H277Q. Their respective PDB accession ID's are also listed in S2 Table.

## Author summary

Pentameric ligand-gated ion channels are an important class of receptors that are involved in many neurological diseases. They have been extensively studied but a full understanding of their mechanism of action has yet to be achieved. In an effort to bypass obstacles in the research of human receptors, bacterial versions have been used to characterize the family's structure-function relationship. One key bacterial receptor, known as GLIC, has
**Funding:** Agence Nationale de la Recherche http://www.agence-nationale-recherche.fr/ (grant number ANR-13-BSV8-0020). Received by ÁN and ZF. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** ASIC, acid-sensing ion channel; DDM, n-dodecyl-β-D-maltoside; ECD, extracellular domain; ESRF, European Synchrotron Radiation Facility; GE, GLIC<sub>ECD</sub>-ELIC<sub>TMD</sub> mutant; GELIC, GLIC<sub>ECD</sub>-ELIC<sub>TMD</sub> mimicking Lily; GLIC, *Gloeobacter violaceus* ligand-gated ion channel; GlyR, glycinergic receptor; MBP, maltose-binding protein; nAChR, nicotinic acetylcholine receptor; OD, optical density; PDB ID, protein data bank identification code; pLGIC, pentameric ligand-gated ion channel; RMSD, root-mean-squared deviation; TBS, tris-buffered saline solution; TMD, transmembrane domain; Wt, wild-type.

lead the way in structural resolution of various mechanistic states along the gating pathway, yet its activation by protons is significantly less understood than its human counterparts. To define the site(s) involved in proton gating, we systematically mutated all titratable residues near the pH<sub>50</sub> of activation: Asp, Glu, and His. We determined that a previously established His residue in the transmembrane domain is structurally important but likely plays little or no role in proton gating. We instead found that proton activation is a complex multiple loci mechanism, with the key contribution stemming from the coupling interface between the extracellular and transmembrane domain, with E35 acting as a key proton-sensing residue.

## Introduction

Pentameric ligand-gated ion channels (pLGICs) are key players of neuronal communication. They promote either cell depolarization or hyperpolarization with the passive permeation of ions through an intrinsic channel, whose opening is stabilized by the binding of specific neurotransmitters. pLGICs are ubiquitously expressed in virtually all neurons, contribute to central nervous system functions, including sensory and motor processing, central autonomous control, memory and attention, sleep and wakefulness, reward, pain, anxiety, emotions, and cognition [1]. As such, they are important drug targets. At least one member of each major subfamily of vertebrate pLGICs has been structurally resolved in the recent years: a serotonergic receptor (5-HT<sub>3A</sub>, [2]) and a nicotinic acetylcholine receptor ([nAChR] α4β2-nAChR, [3]) from the cationic receptors, and a GABAergic receptor (β3-GABA<sub>A</sub>, [4]) and two glycinergic receptors ([GlyRs] α1- and α3-GlyRs, [5,6]) from the anionic receptors. However, relating these 3D structures to the physiologically relevant allosteric states that mediate pLGIC activation and desensitization remains an open and debated question [1].

At present, the best structurally characterized pLGIC is the *G. violaceus* ligand-gated ion channel (GLIC), of prokaryotic origin. Its structure shows a highly conserved fold with the subsequently resolved structures in the pLGIC family. Each subunit consists of an extracellular domain (ECD), predominantly in a β-sandwich fold, and a transmembrane domain (TMD) composed of four helices, labeled M1–M4. The available structures of pLGICs support that the GLIC globally shares a common gating mechanism with its eukaryotic cousins, although molecular details differ [1].

Due to the relative ease for overexpression in *Escherichia coli*, as well as its biochemical robustness toward detergent-solubilization and mutations, the GLIC has been resolved in four distinct conformations. It is the first receptor to be resolved in an apparently open channel confirmation, as well as three various closed channel conformations [7–9]. In addition, the GLIC has been solved as a complex with a variety of allosteric modulators, such as barbiturates, bromoform, lidocaine, propofol, and xenon [10–14]. Membrane-inserted GLIC proteins have been studied by electron paramagnetic resonance spectroscopy and fluorescence-quenching experiments following site-directed labeling, relating some local conformational changes to the activation and desensitization transitions monitored by electrophysiology [8,15–18]. The GLIC has also been studied by computational methods including molecular dynamic simulations [19–23]. This combined set of data support that, upon activation, the subunits’ ECD regions move closer together, which precedes a concerted tilt of the pore-lining M2-α-helices to open the pore gate and activate the receptor.

Chimeras between the ECD of the GLIC and the TMD of the α1-GlyR were additionally shown to fold properly and to be functional, revealing compatibility between prokaryotic and
eukaryotic domains [24,25]. However, the GLIC is a pH-gated channel, activated by lowering the pH, with a maximal activation at pH 4 and a pH₅₀ around 5 [26]. This sharply contrasts with most eukaryotic pLGICs, which are activated by a neurotransmitter binding to a well-described cavity in an intersubunit interface of the ECD. Mammalian pLGICs are not directly activated upon pH changes, but the agonist-elicited responses are modulated by pH, notably for the α3β4-, α3β2-, and α4β2-nAChRs [27], α1-GlyR [28,29], and various GABA<sub>A</sub> receptors [30,31]. So far, a handful of invertebrate pLGICs were found to be directly activated by pH, a nAChR from *Caenorhabditis elegans* by low pH [32], and two insect GABA<sub>A</sub> receptors by high pH [33,34].

Fully understanding the molecular mechanism of GLIC signal transduction thus requires identifying the locus where protons act to activate the channel and which titratable groups are crucial in this process. The observation that a chimera composed of the ECD of the GLIC fused to the TMD of the α1-GlyR (GLIC<sub>ECD</sub>-GlyR<sub>TMD</sub> known as Lily [25]) or the *Erwinia chrysanthemi* ligand-gated ion channel (GLIC<sub>ECD</sub>-ELIC<sub>TMD</sub> [35,36]) is activated by protons indicates that the major proton-sensing motifs are likely situated in the ECD. However, an inverse chimera, the ELIC<sub>ECD</sub>-GLIC<sub>TMD</sub>, yields pH-gated currents when the gain-of-function mutation of I9' A is incorporated, suggesting either proton sensors too weak to activate in the ELIC ECD or their presence in the GLIC TMD [35,37]. Mutation of a His residue located in the middle of the TMD (His235) abolishes GLIC function, raising the possibility that it might be involved in pH sensing [38,39], but combination of this mutation with a gain of function mutation restores at least part of the pH-gated function [8,37]. Altogether, the location of the proton activation site(s) remain(s) essentially unknown thus far.

The present study provides an exhaustive mutational analysis of GLIC titratable residues, with pKa's in the pH 6 to pH 4 range (namely Asp, Glu, and His), given the pH₅₀ of GLIC activation, to search for the proton-sensing site(s). Each GLIC subunit contains 19 Asp, 16 Glu, and 3 His residues that were individually and/or collectively mutated. Mutation of an Asp/Glu/His residue may be expected to affect the GLIC function in two ways: (1) by altering direct proton sensing when the protonated form of this residue stabilizes the active state, as compared to the nonprotonated form, or (2) by altering the gating equilibrium independently of side chain protonation. In an effort to discriminate between these two possibilities, Asp/Glu residues were systematically mutated to Asn/Gln, replacing the carboxylate moiety by a nontitratable amide group, thereby tentatively mimicking a permanently protonated form, as well as to Ala, removing the titratable moiety. Selected mutants were also solved by X-ray crystallography to characterize the effect of the mutations on the local protein structure.

To identify the determinants of proton-modulation/gating, two-electrode voltage clamp electrophysiology using expression in *Xenopus laevis* oocytes was employed; activation was elicited by dropping the pH from neutral (pH 7.3-8) to lower values (minimum pH 3.7). All recordings showed a slow onset of the response, with no apparent desensitization during a 30-90 s pH application. Distinguishing the kinetics of activation and desensitization as compared to the wild-type (Wt) GLIC, proved difficult for the vast majority of mutants, therefore mutants were characterized on the basis of their pH₅₀, defined as the pH eliciting half of the maximal current.

To minimize the influence of the intrinsic variability between oocyte batches, which show some variation in the Wt response to pH changes, mutants were also characterized by a ΔpH₅₀, which corresponds to the variation of pH₅₀ between each cell expressing a mutant, and the Wt cell(s) in the same batch of oocytes. Significance of the results in this report was determined as values larger than 0.5 pH units for the mean ΔpH₅₀, and mean pH₅₀, as compared to Wt, based on the standard deviation of 83 measurements of the Wt (with a standard deviation of 0.4 pH units). All values that fulfilled this criterion of significance also had a *p* value < 0.01
<figure><img src="image_4.png" /></figure>

**Fig 1.** Overview and typical current traces. (A) Three subunit view of the pH 4, apparently open-form, structure of GLIC (PDB ID: 4HFI), with solid black bars representing the membrane level. The central subunit is shown in black ribbon and transparent cartoon representation and the neighboring subunits are in white ribbon. Titratable Asp/Glu/His residues are shown as Cα spheres, color-coded based on mutation results with pH₅₀ decreasing, increasing, and no significant change colored red, blue, and yellow, respectively. The outlined boxes represent the regions of interest (TMD and ECD: apical, basal complementary [–], and principal [+] faces, along with Loop C) present in subsequent figures. (B) Typical traces of proton-elicited currents from prominent and interesting mutants. Traces are aligned and scaled in time to match the Wt trace. Vertical scale bars are identical, whereas horizontal scale bars vary depending on perfusion time required to achieve maximal currents. The dashed line represents 0 µA current. Loss of function mutants (smaller pH₅₀) generally had slower rising phases requiring longer activation perfusions. For example, with the quadruple mutant, the current did not reach a maximum within two min of perfusion at any pH value. Some gain-of-function mutations on the other hand, required longer rinse times to recover the baseline. (C) Normalized dose-response curves for the corresponding mutants represented in (B). Plots are mean ± standard deviation for I/Imax of all recorded values noted in S1 Table for the given mutant, where Imax is the average absolute current for the maximal pH value per oocyte tested. ECD, extracellular domain; GLIC, *G. violaceus* ligand-gated ion channel; I, current (µA); PDB ID, protein data bank identification code; TMD, transmembrane domain; Wt, wild-type.

https://doi.org/10.1371/journal.pbio.2004470.g001

in the Student *t* test against the Wt. A pH₅₀ could not be established for some mutants, often with minimal currents.

For clarity, the presentation of the results is organized according to five regions of a GLIC subunit: the apical top of the ECD, the Loop B and C on the principal (+) face of the interface, the basal ECD principal (+) and complementary (−) faces, and finally the TMD (Fig 1).

## Results

### Apical ECD mutations weakly affect proton sensitivity/modulation

In the apical region, the residues either face the solvent (D13, E14, D55, E67, E69, D97), are located at the subunit interface (D49, E75, D136), or are located at the middle of the ECD β-sandwich (E104). Their mutation produces weak effects. D13N-E14Q, D49N-D55N, E67Q-E69Q, and E67Q-E75Q produced no significant change in pH₅₀, and neither did the corresponding single mutants of D55N, E67Q, E67A, E69Q, and E75Q. The mutants D49N,
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Apical residues do not have an effect on the GLIC pH<sub>50</sub>. Stick representation of the evaluated Asp and Glu residues are labeled and colored yellow for nonsignificant or weak effects. The main subunit is shown in black with ribbon and cartoon representation, whereas neighboring subunits are represented in white ribbon. Table of the pH<sub>50</sub> and ΔpH<sub>50</sub> is shown below with mean values ± standard deviation for evaluated mutants with recorded values noted in S1 Table. Light red and green zones with dashed lines represent the respective ½ log decreases and increases from Wt. Values outside this region are color coded for ease of interpretation with red and green respectively. GLIC, *G. violaceus* ligand-gated ion channel; Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g002

E75A, D97N, D136N, and D136A display small increases in pH<sub>50</sub>, but do not meet the criterion for a significant effect. The effect of D49 replacement is corroborated in a recent study, whereas the same study showed D136A to have a decreased pH<sub>50</sub> with an increased Hill-slope [36].

Finally, E104, which makes internal intrasubunit contacts, showed no significant change when mutated to Q. Mutation to A, meanwhile, produced a significant ΔpH<sub>50</sub> and a noticeable, yet nonsignificant, decrease in pH<sub>50</sub> (Fig 2).

## A complex network in Loop C modulates proton sensitivity

The β9-β10-loop (Loop C) was investigated in conjunction with R133 (Loop B), both of which are located on the principal (+) intersubunit interface, a region that contributes to neurotransmitter binding in eukaryotic pLGICs. The mutation D185N at the base of Loop C did not have a significant effect (Fig 3). This position was also recently reported with Wt-like properties as both Asn and Ala mutations [36]. Removing all remaining titratable moieties in the sextuple mutant R133A-E177Q-D178N-R179Q-E181Q-K183Q produced a marked increase in pH<sub>50</sub>. In contrast, mutating only the three acidic residues
<figure><img src="image_6.png" /><figcaption>**Fig 3. Loop C residues modulate activation through a complex network of interactions.** Stick representation of the evaluated Arg, Asp, Glu, and Lys residues are labeled and color coded based upon effect, with Asp/Glu pink residues producing a strong decrease in pH₅₀ when mutated in combination, and yellow residues producing a nonsignificant or weak effect. R179 is shown in teal. The main subunit is shown in black with ribbon and cartoon representation, whereas the complementary neighboring subunit is represented in white ribbon, with important motifs/β-sheets labeled. Distances of intrasubunit contacts are shown between D136 and R179, as well as R133 and E181. Table of the pH₅₀ and ΔpH₅₀ are shown below with mean values ± standard deviation for evaluated mutants with recorded values noted in S1 Table. Light red and green zones with dashed lines represent the respective ½ log decreases and increases from Wt. Values outside this region are color coded for ease of interpretation with red as decreases and green and dark green as varying degrees of increases. Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g003

(E177Q-D178N-E181Q) produced a marked decrease in pH₅₀, whereas the individual Asp/Glu residue mutations show Wt-like responses with E181Q having a tendency to decrease the pH₅₀. The mutations E177A and E181A were previously reported in a cinnamic acid study with pH₅₀ values of 5.4 ± 0.1 and 5.6 ± 0.1, respectively, as compared to a Wt value of 5.2 ± 0.1 [40]. D178 was also recently reported as displaying a Wt phenotype for both the Asn and Ala mutations [36]. Evaluation of the basic residues shows R179 mutation increases the pH₅₀, with R179Q producing a marked increase, whereas R179A only has a tendency to increase the pH₅₀, showing that the removal of the side chain guanidinium is not solely responsible for the phenotype at this position. Meanwhile, replacing the remaining basic residue of Loop C, K183Q, has no effect. The single mutation R133A on Loop B also has no effect, which was also previously shown in [40]. Despite containing the triple mutation that results in a decrease in pH₅₀, the mutant that removes all titratable residues from Loop C along with R133A has a similar but stronger change in pH₅₀ as compared to the single R179Q mutant, yet their ΔpH₅₀s vary by 0.5 units. This suggests that the main player in this region is R179, along with a complex network of other side chain interactions modulating GLIC activation at this level.
## Mutations of the complementary (–) face of the ECD tend to decrease proton sensitivity

The complementary (–) face of the ECD contains eight Asp/Glu residues at the intersubunit interface and two solvent-exposed residues (D161 and E163) at the bottom of β9 near the pre-M1-π-helical loop [41]. The double mutants of pairs of proximal residues were initially tested: D86N-D88N, D145N-E147Q, D153N-D154N, and D161N-E163Q. The mutant of the vestibular pair on the β5 strand, D86N-D88N, yields a significant decrease in pH₅₀ greater than 1 pH unit. The solvent-facing pair D161N-E163Q, and the pair near the bottom of the β8-β9-loop (Loop F) D153N-D154N, tend to both decrease the pH₅₀, but below the threshold of significance. The pair D145N-E147Q, on Loop F, has no effect.

All residues were also tested individually, with most mutations having no significant effect, although a large majority tend to decrease the pH₅₀ (notably D86N/A, D88N/A, and D91N). Although D91A was found to insignificantly increase the pH₅₀ as compared to Wt, D91N and D91A were recently reported to have a slight decrease in pH₅₀ values [36]. In contrast to the nonsignificant individual mutations, E26, which is found interacting with the bottom of the β4-β5 loop of the principal (+) face, shows a robust decrease of the pH₅₀ when mutated to both Q and A.

A combination of E26Q with mutations of vestibular-facing residues that tend to decrease the pH₅₀ produces an obvious decrease in pH₅₀. The E26Q-D86N-D88N-D91N mutant exhibits strongly reduced maximal currents, which prevented reliable evaluation of the pH₅₀ in most of the cells (Figs 1B and 4A). Of the 11 oocytes tested (over four injections), only three recordings had enough current to properly evaluate (S1 Table).

It is striking that among the 23 mutants tested in this region, most show a propensity to decrease the pH₅₀ (Fig 4A).

## A key mutation at the principal (+) face of the ECD increases proton sensitivity

The basal principal (+) face, below Loop C, contains E82 and H127, along with a cluster of Asp/Glu residues found close to the TMD: D31, D32, and E35 from Loop 2, as well as D115 and D122 from Loop 7.

Consistent with some of the previous studies of H127 [38], H127N and Q both produce Wt-like currents (Fig 4). A structure was resolved for H127N, which resulted in no appreciable difference from the open form of the GLIC (protein data bank identification codes [PDB IDs]: 4HFI/3EAM, S1 Fig).

The effects of D32 and D122, both of which are engaged in salt bridges with R192, have been previously studied in other pLGICs [42], as well as in the GLIC [22,43]. D32A and D122A were shown to be nonfunctional yet weakly expressed [43], whereas D32N shows marked decrease in pH₅₀ and in maximal currents (Fig 4B, [22]). D122N showed no function, and subsequent expression studies showed no expression (Figs 4B and 5). D31, which points towards the vestibule, has no effect when mutated to D31N in this study, whereas a slight loss of function, insignificant with the criterion and to the Wt value of this study, was observed with both Ala and Asn mutations by Alqazzaz et al. [36].

Mutants of the solvent-exposed D115, further away from this triad on Loop 7, exhibit Wt-like properties. Yet, the D115N and D115A mutations show opposite phenotypes, with Asn having a tendency to increase the pH₅₀ and Ala to decrease it. Of the mutants which had both Asn/Gln and Ala mutations tested, only the pair of E82Q and E82A follows the same pattern (Fig 4B).

Finally, mutation of E35 produced the strongest effect; therefore, a more extensive study of this position was performed. The mutations E35Q, E35A, and E35M (the GlyR-position
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Basal ECD mutations have the greatest influence on activation. Stick representation of the evaluated Asp, Glu, and His residues are labeled and color coded based upon effect, with residues in red producing significant decrease in pH₅₀ as a single mutation, in pink a significant decrease in combination, in yellow an insignificant or weak effect, and in blue a significant increase in pH₅₀. The main subunit is shown in black with ribbon and cartoon representation, whereas the neighboring subunits are represented in white ribbon. Tables of the pH₅₀ and ΔpH₅₀ values are shown below with mean values ± standard deviation for evaluated mutants with recorded values noted in S1 Table. Light red and green zones with dashed lines represent the respective ½ log decrease and increase from Wt. Values outside this region are color coded similarly with darker bars representing stronger shifts. Mutants with too few replicates as a result of diminished maximal currents to discern an accurate pH₅₀ have the error bars removed. (A) Complementary inter-subunit residues. Almost all mutated residues decrease with respect to Wt, with only E26 having a significant loss of function shift. D86 and D88, colored in pink, produce a significant decrease in combination. (B) Principal and ECD-TMD interface residues. The mutation of H127, shown in orange, has no effect, whereas E82 is shown in teal to identify the residue. Hydroxyl-containing residues on the complementary subunit that were evaluated are also shown with the only significant residue, Y28, shown in cyan. D122N is bolded and in red as it was nonfunctional. ECD, extracellular domain; TMD, transmembrane domain; Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g004
<figure><img src="image_8.png" /><figcaption>**Fig 5. Protein expression of nonfunctional mutants.** Rabbit anti-HA tag Alexa Fluor-695 immunostaining results of the mutants that produced no detectable current in functional tests, compared to Wt and GFP-alone injected oocytes. Left, colored merge of GFP and Alexa-695 and right, the Alexa-695 imaging alone. Both D122N and H235A show no expression. GLIC, *G. violaceus* ligand-gated ion channel; GFP, green fluorescent protein; HA, human influenza hemagglutinin; Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g005

equivalent) display a significant increase in pH₅₀, whereas E35K shows a weaker nonsignificant effect, and E35H appears to show the inverse phenotype with a marked decrease in pH₅₀. Interestingly, fitting E35H data with a single sigmoidal curve yielded poor fits, preventing reliable calculation of the pH₅₀. Of the 13 recorded oocytes (over four injections), only two recordings could be fit properly (Figs 1B, 1C and 4B, S1 Table). The poor fit of E35H indicates a more complex mechanism at play, which could possibly be elucidated with a more in-depth study of mutation of this position. Altogether, these data show that the loss of charge at the E35 position is important to the gain-of-function effect found.

Three hydroxyl-containing residues that are situated near E35 and E82 on the complementary face (–) of the ECD were subsequently mutated: Y28, S29, and T158. Individually removing each hydroxyl group shows a tendency to increase the pH₅₀, with Y28F (near E82) having a significant increase. The mutation of T158A did not have the same marked effect, although T158 is the residue positioned closest to E35, indicating that the global environmental change in charge at position 35 is more important than the direct residue–residue interactions.

Subsequently, the additivity of the Asp/Glu mutations was also tested; however, neither E82Q-D115N, E35Q-E82Q, nor E35Q-E82Q-D115N showed increased proton sensitivity as compared to E35Q. A combination of all other weak pH₅₀-increasing mutations was performed with the addition of R179Q or E75Q-D97N-D136N to the triple mutant E35Q-E82Q-D115N, with neither of these mutants producing any shift greater than E35Q alone (Fig 4B).

The pH₅₀ evaluation of this region effectively identified key residues involved in pH-sensing, but the results of the combined mutants demonstrate the complexity of the mechanism involved in GLIC-gating.
## Crystal Structures of ECD mutants show no significant change, with the exception of E35A and E82Q

A structural analysis was performed on E26Q, E26A, E35Q, E35A, E67A, E75A, E82Q, E82A, D86A, D88N, D88A, H127N, E181A, and H277Q mutants, by solving their structure at pH 4. The respective PDB IDs and crystallographic statistics are listed in S2 Table. All structures were in the apparently open conformation. Each residue's root-mean-squared deviation (RMSD) and Cα RMSD was evaluated in relation to the intrinsic variability between subunits and across the two Wt structures (3EAM and 4HFI) at pH 4.6 and pH 4, respectively, as a means to control for crystal variability. There are no significant differences, measured as 5-fold over Wt, seen in the backbone Cα residues of any of the resolved structures, with the exception of E35A and E82Q (Fig 6 and S1 Fig). The E82Q side chain takes on a different rotamer, causing a local change in the backbone, as well as a cascading rotamer/conformational change in Y28 and F156, the latter of which flips out towards the aqueous environment. Interestingly, E35A also produces the greatest structural deviations around Y28 along with neighboring residues. Of the few conformational changes seen, it is interesting to note that Y28F had the greatest functional effect of the nontitratable mutations tested.

## Mutations at the TMD decrease proton sensitivity

It has been proposed that H235 is a key residue mediating proton sensing [38,39], but recent data show that the GLIC can still gate when H235 is mutated to nontitratable residues in combination with a strong gain-of-function mutation [8,37]. It has also previously been shown that GLIC expression in bacteria is highly sensitive to mutation at H235 [7]. Both of these findings are confirmed, firstly with the mutation H235A, which abolishes expression in oocytes (Fig 5), indicating a structural importance of the residue, and secondly, with the mutation H235Q, which does actually produce functional receptors, albeit with a strong decrease in pH₅₀ and reduced maximal currents (Fig 7, S1 Table). H235Q does not allow for a charge or a change in protonation state at this position, and therefore one would expect a completely nonfunctional receptor if this residue were solely responsible for the proton-modulated gating of the GLIC. Rather, this finding confirms the structural importance of H235 that has been previously reported, along with the importance of a hydrogen-bonding network, between the M2 α-helix and the neighboring M3 α-helix, with the backbone carbonyl of I262 [7,38].

Among the other titratable residues of the TMD, neither a mutation of E272Q nor E282Q produced an overall significant effect, albeit a significant negative ΔpH₅₀ was observed for E272Q, which was previously reported to be nonfunctional when mutated to A [44]. In contrast, the mutation of either E222 or E243 as Q or A, as well as H277Q, all produced significant loss-of-function effects, pointing to the key role of these residues in activation. E222 (E-2'), at the beginning of the M2 α-helix, has been extensively studied as the key component of the selectivity filter, a charged constriction point in the lower part of the pore, for cationic pLGICs [45,46]. The E222A mutation has previously been resolved crystallographically and has no change in structure as compared to Wt at pH 4 [10]. Therefore, the strong loss-of-function mutation does not appear to affect conformation. E243 flanks the apical section of the TMD pore, making a pentameric ring at the end of the M2 α-helix. A mutation of E243C was shown to have a similar pH₅₀ as Wt [18], whereas the mutation E243P is reported as nonfunctional and found in a locally closed conformation [7]. Mutation of H277, which lies in the adjacent M3 α-helix (nearby E222) with possible electrostatic interaction, has also previously been shown, using noncanonical amino acid substitution to decrease the pH₅₀ [38].
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Crystallographic structural evaluation. (A) X-ray crystallographic resolved structures of E82Q (PDB ID: 6F0R) and E82A (PDB ID: 6F0N) mutants in green and yellow, respectively, overlaid with the apparently open form GLIC-structure (PDB ID: 4HFI, black), with the principal (+) and complementary (−) subunits indicated. Residues that change conformation are labeled and shown as sticks, with the 2mFo-DFc map carved at 1σ around the E82Q structure residues. (B) Normalized Cα RMSD for the two structures that have a ΔCα/ΔCα<sub>ctrl</sub> greater than 5-fold: E35A (cyan) and E82Q (purple), in comparison to the 4HFI reference (Wt green). Both structures have significant deviation around Y28, whereas E35A has additional deviation with the rest of Loop 2. The important principal (orange) subunit loops are indicated along with Loop F and the Pro-Loop (gray) of the ECD, as well as the M2 α-helix (red) of the TMD. (C) Single subunit representation viewed from the vestibule facing outward of E82Q mutant with gradient coloring, ranging from white to blue (0–10 ΔCα/ΔCα<sub>ctrl</sub>), representing the greatest RMSD deviations from the Wt pH 4 structures. The Cα atom is represented as a sphere for orientation. ECD, extracellular domain; GLIC, *G. violaceus* ligand-gated ion channel; PDB ID, protein data bank identification code; RMSD, root-mean-squared deviation; TMD, transmembrane domain; Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g006
<figure><img src="image_10.png" /><figcaption>**Fig 7.** TMD residues decrease sensitivity. Stick representation of the evaluated Glu and His residues are labeled. Glu residues are color coded based upon effect, with residues in red producing a significant decrease in pH<sub>50</sub>, and in yellow, an insignificant or weak effect, whereas the dark purple and magenta for His residues are synonymous to the red of Glu residues. The main subunit is shown in black, with ribbon and cartoon representation with the α-helices labeled, whereas neighboring subunits are represented in white ribbon. The M2 α-helix of the neighboring subunit (labeled M2') is also highlighted in cyan and in cartoon representation. A table of the pH<sub>50</sub> and ΔpH<sub>50</sub> is shown below with mean values ± standard deviation for evaluated mutants with recorded values noted in S1 Table. Light red and green zones with dashed lines represent the respective ½ log decrease and increase from Wt. Values outside this region are color coded for ease of interpretation, with red and darker red as significant decreases and a decrease of greater significance, respectively. TMD, transmembrane domain; Wt, wild-type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.2004470.g007

The structure of H277Q was performed and was also found to have no apparent deviation from the apparently open pH 4 conformation, which further corroborates the several previously published works that show H277 does not appear to play a role in proton gating (S1 Fig).

**Discussion and conclusion**

All of the titratable Asp/Glu/His residues of the GLIC were mutated and evaluated for their impact on the receptor sensitivity to protons. This approach allows for the identification of key regions controlling gating.
Several mutations within Loop C markedly alter the pH<sub>50</sub>. The decrease in proton sensitivity following neutralization of the E177/D178/E181 cluster, which is completely negated and reversed with the neutralization of the basic R133, K183, and R179 residues, as well as the importance of R179 mutations alone, suggests a role in gating rather than in proton sensing. This idea is in line with the observation that replacement of the entire Loop C with other pLGIC Loop C sequences or a polyglycine segment is still compatible with a pH-gated channel [47]. Interestingly, this region was suggested to be the binding site for organic acids that act as negative allosteric modulators, including caffeic [40] and crotonic [48] acids. These data thus further document the key allosteric role of this region, which is homologous to the neurotransmitter binding site of eukaryotic receptors.

Mutants of E35 are unique in strongly increasing the proton sensitivity. Neutralization of the E35 side chain by mutation to A, Q, or M produce a similar increase in proton sensitivity, whereas the mutation of E35K exhibits the same behavior with half the increase. The E35H mutation, however, shows the inverse phenotype. This suggests that a charge at this level impairs activation, and a hydrophobic/polar side chain is preferred for stabilization of the active state. Interestingly, the structures of the Wt GLIC, E35Q, and E35A at pH 4 show that the side chain of residue 35 is located in a hydrophobic environment bordered by the side chains of P113, L114, F116, F156, T158, and P247. In the GLIC structure at pH 7, which is representative of the resting state, the carboxyl moiety of E35 is more solvent exposed toward the inner vestibule because the side chain of P247 is moved away following the key revolving motion of the M2–M3 loop (S2 Fig). Therefore, the resting state would be predicted to better accommodate a charged residue, which is consistent with the functional data. It is noteworthy that E35K, while charged, does not decrease the proton sensitivity. However, the length of the lysine side chain likely places the charged ammonium away from the E35 carboxylate locus. Altogether, the data support that E35 is a bona fide proton sensor, whose charged form stabilizes the resting conformation, whereupon protonation favors the active state. The functional proton gating of E35Q and E35A mutants additionally suggests that if E35 is a bona fide proton sensor, it can't be the only one in the GLIC.

In the lower part of the ECD, the principal (+) face also carries two acidic residues, E82 and D115. E82Q and D115N display significantly higher proton sensitivity than E82A and D115A, respectively. Assuming that Gln and Asn faithfully mimic protonated side chains of Glu and Asp, this would suggest that residues E82 and D115 are also proton sensors. In this case, their protonated side chain would elicit stabilizing interactions specifically in the active state. To investigate this idea, the structures of E82Q and E82A were solved at pH 4. The E82A structure is similar to that of the Wt pH 4 structure, whereas that of E82Q shows a marked local reorganization of the Q82 moiety. Assuming E82 to be in a protonated state at pH 4, Gln is not a good mimic of a protonated Glu residue because Q82 is engaged in a completely different set of interactions with neighboring residues. Therefore, the E82Q mutant phenotype is likely due to local side chain reorganizations of Y28 and F156 rather than being a proton sensor itself.

Finally, D32 mutations significantly decreased the pH<sub>50</sub>. However, this residue is shown in the pH 4 structure to interact through a salt bridge with the already salt-bridged pair of R192 and D122. This interaction interlocking Loop 2, Loop 7, and the pre-M1 loop is often conserved throughout the pLGIC family [23], with notable exception of the ELIC and the α-glutamate-gated chloride channel from *C. elegans*. Thus, it is unlikely that it undergoes side chain protonation at pH 4, and likely contributes in gating rather than in proton sensing [22].

Mutation of Asp/Glu residues at the complementary (–) face, in general, decreased the proton sensitivity of the GLIC, with additive effects in many cases, especially for D88 and D86. However, the pattern of phenotypes does not allow discriminating whether those residues are involved in proton sensing or gating, because mutations to Asn/Gln also produced decreased
proton sensitivity. The key mutants at E26, D86, and D88 were solved at pH 4, each showing a structure similar to that of the Wt (S1 Fig). E26 clearly is an important residue in GLIC function, but other techniques are required to elucidate its mechanistic role.

Mutations at both entrances of the ion channel, E243 (outer ring) and E222/H277 (inner rings), as well as at the middle of the TMD within each subunit four-helix bundle (H235), produced significant decreases in proton sensitivity, but again here, the pattern of phenotypes does not allow discriminating between proton sensing and gating. H235 was previously proposed to be the major proton-sensing residue for GLIC activation, but the mutation H235Q, which is no longer titratable, still produces a functional proton-sensitive receptor, questioning its role as the sole proton sensor. Clearly, H235 is a very important residue for functionality, although its protonated state would be hardly accommodated by the hydrophobic environment of the TMD in both known apparently open and closed forms.

Molecular dynamics simulations using a string method optimization made predictions of the change in protonated state and sensitivity of titratable residues throughout the gating mechanism of the GLIC [23]. The report identified E26 and E177 to have loss-of-function effects when mutated. E26 indeed results in a loss-of-function effect, whereas E177 seems to have no effect. Additionally, the prediction identified E35, E75, and E243 as gain-of-function locations if mutated, whereas E75 and E243, when mutated, have either no effect or the opposite effect. Therefore, these results suggest that molecular dynamics simulations can indeed help identify potentially interesting residues, as the most potent residues were identified, but that the predictions from the simulations still need further experimental testing as the actual functional results do not all agree with predictions of the molecular dynamics data.

## Global robustness of GLIC gating toward Asp/Glu/His mutations

Despite combined mutation attempts, channels displaying constitutive activity were not observed, nor was there a complete abatement of function. This shows that proton sensing is not mediated by a single Asp/Glu/His residue, but rather by several residues located on different parts of the protein. Individual residues in many places may contribute partially to proton-sensitive channel gating, unraveling a mode of proton-controlled activation quite different from that of classical agonists of pLGIC family receptors, which act at a well-characterized orthosteric site.

The results indicated that E35 is an important proton sensor, as well as the existence of a number of other proton-sensing candidates, notably E26 and the D86/D88 pair in the ECD, E243 at the top of the channel, and E222/H277 at the bottom of the channel. E35, E26, and E243 are all located nearby the ECD-TMD interface, suggesting that at least part, if not all, of the pH-elicited activation of the GLIC bypasses the orthosteric site. E222/H277 might also contribute to pH sensing, implying in this case the diffusion of protons from the extracellular compartment, possibly through the ion channel itself. The E222/H277 pair most likely confers the intracellular proton concentration sensitivity of the GLIC, which was reported from inside-out patch clamp experiments [49].

Complementary approaches that would be too cumbersome to perform on all the mutants evaluated in this study, such as reporting on the open probability, channel kinetics, and linking the efficiency of gating to the apparent affinities identified, are necessary to elucidate the specific effects of the key residues in the complex gating mechanism of the GLIC.

The recent studies employing chimeras between ELIC domains and GLIC domains validate the possibility of multiple proton-sensing sites. Interestingly, the pH₅₀ of Lily, the GLIC<sub>ECD</sub>-ELIC<sub>TMD</sub> mimicking Lily (GELIC), and E35Q/A mutants all closely hover around pH 6.5 [25,36]. Altering the TMD residues surrounding E35 could have the same effect as charge
neutralization. However, another GLIC<sub>ECD</sub>-ELIC<sub>TMD</sub> mutant (GE), which lacks the mutation Y119F and the C-terminal conversion of RGITL<sub>ELIC</sub> to LFFGF<sub>GLIC</sub>, has a reported pH<sub>50</sub> of 3.63 [35]. It is important to note that neither the study of Lily nor of GELIC tested below a pH of 5, and both found significantly reduced currents, whereas GE displayed robust currents at pH's below 5. These discrepancies point to the intricate interactions between the ECD and TMD, which influence gating and possibly proton sensitivity. These studies combined with the current results also clearly show that the principal component of proton activation lies within the ECD and not the TMD, as previously assumed. The ELIC<sub>ECD</sub>-GLIC<sub>TMD</sub> chimera is not susceptible to proton activation until further mutation by I9'A, which has previously been shown to destabilize the recovery time for the receptor to return to the resting state [26].

Residues other than Asp/Glu/His may also contribute to proton sensing, notably Arg/Lys side chains that can display pKa's in the 4-6 range when located in very hydrophobic environments [50], or aromatic residues through cation-pi interaction with hydronium ions [51].

## Cases studies of other pH-sensing ion channels

Among pH-gated ion channels, two were studied in molecular details. First, the bacterial potassium channel KcsA was found to be inhibited by a local network of ionic/H-bond interactions between two Glu, two Arg, and a single His residue. A disruption of this network upon protonation allows channel opening [52]. In this case, E to A mutation of the two key residues increased the sensitivity of the channel to protons. Such a "suppression of charge to activate the channel" mechanism on the GLIC is observed for the single E35 residue.

The pattern of phenotypes observed here is reminiscent of acid-sensing ion channels (ASICs). Indeed, ASICs have been extensively studied, but efforts to map the sites for proton binding have so far yielded inconclusive results because mutation of multiple individual Asp and Glu residues independently produces changes in proton sensitivity [53]. Simple mutation of basic Arg/His/Lys residues cannot mimic a protonated state either, therefore further limiting this approach to study ASICs.

## Proton sensitive positions and their involvement, or lack of, in gating

The combined mutagenesis data support that, for several pH-sensitive Asp/Glu positions, notably E26Q, D32N, E82Q, D86N-D88N, D122N, E222Q, and E243Q, the change of charge does not contribute to activation. In these cases, the side chain carboxylic acid may engage in active state stabilizing interactions with neighboring residues or water molecules. Using the site-directed mutagenesis approach may only conclusively evaluate the perturbation of residue titratability, as Asn and Gln residues may in fact be poor mimics of protonated Asp and Glu residues. Either this is the case or the aforementioned residues must maintain their deprotonated state at pH 4 and provide stability to the active state as charged moieties, contradicting pKa predictions. Additionally, it is expected that a good mimic would maintain the same conformation as the protonated version of an Asp/Glu residue, and that this should be witnessed in a crystallographic structure. Yet the E82Q structure differs from Wt GLIC at pH 4 in which E82 would be presumably protonated. As previously mentioned, replacement of His residues may also only evaluate the removal of their titratability, whereas a simple mutagenesis cannot mimic the His protonated state. The functional results indicate either a sensitivity to protons or a structural importance for both H235 and H277, neither of which are crucial for proton-elicited gating of the GLIC.

Overall, the data suggest a complex network of H-bonds and polar interactions, with important positions below the orthosteric site, in the pH-sensitive channel opening
mechanism of the GLIC. A specific importance of the ECD-TMD interface was identified with the position E35 acting as a key sensor next to the D32-R192-D122 triad involved in the signal transduction between domains [43].

## Materials and methods

Residue numbering follows that which was established by Bocquet et al. [54] for consistency. This numbering is shifted by one from other reports [35,49] that use the correct GLIC numbering.

## Two-electrode voltage clamp electrophysiology

X. laevis oocytes were obtained from the Centre de Ressources Biologiques-Rennes, France. Defolliculated oocytes were maintained at 4°C in a modified Earth's saline solution (88 mM NaCl, 1 mM KCl, 1 mM MgSO₄, 2.5 mM NaHCO₃, 5 mM HEPES/Na pH 7.3) with 0.7 mM CaCl₂. After intranucleus injection of approximately 30 nL cDNA (80 ng/µl specified clone cDNA with 20 ng/µl of GFP cDNA), using a compressed air microinjection system, the oocytes were transferred to 18°C. One to two days later they were evaluated for GFP expression, and subsequently maintained at 15°C. Recordings were made 1–5 d after injection using low-resistance (0.2–2 MΩ) electrodes filled with 3 M KCl, with a –40 mV holding potential. The standard solution superfusing the oocyte during recording at room temperature was 100 mM NaCl, 3 mM KCl, 1 mM MgCl₂, 1 mM CaCl₂, 10 mM MES at either pH 7.3, 7.5, or 8 using NaOH. In order to maintain the desired pH and maintain an equivalent Na⁺ concentration in all solutions, the stock solution was adjusted to the indicated pH using NaOH, and lower pH solutions were subsequently obtained using HCl and the addition of the stock solution.

Measurements were performed using pClamp 10.5 software, with a Digidata 1440A digitizer and GeneClamp 500 amplifier (Molecular Devices, LLC, Sunnyvale, CA), using an 8-valve (PS-8H) programmable gravity-driven pinch valve perfusion system (Bioscience Tools, San Diego, CA). pH-dependent responses were elicited by switching from pH 7.3–pH 8 to a series of pH values, with a minimal pH of 3.7, and a 0.5-log-unit increment from either pH 6.5, or pH 7.5 for elevated pH₅₀ mutants. Perfusion times varied from 30 s to 90 s, with equivalent recorded wash periods in the holding buffer. pH series were performed either in an ascending order directly followed by a descent, or a descending order directly followed by an ascent, in order to remove bias. Evaluation of currents was done using Clampfit 10.5 (Molecular Devices, LLC), with Imax (µA) reported as the peak amplitude of negative going current with the holding current subtracted. The average of the two recorded peak values for a given pH was plotted in GraphPad Prism 4 (GraphPad Software, Inc, La Jolla, CA) against the pH and fitted with a nonlinear sigmoidal dose-response fit to obtain 1 (n-unit) value of pH₅₀. The given number of injections and number of recorded oocytes per construct are listed in S1 Table, with generally 2–4 oocytes recorded per injection. Fits with an R² < 0.9, a Hill-slope < 0.6 or > 4, or an absolute Imax < 0.9 µA were excluded from inclusion into mean values, and therefore not counted in the n-unit either. The Imax cutoff was chosen due to endogenous current at pH 4 or lower that appears on occasion. To be sure that this current does not greatly influence the pH₅₀ calculation, a cutoff greater than 3-fold was chosen. The arbitrary cutoffs for Hill-slope and R² were chosen to remain consistent in the removal of data that could not be properly fitted as a result of any circumstance. In order to minimize the influence of the intrinsic variability between oocyte batches, which show some variation in the Wt response to pH changes, mutants were also characterized by a ΔpH₅₀. The ΔpH₅₀ value corresponds to the variation of pH₅₀ between each mutant expressing cell and the Wt cell(s) measured in the same
batch of oocytes, using the same solutions of pH. The pH₅₀ values and ΔpH₅₀ are reported as mean ± standard deviation.

## Protein production and purification

GLIC variants were purified as previously reported [54]. PET20b vectors carrying the GLIC constructs fused with an N-terminal maltose-binding protein (MBP) tag were used to transform *E. coli* C43 cells, cultured at 37°C in the 2YT medium containing 100 mg/ml ampicillin. At an optical density (OD) of 0.8, the cultures were cooled to 20°C and 0.1 mM IPTG was added for an overnight induction. All the purification steps were carried out at 4°C. Proteins were extracted from the cell membrane with a Tris-buffered saline solution (TBS, 300 mM NaCl, 20 mM Tris pH 7.6) containing 2% n-dodecyl-β-D-maltoside (DDM). Solubilized proteins were subsequently isolated by ultracentrifugation, loaded onto an amylose resin, and incubated for approximately 1 h. The resin was extensively washed using a TBS containing 0.1% DDM and subsequently with a TBS containing 0.02% DDM. Thrombin enzyme was added into the MBP-GLIC-bound resin and incubated overnight. The GLIC protein was eluted using a TBS containing 0.02% DDM and 20 mM maltose. A further purification step was carried out by size exclusion chromatography on a Superose 6 10/300 column (GE Healthcare, Little Chalfont, United Kingdom), which was equilibrated with a TBS containing 0.02% DDM. The fractions of the peak corresponding to the molecular weight of the GLIC pentamer were collected and concentrated to around 10 mg/ml for crystallization.

## Crystallography

The concentrated protein was mixed at 1:1 volume ratio with a mother liquor solution containing 12%-14% PEG 4000, 400 mM NaSCN, 15% glycerol, 3% DMSO, and 0.1 M NaAcetate pH 4. The crystallization procedure was performed at 20°C using the hanging drop method. Microseeding was performed after an initial crystallization setup. Crystals appeared overnight and grew to full size in 2 wk. The crystals were flash frozen using liquid nitrogen. The diffraction data sets were collected either on the beamlines Proxima-1 of the SOLEIL Synchrotron or the European Synchrotron Radiation Facility (ESRF) ID29 and ID23A. The data sets were processed with xdsme [55] and further processed by CCP4 programs [56]. The structures were solved by molecular replacement in PHASER [57] using the GLIC (PDB ID: 4HFI) as the initial model. Further refinement was carried out using BUSTER refinement [58]. The quality of the structural models was checked by Molprobity web server [59]. All structures have been deposited in the Research Collaboratory for Structural Bioinformatics protein data bank (https://www.rcsb.org), respective deposition IDs and statistics for all crystal structures are listed in S2 Table.

## Crystallographic analysis

RMSD and Cα distance calculations were performed by aligning a given structure (M) using the 2 subunit pair chains A+B, B+C, C+D, and D+E upon the chain pairs A+B, B+C, C+D, and D+E of a reference Wt structure at either pH 4 or 4.6 (Wt). The alignment of pair E+A was replaced with only an alignment of the E chain of M upon each chain of Wt, and each individual chain of M was aligned upon chain E of Wt as the Pymol structural alignment algorithm had difficulties doing an alignment with nonconsecutive chain pairs. The pairwise alignment method was chosen to include quaternary intersubunit interface interactions that a simple single chain alignment would ignore. The RMSD of each residue, for which alternate side chain conformations were removed, as well as the Cα atom distance, was calculated between the 25 pairwise alignments, and subsequently averaged to yield V<sub>RMSD</sub>(M-Wt) and V<sub>ΔCα</sub>(M-Wt). The
calculated intrinsic variation, V<sub>RMSD</sub>(Wt<sub>2</sub>-Wt<sub>1</sub>) and V<sub>ΔCα</sub>(Wt<sub>2</sub>-Wt<sub>1</sub>), between the two Wt structures at pH 4.6 (Wt<sub>1</sub>) and 4.0 (Wt<sub>2</sub>) in which Wt<sub>1</sub> and Wt<sub>2</sub> are 3EAM and 4HFI, respectively, was used to obtain the reported “normalized” value (RMSD/RMSD<sub>ctrl</sub> and ΔCα/ΔCα<sub>ctrl</sub>) in which RMSD/RMSD<sub>ctrl</sub> = [V<sub>RMSD</sub>(M-Wt<sub>1</sub>)+V<sub>RMSD</sub>(M-Wt<sub>2</sub>)]/2 and ΔCα/ΔCα<sub>ctrl</sub> = [V<sub>ΔCα</sub>(M-Wt<sub>1</sub>)+V<sub>ΔCα</sub>(M-Wt<sub>2</sub>)]/2 for a given mutant, M, in Fig 6 and S1 Fig.

**Supporting information**

**S1 Table. Electrophysiological values.**
(DOCX)

**S2 Table. Crystallographic statistics.**
(DOCX)

**S1 Fig. Normalized crystallographic ΔCα distance and RMSD values.** (A) Normalized Cα distance values (ΔCα/ΔCα<sub>ctrl</sub>) for the remaining resolved structures not shown in Fig 6. The reference values of PDB IDs: 3EAM and 4HFI are shown with each group of structures. (B) Normalized RMSD for each residue (excluding the mutated residue) between all resolved mutant structures and the Wt pH 4, 4.6 structures, with the reference value of PDB IDs: 3EAM and 4HFI shown with each group of structures. PDB ID, protein data bank identification code; RMSD, root-mean-squared deviation; Wt, wild-type.
(TIF)

**S2 Fig. Comparison of E35 structural conformation between pH 7 and pH 4 structures.** (A) Peripheral view facing towards the vestibule and (B) 90° rotation to show the side orientation of the ECD-TMD interface around E35. E35 is shown in stick representation for the GLIC Wt pH 7 (PDB ID: 4NPQ, pink) and pH 4 (PDB ID: 4HFI, black) structures. Residues within 6 Å are shown as lines and labeled next to their pH 7 representations. The large NCS variances of residue conformations for the pH 7 structure are transparently shown as lines and in pale colors. There was no NCS variance for the pH 4 structure. The minimal cavity calculated from the union of all NCS variances for pH 7 is represented as a semitransparent pink surface in order to allow for the visualization of K33 and P247 of the pH 4 structure, which lie behind this cavity. This cavity all but disappears in the pH 4 conformation, with the only remaining region being completely solvent exposed. This is shown in black, again in a semitransparent surface representation. ECD, extracellular domain; GLIC, G. violaceus ligand-gated ion channel; NCS, noncrystallographic symmetry; PDB ID protein data bank identification code; TMD, transmembrane domain; Wt, wild-type.
(TIF)

**Acknowledgments**

The authors would like to thank Marie Prevost for use of GLIC clones; Karima Medjebeur and Nathalie Barilone for molecular biology assistance; Frédéric Poitevin for aid with scripting and discussion about crystallographic analysis; and Marc Gielen and Solène Lefebvre for critical reading of the manuscript.

**Author Contributions**

**Conceptualization:** Marc Delarue, Pierre-Jean Corringer.

**Data curation:** Ákos Nemcz, Haidai Hu, Zaineb Fourati, Catherine Van Renterghem.
**Formal analysis:** Ákos Nemecz, Haidai Hu, Zaineb Fourati, Marc Delarue, Pierre-Jean Corringer.

**Funding acquisition:** Marc Delarue, Pierre-Jean Corringer.

**Investigation:** Ákos Nemecz, Haidai Hu, Zaineb Fourati, Catherine Van Renterghem.

**Methodology:** Ákos Nemecz, Marc Delarue, Pierre-Jean Corringer.

**Project administration:** Ákos Nemecz, Marc Delarue, Pierre-Jean Corringer.

**Resources:** Marc Delarue, Pierre-Jean Corringer.

**Supervision:** Marc Delarue, Pierre-Jean Corringer.

**Validation:** Ákos Nemecz, Haidai Hu, Marc Delarue, Pierre-Jean Corringer.

**Visualization:** Ákos Nemecz, Haidai Hu, Zaineb Fourati, Marc Delarue, Pierre-Jean Corringer.

**Writing – original draft:** Ákos Nemecz, Pierre-Jean Corringer.

**Writing – review & editing:** Ákos Nemecz, Haidai Hu, Catherine Van Renterghem, Marc Delarue, Pierre-Jean Corringer.

## References

1. Nemecz À, Prevost MS, Menny A, Corringer P-J (2016) Emerging Molecular Mechanisms of Signal Transduction in Pentameric Ligand-Gated Ion Channels. Neuron 90: 452–470. https://doi.org/10.1016/j.neuron.2016.03.032 PMID: 27151638

2. Hassaine G, Deluz C, Grasso L, Wyss R, Tol MB, et al. (2014) X-ray structure of the mouse serotonin 5-HT3 receptor. Nature 512: 276–281. https://doi.org/10.1038/nature13552 PMID: 25119048

3. Morales-Perez CL, Noviello CM, Hibbs RE (2016) X-ray structure of the human α4β2 nicotinic receptor. Nature 538: 411–415. https://doi.org/10.1038/nature19785 PMID: 27698419

4. Miller PS, Aricescu AR (2014) Crystal structure of a human GABAA receptor. Nature 512: 270–275. https://doi.org/10.1038/nature13293 PMID: 24909990

5. Du J, Lü W, Wu S, Cheng Y, Gouaux E (2015) Glycine receptor mechanism elucidated by electron cryo-microscopy. Nature 526: 224–229. https://doi.org/10.1038/nature14853 PMID: 26344198

6. Huang X, Chen H, Michelsen K, Schneider S, Shaffer PL (2015) Crystal structure of human glycine receptor-α3 bound to antagonist strychnine. Nature 526: 277–280. https://doi.org/10.1038/nature14972 PMID: 26416729

7. Prevost MS, Sauguet L, Nury H, Van Renterghem C, Huon C, et al. (2012) A locally closed conformation of a bacterial pentameric proton-gated ion channel. Nat Struct Mol Biol 19: 642–649. https://doi.org/10.1038/nsmb.2307 PMID: 22580559

8. Basak S, Schmandt N, Gicheru Y, Chakrapani S (2017) Crystal structure and dynamics of a lipid-induced potential desensitized-state of a pentameric ligand-gated channel. Elife 6. https://doi.org/10.7554/eLife.23886 PMID: 28262093

9. Sauguet L, Poitevin F, Murail S, Van Renterghem C, Moraga-Cid G, et al. (2013) Structural basis for ion permeation mechanism in pentameric ligand-gated ion channels. EMBO J 32: 728–741. https://doi.org/10.1038/emboj.2013.17 PMID: 23403925

10. Hilf RJ, Bertozzi C, Zimmermann I, Reiter A, Trauner D, et al. (2010) Structural basis of open channel block in a prokaryotic pentameric ligand-gated ion channel. Nat Struct Mol Biol 17: 1330–1336. https://doi.org/10.1038/nsmb.1933 PMID: 21037567

11. Fourati Z, Sauguet L, Delarue M (2015) Genuine open form of the pentameric ligand-gated ion channel GLIC. Acta Crystallogr D Biol Crystallogr 71: 454–460. https://doi.org/10.1107/S1399004714026698 PMID: 25760595

12. Fourati Z, Ruza RR, Laverty D, Drège E, Delarue-Cochin S, et al. (2016) Barbiturates Bind in the GLIC Ion Channel Pore and Cause Inhibition by Stabilizing a Closed State. J Biol Chem 292: 1550–1558. https://doi.org/10.1074/jbc.M116.766964 PMID: 27986812
13. Sauguet L, Howard RJ, Malherbe L, Lee US, Corringer P-J, et al. (2013) Structural basis for potentiation by alcohols and anaesthetics in a ligand-gated ion channel. Nat Commun 4: 1697. https://doi.org/10.1038/ncomms2682 PMID: 23591864

14. Nury H, Van Renterghem C, Weng Y, Tran A, Baaden M, et al. (2011) X-ray structures of general anaesthetics bound to a pentameric ligand-gated ion channel. Nature 469: 428–431. https://doi.org/10.1038/nature09647 PMID: 21248852

15. Velisetty P, Chakrapani S (2012) Desensitization mechanism in prokaryotic ligand-gated ion channel. J Biol Chem 287: 18467–18477. https://doi.org/10.1074/jbc.M112.348045 PMID: 22474322

16. Velisetty P, Chalamalasetti SV, Chakrapani S (2012) Conformational transitions underlying pore opening and desensitization in membrane-embedded Gloeobacter violaceus ligand-gated ion channel (GLIC). J Biol Chem 287: 36864–36872. https://doi.org/10.1074/jbc.M112.401067 PMID: 22977232

17. Dellisanti CD, Ghosh B, Hanson SM, Raspanti JM, Grant VA, et al. (2013) Site-directed spin labeling reveals pentameric ligand-gated ion channel gating motions. PLoS Biol 11: e1001714. https://doi.org/10.1371/journal.pbio.1001714 PMID: 24260024

18. Menny A, Lefebvre SN, Schmidpeter PA, Drège E, Fourati Z, et al. (2017) Identification of a pre-active conformation of a pentameric channel receptor. Elife 6. https://doi.org/10.7554/eLife.23955 PMID: 28294942

19. Murail S, Howard RJ, Broemstrup T, Bertaccini EJ, Harris RA, et al. (2012) Molecular mechanism for the dual alcohol modulation of Cys-loop receptors. PLoS Comput Biol 8: e1002710. https://doi.org/10.1371/journal.pcbi.1002710 PMID: 23055913

20. Calimet N, Simoes M, Changeux J-P, Karplus M, Taly A, et al. (2013) A gating mechanism of pentameric ligand-gated ion channels. Proc Natl Acad Sci USA 110: E3987–E3996. https://doi.org/10.1073/pnas.1313785110 PMID: 24043807

21. Arcario MJ, Mayne CG, Tajkhorshid E (2017) A membrane-embedded pathway delivers general anaesthetics to two interacting binding sites in the Gloeobacter violaceus ion channel. J Biol Chem 292: 9480–9492. https://doi.org/10.1074/jbc.M117.780197 PMID: 28420728

22. Sauguet L, Shahsavar A, Poitevin F, Huon C, Menny A, et al. (2014) Crystal structures of a pentameric ligand-gated ion channel provide a mechanism for activation. Proc Natl Acad Sci USA 111: 966–971. https://doi.org/10.1073/pnas.1314997111 PMID: 24367074

23. Lev B, Murail S, Poitevin F, Cromer BA, Baaden M, et al. (2017) String method solution of the gating pathways for a pentameric ligand-gated ion channel. Proc Natl Acad Sci USA 114: E4158–E4167. https://doi.org/10.1073/pnas.1617567114 PMID: 28487483

24. Moraga-Cid G, Sauguet L, Huon C, Malherbe L, Girard-Blanc C, et al. (2015) Allosteric and hyperekplexic mutant phenotypes investigated on an α1 glycine receptor transmembrane structure. Proc Natl Acad Sci USA 112: 2865–2870. https://doi.org/10.1073/pnas.1417864112 PMID: 25730860

25. Duret G, Van Renterghem C, Weng Y, Prevost M, Moraga-Cid G, et al. (2011) Functional prokaryotic-eukaryotic chimera from the pentameric ligand-gated ion channel family. Proc Natl Acad Sci USA 108: 12143–12148. https://doi.org/10.1073/pnas.1104494108 PMID: 21730130

26. Bocquet N, de Carvalho LP, Cartaud J, Neyton J, Le Poupon C, et al. (2007) A prokaryotic proton-gated ion channel from the nicotinic acetylcholine receptor family. Nature 445: 116–119. https://doi.org/10.1038/nature05371 PMID: 17167423

27. Abdrakhmanova G, Cleemann L, Lindstrom J, Morad M (2004) Differential modulation of beta2 and beta4 subunits of human neuronal nicotinic acetylcholine receptors by acidification. Mol Pharmacol 66: 347–355. https://doi.org/10.1124/mol.66.2.347 PMID: 15266026

28. Chen Z, Dillon GH, Huang R (2003) Molecular determinants of proton modulation of glycine receptors. J Biol Chem 279: 876–883. https://doi.org/10.1074/jbc.M307684200 PMID: 14563849

29. Chen Z, Huang R (2007) Identification of residues mediating inhibition of glycine receptors by protons. Neuropharmacology 52: 1606–1615. https://doi.org/10.1016/j.neuropharm.2007.03.005 PMID: 17459427

30. Wilkins ME, Hosie AM, Smart TG (2005) Proton modulation of recombinant GABAA receptors: influence of GABA concentration and the β subunit TM2-TM3 domain. J Physiol 567: 365–377. https://doi.org/10.1113/jphysiol.2005.088823 PMID: 15946973

31. Wilkins ME, Hosie AM, Smart TG (2002) Identification of a beta subunit TM2 residue mediating proton modulation of GABA type A receptors. J Neurosci 22: 5328–5333. PMID: 12097484.

32. Beg AA, Ernst GG, Nix P, Davis MW, Jorgensen EM (2008) Protons act as a transmitter for muscle contraction in C. elegans. Cell 132: 149–160. https://doi.org/10.1016/j.cell.2007.10.058 PMID: 18191228

33. Mounsey KE, Dent JA, Holt DC, McCarthy J, Currie BJ, et al. (2007) Molecular characterisation of a pH-gated chloride channel from Sarcoptes scabiei. Invert Neurosci 7: 149–156. https://doi.org/10.1007/s10158-007-0050-6 PMID: 17602250
34. Schnizler K, Saeger B, Pfeffer C, Gerbaulet A, Ebbinghaus-Kintscher U, et al. (2005) A novel chloride channel in Drosophila melanogaster is inhibited by protons. J Biol Chem 280: 16254-16262. https://doi.org/10.1074/jbc.M411759200 PMID: 15713676

35. Hénault CM, Baenziger JE (2016) Functional characterization of two prokaryotic pentameric ligand-gated ion channel chimeras—role of the GLIC transmembrane domain in proton sensing. Biochim Bio-phys Acta 1859: 218-227. https://doi.org/10.1016/j.bbamem.2016.11.006 PMID: 27845033

36. Alqazzaz MA, Price KL, Lummis SCR (2017) The Proton Responsiveness in the Extracellular Domain of GLIC Differs in the Presence of the ELIC Transmembrane Domain. Biochemistry 56: 2134-2138. https://doi.org/10.1021/acs.biochem.6b00900 PMID: 28383883

37. Schmandt N, Velisetty P, Chalamalasetti SV, Stein RA, Bonner R, et al. (2015) A chimeric prokaryotic pentameric ligand-gated channel reveals distinct pathways of activation. J Gen Physiol 146: 323-340. https://doi.org/10.1085/jgp.201511478 PMID: 26415570

38. Rienzo M, Lummis SCR, Dougherty DA (2014) Structural requirements in the transmembrane domain of GLIC revealed by incorporation of noncanonical histidine analogs. Chem Biol 21: 1700-1706. https://doi.org/10.1016/j.chembiol.2014.10.019 PMID: 25525989

39. Cheng X, Sine SM (2011) Intramembrane proton binding site linked to activation of bacterial pentameric ion channel. J Biol Chem 287: 6482-6489. https://doi.org/10.1074/jbc.M111.305839 PMID: 22084238

40. Prevost MS, Delarue-Cochin S, Marteaux J, Colas C, Van Renterghem C, et al. (2013) Identification of cinnamic acid derivatives as novel antagonists of the prokaryotic proton-gated ion channel GLIC. J Med Chem 56: 4619-4630. https://doi.org/10.1021/jm400374q PMID: 23682762

41. Purohit P, Chakraborty S, Auerbach A (2015) Function of the M1 π-helix in endplate receptor activation and desensitization. J Physiol 593: 2851-2866. https://doi.org/10.1113/JP270223 PMID: 25929452

42. Lee WY, Sine SM (2005) Principal pathway coupling agonist binding to channel gating in nicotinic receptors. Nature 438: 243-247. https://doi.org/10.1038/nature04156 PMID: 16281039

43. Bertozzi C, Zimmermann I, Engeler S, Hilf RJC, Dutzler R (2016) Signal Transduction at the Domain Interface of Prokaryotic Pentameric Ligand-Gated Ion Channels. PLoS Biol 14: e1002393. https://doi.org/10.1371/journal.pbio.1002393 PMID: 26943937

44. Therien JPD, Baenziger JE (2017) Pentameric ligand-gated ion channels exhibit distinct transmembrane domain archetypes for folding/expression and function. Sci Rep 7: 450. https://doi.org/10.1038/s41598-017-00573-2 PMID: 28348412

45. Fritsch S, Ivanov I, Wang H, Cheng X (2011) Ion selectivity mechanism in a bacterial pentameric ligand-gated ion channel. Biophys J 100: 390-398. https://doi.org/10.1016/j.bpj.2010.11.077 PMID: 21244835

46. Cymes GD, Grosman C (2016) Identifying the elusive link between amino acid sequence and charge selectivity in pentameric ligand-gated ion channels. Proc Natl Acad Sci USA 113: E7106. https://doi.org/10.1073/pnas.1608519113 PMID: 27791102

47. Gonzalez-Gutierrez G, Cuello LG, Nair SK, Grosman C (2013) Gating of the proton-gated ion channel from Gloeobacter violaceus at pH 4 as revealed by X-ray crystallography. Proc Natl Acad Sci USA 110: 18716-18721. https://doi.org/10.1073/pnas.1313156110 PMID: 24167270

48. Alqazzaz MA, Price KL, Lummis SCR (2016) Crotonic Acid Blocks the Gloeobacter Ligand-Gated Ion Channel (GLIC) via the Extracellular Domain. Biochemistry 55: 5947-5951. https://doi.org/10.1021/acs.biochem.6b00531 PMID: 27739668

49. Hilf RJC, Dutzler R (2009) Structure of a potentially open state of a proton-activated pentameric ligand-gated ion channel. Nature 457: 115-118. https://doi.org/10.1038/nature07461 PMID: 18987630

50. Cymes GD, Ni Y, Grosman C (2005) Probing ion-channel pores one proton at a time. Nature 438: 975-980. https://doi.org/10.1038/nature04293 PMID: 16355215

51. Wang S, Yan R, Zhang X, Chu Q, Shi Y (2014) Molecular mechanism of pH-dependent substrate transport by an arginine-agmatine antiporter. Proc Natl Acad Sci USA 111: 12734-12739. https://doi.org/10.1073/pnas.1414093111 PMID: 25136114

52. Thompson AN, Posson DJ, Parsa PV, Nimigean CM (2008) Molecular mechanism of pH sensing in KcsA potassium channels. Proc Natl Acad Sci USA 105: 6900-6905. https://doi.org/10.1073/pnas.0800873105 PMID: 18443286

53. Vullo S, Bonifacio G, Roy S, Johner N, Bernèche S, et al. (2017) Conformational dynamics and role of the acidic pocket in ASIC pH-dependent gating. Proc Natl Acad Sci USA 114: 3768-3773. https://doi.org/10.1073/pnas.1620560114 PMID: 28320963

54. Bocquet N, Nury H, Baaden M, Le Poupon C, Changeux J-P, et al. (2008) X-ray structure of a pentameric ligand-gated ion channel in an apparently open conformation. Nature 457: 111-114. https://doi.org/10.1038/nature07462 PMID: 18987633
55. Kabsch W (2010) Integration, scaling, space-group assignment and post-refinement. *Acta Crystallogr D Biol Crystallogr* 66: 133–144. https://doi.org/10.1107/S0907444909047374 PMID: 20124693

56. The CCP4 suite: programs for protein crystallography. (1994) The CCP4 suite: programs for protein crystallography. *Acta Crystallogr D Biol Crystallogr* 50: 760–763. https://doi.org/10.1107/S0907444994003112 PMID: 15299374

57. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, et al. (2007) Phaser crystallographic software. *J Appl Crystallogr* 40: 658–674. https://doi.org/10.1107/S0021889807021206 PMID: 19461840

58. Smart OS, Womack TO, Flensburg C, Keller P, Paciorek W, et al. (2012) Exploiting structure similarity in refinement: automated NCS and target-structure restraints in BUSTER. *Acta Crystallogr D Biol Crystallogr* 68: 368–380. https://doi.org/10.1107/S0907444911056058 PMID: 22505257

59. Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, et al. (2009) MolProbity: all-atom structure validation for macromolecular crystallography. *Acta Crystallogr D Biol Crystallogr* 66: 12–21. https://doi.org/10.1107/S0907444909042073 PMID: 20057044