An Ion-Permeable State of the Glycine Receptor Captured by Molecular Dynamics

Mechanism of hydrophobic gating in the acetylcholine receptor channel pore

Neuromuscular acetylcholine receptors (AChRs) are hetero-pentameric, ligand-gated ion channels. The binding of the neurotransmitter acetylcholine (ACh) to two target sites promotes a global conformational change of the receptor that opens the channel and allows ion conduction through the channel pore. Here, by measuring free-energy changes from single-channel current recordings and using molecular dynamics simulations, we elucidate how a constricted hydrophobic region acts as a "gate" to regulate the channel opening in the pore of AChRs. Mutations of gate residues, including those implicated in congenital myasthenia syndrome, lower the permeation barrier of the channel substantially and increase the unliganded gating equilibrium constant (constitutive channel openings). Correlations between hydrophobicity and the observed free-energy changes, supported by calculations of water densities in the wild-type versus mutant channel pores, provide evidence for hydrophobic wetting-dewetting transition at the gate. The analysis of a coupled interaction network provides insight into the molecular mechanism of closed- versus open-state conformational changes at the gate. Studies of the transition state by "phi"(φ)-value analysis indicate that agonist binding serves to stabilize both the transition and the open state. Intersubunit interaction energy measurements and molecular dynamics simulations suggest that channel opening involves tilting of the pore-lining M2 helices, asymmetric outward rotation of amino acid side chains, and wetting transition of the gate region that lowers the barrier to ion permeation and stabilizes the channel open conformation. Our work provides new insight into the hydrophobic gate opening and shows why the gate mutations result in constitutive AChR channel activity.

Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel

Lipid nanodiscs have become a standard tool for studying membrane proteins, including using single particle cryo-electron microscopy (cryo-EM). We find that reconstituting the pentameric ligand-gated ion channel (pLGIC), Erwinia ligand-gated ion channel (ELIC), in different nanodiscs produces distinct structures by cryo-EM. The effect of the nanodisc on ELIC structure extends to the extracellular domain and agonist binding site. Additionally, molecular dynamic simulations indicate that nanodiscs of different size impact ELIC structure and that the nanodisc scaffold directly interacts with ELIC. These findings suggest that the nanodisc plays a crucial role in determining the structure of pLGICs, and that reconstitution of ion channels in larger nanodiscs may better approximate a lipid membrane environment.

Rational Prediction of PROTAC-Compatible Protein-Protein Interfaces by Molecular Docking

Proteolysis targeting chimeras (PROTACs) are heterobifunctional ligands that mediate the interaction between a protein target and an E3 ligase, resulting in a ternary complex, whose interaction with the ubiquitination machinery leads to target degradation. This technology is emerging as an exciting new avenue for therapeutic development, with several PROTACs currently undergoing clinical trials targeting cancer. Here, we describe a general and computationally efficient methodology combining restraint-based docking, energy-based rescoring, and a filter based on the minimal solvent-accessible surface distance to produce PROTAC-compatible PPIs suitable for when there is no a priori known PROTAC ligand. In a benchmark employing a manually curated data set of 13 ternary complex crystals, we achieved an accuracy of 92% when starting from bound structures and 77% when starting from unbound structures, respectively. Our method only requires that the ligand-bound structures of the monomeric forms of the E3 ligase and target proteins be given to run, making it general, accurate, and highly efficient, with the ability to impact early-stage PROTAC-based drug design campaigns where no structural information about the ternary complex structure is available.

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.

Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

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

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

When is a hydrophobic gate not a hydrophobic gate?

The flux of ions through a channel is most commonly regulated by changes that result in steric occlusion of its pore. However, ion permeation can also be prevented by formation of a desolvation barrier created by hydrophobic residues that line the pore. As a result of relatively minor structural changes, confined hydrophobic regions in channels may undergo transitions between wet and dry states to gate the pore closed without physical constriction of the permeation pathway. This concept is referred to as hydrophobic gating, and many examples of this process have been demonstrated. However, the term is also now being used in a much broader context that often deviates from its original meaning. In this Viewpoint, we explore the formal definition of a hydrophobic gate, discuss examples of this process compared with other gating mechanisms that simply exploit hydrophobic residues and/or lipids in steric closure of the pore, and describe the best practice for identification of a hydrophobic gate.

Lateral fenestrations in the extracellular domain of the glycine receptor contribute to the main chloride permeation pathway

Glycine receptors (GlyRs) are ligand-gated ion channels mediating signal transduction at chemical synapses. Since the early patch-clamp electrophysiology studies, the details of the ion permeation mechanism have remained elusive. Here, we combine molecular dynamics simulations of a zebrafish GlyR-α1 model devoid of the intracellular domain with mutagenesis and single-channel electrophysiology of the full-length human GlyR-α1. We show that lateral fenestrations between subunits in the extracellular domain provide the main translocation pathway for chloride ions to enter/exit a central water-filled vestibule at the entrance of the transmembrane channel. In addition, we provide evidence that these fenestrations are at the origin of current rectification in known anomalous mutants and design de novo two inward-rectifying channels by introducing mutations within them. These results demonstrate the central role of lateral fenestrations on synaptic neurotransmission.

Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations

Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.

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.

Partial agonists go molecular

Single-channel analysis previously revealed a key role for a short-lived 'flipped' state during glycine receptor activation by partial agonists. Structures solved by Yu and colleagues now reveal a surprising mechanism involving a partially activated agonist-bound closed state that is too long-lived to be considered the flipped state.

On the Functional Annotation of Open-Channel Structures in the Glycine Receptor

The glycine receptor (GlyR) is by far the best-characterized pentameric ligand-gated ion channel, with several high-resolution structures from X-ray crystallography, cryoelectron microscopy (cryo-EM), and modeling. Nonetheless, the significance of the currently available open-pore conformations is debated due to their diversity in the pore geometry. Here, we discuss the physiological significance of existing models of the GlyR active state based on conductance and selectivity measurements by computational electrophysiology. The results support the conclusion that the original cryo-EM reconstruction of the active state obtained in detergents as well as its subsequent refinement by molecular dynamics simulations are likely to be non-physiological as they feature artificially dilated ion pores. In addition, the calculations indicate that a physiologically relevant open pore should be constricted within a radius of 2.5 and 2.8 Å, which is consistent with previous modeling, electrophysiology measurements, and the most recent cryo-EM structures obtained in a native lipid membrane environment.

Mechanism of gating and partial agonist action in the glycine receptor

Ligand-gated ion channels mediate signal transduction at chemical synapses and transition between resting, open, and desensitized states in response to neurotransmitter binding. Neurotransmitters that produce maximum open channel probabilities (Po) are full agonists, whereas those that yield lower than maximum Po are partial agonists. Cys-loop receptors are an important class of neurotransmitter receptors, yet a structure-based understanding of the mechanism of partial agonist action has proven elusive. Here, we study the glycine receptor with the full agonist glycine and the partial agonists taurine and γ-amino butyric acid (GABA). We use electrophysiology to show how partial agonists populate agonist-bound, closed channel states and cryo-EM reconstructions to illuminate the structures of intermediate, pre-open states, providing insights into previously unseen conformational states along the receptor reaction pathway. We further correlate agonist-induced conformational changes to Po across members of the receptor family, providing a hypothetical mechanism for partial and full agonist action at Cys-loop receptors.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

The desensitization pathway of GABAA receptors, one subunit at a time

GABAA receptors mediate most inhibitory synaptic transmission in the brain of vertebrates. Following GABA binding and fast activation, these receptors undergo a slower desensitization, the conformational pathway of which remains largely elusive. To explore the mechanism of desensitization, we used concatemeric α1β2γ2 GABAA receptors to selectively introduce gain-of-desensitization mutations one subunit at a time. A library of twenty-six mutant combinations was generated and their bi-exponential macroscopic desensitization rates measured. Introducing mutations at the different subunits shows a strongly asymmetric pattern with a key contribution of the γ2 subunit, and combining mutations results in marked synergistic effects indicating a non-concerted mechanism. Kinetic modelling indeed suggests a pathway where subunits move independently, the desensitization of two subunits being required to occlude the pore. Our work thus hints towards a very diverse and labile conformational landscape during desensitization, with potential implications in physiology and pharmacology.

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.

The general anesthetic etomidate and fenamate mefenamic acid oppositely affect GABAAR and GlyR: a structural explanation

GABA and glycine act as inhibitory neurotransmitters in the CNS. Inhibitory neurotransmission is mediated via activation of ionotropic GABA_A and glycine receptors. We used a modeling approach to explain the opposite effects of the general anesthetic etomidate (ETM) and fenamate mefenamic acid (MFA) on GABA- and glycine-activated currents recorded in isolated cerebellar Purkinje cells and hippocampal pyramidal neurons, respectively. These drugs potentiated GABA_ARs but blocked GlyRs. We built a homology model of α1β GlyR based on the cryo-EM structure of open α1 GlyR, used the α1β3γ2 GABA_AR structure from the PDB, and applied Monte-Carlo energy minimization to optimize models of receptors and ligand-receptor complexes. In silico docking suggests that ETM/MFA bind at the transmembrane β( +)/α( -) intersubunit interface in GABA_AR. Our models predict that the bulky side chain of the highly conserved Arg^19' residue at the plus interface side wedges the interface and maintains the conducting receptor state. We hypothesized that MFA/ETM binding at the β( +)/α( -) interface leads to prolongation of receptor life-time in the open state. Having analyzed different GABA_AR and GlyR structures available in the PDB, we found that mutual arrangement of the Arg^19' and Gln^-26' side chains at the plus and minus interface sides, respectively, plays an important role when the receptor switches from the open to closed state. We show that this process is accompanied by narrowing of the intersubunit interfaces, leading to extrusion of the Arg^19' side chain from the interface. Our models allow us to explain the lack of GlyR potentiation in our electrophysiological experiments.

Mechanisms of activation and desensitization of full-length glycine receptor in lipid nanodiscs

Glycinergic synapses play a central role in motor control and pain processing in the central nervous system. Glycine receptors (GlyRs) are key players in mediating fast inhibitory neurotransmission at these synapses. While previous high-resolution structures have provided insights into the molecular architecture of GlyR, several mechanistic questions pertaining to channel function are still unanswered. Here, we present Cryo-EM structures of the full-length GlyR protein complex reconstituted into lipid nanodiscs that are captured in the unliganded (closed), glycine-bound (open and desensitized), and allosteric modulator-bound conformations. A comparison of these states reveals global conformational changes underlying GlyR channel gating and modulation. The functional state assignments were validated by molecular dynamics simulations, and the observed permeation events are in agreement with the anion selectivity and conductance of GlyR. These studies provide the structural basis for gating, ion selectivity, and single-channel conductance properties of GlyR in a lipid environment.

Glycinergic synapses play a central role in motor control and pain processing in the central nervous system. Here, authors present cryo-EM structures of the full-length glycine receptors (GlyRs) reconstituted into lipid nanodiscs in the unliganded, glycine-bound and allosteric modulator-bound conformations and reveal global conformational changes underlying GlyR channel gating and modulation.

Comment on "On the Functional Annotation of Open-Channel Structures in the Glycine Receptor"

Recently, we reported the simulation of a stable open state of the glycine receptor. Central to the stability of the simulations was the behavior of the highly conserved leucine residues at the 9' gate, which were found to rotate into adjacent pockets, thus providing a structural rationale for decades of biochemical observations. In contrast, a previously reported model from Cerdan et al. (2018) resembled a more collapsed state. However, in support of their model, they draw attention to the agreement between calculated and experimental conductance measurements and argue that our model tends to overestimate ion flow. Here, we argue that there are many pitfalls with this approach and that the apparent agreement most likely reflects a fortuitous cancellation of errors. The computed values are highly sensitive to very small changes in model parameters. It is also likely that polarization effects will be very significant, and these have not yet been considered.

The Glycine Receptor Allosteric Ligands Library (GRALL)

MOTIVATION

Glycine receptors (GlyRs) mediate fast inhibitory neurotransmission in the brain and have been recognized as key pharmacological targets for pain. A large number of chemically diverse compounds that are able to modulate GlyR function both positively and negatively have been reported, which provides useful information for the development of pharmacological strategies and models for the allosteric modulation of these ion channels.

RESULTS

Based on existing literature, we have collected 218 unique chemical entities with documented modulatory activities at homomeric GlyR-α1 and -α3 and built a database named GRALL. This collection includes agonists, antagonists, positive and negative allosteric modulators and a number of experimentally inactive compounds. Most importantly, for a large fraction of them a structural annotation based on their putative binding site on the receptor is provided. This type of annotation, which is currently missing in other drug banks, along with the availability of cooperativity factors from radioligand displacement experiments are expected to improve the predictivity of in silico methodologies for allosteric drug discovery and boost the development of conformation-based pharmacological approaches.

AVAILABILITY AND IMPLEMENTATION

The GRALL library is distributed as a web-accessible database at the following link: https://ifm.chimie.unistra.fr/grall. For each molecular entry, it provides information on the chemical structure, the ligand-binding site, the direction of modulation, the potency, the 3D molecular structure and quantum-mechanical charges as determined by our in-house pipeline.

CONTACT

mcecchini@unistra.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The startle disease mutation α1S270T predicts shortening of glycinergic synaptic currents

KEY POINTS

Loss-of-function mutations in proteins found at glycinergic synapses, most commonly in the α1 subunit of the glycine receptor (GlyR), cause the startle disease/hyperekplexia channelopathy in man. It was recently proposed that the receptors responsible are presynaptic homomeric GlyRs, rather than postsynaptic heteromeric GlyRs (which mediate glycinergic synaptic transmission), because heteromeric GlyRs are less affected by many startle mutations than homomers. We examined the α1 startle mutation S270T, at the extracellular end of the M2 transmembrane helix. Recombinant heteromeric GlyRs were less impaired than homomers by this mutation when we measured their response to equilibrium applications of glycine. However, currents elicited by synaptic-like millisecond applications of glycine to outside-out patches were much shorter (7- to 10-fold) in all mutant receptors, both homomeric and heteromeric. Thus, the synaptic function of heteromeric receptors is likely to be impaired by the mutation.

ABSTRACT

Human startle disease is caused by mutations in glycine receptor (GlyR) subunits or in other proteins associated with glycinergic synapses. Many startle mutations are known, but it is hard to correlate the degree of impairment at molecular level with the severity of symptoms in patients. It was recently proposed that the disease is caused by disruption in the function of presynaptic homomeric GlyRs (rather than postsynaptic heteromeric GlyRs), because homomeric GlyRs are more sensitive to loss-of-function mutations than heteromers. Our patch-clamp recordings from heterologously expressed GlyRs characterised in detail the functional consequences of the α1S270T startle mutation, which is located at the extracellular end of the pore lining M2 transmembrane segment (18'). This mutation profoundly decreased the maximum single-channel open probability of homomeric GlyRs (to 0.16; cf. 0.99 for wild type) but reduced only marginally that of heteromeric GlyRs (0.96; cf. 0.99 for wild type). However, both heteromeric and homomeric mutant GlyRs became less sensitive to the neurotransmitter glycine. Responses evoked by brief, quasi-synaptic pulses of glycine onto outside-out patches were impaired in mutant receptors, as deactivation was approximately 10- and 7-fold faster for homomeric and heteromeric GlyRs, respectively. Our data suggest that the α1S270T mutation is likely to affect the opening step in GlyR activation. The faster decay of synaptic currents mediated by mutant heteromeric GlyRs is expected to reduce charge transfer at the synapse, despite the high equilibrium open probability of these mutant channels.

Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor Channel

Ion channel proteins form water-filled nanoscale pores within lipid bilayers, and their properties are dependent on the complex behavior of water in a nanoconfined environment. Using a simplified model of the pore of the 5-HT3 receptor (5HT3R) which restrains the backbone structure to that of the parent channel protein from which it is derived, we compare additive with polarizable models in describing the behavior of water in nanopores. Molecular dynamics simulations were performed with four conformations of the channel: two closed state structures, an intermediate state, and an open state, each embedded in a phosphatidylcholine bilayer. Water density profiles revealed that for all water models, the closed and intermediate states exhibited strong dewetting within the central hydrophobic gate region of the pore. However, the open state conformation exhibited varying degrees of hydration, ranging from partial wetting for the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. Water dipole moments calculated using polarizable force fields also revealed that water molecules remaining within dewetted sections of the pore resemble gas phase water. Free energy profiles for Na+ and for Cl- ions within the open state pore revealed more rugged energy landscapes using polarizable force fields, and the hydration number profiles of these ions were also sensitive to induced polarization resulting in a substantive reduction of the number of waters within the first hydration shell of Cl- while it permeates the pore. These results demonstrate that induced polarization can influence the complex behavior of water and ions within nanoscale pores and provides important new insights into their chemical properties.

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

A Refined Open State of the Glycine Receptor Obtained via Molecular Dynamics Simulations

Pentameric ligand-gated ion channels are key players in mediating fast neurotransmission. Glycine receptors are chloride-selective members of this receptor family that mediate inhibitory synaptic transmission and are implicated in neurological disorders including autism and hyperekplexia. They have been structurally characterized by both X-ray crystallography and cryoelectron microscopy (cryo-EM) studies, with the latter giving rise to what was proposed as a possible open state. However, recent work has questioned the physiological relevance of this open state structure, since it rapidly collapses in molecular dynamics simulations. Here, we show that the collapse can be avoided by a careful equilibration protocol that reconciles the more problematic regions of the original density map and gives a stable open state that shows frequent selective chloride permeation. The protocol developed in this work provides a means to refine open-like structures of the whole pentameric ligand-gated ion channel superfamily and reconciles the previous issues with the cryo-EM structure.

Agonist Selectivity and Ion Permeation in the α3β4 Ganglionic Nicotinic Receptor

Nicotinic acetylcholine receptors are pentameric ion channels that mediate fast chemical neurotransmission. The α3β4 nicotinic receptor subtype forms the principal relay between the central and peripheral nervous systems in the autonomic ganglia. This receptor is also expressed focally in brain areas that affect reward circuits and addiction. Here, we present structures of the α3β4 nicotinic receptor in lipidic and detergent environments, using functional reconstitution to define lipids appropriate for structural analysis. The structures of the receptor in complex with nicotine, as well as the α3β4-selective ligand AT-1001, complemented by molecular dynamics, suggest principles of agonist selectivity. The structures further reveal much of the architecture of the intracellular domain, where mutagenesis experiments and simulations define residues governing ion conductance.

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'.