Troubleshooting problems with in vitro screening of drugs for QT interval prolongation using HERG K+ channels expressed in mammalian cell lines and Xenopus oocytes

Experimentally validated modeling of dynamic drug-hERG channel interactions reproducing the binding mechanisms and its importance in action potential duration

BACKGROUND AND OBJECTIVE

Assessment of drug cardiotoxicity is critical in the development of new compounds and modeling of drug-binding dynamics to hERG can improve early cardiotoxicity assessment. We previously developed a methodology to generate Markovian models reproducing preferential state-dependent binding properties, trapping dynamics and the onset of IKr block using simple voltage clamp protocols. Here, we test this methodology with real IKr blockers and investigate the impact of drug dynamics on action potential prolongation.

METHODS

Experiments were performed on HEK cells stably transfected with hERG and using the Nanion SyncroPatch 384i. Three protocols, P-80, P0 and P 40, were applied to obtain the experimental data from the drugs and the Markovian models were generated using our pipeline. The corresponding static models were also generated and a modified version of the O´Hara-Rudy action potential model was used to simulate the action potential duration.

RESULTS

The experimental Hill plots and the onset of IKr block of ten compounds were obtained using our voltage clamp protocols and the models generated successfully mimicked these experimental data, unlike the CiPA dynamic models. Marked differences in APD prolongation were observed when drug effects were simulated using the dynamic models and the static models.

CONCLUSIONS

These new dynamic models of ten well-known IKr blockers constitute a validation of our methodology to model dynamic drug-hERG channel interactions and highlight the importance of state-dependent binding, trapping dynamics and the time-course of IKr block to assess drug effects even at the steady-state.

Tricyclic hydrocarbon fluorene attenuates ventricular ionic currents and pressure development in the navaga cod

The release of polycyclic aromatic hydrocarbons (PAHs) into the environment due to oil and diesel fuel spills is a serious threat to Arctic fish populations. PAHs produce multiple toxic effects in fish, but disturbance of electrical and contractile activity of the heart seems to be the most negative effect. Our study focused on the effects of fluorene, a tricyclic PAH resembling the well-investigated tricyclic phenanthrene, on major ionic currents and action potential (AP) waveform in isolated ventricular myocytes and on contractile activity in isolated whole hearts of polar navaga cod (Eleginus nawaga). Among the studied currents, the repolarizing rapid delayed rectifier K+ current IKr demonstrated the highest sensitivity to fluorene with IC50 of 0.54 μM. The depolarizing inward currents, INa and ICaL, were inhibited with 10 μM fluorene by 20.2 ± 2.8 % and 27.9 ± 8.4 %, respectively, thereby being much less sensitive to fluorene than IKr. Inward rectifier IK1 current was insensitive to fluorene (up to 10 μM). While 3 μM fluorene prolonged APs, 10 μM also slowed the AP upstroke. Resting membrane potential was not affected by any tested concentrations. In isolated heart experiments 10 μM fluorene caused modest depression of ventricular contractile activity. Thus, we have demonstrated that fluorene, a tricyclic PAH present in high quantities in crude oil, strongly impacts electrical activity with only slight effects on contractile activity in the heart of the polar fish, the navaga cod.

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.

Automatic modeling of dynamic drug-hERG channel interactions using three voltage protocols and machine learning techniques: A simulation study

BACKGROUND

Assessment of drug cardiac safety is critical in the development of new compounds and is commonly addressed by evaluating the half-maximal blocking concentration of the potassium human ether-à-go-go related gene (hERG) channels. However, recent works have evidenced that the modelling of drug-binding dynamics to hERG can help to improve early cardiac safety assessment. Our goal is to develop a methodology to automatically generate Markovian models of the drug-hERG channel interactions.

METHODS

The training and the test sets consisted of 20800 and 5200 virtual drugs, respectively, distributed into 104 groups with different affinities and kinetics to the conformational states of the channel. In our system, drugs may bind to any state (individually or simultaneously), with different degrees of preference for a conformational state and the change of the conformational state of the drug bound channels may be restricted or allowed. To model such a wide range of possibilities, 12 Markovian chains are considered. Our approach uses the response of the drugs to our three previously developed voltage clamp protocols, which enhance the differences in the probabilities of occupying a certain conformational state of the channel (open, closed and inactivated). The computing tool is comprised of a classifier and a parameter optimizer and uses linear interpolation, support vector machines and a simplex method for function minimization.

RESULTS

We propose a novel methodology that automatically generates dynamic drug models using Markov model formulations and that elucidates the states where the drug binds and unbinds and the preferential binding state using data obtained from simple voltage clamp protocols that captures the preferential state-dependent binding properties, the relative affinities, trapping and non-trapping dynamics and the onset of I_Kr block. Overall, the tool correctly predicted the class of 92.04% of the drugs and the model provided by the tool accurately fitted the response of the target compound, the mean accuracy being 97.53%. Moreover, generation of the dynamic model of an I_Kr blocker from its response to our voltage clamp protocols usually takes less than an hour on a common desktop computer.

CONCLUSION

Our methodology could be very useful to model and simulate dynamic drug-hERG channel interactions. It would contribute to the improvement of the preclinical assessment of the proarrhythmic risk of drugs that inhibit I_Kr and the efficacy of antiarrhythmic I_Kr blockers.

Design of New Potent and Selective Thiophene-Based KV1.3 Inhibitors and Their Potential for Anticancer Activity

Simple Summary

In this article, we describe the discovery of a new class of potent and selective thiophene-based inhibitors of the voltage-gated potassium channel K_V1.3 and their potential to induce apoptosis and inhibit proliferation. The K_V1.3 channel has only recently emerged as a molecular target in cancer therapy. The most potent K_V1.3 inhibitor 44 had an IC₅₀ K_V1.3 value of 470 nM (oocytes) and 950 nM (Ltk⁻ cells) and appropriate selectivity for other K_V channels. New K_V1.3 inhibitors significantly inhibited proliferation of Panc-1 cells and K_V1.3 inhibitor 4 induced significant apoptosis in tumor spheroids of Colo-357 cells.

Abstract

The voltage-gated potassium channel K_V1.3 has been recognized as a tumor marker and represents a promising new target for the discovery of new anticancer drugs. We designed a novel structural class of K_V1.3 inhibitors through structural optimization of benzamide-based hit compounds and structure-activity relationship studies. The potency and selectivity of the new K_V1.3 inhibitors were investigated using whole-cell patch- and voltage-clamp experiments. 2D and 3D cell models were used to determine antiproliferative activity. Structural optimization resulted in the most potent and selective K_V1.3 inhibitor 44 in the series with an IC₅₀ value of 470 nM in oocytes and 950 nM in Ltk⁻ cells. K_V1.3 inhibitor 4 induced significant apoptosis in Colo-357 spheroids, while 14, 37, 43, and 44 significantly inhibited Panc-1 proliferation.

Evaluation of electrocardiographic repolarization parameters after administration of trimethoprim-sulfadiazine, detomidine, or their combination in horses

OBJECTIVE

To determine whether administration of trimethoprim-sulfadiazine (TMS), detomidine (DET), or TMS plus DET would be associated with changes in ECG repolarization parameters in horses.

ANIMALS

9 healthy adult horses.

PROCEDURES

Each horse received 4 treatments in a blinded, randomized, crossover study design as follows: TMS, 16 to 24 mg/kg, IV; DET, 0.015 to 0.02 mg/kg, IV; TMS plus DET; and saline (0.9% NaCl) solution. Surface ECG traces were obtained over 24 hours, and repolarization parameters were measured at predefined time points after each treatment and compared with a 2-way ANOVA for repeated measures.

RESULTS

Heart rate-corrected QT intervals (QTc) were significantly increased after administration of DET (mean ± SD difference in QTc, 36.57 ± 23.07 milliseconds; increase of 7%) and TMS plus DET (44.96 ± 29.16 milliseconds; increase of 9%), compared with baseline (before treatment) values and values after administration of saline solution. Saline solution and TMS alone did not affect QTc.

CONCLUSIONS AND CLINICAL RELEVANCE

Administration of DET or TMS plus DET was associated with a significant and possibly clinically relevant prolongation of QTc, with prolongation of 7% to 9%, a range that is considered as a risk factor for the development of cardiac arrhythmias in people. Results were unexpected because DET is considered to be a safe sedative for horses.

Critical Assessment of Artificial Intelligence Methods for Prediction of hERG Channel Inhibition in the "Big Data" Era

The rise of novel artificial intelligence (AI) methods necessitates their benchmarking against classical machine learning for a typical drug-discovery project. Inhibition of the potassium ion channel, whose alpha subunit is encoded by the human ether-à-go-go-related gene (hERG), leads to a prolonged QT interval of the cardiac action potential and is a significant safety pharmacology target for the development of new medicines. Several computational approaches have been employed to develop prediction models for the assessment of hERG liabilities of small molecules including recent work using deep learning methods. Here, we perform a comprehensive comparison of hERG effect prediction models based on classical approaches (random forests and gradient boosting) and modern AI methods [deep neural networks (DNNs) and recurrent neural networks (RNNs)]. The training set (∼9000 compounds) was compiled by integrating the hERG bioactivity data from the ChEMBL database with experimental data generated from an in-house, high-throughput thallium flux assay. We utilized different molecular descriptors including the latent descriptors, which are real-value continuous vectors derived from chemical autoencoders trained on a large chemical space (>1.5 million compounds). The models were prospectively validated on ∼840 in-house compounds screened in the same thallium flux assay. The best results were obtained with the XGBoost method and RDKit descriptors. The comparison of models based only on latent descriptors revealed that the DNNs performed significantly better than the classical methods. The RNNs that operate on SMILES provided the highest model sensitivity. The best models were merged into a consensus model that offered superior performance compared to reference models from academic and commercial domains. Furthermore, we shed light on the potential of AI methods to exploit the big data in chemistry and generate novel chemical representations useful in predictive modeling and tailoring a new chemical space.

Current development of integrated web servers for preclinical safety and pharmacokinetics assessments in drug development

In drug development, preclinical safety and pharmacokinetics assessments of candidate drugs to ensure the safety profile are a must. While in vivo and in vitro tests are traditionally used, experimental determinations have disadvantages, as they are usually time-consuming and costly. In silico predictions of these preclinical endpoints have each been developed in the past decades. However, only a few web-based tools have integrated different models to provide a simple one-step platform to help researchers thoroughly evaluate potential drug candidates. To efficiently achieve this approach, a platform for preclinical evaluation must not only predict key ADMET (absorption, distribution, metabolism, excretion and toxicity) properties but also provide some guidance on structural modifications to improve the undesired properties. In this review, we organized and compared several existing integrated web servers that can be adopted in preclinical drug development projects to evaluate the subject of interest. We also introduced our new web server, Virtual Rat, as an alternative choice to profile the properties of drug candidates. In Virtual Rat, we provide not only predictions of important ADMET properties but also possible reasons as to why the model made those structural predictions. Multiple models were implemented into Virtual Rat, including models for predicting human ether-a-go-go-related gene (hERG) inhibition, cytochrome P450 (CYP) inhibition, mutagenicity (Ames test), blood-brain barrier penetration, cytotoxicity and Caco-2 permeability. Virtual Rat is free and has been made publicly available at https://virtualrat.cmdm.tw/.

When Does the IC50 Accurately Assess the Blocking Potency of a Drug?

Preclinical assessment of drug-induced proarrhythmicity is typically evaluated by the potency of the drug to block the potassium human ether-à-go-go-related gene (hERG) channels, which is currently quantified by the IC₅₀. However, channel block depends on the experimental conditions. Our aim is to improve the evaluation of the blocking potency of drugs by designing experimental stimulation protocols to measure the IC₅₀ that will help to decide whether the IC₅₀ is representative enough. We used the state-of-the-art mathematical models of the cardiac electrophysiological activity to design three stimulation protocols that enhance the differences in the probabilities to occupy a certain conformational state of the channel and, therefore, the potential differences in the blocking effects of a compound. We simulated an extensive set of 144 in silico I_Kr blockers with different kinetics and affinities to conformational states of the channel and we also experimentally validated our key predictions. Our results show that the IC₅₀ protocol dependency relied on the tested compounds. Some of them showed no differences or small differences on the IC₅₀ value, which suggests that the IC₅₀ could be a good indicator of the blocking potency in these cases. However, others provided highly protocol dependent IC₅₀ values, which could differ by even 2 orders of magnitude. Moreover, the protocols yielding the maximum IC₅₀ and minimum IC₅₀ depended on the drug, which complicates the definition of a "standard" protocol to minimize the influence of the stimulation protocol on the IC₅₀ measurement in safety pharmacology. As a conclusion, we propose the adoption of our three-protocol IC₅₀ assay to estimate the potency to block hERG in vitro. If the IC₅₀ values obtained for a compound are similar, then the IC₅₀ could be used as an indicator of its blocking potency, otherwise kinetics and state-dependent binding properties should be accounted.

The Catch-22 of Predicting hERG Blockade Using Publicly Accessible Bioactivity Data

Drug-induced inhibition of the human ether-à-go-go-related gene (hERG)-encoded potassium ion channels can lead to fatal cardiotoxicity. Several marketed drugs and promising drug candidates were recalled because of this concern. Diverse modeling methods ranging from molecular similarity assessment to quantitative structure-activity relationship analysis employing machine learning techniques have been applied to data sets of varying size and composition (number of blockers and nonblockers). In this study, we highlight the challenges involved in the development of a robust classifier for predicting the hERG end point using bioactivity data extracted from the public domain. To this end, three different modeling methods, nearest neighbors, random forests, and support vector machines, were employed to develop predictive models using different molecular descriptors, activity thresholds, and training set compositions. Our models demonstrated superior performance in external validations in comparison with those reported in the previous studies from which the data sets were extracted. The choice of descriptors had little influence on the model performance, with minor exceptions. The criteria used to filter bioactivity data, the activity threshold settings used to separate blockers from nonblockers, and the structural diversity of blockers in training data set were found to be the crucial indicators of model performance. Training sets based on a binary threshold of 1 μM/10 μM to separate blockers (IC₅₀/ K_i ≤ 1 μM) from nonblockers (IC₅₀/ K_i > 10 μM) provided superior performance in comparison with those defined using a single threshold (1 μM or 10 μM). A major limitation in using the public domain hERG activity data is the abundance of blockers in comparison with nonblockers at usual activity thresholds, since not many studies report the latter.

The tarantula toxin β/δ-TRTX-Pre1a highlights the importance of the S1-S2 voltage-sensor region for sodium channel subtype selectivity

Voltage-gated sodium (Na_V) channels are essential for the transmission of pain signals in humans making them prime targets for the development of new analgesics. Spider venoms are a rich source of peptide modulators useful to study ion channel structure and function. Here we describe β/δ-TRTX-Pre1a, a 35-residue tarantula peptide that selectively interacts with neuronal Na_V channels inhibiting peak current of hNa_V1.1, rNa_V1.2, hNa_V1.6, and hNa_V1.7 while concurrently inhibiting fast inactivation of hNa_V1.1 and rNa_V1.3. The DII and DIV S3-S4 loops of Na_V channel voltage sensors are important for the interaction of Pre1a with Na_V channels but cannot account for its unique subtype selectivity. Through analysis of the binding regions we ascertained that the variability of the S1-S2 loops between Na_V channels contributes substantially to the selectivity profile observed for Pre1a, particularly with regards to fast inactivation. A serine residue on the DIV S2 helix was found to be sufficient to explain Pre1a’s potent and selective inhibitory effect on the fast inactivation process of Na_V1.1 and 1.3. This work highlights that interactions with both S1-S2 and S3-S4 of Na_V channels may be necessary for functional modulation, and that targeting the diverse S1-S2 region within voltage-sensing domains provides an avenue to develop subtype selective tools.

Inhibition of inwardly rectifying Kir2.x channels by the novel anti-cancer agent gambogic acid depends on both pore block and PIP2 interference

The caged xanthone gambogic acid (GA) is a novel anti-cancer agent which exhibits anti-proliferative, anti-inflammatory and cytotoxic effects in many types of cancer tissues. In a recent phase IIa study, GA exhibits a favourable safety profile. However, limited data are available concerning its interaction with cardiac ion channels. Heteromeric assembly of Kir2.x channels underlies the cardiac inwardly rectifying I_K1 current which is responsible for the stabilization of the diastolic resting membrane potential. Inhibition of the cardiac I_K1 current may lead to ventricular arrhythmia due to delayed afterdepolarizations. Compared to Kv2.1, hERG and Kir1.1, a slow, delayed inhibition of Kir2.1 channels by GA in a mammalian cell line was reported before but no data exist in literature concerning action of GA on homomeric Kir2.2 and Kir2.3 and heteromeric Kir2.x channels. Therefore, the aim of this study was to provide comparative data on the effect of GA on homomeric and heteromeric Kir2.x channels. Homomeric and heteromeric Kir2.x channels were heterologously expressed in Xenopus oocytes, and the two-microelectrode voltage-clamp technique was used to record Kir2.x currents. To investigate the mechanism of the channel inhibition by GA, alanine-mutated Kir2.x channels with modifications in the channels pore region or at phosphatidylinositol 4,5-bisphosphate (PIP₂)-binding sites were employed. GA caused a slow inhibition of homomeric and heteromeric Kir2.x channels at low micromolar concentrations (with IC₅₀ Kir2.1/2.2 < Kir2.2 < Kir2.2/2.3 < Kir2.3 < Kir2.1 < Kir2.1/2.3). The effect did not reach saturation within 60 min and was not reversible upon washout for 30 min. The inhibition showed no strong voltage dependence. We provide evidence for a combination of direct channel pore blockade and a PIP₂-dependent mechanism as a molecular basis for the observed effect. We conclude that Kir2.x channel inhibition by GA may be relevant in patients with pre-existing cardiac disorders such as chronic heart failure or certain rhythm disorders and recommend a close cardiac monitoring for those patients when treated with GA.

Utility of hERG Assays as Surrogate Markers of Delayed Cardiac Repolarization and QT Safety

HERG (human-ether-a-go-go-related gene) encodes for a cardiac potassium channel that plays a critical role in defining ventricular repolarization. Noncardiovascular drugs associated with a rare but potentially lethal ventricular arrhythmia (Torsades de Pointes) have been linked to delayed cardiac repolarization and block of hERG current. This brief overview will discuss the role of hERG current in cardiac electrophysiology, its involvement in drug-induced delayed repolarization, and approaches used to define drug effects on hERG current. In addition, examples of hERG blocking drugs acting differently (i.e., overt and covert hERG blockade due to multichannel block) together with the utility and limitations of hERG assays as tools to predict the risk of delayed repolarization and proarrhythmia are discussed.

Differential Inhibition of Inducible T Cell Cytokine Secretion by Potent Iron Chelators

Effector functions and proliferation of T helper (Th) cells are influenced by cytokines in the environment. Th1 cells respond to a synergistic effect of interleukin-12 (IL-12) and interleukin-18 (IL-18) to secrete interferon-gamma (IFN-γ). In contrast, Th2 cells respond to interleukin-4 (IL-4) to secrete IL-4, interleukin-13 (IL-13), interleukin-5 (IL-5), and interleukin-10 (IL-10). The authors were interested in identifying nonpeptide inhibitors of the Th1 response selective for the IL-12/IL-18-mediated secretion of IFN-γ while leaving the IL-4-mediated Th2 cytokine secretion relatively intact. The authors established a screening protocol using human peripheral blood mononuclear cells (PBMCs) and identified the hydrazino anthranilate compound 1 as a potent inhibitor of IL-12/IL-18-mediated IFN-γ secretion from CD3+ cells with an IC50 around 200 nM. The inhibitor was specific because it had virtually no effect on IL-4-mediated IL-13 release from the same population of cells. Further work established that compound 1 was a potent intracellular iron chelator that inhibited both IL-12/IL-18- and IL-4-mediated T cell proliferation. Iron chelation affects multiple cellular pathways in T cells. Thus, the IL-12/IL-18-mediated proliferation and IFN-γ secretion are very sensitive to intracellular iron concentration. However, the IL-4-mediated IL-13 secretion does not correlate with proliferation and is partially resistant to potent iron chelation

Evaluation of a High-Throughput Fluorescence Assay Method for hERG Potassium Channel Inhibition

The number of projects in drug development that fail in late phases because of cardiac side effects such as QT prolongation can impede drug discovery and development of projects. The molecular target responsible for QT prolongation by a wide range of pharmaceutical agents is the myocardial hERG potassium channel. It is therefore desirable to screen for compound interactions with the hERG channel at an early stage of drug development. Here, the authors report a cell-based fluorescence assay using membrane potential-sensitive fluorescent dyes and stably transfected hERG channels from CHO cells. The assay allows semiautomated screening of compounds for hERG activity on 384-well plates and is sufficiently rapid for testing a large number of compounds. The assay is robust as indicated by a Z′ factor larger than 0.6. The throughput is in the range of 10,000 data points per day, which is significantly higher than any other method presently available for hERG. The data obtained with the fluorescence assay were in qualitative agreement with those from patch-clamp electrophysiological analysis. There were no false-positive hits, and the rate of false-negative compounds is currently 12% but might be further reduced by testing compounds at higher concentration. Quantitative differences between fluorescence and electrophysiological methods may be due to the use- or voltage-dependentactivity of the antagonists.

Interactions between amiodarone and the hERG potassium channel pore determined with mutagenesis and in silico docking

The antiarrhythmic drug amiodarone delays cardiac repolarisation through inhibition of hERG-encoded potassium channels responsible for the rapid delayed rectifier potassium current (I_Kr). This study aimed to elucidate molecular determinants of amiodarone binding to the hERG channel. Whole-cell patch-clamp recordings were made at 37 °C of ionic current (I_hERG) carried by wild-type (WT) or mutant hERG channels expressed in HEK293 cells. Alanine mutagenesis and ligand docking were used to investigate the roles of pore cavity amino-acid residues in amiodarone binding. Amiodarone inhibited WT outward I_hERG tails with a half-maximal inhibitory concentration (IC₅₀) of ∼45 nM, whilst inward I_hERG tails in a high K⁺ external solution ([K⁺]_e) of 94 mM were blocked with an IC₅₀ of 117.8 nM. Amiodarone’s inhibitory action was contingent upon channel gating. Alanine-mutagenesis identified multiple residues directly or indirectly involved in amiodarone binding. The IC₅₀ for the S6 aromatic Y652A mutation was increased to ∼20-fold that of WT I_hERG, similar to the pore helical mutant S624A (∼22-fold WT control). The IC₅₀ for F656A mutant I_hERG was ∼17-fold its corresponding WT control. Computational docking using a MthK-based hERG model differentiated residues likely to interact directly with drug and those whose Ala mutation may affect drug block allosterically. The requirements for amiodarone block of aromatic residues F656 and Y652 within the hERG pore cavity are smaller than for other high affinity I_hERG inhibitors, with relative importance to amiodarone binding of the residues investigated being S624A ∼ Y652A > F656A > V659A > G648A > T623A.

Investigation of miscellaneous hERG inhibition in large diverse compound collection using automated patch-clamp assay

AIM

hERG potassium channels display miscellaneous interactions with diverse chemical scaffolds. In this study we assessed the hERG inhibition in a large compound library of diverse chemical entities and provided data for better understanding of the mechanisms underlying promiscuity of hERG inhibition.

METHODS

Approximately 300 000 compounds contained in Molecular Library Small Molecular Repository (MLSMR) library were tested. Compound profiling was conducted on hERG-CHO cells using the automated patch-clamp platform-IonWorks Quattro(™).

RESULTS

The compound library was tested at 1 and 10 μmol/L. IC50 values were predicted using a modified 4-parameter logistic model. Inhibitor hits were binned into three groups based on their potency: high (IC50<1 μmol/L), intermediate (1 μmol/L< IC50<10 μmol/L), and low (IC50>10 μmol/L) with hit rates of 1.64%, 9.17% and 16.63%, respectively. Six physiochemical properties of each compound were acquired and calculated using ACD software to evaluate the correlation between hERG inhibition and the properties: hERG inhibition was positively correlative to the physiochemical properties ALogP, molecular weight and RTB, and negatively correlative to TPSA.

CONCLUSION

Based on a large diverse compound collection, this study provides experimental evidence to understand the promiscuity of hERG inhibition. This study further demonstrates that hERG liability compounds tend to be more hydrophobic, high-molecular, flexible and polarizable.

In silico assessment of kinetics and state dependent binding properties of drugs causing acquired LQTS

The Kv11.1 or hERG potassium channel is responsible for one of the major repolarising currents (IKr) in cardiac myocytes. Drug binding to hERG can result in reduction in IKr, action potential prolongation, acquired long QT syndrome and fatal cardiac arrhythmias. The current guidelines for pre-clinical assessment of drugs in development is based on the measurement of the drug concentration that causes 50% current block, i.e., IC50. However, drugs with the same apparent IC50 may have very different kinetics of binding and unbinding, as well as different affinities for the open and inactivated states of Kv11.1. Therefore, IC50 measurements may not reflect the true risk of drug induced arrhythmias. Here we have used an in silico approach to test the hypothesis that drug binding kinetics and differences in state-dependent affinity will influence the extent of cardiac action potential prolongation independent of apparent IC50 values. We found, in general that drugs with faster overall kinetics and drugs with higher affinity for the open state relative to the inactivated state cause more action potential prolongation. These characteristics of drug-hERG interaction are likely to be more arrhythmogenic but cannot be predicted by IC50 measurement alone. Our results suggest that the pre-clinical assessment of Kv11.1-drug interactions should include descriptions of the kinetics and state dependence of drug binding. Further, incorporation of this information into sophisticated in silico models should be able to better predict arrhythmia risk and therefore more accurately assess safety of new drugs in development.

Human Ether-à-go-go Related Gene (hERG) Channel Blocking Aporphine Alkaloids from Lotus Leaves and Their Quantitative Analysis in Dietary Weight Loss Supplements

Blockage of the human ether-à-go-go related gene (hERG) channel can result in life-threatening ventricular tachyarrhythmia. In an in vitro screening of herbal materials for hERG blockers using an automated two-microelectrode voltage clamp assay on Xenopus oocytes, an alkaloid fraction of Nelumbo nucifera Gaertn. (lotus) leaves induced ∼50% of hERG current inhibition at 100 μg/mL. Chromatographic separation resulted in the isolation and identification of (-)-asimilobine, 1, nuciferine, 2, O-nornuciferine, 3, N-nornuciferine, 4, and liensinine, 5. In agreement with in silico predicted ligand-target interactions, 2, 3, and 4 revealed distinct in vitro hERG blockages measured in HEK293 cells with IC50 values of 2.89, 7.91, and 9.75 μM, respectively. Because lotus leaf dietary weight loss supplements are becoming increasingly popular, the identified hERG-blocking alkaloids were quantitated in five commercially available products. Results showed pronounced differences in the content of hERG-blocking alkaloids ranging up to 992 μg (2) in the daily recommended dose.

Computer-based prediction of the drug proarrhythmic effect: problems, issues, known and suspected challenges

It is likely that computer modelling and simulations will become an element of comprehensive cardiac safety testing. Their role would be primarily the integration and the interpretation of previously gathered data. There are still unanswered questions and issues which we list and describe below. They include sources of data used for the development of the models as well as data utilized as input information, which can come from the in vitro studies and the quantitative structure-activity relationship models. The pharmacokinetics of the drugs in question play a crucial role as their active concentration should be considered, yet the question remains where is the right place to assess it. The pharmacodynamic angle includes complications coming from multiple drugs (i.e. active metabolites) acting in parallel as well as the type of interaction with (potentially) multiple affected channels. Once established, the model and the methodology of its use should be further validated, optimistically against individual data reported at the clinical level as the physiological, anatomical, and genetic parameters play a crucial role in the drug-triggered arrhythmia induction. All the abovementioned issues should be at least considered and-hopefully-resolved, to properly utilize the mathematical models for a cardiac safety assessment.

The Therapeutic Potential of hERG1 K+ Channels for Treating Cancer and Cardiac Arrhythmias

hERG potassium channels present pharmacologists and medicinal chemists with a dilemma. On the one hand hERG is a major reason for drugs being withdrawn from the market because of drug induced long QT syndrome and the associated risk of inducing sudden cardiac death, and yet hERG blockers are still widely used in the clinic to treat cardiac arrhythmias. Moreover, in the last decade overwhelming evidence has been provided that hERG channels are aberrantly expressed in cancer cells and that they contribute to tumour cell proliferation, resistance to apoptosis, and neoangiogenesis. Here we provide an overview of the properties of hERG channels and their role in excitable cells of the heart and nervous system as well as in cancer. We consider the therapeutic potential of hERG, not only with regard to the negative impact due to drug induced long QT syndrome, but also its future potential as a treatment in the fight against cancer.

Enhanced slow inactivation of the human skeletal muscle sodium channel causing normokalemic periodic paralysis

Normokalemic periodic paralysis (normoPP) is a type of skeletal muscle function disorder which is characterized by paralysis attack with concomitant normal serum potassium level. We previously reported that R675Q mutation of human skeletal muscle voltage-gated sodium channel α subunit (SCN4A) may be the novel mutation which caused normoPP in Chinese families. However, it is still not clear how this mutation affects the SCN4A channel function. In this study, we used patch-clamp recording to study the function of wild type (WT) and R675Q mutant of SCN4A channels expressed in human embryonic kidney (HEK293) cells. We found that R675Q mutation did not affect the voltage dependence of sodium channel activation. The fast inactivation was also not significantly affected by R675Q mutation. However, R675Q mutation of SCN4A channels exhibited an 11.1 mV hyperpolarized shift in the voltage dependence of slow inactivation and significantly prolonged the recovery from prolonged inactivation state. Our results thus indicate that SCN4A was functionally affected by R675Q mutation, suggesting a possible reason for causing normoPP in Chinese patients.

UNC1062, a new and potent Mer inhibitor

Abnormal activation of Mer kinase has been implicated in the oncogenesis of many human cancers including acute lymphoblastic and myeloid leukemia, non-small cell lung cancer, and glioblastoma. We have discovered a new family of small molecule Mer inhibitors, pyrazolopyrimidine sulfonamides, that potently inhibit the kinase activity of Mer. Importantly, these compounds do not demonstrate significant hERG activity in the PatchXpress assay. Through structure-activity relationship studies, 35 (UNC1062) was identified as a potent (IC50 = 1.1 nM) and selective Mer inhibitor. When applied to live tumor cells, UNC1062 inhibited Mer phosphorylation and colony formation in soft agar. Given the potential of Mer as a therapeutic target, UNC1062 is a promising candidate for further drug development.

Role of the activation gate in determining the extracellular potassium dependency of block of HERG by trapped drugs

Drug induced long QT syndrome (diLQTS) results primarily from block of the cardiac potassium channel HERG (human-ether-a-go-go related gene). In some cases long QT syndrome can result in the lethal arrhythmia torsade de pointes, an arrhythmia characterized by a rapid heart rate and severely compromised cardiac output. Many patients requiring medication present with serum potassium abnormalities due to a variety of conditions including gastrointestinal dysfunction, renal and endocrine disorders, diuretic use, and aging. Extracellular potassium influences HERG channel inactivation and can alter block of HERG by some drugs. However, block of HERG by a number of drugs is not sensitive to extracellular potassium. In this study, we show that block of WT HERG by bepridil and terfenadine, two drugs previously shown to be trapped inside the HERG channel after the channel closes, is insensitive to extracellular potassium over the range of 0 mM to 20 mM. We also show that bepridil block of the HERG mutant D540K, a mutant channel that is unable to trap drugs, is dependent on extracellular potassium, correlates with the permeant ion, and is independent of HERG inactivation. These results suggest that the lack of extracellular potassium dependency of block of HERG by some drugs may in part be related to the ability of these drugs to be trapped inside the channel after the channel closes.

In vitro cardiovascular effects of dihydroartemisin-piperaquine combination compared with other antimalarials

The in vitro cardiac properties of dihydroartemisinin (DHA) plus piperaquine phosphate (PQP) were compared with those of other antimalarial compounds. Results with antimalarial drugs, chosen on the basis of their free therapeutic maximum concentration in plasma (C(max)), were expressed as the fold of that particular effect with respect to their C(max). The following tests were used at 37 °C: hERG (human ether-à-go-go-related gene) blockade and trafficking, rabbit heart ventricular preparations, and sodium and slow potassium ion current interference (I(Na) and I(Ks), respectively). Chloroquine, halofantrine, mefloquine, and lumefantrine were tested in the hERG studies, but only chloroquine, dofetilide, lumefantrine, and the combination of artemether-lumefantrine were used in the rabbit heart ventricular preparations, hERG trafficking studies, and I(Na) and I(Ks) analyses. A proper reference was used in each test. In hERG studies, the high 50% inhibitory concentration (IC(50)) of halofantrine, which was lower than its C(max), was confirmed. All the other compounds blocked hERG, with IC(50)s ranging from 3- to 30-fold their C(max)s. In hERG trafficking studies, the facilitative effects of chloroquine at about 30-fold its C(max) were confirmed and DHA blocked it at a concentration about 300-fold its C(max). In rabbit heart ventricular preparations, dofetilide, used as a positive control, revealed a high risk of torsades de pointes, whereas chloroquine showed a medium risk. Neither DHA-PQP nor artemether-lumefantrine displayed an in vitro signal for a significant proarrhythmic risk. Only chloroquine blocked the I(Na) ion current and did so at about 30-fold its C(max). No effect on I(Ks) was detected. In conclusion, despite significant hERG blockade, DHA-PQP and artemether-lumefantrine do not appear to induce potential torsadogenic effects in vitro, affect hERG trafficking, or block sodium and slow potassium ion currents.

Strong activation of ether-à-go-go-related gene 1 K+ channel isoforms by NS1643 in human embryonic kidney 293 and Chinese hamster ovary cells

Two different mechanisms leading to increased current have been described for the small-molecule human ether-à-go-go-related gene (herg) activator NS1643 [1,3-bis-(2-hydroxy-5-trifluoromethylphenyl)-urea]. On herg1a channels expressed in Xenopus laevis oocytes, it mainly acts via attenuation of inactivation and for rat (r) erg1b channels expressed in human embryonic kidney (HEK)-293 cells, it strongly shifts the activation curve to the left. We now investigated the NS1643 effects on erg1b channels in more detail and performed comparative experiments with rat and human erg1a in different expression systems. Significant differences were observed between expression systems, but not between the rat and human isoform. In HEK-293 or Chinese hamster ovary (CHO) cells, activation of rat erg1b channels occurred in a dose-dependent manner with a maximum current increase of 300% obtained with 10 μM NS1643. In contrast, the NS1643-induced strong leftward shift in the voltage dependence of activation further increased with higher drug concentration, needed more time to develop, and exhibited use dependence. Coexpression of KCNE1 or KCNE2 did not attenuate this NS1643 effect on erg1 channel activation and did thus not mimic the lower drug potency on this parameter observed in oocytes. NS1643 (10 μM) slowed erg1b channel deactivation and recovery from inactivation without significant changes in activation and inactivation kinetics. With the exception of accelerated activation, NS1643 affected erg1a channels similarly, but the effect was less pronounced than in erg1b or erg1a/1b channels. It is noteworthy that rerg1b and herg1a inactivation estimated from fully activated current voltage relationships were unaltered in the continued presence of 10 μM NS1643 in the mammalian expression systems, indicating qualitative differences from NS1643 effects in X. laevis oocytes.

Searching for new cardiovascular drugs: towards improved systems for drug screening?

INTRODUCTION

The pharmaceutical industry urgently needs new ways of profiling the safety and efficacy of new cardiovascular (CV) drugs and more effectively transitioning these compounds through the stages of CV drug screening. This article reviews new technologies and methodological innovations and assesses whether these frameworks offer improved solutions to the problems facing the contemporary CV drug development.

AREAS COVERED

The article comprises literature derived from a systematic search (from 2000 onwards) using the US patent office and ESP@CENET search engines as well as through multiple Boolean terms. The article focuses on patents relating to technologies and resources and categorises the patents according to their niche in the CV drug screening landscape.

EXPERT OPINION

The CV drug pipeline is stalling due to the inability of many contemporary drug screening frameworks to discriminate between safe, efficacious therapy and hazardous off-target effect. Given the current limitations of drug screening frameworks, there is little scope for expanding the CV drug portfolio with newer, safer drugs with improved mechanisms of action. New screening modalities are urgently needed. Searches reveal that there are few examples of truly new technologies and systems in the patent literature. This apparent failure to revamp facets of the CV drug screening process can only perpetuate the inability of current platforms to improve the CV drug pipeline. Consequently, with few exceptions, there is stagnation in pre-clinical assay design that limits the pharmaceutical industry's ability to search for new drugs in new and more effective ways.

Molecular determinants for activation of human ether-à-go-go-related gene 1 potassium channels by 3-nitro-n-(4-phenoxyphenyl) benzamide

Human ether-à-go-go-related gene 1 (hERG1) channels mediate repolarization of cardiac action potentials. Inherited long QT syndrome (LQTS) caused by loss-of-function mutations, or unintended blockade of hERG1 channels by many drugs, can lead to severe arrhythmia and sudden death. Drugs that activate hERG1 are a novel pharmacological approach to treat LQTS. 3-Nitro-n-(4-phenoxyphenyl) benzamide [ICA-105574 (ICA)] has been discovered to activate hERG1 by strong attenuation of pore-type inactivation. Here, we used scanning mutagenesis of hERG1 to identify the molecular determinants of ICA action. Three mutations abolished the activator effects of 30 μM ICA, including L622C in the pore helix, F557L in the S5 segment, and Y652A in the S6 segment. One mutation in S6 (A653M) switched the activity of ICA from an activator to an inhibitor, revealing its partial agonist activity. This was confirmed by showing that the noninactivating mutant hERG1 channel (G628C/S631C) was inhibited by ICA and that the addition of the F557L mutation rendered the channel drug-insensitive. Simulated molecular docking of ICA to homology models of hERG1 corroborated the scanning mutagenesis findings. Together, our findings indicate that ICA is a mixed agonist of hERG1 channels. Activation or inhibition of currents is mediated by the same or overlapping binding site located in the pore module between two adjacent subunits of the homotetrameric channel.

The open-access dataset for insilico cardiotoxicity prediction system

Drug cardiotoxicity is one of the main reasons of fatal drug related problem events and the subsequent withdrawals. Therefore, its early assessment is a crucial element of the drug development process. For the drug driven hERG inhibition assessment, which is assumed to be the main reason for toxicity, in vitro techniques are used. Gold standards are based on the Patch Clamp method with the use of various cell models but due to its low throughput, insilico models have become more appreciated. To develop a reliable empirical QSAR model, wide dataset containing a variety of cases has to be available. In this article, a freely available for download, set of data is described. It is based on literature peer-reviewed reports and contains hERG inhibition information expressed as IC50 value for 263 molecules described in 642 records. All studies were done with the use of three cell models (XO, CHO, HEK) and other elements describe the electrophysiological settings of the in vitro study. The above mentioned set was used for the successful development of the predictive models.

Ion channels in key marine invertebrates; their diversity and potential for applications in biotechnology

Of the intra-membrane proteins, the class that comprises voltage and ligand-gated ion channels represents the major substrate whereby signals pass between and within cells in all organisms. It has been presumed that vertebrate and particularly mammalian ion channels represent the apex of evolutionary complexity and diversity and much effort has been focused on understanding their function. However, the recent availability of cheap high throughput genome sequencing has massively broadened and deepened the quality of information across phylogeny and is radically changing this view. Here we review current knowledge on such channels in key marine invertebrates where physiological evidence is backed up by molecular sequences and expression/functional studies. As marine invertebrates represent a much greater range of phyla than terrestrial vertebrates and invertebrates together, we argue that these animals represent a highly divergent, though relatively underused source of channel novelty. As ion channels are exquisitely selective sensors for voltage and ligands, their potential and actual applications in biotechnology are manifold.

Role of voltage-gated sodium, potassium and calcium channels in the development of cocaine-associated cardiac arrhythmias

Cocaine is a highly active stimulant that alters dopamine metabolism in the central nervous system resulting in a feeling of euphoria that with time can lead to addictive behaviours. Cocaine has numerous deleterious effects in humans including seizures, vasoconstriction, ischaemia, increased heart rate and blood pressure, cardiac arrhythmias and sudden death. The cardiotoxic effects of cocaine are indirectly mediated by an increase in sympathomimetic stimulation to the heart and coronary vasculature and by a direct effect on the ion channels responsible for maintaining the electrical excitability of the heart. The direct and indirect effects of cocaine work in tandem to disrupt the co-ordinated electrical activity of the heart and have been associated with life-threatening cardiac arrhythmias. This review focuses on the direct effects of cocaine on cardiac ion channels, with particular focus on sodium, potassium and calcium channels, and on the contributions of these channels to cocaine-induced arrhythmias. Companion articles in this edition of the journal examine the epidemiology of cocaine use (Wood & Dargan) and the treatment of cocaine-associated arrhythmias (Hoffmann).

The use of Xenopus oocytes in drug screening

INTRODUCTION

The physiological roles of ion channels are receiving increased interest in both basic research and drug discovery, and a demand for pharmacological approaches that can characterize or screen ion channels and their ligands with higher throughput has emerged. Traditionally, screening of compound libraries at ion channel targets has been performed using assays such as binding assays, fluorescence-based assays and flux assays that allow high-throughput, but sacrifice high data quality. The use of these assays with ion channel targets can also be problematic, emphasizing the usefulness of automated Xenopus oocyte electrophysiological assays in drug screening.

AREAS COVERED

This review summarizes the use of Xenopus oocytes in drug screening, presents the advantages and disadvantages of the use of Xenopus oocytes as expression system, and addresses the options available for automated two-electrode voltage-clamp recordings from Xenopus oocytes.

EXPERT OPINION

Automated and manual Xenopus oocyte two-electrode voltage-clamp recordings are useful and important techniques in drug screening. Although they are not compatible with high-throughput experimentation, these techniques are excellent in combination or as alternatives to fluorescence-based assays for hit validation, screening of focused compound libraries and safety screening on ion channels with their high flexibility for the choice of molecular targets, quality of data and reproducibility.

Development of a selective small-molecule inhibitor of Kir1.1, the renal outer medullary potassium channel

The renal outer medullary potassium (K+) channel, ROMK (Kir1.1), is a putative drug target for a novel class of loop diuretic that would lower blood volume and pressure without causing hypokalemia. However, the lack of selective ROMK inhibitors has hindered efforts to assess its therapeutic potential. In a high-throughput screen for small-molecule modulators of ROMK, we previously identified a potent and moderately selective ROMK antagonist, 7,13-bis(4-nitrobenzyl)-1,4,10-trioxa-7,13-diazacyclopentadecane (VU590), that also inhibits Kir7.1. Because ROMK and Kir7.1 are coexpressed in the nephron, VU590 is not a good probe of ROMK function in the kidney. Here we describe the development of the structurally related inhibitor 2,2'-oxybis(methylene)bis(5-nitro-1H-benzo[d]imidazole) (VU591), which is as potent as VU590 but is selective for ROMK over Kir7.1 and more than 65 other potential off-targets. VU591 seems to block the intracellular pore of the channel. The development of VU591 may enable studies to explore the viability of ROMK as a diuretic target.

Cardiac slices as a predictive tool for arrhythmogenic potential of drugs and chemicals

IMPORTANCE OF THE FIELD

cardiac arrhythmia represents one of the primary safety pharmacological concerns in drug development. The most prominent example is drug induced ventricular tachycardia of the Torsade des Pointes type. The mechanism how this type of arrhythmia develops is a complex multi-cellular phenomenon. It can only be insufficiently reflected by cellular or molecular assays. However, organ models - such as Langendorff hearts - or in vivo experiments are expensive and time consuming and not suitable for assays requiring an increased throughput.

AREAS COVERED IN THIS REVIEW

here, we describe and review an assay bridging the gap between cardiomyocyte based assays and organ based systems - cardiac slices. This assay is reviewed in direct comparison with established safety pharmacological assays.

WHAT THE READER WILL GAIN

while slices have played an important role in brain research for > 2 decades, cardiac slices are experiencing a renaissance due to the novel challenges in safety pharmacology just in the last few years. Cardiac slices can be cultured and recorded over several days. It is possible to access electrophysiological data with a high number of electrodes - up to 256 electrodes - embedded in the surface of a microelectrode array.

TAKE HOME MESSAGE

cardiac slices close the gap between cellular and organ based assays in cardiac safety pharmacology. The tissue properties of a functional cardiac syncytium are more accurately reflected by a slice rather than a single cell.

Exploring chemical substructures essential for HERG k(+) channel blockade by synthesis and biological evaluation of dofetilide analogues

In this study we followed a new approach to analyze molecular substructures required for hERG channel blockade. We designed and synthesized 40 analogues of dofetilide (1), a potent hERG potassium channel blocker, and established structure-activity relationships (SAR) for their interaction with this important cardiotoxicity-related off-target. Structural modifications to dofetilide were made by diversifying the substituents on the phenyl rings and the protonated nitrogen and by varying the carbon chain length. The analogues were evaluated in a radioligand binding assay and SAR data were derived with the aim to specify structural features that give rise to hERG toxicity.

Molecular mechanisms for drug interactions with hERG that cause long QT syndrome

The human ether-a-go-go-related gene (hERG) encodes the pore-forming alpha-subunit of a voltage-gated potassium (K(+)) channel. A variety of unrelated compounds reduce K(+ )current in the heart by blocking the pore or disrupting trafficking of the hERG channel to the membrane surface. This induces a syndrome known as long QT, which arises from abnormalities in action potential repolarisation and can degenerate into lethal cardiac arrhythmias. As a result, this undesirable side effect has severely hindered safe drug development. This review describes progress in understanding the molecular basis for drug binding to hERG, outlines the characteristics of hERG ligands and discusses experimental and in silico approaches for identifying compounds with QT liabilities. Recent developments should enable recognition of hERG-positive compounds at the early stages of their development.

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.