A structure-based computational workflow to predict liability and binding modes of small molecules to hERG

CardioGenAI: a machine learning-based framework for re-engineering drugs for reduced hERG liability

The link between in vitro hERG ion channel inhibition and subsequent in vivo QT interval prolongation, a critical risk factor for the development of arrythmias such as Torsade de Pointes, is so well established that in vitro hERG activity alone is often sufficient to end the development of an otherwise promising drug candidate. It is therefore of tremendous interest to develop advanced methods for identifying hERG-active compounds in the early stages of drug development, as well as for proposing redesigned compounds with reduced hERG liability and preserved primary pharmacology. In this work, we present CardioGenAI, a machine learning-based framework for re-engineering both developmental and commercially available drugs for reduced hERG activity while preserving their pharmacological activity. The framework incorporates novel state-of-the-art discriminative models for predicting hERG channel activity, as well as activity against the voltage-gated NaV1.5 and CaV1.2 channels due to their potential implications in modulating the arrhythmogenic potential induced by hERG channel blockade. We applied the complete framework to pimozide, an FDA-approved antipsychotic agent that demonstrates high affinity to the hERG channel, and generated 100 refined candidates. Remarkably, among the candidates is fluspirilene, a compound which is of the same class of drugs as pimozide (diphenylmethanes) and therefore has similar pharmacological activity, yet exhibits over 700-fold weaker binding to hERG. Furthermore, we demonstrated the framework's ability to optimize hERG, NaV1.5 and CaV1.2 profiles of multiple FDA-approved compounds while maintaining the physicochemical nature of the original drugs. We envision that this method can effectively be applied to developmental compounds exhibiting hERG liabilities to provide a means of rescuing drug development programs that have stalled due to hERG-related safety concerns. Additionally, the discriminative models can also serve independently as effective components of virtual screening pipelines. We have made all of our software open-source at https://github.com/gregory-kyro/CardioGenAI to facilitate integration of the CardioGenAI framework for molecular hypothesis generation into drug discovery workflows.Scientific contributionThis work introduces CardioGenAI, an open-source machine learning-based framework designed to re-engineer drugs for reduced hERG liability while preserving their pharmacological activity. The complete CardioGenAI framework can be applied to developmental compounds exhibiting hERG liabilities to provide a means of rescuing drug discovery programs facing hERG-related challenges. In addition, the framework incorporates novel state-of-the-art discriminative models for predicting hERG, NaV1.5 and CaV1.2 channel activity, which can function independently as effective components of virtual screening pipelines.

Role of CSF1R 550th-tryptophan in kusunokinin and CSF1R inhibitor binding and ligand-induced structural effect

Binding affinity is an important factor in drug design to improve drug-target selectivity and specificity. In this study, in silico techniques based on molecular docking followed by molecular dynamics (MD) simulations were utilized to identify the key residue(s) for CSF1R binding affinity among 14 pan-tyrosine kinase inhibitors and 15 CSF1R-specific inhibitors. We found tryptophan at position 550 (W550) on the CSF1R binding site interacted with the inhibitors' aromatic ring in a π-π way that made the ligands better at binding. Upon W550-Alanine substitution (W550A), the binding affinity of trans-(-)-kusunokinin and imatinib to CSF1R was significantly decreased. However, in terms of structural features, W550 did not significantly affect overall CSF1R structure, but provided destabilizing effect upon mutation. The W550A also did not either cause ligand to change its binding site or conformational changes due to ligand binding. As a result of our findings, the π-π interaction with W550's aromatic ring could be still the choice for increasing binding affinity to CSF1R. Nevertheless, our study showed that the increasing binding to W550 of the design ligand may not ensure CSF1R specificity and inhibition since W550-ligand bound state did not induce significantly conformational change into inactive state.

Applications of Molecular Dynamics Simulations in Drug Discovery

In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.

Synthesis of 1,3,5-triphenyl-1,2,4-triazole derivatives and their neuroprotection by anti-oxidative stress and anti-inflammation and protecting BBB

Acute ischemic stroke (AIS) is a serious cardiovascular and cerebrovascular disease; Oxidative stress and neuroinflammation are important factors which destroy blood-brain barrier (BBB) in AIS. In the study, a series of 1,3,5-triphenyl-1,2,4-triazole derivatives were designed and synthesized; the optimal compound 9 was obtained by screening their anti-oxidant and anti-inflammatory effects; the neuroprotection effect of compound 9 was evaluated with a rat middle cerebral artery occlusion (MCAO) model. Subsequently, the mechanism of neuroprotection were explored via Western blot. The results prompt compound 9 maybe exert anti-AIS neuroprotection by inhibiting oxidative stress and neuroinflammation inhibition by inhibiting Keap1, COX-2 and iNOS. At the same time, it can protect BBB by reducing glycocalyx degradation and matrix metallopeptidase-9 levels. Its LD50 > 1000 mg/kg on mice and hERG channel inhibition IC50 > 30 μM, which lower acute toxicity and hERG channel inhibition would make compound 9 a promising stroke treatment candidate.

Structural modeling of hERG channel-drug interactions using Rosetta

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.

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.

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.

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.

Use of Solvent Mapping for Characterizing the Binding Site and for Predicting the Inhibition of the Human Ether-á-Go-Go-Related K+ Channel

Molecular dynamics was used to optimize the droperidol-hERG complex obtained from docking. To accommodate the inhibitor, residues T623, S624, V625, G648, Y652, and F656 did not move significantly during the simulation, while F627 moved significantly. Binding sites in cryo-EM structures and in structures obtained from molecular dynamics simulations were characterized using solvent mapping and Atlas ligands, which were negative images of the binding site, were generated. Atlas ligands were found to be useful for identifying human ether-á-go-go-related potassium channel (hERG) inhibitors by aligning compounds to them or by guiding the docking of compounds in the binding site. A molecular dynamics optimized structure of hERG led to improved predictions using either compound alignment to the Atlas ligand or docking. The structure was also found to be suitable to define a strategy for lowering inhibition based on the proposed binding mode of compounds in the channel.

TSSF-hERG: A machine-learning-based hERG potassium channel-specific scoring function for chemical cardiotoxicity prediction

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 (I_Kr) 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.

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.

Molecular Dynamics-Derived Pharmacophore Model Explaining the Nonselective Aspect of KV10.1 Pore Blockers

The K_V10.1 voltage-gated potassium channel is highly expressed in 70% of tumors, and thus represents a promising target for anticancer drug discovery. However, only a few ligands are known to inhibit K_V10.1, and almost all also inhibit the very similar cardiac hERG channel, which can lead to undesirable side-effects. In the absence of the structure of the K_V10.1–inhibitor complex, there remains the need for new strategies to identify selective K_V10.1 inhibitors and to understand the binding modes of the known K_V10.1 inhibitors. To investigate these binding modes in the central cavity of K_V10.1, a unique approach was used that allows derivation and analysis of ligand–protein interactions from molecular dynamics trajectories through pharmacophore modeling. The final molecular dynamics-derived structure-based pharmacophore model for the simulated K_V10.1–ligand complexes describes the necessary pharmacophore features for K_V10.1 inhibition and is highly similar to the previously reported ligand-based hERG pharmacophore model used to explain the nonselectivity of K_V10.1 pore blockers. Moreover, analysis of the molecular dynamics trajectories revealed disruption of the π–π network of aromatic residues F359, Y464, and F468 of K_V10.1, which has been reported to be important for binding of various ligands for both K_V10.1 and hERG channels. These data indicate that targeting the K_V10.1 channel pore is also likely to result in undesired hERG inhibition, and other potential binding sites should be explored to develop true K_V10.1-selective inhibitors as new anticancer agents.

Assessing hERG1 Blockade from Bayesian Machine-Learning-Optimized Site Identification by Ligand Competitive Saturation Simulations

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.