Small size ionic heterogeneities in the human heart can attract rotors

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.

Dynamical anchoring of distant arrhythmia sources by fibrotic regions via restructuring of the activation pattern

Rotors are functional reentry sources identified in clinically relevant cardiac arrhythmias, such as ventricular and atrial fibrillation. Ablation targeting rotor sites has resulted in arrhythmia termination. Recent clinical, experimental and modelling studies demonstrate that rotors are often anchored around fibrotic scars or regions with increased fibrosis. However, the mechanisms leading to abundance of rotors at these locations are not clear. The current study explores the hypothesis whether fibrotic scars just serve as anchoring sites for the rotors or whether there are other active processes which drive the rotors to these fibrotic regions. Rotors were induced at different distances from fibrotic scars of various sizes and degree of fibrosis. Simulations were performed in a 2D model of human ventricular tissue and in a patient-specific model of the left ventricle of a patient with remote myocardial infarction. In both the 2D and the patient-specific model we found that without fibrotic scars, the rotors were stable at the site of their initiation. However, in the presence of a scar, rotors were eventually dynamically anchored from large distances by the fibrotic scar via a process of dynamical reorganization of the excitation pattern. This process coalesces with a change from polymorphic to monomorphic ventricular tachycardia.

Rotors are waves of cardiac excitation like a tornado causing cardiac arrhythmia. Recent research shows that they are found in ventricular and atrial fibrillation. Burning (via ablation) the site of a rotor can result in the termination of the arrhythmia. Recent studies showed that rotors are often anchored to regions surrounding scar tissue, where part of the tissue still survived called fibrotic tissue. However, it is unclear why these rotors anchor to these locations. Therefore, in this work, we investigated why rotors are so abundant in fibrotic tissue with the help of computer simulations. We performed simulations in a 2D model of human ventricular tissue and in a patient-specific model of a patient with an infarction. We found that even when rotors are initially at large distances from the fibrotic region, they are attracted by this region, to finally end up at the fibrotic tissue. We called this process dynamical anchoring and explained how the process works.

Myocyte Remodeling Due to Fibro-Fatty Infiltrations Influences Arrhythmogenicity

The onset of cardiac arrhythmias depends on the electrophysiological and structural properties of cardiac tissue. Electrophysiological remodeling of myocytes due to the presence of adipocytes constitutes a possibly important pathway in the pathogenesis of atrial fibrillation. In this paper we perform an in-silico study of the effect of such myocyte remodeling on the onset of atrial arrhythmias and study the dynamics of arrhythmia sources—spiral waves. We use the Courtemanche model for atrial myocytes and modify their electrophysiological properties based on published cellular electrophysiological measurements in myocytes co-cultered with adipocytes (a 69–87 % increase in APD₉₀ and an increase of the RMP by 2.5–5.5 mV). In a generic 2D setup we show that adipose tissue remodeling substantially affects the spiral wave dynamics resulting in complex arrhythmia and such arrhythmia can be initiated under high frequency pacing if the size of the remodeled tissue is sufficiently large. These results are confirmed in simulations with an anatomically accurate model of the human atria.

Optogenetics enables real-time spatiotemporal control over spiral wave dynamics in an excitable cardiac system

Propagation of non-linear waves is key to the functioning of diverse biological systems. Such waves can organize into spirals, rotating around a core, whose properties determine the overall wave dynamics. Theoretically, manipulation of a spiral wave core should lead to full spatiotemporal control over its dynamics. However, this theory lacks supportive evidence (even at a conceptual level), making it thus a long-standing hypothesis. Here, we propose a new phenomenological concept that involves artificially dragging spiral waves by their cores, to prove the aforementioned hypothesis in silico, with subsequent in vitro validation in optogenetically modified monolayers of rat atrial cardiomyocytes. We thereby connect previously established, but unrelated concepts of spiral wave attraction, anchoring and unpinning to demonstrate that core manipulation, through controlled displacement of heterogeneities in excitable media, allows forced movement of spiral waves along pre-defined trajectories. Consequently, we impose real-time spatiotemporal control over spiral wave dynamics in a biological system.

From a spinning galaxy to a swarm of honeybees, rotating spirals are widespread in nature. Even within the muscles of the heart, waves of electrical activity sometimes rotate spirally, leading to irregular heart rhythms or arrhythmia – a condition that can be fatal.

Irrespective of where they occur, spiral waves organize around a center or core with different biophysical properties compared to the rest of the medium. The properties of the core determine the overall dynamics of the spiral. This means that, theoretically, it should be possibly to completely control a spiral wave just by manipulating its core.

Now, Majumder, Feola et al. have tested this long-standing hypothesis using a combination of computer modeling and experiments with single layers of rat heart cells grown in a laboratory. First, the heart cells were genetically modified so that their electrical properties could be altered with light; in other words, the cells were put under optical control. Next, by using of a narrow beam of light, Majumder, Feola et al. precisely controlled the electrical properties of a small number of cells, which then attracted and supported a rotating spiral wave by acting as its new core. Moving the light beam allowed the core of the spiral wave to be shifted too, meaning the spiral wave could now be steered along any desired path in the cell layer.

Majumder, Feola et al. hope that these underlying principles may one day provide the basis of new treatments for irregular heartbeats that are more effective and less damaging to the heart than existing options. Yet first, more work is needed to translate these findings from single layers of cells to actual hearts.

Wavelength and Fibrosis Affect Phase Singularity Locations During Atrial Fibrillation

The mechanisms underlying atrial fibrillation (AF), the most common sustained cardiac rhythm disturbance, remain elusive. Atrial fibrosis plays an important role in the development of AF and rotor dynamics. Both electrical wavelength (WL) and the degree of atrial fibrosis change as AF progresses. However, their combined effect on rotor core location remains unknown. The aim of this study was to analyze the effects of WL change on rotor core location in both fibrotic and non-fibrotic atria. Three patient specific fibrosis distributions (total fibrosis content: 16.6, 22.8, and 19.2%) obtained from clinical imaging data of persistent AF patients were incorporated in a bilayer atrial computational model. Fibrotic effects were modeled as myocyte-fibroblast coupling + conductivity remodeling; structural remodeling; ionic current changes + conductivity remodeling; and combinations of these methods. To change WL, action potential duration (APD) was varied from 120 to 240ms, representing the range of clinically observed AF cycle length, by modifying the inward rectifier potassium current (IK1) conductance between 80 and 140% of the original value. Phase singularities (PSs) were computed to identify rotor core locations. Our results show that IK1 conductance variation resulted in a decrease of APD and WL across the atria. For large WL in the absence of fibrosis, PSs anchored to regions with high APD gradient at the center of the left atrium (LA) anterior wall and near the junctions of the inferior pulmonary veins (PVs) with the LA. Decreasing the WL induced more PSs, whose distribution became less clustered. With fibrosis, PS locations depended on the fibrosis distribution and the fibrosis implementation method. The proportion of PSs in fibrotic areas and along the borders varied with both WL and fibrosis modeling method: for patient one, this was 4.2-14.9% as IK1 varied for the structural remodeling representation, but 12.3-88.4% using the combination of structural remodeling with myocyte-fibroblast coupling. The degree and distribution of fibrosis and the choice of implementation technique had a larger effect on PS locations than the WL variation. Thus, distinguishing the fibrotic mechanisms present in a patient is important for interpreting clinical fibrosis maps to create personalized models.

Localized Optogenetic Targeting of Rotors in Atrial Cardiomyocyte Monolayers

BACKGROUND

Recently, a new ablation strategy for atrial fibrillation has emerged, which involves the identification of rotors (ie, local drivers) followed by the localized targeting of their core region by ablation. However, this concept has been subject to debate because the mode of arrhythmia termination remains poorly understood, as dedicated models and research tools are lacking. We took a unique optogenetic approach to induce and locally target a rotor in atrial monolayers.

METHODS AND RESULTS

Neonatal rat atrial cardiomyocyte monolayers expressing a depolarizing light-gated ion channel (Ca²⁺-translocating channelrhodopsin) were subjected to patterned illumination to induce single, stable, and centralized rotors by optical S1-S2 cross-field stimulation. Next, the core region of these rotors was specifically and precisely targeted by light to induce local conduction blocks of circular or linear shapes. Conduction blocks crossing the core region, but not reaching any unexcitable boundary, did not lead to termination. Instead, electric waves started to propagate along the circumference of block, thereby maintaining reentrant activity, although of lower frequency. If, however, core-spanning lines of block reached at least 1 unexcitable boundary, reentrant activity was consistently terminated by wave collision. Lines of block away from the core region resulted merely in rotor destabilization (ie, drifting).

CONCLUSIONS

Localized optogenetic targeting of rotors in atrial monolayers could lead to both stabilization and destabilization of reentrant activity. For termination, however, a line of block is required reaching from the core region to at least 1 unexcitable boundary. These findings may improve our understanding of the mechanisms involved in rotor-guided ablation.

Impact of number of co-existing rotors and inter-electrode distance on accuracy of rotor localization

BACKGROUND

Conflicting evidence exists on the efficacy of focal impulse and rotor modulation on atrial fibrillation ablation. A potential explanation is inaccurate rotor localization from multiple rotors coexistence and a relatively large (9-11mm) inter-electrode distance (IED) of the multi-electrode basket catheter.

METHODS AND RESULTS

We studied a numerical model of cardiac action potential to reproduce one through seven rotors in a two-dimensional lattice. We estimated rotor location using phase singularity, Shannon entropy and dominant frequency. We then spatially downsampled the time series to create IEDs of 2-30mm. The error of rotor localization was measured with reference to the dynamics of phase singularity at the original spatial resolution (IED=1mm). IED has a significant impact on the error using all the methods. When only one rotor is present, the error increases exponentially as a function of IED. At the clinical IED of 10mm, the error is 3.8mm (phase singularity), 3.7mm (dominant frequency), and 11.8mm (Shannon entropy). When there are more than one rotors, the error of rotor localization increases 10-fold. The error based on the phase singularity method at the clinical IED of 10mm ranges from 30.0mm (two rotors) to 96.1mm (five rotors).

CONCLUSIONS

The magnitude of error of rotor localization using a clinically available basket catheter, in the presence of multiple rotors might be high enough to impact the accuracy of targeting during AF ablation. Improvement of catheter design and development of high-density mapping catheters may improve clinical outcomes of FIRM-guided AF ablation.

Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms

AIMS

Catheter ablation is an effective technique for terminating atrial arrhythmia. However, given a high atrial fibrillation (AF) recurrence rate, optimal ablation strategies have yet to be defined. Computer modelling can be a powerful aid but modelling of fibrosis, a major factor associated with AF, is an open question. Several groups have proposed methodologies based on imaging data, but no comparison to determine which methodology best corroborates clinically observed reentrant behaviour has been performed. We examined several methodologies to determine the best method for capturing fibrillation dynamics.

METHODS AND RESULTS

Patient late gadolinium-enhanced magnetic resonance imaging data were transferred onto a bilayer atrial computer model and used to assign fibrosis distributions. Fibrosis was modelled as conduction disturbances (lower conductivity, edge splitting, or percolation), transforming growth factor-β1 ionic channel effects, myocyte-fibroblast coupling, and combinations of the preceding. Reentry was induced through pulmonary vein ectopy and the ensuing rotor dynamics characterized. Non-invasive electrocardiographic imaging data of the patients in AF was used for comparison. Electrograms were computed and the fractionation durations measured over the surface. Edge splitting produced more phase singularities from wavebreaks than the other representations. The number of phase singularities seen with percolation was closer to the clinical values. Addition of fibroblast coupling had an organizing effect on rotor dynamics. Simple tissue conductivity changes with ionic changes localized rotors over fibrosis which was not observed with clinical data.

CONCLUSION

The specific representation of fibrosis has a large effect on rotor dynamics and needs to be carefully considered for patient specific modelling.

The Cardiac Electrophysiology Web Lab

Computational modeling of cardiac cellular electrophysiology has a long history, and many models are now available for different species, cell types, and experimental preparations. This success brings with it a challenge: how do we assess and compare the underlying hypotheses and emergent behaviors so that we can choose a model as a suitable basis for a new study or to characterize how a particular model behaves in different scenarios? We have created an online resource for the characterization and comparison of electrophysiological cell models in a wide range of experimental scenarios. The details of the mathematical model (quantitative assumptions and hypotheses formulated as ordinary differential equations) are separated from the experimental protocol being simulated. Each model and protocol is then encoded in computer-readable formats. A simulation tool runs virtual experiments on models encoded in CellML, and a website (https://chaste.cs.ox.ac.uk/WebLab) provides a friendly interface, allowing users to store and compare results. The system currently contains a sample of 36 models and 23 protocols, including current-voltage curve generation, action potential properties under steady pacing at different rates, restitution properties, block of particular channels, and hypo-/hyperkalemia. This resource is publicly available, open source, and free, and we invite the community to use it and become involved in future developments. Investigators interested in comparing competing hypotheses using models can make a more informed decision, and those developing new models can upload them for easy evaluation under the existing protocols, and even add their own protocols.

Pro-arrhythmogenic effects of CACNA1C G1911R mutation in human ventricular tachycardia: insights from cardiac multi-scale models

Mutations in the CACNA1C gene are associated with ventricular tachycardia (VT). Although the CACNA1C mutations were well identified in patients with cardiac arrhythmias, mechanisms by which cardiac arrhythmias are generated in such genetic mutation conditions remain unclear. In this study, we identified a novel mechanism of VT resulted from enhanced repolarization dispersion which is a key factor for arrhythmias in the CACNA1C G1911R mutation using multi-scale computational models of the human ventricle. The increased calcium influx in the mutation prolonged action potential duration (APD), produced steepened action potential duration restitution (APDR) curves as well as augmented membrane potential differences among different cell types during repolarization, increasing transmural dispersion of repolarization (DOR) and the spatial and temporal heterogeneity of cardiac electrical activities. Consequentially, the vulnerability to unidirectional conduction block in response to a premature stimulus increased at tissue level in the G1911R mutation. The increased functional repolarization dispersion anchored reentrant excitation waves in tissue and organ models, facilitating the initiation and maintenance of VT due to less meandering rotor tip. Thus, the increased repolarization dispersion caused by the G1911R mutation is a primary factor that may primarily contribute to the genesis of cardiac arrhythmias in Timothy Syndrome.

Global alternans instability and its effect on non-linear wave propagation: dynamical Wenckebach block and self terminating spiral waves

The main mechanism of formation of reentrant cardiac arrhythmias is via formation of waveblocks at heterogeneities of cardiac tissue. We report that heterogeneity and the area of waveblock can extend itself in space and can result formation of new additional sources, or termination of existing sources of arrhythmias. This effect is based on a new form of instability, which we coin as global alternans instability (GAI). GAI is closely related to the so-called (discordant) alternans instability, however its onset is determined by the global properties of the APD-restitution curve and not by its slope. The APD-restitution curve relates the duration of the cardiac pulse (APD) to the time interval between the pulses, and can easily be measured in an experimental or even clinical setting. We formulate the conditions for the onset of GAI, study its manifestation in various 1D and 2D situations and discuss its importance for the onset of cardiac arrhythmias.

A Computational Study of the Factors Influencing the PVC-Triggering Ability of a Cluster of Early Afterdepolarization-Capable Myocytes

Premature ventricular complexes (PVCs), which are abnormal impulse propagations in cardiac tissue, can develop because of various reasons including early afterdepolarizations (EADs). We show how a cluster of EAD-generating cells (EAD clump) can lead to PVCs in a model of cardiac tissue, and also investigate the factors that assist such clumps in triggering PVCs. In particular, we study, through computer simulations, the effects of the following factors on the PVC-triggering ability of an EAD clump: (1) the repolarization reserve (RR) of the EAD cells; (2) the size of the EAD clump; (3) the coupling strength between the EAD cells in the clump; and (4) the presence of fibroblasts in the EAD clump. We find that, although a low value of RR is necessary to generate EADs and hence PVCs, a very low value of RR leads to low-amplitude EAD oscillations that decay with time and do not lead to PVCs. We demonstrate that a certain threshold size of the EAD clump, or a reduction in the coupling strength between the EAD cells, in the clump, is required to trigger PVCs. We illustrate how randomly distributed inexcitable obstacles, which we use to model collagen deposits, affect PVC-triggering by an EAD clump. We show that the gap-junctional coupling of fibroblasts with myocytes can either assist or impede the PVC-triggering ability of an EAD clump, depending on the resting membrane potential of the fibroblasts and the coupling strength between the myocyte and fibroblasts. We also find that the triggering of PVCs by an EAD clump depends sensitively on factors like the pacing cycle length and the distribution pattern of the fibroblasts.

Post-repolarization refractoriness increases vulnerability to block and initiation of reentrant impulses in heterogeneous infarcted myocardium

Myocardial infarction causes remodeling of the tissue structure and the density and kinetics of several ion channels in the cell membrane. Heterogeneities in refractory period (ERP) have been shown to occur in the infarct border zone and have been proposed to lead to initiation of arrhythmias. The purpose of this study is to quantify the window of vulnerability (WV) to block and initiation of reentrant impulses in myocardium with ERP heterogeneities using computer simulations. We found that ERP transitions at the border between normal ventricular cells (NZ) with different ERPs are smooth, whereas ERP transitions between NZ and infarct border zone cells (IZ) are abrupt. The profile of the ERP transitions is a combination of electrotonic interaction between NZ and IZ cells and the characteristic post-repolarization refractoriness (PRR) of IZ cells. ERP heterogeneities between NZ and IZ cells are more vulnerable to block and initiation of reentrant impulses than ERP heterogeneities between NZ cells. The relationship between coupling intervals of premature impulses (V1V2) and coupling intervals between premature and first reentrant impulses (V2T1) at NZ/NZ and NZ/IZ borders is inverse (i.e. the longer the coupling intervals of premature impulses the shorter the coupling interval between the premature and first reentrant impulses); this is in contrast with the reported V1V2/V2T1 relationship measured during initiation of reentrant impulses in canine infarcted hearts which is direct.

IN CONCLUSION

(1) ERP transitions at the NZ-IZ border are abrupt as a consequence of PRR; (2) PRR increases the vulnerability to block and initiation of reentrant impulses in heterogeneous myocardium; (3) V1V2/V2T1 relationships measured at ERP heterogeneities in the computer model and in experimental canine infarcts are not consistent. Therefore, it is likely that other mechanisms like micro and/or macro structural heterogeneities also contribute to initiation of reentrant impulses in infarcted hearts.