Molecular docking and molecular dynamics studies on the structure-activity relationship of fluoroquinolone for the HERG channel

Utilizing Heteroatom Types and Numbers from Extensive Ligand Libraries to Develop Novel hERG Blocker QSAR Models Using Machine Learning-Based Classifiers

The human ether-à-go-go-related gene (hERG) channel plays a crucial role in membrane repolarization. Any disruptions in its function can lead to severe cardiovascular disorders such as long QT syndrome (LQTS), which increases the risk of serious cardiovascular problems such as tachyarrhythmia and sudden cardiac death. Drug-induced LQTS is a significant concern and has resulted in drug withdrawals from the market in the past. The main objective of this study is to pinpoint crucial heteroatoms present in ligands that initiate interactions leading to the effective blocking of the hERG channel. To achieve this aim, ligand-based quantitative structure-activity relationships (QSAR) models were constructed using extensive ligand libraries, considering the heteroatom types and numbers, and their associated hERG channel blockage pIC50 values. Machine learning-assisted QSAR models were developed to analyze the key structural components influencing compound activity. Among the various methods, the KPLS method proved to be the most efficient, allowing the construction of models based on eight distinct fingerprints. The study delved into investigating the influence of heteroatoms on the activity of hERG blockers, revealing their significant role. Furthermore, by quantifying the effect of heteroatom types and numbers on ligand activity at the hERG channel, six compound pairs were selected for molecular docking. Subsequent molecular dynamics simulations and per residue MM/GBSA calculations were performed to comprehensively analyze the interactions of the selected pair compounds.

Machine-learning technique, QSAR and molecular dynamics for hERG-drug interactions

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.

Repurposing drugs as COVID-19 therapies: a toxicity evaluation

Drug repurposing is an appealing method to address the Coronavirus 2019 (COVID-19) pandemic because of the low cost and efficiency. We analyzed our in-house database of approved drug screens and compared their activity profiles with results from a severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) cytopathic effect (CPE) assay. The activity profiles of the human ether-à-go-go-related gene (hERG), phospholipidosis (PLD), and many cytotoxicity screens were found significantly correlated with anti-SARS-CoV-2 activity. hERG inhibition is a nonspecific off-target effect that has contributed to promiscuous drug interactions, whereas drug-induced PLD is an undesirable effect linked to hERG blockers. Thus, this study identifies preferred drug candidates as well as chemical structures that should be avoided because of their potential to induce toxicity. Lastly, we highlight the hERG liability of anti-SARS-CoV-2 drugs currently enrolled in clinical trials.

Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study

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 (AUC_MAX = 0.86 ± 0.01) and negative predictive values (NPV_MAX = 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.

Selection of safe artemisinin derivatives using a machine learning-based cardiotoxicity platform and in vitro and in vivo validation

The majority of drug candidates fails the approval phase due to unwanted toxicities and side effects. Establishment of an effective toxicity prediction platform is of utmost importance, to increase the efficiency of the drug discovery process. For this purpose, we developed a toxicity prediction platform with machine-learning strategies. Cardiotoxicity prediction was performed by establishing a model with five parameters (arrhythmia, cardiac failure, heart block, hypertension, myocardial infarction) and additional toxicity predictions such as hepatotoxicity, reproductive toxicity, mutagenicity, and tumorigenicity are performed by using Data Warrior and Pro-Tox-II software. As a case study, we selected artemisinin derivatives to evaluate the platform and to provide a list of safe artemisinin derivatives. Artemisinin from Artemisia annua was described first as an anti-malarial compound and later its anticancer properties were discovered. Here, random forest feature selection algorithm was used for the establishment of cardiotoxicity models. High AUC scores above 0.830 were achieved for all five cardiotoxicity indications. Using a chemical library of 374 artemisinin derivatives as a case study, 7 compounds (deoxydihydro-artemisinin, 3-hydroxy-deoxy-dihydroartemisinin, 3-desoxy-dihydroartemisinin, dihydroartemisinin-furano acetate-d3, deoxyartemisinin, artemisinin G, artemisinin B) passed the toxicity filtering process for hepatotoxicity, mutagenicity, tumorigenicity, and reproductive toxicity in addition to cardiotoxicity. Experimental validation with the cardiomyocyte cell line AC16 supported the findings from the in silico cardiotoxicity model predictions. Transcriptomic profiling of AC16 cells upon artemisinin B treatment revealed a similar gene expression profile as that of the control compound, dexrazoxane. In vivo experiments with a Zebrafish model further substantiated the in silico and in vitro data, as only slight cardiotoxicity in picomolar range was observed. In conclusion, our machine-learning approach combined with in vitro and in vivo experimentation represents a suitable method to predict cardiotoxicity of drug candidates.

In vivo and in silico evaluations of survival and cardiac developmental toxicity of quinolone antibiotics in zebrafish embryos (Danio rerio)

Quinolones are ranked as the second most commonly used class of antibiotics in China, despite their adverse clinical and environmental effects. However, information on their cardiac developmental toxicity to zebrafish is limited. This study investigates the relationships between different quinolone structures and toxicity in zebrafish embryos using in vivo and in silico methods. All of the experimentally tested quinolones show cardiac developmental toxicity potential and present mortality and teratogenic effects in a dose-dependent manner. Theoretically, the acute toxicity values predicted using quantitative structure-toxicity relationship (QSTR) modeling based on previously reported LC₅₀ values are in good agreement with the in vivo results. Further investigation demonstrates that the hormetic concentration response of some quinolones may be related to methylation on the piperazine ring at the C-7 position. The amino group at the C-5 position, the methylated or ethylated piperazine group at the C-7 position, halogens at the C-8 position and a cyclopropyl ring at N1 position may be responsible for cardiac developmental toxicity. In terms of survival (key ecological endpoint), the naridine ring is more toxic than the quinoline ring. This combined approach can predict the acute and cardiac developmental toxicity of other quinolones and impurities.

In silico toxicity evaluation of Salubrinal and its analogues

This paper reports on a comprehensive in silico toxicity assessment of Salubrinal and its analogues containing a cinnamic acid residue or quinoline ring using the online servers admetSAR, ADMETlab, ProTox, ADVERPred, Pred-hERG and Vienna LiverTox. Apart from rare exceptions, in all 55 studied structures, mild or practical absence of acute toxicity was predicted for rats (III or IV toxicity class). Cardiotoxic, hepatotoxic and immunotoxic effects were predicted for Salubrinal and its analogues. We constructed models of the main predicted anti-targets hERG, BSEP, MRP3, MRP4 and AhR using the principle of homologous modeling. Molecular docking studies were carried out with the obtained models. We carried out molecular docking for all targets using AutoDock Vina, implemented in the PyRx 0.8 software package. According to the results of molecular docking, the compounds analyzed are potential moderate or weak hERG blockers. Induction of cholestasis and, as a consequence, liver damage by these drugs, directly related to inhibition of BSEP, MRP3 and MRP4, most likely will not be observed. Interaction with AhR for the studied compounds is impossible for steric reasons and, as a consequence, toxic effects on the immune and other organ systems associated with the activation of the AhR signaling pathway are excluded.

Computational investigations of hERG channel blockers: New insights and current predictive models

Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.