Rapid Delayed Rectifier Potassium Current Research Articles

Abstract We142: Identification of the Molecular Determinants for the Circadian Regulation of the Human KCNH2 promoter

Background: KCNH2 encodes a voltage-gated potassium channel (Kv11.1), which conducts the rapid delayed rectifier potassium current (IKr) critical for ventricular repolarization in the heart. Circadian clock factors regulate the expression of the KCNH2 mRNA transcript in the heart. However, the molecular mechanisms for the circadian expression of KCNH2 transcripts have yet to be defined. Hypothesis: A conserved cis-regulatory element in the core KCNH2 promoter drives circadian regulation of human KCNH2 (hKCNH2) promoter activity. Objective: Identify the molecular determinants of the circadian regulation of hKCNH2 promoter activity and determine the impact of single nucleotide polymorphisms (SNPs) in this region. Methods: Real-time bioluminescence (LumiCycle, Actimetrics) and Dual Glo luciferase assays were employed to identify critical cis elements that regulate the circadian expression of the human KCNH2 promoter using promoter-luciferase reporter constructs. We studied several deletion and mutant hKCNH2 promoter reporter vectors transiently transfected in a myogenic cell line (C2C12 cells). Bioluminescence was measured at 10-minute intervals for 7-10 days. Data were analyzed for circadian characteristics (period, phase and amplitude). Mutant constructs included both engineered and constructs containing naturally occurring rare and common SNPs. Results: Lumicycle data demonstrated that the hKCNH2 promoter reporter showed robust circadian and overall activity in C2C12 cells. Deletion and mutational analysis of the hKCNH2 promoter-reporter construct identified a highly conserved tandem E-box element within 1 Kb of the exon 1 start site that was critical for circadian and overall regulation of hKCNH2 promoter activity. Cells expressing hKCNH2 promoter-reporter constructs with different SNPs in this tandem E-box element showed SNP-specific effects. Surprisingly, a relatively common SNP decreased hKCNH2 circadian promoter activity by 66% (n = 6, p=0.0049). Conclusion: We identified a conserved tandem E-box in the proximal promoter of hKCNH2 Exon1 necessary for circadian and overall hKCNH2 promoter activity. Several SNPs in this region suggested a functional impact that might lead to decreased levels of IKr.

Recording ten-fold larger IKr conductances with automated patch clamping using equimolar Cs+ solutions.

Background: The rapid delayed rectifier potassium current (IKr) is important for cardiac repolarization and is most often involved in drug-induced arrhythmias. However, accurately measuring this current can be challenging in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes because of its small current density. Interestingly, the ion channel conducting IKr, hERG channel, is not only permeable to K+ ions but also to Cs+ ions when present in equimolar concentrations inside and outside of the cell. Methods: In this study, IhERG was measured from Chinese hamster ovary (CHO)-hERG cells and hiPSC-CM using either Cs+ or K+ as the charge carrier. Equimolar Cs+ has been used in the literature in manual patch-clamp experiments, and here, we apply this approach using automated patch-clamp systems. Four different (pre)clinical drugs were tested to compare their effects on Cs+- and K+-based currents. Results: Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. Comparison of Cs+- and K+-mediated currents upon application of dofetilide, desipramine, moxifloxacin, or LUF7244 revealed many similarities in inhibition or activation properties of the drugs studied. Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. In hiPSC-CM, the Cs+-based conductance is larger compared to the known K+-based conductance, and the Cs+ hERG conductance can be inhibited similarly to the K+-based conductance. Conclusion: Using equimolar Cs+ instead of K+ for IhERG measurements in an automated patch-clamp system gives rise to a new method by which, for example, quick scans can be performed on effects of drugs on hERG currents. This application is specifically relevant when such experiments are performed using cells which express small IKr current densities in combination with small membrane capacitances.

Evaluating pro-arrhythmogenic effects of the T634S-hERG mutation: insights from a simulation study.

A mutation to serine of a conserved threonine (T634S) in the hERG K+ channel S6 pore region has been identified as a variant of uncertain significance, showing a loss-of-function effect. However, its potential consequences for ventricular excitation and arrhythmogenesis have not been reported. This study evaluated possible functional effects of the T634S-hERG mutation on ventricular excitation and arrhythmogenesis by using multi-scale computer models of the human ventricle. A Markov chain model of the rapid delayed rectifier potassium current (IKr) was reconstructed for wild-type and T634S-hERG mutant conditions and incorporated into the ten Tusscher et al. models of human ventricles at cell and tissue (1D, 2D and 3D) levels. Possible functional impacts of the T634S-hERG mutation were evaluated by its effects on action potential durations (APDs) and their rate-dependence (APDr) at the cell level; and on the QT interval of pseudo-ECGs, tissue vulnerability to unidirectional conduction block (VW), spiral wave dynamics and repolarization dispersion at the tissue level. It was found that the T634S-hERG mutation prolonged cellular APDs, steepened APDr, prolonged the QT interval, increased VW, destablized re-entry and augmented repolarization dispersion across the ventricle. Collectively, these results imply potential pro-arrhythmic effects of the T634S-hERG mutation, consistent with LQT2.

Ionic current changes underlying action potential repolarization responses to physiological pacing and adrenergic stimulation in adult rat ventricular myocytes.

This study aimed to simulate ventricular responses to elevations in myocyte pacing and adrenergic stimulation using a novel electrophysiological rat model and investigate ion channel responses underlying action potential (AP) modulations. Peak ion currents and AP repolarization to 50% and 90% of full repolarization (APD50-90 ) were recorded during simulations at 1-10 Hz pacing under control and adrenergic stimulation conditions. Further simulations were performed with incremental ion current block (L-type calcium current, ICa ; transient outward current, Ito ; slow delayed rectifier potassium current, IKs ; rapid delayed rectifier potassium current, IKr ; inward rectifier potassium current, IK1 ) to identify current influence on AP response to exercise. Simulated APD50-90 closely resembled experimental findings. Rate-dependent increases in IKs (6%-101%), IKr (141%-1339%), and ICa (0%-15%) and reductions in Ito (11%-57%) and IK1 (1%-9%) were observed. Meanwhile, adrenergic stimulation triggered moderate increases in all currents (23%-67%) except IK1 . Further analyses suggest AP plateau is most sensitive to modulations in Ito and ICa while late repolarization is most sensitive to IK1 , ICa , and IKs , with alterations in IKs predominantly stimulating the greatest magnitude of influence on late repolarization (35%-846% APD90 prolongation). The modified Leeds rat model (mLR) is capable of accurately modeling APs during physiological stress. This study highlights the importance of ICa , Ito , IK1, and IKs in controlling electrophysiological responses to exercise. This work will benefit the study of cardiac dysfunction, arrythmia, and disease, though future physiologically relevant experimental studies and model development are required.

Importance of modelling hERG binding in predicting drug-induced action potential prolongations for drug safety assessment.

Reduction of the rapid delayed rectifier potassium current (I Kr) via drug binding to the human Ether-à-go-go-Related Gene (hERG) channel is a well recognised mechanism that can contribute to an increased risk of Torsades de Pointes. Mathematical models have been created to replicate the effects of channel blockers, such as reducing the ionic conductance of the channel. Here, we study the impact of including state-dependent drug binding in a mathematical model of hERG when translating hERG inhibition to action potential changes. We show that the difference in action potential predictions when modelling drug binding of hERG using a state-dependent model versus a conductance scaling model depends not only on the properties of the drug and whether the experiment achieves steady state, but also on the experimental protocols. Furthermore, through exploring the model parameter space, we demonstrate that the state-dependent model and the conductance scaling model generally predict different action potential prolongations and are not interchangeable, while at high binding and unbinding rates, the conductance scaling model tends to predict shorter action potential prolongations. Finally, we observe that the difference in simulated action potentials between the models is determined by the binding and unbinding rate, rather than the trapping mechanism. This study demonstrates the importance of modelling drug binding and highlights the need for improved understanding of drug trapping which can have implications for the uses in drug safety assessment.

KCNH2 encodes a nuclear-targeted polypeptide that mediates hERG1 channel gating and expression.

KCNH2 encodes hERG1, the voltage-gated potassium channel that conducts the rapid delayed rectifier potassium current (IKr) in human cardiac tissue. hERG1 is one of the first channels expressed during early cardiac development, and its dysfunction is associated with intrauterine fetal death, sudden infant death syndrome, cardiac arrhythmia, and sudden cardiac death. Here, we identified a novel hERG1 polypeptide (hERG1NP) that is targeted to the nuclei of immature cardiac cells, including hiPSC-CMs and neonatal rat cardiomyocytes. The nuclear hERG1NP immunofluorescent signal is diminished in matured hiPSC-CMs and absent from adult rat cardiomyocytes. Antibodies targeting distinct hERG1 channel epitopes demonstrated that the hERG1NP signal maps to the hERG1 distal C-terminal domain. KCNH2 deletion using CRISPR simultaneously abolished IKr and the hERG1NP signal in hiPSC-CMs. We then identified a putative nuclear localization sequence (NLS) within the distal hERG1 C-terminus, 883-RQRKRKLSFR-892. Interestingly, the distal hERG1 C-terminal domain was targeted almost exclusively to the nuclei when overexpressed HEK293 cells. Conversely, deleting the NLS from the distal peptide abolished nuclear targeting. Similarly, blocking α or β1 karyopherin activity diminished nuclear targeting. Finally, overexpressing the distal hERG1 C-terminus significantly reduced hERG1a current density and slowed deactivation, compared to cells expressing GFP. Deleting the NLS from the distal hERG1 C-terminal domain abolished the effects of the distal C-terminus on hERG1a current amplitude and kinetics. These data identify a developmentally regulated polypeptide encoded by KCNH2, hERGNP, whose presence in the nucleus indirectly modulates hERG1 current magnitude and kinetics.

KCNH2 encodes a nuclear-targeted polypeptide that mediates hERG1 channel gating and expression

KCNH2 encodes hERG1, the voltage-gated potassium channel that conducts the rapid delayed rectifier potassium current (IKr) in human cardiac tissue. hERG1 is one of the first channels expressed during early cardiac development, and its dysfunction is associated with intrauterine fetal death, sudden infant death syndrome, cardiac arrhythmia, and sudden cardiac death. Here, we identified a hERG1 polypeptide (hERG1NP) that is targeted to the nuclei of immature cardiac cells, including human stem cell-derived cardiomyocytes (hiPSC-CMs) and neonatal rat cardiomyocytes. The nuclear hERG1NP immunofluorescent signal is diminished in matured hiPSC-CMs and absent from adult rat cardiomyocytes. Antibodies targeting distinct hERG1 channel epitopes demonstrated that the hERG1NP signal maps to the hERG1 distal C-terminal domain. KCNH2 deletion using CRISPR simultaneously abolished IKr and the hERG1NP signal in hiPSC-CMs. We then identified a putative nuclear localization sequence (NLS) within the distal hERG1 C-terminus, 883-RQRKRKLSFR-892. Interestingly, the distal C-terminal domain was targeted almost exclusively to the nuclei when overexpressed HEK293 cells. Conversely, deleting the NLS from the distal peptide abolished nuclear targeting. Similarly, blocking α or β1 karyopherin activity diminished nuclear targeting. Finally, overexpressing the putative hERG1NP peptide in the nuclei of HEK cells significantly reduced hERG1a current density, compared to cells expressing the NLS-deficient hERG1NP or GFP. These data identify a developmentally regulated polypeptide encoded by KCNH2, hERG1NP, whose presence in the nucleus indirectly modulates hERG1 current magnitude and kinetics.

In silico mechanisms of arsenic trioxide-induced cardiotoxicity.

It has been found that arsenic trioxide (ATO) is effective in treating acute promyelocytic leukemia (APL). However, long QT syndrome was reported in patients receiving therapy using ATO, which even led to sudden cardiac death. The underlying mechanisms of ATO-induced cardiotoxicity have been investigated in some biological experiments, showing that ATO affects human ether-à-go-go-related gene (hERG) channels, coding rapid delayed rectifier potassium current (I Kr ), as well as L-type calcium (I CaL ) channels. Nevertheless, the mechanism by which these channel reconstitutions induced the arrhythmia in ventricular tissue remains unsolved. In this study, a mathematical model was developed to simulate the effect of ATO on ventricular electrical excitation at cellular and tissue levels by considering ATO's effects on I Kr and I CaL . The ATO-dose-dependent pore block model was incorporated into the I Kr model, and the enhanced degree of ATO to I CaL was based on experimental data. Simulation results indicated that ATO extended the action potential duration of three types of ventricular myocytes (VMs), including endocardial cells (ENDO), midmyocardial cells (MCELL), and epicardial cells (EPI), and exacerbated the heterogeneity among them. ATO could also induce alternans in all three kinds of VMs. In a cable model of the intramural ventricular strand, the effects of ATO are reflected in a prolonged QT interval of simulated pseudo-ECG and a wide vulnerable window, thus increasing the possibility of spiral wave formation in ventricular tissue. In addition to showing that ATO prolonged QT, we revealed that the heterogeneity caused by ATO is also an essential hazard factor. Based on this, a pharmacological intervention of ATO toxicity by resveratrol was undertaken. This study provides a further understanding of ATO-induced cardiotoxicity, which may help to improve the treatment for APL patients.

HERG1 channel subunit composition mediates proton inhibition of rapid delayed rectifier potassium current (IKr) in cardiomyocytes derived from hiPSCs

The voltage-gated channel, hERG1, conducts the rapid delayed rectifier potassium current (IKr) and is critical for human cardiac repolarization. Reduced IKr causes long QT syndrome and increases the risk for cardiac arrhythmia and sudden death. At least two subunits form functional hERG1 channels, hERG1a and hERG1b. Changes in hERG1a/1b abundance modulate IKr kinetics, magnitude, and drug sensitivity. Studies from native cardiac tissue suggest that hERG1 subunit abundance is dynamically regulated, but the impact of altered subunit abundance on IKr and its response to external stressors is not well understood. Here, we used a substrate-driven human-induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM) maturation model to investigate how changes in relative hERG1a/1b subunit abundance impact the response of native IKr to extracellular acidosis, a known component of ischemic heart disease and sudden infant death syndrome. IKr recorded from immatured hiPSC-CMs displays a 2-fold greater inhibition by extracellular acidosis (pH 6.3) compared with matured hiPSC-CMs. Quantitative RT-PCR and immunocytochemistry demonstrated that hERG1a subunit mRNA and protein were upregulated and hERG1b subunit mRNA and protein were downregulated in matured hiPSC-CMs compared with immatured hiPSC-CMs. The shift in subunit abundance in matured hiPSC-CMs was accompanied by increased IKr. Silencing hERG1b's impact on native IKr kinetics by overexpressing a polypeptide identical to the hERG1a N-terminal Per-Arnt-Sim domain reduced the magnitude of IKr proton inhibition in immatured hiPSC-CMs to levels comparable to those observed in matured hiPSC-CMs. These data demonstrate that hERG1 subunit abundance is dynamically regulated and determines IKr proton sensitivity in hiPSC-CMs.

Interleukin-6 Elevation Is a Key Pathogenic Factor Underlying COVID-19-Associated Heart Rate-Corrected QT Interval Prolongation.

BackgroundHeart rate-corrected QT interval (QTc) prolongation is prevalent in patients with severe coronavirus disease 2019 (COVID-19) and is associated with poor outcomes. Recent evidence suggests that the exaggerated host immune-inflammatory response characterizing the disease, specifically interleukin-6 (IL-6) increase, may have an important role, possibly via direct effects on cardiac electrophysiology. The aim of this study was to dissect the short-term discrete impact of IL-6 elevation on QTc in patients with severe COVID-19 infection and explore the underlying mechanisms.MethodsWe investigated the following mechanisms: (1) the QTc duration in patients with COVID-19 during the active phase and recovery, and its association with C-reactive protein (CRP) and IL-6 levels; (2) the acute impact of IL-6 administration on QTc in an in vivo guinea pig model; and (3) the electrophysiological effects of IL-6 on ventricular myocytes in vitro.ResultsIn patients with active severe COVID-19 and elevated IL-6 levels, regardless of acute myocardial injury/strain and concomitant QT-prolonging risk factors, QTc was significantly prolonged and rapidly normalized in correlation with IL-6 decrease. The direct administration of IL-6 in an in vivo guinea pig model acutely prolongs QTc duration. Moreover, ventricular myocytes incubated in vitro with IL-6 show evident prolongation in the action potential, along with significant inhibition in the rapid delayed rectifier potassium current (IKr).ConclusionFor the first time, we demonstrated that in severe COVID-19, systemic inflammatory activation can per se promote QTc prolongation via IL-6 elevation, leading to ventricular electric remodeling. Despite being transitory, such modifications may significantly contribute to arrhythmic events and associated poor outcomes in COVID-19. These findings provide a further rationale for current anti-inflammatory treatments for COVID-19, including IL-6-targeted therapies.

Identification and validation of a nuclear localization sequence in the hERG1 potassium channel

KCNH2 encodes hERG1, the voltage-gated potassium channel that conducts the rapid delayed rectifier potassium current (IKr) in cardiac tissue. Reduced IKr causes the cardiac disorder long QT syndrome (LQTS), and increases the risk of cardiac arrhythmia and sudden cardiac death. We observed a hERG1-specific immunofluorescent signal within the nucleus of human stem cell-derived cardiomyocytes and neonatal rat cardiomyocytes. Using antibodies targeting distinct hERG1 channel epitopes, we mapped the nuclear hERG1 peptide to the channel's cytoplasmic C-terminal domain.

Prediction of arrhythmia susceptibility through mathematical modeling and machine learning.

At present, the QT interval on the electrocardiographic (ECG) waveform is the most common metric for assessing an individual's susceptibility to ventricular arrhythmias, with a long QT, or, at the cellular level, a long action potential duration (APD) considered high risk. However, the limitations of this simple approach have long been recognized. Here, we sought to improve prediction of arrhythmia susceptibility by combining mechanistic mathematical modeling with machine learning (ML). Simulations with a model of the ventricular myocyte were performed to develop a large heterogenous population of cardiomyocytes (n = 10,586), and we tested each variant's ability to withstand three arrhythmogenic triggers: 1) block of the rapid delayed rectifier potassium current (IKr Block), 2) augmentation of the L-type calcium current (ICaL Increase), and 3) injection of inward current (Current Injection). Eight ML algorithms were trained to predict, based on simulated AP features in preperturbed cells, whether each cell would develop arrhythmic dynamics in response to each trigger. We found that APD can accurately predict how cells respond to the simple Current Injection trigger but cannot effectively predict the response to IKr Block or ICaL Increase. ML predictive performance could be improved by incorporating additional AP features and simulations of additional experimental protocols. Importantly, we discovered that the most relevant features and experimental protocols were trigger specific, which shed light on the mechanisms that promoted arrhythmia formation in response to the triggers. Overall, our quantitative approach provides a means to understand and predict differences between individuals in arrhythmia susceptibility.

Genetic Algorithm For Fitting Cardiac Cell Biophysical Model Formulations.

Modeling cardiac cell electrophysiology relies on fitting model equations to experimental data obtained under voltage/current clamping conditions. The fitting procedure for these often-nonlinear ionic current equations are mostly executed by trial-and-error by hand or by gradient-based optimization approaches. These methods, though sometimes sufficient at converging at optimal solutions is based on the premise that the characteristic objective function is convex, which often does not apply to cardiac model equations. Meta-heuristic methods, such as evolutionary algorithms and particle swarm algorithms, have proven resilient against early convergence to local optima and saddle-point parameter solutions. This work presents a genetic algorithm-based approach for fitting the adult cardiomyocyte biophysical model formulations to the experimental data obtained in human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM). Specifically, whole-cell patch clamp ionic current data of rapid delayed rectifier potassium current, IKr, transient outward potassium current, Ito and hyperpolarization-activated current, If, was used for fitting. Using a two-point crossover scheme along with initial population and mutation constraints randomly selected from a uniformly distributed constrained parameter space, near-optimal fitting was achieved with R2 values (n = 5) of 0.9960±0.0007, 0.9995±0.0002, and 0.9974±0.0014 for IKr, Ito and If respectively.

Rapid Characterisation of R56Q Mutant hERG Channel Kinetics using Sinusoidal Voltage Protocols

The voltage-gated human-Ether-à-go-go-Related Gene (hERG) potassium channel is the molecular correlate of the rapid delayed rectifier potassium current, IKr, which is an important determinant of repolarisation of the action potential in cardiac myocytes. Loss of function in hERG due to inherited mutations results in Long QT Syndrome Type 2 (LQT2), which can result in the potentially fatal arrhythmia, Torsades de Pointes. However, our ability to predict how such mutations influence cellular electrophysiology and thus behaviour at larger scales depends critically on a faithful characterisation of ion channel kinetics. In this study, we tested the ability of a previously-developed sinusoidal voltage clamp protocol to characterise kinetics of the LQT2-related R56Q hERG mutation. Membrane currents were recorded in Xenopus laevis oocytes injected with cRNA encoding wild type (WT) or R56Q mutant hERG channels at room temperature using the two-electrode voltage clamp technique. We optimised a mathematical model to whole-cell current traces evoked by the sinusoidal voltage clamp protocol, as well as traditional activation and inactivation square step voltage clamp protocols. The models were then used to predict the response to a premature stimulation protocol, where R56Q hERG showed smaller protective currents in response to premature stimuli than WT due to accelerated channel deactivation. Furthermore, the models were able to predict the response to a complex series of physiologically-relevant action potential waveforms with quantitative accuracy. Our study demonstrates the broader applicability of sinusoidal voltage protocols for rapid characterisation of ion channel kinetics, extending the approach to mutations in hERG which underlie congenital rhythm disorders. The resulting predictive models can be used for improved assessment of clinical phenotypes arising from abnormalities in hERG, due to genetic mutation or pharmacological modulation.

High-throughput phenotyping of heteromeric human ether-à-go-go-related gene potassium channel variants can discriminate pathogenic from rare benign variants

KCNH2 encodes the human ether-à-go-go-related gene potassium channel, which passes the rapid delayed rectifier potassium current. Loss-of-function variants in KCNH2 cause long QT syndrome type 2, which is associated with a markedly increased risk of cardiac arrhythmias. The majority of rare KCNH2 variants, however, are likely to be benign. The purpose of this study was to develop a high-throughput assay for discriminating pathogenic from benign KCNH2 variants. Nonsynonymous homozygous KCNH2 variants stably expressed in Flp-In human embryonic kidney 293 cell lines were phenotyped using an automated patch-clamp platform and a cell surface enzyme-linked immunosorbent assay. Functional phenotyping of heterozygous KCNH2 variants stably expressed in Flp-In human embryonic kidney 293 cell lines using a bicistronic vector was performed using an automated patch-clamp platform. In homozygous KCNH2 variant cell lines, discrepancies between current density and cell surface expression levels measured using an enzyme-linked immunosorbent assay can be explained by changes in gating properties of the variant channels. For the 30 heterozygous KCNH2 variant cell lines studied, the assay correctly predicted the ClinVar ascribed classification for 17/17 pathogenic/likely pathogenic/benign variants. Of the 13 pore-domain variants studied, 11 had a dominant-negative expression defect while the remaining 2 had enhanced inactivation gating, resulting in a dominant-negative phenotype. High-throughput electrophysiological phenotyping of heterozygous KCNH2 variants can accurately distinguish between dominant-negative, haploinsufficient loss-of-function, and benign variants. This assay will help with future classification of KCNH2 variants.

In silico Assessment of Pharmacotherapy for Human Atrial Patho-Electrophysiology Associated With hERG-Linked Short QT Syndrome.

Short QT syndrome variant 1 (SQT1) arises due to gain-of-function mutations to the human Ether-à-go-go-Related Gene (hERG), which encodes the α subunit of channels carrying rapid delayed rectifier potassium current, IKr. In addition to QT interval shortening and ventricular arrhythmias, SQT1 is associated with increased risk of atrial fibrillation (AF), which is often the only clinical presentation. However, the underlying basis of AF and its pharmacological treatment remain incompletely understood in the context of SQT1. In this study, computational modeling was used to investigate mechanisms of human atrial arrhythmogenesis consequent to a SQT1 mutation, as well as pharmacotherapeutic effects of selected class I drugs–disopyramide, quinidine, and propafenone. A Markov chain formulation describing wild type (WT) and N588K-hERG mutant IKr was incorporated into a contemporary human atrial action potential (AP) model, which was integrated into one-dimensional (1D) tissue strands, idealized 2D sheets, and a 3D heterogeneous, anatomical human atria model. Multi-channel pharmacological effects of disopyramide, quinidine, and propafenone, including binding kinetics for IKr/hERG and sodium current, INa, were considered. Heterozygous and homozygous formulations of the N588K-hERG mutation shortened the AP duration (APD) by 53 and 86 ms, respectively, which abbreviated the effective refractory period (ERP) and excitation wavelength in tissue, increasing the lifespan and dominant frequency (DF) of scroll waves in the 3D anatomical human atria. At the concentrations tested in this study, quinidine most effectively prolonged the APD and ERP in the setting of SQT1, followed by disopyramide and propafenone. In 2D simulations, disopyramide and quinidine promoted re-entry termination by increasing the re-entry wavelength, whereas propafenone induced secondary waves which destabilized the re-entrant circuit. In 3D simulations, the DF of re-entry was reduced in a dose-dependent manner for disopyramide and quinidine, and propafenone to a lesser extent. All of the anti-arrhythmic agents promoted pharmacological conversion, most frequently terminating re-entry in the order quinidine > propafenone = disopyramide. Our findings provide further insight into mechanisms of SQT1-related AF and a rational basis for the pursuit of combined IKr and INa block based pharmacological strategies in the treatment of SQT1-linked AF.

An exploration of interactions between the antiarrhythmic drug dronedarone and hERG potassium channel pore

Dronedarone is a non-iodinated analogue of the Class III antiarrhythmic agent amiodarone. It exerts potent inhibition of “hERG” potassium channels that underpin the cardiac rapid delayed rectifier potassium current, IKr. This study aimed to extend understanding of interactions between dronedarone and the hERG channel. Whole-cell patch-clamp recordings were made at 37?C of hERG channel current (IhERG) from HEK-293 cells expressing wild-type (WT) hERG or alanine mutants of residues in the channel’s pore-helix/selectivity filter region (T623, S624, V625) or S6 helices (S649, Y652, F656, V659). Molecular docking simulations were performed using a cryo-EM structure of hERG and a MthK-based homology model. The half-maximal inhibitory (IC50) value for WT IhERG inhibition by dronedarone was 42.6 ± 3.9 nM (n= at least 5 cells for each of 6 concentrations). 600 nM dronedarone exerted reduced WT IhERG block when the direction of K+ flux was reversed, consistent with interactions between the drug and permeant ion. In contrast with recently reported data for amiodarone, the S624A mutation did not attenuate IhERG blockade, whilst T623A and V625A channels exhibited modestly attenuated block. The S649A mutation was without significant effect and the Y652A and F656A mutations exhibited modest reductions in block. The V659A mutation produced the most marked effect on dronedarone action. Docking simulations were generally consistent with modest interactions with canonical binding residues and suggested an indirect rather than direct effect of the V659A mutation on the drug’s action. These findings leave open the possibility that as yet unexplored residue(s) could act as key determinants of high affinity hERG channel block by dronedarone.

Molecular determinants of Kv7.1/KCNE1 channel inhibition by amitriptyline

Amitriptyline (AMIT) is a compound widely prescribed for psychiatric and non-psychiatric conditions including depression, migraine, chronic pain, and anorexia. However, AMIT has been associated with risks of cardiac arrhythmia and sudden death since it can induce prolongation of the QT interval on the surface electrocardiogram and torsade de pointes ventricular arrhythmia. These complications have been attributed to the inhibition of the rapid delayed rectifier potassium current (IKr). The slow delayed rectifier potassium current (IKs) is the main repolarizing cardiac current when IKr is compromised and it has an important role in cardiac repolarization at fast heart rates induced by an elevated sympathetic tone. Therefore, we sought to characterize the effects of AMIT on Kv7.1/KCNE1 and homomeric Kv7.1 channels expressed in HEK-293H cells. Homomeric Kv7.1 and Kv7.1/KCNE1 channels were inhibited by AMIT in a concentration-dependent manner with IC50 values of 8.8 ± 2.1 μM and 2.5 ± 0.8 μM, respectively. This effect was voltage-independent for both homomeric Kv7.1 and Kv7.1/KCNE1 channels. Moreover, mutation of residues located on the P-loop and S6 domain along with molecular docking, suggest that T312, I337 and F340 are the most important molecular determinants for AMIT-Kv7.1 channel interaction. Our experimental findings and modeling suggest that AMIT preferentially blocks the open state of Kv7.1/KCNE1 channels by interacting with specific residues that were previously reported to be important for binding of other compounds, such as chromanol 293B and the benzodiazepine L7.

Optimization of an In silico Cardiac Cell Model for Proarrhythmia Risk Assessment.

Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive in vitro Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using in silico models and in vitro multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model). We also developed a novel metric that could separate drugs with different TdP liabilities at high concentrations based on total electronic charge carried by the major inward ionic currents during the action potential. In this study, we further optimized the IKr-dynamic ORd model by refining model parameters using published human cardiomyocyte experimental data under control and drug block conditions. Using this optimized model and manual patch clamp data, we developed an updated version of the metric that quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state action potential, which could classify the level of drug-induced TdP risk across a wide range of concentrations and pacing rates. We also established a framework to quantitatively evaluate a system's robustness against the induction of early afterdepolarizations (EADs), and demonstrated that the new metric is correlated with the cell's robustness to the pro-EAD perturbation of IKr conductance reduction. In summary, in this work we present an optimized model that is more consistent with experimental data, an improved metric that can classify drugs at concentrations both near and higher than clinical exposure, and a physiological framework to check the relationship between a metric and EAD. These findings provide a solid foundation for using in silico models for the regulatory assessment of TdP risk under the CiPA paradigm.

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