Identification of the Molecular Site of Ivabradine Binding to HCN4 Channels

Ivabradine is a specific heart rate-reducing agent approved as a treatment of chronic stable angina. Its mode of action involves a selective and specific block of HCN channels, the molecular components of sinoatrial "funny" (f)-channels. Different studies suggest that the binding site of ivabradine is located in the inner vestibule of HCN channels, but the molecular details of ivabradine binding are unknown. We thus sought to investigate by mutagenesis and in silico analysis which residues of the HCN4 channel, the HCN isoform expressed in the sinoatrial node, are involved in the binding of ivabradine. Using homology modeling, we verified the presence of an inner cavity below the channel pore and identified residues lining the cavity; these residues were replaced with alanine (or valine) either alone or in combination, and WT and mutant channels were expressed in HEK293 cells. Comparison of the block efficiency of mutant vs WT channels, measured by patch-clamp, revealed that residues Y506, F509 and I510 are involved in ivabradine binding. For each mutant channel, docking simulations correctly explain the reduced block efficiency in terms of proportionally reduced affinity for ivabradine binding. In summary our study shows that ivabradine occupies a cavity below the channel pore, and identifies specific residues facing this cavity that interact and stabilize the ivabradine molecule. This study provides an interpretation of known properties of f/HCN4 channel block by ivabradine such as the “open channel block”, the current-dependence of block and the property of "trapping" of drug molecules in the closed configuration.

[1]  Martin Biel,et al.  Exploring HCN channels as novel drug targets , 2011, Nature Reviews Drug Discovery.

[2]  S. Franceschetti,et al.  Recessive Loss-of-Function Mutation in the Pacemaker HCN2 Channel Causing Increased Neuronal Excitability in a Patient with Idiopathic Generalized Epilepsy , 2011, The Journal of Neuroscience.

[3]  P. McNaughton,et al.  HCN2 Ion Channels Play a Central Role in Inflammatory and Neuropathic Pain , 2011, Science.

[4]  Benoît Roux,et al.  Structural basis for the coupling between activation and inactivation gates in K+ channels , 2010, Nature.

[5]  J. Vilaine,et al.  I(f) inhibition in cardiovascular diseases. , 2010, Advances in pharmacology.

[6]  C. Lau,et al.  Probing the bradycardic drug binding receptor of HCN-encoded pacemaker channels , 2009, Pflügers Archiv - European Journal of Physiology.

[7]  C. Wahl-Schott,et al.  Hyperpolarization-activated cation channels: from genes to function. , 2009, Physiological reviews.

[8]  Krista I Kinard,et al.  Molecular Mapping of the Binding Site for a Blocker of Hyperpolarization-Activated, Cyclic Nucleotide-Modulated Pacemaker Channels , 2007, Journal of Pharmacology and Experimental Therapeutics.

[9]  M. Sutcliffe,et al.  Drug block of the hERG potassium channel: Insight from modeling , 2007, Proteins.

[10]  A. Sali,et al.  Statistical potential for assessment and prediction of protein structures , 2006, Protein science : a publication of the Protein Society.

[11]  D. DiFrancesco,et al.  Properties of ivabradine‐induced block of HCN1 and HCN4 pacemaker channels , 2006, The Journal of physiology.

[12]  E. Campbell,et al.  Crystal Structure of a Mammalian Voltage-Dependent Shaker Family K+ Channel , 2005, Science.

[13]  V. Torre,et al.  A homology model of the pore region of HCN channels. , 2005, Biophysical journal.

[14]  Gregory W. Kauffman,et al.  Physicochemical Features of the hERG Channel Drug Binding Site* , 2004, Journal of Biological Chemistry.

[15]  S. Siegelbaum,et al.  Hyperpolarization-activated cation currents: from molecules to physiological function. , 2003, Annual review of physiology.

[16]  G. Yellen,et al.  Movements near the Gate of a Hyperpolarization-activated Cation Channel , 2003, The Journal of general physiology.

[17]  Michael C Sanguinetti,et al.  Voltage-dependent profile of human ether-a-go-go-related gene channel block is influenced by a single residue in the S6 transmembrane domain. , 2003, Molecular pharmacology.

[18]  V. Torre,et al.  Pore topology of the hyperpolarization-activated cyclic nucleotide-gated channel from sea urchin sperm. , 2002, Biophysical journal.

[19]  D. DiFrancesco,et al.  Current-dependent Block of Rabbit Sino-Atrial Node If Channels by Ivabradine , 2002, The Journal of general physiology.

[20]  P. Phale,et al.  Voltage-Controlled Gating at the Intracellular Entrance to a Hyperpolarization-Activated Cation Channel , 2002, The Journal of general physiology.

[21]  István Simon,et al.  The HMMTOP transmembrane topology prediction server , 2001, Bioinform..

[22]  Dario DiFrancesco,et al.  Integrated Allosteric Model of Voltage Gating of Hcn Channels , 2001, The Journal of general physiology.

[23]  G. Yellen,et al.  Blocker State Dependence and Trapping in Hyperpolarization-Activated Cation Channels , 2001, The Journal of general physiology.

[24]  A. Krogh,et al.  Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. , 2001, Journal of molecular biology.

[25]  Jun Chen,et al.  A structural basis for drug-induced long QT syndrome. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[26]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[27]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[28]  G. Murshudov,et al.  Refinement of macromolecular structures by the maximum-likelihood method. , 1997, Acta crystallographica. Section D, Biological crystallography.

[29]  J. Lenfant,et al.  Mode of action of bradycardic agent, S 16257, on ionic currents of rabbit sinoatrial node cells , 1996, British journal of pharmacology.

[30]  Manuel G. Claros,et al.  TopPred II: an improved software for membrane protein structure predictions , 1994, Comput. Appl. Biosci..

[31]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[32]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[33]  D S Moss,et al.  Main-chain bond lengths and bond angles in protein structures. , 1993, Journal of molecular biology.