The pharmacophore hypotheses of IKr potassium channel blockers: novel class III antiarrhythmic agents

Rollover Cyclometalation vs Nitrogen Coordination in Tetrapyridyl Anticancer Gold(III) Complexes: Effect on Protein Interaction and Toxicity

In this work, a pair of gold(III) complexes derived from the analogous tetrapyridyl ligands H₂biqbpy1 and H₂biqbpy2 was prepared: the rollover, bis-cyclometalated [Au(biqbpy1)Cl ([1]Cl) and its isomer [Au(biqbpy2)Cl ([2]Cl). In [1]⁺, two pyridyl rings coordinate to the metal via a Au-C bond (C^∧N^∧N^∧C coordination) and the two noncoordinated amine bridges of the ligand remain protonated, while in [2]⁺ all four pyridyl rings of the ligand coordinate to the metal via a Au-N bond (N^∧N^∧N^∧N coordination), but both amine bridges are deprotonated. As a result, both complexes are monocationic, which allowed comparison of the sole effect of cyclometalation on the chemistry, protein interaction, and anticancer properties of the gold(III) compounds. Due to their identical monocationic charge and similar molecular shape, both complexes [1]Cl and [2]Cl displaced reference radioligand [³H]dofetilide equally well from cell membranes expressing the K_v11.1 (hERG) potassium channel, and more so than the tetrapyridyl ligands H₂biqbpy1 and H₂biqbpy2. By contrast, cyclometalation rendered [1]Cl coordinatively stable in the presence of biological thiols, while [2]Cl was reduced by a millimolar concentration of glutathione into metastable Au(I) species releasing the free ligand H₂biqbpy2 and TrxR-inhibiting Au⁺ ions. The redox stability of [1]Cl dramatically decreased its thioredoxin reductase (TrxR) inhibition properties, compared to [2]Cl. On the other hand, unlike [2]Cl, [1]Cl aggregated into nanoparticles in FCS-containing medium, which resulted in much more efficient gold cellular uptake. [1]Cl had much more selective anticancer properties than [2]Cl and cisplatin, as it was almost 10 times more cytotoxic to human cancer cells (A549, A431, A375, and MCF7) than to noncancerous cells (MRC5). Mechanistic studies highlight the strikingly different mode of action of the two compounds: while for [1]Cl high gold cellular uptake, nuclear DNA damage, and interaction with hERG may contribute to cell killing, for [2]Cl extracellular reduction released TrxR-inhibiting Au⁺ ions that were taken up in minute amounts in the cytosol, and a toxic tetrapyridyl ligand also capable of binding to hERG. These results demonstrate that bis-cyclometalation is an appealing method to improve the redox stability of Au(III) compounds and to develop gold-based cytotoxic compounds that do not rely on TrxR inhibition to kill cancer cells.

Environment-Sensitive Fluorescent Probe for the Human Ether-a-go-go-Related Gene Potassium Channel

A novel environment-sensitive probe S2 with turn-on switch for Human Ether-a-go-go-Related Gene (hERG) potassium channel was developed herein. After careful evaluation, this fluorescent probe showed high binding affinity with hERG potassium channel with an IC₅₀ value of 41.65 nM and can be well applied to hERG channel imaging or cellular distribution study for hERG channel blockers. Compared with other imaging techniques, such as immunofluorescence and fluorescent protein-based approaches, this method is convenient and affordable, especially since a washing procedure is not needed. Meanwhile, this environment-sensitive turn-on design strategy may provide a good example for the probe development for these targets that have no reactive or catalytic activity.

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

Global analysis reveals families of chemical motifs enriched for HERG inhibitors

Promiscuous inhibition of the human ether-à-go-go-related gene (hERG) potassium channel by drugs poses a major risk for life threatening arrhythmia and costly drug withdrawals. Current knowledge of this phenomenon is derived from a limited number of known drugs and tool compounds. However, in a diverse, naïve chemical library, it remains unclear which and to what degree chemical motifs or scaffolds might be enriched for hERG inhibition. Here we report electrophysiology measurements of hERG inhibition and computational analyses of >300,000 diverse small molecules. We identify chemical ‘communities’ with high hERG liability, containing both canonical scaffolds and structurally distinctive molecules. These data enable the development of more effective classifiers to computationally assess hERG risk. The resultant predictive models now accurately classify naïve compound libraries for tendency of hERG inhibition. Together these results provide a more complete reference map of characteristic chemical motifs for hERG liability and advance a systematic approach to rank chemical collections for cardiotoxicity risk.

Strategies to reduce HERG K+ channel blockade. Exploring heteroaromaticity and rigidity in novel pyridine analogues of dofetilide

Drug-induced blockade of the human ether-a-go-go-related gene K(+) channel (hERG) represents one of the major antitarget concerns in pharmaceutical industry. SAR studies of this ion channel have shed light on the structural requirements for hERG interaction but most importantly may reveal drug design principles to reduce hERG affinity. In the present study, a novel library of neutral and positively charged heteroaromatic derivatives of the class III antiarrhythmic agent dofetilide was synthesized and assessed for hERG affinity in radioligand binding and manual patch clamp assays. Structural modifications of the pyridine moiety, side chain, and peripheral aromatic moieties were evaluated, thereby revealing approaches for reducing hERG binding affinity. In particular, we found that the extra rigidity imposed close to the positively charged pyridine moiety can be very efficient in decreasing hERG affinity.

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.

Non-basic ligands for aminergic GPCRs: the discovery and development diaryl sulfones as selective, orally bioavailable 5-HT2A receptor antagonists for the treatment of sleep disorders

Scaffold hopping from a non-basic series of 5-HT(2A) receptor antagonists developed in-house that possessed reduced activity in vivo enabled the discovery of a novel series of diaryl sulfones that gave excellent occupancy on oral dosing. Not only does this work further demonstrate that oral bioavailability of a given series can be enhanced by improving physicochemical parameters such as log P, but it corroborates the growing evidence that a protonated amine is not essential for affinity at aminergic GPCRs.

Highly potent, non-basic 5-HT6 ligands. Site mutagenesis evidence for a second binding mode at 5-HT6 for antagonism

A series of 5-HT(6) ligands derived from (R)-1-(amino)methyl-6-(phenyl)sulfonyltetralin was prepared that yielded several non-basic analogs having sub-nanomolar affinity. Ligand structure-activity relationships, receptor point mutation studies, and molecular modeling of these novel ligands all combined to reveal a new alternative binding mode to 5-HT(6) for antagonism.

Collation, assessment and analysis of literature in vitro data on hERG receptor blocking potency for subsequent modeling of drugs' cardiotoxic properties

The assessment of the torsadogenic potency of a new chemical entity is a crucial issue during lead optimization and the drug development process. It is required by the regulatory agencies during the registration process. In recent years, there has been a considerable interest in developing in silico models, which allow prediction of drug-hERG channel interaction at the early stage of a drug development process. The main mechanism underlying an acquired QT syndrome and a potentially fatal arrhythmia called torsades de pointes is the inhibition of potassium channel encoded by hERG (the human ether-a-go-go-related gene). The concentration producing half-maximal block of the hERG potassium current (IC(50)) is a surrogate marker for proarrhythmic properties of compounds and is considered a test for cardiac safety of drugs or drug candidates. The IC(50) values, obtained from data collected during electrophysiological studies, are highly dependent on experimental conditions (i.e. model, temperature, voltage protocol). For the in silico models' quality and performance, the data quality and consistency is a crucial issue. Therefore the main objective of our work was to collect and assess the hERG IC(50) data available in accessible scientific literature to provide a high-quality data set for further studies.

Targeting the fatty acid biosynthesis enzyme, beta-ketoacyl-acyl carrier protein synthase III (PfKASIII), in the identification of novel antimalarial agents

The importance of fatty acids to the human malaria parasite, Plasmodium falciparum, and differences due to a type I fatty acid synthesis (FAS) pathway in the parasite, make it an attractive drug target. In the present study, we developed and a utilized a pharmacophore to select compounds for testing against PfKASIII, the initiating enzyme of FAS. This effort identified several PfKASIII inhibitors that grouped into various chemical classes of sulfides, sulfonamides, and sulfonyls. Approximately 60% of the submicromolar inhibitors of PfKASIII inhibited in vitro growth of the malaria parasite. These compounds inhibited both drug sensitive and resistant parasites and testing against a mammalian cell line revealed an encouraging in vitro therapeutic index for the most active compounds. Docking studies into the active site of PfKASIII suggest a potential binding mode that exploits amino acid residues at the mouth of the substrate tunnel.

A composite model for HERG blockade

hERG blockade is one of the major toxicological problems in lead structure optimization. Reliable ligand-based in silico models for predicting hERG blockade therefore have considerable potential for saving time and money, as patch-clamp measurements are very expensive and no crystal structures of the hERG-encoded channel are available. Herein we present a predictive QSAR model for hERG blockade that differentiates between specific and nonspecific binding. Specific binders are identified by preliminary pharmacophore scanning. In addition to descriptor-based models for the compounds selected as hitting one of two different pharmacophores, we also use a model for nonspecific binding that reproduces blocking properties of molecules that do not fit either of the two pharmacophores. PLS and SVR models based on interpretable quantum mechanically derived descriptors on a literature dataset of 113 molecules reach overall R(2) values between 0.60 and 0.70 for independent validation sets and R(2) values between 0.39 and 0.76 after partitioning according to the pharmacophore search for the test sets. Our findings suggest that hERG blockade may occur through different types of binding, so that several different models may be necessary to assess hERG toxicity.

Discovery of novel and orally active NR2B-selective N-methyl-D-aspartate (NMDA) antagonists, pyridinol derivatives with reduced HERG binding affinity

Novel NR2B antagonists with an amide tether were found by an approach to avoid pharmacophoric similarity to dofetilide. Structure-activity relationship investigation led to N-[cis-4-hydroxy-4-(5-hydroxypyridin-2-yl)cyclohexyl]-3-henylpropanamide as an orally active NR2B-subtype selective N-methyl-D-aspartate (NMDA) receptor antagonist with very weak HERG (human ether-a-go-go related gene) binding (IC(50)> 30 microM). This compound exhibited potent in vivo anti-allodynic activity in the mouse partial sciatic nerve ligation (PSL) model (minimum effective dose=10 mg/kg, po).

Exploring QSTR and toxicophore of hERG K+ channel blockers using GFA and HypoGen techniques

Predictive quantitative structure-toxicity and toxicophore models were developed for a diverse series of hERG K+ channel blockers, acting as anti-arrhythmic agents using QSAR+ module in Cerius2 and HypoGen module in Catalyst software, respectively. The 2D-QSTR analysis has been performed on a dataset of 68 molecules carefully selected from literature for which IC50 values measured on hERG K+ channels expressed in mammalian cells lines using the voltage patch clamp assay technique were reported. Their biological data, expressed as IC50, spanned from 7.0nM to 1.4mM, with 7 orders difference. Several types of descriptors including electrotopological, thermodynamic, ADMET, graph theoretical (topological and information content) were used to derive a quantitative relationship between the channel blockers and its physico-chemical properties. Statistically significant QSTR model was obtained using genetic function approximation methodology, having seven descriptors, with a correlation coefficient (r2) of 0.837, cross-validated correlation coefficient (q2) of 0.776 and predictive correlation coefficient (r2 pred) of 0.701, indicating the robustness of the model. Toxicophore model generated using HypoGen module in Catalyst, on these datasets, showed three important features for hERG K+ channel blockers, (i) hydrophobic group (HP), (ii) ring aromatic group (RA) and (iii) hydrogen bond acceptor lipid group (HBAl). The most predictive hypothesis (Hypo 1), consisting of these three features had a best correlation coefficient of 0.820, a low rms deviation of 1.740, and a high cost difference of 113.50, which represents a true correlation and a good predictivity. The hypothesis, Hypo 1 was validated by a test set consisting of 12 molecules and by a cross-validation of 95% confidence level. Accordingly, our 2D-QSTR and toxicophore model has strong predictivity to identify structurally diverse hERG K+ channel blockers with desired biological activity. These models provide a useful framework for understanding binding, and gave structural insight into the specific protein-ligand interactions responsible for affinity, and how one might modify any given structure to mitigate binding.

Understanding hERG inhibition with QSAR models based on a one-dimensional molecular representation

Blockage of the potassium channel encoded by the human ether-a-go-go related gene (hERG) is well understood to be the root cause of the cardio-toxicity of numerous approved and investigational drugs. As such, a cascade of in vitro and in vivo assays have been developed to filter compounds with hERG inhibitory activity. Quantitative structure activity relationship (QSAR) models are used at the very earliest part of this cascade to eliminate compounds that are likely to have this undesirable activity prior to synthesis. Here a new QSAR technique based on the one-dimensional representation is described in the context of the development of a model to predict hERG inhibition. The model is shown to perform close to the limits of the quality of the data used for model building. In order to make optimal use of the available data, a general robust mathematical scheme was developed and is described to simultaneously incorporate quantitative data, such as IC50 = 50 nM, and qualitative data, such as inactive or IC50 > 30 microM into QSAR models without discarding any experimental information.

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.

Determination of hERG channel blockers using a decision tree

A decision tree approach for the in silico prediction of Torsade de Pointes (TdP)-causing drugs is presented. As TdP is frequently associated with QT-interval prolongation due to inhibition of the rapid activating delayed rectifier potassium channel in the heart (hERG channel), the properties of such blockers were investigated by molecular modeling and semi-empirical AM1 molecular orbital calculations. In addition, we derived a pharmacophoric SMARTS string using structural information from high affinity compounds. A corresponding search in the PubChem database identified several compounds that exhibit QT-interval prolonging activity that were not among our data set. This SMARTS string furthermore showed to be the most significant descriptor in the decision tree approach from which guidelines for the design of safe compounds are suggested.

A new class of selective, non-basic 5-HT2A receptor antagonists

Based on an existing series of 5-HT2A receptor ligands containing a basic nitrogen, we designed a non-basic lead that had reduced affinity for both the 5-HT2A receptor and the IKr potassium channel. The present paper describes the development of this lead to a novel series of non-basic piperidine sulfonamides and amides that have high affinity for the 5-HT2A receptor, whilst maintaining excellent selectivity over off target activities such as the IKr channel. This work has shown that the proposed pharmacophore model for the 5-HT2A receptor which suggests that a basic nitrogen is required for the binding of ligands is questionable.

Pharmacophore identification of α1A-adrenoceptor antagonists

A chemical feature based pharmacophore model was developed for alpha(1A)-adrenoceptor antagonists by HypoGen module implemented in catalyst software package. The best scoring pharmacophore hypothesis, Hypo1, consisted of four important chemical features (one positive ion, one hydrogen-bond donor, one aromatic ring, and one hydrophobic group). The results of our study provide a valuable tool in designing new leads with desired biological activity by virtual screening.