Stochastic spontaneous calcium release events and sodium channelopathies promote ventricular arrhythmias

Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 μM vs. 22 μM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Patchy fibrosis promotes trigger-substrate interactions that both generate and maintain atrial fibrillation

Fibrosis has been mechanistically linked to arrhythmogenesis in multiple cardiovascular conditions, including atrial fibrillation (AF). Previous studies have demonstrated that fibrosis can create functional barriers to conduction which may promote excitation wavebreak and the generation of re-entry, while also acting to pin re-entrant excitation in stable rotors during AF. However, few studies have investigated the role of fibrosis in the generation of AF triggers in detail. We apply our in-house computational framework to study the impact of fibrosis on the generation of AF triggers and trigger-substrate interactions in two- and three-dimensional atrial tissue models. Our models include a reduced and efficient description of stochastic, spontaneous cellular triggers as well as a simple model of heterogeneous inter-cellular coupling. Our results demonstrate that fibrosis promotes the emergence of focal excitations, primarily through reducing the electrotonic load on individual fibre strands. This enables excitation to robustly initiate within these single strands before spreading to neighbouring strands and inducing a full tissue focal excitation. Enhanced conduction block can allow trigger-substrate interactions that result in the emergence of complex, re-entrant excitation patterns. This study provides new insight into the mechanisms by which fibrosis promotes the triggers and substrate necessary to induce and sustain arrhythmia.

Attenuation of stretch-induced arrhythmias following chemical ablation of Purkinje fibres, in isolated rabbit hearts

Purkinje fibres (PFs) play an important role in some ventricular arrhythmias and acute ventricular stretch can evoke mechanically-induced arrhythmias. We tested whether Purkinje fibres, play a role in these arrhythmias. Pseudo-ECGs were recorded in isolated, Langendorff-perfused, rabbit hearts in which the left ventricular endocardial surface was also irrigated with Tyrode, via an indwelling catheter placed in the left ventricular lumen. The number and period of ectopic activations was measured during left ventricular lumen inflation via an indwelling fluid-filled balloon (500 μL added over 2 s and maintained for 15 s in total). Mechanically-induced arrhythmias occurred in 70% of balloon inflations: they were maximal in the first 5 s and ceased within 15 s. Brief, (10 s) irrigation of the left ventricular lumen with Lugol solution (IK/I₂), via the indwelling catheter, reduced inflation-induced ectopics by 98% (p < 0.05). Ablation of endocardial PFs by Lugol was confirmed by Triphenyltetrazolium Chloride staining. Optical mapping revealed the left ventricular epicardial activation patterns of ectopics could have PF-mediated and focal sources. In silico modelling predicted ectopic sources originating in the endocardial region propagate to and through the Purkinje fibres network. Acute distention-induced ectopics are multi-focal, their attenuation by Lugol, their activation patterns and in silico modelling indicate a participation of Purkinje fibres in these arrhythmias.

Malignant Purkinje Ectopy induced by Sodium-channel Blocker

BACKGROUND

Sodium-Channel Blocker (SCB) infusion is used to unmask the ECG pattern of Brugada syndrome. The test may also induce premature ventricular complexes (PVCs) in individuals without Brugada pattern, the clinical relevance of which is little known.

OBJECTIVE

To describe the prevalence of short-coupled (Sc) PVCs induced by ajmaline or flecainide in patients with suspected or documented severe ventricular arrhythmias.

METHODS

We reviewed the SCB tests performed in 335 patients with suspected ventricular arrhythmias and structurally normal hearts in 9 centers. ScPVCs were defined as frequent and repetitive PVCs with an R-on-T pattern on SCB tests. Repeated SCB tests were performed in seven patients and electrophysiological mapping of ScPVCs in nine.

RESULTS

Sixteen patients (8 males, 36±11 years) showed ScPVCs and were included. ScPVCs appeared 229±118 sec after the start of the infusion and displayed coupling intervals of 288±28 ms. ScPVC patterns were monomorphic in 10 and polymorphic in 5 patients, originating from the Purkinje system in mapped patients. Repetitive PVCs were induced in 15 (94%) patients including polymorphic ventricular tachycardias in 9 (56%). SCB was repeated 45 (IQR 0.6-46) months later and induced identical ScPVC in all. SCB test was the only mean to reveal the malignant arrhythmia in six patients. Catheter ablation was performed in 9 patients, resulting in arrhythmia elimination in 8 with a follow-up of 6 (IQR 2-9) years.

CONCLUSION

SCB can induce ScPVC, mostly from Purkinje tissue, in a small subset of patients with idiopathic ventricular arrhythmias. Its high reproducibility suggests a distinct individual mechanism.

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.

Delayed afterdepolarization-induced triggered activity in cardiac purkinje cells mediated through cytosolic calcium diffusion waves

Cardiac Purkinje cells (PCs) are more susceptible to action potential abnormalities as compared to ventricular myocytes (VMs), which could be associated with their distinct intracellular calcium handling. We developed a detailed biophysical model of a mouse cardiac PC, which importantly reproduces the experimentally observed biphasic cytosolic calcium waves. The model includes a stochastic gating formulation for the opening and closing of ryanodine receptor (RyR) channels, simulated with a Monte Carlo method, to accurately reproduce cytosolic calcium wave propagation and the effects of spontaneous calcium release events. Simulations predict that during an action potential, smaller cytosolic calcium wavelets propagated from the sarcolemma towards the center of the cell and initiated larger magnitude cell‐wide calcium waves via a calcium‐induced‐calcium release mechanism. In the presence of RyR mutations, frequent spontaneous calcium leaks from sarcoplasmic reticulum (SR) initiated calcium waves, which upon reaching the cell periphery produced delayed afterdepolarizations (DADs) via sodium‐calcium exchanger (NCX) and T‐type calcium (I_CaT) channel activation. In the presence of isoproterenol‐mediated effects, DADs induced triggered activity by reactivation of fast sodium channels. Based on our model, we found that the activation of either L‐type calcium channels (I_CaL), I_CaT, sodium‐potassium exchanger (I_NaK) or NCX is sufficient for occurrence of triggered activity; however, a partial blockade of I_CaT or I_NaK is essential for its successful termination. Our modeling study highlights valuable insights into the mechanisms of DAD‐induced triggered activity mediated via cytosolic calcium waves in cardiac PCs and may elucidate the increased arrhythmogeneity in PCs.

Arrhythmia mechanisms and spontaneous calcium release: Bi-directional coupling between re-entrant and focal excitation

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation.

Ventricular Endocardial Tissue Geometry Affects Stimulus Threshold and Effective Refractory Period

BACKGROUND

Understanding the biophysical processes by which electrical stimuli applied to cardiac tissue may result in local activation is important in both the experimental and clinical electrophysiology laboratory environments, as well as for gaining a more in-depth knowledge of the mechanisms of focal-trigger-induced arrhythmias. Previous computational models have predicted that local myocardial tissue architecture alone may significantly modulate tissue excitability, affecting both the local stimulus current required to excite the tissue and the local effective refractory period (ERP). In this work, we present experimental validation of this structural modulation of local tissue excitability on the endocardial tissue surface, use computational models to provide mechanistic understanding of this phenomena in relation to localized changes in electrotonic loading, and demonstrate its implications for the capture of afterdepolarizations.

METHODS AND RESULTS

Experiments on rabbit ventricular wedge preparations showed that endocardial ridges (surfaces of negative mean curvature) had a stimulus capture threshold that was 0.21 ± 0.03 V less than endocardial grooves (surfaces of positive mean curvature) for pairwise comparison (24% reduction, corresponding to 56.2 ± 6.4% of the energy). When stimulated at the minimal stimulus strength for capture, ridge locations showed a shorter ERP than grooves (n = 6, mean pairwise difference 7.4 ± 4.2 ms). When each site was stimulated with identical-strength stimuli, the difference in ERP was further increased (mean pairwise difference 15.8 ± 5.3 ms). Computational bidomain models of highly idealized cylindrical endocardial structures qualitatively agreed with these findings, showing that such changes in excitability are driven by structural modulation in electrotonic loading, quantifying this relationship as a function of surface curvature. Simulations further showed that capture of delayed afterdepolarizations was more likely in trabecular ridges than grooves, driven by this difference in loading.

CONCLUSIONS

We have demonstrated experimentally and explained mechanistically in computer simulations that the ability to capture tissue on the endocardial surface depends upon the local tissue architecture. These findings have important implications for deepening our understanding of excitability differences related to anatomical structure during stimulus application that may have important applications in the translation of novel experimental optogenetics pacing strategies. The uncovered preferential vulnerability to capture of afterdepolarizations of endocardial ridges, compared to grooves, provides important insight for understanding the mechanisms of focal-trigger-induced arrhythmias.

Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data.