Simulating ion channel activation mechanisms using swarms of trajectories

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

[1]  W. Catterall,et al.  THE CRYSTAL STRUCTURE OF A VOLTAGE-GATED SODIUM CHANNEL , 2011, Nature.

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

[3]  Gerhard Hummer,et al.  Characterization of a dynamic string method for the construction of transition pathways in molecular reactions. , 2012, The journal of physical chemistry. B.

[4]  K. Satyshur,et al.  Site-Directed Spin Labeling Reveals Pentameric Ligand-Gated Ion Channel Gating Motions , 2012, PLoS biology.

[5]  Harish Vashisth,et al.  "DFG-flip" in the insulin receptor kinase is facilitated by a helical intermediate state of the activation loop. , 2011, Biophysical journal.

[6]  J. Changeux,et al.  Conformational Dynamics in a Nicotinic Receptor Homologue Probed by Simulations , 2011 .

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

[8]  J. Straub,et al.  The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature , 1997 .

[9]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[10]  J. Galpin,et al.  Contributions of Conserved Residues at the Gating Interface of Glycine Receptors* , 2011, The Journal of Biological Chemistry.

[11]  N Le Novère,et al.  Conformational spread in a ring of proteins: a stochastic approach to allostery. , 2001, Journal of molecular biology.

[12]  Massimo Marchi,et al.  Adiabatic bias molecular dynamics: A method to navigate the conformational space of complex molecular systems , 1999 .

[13]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[14]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[15]  Marco Beato,et al.  Single-Channel Behavior of Heteromeric α1β Glycine Receptors: An Attempt to Detect a Conformational Change before the Channel Opens , 2004, The Journal of Neuroscience.

[16]  Dimitar V. Pachov,et al.  Free energy landscape of activation in a signaling protein at atomic resolution , 2015, Nature Communications.

[17]  I. Vorobyov,et al.  Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation. , 2019, Chemical reviews.

[18]  F. Findeisen,et al.  Structure of a prokaryotic sodium channel pore reveals essential gating elements and an outer ion binding site common to eukaryotic channels. , 2014, Journal of molecular biology.

[19]  Simone Marsili,et al.  Self-healing umbrella sampling: a non-equilibrium approach for quantitative free energy calculations. , 2006, The journal of physical chemistry. B.

[20]  David N. LeBard,et al.  Multiple binding sites for the general anesthetic isoflurane identified in the nicotinic acetylcholine receptor transmembrane domain , 2010, Proceedings of the National Academy of Sciences.

[21]  Jason M. Swails,et al.  Constant pH Replica Exchange Molecular Dynamics in Explicit Solvent Using Discrete Protonation States: Implementation, Testing, and Validation , 2014, Journal of chemical theory and computation.

[22]  C. Brooks,et al.  Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism , 2014, Proteins.

[23]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[24]  A. R. Srinivasan,et al.  Quasi‐harmonic method for studying very low frequency modes in proteins , 1984, Biopolymers.

[25]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[26]  A. Szabó,et al.  Electron transfer reaction dynamics in non-Debye solvents , 1998 .

[27]  B. Roux,et al.  Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data. , 2003, Journal of the American Chemical Society.

[28]  J. Changeux Allostery and the Monod-Wyman-Changeux model after 50 years. , 2012, Annual review of biophysics.

[29]  Francisco Bezanilla,et al.  Intermediate state trapping of a voltage sensor , 2012, The Journal of general physiology.

[30]  P. Biggin,et al.  Agonist and antagonist binding in human glycine receptors. , 2014, Biochemistry.

[31]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[32]  B. Wallace,et al.  The complete structure of an activated open sodium channel , 2017, Nature Communications.

[33]  M. Delarue,et al.  Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the Gloeobacter violaceus pentameric ligand-gated ion channel , 2017, PLoS biology.

[34]  Albert C. Pan,et al.  Finding transition pathways using the string method with swarms of trajectories. , 2008, The journal of physical chemistry. B.

[35]  Andrew Jenkins,et al.  General Anesthetic Actions on GABAA Receptors , 2010, Current neuropharmacology.

[36]  J. Changeux,et al.  Crystal structures of a pentameric ligand-gated ion channel provide a mechanism for activation , 2013, Proceedings of the National Academy of Sciences.

[37]  Albert C. Pan,et al.  Atomic-level characterization of protein–protein association , 2018, Proceedings of the National Academy of Sciences.

[38]  J. Changeux,et al.  X-ray structures of general anaesthetics bound to a pentameric ligand-gated ion channel , 2011, Nature.

[39]  J. Changeux,et al.  Molecular tuning of fast gating in pentameric ligand-gated ion channels. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[40]  B. Roux,et al.  A molecular dynamics study of gating in dioxolane-linked gramicidin A channels. , 1994, Biophysical journal.

[41]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[42]  E. Vanden-Eijnden,et al.  Solvent coarse-graining and the string method applied to the hydrophobic collapse of a hydrated chain , 2007, Proceedings of the National Academy of Sciences.

[43]  R. Dutzler,et al.  A prokaryotic perspective on pentameric ligand-gated ion channel structure. , 2009, Current opinion in structural biology.

[44]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[45]  Christophe Chipot,et al.  Exploring the free-energy landscape of a short peptide using an average force. , 2005, The Journal of chemical physics.

[46]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[47]  Benoît Roux,et al.  Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer. , 2012, Journal of chemical theory and computation.

[48]  N. Yan,et al.  Structures of human Nav1.7 channel in complex with auxiliary subunits and animal toxins , 2019, Science.

[49]  Eric Darve,et al.  Calculating Free Energies Using a Scaled-Force Molecular Dynamics Algorithm , 2002 .

[50]  S. Sine,et al.  Principal pathway coupling agonist binding to channel gating in nicotinic receptors , 2005, Nature.

[51]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[52]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[53]  Erik Lindahl,et al.  Conformational gating dynamics in the GluCl anion-selective chloride channel. , 2015, ACS chemical neuroscience.

[54]  Henry A. Lester,et al.  Cis–trans isomerization at a proline opens the pore of a neurotransmitter-gated ion channel , 2005, Nature.

[55]  M. Delarue,et al.  Sites of Anesthetic Inhibitory Action on a Cationic Ligand-Gated Ion Channel. , 2016, Structure.

[56]  A. Auerbach,et al.  Acetylcholine Receptor Gating at Extracellular Transmembrane Domain Interface: the “Pre-M1” Linker , 2007, The Journal of general physiology.

[57]  Sunhwan Jo,et al.  Cholesterol flip-flop: insights from free energy simulation studies. , 2010, The journal of physical chemistry. B.

[58]  Fabio Pietrucci,et al.  Strategies for the exploration of free energy landscapes: Unity in diversity and challenges ahead , 2017 .

[59]  P. Koehl,et al.  Electrostatics, proton sensor, and networks governing the gating transition in GLIC, a proton-gated pentameric ion channel , 2018, Proceedings of the National Academy of Sciences.

[60]  Emanuele Paci,et al.  Boxed molecular dynamics: a simple and general technique for accelerating rare event kinetics and mapping free energy in large molecular systems. , 2009, The journal of physical chemistry. B.

[61]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[62]  R. Dutzler,et al.  Signal Transduction at the Domain Interface of Prokaryotic Pentameric Ligand-Gated Ion Channels , 2016, PLoS biology.

[63]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[64]  Antoine Taly,et al.  A gating mechanism of pentameric ligand-gated ion channels , 2013, Proceedings of the National Academy of Sciences.

[65]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[66]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[67]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[68]  J. P. Grossman,et al.  Biomolecular simulation: a computational microscope for molecular biology. , 2012, Annual review of biophysics.

[69]  G. Hummer,et al.  Classical Molecular Dynamics with Mobile Protons , 2017, J. Chem. Inf. Model..

[70]  Vijay S Pande,et al.  Improvements in Markov State Model Construction Reveal Many Non-Native Interactions in the Folding of NTL9. , 2013, Journal of chemical theory and computation.

[71]  Spencer T. Stober,et al.  Energetics and mechanism of the normal-to-amyloidogenic isomerization of β2-microglobulin: on-the-fly string method calculations. , 2012, The journal of physical chemistry. B.

[72]  Ron O. Dror,et al.  Mechanism of Voltage Gating in Potassium Channels , 2012, Science.

[73]  C. Kindler,et al.  Anaesthetic mechanisms: update on the challenge of unravelling the mystery of anaesthesia , 2009, European journal of anaesthesiology.

[74]  Fu Kit Sheong,et al.  TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromolecules. , 2019, The Journal of chemical physics.

[75]  T. Tomizaki,et al.  X-ray structure of the mouse serotonin 5-HT3 receptor , 2014, Nature.

[76]  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.

[77]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[78]  Jianping Wu,et al.  Structure of a eukaryotic voltage-gated sodium channel at near-atomic resolution , 2017, Science.

[79]  Patrick G. Blachly,et al.  Measuring the successes and deficiencies of constant pH molecular dynamics: A blind prediction study , 2011, Proteins.

[80]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[81]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[82]  A. Laio,et al.  Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamics. , 2006, The journal of physical chemistry. B.

[83]  J. Changeux,et al.  A prokaryotic proton-gated ion channel from the nicotinic acetylcholine receptor family , 2007, Nature.

[84]  Klaus Schulten,et al.  Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems. , 2017, Journal of chemical theory and computation.

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

[86]  Schuster,et al.  Separation of a mixture of independent signals using time delayed correlations. , 1994, Physical review letters.

[87]  Jianlin Lei,et al.  Structure of the Nav1.4-β1 Complex from Electric Eel , 2017, Cell.

[88]  Nela Durisic,et al.  Probing the Structural Mechanism of Partial Agonism in Glycine Receptors Using the Fluorescent Artificial Amino Acid, ANAP. , 2017, ACS chemical biology.

[89]  A. Auerbach,et al.  Functional anatomy of an allosteric protein , 2013, Nature Communications.

[90]  Srirupa Chakraborty,et al.  A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip , 2017, The Journal of general physiology.

[91]  Victor Ovchinnikov,et al.  Analysis and elimination of a bias in targeted molecular dynamics simulations of conformational transitions: application to calmodulin. , 2012, The journal of physical chemistry. B.

[92]  J. Shepherd,et al.  Climate: More ways to govern geoengineering , 2012, Nature.

[93]  F. Noé,et al.  Transition networks for modeling the kinetics of conformational change in macromolecules. , 2008, Current opinion in structural biology.

[94]  Gerhard Hummer,et al.  Position-dependent diffusion coefficients and free energies from Bayesian analysis of equilibrium and replica molecular dynamics simulations , 2005 .

[95]  Samuel F. Berkovic,et al.  Mechanisms of human inherited epilepsies , 2009, Progress in Neurobiology.

[96]  W. Catterall,et al.  Structures of closed and open states of a voltage-gated sodium channel , 2017, Proceedings of the National Academy of Sciences.

[97]  Eric Vanden-Eijnden,et al.  Free energy of conformational transition paths in biomolecules: the string method and its application to myosin VI. , 2011, The Journal of chemical physics.

[98]  T. Harpole,et al.  Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations. , 2017, Biochimica et biophysica acta. Biomembranes.

[99]  Jana K. Shen,et al.  All-Atom Continuous Constant pH Molecular Dynamics With Particle Mesh Ewald and Titratable Water. , 2016, Journal of chemical theory and computation.

[100]  Andrew E. Torda,et al.  Local elevation: A method for improving the searching properties of molecular dynamics simulation , 1994, J. Comput. Aided Mol. Des..

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

[102]  E. Tajkhorshid,et al.  Atomic-level characterization of transport cycle thermodynamics in the glycerol-3-phosphate:phosphate antiporter , 2015, Nature Communications.

[103]  Frank Noé,et al.  Markov state models of biomolecular conformational dynamics. , 2014, Current opinion in structural biology.

[104]  L. Sivilotti What single‐channel analysis tells us of the activation mechanism of ligand‐gated channels: the case of the glycine receptor , 2010, The Journal of physiology.

[105]  W. Im,et al.  Theoretical and computational models of biological ion channels , 2004, Quarterly Reviews of Biophysics.

[106]  Stephan A. Pless,et al.  Magnitude of a Conformational Change in the Glycine Receptor β1-β2 Loop Is Correlated with Agonist Efficacy* , 2009, The Journal of Biological Chemistry.

[107]  Anthony Auerbach,et al.  A stepwise mechanism for acetylcholine receptor channel gating , 2007, Nature.

[108]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[109]  Eric Darve,et al.  Assessing the efficiency of free energy calculation methods. , 2004, The Journal of chemical physics.

[110]  Sonya M. Hanson,et al.  Packing of the Extracellular Domain Hydrophobic Core Has Evolved to Facilitate Pentameric Ligand-gated Ion Channel Function , 2010, The Journal of Biological Chemistry.

[111]  Hiroshi Fujisaki,et al.  Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase , 2012, PLoS Comput. Biol..

[112]  M. Prevost,et al.  Functional prokaryotic–eukaryotic chimera from the pentameric ligand-gated ion channel family , 2011, Proceedings of the National Academy of Sciences.

[113]  M. Delarue,et al.  Structural basis for potentiation by alcohols and anaesthetics in a ligand-gated ion channel , 2013, Nature Communications.

[114]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[115]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[116]  J. Changeux,et al.  Chimaeric nicotinic–serotonergic receptor combines distinct ligand binding and channel specificities , 1993, Nature.

[117]  A. Voter A method for accelerating the molecular dynamics simulation of infrequent events , 1997 .

[118]  E. Gouaux,et al.  Glycine receptor mechanism elucidated by electron cryo-microscopy , 2015, Nature.

[119]  Xin Huang,et al.  Crystal structure of human glycine receptor-α3 bound to antagonist strychnine , 2015, Nature.

[120]  B. Wallace,et al.  Structure of a bacterial voltage-gated sodium channel pore reveals mechanisms of opening and closing , 2012, Nature Communications.

[121]  Gerhard Hummer,et al.  Pore opening and closing of a pentameric ligand-gated ion channel , 2010, Proceedings of the National Academy of Sciences.

[122]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[123]  Jean-Pierre Changeux,et al.  Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics , 2017, PLoS Comput. Biol..

[124]  V. Pande,et al.  Markov State Models: From an Art to a Science. , 2018, Journal of the American Chemical Society.

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

[126]  M. Delarue,et al.  String method solution of the gating pathways for a pentameric ligand-gated ion channel , 2017, Proceedings of the National Academy of Sciences.

[127]  Surajit Banerjee,et al.  X-ray structures of GluCl in apo states reveal a gating mechanism of Cys-loop receptors , 2014, Nature.

[128]  Christof Schütte,et al.  Metastability and Markov State Models in Molecular Dynamics Modeling, Analysis , 2016 .

[129]  J. Changeux,et al.  Structure and pharmacology of pentameric receptor channels: from bacteria to brain. , 2012, Structure.

[130]  C. Domene,et al.  Capturing the Molecular Mechanism of Anesthetic Action by Simulation Methods. , 2018, Chemical reviews.

[131]  Marc Baaden,et al.  One-microsecond molecular dynamics simulation of channel gating in a nicotinic receptor homologue , 2010, Proceedings of the National Academy of Sciences.

[132]  Christoph Dellago,et al.  On the calculation of reaction rate constants in the transition path ensemble , 1999 .

[133]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

[134]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[135]  Modesto Orozco,et al.  Prediction and validation of protein intermediate states from structurally rich ensembles and coarse-grained simulations , 2016, Nature Communications.

[136]  David Colquhoun,et al.  On the nature of partial agonism in the nicotinic receptor superfamily , 2008, Nature.

[137]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[138]  Frank Noé,et al.  An Introduction to Markov State Models and Their Application to Long Timescale Molecular Simulation , 2014, Advances in Experimental Medicine and Biology.

[139]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[140]  M. Delarue,et al.  A locally closed conformation of a bacterial pentameric proton-gated ion channel , 2012, Nature Structural &Molecular Biology.

[141]  Toni Giorgino,et al.  Identification of slow molecular order parameters for Markov model construction. , 2013, The Journal of chemical physics.

[142]  D. Mowrey,et al.  Conformational Changes Underlying Desensitization of the Pentameric Ligand-Gated Ion Channel ELIC. , 2015, Structure.