1
|
Emigh Cortez AM, DeMarco KR, Furutani K, Bekker S, Sack JT, Wulff H, Clancy CE, Vorobyov I, Yarov-Yarovoy V. Structural modeling of hERG channel-drug interactions using Rosetta. Front Pharmacol 2023; 14:1244166. [PMID: 38035013 PMCID: PMC10682396 DOI: 10.3389/fphar.2023.1244166] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2023] [Accepted: 10/23/2023] [Indexed: 12/02/2023] Open
Abstract
The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.
Collapse
Affiliation(s)
- Aiyana M. Emigh Cortez
- Biophysics Graduate Group, University of California, Davis, Davis, CA, United States
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
| | - Kevin R. DeMarco
- Biophysics Graduate Group, University of California, Davis, Davis, CA, United States
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
| | - Kazuharu Furutani
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- Department of Pharmacology, Tokushima Bunri University, Tokushima, Japan
| | - Slava Bekker
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- American River College, Sacramento, CA, United States
| | - Jon T. Sack
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- Department of Anesthesiology and Pain Medicine, University of California, Davis, Davis, CA, United States
| | - Heike Wulff
- Department of Pharmacology, University of California, Davis, Davis, CA, United States
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- Department of Pharmacology, University of California, Davis, Davis, CA, United States
- Center for Precision Medicine and Data Sciences, University of California, Davis, Davis, CA, United States
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- Department of Pharmacology, University of California, Davis, Davis, CA, United States
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
- Department of Anesthesiology and Pain Medicine, University of California, Davis, Davis, CA, United States
| |
Collapse
|
2
|
Maly J, Emigh AM, DeMarco KR, Furutani K, Sack JT, Clancy CE, Vorobyov I, Yarov-Yarovoy V. Structural modeling of the hERG potassium channel and associated drug interactions. Front Pharmacol 2022; 13:966463. [PMID: 36188564 PMCID: PMC9523588 DOI: 10.3389/fphar.2022.966463] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Accepted: 08/25/2022] [Indexed: 11/13/2022] Open
Abstract
The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.
Collapse
Affiliation(s)
- Jan Maly
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,Biophysics Graduate Group, University of California, Davis, Davis, CA, United States
| | - Aiyana M. Emigh
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,Biophysics Graduate Group, University of California, Davis, Davis, CA, United States
| | - Kevin R. DeMarco
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,Biophysics Graduate Group, University of California, Davis, Davis, CA, United States
| | - Kazuharu Furutani
- Department of Pharmacology, Tokushima Bunri University, Tokushima, Japan
| | - Jon T. Sack
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,Department of Pharmacology, University of California, Davis, Davis, CA, United States
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,Department of Pharmacology, University of California, Davis, Davis, CA, United States,*Correspondence: Igor Vorobyov, ; Vladimir Yarov-Yarovoy,
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, United States,*Correspondence: Igor Vorobyov, ; Vladimir Yarov-Yarovoy,
| |
Collapse
|
3
|
Kudaibergenova M, Guo J, Khan HM, Zahid F, Lees-Miller J, Noskov SY, Duff HJ. Allosteric Coupling Between Drug Binding and the Aromatic Cassette in the Pore Domain of the hERG1 Channel: Implications for a State-Dependent Blockade. Front Pharmacol 2020; 11:914. [PMID: 32694995 PMCID: PMC7338687 DOI: 10.3389/fphar.2020.00914] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2019] [Accepted: 06/04/2020] [Indexed: 12/18/2022] Open
Abstract
Human-ether-a-go-go-related channel (hERG1) is the pore-forming domain of the delayed rectifier K+ channel in the heart which underlies the IKr current. The channel has been extensively studied due to its propensity to bind chemically diverse group of drugs. The subsequent hERG1 block can lead to a prolongation of the QT interval potentially leading to an abnormal cardiac electrical activity. The recently solved cryo-EM structure featured a striking non-swapped topology of the Voltage-Sensor Domain (VSD) which is packed against the pore-domain as well as a small and hydrophobic intra-cavity space. The small size and hydrophobicity of the cavity was unexpected and challenges the already-established hypothesis of drugs binding to the wide cavity. Recently, we showed that an amphipathic drug, ivabradine, may favorably bind the channel from the lipid-facing surface and we discovered a mutant (M651T) on the lipid facing domain between the VSD and the PD which inhibited the blocking capacity of the drug. Using multi-microseconds Molecular Dynamics (MD) simulations of wild-type and M651T mutant hERG1, we suggested the block of the channel through the lipid mediated pathway, the opening of which is facilitated by the flexible phenylalanine ring (F656). In this study, we characterize the dynamic interaction of the methionine-aromatic cassette in the S5-S6 helices by combining data from electrophysiological experiments with MD simulations and molecular docking to elucidate the complex allosteric coupling between drug binding to lipid-facing and intra-cavity sites and aromatic cassette dynamics. We investigated two well-established hERG1 blockers (ivabradine and dofetilide) for M651 sensitivity through electrophysiology and mutagenesis techniques. Our electrophysiology data reveal insensitivity of dofetilide to the mutations at site M651 on the lipid facing side of the channel, mirroring our results obtained from docking experiments. Moreover, we show that the dofetilide-induced block of hERG1 occurs through the intracellular space, whereas little to no block of ivabradine is observed during the intracellular application of the drug. The dynamic conformational rearrangement of the F656 appears to regulate the translocation of ivabradine into the central cavity. M651T mutation appears to disrupt this entry pathway by altering the molecular conformation of F656.
Collapse
Affiliation(s)
- Meruyert Kudaibergenova
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada.,Cumming School of Medicine, Libin Cardiovascular Institute of Alberta, University of Calgary, Calgary, AB, Canada
| | - Jiqing Guo
- Cumming School of Medicine, Libin Cardiovascular Institute of Alberta, University of Calgary, Calgary, AB, Canada
| | - Hanif M Khan
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Farhan Zahid
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - James Lees-Miller
- Cumming School of Medicine, Libin Cardiovascular Institute of Alberta, University of Calgary, Calgary, AB, Canada
| | - Sergei Yu Noskov
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Henry J Duff
- Cumming School of Medicine, Libin Cardiovascular Institute of Alberta, University of Calgary, Calgary, AB, Canada
| |
Collapse
|
4
|
Abstract
The human voltage-sensitive K+ channel hERG plays a fundamental role in cardiac action potential repolarization, effectively controlling the QT interval of the electrocardiogram. Inherited loss- or gain-of-function mutations in hERG can result in dangerous “long” (LQTS) or “short” QT syndromes (SQTS), respectively, and the anomalous susceptibility of hERG to block by a diverse range of drugs underlies an acquired LQTS. A recent open channel cryo-EM structure of hERG should greatly advance understanding of the molecular basis of hERG channelopathies and drug-induced LQTS. Here we describe an update of recent research that addresses the nature of the particular gated state of hERG captured in the new structure, and the insight afforded by the structure into the molecular basis for high affinity drug block of hERG, the binding of hERG activators and the molecular basis of hERG's peculiar gating properties. Interpretation of the pharmacology of natural SQTS mutants in the context of the structure is a promising approach to understanding the molecular basis of hERG inactivation, and the structure suggests how voltage-dependent changes in the membrane domain may be transmitted to an extracellular “turret” to effect inactivation through aromatic side chain motifs that are conserved throughout the KCNH family of channels.
Collapse
Affiliation(s)
- Andrew Butler
- School of Physiology, Pharmacology and Neuroscience, Medical Sciences Building, Bristol, United Kingdom
| | - Matthew V Helliwell
- School of Physiology, Pharmacology and Neuroscience, Medical Sciences Building, Bristol, United Kingdom
| | - Yihong Zhang
- School of Physiology, Pharmacology and Neuroscience, Medical Sciences Building, Bristol, United Kingdom
| | - Jules C Hancox
- School of Physiology, Pharmacology and Neuroscience, Medical Sciences Building, Bristol, United Kingdom
| | | |
Collapse
|
5
|
Dutta S, Chang KC, Beattie KA, Sheng J, Tran PN, Wu WW, Wu M, Strauss DG, Colatsky T, Li Z. Corrigendum: Optimization of an In silico Cardiac Cell Model for Proarrhythmia Risk Assessment. Front Physiol 2017; 8:1025. [PMID: 29230183 PMCID: PMC5723676 DOI: 10.3389/fphys.2017.01025] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2017] [Accepted: 11/27/2017] [Indexed: 11/30/2022] Open
Affiliation(s)
- Sara Dutta
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Kelly C Chang
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Kylie A Beattie
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Jiansong Sheng
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Phu N Tran
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Wendy W Wu
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Min Wu
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - David G Strauss
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | - Thomas Colatsky
- Marshview Life Science Advisors, Seabrook Island, SC, United States
| | - Zhihua Li
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, MD, United States
| |
Collapse
|
6
|
Dutta S, Chang KC, Beattie KA, Sheng J, Tran PN, Wu WW, Wu M, Strauss DG, Colatsky T, Li Z. Optimization of an In silico Cardiac Cell Model for Proarrhythmia Risk Assessment. Front Physiol 2017; 8:616. [PMID: 28878692 PMCID: PMC5572155 DOI: 10.3389/fphys.2017.00616] [Citation(s) in RCA: 100] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Accepted: 08/09/2017] [Indexed: 12/19/2022] Open
Abstract
Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive in vitro Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using in silico models and in vitro multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model). We also developed a novel metric that could separate drugs with different TdP liabilities at high concentrations based on total electronic charge carried by the major inward ionic currents during the action potential. In this study, we further optimized the IKr-dynamic ORd model by refining model parameters using published human cardiomyocyte experimental data under control and drug block conditions. Using this optimized model and manual patch clamp data, we developed an updated version of the metric that quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state action potential, which could classify the level of drug-induced TdP risk across a wide range of concentrations and pacing rates. We also established a framework to quantitatively evaluate a system's robustness against the induction of early afterdepolarizations (EADs), and demonstrated that the new metric is correlated with the cell's robustness to the pro-EAD perturbation of IKr conductance reduction. In summary, in this work we present an optimized model that is more consistent with experimental data, an improved metric that can classify drugs at concentrations both near and higher than clinical exposure, and a physiological framework to check the relationship between a metric and EAD. These findings provide a solid foundation for using in silico models for the regulatory assessment of TdP risk under the CiPA paradigm.
Collapse
Affiliation(s)
- Sara Dutta
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Kelly C Chang
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Kylie A Beattie
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Jiansong Sheng
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Phu N Tran
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Wendy W Wu
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Min Wu
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - David G Strauss
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Thomas Colatsky
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| | - Zhihua Li
- Division of Applied Regulatory Science, Office of Clinical Pharmacology, Office of Translational Sciences, Center for Drug Evaluation and Research, U.S. Food and Drug AdministrationSilver Spring, MD, United States
| |
Collapse
|