Use of topological charge to determine filament location and dynamics in a numerical model of scroll wave activity

Human atrial fibrillation and genetic defects in transient outward currents: mechanistic insights from multi-scale computational models

Previous studies have linked dysfunctional I_to arising from mutations to KCND3-encoded Kv4.3 and KCND2-encoded Kv4.2 to atrial fibrillation. Using computational models, this study aimed to investigate the mechanisms underlying pro-arrhythmic effects of the gain-of-function Kv4.3 (T361S, A545P) and Kv4.2 (S447R) mutations. Wild-type and mutant I_to formulations were developed from and validated against experimental data and incorporated into the Colman et al. model of human atrial cells. Single-cell models were incorporated into one- (1D) and two-dimensional (2D) models of atrial tissue, and a three-dimensional (3D) realistic model of the human atria. The three gain-of-function mutations had similar, albeit quantitatively different, effects: shortening of the action potential duration; lowering the plateau membrane potential, abbreviating the effective refractory period (ERP) and the wavelength (WL) of atrial excitation at the tissue level. Restitution curves for the WL, the ERP and the conduction velocity were leftward shifted, facilitating the conduction of atrial excitation waves at high excitation rates. The mutations also increased lifespan and stationarity of re-entry in both 2D and 3D simulations, which further highlighted a mutation-induced increase in spatial dispersion of repolarization. Collectively, these changes account for pro-arrhythmic effects of these Kv4.3 and Kv4.2 mutations in facilitating AF. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Electrocardiographic imaging in the atria

The inverse problem of electrocardiography or electrocardiographic imaging (ECGI) is a technique for reconstructing electrical information about cardiac surfaces from noninvasive or non-contact recordings. ECGI has been used to characterize atrial and ventricular arrhythmias. Although it is a technology with years of progress, its development to characterize atrial arrhythmias is challenging. Complications can arise when trying to describe the atrial mechanisms that lead to abnormal propagation patterns, premature or tachycardic beats, and reentrant arrhythmias. This review addresses the various ECGI methodologies, regularization methods, and post-processing techniques used in the atria, as well as the context in which they are used. The current advantages and limitations of ECGI in the fields of research and clinical diagnosis of atrial arrhythmias are outlined. In addition, areas where ECGI efforts should be concentrated to address the associated unsatisfied needs from the atrial perspective are discussed.

AF driver detection in pulmonary vein area by electropcardiographic imaging: Relation with a favorable outcome of pulmonary vein isolation

Pulmonary vein isolation (PVI) is the most successful treatment for atrial fibrillation (AF) nowadays. However, not all AF patients benefit from PVI. In this study, we evaluate the use of ECGI to identify reentries and relate rotor density in the pulmonary vein (PV) area as an indicator of PVI outcome. Rotor maps were computed in a set of 29 AF patients using a new rotor detection algorithm. The relationship between the distribution of reentrant activity and the clinical outcome after PVI was studied. The number of rotors and proportion of PSs in different atrial regions were computed and compared retrospectively in two groups of patients: patients that remained in sinus rhythm 6 months after PVI and patients with arrhythmia recurrence. The total number of rotors obtained was higher in patients returning to arrhythmia after the ablation (4.31 ± 2.77 vs. 3.58 ± 2.67%, p = 0.018). However, a significantly higher concentration of PSs in the pulmonary veins was found in patients that remained in sinus rhythm (10.20 ± 12.40% vs. 5.19 ± 9.13%, p = 0.011) 6 months after PVI. The results obtained show a direct relationship between the expected AF mechanism and the electrophysiological parameters provided by ECGI, suggesting that this technology offers relevant information to predict the clinical outcome after PVI in AF patients.

Topological charge-density-vector method of identifying filaments of scroll waves

Scroll waves have been found in a variety of three-dimensional excitable media, including physical, chemical, and biological origins. Scroll waves in cardiac tissue are of particular significance as they underlie ventricular fibrillation that can cause sudden death. The behavior of a scroll wave is characterized by a line of phase singularity at its organizing center, known as a filament. A thorough investigation into the filament dynamics is the key to further exploration of the general theory of scroll waves in excitable media and the mechanisms of ventricular fibrillation. In this paper, we propose a method to identify filaments of scroll waves in excitable media. From the definition of the topological charge of filaments, we obtain the discrete expression of the topological charge-density vector, which is useful in calculating the topological charge vectors at each grid in the space directly. The set of starting points of these topological charge vectors represents a set of phase singularities, thereby forming a line of phase singularity, that is, a filament of a scroll wave.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Understanding the origins of the basic equations of statistical fibrillatory dynamics

The mechanisms governing cardiac fibrillation remain unclear; however, it most likely represents a form of spatiotemporal chaos with conservative system dynamics. Renewal theory has recently been suggested as a statistical formulation with governing equations to quantify the formation and destruction of wavelets and rotors in fibrillatory dynamics. In this perspective Review, we aim to explain the origin of the renewal theory paradigm in spatiotemporal chaos. The ergodic nature of pattern formation in spatiotemporal chaos is demonstrated through the use of three chaotic systems: two classical systems and a simulation of cardiac fibrillation. The logistic map and the baker's transformation are used to demonstrate how the apparently random appearance of patterns in classical chaotic systems has macroscopic parameters that are predictable in a statistical sense. We demonstrate that the renewal theory approach developed for cardiac fibrillation statistically predicts pattern formation in these classical chaotic systems. Renewal theory provides governing equations to describe the apparently random formation and destruction of wavelets and rotors in atrial fibrillation (AF) and ventricular fibrillation (VF). This statistical framework for fibrillatory dynamics provides a holistic understanding of observed rotor and wavelet dynamics and is of conceptual significance in informing the clinical and mechanistic research of the rotor and multiple-wavelet mechanisms of AF and VF.

Topological charge-density method of identifying phase singularities in cardiac fibrillation

Spiral waves represent the key motifs of typical self-sustained dynamical patterns in excitable systems such as cardiac tissue. The motion of phase singularities (PSs) that lies at the center of spiral waves captures many qualitative and, in some cases, quantitative features of their complex dynamics. Recent clinical studies suggested that ablating the tissue at PS locations may cure atrial fibrillation. Here, we propose a different method to determine the location of PSs. Starting from the definition of the topological charge of spiral waves, we obtain the expression of the topological charge density in a discrete case. With this expression, we can calculate the topological charge at each grid in the space directly, so as to accurately identify the position of PSs.

Spatial phase discontinuity at the center of moving cardiac spiral waves

BACKGROUND

Precise analysis of cardiac spiral wave (SW) dynamics is essential for effective arrhythmia treatment. Although the phase singularity (PS) point in the spatial phase map has been used to determine the cardiac SW center for decades, quantitative detection algorithms that assume PS as a point fail to trace complex and rapid PS dynamics. Through a detailed analysis of numerical simulations, we examined our hypothesis that a boundary of spatial phase discontinuity induced by a focal conduction block exists around the moving SW center in the phase map.

METHOD

In a numerical simulation model of a 2D cardiac sheet, three different types of SWs (short wavelength; long wavelength; and low excitability) were induced by regulating ion channels. Discontinuities of all boundaries among adjacent cells at each instance were evaluated by calculating the phase bipolarity (PB). The total amount of phase transition (PTA) in each cell during the study period was evaluated.

RESULTS

Pivoting, drifting, and shifting SWs were observed in the short-wavelength, low-excitability, and long-wavelength models, respectively. For both the drifting and shifting cases, long high-PB edges were observed on the SW trajectories. In all cases, the conduction block (CB) was observed at the same boundaries. These were also identical to the boundaries in the PTA maps.

CONCLUSIONS

The analysis of the simulations revealed that the conduction block at the center of a moving SW induces discontinuous boundaries in spatial phase maps that represent a more appropriate model of the SW center than the PS point.

M/M/Infinity Birth-Death Processes - A Quantitative Representational Framework to Summarize and Explain Phase Singularity and Wavelet Dynamics in Atrial Fibrillation

Rationale

A quantitative framework to summarize and explain the quasi-stationary population dynamics of unstable phase singularities (PS) and wavelets in human atrial fibrillation (AF) is at present lacking. Building on recent evidence showing that the formation and destruction of PS and wavelets in AF can be represented as renewal processes, we sought to establish such a quantitative framework, which could also potentially provide insight into the mechanisms of spontaneous AF termination.

Objectives

Here, we hypothesized that the observed number of PS or wavelets in AF could be governed by a common set of renewal rate constants λ _f (for PS or wavelet formation) and λ _d (PS or wavelet destruction), with steady-state population dynamics modeled as an M/M/∞ birth-death process. We further hypothesized that changes to the M/M/∞ birth-death matrix would explain spontaneous AF termination.

Methods and Results

AF was studied in in a multimodality, multispecies study in humans, animal experimental models (rats and sheep) and Ramirez-Nattel-Courtemanche model computer simulations. We demonstrated: (i) that λ _f and λ _d can be combined in a Markov M/M/∞ process to accurately model the observed average number and population distribution of PS and wavelets in all systems at different scales of mapping; and (ii) that slowing of the rate constants λ _f and λ _d is associated with slower mixing rates of the M/M/∞ birth-death matrix, providing an explanation for spontaneous AF termination.

Conclusion

M/M/∞ birth-death processes provide an accurate quantitative representational architecture to characterize PS and wavelet population dynamics in AF, by providing governing equations to understand the regeneration of PS and wavelets during sustained AF, as well as providing insight into the mechanism of spontaneous AF termination.

EEG microstate periodicity explained by rotating phase patterns of resting-state alpha oscillations

Spatio-temporal patterns in electroencephalography (EEG) can be described by microstate analysis, a discrete approximation of the continuous electric field patterns produced by the cerebral cortex. Resting-state EEG microstates are largely determined by alpha frequencies (8-12 Hz) and we recently demonstrated that microstates occur periodically with twice the alpha frequency. To understand the origin of microstate periodicity, we analyzed the analytic amplitude and the analytic phase of resting-state alpha oscillations independently. In continuous EEG data we found rotating phase patterns organized around a small number of phase singularities which varied in number and location. The spatial rotation of phase patterns occurred with the underlying alpha frequency. Phase rotors coincided with periodic microstate motifs involving the four canonical microstate maps. The analytic amplitude showed no oscillatory behaviour and was almost static across time intervals of 1-2 alpha cycles, resulting in the global pattern of a standing wave. In n=23 healthy adults, time-lagged mutual information analysis of microstate sequences derived from amplitude and phase signals of awake eyes-closed EEG records showed that only the phase component contributed to the periodicity of microstate sequences. Phase sequences showed mutual information peaks at multiples of 50 ms and the group average had a main peak at 100 ms (10 Hz), whereas amplitude sequences had a slow and monotonous information decay. This result was confirmed by an independent approach combining temporal principal component analysis (tPCA) and autocorrelation analysis. We reproduced our observations in a generic model of EEG oscillations composed of coupled non-linear oscillators (Stuart-Landau model). Phase-amplitude dynamics similar to experimental EEG occurred when the oscillators underwent a supercritical Hopf bifurcation, a common feature of many computational models of the alpha rhythm. These findings explain our previous description of periodic microstate recurrence and its relation to the time scale of alpha oscillations. Moreover, our results corroborate the predictions of computational models and connect experimentally observed EEG patterns to properties of critical oscillator networks.

Standardizing Single-Frame Phase Singularity Identification Algorithms and Parameters in Phase Mapping During Human Atrial Fibrillation

PURPOSE

Recent investigations failed to reproduce the positive rotor-guided ablation outcomes shown by initial studies for treating persistent atrial fibrillation (persAF). Phase singularity (PS) is an important feature for AF driver detection, but algorithms for automated PS identification differ. We aim to investigate the performance of four different techniques for automated PS detection.

METHODS

2048-channel virtual electrogram (VEGM) and electrocardiogram signals were collected for 30 s from 10 patients undergoing persAF ablation. QRST-subtraction was performed and VEGMs were processed using sinusoidal wavelet reconstruction. The phase was obtained using Hilbert transform. PSs were detected using four algorithms: (1) 2D image processing based and neighbor-indexing algorithm; (2) 3D neighbor-indexing algorithm; (3) 2D kernel convolutional algorithm estimating topological charge; (4) topological charge estimation on 3D mesh. PS annotations were compared using the structural similarity index (SSIM) and Pearson's correlation coefficient (CORR). Optimized parameters to improve detection accuracy were found for all four algorithms using F β score and 10-fold cross-validation compared with manual annotation. Local clustering with density-based spatial clustering of applications with noise (DBSCAN) was proposed to improve algorithms 3 and 4.

RESULTS

The PS density maps created by each algorithm with default parameters were poorly correlated. Phase gradient threshold and search radius (or kernels) were shown to affect PS detections. The processing times for the algorithms were significantly different (p < 0.0001). The F β scores for algorithms 1, 2, 3, 3 + DBSCAN, 4 and 4 + DBSCAN were 0.547, 0.645, 0.742, 0.828, 0.656, and 0.831. Algorithm 4 + DBSCAN achieved the best classification performance with acceptable processing time (2.0 ± 0.3 s).

CONCLUSION

AF driver identification is dependent on the PS detection algorithms and their parameters, which could explain some of the inconsistencies in rotor-guided ablation outcomes in different studies. For 3D triangulated meshes, algorithm 4 + DBSCAN with optimal parameters was the best solution for real-time, automated PS detection due to accuracy and speed. Similarly, algorithm 3 + DBSCAN with optimal parameters is preferred for uniform 2D meshes. Such algorithms - and parameters - should be preferred in future clinical studies for identifying AF drivers and minimizing methodological heterogeneities. This would facilitate comparisons in rotor-guided ablation outcomes in future works.

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.

Investigation of the Role of Myocyte Orientations in Cardiac Arrhythmia Using Image-Based Models

Cardiac electrical excitation-propagation is influenced by myocyte orientations (cellular organization). Quantitatively understanding this relationship presents a significant research challenge, especially during arrhythmias in which excitation patterns become complex. Tissue-scale simulations of cardiac electrophysiology, incorporating both dynamic action potential behavior and image-based myocardial architecture, provide an approach to investigate three-dimensional (3D) propagation of excitation waves in the heart. In this study, we aimed to assess the importance of natural variation in myocyte orientations on cardiac arrhythmogenesis using 3D tissue electrophysiology simulations. Three anatomical models (i.e., describing myocyte orientations) of healthy rat ventricles—obtained using diffusion tensor imaging at 100 μm resolution—were registered to a single biventricular geometry (i.e., a single cardiac shape), in which the myocyte orientations could be represented by each of the diffusion tensor imaging data sets or by an idealized rule-based description. The Fenton-Karma cellular excitation model was modified to reproduce rat ventricular action potential duration restitution to create reaction-diffusion cardiac electrophysiology models. Over 250 3D simulations were performed to investigate the effects of myocyte orientations on the following: 1) ventricular activation, 2) location-dependent arrhythmia induction via rapid pacing, and 3) dynamics of re-entry averaged over multiple episodes. It was shown that 1) myocyte orientation differences manifested themselves in local activation times, but the influence on total activation time was small; 2) differences in myocyte orientations could critically affect the inducibility and persistence of arrhythmias for specific stimulus-location/cycle-length combinations; and 3) myocyte orientations alone could be an important determinant of scroll wave break, although no significant differences were observed in averaged arrhythmia dynamics between the four myocyte orientation scenarios considered. Our results show that myocyte orientations are an important determinant of arrhythmia inducibility, persistence, and scroll wave break. These findings suggest that where specificity is desired (for example, when predicting location-dependent, patient-specific arrhythmia inducibility), subject-specific myocyte orientations may be important.

Directed Networks as a Novel Way to Describe and Analyze Cardiac Excitation: Directed Graph Mapping

Networks provide a powerful methodology with applications in a variety of biological, technological and social systems such as analysis of brain data, social networks, internet search engine algorithms, etc. To date, directed networks have not yet been applied to characterize the excitation of the human heart. In clinical practice, cardiac excitation is recorded by multiple discrete electrodes. During (normal) sinus rhythm or during cardiac arrhythmias, successive excitation connects neighboring electrodes, resulting in their own unique directed network. This in theory makes it a perfect fit for directed network analysis. In this study, we applied directed networks to the heart in order to describe and characterize cardiac arrhythmias. Proof-of-principle was established using in-silico and clinical data. We demonstrated that tools used in network theory analysis allow determination of the mechanism and location of certain cardiac arrhythmias. We show that the robustness of this approach can potentially exceed the existing state-of-the art methodology used in clinics. Furthermore, implementation of these techniques in daily practice can improve the accuracy and speed of cardiac arrhythmia analysis. It may also provide novel insights in arrhythmias that are still incompletely understood.

Patient-specific simulations predict efficacy of ablation of interatrial connections for treatment of persistent atrial fibrillation

Aims

Treatments for persistent atrial fibrillation (AF) offer limited efficacy. One potential strategy aims to return the right atrium (RA) to sinus rhythm (SR) by ablating interatrial connections (IAC) to isolate the atria, but there is limited clinical data to evaluate this ablation approach. We aimed to use simulation to evaluate and predict patient-specific suitability for ablation of IAC to treat AF.

Methods and results

Persistent AF was simulated in 12 patient-specific geometries, incorporating electrophysiological heterogeneity and fibres, with IAC at Bachmann’s bundle, the coronary sinus, and fossa ovalis. Simulations were performed to test the effect of left atrial (LA)-to-RA frequency gradient and fibrotic remodelling on IAC ablation efficacy. During AF, we simulated ablation of one, two, or all three IAC, with or without pulmonary vein isolation and determined if this altered or terminated the arrhythmia. For models without structural remodelling, ablating all IAC terminated RA arrhythmia in 83% of cases. Models with the LA-to-RA frequency gradient removed had an increased success rate (100% success). Ablation of IACs is less effective in cases with fibrotic remodelling (interstitial fibrosis 50% success rate; combination remodelling 67%). Mean number of phase singularities in the RA was higher pre-ablation for IAC failure (success 0.6 ± 0.8 vs. failure 3.2 ± 2.5, P < 0.001).

Conclusion

This simulation study predicts that IAC ablation is effective in returning the RA to SR for many cases. Patient-specific modelling approaches have the potential to stratify patients prior to ablation by predicting if drivers are located in the LA or RA. We present a platform for predicting efficacy and informing patient selection for speculative treatments.

In silico Assessment of Pharmacotherapy for Human Atrial Patho-Electrophysiology Associated With hERG-Linked Short QT Syndrome

Short QT syndrome variant 1 (SQT1) arises due to gain-of-function mutations to the human Ether-à-go-go-Related Gene (hERG), which encodes the α subunit of channels carrying rapid delayed rectifier potassium current, I_Kr. In addition to QT interval shortening and ventricular arrhythmias, SQT1 is associated with increased risk of atrial fibrillation (AF), which is often the only clinical presentation. However, the underlying basis of AF and its pharmacological treatment remain incompletely understood in the context of SQT1. In this study, computational modeling was used to investigate mechanisms of human atrial arrhythmogenesis consequent to a SQT1 mutation, as well as pharmacotherapeutic effects of selected class I drugs–disopyramide, quinidine, and propafenone. A Markov chain formulation describing wild type (WT) and N588K-hERG mutant I_Kr was incorporated into a contemporary human atrial action potential (AP) model, which was integrated into one-dimensional (1D) tissue strands, idealized 2D sheets, and a 3D heterogeneous, anatomical human atria model. Multi-channel pharmacological effects of disopyramide, quinidine, and propafenone, including binding kinetics for I_Kr/hERG and sodium current, I_Na, were considered. Heterozygous and homozygous formulations of the N588K-hERG mutation shortened the AP duration (APD) by 53 and 86 ms, respectively, which abbreviated the effective refractory period (ERP) and excitation wavelength in tissue, increasing the lifespan and dominant frequency (DF) of scroll waves in the 3D anatomical human atria. At the concentrations tested in this study, quinidine most effectively prolonged the APD and ERP in the setting of SQT1, followed by disopyramide and propafenone. In 2D simulations, disopyramide and quinidine promoted re-entry termination by increasing the re-entry wavelength, whereas propafenone induced secondary waves which destabilized the re-entrant circuit. In 3D simulations, the DF of re-entry was reduced in a dose-dependent manner for disopyramide and quinidine, and propafenone to a lesser extent. All of the anti-arrhythmic agents promoted pharmacological conversion, most frequently terminating re-entry in the order quinidine > propafenone = disopyramide. Our findings provide further insight into mechanisms of SQT1-related AF and a rational basis for the pursuit of combined I_Kr and I_Na block based pharmacological strategies in the treatment of SQT1-linked AF.

Analytical approaches for myocardial fibrillation signals

Atrial and ventricular fibrillation are complex arrhythmias, and their underlying mechanisms remain widely debated and incompletely understood. This is partly because the electrical signals recorded during myocardial fibrillation are themselves complex and difficult to interpret with simple analytical tools. There are currently a number of analytical approaches to handle fibrillation data. Some of these techniques focus on mapping putative drivers of myocardial fibrillation, such as dominant frequency, organizational index, Shannon entropy and phase mapping. Other techniques focus on mapping the underlying myocardial substrate sustaining fibrillation, such as voltage mapping and complex fractionated electrogram mapping. In this review, we discuss these techniques, their application and their limitations, with reference to our experimental and clinical data. We also describe novel tools including a new algorithm to map microreentrant circuits sustaining fibrillation.

Determinants of new wavefront locations in cholinergic atrial fibrillation

AIMS

Atrial fibrillation (AF) wavefront dynamics are complex and difficult to interpret, contributing to uncertainty about the mechanisms that maintain AF. We aimed to investigate the interplay between rotors, wavelets, and focal sources during fibrillation.

METHODS AND RESULTS

Arrhythmia wavefront dynamics were analysed for four optically mapped canine cholinergic AF preparations. A bilayer computer model was tuned to experimental preparations, and varied to have (i) fibrosis in both layers or the epicardium only, (ii) different spatial acetylcholine distributions, (iii) different intrinsic action potential duration between layers, and (iv) varied interlayer connectivity. Phase singularities (PSs) were identified and tracked over time to identify rotational drivers. New focal wavefronts were identified using phase contours. Phase singularity density and new wavefront locations were calculated during AF. There was a single dominant mechanism for sustaining AF in each of the preparations, either a rotational driver or repetitive new focal wavefronts. High-density PS sites existed preferentially around the pulmonary vein junctions. Three of the four preparations exhibited stable preferential sites of new wavefronts. Computational simulations predict that only a small number of connections are functionally important in sustaining AF, with new wavefront locations determined by the interplay between fibrosis distribution, acetylcholine concentration, and heterogeneity in repolarization within layers.

CONCLUSION

We were able to identify preferential sites of new wavefront initiation and rotational activity, in order to determine the mechanisms sustaining AF. Electrical measurements should be interpreted differently according to whether they are endocardial or epicardial recordings.

Inter-scale information flow as a surrogate for downward causation that maintains spiral waves

A rotor, the rotation center of spiral waves, has been proposed as a causal mechanism to maintain atrial fibrillation (AF) in human. However, our current understanding of the causality between rotors and spiral waves remains incomplete. One approach to improving our understanding is to determine the relationship between rotors and downward causation from the macro-scale collective behavior of spiral waves to the micro-scale behavior of individual components in a cardiac system. This downward causation is quantifiable as inter-scale information flow that can be used as a surrogate for the mechanism that maintains spiral waves. We used a numerical model of a cardiac system and generated a renormalization group with system descriptions at multiple scales. We found that transfer entropy quantified the upward and downward inter-scale information flow between micro- and macro-scale descriptions of the cardiac system with spiral waves. In addition, because the spatial profile of transfer entropy and intrinsic transfer entropy was identical, there were no synergistic effects in the system. Furthermore, inter-scale information flow significantly decreased as the description of the system became more macro-scale. Finally, downward information flow was significantly correlated with the number of rotors, but the higher numbers of rotors were not necessarily associated with higher downward information flow. This finding contradicts the concept that the rotors are the causal mechanism that maintains spiral waves, and may account for the conflicting evidence from clinical studies targeting rotors to eliminate AF.

Effects of early afterdepolarizations on excitation patterns in an accurate model of the human ventricles

Early Afterdepolarizations, EADs, are defined as the reversal of the action potential before completion of the repolarization phase, which can result in ectopic beats. However, the series of mechanisms of EADs leading to these ectopic beats and related cardiac arrhythmias are not well understood. Therefore, we aimed to investigate the influence of this single cell behavior on the whole heart level. For this study we used a modified version of the Ten Tusscher-Panfilov model of human ventricular cells (TP06) which we implemented in a 3D ventricle model including realistic fiber orientations. To increase the likelihood of EAD formation at the single cell level, we reduced the repolarization reserve (RR) by reducing the rapid delayed rectifier Potassium current and raising the L-type Calcium current. Varying these parameters defined a 2D parametric space where different excitation patterns could be classified. Depending on the initial conditions, by either exciting the ventricles with a spiral formation or burst pacing protocol, we found multiple different spatio-temporal excitation patterns. The spiral formation protocol resulted in the categorization of a stable spiral (S), a meandering spiral (MS), a spiral break-up regime (SB), spiral fibrillation type B (B), spiral fibrillation type A (A) and an oscillatory excitation type (O). The last three patterns are a 3D generalization of previously found patterns in 2D. First, the spiral fibrillation type B showed waves determined by a chaotic bi-excitable regime, i.e. mediated by both Sodium and Calcium waves at the same time and in same tissue settings. In the parameter region governed by the B pattern, single cells were able to repolarize completely and different (spiral) waves chaotically burst into each other without finishing a 360 degree rotation. Second, spiral fibrillation type A patterns consisted of multiple small rotating spirals. Single cells failed to repolarize to the resting membrane potential hence prohibiting the Sodium channel gates to recover. Accordingly, we found that Calcium waves mediated these patterns. Third, a further reduction of the RR resulted in a more exotic parameter regime whereby the individual cells behaved independently as oscillators. The patterns arose due to a phase-shift of different oscillators as disconnection of the cells resulted in continuation of the patterns. For all patterns, we computed realistic 9 lead ECGs by including a torso model. The B and A type pattern exposed the behavior of Ventricular Tachycardia (VT). We conclude that EADs at the single cell level can result in different types of cardiac fibrillation at the tissue and 3D ventricle level.

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.

A statistical model of false negative and false positive detection of phase singularities

The complexity of cardiac fibrillation dynamics can be assessed by analyzing the distribution of phase singularities (PSs) observed using mapping systems. Interelectrode distance, however, limits the accuracy of PS detection. To investigate in a theoretical framework the PS false negative and false positive rates in relation to the characteristics of the mapping system and fibrillation dynamics, we propose a statistical model of phase maps with controllable number and locations of PSs. In this model, phase maps are generated from randomly distributed PSs with physiologically-plausible directions of rotation. Noise and distortion of the phase are added. PSs are detected using topological charge contour integrals on regular grids of varying resolutions. Over 100 × 10⁶ realizations of the random field process are used to estimate average false negative and false positive rates using a Monte-Carlo approach. The false detection rates are shown to depend on the average distance between neighboring PSs expressed in units of interelectrode distance, following approximately a power law with exponents in the range of 1.14 to 2 for false negatives and around 2.8 for false positives. In the presence of noise or distortion of phase, false detection rates at high resolution tend to a non-zero noise-dependent lower bound. This model provides an easy-to-implement tool for benchmarking PS detection algorithms over a broad range of configurations with multiple PSs.

Computational assessment of the functional role of sinoatrial node exit pathways in the human heart

AIM

The human right atrium and sinoatrial node (SAN) anatomy is complex. Optical mapping experiments suggest that the SAN is functionally insulated from atrial tissue except at discrete SAN-atrial electrical junctions called SAN exit pathways, SEPs. Additionally, histological imaging suggests the presence of a secondary pacemaker close to the SAN. We hypothesise that a) an insulating border-SEP anatomical configuration is related to SAN arrhythmia; and b) a secondary pacemaker, the paranodal area, is an alternate pacemaker but accentuates tachycardia. A 3D electro-anatomical computational model was used to test these hypotheses.

METHODS

A detailed 3D human SAN electro-anatomical mathematical model was developed based on our previous anatomical reconstruction. Electrical activity was simulated using tissue specific variants of the Fenton-Karma action potential equations. Simulation experiments were designed to deploy this complex electro-anatomical system to assess the roles of border-SEPs and paranodal area by mimicking experimentally observed SAN arrhythmia. Robust and accurate numerical algorithms were implemented for solving the mono domain reaction-diffusion equation implicitly, calculating 3D filament traces, and computing dominant frequency among other quantitative measurements.

RESULTS

A centre to periphery gradient of increasing diffusion was sufficient to permit initiation of pacemaking at the centre of the 3D SAN. Re-entry within the SAN, micro re-entry, was possible by imposing significant SAN fibrosis in the presence of the insulating border. SEPs promoted the micro re-entry to generate more complex SAN-atrial tachycardia. Simulation of macro re-entry, i.e. re-entry around the SAN, was possible by inclusion of atrial fibrosis in the presence of the insulating border. The border shielded the SAN from atrial tachycardia. However, SAN micro-structure intercellular gap junctional coupling and the paranodal area contributed to prolonged atrial fibrillation. Finally, the micro-structure was found to be sufficient to explain shifts of leading pacemaker site location.

CONCLUSIONS

The simulations establish a relationship between anatomy and SAN electrical function. Microstructure, in the form of intercellular gap junction coupling, was found to regulate SAN function and arrhythmia.

Spatial Resolution Requirements for Accurate Identification of Drivers of Atrial Fibrillation

Supplemental Digital Content is available in the text.

Background—

Recent studies have demonstrated conflicting mechanisms underlying atrial fibrillation (AF), with the spatial resolution of data often cited as a potential reason for the disagreement. The purpose of this study was to investigate whether the variation in spatial resolution of mapping may lead to misinterpretation of the underlying mechanism in persistent AF.

Methods and Results—

Simulations of rotors and focal sources were performed to estimate the minimum number of recording points required to correctly identify the underlying AF mechanism. The effects of different data types (action potentials and unipolar or bipolar electrograms) and rotor stability on resolution requirements were investigated. We also determined the ability of clinically used endocardial catheters to identify AF mechanisms using clinically recorded and simulated data. The spatial resolution required for correct identification of rotors and focal sources is a linear function of spatial wavelength (the distance between wavefronts) of the arrhythmia. Rotor localization errors are larger for electrogram data than for action potential data. Stationary rotors are more reliably identified compared with meandering trajectories, for any given spatial resolution. All clinical high-resolution multipolar catheters are of sufficient resolution to accurately detect and track rotors when placed over the rotor core although the low-resolution basket catheter is prone to false detections and may incorrectly identify rotors that are not present.

Conclusions—

The spatial resolution of AF data can significantly affect the interpretation of the underlying AF mechanism. Therefore, the interpretation of human AF data must be taken in the context of the spatial resolution of the recordings.

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in I_K1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on I_K1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of I_K1 and I_Kr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas I_Kr block alone may be a safer alternative for the E299V mutation. Combined inhibition of I_Kr and I_Kur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, and is characterised by complex and irregular electrical activation of the upper chambers of the heart. One rare, genetic condition associated with increased risk of AF is the short QT syndrome (SQTS), which is caused by mutations in genes involved in normal electrical function of the heart. Underlying mechanisms by which SQTS-related gene mutations facilitate development of arrhythmias in the human atria are not well understood. In this study, sophisticated computer models representing ‘virtual’ human atria, incorporating detailed electrophysiological data at the ‘ion channel’ protein level into both idealised and realistic multi-scale tissue geometries, were used to dissect mechanisms by which two mutations in the KCNJ2 gene responsible for SQTS variant 3 (SQT3) promote initiation and sustenance of arrhythmias. It was found that the D172N and E299V mutations to KCNJ2 accelerated the repolarisation process at the cellular level through distinct mechanisms. This, along with the way the mutations affected heterogeneity in electrical behaviour at the organ level, mediated stability of arrhythmias and response to simulated ion channel block. This study improves understanding of mechanisms underlying increased AF risk associated with D172N and E299V KCNJ2 mutations, and outlines potential therapeutic strategies.

Rotor Tracking Using Phase of Electrograms Recorded During Atrial Fibrillation

Extracellular electrograms recorded during atrial fibrillation (AF) are challenging to interpret due to the inherent beat-to-beat variability in amplitude and duration. Phase mapping represents these voltage signals in terms of relative position within the cycle, and has been widely applied to action potential and unipolar electrogram data of myocardial fibrillation. To date, however, it has not been applied to bipolar recordings, which are commonly acquired clinically. The purpose of this study is to present a novel algorithm for calculating phase from both unipolar and bipolar electrograms recorded during AF. A sequence of signal filters and processing steps are used to calculate phase from simulated, experimental, and clinical, unipolar and bipolar electrograms. The algorithm is validated against action potential phase using simulated data (trajectory centre error <0.8 mm); between experimental multi-electrode array unipolar and bipolar phase; and for wavefront identification in clinical atrial tachycardia. For clinical AF, similar rotational content (R 2 = 0.79) and propagation maps (median correlation 0.73) were measured using either unipolar or bipolar recordings. The algorithm is robust, uses standard signal processing techniques, and accurately quantifies AF wavefronts and sources. Identifying critical sources, such as rotors, in AF, may allow for more accurate targeting of ablation therapy and improved patient outcomes.

A technique for visualising three-dimensional left atrial cardiac activation data in two dimensions with minimal distance distortion

Electro-anatomic mapping and medical imaging systems, used during clinical procedures for treatment of atrial arrhythmias, frequently record and display measurements on an anatomical surface of the left atrium. As such, obtaining a complete picture of activation necessitates simultaneous views from multiple angles. In addition, post-processing of three-dimensional surface data is challenging, since algorithms are typically applicable to planar or volumetric data. We applied a surface flattening methodology to medical imaging data and electro-anatomic mapping data to generate a two-dimensional representation that best preserves distances, since the calculation of many clinically relevant metrics, including conduction velocity and rotor trajectory identification require an accurate representation of distance. Distance distortions were small and improved upon exclusion of the pulmonary veins. The technique is demonstrated using maps of local activation time, based on clinical data, and plotting rotor-core trajectories, using simulated data.

Models of ventricular arrhythmia mechanisms

The mechanisms that initiate and sustain ventricular arrhythmias in the human heart are clinically important, but hard to study experimentally. In this study, a monodomain model of electrical activation was used to examine how dynamics of electrophysiology at the cell scale influence the surface activation patterns of VF at the tissue scale. Cellular electrophysiology was described with two variants of a phenomenological model of the human ventricular epicardial action potential. The tissue geometry was an 8.0 × 8.0 × 1.2 cm 3D tissue slab with axially symmetric anisotropy. In both cases an initial re-entrant wave fragmented into multiple wavelets of activation. The model variant with steep action potential duration restitution produced much more complex activation, with a greater average number of filaments (13.79) than the variant with less steep restitution (3.08). More complex activation was associated with proportionally fewer transmural filaments, and so the average number of epicardial wavefronts and phase singularities per filament was lower. The average number of epicardial phase singularities and wavefronts for the model variant with less steep restitution were consistent with experimental observations in the human heart. This study shows that small changes in cell scale dynamics can have a large influence on the complexity of re-entrant activation in simulated 3D tissue, as well as on the features observed on the epicardial surface.

Analysis of cardiac fibrillation using phase mapping

The sequence of myocardial electrical activation during fibrillation is complex and changes with each cycle. Phase analysis represents the electrical activation-recovery process as an angle. Lines of equal phase converge at a phase singularity at the center of rotation of a reentrant wave, and the identification of reentry and tracking of reentrant wavefronts can be automated. We examine the basic ideas behind phase analysis. With the exciting prospect of using phase analysis of atrial electrograms to guide ablation in the human heart, we highlight several recent developments in preprocessing electrograms so that phase can be estimated reliably.

Human ventricular fibrillation during global ischemia and reperfusion: paradoxical changes in activation rate and wavefront complexity

BACKGROUND

Ischemic ventricular fibrillation in experimental models has been shown to progress through a series of stages. Progression of ischemic VF in the in vivo human heart has not been determined.

METHODS AND RESULTS

We studied 10 patients undergoing cardiac surgery. Ventricular fibrillation was induced by burst pacing. After 30 seconds, global myocardial ischemia was induced by aortic cross-clamp and maintained for 2.5 minutes, followed by coronary reflow. Epicardial activity was sampled (1 kHz) with a sock that contained 256 unipolar contact electrodes. Dominant frequencies were calculated with a fast Fourier transform with a moving window. The locations of phase singularities and activation wavefronts were identified at 10-ms intervals. Preischemic (perfused) ventricular fibrillation was maintained by a disorganized mix of large and small wavefronts. During global myocardial ischemia, mean dominant frequencies decreased from 6.4 to 4.7 Hz at a rate of -0.011±0.002 Hz s(-1) (P<0.001) and then increased rapidly to 7.4 Hz within 30 seconds of reflow. In contrast, the average number of epicardial phase singularities increased during ischemia from 7.7 to 9.7 at a rate of 0.013±0.005 phase singularities per second (P<0.01) and remained unchanged during reflow, at 10.3. The number of wavefronts showed a similar time course to the number of phase singularities.

CONCLUSIONS

In human ventricular fibrillation, we found an increase in complexity of electric activation patterns during global myocardial ischemia, and this was not reversed during reflow despite an increase in activation rate.

Experiment-model interaction for analysis of epicardial activation during human ventricular fibrillation with global myocardial ischaemia

We describe a combined experiment-modelling framework to investigate the effects of ischaemia on the organisation of ventricular fibrillation in the human heart. In a series of experimental studies epicardial activity was recorded from 10 patients undergoing routine cardiac surgery. Ventricular fibrillation was induced by burst pacing, and recording continued during 2.5 min of global cardiac ischaemia followed by 30 s of coronary reflow. Modelling used a 2D description of human ventricular tissue. Global cardiac ischaemia was simulated by (i) decreased intracellular ATP concentration and subsequent activation of an ATP sensitive K⁺ current, (ii) elevated extracellular K⁺ concentration, and (iii) acidosis resulting in reduced magnitude of the L-type Ca²⁺ current I(Ca,L). Simulated ischaemia acted to shorten action potential duration, reduce conduction velocity, increase effective refractory period, and flatten restitution. In the model, these effects resulted in slower re-entrant activity that was qualitatively consistent with our observations in the human heart. However, the flattening of restitution also resulted in the collapse of many re-entrant waves to several stable re-entrant waves, which was different to the overall trend we observed in the experimental data. These findings highlight a potential role for other factors, such as structural or functional heterogeneity in sustaining wavebreak during human ventricular fibrillation with global myocardial ischaemia.

Regional increase of extracellular potassium leads to electrical instability and reentry occurrence through the spatial heterogeneity of APD restitution

The heterogeneities of electrophysiological properties of cardiac tissue are the main factors that control both arrhythmia induction and maintenance. Although the local increase of extracellular potassium ([K(+)](o)) due to coronary occlusion is a well-established metabolic response to acute ischemia, the role of local [K(+)](o) heterogeneity in phase 1a arrhythmias has yet to be determined. In this work, we created local [K(+)](o) heterogeneity and investigated its role in fast pacing response and arrhythmia induction. The left marginal vein of a Langendorff-perfused rabbit heart was cannulated and perfused separately with solutions containing 4, 6, 8, 10, and 12 mM of K(+). The fluorescence dye was utilized to map the voltage distribution. We tested stimulation rates, starting from 400 ms down to 120 ms, with steps of 5-50 ms. We found that local [K(+)](o) heterogeneity causes action potential (AP) alternans, 2:1 conduction block, and wave breaks. The effect of [K(+)](o) heterogeneity on electrical stability and vulnerability to arrhythmia induction was largest during regional perfusion with 10 mM of K(+). We detected three concurrent dynamics: normally propagating activation when excitation waves spread over tissue perfused with normal K(+), alternating 2:2 rhythm near the border of [K(+)](o) heterogeneity, and 2:1 aperiodicity when propagation was within the high [K(+)](o) area. [K(+)](o) elevation changed the AP duration (APD) restitution and shifted the restitution curve toward longer diastolic intervals and shorter APD. We conclude that spatial heterogeneity of the APD restitution, created with regional elevation of [K(+)](o), can lead to AP instability, 2:1 block, and reentry induction.

Frequency clusters and defect structures in nonlinear dust-density waves under microgravity conditions

Density waves in a dusty plasma emerge spontaneously at low gas pressures and high dust densities. These acousticlike wave modes were studied in a radio-frequency discharge under microgravity conditions. The complex three-dimensional wave pattern shows a spatially varying wavelength that leads to bifurcations, i.e., topological defects, where wave fronts split or merge. The calculation of instantaneous wave attributes from the spatiotemporal evolution of the dust density allows a precise analysis of those structures. Investigations of the spatial frequency distribution inside the wave field revealed that the wave frequency decreases from the bulk to the edge of the cloud in terms of frequency jumps. Between those jumps, regions of almost constant frequency appear. The formation of frequency clusters is strongly correlated with defects that occur exclusively at the cluster boundaries. It is shown that the nonlinearity of the waves has a significant influence on the topology of the wave pattern.

The potential of dual camera systems for multimodal imaging of cardiac electrophysiology and metabolism

Fluorescence imaging has become a common modality in cardiac electrodynamics. A single fluorescent parameter is typically measured. Given the growing emphasis on simultaneous imaging of more than one cardiac variable, we present an analysis of the potential of dual camera imaging, using as an example our straightforward dual camera system that allows simultaneous measurement of two dynamic quantities from the same region of the heart. The advantages of our system over others include an optional software camera calibration routine that eliminates the need for precise camera alignment. The system allows for rapid setup, dichroic image separation, dual-rate imaging, and high spatial resolution, and it is generally applicable to any two-camera measurement. This type of imaging system offers the potential for recording simultaneously not only transmembrane potential and intracellular calcium, two frequently measured quantities, but also other signals more directly related to myocardial metabolism, such as [K(+)](e), NADH, and reactive oxygen species, leading to the possibility of correlative multimodal cardiac imaging. We provide a compilation of dye and camera information critical to the design of dual camera systems and experiments.

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.

Quantitative panoramic imaging of epicardial electrical activity

Fluorescent imaging with voltage- and/or calcium-sensitive dyes has revolutionized cardiac physiology research. Here we present improved panoramic imaging for optically mapping electrical activity from the entire epicardium of the Langendorff-perfused rabbit heart. Combined with reconstruction of the 3D heart surface, the functional data can be conveniently visualized on the realistic heart geometry. Methods to quantify the panoramic data set are introduced by first describing a simple approach to mesh the heart in regular grid form. The regular grid mesh provides substrate for easy translation of previously available non-linear dynamics methods for 2D array data. It also simplifies the unwrapping of curved three-dimensional surface to 2D surface for global epicardial visualization of the functional data. The translated quantification methods include activation maps (isochrones), phase maps, phase singularity, and electric stimulus-induced virtual electrode polarization (VEP) maps. We also adapt a method to calculate the conduction velocities on the global epicardial surface by taking the curvature of the heart surface into account.

Twisted vortex filaments in the three-dimensional complex Ginzburg-Landau equation

The structure and dynamics of vortex filaments that form the cores of scroll waves in three-dimensional oscillatory media described by the complex Ginzburg-Landau equation are investigated. The study focuses on the role that twist plays in determining the bifurcation structure in various regions of the (alpha,beta) parameter space of this equation. As the degree of twist increases, initially straight filaments first undergo a Hopf bifurcation to helical filaments; further increase in the twist leads to a secondary Hopf bifurcation that results in supercoiled helices. In addition, localized states composed of superhelical segments interspersed with helical segments are found. If the twist is zero, zigzag filaments are found in certain regions of the parameter space. In very large systems disordered states comprising zigzag and helical segments with positive and negative senses exist. The behavior of vortex filaments in different regions of the parameter space is explored in some detail. In particular, an instability for nonzero twist near the alpha=beta line suggests the existence of a nonsaturating state that reduces the stability domain of straight filaments. The results are obtained through extensive simulations of the complex Ginzburg-Landau equation on large domains for long times, in conjunction with simulations on equivalent two-dimensional reductions of this equation and analytical considerations based on topological concepts.

Mother rotor anchoring in branching tissue with heterogeneous membrane properties

Abstract Current understanding of atrial fibrillation is based on the co-existence of multiple re-entrant waves propagating randomly throughout the tissue. However, recent experimental results indicate that in many cases one or a small number of periodic, high-frequency re-entrant sources (mother rotors) can drive the arrhythmia. Owing to the high activation rate, mother rotors seem to be located in regions of shortened action potential duration. In this study a computer model of cardiac propagation was applied to investigate mechanisms leading to the formation and maintenance of such mother rotors. For this purpose, a region of short action potential duration was generated by varying the acetylcholine concentration across the tissue. A mother rotor initiated in the center of this region drifts away, and the activation terminates. If an additional heterogeneity such as a bundle is included into the model, a further drift mechanism directed to the bundle is observed and the rotor can be stabilized. Therefore, bundle insertions may play an important role in the maintenance of mother rotors. The influence of the driving rotor on the activation pattern was studied in a three-dimensional model of rectangular shape and a monolayer model of anatomically correct atrial geometry.

Fibrillatory conduction in branching atrial tissue--Insight from volumetric and monolayer computer models

Increased local load in branching atrial tissue (muscle fibers and bundle insertions) influences wave propagation during atrial fibrillation (AF). This computer model study reveals two principal phenomena: if the branching is distant from the driving rotor (>19 mm), the load causes local slowing of conduction or wavebreaks. If the driving rotor is close to the branching, the increased load causes first a slow drift of the rotor towards the branching. Finally, the rotor anchors, and a stable, repeatable pattern of activation can be observed. Variation of the bundle geometry from a cylindrical, volumetric structure to a flat strip of a comparable load in a monolayer model changed the local activation sequence in the proximity of the bundle. However, the global behavior and the basic effects are similar in all models. Wavebreaks in branching tissue contribute to the chaotic nature of AF (fibrillatory conduction). The stabilization (anchoring) of driving rotors by branching tissue might contribute to maintain sustained AF.

Spiral wave dynamics in excitable media with spherical geometries

We describe the spatial and temporal organization of spiral and scroll waves in spherical shells of different sizes and solid spheres. We present simulation results for the evolution of the dynamics and clustering of spiral waves as a function of the excitability of the medium. The excitability, topology, and size of the domain places restrictions on how single and multiarmed spiral waves are organized in space. The results in spherical geometries are compared with those in planar two-dimensional media. These studies are relevant to the dynamics of spiral waves in a variety of media including the heart, and chemical reactions on spherical surfaces.

Evidence for multiple mechanisms in human ventricular fibrillation

BACKGROUND

The mechanisms that sustain ventricular fibrillation (VF) in the human heart remain unclear. Experimental models have demonstrated either a periodic source (mother rotor) or multiple wavelets as the mechanism underlying VF. The aim of this study was to map electrical activity from the entire ventricular epicardium of human hearts to establish the relative roles of these mechanisms in sustaining early human VF.

METHODS AND RESULTS

In 10 patients undergoing cardiac surgery, VF was induced by burst pacing, and 20 to 40 seconds of epicardial activity was sampled (1 kHz) with a sock containing 256 unipolar contact electrodes connected to a UnEmap system. Signals were interpolated from the electrode sites to a fine regular grid (100x100 points), and dominant frequencies (DFs) were calculated with a fast Fourier transform with a moving 4096-ms window (10-ms increments). Epicardial phase was calculated at each grid point with the Hilbert transform, and phase singularities and activation wavefronts were identified at 10-ms intervals. Early human VF was sustained by large coherent wavefronts punctuated by periods of disorganized wavelet behavior. The initial fitted DF intercept was 5.11 +/- 0.25 (mean +/- SE) Hz (P < 0.0001), and DF increased at a rate of 0.018 +/- 0.005 Hz/s (P < 0.01) during VF, whereas combinations of homogeneous, heterogeneous, static, and mobile DF domains were observed for each of the patients. Epicardial reentry was present in all fibrillating hearts, typically with low numbers of phase singularities. In some cases, persistent phase singularities interacted with multiple complex wavelets; in other cases, VF was driven at times by a single reentrant wave that swept the entire epicardium for several cycles.

CONCLUSIONS

Our data support both the mother rotor and multiple wavelet mechanisms of VF, which do not appear to be mutually exclusive in the human heart.

Mechanisms of unpinning and termination of ventricular tachycardia

High-energy defibrillation shock is the only therapy for ventricular tachyarrhythmias. However, because of adverse side effects, lowering defibrillation energy is desirable. We investigated mechanisms of unpinning, destabilization, and termination of ventricular tachycardia (VT) by low-energy shocks in isolated rabbit right ventricular preparations (n = 22). Stable VT was initiated with burst pacing and was optically mapped. Monophasic "unpinning" shocks (10 ms) of different strengths were applied at various phases throughout the reentry cycle. In 8 of 22 preparations, antitachycardia pacing (ATP: 8-20 pulses, 50-105% of period, 0.8-10 mA) was also applied. Termination of reentry by ATP was achieved in only 5 of 8 preparations. Termination by unpinning occurred in all 22 preparations. Rayleigh's test showed a statistically significant unpinning phase window, during which reentry could be unpinned and subsequently terminated with E80 (magnitude at which 80% of reentries were unpinned) = 1.2 V/cm. All reentries were unpinned with field strengths < or = 2.4 V/cm. Unpinning was achieved by inducing virtual electrode polarization and secondary sources of excitation at the core of reentry. Optical mapping revealed the mechanisms of phase-dependent unpinning of reentry. These results suggest that a 20-fold reduction in energy could be achieved compared with conventional high-energy defibrillation and that the unpinning method may be more effective than ATP for terminating stable, pinned reentry in this experimental model.

Mechanisms of enhanced shock-induced arrhythmogenesis in the rabbit heart with healed myocardial infarction

Shock-induced vulnerability and defibrillation have been mostly studied in structurally normal hearts. However, defibrillation therapy is normally applied to patients with diseased hearts, frequently those with prior myocardial infarction (MI). Shock-induced vulnerability and defibrillation have not been well studied under this condition. We sought to examine the mechanisms of shock-induced arrhythmogenesis and arrhythmia maintenance in a rabbit model of healed MI (4 wk or more postinfarction). Ligation of the lateral division or posterolateral division of the left coronary artery at a level of 40-70% from the apex was performed 53 +/- 21 days before acute experiments. Shock-induced vulnerability was assessed in infarcted (n = 8) and structurally normal (n = 8) hearts by delivering internal monophasic shocks at different shock strengths and delivery phases. Electrical activities from the anterior epicardium during shock application and during shock-induced arrhythmias were optically recorded and quantitatively analyzed. Ligation resulted in a transmural left ventricular free wall infarction mainly located at the apical region with a consistent endocardial border zone (BZ) as confirmed by histological studies. There were significant increases in the incidence, severity, and duration of shock-induced arrhythmias in the infarcted hearts versus controls due to 1) postshock break-excitation wavefronts that frequently originated near the infarction BZ and 2) the existence of an infarction BZ that created an anatomic reentry pathway and facilitated arrhythmia maintenance. In conclusion, the infarction BZ contributes to both increased shock-induced arrhythmogenesis and arrhythmia maintenance in the rabbit model of healed MI.

Dynamics of virtual electrode-induced scroll-wave reentry in a 3D bidomain model

Functional reentry in the heart can be caused by a wave front of excitation rotating around its edge. Previous simulations on the basis of monodomain cable equations predicted the existence of self-sustained, vortex-like wave fronts (scroll waves) rotating around a filament in three dimensions. In our simulations, we used the more accurate bidomain model with modified Beeler-Reuter ionic kinetics to study the dynamics of scroll-wave filaments in a 16 x 8 x 1.5-mm slab of ventricular tissue with straight fibers. Wave fronts were identified as the areas with inward current. Their edges represented the filaments. Both transmural and intramural reentries with I- and U-shaped filaments, respectively, were obtained by the S1-S2 point stimulation protocol through the virtual electrode-induced phase singularity mechanism. The filaments meandered along elongated trajectories and tended to attach to the tissue boundaries exposed to air (no current flow) rather than to the bath (zero extracellular potential). They completely detached from electroporated (zero transmembrane potential) boundaries. In our simulations, the presence of the bath led to generation of only U-shaped filaments, which survived for the 1.5-mm-thick slab but not for the slabs of 0.5- or 3-mm thicknesses. Thus boundary conditions may be another determinant of the type and dynamics of reentry.

Examination of optical depth effects on fluorescence imaging of cardiac propagation

Optical mapping with voltage-sensitive dyes provides a high-resolution technique to observe cardiac electrodynamic behavior. Although most studies assume that the fluorescent signal is emitted from the surface layer of cells, the effects of signal attenuation with depth on signal interpretation are still unclear. This simulation study examines the effects of a depth-weighted signal on epicardial activation patterns and filament localization. We simulated filament behavior using a detailed cardiac model, and compared the signal obtained from the top (epicardial) layer of the spatial domain with the calculated weighted signal. General observations included a prolongation of the action upstroke duration, early upstroke initiation, and reduction in signal amplitude in the weighted signal. A shallow filament was found to produce a dual-humped action potential morphology consistent with previously reported observations. Simulated scroll wave breakup exhibited effects such as the false appearance of graded potentials, apparent supramaximal conduction velocities, and a spatially blurred signal with the local amplitude dependent upon the immediate subepicardial activity; the combination of these effects produced a corresponding change in the accuracy of filament localization. Our results indicate that the depth-dependent optical signal has significant consequences on the interpretation of epicardial activation dynamics.

Combined phase singularity and wavefront analysis for optical maps of ventricular fibrillation

Much of the research into the mechanisms of ventricular fibrillation (VF) employs high-resolution mapping of electrical activation and recovery patterns. We previously developed a method for analyzing electrically mapped VF patterns that was based on identifying individual VF wavefronts. We now introduce a related method designed to take into account the information on repolarization that is present in optically mapped VF data. The new method first converts raw fluorescence data to an angular variable that tracks the phase of the mapped tissue through the depolarization-repolarization cycle. We define wavefronts in this context as isolines of phase that terminate either at boundaries or at singular points within the phase field. These singularities are the pivots of functional reentry and are important determinants of VF patterns. We parameterize VF by constructing data structures that describe wavefronts and singularities and also maintain wavefront-wavefront, wavefront-singularity, and singularity-singularity relationships. We describe one important application of this parameterization, which is to identify, localize, and characterize the importance of occurrences of propagation block during VF.

Effects of lidocaine on shock-induced vulnerability

INTRODUCTION

Lidocaine is known to increase the defibrillation threshold (DFT) of monophasic shocks (MS) and have no effect on DFT of biphasic shocks (BS). The aim of this study was to enhance our understanding of the mechanisms of vulnerability and defibrillation through the investigation of this difference.

METHODS AND RESULTS

We studied the effect of 15 microM lidocaine on shock-induced vulnerability using fluorescent imaging of Langendorff-perfused rabbit hearts. Vulnerability was assessed as vulnerable window with shock strengths of 15 to 150 V and vulnerable period (VP) with shock delivery phase of 0% to 100% of action potential duration (% APD). With MS, lidocaine caused a significant increase in both the upper limit of vulnerability (ULV, 71 +/- 17 V vs 120 +/- 1.5 V, P < 0.01) and upper limit of VP (91 +/- 8.0% APD vs 110 +/- 4.2% APD, P < 0.01). With BS, lidocaine had no effect on ULV (40 +/- 3.4 V vs 45 +/- 4.5 V) and did not increase the upper limit of VP (78 +/- 8.9% APD vs 96 +/- 12% APD, P < 0.01). Lidocaine caused reduction of the conduction velocity during pacing (0.58 +/- 0.08 m/s vs 0.44 +/- 0.05 m/s, P < 0.01), shock-induced break excitation (0.82 +/- 0.17 m/s vs 0.30 +/- 0.07 m/s, P < 0.01), and postshock reentry (0.34 +/- 0.07 m/s vs 0.19 +/- 0.08 m/s, P < 0.01). Lidocaine had no effect on shock-induced virtual electrode polarization.

CONCLUSION

Lidocaine increased MS ULV due to slowing of shock-induced break-excitation wavefronts, which resulted in enhanced probability of survival of virtual electrode induced phase singularity. Lidocaine had no effect on BS ULV because no break excitation was induced by BS. Reduction of conduction velocity by lidocaine resulted in increased dispersion of repolarization and led to upper limit of VP increase for both MS and BS.