1
|
Alberini G, Alexis Paz S, Corradi B, Abrams CF, Benfenati F, Maragliano L. Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels. J Chem Theory Comput 2023; 19:2953-2972. [PMID: 37116214 DOI: 10.1021/acs.jctc.2c00990] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na+ translocation through the SF of the homology-based Nav1.2 model and the cryo-EM Nav1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na+ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav channels.
Collapse
Affiliation(s)
- Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi 10, 16132 Genova, Italy
- IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi 10, 16132 Genova, Italy
| | - Sergio Alexis Paz
- Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, X5000HUA Córdoba, Argentina
- Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto de Fisicoquímica de Córdoba (INFIQC), X5000HUA Córdoba, Argentina
| | - Beatrice Corradi
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi 10, 16132 Genova, Italy
- Department of Experimental Medicine, Università degli Studi di Genova, Viale Benedetto XV 3, 16132 Genova, Italy
| | - Cameron F Abrams
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi 10, 16132 Genova, Italy
- IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi 10, 16132 Genova, Italy
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi 10, 16132 Genova, Italy
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| |
Collapse
|
2
|
Wang G, Xu L, Chen H, Liu Y, Pan P, Hou T. Recent advances in computational studies on voltage‐gated sodium channels: Drug design and mechanism studies. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1641] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Affiliation(s)
- Gaoang Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University College of Pharmaceutical Sciences, Zhejiang University Hangzhou Zhejiang China
| | - Lei Xu
- Institute of Bioinformatics and Medical Engineering School of Electrical and Information Engineering, Jiangsu University of Technology Changzhou Jiangsu China
| | - Haiyi Chen
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University College of Pharmaceutical Sciences, Zhejiang University Hangzhou Zhejiang China
| | - Yifei Liu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University College of Pharmaceutical Sciences, Zhejiang University Hangzhou Zhejiang China
| | - Peichen Pan
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University College of Pharmaceutical Sciences, Zhejiang University Hangzhou Zhejiang China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University College of Pharmaceutical Sciences, Zhejiang University Hangzhou Zhejiang China
| |
Collapse
|
3
|
Structure of human Na V1.6 channel reveals Na + selectivity and pore blockade by 4,9-anhydro-tetrodotoxin. Nat Commun 2023; 14:1030. [PMID: 36823201 PMCID: PMC9950489 DOI: 10.1038/s41467-023-36766-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 02/15/2023] [Indexed: 02/25/2023] Open
Abstract
The sodium channel NaV1.6 is widely expressed in neurons of the central and peripheral nervous systems, which plays a critical role in regulating neuronal excitability. Dysfunction of NaV1.6 has been linked to epileptic encephalopathy, intellectual disability and movement disorders. Here we present cryo-EM structures of human NaV1.6/β1/β2 alone and complexed with a guanidinium neurotoxin 4,9-anhydro-tetrodotoxin (4,9-ah-TTX), revealing molecular mechanism of NaV1.6 inhibition by the blocker. The apo-form structure reveals two potential Na+ binding sites within the selectivity filter, suggesting a possible mechanism for Na+ selectivity and conductance. In the 4,9-ah-TTX bound structure, 4,9-ah-TTX binds to a pocket similar to the tetrodotoxin (TTX) binding site, which occupies the Na+ binding sites and completely blocks the channel. Molecular dynamics simulation results show that subtle conformational differences in the selectivity filter affect the affinity of TTX analogues. Taken together, our results provide important insights into NaV1.6 structure, ion conductance, and inhibition.
Collapse
|
4
|
Berselli A, Alberini G, Benfenati F, Maragliano L. Computational study of ion permeation through claudin-4 paracellular channels. Ann N Y Acad Sci 2022; 1516:162-174. [PMID: 35811406 PMCID: PMC9796105 DOI: 10.1111/nyas.14856] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Claudins (Cldns) form a large family of protein homologs that are essential for the assembly of paracellular tight junctions (TJs), where they form channels or barriers with tissue-specific selectivity for permeants. In contrast to several family members whose physiological role has been identified, the function of claudin 4 (Cldn4) remains elusive, despite experimental evidence suggesting that it can form anion-selective TJ channels in the renal epithelium. Computational approaches have recently been employed to elucidate the molecular basis of Cldns' function, and hence could help in clarifying the role of Cldn4. In this work, we use structural modeling and all-atom molecular dynamics simulations to transfer two previously introduced structural models of Cldn-based paracellular complexes to Cldn4 to reproduce a paracellular anion channel. Free energy calculations for ionic transport through the pores allow us to establish the thermodynamic properties driving the ion-selectivity of the structures. While one model shows a cavity permeable to chloride and repulsive to cations, the other forms barrier to the passage of all the major physiological ions. Furthermore, our results confirm the charge selectivity role of the residue Lys65 in the first extracellular loop of the protein, rationalizing Cldn4 control of paracellular permeability.
Collapse
Affiliation(s)
- Alessandro Berselli
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe)Istituto Italiano di TecnologiaGenovaItaly,Department of Experimental MedicineUniversità degli Studi di GenovaGenovaItaly
| | - Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe)Istituto Italiano di TecnologiaGenovaItaly,IRCCS Ospedale Policlinico San MartinoGenovaItaly
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe)Istituto Italiano di TecnologiaGenovaItaly,IRCCS Ospedale Policlinico San MartinoGenovaItaly
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe)Istituto Italiano di TecnologiaGenovaItaly,Department of Life and Environmental SciencesPolytechnic University of MarcheAnconaItaly
| |
Collapse
|
5
|
Zhu Z, Deng Z, Wang Q, Wang Y, Zhang D, Xu R, Guo L, Wen H. Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design. Front Pharmacol 2022; 13:939555. [PMID: 35837274 PMCID: PMC9275593 DOI: 10.3389/fphar.2022.939555] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 06/07/2022] [Indexed: 11/13/2022] Open
Abstract
Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.
Collapse
Affiliation(s)
- Zhengdan Zhu
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- Beijing Institute of Big Data Research, Beijing, China
| | - Zhenfeng Deng
- DP Technology, Beijing, China
- School of Pharmaceutical Sciences, Peking University, Beijing, China
| | | | | | - Duo Zhang
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- DP Technology, Beijing, China
| | - Ruihan Xu
- DP Technology, Beijing, China
- National Engineering Research Center of Visual Technology, Peking University, Beijing, China
| | | | - Han Wen
- DP Technology, Beijing, China
| |
Collapse
|
6
|
Guardiani C, Cecconi F, Chiodo L, Cottone G, Malgaretti P, Maragliano L, Barabash ML, Camisasca G, Ceccarelli M, Corry B, Roth R, Giacomello A, Roux B. Computational methods and theory for ion channel research. ADVANCES IN PHYSICS: X 2022; 7:2080587. [PMID: 35874965 PMCID: PMC9302924 DOI: 10.1080/23746149.2022.2080587] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 05/15/2022] [Indexed: 06/15/2023] Open
Abstract
Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.
Collapse
Affiliation(s)
- C. Guardiani
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
| | - F. Cecconi
- CNR - Istituto dei Sistemi Complessi, Rome, Italy and Istituto Nazionale di Fisica Nucleare, INFN, Roma1 section. 00185, Roma, Italy
| | - L. Chiodo
- Department of Engineering, Campus Bio-Medico University, Rome, Italy
| | - G. Cottone
- Department of Physics and Chemistry-Emilio Segrè, University of Palermo, Palermo, Italy
| | - P. Malgaretti
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Forschungszentrum Jülich, Erlangen, Germany
| | - L. Maragliano
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy, and Center for Synaptic Neuroscience and Technology, Istituto Italiano di Tecnologia, Genova, Italy
| | - M. L. Barabash
- Department of Materials Science and Nanoengineering, Rice University, Houston, TX 77005, USA
| | - G. Camisasca
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
- Dipartimento di Fisica, Università Roma Tre, Rome, Italy
| | - M. Ceccarelli
- Department of Physics and CNR-IOM, University of Cagliari, Monserrato 09042-IT, Italy
| | - B. Corry
- Research School of Biology, The Australian National University, Canberra, ACT 2600, Australia
| | - R. Roth
- Institut Für Theoretische Physik, Eberhard Karls Universität Tübingen, Tübingen, Germany
| | - A. Giacomello
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
| | - B. Roux
- Department of Biochemistry & Molecular Biology, University of Chicago, Chicago IL, USA
| |
Collapse
|
7
|
Gibby WAT, Barabash ML, Guardiani C, Luchinsky DG, McClintock PVE. Physics of Selective Conduction and Point Mutation in Biological Ion Channels. PHYSICAL REVIEW LETTERS 2021; 126:218102. [PMID: 34114848 DOI: 10.1103/physrevlett.126.218102] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 04/23/2021] [Indexed: 06/12/2023]
Abstract
We introduce a statistical and linear response theory of selective conduction in biological ion channels with multiple binding sites and possible point mutation. We derive an effective grand-canonical ensemble and generalized Einstein relations for the selectivity filter, assuming strongly coordinated ionic motion, and allowing for ionic Coulomb blockade. The theory agrees well with data from the KcsA K^{+} channel and a mutant. We show that the Eisenman relations for thermodynamic selectivity follow from the condition for fast conduction and find that maximum conduction requires the binding sites to be nearly identical.
Collapse
Affiliation(s)
- W A T Gibby
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| | - M L Barabash
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| | - C Guardiani
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
- Department of Mechanical and Aerospace Engineering, Sapienza University, Rome 00184, Italy
| | - D G Luchinsky
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
- KBR Inc., Ames Research Center, Moffett Field, Mountain View, California 94035, USA
| | - P V E McClintock
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| |
Collapse
|
8
|
Gibby WAT, Fedorenko OA, Guardiani C, Barabash ML, Mumby T, Roberts SK, Luchinsky DG, McClintock PVE. Application of a Statistical and Linear Response Theory to Multi-Ion Na + Conduction in NaChBac. ENTROPY (BASEL, SWITZERLAND) 2021; 23:249. [PMID: 33670053 PMCID: PMC7926348 DOI: 10.3390/e23020249] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Revised: 02/03/2021] [Accepted: 02/11/2021] [Indexed: 01/04/2023]
Abstract
Biological ion channels are fundamental to maintaining life. In this manuscript we apply our recently developed statistical and linear response theory to investigate Na+ conduction through the prokaryotic Na+ channel NaChBac. This work is extended theoretically by the derivation of ionic conductivity and current in an electrochemical gradient, thus enabling us to compare to a range of whole-cell data sets performed on this channel. Furthermore, we also compare the magnitudes of the currents and populations at each binding site to previously published single-channel recordings and molecular dynamics simulations respectively. In doing so, we find excellent agreement between theory and data, with predicted energy barriers at each of the four binding sites of ∼4,2.9,3.6, and 4kT.
Collapse
Affiliation(s)
- William A. T. Gibby
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
| | - Olena A. Fedorenko
- School of Life Sciences, University of Nottingham, Nottingham NG7 2UH, UK;
- Division of Biomedical and Life Sciences, Lancaster University, Lancaster LA1 4YQ, UK;
| | - Carlo Guardiani
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
- Department of Mechanical and Aerospace Engineering, Sapienza University, 00185 Rome, Italy
| | - Miraslau L. Barabash
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
| | - Thomas Mumby
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
| | - Stephen K. Roberts
- Division of Biomedical and Life Sciences, Lancaster University, Lancaster LA1 4YQ, UK;
| | - Dmitry G. Luchinsky
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
- KBR Inc., Ames Research Center, Mountain View, CA 94035, USA
| | - Peter V. E. McClintock
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; (C.G.); (M.L.B.); (T.M.); (D.G.L.)
| |
Collapse
|
9
|
Changes in Ion Selectivity Following the Asymmetrical Addition of Charge to the Selectivity Filter of Bacterial Sodium Channels. ENTROPY 2020; 22:e22121390. [PMID: 33316962 PMCID: PMC7764494 DOI: 10.3390/e22121390] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 11/30/2020] [Accepted: 12/07/2020] [Indexed: 12/30/2022]
Abstract
Voltage-gated sodium channels (NaVs) play fundamental roles in eukaryotes, but their exceptional size hinders their structural resolution. Bacterial NaVs are simplified homologues of their eukaryotic counterparts, but their use as models of eukaryotic Na+ channels is limited by their homotetrameric structure at odds with the asymmetric Selectivity Filter (SF) of eukaryotic NaVs. This work aims at mimicking the SF of eukaryotic NaVs by engineering radial asymmetry into the SF of bacterial channels. This goal was pursued with two approaches: the co-expression of different monomers of the NaChBac bacterial channel to induce the random assembly of heterotetramers, and the concatenation of four bacterial monomers to form a concatemer that can be targeted by site-specific mutagenesis. Patch-clamp measurements and Molecular Dynamics simulations showed that an additional gating charge in the SF leads to a significant increase in Na+ and a modest increase in the Ca2+ conductance in the NavMs concatemer in agreement with the behavior of the population of random heterotetramers with the highest proportion of channels with charge -5e. We thus showed that charge, despite being important, is not the only determinant of conduction and selectivity, and we created new tools extending the use of bacterial channels as models of eukaryotic counterparts.
Collapse
|
10
|
Fedorenko OA, Kaufman IK, Gibby WAT, Barabash ML, Luchinsky DG, Roberts SK, McClintock PVE. Ionic Coulomb blockade and the determinants of selectivity in the NaChBac bacterial sodium channel. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2020; 1862:183301. [PMID: 32360369 DOI: 10.1016/j.bbamem.2020.183301] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2019] [Revised: 01/30/2020] [Accepted: 03/29/2020] [Indexed: 10/24/2022]
Abstract
Mutation-induced transformations of conductivity and selectivity in NaChBac bacterial channels are studied experimentally and interpreted within the framework of ionic Coulomb blockade (ICB), while also taking account of resonant quantised dehydration (QD) and site protonation. Site-directed mutagenesis and whole-cell patch-clamp experiments are used to investigate how the fixed charge Qf at the selectivity filter (SF) affects both valence selectivity and same-charge selectivity. The new ICB/QD model predicts that increasing ∣Qf∣ should lead to a shift in selectivity sequences toward larger ion sizes, in agreement with the present experiments and with earlier work. Comparison of the model with experimental data leads to the introduction of an effective charge Qf∗ at the SF, which was found to differ between Aspartate and Glutamate charged rings, and also to depend on position within the SF. It is suggested that protonation of the residues within the restricted space of the SF is important in significantly reducing the effective charge of the EEEE ring. Values of Qf∗ derived from experiments on divalent blockade agree well with expectations based on the ICB/QD model and have led to the first demonstration of ICB oscillations in Ca2+ conduction as a function of the fixed charge. Preliminary studies of the dependence of Ca2+ conduction on pH are qualitatively consistent with the predictions of the model.
Collapse
Affiliation(s)
- O A Fedorenko
- Division of Biomedical and Life Sciences, Lancaster University, Lancaster LA1 4YQ, UK; School of Life Sciences, University of Nottingham, Nottingham NG7 2UH, UK.
| | - I Kh Kaufman
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
| | - W A T Gibby
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK.
| | - M L Barabash
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK.
| | - D G Luchinsky
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK; SGT, Inc., Greenbelt, MD 20770, USA.
| | - S K Roberts
- Division of Biomedical and Life Sciences, Lancaster University, Lancaster LA1 4YQ, UK.
| | - P V E McClintock
- Department of Physics, Lancaster University, Lancaster LA1 4YB, UK.
| |
Collapse
|
11
|
Damjanovic A, Chen AY, Rosenberg RL, Roe DR, Wu X, Brooks BR. Protonation state of the selectivity filter of bacterial voltage‐gated sodium channels is modulated by ions. Proteins 2019; 88:527-539. [DOI: 10.1002/prot.25831] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Revised: 09/03/2019] [Accepted: 09/17/2019] [Indexed: 01/28/2023]
Affiliation(s)
- Ana Damjanovic
- Department of BiophysicsJohns Hopkins University Baltimore Maryland
| | - Ada Y. Chen
- Department of PhysicsJohns Hopkins University Baltimore Maryland
| | | | - Daniel R. Roe
- Laboratory of Computational Biology, National Heart, Lung and Blood InstituteNational Institutes of Health Bethesda Maryland
| | - Xiongwu Wu
- Laboratory of Computational Biology, National Heart, Lung and Blood InstituteNational Institutes of Health Bethesda Maryland
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood InstituteNational Institutes of Health Bethesda Maryland
| |
Collapse
|
12
|
Flood E, Boiteux C, Lev B, Vorobyov I, Allen TW. Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation. Chem Rev 2019; 119:7737-7832. [DOI: 10.1021/acs.chemrev.8b00630] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Emelie Flood
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Céline Boiteux
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Bogdan Lev
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Igor Vorobyov
- Department of Physiology & Membrane Biology/Department of Pharmacology, University of California, Davis, 95616, United States
| | - Toby W. Allen
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| |
Collapse
|
13
|
Ngo ST, Vu KB, Bui LM, Vu VV. Effective Estimation of Ligand-Binding Affinity Using Biased Sampling Method. ACS OMEGA 2019; 4:3887-3893. [PMID: 31459599 PMCID: PMC6648447 DOI: 10.1021/acsomega.8b03258] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Accepted: 02/08/2019] [Indexed: 05/09/2023]
Abstract
The binding between two biomolecules is one of the most critical factors controlling many bioprocesses. Therefore, it is of great interest to derive a reliable method to calculate the free binding energy between two biomolecules. In this work, we have demonstrated that the binding affinity of ligands to proteins can be determined through biased sampling simulations. The umbrella sampling (US) method was applied on 20 protein-ligand complexes, including the cathepsin K (CTSK), type II dehydroquinase (DHQase), heat shock protein 90 (HSP90), and factor Xa (FXa) systems. The ligand-binding affinity was evaluated as the difference between the largest and smallest values of the free-energy curve, which was obtained via a potential of mean force analysis. The calculated affinities differ sizably from the previously reported experimental values, with an average difference of ∼3.14 kcal/mol. However, the calculated results are in good correlation with the experimental data, with correlation coefficients of 0.76, 0.87, 0.96, and 0.97 for CTSK, DHQase, HSP90, and FXa, respectively. Thus, the binding free energy of a new ligand can be reliably estimated using our US approach. Furthermore, the root-mean-square errors (RMSEs) of binding affinity of these systems are 1.13, 0.90, 0.37, and 0.25 kcal/mol, for CTSK, DHQase, HSP90, and FXa, respectively. The small RMSE values indicate the good precision of the biased sampling method that can distinguish the ligands exhibiting similar binding affinities.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical
and Computational Biophysics, Ton Duc Thang
University, Ho Chi Minh City 7000000, Vietnam
- Faculty
of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 7000000, Vietnam
| | - Khanh B. Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam
| | - Le Minh Bui
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam
| | - Van V. Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam
| |
Collapse
|
14
|
Callahan KM, Roux B. Molecular Dynamics of Ion Conduction through the Selectivity Filter of the Na VAb Sodium Channel. J Phys Chem B 2018; 122:10126-10142. [PMID: 30351118 DOI: 10.1021/acs.jpcb.8b09678] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The determination of the atomic structures of voltage-gated bacterial sodium channels using X-ray crystallography has provided a first view of this family of membrane proteins. Molecular dynamics simulations offer one approach to clarify the underlying mechanism of permeation and selectivity in these channels. However, it appears that the intracellular gate of the pore domain is either closed or only open partially in the available X-ray structures. The lack of structure with a fully open intracellular gate poses a special challenge to computational studies aimed at simulating ion conduction. To circumvent this problem, we simulated a model of the NaVAb channel in which the transmembrane S5 and S6 helices of the pore domain have been truncated to provide direct open access to the intracellular entryway to the pore. Molecular dynamics simulations were carried out over a range of membrane potential and ion concentration of sodium and potassium. The simulations show that the NaVAb selectivity filter is essentially a cationic pore supporting the conduction of ions at a rate comparable to aqueous diffusion with no significant selectivity for sodium. Conductance and selectivity vary as a function of ion concentration for both cations. Permeation occurs primarily via a knock-on mechanism for both sodium and potassium, although the ion ordering in single file along the pore is not strictly maintained. The character of the outward current appears quite different from the inward current, with a buildup on ions in the selectivity filter prior to escape toward the extracellular side, indicating the presence of a rectification effect that is overcome by nonphysiological applied voltages.
Collapse
Affiliation(s)
- Karen M Callahan
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
15
|
Alberini G, Benfenati F, Maragliano L. Molecular Dynamics Simulations of Ion Selectivity in a Claudin-15 Paracellular Channel. J Phys Chem B 2018; 122:10783-10792. [DOI: 10.1021/acs.jpcb.8b06484] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy
- Department of Experimental Medicine, Università degli Studi di Genova, Viale Benedetto XV, 3, 16132 Genova, Italy
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy
- IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi, 10, 16132 Genova, Italy
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy
- IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi, 10, 16132 Genova, Italy
| |
Collapse
|
16
|
Boiteux C, Flood E, Allen TW. Comparison of permeation mechanisms in sodium-selective ion channels. Neurosci Lett 2018; 700:3-8. [PMID: 29807068 DOI: 10.1016/j.neulet.2018.05.036] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2018] [Revised: 05/24/2018] [Accepted: 05/25/2018] [Indexed: 12/27/2022]
Abstract
Voltage-gated sodium channels are the molecular components of electrical signaling in the body, yet the molecular origins of Na+-selective transport remain obscured by diverse protein chemistries within this family of ion channels. In particular, bacterial and mammalian sodium channels are known to exhibit similar relative ion permeabilities for Na+ over K+ ions, despite their distinct signature EEEE and DEKA sequences. Atomic-level molecular dynamics simulations using high-resolution bacterial channel structures and mammalian channel models have begun to describe how these sequences lead to analogous high field strength ion binding sites that drive Na+ conduction. Similar complexes have also been identified in unrelated acid sensing ion channels involving glutamate and aspartate side chains that control their selectivity. These studies suggest the possibility of a common origin for Na+ selective binding and transport.
Collapse
Affiliation(s)
- Céline Boiteux
- School of Science, RMIT University, Melbourne, Australia
| | - Emelie Flood
- School of Science, RMIT University, Melbourne, Australia
| | - Toby W Allen
- School of Science, RMIT University, Melbourne, Australia.
| |
Collapse
|
17
|
Sakipov S, Sobolevsky AI, Kurnikova MG. Ion Permeation Mechanism in Epithelial Calcium Channel TRVP6. Sci Rep 2018; 8:5715. [PMID: 29632318 PMCID: PMC5890290 DOI: 10.1038/s41598-018-23972-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Accepted: 03/20/2018] [Indexed: 12/14/2022] Open
Abstract
Calcium is the most abundant metal in the human body that plays vital roles as a cellular electrolyte as well as the smallest and most frequently used signaling molecule. Calcium uptake in epithelial tissues is mediated by tetrameric calcium-selective transient receptor potential (TRP) channels TRPV6 that are implicated in a variety of human diseases, including numerous forms of cancer. We used TRPV6 crystal structures as templates for molecular dynamics simulations to identify ion binding sites and to study the permeation mechanism of calcium and other ions through TRPV6 channels. We found that at low Ca2+ concentrations, a single calcium ion binds at the selectivity filter narrow constriction formed by aspartates D541 and allows Na+ permeation. In the presence of ions, no water binds to or crosses the pore constriction. At high Ca2+ concentrations, calcium permeates the pore according to the knock-off mechanism that includes formation of a short-lived transition state with three calcium ions bound near D541. For Ba2+, the transition state lives longer and the knock-off permeation occurs slower. Gd3+ binds at D541 tightly, blocks the channel and prevents Na+ from permeating the pore. Our results provide structural foundations for understanding permeation and block in tetrameric calcium-selective ion channels.
Collapse
Affiliation(s)
- Serzhan Sakipov
- Chemistry Department, Carnegie Mellon University, 4400 Fifth Ave., Pittsburgh, PA, 15213, USA
| | - Alexander I Sobolevsky
- Department of Biochemistry and Molecular Biophysics, Columbia University, 650 West 168th St., New York, NY, 10032, USA
| | - Maria G Kurnikova
- Chemistry Department, Carnegie Mellon University, 4400 Fifth Ave., Pittsburgh, PA, 15213, USA.
| |
Collapse
|
18
|
Guardiani C, Fedorenko OA, Roberts SK, Khovanov IA. On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation. Phys Chem Chem Phys 2018; 19:29840-29854. [PMID: 29090695 DOI: 10.1039/c7cp05928k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Ion channel selectivity is essential for their function, yet the molecular basis of a channel's ability to select between ions is still rather controversial. In this work, using a combination of molecular dynamics simulations and electrophysiological current measurements we analyze the ability of the NaChBac channel to discriminate between calcium and sodium. Our simulations show that a single calcium ion can access the Selectivity Filter (SF) interacting so strongly with the glutamate ring so as to remain blocked inside. This is consistent with the tiny calcium currents recorded in our patch-clamp experiments. Two reasons explain this scenario. The first is the higher free energy of ion/SF binding of Ca2+ with respect to Na+. The second is the strong electrostatic repulsion exerted by the resident ion that turns back a second potentially incoming Ca2+, preventing the knock-on permeation mechanism. Finally, we analyzed the possibility of the Anomalous Mole Fraction Effect (AMFE), i.e. the ability of micromolar Ca2+ concentrations to block Na+ currents. Current measurements in Na+/Ca2+ mixed solutions excluded the AMFE, in agreement with metadynamics simulations showing the ability of a sodium ion to by-pass and partially displace the resident calcium. Our work supports a new scenario for Na+/Ca2+ selectivity in the bacterial sodium channel, challenging the traditional notion of an exclusion mechanism strictly confining Ca2+ ions outside the channel.
Collapse
Affiliation(s)
- Carlo Guardiani
- School of Engineering, University of Warwick, Coventry CV4 7AL, UK.
| | | | | | | |
Collapse
|
19
|
Kaufman IK, Fedorenko OA, Luchinsky DG, Gibby WA, Roberts SK, McClintock PV, Eisenberg RS. Ionic Coulomb blockade and anomalous mole fraction effect in the NaChBac bacterial ion channel and its charge-varied mutants. ACTA ACUST UNITED AC 2017. [DOI: 10.1051/epjnbp/2017003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|