Representing variability and transmural differences in a model of human heart failure

Facilitation of hERG Activation by Its Blocker: A Mechanism to Reduce Drug-Induced Proarrhythmic Risk

Modulation of the human Ether-à-go-go-Related Gene (hERG) channel, a crucial voltage-gated potassium channel in the repolarization of action potentials in ventricular myocytes of the heart, has significant implications on cardiac electrophysiology and can be either antiarrhythmic or proarrhythmic. For example, hERG channel blockade is a leading cause of long QT syndrome and potentially life-threatening arrhythmias, such as torsades de pointes. Conversely, hERG channel blockade is the mechanism of action of Class III antiarrhythmic agents in terminating ventricular tachycardia and fibrillation. In recent years, it has been recognized that less proarrhythmic hERG blockers with clinical potential or Class III antiarrhythmic agents exhibit, in addition to their hERG-blocking activity, a second action that facilitates the voltage-dependent activation of the hERG channel. This facilitation is believed to reduce the proarrhythmic potential by supporting the final repolarizing of action potentials. This review covers the pharmacological characteristics of hERG blockers/facilitators, the molecular mechanisms underlying facilitation, and their clinical significance, as well as unresolved issues and requirements for research in the fields of ion channel pharmacology and drug-induced arrhythmias.

In silico assessment of pharmacotherapy for carbon monoxide induced arrhythmias in healthy and failing human hearts

Background: Carbon monoxide (CO) is gaining increased attention in air pollution-induced arrhythmias. The severe cardiotoxic consequences of CO urgently require effective pharmacotherapy to treat it. However, existing evidence demonstrates that CO can induce arrhythmias by directly affecting multiple ion channels, which is a pathway distinct from heart ischemia and has received less concern in clinical treatment.

Objective: To evaluate the efficacy of some common clinical antiarrhythmic drugs for CO-induced arrhythmias, and to propose a potential pharmacotherapy for CO-induced arrhythmias through the virtual pathological cell and tissue models.

Methods: Two pathological models describing CO effects on healthy and failing hearts were constructed as control baseline models. After this, we first assessed the efficacy of some common antiarrhythmic drugs like ranolazine, amiodarone, nifedipine, etc., by incorporating their ion channel-level effects into the cell model. Cellular biomarkers like action potential duration and tissue-level biomarkers such as the QT interval from pseudo-ECGs were obtained to assess the drug efficacy. In addition, we also evaluated multiple specific IKr activators in a similar way to multi-channel blocking drugs, as the IKr activator showed great potency in dealing with CO-induced pathological changes.

Results: Simulation results showed that the tested seven antiarrhythmic drugs failed to rescue the heart from CO-induced arrhythmias in terms of the action potential and the ECG manifestation. Some of them even worsened the condition of arrhythmogenesis. In contrast, IKr activators like HW-0168 effectively alleviated the proarrhythmic effects of CO.

Conclusion: Current antiarrhythmic drugs including the ranolazine suggested in previous studies did not achieve therapeutic effects for the cardiotoxicity of CO, and we showed that the specific IKr activator is a promising pharmacotherapy for the treatment of CO-induced arrhythmias.

Deterministic and Stochastic Cellular Mechanisms Contributing to Carbon Monoxide Induced Ventricular Arrhythmias

Chronic exposure to low levels of Carbon Monoxide is associated with an increased risk of cardiac arrhythmia. Microelectrode recordings from rat and guinea pig single isolated ventricular myocytes exposed to CO releasing molecule CORM-2 and excited at 0.2/s show repolarisation changes that develop over hundreds of seconds: action potential prolongation by delayed repolarisation, EADs, multiple EADs and oscillations around the plateau, leading to irreversible repolarisation failure. The measured direct effects of CO on currents in these cells, and ion channels expressed in mammalian systems showed an increase in prolonged late Na+, and a decrease in the maximal T- and L-type Ca++. peak and late Na+, ultra-rapid delayed, delayed rectifier, and the inward rectifier K+ currents. Incorporation of these CO induced changes in maximal currents in ventricular cell models; (Gattoni et al., J. Physiol., 2016, 594, 4193-4224) (rat) and (Luo and Rudy, Circ. Res., 1994, 74, 1071-1096) (guinea-pig) and human endo-, mid-myo- and epi-cardial (O'Hara et al., PLoS Comput. Biol., 2011, 7, e1002061) models, by changes in maximal ionic conductance reproduces these repolarisation abnormalities. Simulations of cell populations with Gaussian distributions of maximal conductance parameters predict a CO induced increase in APD and its variability. Incorporation of these predicted CO induced conductance changes in human ventricular cell electrophysiology into ventricular tissue and wall models give changes in indices for the probability of the initiation of re-entrant arrhythmia.

Nanoscale regulation of L-type calcium channels differentiates between ischemic and dilated cardiomyopathies

BACKGROUND

Subcellular localization and function of L-type calcium channels (LTCCs) play an important role in regulating contraction of cardiomyocytes. Understanding how this is affected by the disruption of transverse tubules during heart failure could lead to new insights into the disease.

METHODS

Cardiomyocytes were isolated from healthy donor hearts, as well as from patients with cardiomyopathies and with left ventricular assist devices. Scanning ion conductance and confocal microscopy was used to study membrane structures in the cells. Super-resolution scanning patch-clamp was used to examine LTCC function in different microdomains. Computational modeling predicted the impact of these changes to arrhythmogenesis at the whole-heart level.

FINDINGS

We showed that loss of structural organization in failing myocytes leads to re-distribution of functional LTCCs from the T-tubules to the sarcolemma. In ischemic cardiomyopathy, the increased LTCC open probability in the T-tubules depends on the phosphorylation by protein kinase A, whereas in dilated cardiomyopathy, the increased LTCC opening probability in the sarcolemma results from enhanced phosphorylation by calcium-calmodulin kinase II. LVAD implantation corrected LTCCs pathophysiological activity, although it did not improve their distribution. Using computational modeling in a 3D anatomically-realistic human ventricular model, we showed how LTCC location and activity can trigger heart rhythm disorders of different severity.

INTERPRETATION

Our findings demonstrate that LTCC redistribution and function differentiate between disease aetiologies. The subcellular changes observed in specific microdomains could be the consequence of the action of distinct protein kinases.

FUNDING

This work was supported by NIH grant (ROI-HL 126802 to NT-JG) and British Heart Foundation (grant RG/17/13/33173 to JG, project grant PG/16/17/32069 to RAC). Funders had no role in study design, data collection, data analysis, interpretation, writing of the report.

Anisotropic shortening in the wavelength of electrical waves promotes onset of electrical turbulence in cardiac tissue: An in silico study

Several pathological conditions introduce spatial variations in the electrical properties of cardiac tissue. These variations occur as localized or distributed gradients in ion-channel functionality over extended tissue media. Electrical waves, propagating through such affected tissue, demonstrate distortions, depending on the nature of the ionic gradient in the diseased substrate. If the degree of distortion is large, reentrant activity may develop, in the form of rotating spiral (2d) and scroll (3d) waves of electrical activity. These reentrant waves are associated with the occurrence of lethal cardiac rhythm disorders, known as arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), which are believed to be common precursors of sudden cardiac arrest. By using state-of-the-art mathematical models for generic, and ionically-realistic (human) cardiac tissue, we study the detrimental effects of these ionic gradients on electrical wave propagation. We propose a possible mechanism for the development of instabilities in reentrant wave patterns, in the presence of ionic gradients in cardiac tissue, which may explain how one type of arrhythmia (VT) can degenerate into another (VF). Our proposed mechanism entails anisotropic reduction in the wavelength of the excitation waves because of anisotropic variation in its electrical properties, in particular the action potential duration (APD). We find that the variation in the APD, which we induce by varying ion-channel conductances, imposes a spatial variation in the spiral- or scroll-wave frequency ω. Such gradients in ω induce anisotropic shortening of wavelength of the spiral or scroll arms and eventually leads to instabilitites.

Facilitation of I Kr current by some hERG channel blockers suppresses early afterdepolarizations

Some hERG channel blockers are clinically safe, but others cause fatal cardiac arrhythmias. Furutani et al. show that safe blockers facilitate channel opening in ventricular myocytes and provide a repolarization reserve at precisely the voltages and times needed to suppress arrhythmias.

Drug-induced block of the cardiac rapid delayed rectifying potassium current (I_Kr), carried by the human ether-a-go-go-related gene (hERG) channel, is the most common cause of acquired long QT syndrome. Indeed, some, but not all, drugs that block hERG channels cause fatal cardiac arrhythmias. However, there is no clear method to distinguish between drugs that cause deadly arrhythmias and those that are clinically safe. Here we propose a mechanism that could explain why certain clinically used hERG blockers are less proarrhythmic than others. We demonstrate that several drugs that block hERG channels, but have favorable cardiac safety profiles, also evoke another effect; they facilitate the hERG current amplitude in response to low-voltage depolarization. To investigate how hERG facilitation impacts cardiac safety, we develop computational models of I_Kr block with and without this facilitation. We constrain the models using data from voltage clamp recordings of hERG block and facilitation by nifekalant, a safe class III antiarrhythmic agent. Human ventricular action potential simulations demonstrate the ability of nifekalant to suppress ectopic excitations, with or without facilitation. Without facilitation, excessive I_Kr block evokes early afterdepolarizations, which cause lethal arrhythmias. When facilitation is introduced, early afterdepolarizations are prevented at the same degree of block. Facilitation appears to prevent early afterdepolarizations by increasing I_Kr during the repolarization phase of action potentials. We empirically test this prediction in isolated rabbit ventricular myocytes and find that action potential prolongation with nifekalant is less likely to induce early afterdepolarization than action potential prolongation with dofetilide, a hERG channel blocker that does not induce facilitation. Our data suggest that hERG channel blockers that induce facilitation increase the repolarization reserve of cardiac myocytes, rendering them less likely to trigger lethal ventricular arrhythmias.

Degradation of T-Tubular Microdomains and Altered cAMP Compartmentation Lead to Emergence of Arrhythmogenic Triggers in Heart Failure Myocytes: An in silico Study

Heart failure (HF) is one of the most common causes of morbidity and mortality worldwide. Although many patients suffering from HF die from sudden cardiac death caused by arrhythmias, the mechanism linking HF remodeling to an increased arrhythmogenic propensity remains incomplete. HF is typically characterized by a progressive loss of transverse tubule (T-tubule) domains, which leads to an altered distribution of L-type calcium channels (LTCCs). Microdomain degradation also causes the disruption of the β₂ adrenergic receptor (β₂AR) and phosphodiesterase (PDE) signaling localization, normally confined to the dyadic space. The goal of this study was to analyze how these subcellular changes affect the function of LTCCs and lead to the emergence of ventricular cell-level triggers of arrhythmias. To accomplish this, we developed a novel computational model of a human ventricular HF myocyte in which LTCCs were divided into six different populations, based on their location and signaling environment they experience. To do so, we included T-tubular microdomain remodeling which led to a subset of LTCCs to be redistributed from the T-tubular to the surface membrane and allowed for different levels of phosphorylation of LTCCs by PKA, based on the presence of β₂ARs and PDEs. The model was used to study the behavior of the LTCC current (I_CaL) under basal and sympathetic stimulation and its effect on cellular action potential. Our results showed that channels redistributed from the T-tubular membrane to the bulk of the sarcolemma displayed an altered function in their new, non-native signaling domain. Incomplete calcium dependent inactivation, which resulted in a longer-lasting and larger-in-magnitude LTCC current, was observed when we decoupled LTCCs from ryanodine receptors and removed them from the dyadic space. The magnitude of the LTCC current, especially in the surface sarcolemma, was also increased via phosphorylation by the redistributed β₂ARs and PDEs. These changes in LTCC current led to the development of early afterdepolarizations. Thus, our study shows that altered LTCC function is a potential cause for the emergence of cell-level triggers of arrhythmia, and that β₂ARs and PDEs present useful therapeutic targets for treatment of HF and prevention of sudden cardiac death.

Sensitivity analysis revealing the effect of modulating ionic mechanisms on calcium dynamics in simulated human heart failure

Abnormal intracellular Ca²⁺ handling is the major contributor to the depressed cardiac contractility observed in heart failure. The electrophysiological remodeling associated with this pathology alters both the action potential and the Ca²⁺ dynamics, leading to a defective excitation-contraction coupling that ends in mechanical dysfunction. The importance of maintaining a correct intracellular Ca²⁺ concentration requires a better understanding of its regulation by ionic mechanisms. To study the electrical activity and ionic homeostasis of failing myocytes, a modified version of the O’Hara et al. human action potential model was used, including electrophysiological remodeling. The impact of the main ionic transport mechanisms was analyzed using single-parameter sensitivity analyses, the first of which explored the modulation of electrophysiological characteristics related to Ca²⁺ exerted by the remodeled parameters. The second sensitivity analysis compared the potential consequences of modulating individual channel conductivities, as one of the main effects of potential drugs, on Ca²⁺ dynamic properties under both normal conditions and in heart failure. The first analysis revealed the important contribution of the sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) dysfunction to the altered Ca²⁺ homeostasis, with the Na⁺/Ca²⁺ exchanger (NCX) and other Ca²⁺ cycling proteins also playing a significant role. Our results highlight the importance of improving the SR uptake function to increase Ca²⁺ content and restore Ca²⁺ homeostasis and contractility. The second sensitivity analysis highlights the different response of the failing myocyte versus the healthy myocyte to potential pharmacological actions on single channels. The result of modifying the conductances of the remodeled proteins such as SERCA and NCX in heart failure has less impact on Ca²⁺ modulation. These differences should be taken into account when designing drug therapies.

Real-Time ECG Simulation for Hybrid Mock Circulatory Loops

Classically, mock circulatory loops only simulate mechanical properties of the circulation. To connect the hydraulic world with electrophysiology, we present a real-time electrical activity model of the heart and show how to integrate this model into a real-time mock loop simulation. The model incorporates a predefined conduction pathway and a simplified volume conductor to solve the bidomain equations and the forward problem of electrocardiography, resulting in a physiological simulation of the electrocardiogram (ECG) at arbitrary electrode positions. A complete physiological simulation of the heart's excitation would be too CPU intensive. Thus, in our model, complexity was reduced to allow real-time simulation of ECG-triggered medical systems in vitro; this decreases time and cost in the development process. Conversely, the presented model can still be adapted to various pathologies by locally changing the properties of the heart's conduction pathway. To simulate the ECG, the heart is divided into suitable areas, which are innervated by the hierarchically structured conduction system. To distinguish different cardiac regions, a segmentation of the heart was performed. In these regions, Prim's algorithm was applied to identify the directed minimal spanning trees for conduction orientation. Each node of the tree was assigned to a cardiac action potential generated by its hybrid automaton to represent the heart's conduction system by the spatial distribution of action potentials. To generate the ECG output, the bidomain equations were implemented and a simple model of the volume conductor of the body was used to solve the forward problem of electrocardiography. As a result, the model simulates potentials at arbitrary electrode positions in real-time. To verify the developed real-time ECG model, measurements were made within a hybrid mock circulatory loop, including a simple ECG-triggered ventricular assist device control. The model's potential value is to simulate physiological and pathological behavior for hardware-in-the-loop testing of medical devices in an ECG-triggered scenario.

Human In Silico Drug Trials Demonstrate Higher Accuracy than Animal Models in Predicting Clinical Pro-Arrhythmic Cardiotoxicity

Early prediction of cardiotoxicity is critical for drug development. Current animal models raise ethical and translational questions, and have limited accuracy in clinical risk prediction. Human-based computer models constitute a fast, cheap and potentially effective alternative to experimental assays, also facilitating translation to human. Key challenges include consideration of inter-cellular variability in drug responses and integration of computational and experimental methods in safety pharmacology. Our aim is to evaluate the ability of in silico drug trials in populations of human action potential (AP) models to predict clinical risk of drug-induced arrhythmias based on ion channel information, and to compare simulation results against experimental assays commonly used for drug testing. A control population of 1,213 human ventricular AP models in agreement with experimental recordings was constructed. In silico drug trials were performed for 62 reference compounds at multiple concentrations, using pore-block drug models (IC₅₀/Hill coefficient). Drug-induced changes in AP biomarkers were quantified, together with occurrence of repolarization/depolarization abnormalities. Simulation results were used to predict clinical risk based on reports of Torsade de Pointes arrhythmias, and further evaluated in a subset of compounds through comparison with electrocardiograms from rabbit wedge preparations and Ca²⁺-transient recordings in human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs). Drug-induced changes in silico vary in magnitude depending on the specific ionic profile of each model in the population, thus allowing to identify cell sub-populations at higher risk of developing abnormal AP phenotypes. Models with low repolarization reserve (increased Ca²⁺/late Na⁺ currents and Na⁺/Ca²⁺-exchanger, reduced Na⁺/K⁺-pump) are highly vulnerable to drug-induced repolarization abnormalities, while those with reduced inward current density (fast/late Na⁺ and Ca²⁺ currents) exhibit high susceptibility to depolarization abnormalities. Repolarization abnormalities in silico predict clinical risk for all compounds with 89% accuracy. Drug-induced changes in biomarkers are in overall agreement across different assays: in silico AP duration changes reflect the ones observed in rabbit QT interval and hiPS-CMs Ca²⁺-transient, and simulated upstroke velocity captures variations in rabbit QRS complex. Our results demonstrate that human in silico drug trials constitute a powerful methodology for prediction of clinical pro-arrhythmic cardiotoxicity, ready for integration in the existing drug safety assessment pipelines.

Microdomain-Specific Modulation of L-Type Calcium Channels Leads to Triggered Ventricular Arrhythmia in Heart Failure

Supplemental Digital Content is available in the text.

Rationale:

Disruption in subcellular targeting of Ca²⁺ signaling complexes secondary to changes in cardiac myocyte structure may contribute to the pathophysiology of a variety of cardiac diseases, including heart failure (HF) and certain arrhythmias.

Objective:

To explore microdomain-targeted remodeling of ventricular L-type Ca²⁺ channels (LTCCs) in HF.

Methods and Results:

Super-resolution scanning patch-clamp, confocal and fluorescence microscopy were used to explore the distribution of single LTCCs in different membrane microdomains of nonfailing and failing human and rat ventricular myocytes. Disruption of membrane structure in both species led to the redistribution of functional LTCCs from their canonical location in transversal tubules (T-tubules) to the non-native crest of the sarcolemma, where their open probability was dramatically increased (0.034±0.011 versus 0.154±0.027, P<0.001). High open probability was linked to enhance calcium–calmodulin kinase II–mediated phosphorylation in non-native microdomains and resulted in an elevated I_Ca,L window current, which contributed to the development of early afterdepolarizations. A novel model of LTCC function in HF was developed; after its validation with experimental data, the model was used to ascertain how HF-induced T-tubule loss led to altered LTCC function and early afterdepolarizations. The HF myocyte model was then implemented in a 3-dimensional left ventricle model, demonstrating that such early afterdepolarizations can propagate and initiate reentrant arrhythmias.

Conclusions:

Microdomain-targeted remodeling of LTCC properties is an important event in pathways that may contribute to ventricular arrhythmogenesis in the settings of HF-associated remodeling. This extends beyond the classical concept of electric remodeling in HF and adds a new dimension to cardiovascular disease.

Electrophysiology of Heart Failure Using a Rabbit Model: From the Failing Myocyte to Ventricular Fibrillation

Heart failure is a leading cause of death, yet its underlying electrophysiological (EP) mechanisms are not well understood. In this study, we use a multiscale approach to analyze a model of heart failure and connect its results to features of the electrocardiogram (ECG). The heart failure model is derived by modifying a previously validated electrophysiology model for a healthy rabbit heart. Specifically, in accordance with the heart failure literature, we modified the cell EP by changing both membrane currents and calcium handling. At the tissue level, we modeled the increased gap junction lateralization and lower conduction velocity due to downregulation of Connexin 43. At the biventricular level, we reduced the apex-to-base and transmural gradients of action potential duration (APD). The failing cell model was first validated by reproducing the longer action potential, slower and lower calcium transient, and earlier alternans characteristic of heart failure EP. Subsequently, we compared the electrical wave propagation in one dimensional cables of healthy and failing cells. The validated cell model was then used to simulate the EP of heart failure in an anatomically accurate biventricular rabbit model. As pacing cycle length decreases, both the normal and failing heart develop T-wave alternans, but only the failing heart shows QRS alternans (although moderate) at rapid pacing. Moreover, T-wave alternans is significantly more pronounced in the failing heart. At rapid pacing, APD maps show areas of conduction block in the failing heart. Finally, accelerated pacing initiated wave reentry and breakup in the failing heart. Further, the onset of VF was not observed with an upregulation of SERCA, a potential drug therapy, using the same protocol. The changes introduced at the cell and tissue level have increased the failing heart’s susceptibility to dynamic instabilities and arrhythmias under rapid pacing. However, the observed increase in arrhythmogenic potential is not due to a steepening of the restitution curve (not present in our model), but rather to a novel blocking mechanism.

Ventricular fibrillation (VF) is one of the leading causes of sudden death. During VF, the electrical wave of activation in the heart breaks up chaotically. Consequently, the heart is unable to contract synchronously and pump blood to the rest of the body. In our work we formulate and validate a model of heart failure (HF) that allows us to evaluate the arrhythmogenic potential of individual and combined electrophysiological changes. In diagnostic cardiology, the electrocardiogram (ECG) is one of the most commonly used tools for detecting abnormalities in the heart electrophysiology. One of our goals is to use our numerical model to link changes at the cellular and tissue level in a failing heart to a numerically computed ECG. This allows us to characterize the precursor to and the risk of VF. In order to understand the mechanisms underlying VF in HF, we design a test that simulates a HF patient performing physical exercise. We show that under fast heart rates with changes in pacing, HF patients are more prone to VF due to a new conduction blocking mechanism. In the long term, our mathematical model is suitable for investigating the effect of drug therapies in HF.