Effects of IKr and IKs heterogeneity on action potential duration and its rate dependence: a simulation study

ST segment and T wave abnormalities: A narrative review

INTRODUCTION

The electrocardiogram (ECG) is a valuable tool for interpreting ventricular repolarization. This article aims to broaden the diagnostic scope beyond the conventional ischemia-centric approach, integrating an understanding of pathophisiological influences on ST-T wave changes.

METHODS

A review was conducted on the physiological underpinnings of ventricular repolarization and the pathophisiological processes that can change ECG patterns. The research encompassed primary repolarization abnormalities due to uniform variations in ventricular action potential, secondary changes from electrical or mechanical alterations, and non-ischemic conditions influencing ST-T segments.

RESULTS

Primary T waves are characterized by symmetrical waves with broad bases and variable QT intervals, indicative of direct myocardial action potential modifications due to ischemia, electrolyte imbalances, and channelopathies. Secondary T waves are asymmetric and often unassociated with significant QT interval changes, suggesting depolarization alterations or changes in cardiac geometry and contractility.

CONCLUSION

We advocate for a unified ECG analysis, recognizing primary and secondary ST-T changes, and their clinical implications. Our proposed analytical framework enhances the clinician's ability to discern a wide array of cardiac conditions, extending diagnostic accuracy beyond myocardial ischemia.

The electrical restitution of the non-propagated cardiac ventricular action potential

Sudden changes in pacing cycle length are frequently associated with repolarization abnormalities initiating cardiac arrhythmias, and physiologists have long been interested in measuring the likelihood of these events before their manifestation. A marker of repolarization stability has been found in the electrical restitution (ER), the response of the ventricular action potential duration to a pre- or post-mature stimulation, graphically represented by the so-called ER curve. According to the restitution hypothesis (ERH), the slope of this curve provides a quantitative discrimination between stable repolarization and proneness to arrhythmias. ER has been studied at the body surface, whole organ, and tissue level, and ERH has soon become a key reference point in theoretical, clinical, and pharmacological studies concerning arrhythmia development, and, despite criticisms, it is still widely adopted. The ionic mechanism of ER and cellular applications of ERH are covered in the present review. The main criticism on ERH concerns its dependence from the way ER is measured. Over the years, in fact, several different experimental protocols have been established to measure ER, which are also described in this article. In reviewing the state-of-the art on cardiac cellular ER, I have introduced a notation specifying protocols and graphical representations, with the aim of unifying a sometime confusing nomenclature, and providing a physiological tool, better defined in its scope and limitations, to meet the growing expectations of clinical and pharmacological research.

Building a search tool for compositely annotated entities using Transformer-based approach: Case study in Biosimulation Model Search Engine (BMSE)

The Transformer-based approaches to solving natural language processing (NLP) tasks such as BERT and GPT are gaining popularity due to their ability to achieve high performance. These approaches benefit from using enormous data sizes to create pre-trained models and the ability to understand the context of words in a sentence. Their use in the information retrieval domain is thought to increase effectiveness and efficiency. This paper demonstrates a BERT-based method (CASBERT) implementation to build a search tool over data annotated compositely using ontologies. The data was a collection of biosimulation models written using the CellML standard in the Physiome Model Repository (PMR). A biosimulation model structurally consists of basic entities of constants and variables that construct higher-level entities such as components, reactions, and the model. Finding these entities specific to their level is beneficial for various purposes regarding variable reuse, experiment setup, and model audit. Initially, we created embeddings representing compositely-annotated entities for constant and variable search (lowest level entity). Then, these low-level entity embeddings were vertically and efficiently combined to create higher-level entity embeddings to search components, models, images, and simulation setups. Our approach was general, so it can be used to create search tools with other data semantically annotated with ontologies - biosimulation models encoded in the SBML format, for example. Our tool is named Biosimulation Model Search Engine (BMSE).

Kv Channel Ancillary Subunits: Where Do We Go from Here?

Voltage-gated potassium (Kv) channels each comprise four pore-forming α-subunits that orchestrate essential duties such as voltage sensing and K+ selectivity and conductance. In vivo, however, Kv channels also incorporate regulatory subunits-some Kv channel specific, others more general modifiers of protein folding, trafficking, and function. Understanding all the above is essential for a complete picture of the role of Kv channels in physiology and disease.

Downregulation of Vascular Hemeoxygenase-1 Leads to Vasculopathy in Systemic Sclerosis

Systemic sclerosis (SSc) is a terminal disease characterized by vasculopathy, tissue fibrosis, and autoimmunity. Although the exact etiology of SSc remains unknown, endothelial dysfunction, oxidative stress, and calcium handling dysregulation have been associated with a large number of SSc-related complications such as neointima formation, vasculogenesis, pulmonary arterial hypertension, impaired angiogenesis, and cardiac arrhythmias. Hemeoxygenase-1 (HO-1) is an antioxidant enzyme involved in multiple biological actions in the cardiovascular system including vascular tone, angiogenesis, cellular proliferation, apoptosis, and oxidative stress. The aim of this work was to investigate the physiological role of HO-1 and its relevance in the cardiovascular complications occurring in SSc. We found that, in early phases of SSc, the expression of HO-1 in dermal fibroblast is lower compared to those isolated from healthy control individuals. This is particularly relevant as reduction of the HO-1/CO signaling pathway is associated with endothelial dysfunction and vasculopathy. We show evidence of the role of HO-1/carbon monoxide (CO) signaling pathway in calcium handling. Using an in vitro model of pulmonary arterial hypertension (PAH) we investigated the role of HO-1 in Ca2+ mobilization from intracellular stores. Our results indicate that HO-1 regulates calcium release from intracellular stores of human pulmonary arterial endothelial cells. We interrogated the activity of HO-1 in angiogenesis using an organotypic co-culture of fibroblast-endothelial cell. Inhibition of HO-1 significantly reduced the ability of endothelial cells to form tubules. We further investigated if this could be associated with cell motility or migration of endothelial cells into the extracellular matrix synthesized by fibroblasts. By mean of holographic imaging, we studied the morphological and functional features of endothelial cells in the presence of an HO-1 activator and selective inhibitors. Our results demonstrate that inhibition of HO-1 significantly reduces cell proliferation and cell motility (migration) of cultured endothelial cells, whilst activation of HO-1 does not modify either morphology, proliferation or motility. In addition, we investigated the actions of CO on the Kv7.1 (KCQN1) channel current, an important component of the cardiac action potential repolarization. Using electrophysiology (whole-cell patch-clamp in a recombinant system overexpressing the KCQN1 channel), we assessed the regulation of KCQN1 by CO. CORM-2, a CO donor, significantly reduced the Kv7.1 current, suggesting that HO-1/CO signaling may play a role in the modulation of the cardiac action potential via regulation of this ion channel. In summary, our results indicate a clear link between: 1) downregulation of HO-1/CO signaling; and 2) pathophysiological processes occurring in early phases of SSc, such as calcium homeostasis dysregulation, impaired angiogenesis and cardiac arrhythmias. A better understanding of the canonical actions (mainly due to the biological actions of CO), and non-canonical actions of HO-1, as well as the interaction of HO-1/CO signaling with other gasotransmitters in SSc will contribute to the development of novel therapeutic approaches.

Mechanisms of gefitinib-induced QT prolongation

Gefitinib, a tyrosine kinase inhibitor, was the first targeted therapy for non-small cell lung cancer (NSCLC). Gefitinib could block human Ether-à-go-go-Related Gene (hERG) channel, an important target in drug-induced long QT syndrome. However, it is unclear whether gefitinib could induce QT interval prolongation. Here, whole-cell patch-clamp technique was used for evaluating the effect of gefitinib on rapidly-activating delayed rectifier K⁺ current (I_Kr), slowly-activating delayed rectifier K⁺ current (I_Ks), transient outward potassium current (I_to), inward rectifier K⁺ current (I_K1) and on action potentials in guinea pig ventricular myocytes. The Langendorff heart perfusion technique was used to determine drug effect on the ECG. Gefitinib depressed I_Kr by binding to open and closed hERG channels in a concentration-dependent way (IC_50: 1.91 μM). The inhibitory effect of gefitinib on wildtype hERG channels was reduced at the hERG mutants Y652A, S636A, F656V and S631A (IC_50: 8.51, 13.97, 18.86, 32.99 μM), indicating that gefitinib is a pore inhibitor of hERG channels. In addition, gefitinib accelerated hERG channel inactivation and decreased channel steady-state inactivation. Gefitinib also decreased I_Ks with IC₅₀ of 23.8 μM. Moreover, gefitinib increased action potential duration (APD) in guinea pig ventricular myocytes and the corrected QT interval (QTc) in isolated perfused guinea pig hearts in a concentration-dependent way (1-30 μM). These findings indicate that gefitinib could prolong QTc interval by potently blocking hERG channel, modulating kinetic properties of hERG channel. Partial block of KCNQ1/KCNE1 could also contribute to delayed repolarization and prolonged QT interval. Thus, caution should be taken when gefitinib is used for NSCLC treatment.

Initiation and entrainment of multicellular automaticity via diffusion limited extracellular domains

Electrically excitable cells often spontaneously and synchronously depolarize in vitro and in vivo preparations. It remains unclear how cells entrain and autorhythmically activate above the intrinsic mean activation frequency of isolated cells with or without pacemaking mechanisms. Recent studies suggest that cyclic ion accumulation and depletion in diffusion-limited extracellular volumes modulate electrophysiology by ephaptic mechanisms (nongap junction or synaptic coupling). This report explores how potassium accumulation and depletion in a restricted extracellular domain induces spontaneous action potentials in two different computational models of excitable cells without gap junctional coupling: Hodgkin-Huxley and Luo-Rudy. Importantly, neither model will spontaneously activate on its own without external stimuli. Simulations demonstrate that cells sharing a diffusion-limited extracellular compartment can become autorhythmic and entrained despite intercellular electrical heterogeneity. Autorhythmic frequency is modulated by the cleft volume and potassium fluxes through the cleft. Additionally, inexcitable cells can suppress or induce autorhythmic activity in an excitable cell via a shared cleft. Diffusion-limited shared clefts can also entrain repolarization. Critically, this model predicts a mechanism by which diffusion-limited shared clefts can initiate, entrain, and modulate multicellular automaticity in the absence of gap junctions.

Not all Long-QTs Are The Same, Proarrhytmic Quantification with Action Potential Triangulation and Alternans

Long-QT is commonly associated with an increased risk of polymorphic ventricular tachycardia from drug therapy. However, not all drugs prolonging QT interval are proarrhythmic. This study aimed to characterize cellular and tissue mechanisms under which QT-interval prolonging drugs and their combination are proarrhythmic, examining arrhythmia susceptibility due to action potential (AP) triangulation and spatial dispersion of action potential duration (APD). Additionally, we aimed to elucidate that Torsades de Pointe (TdP) associated with long-QT are not necessarily caused by early-after-depolarization (EADs) but are related to the presence of AP alternans in both time and space. Isolated Guinea Pig hearts were Langendorff perfused, and optical mapping was done with a voltage dye-sensitive dye. Two commonly used drugs at the beginning of the COVID-19 pandemic, hydroxychloroquine (HCQ) and Azithromycin (AZM), were added to study the effects of QT interval prolongation. Alternans in time and space were characterized by performing restitution pacing protocols. Comparing APs, HCQ prolongs APD during phase-III repolarization, resulting in a higher triangulation ratio than AZM alone or AZM combined with HCQ. Lower triangulation ratios with AZM are associated with phase-II prolongation, lower arrhythmia, and lower incidence of spatially discordant alternans.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Circadian Rhythms of Early Afterdepolarizations and Ventricular Arrhythmias in a Cardiomyocyte Model

Sudden cardiac arrest is a malfunction of the heart's electrical system, typically caused by ventricular arrhythmias, that can lead to sudden cardiac death (SCD) within minutes. Epidemiological studies have shown that SCD and ventricular arrhythmias are more likely to occur in the morning than in the evening, and laboratory studies indicate that these daily rhythms in adverse cardiovascular events are at least partially under the control of the endogenous circadian timekeeping system. However, the biophysical mechanisms linking molecular circadian clocks to cardiac arrhythmogenesis are not fully understood. Recent experiments have shown that L-type calcium channels exhibit circadian rhythms in both expression and function in guinea pig ventricular cardiomyocytes. We developed an electrophysiological model of these cells to simulate the effect of circadian variation in L-type calcium conductance. In our simulations, we found that there is a circadian pattern in the occurrence of early afterdepolarizations (EADs), which are abnormal depolarizations during the repolarization phase of a cardiac action potential that can trigger fatal ventricular arrhythmias. Specifically, the model produces EADs in the morning, but not at other times of day. We show that the model exhibits a codimension-2 Takens-Bogdanov bifurcation that serves as an organizing center for different types of EAD dynamics. We also simulated a two-dimensional spatial version of this model across a circadian cycle. We found that there is a circadian pattern in the breakup of spiral waves, which represents ventricular fibrillation in cardiac tissue. Specifically, the model produces spiral wave breakup in the morning, but not in the evening. Our computational study is the first, to our knowledge, to propose a link between circadian rhythms and EAD formation and suggests that the efficacy of drugs targeting EAD-mediated arrhythmias may depend on the time of day that they are administered.

Mechanisms of Arrhythmias in the Brugada Syndrome

Arrhythmias in Brugada syndrome patients originate in the right ventricular outflow tract (RVOT). Over the past few decades, the characterization of the unique anatomy and electrophysiology of the RVOT has revealed the arrhythmogenic nature of this region. However, the mechanisms that drive arrhythmias in Brugada syndrome patients remain debated as well as the exact site of their occurrence in the RVOT. Identifying the site of origin and mechanism of Brugada syndrome would greatly benefit the development of mechanism-driven treatment strategies.

Structures Illuminate Cardiac Ion Channel Functions in Health and in Long QT Syndrome

The cardiac action potential is critical to the production of a synchronized heartbeat. This electrical impulse is governed by the intricate activity of cardiac ion channels, among them the cardiac voltage-gated potassium (Kv) channels KCNQ1 and hERG as well as the voltage-gated sodium (Nav) channel encoded by SCN5A. Each channel performs a highly distinct function, despite sharing a common topology and structural components. These three channels are also the primary proteins mutated in congenital long QT syndrome (LQTS), a genetic condition that predisposes to cardiac arrhythmia and sudden cardiac death due to impaired repolarization of the action potential and has a particular proclivity for reentrant ventricular arrhythmias. Recent cryo-electron microscopy structures of human KCNQ1 and hERG, along with the rat homolog of SCN5A and other mammalian sodium channels, provide atomic-level insight into the structure and function of these proteins that advance our understanding of their distinct functions in the cardiac action potential, as well as the molecular basis of LQTS. In this review, the gating, regulation, LQTS mechanisms, and pharmacological properties of KCNQ1, hERG, and SCN5A are discussed in light of these recent structural findings.

Short-Long Heart Rate Variation Increases Dispersion of Action Potential Duration in Long QT Type 2 Transgenic Rabbit Model

The initiation of polymorphic ventricular tachycardia in long QT syndrome type 2 (LQT2) has been associated with a characteristic ECG pattern of short-long RR intervals. We hypothesize that this characteristic pattern increases APD dispersion in LQT2, thereby promoting arrhythmia. We investigated APD dispersion and its dependence on two previous cycle lengths (CLs) in transgenic rabbit models of LQT2, LQT1, and their littermate controls (LMC) using random stimulation protocols. The results show that the short-long RR pattern was associated with a larger APD dispersion in LQT2 but not in LQT1 rabbits. The multivariate analyses of APD as a function of two previous CLs (APD_n = C + α₁CL_n−1 + α₂CL_n−2) showed that α₁ (APD restitution slope) is largest and heterogeneous in LQT2 but uniform in LQT1, enhancing APD dispersion under long CL_n−1 in LQT2. The α₂ (short-term memory) was negative in LQT2 while positive in LQT1, and the spatial pattern of α₁ was inversely correlated to α₂ in LQT2, which explains why a short-long combination causes a larger APD dispersion in LQT2 but not in LQT1 rabbits. In conclusion, short-long RR pattern increased APD dispersion only in LQT2 rabbits through heterogeneous APD restitution and the short-term memory, underscoring the genotype-specific triggering of arrhythmias in LQT syndrome.

Human-induced pluripotent stem cell-derived cardiomyocytes have limited IKs for repolarization reserve as revealed by specific KCNQ1/KCNE1 blocker

Objective

We investigated if there is I_Ks, and if there is repolarization reserve by I_Ks in human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs).

Design

We used a specific KCNQ1/KCNE1 channel blocker, L-000768673, with an IC₅₀ of 9 nM, and four hERG-specific blockers, astemizole, cisapride, dofetilide, and E-4031 to investigate the issue.

Results

L-000768673 concentration-dependently prolonged feature point duration (FPD)―a surrogate signal of action potential duration―from 1 to 30 nM without pacing or paced at 1.2 Hz, resulting from I_Ks blockade in hiPSC-CMs. At higher concentrations, the effect of L-000768673 on I_Ks was mitigated by its effect on I_Ca-L, resulting in shortened FPD, reduced impedance amplitude, and increased beating rate at 1 µM and above, recapitulating the self-limiting properties of L-000768673 on action potentials. All four hERG-specific blockers prolonged FPD as expected. Co-application of L-000768673 at sub-threshold (0.1 and 0.3 nM) and threshold (1 nM) concentrations failed to synergistically enhance the effects of hERG blockers on FPD prolongation, rather it showed additive effects, inconsistent with the repolarization reserve role of I_Ks in mature human myocytes that enhanced I_Kr response, implying a difference between hiPSC-CMs used in this study and mature human cardiomyocytes.

Conclusion

There was I_Ks current in hiPSC-CMs, and blockade of I_Ks current caused prolongation of action potential of hiPSC-CMs. However, we could not demonstrate any synergistic effects on action potential duration prolongation of hiPSC-CMs by blocking hERG current and I_Ks current simultaneously, implying little or no repolarization reserve by I_Ks current in hiPSC-CMs used in this study.

Investigating the Complex Arrhythmic Phenotype Caused by the Gain-of-Function Mutation KCNQ1-G229D

The congenital long QT syndrome (LQTS) is a cardiac electrophysiological disorder that can cause sudden cardiac death. LQT1 is a subtype of LQTS caused by mutations in KCNQ1, affecting the slow delayed-rectifier potassium current (I Ks), which is essential for cardiac repolarization. Paradoxically, gain-of-function mutations in KCNQ1 have been reported to cause borderline QT prolongation, atrial fibrillation (AF), sinus bradycardia, and sudden death, however, the mechanisms are not well understood. The goal of the study is to investigate the ionic, cellular and tissue mechanisms underlying the complex phenotype of a gain-of-function mutation in KCNQ1, c.686G > A (p.G229D) using computer modeling and simulations informed by in vitro measurements. Previous studies have shown this mutation to cause AF and borderline QT prolongation. We report a clinical description of a family that carry this mutation and that a member of the family died suddenly during sleep at 21 years old. Using patch-clamp experiments, we confirm that KCNQ1-G229D causes a significant gain in channel function. We introduce the effect of the mutation in populations of atrial, ventricular and sinus node (SN) cell models to investigate mechanisms underlying phenotypic variability. In a population of human atrial and ventricular cell models and tissue, the presence of KCNQ1-G229D predominantly shortens atrial action potential duration (APD). However, in a subset of models, KCNQ1-G229D can act to prolong ventricular APD by up to 7% (19 ms) and underlie depolarization abnormalities, which could promote QT prolongation and conduction delays. Interestingly, APD prolongations were predominantly seen at slow pacing cycle lengths (CL > 1,000 ms), which suggests a greater arrhythmic risk during bradycardia, and is consistent with the observed sudden death during sleep. In a population of human SN cell models, the KCNQ1-G229D mutation results in slow/abnormal sinus rhythm, and we identify that a stronger L-type calcium current enables the SN to be more robust to the mutation. In conclusion, our computational modeling experiments provide novel mechanistic explanations for the observed borderline QT prolongation, and predict that KCNQ1-G229D could underlie SN dysfunction and conduction delays. The mechanisms revealed in the study can potentially inform management and treatment of KCNQ1 gain-of-function mutation carriers.

Electrocardiographic manifestations in severe hypokalemia

Hypokalemia is one of the most common electrolyte disturbances in the clinic and it can increase the risk of life-threatening arrhythmias. Electrocardiographic characteristics associated with hypokalemia include dynamic changes in T-wave morphology, ST-segment depression, and U waves, which are often best seen in the mid-precordial leads (V2–V4). The PR interval can also be prolonged along with an increase in the amplitude of the P wave. We report a case of a patient with hypokalemia (1.31 mmol/L) who showed typical electrocardiographic characteristics of hypokalemia.

Calcium in the Pathophysiology of Atrial Fibrillation and Heart Failure

Atrial fibrillation (AF) is commonly associated with heart failure. A bidirectional relationship exists between the two-AF exacerbates heart failure causing a significant increase in heart failure symptoms, admissions to hospital and cardiovascular death, while pathological remodeling of the atria as a result of heart failure increases the risk of AF. A comprehensive understanding of the pathophysiology of AF is essential if we are to break this vicious circle. In this review, the latest evidence will be presented showing a fundamental role for calcium in both the induction and maintenance of AF. After outlining atrial electrophysiology and calcium handling, the role of calcium-dependent afterdepolarizations and atrial repolarization alternans in triggering AF will be considered. The atrial response to rapid stimulation will be discussed, including the short-term protection from calcium overload in the form of calcium signaling silencing and the eventual progression to diastolic calcium leak causing afterdepolarizations and the development of an electrical substrate that perpetuates AF. The role of calcium in the bidirectional relationship between heart failure and AF will then be covered. The effects of heart failure on atrial calcium handling that promote AF will be reviewed, including effects on both atrial myocytes and the pulmonary veins, before the aspects of AF which exacerbate heart failure are discussed. Finally, the limitations of human and animal studies will be explored allowing contextualization of what are sometimes discordant results.

Transcripts of Kv7.1 and MinK channels and slow delayed rectifier K+ current (IKs) are expressed in zebrafish (Danio rerio) heart

Zebrafish are increasingly used as a model for human cardiac electrophysiology, arrhythmias, and drug screening. However, K⁺ ion channels of the zebrafish heart, which determine the rate of repolarization and duration of cardiac action potential (AP) are still incompletely known and characterized. Here, we provide the first evidence for the presence of the slow component of the delayed rectifier K⁺channels in the zebrafish heart and characterize electrophysiological properties of the slow component of the delayed rectifier K⁺current, I_Ks. Zebrafish atrium and ventricle showed strong transcript expression of the kcnq1 gene, which encodes the Kv7.1 α-subunit of the slow delayed rectifier K⁺ channel. In contrast, the kcne1 gene, encoding the MinK β-subunit of the delayed rectifier, was expressed at 21 and 17 times lower level in ventricle and atrium, respectively, in comparison to the kcnq1. I_Ks was observed in 62% of ventricular myocytes with mean (± SEM) density of 1.23 ± 0.37 pA/pF at + 30 mV. Activation rate of I_Ks was 38% faster (τ₅₀ = 1248 ± 215 ms) than kcnq1:kcne1 channels (1725 ± 792 ms) expressed in 3:1 ratio in Chinese hamster ovary cells. Microelectrode experiments demonstrated the functional relevance of I_Ks in the zebrafish heart, since 100 μM chromanol 293B produced a significant prolongation of AP in zebrafish ventricle. We conclude that AP repolarization in zebrafish ventricle is contributed by I_Ks, which is mainly generated by homotetrameric Kv7.1 channels not coupled to MinK ancillary β-subunits. This is a clear difference to the human heart, where MinK is an essential component of the slow delayed rectifier K⁺channel.

Determinants of Isoform-Specific Gating Kinetics of hERG1 Channel: Combined Experimental and Simulation Study

I_Kr is the rapidly activating component of the delayed rectifier potassium current, the ion current largely responsible for the repolarization of the cardiac action potential. Inherited forms of long QT syndrome (LQTS) (Lees-Miller et al., 1997) in humans are linked to functional modifications in the Kv11.1 (hERG) ion channel and potentially life threatening arrhythmias. There is little doubt now that hERG-related component of I_Kr in the heart depends on the tetrameric (homo- or hetero-) channels formed by two alternatively processed isoforms of hERG, termed hERG1a and hERG1b. Isoform composition (hERG1a- vs. the b-isoform) has recently been reported to alter pharmacologic responses to some hERG blockers and was proposed to be an essential factor pre-disposing patients for drug-induced QT prolongation. Very little is known about the gating and pharmacological properties of two isoforms in heart membranes. For example, how gating mechanisms of the hERG1a channels differ from that of hERG1b is still unknown. The mechanisms by which hERG 1a/1b hetero-tetramers contribute to function in the heart, or what role hERG1b might play in disease are all questions to be answered. Structurally, the two isoforms differ only in the N-terminal region located in the cytoplasm: hERG1b is 340 residues shorter than hERG1a and the initial 36 residues of hERG1b are unique to this isoform. In this study, we combined electrophysiological measurements for HEK cells, kinetics and structural modeling to tease out the individual contributions of each isoform to Action Potential formation and then make predictions about the effects of having various mixture ratios of the two isoforms. By coupling electrophysiological data with computational kinetic modeling, two proposed mechanisms of hERG gating in two homo-tetramers were examined. Sets of data from various experimental stimulation protocols (HEK cells) were analyzed simultaneously and fitted to Markov-chain models (M-models). The minimization procedure presented here, allowed assessment of suitability of different Markov model topologies and the corresponding parameters that describe the channel kinetics. The kinetics modeling pointed to key differences in the gating kinetics that were linked to the full channel structure. Interactions between soluble domains and the transmembrane part of the channel appeared to be critical determinants of the gating kinetics. The structures of the full channel in the open and closed states were compared for the first time using the recent Cryo-EM resolved structure for full open hERG channel and an homology model for the closed state, based on the highly homolog EAG1 channel. Key potential interactions which emphasize the importance of electrostatic interactions between N-PAS cap, S4-S5, and C-linker are suggested based on the structural analysis. The derived kinetic parameters were later used in higher order models of cells and tissue to track down the effect of varying the ratios of hERG1a and hERG1b on cardiac action potentials and computed electrocardiograms. Simulations suggest that the recovery from inactivation of hERG1b may contribute to its physiologic role of this isoform in the action potential. Finally, the results presented here contribute to the growing body of evidence that hERG1b significantly affects the generation of the cardiac I_kr and plays an important role in cardiac electrophysiology. We highlight the importance of carefully revisiting the Markov models previously proposed in order to properly account for the relative abundance of the hERG1 a- and b- isoforms.

Restitution slope is principally determined by steady-state action potential duration

Aims

The steepness of the action potential duration (APD) restitution curve and local tissue refractoriness are both thought to play important roles in arrhythmogenesis. Despite this, there has been little recognition of the apparent association between steady-state APD and the slope of the restitution curve. The objective of this study was to test the hypothesis that restitution slope is determined by APD and to examine the relationship between restitution slope, refractoriness and susceptibility to VF.

Methods and results

Experiments were conducted in isolated hearts and ventricular myocytes from adult guinea pigs and rabbits. Restitution curves were measured under control conditions and following intervention to prolong (clofilium, veratridine, bretylium, low [Ca]e, chronic transverse aortic constriction) or shorten (catecholamines, rapid pacing) ventricular APD. Despite markedly differing mechanisms of action, all interventions that prolonged the action potential led to a steepening of the restitution curve (and vice versa). Normalizing the restitution curve as a % of steady-state APD abolished the difference in restitution curves with all interventions. Effects on restitution were preserved when APD was modulated by current injection in myocytes pre-treated with the calcium chelator BAPTA-AM – to abolish the intracellular calcium transient. The non-linear relation between APD and the rate of repolarization of the action potential is shown to underpin the common influence of APD on the slope of the restitution curve. Susceptibility to VF was found to parallel changes in APD/refractoriness, rather than restitution slope.

Conclusion(s)

Steady-state APD is the principal determinant of the slope of the ventricular electrical restitution curve. In the absence of post-repolarization refractoriness, factors that prolong the action potential would be expected to steepen the restitution curve. However, concomitant changes in tissue refractoriness act to reduce susceptibility to sustained VF. Dependence on steady-state APD may contribute to the failure of restitution slope to predict sudden cardiac death.

A hERG mutation E1039X produced a synergistic lesion on IKs together with KCNQ1-R174C mutation in a LQTS family with three compound mutations

Congenital long QT syndrome (LQTS) caused by compound mutations is usually associated with more severe clinical phenotypes. We identified a LQTS family harboring three compound mutations in different genes (KCNQ1-R174C, hERG-E1039X and SCN5A-E428K). KCNQ1-R174C, hERG-E1039X and SCN5A-E428K mutations and/or relevant wild-type (WT) cDNAs were respectively expressed in mammalian cells. I_Ks-like, I_Kr-like, I_Na-like currents and the functional interaction between KCNQ1-R174C and hERG-E1039X channels were studied using patch-clamp and immunocytochemistry techniques. (1) Expression of KCNQ1-R174C alone showed no I_Ks. Co-expression of KCNQ1-WT + KCNQ1-R174C caused a loss-of-function in I_Ks and blunted the activation of I_Ks in response to isoproterenol. (2) Expression of hERG-E1039X alone and co-expression of hERG-WT + hERG-E1039X negatively shifted inactivation curves and decelerated the recovery time from inactivation. (3) Expression of SCN5A-E428K increased peak I_Na, but had no effect on late I_Na. (4) I_Ks and I_Kr interact, and hERG-E1039X caused a loss-of-function in I_Ks. (5) Immunocytochemical studies indicated that KCNQ1-R174C is trafficking defective and hERG-E1039X is defective in biosynthesis/degradation, but the abnormities were rescued by co-expression with WT. Thus, KCNQ1-R174C and hERG-E1039X disrupted I_Ks and I_Kr functions, respectively. The synergistic lesion, caused by KCNQ1-R174C and hERG-E1039X in I_Ks, is very likely why patients showed more severe phenotypes in the compound mutation case.

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.

D242N, a KV7.1 LQTS mutation uncovers a key residue for IKs voltage dependence

K_V7.1 and KCNE1 co-assemble to give rise to the I_Ks current, one of the most important repolarizing currents of the cardiac action potential. Its relevance is underscored by the identification of >500 mutations in K_V7.1 and, at least, 36 in KCNE1, that cause Long QT Syndrome (LQTS). The aim of this study was to characterize the biophysical and cellular consequences of the D242N K_V7.1 mutation associated with the LQTS. The mutation is located in the S4 transmembrane segment, within the voltage sensor of the K_V7.1 channel, disrupting the conserved charge balance of this region. Perforated patch-clamp experiments show that, unexpectedly, the mutation did not disrupt the voltage-dependent activation but it removed the inactivation and slowed the activation kinetics of D242N K_V7.1 channels. Biotinylation of cell-surface protein and co-immunoprecipitation experiments revealed that neither plasma membrane targeting nor co-assembly between K_V7.1 and KCNE1 was altered by the mutation. However, the association of D242N K_V7.1 with KCNE1 strongly shifted the voltage dependence of activation to more depolarized potentials (+50mV), hindering I_Ks current at physiologically relevant membrane potentials. Both functional and computational analysis suggest that the clinical phenotype of the LQTS patients carrying the D242N mutation is due to impaired action potential adaptation to exercise and, in particular, to increase in heart rate. Moreover, our data identify D242 aminoacidic position as a potential residue involved in the KCNE1-mediated regulation of the voltage dependence of activation of the K_V7.1 channel.

Electrical and Structural Substrate of Arrhythmogenic Right Ventricular Cardiomyopathy Determined Using Noninvasive Electrocardiographic Imaging and Late Gadolinium Magnetic Resonance Imaging

BACKGROUND

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a significant cause of sudden cardiac death in the young. Improved noninvasive assessment of ARVC and better understanding of the disease substrate are important for improving patient outcomes.

METHODS AND RESULTS

We studied 20 genotyped ARVC patients with a broad spectrum of disease using electrocardiographic imaging (a method for noninvasive cardiac electrophysiology mapping) and advanced late gadolinium enhancement cardiac magnetic resonance scar imaging. Compared with 20 healthy controls, ARVC patients had longer ventricular activation duration (median, 52 versus 42 ms; P=0.007) and prolonged mean epicardial activation-recovery intervals (a surrogate for local action potential duration; median, 275 versus 241 ms; P=0.014). In these patients, we observed abnormal and varied epicardial activation breakthrough locations and regions of nonuniform conduction and fractionated electrograms. Nonuniform conduction and fractionated electrograms were present in the early concealed phase of ARVC. Electrophysiological abnormalities colocalized with late gadolinium enhancement scar, indicating a relationship with structural disease. Premature ventricular contractions were common in ARVC patients with variable initiation sites in both ventricles. Premature ventricular contraction rate increased with exercise, and within anatomic segments, it correlated with prolonged repolarization, electric markers of scar, and late gadolinium enhancement (all P<0.001).

CONCLUSIONS

Electrocardiographic imaging reveals electrophysiological substrate properties that differ in ARVC patients compared with healthy controls. A novel mechanistic finding is the presence of repolarization abnormalities in regions where ventricular ectopy originates. The results suggest a potential role for electrocardiographic imaging and late gadolinium enhancement in early diagnosis and noninvasive follow-up of ARVC patients.

Sex and regional differences in rabbit right ventricular L-type calcium current levels and mathematical modelling of arrhythmia vulnerability

NEW FINDINGS

What is the central question of this study? Regional variations of ventricular L-type calcium current (I_Ca-L ) amplitude may underlie the increased arrhythmia risk in adult females. Current amplitude variations have been described for the left ventricle but not for the right ventricle. What is the main finding and its importance? Adult female rabbit right ventricular base myocytes exhibit elevated I_Ca-L compared with female apex or male myocytes. Oestrogen upregulated I_Ca-L in cultured female myocytes. Mathematical simulations modelling long QT syndrome type 2 demonstrated that elevated I_Ca-L prolonged action potentials and induced early after-depolarizations. Thus, ventricular arrhythmias in adult females may be associated with an oestrogen-induced upregulation of I_Ca-L . Previous studies have shown that adult rabbit left ventricular myocytes exhibit sex and regional differences in L-type calcium current (I_Ca-L ) levels that contribute to increased female susceptibility to arrhythmogenic early after-depolarizations (EADs). We used patch-clamp recordings from isolated adult male and female rabbit right ventricular myocytes to determine apex-base differences in I_Ca-L density and used mathematical modelling to examine the contribution of I_Ca-L to EAD formation. Current density measured at 0 mV in female base myocytes was 67% higher than in male base myocytes and 55% higher than in female apex myocytes. No differences were observed between male and female apex myocytes, between male apex and base myocytes, or in the voltage dependences of I_Ca-L activation or inactivation. The role of oestrogen was investigated using cultured adult female right ventricular base myocytes. After 2 days, 17β-estradiol (1 nm) produced a 65% increase in I_Ca-L density compared with untreated control myocytes, suggesting an oestrogen-induced upregulation of I_Ca-L . Action potential simulations using a modified Luo-Rudy cardiomyocyte model showed that increased I_Ca-L density, at the level observed in female base myocytes, resulted in longer duration action potentials, and when combined with a 50% reduction of the rapidly inactivating delayed rectifier potassium current conductance to model long QT syndrome type 2, the action potential was accompanied by one or more EADs. Thus, we found higher levels of I_Ca-L in adult female right ventricle base myocytes and the upregulation of this current by oestrogen. Simulations of long QT syndrome type 2 showed that elevated I_Ca-L contributed to genesis of EADs.

Recent heart rate history affects QT interval duration in atrial fibrillation

QT interval prolongation is associated with a risk of polymorphic ventricular tachycardia. QT interval shortens with increasing heart rate and correction for this effect is necessary for meaningful QT interval assessment. We aim to improve current methods of correcting the QT interval during atrial fibrillation (AF). Digitized Holter recordings were analyzed from patients with AF. Models of QT interval dependence on RR intervals were tested by sorting the beats into 20 bins based on corrected RR interval and assessing ST-T variability within the bins. Signal-averaging within bins was performed to determine QT/RR dependence. Data from 30 patients (29 men, 69.3±7.3 years) were evaluated. QT behavior in AF is well described by a linear function (slope ~0.19) of steady-state corrected RR interval. Corrected RR is calculated as a combination of an exponential weight function with time-constant of 2 minutes and a smaller “immediate response” component (weight ~ 0.18). This model performs significantly (p<0.0001) better than models based on instantaneous RR interval only including Bazett and Fridericia. It also outperforms models based on shorter time-constants and other previously proposed models. This model may improve detection of repolarization delay in AF. QT response to heart rate changes in AF is similar to previously published QT dynamics during atrial pacing and in sinus rhythm.

Murine Electrophysiological Models of Cardiac Arrhythmogenesis

Cardiac arrhythmias can follow disruption of the normal cellular electrophysiological processes underlying excitable activity and their tissue propagation as coherent wavefronts from the primary sinoatrial node pacemaker, through the atria, conducting structures and ventricular myocardium. These physiological events are driven by interacting, voltage-dependent, processes of activation, inactivation, and recovery in the ion channels present in cardiomyocyte membranes. Generation and conduction of these events are further modulated by intracellular Ca2+ homeostasis, and metabolic and structural change. This review describes experimental studies on murine models for known clinical arrhythmic conditions in which these mechanisms were modified by genetic, physiological, or pharmacological manipulation. These exemplars yielded molecular, physiological, and structural phenotypes often directly translatable to their corresponding clinical conditions, which could be investigated at the molecular, cellular, tissue, organ, and whole animal levels. Arrhythmogenesis could be explored during normal pacing activity, regular stimulation, following imposed extra-stimuli, or during progressively incremented steady pacing frequencies. Arrhythmic substrate was identified with temporal and spatial functional heterogeneities predisposing to reentrant excitation phenomena. These could arise from abnormalities in cardiac pacing function, tissue electrical connectivity, and cellular excitation and recovery. Triggering events during or following recovery from action potential excitation could thereby lead to sustained arrhythmia. These surface membrane processes were modified by alterations in cellular Ca2+ homeostasis and energetics, as well as cellular and tissue structural change. Study of murine systems thus offers major insights into both our understanding of normal cardiac activity and its propagation, and their relationship to mechanisms generating clinical arrhythmias.

The Major Role of IK1 in Mechanisms of Rotor Drift in the Atria: A Computational Study

Maintenance of paroxysmal atrial fibrillation (AF) by fast rotors in the left atrium (LA) or at the pulmonary veins (PVs) is not fully understood. This review describes the role of the heterogeneous distribution of transmembrane currents in the PVs and LA junction (PV-LAJ) in the localization of rotors in the PVs. Experimentally observed heterogeneities in IK1, IKs, IKr, Ito, and ICaL in the PV-LAJ were incorporated into models of human atrial kinetics to simulate various conditions and investigate rotor drifting mechanisms. Spatial gradients in the currents resulted in shorter action potential duration, less negative minimum diastolic potential, slower upstroke and conduction velocity for rotors in the PV region than in the LA. Rotors under such conditions drifted toward the PV and stabilized at the less excitable region. Our simulations suggest that IK1 heterogeneity is dominant in determining the drift direction through its impact on the excitability gradient. These results provide a novel framework for understanding the complex dynamics of rotors in AF.

Molecular pathogenesis of long QT syndrome type 1

Long QT syndrome type 1 (LQT1) is a subtype of a congenital cardiac syndrome caused by mutation in the KCNQ1 gene, which encodes the α-subunit of the slow component of delayed rectifier K+ current (IKs) channel. Arrhythmias in LQT1 are characterized by prolongation of the QT interval on ECG, as well as the occurrence of life-threatening cardiac events, frequently triggered by adrenergic stimuli (e.g., physical or emotional stress). During the past two decades, much advancement has been made in understanding the molecular pathogenesis underlying LQT1. Uncovering the genotype-phenotype correlations in LQT1 is of clinical importance to better understand the gene-specific differences that may influence the propensity for developing life-threatening arrhythmias under specific conditions. Elucidation of these mechanisms will also help to improve the diagnosis and management of this cardiac disorder based on gene-specific considerations. This review describes the current medical consensus and recent developments regarding the molecular pathogenesis of LQT1 and provides a novel insight into the adrenergic regulation of this disease.

Assessment of drug-induced proarrhythmia: The importance of study design in the rabbit left ventricular wedge model

In the present study, we investigated an impact of the stimulation rate on the detection of the proarrhythmic potential of 10 reference compounds with effects on different cardiac ion channels in the isolated arterially-perfused rabbit left ventricular wedge preparation. The compounds were tested in the wedge model using two distinct protocols; including baseline stimulation at 1-Hz followed by a brief period at 0.5-Hz, either without an additional brief period of 2-Hz stimulation (i.e. Protocol 1) or with 2-Hz stimulation (i.e. Protocol 2). As expected, QT-prolonging drugs (ibutilide and quinidine) prolonged the QT interval, similarly increased the Torsades de Pointes (TdP) score, and elicited early afterdepolarizations (EADs) in both protocols. HMR1556 and JNJ-303 (IKs blockers) also prolonged the QT interval up to 1μM similarly in both protocols. Nifedipine (Ca(2+) antagonist) shortened the QT interval, and reduced force of contraction similarly in both protocols. However, Na(+) channel blockers (Ia, Ib, Ic) widened the QRS duration more in Protocol 2 than in Protocol 1. Furthermore, it was only possible to detect non-TdP-like ventricular tachycardia/fibrillation (VT/VF) induced by Na(+) blockers and by QT-shortening drugs (levcromakalim and mallotoxin) using the 2-Hz stimulation (Protocol 2). Our data suggest that the inclusion of a brief period of fast stimulation at 2Hz is critical for detecting drug-induced slowing of conduction (QRS widening), QT shortening and associated (non-TdP-like) VT/VF, which are distinct from the QT prolongation/TdP proarrhythmia in isolated, arterially-perfused rabbit left ventricular wedges.

Computational modeling of cardiac fibroblasts and fibrosis

Altered fibroblast behavior can lead to pathologic changes in the heart such as arrhythmia, diastolic dysfunction, and systolic dysfunction. Computational models are increasingly used as a tool to identify potential mechanisms driving a phenotype or potential therapeutic targets against an unwanted phenotype. Here we review how computational models incorporating cardiac fibroblasts have clarified the role for these cells in electrical conduction and tissue remodeling in the heart. Models of fibroblast signaling networks have primarily focused on fibroblast cell lines or fibroblasts from other tissues rather than cardiac fibroblasts, specifically, but they are useful for understanding how fundamental signaling pathways control fibroblast phenotype. In the future, modeling cardiac fibroblast signaling, incorporating -omics and drug-interaction data into signaling network models, and utilizing multi-scale models will improve the ability of in silico studies to predict potential therapeutic targets against adverse cardiac fibroblast activity.

The Contribution of Ionic Currents to Rate-Dependent Action Potential Duration and Pattern of Reentry in a Mathematical Model of Human Atrial Fibrillation

Persistent atrial fibrillation (PeAF) in humans is characterized by shortening of action potential duration (APD) and attenuation of APD rate-adaptation. However, the quantitative influences of particular ionic current alterations on rate-dependent APD changes, and effects on patterns of reentry in atrial tissue, have not been systematically investigated. Using mathematical models of human atrial cells and tissue and performing parameter sensitivity analysis, we evaluated the quantitative contributions to action potential (AP) shortening and APD rate-adaptation of ionic current remodeling seen with PeAF. Ionic remodeling in PeAF was simulated by reducing L-type Ca²⁺ channel current (I_CaL), increasing inward rectifier K⁺ current (I_K1) and modulating five other ionic currents. Parameter sensitivity analysis, which quantified how each ionic current influenced APD in control and PeAF conditions, identified interesting results, including a negative effect of Na⁺/Ca²⁺ exchange on APD only in the PeAF condition. At high pacing rate (2 Hz), electrical remodeling in I_K1 alone accounts for the APD reduction of PeAF, but at slow pacing rate (0.5 Hz) both electrical remodeling in I_CaL alone (-70%) and I_K1 alone (+100%) contribute equally to the APD reduction. Furthermore, AP rate-adaptation was affected by I_Kur in control and by I_NaCa in the PeAF condition. In a 2D tissue model, a large reduction (-70%) of I_CaL becomes a dominant factor leading to a stable spiral wave in PeAF. Our study provides a quantitative and unifying understanding of the roles of ionic current remodeling in determining rate-dependent APD changes at the cellular level and spatial reentry patterns in tissue.

Conformational changes of an ion-channel during gating and emerging electrophysiologic properties: Application of a computational approach to cardiac Kv7.1

Ion channels are the "building blocks" of the excitation process in excitable tissues. Despite advances in determining their molecular structure, understanding the relationship between channel protein structure and electrical excitation remains a challenge. The Kv7.1 potassium channel is an important determinant of the cardiac action potential and its adaptation to rate changes. It is subject to beta adrenergic regulation, and many mutations in the channel protein are associated with the arrhythmic long QT syndrome. In this theoretical study, we use a novel computational approach to simulate the conformational changes that Kv7.1 undergoes during activation gating and compute the resulting electrophysiologic function in terms of single-channel and macroscopic currents. We generated all possible conformations of the S4-S5 linker that couples the S3-S4 complex (voltage sensor domain, VSD) to the pore, and all associated conformations of VSD and the pore (S6). Analysis of these conformations revealed that VSD-to-pore mechanical coupling during activation gating involves outward translation of the voltage sensor, accompanied by a translation away from the pore and clockwise twist. These motions cause pore opening by moving the S4-S5 linker upward and away from the pore, providing space for the S6 tails to move away from each other. Single channel records, computed from the simulated motion trajectories during gating, have stochastic properties similar to experimentally recorded traces. Macroscopic current through an ensemble of channels displays two key properties of Kv7.1: an initial delay of activation and fast inactivation. The simulations suggest a molecular mechanism for fast inactivation; a large twist of the VSD following its outward translation results in movement of the base of the S4-S5 linker toward the pore, eliminating open pore conformations to cause inactivation.

QT Adaptation and Intrinsic QT Variability in Congenital Long QT Syndrome

BACKGROUND

Increased variability of QT interval (QTV) has been linked to arrhythmias in animal experiments and multiple clinical situations. Congenital long QT syndrome (LQTS), a pure repolarization disease, may provide important information on the relationship between delayed repolarization and QTV.

METHODS AND RESULTS

Twenty-four-hour Holter monitor tracings from 78 genotyped congenital LQTS patients (52 females; 51 LQT1, 23 LQT2, 2 LQT5, 2 JLN, 27 symptomatic; age, 35.2±12.3 years) were evaluated with computer-assisted annotation of RR and QT intervals. Several models of RR-QT relationship were tested in all patients. A model assuming exponential decrease of past RR interval contributions to QT duration with 60-second time constant provided the best data fit. This model was used to calculate QTc and residual "intrinsic" QTV, which cannot be explained by heart rate change. The intrinsic QTV was higher in patients with long QTc (r=0.68; P<10(-4)), and in LQT2 than in LQT1/5 patients (5.65±1.28 vs 4.46±0.82; P<0.0002). Both QTc and intrinsic QTV were similar in symptomatic and asymptomatic patients (467±52 vs 459±53 ms and 5.10±1.19 vs 4.74±1.09, respectively).

CONCLUSIONS

In LQTS patients, QT interval adaptation to heart rate changes occurs with time constant ≈60 seconds, similar to results reported in control subjects. Intrinsic QTV correlates with the degree of repolarization delay and might reflect action potential instability observed in animal models of LQTS.

A computational modelling approach combined with cellular electrophysiology data provides insights into the therapeutic benefit of targeting the late Na+ current

KEY POINTS

The ventricular action potential plateau is a phase of high resistance, which makes ventricular myocytes vulnerable to small electrical perturbations. We developed a computationally based model of GS-458967 interaction with the cardiac Na+ channel, informed by experimental data recorded from guinea pig isolated single ventricular myocytes. The model predicts that the therapeutic potential of GS-458967 derives largely from the designed property of significant potent selectivity for INaL.

ABSTRACT

Selective inhibition of the slowly inactivating or late Na(+) current (INaL) in patients with inherited or acquired arrhythmia syndrome may confer therapeutic benefit by reducing the incidence of triggers for arrhythmia and suppressing one component of arrhythmia-promoting cardiac substrates (e.g. prolonged refractoriness and spatiotemporal dispersion of action potential duration). Recently, a novel compound that preferentially and potently reduces INaL, GS-458967 (IC50 for block of INaL = 130 nM) has been studied. Experimental measurements of the effects of GS-458967 on endogenous INaL in guinea pig ventricular myocytes demonstrate a robust concentration-dependent reduction in action potential duration (APD). Using experimental data to calibrate INaL and the rapidly activating delayed rectifier K(+) current, IKr, in the Faber-Rudy computationally based model of the guinea pig ventricular action potential, we simulated effects of GS-458967 on guinea pig ventricular APD. GS-458967 (0.1 μM) caused a 28.67% block of INaL and 12.57% APD shortening in experiments, while the model predicted 10.06% APD shortening with 29.33% block of INaL. An additional effect of INaL block is to reduce the time during which the membrane potential is in a high resistance state (i.e. the action potential plateau). To test the hypothesis that targeted block of INaL would make ventricular myocytes less susceptible to small electrical perturbations, we used the computational model to test the degree of APD prolongation induced by small electrical perturbations in normal cells and in cells with simulated long QT syndrome. The model predicted a substantial dose-dependent reduction in sensitivity to small electrical perturbations as evidenced by action potential duration at 90% repolarization variability in the presence of GS-458967-induced INaL block. This effect was especially potent in the 'disease setting' of inherited long QT syndrome. Using a combined experimental and theoretical approach, our results suggest that INaL block is a potent therapeutic strategy. This is because reduction of INaL stabilizes the action potential waveform by reducing depolarizing current during the plateau phase of the action potential. This reduces the most vulnerable phase of the action potential with high membrane resistance. In summary, by reducing the sensitivity of the myocardial substrate to small electrical perturbations that promote arrhythmia triggers, agents such as GS-458967 may constitute an effective antiarrhythmic pharmacological strategy.

An electromechanical left ventricular wedge model to study the effects of deformation on repolarization during heart failure

Heart failure is a major and costly problem in public health, which, in certain cases, may lead to death. The failing heart undergo a series of electrical and structural changes that provide the underlying basis for disturbances like arrhythmias. Computer models of coupled electrical and mechanical activities of the heart can be used to advance our understanding of the complex feedback mechanisms involved. In this context, there is a lack of studies that consider heart failure remodeling using strongly coupled electromechanics. We present a strongly coupled electromechanical model to study the effects of deformation on a human left ventricle wedge considering normal and hypertrophic heart failure conditions. We demonstrate through a series of simulations that when a strongly coupled electromechanical model is used, deformation results in the thickening of the ventricular wall that in turn increases transmural dispersion of repolarization. These effects were analyzed in both normal and failing heart conditions. We also present transmural electrograms obtained from these simulations. Our results suggest that the waveform of electrograms, particularly the T-wave, is influenced by cardiac contraction on both normal and pathological conditions.

A new KCNQ1 mutation at the S5 segment that impairs its association with KCNE1 is responsible for short QT syndrome

AIMS

KCNQ1 and KCNE1 encode Kv7.1 and KCNE1, respectively, the pore-forming and the accessory subunits of the slow delayed rectifier potassium current, IKs. KCNQ1 mutations are associated with long and short QT syndrome. The aim of this study was to characterize the biophysical and cellular phenotype of a KCNQ1 missense mutation, F279I, found in a 23-year-old man with a corrected QT interval (QTc) of 356 ms and a family history of sudden cardiac death.

METHODS AND RESULTS

Experiments were performed using perforated patch-clamp, western blot, co-immunoprecipitation, biotinylation, and immunocytochemistry techniques in HEK293, COS7 cells and in cardiomyocytes transfected with WT Kv7.1/KCNE1 or F279I Kv7.1/KCNE1 channels. In the absence of KCNE1, F279I Kv7.1 current exhibited a lesser degree of inactivation than WT Kv7.1. Also, functional analysis of F279I Kv7.1 in the presence of KCNE1 revealed a negative shift in the activation curve and an acceleration of the activation kinetics leading to a gain of function in IKs. The co-assembly between F279I Kv7.1 channels and KCNE1 was markedly decreased compared with WT Kv7.1 channels, as revealed by co-immunoprecipitation and Föster Resonance Energy Transfer experiments. All these effects contribute to the increase of IKs when channels incorporate F279I Kv7.1 subunits, as shown by a computer model simulation of these data that predicts a shortening of the action potential (AP) consistent with the patient phenotype.

CONCLUSION

The F279I mutation induces a gain of function of IKs due to an impaired gating modulation of Kv7.1 induced by KCNE1, leading to a shortening of the cardiac AP.

Spatially Discordant Alternans and Arrhythmias in Tachypacing-Induced Cardiac Myopathy in Transgenic LQT1 Rabbits: The Importance of IKs and Ca2+ Cycling

BACKGROUND

Remodeling of cardiac repolarizing currents, such as the downregulation of slowly activating K+ channels (IKs), could underlie ventricular fibrillation (VF) in heart failure (HF). We evaluated the role of Iks remodeling in VF susceptibility using a tachypacing HF model of transgenic rabbits with Long QT Type 1 (LQT1) syndrome.

METHODS AND RESULTS

LQT1 and littermate control (LMC) rabbits underwent three weeks of tachypacing to induce cardiac myopathy (TICM). In vivo telemetry demonstrated steepening of the QT/RR slope in LQT1 with TICM (LQT1-TICM; pre: 0.26±0.04, post: 0.52±0.01, P<0.05). In vivo electrophysiology showed that LQT1-TICM had higher incidence of VF than LMC-TICM (6 of 11 vs. 3 of 11, respectively). Optical mapping revealed larger APD dispersion (16±4 vs. 38±6 ms, p<0.05) and steep APD restitution in LQT1-TICM compared to LQT1-sham (0.53±0.12 vs. 1.17±0.13, p<0.05). LQT1-TICM developed spatially discordant alternans (DA), which caused conduction block and higher-frequency VF (15±1 Hz in LQT1-TICM vs. 13±1 Hz in LMC-TICM, p<0.05). Ca2+ DA was highly dynamic and preceded voltage DA in LQT1-TICM. Ryanodine abolished DA in 5 out of 8 LQT1-TICM rabbits, demonstrating the importance of Ca2+ in complex DA formation. Computer simulations suggested that HF remodeling caused Ca2+-driven alternans, which was further potentiated in LQT1-TICM due to the lack of IKs.

CONCLUSIONS

Compared with LMC-TICM, LQT1-TICM rabbits exhibit steepened APD restitution and complex DA modulated by Ca2+. Our results strongly support the contention that the downregulation of IKs in HF increases Ca2+ dependent alternans and thereby the risk of VF.

Attraction of rotors to the pulmonary veins in paroxysmal atrial fibrillation: a modeling study

Maintenance of paroxysmal atrial fibrillation (AF) by fast rotors in the left atrium (LA) or at the pulmonary veins (PVs) is not fully understood. To gain insight into this dynamic and complex process, we studied the role of the heterogeneous distribution of transmembrane currents in the PVs and LA junction (PV-LAJ) in the localization of rotors in the PVs. We also investigated whether simple pacing protocols could be used to predict rotor drift in the PV-LAJ. Experimentally observed heterogeneities in IK1, IKs, IKr, Ito, and ICaL in the PV-LAJ were incorporated into two- and pseudo three-dimensional models of Courtemanche-Ramirez-Nattel-Kneller human atrial kinetics to simulate various conditions and investigate rotor drifting mechanisms. Spatial gradients in the currents resulted in shorter action potential duration, minimum diastolic potential that was less negative, and slower upstroke and conduction velocity for rotors in the PV region than in the LA. Rotors under such conditions drifted toward the PV and stabilized at the shortest action potential duration and less-excitable region, consistent with drift direction under intercellular coupling heterogeneities and regardless of the geometrical constraint in the PVs. Simulations with various IK1 gradient conditions and current-voltage relationships substantiated its major role in the rotor drift. In our 1:1 pacing protocol, we found that among various action potential properties, only the minimum diastolic potential gradient was a rate-independent predictor of rotor drift direction. Consistent with experimental and clinical AF studies, simulations in an electrophysiologically heterogeneous model of the PV-LAJ showed rotor attraction toward the PV. Our simulations suggest that IK1 heterogeneity is dominant compared to other currents in determining the drift direction through its impact on the excitability gradient. These results provide a believed novel framework for understanding the complex dynamics of rotors in AF.

Three-dimensional computer model of the right atrium including the sinoatrial and atrioventricular nodes predicts classical nodal behaviours

The aim of the study was to develop a three-dimensional (3D) anatomically-detailed model of the rabbit right atrium containing the sinoatrial and atrioventricular nodes to study the electrophysiology of the nodes. A model was generated based on 3D images of a rabbit heart (atria and part of ventricles), obtained using high-resolution magnetic resonance imaging. Segmentation was carried out semi-manually. A 3D right atrium array model (∼3.16 million elements), including eighteen objects, was constructed. For description of cellular electrophysiology, the Rogers-modified FitzHugh-Nagumo model was further modified to allow control of the major characteristics of the action potential with relatively low computational resource requirements. Model parameters were chosen to simulate the action potentials in the sinoatrial node, atrial muscle, inferior nodal extension and penetrating bundle. The block zone was simulated as passive tissue. The sinoatrial node, crista terminalis, main branch and roof bundle were considered as anisotropic. We have simulated normal and abnormal electrophysiology of the two nodes. In accordance with experimental findings: (i) during sinus rhythm, conduction occurs down the interatrial septum and into the atrioventricular node via the fast pathway (conduction down the crista terminalis and into the atrioventricular node via the slow pathway is slower); (ii) during atrial fibrillation, the sinoatrial node is protected from overdrive by its long refractory period; and (iii) during atrial fibrillation, the atrioventricular node reduces the frequency of action potentials reaching the ventricles. The model is able to simulate ventricular echo beats. In summary, a 3D anatomical model of the right atrium containing the cardiac conduction system is able to simulate a wide range of classical nodal behaviours.

Exponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes

The modern Markov chain models of ionic channels in excitable membranes are numerically stiff. The popular numerical methods for these models require very small time steps to ensure stability. Our objective is to formulate and test two methods addressing this issue, so that the timestep can be chosen based on accuracy rather than stability. Both proposed methods extend Rush-Larsen technique, which was originally developed to Hogdkin-Huxley type gate models. One method, "matrix Rush-Larsen" (MRL) uses a matrix reformulation of the Rush-Larsen scheme, where the matrix exponentials are calculated using precomputed tables of eigenvalues and eigenvectors. The other, "hybrid operator splitting" (HOS) method exploits asymptotic properties of a particular Markov chain model, allowing explicit analytical expressions for the substeps. We test both methods on the Clancy and Rudy (2002) I(Na)Markov chain model. With precomputed tables for functions of the transmembrane voltage, both methods are comparable to the forward Euler method in accuracy and computational cost, but allow longer time steps without numerical instability. We conclude that both methods are of practical interest. MRL requires more computations than HOS, but is formulated in general terms which can be readily extended to other Markov chain channel models, whereas the utility of HOS depends on the asymptotic properties of a particular model. The significance of the methods is that they allow a considerable speed-up of large-scale computations of cardiac excitation models by increasing the time step, while maintaining acceptable accuracy and preserving numerical stability.

Application of stochastic phenomenological modelling to cell-to-cell and beat-to-beat electrophysiological variability in cardiac tissue

Variability in the action potential of isolated myocytes and tissue samples is observed in experimental studies. Variability is manifested as both differences in the action potential (AP) morphology between cells (extrinsic variability), and also ‘intrinsic’ or beat-to-beat variability of repolarization (BVR) in the AP duration of each cell. We studied the relative contributions of experimentally recorded intrinsic and extrinsic variability to dispersion of repolarization in tissue. We developed four cell-specific parameterizations of a phenomenological stochastic differential equation AP model exhibiting intrinsic variability using APs recorded from isolated guinea pig ventricular myocytes exhibiting BVR. We performed simulations in tissue using the four different model parameterizations in the presence and the absence of both intrinsic and extrinsic variability. We altered the coupling of the tissue to determine how inter-cellular coupling affected the dispersion of the AP duration in tissue. Both intrinsic and extrinsic variability were gradually revealed by reduction of tissue coupling. However, the recorded extrinsic variability between individual myocytes produced a greater degree of dispersion in repolarization in tissue than the intrinsic variability of each myocyte.

•

We modelled inter-cell and beat-to-beat repolarization variability in cardiomyocytes.

•

We coupled the cells together into cardiac tissue.

•

Reducing tissue coupling increased repolarization dispersion in tissue.

•

Inter-cell variability had a greater effect on repolarization dispersion.

In silico ischaemia-induced reentry at the Purkinje-ventricle interface

AIMS

This computational modelling work illustrates the influence of hyperkalaemia and electrical uncoupling induced by defined ischaemia on action potential (AP) propagation and the incidence of reentry at the Purkinje-ventricle interface in mammalian hearts.

METHODS AND RESULTS

Unidimensional and bidimensional models of the Purkinje-ventricle subsystem, including ischaemic conditions (defined as phase 1B) in the ventricle and an ischaemic border zone, were developed by altering several important electrophysiological parameters of the Luo-Rudy AP model of the ventricular myocyte. Purkinje electrical activity was modelled using the equations of DiFrancesco and Noble. Our study suggests that an extracellular potassium concentration [K(+)]o >14 mM and a slight decrease in intercellular coupling induced by ischaemia in ventricle can cause conduction block from Purkinje to ventricle. Under these conditions, propagation from ventricle to Purkinje is possible. Thus, unidirectional block (UDB) and reentry can result. When conditions of UDB are met, retrograde propagation with a long delay (320 ms) may re-excite Purkinje cells, and give rise to a reentrant pathway. This induced reentry may be the origin of arrhythmias observed in phase 1B ischaemia.

CONCLUSION

In a defined setting of ischaemia (phase 1B), a small amount of uncoupling between ventricular cells, as well as between Purkinje and ventricular tissue, may induce UDBs and reentry. Hyperkalaemia is also confirmed to be an important factor in the genesis of reentrant rhythms, since it regulates the range of coupling in which UDBs may be induced.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Electrophysiologic substrate in congenital Long QT syndrome: noninvasive mapping with electrocardiographic imaging (ECGI)

BACKGROUND

Congenital Long QT syndrome (LQTS) is an arrhythmogenic disorder that causes syncope and sudden death. Although its genetic basis has become well-understood, the mechanisms whereby mutations translate to arrhythmia susceptibility in the in situ human heart have not been fully defined. We used noninvasive ECG imaging to map the cardiac electrophysiological substrate and examine whether LQTS patients display regional heterogeneities in repolarization, a substrate that promotes arrhythmogenesis.

METHODS AND RESULTS

Twenty-five subjects (9 LQT1, 9 LQT2, 5 LQT3, and 2 LQT5) with genotype and phenotype positive LQTS underwent ECG imaging. Seven normal subjects provided control. Epicardial maps of activation, recovery times, activation-recovery intervals, and repolarization dispersion were constructed. Activation was normal in all patients. However, recovery times and activation-recovery intervals were prolonged relative to control, indicating delayed repolarization and abnormally long action potential duration (312±30 ms versus 235±21 ms in control). Activation-recovery interval prolongation was spatially heterogeneous, with repolarization gradients much steeper than control (119±19 ms/cm versus 2.0±2.0 ms/cm). There was variability in steepness and distribution of repolarization gradients between and within LQTS types. Repolarization gradients were steeper in symptomatic patients (130±27 ms/cm in 12 symptomatic patients versus 98±19 ms/cm in 13 asymptomatic patients; P<0.05).

CONCLUSIONS

LQTS patients display regions with steep repolarization dispersion caused by localized action potential duration prolongation. This defines a substrate for reentrant arrhythmias, not detectable by surface ECG. Steeper dispersion in symptomatic patients suggests a possible role for ECG imaging in risk stratification.

Differential modulations of KCNQ1 by auxiliary proteins KCNE1 and KCNE2

KCNQ1 channels play vital roles in cardiovascular, gastric and other systems. The conductance and dynamics of KCNQ1 could be modulated by different single transmembrane helical auxiliary proteins (such as KCNE1, KCNE2 and others). In this study, detail KCNQ1 function modulations by different regions of KCNE1 or KCNE2 were examined using combinational methods of electrophysiology, immunofluorescence, solution NMR and related backbone flexibility analysis. In the presence of KCNE2 N-terminus, decreased surface expression and consequent low activities of KCNQ1 were observed. The transmembrane domains (TMDs) of KCNE1 and KCNE2 were illustrated to associate with the KCNQ1 channel in different modes: Ile64 in KCNE2-TMD interacting with Phe340 and Phe275 in KCNQ1, while two pairs of interacting residues (Phe340-Thr58 and Ala244-Tyr65) in the KCNQ1/KCNE1 complex. The KCNE1 C-terminus could modulate gating property of KCNQ1, whereas KCNE2 C-terminus had only minimal influences on KCNQ1. All of the results demonstrated different KCNQ1 function modulations by different regions of the two auxiliary proteins.