Molecular basis of hERG potassium channel blockade by the class Ic antiarrhythmic flecainide

Engineered cocultures of iPSC-derived atrial cardiomyocytes and atrial fibroblasts for modeling atrial fibrillation

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia treatable with antiarrhythmic drugs; however, patient responses remain highly variable. Human induced pluripotent stem cell-derived atrial cardiomyocytes (iPSC-aCMs) are useful for discovering precision therapeutics, but current platforms yield phenotypically immature cells and are not easily scalable for high-throughput screening. Here, primary adult atrial, but not ventricular, fibroblasts induced greater functional iPSC-aCM maturation, partly through connexin-40 and ephrin-B1 signaling. We developed a protein patterning process within multiwell plates to engineer patterned iPSC-aCM and atrial fibroblast coculture (PC) that significantly enhanced iPSC-aCM structural, electrical, contractile, and metabolic maturation for 6+ weeks compared to conventional mono-/coculture. PC displayed greater sensitivity for detecting drug efficacy than monoculture and enabled the modeling and pharmacological or gene editing treatment of an AF-like electrophysiological phenotype due to a mutated sodium channel. Overall, PC is useful for elucidating cell signaling in the atria, drug screening, and modeling AF.

Structural modeling of hERG channel-drug interactions using Rosetta

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

Utilizing Heteroatom Types and Numbers from Extensive Ligand Libraries to Develop Novel hERG Blocker QSAR Models Using Machine Learning-Based Classifiers

The human ether-à-go-go-related gene (hERG) channel plays a crucial role in membrane repolarization. Any disruptions in its function can lead to severe cardiovascular disorders such as long QT syndrome (LQTS), which increases the risk of serious cardiovascular problems such as tachyarrhythmia and sudden cardiac death. Drug-induced LQTS is a significant concern and has resulted in drug withdrawals from the market in the past. The main objective of this study is to pinpoint crucial heteroatoms present in ligands that initiate interactions leading to the effective blocking of the hERG channel. To achieve this aim, ligand-based quantitative structure-activity relationships (QSAR) models were constructed using extensive ligand libraries, considering the heteroatom types and numbers, and their associated hERG channel blockage pIC50 values. Machine learning-assisted QSAR models were developed to analyze the key structural components influencing compound activity. Among the various methods, the KPLS method proved to be the most efficient, allowing the construction of models based on eight distinct fingerprints. The study delved into investigating the influence of heteroatoms on the activity of hERG blockers, revealing their significant role. Furthermore, by quantifying the effect of heteroatom types and numbers on ligand activity at the hERG channel, six compound pairs were selected for molecular docking. Subsequent molecular dynamics simulations and per residue MM/GBSA calculations were performed to comprehensively analyze the interactions of the selected pair compounds.

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.

Store-Operated Ca2+ Entry as a Putative Target of Flecainide for the Treatment of Arrhythmogenic Cardiomyopathy

Arrhythmogenic cardiomyopathy (ACM) is a genetic disorder that may lead patients to sudden cell death through the occurrence of ventricular arrhythmias. ACM is characterised by the progressive substitution of cardiomyocytes with fibrofatty scar tissue that predisposes the heart to life-threatening arrhythmic events. Cardiac mesenchymal stromal cells (C-MSCs) contribute to the ACM by differentiating into fibroblasts and adipocytes, thereby supporting aberrant remodelling of the cardiac structure. Flecainide is an Ic antiarrhythmic drug that can be administered in combination with β-adrenergic blockers to treat ACM due to its ability to target both Nav1.5 and type 2 ryanodine receptors (RyR2). However, a recent study showed that flecainide may also prevent fibro-adipogenic differentiation by inhibiting store-operated Ca2+ entry (SOCE) and thereby suppressing spontaneous Ca2+ oscillations in C-MSCs isolated from human ACM patients (ACM C-hMSCs). Herein, we briefly survey ACM pathogenesis and therapies and then recapitulate the main molecular mechanisms targeted by flecainide to mitigate arrhythmic events, including Nav1.5 and RyR2. Subsequently, we describe the role of spontaneous Ca2+ oscillations in determining MSC fate. Next, we discuss recent work showing that spontaneous Ca2+ oscillations in ACM C-hMSCs are accelerated to stimulate their fibro-adipogenic differentiation. Finally, we describe the evidence that flecainide suppresses spontaneous Ca2+ oscillations and fibro-adipogenic differentiation in ACM C-hMSCs by inhibiting constitutive SOCE.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

Intracellular Binding of Terfenadine Competes with Its Access to Pancreatic ß-cell ATP-Sensitive K+ Channels and Human ether-à-go-go-Related Gene Channels

Most blockers of both hERG (human ether-à-go-go-related gene) channels and pancreatic ß-cell ATP-sensitive K⁺ (K_ATP) channels access their binding sites from the cytoplasmic side of the plasma membrane. It is unknown whether binding to intracellular components competes with binding of these substances to K⁺ channels. The whole-cell configuration of the patch-clamp technique, a laser-scanning confocal microscope, and fluorescence correlation spectroscopy (FCS) were used to study hERG channels expressed in HEK (human embryonic kidney) 293 cells and K_ATP channels from the clonal insulinoma cell line RINm5F. When applied via the pipette solution in the whole-cell configuration, terfenadine blocked both hERG and K_ATP currents with much lower potency than after application via the bath solution, which was not due to P-glycoprotein-mediated efflux of terfenadine. Such a difference was not observed with dofetilide and tolbutamide. 37-68% of hERG/EGFP (enhanced green-fluorescent protein) fusion proteins expressed in HEK 293 cells were slowly diffusible as determined by laser-scanning microscopy in the whole-cell configuration and by FCS in intact cells. Bath application of a green-fluorescent sulphonylurea derivative (Bodipy-glibenclamide) induced a diffuse fluorescence in the cytosol of RINm5F cells under whole-cell patch-clamp conditions. These observations demonstrate the presence of intracellular binding sites for hERG and K_ATP channel blockers not dialyzable by the patch-pipette solution. Intracellular binding of terfenadine was not influenced by a mutated hERG (Y652A) channel. In conclusion, substances with high lipophilicity are not freely diffusible inside the cell but steep concentration gradients might exist within the cell and in the sub-membrane space.

New Insights into Ion Channels: Predicting hERG-Drug Interactions

Drug-induced long QT syndrome can be a very dangerous side effect of existing and developmental drugs. In this work, a model proposed two decades ago addressing the ion specificity of potassium channels is extended to the human ether-à-gogo gene (hERG). hERG encodes the protein that assembles into the potassium channel responsible for the delayed rectifier current in ventricular cardiac myocytes that is often targeted by drugs associated with QT prolongation. The predictive value of this model can guide a rational drug design decision early in the drug development process and enhance NCE (New Chemical Entity) retention. Small molecule drugs containing a nitrogen that can be protonated to afford a formal +1 charge can interact with hERG to prevent the repolarization of outward rectifier currents. Low-level ab initio calculations are employed to generate electronic features of the drug molecules that are known to interact with hERG. These calculations were employed to generate structure–activity relationships (SAR) that predict whether a small molecule drug containing a protonated nitrogen has the potential to interact with and inhibit the activity of the hERG potassium channels of the heart. The model of the mechanism underlying the ion specificity of potassium channels offers predictive value toward optimizing drug design and, therefore, minimizes the effort and expense invested in compounds with the potential for life-threatening inhibitory activity of the hERG potassium channel.

How does flecainide impact RyR2 channel function?

Flecainide, a cardiac class 1C blocker of the surface membrane sodium channel (NaV1.5), has also been reported to reduce cardiac ryanodine receptor (RyR2)-mediated sarcoplasmic reticulum (SR) Ca2+ release. It has been introduced as a clinical antiarrhythmic agent for catecholaminergic polymorphic ventricular tachycardia (CPVT), a condition most commonly associated with gain-of-function RyR2 mutations. Current debate concerns both cellular mechanisms of its antiarrhythmic action and molecular mechanisms of its RyR2 actions. At the cellular level, it targets NaV1.5, RyR2, Na+/Ca2+ exchange (NCX), and additional proteins involved in excitation-contraction (EC) coupling and potentially contribute to the CPVT phenotype. This Viewpoint primarily addresses the various direct molecular actions of flecainide on isolated RyR2 channels in artificial lipid bilayers. Such studies demonstrate different, multifarious, flecainide binding sites on RyR2, with voltage-dependent binding in the channel pore or voltage-independent binding at distant peripheral sites. In contrast to its single NaV1.5 pore binding site, flecainide may bind to at least four separate inhibitory sites on RyR2 and one activation site. None of these binding sites have been specifically located in the linear RyR2 sequence or high-resolution structure. Furthermore, it is not clear which of the inhibitory sites contribute to flecainide's reduction of spontaneous Ca2+ release in cellular studies. A confounding observation is that flecainide binding to voltage-dependent inhibition sites reduces cation fluxes in a direction opposite to physiological Ca2+ flow from SR lumen to cytosol. This may suggest that, rather than directly blocking Ca2+ efflux, flecainide can reduce Ca2+ efflux by blocking counter currents through the pore which otherwise limit SR membrane potential change during systolic Ca2+ efflux. In summary, the antiarrhythmic effects of flecainide in CPVT seem to involve multiple components of EC coupling and multiple actions on RyR2. Their clarification may identify novel specific drug targets and facilitate flecainide's clinical utilization in CPVT.

Molecular stratification of arrhythmogenic mechanisms in the Andersen Tawil Syndrome

Andersen Tawil Syndrome (ATS) is a rare inheritable disease associated with loss-of-function mutations in KCNJ2, the gene coding the strong inward rectifier potassium channel Kir2.1, which forms an essential membrane protein controlling cardiac excitability. ATS is usually marked by a triad of periodic paralysis, life-threatening cardiac arrhythmias and dysmorphic features, but its expression is variable and not all patients with a phenotype linked to ATS have a known genetic alteration. The mechanisms underlying this arrhythmogenic syndrome are poorly understood. Knowing such mechanisms would be essential to distinguish ATS from other channelopathies with overlapping phenotypes and to develop individualized therapies. For example, the recently suggested role of Kir2.1 as a countercurrent to sarcoplasmic calcium reuptake might explain the arrhythmogenic mechanisms of ATS and its overlap with catecholaminergic polymorphic ventricular tachycardia (CPVT). Here we summarize current knowledge on the mechanisms of arrhythmias leading to sudden cardiac death in ATS. We first provide an overview of the syndrome and its pathophysiology, from the patient´s bedside to the protein, and discuss the role of essential regulators and interactors that could play a role in cases of ATS. The review highlights novel ideas related to some post-translational channel interactions with partner proteins that might help define the molecular bases of the arrhythmia phenotype. We then propose a new all-embracing classification of the currently known ATS loss-of-function mutations according to their position in the Kir2.1 channel structure and their functional implications. We also discuss specific ATS pathogenic variants, their clinical manifestations and treatment stratification. The goal is to provide a deeper mechanistic understanding of the syndrome toward the development of novel targets and personalized treatment strategies.

Common Structural Pattern for Flecainide Binding in Atrial-Selective Kv1.5 and Nav1.5 Channels: A Computational Approach

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Its treatment includes antiarrhythmic drugs (AADs) to modulate the function of cardiac ion channels. However, AADs have been limited by proarrhythmic effects, non-cardiovascular toxicities as well as often modest antiarrhythmic efficacy. Theoretical models showed that a combined blockade of Na_v1.5 (and its current, I_Na) and K_v1.5 (and its current, I_Kur) ion channels yield a synergistic anti-arrhythmic effect without alterations in ventricles. We focused on K_v1.5 and Na_v1.5 to search for structural similarities in their binding site (BS) for flecainide (a common blocker and widely prescribed AAD) as a first step for prospective rational multi-target directed ligand (MTDL) design strategies. We present a computational workflow for a flecainide BS comparison in a flecainide-K_v1.5 docking model and a solved structure of the flecainide-Na_v1.5 complex. The workflow includes docking, molecular dynamics, BS characterization and pattern matching. We identified a common structural pattern in flecainide BS for these channels. The latter belongs to the central cavity and consists of a hydrophobic patch and a polar region, involving residues from the S6 helix and P-loop. Since the rational MTDL design for AF is still incipient, our findings could advance multi-target atrial-selective strategies for AF treatment.

Moving in the right direction: elucidating the mechanisms of interaction between flecainide and the cardiac ryanodine receptor

Flecainide is used to treat catecholaminergic polymorphic ventricular tachycardia (CPVT), an arrhythmia caused by disrupted cellular Ca²⁺ handling following β-adrenergic stimulation. The clinical efficacy of flecainide in this context involves complex effects on multiple ion channels that may be influenced by the disease state. A compelling narrative has been constructed around flecainide's nonselective block of sarcoplasmic reticulum (SR) lumen-to-cytoplasm Ca²⁺ release through intracellular calcium release channels (RyR2). However, ion fluxes across the SR membrane during heart contraction are bidirectional, and here, we review experimental evidence that flecainide's principal action on RyR2 involves the partial block of ion flow in the cytoplasm-to-lumen direction (i.e., flecainide inhibits RyR2-mediated SR 'countercurrent'). Experimental approaches that could advance new knowledge on the mechanism of RyR2 block by flecainide are proposed. Some impediments to progress in this area, that must be overcome to enable the development of superior drugs to treat CPVT, are also considered.

Electrocardiographic marker of the cardiac action potential triangulation induced by antiarrhythmic drugs in perfused guinea-pig heart

NEW FINDINGS

What is the central question of this study? Can triangular appearance of ventricular action potential, indicating proarrhythmic profile of antiarrhythmic agent, be approximated by specific changes on ECG? What is the main finding and its importance? The triangulation of the ventricular action potential seen when antiarrhythmic drugs induce a greater lengthening of the late repolarization compared to the initial repolarization in epicardium, is closely approximated by a greater prolongation of the T wave upslope relative to the interval between the J point and the start of the T wave (the JT_start interval) on ECG. These findings may improve the power of ECG assessments in predicting the drug-induced arrhythmia resulting from slowed phase 3 repolarization.

ABSTRACT

Antiarrhythmic drugs prescribed to treat atrial fibrillation can occasionally precipitate ventricular tachyarrhythmia through a prominent slowing of the phase 3 repolarization. The latter results in the triangular shape of ventricular action potential, indicating high arrhythmic risks. However, clinically, the utilility of triangulation assessments for predicting arrhythmia is limited owing to the invasive nature of the ventricular action potential recordings. This study examined whether the triangulation effect can be detected indirectly from ECG analysis. Epicardial monophasic action potentials and ECG were simultaneously recorded in perfused guinea-pig hearts. With antiarrhythmics (dofetilide, quinidine, procainamide and flecainide), a prolongation of the initial repolarization seen in the action potential recordings was closely approximated by lengthening of the interval bewteen the J point and the start of the T wave (the JT_start interval) on ECG, whereas a prolongation of the late repolarization was paralleled by widening of the T wave upslope. Dofetilide, quinidine and procainamide induced a prominent slowing of the phase 3 repolarization in epicardium, leading to triangulation of the action potential. These effects were accompanied by a greater prolongation of the T wave upslope compared to the JT_start interval. Flecainide elicited a proportional prolongation of the initial and the late ventricular repolarization, and therefore failed to induce triangulation, based on analysis of both epicardial action potential and ECG profiles. Collectively, these findings suggest that the ratio between the durations of the T wave upslope and the JT_start interval may represent ECG metric of the ventricular action potential triangulation induced by antiarrhythmic drugs. This article is protected by copyright. All rights reserved.

Metabolically driven maturation of human-induced-pluripotent-stem-cell-derived cardiac microtissues on microfluidic chips

The immature physiology of cardiomyocytes derived from human induced pluripotent stem cells (hiPSCs) limits their utility for drug screening and disease modelling. Here we show that suitable combinations of mechanical stimuli and metabolic cues can enhance the maturation of hiPSC-derived cardiomyocytes, and that the maturation-inducing cues have phenotype-dependent effects on the cells' action-potential morphology and calcium handling. By using microfluidic chips that enhanced the alignment and extracellular-matrix production of cardiac microtissues derived from genetically distinct sources of hiPSC-derived cardiomyocytes, we identified fatty-acid-enriched maturation media that improved the cells' mitochondrial structure and calcium handling, and observed divergent cell-source-dependent effects on action-potential duration (APD). Specifically, in the presence of maturation media, tissues with abnormally prolonged APDs exhibited shorter APDs, and tissues with aberrantly short APDs displayed prolonged APDs. Regardless of cell source, tissue maturation reduced variabilities in spontaneous beat rate and in APD, and led to converging cell phenotypes (with APDs within the 300-450 ms range characteristic of human left ventricular cardiomyocytes) that improved the modelling of the effects of pro-arrhythmic drugs on cardiac tissue.

The Antiarrhythmic Mechanisms of Flecainide in Catecholaminergic Polymorphic Ventricular Tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a severe yet rare inherited arrhythmia disorder. The cornerstone of CPVT medical therapy is the use of β-blockers; 30% of patients with CPVT do not respond well to optimal β-blocker treatment. Studies have shown that flecainide effectively prevents life-threatening arrhythmias in CPVT. Flecainide is a class IC antiarrhythmic drug blocking cardiac sodium channels. RyR2 inhibition is proposed as the principal mechanism of antiarrhythmic action of flecainide in CPVT, while it is highly debated. In this article, we review the current progress of this issue.

Conditional immortalization of human atrial myocytes for the generation of in vitro models of atrial fibrillation

The lack of a scalable and robust source of well-differentiated human atrial myocytes constrains the development of in vitro models of atrial fibrillation (AF). Here we show that fully functional atrial myocytes can be generated and expanded one-quadrillion-fold via a conditional cell-immortalization method relying on lentiviral vectors and the doxycycline-controlled expression of a recombinant viral oncogene in human foetal atrial myocytes, and that the immortalized cells can be used to generate in vitro models of AF. The method generated 15 monoclonal cell lines with molecular, cellular and electrophysiological properties resembling those of primary atrial myocytes. Multicellular in vitro models of AF generated using the immortalized atrial myocytes displayed fibrillatory activity (with activation frequencies of 6-8 Hz, consistent with the clinical manifestation of AF), which could be terminated by the administration of clinically approved antiarrhythmic drugs. The conditional cell-immortalization method could be used to generate functional cell lines from other human parenchymal cells, for the development of in vitro models of human disease.

Electrophysiological and Proarrhythmic Effects of Hydroxychloroquine Challenge in Guinea-Pig Hearts

Hydroxychloroquine (HCQ), clinically established in antimalarial and autoimmune therapy, recently raised cardiac arrhythmogenic concerns when used alone or with azithromycin (HCQ+AZM) in Covid-19. We report complementary, experimental, studies of its electrophysiological effects. In patch clamped HEK293 cells expressing human cardiac ion channels, HCQ inhibited IKr and IK1 at a therapeutic concentrations (IC50s: 10 ± 0.6 and 34 ± 5.0 μM). INa and ICaL showed higher IC50s; Ito and IKs were unaffected. AZM slightly inhibited INa, ICaL, IKs, and IKr, sparing IK1 and Ito. (HCQ+AZM) inhibited IKr and IK1 (IC50s: 7.7 ± 0.8 and 30.4 ± 3.0 μM), sparing INa, ICaL, and Ito. Molecular induced-fit docking modeling confirmed potential HCQ-hERG but weak AZM-hERG binding. Effects of μM-HCQ were studied in isolated perfused guinea-pig hearts by multielectrode, optical RH237 voltage, and Rhod-2 mapping. These revealed reversibly reduced left atrial and ventricular action potential (AP) conduction velocities increasing their heterogeneities, increased AP durations (APDs), and increased durations and dispersions of intracellular [Ca2+] transients, respectively. Hearts also became bradycardic with increased electrocardiographic PR and QRS durations. The (HCQ+AZM) combination accentuated these effects. Contrastingly, (HCQ+AZM) and not HCQ alone disrupted AP propagation, inducing alternans and torsadogenic-like episodes on voltage mapping during forced pacing. O'Hara-Rudy modeling showed that the observed IKr and IK1 effects explained the APD alterations and the consequently prolonged Ca2+ transients. The latter might then downregulate INa, reducing AP conduction velocity through recently reported INa downregulation by cytosolic [Ca2+] in a novel scheme for drug action. The findings may thus prompt future investigations of HCQ's cardiac safety under particular, chronic and acute, clinical situations.

Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study

Drug-induced blockade of the human ether-à-go-go-related gene (hERG) channel is today considered the main cause of cardiotoxicity in postmarketing surveillance. Hence, several ligand-based approaches were developed in the last years and are currently employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates. Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support vector machines were applied to integrate docking scores and protein–ligand interaction fingerprints. A total of 396 models were trained and validated based on: (i) high-quality experimental bioactivity information returned by 8337 curated compounds extracted from ChEMBL (version 25) and (ii) structural predictor data. Molecular docking simulations were performed using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and two published homology models selected for comparison. Interestingly, some classifiers return performances comparable to ligand-based models in terms of area under the ROC curve (AUC_MAX = 0.86 ± 0.01) and negative predictive values (NPV_MAX = 0.81 ± 0.01), thus putting forward the herein proposed computational workflow as a valuable tool for predicting hERG-related cardiotoxicity without the limitations of ligand-based models, typically affected by low interpretability and a limited applicability domain. From a methodological point of view, our study represents the first example of a successful integration of docking scores and protein–ligand interaction fingerprints (IFs) through a support vector machine (SVM) LASSO regularized strategy. Finally, the study highlights the importance of using hERG structural models accounting for ligand-induced fit effects and allowed us to select the best-performing protein conformation (made available in the Supporting Information, SI) to be employed for a reliable structure-based prediction of hERG-related cardiotoxicity.

Effects of Antiarrhythmic Drugs on hERG Gating in Human-Induced Pluripotent Stem Cell-Derived Cardiomyocytes From a Patient With Short QT Syndrome Type 1

Aims: The short QT syndrome type 1 (SQT1) is linked to hERG channel mutations (e.g., N588K). Drug effects on hERG channel gating kinetics in SQT1-cells have not been investigated. Methods: This study used hiPSC-CMs of a healthy donor and a SQT1-patient carrying the N588K mutation and patch clamp to examine the drug effects on hERG channel gating kinetics. Results: Ajmaline, amiodarone, ivabradine, flecainide, quinidine, mexiletine and ranolazine inhibited the hERG channel current (I_Kr) less strongly in hiPSC-CMs from the SQTS1-patient (SQT1-hiPSC-CMs) comparing with cells from the healthy donor (donor-hiPSC-CMs). Quinidine and mexiletine reduced, but ajmaline, amiodarone, ivabradine and ranolazine increased the time to peak of I_Kr similarly in SQT1-hiPSC-CMs and donor-hiPSC-CMs. Although regarding the shift of activation and inactivation curves, tested drugs showed differential effects in donor- and SQT1-hiPSC-CMs, quinidine, ajmaline, ivabradine and mexiletine but not amiodarone, flecainide and ranolazine reduced the window current in SQT1-hiPSC-CMs. Quinidine, ajmaline, ivabradine and mexiletine differentially changed the time constant of recovery from inactivation, but all of them increased the time constant of deactivation in SQT1-hiPSC-CMs. Conclusion: The window current-reducing and deactivation-slowing effects may be important for the antiarrhythmic effect of ajmaline, ivabradine, quinidine and mexiletine in SQT1-cells. This information may be helpful for selecting drugs for treating SQT1-patients with hERG channel mutation.

In silico toxicity evaluation of Salubrinal and its analogues

This paper reports on a comprehensive in silico toxicity assessment of Salubrinal and its analogues containing a cinnamic acid residue or quinoline ring using the online servers admetSAR, ADMETlab, ProTox, ADVERPred, Pred-hERG and Vienna LiverTox. Apart from rare exceptions, in all 55 studied structures, mild or practical absence of acute toxicity was predicted for rats (III or IV toxicity class). Cardiotoxic, hepatotoxic and immunotoxic effects were predicted for Salubrinal and its analogues. We constructed models of the main predicted anti-targets hERG, BSEP, MRP3, MRP4 and AhR using the principle of homologous modeling. Molecular docking studies were carried out with the obtained models. We carried out molecular docking for all targets using AutoDock Vina, implemented in the PyRx 0.8 software package. According to the results of molecular docking, the compounds analyzed are potential moderate or weak hERG blockers. Induction of cholestasis and, as a consequence, liver damage by these drugs, directly related to inhibition of BSEP, MRP3 and MRP4, most likely will not be observed. Interaction with AhR for the studied compounds is impossible for steric reasons and, as a consequence, toxic effects on the immune and other organ systems associated with the activation of the AhR signaling pathway are excluded.

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.

Use of Flecainide for the Treatment of Atrial Fibrillation

Atrial fibrillation (AF) is the most common sustained arrhythmia and is associated with substantial morbidity and impairment of quality of life. Restoration and maintenance of normal sinus rhythm is a desirable goal for many patients with AF; however, this strategy is limited by the relatively small number of antiarrhythmic drugs (AADs) available for AF rhythm control. Although it is recommended in current medical guidelines as first-line therapy for patients without structural heart disease, the use of flecainide has been curtailed since the completion of the Cardiac Arrhythmia Suppression Trial. In clinical trials and real-world use, flecainide has proven to be more effective than other AADs for the acute termination of recent onset AF. Flecainide is also moderately effective and, with the exception of amiodarone, equivalent to other AADs for the chronic suppression of paroxysmal and persistent AF. In patients without structural heart disease, flecainide has been demonstrated to be safe and well tolerated relative to other AADs. Despite this favorable profile, flecainide is underutilized, likely due to a perceived risk of ventricular proarrhythmia, a concern that has not been borne out in patients without underlying structural heart disease. Guidelines for administration and use of flecainide are summarized in this review.

Synthesis, Antitumor Activity and Molecular Docking Studies on Seven Novel Thiazacridine Derivatives

AIM AND OBJECTIVE

In the last decades, cancer has become a major problem in public health all around the globe. Chimeric chemical structures have been established as an important trend on medicinal chemistry in the last years. Thiazacridines are hybrid molecules composed of a thiazolidine and acridine nucleus, both pharmacophores that act on important biological targets for cancer. By the fact it is a serious disease, seven new 3-acridin-9-ylmethyl-thiazolidine-2,4-dione derivatives were synthesized, characterized, analyzed by computer simulation and tested in tumor cells. In order to find out if the compounds have therapeutic potential.

MATERIALS AND METHODS

Seven new 3-acridin-9-ylmethyl-thiazolidine-2,4-dione derivatives were synthesized through Michael addition and Knoevenagel condensation strategies. Characterization was performed by NMR and Infrared spectroscopy techniques. Regarding biological activity, thiazacridines were tested against solid and hematopoietic tumoral cell lines, namely Jurkat (acute T-cell leukemia); HL-60 (acute promyelocytic leukemia); DU 145 (prostate cancer); MOLT-4 (acute lymphoblastic leukemia); RAJI (Burkitt's lymphoma); K562 (chronic myelogenous leukemia) and normal cells PBMC (healthy volunteers). Molecular docking analysis was also performed in order to assess major targets of these new compounds. Cell cycle and clonogenic assay were also performed.

RESULTS

Compound LPSF/AA-62 (9f) exhibited the most potent anticancer activity against HL-60 (IC50 3,7±1,7 μM), MOLT-4 (IC50 5,7±1,1 μM), Jurkat (IC50 18,6 μM), Du-145 (IC50 20±5 μM) and Raji (IC50 52,3±9,2 μM). While the compound LPSF/AA-57 (9b) exhibited anticancer activity against the K562 cell line (IC50 51,8±7,8 μM). Derivative LPSF/AA-62 (9f) did not interfere in the cell cycle phases of the Molt-4 lineage. However, the LPSF/AA-62 (9f) derivative significantly reduced the formation of prostate cancer cell clones. The compound LPSF/AA-62 (9f) has shown strong anchorage stability with enzymes topoisomerases 1 and 2, in particular due the presence of chlorine favored hydrogen bonds with topoisomerase 1.

CONCLUSION

The 3-(acridin-9-ylmethyl)-5-((10-chloroanthracen-9-yl)methylene)thiazolidine-2,4-dione (LPSF/AA-62) presented the most promising results, showing anti-tumor activity in 5 of the 6 cell types tested, especially inhibiting the formation of colonies of prostate tumor cells (DU-145).

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.

Mitragynine, an euphoric compound inhibits hERG1a/1b channel current and upregulates the complexation of hERG1a-Hsp90 in HEK293-hERG1a/1b cells

Mitragyna speciosa Korth (M. speciosa) has been widely used as a recreational product, however, there are growing concerns on the abuse potentials and toxicity of the plant. Several poisoning and fatal cases involving kratom and mitragynine have been reported but the underlying causes remain unclear. The human ether-a-go-go-related gene 1 (hERG1) encodes the pore-forming subunit underlying cardiac rapidly delayed rectifier potassium current (I_Kr). Pharmacological blockade of the I_Kr can cause acquired long QT syndrome, leading to lethal cardiac arrhythmias. This study aims to elucidate the mechanisms of mitragynine-induced inhibition on hERG1a/1b current. Electrophysiology experiments were carried out using Port-a-Patch system. Quantitative RT-PCR, Western blot analysis, immunofluorescence and co-immunoprecipitation methods were used to determine the effects of mitragynine on hERG1a/1b expression and hERG1-cytosolic chaperones interaction. Mitragynine was found to inhibit the I_Kr current with an IC₅₀ value of 332.70 nM. It causes a significant reduction of the fully-glycosylated (fg) hERG1a protein expression but upregulates both core-glycosylated (cg) expression and hERG1a-Hsp90 complexes, suggesting possible impaired hERG1a trafficking. In conclusion, mitragynine inhibits hERG1a/1b current through direct channel blockade at lower concentration, but at higher concentration, it upregulates the complexation of hERG1a-Hsp90 which may be inhibitory towards channel trafficking.

Cardiac Na+/Ca2+ exchange stimulators among cardioprotective drugs

We previously reviewed our study of the pharmacological properties of cardiac Na+/Ca2+ exchange (NCX1) inhibitors among cardioprotective drugs, such as amiodarone, bepridil, dronedarone, cibenzoline, azimilide, aprindine, and benzyl-oxyphenyl derivatives (Watanabe et al. in J Pharmacol Sci 102:7-16, 2006). Since then we have continued our studies further and found that some cardioprotective drugs are NCX1 stimulators. Cardiac Na+/Ca2+ exchange current (INCX1) was stimulated by nicorandil (a hybrid ATP-sensitive K+ channel opener), pinacidil (a non-selective ATP-sensitive K+ channel opener), flecainide (an antiarrhythmic drug), and sodium nitroprusside (SNP) (an NO donor). Sildenafil (a phosphodiesterase-5 inhibitor) further increased the pinacidil-induced augmentation of INCX1. In paper, here I review the NCX stimulants that enhance NCX function among the cardioprotective agents we examined such as nicorandil, pinacidil, SNP, sildenafil and flecainide, in addition to atrial natriuretic (ANP) and dofetilide, which were reported by other investigators.

Computational Modeling of Electrophysiology and Pharmacotherapy of Atrial Fibrillation: Recent Advances and Future Challenges

The pathophysiology of atrial fibrillation (AF) is broad, with components related to the unique and diverse cellular electrophysiology of atrial myocytes, structural complexity, and heterogeneity of atrial tissue, and pronounced disease-associated remodeling of both cells and tissue. A major challenge for rational design of AF therapy, particularly pharmacotherapy, is integrating these multiscale characteristics to identify approaches that are both efficacious and independent of ventricular contraindications. Computational modeling has long been touted as a basis for achieving such integration in a rapid, economical, and scalable manner. However, computational pipelines for AF-specific drug screening are in their infancy, and while the field is progressing quite rapidly, major challenges remain before computational approaches can fill the role of workhorse in rational design of AF pharmacotherapies. In this review, we briefly detail the unique aspects of AF pathophysiology that determine requirements for compounds targeting AF rhythm control, with emphasis on delimiting mechanisms that promote AF triggers from those providing substrate or supporting reentry. We then describe modeling approaches that have been used to assess the outcomes of drugs acting on established AF targets, as well as on novel promising targets including the ultra-rapidly activating delayed rectifier potassium current, the acetylcholine-activated potassium current and the small conductance calcium-activated potassium channel. Finally, we describe how heterogeneity and variability are being incorporated into AF-specific models, and how these approaches are yielding novel insights into the basic physiology of disease, as well as aiding identification of the important molecular players in the complex AF etiology.

Action potential clamp characterization of the S631A hERG mutation associated with short QT syndrome

The hERG potassium channel is critical to normal repolarization of cardiac action potentials (APs) and loss- and gain-of-function hERG mutations are associated, respectively, with long and short QT syndromes, pathological conditions that can lead to arrhythmias and sudden death. hERG current (IhERG ) exhibits uniquely fast inactivation involving conformational changes to the channel pore. The S631A hERG pore mutation was originally engineered to interrogate hERG channel inactivation, but has very recently been found in a family with short QT syndrome (SQTS). Accordingly, this study characterized the effects of the S631A mutation on IhERG profile during ventricular, atrial, and Purkinje fiber (PF) AP waveforms, using patch clamp recording from hERG expressing HEK 293 cells at 37°C. Under conventional voltage clamp, the current-voltage (I-V) relation for IhERG exhibited a marked right-ward shift in the region of negative slope at positive membrane potentials. Under ventricular AP clamp, the S631A mutation resulted in augmented IhERG , which also peaked much earlier during the AP plateau than did wild-type (WT) IhERG . Instantaneous I-V relations showed a marked positive shift in peak repolarizing current during the ventricular AP in the S631A setting, while the instantaneous conductance-voltage relation showed an earlier and more sustained rise in S631A compared to WT IhERG conductance during ventricular repolarization. Experiments with atrial and PF APs in each case also showed augmented and positively shifted IhERG in the S631A setting, indicating that the S631A mutation is likely to accelerate repolarization in all three cardiac regions. Ventricular AP clamp experiments showed retained effectiveness of the class Ia antiarrhythmic drug quinidine (1 μmol/L) against S631A IhERG . Quinidine is thus likely to be effective in reducing excessively fast repolarization in SQTS resulting from the S631A hERG mutation.

Slow Recovery of Excitability Increases Ventricular Fibrillation Risk as Identified by Emulation

Purpose: Rotor stability and meandering are key mechanisms determining and sustaining cardiac fibrillation, with important implications for anti-arrhythmic drug development. However, little is yet known on how rotor dynamics are modulated by variability in cellular electrophysiology, particularly on kinetic properties of ion channel recovery.

Methods: We propose a novel emulation approach, based on Gaussian process regression augmented with machine learning, for data enrichment, automatic detection, classification, and analysis of re-entrant biomarkers in cardiac tissue. More than 5,000 monodomain simulations of long-lasting arrhythmic episodes with Fenton-Karma ionic dynamics, further enriched by emulation to 80 million electrophysiological scenarios, were conducted to investigate the role of variability in ion channel densities and kinetics in modulating rotor-driven arrhythmic behavior.

Results: Our methods predicted the class of excitation behavior with classification accuracy up to 96%, and emulation effectively predicted frequency, stability, and spatial biomarkers of functional re-entry. We demonstrate that the excitation wavelength interpretation of re-entrant behavior hides critical information about rotor persistence and devolution into fibrillation. In particular, whereas action potential duration directly modulates rotor frequency and meandering, critical windows of excitability are identified as the main determinants of breakup. Further novel electrophysiological insights of particular relevance for ventricular arrhythmias arise from our multivariate analysis, including the role of incomplete activation of slow inward currents in mediating tissue rate-dependence and dispersion of repolarization, and the emergence of slow recovery of excitability as a significant promoter of this mechanism of dispersion and increased arrhythmic risk.

Conclusions: Our results mechanistically explain pro-arrhythmic effects of class Ic anti-arrhythmics in the ventricles despite their established role in the pharmacological management of atrial fibrillation. This is mediated by their slow recovery of excitability mode of action, promoting incomplete activation of slow inward currents and therefore increased dispersion of repolarization, given the larger influence of these currents in modulating the action potential in the ventricles compared to the atria. These results exemplify the potential of emulation techniques in elucidating novel mechanisms of arrhythmia and further application to cardiac electrophysiology.

Determinants of slowed conduction in premature ventricular beats induced during programmed stimulations in perfused guinea-pig heart

NEW FINDINGS

What is the central question of this study? Is the slowed conduction upon premature ventricular activations during clinical electrophysiological testing attributable to the prolonged activation latency, or increased impulse propagation time, or both? What is the main finding and its importance? Prolonged activation latency at the stimulation site is the critical determinant of conduction slowing and associated changes in the ventricular response intervals in premature beats initiated during phase 3 repolarization in perfused guinea-pig heart. These relations are likely to have an effect on arrhythmia induction and termination independently of the presence of ventricular conduction defects or the proximity of the stimulation site to the re-entrant circuit.

ABSTRACT

During cardiac electrophysiological testing, slowed conduction upon premature ventricular activation can limit the delivery of the closely coupled impulses from the stimulation site to the region of tachycardia origin. In order to examine the contributing factors, in this study, cardiac conduction intervals and refractory periods were determined from left ventricular (LV) and the right ventricular (RV) monophasic action potential recordings obtained in perfused guinea-pig hearts. A premature activation induced immediately after the termination of the refractory period was associated with conduction slowing. The latter was primarily accounted for by the markedly increased (+54%) activation latency at the LV stimulation site, with only negligible changes (+12%) noted in the LV-to-RV delay. The prolonged activation latency was acting to limit the shortest interval at which two successive action potentials can be induced in the LV and RV chambers. The prolongation of the activation latency in premature beats was accentuated upon an increase in the stimulating current intensity, or during hypokalaemia. This change was related to the reduced ratio of the refractory period to the action potential duration, which allowed extrastimulus capture to occur earlier during phase 3 repolarization. Flecainide, a Na⁺ channel blocker, prolonged both the activation latency and the LV-to-RV delay, without changing their relative contributions to conduction slowing. In summary, these findings suggest that the activation latency is the critical determinant of conduction slowing and associated changes in the ventricular response intervals upon extrastimulus application during phase 3 of the action potential.

Antiarrhythmic Effects of Carvedilol and Flecainide in Cardiomyocytes Derived from Catecholaminergic Polymorphic Ventricular Tachycardia Patients

Mutations in the cardiac ryanodine receptor (RYR2) are the leading cause for catecholaminergic polymorphic ventricular tachycardia (CPVT). In this study, we evaluated antiarrhythmic efficacy of carvedilol and flecainide in CPVT patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) carrying different mutations in RYR2. iPSC-CMs were generated from skin biopsies of CPVT patients carrying exon 3 deletion and L4115 or V4653F mutation in RYR2 and of a healthy individual. Ca²⁺ kinetics and drug effects were studied with Fluo-4 AM indicator. Carvedilol abolished Ca²⁺ abnormalities in 31% of L4115F, 36% of V4653F, and 46% of exon 3 deletion carrying CPVT cardiomyocytes and flecainide 33%, 30%, and 52%, respectively. Both drugs lowered the intracellular Ca²⁺ level and beating rate of the cardiomyocytes significantly. Moreover, flecainide caused abnormal Ca²⁺ transients in 61% of controls compared to 26% of those with carvedilol. Carvedilol and flecainide were equally effective in CPVT iPSC-CMs. However, flecainide induced arrhythmias in 61% of control cells. CPVT cardiomyocytes carrying the exon 3 deletion had the most severe Ca²⁺ abnormalities, but they had the best response to drug therapies. According to this study, the arrhythmia-abolishing effect of neither of the drugs is optimal. iPSC-CMs provide a unique platform for testing drugs for CPVT.

Multiple targets for flecainide action: implications for cardiac arrhythmogenesis

Flecainide suppresses cardiac tachyarrhythmias including paroxysmal atrial fibrillation, supraventricular tachycardia and arrhythmic long QT syndromes (LQTS), as well as the Ca2+ -mediated, catecholaminergic polymorphic ventricular tachycardia (CPVT). However, flecainide can also exert pro-arrhythmic effects most notably following myocardial infarction and when used to diagnose Brugada syndrome (BrS). These divergent actions result from its physiological and pharmacological actions at multiple, interacting levels of cellular organization. These were studied in murine genetic models with modified Nav channel or intracellular ryanodine receptor (RyR2)-Ca2+ channel function. Flecainide accesses its transmembrane Nav 1.5 channel binding site during activated, open, states producing a use-dependent antagonism. Closing either activation or inactivation gates traps flecainide within the pore. An early peak INa related to activation of Nav channels followed by rapid de-activation, drives action potential (AP) upstrokes and their propagation. This is diminished in pro-arrhythmic conditions reflecting loss of function of Nav 1.5 channels, such as BrS, accordingly exacerbated by flecainide challenge. Contrastingly, pro-arrhythmic effects attributed to prolonged AP recovery by abnormal late INaL following gain-of-function modifications of Nav 1.5 channels in LQTS3 are reduced by flecainide. Anti-arrhythmic effects of flecainide that reduce triggering in CPVT models mediated by sarcoplasmic reticular Ca2+ release could arise from its primary actions on Nav channels indirectly decreasing [Ca2+ ]i through a reduced [Na+ ]i and/or direct open-state RyR2-Ca2+ channel antagonism. The consequent [Ca2+ ]i alterations could also modify AP propagation velocity and therefore arrhythmic substrate through its actions on Nav 1.5 channel function. This is consistent with the paradoxical differences between flecainide actions upon Na+ currents, AP conduction and arrhythmogenesis under circumstances of normal and increased RyR2 function.

LINKED ARTICLES

This article is part of a themed section on Spotlight on Small Molecules in Cardiovascular Diseases. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.8/issuetoc.

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.

Arrhythmogenic drugs can amplify spatial heterogeneities in the electrical restitution in perfused guinea-pig heart: An evidence from assessments of monophasic action potential durations and JT intervals

Non-uniform shortening of the action potential duration (APD90) in different myocardial regions upon heart rate acceleration can set abnormal repolarization gradients and promote arrhythmia. This study examined whether spatial heterogeneities in APD90 restitution can be amplified by drugs with clinically proved proarrhythmic potential (dofetilide, quinidine, procainamide, and flecainide) and, if so, whether these effects can translate to the appropriate changes of the ECG metrics of ventricular repolarization, such as JT intervals. In isolated, perfused guinea-pig heart preparations, monophasic action potentials and volume-conducted ECG were recorded at progressively increased pacing rates. The APD90 measured at distinct ventricular sites, as well as the JTpeak and JTend values were plotted as a function of preceding diastolic interval, and the maximum slopes of the restitution curves were determined at baseline and upon drug administration. Dofetilide, quinidine, and procainamide reverse rate-dependently prolonged APD90 and steepened the restitution curve, with effects being greater at the endocardium than epicardium, and in the right ventricular (RV) vs. the left ventricular (LV) chamber. The restitution slope was increased to a greater extent for the JTend vs. the JTpeak interval. In contrast, flecainide reduced the APD90 restitution slope at LV epicardium without producing effect at LV endocardium and RV epicardium, and reduced the JTpeak restitution slope without changing the JTend restitution. Nevertheless, with all agents, these effects translated to the amplified epicardial-to-endocardial and the LV-to-RV non-uniformities in APD90 restitution, paralleled by the increased JTend vs. JTpeak difference in the restitution slope. In summary, these findings suggest that arrhythmic drug profiles are partly attributable to the accentuated regional heterogeneities in APD90 restitution, which can be indirectly determined through ECG assessments of the JTend vs. JTpeak dynamics at variable pacing rates.

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

Repolarization of the cardiac action potential is primarily mediated by two voltage-dependent potassium currents: I_Kr and I_Ks. The voltage-gated potassium channel that gives rise to I_Kr, K_v11.1 (hERG), is uniquely susceptible to high-affinity block by a wide range of drug classes. Pore residues Tyr652 and Phe656 are critical to potent drug interaction with hERG. It is considered that the molecular basis of this broad-spectrum drug block phenomenon occurs through interactions specific to the aromatic nature of the side chains at Tyr652 and Phe656. In this study, we used nonsense suppression to incorporate singly and doubly fluorinated phenylalanine residues at Tyr652 and Phe656 to assess cation-π interactions in hERG terfenadine, quinidine, and dofetilide block. Incorporation of these unnatural amino acids was achieved with minimal alteration to channel activation or inactivation gating. Our assessment of terfenadine, quinidine, and dofetilide block did not reveal evidence of a cation-π interaction at either aromatic residue, but, interestingly, shows that certain fluoro-Phe substitutions at position 652 result in weaker  drug potency.

Frequency-Dependent Multi-Well Cardiotoxicity Screening Enabled by Optogenetic Stimulation

Side effects on cardiac ion channels causing lethal arrhythmias are one major reason for drug withdrawals from the market. Field potential (FP) recording from cardiomyocytes, is a well-suited tool to assess such cardiotoxic effects of drug candidates in preclinical drug development, but it is currently limited to the spontaneous beating of the cardiomyocytes and manual analysis. Herein, we present a novel optogenetic cardiotoxicity screening system suited for the parallel automated frequency-dependent analysis of drug effects on FP recorded from human-induced pluripotent stem cell-derived cardiomyocytes. For the expression of the light-sensitive cation channel Channelrhodopsin-2, we optimised protocols using virus transduction or transient mRNA transfection. Optical stimulation was performed with a new light-emitting diode lid for a 96-well FP recording system. This enabled reliable pacing at physiologically relevant heart rates and robust recording of FP. Thereby we detected rate-dependent effects of drugs on Na⁺, Ca2+ and K⁺ channel function indicated by FP prolongation, FP shortening and the slowing of the FP downstroke component, as well as generation of afterdepolarisations. Taken together, we present a scalable approach for preclinical frequency-dependent screening of drug effects on cardiac electrophysiology. Importantly, we show that the recording and analysis can be fully automated and the technology is readily available using commercial products.

Inhibitory effects and mechanism of dihydroberberine on hERG channels expressed in HEK293 cells

The human ether-a-go-go-related gene (hERG) potassium channel conducts rapid delayed rectifier potassium currents (I_Kr) and contributes to phase III cardiac action potential repolarization. Drugs inhibit hERG channels by binding to aromatic residues in hERG helixes. Berberine (BBR) has multiple actions, and its hydrogenated derivative dihydroberberine (DHB) is a potential candidate for developing new drugs. Previous studies have demonstrated that BBR blocks hERG channels and prolongs action potential duration (APD). Our present study aimed to investigate the effects and mechanism of DHB on hERG channels. Protein expression and the hERG current were analyzed using western blotting and patch-clamp, respectively. DHB inhibited the hERG current concentration-dependently after instantaneous perfusion, accelerated channel inactivation by directly binding tyrosine (Tyr652) and phenylalanine (Phe656), and decreased mature (155-kDa) and simultaneously increased immature (135-kDa) hERG expression, respectively. This suggests disruption of forward trafficking of hERG channels. Besides, DHB remarkably reduced heat shock protein 90 (Hsp90) expression and its interaction with hERG, indicating that DHB disrupted hERG trafficking by impairing channel folding. Meanwhie, DHB enhanced the expression of cleaved activating transcription factor-6 (ATF-6), a biomarker of unfolded protein response (UPR). Expression of calnexin and calreticulin, chaperones activated by ATF-6 to facilitate channel folding, were also increased, which indicating UPR activation. Additionally, the degradation rate of mature 155-kDa hERG increased following DHB exposure. In conclusion, we demonstrated that DHB acutely blocked hERG channels by binding the aromatic Tyr652 and Phe656. DHB may decrease hERG plasma membrane expression through two pathways involving disruption of forward trafficking of immature hERG channels and enhanced degradation of mature hERG channels. Furthermore, forward trafficking was disrupted by impaired channel folding associated with altered interactions between hERG proteins and chaperones. Finally, trafficking inhibition activated UPR, and mature hERG channel degradation was increased by DHB.

Role of abnormal repolarization in the mechanism of cardiac arrhythmia

In cardiac patients, life-threatening tachyarrhythmia is often precipitated by abnormal changes in ventricular repolarization and refractoriness. Repolarization abnormalities typically evolve as a consequence of impaired function of outward K⁺ currents in cardiac myocytes, which may be caused by genetic defects or result from various acquired pathophysiological conditions, including electrical remodelling in cardiac disease, ion channel modulation by clinically used pharmacological agents, and systemic electrolyte disorders seen in heart failure, such as hypokalaemia. Cardiac electrical instability attributed to abnormal repolarization relies on the complex interplay between a provocative arrhythmic trigger and vulnerable arrhythmic substrate, with a central role played by the excessive prolongation of ventricular action potential duration, impaired intracellular Ca²⁺ handling, and slowed impulse conduction. This review outlines the electrical activity of ventricular myocytes in normal conditions and cardiac disease, describes classical electrophysiological mechanisms of cardiac arrhythmia, and provides an update on repolarization-related surrogates currently used to assess arrhythmic propensity, including spatial dispersion of repolarization, activation-repolarization coupling, electrical restitution, TRIaD (triangulation, reverse use dependence, instability, and dispersion), and the electromechanical window. This is followed by a discussion of the mechanisms that account for the dependence of arrhythmic vulnerability on the location of the ventricular pacing site. Finally, the review clarifies the electrophysiological basis for cardiac arrhythmia produced by hypokalaemia, and gives insight into the clinical importance and pathophysiology of drug-induced arrhythmia, with particular focus on class Ia (quinidine, procainamide) and Ic (flecainide) Na⁺ channel blockers, and class III antiarrhythmic agents that block the delayed rectifier K⁺ channel (dofetilide).

Assessments of the QT/QRS restitution in perfused guinea-pig heart can discriminate safe and arrhythmogenic drugs

INTRODUCTION

Drug-induced arrhythmia remains a matter of serious clinical concern, partly due to low prognostic value of currently available arrhythmic biomarkers.

METHODS

This study examined whether arrhythmogenic risks can be predicted through assessments of the rate adaptation of QT interval, ventricular effective refractory period (ERP), or the QT/QRS ratio, in perfused guinea-pig hearts.

RESULTS

When the maximum restitution slope was taken as a metric of proarrhythmia, neither QT interval nor ERP measurements at progressively increased pacing rates were found to fully discriminate arrhythmogenic drugs (dofetilide, quinidine, flecainide, and procainamide) from those recognized as safe antiarrhythmics (lidocaine and mexiletine). For example, the slope of QT restitution was increased by dofetilide and quinidine, but remained unchanged by flecainide, procainamide, lidocaine, and mexiletine. With ERP rate adaptation, even though the restitution slope was increased by dofetilide, all class I agents reduced the slope value independently of their safety profile. The QRS measurements revealed variable drug effects, ranging from significant use-dependent conduction slowing (flecainide, quinidine, and procainamide) to only modest increase in QRS (lidocaine and mexiletine), or no change at all (dofetilide). However, with the QT/QRS rate adaptation, the restitution slope was significantly increased by all agents which have been reported to produce proarrhythmic effects (dofetilide, quinidine, flecainide, and procainamide), but not changed by lidocaine and mexiletine.

DISCUSSION

These findings suggest that the slope of the QT/QRS rate adaptation can be considered as a novel electrophysiological biomarker in predicting potential arrhythmic risks associated with pharmacotherapy in cardiac patients.

Non-invasive phenotyping and drug testing in single cardiomyocytes or beta-cells by calcium imaging and optogenetics

Identification of drug induced electrical instability of the heart curtails development, and introduction, of potentially proarrhythmic drugs. This problem usually requires complimentary contact based approaches such as patch-clamp electrophysiology combined with field stimulation electrodes to observe and control the cell. This produces data with high signal to noise but requires direct physical contact generally preventing high-throughput, or prolonged, phenotyping of single cells or tissues. Combining genetically encoded optogenetic control and spectrally compatible calcium indicator tools into a single adenoviral vector allows the analogous capability for cell control with simultaneous cellular phenotyping without the need for contact. This combination can be applied to single rodent primary adult cardiomyocytes, and human stem cell derived cardiomyocytes, enabling contactless small molecule evaluation for inhibitors of sodium, potassium and calcium channels suggesting it may be useful for early toxicity work. In pancreatic beta-cells it reveals the effects of glucose and the KATP inhibitor gliclazide.

Effect of flecainide derivatives on sarcoplasmic reticulum calcium release suggests a lack of direct action on the cardiac ryanodine receptor

BACKGROUND AND PURPOSE

Flecainide is a use-dependent blocker of cardiac Na(+) channels. Mechanistic analysis of this block showed that the cationic form of flecainide enters the cytosolic vestibule of the open Na(+) channel. Flecainide is also effective in the treatment of catecholaminergic polymorphic ventricular tachycardia but, in this condition, its mechanism of action is contentious. We investigated how flecainide derivatives influence Ca(2) (+) -release from the sarcoplasmic reticulum through the ryanodine receptor channel (RyR2) and whether this correlates with their effectiveness as blockers of Na(+) and/or RyR2 channels.

EXPERIMENTAL APPROACH

We compared the ability of fully charged (QX-FL) and neutral (NU-FL) derivatives of flecainide to block individual recombinant human RyR2 channels incorporated into planar phospholipid bilayers, and their effects on the properties of Ca(2) (+) sparks in intact adult rat cardiac myocytes.

KEY RESULTS

Both QX-FL and NU-FL were partial blockers of the non-physiological cytosolic to luminal flux of cations through RyR2 channels but were significantly less effective than flecainide. None of the compounds influenced the physiologically relevant luminal to cytosol cation flux through RyR2 channels. Intracellular flecainide or QX-FL, but not NU-FL, reduced Ca(2) (+) spark frequency.

CONCLUSIONS AND IMPLICATIONS

Given its inability to block physiologically relevant cation flux through RyR2 channels, and its lack of efficacy in blocking the cytosolic-to-luminal current, the effect of QX-FL on Ca(2) (+) sparks is likely, by analogy with flecainide, to result from Na(+) channel block. Our data reveal important differences in the interaction of flecainide with sites in the cytosolic vestibules of Na(+) and RyR2 channels.

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.