Quantitative Structure−Activity Relationship Studies on Inhibition of HERG Potassium Channels

hERG agonists pose challenges to web-based machine learning methods for prediction of drug-hERG channel interaction

Pharmacological blockade of the IKr channel (hERG) by diverse drugs in clinical use is associated with the Long QT Syndrome that can lead to life threatening arrhythmia. Various computational tools including machine learning models (MLM) for the prediction of hERG inhibition have been developed to facilitate the throughput screening of drugs in development and optimise thus the prediction of hERG liabilities. The use of MLM relies on large libraries of training compounds for the quantitative structure-activity relationship (QSAR) modelling of hERG inhibition. The focus on inhibition omits potential effects of hERG channel agonist molecules and their associated QT shortening risk. It is instructive, therefore, to consider how known hERG agonists are handled by MLM. Here, two highly developed online computational tools for the prediction of hERG liability, Pred-hERG and HergSPred were probed for their ability to detect hERG activator drug molecules as hERG interactors. In total, 73 hERG blockers were tested with both computational tools giving overall good predictions for hERG blockers with reported IC50s below Pred-hERG and HergSPred cut-off threshold for hERG inhibition. However, for compounds with reported IC50s above this threshold such as disopyramide or sotalol discrepancies were observed. HergSPred identified all 20 hERG agonists selected as interacting with the hERG channel. Further studies are warranted to improve online MLM prediction of hERG related cardiotoxicity, by explicitly taking into account channel agonism as well as inhibition.

Discovery of Investigational Drug CT1812, an Antagonist of the Sigma-2 Receptor Complex for Alzheimer's Disease

An unbiased phenotypic neuronal assay was developed to measure the synaptotoxic effects of soluble Aβ oligomers. A collection of CNS druglike small molecules prepared by conditioned extraction was screened. Compounds that prevented and reversed synaptotoxic effects of Aβ oligomers in neurons were discovered to bind to the sigma-2 receptor complex. Select development compounds displaced receptor-bound Aβ oligomers, rescued synapses, and restored cognitive function in transgenic hAPP Swe/Ldn mice. Our first-in-class orally administered small molecule investigational drug 7 (CT1812) has been advanced to Phase II clinical studies for Alzheimer's disease.

Capsule Networks Showed Excellent Performance in the Classification of hERG Blockers/Nonblockers

Capsule networks (CapsNets), a new class of deep neural network architectures proposed recently by Hinton et al., have shown a great performance in many fields, particularly in image recognition and natural language processing. However, CapsNets have not yet been applied to drug discovery-related studies. As the first attempt, we in this investigation adopted CapsNets to develop classification models of hERG blockers/nonblockers; drugs with hERG blockade activity are thought to have a potential risk of cardiotoxicity. Two capsule network architectures were established: convolution-capsule network (Conv-CapsNet) and restricted Boltzmann machine-capsule networks (RBM-CapsNet), in which convolution and a restricted Boltzmann machine (RBM) were used as feature extractors, respectively. Two prediction models of hERG blockers/nonblockers were then developed by Conv-CapsNet and RBM-CapsNet with the Doddareddy's training set composed of 2,389 compounds. The established models showed excellent performance in an independent test set comprising 255 compounds, with prediction accuracies of 91.8 and 92.2% for Conv-CapsNet and RBM-CapsNet models, respectively. Various comparisons were also made between our models and those developed by other machine learning methods including deep belief network (DBN), convolutional neural network (CNN), multilayer perceptron (MLP), support vector machine (SVM), k-nearest neighbors (kNN), logistic regression (LR), and LightGBM, and with different training sets. All the results showed that the models by Conv-CapsNet and RBM-CapsNet are among the best classification models. Overall, the excellent performance of capsule networks achieved in this investigation highlights their potential in drug discovery-related studies.

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

The hERG (human ether-a-go-go-related gene) encoded potassium ion (K⁺) channel plays a major role in cardiac repolarization. Drug-induced blockade of hERG has been a major cause of potentially lethal ventricular tachycardia termed Torsades de Pointes (TdPs). Therefore, we presented a pharmacoinformatics strategy using combined ligand and structure based models for the prediction of hERG inhibition potential (IC₅₀) of new chemical entities (NCEs) during early stages of drug design and development. Integrated GRid-INdependent Descriptor (GRIND) models, and lipophilic efficiency (LipE), ligand efficiency (LE) guided template selection for the structure based pharmacophore models have been used for virtual screening and subsequent hERG activity (pIC₅₀) prediction of identified hits. Finally selected two hits were experimentally evaluated for hERG inhibition potential (pIC₅₀) using whole cell patch clamp assay. Overall, our results demonstrate a difference of less than ±1.6 log unit between experimentally determined and predicted hERG inhibition potential (IC₅₀) of the selected hits. This revealed predictive ability and robustness of our models and could help in correctly rank the potency order (lower μM to higher nM range) against hERG.

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.

In silico prediction of hERG inhibition

The voltage-gated potassium channel encoded by hERG carries a delayed rectifying potassium current (IKr) underlying repolarization of the cardiac action potential. Pharmacological blockade of the hERG channel results in slowed repolarization and therefore prolongation of action potential duration and an increase in the QT interval as measured on an electrocardiogram. Those are possible to cause sudden death, leading to the withdrawals of many drugs, which is the reason for hERG screening. Computational in silico prediction models provide a rapid, economic way to screen compounds during early drug discovery. In this review, hERG prediction models are classified as 2D and 3D quantitative structure–activity relationship models, pharmacophore models, classification models, and structure based models (using homology models of hERG).

A human ether-á-go-go-related (hERG) ion channel atomistic model generated by long supercomputer molecular dynamics simulations and its use in predicting drug cardiotoxicity

Acquired cardiac long QT syndrome (LQTS) is a frequent drug-induced toxic event that is often caused through blocking of the human ether-á-go-go-related (hERG) K(+) ion channel. This has led to the removal of several major drugs post-approval and is a frequent cause of termination of clinical trials. We report here a computational atomistic model derived using long molecular dynamics that allows sensitive prediction of hERG blockage. It identified drug-mediated hERG blocking activity of a test panel of 18 compounds with high sensitivity and specificity and was experimentally validated using hERG binding assays and patch clamp electrophysiological assays. The model discriminates between potent, weak, and non-hERG blockers and is superior to previous computational methods. This computational model serves as a powerful new tool to predict hERG blocking thus rendering drug development safer and more efficient. As an example, we show that a drug that was halted recently in clinical development because of severe cardiotoxicity is a potent inhibitor of hERG in two different biological assays which could have been predicted using our new computational model.

Novel Bayesian classification models for predicting compounds blocking hERG potassium channels

AIM

A large number of drug-induced long QT syndromes are ascribed to blockage of hERG potassium channels. The aim of this study was to construct novel computational models to predict compounds blocking hERG channels.

METHODS

Doddareddy's hERG blockage data containing 2644 compounds were used, which divided into training (2389) and test (255) sets. Laplacian-corrected Bayesian classification models were constructed using Discovery Studio. The models were internally validated with the training set of compounds, and then applied to the test set for validation. Doddareddy's experimentally validated dataset with 60 compounds was used for external test set validation.

RESULTS

A Bayesian classification model considering the effects of four molecular properties (Mw, PPSA, ALogP and pKa_basic) as well as extended-connectivity fingerprints (ECFP_14) exhibited a global accuracy (91%), parameter sensitivity (90%) and specificity (92%) in the test set validation, and a global accuracy (58%), parameter sensitivity (61%) and specificity (57%) in the external test set validation.

CONCLUSION

The novel model is better than those in the literatures for predicting compounds blocking hERG channels, and can be used for large-scale prediction.

Recent approaches to medicinal chemistry and therapeutic potential of dipeptidyl peptidase-4 (DPP-4) inhibitors

Dipeptidyl peptidase-4 (DPP-4) is one of the widely explored novel targets for Type 2 diabetes mellitus (T2DM) currently. Research has been focused on the strategy to preserve the endogenous glucagon like peptide (GLP)-1 activity by inhibiting the DPP-4 action. The DPP-4 inhibitors are weight neutral, well tolerated and give better glycaemic control over a longer duration of time compared to existing conventional therapies. The journey of DPP-4 inhibitors in the market started from the launch of sitagliptin in 2006 to latest drug teneligliptin in 2012. This review is mainly focusing on the recent medicinal aspects and advancements in the designing of DPP-4 inhibitors with the therapeutic potential of DPP-4 as a target to convey more clarity in the diffused data.

hERG K+ Channels: Structure, Function, and Clinical Significance

The human ether-a-go-go related gene (hERG) encodes the pore-forming subunit of the rapid component of the delayed rectifier K+ channel, Kv11.1, which are expressed in the heart, various brain regions, smooth muscle cells, endocrine cells, and a wide range of tumor cell lines. However, it is the role that Kv11.1 channels play in the heart that has been best characterized, for two main reasons. First, it is the gene product involved in chromosome 7-associated long QT syndrome (LQTS), an inherited disorder associated with a markedly increased risk of ventricular arrhythmias and sudden cardiac death. Second, blockade of Kv11.1, by a wide range of prescription medications, causes drug-induced QT prolongation with an increase in risk of sudden cardiac arrest. In the first part of this review, the properties of Kv11.1 channels, including biogenesis, trafficking, gating, and pharmacology are discussed, while the second part focuses on the pathophysiology of Kv11.1 channels.

Analysis of van der Waals surface area properties for human ether-a-go-go-related gene blocking activity: computational study on structurally diverse compounds

In the present investigation, a computational analysis was performed on a data set comprised of human ether-a-go-go-related gene (hERG) blockers (triethanolamine, 1,3-thiazol-2-yl and tetrasubstituted imidazoline derivatives) in order to investigate the structural features required to reduce the hERG-induced cardiotoxicity problems in an early stage of drug discovery. The results derived from the quantitative structure-activity relationship (QSAR) analysis showed that the volume, surface area and shape descriptors (vsurf_) contributed significantly in all the models. This reveals that the hydrogen-bonding and hydrophilicity properties (vsurf_HB1, vsurf_CW4 and a_acc) on the van der Waals (vdW) surface of the molecule is negatively contributed for the hERG blocking activity and the hydrophobic property (vsurf_D6) and the total polar volume (vsurf_Wp2) on the vdW surface of the molecule are favourable for the activity. Further, the pharmacophore analysis also shows that the Aro/Hyd/Acc contour is one of the important biophore sites for the hERG blocking activity. This suggests that the presence of aromatic, hydrophobic and hydrogen-bonding groups in the molecules is favourable for interaction. In comparison with our earlier works (explaining the role of topological and hydrophobicity properties for the hERG blocking activity), these studies provided additional information on the importance of vdW surface area properties for the hERG blocking activity. These results can be used with other molecular modelling studies for the design of novel molecules that are free of cardiotoxicity.

A critical assessment of combined ligand- and structure-based approaches to HERG channel blocker modeling

Blockade of human ether-à-go-go related gene (hERG) channel prolongs the duration of the cardiac action potential and is a common reason for drug failure in preclinical safety trials. Therefore, it is of great importance to develop robust in silico tools to predict potential hERG blockers in the early stages of drug discovery and development. Herein we described comprehensive approaches to assess the discrimination of hERG-active and -inactive compounds by combining quantitative structure-activity relationship (QSAR) modeling, pharmacophore analysis, and molecular docking. Our consensus models demonstrated high-predictive capacity and improved enrichment and could correctly classify 91.8% of 147 hERG blockers from 351 inactives. To further enhance our modeling effort, hERG homology models were constructed, and molecular docking studies were conducted, resulting in high correlations (R² = 0.81) between predicted and experimental pIC₅₀s. We expect our unique models can be applied to efficient screening for hERG blockades, and our extensive understanding of the hERG-inhibitor interactions will facilitate the rational design of drugs devoid of hERG channel activity and hence with reduced cardiac toxicities.

Rank order entropy: why one metric is not enough

The use of Quantitative Structure-Activity Relationship models to address problems in drug discovery has a mixed history, generally resulting from the misapplication of QSAR models that were either poorly constructed or used outside of their domains of applicability. This situation has motivated the development of a variety of model performance metrics (r(2), PRESS r(2), F-tests, etc.) designed to increase user confidence in the validity of QSAR predictions. In a typical workflow scenario, QSAR models are created and validated on training sets of molecules using metrics such as Leave-One-Out or many-fold cross-validation methods that attempt to assess their internal consistency. However, few current validation methods are designed to directly address the stability of QSAR predictions in response to changes in the information content of the training set. Since the main purpose of QSAR is to quickly and accurately estimate a property of interest for an untested set of molecules, it makes sense to have a means at hand to correctly set user expectations of model performance. In fact, the numerical value of a molecular prediction is often less important to the end user than knowing the rank order of that set of molecules according to their predicted end point values. Consequently, a means for characterizing the stability of predicted rank order is an important component of predictive QSAR. Unfortunately, none of the many validation metrics currently available directly measure the stability of rank order prediction, making the development of an additional metric that can quantify model stability a high priority. To address this need, this work examines the stabilities of QSAR rank order models created from representative data sets, descriptor sets, and modeling methods that were then assessed using Kendall Tau as a rank order metric, upon which the Shannon entropy was evaluated as a means of quantifying rank-order stability. Random removal of data from the training set, also known as Data Truncation Analysis (DTA), was used as a means for systematically reducing the information content of each training set while examining both rank order performance and rank order stability in the face of training set data loss. The premise for DTA ROE model evaluation is that the response of a model to incremental loss of training information will be indicative of the quality and sufficiency of its training set, learning method, and descriptor types to cover a particular domain of applicability. This process is termed a "rank order entropy" evaluation or ROE. By analogy with information theory, an unstable rank order model displays a high level of implicit entropy, while a QSAR rank order model which remains nearly unchanged during training set reductions would show low entropy. In this work, the ROE metric was applied to 71 data sets of different sizes and was found to reveal more information about the behavior of the models than traditional metrics alone. Stable, or consistently performing models, did not necessarily predict rank order well. Models that performed well in rank order did not necessarily perform well in traditional metrics. In the end, it was shown that ROE metrics suggested that some QSAR models that are typically used should be discarded. ROE evaluation helps to discern which combinations of data set, descriptor set, and modeling methods lead to usable models in prioritization schemes and provides confidence in the use of a particular model within a specific domain of applicability.

In vivo and in vitro characterization of a first-in-class novel azole analog that targets pregnane X receptor activation

The pregnane X receptor (PXR) is a master regulator of xenobiotic clearance and is implicated in deleterious drug interactions (e.g., acetaminophen hepatotoxicity) and cancer drug resistance. However, small-molecule targeting of this receptor has been difficult; to date, directed synthesis of a relatively specific PXR inhibitor has remained elusive. Here we report the development and characterization of a first-in-class novel azole analog [1-(4-(4-(((2R,4S)-2-(2,4-difluorophenyl)-2-methyl-1,3-dioxolan-4-yl)methoxy)phenyl)piperazin-1-yl)ethanone (FLB-12)] that antagonizes the activated state of PXR with limited effects on other related nuclear receptors (i.e., liver X receptor, farnesoid X receptor, estrogen receptor α, peroxisome proliferator-activated receptor γ, and mouse constitutive androstane receptor). We investigated the toxicity and PXR antagonist effect of FLB-12 in vivo. Compared with ketoconazole, a prototypical PXR antagonist, FLB-12 is significantly less toxic to hepatocytes. FLB-12 significantly inhibits the PXR-activated loss of righting reflex to 2,2,2-tribromoethanol (Avertin) in vivo, abrogates PXR-mediated resistance to 7-ethyl-10-hydroxycamptothecin (SN-38) in colon cancer cells in vitro, and attenuates PXR-mediated acetaminophen hepatotoxicity in vivo. Thus, relatively selective targeting of PXR by antagonists is feasible and warrants further investigation. This class of agents is suitable for development as chemical probes of PXR function as well as potential PXR-directed therapeutics.

Development, validation, and use of quantitative structure-activity relationship models of 5-hydroxytryptamine (2B) receptor ligands to identify novel receptor binders and putative valvulopathic compounds among common drugs

Some antipsychotic drugs are known to cause valvular heart disease by activating serotonin 5-HT(2B) receptors. We have developed and validated binary classification QSAR models capable of predicting potential 5-HT(2B) actives. The classification accuracies of the models built to discriminate 5-HT(2B) actives from the inactives were as high as 80% for the external test set. These models were used to screen in silico 59,000 compounds included in the World Drug Index, and 122 compounds were predicted as actives with high confidence. Ten of them were tested in radioligand binding assays and nine were found active, suggesting a success rate of 90%. All validated actives were then tested in functional assays, and one compound was identified as a true 5-HT(2B) agonist. We suggest that the QSAR models developed in this study could be used as reliable predictors to flag drug candidates that are likely to cause valvulopathy.

hERG-related drug toxicity and models for predicting hERG liability and QT prolongation

BACKGROUND

hERG K(+) channels have been recognized as a primary antitarget in safety pharmacology. Their blockade, caused by several drugs with different therapeutic indications, may lead to QT prolongation and, eventually, to potentially fatal arrhythmia, namely torsade de pointes. Therefore, a number of preclinical models have been developed to predict hERG liability early in the drug development process.

OBJECTIVE

The aim of this review is to outline the present state of the art on drug-induced hERG blockade, providing insights on the predictive value of in vitro and in silico models for hERG liability.

METHODS

On the basis of latest reports, high-throughput preclinical models have been discussed outlining advantages and limitations.

CONCLUSION

Although no single model has an absolute value, an integrated risk assessment is recommended to predict the pro-arrhythmic risk of a given drug. This prediction requires expertise from different areas and should encompass emerging issues such as interference with hERG trafficking and QT shortening.

How not to develop a quantitative structure-activity or structure-property relationship (QSAR/QSPR)

Although thousands of quantitative structure-activity and structure-property relationships (QSARs/QSPRs) have been published, as well as numerous papers on the correct procedures for QSAR/QSPR analysis, many analyses are still carried out incorrectly, or in a less than satisfactory manner. We have identified 21 types of error that continue to be perpetrated in the QSAR/QSPR literature, and each of these is discussed, with examples (including some of our own). Where appropriate, we make recommendations for avoiding errors and for improving and enhancing QSAR/QSPR analyses.

Topological polar surface area: a useful descriptor in 2D-QSAR

Topological polar surface area (TPSA), which makes use of functional group contributions based on a large database of structures, is a convenient measure of the polar surface area that avoids the need to calculate ligand 3D structure or to decide which is the relevant biological conformation or conformations. We demonstrate the utility of TPSA in 2D-QSAR for 14 sets of diverse pharmacological activity data. Even though a large pool of reports showing the importance of the classic 2D descriptors such as calculated logP (ClogP) and calculated molar refractivity (CMR) exists in the 2D-QSAR literature, this is the first report to demonstrate the value of TPSA as a relevant descriptor applicable to a large, structurally and pharmacologically diverse set of classes of compounds. We also address the limitations of applicability of this descriptor for 2D-QSAR analysis. We observed a negative correlation of TPSA with activity data for anticancer alkaloids, MT1 and MT2 agonists, MAO-B and tumor necrosis factor-alpha inhibitors and a positive correlation with inhibitory activity data for telomerase, PDE-5, GSK-3, DNA-PK, aromatase, malaria, trypanosomatids and CB2 agonists.

Pharmacokinetics in drug discovery

The aim of this current review is to summarize the present status of pharmacokinetics in Drug Discovery. The review is structured into four sections. The first section is a general overview of what we understand by pharmacokinetics and the different LADMET aspects: Liberation, Absorption, Distribution, Metabolism, Excretion, and Toxicity. The second section highlights the different computational or in silico approaches to estimate/predict one or several aspects of the pharmacokinetic profile of a discovery lead compound. The third section discusses the most commonly used in vitro methodologies. The fourth and last section examines the various approaches employed towards the pharmacokinetic assessment of discovery molecules; including all the LADME processes, discussing the different mathematical methodologies available to establish the PK profile of a test compound; what the main differences are and what should be the criteria for using one or another mathematical approach. The major conclusion of this review is that the use of the appropriate preclinical assays has a key role in the long-term viability of a pharmaceutical company since applying the right tools early in discovery will play a key role in determining the company's ability to discover novel safe and effective therapeutics to patients as quickly as possible.

A novel structure-based virtual screening model for the hERG channel blockers

The hERG potassium channel is a key effector of cardiac repolarization and the blockade of this channel could cause arrhythmia. Thus, hERG channel blockade plays an important role for the potential pro-arrhythmic liability. In this report, binding of blockers to the hERG potassium channel is investigated using a combination of homology modeling, molecular docking, and molecular simulations, where blockade activities are evaluated using the linear regression model of GoldScore fitness. This structure-based virtual screening model is able to estimate the pIC(50) value of a wide range of ligands for the hERG potassium channel. The docked poses for ligands are also consistent with published mutation. Therefore, this model for the prediction of hERG channel blockade has the potential to provide cost-effective virtual screening tools for the evaluation of the cardiac liability of new chemical entities.