A Dynamical Threshold for Cardiac Delayed Afterdepolarization-Mediated Triggered Activity

Integrative human atrial modelling unravels interactive protein kinase A and Ca2+/calmodulin-dependent protein kinase II signalling as key determinants of atrial arrhythmogenesis

AIMS

Atrial fibrillation (AF), the most prevalent clinical arrhythmia, is associated with atrial remodelling manifesting as acute and chronic alterations in expression, function, and regulation of atrial electrophysiological and Ca2+-handling processes. These AF-induced modifications crosstalk and propagate across spatial scales creating a complex pathophysiological network, which renders AF resistant to existing pharmacotherapies that predominantly target transmembrane ion channels. Developing innovative therapeutic strategies requires a systems approach to disentangle quantitatively the pro-arrhythmic contributions of individual AF-induced alterations.

METHODS AND RESULTS

Here, we built a novel computational framework for simulating electrophysiology and Ca2+-handling in human atrial cardiomyocytes and tissues, and their regulation by key upstream signalling pathways [i.e. protein kinase A (PKA), and Ca2+/calmodulin-dependent protein kinase II (CaMKII)] involved in AF-pathogenesis. Populations of atrial cardiomyocyte models were constructed to determine the influence of subcellular ionic processes, signalling components, and regulatory networks on atrial arrhythmogenesis. Our results reveal a novel synergistic crosstalk between PKA and CaMKII that promotes atrial cardiomyocyte electrical instability and arrhythmogenic triggered activity. Simulations of heterogeneous tissue demonstrate that this cellular triggered activity is further amplified by CaMKII- and PKA-dependent alterations of tissue properties, further exacerbating atrial arrhythmogenesis.

CONCLUSIONS

Our analysis reveals potential mechanisms by which the stress-associated adaptive changes turn into maladaptive pro-arrhythmic triggers at the cellular and tissue levels and identifies potential anti-AF targets. Collectively, our integrative approach is powerful and instrumental to assemble and reconcile existing knowledge into a systems network for identifying novel anti-AF targets and innovative approaches moving beyond the traditional ion channel-based strategy.

EDEN (Electrocardiographic Radial Depth Navigation): A Novel Approach to Navigate Inside Heart Muscle

BACKGROUND

Intramyocardial guidewire navigation is a novel technique that allows free transcatheter movement within ventricular muscle. Guidewire radial depth, between endocardial and epicardial surfaces, is ambiguous by x-ray and echocardiography.

OBJECTIVES

The aim of this study was to develop a simple tool, EDEN (Electrocardiographic Radial Depth Navigation), to indicate radial depth during intramyocardial guidewire navigation. Combined with routine imaging, EDEN facilitates a new family of intramyocardial catheter procedures to slice, reshape, pace, and ablate the heart.

METHODS

We mapped intramyocardial electrograms of left and right ventricular walls and septum during open- and closed-chest swine procedures (N = 53), including MIRTH (Myocardial Intramural Remodeling by Transvenous Tether) ventriculoplasty. We identified radial depth-dependent features on unipolar electrograms. We developed a machine learning-based classifier to indicate categorical position, and modeled the findings in silico to test understanding of the physiology.

RESULTS

EDEN signatures distinguished 5 depth zones throughout left and right ventricular free walls and interventricular septum. Relative ST-segment elevation magnitude best discriminated position and was maximum (40.1 ± 6.5 mV) in the midmyocardium. Subendocardial positions exhibited dominant Q waves with lower-amplitude ST segments (16.8 ± 5.8 mV), whereas subepicardial positions exhibited dominant R waves with lower-amplitude ST segments (15.7 ± 4.8 mV). EDEN was unaffected by pacing-induced left bundle branch block. ST-segment elevation declined over minutes and reappeared after submillimeter guidewire manipulation. Modeling recapitulated EDEN features. The machine learning-based classifier was 97% accurate. EDEN successfully guided MIRTH ventriculoplasty.

CONCLUSIONS

EDEN provides a simple and reproducible real-time reflection of categorical guidewire-tip radial depth during intramyocardial guidewire navigation. Used in tandem with x-ray, EDEN enables novel, transcatheter, intramyocardial therapies such as MIRTH, SESAME (Septal Surfing Along Midline Endocardium), and cerclage ventriculoplasty.

R-on-T and the initiation of reentry revisited: Integrating old and new concepts

Initiation of reentry requires 2 factors: (1) a triggering event, most commonly focal excitations such as premature ventricular complexes (PVCs); and (2) a vulnerable substrate with regional dispersion of refractoriness and/or excitability, such as occurs during the T wave of the electrocardiogram when some areas of the ventricle have repolarized and recovered excitability but others have not. When the R wave of a PVC coincides in time with the T wave of the previous beat, this timing can lead to unidirectional block and initiation of reentry, known as the R-on-T phenomenon. Classically, the PVC triggering reentry has been viewed as arising focally from 1 region and propagating into another region whose recovery is delayed, resulting in unidirectional conduction block and reentry initiation. However, more recent evidence indicates that PVCs also can arise from the T wave itself. In the latter case, the PVC initiating reentry is not a separate event from the T wave but rather is causally generated from the repolarization gradient that manifests as the T wave. We call the former an "R-to-T" mechanism and the latter an "R-from-T" mechanism, which are initiation mechanisms distinct from each other. Both are important components of the R-on-T phenomenon and need to be taken into account when designing antiarrhythmic strategies. Strategies targeting suppression of triggers alone or vulnerable substrate alone may be appropriate in some instances but not in others. Preventing R-from-T arrhythmias requires suppressing the underlying dynamic tissue instabilities responsible for producing both triggers and substrate vulnerability simultaneously. The same principles are likely to apply to supraventricular arrhythmias.

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

Quantitative roles of ion channel dynamics on ventricular action potential

Mathematical models for the action potential (AP) generation of the electrically excitable cells including the heart are involved different mechanisms including the voltage-dependent currents with nonlinear time- and voltage-gating properties. From the shape of the AP waveforms to the duration of the refractory periods or heart rhythms are greatly affected by the functions describing the features or the quantities of these ion channels. In this work, a mathematical measure to analyze the regional contributions of voltage-gated channels is defined by dividing the AP into phases, epochs, and intervals of interest. The contribution of each time-dependent current for the newly defined cardiomyocyte model is successfully calculated and it is found that the contribution of dominant ion channels changes substantially not only for each phase but also for different regions of the cardiac AP. Besides, the defined method can also be applied in all Hodgkin-Huxley types of electrically excitable cell models to be able to understand the underlying dynamics better.

Estimating ectopic beat probability with simplified statistical models that account for experimental uncertainty

Ectopic beats (EBs) are cellular arrhythmias that can trigger lethal arrhythmias. Simulations using biophysically-detailed cardiac myocyte models can reveal how model parameters influence the probability of these cellular arrhythmias, however such analyses can pose a huge computational burden. Here, we develop a simplified approach in which logistic regression models (LRMs) are used to define a mapping between the parameters of complex cell models and the probability of EBs (P(EB)). As an example, in this study, we build an LRM for P(EB) as a function of the initial value of diastolic cytosolic Ca²⁺ concentration ([Ca²⁺]_iⁱⁿⁱ), the initial state of sarcoplasmic reticulum (SR) Ca²⁺ load ([Ca²⁺]_SRⁱⁿⁱ), and kinetic parameters of the inward rectifier K⁺ current (I_K1) and ryanodine receptor (RyR). This approach, which we refer to as arrhythmia sensitivity analysis, allows for evaluation of the relationship between these arrhythmic event probabilities and their associated parameters. This LRM is also used to demonstrate how uncertainties in experimentally measured values determine the uncertainty in P(EB). In a study of the role of [Ca²⁺]_SRⁱⁿⁱ uncertainty, we show a special property of the uncertainty in P(EB), where with increasing [Ca²⁺]_SRⁱⁿⁱ uncertainty, P(EB) uncertainty first increases and then decreases. Lastly, we demonstrate that I_K1 suppression, at the level that occurs in heart failure myocytes, increases P(EB).

An ectopic beat is an abnormal cellular electrical event which can trigger dangerous arrhythmias in the heart. Complex biophysical models of the cardiac myocyte can be used to reveal how cell properties affect the probability of ectopic beats. However, such analyses can pose a huge computational burden. We develop a simplified approach that enables a highly complex biophysical model to be reduced to a rather simple statistical model from which the functional relationship between myocyte model parameters and the probability of an ectopic beat is determined. We refer to this approach as arrhythmia sensitivity analysis. Given the efficiency of our approach, we also use it to demonstrate how uncertainties in experimentally measured myocyte model parameters determine the uncertainty in ectopic beat probability. We find that, with increasing model parameter uncertainty, the uncertainty in probability of ectopic beat first increases and then decreases. In general, our approach can efficiently analyze the relationship between cardiac myocyte parameters and the probability of ectopic beats and can be used to study how uncertainty of these cardiac myocyte parameters influences the ectopic beat probability.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Stabilizer Cell Gene Therapy: A Less-Is-More Strategy to Prevent Cardiac Arrhythmias

BACKGROUND

In cardiac gene therapy to improve contractile function, achieving gene expression in the majority of cardiac myocytes is essential. In preventing cardiac arrhythmias, however, this goal may not be as important since transduction efficiencies as low as 40% suppressed ventricular arrhythmias in genetically modified mice with catecholaminergic polymorphic ventricular tachycardia.

METHODS

Using computational modeling, we simulated 1-, 2-, and 3-dimensional tissue under a variety of conditions to test the ability of genetically engineered nonarrhythmogenic stabilizer cells to suppress triggered activity due to delayed or early afterdepolarizations.

RESULTS

Due to source-sink relationships in cardiac tissue, a minority (20%-50%) of randomly distributed stabilizer cells engineered to be nonarrhythmogenic can suppress the ability of arrhythmogenic cells to generate delayed and early afterdepolarizations-related arrhythmias. Stabilizer cell gene therapy strategy can be designed to correct a specific arrhythmogenic mutation, as in the catecholaminergic polymorphic ventricular tachycardia mice studies, or more generally to suppress delayed or early afterdepolarizations from any cause by overexpressing the inward rectifier K channel Kir2.1 in stabilizer cells.

CONCLUSIONS

This promising antiarrhythmic strategy warrants further testing in experimental models to evaluate its clinical potential.

Photoelectric Cardiac Pacing by Flexible and Degradable Amorphous Si Radial Junction Stimulators

Implanted pacemakers are usually bulky and rigid electronics that are constraint by limited battery lifetimes, and need to be installed and repaired via surgeries that risk secondary infection and injury. In this work, a flexible self-powered photoelectric cardiac stimulator is demonstrated based on hydrogenated amorphous Si (a-Si:H) radial p-i-n junctions (RJs), constructed upon standing Si nanowires grown directly on aluminum thin foils. The flexible RJ stimulators, with an open-circuit voltage of 0.67 V and short-circuit current density of 12.7 mA cm^-2 under standard AM1.5G illumination, can be conformally attached to the uneven tissue surface to pace heart-beating under modulated 650 nm laser illumination. In vivo pacing evaluations on porcine hearts show that the heart rate can be effectively controlled by the external photoelectric stimulations, to increase from the normal rate of 101-128 beating min^-1 . Importantly, the a-Si:H RJ units are highly biofriendly and biodegradable, with tunable lifetimes in phosphate-buffered saline environment controlled by surface coating and passivation, catering to the needs of short term or lasting cardiac pacing applications. This implantable a-Si:H RJ photoelectric stimulation strategy has the potential to establish eventually a self-powered, biocompatible, and conformable cardiac pacing technology for clinical therapy.

Transverse tubular network structures in the genesis of intracellular calcium alternans and triggered activity in cardiac cells

RATIONALE

The major role of a transverse-tubular (TT) network in a cardiac cell is to facilitate effective excitation-contraction coupling and signaling. The TT network structures are heterogeneous within a single cell, and vary between different types of cells and species. They are also remodeled in cardiac diseases. However, how different TT network structures predispose cardiac cells to arrhythmogenesis remains to be revealed.

OBJECTIVE

To systematically investigate the roles of TT network structure and the underlying mechanisms in the genesis of intracellular calcium (Ca2+) alternans and triggered activity (TA).

METHODS AND RESULTS

Based on recent experimental observations, different TT network structures, including uniformly and non-uniformly random TT distributions, were modeled in a cardiac cell model consisting of a three-dimensional network of Ca2+ release units (CRUs). Our simulations showed that both Ca2+ alternans and Ca2+ wave-mediated TA were promoted when the fraction of orphaned CRUs was in an intermediate range, but suppressed in cells exhibiting either well-organized TT networks or low TT densities. Ca2+ alternans and TA could be promoted by low TT densities when the cells were small or the CRU coupling was strong. Both alternans and TA occurred more easily in uniformly random TT networks than in non-uniformly random TT networks. Subcellular spatially discordant Ca2+ alternans was promoted by non-uniformly random TT networks but suppressed by increasing CRU coupling strength. These mechanistic insights provide a holistic understanding of the effects of TT network structure on the susceptibility to arrhythmogenesis.

CONCLUSIONS

The TT network plays important roles in promoting Ca2+ alternans and TA, and different TT network structures may predispose cardiac cells differently to arrhythmogenesis.

Cardiac action potential repolarization revisited: early repolarization shows all-or-none behaviour

In healthy mammalian hearts the action potential (AP) waveform initiates and modulates each contraction, or heartbeat. As a result, AP height and duration are key physiological variables. In addition, rate-dependent changes in ventricular AP duration (APD), and variations in APD at a fixed heart rate are both reliable biomarkers of electrophysiological stability. Present guidelines for the likelihood that candidate drugs will increase arrhythmias rely on small changes in APD and Q-T intervals as criteria for safety pharmacology decisions. However, both of these measurements correspond to the final repolarization of the AP. Emerging clinical evidence draws attention to the early repolarization phase of the action potential (and the J-wave of the ECG) as an additional important biomarker for arrhythmogenesis. Here we provide a mechanistic background to this early repolarization syndrome by summarizing the evidence that both the initial depolarization and repolarization phases of the cardiac action potential can exhibit distinct time- and voltage-dependent thresholds, and also demonstrating that both can show regenerative all-or-none behaviour. An important consequence of this is that not all of the dynamics of action potential repolarization in human ventricle can be captured by data from single myocytes when these results are expressed as 'repolarization reserve'. For example, the complex pattern of cell-to-cell current flow that is responsible for AP conduction (propagation) within the mammalian myocardium can change APD and the Q-T interval of the electrocardiogram alter APD stability, and modulate responsiveness to pharmacological agents (such as Class III anti-arrhythmic drugs).

Multiscale Determinants of Delayed Afterdepolarization Amplitude in Cardiac Tissue

Spontaneous calcium (Ca) waves in cardiac myocytes underlie delayed afterdepolarizations (DADs) that trigger cardiac arrhythmias. How these subcellular/cellular events overcome source-sink factors in cardiac tissue to generate DADs of sufficient amplitude to trigger action potentials is not fully understood. Here, we evaluate quantitatively how factors at the subcellular scale (number of Ca wave initiation sites), cellular scale (sarcoplasmic reticulum (SR) Ca load), and tissue scale (synchrony of Ca release in populations of myocytes) determine DAD features in cardiac tissue using a combined experimental and computational modeling approach. Isolated patch-clamped rabbit ventricular myocytes loaded with Fluo-4 to image intracellular Ca were rapidly paced during exposure to elevated extracellular Ca (2.7 mmol/L) and isoproterenol (0.25 μmol/L) to induce diastolic Ca waves and subthreshold DADs. As the number of paced beats increased from 1 to 5, SR Ca content (assessed with caffeine pulses) increased, the number of Ca wave initiation sites increased, integrated Ca transients and DADs became larger and shorter in duration, and the latency period to the onset of Ca waves shortened with reduced variance. In silico analysis using a computer model of ventricular tissue incorporating these experimental measurements revealed that whereas all of these factors promoted larger DADs with higher probability of generating triggered activity, the latency period variance and SR Ca load had the greatest influences. Therefore, incorporating quantitative experimental data into tissue level simulations reveals that increased intracellular Ca promotes DAD-mediated triggered activity in tissue predominantly by increasing both the synchrony (decreasing latency variance) of Ca waves in nearby myocytes and SR Ca load, whereas the number of Ca wave initiation sites per myocyte is less important.