A method to quantify the dynamics and complexity of re-entry in computational models of ventricular fibrillation

Time resolution for wavefront and phase singularity tracking using activation maps in cardiac propagation models

The dynamics of cardiac fibrillation can be described by the number, the trajectory, the stability, and the lifespan of phase singularities (PSs). Accurate PS tracking is straightforward in simple uniform tissues but becomes more challenging as fibrosis, structural heterogeneity, and strong anisotropy are combined. In this paper, we derive a mathematical formulation for PS tracking in two-dimensional reaction-diffusion models. The method simultaneously tracks wavefronts and PS based on activation maps at full spatiotemporal resolution. PS tracking is formulated as a linear assignment problem solved by the Hungarian algorithm. The cost matrix incorporates information about distances between PS, chirality, and wavefronts. A graph of PS trajectories is generated to represent the creations and annihilations of PS pairs. Structure-preserving graph transformations are applied to provide a simplified description at longer observation time scales. The approach is validated in 180 simulations of fibrillation in four different types of substrates featuring, respectively, wavebreaks, ionic heterogeneities, fibrosis, and breakthrough patterns. The time step of PS tracking is studied in the range from 0.1 to 10 ms. The results show the benefits of improving time resolution from 1 to 0.1 ms. The tracking error rate decreases by an order of magnitude because the occurrence of simultaneous events becomes less likely. As observed on PS survival curves, the graph-based analysis facilitates the identification of macroscopically stable rotors despite wavefront fragmentation by fibrosis.

Effect of Heart Structure on Ventricular Fibrillation in the Rabbit: A Simulation Study

Ventricular fibrillation (VF) is a lethal condition that affects millions worldwide. The mechanism underlying VF is unstable reentrant electrical waves rotating around lines called filaments. These complex spatio-temporal patterns can be studied using both experimental and numerical methods. Computer simulations provide unique insights including high resolution dynamics throughout the heart and systematic control of quantities such as fiber orientation and cellular kinetics that are not feasible experimentally. Here we study filament dynamics using two bi-ventricular 3-D high-resolution rabbit heart geometries, one with detailed fine structure and another without fine structure. We studied filament dynamics using anisotropic and isotropic conductivities, and with four cellular action potential models with different recovery kinetics. Spiral wave dynamics observed in isotropic two-dimensional sheets were not predictive of the behavior in the whole heart. In 2-D the four cell models exhibited stable reentry, meandering spiral waves, and spiral-wave breakup. In the whole heart with fine structure, all simulation results exhibited complex dynamics reminiscent of fibrillation observed experimentally. In the whole heart without fine structure, anisotropy acted to destabilize filament dynamics although the number of filaments was reduced compared to the heart with structure. In addition, in isotropic hearts without structure the two cell models that exhibited meandering spiral waves in 2-D, stabilized into figure-of-eight surface patterns. We also studied the sensitivity of filament dynamics to computer system configuration and initial conditions. After large simulation times, different macroscopic results sometimes occurred across different system configurations, likely due to a lack of bitwise reproducibility. The study conclusions were insensitive to initial condition perturbations, however, the exact number of filaments over time and their trends were altered by these changes. In summary, we present the following new results. First, we provide a new cell model that resembles the surface patterns of VF in the rabbit heart both qualitatively and quantitatively. Second, filament dynamics in the whole heart cannot be predicted from spiral wave dynamics in 2-D and we identified anisotropy as one destabilizing factor. Third, the exact dynamics of filaments are sensitive to a variety of factors, so we suggest caution in their interpretation and their quantitative analyses.

Role of Apamin-Sensitive Calcium-Activated Small-Conductance Potassium Currents on the Mechanisms of Ventricular Fibrillation in Pacing-Induced Failing Rabbit Hearts

BACKGROUND

Ventricular fibrillation (VF) during heart failure is characterized by stable reentrant spiral waves (rotors). Apamin-sensitive small-conductance calcium-activated potassium currents (I_KAS) are heterogeneously upregulated in failing hearts. We hypothesized that I_KAS influences the location and stability of rotors during VF.

METHODS AND RESULTS

Optical mapping was performed on 9 rabbit hearts with pacing-induced heart failure. The epicardial right ventricular and left ventricular surfaces were simultaneously mapped in a Langendorff preparation. At baseline and after apamin (100 nmol/L) infusion, the action potential duration (APD₈₀) was determined, and VF was induced. Areas with a >50% increase in the maximum action potential duration (ΔAPD) after apamin infusion were considered to have a high I_KAS distribution. At baseline, the distribution density of phase singularities during VF in high I_KAS distribution areas was higher than in other areas (0.0035±0.0011 versus 0.0014±0.0010 phase singularities/pixel; P=0.004). In addition, high dominant frequencies also colocalized to high I_KAS distribution areas (26.0 versus 17.9 Hz; P=0.003). These correlations were eliminated during VF after apamin infusion, as the number of phase singularities (17.2 versus 11.0; P=0.009) and dominant frequencies (22.1 versus 16.2 Hz; P=0.022) were all significantly decreased. In addition, reentrant spiral waves became unstable after apamin infusion, and the duration of VF decreased.

CONCLUSIONS

The I_KAS current influences the mechanism of VF in failing hearts as phase singularities, high dominant frequencies, and reentrant spiral waves all correlated to areas of high I_KAS. Apamin eliminated this relationship and reduced VF vulnerability.

A Mathematical Model of Neonatal Rat Atrial Monolayers with Constitutively Active Acetylcholine-Mediated K+ Current

Atrial fibrillation (AF) is the most frequent form of arrhythmia occurring in the industrialized world. Because of its complex nature, each identified form of AF requires specialized treatment. Thus, an in-depth understanding of the bases of these arrhythmias is essential for therapeutic development. A variety of experimental studies aimed at understanding the mechanisms of AF are performed using primary cultures of neonatal rat atrial cardiomyocytes (NRAMs). Previously, we have shown that the distinct advantage of NRAM cultures is that they allow standardized, systematic, robust re-entry induction in the presence of a constitutively-active acetylcholine-mediated K⁺ current (I_KACh-c). Experimental studies dedicated to mechanistic explorations of AF, using these cultures, often use computer models for detailed electrophysiological investigations. However, currently, no mathematical model for NRAMs is available. Therefore, in the present study we propose the first model for the action potential (AP) of a NRAM with constitutively-active acetylcholine-mediated K⁺ current (I_KACh-c). The descriptions of the ionic currents were based on patch-clamp data obtained from neonatal rats. Our monolayer model closely mimics the action potential duration (APD) restitution and conduction velocity (CV) restitution curves presented in our previous in vitro studies. In addition, the model reproduces the experimentally observed dynamics of spiral wave rotation, in the absence and in the presence of drug interventions, and in the presence of localized myofibroblast heterogeneities.

A fundamentally important element in cardiac in silico research is a model for the cardiac cell. It provides a link between measurable characteristics at the subcellular level and biological processes at the whole cell level, thereby allowing the researcher to study mechanisms of cardiac arrhythmias from a molecular cell biological perspective. Such studies are of vast importance for the advancement of understanding of living systems from cells to patient populations. This paper is a joint in silico-experimental study in which we propose the first model for the action potential of an NRAM. To develop this model, we fitted patch-clamp data from recent literature, while additionally performing specific measurements of I_KACh-c in NRAMs. I_KACh-c is an important factor in atrial arrhythmogenesis and a promising target for pharmacological AF-management. The model reproduces in vitro results such as standard characteristics of AP morphology, restitution, and spiral wave dynamics in monolayers, with effects of a subsequent drug-intervention and in the presence of localized myofibroblast heterogeneities. Thus it can be used as a tool to provide computational support to a variety of systematic experimental studies that investigate the mechanisms underlying atrial fibrillation (AF) in NRAM cultures.

Biophysical Modeling Predicts Ventricular Tachycardia Inducibility and Circuit Morphology: A Combined Clinical Validation and Computer Modeling Approach

INTRODUCTION

Computational modeling of cardiac arrhythmogenesis and arrhythmia maintenance has made a significant contribution to the understanding of the underlying mechanisms of arrhythmia. We hypothesized that a cardiac model using personalized electro-anatomical parameters could define the underlying ventricular tachycardia (VT) substrate and predict reentrant VT circuits. We used a combined modeling and clinical approach in order to validate the concept.

METHODS AND RESULTS

Non-contact electroanatomic mapping studies were performed in 7 patients (5 ischemics, 2 non-ischemics). Three ischemic cardiomyopathy patients underwent a clinical VT stimulation study. Anatomical information was obtained from cardiac magnetic resonance imaging (CMR) including high-resolution scar imaging. A simplified biophysical mono-domain action potential model personalized with the patients' anatomical and electrical information was used to perform in silico VT stimulation studies for comparison. The personalized in silico VT stimulations were able to predict VT inducibility as well as the macroscopic characteristics of the VT circuits in patients who had clinical VT stimulation studies. The patients with positive clinical VT stimulation studies had wider distribution of action potential duration restitution curve (APD-RC) slopes and APDs than the patient with a negative VT stimulation study. The exit points of reentrant VT circuits encompassed a higher percentage of the maximum APD-RC slope compared to the scar and non-scar areas, 32%, 4%, and 0.2%, respectively.

CONCLUSIONS

VT stimulation studies can be simulated in silico using a personalized biophysical cardiac model. Myocardial spatial heterogeneity of APD restitution properties and conductivity may help predict the location of crucial entry/exit points of reentrant VT circuits.

Virtual ablation for atrial fibrillation in personalized in-silico three-dimensional left atrial modeling: comparison with clinical catheter ablation

BACKGROUND

Although catheter ablation is an effective rhythm control strategy for atrial fibrillation (AF), empirically-based ablation has a substantial recurrence rate. The purposes of this study were to develop a computational platform for patient-specific virtual AF ablation and to compare the anti-fibrillatory effects of 5 different virtual ablation protocols with empirically chosen clinical ablations.

METHODS

We included 20 patients with AF (65% male, 60.1 ± 10.5 years old, 80% persistent AF [PeAF]) who had undergone empirically-based catheter ablation: circumferential pulmonary vein isolation (CPVI) for paroxysmal AF (PAF) and additional posterior box lesion (L1) and anterior line (L2) for PeAF. Using patient-specific three-dimensional left atrial (LA) geometry, we generated a finite element model and tested the AF termination rate after 5 different virtual ablations: CPVI alone, CPVI + L1, CPVI + L1,2, CPVI with complex fractionated atrial electrogram (CFAE) ablation, and CFAE ablation alone.

RESULTS

1. Virtual CPVI + L1,2 ablation showed the highest AF termination rate in overall patients (55%) and PeAF patients (n = 16, 62.5%). 2. The virtual AF maintenance duration was shortest in the case of virtual CPVI + L1,2 ablation in overall patients (2.19 ± 1.28 vs. 2.91 ± 1.04 s, p = 0.009) and in patients with PeAF (2.05 ± 1.23 vs. 2.93 ± 10.2 s, p = 0.004) compared with other protocols.

CONCLUSION

Virtual AF ablation using personalized in-silico model of LA is feasible. Virtual ablation with CPVI + L1,2 shows the highest antifibrillatory effect, concordant with the empirical ablation protocol in patients with PeAF.

The relationship among complex fractionated electrograms, wavebreak, phase singularity, and local dominant frequency in fibrillation wave-dynamics: a modeling comparison study

Although complex fractionated electrogram (CFE) is known to be a target for catheter ablation of fibrillation, its physiological meaning in fibrillation wave-dynamics remains to be clarified. We evaluated the spatiotemporal relationships among the parameters of fibrillation wave-dynamics by simulation modeling. We generated maps of CFE-cycle length (CFE-CL), local dominant frequency (LDF), wave break (WB), and phase singularity (PS) of fibrillation in 2-dimensional homogeneous bidomain cardiac modeling (1,000 × 1,000 cells ten Tusscher model). We compared spatiotemporal correlations by dichotomizing each maps into 10 × 10 lattice zones. In spatial distribution, WB and PS showed excellent correlation (R = 0.963, P < 0.001). CFE-CL had weak correlations with WB (R = 0.288, P < 0.001), PS (R = 0.313, P < 0.001), and LDF (R = -0.411, P < 0.001). However, LDF did not show correlation with PS or WB. PSs were mostly distributed at the periphery of low CFE-CL area. Virtual ablation (5% of critical mass) of CFE-CL < 100 ms terminated fibrillation at 14.3 sec, and high LDF ablation (5% of critical mass) changed fibrillation to organized tachycardia, respectively. In homogeneous 2D fibrillation modeling, CFE-CL was weakly correlated with WB, PS, and LDF, spatiotemporally. PSs are mostly positioned at the periphery of low CFE-CL areas, and virtual ablation targeting low CFE-CL regions terminated fibrillation successfully.

Effective dynamics of twisted and curved scroll waves using virtual filaments

Scroll waves are three-dimensional excitation patterns that rotate around a central filament curve; they occur in many physical, biological, and chemical systems. We explicitly derive the equations of motion for scroll wave filaments in reaction-diffusion systems with isotropic diffusion up to third order in the filament's twist and curvature. The net drift components define at every instance of time a virtual filament which lies close to the instantaneous filament. Importantly, virtual filaments obey simpler, time-independent laws of motion which we analytically derive here and illustrate with numerical examples. Stability analysis of scroll waves is performed using virtual filaments, showing that filament curvature and twist add as quadratic terms to the nominal filament tension. Applications to oscillating chemical reactions and cardiac tissue are discussed.

The functional role of electrophysiological heterogeneity in the rabbit ventricle during rapid pacing and arrhythmias

Electrophysiological heterogeneity in action potential recordings from healthy intact hearts remains highly variable and, where present, is almost entirely abolished at fast pacing rates. Consequently, the functional importance of intrinsic action potential duration (APD) heterogeneity in healthy ventricles, and particularly its role during rapidly activating reentrant arrhythmias, remain poorly understood. By incorporating both transmural and apicobasal APD heterogeneity within a biventricular rabbit computational model and comparing with an equivalent homogeneous model, we directly investigated the functional importance of intrinsic APD heterogeneity under fast pacing and arrhythmogenic protocols. Although differences in APD were significantly modulated at the tissue level during pacing and further reduced as pacing frequency increased, small differences were still noticeable. Such differences were further marginally accentuated/attenuated via electrotonic effects relative to wavefront propagation directions. The remaining small levels of APD heterogeneity under the fastest pacing frequencies resulted in arrhythmia initiation via heterogeneous conduction block, in contrast to complete block in the homogeneous model. Such induction mechanisms were more evident during premature stimuli at slower paced rhythms where intrinsic heterogeneity remained to a greater degree. During sustained arrhythmias, however, intrinsic heterogeneity made little difference to overall reentrant behavior, either visually, or in terms of duration, metrics quantifying filament/phase singularity dynamics, and global electrocardiogram characteristics. These findings suggest that, despite being important during arrhythmia initiation, intrinsic electrophysiological heterogeneity plays little functional role during rapid pacing and sustained arrhythmia dynamics in the healthy ventricle and thus questions the need to incorporate such detail in computational models when simulating rapid arrhythmias.

The role of fine-scale anatomical structure in the dynamics of reentry in computational models of the rabbit ventricles

Fine-scale anatomical structures in the heart may play an important role in sustaining cardiac arrhythmias. However, the extent of this role and how it may differ between species are not fully understood. In this study we used computational modelling to assess the impact of anatomy upon arrhythmia maintenance in the rabbit ventricles. Specifically, we quantified the dynamics of excitation wavefronts during episodes of simulated tachyarrhythmias and fibrillatory arrhythmias, defined as being respectively characterised by relatively low and high spatio-temporal disorganisation.Two computational models were used: a highly anatomically detailed MR-derived rabbit ventricular model (representing vasculature, endocardial structures) and a simplified equivalent model, constructed from the same MR-data but lacking such fine-scale anatomical features. During tachyarrhythmias, anatomically complex and simplified models showed very similar dynamics; however, during fibrillatory arrhythmias, as activation wavelength decreased, the presence of fine-scale anatomical details appeared to marginally increase disorganisation of wavefronts during arrhythmias in the complex model. Although a small amount of clustering of reentrant rotor centres (filaments) around endocardial structures was witnessed in follow-up analysis (which slightly increased during fibrillation as rotor size decreased), this was significantly less than previously reported in large animals. Importantly, no anchoring of reentrant rotors was visibly identifiable in arrhythmia movies. These differences between tachy- and fibrillatory arrhythmias suggest that the relative size of reentrant rotors with respect to anatomical obstacles governs the influence of fine-scale anatomy in the maintenance of ventricular arrhythmias in the rabbit. In conclusion, our simulations suggest that fine-scale anatomical features play little apparent role in the maintenance of tachyarrhythmias in the rabbit ventricles and, contrary to experimental reports in larger animals, appear to play only a minor role in the maintenance of fibrillatory arrhythmias. These findings also have important implications in optimising the level of detail required in anatomical computational meshes frequently used in arrhythmia investigations.

A computational study of mother rotor VF in the human ventricles

Sudden cardiac death is one of the major causes of death in the industrialized world. It is most often caused by a cardiac arrhythmia called ventricular fibrillation (VF). Despite its large social and economical impact, the mechanisms for VF in the human heart yet remain to be identified. Two of the most frequently discussed mechanisms observed in experiments with animal hearts are the multiple wavelet and mother rotor hypotheses. Most recordings of VF in animal hearts are consistent with the multiple wavelet mechanism. However, in animal hearts, mother rotor fibrillation has also been observed. For both multiple wavelet and mother rotor VF, cardiac heterogeneity plays an important role. Clinical data of action potential restitution measured from the surface of human hearts have been recently published. These in vivo data show a substantial degree of spatial heterogeneity. Using these clinical restitution data, we studied the dynamics of VF in the human heart using a heterogeneous computational model of human ventricles. We hypothesized that this observed heterogeneity can serve as a substrate for mother rotor fibrillation. We found that, based on these data, mother rotor VF can occur in the human heart and that ablation of the mother rotor terminates VF. Furthermore, we found that both mother rotor and multiple wavelet VF can occur in the same heart depending on the initial conditions at the onset of VF. We studied the organization of these two types of VF in terms of filament numbers, excitation periods, and frequency domains. We conclude that mother rotor fibrillation is a possible mechanism in the human heart.

Vortex filament dynamics in computational models of ventricular fibrillation in the heart

In three-dimensional cardiac tissue, the re-entrant waves that sustain ventricular fibrillation rotate around a line of phase singularity or vortex filament. The aim of this study was to investigate how the behavior of these vortex filaments is influenced by membrane kinetics, initial conditions, and tissue geometry in computational models of excitable tissue. A monodomain model of cardiac tissue was used, with kinetics described by a three-variable simplified ionic model (3V-SIM). Two versions of 3V-SIM were used, one with steep action potential duration restitution, and one with reduced excitability. Re-entrant fibrillation was then simulated in three tissue geometries: a cube, a slab, and an anatomically detailed model of rabbit ventricles. Filaments were identified using a phase-based method, and the number, size, origin, and orientation of filaments was tracked throughout each simulation. The main finding of this study is that kinetics, initial conditions, geometry, and anisotropy all affected the number, proliferation, and orientation of vortex filaments in re-entrant fibrillation. An important finding of this study was that the behavior of vortex filaments in simplified slab geometry representing part of the ventricular wall did not necessarily predict behavior in an anatomically detailed model of the rabbit ventricles.

Effect of heterogeneous APD restitution on VF organization in a model of the human ventricles

The onset of ventricular fibrillation (VF) has been associated with steep action potential duration restitution in both clinical and computational studies. Recently, detailed clinical restitution properties in cardiac patients were reported showing a substantial degree of heterogeneity in restitution slopes at the epicardium of the ventricles. The aim of the present study was to investigate the effect of heterogeneous restitution properties in a three-dimensional model of the ventricles using these clinically measured restitution data. We used a realistic model of the human ventricles, including detailed descriptions of cell electrophysiology, ventricular anatomy, and fiber direction anisotropy. We extended this model by mapping the clinically observed epicardial restitution data to our anatomic representation using a diffusion-based algorithm. Restitution properties were then fitted by regionally varying parameters of the electrophysiological model. We studied the effects of restitution heterogeneity on the organization of VF by analyzing filaments and the distributions of excitation periods. We found that the number of filaments and the excitation periods were both dependent on the extent of heterogeneity. An increased level of heterogeneity leads to a greater number of filaments and a broader distribution of excitation periods, thereby increasing the complexity and dynamics of VF. Restitution heterogeneity may play an important role in providing a substrate for cardiac arrhythmias.

A guide to modelling cardiac electrical activity in anatomically detailed ventricles

One of the most recent trends in cardiac electrophysiology is the development of integrative anatomically accurate models of the heart, which include description of cardiac activity from sub-cellular and cellular level to the level of the whole organ. In order to construct this type of model, a researcher needs to collect a wide range of information from books and journal articles on various aspects of biology, physiology, electrophysiology, numerical mathematics and computer programming. The aim of this methodological article is to survey recent developments in integrative modelling of electrical activity in the ventricles of the heart, and to provide a practical guide to the resources and tools that are available for work in this exciting and challenging area.

Organization of ventricular fibrillation in the human heart

Sudden cardiac death is a major cause of death in the industrialized world, claiming approximately 300,000 victims annually in the United States alone. In most cases, sudden cardiac death is caused by ventricular fibrillation (VF). Experimental studies in large animal hearts have shown that the uncoordinated contractions during VF are caused by large numbers of chaotically wandering reentrant waves of electrical activity. However, recent clinical data on VF in the human heart seem to suggest that human VF may have a markedly different organization. Here, we use a detailed model of the human ventricles, including a detailed description of cell electrophysiology, ventricular anatomy, and fiber direction anisotropy, to study the organization of human VF. We show that characteristics of our simulated VF are qualitatively similar to the clinical data. Furthermore, we find that human VF is driven by only approximately 10 reentrant sources and thus is much more organized than VF in animal hearts of comparable size, where VF is driven by approximately 50 sources. We investigate the influence of anisotropy ratio, tissue excitability, and restitution properties on the number of reentrant sources driving VF. We find that the number of rotors depends strongest on minimum action potential duration, a property that differs significantly between human and large animal hearts. Based on these findings, we suggest that the simpler spatial organization of human VF relative to VF in large animal hearts may be caused by differences in minimum action potential duration. Both the simpler spatial organization of human VF and its suggested cause may have important implications for treating and preventing this dangerous arrhythmia in humans.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram (ECG) (mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573-89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (I(CaL)), fast sodium current (I(Na)), transient outward potassium current (I(to)), rapidly and slowly delayed rectifier potassium currents (I(Kr) and I(Ks)) and inward rectifier potassium current (I(Ki)). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of I(Ks) and I(to). The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart-torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Whole heart action potential duration restitution properties in cardiac patients: a combined clinical and modelling study

Steep action potential duration (APD) restitution has been shown to facilitate wavebreak and ventricular fibrillation. The global APD restitution properties in cardiac patients are unknown. We report a combined clinical electrophysiology and computer modelling study to: (1) determine global APD restitution properties in cardiac patients; and (2) examine the interaction of the observed APD restitution with known arrhythmia mechanisms. In 14 patients aged 52-85 years undergoing routine cardiac surgery, 256 electrode epicardial mapping was performed. Activation-recovery intervals (ARI; a surrogate for APD) were recorded over the entire ventricular surface. Mono-exponential restitution curves were constructed for each electrode site using a standard S1-S2 pacing protocol. The median maximum restitution slope was 0.91, with 27% of all electrode sites with slopes<0.5, 29% between 0.5 and 1.0, and 20% between 1.0 and 1.5. Eleven per cent of restitution curves maintained slope>1 over a range of diastolic intervals of at least 30 ms; and 0.3% for at least 50 ms. Activation-recovery interval restitution was spatially heterogeneous, showing regional organization with multiple discrete areas of steep and shallow slope. We used a simplified computer model of 2-D cardiac tissue to investigate how heterogeneous APD restitution can influence vulnerability to, and stability of re-entry. Our model showed that heterogeneity of restitution can act as a potent arrhythmogenic substrate, as well as influencing the stability of re-entrant arrhythmias. Global epicardial mapping in humans showed that APD restitution slopes were organized into regions of shallow and steep slopes. This heterogeneous organization of restitution may provide a substrate for arrhythmia.

Regional differences in APD restitution can initiate wavebreak and re-entry in cardiac tissue: a computational study

Background

Regional differences in action potential duration (APD) restitution in the heart favour arrhythmias, but the mechanism is not well understood.

Methods

We simulated a 150 × 150 mm 2D sheet of cardiac ventricular tissue using a simplified computational model. We investigated wavebreak and re-entry initiated by an S1S2S3 stimulus protocol in tissue sheets with two regions, each with different APD restitution. The two regions had a different APD at short diastolic interval (DI), but similar APD at long DI. Simulations were performed twice; once with both regions having steep (slope > 1), and once with both regions having flat (slope < 1) APD restitution.

Results

Wavebreak and re-entry were readily initiated using the S1S2S3 protocol in tissue sheets with two regions having different APD restitution properties. Initiation occurred irrespective of whether the APD restitution slopes were steep or flat. With steep APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms with S1S2 of 250 ms, to 75 ms (S1S2 180 ms). With flat APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms (S1S2 250 ms), to 21 ms (S1S2 340 ms) and then 11 ms (S1S2 400 ms).

Conclusion

Regional differences in APD restitution are an arrhythmogenic substrate that can be concealed at normal heart rates. A premature stimulus produces regional differences in repolarisation, and a further premature stimulus can then result in wavebreak and initiate re-entry. This mechanism for initiating re-entry is independent of the steepness of the APD restitution curve.

Propagation of normal beats and re-entry in a computational model of ventricular cardiac tissue with regional differences in action potential shape and duration

There is substantial experimental evidence from studies using both intact tissue and isolated single cells to support the existence of different cell types within the ventricular wall of the heart, each possessing different electrical properties. However other studies have failed to find these differences, and instead support the idea that electrical coupling in vivo between regions with different cell types smoothes out differences in action potential shape and duration. In this study we have used a computational model of electrical activation in heterogenous 2D and 3D cardiac tissue to investigate the propagation of both normal beats and arrhythmias. We used the Luo-Rudy dynamic model for guinea pig ventricular cells, with simplified Ca2+ handling and transmural heterogeneity in IKs and Ito. With normal cell-to-cell coupling, a layer of M cells was not necessary for the formation of an upright T wave in the simulated electrocardiogram, and the amplitude and configuration of the T wave was not greatly affected by the thickness and configuration of the M cell layer. Transmural gradients in repolarisation pushed re-entrant waves with an intramural filament towards either the base or the apex of the ventricles, and caused transient break up of re-entry with a transmural filament.

Spatiotemporal correlation between phase singularities and wavebreaks during ventricular fibrillation

Phase Singularity and Wavebreak.

INTRODUCTION

Phase maps and the detection of phase singularities (PSs) have become a well-developed method for characterizing the organization of ventricular fibrillation (VF). How precisely PS colocalizes with wavebreak (WB) during VF, however, is unknown.

METHODS AND RESULTS

We performed optical mapping of 27 episodes of VF in nine Langendorff-perfused rabbit hearts. A WB is a point where the activation wavefront and the repolarization waveback meet. A PS is a site where its phase is ambiguous and its neighboring pixels exhibit a continuous phase progression from -pi to +pi. The correlation coefficient between the number of WBs and PSs was 0.78 +/- 0.09 for each heart and 0.81 for all VF episodes (P < 0.001), indicating a significant temporal correlation. We then superimposed the WBs and PSs for every 100 frames of each episode. These maps showed a high degree of spatial colocalization. To quantify spatial colocalization, the spatial shifts between the cumulative maps of WBs and PSs in corresponding frames were calculated by automatic alignment to obtain maximum overlap between these two maps. The spatial shifts were 0.04 +/- 0.31 mm on the x-axis and 0.06 +/- 0.27 mm on the y-axis over a 20 x 20 mm2 mapped field, indicating highly significant spatial correlation.

CONCLUSION

Phase mapping provides a convenient and robust approach to quantitatively describe wave propagation and organization during VF. The close spatiotemporal correlation between PSs and WBs establishes that PSs are a valid alternate representation of WB during VF and further validated the use of phase mapping in the study of VF dynamics.

Filament behavior in a computational model of ventricular fibrillation in the canine heart

The aim of this paper was to quantify the behavior of filaments in a computational model of re-entrant ventricular fibrillation. We simulated cardiac activation in an anisotropic monodomain with excitation described by the Fenton-Karma model with Beeler-Reuter restitution, and geometry by the Auckland canine ventricle. We initiated re-entry in the left and right ventricular free walls, as well as the septum. The number of filaments increased during the first 1.5 s before reaching a plateau with a mean value of about 36 in each simulation. Most re-entrant filaments were between 10 and 20 mm long. The proportion of filaments touching the epicardial surface was 65%, but most of these were visible for much less than one period of re-entry. This paper shows that useful information about filament dynamics can be gleaned from models of fibrillation in complex geometries, and suggests that the interplay of filament creation and destruction may offer a target for antifibrillatory therapy.

Effect of regional differences in cardiac cellular electrophysiology on the stability of ventricular arrhythmias: a computational study

Re-entry is an important mechanism of cardiac arrhythmias. During re-entry a wave of electrical activation repeatedly propagates into recovered tissue, rotating around a rod-like filament. Breakdown of a single re-entrant wave into multiple waves is believed to underlie the transition from ventricular tachycardia to ventricular fibrillation. Several mechanisms of breakup have been identified including the effect of anisotropic conduction in the ventricular wall. Cells in the inner and outer layers of the ventricular wall have different action potential durations (APD), and support re-entrant waves with different periods. The aim of this study was to use a computational approach to study twisting and breakdown in a transmural re-entrant wave spanning these regions, and examine the relative role of this effect and anisotropic conduction. We used a simplified model of action potential conduction in the ventricular wall that we modified so that it supported stable re-entry in an anisotropic model with uniform APD. We first examined the effect of regional differences on breakdown in an isotropic model with transmural differences in APD, and found that twisting of the re-entrant filament resulted in buckling and breakdown during the second cycle of re-entry. We found that breakdown was amplified in the anisotropic model, resulting in complex activation in the region of longest APD. This study shows that regional differences in cardiac electrophysiology are a potentially important mechanism for destabilizing re-entry and may act synergistically with other mechanisms to mediate the transition from ventricular tachycardia to ventricular fibrillation.

Computational framework for simulating the mechanisms and ECG of re-entrant ventricular fibrillation

Ventricular arrhythmias remain an important cause of morbidity and mortality in the Western world. Although the underlying mechanisms of these arrhythmias can be studied experimentally, these investigations are in general limited to mapping electrical activity on the heart surface. Computational models of action potential propagation offer a potentially powerful way to study electrical activation and arrhythmias, but current models are not easy to link to the clinical environment. In this paper, we describe a framework for computing action potential propagation in which the geometry, electrophysiology and regional properties of ventricular myocardium can be specified so that, for example, different models for cardiac cellular electrophysiology can be used. We have computed action potential propagation during both normal beats and re-entry in an anatomically accurate model of ventricular geometry. By computing the resultant electric current flow in the torso we have also generated simulated ECG signals that result from specific activation patterns in the ventricular model. Models can be powerful tools for explaining observations, and this approach is able to provide a direct link between the different configurations of re-entry observed in computational and experimental studies, and the ECG signals observed in patients.

Enhanced self-termination of re-entrant arrhythmias as a pharmacological strategy for antiarrhythmic action

Ventricular tachycardia and fibrillation are potentially lethal cardiac arrhythmias generated by high frequency, irregular spatio-temporal electrical activity. Re-entrant propagation has been demonstrated as a mechanism generating these arrhythmias in computational and in vitro animal models of these arrhythmias. Re-entry can be idealised in homogenous isotropic virtual cardiac tissues as spiral and scroll wave solutions of reaction-diffusion equations. A spiral wave in a bounded medium can be terminated if its core reaches a boundary. Ventricular tachyarrhythmias in patients are sometimes observed to spontaneously self-terminate. One possible mechanism for self-termination of a spiral wave is meander of its core to an inexcitable boundary. We have previously proposed the hypothesis that the spatial extent of meander of a re-entrant wave in the heart can be directly related to its probability of self-termination, and so inversely related to its lethality. Meander in two-dimensional virtual ventricular tissues based on the Oxsoft family of cell models, with membrane excitation parameters simulating the inherited long Q-T syndromes has been shown to be consistent with this hypothesis: the largest meander is seen in the syndrome with the lowest probability of death per arrhythmic episode. Here we extend our previous results to virtual tissues based on the Luo-Rudy family of models. Consistent with our hypothesis, for both families of models, whose different ionic mechanisms produce different patterns of meander, the LQT virtual tissue with the larger meander simulates the syndrome with the lower probability of death per episode. Further, we search the parameter space of the repolarizing currents to find their conductance parameter values that give increased meander of spiral waves. These parameters may provide targets for antiarrhythmic drugs designed to act by increasing the likelihood of self-termination of re-entrant arrhythmias. (c) 2002 American Institute of Physics.

Modified ionic models of cardiac tissue for efficient large scale computations

Recirculation of excitation, or re-entry, is one of the most important mechanisms of life-threatening cardiac arrhythmias and fibrillation. Modelling these phenomena requires large scale computations in two and three-dimensional slabs of cardiac tissue. Because of computational constraints, most of the studies use simplified (non-ionic) models of cardiac tissue, which are electrophysiologically less accurate than the detailed ionic models. In this paper, we propose a method to modify ionic models of cardiac tissue into an intermediate class of models, which are almost as efficient for computations as simplified models, and retain most of the properties of the original ionic models, such as the shape of the action potential, the restitution of action potential duration and of the conduction velocity, as well as unchanged description of most of the ionic currents.

Dynamics and interaction of filaments in a computational model of re-entrant ventricular fibrillation

Ventricular fibrillation (VF) is a lethal cardiac arrhythmia. Re-entry, in which action potential wavefronts rotate around filaments, is believed to sustain VF. In this study we used a computational model of multiple wavelet fibrillation in the thin-walled right ventricle (10 mm thick) and the thicker walled left ventricle (16 mm thick) to investigate the effect of tissue thickness and initiation protocol on re-entry, and to examine whether filament dynamics and interaction in the model could explain why re-entry is both rarely observed and short-lived in experimental studies that map electrical activation on the heart surface. We found (i) that the density of filaments, the proportion of transmural filaments and the proportion of filaments visible on the model surface were all higher in the 10 mm simulation, (ii) that the initiation protocol influences the rate of filament breakdown but not the number of filaments present after 1 s, and (iii) that although many filaments are visible on the surface of the model, the majority are visible for less than one rotation. This study shows that tissue thickness, geometry and initiation protocol influence electrical activation during VF, and that the rapid motion and interaction of filaments result in transient appearance of surface re-entry.