Reversal of ion-charge selectivity renders the pentameric ligand-gated ion channel GLIC insensitive to anaesthetics

Resonance-driven ion transport and selectivity in prokaryotic ion channels

Ion channels exhibit a remarkably high accuracy in selecting uniquely its associated type of ion. The mechanisms behind ion selectivity are not well understood. Current explanations build mainly on molecular biology and bioinformatics. Here we propose a simple physical model for ion selectivity based on the driven damped harmonic oscillator (DDHO). The driving force for this oscillator is provided by self-organizing harmonic turbulent structures in the dehydrating ionic flow through the ion channel, namely, oscillating pressure waves in one dimension, and toroidal vortices in two and three dimensions. Density fluctuations caused by these turbulences efficiently transmit their energy to aqua ions that resonate with the driving frequency. Consequently, these release their hydration shell and exit the ion channel as free ions. Existing modeling frameworks do not express the required complex spatiotemporal dynamics, hence we introduce a macroscopic continuum model for ionic dehydration and transport, based on the hydrodynamics of a dissipative ionic flow through an ion channel, subject to electrostatic and amphiphilic interactions. This model combines three classical physical fields: Navier-Stokes equations from hydrodynamics, Gauss's law from Maxwell theory, and the convection-diffusion equation from continuum physics. Numerical experiments with mixtures of chemical species of ions in various degrees of hydration indeed reveal the emergence of strong oscillations in the ionic flow that are instrumental in the efficient dehydration and cause a strong ionic jet into the cell. As such, they provide a powerful engine for the DDHO mechanism. Theoretical predictions of the modeling framework match significantly with empirical patch-clamp data. The DDHO standard response curve defines a unique resonance frequency that depends on the mass and charge of the ion. In this way, the driving oscillations act as a selection mechanism that filters out one specific ion. Application of the DDHO model to real ion data shows that this mechanism indeed clearly distinguishes between chemical species and between aqua and bare ions with a large Mahalanobis distance and high oscillator quality. The DDHO framework helps to understand how SNP mutations can cause severe genetic pathologies as they destroy the geometry of the channel and so alter the resonance peaks of the required ion type.

Structural Basis for a Bimodal Allosteric Mechanism of General Anesthetic Modulation in Pentameric Ligand-Gated Ion Channels

Ion channel modulation by general anesthetics is a vital pharmacological process with implications for receptor biophysics and drug development. Functional studies have implicated conserved sites of both potentiation and inhibition in pentameric ligand-gated ion channels, but a detailed structural mechanism for these bimodal effects is lacking. The prokaryotic model protein GLIC recapitulates anesthetic modulation of human ion channels, and it is accessible to structure determination in both apparent open and closed states. Here, we report ten X-ray structures and electrophysiological characterization of GLIC variants in the presence and absence of general anesthetics, including the surgical agent propofol. We show that general anesthetics can allosterically favor closed channels by binding in the pore or favor open channels via various subsites in the transmembrane domain. Our results support an integrated, multi-site mechanism for allosteric modulation, and they provide atomic details of both potentiation and inhibition by one of the most common general anesthetics.

Propofol inhibits the voltage-gated sodium channel NaChBac at multiple sites

Voltage-gated sodium (NaV) channels are important targets of general anesthetics, including the intravenous anesthetic propofol. Electrophysiology studies on the prokaryotic NaV channel NaChBac have demonstrated that propofol promotes channel activation and accelerates activation-coupled inactivation, but the molecular mechanisms of these effects are unclear. Here, guided by computational docking and molecular dynamics simulations, we predict several propofol-binding sites in NaChBac. We then strategically place small fluorinated probes at these putative binding sites and experimentally quantify the interaction strengths with a fluorinated propofol analogue, 4-fluoropropofol. In vitro and in vivo measurements show that 4-fluoropropofol and propofol have similar effects on NaChBac function and nearly identical anesthetizing effects on tadpole mobility. Using quantitative analysis by 19F-NMR saturation transfer difference spectroscopy, we reveal strong intermolecular cross-relaxation rate constants between 4-fluoropropofol and four different regions of NaChBac, including the activation gate and selectivity filter in the pore, the voltage sensing domain, and the S4-S5 linker. Unlike volatile anesthetics, 4-fluoropropofol does not bind to the extracellular interface of the pore domain. Collectively, our results show that propofol inhibits NaChBac at multiple sites, likely with distinct modes of action. This study provides a molecular basis for understanding the net inhibitory action of propofol on NaV channels.

Recent progress on the molecular pharmacology of propofol

The precise mechanism by which propofol enhances GABAergic transmission remains unclear, but much progress has been made regarding the underlying structural and dynamic mechanisms. Furthermore, it is now clear that propofol has additional molecular targets, many of which are functionally influenced at concentrations achieved clinically. Focusing primarily on molecular targets, this brief review attempts to summarize some of this recent progress while pointing out knowledge gaps and controversies. It is not intended to be comprehensive but rather to stimulate further thought, discussion, and study on the mechanisms by which propofol produces its pleiotropic effects.

Structural Basis of Alcohol Inhibition of the Pentameric Ligand-Gated Ion Channel ELIC

The structural basis for alcohol modulation of neuronal pentameric ligand-gated ion channels (pLGICs) remains elusive. We determined an inhibitory mechanism of alcohol on the pLGIC Erwinia chrysanthemi (ELIC) through direct binding to the pore. X-ray structures of ELIC co-crystallized with 2-bromoethanol, in both the absence and presence of agonist, reveal 2-bromoethanol binding in the pore near T237(6') and the extracellular domain (ECD) of each subunit at three different locations. Binding to the ECD does not appear to contribute to the inhibitory action of 2-bromoethanol and ethanol as indicated by the same functional responses of wild-type ELIC and mutants. In contrast, the ELIC-α1β3GABA_AR chimera, replacing the ELIC transmembrane domain (TMD) with the TMD of α1β3GABA_AR, is potentiated by 2-bromoethanol and ethanol. The results suggest a dominant role of the TMD in modulating alcohol effects. The X-ray structures and functional measurements support a pore-blocking mechanism for inhibitory action of short-chain alcohols.

Common Anesthetic-binding Site for Inhibition of Pentameric Ligand-gated Ion Channels

BACKGROUND

Identifying functionally relevant anesthetic-binding sites in pentameric ligand-gated ion channels (pLGICs) is an important step toward understanding the molecular mechanisms underlying anesthetic action. The anesthetic propofol is known to inhibit cation-conducting pLGICs, including a prokaryotic pLGIC from Erwinia chrysanthemi (ELIC), but the sites responsible for functional inhibition remain undetermined.

METHODS

We photolabeled ELIC with a light-activated derivative of propofol (AziPm) and performed fluorine-19 nuclear magnetic resonance experiments to support propofol binding to a transmembrane domain (TMD) intrasubunit pocket. To differentiate sites responsible for propofol inhibition from those that are functionally irrelevant, we made an ELIC-γ-aminobutyric acid receptor (GABAAR) chimera that replaced the ELIC-TMD with the α1β3GABAAR-TMD and compared functional responses of ELIC-GABAAR and ELIC with propofol modulations.

RESULTS

Photolabeling showed multiple AziPm-binding sites in the extracellular domain (ECD) but only one site in the TMD with labeled residues M265 and F308 in the resting state of ELIC. Notably, this TMD site is an intrasubunit pocket that overlaps with binding sites for anesthetics, including propofol, found previously in other pLGICs. Fluorine-19 nuclear magnetic resonance experiments supported propofol binding to this TMD intrasubunit pocket only in the absence of agonist. Functional measurements of ELIC-GABAAR showed propofol potentiation of the agonist-elicited current instead of inhibition observed on ELIC.

CONCLUSIONS

The distinctly different responses of ELIC and ELIC-GABAAR to propofol support the functional relevance of propofol binding to the TMD. Combining the newly identified TMD intrasubunit pocket in ELIC with equivalent TMD anesthetic sites found previously in other cationic pLGICs, we propose this TMD pocket as a common site for anesthetic inhibition of pLGICs.

Evolution of Pentameric Ligand-Gated Ion Channels: Pro-Loop Receptors

Pentameric ligand-gated ion channels (pLGICs) are ubiquitous neurotransmitter receptors in Bilateria, with a small number of known prokaryotic homologues. Here we describe a new inventory and phylogenetic analysis of pLGIC genes across all kingdoms of life. Our main finding is a set of pLGIC genes in unicellular eukaryotes, some of which are metazoan-like Cys-loop receptors, and others devoid of Cys-loop cysteines, like their prokaryotic relatives. A number of such “Cys-less” receptors also appears in invertebrate metazoans. Together, those findings draw a new distribution of pLGICs in eukaryotes. A broader distribution of prokaryotic channels also emerges, including a major new archaeal taxon, Thaumarchaeota. More generally, pLGICs now appear nearly ubiquitous in major taxonomic groups except multicellular plants and fungi. However, pLGICs are sparsely present in unicellular taxa, suggesting a high rate of gene loss and a non-essential character, contrasting with their essential role as synaptic receptors of the bilaterian nervous system. Multiple alignments of these highly divergent sequences reveal a small number of conserved residues clustered at the interface between the extracellular and transmembrane domains. Only the “Cys-loop” proline is absolutely conserved, suggesting the more fitting name “Pro loop” for that motif, and “Pro-loop receptors” for the superfamily. The infered molecular phylogeny shows a Cys-loop and a Cys-less clade in eukaryotes, both containing metazoans and unicellular members. This suggests new hypotheses on the evolutionary history of the superfamily, such as a possible origin of the Cys-loop cysteines in an ancient unicellular eukaryote. Deeper phylogenetic relationships remain uncertain, particularly around the split between bacteria, archaea, and eukaryotes.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Direct Pore Binding as a Mechanism for Isoflurane Inhibition of the Pentameric Ligand-gated Ion Channel ELIC

Pentameric ligand-gated ion channels (pLGICs) are targets of general anesthetics, but molecular mechanisms underlying anesthetic action remain debatable. We found that ELIC, a pLGIC from Erwinia chrysanthemi, can be functionally inhibited by isoflurane and other anesthetics. Structures of ELIC co-crystallized with isoflurane in the absence or presence of an agonist revealed double isoflurane occupancies inside the pore near T237(6′) and A244(13′). A pore-radius contraction near the extracellular entrance was observed upon isoflurane binding. Electrophysiology measurements with a single-point mutation at position 6′ or 13′ support the notion that binding at these sites renders isoflurane inhibition. Molecular dynamics simulations suggested that isoflurane binding was more stable in the resting than in a desensitized pore conformation. This study presents compelling evidence for a direct pore-binding mechanism of isoflurane inhibition, which has a general implication for inhibitory action of general anesthetics on pLGICs.

Conformational Changes Underlying Desensitization of the Pentameric Ligand-Gated Ion Channel ELIC

Structural rearrangements underlying functional transitions of pentameric ligand-gated ion channels (pLGICs) are not fully understood. Using (19)F nuclear magnetic resonance and electron spin resonance spectroscopy, we found that ELIC, a pLGIC from Erwinia chrysanthemi, expanded the extracellular end and contracted the intracellular end of its pore during transition from the resting to an apparent desensitized state. Importantly, the contraction at the intracellular end of the pore likely forms a gate to restrict ion transport in the desensitized state. This gate differs from the hydrophobic gate present in the resting state. Conformational changes of the TM2-TM3 loop were limited to the N-terminal end. The TM4 helices and the TM3-TM4 loop appeared relatively insensitive to agonist-mediated structural rearrangement. These results indicate that conformational changes accompanying functional transitions are not uniform among different ELIC regions. This work also revealed the co-existence of multiple conformations for a given state and suggested asymmetric conformational arrangements in a homomeric pLGIC.