1
|
Arab I, Laukens K, Bittremieux W. Semisupervised Learning to Boost hERG, Nav1.5, and Cav1.2 Cardiac Ion Channel Toxicity Prediction by Mining a Large Unlabeled Small Molecule Data Set. J Chem Inf Model 2024; 64:6410-6420. [PMID: 39110924 DOI: 10.1021/acs.jcim.4c01102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/27/2024]
Abstract
Predicting drug toxicity is a critical aspect of ensuring patient safety during the drug design process. Although conventional machine learning techniques have shown some success in this field, the scarcity of annotated toxicity data poses a significant challenge in enhancing models' performance. In this study, we explore the potential of leveraging large unlabeled small molecule data sets using semisupervised learning to improve drug cardiotoxicity predictive performance across three cardiac ion channel targets: the voltage-gated potassium channel (hERG), the voltage-gated sodium channel (Nav1.5), and the voltage-gated calcium channel (Cav1.2). We extensively mined the ChEMBL database, comprising approximately 2 million small molecules, and then employed semisupervised learning to construct robust classification models for this purpose. We achieved a performance boost on highly diverse (i.e., structurally dissimilar) test data sets across all three targets. Using our built models, we screened the whole ChEMBL database and a large set of FDA-approved drugs, identifying several compounds with potential cardiac ion channel activity. To ensure broad accessibility and usability for both technical and nontechnical users, we developed a cross-platform graphical user interface that allows users to make predictions and gain insights into the cardiotoxicity of drugs and other small molecules. The software is made available as open source under the permissive MIT license at https://github.com/issararab/CToxPred2.
Collapse
Affiliation(s)
- Issar Arab
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (biomina), 2020 Antwerp, Belgium
| | - Kris Laukens
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (biomina), 2020 Antwerp, Belgium
| | - Wout Bittremieux
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (biomina), 2020 Antwerp, Belgium
| |
Collapse
|
2
|
Arab I, Egghe K, Laukens K, Chen K, Barakat K, Bittremieux W. Benchmarking of Small Molecule Feature Representations for hERG, Nav1.5, and Cav1.2 Cardiotoxicity Prediction. J Chem Inf Model 2024; 64:2515-2527. [PMID: 37870574 DOI: 10.1021/acs.jcim.3c01301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
In the field of drug discovery, there is a substantial challenge in seeking out chemical structures that possess desirable pharmacological, toxicological, and pharmacokinetic properties. Complications arise when drugs interfere with the functioning of cardiac ion channels, leading to serious cardiovascular consequences. The discontinuation and removal of numerous approved drugs from the market or at late development stages in the pipeline due to such inhibitory effects further highlight the urgency of addressing this issue. Consequently, the early prediction of potential blockers targeting cardiac ion channels during the drug discovery process is of paramount importance. This study introduces a deep learning framework that computationally determines the cardiotoxicity associated with the voltage-gated potassium channel (hERG), the voltage-gated calcium channel (Cav1.2), and the voltage-gated sodium channel (Nav1.5) for drug candidates. The predictive capabilities of three feature representations─molecular fingerprints, descriptors, and graph-based numerical representations─are rigorously benchmarked. Additionally, a novel training and evaluation data set framework is presented, enabling predictive model training of drug off-target cardiotoxicity using a comprehensive and large curated data set covering these three cardiac ion channels. To facilitate these predictions, a robust and comprehensive small molecule cardiotoxicity prediction tool named CToxPred has been developed. It is made available as open source under the permissive MIT license at https://github.com/issararab/CToxPred.
Collapse
Affiliation(s)
- Issar Arab
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (Biomina), 2020 Antwerp, Belgium
| | - Kristof Egghe
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
| | - Kris Laukens
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (Biomina), 2020 Antwerp, Belgium
| | - Ke Chen
- Chair for Theoretical Chemistry, Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Khaled Barakat
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, Alberta 8613, Canada
| | - Wout Bittremieux
- Department of Computer Science, University of Antwerp, 2020 Antwerp, Belgium
- Biomedical Informatics Network Antwerpen (Biomina), 2020 Antwerp, Belgium
| |
Collapse
|
3
|
Lei CL, Whittaker DG, Mirams GR. The impact of uncertainty in hERG binding mechanism on in silico predictions of drug-induced proarrhythmic risk. Br J Pharmacol 2024; 181:987-1004. [PMID: 37740435 DOI: 10.1111/bph.16250] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Revised: 08/23/2023] [Accepted: 08/28/2023] [Indexed: 09/24/2023] Open
Abstract
BACKGROUND AND PURPOSE Drug-induced reduction of the rapid delayed rectifier potassium current carried by the human Ether-à-go-go-Related Gene (hERG) channel is associated with increased risk of arrhythmias. Recent updates to drug safety regulatory guidelines attempt to capture each drug's hERG binding mechanism by combining in vitro assays with in silico simulations. In this study, we investigate the impact on in silico proarrhythmic risk predictions due to uncertainty in the hERG binding mechanism and physiological hERG current model. EXPERIMENTAL APPROACH Possible pharmacological binding models were designed for the hERG channel to account for known and postulated small molecule binding mechanisms. After selecting a subset of plausible binding models for each compound through calibration to available voltage-clamp electrophysiology data, we assessed their effects, and the effects of different physiological models, on proarrhythmic risk predictions. KEY RESULTS For some compounds, multiple binding mechanisms can explain the same data produced under the safety testing guidelines, which results in different inferred binding rates. This can result in substantial uncertainty in the predicted torsade risk, which often spans more than one risk category. By comparison, we found that the effect of a different hERG physiological current model on risk classification was subtle. CONCLUSION AND IMPLICATIONS The approach developed in this study assesses the impact of uncertainty in hERG binding mechanisms on predictions of drug-induced proarrhythmic risk. For some compounds, these results imply the need for additional binding data to decrease uncertainty in safety-critical applications.
Collapse
Affiliation(s)
- Chon Lok Lei
- Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Macau, China
- Department of Biomedical Sciences, Faculty of Health Sciences, University of Macau, Macau, China
| | - Dominic G Whittaker
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| |
Collapse
|
4
|
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: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [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
|
5
|
Furutani K. Facilitation of hERG Activation by Its Blocker: A Mechanism to Reduce Drug-Induced Proarrhythmic Risk. Int J Mol Sci 2023; 24:16261. [PMID: 38003453 PMCID: PMC10671758 DOI: 10.3390/ijms242216261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 11/08/2023] [Accepted: 11/12/2023] [Indexed: 11/26/2023] Open
Abstract
Modulation of the human Ether-à-go-go-Related Gene (hERG) channel, a crucial voltage-gated potassium channel in the repolarization of action potentials in ventricular myocytes of the heart, has significant implications on cardiac electrophysiology and can be either antiarrhythmic or proarrhythmic. For example, hERG channel blockade is a leading cause of long QT syndrome and potentially life-threatening arrhythmias, such as torsades de pointes. Conversely, hERG channel blockade is the mechanism of action of Class III antiarrhythmic agents in terminating ventricular tachycardia and fibrillation. In recent years, it has been recognized that less proarrhythmic hERG blockers with clinical potential or Class III antiarrhythmic agents exhibit, in addition to their hERG-blocking activity, a second action that facilitates the voltage-dependent activation of the hERG channel. This facilitation is believed to reduce the proarrhythmic potential by supporting the final repolarizing of action potentials. This review covers the pharmacological characteristics of hERG blockers/facilitators, the molecular mechanisms underlying facilitation, and their clinical significance, as well as unresolved issues and requirements for research in the fields of ion channel pharmacology and drug-induced arrhythmias.
Collapse
Affiliation(s)
- Kazuharu Furutani
- Department of Pharmacology, Faculty of Pharmaceutical Sciences, Tokushima Bunri University, 180 Nishihama-Boji, Yamashiro-cho, Tokushima 770-8514, Japan
| |
Collapse
|
6
|
Das NR, Sharma T, Toropov AA, Toropova AP, Tripathi MK, Achary PGR. Machine-learning technique, QSAR and molecular dynamics for hERG-drug interactions. J Biomol Struct Dyn 2023; 41:13766-13791. [PMID: 37021352 DOI: 10.1080/07391102.2023.2193641] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 02/06/2023] [Indexed: 04/07/2023]
Abstract
One of the most well-known anti-targets defining medication cardiotoxicity is the voltage-dependent hERG K + channel, which is well-known for its crucial involvement in cardiac action potential repolarization. Torsades de Pointes, QT prolongation, and sudden death are all caused by hERG (the human Ether-à-go-go-Related Gene) inhibition. There is great interest in creating predictive computational (in silico) tools to identify and weed out potential hERG blockers early in the drug discovery process because testing for hERG liability and the traditional experimental screening are complicated, expensive and time-consuming. This study used 2D descriptors of a large curated dataset of 6766 compounds and machine learning approaches to build robust descriptor-based QSAR and predictive classification models for KCNH2 liability. Decision Tree, Random Forest, Logistic Regression, Ada Boosting, kNN, SVM, Naïve Bayes, neural network and stochastic gradient classification classifier algorithms were used to build classification models. If a compound's IC50 value was between 10 μM and less, it was classified as a blocker (hERG-positive), and if it was more, it was classified as a non-blocker (hERG-negative). Matthew's correlation coefficient formula and F1score were applied to compare and track the developed models' performance. Molecular docking and dynamics studies were performed to understand the cardiotoxicity relating to the hERG-gene. The hERG residues interacting after 100 ns are LEU:697, THR:708, PHE:656, HIS:674, HIS:703, TRP:705 and ASN:709 and the hERG-ligand-16 complex trajectory showed stable behaviour with lesser fluctuations in the entire simulation of 200 ns.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Nilima Rani Das
- Department of CA, Siksha 'O' Anusandhan Deemed to be University, Bhubaneswar, Odisha, India
| | - Tripti Sharma
- School of Pharmaceutical Sciences, Siksha 'O' Anusandhan Deemed to be University, Bhubaneswar, Odisha, India
| | - Andrey A Toropov
- Department of Environmental Health Science, Laboratory of Environmental Chemistry and Toxicology, Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Milan, Italy
| | - Alla P Toropova
- Department of Environmental Health Science, Laboratory of Environmental Chemistry and Toxicology, Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Milan, Italy
| | | | - P Ganga Raju Achary
- Department of Chemistry, Siksha 'O' Anusandhan Deemed to be University, Bhubaneswar, Odisha, India
| |
Collapse
|
7
|
Wang T, Sun J, Zhao Q. Investigating cardiotoxicity related with hERG channel blockers using molecular fingerprints and graph attention mechanism. Comput Biol Med 2023; 153:106464. [PMID: 36584603 DOI: 10.1016/j.compbiomed.2022.106464] [Citation(s) in RCA: 112] [Impact Index Per Article: 112.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Revised: 12/12/2022] [Accepted: 12/19/2022] [Indexed: 12/24/2022]
Abstract
Human ether-a-go-go-related gene (hERG) channel blockade by small molecules is a big concern during drug development in the pharmaceutical industry. Failure or inhibition of hERG channel activity caused by drug molecules can lead to prolonging QT interval, which will result in serious cardiotoxicity. Thus, evaluating the hERG blocking activity of all these small molecular compounds is technically challenging, and the relevant procedures are expensive and time-consuming. In this study, we develop a novel deep learning predictive model named DMFGAM for predicting hERG blockers. In order to characterize the molecule more comprehensively, we first consider the fusion of multiple molecular fingerprint features to characterize its final molecular fingerprint features. Then, we use the multi-head attention mechanism to extract the molecular graph features. Both molecular fingerprint features and molecular graph features are fused as the final features of the compounds to make the feature expression of compounds more comprehensive. Finally, the molecules are classified into hERG blockers or hERG non-blockers through the fully connected neural network. We conduct 5-fold cross-validation experiment to evaluate the performance of DMFGAM, and verify the robustness of DMFGAM on external validation datasets. We believe DMFGAM can serve as a powerful tool to predict hERG channel blockers in the early stages of drug discovery and development.
Collapse
Affiliation(s)
- Tianyi Wang
- School of Computer Science and Software Engineering, University of Science and Technology Liaoning, Anshan, 114051, China
| | - Jianqiang Sun
- School of Automation and Electrical Engineering, Linyi University, Linyi, 276000, China
| | - Qi Zhao
- School of Computer Science and Software Engineering, University of Science and Technology Liaoning, Anshan, 114051, China.
| |
Collapse
|
8
|
Goel H, Yu W, MacKerell AD. hERG Blockade Prediction by Combining Site Identification by Ligand Competitive Saturation and Physicochemical Properties. CHEMISTRY (BASEL, SWITZERLAND) 2022; 4:630-646. [PMID: 36712295 PMCID: PMC9881610 DOI: 10.3390/chemistry4030045] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Human ether-a-go-go-related gene (hERG) potassium channel is well-known contributor to drug-induced cardiotoxicity and therefore an extremely important target when performing safety assessments of drug candidates. Ligand-based approaches in connection with quantitative structure active relationships (QSAR) analyses have been developed to predict hERG toxicity. Availability of the recent published cryogenic electron microscopy (cryo-EM) structure for the hERG channel opened the prospect for using structure-based simulation and docking approaches for hERG drug liability predictions. In recent time, the idea of combining structure- and ligand-based approaches for modeling hERG drug liability has gained momentum offering improvements in predictability when compared to ligand-based QSAR practices alone. The present article demonstrates uniting the structure-based SILCS (site-identification by ligand competitive saturation) approach in conjunction with physicochemical properties to develop predictive models for hERG blockade. This combination leads to improved model predictability based on Pearson's R and percent correct (represents rank-ordering of ligands) metric for different validation sets of hERG blockers involving diverse chemical scaffold and wide range of pIC50 values. The inclusion of the SILCS structure-based approach allows determination of the hERG region to which compounds bind and the contribution of different chemical moieties in the compounds to blockade, thereby facilitating the rational ligand design to minimize hERG liability.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn St. Baltimore, MD 21201, United States
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn St. Baltimore, MD 21201, United States
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn St. Baltimore, MD 21201, United States
| |
Collapse
|
9
|
Ding W, Nan Y, Wu J, Han C, Xin X, Li S, Liu H, Zhang L. Combining multi-dimensional molecular fingerprints to predict the hERG cardiotoxicity of compounds. Comput Biol Med 2022; 144:105390. [DOI: 10.1016/j.compbiomed.2022.105390] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Revised: 03/06/2022] [Accepted: 03/07/2022] [Indexed: 01/28/2023]
|
10
|
Combined Pharmacophore and Grid-Independent Molecular Descriptors (GRIND) Analysis to Probe 3D Features of Inositol 1,4,5-Trisphosphate Receptor (IP 3R) Inhibitors in Cancer. Int J Mol Sci 2021; 22:ijms222312993. [PMID: 34884798 PMCID: PMC8657927 DOI: 10.3390/ijms222312993] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 11/18/2021] [Accepted: 11/24/2021] [Indexed: 12/11/2022] Open
Abstract
Inositol 1, 4, 5-trisphosphate receptor (IP3R)-mediated Ca2+ signaling plays a pivotal role in different cellular processes, including cell proliferation and cell death. Remodeling Ca2+ signals by targeting the downstream effectors is considered an important hallmark in cancer progression. Despite recent structural analyses, no binding hypothesis for antagonists within the IP3-binding core (IBC) has been proposed yet. Therefore, to elucidate the 3D structural features of IP3R modulators, we used combined pharmacoinformatic approaches, including ligand-based pharmacophore models and grid-independent molecular descriptor (GRIND)-based models. Our pharmacophore model illuminates the existence of two hydrogen-bond acceptors (2.62 Å and 4.79 Å) and two hydrogen-bond donors (5.56 Å and 7.68 Å), respectively, from a hydrophobic group within the chemical scaffold, which may enhance the liability (IC50) of a compound for IP3R inhibition. Moreover, our GRIND model (PLS: Q2 = 0.70 and R2 = 0.72) further strengthens the identified pharmacophore features of IP3R modulators by probing the presence of complementary hydrogen-bond donor and hydrogen-bond acceptor hotspots at a distance of 7.6-8.0 Å and 6.8-7.2 Å, respectively, from a hydrophobic hotspot at the virtual receptor site (VRS). The identified 3D structural features of IP3R modulators were used to screen (virtual screening) 735,735 compounds from the ChemBridge database, 265,242 compounds from the National Cancer Institute (NCI) database, and 885 natural compounds from the ZINC database. After the application of filters, four compounds from ChemBridge, one compound from ZINC, and three compounds from NCI were shortlisted as potential hits (antagonists) against IP3R. The identified hits could further assist in the design and optimization of lead structures for the targeting and remodeling of Ca2+ signals in cancer.
Collapse
|
11
|
Meng J, Zhang L, Wang L, Li S, Xie D, Zhang Y, Liu H. TSSF-hERG: A machine-learning-based hERG potassium channel-specific scoring function for chemical cardiotoxicity prediction. Toxicology 2021; 464:153018. [PMID: 34757159 DOI: 10.1016/j.tox.2021.153018] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Revised: 10/15/2021] [Accepted: 10/26/2021] [Indexed: 11/27/2022]
Abstract
The human ether-à-go-go-related gene (hERG) encodes the Kv11.1 voltage-gated potassium ion (K+) channel that conducts the rapidly activating delayed rectifier current (IKr) in cardiomyocytes to regulate the repolarization process. Some drugs, as blockers of hERG potassium channels, cannot be marketed due to prolonged QT intervals, as well known as cardiotoxicity. Predetermining the binding affinity values between drugs and hERG through in silico methods can greatly reduce the time and cost required for experimental verification. In this study, we collected 9,215 compounds with AutoDock Vina's docking structures as training set, and collected compounds from four references as test sets. A series of models for predicting the binding affinities of hERG blockers were built based on five machine learning algorithms and combinations of interaction features and ligand features. The model built by support vector regression (SVR) using the combination of all features achieved the best performance on both tenfold cross-validation and external verification, which was selected and named as TSSF-hERG (target-specific scoring function for hERG). TSSF-hERG is more accurate than the classic scoring function of AutoDock Vina and the machine-learning-based generic scoring function RF-Score, with a Pearson's correlation coefficient (Rp) of 0.765, a Spearman's rank correlation coefficient (Rs) of 0.757, a root-mean-square error (RMSE) of 0.585 in a tenfold cross-validation study. All results demonstrated that TSSF-hERG would be useful for improving the power of binding affinity prediction between hERG and compounds, which can be further used for prediction or virtual screening of the hERG-related cardiotoxicity of drug candidates.
Collapse
Affiliation(s)
- Jinhui Meng
- School of Life Science, Liaoning University, Shenyang, 110036, China
| | - Li Zhang
- School of Life Science, Liaoning University, Shenyang, 110036, China; Technology Innovation Center for Computer Simulating and Information Processing of Bio-macromolecules of Liaoning Province, Shenyang, 110036, China; Engineering Laboratory for Molecular Simulation and Designing of Drug Molecules of Liaoning, Shenyang, 110036, China
| | - Lianxin Wang
- School of Life Science, Liaoning University, Shenyang, 110036, China
| | - Shimeng Li
- School of Life Science, Liaoning University, Shenyang, 110036, China
| | - Di Xie
- School of Life Science, Liaoning University, Shenyang, 110036, China
| | - Yuxi Zhang
- School of Life Science, Liaoning University, Shenyang, 110036, China
| | - Hongsheng Liu
- Technology Innovation Center for Computer Simulating and Information Processing of Bio-macromolecules of Liaoning Province, Shenyang, 110036, China; Engineering Laboratory for Molecular Simulation and Designing of Drug Molecules of Liaoning, Shenyang, 110036, China; School of Pharmacy, Liaoning University, Shenyang, 110036, China.
| |
Collapse
|
12
|
Creanza TM, Delre P, Ancona N, Lentini G, Saviano M, Mangiatordi GF. Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study. J Chem Inf Model 2021; 61:4758-4770. [PMID: 34506150 PMCID: PMC9282647 DOI: 10.1021/acs.jcim.1c00744] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Drug-induced blockade of the human
ether-à-go-go-related
gene (hERG) channel is today considered the main
cause of cardiotoxicity in postmarketing surveillance. Hence, several
ligand-based approaches were developed in the last years and are currently
employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates.
Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support
vector machines were applied to integrate docking scores and protein–ligand
interaction fingerprints. A total of 396 models were trained and validated
based on: (i) high-quality experimental bioactivity information returned
by 8337 curated compounds extracted from ChEMBL (version 25) and (ii)
structural predictor data. Molecular docking simulations were performed
using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures
obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and
two published homology models selected for comparison. Interestingly,
some classifiers return performances comparable to ligand-based models
in terms of area under the ROC curve (AUCMAX = 0.86 ±
0.01) and negative predictive values (NPVMAX = 0.81 ±
0.01), thus putting forward the herein proposed computational workflow
as a valuable tool for predicting hERG-related cardiotoxicity
without the limitations of ligand-based models, typically affected
by low interpretability and a limited applicability domain. From a
methodological point of view, our study represents the first example
of a successful integration of docking scores and protein–ligand
interaction fingerprints (IFs) through a support vector machine (SVM)
LASSO regularized strategy. Finally, the study highlights the importance
of using hERG structural models accounting for ligand-induced
fit effects and allowed us to select the best-performing protein conformation
(made available in the Supporting Information, SI) to be employed
for a reliable structure-based prediction of hERG-related cardiotoxicity.
Collapse
Affiliation(s)
- Teresa Maria Creanza
- CNR-Institute of Intelligent Industrial Technologies and Systems for Advanced Manufacturing, Via Amendola 122/o, 70126 Bari, Italy
| | - Pietro Delre
- Chemistry Department, University of Bari "Aldo Moro", via E. Orabona, 4, I-70125 Bari, Italy.,CNR-Institute of Crystallography, Via Amendola 122/o, 70126 Bari, Italy
| | - Nicola Ancona
- CNR-Institute of Intelligent Industrial Technologies and Systems for Advanced Manufacturing, Via Amendola 122/o, 70126 Bari, Italy
| | - Giovanni Lentini
- Department of Pharmacy-Pharmaceutical Sciences, University of Bari "Aldo Moro", via E. Orabona, 4, I-70125 Bari, Italy
| | - Michele Saviano
- CNR-Institute of Crystallography, Via Amendola 122/o, 70126 Bari, Italy
| | | |
Collapse
|
13
|
Mousaei M, Kudaibergenova M, MacKerell AD, Noskov S. Assessing hERG1 Blockade from Bayesian Machine-Learning-Optimized Site Identification by Ligand Competitive Saturation Simulations. J Chem Inf Model 2020; 60:6489-6501. [PMID: 33196188 PMCID: PMC7839320 DOI: 10.1021/acs.jcim.0c01065] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Drug-induced cardiotoxicity is a potentially lethal and yet one of the most common side effects with the drugs in clinical use. Most of the drug-induced cardiotoxicity is associated with an off-target pharmacological blockade of K+ currents carried out by the cardiac Human-Ether-a-go-go-Related (hERG1) potassium channel. There is a compulsory preclinical stage safety assessment for the hERG1 blockade for all classes of drugs, which adds substantially to the cost of drug development. The availability of a high-resolution cryogenic electron microscopy (cryo-EM) structure for the channel in its open/depolarized state solved in 2017 enabled the application of molecular modeling for rapid assessment of drug blockade by molecular docking and simulation techniques. More importantly, if successful, in silico methods may allow a path to lead-compound salvaging by mapping out key block determinants. Here, we report the blind application of the site identification by the ligand competitive saturation (SILCS) protocol to map out druggable/regulatory hotspots in the hERG1 channel available for blockers and activators. The SILCS simulations use small solutes representative of common functional groups to sample the chemical space for the entire protein and its environment using all-atom simulations. The resulting chemical maps, FragMaps, explicitly account for receptor flexibility, protein-fragment interactions, and fragment desolvation penalty allowing for rapid ranking of potential ligands as blockers or nonblockers of hERG1. To illustrate the power of the approach, SILCS was applied to a test set of 55 blockers with diverse chemical scaffolds and pIC50 values measured under uniform conditions. The original SILCS model was based on the all-atom modeling of the hERG1 channel in an explicit lipid bilayer and was further augmented with a Bayesian-optimization/machine-learning (BML) stage employing an independent literature-derived training set of 163 molecules. BML approach was used to determine weighting factors for the FragMaps contributions to the scoring function. pIC50 predictions from the combined SILCS/BML approach to the 55 blockers showed a Pearson correlation (PC) coefficient of >0.535 relative to the experimental data. SILCS/BML model was shown to yield substantially improved performance as compared to commonly used rigid and flexible molecular docking methods for a well-established cohort of hERG1 blockers, where no correlation with experimental data was recorded. SILCS/BML results also suggest that a proper weighting of protonation states of common blockers present at physiological pH is essential for accurate predictions of blocker potency. The precalculated and optimized SILCS FragMaps can now be used for the rapid screening of small molecules for their cardiotoxic potential as well as for exploring alternative binding pockets in the hERG1 channel with applications to the rational design of activators.
Collapse
Affiliation(s)
- Mahdi Mousaei
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Meruyert Kudaibergenova
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Alexander D. MacKerell
- Computer-Aided Drug Design Center, Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| | - Sergei Noskov
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| |
Collapse
|
14
|
Comparing AutoDock and Vina in Ligand/Decoy Discrimination for Virtual Screening. APPLIED SCIENCES-BASEL 2019. [DOI: 10.3390/app9214538] [Citation(s) in RCA: 63] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
AutoDock and Vina are two of the most widely used protein–ligand docking programs. The fact that these programs are free and available under an open source license, also makes them a very popular first choice for many users and a common starting point for many virtual screening campaigns, particularly in academia. Here, we evaluated the performance of AutoDock and Vina against an unbiased dataset containing 102 protein targets, 22,432 active compounds and 1,380,513 decoy molecules. In general, the results showed that the overall performance of Vina and AutoDock was comparable in discriminating between actives and decoys. However, the results varied significantly with the type of target. AutoDock was better in discriminating ligands and decoys in more hydrophobic, poorly polar and poorly charged pockets, while Vina tended to give better results for polar and charged binding pockets. For the type of ligand, the tendency was the same for both Vina and AutoDock. Bigger and more flexible ligands still presented a bigger challenge for these docking programs. A set of guidelines was formulated, based on the strengths and weaknesses of both docking program and their limits of validation.
Collapse
|