Side-chain conformation at the selectivity filter shapes the permeation free-energy landscape of an ion channel

Significance Our mutagenesis work on the ring of glutamates in the charge-selectivity filter region of the muscle nicotinic acetylcholine receptor led us to propose that the rate at which ions permeate through the open channel depends not only on the number of these glutamates, but also on the conformation of their side chains. Because our inferences were made on the basis of electrophysiological observations, however, we decided to test the plausibility of these ideas using computer simulations. Remarkably, the simulations gave ample credence to all aspects of our proposal and allowed us to gain insight into the effect of specific glutamate rotamers on single-channel conductance. On the basis of single-channel currents recorded from the muscle nicotinic acetylcholine receptor (AChR), we have recently hypothesized that the conformation adopted by the glutamate side chains at the first turn of the pore-lining α-helices is a key determinant of the rate of ion permeation. In this paper, we set out to test these ideas within a framework of atomic detail and stereochemical rigor by conducting all-atom molecular dynamics and Brownian dynamics simulations on an extensively validated model of the open-channel muscle AChR. Our simulations provided ample support to the notion that the different rotamers of these glutamates partition into two classes that differ markedly in their ability to catalyze ion conduction, and that the conformations of the four wild-type glutamates are such that two of them “fall” in each rotamer class. Moreover, the simulations allowed us to identify the mm (χ1 ≅ –60°; χ2 ≅ –60°) and tp (χ1 ≅ 180°; χ2 ≅ +60°) rotamers as the likely conduction-catalyzing conformations of the AChR’s selectivity-filter glutamates. More generally, our work shows an example of how experimental benchmarks can guide molecular simulations into providing a type of structural and mechanistic insight that seems otherwise unattainable.

[1]  J. Richardson,et al.  Amino acid preferences for specific locations at the ends of alpha helices. , 1988, Science.

[2]  G. Rose,et al.  Helix signals in proteins. , 1988, Science.

[3]  K. Fatima-Shad,et al.  Anion permeation in GABA- and glycine-gated channels of mammalian cultured hippocampal neurons , 1993, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[4]  G. Rose,et al.  Helix stop signals in proteins and peptides: the capping box. , 1993, Biochemistry.

[5]  G. Cymes,et al.  Decremental Response to High-Frequency Trains of Acetylcholine Pulses but Unaltered Fractional Ca2+ Currents in a Panel of “Slow-Channel Syndrome” Nicotinic Receptor Mutants , 2009, The Journal of general physiology.

[6]  R. MacKinnon,et al.  Two identical noninteracting sites in an ion channel revealed by proton transfer. , 1994, Science.

[7]  B. Sakmann,et al.  Rings of negatively charged amino acids determine the acetylcholine receptor channel conductance , 1988, Nature.

[8]  P. Argos,et al.  Knowledge‐based protein secondary structure assignment , 1995, Proteins.

[9]  Sunhwan Jo,et al.  Web interface for brownian dynamics simulation of ion transport and its applications to beta‐barrel pores , 2012, J. Comput. Chem..

[10]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[11]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[12]  S. Nair,et al.  Mutations that stabilize the open state of the Erwinia chrisanthemi ligand-gated ion channel fail to change the conformation of the pore domain in crystals , 2012, Proceedings of the National Academy of Sciences.

[13]  B. Wallace,et al.  HOLE: a program for the analysis of the pore dimensions of ion channel structural models. , 1996, Journal of molecular graphics.

[14]  P. Emsley,et al.  Features and development of Coot , 2010, Acta crystallographica. Section D, Biological crystallography.

[15]  Michael P. Eastwood,et al.  Atomic-level simulation of current–voltage relationships in single-file ion channels , 2013, The Journal of general physiology.

[16]  V. Bondarenko,et al.  NMR structures of the transmembrane domains of the α4β2 nAChR. , 2012, Biochimica et biophysica acta.

[17]  U. Sauer,et al.  Dissection of helix capping in T4 lysozyme by structural and thermodynamic analysis of six amino acid substitutions at Thr 59. , 1992, Biochemistry.

[18]  R. Dutzler,et al.  X-ray structure of a prokaryotic pentameric ligand-gated ion channel , 2008, Nature.

[19]  A. Doig,et al.  Side-chain structures in the first turn of the alpha-helix. , 1999, Journal of molecular biology.

[20]  G. Cymes,et al.  The unanticipated complexity of the selectivity-filter glutamates of nicotinic receptors , 2012, Nature chemical biology.

[21]  S. Dasgupta,et al.  Design of helix ends. Amino acid preferences, hydrogen bonding and electrostatic interactions. , 2009, International journal of peptide and protein research.

[22]  P. Taylor,et al.  An Ion Selectivity Filter in the Extracellular Domain of Cys-loop Receptors Reveals Determinants for Ion Conductance* , 2008, Journal of Biological Chemistry.

[23]  W. Im,et al.  A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. , 2000, Biophysical journal.

[24]  R. Dutzler,et al.  Structure of a potentially open state of a proton-activated pentameric ligand-gated ion channel , 2009, Nature.

[25]  J. Richardson,et al.  Corrections: Amino Acid Preferences for Specific Locations at the Ends of α Helices , 1988 .

[26]  G. Cymes,et al.  Probing ion-channel pores one proton at a time , 2005, Nature.

[27]  B. Sakmann,et al.  Mechanism of anion permeation through channels gated by glycine and gamma‐aminobutyric acid in mouse cultured spinal neurones. , 1987, The Journal of physiology.

[28]  J. Richardson,et al.  The penultimate rotamer library , 2000, Proteins.

[29]  Y. Fujiyoshi,et al.  Gating Movement of Acetylcholine Receptor Caught by Plunge-Freezing , 2012, Journal of molecular biology.

[30]  D. Bertrand,et al.  Pentameric ligand-gated ion channel ELIC is activated by GABA and modulated by benzodiazepines , 2012, Proceedings of the National Academy of Sciences.

[31]  Toby W Allen,et al.  Ion conduction and conformational flexibility of a bacterial voltage-gated sodium channel , 2014, Proceedings of the National Academy of Sciences.

[32]  Alexander D. MacKerell,et al.  Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. , 2012, Journal of chemical theory and computation.

[33]  R. Tsien,et al.  Molecular basis of proton block of L-type Ca2+ channels , 1996, The Journal of general physiology.

[34]  R. Aurora,et al.  Helix capping , 1998, Protein science : a publication of the Protein Society.

[35]  J. Yang,et al.  Ion permeation through 5-hydroxytryptamine-gated channels in neuroblastoma N18 cells , 1990, The Journal of general physiology.

[36]  Roland L. Dunbrack,et al.  proteins STRUCTURE O FUNCTION O BIOINFORMATICS Improved prediction of protein side-chain conformations with SCWRL4 , 2022 .

[37]  M. Jansen,et al.  Experimental determination of the vertical alignment between the second and third transmembrane segments of muscle nicotinic acetylcholine receptors , 2013, Journal of neurochemistry.

[38]  J. A. Peters,et al.  Novel structural determinants of single-channel conductance in nicotinic acetylcholine and 5-hydroxytryptamine type-3 receptors. , 2006, Biochemical Society transactions.

[39]  L. Presta,et al.  Contribution to global protein stabilization of the N-capping box in human growth hormone. , 1994, Biochemistry.

[40]  Eric Gouaux,et al.  Principles of activation and permeation in an anion-selective Cys-loop receptor , 2011, Nature.

[41]  Christopher Ing,et al.  Catalysis of Na+ permeation in the bacterial sodium channel NaVAb , 2013, Proceedings of the National Academy of Sciences.

[42]  G. Hummer,et al.  Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides. , 2009, The journal of physical chemistry. B.

[43]  K. A. Jones,et al.  Functional properties of a cloned 5‐hydroxytryptamine ionotropic receptor subunit: comparison with native mouse receptors. , 1994, The Journal of physiology.

[44]  B. Roux The membrane potential and its representation by a constant electric field in computer simulations. , 2008, Biophysical journal.

[45]  A. R. Aricescu,et al.  Crystal structure of a human GABAA receptor , 2014, Nature.

[46]  Uwe Sauer,et al.  Dissection of helix capping in T4 lysozyme by structural and thermodynamic analysis of six amino acid substitutions at Thr 59 , 1993 .

[47]  George D. Rose,et al.  Sequence determinants of the capping box, a stabilizing motif at the N‐termini of α‐helices , 1994 .

[48]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[49]  T. Dwyer,et al.  The permeability of the endplate channel to organic cations in frog muscle , 1980, The Journal of general physiology.

[50]  Leonardo G. Trabuco,et al.  Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics. , 2008, Structure.

[51]  J. Changeux,et al.  X-ray structure of a pentameric ligand-gated ion channel in an apparently open conformation , 2009, Nature.