1
|
Li Z, Chan KC, Nickels JD, Cheng X. Molecular Dynamics Refinement of Open State Serotonin 5-HT 3A Receptor Structures. J Chem Inf Model 2023; 63:1196-1207. [PMID: 36757760 DOI: 10.1021/acs.jcim.2c01441] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT3 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-HT3 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-HT3 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-HT3 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-HT3 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-HT3 receptor.
Collapse
Affiliation(s)
- Zoe Li
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio 43210, United States
| | - Kevin C Chan
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio 43210, United States
| | - Jonathan D Nickels
- Department of Chemical and Environmental Engineering, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Xiaolin Cheng
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio 43210, United States.,Translational Data Analytics Institute (TDAI) at The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
2
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
3
|
Westra RL. Resonance-driven ion transport and selectivity in prokaryotic ion channels. Phys Rev E 2019; 100:062410. [PMID: 31962411 DOI: 10.1103/physreve.100.062410] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Indexed: 06/10/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ronald L Westra
- Department of Data Science and Knowledge Engineering, Faculty of Science and Engineering, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
4
|
Oliveira ASF, Shoemark DK, Campello HR, Wonnacott S, Gallagher T, Sessions RB, Mulholland AJ. Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations. Structure 2019; 27:1171-1183.e3. [PMID: 31130483 DOI: 10.1016/j.str.2019.04.008] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 01/28/2019] [Accepted: 04/10/2019] [Indexed: 02/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- School of Biochemistry, University of Bristol, Bristol BS8 1DT, UK; Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Hugo Rego Campello
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Susan Wonnacott
- Department of Biology and Biochemistry, University of Bath, Bath BA2 7AY, UK
| | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| |
Collapse
|
5
|
Crnjar A, Comitani F, Melis C, Molteni C. Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution. Interface Focus 2019; 9:20180067. [PMID: 31065340 PMCID: PMC6501341 DOI: 10.1098/rsfs.2018.0067] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/26/2019] [Indexed: 12/21/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Alessandro Crnjar
- King’s College London, Department of Physics, Strand, London WC2R 2LS, UK
| | - Federico Comitani
- Program in Genetics and Genome Biology, The Hospital for Sick Children, Toronto, Ontario, Canada
| | - Claudio Melis
- Universitá degli Studi di Cagliari, Complesso Universitario di Monserrato, Dipartimento di Fisica, S.P. Monserrato-Sestu Km 0,700, Monserrato (CA) 09042, Italy
| | - Carla Molteni
- King’s College London, Department of Physics, Strand, London WC2R 2LS, UK
| |
Collapse
|
6
|
Changeux JP. The nicotinic acetylcholine receptor: a typical 'allosteric machine'. Philos Trans R Soc Lond B Biol Sci 2019; 373:rstb.2017.0174. [PMID: 29735728 DOI: 10.1098/rstb.2017.0174] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/06/2017] [Indexed: 12/26/2022] Open
Abstract
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 Ca2+ 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'.
Collapse
Affiliation(s)
- Jean-Pierre Changeux
- CNRS UMR 3571, Institut Pasteur, Paris 75724, France .,Communications Cellulaires, Collège de France, Paris 75005, France
| |
Collapse
|
7
|
Competitive docking model for prediction of the human nicotinic acetylcholine receptor α7 binding of tobacco constituents. Oncotarget 2018; 9:16899-16916. [PMID: 29682193 PMCID: PMC5908294 DOI: 10.18632/oncotarget.24458] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2017] [Accepted: 02/01/2018] [Indexed: 12/21/2022] Open
Abstract
The detrimental health effects associated with tobacco use constitute a major public health concern. The addiction associated with nicotine found in tobacco products has led to difficulty in quitting among users. Nicotinic acetylcholine receptors (nAChRs) are the targets of nicotine and are responsible for addiction to tobacco products. However, it is unknown if the other >8000 tobacco constituents are addictive. Since it is time-consuming and costly to experimentally assess addictive potential of such larger number of chemicals, computationally predicting human nAChRs binding is important for in silico evaluation of addiction potential of tobacco constituents and needs structures of human nAChRs. Therefore, we constructed three-dimensional structures of the ligand binding domain of human nAChR α7 subtype and then developed a predictive model based on the constructed structures to predict human nAChR α7 binding activity of tobacco constituents. The predictive model correctly predicted 11 out of 12 test compounds to be binders of nAChR α7. The model is a useful tool for high-throughput screening of potential addictive tobacco constituents. These results could inform regulatory science research by providing a new validated predictive tool using cutting-edge computational methodology to high-throughput screen tobacco additives and constituents for their binding interaction with the human α7 nicotinic receptor. The tool represents a prediction model capable of screening thousands of chemicals found in tobacco products for addiction potential, which improves the understanding of the potential effects of additives.
Collapse
|
8
|
Trick JL, Chelvaniththilan S, Klesse G, Aryal P, Wallace EJ, Tucker SJ, Sansom MSP. Functional Annotation of Ion Channel Structures by Molecular Simulation. Structure 2016; 24:2207-2216. [PMID: 27866853 PMCID: PMC5145807 DOI: 10.1016/j.str.2016.10.005] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2016] [Revised: 08/12/2016] [Accepted: 10/12/2016] [Indexed: 01/30/2023]
Abstract
Ion channels play key roles in cell membranes, and recent advances are yielding an increasing number of structures. However, their functional relevance is often unclear and better tools are required for their functional annotation. In sub-nanometer pores such as ion channels, hydrophobic gating has been shown to promote dewetting to produce a functionally closed (i.e., non-conductive) state. Using the serotonin receptor (5-HT3R) structure as an example, we demonstrate the use of molecular dynamics to aid the functional annotation of channel structures via simulation of the behavior of water within the pore. Three increasingly complex simulation analyses are described: water equilibrium densities; single-ion free-energy profiles; and computational electrophysiology. All three approaches correctly predict the 5-HT3R crystal structure to represent a functionally closed (i.e., non-conductive) state. We also illustrate the application of water equilibrium density simulations to annotate different conformational states of a glycine receptor.
Collapse
Affiliation(s)
- Jemma L Trick
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK
| | - Sivapalan Chelvaniththilan
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK; Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Gianni Klesse
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK; Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Prafulla Aryal
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK; Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK; OXION Initiative in Ion Channels and Disease, University of Oxford, Oxford, UK
| | | | - Stephen J Tucker
- Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK; OXION Initiative in Ion Channels and Disease, University of Oxford, Oxford, UK
| | - Mark S P Sansom
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK; OXION Initiative in Ion Channels and Disease, University of Oxford, Oxford, UK.
| |
Collapse
|
9
|
Yonkunas M, Kurnikova M. The Hydrophobic Effect Contributes to the Closed State of a Simplified Ion Channel through a Conserved Hydrophobic Patch at the Pore-Helix Crossing. Front Pharmacol 2015; 6:284. [PMID: 26640439 PMCID: PMC4661268 DOI: 10.3389/fphar.2015.00284] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2015] [Accepted: 11/10/2015] [Indexed: 11/13/2022] Open
Abstract
Ion selectivity-filter structures are strikingly similar throughout the large family of K++ channels and other p-loop-like receptors (i.e., glutamate receptors). At the same time, the triggers for opening these channels, or gating, are diverse. Two questions that remain unanswered regarding these channels are: (1) what force(s) stabilize the closed non-conducting channel-pore conformation? And (2) what is the free energy associated with transitioning from a closed (non-conducting) to an open (conducting) channel-pore conformation? The effects of charge and hydrophobicity on the conformational states of a model tetrameric biological ion channel are shown utilizing the amino acid sequence from the K+ channel KcsA as the model “channel”. Its widely conserved hydrophobic bundle crossing located adjacent to the lipid head-groups at the intracellular side of the membrane was calculated to have a 5 kcal/mol free energy difference between modeled open and closed conformations. Simulated mutants of amino acids within the hydrophobic region significantly contribute to the size of this difference. Specifically for KcsA, these residues are part of the pH sensor important for channel gating and our results are in agreement with published electrophysiology data. Our simulations support the idea that the hydrophobic effect contributes significantly to the stability of the closed conformation in tetrameric ion channels with a hydrophobic bundle crossing positioned in proximity to the lipid head groups of the biological membrane.
Collapse
Affiliation(s)
- Michael Yonkunas
- Department of Chemistry, Carnegie Mellon University, Pittsburgh PA, USA
| | - Maria Kurnikova
- Department of Chemistry, Carnegie Mellon University, Pittsburgh PA, USA
| |
Collapse
|
10
|
Grazioso G, Sgrignani J, Capelli R, Matera C, Dallanoce C, De Amici M, Cavalli A. Allosteric Modulation of Alpha7 Nicotinic Receptors: Mechanistic Insight through Metadynamics and Essential Dynamics. J Chem Inf Model 2015; 55:2528-39. [PMID: 26569022 DOI: 10.1021/acs.jcim.5b00459] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Increasing attention has recently been devoted to allosteric modulators, as they can provide inherent advantages over classic receptor agonists. In the field of nicotinic receptors (nAChRs), the main advantage is that allosteric modulators can trigger pharmacological responses, limiting receptor desensitization. Most of the known allosteric ligands are "positive allosteric modulators" (PAMs), which increase both sensitivity to receptor agonists and current amplitude. Intriguingly, some allosteric modulators are also able to activate the α7 receptor (α7-nAChR) even in the absence of orthosteric agonists. These compounds have been named "ago-allosteric modulators" and GAT107 has been studied in depth because of its unique mechanism of action. We here investigate by molecular dynamics simulations, metadynamics, and essential dynamics the activation mechanism of α7-nAChR, in the presence of different nicotinic modulators. We determine the free energy profiles associated with the closed-to-open motion of the loop C, and we highlight mechanistic differences observed in the presence of different modulators. In particular, we demonstrate that GAT107 triggers conformational motions and cross-talk similar to those observed when the α7-nACh receptor is in complex with both an agonist and an allosteric modulator.
Collapse
Affiliation(s)
- Giovanni Grazioso
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Farmaceutica "Pietro Pratesi", Università degli Studi di Milano , Via L. Mangiagalli 25, 20133 Milan, Italy
| | - Jacopo Sgrignani
- Institute of Research in Biomedicine (IRB) , Via Vincenzo Vela 6, 6500 Bellinzona, Switzerland
| | - Romina Capelli
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Farmaceutica "Pietro Pratesi", Università degli Studi di Milano , Via L. Mangiagalli 25, 20133 Milan, Italy
| | - Carlo Matera
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Farmaceutica "Pietro Pratesi", Università degli Studi di Milano , Via L. Mangiagalli 25, 20133 Milan, Italy
| | - Clelia Dallanoce
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Farmaceutica "Pietro Pratesi", Università degli Studi di Milano , Via L. Mangiagalli 25, 20133 Milan, Italy
| | - Marco De Amici
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Farmaceutica "Pietro Pratesi", Università degli Studi di Milano , Via L. Mangiagalli 25, 20133 Milan, Italy
| | - Andrea Cavalli
- Drug Discovery and Development-Computation, Istituto Italiano di Tecnologia , Via Morego 30, 16163 Genoa, Italy.,Department of Pharmacy and Biotecnology, University of Bologna , Via Belmeloro 6, 40126 Bologna, Italy
| |
Collapse
|
11
|
Yoluk Ö, Lindahl E, Andersson M. Conformational gating dynamics in the GluCl anion-selective chloride channel. ACS Chem Neurosci 2015; 6:1459-67. [PMID: 25992588 DOI: 10.1021/acschemneuro.5b00111] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Cys-loop receptors are central to propagation of signals in the nervous system. The gating of the membrane-spanning pore is triggered by structural rearrangements in the agonist-binding site, located some 50 Å away from the pore. A sequential conformational change, propagating from the ligand-binding site to the pore, has been proposed to govern gating in all Cys-loop receptors. Here, we identify structural and dynamic components of the conformational gating in the eukaryotic glutamate-gated chloride channel (GluCl) by means of molecular dynamics (MD) simulations with and without the l-glutamate agonist bound. A significant increase in pore opening and accompanying hydration is observed in the presence of glutamate. Potential of mean force calculations reveal that the barrier for ion passage drops from 15 kcal/mol to 5-10 kcal/mol with the agonist bound. This appears to be explained by agonist binding that leads to significant changes in the intersubunit hydrogen-bonding pattern, which induce a slight tilt of the extracellular domain relative to the transmembrane domain in the simulations. This rearrangement is subtle, but correspond to the direction of the quaternary twist observed as a key difference between open and closed X-ray structures. While the full reversible gating is still a much slower process, the observed structural dynamics sheds new light on the early stages of how the agonist influences the extracellular domain, how the extracellular domain interacts with the transmembrane domain, and how changes in the transmembrane domain alter the free energy of ion passage.
Collapse
Affiliation(s)
- Özge Yoluk
- Science for Life Laboratory, Stockholm and Uppsala, 171 21 Stockholm, Sweden
- Theoretical and Computational Biophysics, Department of Theoretical Physics, KTH Royal Institute of Technology, 106 91 Stockholm, Sweden
| | - Erik Lindahl
- Science for Life Laboratory, Stockholm and Uppsala, 171 21 Stockholm, Sweden
- Theoretical and Computational Biophysics, Department of Theoretical Physics, KTH Royal Institute of Technology, 106 91 Stockholm, Sweden
- Department of Biochemistry and Biophysics, Center for Biomembrane Research, Stockholm University, 114 18 Stockholm, Sweden
| | - Magnus Andersson
- Science for Life Laboratory, Stockholm and Uppsala, 171 21 Stockholm, Sweden
- Theoretical and Computational Biophysics, Department of Theoretical Physics, KTH Royal Institute of Technology, 106 91 Stockholm, Sweden
| |
Collapse
|
12
|
Neale C, Chakrabarti N, Pomorski P, Pai EF, Pomès R. Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA. PLoS Comput Biol 2015; 11:e1004303. [PMID: 26181442 PMCID: PMC4504495 DOI: 10.1371/journal.pcbi.1004303] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2014] [Accepted: 04/28/2015] [Indexed: 12/17/2022] Open
Abstract
Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg2+ into cells, we examine early events in the relaxation of Mg2+ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg2+ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction. This study shows how rapid wetting/dewetting transitions in the pores of ion channels participate in the control of biological ion permeation. Ion channels catalyze ionic permeation across non-polar membranes via water-filled pores. However, non-polar stretches or hydrophobic bottlenecks are present in the pores of many ion channels. To clarify the relationship between channel regulation, pore hydration, and ion permeation, we examine how the slow relaxation of magnesium channel CorA from its closed state towards its open state modulates wetting of its hydrophobic bottleneck. Results provide a quantitative description of wetting and dewetting probabilities and kinetics and a quantitative relationship between the extent of pore hydration and the energetics of ion permeation, consistent with a mechanism of hydrophobic gating.
Collapse
Affiliation(s)
- Chris Neale
- Molecular Structure and Function, The Hospital for Sick Children, Toronto, Ontario, Canada
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
| | - Nilmadhab Chakrabarti
- Molecular Structure and Function, The Hospital for Sick Children, Toronto, Ontario, Canada
| | - Pawel Pomorski
- Shared Hierarchical Academic Research Computing Network, Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario, Canada
| | - Emil F. Pai
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
- Ontario Cancer Institute/Princess Margaret Cancer Centre, Campbell Family Institute for Cancer Research, Toronto, Ontario, Canada
| | - Régis Pomès
- Molecular Structure and Function, The Hospital for Sick Children, Toronto, Ontario, Canada
- Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
- * E-mail:
| |
Collapse
|
13
|
Masood TB, Sandhya S, Chandra N, Natarajan V. CHEXVIS: a tool for molecular channel extraction and visualization. BMC Bioinformatics 2015; 16:119. [PMID: 25888118 PMCID: PMC4411761 DOI: 10.1186/s12859-015-0545-9] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2014] [Accepted: 03/20/2015] [Indexed: 12/11/2022] Open
Abstract
Background Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties. Results We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, ChExVis. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that ChExVis performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how ChExVis can be used to study channels, extract their properties and gain insights into molecular function. Conclusion ChExVis supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The ChExVis web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0545-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Talha Bin Masood
- Department of Computer Science and Automation, Indian Institute of Science, Bangalore, 560012, India.
| | - Sankaran Sandhya
- Department of Biochemistry, Indian Institute of Science, Bangalore, 560012, India.
| | - Nagasuma Chandra
- Department of Biochemistry, Indian Institute of Science, Bangalore, 560012, India.
| | - Vijay Natarajan
- Department of Computer Science and Automation, Indian Institute of Science, Bangalore, 560012, India. .,Supercomputer Education and Research Centre, Indian Institute of Science, Bangalore, 560012, India.
| |
Collapse
|
14
|
Changeux JP. Protein dynamics and the allosteric transitions of pentameric receptor channels. Biophys Rev 2014; 6:311-321. [PMID: 25505495 PMCID: PMC4256460 DOI: 10.1007/s12551-014-0149-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2014] [Accepted: 11/13/2014] [Indexed: 01/30/2023] Open
Abstract
The recent application of molecular dynamics (MD) methodology to investigate the allosteric transitions of the acetylcholine receptor and its prokaryotic and eukaryotic pentameric homologs has yielded new insights into the mechanisms of signal transduction by these receptors. Combined with available data on X-ray structures, MD techniques enable description of the dynamics of the conformational change at the atomic level, intra-molecular propagation of this signal transduction mechanism as a concerted stepwise process at physiological timescales and the control of this process by allosteric modulators, thereby offering new perspectives for drug design.
Collapse
Affiliation(s)
- Jean-Pierre Changeux
- UMR 3571 CNRS, Institut Pasteur, 25 rue du Docteur Roux, 75015 Paris, France
- Collège de France, 11 Place Marcelin Berthelot, 75005 Paris, France
- Kavli Brain-Mind Institute University of California, San Diego, CA USA
| |
Collapse
|
15
|
Thompson AJ, Alqazzaz M, Price KL, Weston DA, Lummis SCR. Phenylalanine in the pore of the Erwinia ligand-gated ion channel modulates picrotoxinin potency but not receptor function. Biochemistry 2014; 53:6183-8. [PMID: 25238029 PMCID: PMC4312132 DOI: 10.1021/bi5008035] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
![]()
The Erwinia ligand-gated ion channel (ELIC) is
a bacterial homologue of eukaryotic Cys-loop ligand-gated ion channels.
This protein has the potential to be a useful model for Cys-loop receptors
but is unusual in that it has an aromatic residue (Phe) facing into
the pore, leading to some predictions that this protein is incapable
of ion flux. Subsequent studies have shown this is not the case, so
here we probe the role of this residue by examining the function of
the ELIC in cases in which the Phe has been substituted with a range
of alternative amino acids, expressed in Xenopus oocytes
and functionally examined. Most of the mutations have little effect
on the GABA EC50, but the potency of the weak pore-blocking
antagonist picrotoxinin at F16′A-, F16′D-, F16′S-,
and F16′T-containing receptors was increased to levels comparable
with those of Cys-loop receptors, suggesting that this antagonist
can enter the pore only when residue 16′ is small. T6′S
has no effect on picrotoxinin potency when expressed alone but abolishes
the increased potency when combined with F16′S, indicating
that the inhibitor binds at position 6′, as in Cys-loop receptors,
if it can enter the pore. Overall, the data support the proposal that
the ELIC pore is a good model for Cys-loop receptor pores if the role
of F16′ is taken into consideration.
Collapse
Affiliation(s)
- Andrew J Thompson
- Department of Biochemistry, University of Cambridge , Cambridge CB2 1QW, U.K
| | | | | | | | | |
Collapse
|
16
|
Structural dynamics of the cell wall precursor lipid II in the presence and absence of the lantibiotic nisin. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:3061-8. [PMID: 25128154 DOI: 10.1016/j.bbamem.2014.07.024] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2014] [Revised: 07/23/2014] [Accepted: 07/25/2014] [Indexed: 01/06/2023]
Abstract
Representing a physiological "Achilles' heel", the cell wall precursor lipid II (LII) is a prime target for various classes of antibiotics. Over the years LII-binding agents have been recognized as promising candidates and templates in the search for new antibacterial compounds to complement or replace existing drugs. To elucidate the molecular structural basis underlying LII functional mechanism and to better understand if and how lantibiotic binding alters the molecular behavior of LII, we performed molecular dynamics (MD) simulations of phospholipid membrane-embedded LII in the absence and presence of the LII-binding lantibiotic nisin. In a series of 2×4 independent, unbiased 100ns MD simulations we sampled the conformational dynamics of nine LII as well as nine LII-nisin complexes embedded in an aqueous 150mM NaCl/POPC phospholipid membrane environment. We found that nisin binding to LII induces a reduction of LII mobility and flexibility, an outward shift of the LII pentapeptide, an inward movement of the LII disaccharide section, and an overall deeper insertion of the LII tail group into the membrane. The latter effect might indicate an initial step in adopting a stabilizing, scaffold-like structure in the process of nisin-induced membrane leakage. At the same time nisin conformation and LII interaction remain similar to the 1WCO LII-nisin NMR solution structure.
Collapse
|
17
|
Crystallographic studies of pharmacological sites in pentameric ligand-gated ion channels. Biochim Biophys Acta Gen Subj 2014; 1850:511-23. [PMID: 24836522 DOI: 10.1016/j.bbagen.2014.05.007] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2014] [Revised: 05/05/2014] [Accepted: 05/06/2014] [Indexed: 12/18/2022]
Abstract
BACKGROUND Pentameric ligand-gated ion channels (pLGICs) mediate fast chemical transmission of nerve signals in the central and peripheral nervous system. On the functional side, these molecules respond to the binding of a neurotransmitter (glycine, GABA, acetylcholine or 5HT3) in the extracellular domain (ECD) by opening their ionotropic pore in the transmembrane domain (TMD). The response to the neurotransmitter binding can be modulated by several chemical compounds acting at topographically distinct sites, as documented by a large body of literature. Notably, these receptors are the target of several classes of world-wide prescribed drugs, including general anesthetics, smoking cessation aids, anxiolytics, anticonvulsants, muscle relaxants, hypnotics and anti-emetics. On the structural side recent progress has been made on the crystallization of pLGICs in its different allosteric states, especially pLGICs of bacterial origin. Therefore, structure-function relationships can now be discussed at the atomic level for pLGICs. SCOPE OF REVIEW This review focuses on the crystallographic structure of complexes of pLGICs with a number of ligands of pharmacological interest. First, we review structural data on two key functional aspects of these receptors: the agonist-induced activation and ion transport itself. The molecular understanding of both these functional aspects is important, as they are those that most pharmacological compounds target. Next, we describe modulation sites that have recently been documented by X-ray crystallography. Finally, we propose a simple geometric classification of all these pharmacological sites in pLGICs, based on icosahedrons. MAJOR CONCLUSIONS This review illustrates the wealth of structural insight gained by comparing all available structures of members of the pLGIC family to rationalize the pharmacology of structurally diverse drugs acting at topographically distinct sites. It will be highlighted how sites that had been described earlier using biochemical techniques can be rationalized using structural data. Surprisingly, the use of icosahedral symmetry allows to link together several modulation sites, in a way that was totally unanticipated. GENERAL SIGNIFICANCE Overall, understanding the interplay between the different modulation sites at the structural level should help the design of future drugs targeting pLGICs. This article is part of a Special Issue entitled structural biochemistry and biophysics of membrane proteins.
Collapse
|
18
|
Taly A, Hénin J, Changeux JP, Cecchini M. Allosteric regulation of pentameric ligand-gated ion channels: an emerging mechanistic perspective. Channels (Austin) 2014; 8:350-60. [PMID: 25478624 PMCID: PMC4203737 DOI: 10.4161/chan.29444] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2014] [Revised: 06/03/2014] [Accepted: 06/03/2014] [Indexed: 12/22/2022] Open
Abstract
Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communications in the nervous system by converting the binding of a chemical messenger—a neurotransmitter—into an ion flux through the postsynaptic membrane. They are oligomeric assemblies that provide prototypical examples of allosterically regulated integral membrane proteins. Here, we present an overview of the most recent advances on the signal transduction mechanism based on the X-ray structures of both prokaryotic and invertebrate eukaryotic pLGICs and atomistic Molecular Dynamics simulations. The present results suggest that ion gating involves a large structural reorganization of the molecule mediated by two distinct quaternary transitions, a global twisting and the blooming of the extracellular domain, which can be modulated by ligand binding at the topographically distinct orthosteric and allosteric sites. The emerging model of gating is consistent with a wealth of functional studies and will boost the development of novel pharmacological strategies.
Collapse
Affiliation(s)
- Antoine Taly
- Laboratoire de Biochimie Théorique; IBPC; CNRS and Université Paris Diderot; Paris, France
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique; IBPC; CNRS and Université Paris Diderot; Paris, France
| | - Jean-Pierre Changeux
- CNRS; URA 2182; F-75015 & Collège de France; Paris, France
- Kavli Institute for Brain & Mind University of California; San Diego La Jolla, CA USA
| | - Marco Cecchini
- ISIS; UMR 7006 CNRS; Université de Strasbourg; F-67083 Strasbourg Cedex, France
| |
Collapse
|
19
|
Abstract
Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communication in the nervous system and are involved in fundamental processes such as attention, learning, and memory. They are oligomeric protein assemblies that convert a chemical signal into an ion flux through the postsynaptic membrane, but the molecular mechanism of gating ions has remained elusive. Here, we present atomistic molecular dynamics simulations of the prokaryotic channels from Gloeobacter violaceus (GLIC) and Erwinia chrysanthemi (ELIC), whose crystal structures are thought to represent the active and the resting states of pLGICs, respectively, and of the eukaryotic glutamate-gated chloride channel from Caenorhabditis elegans (GluCl), whose open-channel structure was determined complexed with the positive allosteric modulator ivermectin. Structural observables extracted from the trajectories of GLIC and ELIC are used as progress variables to analyze the time evolution of GluCl, which was simulated in the absence of ivermectin starting from the structure with bound ivermectin. The trajectory of GluCl with ivermectin removed shows a sequence of structural events that couple agonist unbinding from the extracellular domain to ion-pore closing in the transmembrane domain. Based on these results, we propose a structural mechanism for the allosteric communication leading to deactivation/activation of the GluCl channel. This model of gating emphasizes the coupling between the quaternary twisting and the opening/closing of the ion pore and is likely to apply to other members of the pLGIC family.
Collapse
|
20
|
Liu Y, Zhu F. Collective diffusion model for ion conduction through microscopic channels. Biophys J 2013; 104:368-76. [PMID: 23442858 DOI: 10.1016/j.bpj.2012.11.3826] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2012] [Revised: 11/21/2012] [Accepted: 11/29/2012] [Indexed: 10/27/2022] Open
Abstract
Ion conduction through microscopic channels is of central importance in both biology and nanotechnology. To better understand the current-voltage (I-V) dependence of ion channels, here we describe and prove a collective diffusion model that quantitatively relates the spontaneous ion permeation at equilibrium to the stationary ionic fluxes driven by small voltages. The model makes it possible to determine the channel conductance in the linear I-V range from equilibrium simulations without the application of a voltage. To validate the theory, we perform molecular-dynamics simulations on two channels-a conical-shaped nanopore and the transmembrane pore of an α-hemolysin-under both equilibrium and nonequilibrium conditions. The simulations reveal substantial couplings between the motions of cations and anions, which are effectively captured by the collective coordinate in the model. Although the two channels exhibit very different linear ranges in the I-V curves, in both cases the channel conductance at small voltages is in reasonable agreement with the prediction from the equilibrium simulation. The simulations also suggest that channel charges, rather than geometric asymmetry, play a more prominent role in current rectification.
Collapse
Affiliation(s)
- Yingting Liu
- Department of Physics, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana, USA
| | | |
Collapse
|
21
|
Ashby JA, McGonigle IV, Price KL, Cohen N, Comitani F, Dougherty DA, Molteni C, Lummis SCR. GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study. Biophys J 2012. [PMID: 23200041 DOI: 10.1016/j.bpj.2012.10.016] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.
Collapse
|
22
|
Zhu F, Hummer G. Drying transition in the hydrophobic gate of the GLIC channel blocks ion conduction. Biophys J 2012; 103:219-27. [PMID: 22853899 DOI: 10.1016/j.bpj.2012.06.003] [Citation(s) in RCA: 79] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2012] [Revised: 05/31/2012] [Accepted: 06/04/2012] [Indexed: 11/19/2022] Open
Abstract
The theoretical prediction of water drying transitions near nonpolar surfaces has stimulated an intensive search for biological processes exploiting this extreme form of hydrophobicity. Here we quantitatively demonstrate that drying of a hydrophobic constriction is the major determinant of ion conductance in the GLIC pentameric ion channel. Molecular-dynamics simulations show that in the closed state, the channel conductance is ∼12 orders-of-magnitude lower than in the open state. This large drop in conductance is remarkable because even in the functionally closed conformation the pore constriction remains wide enough for the passage of sodium ions, aided by a continuous bridge of ∼12 water molecules. However, we find that the free energy cost of hydrating the hydrophobic gate is large, accounting almost entirely for the energetic barrier blocking ion passage. The free energies of transferring a sodium ion into a prehydrated gate in functionally closed and open states differ by only 1.2 kcal/mol, compared to an 11 kcal/mol difference in the costs of hydrating the hydrophobic gate. Conversely, ion desolvation effects play only minor roles in GLIC ion channel gating. Our simulations help rationalize experiments probing the gating kinetics of the nicotinic acetylcholine receptor in response to mutations of pore-lining residues. The molecular character and phase behavior of water should thus be included in quantitative descriptions of ion channel gating.
Collapse
Affiliation(s)
- Fangqiang Zhu
- Department of Physics, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana, USA.
| | | |
Collapse
|
23
|
Martínez-Delgado G, Estrada-Mondragón A, Miledi R, Martínez-Torres A. An Update on GABAρ Receptors. Curr Neuropharmacol 2011; 8:422-33. [PMID: 21629448 PMCID: PMC3080597 DOI: 10.2174/157015910793358141] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2009] [Revised: 04/08/2010] [Accepted: 06/21/2010] [Indexed: 01/29/2023] Open
Abstract
The present review discusses the functional and molecular diversity of GABAρ receptors. These receptors were originally described in the mammalian retina, and their functional role in the visual pathway has been recently elucidated; however new studies on their distribution in the brain and spinal cord have revealed that they are more spread than originally thought, and thus it will be important to determine their physiological contribution to the GABAergic transmission in other areas of the central nervous system. In addition, molecular modeling has revealed peculiar traits of these receptors that have impacted on the interpretations of the latest pharmacolgical and biophysical findings. Finally, sequencing of several vertebrate genomes has permitted a comparative analysis of the organization of the GABAρ genes.
Collapse
Affiliation(s)
- Gustavo Martínez-Delgado
- Instituto de Neurbiología, Departamento de Neurobiología Celular y Molecular, Laboratorio D15, Campus UNAM Juriquilla. Querétaro 76230, México
| | | | | | | |
Collapse
|
24
|
Zimmermann I, Dutzler R. Ligand activation of the prokaryotic pentameric ligand-gated ion channel ELIC. PLoS Biol 2011; 9:e1001101. [PMID: 21713033 PMCID: PMC3119659 DOI: 10.1371/journal.pbio.1001101] [Citation(s) in RCA: 89] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2011] [Accepted: 05/31/2011] [Indexed: 11/18/2022] Open
Abstract
While the pentameric ligand-gated ion channel ELIC has recently provided first insight into the architecture of the family at high resolution, its detailed investigation was so far prevented by the fact that activating ligands were unknown. Here we describe a study on the functional characterization of ELIC by electrophysiology and X-ray crystallography. ELIC is activated by a class of primary amines that include the neurotransmitter GABA at high micro- to millimolar concentrations. The ligands bind to a conserved site and evoke currents that slowly desensitize over time. The protein forms cation selective channels with properties that resemble the nicotinic acetylcholine receptor. The high single channel conductance and the comparably simple functional behavior make ELIC an attractive model system to study general mechanisms of ion conduction and gating in this important family of neurotransmitter receptors.
Collapse
Affiliation(s)
- Iwan Zimmermann
- Department of Biochemistry, University of Zurich, Zurich, Switzerland
| | | |
Collapse
|
25
|
Zheng W, Auerbach A. Decrypting the sequence of structural events during the gating transition of pentameric ligand-gated ion channels based on an interpolated elastic network model. PLoS Comput Biol 2011; 7:e1001046. [PMID: 21253563 PMCID: PMC3017109 DOI: 10.1371/journal.pcbi.1001046] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2010] [Accepted: 12/01/2010] [Indexed: 01/22/2023] Open
Abstract
Despite many experimental and computational studies of the gating transition of pentameric ligand-gated ion channels (pLGICs), the structural basis of how ligand binding couples to channel gating remains unknown. By using a newly developed interpolated elastic network model (iENM), we have attempted to compute a likely transition pathway from the closed- to the open-channel conformation of pLGICs as captured by the crystal structures of two prokaryotic pLGICs. The iENM pathway predicts a sequence of structural events that begins at the ligand-binding loops and is followed by the displacements of two key loops (loop 2 and loop 7) at the interface between the extracellular and transmembrane domain, the tilting/bending of the pore-lining M2 helix, and subsequent movements of M4, M3 and M1 helices in the transmembrane domain. The predicted order of structural events is in broad agreement with the Φ-value analysis of α subunit of nicotinic acetylcholine receptor mutants, which supports a conserved core mechanism for ligand-gated channel opening in pLGICs. Further perturbation analysis has supported the critical role of certain intra-subunit and inter-subunit interactions in dictating the above sequence of events.
Collapse
Affiliation(s)
- Wenjun Zheng
- Department of Physics, University at Buffalo, Buffalo, New York, United States of America
- * E-mail:
| | - Anthony Auerbach
- Department of Physiology and Biophysics, University at Buffalo, Buffalo, New York, United States of America
| |
Collapse
|
26
|
Dellisanti CD, Hanson SM, Chen L, Czajkowski C. Packing of the extracellular domain hydrophobic core has evolved to facilitate pentameric ligand-gated ion channel function. J Biol Chem 2010; 286:3658-70. [PMID: 21098036 DOI: 10.1074/jbc.m110.156851] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Protein function depends on conformational flexibility and folding stability. Loose packing of hydrophobic cores is not infrequent in proteins, as the enhanced flexibility likely contributes to their biological function. Here, using experimental and computational approaches, we show that eukaryotic pentameric ligand-gated ion channels are characterized by loose packing of their extracellular domain β-sandwich cores, and that loose packing contributes to their ability to rapidly switch from closed to open channel states in the presence of ligand. Functional analyses of GABA(A) receptors show that increasing the β-core packing disrupted GABA-mediated currents, with impaired GABA efficacy and slowed GABA current activation and desensitization. We propose that loose packing of the hydrophobic β-core developed as an evolutionary strategy aimed to facilitate the allosteric mechanisms of eukaryotic pentameric ligand-gated ion channels.
Collapse
Affiliation(s)
- Cosma D Dellisanti
- Department of Physiology, University of Wisconsin, Madison, Wisconsin 53711, USA
| | | | | | | |
Collapse
|
27
|
Abstract
Nerve signaling in humans and chemical sensing in bacteria both rely on the controlled opening and closing of the ion-conducting pore in pentameric ligand-gated ion channels. With the help of a multiscale simulation approach that combines mixed elastic network model calculations with molecular dynamics simulations, we study the opening and closing of the pore in Gloeobacter violaceus channel GLIC at atomic resolution. In our simulations of the GLIC transmembrane domain, we first verify that the two endpoints of the transition are open and closed to sodium ion conduction, respectively. We then show that a two-stage tilting of the pore-lining helices induces cooperative drying and iris-like closing of the channel pore. From the free energy profile of the gating transition and from unrestrained simulations, we conclude that the pore of the isolated GLIC transmembrane domain closes spontaneously. The mechanical work of opening the pore is performed primarily on the M2-M3 loop. Strong interactions of this short and conserved loop with the extracellular domain are therefore crucial to couple ligand binding to channel opening.
Collapse
Affiliation(s)
- Fangqiang Zhu
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
| | - Gerhard Hummer
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
| |
Collapse
|
28
|
Majd S, Yusko EC, Billeh YN, Macrae MX, Yang J, Mayer M. Applications of biological pores in nanomedicine, sensing, and nanoelectronics. Curr Opin Biotechnol 2010; 21:439-76. [PMID: 20561776 PMCID: PMC3121537 DOI: 10.1016/j.copbio.2010.05.002] [Citation(s) in RCA: 237] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2010] [Revised: 05/03/2010] [Accepted: 05/06/2010] [Indexed: 12/29/2022]
Abstract
Biological protein pores and pore-forming peptides can generate a pathway for the flux of ions and other charged or polar molecules across cellular membranes. In nature, these nanopores have diverse and essential functions that range from maintaining cell homeostasis and participating in cell signaling to activating or killing cells. The combination of the nanoscale dimensions and sophisticated - often regulated - functionality of these biological pores make them particularly attractive for the growing field of nanobiotechnology. Applications range from single-molecule sensing to drug delivery and targeted killing of malignant cells. Potential future applications may include the use of nanopores for single strand DNA sequencing and for generating bio-inspired, and possibly, biocompatible visual detection systems and batteries. This article reviews the current state of applications of pore-forming peptides and proteins in nanomedicine, sensing, and nanoelectronics.
Collapse
Affiliation(s)
- Sheereen Majd
- Department of Biomedical Engineering, University of Michigan, 1101 Beal Avenue, Ann Arbor, Michigan 48109-2110, USA
| | | | | | | | | | | |
Collapse
|
29
|
Abstract
Large chemical and biological systems such as fuel cells, ion channels, molecular motors, and viruses are of great importance to the scientific community and public health. Typically, these complex systems in conjunction with their aquatic environment pose a fabulous challenge to theoretical description, simulation, and prediction. In this work, we propose a differential geometry based multiscale paradigm to model complex macromolecular systems, and to put macroscopic and microscopic descriptions on an equal footing. In our approach, the differential geometry theory of surfaces and geometric measure theory are employed as a natural means to couple the macroscopic continuum mechanical description of the aquatic environment with the microscopic discrete atomistic description of the macromolecule. Multiscale free energy functionals, or multiscale action functionals are constructed as a unified framework to derive the governing equations for the dynamics of different scales and different descriptions. Two types of aqueous macromolecular complexes, ones that are near equilibrium and others that are far from equilibrium, are considered in our formulations. We show that generalized Navier-Stokes equations for the fluid dynamics, generalized Poisson equations or generalized Poisson-Boltzmann equations for electrostatic interactions, and Newton's equation for the molecular dynamics can be derived by the least action principle. These equations are coupled through the continuum-discrete interface whose dynamics is governed by potential driven geometric flows. Comparison is given to classical descriptions of the fluid and electrostatic interactions without geometric flow based micro-macro interfaces. The detailed balance of forces is emphasized in the present work. We further extend the proposed multiscale paradigm to micro-macro analysis of electrohydrodynamics, electrophoresis, fuel cells, and ion channels. We derive generalized Poisson-Nernst-Planck equations that are coupled to generalized Navier-Stokes equations for fluid dynamics, Newton's equation for molecular dynamics, and potential and surface driving geometric flows for the micro-macro interface. For excessively large aqueous macromolecular complexes in chemistry and biology, we further develop differential geometry based multiscale fluid-electro-elastic models to replace the expensive molecular dynamics description with an alternative elasticity formulation.
Collapse
Affiliation(s)
- Guo-Wei Wei
- Department of Mathematics, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
30
|
One-microsecond molecular dynamics simulation of channel gating in a nicotinic receptor homologue. Proc Natl Acad Sci U S A 2010; 107:6275-80. [PMID: 20308576 DOI: 10.1073/pnas.1001832107] [Citation(s) in RCA: 137] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Recently discovered bacterial homologues of eukaryotic pentameric ligand-gated ion channels, such as the Gloeobacter violaceus receptor (GLIC), are increasingly used as structural and functional models of signal transduction in the nervous system. Here we present a one-microsecond-long molecular dynamics simulation of the GLIC channel pH stimulated gating mechanism. The crystal structure of GLIC obtained at acidic pH in an open-channel form is equilibrated in a membrane environment and then instantly set to neutral pH. The simulation shows a channel closure that rapidly takes place at the level of the hydrophobic furrow and a progressively increasing quaternary twist. Two major events are captured during the simulation. They are initiated by local but large fluctuations in the pore, taking place at the top of the M2 helix, followed by a global tertiary relaxation. The two-step transition of the first subunit starts within the first 50 ns of the simulation and is followed at 450 ns by its immediate neighbor in the pentamer, which proceeds with a similar scenario. This observation suggests a possible two-step domino-like tertiary mechanism that takes place between adjacent subunits. In addition, the dynamical properties of GLIC described here offer an interpretation of the paradoxical properties of a permeable A13'F mutant whose crystal structure determined at 3.15 A shows a pore too narrow to conduct ions.
Collapse
|
31
|
Zhu F, Hummer G. Gating transition of pentameric ligand-gated ion channels. Biophys J 2010; 97:2456-63. [PMID: 19883588 DOI: 10.1016/j.bpj.2009.08.020] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2009] [Revised: 08/12/2009] [Accepted: 08/13/2009] [Indexed: 02/03/2023] Open
Abstract
Pentameric ligand-gated ion channels are an important family of membrane proteins and play key roles in physiological processes, including signal transduction at chemical synapses. Here, we study the conformational changes associated with the opening and closing of the channel pore. Based on recent crystal structures of two prokaryotic members of the family in open and closed states, respectively, mixed elastic network models are constructed for the transmembrane domain. To explore the conformational changes in the gating transition, a coarse-grained transition path is computed that smoothly connects the closed and open conformations of the channel. We find that the conformational transition involves no major rotations of the transmembrane helices, and is instead characterized by a concerted tilting of helices M2 and M3. In addition, helix M2 changes its bending state, which results in an early closure of the pore during the open-to-closed transition.
Collapse
Affiliation(s)
- Fangqiang Zhu
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, USA
| | | |
Collapse
|
32
|
Song C, Corry B. Ion conduction in ligand-gated ion channels: Brownian dynamics studies of four recent crystal structures. Biophys J 2010; 98:404-11. [PMID: 20141753 PMCID: PMC2814205 DOI: 10.1016/j.bpj.2009.10.032] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2009] [Revised: 10/21/2009] [Accepted: 10/23/2009] [Indexed: 01/12/2023] Open
Abstract
Four x-ray crystal structures of prokaryotic homologs of ligand-gated ion channels have recently been determined: ELIC from Erwinia chrysanthemi, two structures of a proton-activated channel from Gloebacter violaceus (GLIC1 and GLIC2) and that of the E221A mutant (GLIC1M). The availability of numerous structures of channels in this family allows for aspects of channel gating and ion conduction to be examined. Here, we determine the likely conduction states of the four structures as well as IV curves, ion selectivity, and steps involved in ion permeation by performing extensive Brownian dynamics simulations. Our results show that the ELIC structure is indeed nonconductive, but that GLIC1 and GLIC1M are both conductive of ions with properties different from those seen in experimental studies of the channel. GLIC2 appears to reflect an open state of the channel with a predicted conductance of 10.8-12.4 pS in 140 mM NaCl solution, which is comparable to the experimental value 8 +/- 2 pS. The extracellular domain of the channel is shown to have an important influence on the channel current, but a less significant role in ion selectivity.
Collapse
Affiliation(s)
| | - Ben Corry
- School of Biomedical, Biomolecular and Chemical Sciences, The University of Western Australia, Crawley, Western Australia
| |
Collapse
|
33
|
Chen L. In pursuit of the high-resolution structure of nicotinic acetylcholine receptors. J Physiol 2009; 588:557-64. [PMID: 19995851 DOI: 10.1113/jphysiol.2009.184085] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
The nicotinic acetylcholine receptor (nAChR) has been studied extensively for well over four decades because of its important physiological roles and medical relevance. A large body of data from biochemical and biophysical studies are now available. The structural information, which is needed to integrate existing data to address the mechanism and function of nAChRs, started to emerge in recent years. Structural studies of acetylcholine binding proteins (AChBPs) have greatly facilitated the study of nAChRs. The recently determined crystal structures of the prokaryotic homologues of nAChRs will probably have similar impact over time. However, a direct structural model of nAChRs at high resolution will be important for mechanistic studies and drug development. Here we will review some of the recent efforts in this area and use the high-resolution structure of the extracellular domains of nAChR alpha1 to illustrate the potential insights one may gain at higher resolution.
Collapse
Affiliation(s)
- Lin Chen
- Molecular and Computational Biology, Departments of Biological Sciences and Chemistry, Norris Comprehensive Cancer Center, Keck School of Medicine, University of Southern California, Los Angeles, CA 90089, USA.
| |
Collapse
|