Structure of a potentially open state of a proton-activated pentameric ligand-gated ion channel

A bupropion modulatory site in the Gloeobacter violaceus ligand-gated ion channel

Bupropion is an atypical antidepressant and smoking cessation drug that causes adverse effects such as insomnia, irritability, and anxiety. Bupropion inhibits dopamine and norepinephrine reuptake transporters and eukaryotic cation-conducting pentameric ligand-gated ion channels, such as nicotinic acetylcholine and serotonin type 3A receptors, at clinically relevant concentrations. Here, we demonstrate that bupropion also inhibits a prokaryotic homolog of pentameric ligand-gated ion channels, the Gloeobacter violaceus ligand-gated ion channel (GLIC). Using the GLIC as a model, we used molecular docking to predict binding sites for bupropion. Bupropion was found to bind to several sites within the transmembrane domain, with the predominant site being localized to the interface between transmembrane segments M1 and M3 of two adjacent subunits. Residues W213, T214, and W217 in the first transmembrane segment, M1, and F267 and I271 in the third transmembrane segment, M3, most frequently reside within a 4 Å distance from bupropion. We then used single amino acid substitutions at these positions and two-electrode voltage-clamp recordings to determine their impact on bupropion inhibitory effects. The substitution T214F alters bupropion potency by shifting the half-maximal inhibitory concentration to a 13-fold higher value compared to wild-type GLIC. Residue T214 is found within a previously identified binding pocket for neurosteroids and lipids in the GLIC. This intersubunit binding pocket is structurally conserved and almost identical to a binding pocket described for neurosteroids in γ-aminobutyric acid type A receptors. Our data thus suggest that the T214 that lines a previously identified lipophilic binding pocket in GLIC and γ-aminobutyric acid type A receptors is also a modulatory site for bupropion interaction with the GLIC.

Cryo-EM structures of prokaryotic ligand-gated ion channel GLIC provide insights into gating in a lipid environment

GLIC, a proton-activated prokaryotic ligand-gated ion channel, served as a model system for understanding the eukaryotic counterparts due to their structural and functional similarities. Despite extensive studies conducted on GLIC, the molecular mechanism of channel gating in the lipid environment requires further investigation. Here, we present the cryo-EM structures of nanodisc-reconstituted GLIC at neutral and acidic pH in the resolution range of 2.6 - 3.4 Å. In our apo state at pH 7.5, the extracellular domain (ECD) displays conformational variations compared to the existing apo structures. At pH 4.0, three distinct conformational states (C1, C2 and O states) are identified. The protonated structures exhibit a compacted and counter-clockwise rotated ECD compared with our apo state. A gradual widening of the pore in the TMD is observed upon reducing the pH, with the widest pore in O state, accompanied by several layers of water pentagons. The pore radius and molecular dynamics (MD) simulations suggest that the O state represents an open conductive state. We also observe state-dependent interactions between several lipids and proteins that may be involved in the regulation of channel gating. Our results provide comprehensive insights into the importance of lipids impact on gating.

Short-chain mono-carboxylates as negative modulators of allosteric transitions in Gloeobacter violaceus ligand-gated ion channel, and impact of a pre-β5 strand (Loop Ω) double mutation on crotonate, not butyrate effect

Using the bacterial proton-activated pentameric receptor-channel Gloeobacter violaceus ligand-gated ion channel (GLIC): (1) We characterize saturated, mono-carboxylates as negative modulators of GLIC (as previously shown for crotonate; Alqazzaz et al., Biochemistry, 2016, 55, 5947). Butyrate and crotonate have indistinguishable properties regarding negative modulation of wt GLIC. (2) We identify a locus in the pre-β5 strand (Loop Ω) whose mutation inverses the effect of the mono-carboxylate crotonate from negative to positive modulation of the allosteric transitions, suggesting an involvement of the pre-β5 strand in coupling the extracellular orthotopic receptor to pore gating. (3) As an extension to the previously proposed "in series" mechanism, we suggest that a orthotopic/orthosteric site-vestibular site-Loop Ω-β5-β6 "sandwich"-Pro-Loop/Cys-Loop series may be an essential component of orthotopic/orthosteric compound-elicited gating control in this pentameric ligand-gated ion channel, on top of which compounds targeting the vestibular site may provide modulation.

Identification of a binding site for bupropion in Gloeobacter violaceus ligand-gated ion channel

Bupropion is an atypical antidepressant and smoking cessation drug which causes adverse effects such as insomnia, irritability, and anxiety. Bupropion inhibits dopamine and norepinephrine reuptake transporters and eukaryotic cation-conducting pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine (nACh) and serotonin type 3A (5-HT3A) receptors, at clinically relevant concentrations. However, the binding sites and binding mechanisms of bupropion are still elusive. To further understand the inhibition of pLGICs by bupropion, in this work, using a prokaryotic homologue of pLGICs as a model, we examined the inhibitory potency of bupropion in Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated ion channel. Bupropion inhibited proton-induced currents in GLIC with an inhibitory potency of 14.9 ± 2.0 μM, comparable to clinically attainable concentrations previously shown to also modulate eukaryotic pLGICs. Using single amino acid substitutions in GLIC and two-electrode voltage-clamp recordings, we further determined a binding site for bupropion in the lower third of the first transmembrane segment M1 at residue T214. The sidechain of M1 T214 together with additional residues of M1 and also of M3 of the adjacent subunit have previously been shown to contribute to binding of other lipophilic molecules like allopregnanolone and pregnanolone.

Structure-function analysis suggests that the photoreceptor LITE-1 is a light-activated ion channel

Sensation of light is essential for all organisms. The eye-less nematode Caenorhabditis elegans detects UV and blue light to evoke escape behavior. The photosensor LITE-1 absorbs UV photons with an unusually high extinction coefficient, involving essential tryptophans. Here, we modeled the structure and dynamics of LITE-1 using AlphaFold2-multimer and molecular dynamics (MD) simulations and performed mutational and behavioral assays in C. elegans to characterize its function. LITE-1 resembles olfactory and gustatory receptors from insects, recently shown to be tetrameric ion channels. We identified residues required for channel gating, light absorption, and mechanisms of photo-oxidation, involving a likely binding site for the peroxiredoxin PRDX-2. Furthermore, we identified the binding pocket for a putative chromophore. Several residues lining this pocket have previously been established as essential for LITE-1 function. A newly identified critical cysteine pointing into the pocket represents a likely chromophore attachment site. We derived a model for how photon absorption, via a network of tryptophans and other aromatic amino acids, induces an excited state that is transferred to the chromophore. This evokes conformational changes in the protein, possibly leading to a state receptive to oxidation of cysteines and, jointly, to channel gating. Electrophysiological data support the idea that LITE-1 is a photon and H2O2-coincidence detector. Other proteins with similarity to LITE-1, specifically C. elegans GUR-3, likely use a similar mechanism for photon detection. Thus, a common protein fold and assembly, used for chemoreception in insects, possibly by binding of a particular compound, may have evolved into a light-activated ion channel.

Ligand-Dependent Intramolecular Motion of Native Nicotinic Acetylcholine Receptors Determined in Living Myotube Cells via Diffracted X-ray Tracking

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that play an important role in signal transduction at the neuromuscular junction (NMJ). Movement of the nAChR extracellular domain following agonist binding induces conformational changes in the extracellular domain, which in turn affects the transmembrane domain and opens the ion channel. It is known that the surrounding environment, such as the presence of specific lipids and proteins, affects nAChR function. Diffracted X-ray tracking (DXT) facilitates measurement of the intermolecular motions of receptors on the cell membranes of living cells, including all the components involved in receptor function. In this study, the intramolecular motion of the extracellular domain of native nAChR proteins in living myotube cells was analyzed using DXT for the first time. We revealed that the motion of the extracellular domain in the presence of an agonist (e.g., carbamylcholine, CCh) was restricted by an antagonist (i.e., alpha-bungarotoxin, BGT).

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Tutorial for Stopped-Flow Water Flux Measurements: Why a Report about “Ultrafast Water Permeation through Nanochannels with a Densely Fluorous Interior Surface” Is Flawed

Millions of years of evolution have produced proteinaceous water channels (aquaporins) that combine perfect selectivity with a transport rate at the edge of the diffusion limit. However, Itoh et al. recently claimed in Science that artificial channels are 100 times faster and almost as selective. The published deflation kinetics of vesicles containing channels or channel elements indicate otherwise, since they do not demonstrate the facilitation of water transport. In an illustrated tutorial on the experimental basis of stopped-flow measurements, we point out flaws in data processing. In contrast to the assumption voiced in Science, individual vesicles cannot simultaneously shrink with two different kinetics. Moreover, vesicle deflation within the dead time of the instrument cannot be detected. Since flawed reports of ultrafast water channels in Science are not a one-hit-wonder as evidenced by a 2018 commentary by Horner and Pohl in Science, we further discuss the achievable limits of single-channel water permeability. After analyzing (i) diffusion limits for permeation through narrow channels and (ii) hydrodynamics in the surrounding reservoirs, we conclude that it is unlikely to fundamentally exceed the evolutionarily optimized water-channeling performance of the fastest aquaporins while maintaining near-perfect selectivity.

Origin of acetylcholine antagonism in ELIC, a bacterial pentameric ligand-gated ion channel

ELIC is a prokaryotic homopentameric ligand-gated ion channel that is homologous to vertebrate nicotinic acetylcholine receptors. Acetylcholine binds to ELIC but fails to activate it, despite bringing about conformational changes indicative of activation. Instead, acetylcholine competitively inhibits agonist-activated ELIC currents. What makes acetylcholine an agonist in an acetylcholine receptor context, and an antagonist in an ELIC context, is not known. Here we use available structures and statistical coupling analysis to identify residues in the ELIC agonist-binding site that contribute to agonism. Substitution of these ELIC residues for their acetylcholine receptor counterparts does not convert acetylcholine into an ELIC agonist, but in some cases reduces the sensitivity of ELIC to acetylcholine antagonism. Acetylcholine antagonism can be abolished by combining two substitutions that together appear to knock out acetylcholine binding. Thus, making the ELIC agonist-binding site more acetylcholine receptor-like, paradoxically reduces the apparent affinity for acetylcholine, demonstrating that residues important for agonist binding in one context can be deleterious in another. These findings reinforce the notion that although agonism originates from local interactions within the agonist-binding site, it is a global property with cryptic contributions from distant residues. Finally, our results highlight an underappreciated mechanism of antagonism, where agonists with appreciable affinity, but negligible efficacy, present as competitive antagonists.

A structural and functional study of the prokaryotic ligand-gated ion channel, ELIC, provides insight into the origin of agonism and antagonism at nicotinic acetylcholine receptors.

Expression of Mutant Glycine Receptors in Xenopus Oocytes Using Canonical and Non-Canonical Amino Acids Reveals Distinct Roles of Conserved Proline Residues

Pentameric ligand-gated ion channels (pLGIC) play important roles in fast neuronal signal transmission. Functional receptors are pentamers, with each subunit having an extracellular domain (ECD), a transmembrane domain (TMD) and an intracellular domain. The binding of the agonist to the ECD induces a structural change that is transduced to the TMD to open the channel. Molecular details of this process are emerging, but a comprehensive understanding is still lacking. Proline (Pro) is one amino acid that has attracted much interest; its unusual features generate bends in loops and kinks and bulges in helices, which can be essential for function in some pLGICs. Here, we explore the roles of four conserved Pros in the glycine receptor (GlyR), creating substitutions with canonical and noncanonical amino acids, characterizing them using two electrode voltage clamp electrophysiology in Xenopus oocytes, and interpreting changes in receptor parameters using structural data from the open and closed states of the receptor. The data reveal that for efficient function, the Pro in the α1β1 loop is needed to create a turn and to be the correct size and shape to interact with nearby residues; the peptide bond of the Pro in the Cys-loop requires the cis conformation; and the Pros in loop A and M1 allow efficient function because of their reduced hydrogen bonding capacity. These data are broadly consistent with data from other pLGICs, and therefore likely represent the important features of these Pros in all members of the family.

Carboxylate binding prefers two cations to one

Almost all studies of specific ion binding by carboxylates (-COO^-) have considered only a single cation, but clustering of ions and ligands is a common phenomenon. We apply density functional theory to investigate how variations in the number of acetate ligands in binding to two monovalent cations affects ion binding preferences. We study a series of monovalent (Li⁺, Na⁺, K⁺, Cs⁺) ions relevant to experimental work on many topics, including ion channels, battery storage, water purification and solar cells. We find that the preferred optimal structure has 3 acetates except for Cs⁺, which has 2 acetates. The optimal coordination of the cation by the carboxylate O atoms is 4 for both Na⁺ and K⁺, and 3 for Li⁺ and Cs⁺. There is a 4-fold coordination minimum just a few kcal mol^-1 higher than the optimal 3-fold structure for Li⁺. For two cations, multiple minima occur in the vicinity of the lowest free energy state. We find that, for Li, Na and K, the preferred optimal structure with two cations is favored over a mixture of single cation complexes, providing a basis for understanding ionic cluster formation that is relevant for engineering proteins and other materials for rapid, selective ion transport.

Potassium Binding Interactions with Aliphatic Amino Acids: Thermodynamic and Entropic Effects Analyzed via a Guided Ion Beam and Computational Study

Noncovalent interactions between alkali metals and amino acids are critical for many biological processes, especially for proper function of protein ion channels; however, many precise binding affinities between alkali metals and amino acids still need to be measured. This study addresses this need by using threshold collision-induced dissociation with a guided ion beam tandem mass spectrometer to measure binding affinities between potassium cations and the aliphatic amino acids: Gly, Ala, hAla, Val, Leu, and Ile. These measurements are supplemented by theoretical calculations and include commentary on effects of enthalpy, entropy, and structural preference. Notably, all levels of theory indicate that the lowest-lying isomers at 298 K have K⁺ binding to the carbonyl oxygen in either a monodentate ([CO]) or bidentate ([CO,OH]) fashion, isomers that are linked in a double-well potential. This complicates the analysis of the data, although does not greatly influence the final results. Analysis of the resulting cross sections includes accounting for multiple ion-molecule collisions, internal energy of reactant ions, and unimolecular decay rates. The resulting experimental bond dissociation energies generally increase as the polarizability of the amino acid increases, results that agree well with quantum chemical calculations done at the B3LYP, B3P86, and MP2(full) levels of theory, with B3LYP-GD3BJ predicting systematically larger values.

Speculation on How RIC-3 and Other Chaperones Facilitate α7 Nicotinic Receptor Folding and Assembly

The process of how multimeric transmembrane proteins fold and assemble in the endoplasmic reticulum is not well understood. The alpha7 nicotinic receptor (α7 nAChR) is a good model for multimeric protein assembly since it has at least two independent and specialized chaperones: Resistance to Inhibitors of Cholinesterase 3 (RIC-3) and Nicotinic Acetylcholine Receptor Regulator (NACHO). Recent cryo-EM and NMR data revealed structural features of α7 nAChRs. A ser-ala-pro (SAP) motif precedes a structurally important but unique "latch" helix in α7 nAChRs. A sampling of α7 sequences suggests the SAP motif is conserved from C. elegans to humans, but the latch sequence is only conserved in vertebrates. How RIC-3 and NACHO facilitate receptor subunits folding into their final pentameric configuration is not known. The artificial intelligence program AlphaFold2 recently predicted structures for NACHO and RIC-3. NACHO is highly conserved in sequence and structure across species, but RIC-3 is not. This review ponders how different intrinsically disordered RIC-3 isoforms from C. elegans to humans interact with α7 nAChR subunits despite having little sequence homology across RIC-3 species. Two models from the literature about how RIC-3 assists α7 nAChR assembly are evaluated considering recent structural information about the receptor and its chaperones.

General anesthesia bullies the gut: a toxic relationship with dysbiosis and cognitive dysfunction

Perioperative neurocognitive disorder (PND) is a common surgery outcome affecting up to a third of the elderly patients, and it is associated with high morbidity and increased risk for Alzheimer's disease development. PND is characterized by cognitive impairment that can manifest acutely in the form of postoperative delirium (POD) or after hospital discharge as postoperative cognitive dysfunction (POCD). Although POD and POCD are clinically distinct, their development seems to be mediated by a systemic inflammatory reaction triggered by surgical trauma that leads to dysfunction of the blood-brain barrier and facilitates the occurrence of neuroinflammation. Recent studies have suggested that the gut microbiota composition may play a pivotal role in the PND development by modulating the risk of neuroinflammation establishment. In fact, modulation of gut microbiome composition with pre- and probiotics seems to be effective for the prevention and treatment of PND in animals. Interestingly, general anesthetics seem to have major responsibility on the gut microbiota composition changes following surgery and, consequently, can be an important element in the process of PND initiation. This concept represents an important milestone for the understanding of PND pathogenesis and may unveil new opportunities for the development of preventive or mitigatory strategies against the development of these conditions. The aim of this review is to discuss how anesthetics used in general anesthesia can interact and alter the gut microbiome composition and contribute to PND development by favoring the emergence of neuroinflammation.

The Hair Cell α9α10 Nicotinic Acetylcholine Receptor: Odd Cousin in an Old Family

Nicotinic acetylcholine receptors (nAChRs) are a subfamily of pentameric ligand-gated ion channels with members identified in most eumetazoan clades. In vertebrates, they are divided into three subgroups, according to their main tissue of expression: neuronal, muscle and hair cell nAChRs. Each receptor subtype is composed of different subunits, encoded by paralogous genes. The latest to be identified are the α9 and α10 subunits, expressed in the mechanosensory hair cells of the inner ear and the lateral line, where they mediate efferent modulation. α9α10 nAChRs are the most divergent amongst all nicotinic receptors, showing marked differences in their degree of sequence conservation, their expression pattern, their subunit co-assembly rules and, most importantly, their functional properties. Here, we review recent advances in the understanding of the structure and evolution of nAChRs. We discuss the functional consequences of sequence divergence and conservation, with special emphasis on the hair cell α9α10 receptor, a seemingly distant cousin of neuronal and muscle nicotinic receptors. Finally, we highlight potential links between the evolution of the octavolateral system and the extreme divergence of vertebrate α9α10 receptors.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Probing solution structure of the pentameric ligand-gated ion channel GLIC by small-angle neutron scattering

Pentameric ligand-gated ion channels undergo subtle conformational cycling to control electrochemical signal transduction in many kingdoms of life. Several crystal structures have now been reported in this family, but the functional relevance of such models remains unclear. Here, we used small-angle neutron scattering (SANS) to probe ambient solution-phase properties of the pH-gated bacterial ion channel GLIC under resting and activating conditions. Data collection was optimized by inline paused-flow size-exclusion chromatography, and exchanging into deuterated detergent to hide the micelle contribution. Resting-state GLIC was the best-fit crystal structure to SANS curves, with no evidence for divergent mechanisms. Moreover, enhanced-sampling molecular-dynamics simulations enabled differential modeling in resting versus activating conditions, with the latter corresponding to an intermediate ensemble of both the extracellular and transmembrane domains. This work demonstrates state-dependent changes in a pentameric ion channel by SANS, an increasingly accessible method for macromolecular characterization with the coming generation of neutron sources.

Regulation of a pentameric ligand-gated ion channel by a semiconserved cationic lipid-binding site

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of noncanonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid-binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABAA receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABAA receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Dynamic closed states of a ligand-gated ion channel captured by cryo-EM and simulations

Ligand-gated ion channels are critical mediators of electrochemical signal transduction across evolution. Biophysical and pharmacological characterization of these receptor proteins relies on high-quality structures in multiple, subtly distinct functional states. However, structural data in this family remain limited, particularly for resting and intermediate states on the activation pathway. Here, we report cryo-electron microscopy (cryo-EM) structures of the proton-activated Gloeobacter violaceus ligand-gated ion channel (GLIC) under three pH conditions. Decreased pH was associated with improved resolution and side chain rearrangements at the subunit/domain interface, particularly involving functionally important residues in the β1-β2 and M2-M3 loops. Molecular dynamics simulations substantiated flexibility in the closed-channel extracellular domains relative to the transmembrane ones and supported electrostatic remodeling around E35 and E243 in proton-induced gating. Exploration of secondary cryo-EM classes further indicated a low-pH population with an expanded pore. These results allow us to define distinct protonation and activation steps in pH-stimulated conformational cycling in GLIC, including interfacial rearrangements largely conserved in the pentameric channel family.

Important amino acids for function of the insect Rdl GABA receptor

BACKGROUND

Insect Rdl GABA receptor is an important insecticide target. To design a novel insecticide, studies on the structures of homologous pentameric ligand-gated ion channels provide information about important amino acids that are necessary for the function of insect Rdl GABA receptors.

RESULTS

L9'A, T12'A, T13'A, T13'S, M15'S, and M15'N mutations in the Drosophila Rdl GABA receptor subunit caused the protein to spontaneously adopt the open state conformation. In contrast, the S16'A, S16'T, S17'A, and S17'H mutant homomers showed the same levels of agonist and antagonist sensitivity as the wild-type receptor. The G336M mutation in the Drosophila Rdl GABA receptor abolished the agonist activities of ivermectin and milbemectin, but the F339M mutation did not. Additionally, the F339M mutation caused spontaneous opening of the receptor. In the Drosophila Rdl model, the hydrophobic girdle plays an important role in stabilization of the closed state. Mutations which decrease hydrophobic interactions resulted in spontaneous opening, supporting the importance of the hydrophobic girdle for keeping the channel closed. Through a mutational study of transmembrane 3 (TM3) cytoplasmic domain and Rdl GABA receptor modeling, hydrophobic interactions between TM3 and TM4 and intersubunit interaction were demonstrated to be important for channel gating. Alternatively, the intrasubunit interaction between TM2 and TM3 domains were less important for channel gating in case of Drosophila Rdl GABA receptor.

CONCLUSION

This study demonstrates important amino acids critical to the function of the Drosophila Rdl GABA receptor based on the mutational studies and Drosophila Rdl GABA receptor modeling approach. © 2020 Society of Chemical Industry.

Elephants in the Dark: Insights and Incongruities in Pentameric Ligand-gated Ion Channel Models

The superfamily of pentameric ligand-gated ion channels (pLGICs) comprises key players in electrochemical signal transduction across evolution, including historic model systems for receptor allostery and targets for drug development. Accordingly, structural studies of these channels have steadily increased, and now approach 250 depositions in the protein data bank. This review contextualizes currently available structures in the pLGIC family, focusing on morphology, ligand binding, and gating in three model subfamilies: the prokaryotic channel GLIC, the cation-selective nicotinic acetylcholine receptor, and the anion-selective glycine receptor. Common themes include the challenging process of capturing and annotating channels in distinct functional states; partially conserved gating mechanisms, including remodeling at the extracellular/transmembrane-domain interface; and diversity beyond the protein level, arising from posttranslational modifications, ligands, lipids, and signaling partners. Interpreting pLGIC structures can be compared to describing an elephant in the dark, relying on touch alone to comprehend the many parts of a monumental beast: each structure represents a snapshot in time under specific experimental conditions, which must be integrated with further structure, function, and simulations data to build a comprehensive model, and understand how one channel may fundamentally differ from another.

Distinct Evolutionary Trajectories of Neuronal and Hair Cell Nicotinic Acetylcholine Receptors

The expansion and pruning of ion channel families has played a crucial role in the evolution of nervous systems. Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels with distinct roles in synaptic transmission at the neuromuscular junction, the central and peripheral nervous system, and the inner ear. Remarkably, the complement of nAChR subunits has been highly conserved along vertebrate phylogeny. To ask whether the different subtypes of receptors underwent different evolutionary trajectories, we performed a comprehensive analysis of vertebrate nAChRs coding sequences, mouse single-cell expression patterns, and comparative functional properties of receptors from three representative tetrapod species. We found significant differences between hair cell and neuronal receptors that were most likely shaped by the differences in coexpression patterns and coassembly rules of component subunits. Thus, neuronal nAChRs showed high degree of coding sequence conservation, coupled to greater coexpression variance and conservation of functional properties across tetrapod clades. In contrast, hair cell α9α10 nAChRs exhibited greater sequence divergence, narrow coexpression pattern, and great variability of functional properties across species. These results point to differential substrates for random change within the family of gene paralogs that relate to the segregated roles of nAChRs in synaptic transmission.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Sequential purification and characterization of Torpedo californica nAChR-DC supplemented with CHS for high-resolution crystallization studies

Over the past 10 years we have been developing a multi-attribute analytical platform that allows for the preparation of milligram amounts of functional, high-pure, and stable Torpedo (muscle-type) nAChR detergent complexes for crystallization purpose. In the present work, we have been able to significantly improve and optimize the purity and yield of nicotinic acetylcholine receptors in detergent complexes (nAChR-DC) without compromising stability and functionality. We implemented new methods in the process, such as analysis and rapid production of samples for future crystallization preparations. Native nAChR was extracted from the electric organ of Torpedo californica using the lipid-like detergent LysoFos Choline 16 (LFC-16), followed by three consecutive steps of chromatography purification. We evaluated the effect of cholesteryl hemisuccinate (CHS) supplementation during the affinity purification steps of nAChR-LFC-16 in terms of receptor secondary structure, stability and functionality. CHS produced significant changes in the degree of β-secondary structure, these changes compromise the diffusion of the nAChR-LFC-16 in lipid cubic phase. The behavior was reversed by Methyl-β-Cyclodextrin treatment. Also, CHS decreased acetylcholine evoked currents of Xenopus leavis oocyte injected with nAChR-LFC-16 in a concentration-dependent manner. Methyl-β-Cyclodextrin treatment do not reverse functionality, however column delipidation produced a functional protein similar to nAChR-LFC-16 without CHS treatment.

Mobility of Lower MA-Helices for Ion Conduction through Lateral Portals in 5-HT3A Receptors

The intracellular domain of the serotonin type 3A receptor, a pentameric ligand-gated ion channel, is crucial for regulating conductance. Ion permeation through the extracellular vestibule and the transmembrane channel is well understood, whereas the specific ion conduction pathway through the intracellular domain is less clear. The intracellular domain starts with a short loop after the third transmembrane segment, followed by a short α-helical segment, a large unstructured loop, and finally, the membrane-associated MA-helix that continues into the last transmembrane segment. The MA-helices from all five subunits form the extension of the transmembrane ion channel and shape what has been described as a “closed vestibule,” with their lateral portals obstructed by loops and their cytosolic ends forming a tight hydrophobic constriction. The question remains whether the lateral portals or cytosolic constriction conduct ions upon channel opening. In our study, we used disulfide bond formation between pairs of engineered cysteines to probe the proximity and mobility of segments of the MA-helices most distal to the membrane bilayer. Our results indicate that the proximity and orientation for cysteine pairs at I409C/R410C, in close proximity to the lateral windows, and L402C/L403C, at the cytosolic ends of the MA-helices, are conducive for disulfide bond formation. Although conformational changes associated with gating promote cross-linking for I409C/R410C, which in turn decreases channel currents, cross-linking of L402C/L403C is functionally silent in macroscopic currents. These results support the hypothesis that concerted conformational changes open the lateral portals for ion conduction, rendering ion conduction through the vertical portal unlikely.

Pair of Residue Substitutions at the Outer Mouth of the Channel Pore Act as Inputs for a Boolean Logic "OR" Gate Based on the Glycine Receptor

The glycine receptor (GlyR) is a ligand-activated chloride channel, whose mutations are the major cause of hereditary hyperekplexia. The hyperekplexia-causing R271Q mutation, which is located at the extracellular outer mouth of the channel pore, dramatically impairs the GlyR function manifesting a reduced sensitivity toward glycine. This study reports that a second mutation, S273D, rescues the function of the R271Q GlyR to that of the wild-type (WT) GlyR. Surprisingly, the S273D mutation, when introduced to the WT GlyR, does not further increase the receptor function. In other words, the compromised function of the 271Q 273S GlyR (i.e., the R271Q GlyR) can be rescued to WT levels by the introduction of either, or both, of the Q271R and S273D substitutions. From the perspective of Boolean logic gates, the Q271R and S273D substitutions act as inputs for an OR gate based on the GlyR. Further experiments revealed that the negative-charge carried by the 273 residue is essential for the expression of the OR gate and that the expression of the OR gate is residue-position-specific. In addition, mechanistic investigation implied that the 273 residue influences the 271 residue, which might underpin the unique nonadditive OR gate relationship between these two residues. Such an ion-channel-based OR gate, expressing output in the form of electrical current, could potentially be developed to digitally manipulate neuronal activity.

Charged residues at the pore extracellular half of the glycine receptor facilitate channel gating: a potential role played by electrostatic repulsion

KEY POINTS

The Arg271Gln mutation of the glycine receptor (GlyR) causes hereditary hyperekplexia. This mutation dramatically compromises GlyR function; however, the underlying mechanism is not yet known. This study, by employing function and computation methods, proposes that charged residues (including the Arg residue) at the pore extracellular half from each of the five subunits of the homomeric α1 GlyR, create an electrostatic repulsive potential to widen the pore, thereby facilitating channel opening. This mechanism explains how the Arg271Gln mutation, in which the positively charged Arg residue is substituted by the neutral Gln residue, compromises GlyR function. This study furthers our understanding of the biophysical mechanism underlying the Arg271Gln mutation compromising GlyR function.

ABSTRACT

The R271(19')Q mutation in the α1 subunit of the glycine receptor (GlyR) chloride channel causes hereditary hyperekplexia. This mutation dramatically compromises channel function; however, the underlying mechanism is not yet known. The R271 residue is located at the extracellular half of the channel pore. In this study, an Arg-scanning mutagenesis was performed at the pore extracellular half from the 262(10') to the 272(20') position on the background of the α1 GlyR carrying the hyperekplexia-causing mutation R271(19')Q. It was found that the placement of the Arg residue rescued channel function to an extent inversely correlated with the distance between the residue and the pore central axis (perpendicular to the plane of the lipid bilayer). Accordingly, it was hypothesized that the placed Arg residues from each of the five subunits of the homomeric α1 GlyR create an electrostatic repulsive potential to widen the pore, thereby facilitating channel opening. This hypothesis was quantitatively verified by theoretical computation via exploiting basic laws of electrostatics and thermodynamics, and further supported by more experimental findings that the placement of another positively charged Lys residue or even a negatively charged Asp residue also rescued channel function in the same manner. This study provides a novel mechanism via which charged residues in the pore region facilitate channel gating, not only for the disease-causing 19'R residue in the GlyR, but also potentially for charged residues in the same region of other ion channels.

In silico studies of the interactions between propofol and fentanyl using Gaussian accelerated molecular dynamics

Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.

Many Proline Residues in the Extracellular Domain Contribute to Glycine Receptor Function

Prolines in signaling proteins are of particular interest because they have a range of unique properties that may be critical for function. Here we show that many proline residues in the extracellular domain (ECD) of the glycine receptor are involved in the correct functioning of this ligand-gated ion channel. We explore their role by creating mutant receptors, expressing them in cells, and using fluorescent membrane potential sensitive dye to monitor receptor activity. We then interpret the changes in receptor parameters using structural information from the open and closed states of the receptor. The data reveal that substitution with alanine of ten of the 13 Pro residues in the ECD alters the function of the receptor: one substitution ablates function, six cause a decrease in the EC₅₀, and three cause an increase. Only three of these mutants result in EC₅₀ values similar to WT. The nonfunctional mutant, Pro30Ala, was further probed in oocytes, and the data suggest a role in both expression and function. Examination of the locations of sensitive Pro residues in the receptor and identification of potential interactions with nearby residues reveal how these residues could contribute to the correct functioning of this typical pentameric ligand-gated ion channel.

Structural basis for allosteric transitions of a multidomain pentameric ligand-gated ion channel

Pentameric ligand-gated ion channels (pLGICs) are allosteric receptors that mediate rapid electrochemical signal transduction in the animal nervous system through the opening of an ion pore upon binding of neurotransmitters. Orthologs have been found and characterized in prokaryotes and they display highly similar structure-function relationships to eukaryotic pLGICs; however, they often encode greater architectural diversity involving additional amino-terminal domains (NTDs). Here we report structural, functional, and normal-mode analysis of two conformational states of a multidomain pLGIC, called DeCLIC, from a Desulfofustis deltaproteobacterium, including a periplasmic NTD fused to the conventional ligand-binding domain (LBD). X-ray structure determination revealed an NTD consisting of two jelly-roll domains interacting across each subunit interface. Binding of Ca²⁺ at the LBD subunit interface was associated with a closed transmembrane pore, with resolved monovalent cations intracellular to the hydrophobic gate. Accordingly, DeCLIC-injected oocytes conducted currents only upon depletion of extracellular Ca²⁺; these were insensitive to quaternary ammonium block. Furthermore, DeCLIC crystallized in the absence of Ca²⁺ with a wide-open pore and remodeled periplasmic domains, including increased contacts between the NTD and classic LBD agonist-binding sites. Functional, structural, and dynamical properties of DeCLIC paralleled those of sTeLIC, a pLGIC from another symbiotic prokaryote. Based on these DeCLIC structures, we would reclassify the previous structure of bacterial ELIC (the first high-resolution structure of a pLGIC) as a "locally closed" conformation. Taken together, structures of DeCLIC in multiple conformations illustrate dramatic conformational state transitions and diverse regulatory mechanisms available to ion channels in pLGICs, particularly involving Ca²⁺ modulation and periplasmic NTDs.

The M4 Helix Is Involved in α7 nACh Receptor Function

Nicotinic acetylcholine receptors (nAChR) are the archetypal members of the pentameric ligand-gated ion channel (pLGIC) family, an important class of cell signaling proteins. In all members of this family, each of the five subunits has four transmembrane α-helices (M1-M4) with M2 lining the pore and then M1 and M3, with M4 outermost and adjacent to the membrane lipids. M4 has a variety of roles: its interaction with neighboring M1 and M3 helices is important for receptor assembly, it can a transmit information on the lipid content of the membrane to the gating mechanism, and it may form a vital link to the extracellular domain via the Cys-loop. This study examines the role of M4 receptor residues in the α7 nAChR using site-directed mutagenesis and subsequent expression in Xenopus oocytes. The data indicate that many of the residues in M4 play a role in receptor function, as substitution with Ala can modify functional parameters; 11 of 24 mutants showed a small gain of function (<10-fold decrease in EC₅₀), and 1 (D446A) did not respond to the agonist; it was also not expressed at the cell surface. Removal or addition of aromatic residues had small or no effects. These results demonstrate the α7 nAChR M4 has a role in receptor function, and a structural model suggests possible interactions of some of these residues with their neighbors.

A Novel Glycine Receptor Variant with Startle Disease Affects Syndapin I and Glycinergic Inhibition

Glycine receptors (GlyRs) are the major mediators of fast synaptic inhibition in the adult human spinal cord and brainstem. Hereditary mutations to GlyRs can lead to the rare, but potentially fatal, neuromotor disorder hyperekplexia. Most mutations located in the large intracellular domain (TM3-4 loop) of the GlyRα1 impair surface expression levels of the receptors. The novel GLRA1 mutation P366L, located in the TM3-4 loop, showed normal surface expression but reduced chloride currents, and accelerated whole-cell desensitization observed in whole-cell recordings. At the single-channel level, we observed reduced unitary conductance accompanied by spontaneous opening events in the absence of extracellular glycine. Using peptide microarrays and tandem MS-based analysis methods, we show that the proline-rich stretch surrounding P366 mediates binding to syndapin I, an F-BAR domain protein involved in membrane remodeling. The disruption of the noncanonical Src homology 3 recognition motif by P366L reduces syndapin I binding. These data suggest that the GlyRα1 subunit interacts with intracellular binding partners and may therefore play a role in receptor trafficking or synaptic anchoring, a function thus far only ascribed to the GlyRβ subunit. Hence, the P366L GlyRα1 variant exhibits a unique set of properties that cumulatively affect GlyR functionality and thus might explain the neuropathological mechanism underlying hyperekplexia in the mutant carriers. P366L is the first dominant GLRA1 mutation identified within the GlyRα1 TM3-4 loop that affects GlyR physiology without altering protein expression at the whole-cell and surface levels.SIGNIFICANCE STATEMENT We show that the intracellular domain of the inhibitory glycine receptor α1 subunit contributes to trafficking and synaptic anchoring. A proline-rich stretch in this receptor domain forms a noncanonical recognition motif important for the interaction with syndapin I (PACSIN1). The disruption of this motif, as present in a human patient with hyperekplexia led to impaired syndapin I binding. Functional analysis revealed that the altered proline-rich stretch determines several functional physiological parameters of the ion channel (e.g., faster whole-cell desensitization) reduced unitary conductance and spontaneous opening events. Thus, the proline-rich stretch from the glycine receptor α1 subunit represents a multifunctional intracellular protein motif.

Progress in nicotinic receptor structural biology

Here we begin by briefly reviewing landmark structural studies on the nicotinic acetylcholine receptor. We highlight challenges that had to be overcome to push through resolution barriers, then focus on what has been gleaned in the past few years from crystallographic and single particle cryo-EM studies of different nicotinic receptor subunit assemblies and ligand complexes. We discuss insights into ligand recognition, ion permeation, and allosteric gating. We then highlight some foundational aspects of nicotinic receptor structural biology that remain unresolved and are areas ripe for future exploration. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Simulating ion channel activation mechanisms using swarms of trajectories

Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.

Modulation of the Erwinia ligand-gated ion channel (ELIC) and the 5-HT3 receptor via a common vestibule site

Pentameric ligand-gated ion channels (pLGICs) or Cys-loop receptors are involved in fast synaptic signaling in the nervous system. Allosteric modulators bind to sites that are remote from the neurotransmitter binding site, but modify coupling of ligand binding to channel opening. In this study, we developed nanobodies (single domain antibodies), which are functionally active as allosteric modulators, and solved co-crystal structures of the prokaryote (Erwinia) channel ELIC bound either to a positive or a negative allosteric modulator. The allosteric nanobody binding sites partially overlap with those of small molecule modulators, including a vestibule binding site that is not accessible in some pLGICs. Using mutagenesis, we extrapolate the functional importance of the vestibule binding site to the human 5-HT3 receptor, suggesting a common mechanism of modulation in this protein and ELIC. Thus we identify key elements of allosteric binding sites, and extend drug design possibilities in pLGICs with an accessible vestibule site.

A Review of the Cholinergic System and Therapeutic Approaches to Treat Brain Disorders

Since its identification over a hundred years ago, the neurotransmitter acetylcholine (ACh) has proven to play an essential role in supporting many diverse functions. Some well-characterized functions include: chemical transmission at the neuromuscular junction; autonomic function in the peripheral nervous system; and, sustained attention, sleep/wake regulation, and learning and memory within the central nervous system. Within the brain, major cholinergic projection pathways from the basal forebrain and the brainstem support these centrally mediated processes, and dysregulation of the cholinergic system is implicated in cognitive decline associated with aging and dementias including Alzheimer's disease. ACh exerts its effects by binding to two different membrane-bound receptor classes: (1) G‑protein coupled muscarinic acetylcholine receptors (mAChRs), and (2) ligand-gated nicotinic acetylcholine receptors (nAChRs). These receptor systems are described in detail within this chapter along with discussion on the successes and failures of synthetic ligands designed to selectively target receptor subtypes for treating brain disorders. New molecular approaches and advances in our understanding of the target biology combined with opportunities to re-purpose existing cholinergic drugs for new indications continue to highlight the exciting opportunities for modulating this system for therapeutic purposes.

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.

Membrane directed expression in Escherichia coli of BBA57 and other virulence factors from the Lyme disease agent Borrelia burgdorferi

Membrane-embedded proteins are critical to the establishment, survival and persistence in the host of the Lyme disease bacterium Borrelia burgdorferi (Bb), but to date, there are no solved structures of transmembrane proteins representing these attractive therapeutic targets. All available structures from the genus Borrelia represent proteins expressed without a membrane-targeting signal peptide, thus avoiding conserved pathways that modify, fold and assemble membrane protein complexes. Towards elucidating structure and function of these critical proteins, we directed translocation of eleven expression-optimized Bb virulence factors, including the signal sequence, to the Escherichia coli membrane, of which five, BBA57, HtrA, BB0238, BB0323, and DipA, were expressed with C-terminal His-tags. P66 was also expressed using the PelB signal sequence fused to maltose binding protein. Membrane-associated BBA57 lipoprotein was solubilized by non-ionic and zwitterionic detergents. We show BBA57 translocation to the outer membrane, purification at a level sufficient for structural studies, and evidence for an α-helical multimer. Previous studies showed multiple critical roles of BBA57 in transmission, joint arthritis, carditis, weakening immune responses, and regulating other Bb outer surface proteins. In describing the first purification of membrane-translocated BBA57, this work will support subsequent studies that reveal the precise mechanisms of this important Lyme disease virulence factor.

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.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP₂ plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP₂ and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABA_A receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP₂ also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.

A structural perspective on GABAA receptor pharmacology

GABA_A receptors are pentameric ligand-gated chloride channels of crucial importance for the vertebrate nervous system physiology. They typically modulate the fast inhibitory neurotransmission, and represent the target receptors for major classes of drugs used in the clinic, such as benzodiazepines and general anesthetics. Recent technological progress in structural biology, in particular single-particle cryo-electron microscopy, has led to fundamental advances in understanding the detailed organization and signalling mechanisms of major GABA_A receptor subtypes. This effort culminated with the high-resolution structural analysis of an intact, full-length human heteropentameric receptor, α1β3γ2, in a lipid bilayer and in complex with small molecule ligands including the commonly used benzodiazepines diazepam (Valium) and alprazolam (Xanax). These structures reveal multiple aspects of receptor activation and provide a path for rational design of subunit-specific GABA_A receptor modulators.

The nicotinic acetylcholine receptor: a typical 'allosteric machine'

The concept of allosteric interaction was initially proposed to account for the inhibitory feedback mechanism mediated by bacterial regulatory enzymes. In contrast with the classical mechanism of competitive, steric, interaction between ligands for a common site, allosteric interactions take place between topographically distinct sites and are mediated by a discrete and reversible conformational change of the protein. The concept was soon extended to membrane receptors for neurotransmitters and shown to apply to the signal transduction process which, in the case of the acetylcholine nicotinic receptor (nAChR), links the ACh binding site to the ion channel. Pharmacological effectors, referred to as allosteric modulators, such as Ca²⁺ ions and ivermectin, were discovered that enhance the transduction process when they bind to sites distinct from the orthosteric ACh site and the ion channel. The recent X-ray and electron microscopy structures, at atomic resolution, of the resting and active conformations of several homologues of the nAChR, in combination with atomistic molecular dynamics simulations reveal a stepwise quaternary transition in the transduction process with tertiary changes modifying the boundaries between subunits. These interfaces host orthosteric and allosteric modulatory sites which structural organization changes in the course of the transition. The nAChR appears as a typical allosteric machine. The model emerging from these studies has led to the conception and development of several new pharmacological agents.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Effects of C-Terminal Carboxylation on α-Conotoxin LsIA Interactions with Human α7 Nicotinic Acetylcholine Receptor: Molecular Simulation Studies

α-Conotoxins selectively bind to nicotinic acetylcholine receptors (nAChRs), which are therapeutic targets due to their important role in signaling transmission in excitable cells. A previous experimental study has demonstrated that carboxylation of the C-terminal of α-conotoxin LsIA reduces its potency to inhibit human α7 nAChR relative to naturally amidated LsIA. However, little is known about the contribution of conformational changes in the receptor and interactions, induced by C-terminal amidation/carboxylation of conotoxins, to selective binding to nAChRs, since most conotoxins and some disulfide-rich peptides from other conotoxin subfamilies possess a naturally amidated C-terminal. In this study, we employ homology modeling and molecular dynamics (MD) simulations to propose the determinants for differential interactions between amidated and carboxylated LsIAs with α7 nAChR. Our findings indicate an overall increased number of contacts favored by binding of amidated LsIA versus its carboxylated counterpart. Toxin-receptor pairwise interactions, which may play a role in enhancing the potency of the former, include ARG10-TRP77, LEU141 and CYS17-GLN79 via persistent hydrogen bonds and cation-π interactions, which are weakened in the carboxylated form due to a strong intramolecular salt-bridge formed by ARG10 and carboxylated C-terminus. The binding of amidated LsIA also induces enhanced movements in loop C and the juxtamembrane Cys-loop that are closely associated with receptor function. Additionally, the impacts of binding of LsIA on the overall structure and inter-subunit contacts were examined using inter-residue network analysis, suggesting a clockwise tilting of the α7 C and F loops upon binding to carboxylated LsIA, which is absent for amidated LsIA binding. The predicted molecular mechanism of LsIA binding to the α7 receptor may provide new insights into the important role of the C-terminal in the binding potency of conotoxins at neuronal nAChRs for pharmacological purposes.

Trans- Cis Proline Switches in a Pentameric Ligand-Gated Ion Channel: How They Are Affected by and How They Affect the Biomolecular Environment

Pentameric ligand-gated ion channels (pLGICs) are important neuroreceptors, embedded in neuronal membranes, that mediate fast synaptic transmission. The molecular details of their working mechanisms have still to be fully unravelled due to their complexity and limited structural information available. Here we focus on a potential molecular switch in a prototypical pLGIC, the serotonin-activated 5-HT₃ receptor, consisting of the trans- cis isomerization of a proline at the interface between the extracellular and transmembrane domain. Mutagenesis electrophysiology experiments previously showed that if such isomerization could not take place, the channel would not open, but the hypothetical role of this mechanism as key to channel gating is still debated. We investigate this switch within the receptor with molecular dynamics and enhanced sampling simulations. We analyze how the isomerization free energy landscape is affected by the receptor environment in comparison to simplified models. Moreover, we reveal how the isomerization, in turn, affects the structural and electrostatic properties of the receptor at the extracellular-transmembrane domain interface, e.g., by tuning the ion selectivity filter.