Glycine receptors (GlyR) mediate fast inhibitory neurotransmission by switching between discrete states in response to ligand-binding events. Recent high-resolution structures from cryoelectron microscopy (cryo-EM) and X-ray crystallography have provided atomistic models for the open and closed states. Notably, the cryo-EM structure in complex with glycine illuminated a previously unreported wide-open state, whose physiological significance is debated. Here, we present the structure of an ion-conducting state of GlyR α1 captured by molecular dynamics and validate its physiological relevance with computational electrophysiology and polyatomic anion permeation simulations. Our analysis suggests that none of the experimental structures is a true representation of the physiologically active state, although previously characterized open channels in GLIC at pH 4, or GluCl/GlyR with ivermectin bound, provide reasonable models. These results open the door to an original functional annotation and support the conclusion that pore closing by desensitization versus deactivation involves the reorientation of the pore-lining helices in opposite directions.

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter -glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The μs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the β1-β2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with -glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.