Reversal of repolarization gradient does not reverse the chirality of shock-induced reentry in the rabbit heart

Bidomain modeling of electrical and mechanical properties of cardiac tissue

Throughout the history of cardiac research, there has been a clear need to establish mathematical models to complement experimental studies. In an effort to create a more complete picture of cardiac phenomena, the bidomain model was established in the late 1970s to better understand pacing and defibrillation in the heart. This mathematical model has seen ongoing use in cardiac research, offering mechanistic insight that could not be obtained from experimental pursuits. Introduced from a historical perspective, the origins of the bidomain model are reviewed to provide a foundation for researchers new to the field and those conducting interdisciplinary research. The interplay of theory and experiment with the bidomain model is explored, and the contributions of this model to cardiac biophysics are critically evaluated. Also discussed is the mechanical bidomain model, which is employed to describe mechanotransduction. Current challenges and outstanding questions in the use of the bidomain model are addressed to give a forward-facing perspective of the model in future studies.

Imaging of Ventricular Fibrillation and Defibrillation: The Virtual Electrode Hypothesis

Ventricular fibrillation is the major underlying cause of sudden cardiac death. Understanding the complex activation patterns that give rise to ventricular fibrillation requires high resolution mapping of localized activation. The use of multi-electrode mapping unraveled re-entrant activation patterns that underlie ventricular fibrillation. However, optical mapping contributed critically to understanding the mechanism of defibrillation, where multi-electrode recordings could not measure activation patterns during and immediately after a shock. In addition, optical mapping visualizes the virtual electrodes that are generated during stimulation and defibrillation pulses, which contributed to the formulation of the virtual electrode hypothesis. The generation of virtual electrode induced phase singularities during defibrillation is arrhythmogenic and may lead to the induction of fibrillation subsequent to defibrillation. Defibrillating with low energy may circumvent this problem. Therefore, the current challenge is to use the knowledge provided by optical mapping to develop a low energy approach of defibrillation, which may lead to more successful defibrillation.

How hyperpolarization and the recovery of excitability affect propagation through a virtual anode in the heart

Researchers have suggested that the fate of a shock-induced wave front at the edge of a "virtual anode" (a region hyperpolarized by the shock) is a key factor determining success or failure during defibrillation of the heart. In this paper, we use a simple one-dimensional computer model to examine propagation speed through a hyperpolarized region. Our goal is to test the hypothesis that rapid propagation through a virtual anode can cause failure of propagation at the edge of the virtual anode. The calculations support this hypothesis and suggest that the time constant of the sodium inactivation gate is an important parameter. These results may be significant in understanding the mechanism of the upper limit of vulnerability.

A mechanism for the upper limit of vulnerability

BACKGROUND

The strongest shock that induces reentry in the heart is the upper limit of vulnerability (ULV). In order to understand defibrillation, one must know what causes the ULV.

OBJECTIVE

The goal of this study was to examine the mechanism of the upper limit of vulnerability.

METHODS

Numerical simulations of cardiac tissue were performed using the bidomain model. An S2 shock was applied during the refractory period of the S1 action potential, and results using a smooth curving fiber geometry were compared with results using a smooth plus random fiber geometry.

RESULTS

When using a smooth fiber geometry only, no ULV was observed. However, when a random fiber geometry was included, the ULV was present. The difference arises from the fate of the shock-induced break wave front when it reaches the edge of the tissue hyperpolarized by the shock (the virtual anode).

CONCLUSION

Our numerical simulations suggest that local heterogeneities throughout the tissue may be crucial for determining the fate of the shock-induced wave front at the edge of the virtual anode, and therefore play an important role in the mechanism underlying the ULV.

The role of mechanoelectric feedback in vulnerability to electric shock

Experimental and clinical studies have shown that ventricular dilatation is associated with increased arrhythmogenesis and elevated defibrillation threshold; however, the underlying mechanisms remain poorly understood. The goal of the present study was to test the hypothesis that (1) stretch-activated channel (SAC) recruitment and (2) geometrical deformations in organ shape and fiber architecture lead to increased arrhythmogenesis by electric shocks following acute ventricular dilatation. To elucidate the contribution of these two factors, the study employed, for the first time, a combined electro-mechanical simulation approach. Acute dilatation was simulated in a model of rabbit ventricular mechanics by raising the LV end-diastolic pressure from 0.6 (control) to 4.2 kPa (dilated). The output of the mechanics model was used in the electrophysiological model. Vulnerability to shocks was examined in the control, the dilated ventricles, and in the dilated ventricles that also incorporated currents through SAC as a function of local strain, by constructing vulnerability grids. Results showed that dilatation-induced deformation alone decreased upper limit of vulnerability (ULV) slightly and did not result in increased vulnerability. With SAC recruitment in the dilated ventricles, the number of shock-induced arrhythmia episodes increased by 37% (from 41 to 56) and the lower limit of vulnerability (LLV) decreased from 9 to 7 V/cm, while ULV did not change. The heterogeneous activation of SAC caused by the heterogeneous fiber strain in the ventricular walls was the main reason for increased vulnerability to electric shocks since it caused dispersion of electrophysiological properties in the tissue, resulting in postshock unidirectional block and establishment of reentry.

The pinwheel experiment revisited: effects of cellular electrophysiological properties on vulnerability to cardiac reentry

In normal heart, ventricular fibrillation can be induced by a single properly timed strong electrical or mechanical stimulus. A mechanism first proposed by Winfree and coined the “pinwheel experiment” emphasizes the timing and strength of the stimulus in inducing figure-of-eight reentry. However, the effects of cellular electrophysiological properties on vulnerability to reentry in the pinwheel scenario have not been investigated. In this study, we extend Winfree's pinwheel experiment to show how the vulnerability to reentry is affected by the graded action potential responses induced by a strong premature stimulus, action potential duration (APD), and APD restitution in simulated monodomain homogeneous two-dimensional tissue. We find that a larger graded response, longer APD, or steeper APD restitution slope reduces the vulnerable window of reentry. Strong graded responses and long APD promote tip-tip interactions at long coupling intervals, causing the two initiated spiral wave tips to annihilate. Steep APD restitution promotes wave front-wave back interaction, causing conduction block in the central common pathway of figure-of-eight reentry. We derive an analytical treatment that shows good agreement with numerical simulation results.

Activation and repolarization patterns in the ventricular epicardium under sinus rhythm in frog and rabbit hearts

Our study compared the contributions of activation sequence and local repolarization durations distribution in the organization of epicardial repolarization in animals with fast (rabbit) and slow (frog) myocardial activation under sinus rhythm. Activation times, repolarization times and activation-recovery intervals (ARI) were obtained from ventricular epicardial unipolar electrograms recorded in 13 Chinchilla rabbits (Oryctolagus cuniculus) and 10 frogs (Rana temporaria). In frogs, depolarization travels from the atrioventricular ring radially. ARIs increased progressively from the apex to the middle portion and finally to the base (502+/-75, 557+/-73, 606+/-79 ms, respectively; P<0.01). In rabbits, depolarization spread from two epicardial breakthroughs with the duration of epicardial activation being lower than that in frogs (17+/-3 vs. 44+/-18 ms; P<0.001). ARI durations were 120+/-37, 143+/-45, and 163+/-40 ms in the left ventricular apex, left, and right ventricular bases, respectively (P<0.05). In both species, repolarization sequence was directed from apex to base according to the ARI distribution with dispersion of repolarization being higher than that of activation (P<0.001). Thus, excitation spread sequence and velocity per se do not play a crucial role in the formation of ventricular epicardial repolarization pattern, but the chief factor governing repolarization sequences is the distribution of local repolarization durations.

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.

Refractory Gradient is Responsible for the Increase in Ventricular Vulnerability Under Sodium Channel Blockade

BACKGROUND

Previous studies have shown that sodium channel (I(Na)) blockade increases ventricular vulnerability; however, there were differences in the degree of the increase. Because the vulnerable window (VW) is altered by the type of preshock refractory gradient (RG), the hypothesis was that the differences in the arrhythmogenesis of I(Na) blockade result from the different types of preshock RG employed.

METHODS AND RESULTS

Simulations of regio(Na)l electric shock following constant pacing stimuli in 2-dimensional bidomain myocardial sheets under I(Na) blockade were conducted using 3 types of preshock RG: longitudinally tilted (LRG), transversely tilted (TRG), and non-tilted RG (NRG). The increase in the degree of I(Na) blockade almost linearly decreased the conduction velocity. The action potential duration in the LRG and TRG cases was non-linearly shortened with the increase in INa blockade because of electrotonic influences, whereas in the case of NRG it was slightly prolonged. In both LRG and TRG cases, the VW for reentry induction by electric shock was considerably widened by the INa blockade; however, this was not the case for NRG in which the VW was rather narrowed by the INa blockade.

CONCLUSION

The type of preshock RG alters the degree of the increase in ventricular vulnerability under INa blockade.

Art Winfree and the bidomain model of cardiac tissue

This paper reviews Art Winfree's contributions to the bidomain model of cardiac tissue. Specifically, he first predicted quatrefoil reentry, he showed that an S1 refractory gradient is not required for an S2 stimulus to induce reentry, and his work on spiral wave meandering led to studies on how the path of the tip of a spiral wave is influenced by tissue anisotropy.

Postshock arrhythmogenesis in a slice of the canine heart

INTRODUCTION

Recent evidence has demonstrated that defibrillation shocks terminate or reset reentrant activity in the myocardium through the generation of virtual electrode polarization (VEP). Previous research has revealed that the shock establishes phase singularities (PSs) in the tissue via the VEP mechanism. The aim of this study was to examine, as a function of shock strength and electrode configuration, the relationship between end-shock PSs and the reentrant circuits established after failed defibrillation attempts.

METHODS AND RESULTS

The study uses a complex three-dimensional finite-element bidomain model of a slice of the canine heart characterized by realistic geometry and fiber architecture and undergoing a single scroll wave. Defibrillation shocks of increasing strength are delivered through three different electrode configurations. The results demonstrated that >98% of all PSs have a lifetime of half a reentrant cycle or less. Stronger shocks result in a faster rate of annihilation of postshock PSs. For failed shocks, the surviving PSs underlie the activity of one or more scroll waves, which remain stationary in the slice. For all electrode configurations tested, the increase in shock strength leads to a rapid initial increase in the number of postshock reentries followed by a slower decrease; similar behavior is observed with regard to end-shock PSs.

CONCLUSION

These results present new evidence regarding the mechanisms underlying failure of defibrillation shocks.

Shock-induced arrhythmogenesis is enhanced by 2,3-butanedione monoxime compared with cytochalasin D

Investigation of the mechanisms of arrhythmia genesis and maintenance has benefited from the use of optical mapping techniques that employ excitation-contraction uncouplers. We investigated the effects of the excitation-contraction uncouplers 2,3-butanedione monoxime (BDM) and cytochalasin D (Cyto D) on the induction and maintenance of arrhythmia by electric shocks. Electrical activity was optically mapped from anterior epicardium of rabbit hearts (n = 9) during shocks (-100 V, 8 ms) applied from a ventricular lead at various phases of action potential duration (APD). Restitution curves were obtained using S1-S2 protocol and measurement of APD values at 70% of repolarization. Compared with Cyto D, BDM significantly shortened APD at 90% of repolarization, although no significant difference in dispersion of repolarization was observed. Wavelength was also shortened with BDM. In general, shock-induced arrhythmias with BDM and Cyto D were ventricular tachycardic in nature. With respect to shock-induced sustained arrhythmias, the vulnerable window was wider and the incidence was higher with BDM than with Cyto D. There was also a difference in the morphology of ventricular tachycardia (VT) between the two agents. The arrhythmias with BDM usually resembled monomorphic VT, especially those that lasted >30 s. In contrast, arrhythmias with Cyto D more resembled polymorphic VT. However, the average number of phase singularities increased under Cyto D vs. BDM, whereas no significant difference in the dominant frequency of shock-induced sustained arrhythmia was observed. BDM reduced the slope of the restitution curve compared with Cyto D, but duration of arrhythmia under BDM was significantly increased compared with Cyto D. In conclusion, BDM increased arrhythmia genesis and maintenance relative to Cyto D.

Mechanistic inquiry into decrease in probability of defibrillation success with increase in complexity of preshock reentrant activity

Energy requirements for successful antiarrhythmia shocks are arrhythmia specific. However, it remains unclear why the probability of shock success decreases with increasing arrhythmia complexity. The goal of this research was to determine whether a diminished probability of shock success results from an increased number of functional reentrant circuits in the myocardium, and if so, to identify the responsible mechanisms. To achieve this goal, we assessed shock efficacy in a bidomain defibrillation model of a 4-mm-thick slice of canine ventricles. Shocks were applied between a right ventricular cathode and a distant anode to terminate either a single scroll wave (SSW) or multiple scroll waves (MSWs). From the 160 simulations conducted, dose-response curves were constructed for shocks given to SSWs and MSWs. The shock strength that yielded a 50% probability of success (ED(50)) for SSWs was found to be 13% less than that for MSWs, which indicates that a larger number of functional reentries results in an increased defibrillation threshold. The results also demonstrate that an isoelectric window exists after both failed and successful shocks; however, shocks of strength near the ED(50) value that were given to SSWs resulted in 16.3% longer isoelectric window durations than the same shocks delivered to MSWs. Mechanistic inquiry into these findings reveals that the two main factors underlying the observed relationships are 1) smaller virtual electrode polarizations in the tissue depth, and 2) differences in preshock tissue state. As a result of these factors, intramural excitable pathways leading to delayed breakthrough on the surface were formed earlier after shocks given to MSWs compared with SSWs and thus resulted in a lower defibrillation threshold for shocks given to SSWs.

Development of a computer algorithm for the detection of phase singularities and initial application to analyze simulations of atrial fibrillation

Atrial fibrillation (AF) is a common cardiac arrhythmia, but its mechanisms are incompletely understood. The identification of phase singularities (PSs) has been used to define spiral waves involved in maintaining the arrhythmia, as well as daughter wavelets. In the past, PSs have often been identified manually. Automated PS detection algorithms have been described previously, but when we attempted to apply a previously developed algorithm we experienced problems with false positives that made the results difficult to use directly. We therefore developed a tool for PS identification that uses multiple strategies incorporating both image analysis and mathematical convolution for automated detection with optimized sensitivity and specificity, followed by manual verification. The tool was then applied to analyze PS behavior in simulations of AF maintained in the presence of spatially distributed acetylcholine effects in cell grids of varying size. These analyses indicated that in almost all cases, a single PS lasted throughout the simulation, corresponding to the central-core tip of a single spiral wave that maintained AF. The sustained PS always localized to an area of low acetylcholine concentration. When the grid became very small and no area of low acetylcholine concentration was surrounded by zones of higher concentration, AF could not be sustained. The behavior of PSs and the mechanisms of AF were qualitatively constant over an 11.1-fold range of atrial grid size, suggesting that the classical emphasis on tissue size as a primary determinant of fibrillatory behavior may be overstated. (c) 2002 American Institute of Physics.

Mechanisms of shock-induced arrhythmogenesis during acute global ischemia

Little is known about the mechanisms of vulnerability and defibrillation under ischemic conditions. We investigated these mechanisms in 18 Langendorff-perfused rabbit hearts during 75% reduced-flow ischemia. Electrical activity was optically mapped from the anterior epicardium during right ventricular shocks applied at various phases of the cardiac cycle while the excitation-contraction decoupler 2,3-butanedione monoxime (BDM; 15 mM) was used to suppress motion artifacts caused by contraction of the heart. During ischemia, vulnerable window width increased [from 30-90% of the action potential duration (APD) in the control to -10 to 100% of the APD in ischemia]. Moreover, arrhythmia severity increased along with the reduction of APD (176 +/- 9 ms in control and 129 +/- 26 ms in ischemia, P < 0.01) and increased dispersion of repolarization (45 +/- 17 ms in control and 73 +/- 28 ms in ischemia, P < 0.01). Shock-induced virtual electrode polarization was preserved. Depolarizing (contrary to hyperpolarizing) response time constants increased. Virtual electrode-induced wavefronts of excitation had much more tortuous pathways leading to wavefront fractionation. Defibrillation failure at all shock strengths was observed in four hearts. Optical mapping revealed that the shock extinguished the arrhythmia; however, the arrhythmia self-originated after an isoelectric window of 339 +/- 189 ms. In conclusion, in most cases, virtual electrode-induced phase singularity (VEIPS) was responsible for shock-induced arrhythmogenesis during acute global ischemia. Enhancement of arrhythmogenesis was associated with an increased dispersion of repolarization and altered deexcitation. In four hearts, arrhythmogenesis could not be explained by VEIPS.

The pinwheel experiment re-revisited

Virtual electrode induced phase singularity hypothesis explains the origin of cardiac arrhythmias caused by artificial electrical induction of rotors, i.e. vortex-like self-sustained sources of activity. This mechanism is thought to underlie both stimulus-induced arrhythmias and shock defibrillation therapy. In this paper, we extend this hypothesis to three dimensions using the bidomain model of cardiac tissue. We predict that virtual electrode polarization can produce three topologically distinct types of rotors anchored to: (1) transmural I-shaped scroll wave filaments; (2) near-surface U-shaped scroll wave filaments; and (3) intramural O-shaped scroll wave filaments.