Probing the molecular basis of hERG drug block with unnatural amino acids

An insight into the mechanisms of action of selected bioactive compounds against epigenetic targets of prostate cancer: implications on histones modifications

Prostate cancer is a leading cause of morbidity and mortality among men globally. In this study, we employed an in silico approach to predict the possible mechanisms of action of selected novel compounds reported against prostate cancer epigenetic targets and their derivatives, exhausting through ADMET profiling, drug-likeness, and molecular docking analyses. The selected compounds: sulforaphane, silibinin, 3, 3'-diindolylmethane (DIM), and genistein largely conformed to ADMET and drug-likeness rules including Lipinski's. Docking studies revealed strong binding energy of sulforaphane with HDAC6 (- 4.2 kcal/ mol), DIM versus HDAC2 (- 5.2 kcal/mol), genistein versus HDAC6 (- 4.1 kcal/mol), and silibinin against HDAC1 (- 7.0 kcal/mol) coupled with improved binding affinities and biochemical stabilities after derivatization. Findings from this study may provide insight into the potential epigenetic reprogramming mechanisms of these compounds against prostate cancer and could pave the way toward more success in prostate cancer phytotherapy.

Genetic and Molecular Aspects of Drug-Induced QT Interval Prolongation

Long QT syndromes can be either acquired or congenital. Drugs are one of the many etiologies that may induce acquired long QT syndrome. In fact, many drugs frequently used in the clinical setting are a known risk factor for a prolonged QT interval, thus increasing the chances of developing torsade de pointes. The molecular mechanisms involved in the prolongation of the QT interval are common to most medications. However, there is considerable inter-individual variability in drug response, thus making the application of personalized medicine a relevant aspect in long QT syndrome, in order to evaluate the risk of every individual from a pharmacogenetic standpoint.

Targeting KNa1.1 channels in KCNT1-associated epilepsy

Gain-of-function (GOF) pathogenic variants of KCNT1, the gene encoding the largest known potassium channel subunit, K_Na1.1, are associated with developmental and epileptic encephalopathies accompanied by severe psychomotor and intellectual disabilities. Blocking hyperexcitable K_Na1.1 channels with quinidine, a class I antiarrhythmic drug, has shown variable success in patients in part because of dose-limiting off-target effects, poor blood-brain barrier (BBB) penetration, and low potency. In recent years, high-resolution cryogenic electron microscopy (cryo-EM) structures of the chicken K_Na1.1 channel in different activation states have been determined, and animal models of the diseases have been generated. Alongside increasing information about the functional effects of GOF pathogenic variants on K_Na1.1 channel behaviour and how they lead to hyperexcitability, these tools will facilitate the development of more effective treatment strategies. We review the range of KCNT1 variants and their functional effects, the challenges posed by current treatment strategies, and recent advances in finding more potent and selective therapeutic interventions for KCNT1-related epilepsies.

Flavonoids and hERG channels: Friends or foes?

The cardiac action potential is regulated by several ion channels. Drugs capable to block these channels, in particular the human ether-à-go-go-related gene (hERG) channel, also known as K_V11.1 channel, may lead to a potentially lethal ventricular tachyarrhythmia called "Torsades de Pointes". Thus, evaluation of the hERG channel off-target activity of novel chemical entities is nowadays required to safeguard patients as well as to avoid attrition in drug development. Flavonoids, a large class of natural compounds abundantly present in food, beverages, herbal medicines, and dietary food supplements, generally escape this assessment, though consumed in consistent amounts. Continuously growing evidence indicates that these compounds may interact with the hERG channel and block it. The present review, by examining numerous studies, summarizes the state-of-the-art in this field, describing the most significant examples of direct and indirect inhibition of the hERG channel current operated by flavonoids. A description of the molecular interactions between a few of these natural molecules and the Rattus norvegicus channel protein, achieved by an in silico approach, is also presented.

Allosteric Coupling Between Drug Binding and the Aromatic Cassette in the Pore Domain of the hERG1 Channel: Implications for a State-Dependent Blockade

Human-ether-a-go-go-related channel (hERG1) is the pore-forming domain of the delayed rectifier K⁺ channel in the heart which underlies the I_Kr current. The channel has been extensively studied due to its propensity to bind chemically diverse group of drugs. The subsequent hERG1 block can lead to a prolongation of the QT interval potentially leading to an abnormal cardiac electrical activity. The recently solved cryo-EM structure featured a striking non-swapped topology of the Voltage-Sensor Domain (VSD) which is packed against the pore-domain as well as a small and hydrophobic intra-cavity space. The small size and hydrophobicity of the cavity was unexpected and challenges the already-established hypothesis of drugs binding to the wide cavity. Recently, we showed that an amphipathic drug, ivabradine, may favorably bind the channel from the lipid-facing surface and we discovered a mutant (M651T) on the lipid facing domain between the VSD and the PD which inhibited the blocking capacity of the drug. Using multi-microseconds Molecular Dynamics (MD) simulations of wild-type and M651T mutant hERG1, we suggested the block of the channel through the lipid mediated pathway, the opening of which is facilitated by the flexible phenylalanine ring (F656). In this study, we characterize the dynamic interaction of the methionine-aromatic cassette in the S5-S6 helices by combining data from electrophysiological experiments with MD simulations and molecular docking to elucidate the complex allosteric coupling between drug binding to lipid-facing and intra-cavity sites and aromatic cassette dynamics. We investigated two well-established hERG1 blockers (ivabradine and dofetilide) for M651 sensitivity through electrophysiology and mutagenesis techniques. Our electrophysiology data reveal insensitivity of dofetilide to the mutations at site M651 on the lipid facing side of the channel, mirroring our results obtained from docking experiments. Moreover, we show that the dofetilide-induced block of hERG1 occurs through the intracellular space, whereas little to no block of ivabradine is observed during the intracellular application of the drug. The dynamic conformational rearrangement of the F656 appears to regulate the translocation of ivabradine into the central cavity. M651T mutation appears to disrupt this entry pathway by altering the molecular conformation of F656.

Structure-Based Identification and Characterization of Inhibitors of the Epilepsy-Associated KNa1.1 (KCNT1) Potassium Channel

Drug-resistant epileptic encephalopathies of infancy have been associated with KCNT1 gain-of-function mutations, which increase the activity of KNa1.1 sodium-activated potassium channels. Pharmacological inhibition of hyperactive KNa1.1 channels by quinidine has been proposed as a stratified treatment, but mostly this has not been successful, being linked to the low potency and lack of specificity of the drug. Here we describe the use of a previously determined cryo-electron microscopy-derived KNa1.1 structure and mutational analysis to identify how quinidine binds to the channel pore and, using computational methods, screened for compounds predicated to bind to this site. We describe six compounds that inhibited KNa1.1 channels with low- and sub-micromolar potencies, likely also through binding in the intracellular pore vestibule. In hERG inhibition and cytotoxicity assays, two compounds were ineffective. These may provide starting points for the development of new pharmacophores and could become tool compounds to study this channel further.

An Update on the Structure of hERG

The human voltage-sensitive K⁺ channel hERG plays a fundamental role in cardiac action potential repolarization, effectively controlling the QT interval of the electrocardiogram. Inherited loss- or gain-of-function mutations in hERG can result in dangerous “long” (LQTS) or “short” QT syndromes (SQTS), respectively, and the anomalous susceptibility of hERG to block by a diverse range of drugs underlies an acquired LQTS. A recent open channel cryo-EM structure of hERG should greatly advance understanding of the molecular basis of hERG channelopathies and drug-induced LQTS. Here we describe an update of recent research that addresses the nature of the particular gated state of hERG captured in the new structure, and the insight afforded by the structure into the molecular basis for high affinity drug block of hERG, the binding of hERG activators and the molecular basis of hERG's peculiar gating properties. Interpretation of the pharmacology of natural SQTS mutants in the context of the structure is a promising approach to understanding the molecular basis of hERG inactivation, and the structure suggests how voltage-dependent changes in the membrane domain may be transmitted to an extracellular “turret” to effect inactivation through aromatic side chain motifs that are conserved throughout the KCNH family of channels.

Histidine at position 462 determines the low quinine sensitivity of ether-à-go-go channel superfamily member Kv 12.1

BACKGROUND AND PURPOSE

The ether-à-go-go (Eag) Kv superfamily comprises closely related Kv 10, Kv 11, and Kv 12 subunits. Kv 11.1 (termed hERG in humans) gained much attention, as drug-induced inhibition of these channels is a frequent cause of sudden death in humans. The exclusive drug sensitivity of Kv 11.1 can be explained by central drug-binding pockets that are absent in most other channels. Currently, it is unknown whether Kv 12 channels are equipped with an analogous drug-binding pocket and whether drug-binding properties are conserved in all Eag superfamily members.

EXPERIMENTAL APPROACH

We analysed sensitivity of recombinant Kv 12.1 channels to quinine, a substituted quinoline that blocks Kv 10.1 and Kv 11.1 at low micromolar concentrations.

KEY RESULTS

Quinine inhibited Kv 12.1, but its affinity was 10-fold lower than for Kv 11.1. Contrary to Kv 11.1, quinine inhibited Kv 12.1 in a largely voltage-independent manner and induced channel opening at more depolarised potentials. Low sensitivity of Kv 12.1 and characteristics of quinine-dependent inhibition were determined by histidine 462, as site-directed mutagenesis of this residue into the homologous tyrosine of Kv 11.1 conferred Kv 11.1-like quinine block to Kv 12.1(H462Y). Molecular modelling demonstrated that the low affinity of Kv 12.1 was determined by only weak interactions of residues in the central cavity with quinine. In contrast, more favourable interactions can explain the higher quinine sensitivity of Kv 12.1(H462Y) and Kv 11.1 channels.

CONCLUSIONS AND IMPLICATIONS

The quinoline-binding "motif" is not conserved within the Eag superfamily, although the overall architecture of these channels is apparently similar. Our findings highlight functional and pharmacological diversity in this group of evolutionary-conserved channels.

Predicting critical drug concentrations and torsadogenic risk using a multiscale exposure-response simulator

Torsades de pointes is a serious side effect of many drugs that can trigger sudden cardiac death, even in patients with structurally normal hearts. Torsadogenic risk has traditionally been correlated with the blockage of a specific potassium channel and a prolonged recovery period in the electrocardiogram. However, the precise mechanisms by which single channel block translates into heart rhythm disorders remain incompletely understood. Here we establish a multiscale exposure-response simulator that converts block-concentration characteristics from single cell recordings into three-dimensional excitation profiles and electrocardiograms to rapidly assess torsadogenic risk. For the drug dofetilide, we characterize the QT interval and heart rate at different drug concentrations and identify the critical concentration at the onset of torsades de pointes: For dofetilide concentrations of 2x, 3x, and 4x, as multiples of the free plasma concentration Cmax = 2.1 nM, the QT interval increased by +62.0%, +71.2%, and +82.3% compared to baseline, and the heart rate changed by -21.7%, -23.3%, and +88.3%. The last number indicates that, at the critical concentration of 4x, the heart spontaneously developed an episode of a torsades-like arrhythmia. Strikingly, this critical drug concentration is higher than the concentration estimated from early afterdepolarizations in single cells and lower than in one-dimensional cable models. Our results highlight the importance of whole heart modeling and explain, at least in part, why current regulatory paradigms often fail to accurately quantify the pro-arrhythmic potential of a drug. Our exposure-response simulator could provide a more mechanistic assessment of pro-arrhythmic risk and help establish science-based guidelines to reduce rhythm disorders, design safer drugs, and accelerate drug development.

Drug-likeness of linear pentamidine analogues and their impact on the hERG K+channel – correlation with structural features

The pentamidines with S atoms or sulfanilide groups in the linker have favorable drug-likeness parameters and low toxicity.

Investigating the state dependence of drug binding in hERG channels using a trapped-open channel phenotype

The hERG channel is a key player in repolarization of the cardiac action potential. Pharmacological blockade of hERG channels depletes the cardiac repolarization reserve, increasing the risk of cardiac arrhythmias. The promiscuous nature of drug interactions with hERG presents a therapeutic challenge for drug design and development. Despite considerable effort, the mechanisms of drug binding remain incompletely understood. One proposed mechanism is that high-affinity drug binding preferentially occurs when channels are in the inactivated state. However, this has been difficult to test, since inactivation is rapid in hERG and access to the drug binding site is limited by slower opening of the activation gate. Here, we have directly assessed the role of inactivation in cisparide and terfenadine drug binding in mutant (I663P) hERG channels where the activation gate is trapped-open. We firstly demonstrate the utility of this approach by showing that inactivation, ion selectivity and high affinity drug binding are preserved in I663P mutant channels. We then assess the role of inactivation by applying cisapride and terfenadine at different membrane voltages, which induce varying degrees of inactivation. We show that the extent of block does not correlate with the extent of inactivation. These data suggest that inactivation is not a major determinant of cisapride or terfenadine binding in hERG channels.

Structural implications of hERG K+ channel block by a high-affinity minimally structured blocker

Cardiac potassium channels encoded by human ether-à-go-go–related gene (hERG) are major targets for structurally diverse drugs associated with acquired long QT syndrome. This study characterized hERG channel inhibition by a minimally structured high-affinity hERG inhibitor, Cavalli-2, composed of three phenyl groups linked by polymethylene spacers around a central amino group, chosen to probe the spatial arrangement of side chain groups in the high-affinity drug-binding site of the hERG pore. hERG current (I_hERG) recorded at physiological temperature from HEK293 cells was inhibited with an IC₅₀ of 35.6 nm with time and voltage dependence characteristic of blockade contingent upon channel gating. Potency of Cavalli-2 action was markedly reduced for attenuated inactivation mutants located near (S620T; 54-fold) and remote from (N588K; 15-fold) the channel pore. The S6 Y652A and F656A mutations decreased inhibitory potency 17- and 75-fold, respectively, whereas T623A and S624A at the base of the selectivity filter also decreased potency (16- and 7-fold, respectively). The S5 helix F557L mutation decreased potency 10-fold, and both F557L and Y652A mutations eliminated voltage dependence of inhibition. Computational docking using the recent cryo-EM structure of an open channel hERG construct could only partially recapitulate experimental data, and the high dependence of Cavalli-2 block on Phe-656 is not readily explainable in that structure. A small clockwise rotation of the inner (S6) helix of the hERG pore from its configuration in the cryo-EM structure may be required to optimize Phe-656 side chain orientations compatible with high-affinity block.