Experimentally Validated Pharmacoinformatics Approach to Predict hERG Inhibition Potential of New Chemical Entities

In Silico Human Cardiomyocyte Action Potential Modeling: Exploring Ion Channel Input Combinations

In silico modeling offers an opportunity to supplement and accelerate cardiac safety testing. With in silico modeling, computational simulation methods are used to predict electrophysiological interactions and pharmacological effects of novel drugs on critical physiological processes. The O'Hara-Rudy's model was developed to predict the response to different ion channel inhibition levels on cardiac action potential duration (APD) which is known to directly correlate with the QT interval. APD data at 30% 60% and 90% inhibition were derived from the model to delineate possible ventricular arrhythmia scenarios and the marginal contribution of each ion channel to the model. Action potential values were calculated for epicardial, myocardial, and endocardial cells, with action potential curve modeling. This study assessed cardiac ion channel inhibition data combinations to consider when undertaking in silico modeling of proarrhythmic effects as stipulated in the Comprehensive in Vitro Proarrhythmia Assay (CiPA). As expected, our data highlight the importance of the delayed rectifier potassium channel (IKr) as the most impactful channel for APD prolongation. The impact of the transient outward potassium channel (Ito) inhibition on APD was minimal while the inward rectifier (IK1) and slow component of the delayed rectifier potassium channel (IKs) also had limited APD effects. In contrast, the contribution of fast sodium channel (INa) and/or L-type calcium channel (ICa) inhibition resulted in substantial APD alterations supporting the pharmacological relevance of in silico modeling using input from a limited number of cardiac ion channels including IKr, INa, and ICa, at least at an early stage of drug development.

Enhancing hERG Risk Assessment with Interpretable Classificatory and Regression Models

The human Ether-à-go-go-Related Gene (hERG) is a transmembrane protein that regulates cardiac action potential, and its inhibition can induce a potentially deadly cardiac syndrome. In vitro tests help identify hERG blockers at early stages; however, the high cost motivates searching for alternative, cost-effective methods. The primary goal of this study was to enhance the Pred-hERG tool for predicting hERG blockage. To achieve this, we developed new QSAR models that incorporated additional data, updated existing classificatory and multiclassificatory models, and introduced new regression models. Notably, we integrated SHAP (SHapley Additive exPlanations) values to offer a visual interpretation of these models. Utilizing the latest data from ChEMBL v30, encompassing over 14,364 compounds with hERG data, our binary and multiclassification models outperformed both the previous iteration of Pred-hERG and all publicly available models. Notably, the new version of our tool introduces a regression model for predicting hERG activity (pIC50). The optimal model demonstrated an R2 of 0.61 and an RMSE of 0.48, surpassing the only available regression model in the literature. Pred-hERG 5.0 now offers users a swift, reliable, and user-friendly platform for the early assessment of chemically induced cardiotoxicity through hERG blockage. The tool provides versatile outcomes, including (i) classificatory predictions of hERG blockage with prediction reliability, (ii) multiclassificatory predictions of hERG blockage with reliability, (iii) regression predictions with estimated pIC50 values, and (iv) probability maps illustrating the contribution of chemical fragments for each prediction. Furthermore, we implemented explainable AI analysis (XAI) to visualize SHAP values, providing insights into the contribution of each feature to binary classification predictions. A consensus prediction calculated based on the predictions of the three developed models is also present to assist the user's decision-making process. Pred-hERG 5.0 has been designed to be user-friendly, making it accessible to users without computational or programming expertise. The tool is freely available at http://predherg.labmol.com.br.

Computational Screening Using a Combination of Ligand-Based Machine Learning and Molecular Docking Methods for the Repurposing of Antivirals Targeting the SARS-CoV-2 Main Protease

BACKGROUND

COVID-19 is an infectious disease caused by SARS-CoV-2, a close relative of SARS-CoV. Several studies have searched for COVID-19 therapies. The topics of these works ranged from vaccine discovery to natural products targeting the SARS-CoV-2 main protease (Mpro), a potential therapeutic target due to its essential role in replication and conserved sequences. However, published research on this target is limited, presenting an opportunity for drug discovery and development.

METHOD

This study aims to repurpose 10692 drugs in DrugBank by using ligand-based virtual screening (LBVS) machine learning (ML) with Konstanz Information Miner (KNIME) to seek potential therapeutics based on Mpro inhibitors. The top candidate compounds, the native ligand (GC-376) of the Mpro inhibitor, and the positive control boceprevir were then subjected to absorption, distribution, metabolism, excretion, and toxicity (ADMET) characterization, drug-likeness prediction, and molecular docking (MD). Protein-protein interaction (PPI) network analysis was added to provide accurate information about the Mpro regulatory network.

RESULTS

This study identified 3,166 compound candidates inhibiting Mpro. The random forest (RF) molecular access system ML model provided the highest confidence score of 0.95 (bromo-7-nitroindazole) and identified the top 22 candidate compounds. Subjecting the 22 candidate compounds, the native ligand GC-376, and boceprevir to further ADMET property characterization and drug-likeness predictions revealed that one compound had two violations of Lipinski's rule. Additional MD results showed that only five compounds had more negative binding energies than the native ligand (- 12.25 kcal/mol). Among these compounds, CCX-140 exhibited the lowest score of - 13.64 kcal/mol. Through literature analysis, six compound classes with potential activity for Mpro were discovered. They included benzopyrazole, azole, pyrazolopyrimidine, carboxylic acids and derivatives, benzene and substituted derivatives, and diazine. Four pathologies were also discovered on the basis of the Mpro PPI network.

CONCLUSION

Results demonstrated the efficiency of LBVS combined with MD. This combined strategy provided positive evidence showing that the top screened drugs, including CCX-140, which had the lowest MD score, can be reasonably advanced to the in vitro phase. This combined method may accelerate the discovery of therapies for novel or orphan diseases from existing drugs.

Predicting Cardiotoxicity of Molecules Using Attention-Based Graph Neural Networks

In drug discovery, the search for new and effective medications is often hindered by concerns about toxicity. Numerous promising molecules fail to pass the later phases of drug development due to strict toxicity assessments. This challenge significantly increases the cost, time, and human effort needed to discover new therapeutic molecules. Additionally, a considerable number of drugs already on the market have been withdrawn or re-evaluated because of their unwanted side effects. Among the various types of toxicity, drug-induced heart damage is a severe adverse effect commonly associated with several medications, especially those used in cancer treatments. Although a number of computational approaches have been proposed to identify the cardiotoxicity of molecules, the performance and interpretability of the existing approaches are limited. In our study, we proposed a more effective computational framework to predict the cardiotoxicity of molecules using an attention-based graph neural network. Experimental results indicated that the proposed framework outperformed the other methods. The stability of the model was also confirmed by our experiments. To assist researchers in evaluating the cardiotoxicity of molecules, we have developed an easy-to-use online web server that incorporates our model.

Insights into Drug Cardiotoxicity from Biological and Chemical Data: The First Public Classifiers for FDA Drug-Induced Cardiotoxicity Rank

Drug-induced cardiotoxicity (DICT) is a major concern in drug development, accounting for 10-14% of postmarket withdrawals. In this study, we explored the capabilities of chemical and biological data to predict cardiotoxicity, using the recently released DICTrank data set from the United States FDA. We found that such data, including protein targets, especially those related to ion channels (e.g., hERG), physicochemical properties (e.g., electrotopological state), and peak concentration in plasma offer strong predictive ability for DICT. Compounds annotated with mechanisms of action such as cyclooxygenase inhibition could distinguish between most-concern and no-concern DICT. Cell Painting features for ER stress discerned most-concern cardiotoxic from nontoxic compounds. Models based on physicochemical properties provided substantial predictive accuracy (AUCPR = 0.93). With the availability of omics data in the future, using biological data promises enhanced predictability and deeper mechanistic insights, paving the way for safer drug development. All models from this study are available at https://broad.io/DICTrank_Predictor.

Derivation of Highly Predictive 3D-QSAR Models for hERG Channel Blockers Based on the Quantum Artificial Neural Network Algorithm

The hERG potassium channel serves as an annexed target for drug discovery because the associated off-target inhibitory activity may cause serious cardiotoxicity. Quantitative structure-activity relationship (QSAR) models were developed to predict inhibitory activities against the hERG potassium channel, utilizing the three-dimensional (3D) distribution of quantum mechanical electrostatic potential (ESP) as the molecular descriptor. To prepare the optimal atomic coordinates of dataset molecules, pairwise 3D structural alignments were carried out in order for the quantum mechanical cross correlation between the template and other molecules to be maximized. This alignment method stands out from the common atom-by-atom matching technique, as it can handle structurally diverse molecules as effectively as chemical derivatives that share an identical scaffold. The alignment problem prevalent in 3D-QSAR methods was ameliorated substantially by dividing the dataset molecules into seven subsets, each of which contained molecules with similar molecular weights. Using an artificial neural network algorithm to find the functional relationship between the quantum mechanical ESP descriptors and the experimental hERG inhibitory activities, highly predictive 3D-QSAR models were derived for all seven molecular subsets to the extent that the squared correlation coefficients exceeded 0.79. Given their simplicity in model development and strong predictability, the 3D-QSAR models developed in this study are expected to function as an effective virtual screening tool for assessing the potential cardiotoxicity of drug candidate molecules.

Insights into Drug Cardiotoxicity from Biological and Chemical Data: The First Public Classifiers for FDA DICTrank

Drug-induced cardiotoxicity (DICT) is a major concern in drug development, accounting for 10-14% of postmarket withdrawals. In this study, we explored the capabilities of various chemical and biological data to predict cardiotoxicity, using the recently released Drug-Induced Cardiotoxicity Rank (DICTrank) dataset from the United States FDA. We analyzed a diverse set of data sources, including physicochemical properties, annotated mechanisms of action (MOA), Cell Painting, Gene Expression, and more, to identify indications of cardiotoxicity. We found that such data, including protein targets, especially those related to ion channels (such as hERG), physicochemical properties (such as electrotopological state) as well as peak concentration in plasma offer strong predictive ability as well as valuable insights into DICT. We also found compounds annotated with particular mechanisms of action, such as cyclooxygenase inhibition, could distinguish between most-concern and no-concern DICT compounds. Cell Painting features related to ER stress discern the most-concern cardiotoxic compounds from non-toxic compounds. While models based on physicochemical properties currently provide substantial predictive accuracy (AUCPR = 0.93), this study also underscores the potential benefits of incorporating more comprehensive biological data in future DICT predictive models. With the availability of - omics data in the future, using biological data promises enhanced predictability and delivers deeper mechanistic insights, paving the way for safer therapeutic drug development. All models and data used in this study are publicly released at https://broad.io/DICTrank_Predictor.

Mitigating hERG Liability of Toll-Like Receptor 9 and 7 Antagonists through Structure-Based Design

hERG is considered to be a primary anti-target in the drug development process, as the K+ channel encoded by hERG plays an important role in cardiac re-polarization. It is desirable to address the hERG safety liability during early-stage development to avoid the expenses of validating leads that will eventually fail at a later stage. We have previously reported the development of highly potent quinazoline-based TLR7 and TLR9 antagonists for possible application against autoimmune disease. Initial experimental hERG assessment showed that most of the lead TLR7 and TLR9 antagonists suffer from hERG liability rendering them ineffective for further development. The present study herein describes a coordinated strategy to integrate the understanding from structure-based protein-ligand interaction to develop non- hERG binders with IC50 >30 μM with retention of TLR7/9 antagonism through a single point change in the scaffold. This structure-guided strategy can serve as a prototype for abolishing hERG liability during lead optimization.

Structural analysis of hERG channel blockers and the implications for drug design

The repolarizing current (I_kr) produced by the hERG potassium channel forms a major component of the cardiac action potential and blocking this current by small molecule drugs can lead to life-threatening cardiotoxicity. Understanding the mechanisms of drug-mediated hERG inhibition is essential to develop a second generation of safe drugs, with minimal cardiotoxic effects. Although various computational tools and drug design guidelines have been developed to avoid binding of drugs to the hERG pore domain, there are many other aspects that are still open for investigation. This includes the use computational modelling to study the implications of hERG mutations on hERG structure and trafficking, the interactions of hERG with hERG chaperone proteins and with membrane-soluble molecules, the mechanisms of drugs that inhibit hERG trafficking and drugs that rescue hERG mutations. The plethora of available experimental data regarding all these aspects can guide the construction of much needed robust computational structural models to study these mechanisms for the rational design of safe drugs.

Identification and Empiric Evaluation of New Inhibitors of the Multidrug Transporter P-Glycoprotein (ABCB1)

The expression of the drug efflux pump ABCB1 correlates negatively with cancer survival, making the transporter an attractive target for therapeutic inhibition. In order to identify new inhibitors of ABCB1, we have exploited the cryo-EM structure of the protein to develop a pharmacophore model derived from the best docked conformations of a structurally diverse range of known inhibitors. The pharmacophore model was used to screen the Chembridge compound library. We identified six new potential inhibitors with distinct chemistry compared to the third-generation inhibitor tariquidar and with favourable lipophilic efficiency (LipE) and lipophilicity (CLogP) characteristics, suggesting oral bioavailability. These were evaluated experimentally for efficacy and potency using a fluorescent drug transport assay in live cells. The half-maximal inhibitory concentrations (IC50) of four of the compounds were in the low nanomolar range (1.35 to 26.4 nM). The two most promising compounds were also able to resensitise ABCB1-expressing cells to taxol. This study demonstrates the utility of cryo-electron microscopy structure determination for drug identification and design.

Triad pyrazole-thiazole-coumarin heterocyclic core effectively inhibit HSP and drive cancer cells to apoptosis

Intensive studies on hepatocellular carcinoma (HCC), which is spreading rapidly around the world and has a high mortality rate, is due to the lack of adequate preventive or curative treatment methods. Treating patients with HCC has become very challenging because of the heterogeneity in the patient population lead activation of different signaling pathways, and pathway crosstalk for patients. Therefore, understanding these molecular mechanisms and combining drugs with molecular therapies to overcome these drawbacks has become an area of utmost importance. In this study, the biological activities of the designed and characterized triad Pyrazole-Thiazol-Coumarin (PTC) compounds were determined by performing cell viability, qPCR array, apoptosis and cell cycle assays. One of the compounds (PTC10) implicitly suppresses multiple pathways (RAS/MAP kinase and PI3K-AKT) simultaneously. This action is provided by (i) arresting cancer cells at G2 phase, (ii) driving cancer cells to apoptosis and (iii) inhibiting HSP network. Remarkably, HSP is an apoptotic factor and help cancer cell to survive. HSP90 also coordinates with Cdk4/Cdc37, therefore inhibiting HSP both drives cells to arrest and apoptosis. ATP hydrolysis and aggregation assay further displayed specific HSP inhibition. Therefore, PTC provides a unique drug template for HCC treatment.

cardioToxCSM: A Web Server for Predicting Cardiotoxicity of Small Molecules

The design of novel, safe, and effective drugs to treat human diseases is a challenging venture, with toxicity being one of the main sources of attrition at later stages of development. Failure due to toxicity incurs a significant increase in costs and time to market, with multiple drugs being withdrawn from the market due to their adverse effects. Cardiotoxicity, for instance, was responsible for the failure of drugs such as fenspiride, propoxyphene, and valdecoxib. While significant effort has been dedicated to mitigate this issue by developing computational approaches that aim to identify molecules likely to be toxic, including quantitative structure-activity relationship models and machine learning methods, current approaches present limited performance and interpretability. To overcome these, we propose a new web-based computational method, cardioToxCSM, which can predict six types of cardiac toxicity outcomes, including arrhythmia, cardiac failure, heart block, hERG toxicity, hypertension, and myocardial infarction, efficiently and accurately. cardioToxCSM was developed using the concept of graph-based signatures, molecular descriptors, toxicophore matchings, and molecular fingerprints, leveraging explainable machine learning, and was validated internally via different cross validation schemes and externally via low-redundancy blind sets. The models presented robust performances with areas under ROC curves of up to 0.898 on 5-fold cross-validation, consistent with metrics on blind tests. Additionally, our models provide interpretation of the predictions by identifying whether substructures that are commonly enriched in toxic compounds were present. We believe cardioToxCSM will provide valuable insight into the potential cardiotoxicity of small molecules early on drug screening efforts. The method is made freely available as a web server at https://biosig.lab.uq.edu.au/cardiotoxcsm.

Ligand-based prediction of hERG-mediated cardiotoxicity based on the integration of different machine learning techniques

Drug-induced cardiotoxicity is a common side effect of drugs in clinical use or under postmarket surveillance and is commonly due to off-target interactions with the cardiac human-ether-a-go-go-related (hERG) potassium channel. Therefore, prioritizing drug candidates based on their hERG blocking potential is a mandatory step in the early preclinical stage of a drug discovery program. Herein, we trained and properly validated 30 ligand-based classifiers of hERG-related cardiotoxicity based on 7,963 curated compounds extracted by the freely accessible repository ChEMBL (version 25). Different machine learning algorithms were tested, namely, random forest, K-nearest neighbors, gradient boosting, extreme gradient boosting, multilayer perceptron, and support vector machine. The application of 1) the best practices for data curation, 2) the feature selection method VSURF, and 3) the synthetic minority oversampling technique (SMOTE) to properly handle the unbalanced data, allowed for the development of highly predictive models (BAMAX = 0.91, AUCMAX = 0.95). Remarkably, the undertaken temporal validation approach not only supported the predictivity of the herein presented classifiers but also suggested their ability to outperform those models commonly used in the literature. From a more methodological point of view, the study put forward a new computational workflow, freely available in the GitHub repository (https://github.com/PDelre93/hERG-QSAR), as valuable for building highly predictive models of hERG-mediated cardiotoxicity.

hERG Blockade Prediction by Combining Site Identification by Ligand Competitive Saturation and Physicochemical Properties

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.

The utility of hERG channel inhibition data in the derivation of occupational exposure limits

Inhibition of the human ether-à-go-go (hERG) channel may lead to QT prolongation and fatal arrhythmia. While pharmaceutical drug candidates that exhibit potent hERG channel inhibition often fail early in development, many drugs with both cardiac and non-cardiac indications proceed to market. In this study, the relationship between in vitro hERG channel inhibition and published occupational exposure limit (OEL) was evaluated. A total of 23 cardiac drugs and 44 drugs with non-cardiac indications with published hERG channel IC₅₀ and published OELs were identified. There was an apparent relationship between hERG IC₅₀ potency and the OEL for cardiac and non-cardiac drugs. Twenty cardiac and non-cardiac drugs were identified that had a potent hERG IC₅₀ (≤25 μM) and a contrastingly large OEL value (≥100 μg/m³). OELs or hazard banding corresponding to ≤100 μg/m³ should be sufficiently protective of effects following occupational exposure to the majority of APIs with hERG IC₅₀ values ≤ 100 μM. It is important to consider hERG IC₅₀ values and possible cardiac effects when deriving OEL values for drugs, regardless of indication. These considerations may be particularly important early in the drug development process for establishing exposure control bands for drugs that do not yet have full clinical safety data.

Molecular Structure-Based Screening of the Constituents of Calotropis procera Identifies Potential Inhibitors of Diabetes Mellitus Target Alpha Glucosidase

Diabetes mellitus is a disorder characterized by higher levels of blood glucose due to impaired insulin mechanisms. Alpha glucosidase is a critical drug target implicated in the mechanisms of diabetes mellitus and its inhibition controls hyperglycemia. Since the existing standard synthetic drugs have therapeutic limitations, it is imperative to identify new potent inhibitors of natural product origin which may slow carbohydrate digestion and absorption via alpha glucosidase. Since plant extracts from Calotropis procera have been extensively used in the treatment of diabetes mellitus, the present study used molecular docking and dynamics simulation techniques to screen its constituents against the receptor alpha glucosidase. Taraxasterol, syriogenin, isorhamnetin-3-O-robinobioside and calotoxin were identified as potential novel lead compounds with plausible binding energies of −40.2, −35.1, −34.3 and −34.3 kJ/mol against alpha glucosidase, respectively. The residues Trp⁴⁸¹, Asp⁵¹⁸, Leu⁶⁷⁷, Leu⁶⁷⁸ and Leu⁶⁸⁰ were identified as critical for binding and the compounds were predicted as alpha glucosidase inhibitors. Structurally similar compounds with Tanimoto coefficients greater than 0.7 were reported experimentally to be inhibitors of alpha glucosidase or antidiabetic. The structures of the molecules may serve as templates for the design of novel inhibitors and warrant in vitro assaying to corroborate their antidiabetic potential.

Combined Pharmacophore and Grid-Independent Molecular Descriptors (GRIND) Analysis to Probe 3D Features of Inositol 1,4,5-Trisphosphate Receptor (IP3R) Inhibitors in Cancer

Inositol 1, 4, 5-trisphosphate receptor (IP₃R)-mediated Ca²⁺ signaling plays a pivotal role in different cellular processes, including cell proliferation and cell death. Remodeling Ca²⁺ 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 IP₃-binding core (IBC) has been proposed yet. Therefore, to elucidate the 3D structural features of IP₃R 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 (IC₅₀) of a compound for IP₃R inhibition. Moreover, our GRIND model (PLS: Q² = 0.70 and R² = 0.72) further strengthens the identified pharmacophore features of IP₃R 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 IP₃R 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 IP₃R. The identified hits could further assist in the design and optimization of lead structures for the targeting and remodeling of Ca²⁺ signals in cancer.

In silico approaches in organ toxicity hazard assessment: Current status and future needs for predicting heart, kidney and lung toxicities

The kidneys, heart and lungs are vital organ systems evaluated as part of acute or chronic toxicity assessments. New methodologies are being developed to predict these adverse effects based on in vitro and in silico approaches. This paper reviews the current state of the art in predicting these organ toxicities. It outlines the biological basis, processes and endpoints for kidney toxicity, pulmonary toxicity, respiratory irritation and sensitization as well as functional and structural cardiac toxicities. The review also covers current experimental approaches, including off-target panels from secondary pharmacology batteries. Current in silico approaches for prediction of these effects and mechanisms are described as well as obstacles to the use of in silico methods. Ultimately, a commonly accepted protocol for performing such assessment would be a valuable resource to expand the use of such approaches across different regulatory and industrial applications. However, a number of factors impede their widespread deployment including a lack of a comprehensive mechanistic understanding, limited in vitro testing approaches and limited in vivo databases suitable for modeling, a limited understanding of how to incorporate absorption, distribution, metabolism, and excretion (ADME) considerations into the overall process, a lack of in silico models designed to predict a safe dose and an accepted framework for organizing the key characteristics of these organ toxicants.

Prediction of hERG inhibition of drug discovery compounds using biomimetic HPLC measurements

The major causes of failure of drug discovery compounds in clinics are the lack of efficacy and toxicity. To reduce late-stage failures in the drug discovery process, it is essential to estimate early the probability of adverse effects and potential toxicity. Cardiotoxicity is one of the most often observed problems related to a compound's inhibition of the hERG channel responsible for the potassium cation flux. Biomimetic HPLC methods can be used for the early screening of a compound's lipophilicity, protein binding and phospholipid partition. Based on the published hERG pIC50 data of 90 marketed drugs and their measured biomimetic properties, a model has been developed to predict the hERG inhibition using the measured binding of compounds to alpha-1-acid-glycoprotein (AGP) and immobilised artificial membrane (IAM). A representative test set of 16 compounds was carefully selected. The training set, involving the remaining compounds, served to establish the linear model. The mechanistic model supports the hypothesis that compounds have to traverse the cell membrane and bind to the hERG ion channel to cause the inhibition. The AGP and the hERG ion channel show structural similarity, as both bind positively charged compounds with strong shape selectivity. In contrast, a good IAM partition is a prerequisite for cell membrane traversal. For reasons of comparison, a corresponding model was derived by replacing the measured biomimetic properties with calculated physicochemical properties. The model established with the measured biomimetic binding properties proved to be superior and can explain over 70% of the variance of the hERG pIC50 values.

Molecular docking and pharmacophore models to probe binding hypothesis of inhibitors of hypoxia inducible factor-1

Hypoxia inducible factor-1 is a heterodimeric transcription factor that regulates cellular responses to hypoxia and is involved in tumor progression and resistance to chemotherapy. Dimerization between HIF-1α and β subunits has been recognized crucial for DNA binding and transcriptional activity of HIF-1. Therefore, inhibitors of α and β dimerization subunits of HIF-1 may potentially evade HIF-1-mediated chemotherapy resistance. In the current study, ligand-based pharmacophore model was developed to determine 3 D binding features of HIF-1 inhibitors. The selected pharmacophore model comprises of one hydrogen bond donor, one hydrogen bond acceptor and one hydrophobic feature. The selected model was used for virtual screening of publically available data base by ChemBridge Corporation. Overall, six potential hits against HIF-1α and β dimerization have been identified. These include, Hit 1 (4-(4-chlorophenyl)-2,6-dimethyl-3,5-pyridinedicarboxylic acid), 3 (2-[2-(2-hydroxybenzoyl)carbonohydrazonoyl]benzoic acid) and 5 (3-(4-methoxyphenyl)-2,4-quinolinedicarboxylic acid) nicotonic acid derivatives, Hit 2 (3-[(1-adamantylamino)sulfonyl]benzoic acid), 4 (5-{[(2-fluorophenyl)amino]sulfonyl}-2-methylbenzoic acid), and 6 (4-({[2-(trifluoromethyl)phenyl]sulfonyl}amino)benzoic acid) sulfonamide derivatives. Additionally, adamantyl moiety of compound 2 shows interactions with the experimentally known hydrophobic amino acid residues (V336, C334, E245) of HIF-1α and β dimerization site. The identified hits showed lower to higher µM biological activity (IC₅₀) values and thus, after further structure optimization may serve as potential inhibitor of HIF-1 dimerization in cancer chemotherapy.Communicated by Ramaswamy H. Sarma.

Effects of Cryopreservation on Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes for Assessing Drug Safety Response Profiles

Burgeoning applications of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in disease modeling, regenerative medicine, and drug screening have broadened the usage of hiPSC-CMs and entailed their long-term storage. Cryopreservation is the most common approach to store hiPSC-CMs. However, the effects of cryopreservation and recovery on hiPSC-CMs remain poorly understood. Here, we characterized the transcriptome, electro-mechanical function, and drug response of fresh hiPSC-CMs without cryopreservation and recovered hiPSC-CMs from cryopreservation. We found that recovered hiPSC-CMs showed upregulation of cell cycle genes, similar or reduced contractility, Ca2+ transients, and field potential duration. When subjected to treatment of drugs that affect electrophysiological properties, recovered hiPSC-CMs showed an altered drug response and enhanced propensity for drug-induced cardiac arrhythmic events. In conclusion, fresh and recovered hiPSC-CMs do not always show comparable molecular and physiological properties. When cryopreserved hiPSC-CMs are used for assessing drug-induced cardiac liabilities, the altered drug sensitivity needs to be considered.

Artificial Intelligence in Drug Discovery: A Comprehensive Review of Data-driven and Machine Learning Approaches

As expenditure on drug development increases exponentially, the overall drug discovery process requires a sustainable revolution. Since artificial intelligence (AI) is leading the fourth industrial revolution, AI can be considered as a viable solution for unstable drug research and development. Generally, AI is applied to fields with sufficient data such as computer vision and natural language processing, but there are many efforts to revolutionize the existing drug discovery process by applying AI. This review provides a comprehensive, organized summary of the recent research trends in AI-guided drug discovery process including target identification, hit identification, ADMET prediction, lead optimization, and drug repositioning. The main data sources in each field are also summarized in this review. In addition, an in-depth analysis of the remaining challenges and limitations will be provided, and proposals for promising future directions in each of the aforementioned areas.

A current understanding of drug-induced QT prolongation and its implications for anticancer therapy

The QT interval, a global index of ventricular repolarization, varies among individuals and is influenced by diverse physiologic and pathophysiologic stimuli such as gender, age, heart rate, electrolyte concentrations, concomitant cardiac disease, and other diseases such as diabetes. Many drugs produce a small but reproducible effect on QT interval but in rare instances this is exaggerated and marked QT prolongation can provoke the polymorphic ventricular tachycardia 'torsades de pointes', which can cause syncope or sudden cardiac death. The generally accepted common mechanism whereby drugs prolong QT is block of a key repolarizing potassium current in heart, IKr, generated by expression of KCNH2, also known as HERG. Thus, evaluation of the potential that a new drug entity may cause torsades de pointes has relied on exposure of normal volunteers or patients to drug at usual and high concentrations, and on assessment of IKr block in vitro. More recent work, focusing on anticancer drugs with QT prolonging liability, is defining new pathways whereby drugs can prolong QT. Notably, the in vitro effects of some tyrosine kinase inhibitors to prolong cardiac action potentials (the cellular correlate of QT) can be rescued by intracellular phosphatidylinositol 3,4,5-trisphosphate, the downstream effector of phosphoinositide 3-kinase. This finding supports a role for inhibition of this enzyme, either directly or by inhibition of upstream kinases, to prolong QT through mechanisms that are being worked out, but include enhanced inward 'late' sodium current during the plateau of the action potential. The definition of non-IKr-dependent pathways to QT prolongation will be important for assessing risk, not only with anticancer therapies but also with other QT prolonging drugs and for generating a refined understanding how variable activity of intracellular signalling systems can modulate QT and associated arrhythmia risk.

DeepHIT: a deep learning framework for prediction of hERG-induced cardiotoxicity

Motivation

Blockade of the human ether-à-go-go-related gene (hERG) channel by small compounds causes a prolonged QT interval that can lead to severe cardiotoxicity and is a major cause of the many failures in drug development. Thus, evaluating the hERG-blocking activity of small compounds is important for successful drug development. To this end, various computational prediction tools have been developed, but their prediction performances in terms of sensitivity and negative predictive value (NPV) need to be improved to reduce false negative predictions.

Results

We propose a computational framework, DeepHIT, which predicts hERG blockers and non-blockers for input compounds. For the development of DeepHIT, we generated a large-scale gold-standard dataset, which includes 6632 hERG blockers and 7808 hERG non-blockers. DeepHIT is designed to contain three deep learning models to improve sensitivity and NPV, which, in turn, produce fewer false negative predictions. DeepHIT outperforms currently available tools in terms of accuracy (0.773), MCC (0.476), sensitivity (0.833) and NPV (0.643) on an external test dataset. We also developed an in silico chemical transformation module that generates virtual compounds from a seed compound, based on the known chemical transformation patterns. As a proof-of-concept study, we identified novel urotensin II receptor (UT) antagonists without hERG-blocking activity derived from a seed compound of a previously reported UT antagonist (KR-36676) with a strong hERG-blocking activity. In summary, DeepHIT will serve as a useful tool to predict hERG-induced cardiotoxicity of small compounds in the early stages of drug discovery and development.

Availability and implementation

https://bitbucket.org/krictai/deephit and https://bitbucket.org/krictai/chemtrans

Supplementary information

Supplementary data are available at Bioinformatics online.

Revealing Molecular Determinants of hERG Blocker and Activator Binding

The Kv11.1 potassium channel, encoded by the human ether-a-go-go-related gene (hERG), plays an essential role in the cardiac action potential. hERG blockade by small molecules can induce "torsade de pointes" arrhythmias and sudden death; as such, it is an important off-target to avoid during drug discovery. Recently, a cryo-EM structure of the open channel state of hERG was reported, opening the door to in silico docking analyses and interpretation of hERG structure-activity relationships, with a view to avoiding blocking activity. Despite this, docking directly to this cryo-EM structure has been reported to yield binding modes that are unable to explain known mutagenesis data. In this work, we use molecular dynamics simulations to sample a range of channel conformations and run ensemble docking campaigns at the known hERG binding site below the selectivity filter, composed of the central cavity and the four deep hydrophobic pockets. We identify a hERG conformational state allowing discrimination of blockers vs nonblockers from docking; furthermore, the binding pocket agrees with mutagenesis data, and blocker binding modes fit the hERG blocker pharmacophore. We then use the same protocol to identify a binding pocket in the hERG channel pore for hERG activators, again agreeing with the reported mutagenesis. Our approach may be useful in drug discovery campaigns to prioritize candidate compounds based on hERG liability via virtual docking screens.

Molecular Docking Guided Grid-Independent Descriptor Analysis to Probe the Impact of Water Molecules on Conformational Changes of hERG Inhibitors in Drug Trapping Phenomenon

Human ether a-go-go related gene (hERG) or KV11.1 potassium channels mediate the rapid delayed rectifier current (I_Kr) in cardiac myocytes. Drug-induced inhibition of hERG channels has been implicated in the development of acquired long QT syndrome type (aLQTS) and fatal arrhythmias. Several marketed drugs have been withdrawn for this reason. Therefore, there is considerable interest in developing better tests for predicting drugs which can block the hERG channel. The drug-binding pocket in hERG channels, which lies below the selectivity filter, normally contains K⁺ ions and water molecules. In this study, we test the hypothesis that these water molecules impact drug binding to hERG. We developed 3D QSAR models based on alignment independent descriptors (GRIND) using docked ligands in open and closed conformations of hERG in the presence (solvated) and absence (non-solvated) of water molecules. The ligand–protein interaction fingerprints (PLIF) scheme was used to summarize and compare the interactions. All models delineated similar 3D hERG binding features, however, small deviations of about ~0.4 Å were observed between important hotspots of molecular interaction fields (MIFs) between solvated and non-solvated hERG models. These small changes in conformations do not affect the performance and predictive power of the model to any significant extent. The model that exhibits the best statistical values was attained with a cryo_EM structure of the hERG channel in open state without water. This model also showed the best R² of 0.58 and 0.51 for the internal and external validation test sets respectively. Our results suggest that the inclusion of water molecules during the docking process has little effect on conformations and this conformational change does not impact the predictive ability of the 3D QSAR models.

Computational determination of hERG-related cardiotoxicity of drug candidates

Background

Drug candidates often cause an unwanted blockage of the potassium ion channel of the human ether-a-go-go-related gene (hERG). The blockage leads to long QT syndrome (LQTS), which is a severe life-threatening cardiac side effect. Therefore, a virtual screening method to predict drug-induced hERG-related cardiotoxicity could facilitate drug discovery by filtering out toxic drug candidates.

Result

In this study, we generated a reliable hERG-related cardiotoxicity dataset composed of 2130 compounds, which were carried out under constant conditions. Based on our dataset, we developed a computational hERG-related cardiotoxicity prediction model. The neural network model achieved an area under the receiver operating characteristic curve (AUC) of 0.764, with an accuracy of 90.1%, a Matthews correlation coefficient (MCC) of 0.368, a sensitivity of 0.321, and a specificity of 0.967, when ten-fold cross-validation was performed. The model was further evaluated using ten drug compounds tested on guinea pigs and showed an accuracy of 80.0%, an MCC of 0.655, a sensitivity of 0.600, and a specificity of 1.000, which were better than the performances of existing hERG-toxicity prediction models.

Conclusion

The neural network model can predict hERG-related cardiotoxicity of chemical compounds with a high accuracy. Therefore, the model can be applied to virtual high-throughput screening for drug candidates that do not cause cardiotoxicity. The prediction tool is available as a web-tool at http://ssbio.cau.ac.kr/CardPred.

Avoiding hERG-liability in drug design via synergetic combinations of different (Q)SAR methodologies and data sources: a case study in an industrial setting

In this paper, we explore the impact of combining different in silico prediction approaches and data sources on the predictive performance of the resulting system. We use inhibition of the hERG ion channel target as the endpoint for this study as it constitutes a key safety concern in drug development and a potential cause of attrition. We will show that combining data sources can improve the relevance of the training set in regard of the target chemical space, leading to improved performance. Similarly we will demonstrate that combining multiple statistical models together, and with expert systems, can lead to positive synergistic effects when taking into account the confidence in the predictions of the merged systems. The best combinations analyzed display a good hERG predictivity. Finally, this work demonstrates the suitability of the SOHN methodology for building models in the context of receptor based endpoints like hERG inhibition when using the appropriate pharmacophoric descriptors.

Open-Access Activity Prediction Tools for Natural Products. Case Study: hERG Blockers

Interference with the hERG potassium ion channel may cause cardiac arrhythmia and can even lead to death. Over the last few decades, several drugs, already on the market, and many more investigational drugs in various development stages, have had to be discontinued because of their hERG-associated toxicity. To recognize potential hERG activity in the early stages of drug development, a wide array of computational tools, based on different principles, such as 3D QSAR, 2D and 3D similarity, and machine learning, have been developed and are reviewed in this chapter. The various available prediction tools Similarity Ensemble Approach, SuperPred, SwissTargetPrediction, HitPick, admetSAR, PASSonline, Pred-hERG, and VirtualToxLab™ were used to screen a dataset of known hERG synthetic and natural product actives and inactives to quantify and compare their predictive power. This contribution will allow the reader to evaluate the suitability of these computational methods for their own related projects. There is an unmet need for natural product-specific prediction tools in this field.