# Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

Sophie Shi, Solène N Lefebvre, Laurie Peverini, Adrien Cerdan, Paula Milán Rodríguez, Marc Gielen, Jean-Pierre Changeux, Marco Cecchini, Pierre-Jean Corringer

► To cite this version:

Sophie Shi, Solène N Lefebvre, Laurie Peverini, Adrien Cerdan, Paula Milán Rodríguez, et al. Illumination of a progressive allosteric mechanism mediating the glycine receptor activation. Nature Communications, 2023, 14 (1), pp.795. 10.1038/s41467-023-36471-7. pasteur-03993599v2

HAL Id: pasteur-03993599

https://pasteur.hal.science/pasteur-03993599v2

Submitted on 17 Feb 2023

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
<figure><img src="image_2.png" /></figure>

# Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

Received: 24 May 2022

Accepted: 2 February 2023

Published online: 13 February 2023

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

Check for updates

Sophie Shi<fnref n="1" />,²,⁶, Solène N. Lefebvre ©<fnref n="1" />,²,⁶, Laurie Peverini ©<fnref n="1" />, Adrien H. Cerdan ©³, Paula Milán Rodríguez ©<fnref n="1" />,³, Marc Gielen ©<fnref n="1" />,⁴, Jean-Pierre Changeux⁵, Marco Cecchini ©³✉ & Pierre-Jean Corringer<fnref n="1" />✉

Pentameric ligand-gated ion channel mediate signal transduction at chemical synapses by transiting between resting and open states upon neurotransmitter binding. Here, we investigate the gating mechanism of the glycine receptor fluorescently labeled at the extracellular-transmembrane interface by voltage-clamp fluorometry (VCF). Fluorescence reports a glycine-elicited conformational change that precedes pore opening. Low concentrations of glycine, partial agonists or specific mixtures of glycine and strychnine trigger the full fluorescence signal while weakly activating the channel. Molecular dynamic simulations of a partial agonist bound-closed Cryo-EM structure show a highly dynamic nature: a marked structural flexibility at both the extracellular-transmembrane interface and the orthosteric site, generating docking properties that recapitulate VCF data. This work illuminates a progressive propagating transition towards channel opening, highlighting structural plasticity within the mechanism of action of allosteric effectors.

Pentameric ligand-gated ion channels (pLGICs), including nAChRs, 5-HT₃Rs, GlyRs, and GABAARs constitute a superfamily of trans-membrane receptors mediating intercellular communications in the nervous system<fnref n="1" />,². They transduce the binding of a neurotransmitter at their orthosteric site into the opening of an intrinsic ion channel, leading to ion fluxes that promote cell excitation or inhibition. The orthosteric binding site also binds partial agonists, with less efficacy than full agonists to elicit channel opening, and competitive antagonists. These distinct pharmacological profiles were interpreted in terms of the concerted two-state model (Monod et al., 1965) based upon a pre-existing equilibrium between the resting and active states that is differentially shifted depending on the nature of the effector³⁻⁵.

In between these two equilibrium states, the conformational pathway that the protein follows during activation remains elusive. Single-channel recordings of numerous mutants of the muscle

nAChR analyzed by REFERs (rate equilibrium linear free energy relationships) suggested an early motion of the extracellular domain (ECD) and a late motion of the transmembrane domain (TMD) during activation⁶. Likewise, single-channel kinetics analysis on GlyRs and nAChRs were interpreted via a model involving a multistep process since the accurate fitting of the experimental data requires the introduction of intermediate states “flip”⁷ or “primed”⁸ in the gating transition. Although capturing structural reorganizations only indirectly, the electrophysiological data suggest that the transition is progressive, with the conformational changes starting from the ECD where the orthosteric site is located and propagating to the TMD to reach the channel gate. Moreover, computational studies of various pLGICs by Molecular Dynamics highlight their passage through a complex conformational landscape with the gating transition being composed of a progressive reorganization of the multimeric receptor architecture⁹⁻¹³.

<footnote n="1">¹Institut Pasteur, Université Paris Cité, CNRS UMR 3571, Channel-Receptors Unit, Paris, France. ²Sorbonne Université, Collège doctoral, Paris, France. ³Institut de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg, Cedex, France. ⁴Sorbonne Université, 21 rue de l'École de Médecine, 75006 Paris, France. ⁵Neuroscience Department, Institut Pasteur, Collège de France, Paris, France. ⁶These authors contributed equally: Sophie Shi, Solène N. Lefebvre. ✉ e-mail: mcecchini@unistra.fr; pierre-jean.corringer@pasteur.fr</footnote>
At the structural level, fluorescence experiments provided evidence that the bacterial pLGIC homolog GLIC undergoes a progressive reorganization, with a first pre-activation step involving quaternary compaction of the ECD¹⁴ followed by pore opening associated with a structural reorganization of the TMD and ECD-TMD interface¹⁵. However, direct evidence for progressive signal transduction in eukaryotic receptors is essentially lacking, despite recent cryo-EM structures of pLGICs revealing that the binding of orthosteric ligands is often accompanied by a significant reorganization of the extracellular domain (ECD), while the channel remains closed at the transmembrane domain (TMD), e.g. in 5-HT₃R¹⁶,¹⁷, GABA<sub>A</sub>R¹⁸, GlyR¹⁹, and nAChR²⁰. Thus, pLGIC structures appear to be more flexible than previously anticipated, suggesting that ligands or classes of molecules could stabilize specific conformation(s) in the conformational landscape accessible to the protein. Whether the captured conformations solved by cryo-EM correspond or not to intermediates during transitions on the pathway to channel activation remains an open question. Addressing these questions requires monitoring the conformational dynamics of a fully functional receptor embedded in a plasma membrane and probing its relationship with both activation and desensitization. The voltage-clamp fluorometry (VCF) technique is perfectly suited to this purpose, since it allows the simultaneous recording of local conformational changes by fluorescence together with transmembrane current fluxes by electrophysiology.

Here, we used the α1 homeric GlyR, which is surely among the best characterized pLGIC by electrophysiology²¹,²², cryo-EM alone or in complex with ligands¹⁶,¹⁹,²³, VCF²⁴⁻³¹, and MD simulations³²,³³. We report the development and characterization of a mutant GlyRα1 bearing a fluorescent sensor at the interface between ECD and TMD. Its pharmacological characterization using different allosteric effectors by VCF coupled to MD simulations and docking illuminates the occurrence of intermediates, displaying dynamic properties of the ECD, on the path to channel opening consistent with a progressive mechanism of gating.

<figure><img src="image_4.png" /><figcaption>**Fig. 1 |** Localization of K143W/Q219C sensor in the α1 glycine receptor structure. A Side view of the zebrafish α1 glycine receptor structure in Apo state (PDB:6PXD) showing the position of the mutated residues K143W in Cys-loop and Q219C in Pre-M1 loop. Right panel: representation of the two isomers of 5(6)-carboxytetramethylrhodamine methanethiosulfonate (MTS-TAMRA) used for the labeling of the mutated cysteine. B Structural comparison of the quenching sensor K143W/Q219C between the Apo state in gray (PDB:6PXD) and Gly-bound open state in green (PDB:6PM6), residues are presented in cyan for apo state and magenta for Gly-bound open state. The distance calculated between the Cα of the two residues varied from 12.9 Å for the Apo state to 15.6 Å for the Open state.</figcaption></figure>

## Results

### Design of a fluorescent-quenching sensor reporting conformational changes at the ECD-TMD interface

To generate a fluorescent sensor related to the activation process, we compared the structure of the zebrafish GlyRα1 in the resting-like apo state (PDB:6PXD) with the glycine-bound open-channel structure (PDB:6PM6)¹⁹. Besides the well-documented contraction of the orthosteric site at the subunit-subunit interfaces in the ECD and the opening of the channel in the TMD, visual inspection also shows an expansion of the ECD-TMD interface, where the Cys-loop of each subunit moves away by nearly 3 Å from the pre-M1 region of the adjacent subunit (Fig. 1A, B). In order to monitor this reorganization, we used the tryptophan quenching technique that can sense the relative place of two residues³⁴,³⁵. This technique is particularly relevant to monitor distance changes of 5–15 Å range¹⁴,¹⁵. We introduced a cysteine at position Q219 in the pre-M1 region previously used for labeling with fluorescent dyes²⁸, together with a tryptophan in Cys-loop at position K143, in order to quench the fluorophore in a conformation-dependent manner (Fig. 1A, B).

### The quenching sensor K143W/Q219C reveals glycine-elicited fluorescence variations at much lower glycine concentrations than currents

The endogenous extracellular single cysteine of the human GlyRα1 was mutated to avoid non-specific labeling, yielding WT-C41V displaying wild-type-like properties in electrophysiology (Supplementary Fig. 1). All mutations shown here present this background C41V mutation. Upon introduction of the quenching sensor mutations, the unlabeled K143W/Q219C was analyzed by two-electrode voltage-clamp electrophysiology (TEVC). In all dose-response curves presented here, we used brief perfusion (10–30 s) of the glycine solutions to reach a current plateau with no apparent desensitization. The unlabeled K143W/Q219C displays an EC₅₀<sup>current</sup> that is 16-fold higher than that of WT-C41V (Supplementary Fig. 3A, B and Table 1), indicating a loss-of-function phenotype. We controlled that the mutations do not impair ion channel permeation by performing outside-out single-channel recordings (Fig. 2C) and measuring a unitary conductance of 86.7 ± 4.6 pS that is identical to that of WT-C41V (85.7 ± 4.5 pS) (Supplementary Fig. 1).

After labeling with MTS-TAMRA (Fig. 1A), K143W/Q219C was analyzed by VCF, using a dedicated recording chamber (Supplementary Fig. 2). Its EC₅₀<sup>current</sup> decreases by 1.6-fold, indicating that the labeling produces a slight gain-of-function that partially counter-balances the effect of the K143W/Q219C mutations. Robust glycine-elicited fluorescence variations were also observed (Fig. 2A, B). They start at non-activating glycine concentrations (5–25 µM) and reach a plateau at the beginning of the current dose-response curve, with a maximal fluorescence amplitude ΔF/F corresponding to a 12.9 ± 1.0% decrease. Therefore, the EC₅₀<sup>fluo</sup> is 73-fold lower than the EC₅₀<sup>current</sup>. As a control, we also recorded the labeled Q219C/K143F and Q219C, that show weak ΔF/F (no variation and a 2 ± 0.2% decrease, respectively) (Supplementary Fig. 3B, C), indicating that the ΔF/F of the K143W/Q219C sensor is dominated by the quenching provided by the introduced tryptophan. Of note, the MTS-TAMRA used here corresponds to a commercial mixture of isomers (Fig. 1A), each of which was verified to produce the same VCF pattern separately (Supplementary Fig. 3D).

### The quenching sensor K143W/Q219C reports an "intermediate" conformational reorganization occurring before the onset of currents

Inspection of the current/fluorescence time courses and quantification by single exponential fitting reveals important features. First, the onset of both the fluorescence and current traces is faster with increased glycine concentrations whereas their offset during agonist rinsing is slower (Fig. 2D, F, Supplementary Table 3, and Supplementary Fig. 4).
<table><caption>Table 1 | EC₅₀ values for current and fluorescence responses to glycine at labeled and unlabeled C41V-GlyRα1 mutants</caption><thead><tr><th></th><th>EC<sub>50</sub><sup>current</sup> (µM)</th><th>n<sub>H</sub></th><th>EC<sub>50</sub><sup>fluo</sup> (µM)</th><th>n<sub>H</sub></th><th>n</th></tr></thead><tbody><tr><td>WT without C41V MTS-TAMRA</td><td>165.84 ± 14.06</td><td>1.55 ± 0.08</td><td></td><td></td><td>5</td></tr><tr><td>WT without C41V unlabeled</td><td>79.30 ± 11.34</td><td>2.16 ± 0.36</td><td></td><td></td><td>8</td></tr><tr><td>WT MTS-TAMRA</td><td>162.87 ± 43.21</td><td>1.18 ± 0.13</td><td></td><td></td><td>5</td></tr><tr><td>WT unlabeled</td><td>152.40 ± 28.87</td><td>1.33 ± 0.26</td><td></td><td></td><td>4</td></tr><tr><td>Q219C MTS-TAMRA</td><td>56.80 ± 6.65</td><td>1.84 ± 0.13</td><td>550.95 ± 67.55</td><td>1.71 ± 0.15</td><td>6</td></tr><tr><td>Q219C unlabeled</td><td>70.52 ± 3.47</td><td>2.23 ± 0.46</td><td></td><td></td><td>5</td></tr><tr><td>K143W/Q219C MTS-TAMRA</td><td>1541.12 ± 338.71</td><td>1.06 ± 0.11</td><td>21.00 ± 2.26</td><td>0.96 ± 0.10</td><td>5</td></tr><tr><td>K143W/Q219C unlabeled</td><td>2485.20 ± 145.08</td><td>1.04 ± 0.03</td><td></td><td></td><td>5</td></tr><tr><td>K143W/Q219C/N46D/N61D MTS-TAMRA</td><td>23029.40 ± 4281.72</td><td>1.41 ± 0.04</td><td>1428.52 ± 396.81</td><td>1.12 ± 0.06</td><td>5</td></tr><tr><td>K143W/Q219C V280M MTS-TAMRA</td><td>22.79 ± 11.24</td><td>0.544 ± 0.06</td><td></td><td></td><td>5</td></tr><tr><td>K143W/Q219C R271Q MTS-TAMRA</td><td></td><td></td><td>118.09 ± 15.11</td><td>0.81 ± 0.06</td><td>6</td></tr><tr><td>R271C MTSR</td><td>2863.40 ± 766.67</td><td>1.36 ± 0.11</td><td>4517.60 ± 1335.87</td><td>1.12 ± 0.11</td><td>5</td></tr></tbody></table>

Second, at each glycine concentration, the rise time of the fluorescence change is systematically faster than that of the currents whereas the offset tends to be slower, particularly at glycine concentrations >200 µM. It is noteworthy that these kinetic values are limited by the rate of agonist perfusion in the recording chamber, and do not reflect the intrinsic kinetics of the receptor. For instance, upon perfusion of 200 µM Gly (Fig. 2D), the glycine concentration at the oocyte surface will progressively rise, first reaching the concentration eliciting full fluorescence variation (100 µM), then reaching after a delay in the concentration eliciting currents (200 µM). Conversely, upon glycine washing, its concentration will rapidly drop below 100 µM, causing deactivation, but will take more time to drop below the concentrations causing fluorescence variations. The relatively slow kinetics of perfusion of the agonist thus separate in time the molecular events causing the fluorescence variation and the current, the former preceding the latter.

Figure 2A, B shows that 100 µM glycine, at steady state, elicits most of the fluorescence variation but no significant current. This indicates that the receptor has completed the full movement reported by the quenching sensor with a closed channel. Hence, VCF data provide direct evidence for a structural conformational transition toward one or a cluster of conformations that are structurally distinct from both resting and active. Importantly, kinetic data show that these conformations appear before channel opening in our recording conditions (Fig. 2D, F). This suggests that they correspond to early “intermediate conformation(s)” along the allosteric transition pathway to activation that are not related to desensitization. To further investigate this possibility, we applied glycine at a high concentration (30 mM) for a prolonged period (1 min), leading to a rapid onset of the current followed by a slower decrease (34.6 ± 4.9% decrease after 1 min) caused by desensitization. The fluorescence strongly decreases during the activation phase, then reaches a plateau and remains stable during the second phase, providing evidence that the quenching sensor does not report movements related to desensitization (Fig. 2E).

## The Gly-elicited fluorescent and current changes are linked to glycine binding to the orthosteric site

The marked separation of the current (ΔI) and fluorescence (ΔF) dose-response curves raises the possibility that both processes could be mediated by different classes of glycine binding sites within and outside the orthosteric site. To challenge this possibility, we introduced two mutations into the orthosteric binding site (N61D/N46D, Fig. 3A) that strongly decrease the affinity for glycine³⁶. The K143W/Q219C/N61D/N46D shows a parallel rightward shift of both ΔI and ΔF curves (15-fold and 68-fold increase in EC₅₀<sup>current</sup> and EC₅₀<sup>fluo</sup>, respectively) (Fig. 3B, C and Table 1), indicating that both processes are strongly linked to glycine binding to the orthosteric site. We also found that ΔF

changes are sensitive to known allosteric hyperekplexia mutations at the ECD-TMD interface. Combined with K143W/Q219C, the loss of function R271Q³⁷, located at the top of the M2 helix, fails to generate currents (certainly due to the combination with K143W) but a ΔF is observed and its curve is shifted to higher glycine concentrations (Supplementary Fig. 5A). Conversely, combined with K143W/Q219C, the gain of function V280M³⁸, located within the loop linking the M2 and M3 helices (M2-M3 loop) does not show fluorescence variation, but currents with a ΔI curve shifted to lower concentrations. In addition, it shows high leak currents in the absence of glycine that are inhibited by strychnine, indicating constitutive channel openings (Supplementary Fig. 5B, C).

## Partial agonists, strychnine, and propofol differentially affect the current and fluorescence variations

First, we investigated the effect of agonists that bind at the orthosteric site but are less effective than glycine at activating the receptor, i.e., partial agonists. We selected β-alanine and taurine. They trigger respectively 54–98% and 20–42% of the glycine-elicited maximum current³⁹⁻⁴¹ in TEVC³⁹,⁴¹⁻⁴³. On K143W/Q219C, we found that their efficacy decreases, eliciting only 6.1 ± 2.2% and 2.0 ± 1.4% of the glycine-elicited maximum current respectively (Fig. 4A, B). Such a decrease in efficacy is commonly observed for loss-of-function mutants³¹. Surprisingly, both partial agonists elicit the full fluorescence variation recorded with glycine. β-Alanine and taurine show a large difference in the ΔF and ΔI dose-response curve, their EC₅₀<sup>fluo</sup> being respectively 66- and 70-fold lower than their EC₅₀<sup>current</sup> (Table 1). β-Alanine and taurine thus trigger the full transition mediating the fluorescence change, but they are weakly effective to trigger the “downward” process of channel opening.

Next, we investigated the competitive antagonist strychnine. When strychnine is applied alone at 5 µM, no current is observed as expected, but the fluorescence increases over the baseline in a direction opposite to the glycine-elicited quenching (dequenching). It then reaches a plateau corresponding to 38 ± 6% in the absolute amplitude of the maximal glycine-elicited fluorescence variation. We then applied strychnine (5 µM) during the perfusion of glycine at different concentrations (Fig. 5). At 100 µM glycine, strychnine inhibits both the current and fluorescence variations elicited by glycine. As before, strychnine tends to increase the fluorescence over the baseline, but the effect is not significant. At 1 mM of glycine, strychnine still totally inhibits the current, but surprisingly the glycine-elicited fluorescence quenching is decreased by only 48 ± 4%. Finally, at 10 mM of glycine, strychnine inhibits only 77 ± 2% of the current, while the fluorescence is almost not affected. Therefore, the co-application of strychnine with various concentrations of glycine differentially affects the fluorescence variations and currents.
<figure><img src="image_5.png" /><figcaption>**Fig. 2 | Electrophysiological and fluorescence characterization of K143W/Q219C sensor by VCF on C41V GlyR.** A Representative VCF recordings on oocytes of the mutant labeled with MTS-TAMRA (black for current and cyan for fluorescence). Glycine application triggers a current variation and a fluorescence quenching phenomenon with a maximum fluorescence variation that reaches 12.1 ± 1.1% of ΔF/F. At low concentrations of glycine (under 25 µM), only fluorescence quenching is observed without any current. B ΔI (black) and ΔF (cyan) dose-response curves with mean ± SEM values (n = 5). C Left panel: representative trace of single-channel recordings obtained in outside-out configuration recorded at -100 mV with concentrations of glycine at 100 µM. Right panel: histograms of current amplitude representing the closed state (c) and the open state (o). D Superimposition of the current and fluorescence recordings evoked by 200 µM</figcaption></figure>

glycine application shows that the onset of the fluorescence is faster than the onset of the current and that the fluorescence offset is slower than the current offset. E Representative recording with a high glycine concentration (30 mM) application shows that the desensitization triggered by a prolonged glycine application does not impact the fluorescence variation. A mean decrease of 30.63 ± 4.90% of the current elicited by 30 mM of glycine is observed (n = 4). F Left panel: examples of single exponential fitting (red line) of the current (left) and fluorescence (right) traces onset. Right panel: time constants τ (onset) values obtained via single exponential fitting with mean and SEM error bars at different glycine concentrations (n = 3 and 4 for glycine concentrations under 25 µM and n = 5 for other concentrations). Unpaired two-sided student t-test indicates the significance of the difference between fluorescence and current onset (\*\*P < 0.005; \*\*\*P < 0.0005).

Last, we investigated the general anesthetic propofol, a positive allosteric modulator of GlyRα1 that binds in the TMD of pLGICs⁴⁴⁻⁴⁷. At 200 and 300 µM, propofol alone elicits a fluorescence variation of about 3–5% ΔF/F without activating the receptor, while at higher

concentrations it also activates the currents in the absence of glycine (Supplementary Fig. 6). At a lower 100 µM concentration, propofol alone does not elicit significant ΔF. Interestingly, the co-application of propofol with a range of glycine concentrations shows that it
<figure><img src="image_6.png" /><figcaption>**Fig. 3 |** Effect of N46D/N61D mutations on ΔI and ΔF curves of K143W/Q219C on C41V GlyR. A Side view of zebrafish α1 glycine receptor structure in Apo state (PDB:6PXD) with details of loops and β-sheets forming the orthosteric site and the localization of N46D/N61D mutations. The glycine molecule is represented in green and N46D/N61D mutated residues (spheres) in black. B Representative VCF</figcaption></figure>

recordings of the mutant K143W/Q219C + N46D/N61D (current in black and fluorescence in cyan). C ΔI (cyan) and ΔF (black) dose-response curves with mean ± SEM show a parallel shift of the current and fluorescence curves of the K143W/Q219C/N46D/N64D mutant (solid line; *n* = 5) compared to K143W/Q219C alone (dotted line).

potentiates by 21 ± 7% the maximal currents and shifts by 8-fold the ΔI curve to lower concentrations. In contrast, propofol does not affect significantly the ΔF curve (Fig. 4C and Supplementary Table 1). Propofol thus does not alter the transition mediating the fluorescence change, but specifically facilitates the “downward” process of channel opening.

In conclusion, various pharmacological conditions yield a common phenotype with a change in ΔF with comparatively weaker or no change in ΔI: they are low glycine concentrations, partial agonists, specific mixtures of glycine and strychnine and 200–400 µM of propofol alone. This phenotype suggests that a substantial fraction of the receptors has underwent the local conformational motion causing quenching of K143W/Q219C in the lower part of the ECD but with a closed channel. However, it is possible that the glycine-strychnine combination and propofol alone change the environment of the sensor in a similar way as partial agonists and agonists do, giving similar ΔF signature, while the rest of the ECD could present different conformations. This is especially possible for propofol alone as it binds to the TMD of GlyRs and might thus present an alternative mechanism.

## The R271C mutant also shows a relatively increased fluorescence variation when activated by taurine as compared to glycine

Among pLGICs, the GlyR has been extensively studied by VCF by the team of Lynch et al., with the labeling of the upper M2 helix at position R271$^{25,31}$, the orthosteric site, the loop 2 and pre-M1 regions$^{28,30}$, and the TMD$^{26,27}$. In all cases, the fluorescence dose-response curves were superimposed or right-shifted as compared to the current dose-response curve under glycine application (Supplementary Fig. 7). Thus, the leftward shift of the ΔF dose-response curve shown here seems unique.

However, another landmark of K143W/Q219C concerns the action of partial agonists, that produce much higher maximal ΔF/ΔI ratios

than glycine. Interestingly, R271C$^{31}$ was also reported to display such a phenotype with taurine. On our VCF setup, with glycine, R271C produces an 18-fold shift of the ΔI curve to higher concentrations as compared to WT, and a near superimposition of the ΔI and ΔF curves. Taurine at saturation elicits significantly higher ΔF than ΔI, respectively 21 ± 3.4% and 5.4 ± 2.3% of glycine values. In addition, for taurine, the ΔF curve is shifted by 2-fold to lower concentrations as compared to the ΔI curve (Fig. 6A, B). These data, that are in qualitative agreement with those of Pless et al. (2007), indicate that a fraction of the taurine-elicited states of R271C presents a fluorescence variation with a closed channel.

## Molecular dynamics simulations of a cryo-EM partial-agonist bound state recapitulate the pharmacological profile of the intermediate state(s)

A recent cryo-EM analysis of the zebrafish GlyR$^{19}$ has revealed that the partial agonists taurine and GABA may stabilize a previously unseen closed-channel state (referred to as “tau-closed”) characterized by a significant reorganization of the ECD compared to the apo state. We hypothesized that the intermediate conformation(s) illuminated by VCF is structurally similar to tau-closed. To explore this hypothesis, we carried out all-atom MD simulations of the zebrafish apo-closed (PDB:6PXD), tau-closed (PDB:6PM3), and tau-open (PDB:6PM2) structures embedded in a POPC bilayer. In addition, since strychnine nicely fits the orthosteric site of apo-closed, we performed MD simulations of the strychnine-closed (stry-closed) complex as well. For each system, five independent simulations of 100 ns were carried out.

The structural stability of the protein was analyzed by monitoring the root-mean-square fluctuations (RMSF) from the average structure extracted from MD (Fig. 7A). In addition, we calculated the evolution of the twisting and blooming angles of the ECD during each simulation (Supplementary Fig. 8 and Supplementary Table 4). The simulations
<figure><img src="image_7.png" /><figcaption>**Fig. 4** | Electrophysiological and fluorescence characterization of partial-agonist (β-alanine and taurine) and propofol effects on K143W/Q219C sensor on C41V GlyR. A Right panel: Representative VCF recordings on an oocyte challenged with both β-alanine and glycine. Left panel: representative VCF recording on oocytes under the β-alanine application. Right panels: the ΔI (cyan) and ΔF (black) curves (normalized to the glycine maximal response recorded on each individual oocyte) with mean and SEM for K143W/Q219C under glycine (solid line) and β-alanine (dotted line) application. The ΔF curves for both molecules are superimposed but the efficacy of β-alanine to activate the receptor is lower than that of glycine (6.1 ± 2.2% of glycine maximum response in current, *n* = 6). B Right panel: Representative VCF recordings on an oocyte challenged with both taurine and</figcaption></figure>

glycine. Left panel: representative VCF recording in oocytes under the taurine application. Right panel: ΔI (cyan) and ΔF (black) curves (normalized to the glycine maximal response recorded on each individual oocyte) with mean and SEM for K143W/Q219C under glycine (solid line) and taurine (dotted line) application. Taurine elicits only 2.0 ± 1.4% of glycine maximum current, *n* = 6. C Left panels: Representative VCF recording on an oocyte challenged with maximal glycine concentration in the presence and absence of propofol at 500 µM. Middle panel: representative VCF recording on an oocyte under glycine and propofol co-application. Right panel: ΔI (cyan) and ΔF (black) dose-response curves with mean and SEM for K143W/Q219C under glycine and propofol (solid line, *n* = 6) and glycine alone (dotted line) application.

show that: i. apo-closed is the most flexible conformational state of the receptor (Fig. 7A) and strychnine binding stabilizes the bloomed conformation (expended ECD as seen in the Apo-closed cryo-EM structure), rigidifying the ECD; ii. the ECD of tau-closed is remarkably more flexible than that in tau-open despite the apparent structural similarity in the cryo-EM coordinates that show an un-bloomed (compaction of the ECDs) structure; iii. tau-closed has a weaker

affinity for taurine than tau-open as evidenced by four spontaneous unbinding events in the tau-closed simulations, that are associated with significant blooming of the structure (Supplementary Fig. 10); iv. the C loop that is critical for orthosteric-ligand binding is more flexible in tau-closed (RMSF of 3.5 Å) than in tau-open (RMSF of 2.5 Å) or stry-closed (RMSF of 2.0 Å); v. the conformational dynamics of the ECD-TMD interface (i.e., β1-β2 loop, Cys-loop, and the β8-β9 loop) is
<figure><img src="image_8.png" /><figcaption>**Fig. 5 | Differential inhibition of strychnine on fluorescence and current on K143W/Q219C sensor on C41V GlyR.** Upper left panel: Representative VCF recording where the application of strychnine alone at 5 µM triggers a fluorescence dequenching. Upper right panels: glycine is applied at different concentrations first alone (green bar) and then in a mixture with 5 µM strychnine (red bar plus green bar), showing differential strychnine-elicited inhibition of current and fluorescence</figcaption></figure>

depending on the glycine concentration (n = 5). Lower panel: fluorescence (cyan) and current (black) variations normalized to the fluorescence and current variations under 10 mM of glycine (mean and SEM (n = 5)). “ns” denotes that the ΔF (under glycine) and ΔF (under strychnine inhibition) are not significantly different; *P* < 0.05 (Unpaired two-sided student t-test).

<figure><img src="image_9.png" /><figcaption>**Fig. 6 | Characterization of R271C by VCF on C41V GlyR.** A Localization of the R271C mutation on the GlyR compared to the sensor K143W/Q219C (PDB:6PXD). B Upper panel: ΔI (cyan) and ΔF (black) dose-response curves (normalized to the glycine maximal response recorded on each individual oocyte) with mean and</figcaption></figure>

SEM error bars for R271C under glycine (solid line) and taurine (dotted line) application. Taurine elicits only 5.4 ± 2.3% of glycine maximum current, n = 5. Lower panel: ΔI (cyan) and ΔF (black) dose-response curves normalized to the maximum value with mean and SEM error bars for R271C under the taurine application.
<figure><img src="image_10.png" /><figcaption>**Fig. 7** | *Dynamic "personality" of the tau-closed state.* A Root-mean-square fluctuations (RMSF) from the average structure extracted from 0.5 µs MD simulations in explicit solvent/membrane are shown for three conformational states of GlyR captured by cryo-EM. Apo-closed corresponds to the resting state and tau-open to the active state of the receptor. For each state, RMSF values ranging from 0.5 to 4.5 Å are shown by the color (from blue to red) and the thickness (from thin to thick) of the sausage representation (B). State-based docking of strychnine. The boxplot defines the center at the median (red line) and mean (yellow dashed line),</figcaption></figure>

minima and maxima are indicated by whiskers, and the box represents data between the first and third quartiles. Relative strychnine-binding affinities for apo-closed, tau-closed, and tau-open was probed by docking strychnine to an ensemble of 500 protein snapshots sampled by MD and comparing the success rate of docking; i.e., a docking experiment was considered as successful if the docking score was within 10% of the score obtained in the strychnine-bound X-ray structure (stry-closed). The mean (yellow) or median (red) dashed lines show that strychnine binding in tau-closed is significantly stronger than in tau-open.

enhanced in tau-closed relative to apo-closed and tau-open. Taken together, the simulation results highlight that the tau-closed conformation features a surprising dynamic character, which is remarkably different from that of tau-open and apo-closed and could not be deduced from the comparison of static cryo-EM structures only.

To explore the correlation between MD and VCF data, we quantified the reorganization of the fluorescence quenching sensor, as well as the binding of strychnine.

First, we calculated the average distance between the Cα carbons of the VCF quenching sensor (corresponding to K143W/Q219C). It is 13.0 ± 0.7 Å in stry-closed and apo-closed and this value increases to 14.7 ± 1.4 Å and 16.4 ± 1.0 Å in tau-closed and tau-open, respectively. Therefore, the significant reorganization seen in the simulations of the sensor from the apo to the tau-closed conformation is compatible with a change in fluorescence during this transition. We speculate that, during the separation of the Cα carbons, the rather long and flexible TAMRA "side chain" anchored at C219 reorganizes to interact more frequently with W143, generating fluorescence quenching.

Second, the binding of strychnine to the various conformations of GlyR was explored by docking. As expected, favorable docking scores were obtained in stry-closed (–12.1 kcal/mol) and apo-closed (–10.3 kcal/mol) from cryo-EM. By contrast, both tau-open and tau-closed from cryo-EM feature an orthosteric site that is too small to accommodate strychnine and the docking scores were largely unfavorable. To account for the intrinsic flexibility of the protein, the docking experiments were repeated using protein conformations extracted from the MD trajectories. Assuming that a docking experiment was successful when the docking score was <–9.3 kcal/mol (i.e., within 10% of the score obtained in stry-closed), the success rate for docking strychnine was determined by counting the fraction of successful docking over 500 protein snapshots from MD of apo-closed, tau-closed, and tau-open. As shown in Fig. 7B, docking of strychnine was successful at 26% (std dev 6.5%) in apo-closed, 6.6% (std dev 7.0%), and with large deviations among replicas in tau-closed, and only 2.2% (std dev 1.4%) in tau-open. We conclude that the enhanced flexibility of the tau-closed structure makes it very different from tau-open and compatible with strychnine binding.

Although this analysis was carried out on the WT zebrafish receptor, we performed control MD simulations of the K143W/Q219C/C41V triple mutant in the tau-open and tau-closed states, giving similar results (Supplementary Fig. 9). The simulation results support the conclusion that the intermediate conformation(s)

revealed by VCF, concerning both (partial) agonists and glycine-strychnine combinations, is consistent with a highly dynamic cluster of states as described by MD simulation started from the cryo-EM tau-closed structure.

## Discussion

The VCF data on K143W/Q219C provide structural evidence for an early intermediate step in the global activation path of GlyR. Indeed, upon glycine application, they reveal the emergence of one or more intermediate conformation(s) characterized by a change in fluorescence but no current, followed by channel opening with no further change in fluorescence. Two lines of evidence support that the intermediate conformation(s) illuminated by VCF is on-path to the global activation process: (i) Receptors in this intermediate conformation(s) are activatable since they appear before the onset of the currents, and are thus not desensitized, (ii) Receptors in the intermediate conformation(s) are stabilized by partial agonists, low glycine concentrations, and mixtures of glycine and strychnine, showing a pharmacological profile "in-between" those stabilizing the resting and active states.

It is noteworthy that the mutations introduced here to anchor the fluorescence quenching sensor (K143W/Q219C) cause an 8-fold increase in EC₅₀<sup>current</sup> of glycine, and a decreased efficacy of the partial agonists to activate the receptor, yet with an intact conductance of the channel. This indicates a destabilization of the active state. Concerning fluorescence, the ΔF curve displays higher sensitivity than the ΔI curve, indicating that partial binding of glycine is sufficient to trigger the resting to intermediate(s) conformational transition. In addition, strychnine alone produces a fluorescence variation in the opposite direction to agonists, suggesting that the intermediate(s) conformation(s) are already substantially populated in the absence of ligand, with strychnine shifting the spontaneous equilibrium back toward the resting state⁴⁸. For partial agonists, they even mainly stabilize the intermediate(s) conformation at saturation, with a very minor representation of the active state.

Interestingly, we found that a non-activating concentration of propofol potentiates the currents without change of the ΔF curve, suggesting a specific action of propofol on the intermediate(s) to active transition. We also report that the gain of function mutant V280M does not show fluorescence changes. We suggest that this mutant could populate mainly the intermediate(s) state in the absence of an agonist.
<figure><img src="image_11.png" /><figcaption>**Fig. 8** | Hypothetical transition pathway of the GlyR. Left panel: side-view cartoon representation of two subunits of the GlyR. C loop contributing to the orthosteric site, the quenching sensor, the TAMRA (hexagon), and tryptophane (indole) are represented. Upper panel: top view representation of the receptor in different conformations with the ECD of each subunit shown in a circle and C loop as a line. Lower panel: top view representation of the TMD showing the M2 helices and the L9' that closes the pore in the resting and intermediate(s) conformations. The black circle represents positions where the fluorescent-quenching sensor is grafted. The indole goes closer to the TAMRA from the resting to the intermediate(s) conformations, generating a quenching of fluorescence.</figcaption></figure>

An important observation from VCF is that taurine binding stabilizes primarily the intermediate state(s). The structure of the “tau-closed” state of the zebrafish GlyR recently solved by cryo-EM is a strong candidate for an early on-pathway intermediate during the global activation process¹⁹. However, docking experiments show that this structure is unable to accommodate the bulky antagonist strychnine, a feature that is inconsistent with our VCF data. Interestingly, MD simulations of the tau-closed structure provide a way to explain this discrepancy. In fact, the simulations reveal that the ECD is surprisingly dynamic in the tau-closed state and samples the spontaneous opening of the orthosteric site on the sub-µs time scale, which makes it compatible with strychnine binding. Combined with the VCF data, these observations suggest that receptors in the presence of mixtures of strychnine and agonists would keep an overall intermediate-like conformation with significant asymmetry at the ECD and at the orthosteric sites. Some sites would bind agonists/partial agonists in an active-like compact state, while others would bind strychnine in a resting-like expanded state. We speculate that the same mechanism applies in the case of partial occupation of the orthosteric sites by glycine, accounting for a maximal fluorescence variation at low glycine concentration with little current, which explains the large separation of the ΔI and ΔF dose-response curves.

The simulations reveal that the tau-closed state presents a unique dynamic personality that is remarkably different from those at the resting and active states. In fact, a comparison of the atomic fluctuations in the different conformations reveals that the flexibility of the ECD-TMD interface is surprisingly enhanced in tau-closed relative to both apo-closed (resting-like) and tau-open (active-like). We speculate that the unique dynamic signature of tau-closed originates from the hybrid nature of this state, which features an ECD in a nearly active conformation along with a resting-like TMD endowed with a closed channel. This apparent mismatch at the ECD-TMD interface, which enhances the flexibility of the ECD domain and its blooming in tau-closed, favors the structural plasticity for promiscuous orthosteric-ligand binding. This hypothesis is consistent with the large reorganization of the pre-M1 and Cys-loop regions probed by VCF in the absence of currents and suggests that the tau-closed structure and its

dynamic character provide a reasonable representation of the early intermediate conformations found not only with partial agonists, but also with glycine at sub-saturation or with glycine and strychnine mixtures.

From combined VCF and MD data, we thus propose a speculative model of receptor activation depicted in Fig. 8, where the receptor transits through a cluster of highly dynamic intermediate conformations, allowing significant local asymmetry at the level of the orthosteric site (and also potentially in the upper ECD), but a more symmetric organization at the ECD-TMD interface where the fluorescent sensors are anchored. Although we cannot exclude the existence of alternative intermediate states in different pharmacological conditions e.g. involving changes in the local environment of the fluorescence sensors, the MD simulations of the cryo-EM Tau-closed state provide a reasonable and unifying explanation of the intermediate states recorded with (partial) agonists and the glycine-strychnine mixture.

Our data also suggest that a fraction of R271C receptors, when activated by taurine, display an intermediate-like conformation. It is noteworthy that R271C and K143W/Q219C cause a similar 10-fold increase in EC₅₀<sup>current</sup> of glycine (Supplementary Fig. 7). We thus speculate that their major difference relies on the resting to intermediate equilibrium, which is strongly shifted toward the intermediate(s) in K143W/Q219C and much less in R271C. Interestingly, K143W/Q219C and R271C are located at both extremities of the M2-M3 loop which undergoes a large outward movement during gating (Fig. 6A).

In regards to the WT zebrafish GlyR studied by cryo-EM, literature data suggest a contribution of around 7% of Tau-closed states at taurine saturation¹⁹. It is thus likely that, in the WT human GlyR context, the intermediate(s) state would also be sparsely populated in the presence of partial agonists. We speculate that the K143W mutation, strongly stabilizes the intermediate(s) conformation, thus allowing its robust monitoring by VCF. In retrospect, this might explain why the “intermediate” phenotype has rarely been observed in the previous VCF studies of the GlyR. Of note, previous single-channel works on nAChRs and GlyRs suggested the contribution of flipped/primed intermediate states on receptor activation. Flipped/primed are incompletely stabilized by partial agonists, underlying their limited efficacy to fully activate the receptor⁷,⁸. This is in contrast with the intermediate conformation(s) unraveled here, which is fully stabilized by partial agonists and consequently must occur earlier than the flip/prime transition.

Phenotypes resembling the one presented here have been already described for the bacterial channel GLIC and cationic pLGICs. With GLIC, we reported a “pre-activation” phenotype for a series of quenching sensors at the ECD interface and M2-M3 loop¹⁴,¹⁵. They all report fluorescence variations at non-activating agonist (proton) concentrations and the fluorescence variation was faster (ms range measured by stopped-flow) than the onset of activation (10–100 of ms measured by patch-clamp)⁴⁹. It is noteworthy that the M2-M3 loop is located just in between the quencher (W143) and the fluorophore (C219) in the GlyR structure (Fig. 1).

In the case of cationic pLGICs using VCF, labeling of the muscle nAChR (αγαδβ subtype) at the top of TMD (position 19' located at the upper end of M2 helix and facing the ECD) reports a 100-fold lower EC₅₀<sup>fluo</sup> as compared to EC₅₀<sup>current</sup> for ACh (binding the αδ interface), while no difference was observed for epibatidine (binding the αγ interface). Moreover, the change in fluorescence occurs with fast kinetics and pharmacological considerations support the occurrence of a conformational transition toward a singly liganded closed-channel state⁵⁰. Labeling at the 19' position has also been performed on the 5-HT₃R, consistent with monitoring a pre-active reorganization¹⁷. Finally, labeling of the α4β4 receptor in the extracellular loop 5 of α4 shows a fivefold lower EC₅₀<sup>fluo</sup> as compared to EC₅₀<sup>current</sup> for ACh suggesting that loop 5 moves before channel activation. Additionally, the antagonist DHBE elicits also a robust fluorescent change⁵¹ which
suggests a highly dynamic nature of loop 5 in the activation and inhi-
bition process. Altogether, data on the GlyR, GLIC, and possibly
nAChR/SHT3R provide evidence for intermediate transitions involving
structural reorganizations mainly at the ECD and/or ECD-TMD inter-
face. This idea is in line with REFERs analysis of the muscle nAChR, that
suggests that the orthosteric site and M2-M3 loop move very early
during the activation transition⁹. It is also consistent with simulations
describing the gating isomerization in terms of a combination of two
consecutive quaternary transitions named twisting and blooming². In
this view, the intermediate conformation(s) illuminated by VCF would
correspond to the completion of the un-blooming isomerization (i.e.,
compaction of the ECD) in a still closed-channel receptor. The appar-
ent absence of significant twisting during GlyR activation (Fig. S7) and
its role in gating remains to be understood. Altogether, while the
detailed activation mechanism may differ between different pLGIC
subtypes, the collected data highlight a general activation mechanism
in pLGICs characterized by progressive conformational propagation
from the ECD to the TMD.

In conclusion, this work illustrates the interest in combining VCF and MD approaches to identify functional intermediate conformations, characterize their flexibility, and annotate Cryo-EM structures to functional states. Noteworthy, such intermediate conformations may contribute to signal transduction in the postsynaptic membrane. The progressive transition also illuminates the mode of action of allosteric effectors and may be valuable for drug-design purposes.

Methods

Site-directed mutagenesis

The gene coding the full-length human glycine α1-subunit was cloned into the pMT3 plasmid containing an Ampicillin resistance gene. The PCR reaction was done with CloneAmp Hifi premix from Takara.

α1 homomeric glycine receptor expression in oocytes

*Xenopus Laevis* oocytes at stage VI are ordered from Portsmouth European Xenopus resource center and Ecocyte Biosciences and kept in Barth's solution (87.34 mM NaCl, 1 mM KCl, 0.66 mM CaNO3, 0.72 mM CaCl2, 0.82 mM MgSO4, 2.4 mM NaHCO3, 10 mM HEPES, and pH adjusted at 7.6 with NaOH). Ovary fragments obtained from Portsmouth European Xenopus resource center are treated as previously described⁵². cDNA coding the α1-subunit at 80 ng/µL is co-injected with a cDNA coding for GFP at 25 ng/µL into the oocyte nucleus by air injection. The injection pipettes are done with the capillary pipette (Nanoject II, Drummond) and PC-10 dual-stage glass micropipette puller. Injected oocytes are incubated at 18 °C for 3–4 days for expression.

Labeling of the cysteine

MTS-TAMRA and MTSR (Clinisciences) are dissolved in DMSO to obtain a final stock concentration of 10 mM and conserved at –20 °C. For labeling, oocytes expressing the receptors are incubated in 10 µM MTS-TAMRA (diluted in Barth's solution to obtain a final DMSO concentration of 0.1%) at 16–18 °C for 20 min before the recording. The procedure is similar for other fluorescent dyes (MTS-5-TAMRA, MTS-6-TAMRA). For MTSR, oocytes are incubated in 10 µM MTSR at 16–18 °C for less than 1 min.

Two-electrode voltage-clamp on oocytes

Oocytes expressing GlyR constructs are recorded in ND96 solution (96 mM NaCl, 2 mM KCl, 5 mM HEPES, 1 mM MgCl2, 1.8 mM CaCl2, and pH adjusted at 7.6 with NaOH). Molecules that are not fully soluble in an aqueous solution are first dissolved in DMSO and then diluted in ND96 to have a DMSO final concentration <0.1%. Currents were recorded by a Warner OC-725C amplifier and digitized by Digidata 1550 A and Clampex 10 (Molecular Devices). Currents were sampled at 500 Hz and filtered at 100 Hz. The voltage-clamp is maintained at

–60 mV at room temperature. Recording pipettes were made with Borosilicate glass with filament (BF150-110-7.5, WPI) with a PC-10 dual-stage glass micropipette puller to obtain pipettes of resistances comprised between 0.2 and 2 MΩ.

Voltage-clamp fluorometry on oocytes

Injected oocytes are placed in a home-designed chamber with the animal pole turning toward the microscope objective. Currents were recorded by GeneClamp 500 voltage patch-clamp amplifier (Axon Instruments) and digitized by Axon Digidata 1400 A digitizer and Clampex 10.6 software (Molecular Devices). TAMRA excitation is done by illumination with a 550 nm LED (pE-4000 CoolLED, 15–20% intensity) through an FF01-543_22 bandpass filter (Semrock). The fluorescence emission goes through an FF01-593_40 bandpass filter (Semrock) and is collected by a photo-multiplicator H10722 series (HAMAMATSU). Recording pipettes were the same as those prepared for TEVC. Currents were sampled at 2 kHz and filtered at 500 Hz. The voltage-clamp is maintained at –60 mV and at room temperature. The recording chamber has been designed to perfuse only the part of the oocyte from which the fluorescence emission is simultaneously collected (Fig. S3). This chamber has several advantages over chambers that allow to record the current of the entire membrane of the oocytes: (1) approximately the same population of receptors is recorded in current and fluorescence simultaneously, (2) the chamber allows for overexpression of the receptor on the oocyte membrane to increase the fluorescence signal without generating huge currents that make the clamp impossible to maintain during a prolonged period. The oocyte is placed in a hole and the geometry of the perfusion channels has been designed to “suck” the oocyte and stick it to the hole through a Venturi effect, ensuring efficient sealing between the two compartments. A limitation of the chamber is that the design is compatible with only relatively slow perfusion of the agonist, with rapid perfusion often expelling the oocyte out of the hole.

Voltage-clamp fluorometry data analysis

For each construct, we performed the experiments at least on two different batches of oocytes (oocytes from different ovary fragments and different animals) to obtain at least *n* = 5. Data were analyzed by Clampfit software (Molecular Devices) and are filtered by a boxcar filter. The baseline was corrected from leaking currents and measurements were done to the peak of each response. Dose-response curves were fitted to the Hill equation by Prism GraphPad software, and the error bars represent the SEM values. Rise time and decaytime analysis were calculated in Clampfit by monoexponential fits of the onset and the offset of the recordings. Statistical analyses were done using a Student *t*-test in Prism Graphpad.

Expression in cultured cells

Human Embryonic Kidney 293 (HEK-293, ATCC CRL-1573) cells were cultured in Dulbecco's modified Eagle's medium (DMEM) with 10% FBS (Invitrogen) in an incubator at 37°C and 5% CO2. After being PBS washed, trypsin-treated (Trypsine-EDTA; ThermoFisher Scientific), and seeded on Petri dishes, cells were transiently transfected using calcium phosphate-DNA co-precipitation with glycine receptor constructs (2 µg DNA) and a construct coding for a green fluorescent protein (0.2 µg). One day after transfection, cells were washed with fresh medium and recordings were carried out within 24 h.

Outside-out recordings and analysis

Recording currents are obtained with an RK-400 amplifier (BioLogic) using pClamp 10.5 software, digitized with a 1550 digidata (Axon instruments). Recording pipettes were obtained from thick-wall borosilicate glass (1.5 × 0.75 mm × 7.5 cm, Sutter Instrument) using a micropipettes puller (P-1000, Sutter Instrument) and fire-polished with a micro-forge (MF-830, Narishige) to be used at resistances
between 7 and 15 MΩ. Micropipettes were filled with internal solution (that contain in mM: 152 NaCl, 1 MgCl₂, 10 BAPTA, 10 HEPES; pH adjusted to 7.3 with NaOH solution, osmolarity measured at 335 mOsm). Extracellular solution (in mM: 152 NaCl, 1 MgCl₂, 10 HEPES; pH adjusted to 7.3 with NaOH solution and osmolarity was adjusted to 340 mOsm with glucose) was delivered by an automated perfusion system (RSC-200, BioLogic). Agonists' solutions are freshly made before sessions of recordings and are obtained with extracellular solution added with 1 to 100 µM of glycine (dissolved from a stock solution of 1 M in water). Acquisition of recordings was performed at the sampling of 20 kHz and low-pass filtered at 1 kHz (using the amplifier 5-pole Bessel filter). Openings are analyzed using Clampfit 10 software and currents were calculated by fitting the all-points histogram distributions of current amplitudes with the sum of two Gaussian curves. No further filtering is performed for the analysis.

## Molecular dynamics simulations

Molecular dynamics simulations of the zebrafish GlyR-α1 in four different conformational states were carried out: apo-closed (PDB:6PXD), Sy9-closed (PDB:6PXD with Sy9 docked), tau-closed (PDB:6PM3), and tau-open (PDB:6PM2). The zebrafish receptor shares more than 87% of sequence homology with the human receptor according to UNIPROT BLAST alignment. For all the studied systems, the protonation state at pH 7.0 was predicted using the Adaptive Poisson-Boltzmann Solver⁵³ (PropKa) via the web server ("https://server.poissonboltzmann.org/pdb2pqr"). Predicted protonation probabilities differed depending on the structure. To homogenize the systems, the most frequent predictions were selected and imposed on all systems: GLU316 was unprotonated, HIS125 protonated in Nδ, HIS217 protonated in Nδ, HIS327 protonated in Nδ, and HIS231 protonated in Nε. The structures were then embedded in a POPC bilayer and solvated with TIP3P water molecules using CHARMM-GUI Membrane Builder⁵⁴. The CHARMM-36m forcefield⁵⁵ for the protein and CGenFF⁵⁶ for the taurine and strychnine were used for the all-atom molecular dynamics simulations. The simulations were performed using GROMACS 2021.4⁷ with input files generated by CHARMM-GUI Membrane Builder for the minimization and equilibration. In short, 5000 steps of steepest descent were performed, then the system was equilibrated over 6 short simulations in the presence of atomic position restraints of decreasing strength, first in NVT, then in NPT ensemble, with the Berendsen thermostat and barostat. The temperature and the pressure were set to 300 K and 1 bar respectively. Finally, the equilibrated systems were carried on in production for 100 ns at a 2 fs timestep, in the NPT ensemble using V-rescale and Parrinello-Rahman as thermostat and barostat, respectively. In all cases, the LINCS algorithm was used to constrain the H bonds and the PME served to treat the long-range Coulomb interactions. For each system, we produced 5 replicas generated with different initial velocities from the equilibration.

MD simulations of the triple mutant used for the VCF experiments (K159W/Q235C/C57V, following the residue order in the PDB files) were also performed. Initial coordinates of the protein were extracted from the tau-closed (PDB: 6PM3) and tau-open (PDB: 6PM2) cryo-EM structures and the mutations were introduced using the CHARMM-GUI web server during preparation. The protocol followed to perform the simulations was the same as applied to the WT (see above).

## Trajectory analysis

The generated trajectories were analyzed using MDanalysis⁵⁸ through Python scripts. Notably, we computed the RMSD of the backbone atoms in the ECD, the TMD, and the structural core (i.e., the inner and outer beta-sheets in ECD, and the M1-M3 helices in TMD) relative to the initial coordinates from cryo-EM. Additionally, we computed the RMSF decomposed per residue and averaged over the 5 subunits of the

receptor. Analyses of the blooming and twisting angles¹⁰ were carried out using WORDOM software (https://sourceforge.net/p/wordom/codeh/ci/default/tree/). For these calculations, only a subset of residues was considered, referred to as the core, which is more relevant to monitor the quaternary structure of the receptor. The selected core residues are 54–69, 73–86, 113–125, 142–154, 165–175, and 207–214 for the ECD and residues 225–234, 241–258, 265–288, and 297–322 for the TMD. Plots were produced using Matplotlib⁵⁹, and the visual representations of the protein with the open-source version of PyMOL⁶⁰ (Schrödinger).

## Docking

To explore strychnine binding to the various conformational states of the receptor, i.e., apo-closed, tau-closed, and tau-open, we extracted 1 snapshot per nanosecond from the trajectory generated by molecular dynamics. By doing so, we ended up with 5 (replica) × 100 (ns) = 500 (snapshots) for each system. The molecular docking of strychnine was then executed using QuickVina 2.1⁶¹ with an exhaustiveness of 8 and a box of dimension 15 × 15 × 15 Å centered on the previously defined binding sites. For this purpose, we prepared the PDBQT inputs of the receptor using prepare_receptort4.py script from MGLTools⁶². Then, the positions of the 5 orthosteric binding sites of a given structure were computed on the fly as the center-of-mass of residues PHE223 of the (+)-subunit and PHE79 of the (-)-subunit using MDanalysis. Of note, for both tau-open and tau-closed, the receptor was simulated by Molecular Dynamics in the presence of taurine, but the docking of Strychnine was carried out without taurine.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

Data supporting the findings of this manuscript are available from the corresponding author upon reasonable request. PDB structures referred to in the course of this work: 6PXD, 6PM2, 6PM3, 6PM6. Source data are provided with this paper.

## References

1. Nemecz, Å., Prevost, M. S., Menny, A. & Corringer, P.-J. Emerging molecular mechanisms of signal transduction in pentameric ligand-gated ion channels. *Neuron* **90**, 452–470 (2016).

2. Cecchini, M. & Changeux, J.-P. The nicotinic acetylcholine receptor and its prokaryotic homologues: structure, conformational transitions & allosteric modulation. *Neuropharmacology* **96**, 137–149 (2015).

3. Cecchini, M. & Changeux, J.-P. Nicotinic receptors: from protein allostery to computational neuropharmacology. *Mol. Asp. Med.* **84**, 101044 (2022).

4. Indurthi, D. C. & Auerbach, A. Agonist efficiency from concentration-response curves: structural implications and applications. *Biophys. J.* **120**, 1800–1813 (2021).

5. Edelstein, S. J. & Changeux, J.-P. Relationships between structural dynamics and functional kinetics in oligomeric membrane receptors. *Biophys. J.* **98**, 2045–2052 (2010).

6. Gupta, S., Chakraborty, S., Vij, R. & Auerbach, A. A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip. *J. Gen. Physiol.* **149**, 85–103 (2017).

7. Lape, R., Colquhoun, D. & Sivilotti, L. G. On the nature of partial agonism in the nicotinic receptor superfamily. *Nature* **454**, 722–727 (2008).

8. Lee, W. Y. & Sine, S. M. Principal pathway coupling agonist binding to channel gating in nicotinic receptors. *Nature* **438**, 243–247 (2005).
9. Rovšnik, U. et al. Dynamic closed states of a ligand-gated ion channel captured by cryo-EM and simulations. *Life Sci. Alliance* **4**, e202101011 (2021).

10. Martin, N. E., Malik, S., Calimet, N., Changeux, J.-P. & Cecchini, M. Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics. *PLoS Comput. Biol.* **13**, e1005784 (2017).

11. Calimet, N. et al. A gating mechanism of pentameric ligand-gated ion channels. *Proc. Natl Acad. Sci. USA* **110**, E3987-E3996 (2013).

12. Bar-Lev, D. D., Degani-Katzav, N., Perelman, A. & Paas, Y. Molecular dissection of Cl--selective Cys-loop receptor points to components that are dispensable or essential for channel activity. *J. Biol. Chem.* **286**, 43830-43841 (2011).

13. Nury, H. et al. One-microsecond molecular dynamics simulation of channel gating in a nicotinic receptor homologue. *Proc. Natl Acad. Sci. USA* **107**, 6275-6280 (2010).

14. Menny, A. et al. Identification of a pre-active conformation of a pentameric channel receptor. *eLife* **6**, e23955 (2017).

15. Lefebvre, S. N., Taly, A., Menny, A., Medjebeur, K. & Corringer, P.-J. Mutational analysis to explore long-range allosteric couplings involved in a pentameric channel receptor pre-activation and activation. *eLife* **10**, e60682 (2021).

16. Basak, S. et al. High-resolution structures of multiple 5-HT3AR-setron complexes reveal a novel mechanism of competitive inhibition. *eLife* **9**, e57870 (2020).

17. Polovinkin, L. et al. Conformational transitions of the serotonin 5-HT3 receptor. *Nature* **563**, 275-279 (2018).

18. Masiulis, S. et al. GABAA receptor signalling mechanisms revealed by structural pharmacology. *Nature* **565**, 454-459 (2019).

19. Yu, J. et al. Mechanism of gating and partial agonist action in the glycine receptor. *Cell* **184**, 957-968.e21 (2021).

20. Rahman, M. M. et al. Structural mechanism of muscle nicotinic receptor desensitization and block by curare. *Nat. Struct. Mol. Biol.* **29**, 386-394 (2022).

21. Lape, R., Plested, A. J. R., Moroni, M., Colquhoun, D. & Sivilotti, L. G. The 1K276E startle disease mutation reveals multiple intermediate states in the gating of glycine receptors. *J. Neurosci.* **32**, 1336-1352 (2012).

22. Bormann, J., Hamill, O. P. & Sakmann, B. Mechanism of anion permeation through channels gated by glycine and gamma-aminobutyric acid in mouse cultured spinal neurones. *J. Physiol.* **385**, 243-286 (1987).

23. Du, J., Lü, W., Wu, S., Cheng, Y. & Gouaux, E. Glycine receptor mechanism elucidated by electron cryo-microscopy. *Nature* **526**, 224-229 (2015).

24. Talwar, S. & Lynch, J. W. Investigating ion channel conformational changes using voltage clamp fluorometry. *Neuropharmacology* **98**, 3-12 (2015).

25. Bode, A. et al. New hyperekplexia mutations provide insight into glycine receptor assembly, trafficking, and activation mechanisms. *J. Biol. Chem.* **288**, 33745-33759 (2013).

26. Han, L., Talwar, S., Wang, Q., Shan, Q. & Lynch, J. W. Phosphorylation of α3 glycine receptors induces a conformational change in the glycine-binding site. *ACS Chem. Neurosci.* **4**, 1361-1370 (2013).

27. Shan, Q., Han, L. & Lynch, J. W. β subunit M2-M3 loop conformational changes are uncoupled from α1 β glycine receptor channel gating: implications for human hereditary hyperekplexia. *PLoS ONE* **6**, e28105 (2011).

28. Wang, Q. & Lynch, J. W. Activation and desensitization induce distinct conformational changes at the extracellular-transmembrane domain interface of the glycine receptor. *J. Biol. Chem.* **286**, 38814-38824 (2011).

29. Pless, S. A. & Lynch, J. W. Ligand-specific conformational changes in the α1 glycine receptor ligand-binding domain. *J. Biol. Chem.* **284**, 15847-15856 (2009).

30. Pless, S. A. & Lynch, J. W. Magnitude of a conformational change in the glycine receptor β1-β2 loop is correlated with agonist efficacy. *J. Biol. Chem.* **284**, 27370-27376 (2009).

31. Pless, S. A., Dibas, M. I., Lester, H. A. & Lynch, J. W. Conformational variability of the glycine receptor M2 domain in response to activation by different agonists. *J. Biol. Chem.* **282**, 36057-36067 (2007).

32. Kumar, A. et al. Mechanisms of activation and desensitization of full-length glycine receptor in lipid nanodiscs. *Nat. Commun.* **11**, 3752 (2020).

33. Cerdan, A. H., Martin, N. É. & Cecchini, M. An ion-permeable state of the glycine receptor captured by molecular dynamics. *Structure* **26**, 1555-1562.e4 (2018).

34. Mansoor, S. E., DeWitt, M. A. & Farrens, D. L. Distance mapping in proteins using fluorescence spectroscopy: the tryptophan-induced quenching (TrlQ) method. *Biochemistry* **49**, 9722-9731 (2010).

35. Doose, S., Neuweiler, H. & Sauer, M. A close look at fluorescence quenching of organic dyes by tryptophan. *ChemPhysChem* **6**, 2277-2285 (2005).

36. Wilkins, M. E., Caley, A., Gielen, M. C., Harvey, R. J. & Smart, T. G. Murine startle mutant *Nmfr11* affects the structural stability of the glycine receptor and increases deactivation: Glycine receptor α1N46K and startle syndrome. *J. Physiol.* **594**, 3589-3607 (2016).

37. Langosch, D., Herbold, A., Schmieden, V., Borman, J. & Kirsch, J. Importance of Arg-219 for correct biogenesis of alpha 1 homopolymeric glycine receptors. *FEBS Lett.* **336**, 540-544 (1993).

38. Bode, A. & Lynch, J. W. Analysis of hyperekplexia mutations identifies transmembrane domain rearrangements that mediate glycine receptor activation. *J. Biol. Chem.* **288**, 33760-33771 (2013).

39. O'Shea, S. M. Propofol restores the function of 'hyperekplexic' mutant glycine receptors in xenopus oocytes and mice. *J. Neurosci.* **24**, 2322-2327 (2004).

40. Laube, B., Langosch, D., Betz, H. & Schmieden, V. Hyperekplexia mutations of the glycine receptor unmask the inhibitory subsite for β-amino-acids. *NeuroReport* **6**, 897-900 (1995).

41. Schmieden, V., Kuhse, J. & Betz, H. Agonist pharmacology of neonatal and adult glycine receptor alpha subunits: identification of amino acid residues involved in taurine activation. *EMBO J.* **11**, 2025-2032 (1992).

42. Jan, D. D. S., David-Watine, B., Korn, H. & Bregestovski, P. Activation of human α1 and α2 homomeric glycine receptors by taurine and GABA. *J. Physiol.* **535**, 741-755 (2001).

43. Schmieden, V. & Betz, H. Pharmacology of the inhibitory glycine receptor: agonist and antagonist actions of amino acids and piperidine carboxylic acid compounds. *Mol. Pharm.* **48**, 919-927 (1995).

44. Kim, J. J. & Hibbs, R. E. Direct structural insights into GABAA receptor pharmacology. *Trends Biochem. Sci.* **46**, 502-517 (2021).

45. Cerdan, A. H. et al. The glycine receptor allosteric ligands library (GRALL). *Bioinformatics* **36**, 3379-3384 (2020).

46. Fourati, Z. et al. Structural basis for a bimodal allosteric mechanism of general anesthetic modulation in pentameric ligand-gated ion channels. *Cell Rep.* **23**, 993-1004 (2018).

47. Nury, H. et al. X-ray structures of general anaesthetics bound to a pentameric ligand-gated ion channel. *Nature* **469**, 428-431 (2011).

48. Sauguet, L. et al. Crystal structures of a pentameric ligand-gated ion channel provide a mechanism for activation. *Proc. Natl Acad. Sci. USA* **111**, 966-971 (2014).

49. Laha, K. T., Ghosh, B. & Czajkowski, C. Macroscopic kinetics of pentameric ligand gated ion channels: comparisons between two prokaryotic channels and one eukaryotic channel. *PLoS ONE* **8**, e80322 (2013).
50. Dahan, D. S. et al. A fluorophore attached to nicotinic acetylcholine receptor M2 detects productive binding of agonist to the site. *Proc. Natl Acad. Sci. USA* **101**, 10195–10200 (2004).

51. Mourot, A., Bamberg, E. & Rettinger, J. Agonist- and competitive antagonist-induced movement of loop 5 on the α subunit of the neuronal α4β4 nicotinic acetylcholine receptor. *J. Neurochem.* **105**, 413–424 (2008).

52. Gielen, M., Barilone, N. & Corringer, P.-J. The desensitization pathway of GABAA receptors, one subunit at a time. *Nat. Commun.* **11**, 5369 (2020).

53. Jurrus, E. et al. Improvements to the *APBS* biomolecular solvation software suite. *Protein Sci.* **27**, 112–128 (2018).

54. Wu, E. L. et al. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. *J. Comput. Chem.* **35**, 1997–2004 (2014).

55. Huang, J. et al. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. *Nat. Methods* **14**, 71–73 (2017).

56. Vanommeslaeghe, K., Raman, E. P. & MacKerell, A. D. Automation of the CHARMM General Force Field (CGenFF) II: assignment of bonded parameters and partial atomic charges. *J. Chem. Inf. Model.* **52**, 3155–3168 (2012).

57. Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. *SoftwareX* **1–2**, 19–25 (2015).

58. Michaud-Agrawal, N., Denning, E. J., Woolf, T. B. & Beckstein, O. MDAnalysis: a toolkit for the analysis of molecular dynamics simulations. *J. Comput. Chem.* **32**, 2319–2327 (2011).

59. Hunter, J. D. Matplotlib: a 2D graphics environment. *Comput. Sci. Eng.* **9**, 90–95 (2007).

60. Schrödinger, L. L. C. The PyMOL molecular graphics system, version 2.5. (2015).

61. Alhossary, A., Handoko, S. D., Mu, Y. & Kwoh, C.-K. Fast, accurate, and reliablemolecular docking with QuickVina 2. *Bioinformatics* **31**, 2214–2216 (2015).

62. Morris, G. M. et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. *J. Comput. Chem.* **30**, 2785–2791 (2009).

## Acknowledgements

The work was supported by the ERC (Grant no. 788974, Dynacotine), by the "Agence Nationale de la Recherche" (Grant ANR-18-CE11-0015-01, Pentacontrol), Specific Grant Agreement No. 945539 (Human Brain Project SGA3), the doctoral school ED3C and the Foundation pour la Recherche Médicale (to Solène N. Lefebvre). The authors would like to thank Antoine Taly for their help in data interpretation, and Alexandre Mourot for critical reading of the manuscript.

## Author contributions

L.P. and A.C. contributed equally to the work. S.S. and S.N.L. set up and performed VCF experiments. L.P. performed single-channel experiments and some TEVC experiments and VCF experiments. M.G. designed and constructed the recording chamber. A.C. and P.M.R. performed MD and docking experiments. S.S., S.N.L., L.P., A.C., P.M.R., M.G., M.C., and P.-J.C. designed the study and analyzed the data. S.S., L.P., A.C., P.M.R., J.-P.C., M.C., and P.-J.C. wrote the manuscript. M.C., P.-J.C., and J.-P.C. acquired funding.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-36471-7.

**Correspondence** and requests for materials should be addressed to Marco Cecchini or Pierre-Jean Corringer.

**Peer review information** *Nature Communications* thanks Erik Lindahl and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023