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

Low-energy, single-pulse surface stimulation defibrillates large mammalian ventricles

BACKGROUND

Strong electric shocks are the gold standard for ventricular defibrillation but are associated with pain and tissue damage. We hypothesized that targeting the excitable gap (EG) of reentry with low-energy surface stimulation is a less damaging and painless alternative for ventricular defibrillation.

OBJECTIVE

The purpose of this study was to determine the conditions under which low-energy surface stimulation defibrillates large mammalian ventricles.

METHODS

Low-energy surface stimulation was delivered with five electrodes that were 7 cm long and placed 1-2 cm apart on the endocardial and epicardial surfaces of perfused pig left ventricle (LV). Rapid pacing (>4 Hz) was used to induce reentry from a single electrode. A 2 ms defibrillation pulse ≤0.5 A was delivered from all electrodes with a varied time delay from the end of the induction protocol (0.1-5 seconds). Optical mapping was performed and arrhythmia dynamics analyzed. For mechanistic insight, simulations of the VF induction and defibrillation protocols were performed in silico with an LV model emulating the experimental conditions and electrodes placed 0.25-2 cm apart.

RESULTS

In living LV, reentry was induced with varying complexity and dominant frequencies ranging between 3.5 to 6.2 Hz over 8 seconds postinitiation. Low-energy defibrillation was achieved with energy <60 mJ and electrode separations up to 2 cm for less complex arrhythmia. In simulations, defibrillation consistently occurred when stimulation captured >75% of the EG, which blocked reentry <2.9 mm in front of the leading reentrant wavefront.

CONCLUSION

Defibrillation with low-energy, single-pulse surface stimulation is feasible with energies below the human pain threshold (100 mJ). Optimal defibrillation occurs when arrhythmia complexity is minimal and electrodes capture >75% of the EG.

Therapy-Resistant Ventricular Arrhythmias Developed More Often in Advanced Than in Therapeutic Mild Hypothermic Condition

Therapy-resistant ventricular arrhythmias can occur during accidental advanced hypothermic conditions. On the other hand, hypothermic therapy using mild cooling has been successfully accomplished with infrequent ventricular arrhythmia events.To further clarify the therapeutic-resistant arrhythmogenic substrate which develops in hypothermic conditions, an experimental study was performed using a perfusion wedge preparation model of porcine ventricle, and electrophysiological characteristics, inducibility of ventricular arrhythmias, and effects of therapeutic interventions were assessed at 3 target temperatures (37, 32 and 28°C).As the myocardial temperature decreased, myocardial contractions and the number of spontaneous beats deceased. Depolarization (QRS width, stimulus-QRS interval) and repolarization (QT interval, ERP) parameters progressively increased, and dispersion of the ventricular repolarization increased. At 28°C, VF tended to be inducible more frequently (1/11 at 37°C, 1/11 at 32°C, and 5/11 hearts at 28°C), and some VFs at 28°C required greater defibrillation energy to resume basic rhythm.An advanced but not a mild hypothermic condition had an enhanced arrhythmogenic potential in our model. In the advanced hypothermic condition, VF with relatively prolonged F-F intervals and a greater defibrillation energy were occasionally inducible based on the arrhythmogenic substrate characterized as slowed conduction and prolonged repolarization of the ventricle.

Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia

Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.

PDE constrained optimization of electrical defibrillation in a 3D ventricular slice geometry

A computational study of an optimal control approach for cardiac defibrillation in a 3D geometry is presented. The cardiac bioelectric activity at the tissue and bath volumes is modeled by the bidomain model equations. The model includes intramural fiber rotation, axially symmetric around the fiber direction, and anisotropic conductivity coefficients, which are extracted from a histological image. The dynamics of the ionic currents are based on the regularized Mitchell-Schaeffer model. The controls enter in the form of electrodes, which are placed at the boundary of the bath volume with the goal of dampening undesired arrhythmias. The numerical optimization is based on Newton techniques. We demonstrated the parallel architecture environment for the computation of potentials on multidomains and for the higher order optimization techniques.

Constitutively Active Acetylcholine-Dependent Potassium Current Increases Atrial Defibrillation Threshold by Favoring Post-Shock Re-Initiation

Electrical cardioversion (ECV), a mainstay in atrial fibrillation (AF) treatment, is unsuccessful in up to 10–20% of patients. An important aspect of the remodeling process caused by AF is the constitutive activition of the atrium-specific acetylcholine-dependent potassium current (I_K,ACh → I_K,ACh-c), which is associated with ECV failure. This study investigated the role of I_K,ACh-c in ECV failure and setting the atrial defibrillation threshold (aDFT) in optically mapped neonatal rat cardiomyocyte monolayers. AF was induced by burst pacing followed by application of biphasic shocks of 25–100 V to determine aDFT. Blocking I_K,ACh-c by tertiapin significantly decreased DFT, which correlated with a significant increase in wavelength during reentry. Genetic knockdown experiments, using lentiviral vectors encoding a Kcnj5-specific shRNA to modulate I_K,ACh-c, yielded similar results. Mechanistically, failed ECV was attributed to incomplete phase singularity (PS) removal or reemergence of PSs (i.e. re-initiation) through unidirectional propagation of shock-induced action potentials. Re-initiation occurred at significantly higher voltages than incomplete PS-removal and was inhibited by I_K,ACh-c blockade. Whole-heart mapping confirmed our findings showing a 60% increase in ECV success rate after I_K,ACh-c blockade. This study provides new mechanistic insight into failing ECV of AF and identifies I_K,ACh-c as possible atrium-specific target to increase ECV effectiveness, while decreasing its harmfulness.

Frequency characteristics and associations with the defibrillation threshold of ventricular fibrillation in patients with implantable cardioverter defibrillators

OBJECTIVE

The dominant frequency (DF) in frequency analyses is considered to represent the objective cycle length and complexity of activation under conditions of ventricular fibrillation (VF). However, knowledge regarding the mechanisms determining the DF in human VF is limited. We studied the characteristics of the DF of human VF and relationship between DF and the defibrillation threshold.

METHODS

Seventy-two implantable cardioverter-defibrillator patients and 211 VF were studied. Using defibrillation tests, we performed a frequency analysis with fast Fourier transformation. The correlations between DF and clinical characteristics, including the defibrillation threshold, were assessed.

RESULTS

The mean DF of all induced VFs was 5.2±0.8 Hz. The patients were divided into two groups according to DF: the low-DF (DF <5.2 Hz, n=32) and high-DF (DF ≥5.2 Hz, n=40) groups. The frequency of structural heart disease was significantly higher in the low-DF group. In addition, the QRS duration, QT interval and effective refractory period of the right ventricle (RV-ERP) were significantly longer in the low-DF group. A multivariate analysis showed RV-ERP to be the only independent predictor of DF. Excluding patients receiving group III anti-arrhythmic drugs, which are known to have potent defibrillation threshold effects, the defibrillation threshold was significantly lower in the low-DF group (p=0.026).

CONCLUSION

We found that the DF of human VF is associated with underlying heart disease, the cardiac function, cardiac conduction, ventricular refractoriness and defibrillation threshold. Our findings may be useful for identifying and managing patients with a high defibrillation threshold.

New insights into defibrillation of the heart from realistic simulation studies

Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices, constitutes the most important means of combating sudden cardiac death. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart, particularly under diseased conditions, is a promising approach to achieve an optimal therapy. The aim of this article is to assess the current state-of-the-art in whole-heart defibrillation modelling, focusing on major insights that have been obtained using defibrillation models, primarily those of realistic heart geometry and disease remodelling. The article showcases the contributions that modelling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modelling of the heart (i.e. cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level, and has the potential to contribute to the betterment of the clinical practice of defibrillation.

Prolongation of minimal action potential duration in sustained fibrillation decreases complexity by transient destabilization

AIMS

Sustained ventricular fibrillation (VF) is maintained by multiple stable rotors. Destabilization of sustained VF could be beneficial by affecting VF complexity (defined by the number of rotors). However, underlying mechanisms affecting VF stability are poorly understood. Therefore, the aim of this study was to correlate changes in arrhythmia complexity with changes in specific electrophysiological parameters, allowing a search for novel factors and underlying mechanisms affecting stability of sustained VF.

METHODS AND RESULTS

Neonatal rat ventricular cardiomyocyte monolayers and Langendorff-perfused adult rat hearts were exposed to increasing dosages of the gap junctional uncoupler 2-aminoethoxydiphenyl borate (2-APB) to induce arrhythmias. Ion channel blockers/openers were added to study effects on VF stability. Electrophysiological parameters were assessed by optical mapping and patch-clamp techniques. Arrhythmia complexity in cardiomyocyte cultures increased with increasing dosages of 2-APB (n > 38), leading to sustained VF: 0.0 ± 0.1 phase singularities/cm(2) in controls vs. 0.0 ± 0.1, 1.0 ± 0.9, 3.3 ± 3.2, 11.0 ± 10.1, and 54.3 ± 21.7 in 5, 10, 15, 20, and 25 µmol/L 2-APB, respectively. Arrhythmia complexity inversely correlated with wavelength. Lengthening of wavelength during fibrillation could only be induced by agents (BaCl(2)/BayK8644) increasing the action potential duration (APD) at maximal activation frequencies (minimal APD); 123 ± 32%/117 ± 24% of control. Minimal APD prolongation led to transient VF destabilization, shown by critical wavefront collision leading to rotor termination, followed by significant decreases in VF complexity and activation frequency (52%/37%). These key findings were reproduced ex vivo in rat hearts (n = 6 per group).

CONCLUSION

These results show that stability of sustained fibrillation is regulated by minimal APD. Minimal APD prolongation leads to transient destabilization of fibrillation, ultimately decreasing VF complexity, thereby providing novel insights into anti-fibrillatory mechanisms.

Slowing of Electrical Activity in Ventricular Fibrillation is Not Associated with Increased Defibrillation Energies in the Isolated Rabbit Heart

Prolonged out-of-hospital ventricular fibrillation (VF) arrests are associated with reduced ECG dominant frequency (DF) and diminished defibrillation success. Partial reversal of ischemia increases ECG DF and improves defibrillation outcome. We have investigated the metabolic components of ischemia responsible for the decline in ECG DF and defibrillation success. Isolated Langendorff-perfused rabbit hearts were loaded with the voltage-sensitive dye RH237. Using a photodiode array, epicardial membrane potentials were recorded at 252 sites (15 mm × 15 mm) on the anterior surface of the left and right ventricles. Simultaneously, a global ECG was recorded. VF was induced by burst pacing, and after 60s, perfusion was either reduced to 6 ml/min or the perfusate composition changed to impose hypoxia (95% N(2)/5% CO(2)), pH 6.7 (80% O(2)/20% CO(2)), or hyperkalemia (8 mM). Using fast Fourier transform, power spectra were created from the optical signals and the global ECG. The optical power spectra were summated to give a global power spectrum (pseudoECG). At 600 s the minimum defibrillation voltage (MDV) was determined by step-up protocol. During VF, the ECG and pseudoECG DF were reduced by low-flow ischemia (9.0 ± 1.0 Hz, p < 0.01, n = 5) and raised [K(+)](o) (12.2 ± 1.3 Hz, p < 0.05, n = 7) compared to control (19.2 ± 1.5 Hz, n = 20), but were unaffected by acidic pH(o) (16.7 ± 1.1 Hz, n = 11) and hypoxia (14.0 ± 1.2 Hz, n = 10). In contrast, the MDV was raised by acidic pH (156.1 ± 26.4 V, p < 0.001) and hypoxia (154.1 ± 22.1 V, p < 0.01) compared to control (65.6 ± 2.3 V), but comparable changes were not observed in low-flow ischemia (61.0 ± 0.5 V) or raised [K(+)](o) (56 ± 3 V). In summary, different metabolites are responsible for the reduction in DF and the increase in defibrillation energy during ischemic VF.

Tunnel propagation following defibrillation with ICD shocks: hidden postshock activations in the left ventricular wall underlie isoelectric window

BACKGROUND

After near-defibrillation threshold (DFT) shocks from an implantable cardioverter-defibrillator (ICD), the first postshock activation that leads to defibrillation failure arises focally after an isoelectric window (IW). The mechanisms underlying the IW remain incompletely understood.

OBJECTIVE

The goal of this study was to provide mechanistic insight into the origins of postshock activations and IW after ICD shocks, and to link shock outcome to the preshock state of the ventricles. We hypothesized that the nonuniform ICD field results in the formation of an intramural excitable area (tunnel) only in the left ventricular (LV) free wall, through which both pre-existing and new shock-induced wavefronts propagate during the IW.

METHODS

Simulations were conducted using a realistic three dimensional (3D) model of defibrillation in the rabbit ventricles. Biphasic ICD shocks of varying strengths were delivered to 27 different fibrillatory states.

RESULTS

After near-DFT shocks, regardless of preshock state, the main postshock excitable area was always located within LV free wall, creating an intramural tunnel. Either pre-existing fibrillatory or shock-induced wavefronts propagated during the IW (duration of up to 74 ms) in this tunnel and emerged as breakthroughs on LV epicardium. Preshock activity within the LV played a significant role in shock outcome: a large number of preshock filaments resulted in an IW associated with tunnel propagation of pre-existing rather than shock-induced wavefronts. Furthermore, shocks were more likely to succeed if the LV excitable area was smaller.

CONCLUSION

The LV intramural excitable area is the primary reason for near-DFT failure. Any intervention that decreases the extent of this area will improve the likelihood of defibrillation success.

Generation of histo-anatomically representative models of the individual heart: tools and application

This paper presents methods to build histo-anatomically detailed individualized cardiac models. The models are based on high-resolution three-dimensional anatomical and/or diffusion tensor magnetic resonance images, combined with serial histological sectioning data, and are used to investigate individualized cardiac function. The current state of the art is reviewed, and its limitations are discussed. We assess the challenges associated with the generation of histo-anatomically representative individualized in silico models of the heart. The entire processing pipeline including image acquisition, image processing, mesh generation, model set-up and execution of computer simulations, and the underlying methods are described. The multifaceted challenges associated with these goals are highlighted, suitable solutions are proposed, and an important application of developed high-resolution structure-function models in elucidating the effect of individual structural heterogeneity upon wavefront dynamics is demonstrated.

Organization of ventricular fibrillation in the human heart: experiments and models

Sudden cardiac death is a major health problem in the industrialized world. The lethal event is typically ventricular fibrillation (VF), during which the co-ordinated regular contraction of the heart is overthrown by a state of mechanical and electrical anarchy. Understanding the excitation patterns that sustain VF is important in order to identify potential therapeutic targets. In this paper, we studied the organization of human VF by combining clinical recordings of electrical excitation patterns on the epicardial surface during in vivo human VF with simulations of VF in an anatomically and electrophysiologically detailed computational model of the human ventricles. We find both in the computational studies and in the clinical recordings that epicardial surface excitation patterns during VF contain around six rotors. Based on results from the simulated three-dimensional excitation patterns during VF, which show that the total number of electrical sources is 1.4 +/- 0.12 times greater than the number of epicardial rotors, we estimate that the total number of sources present during clinically recorded VF is 9.0 +/- 2.6. This number is approximately fivefold fewer compared with that observed during VF in dog and pig hearts, which are of comparable size to human hearts. We explain this difference by considering differences in action potential duration dynamics across these species. The simpler spatial organization of human VF has important implications for treatment and prevention of this dangerous arrhythmia. Moreover, our findings underline the need for integrated research, in which human-based clinical and computational studies complement animal research.

From mitochondrial ion channels to arrhythmias in the heart: computational techniques to bridge the spatio-temporal scales

Computer simulations of electrical behaviour in the whole ventricles have become commonplace during the last few years. The goals of this article are (i) to review the techniques that are currently employed to model cardiac electrical activity in the heart, discussing the strengths and weaknesses of the various approaches, and (ii) to implement a novel modelling approach, based on physiological reasoning, that lifts some of the restrictions imposed by current state-of-the-art ionic models. To illustrate the latter approach, the present study uses a recently developed ionic model of the ventricular myocyte that incorporates an excitation-contraction coupling and mitochondrial energetics model. A paradigm to bridge the vastly disparate spatial and temporal scales, from subcellular processes to the entire organ, and from sub-microseconds to minutes, is presented. Achieving sufficient computational efficiency is the key to success in the quest to develop multiscale realistic models that are expected to lead to better understanding of the mechanisms of arrhythmia induction following failure at the organelle level, and ultimately to the development of novel therapeutic applications.

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.

Three distinct phases of VF during global ischemia in the isolated blood-perfused pig heart

Changes in ventricular fibrillation (VF) organization occurring after the onset of global ischemia are relevant to defibrillation and survival but remain poorly understood. We hypothesized that ischemia-specific dynamic instability of the action potential (AP) causes a loss of spatiotemporal periodicity of propagation and broadening of the electrocardiogram (ECG) frequency spectrum during VF in the ischemic myocardium. We recorded voltage-sensitive fluorescence of di-4-ANEPPS (anterior left ventricle, 35 x 35 mm, 64 x 64 pixels) and the volume-conducted ECG in six blood-perfused hearts during 10 min of VF and global ischemia. We used coefficient of variation (CV) to estimate variability of AP amplitude, AP duration, and diastolic interval (CV-APA, CV-APD, and CV-DI, respectively). We computed excitation median frequency (Median_F), spectral width of the AP and ECG (SpW-AP and SpW-ECG, respectively), wavebreak incidence (WBI), and recurrence of propagation direction (RPD). We found three distinct phases of local VF dynamics: "relatively periodic" (<or=1 min, high Median_F, moderate AP variability, high WBI, low RPD), "highly periodic" (1-2 min, reduced Median_F, low AP variability, low WBI, high RPD), and "aperiodic" (3-10 min, low Median_F, high AP variability, high WBI, low RPD). In one experiment, spontaneous conversion from the aperiodic to the highly periodic phase occurred after 5 min of ischemia. The SpW-ECG was correlated with SpW-AP, CV-APD, and CV-APA. We conclude that 1) at least three distinct phases of VF dynamics are present in our model, and 2) the newly described aperiodic phase is related to ischemia-specific dynamic instability of the AP shape, which underlies broadening of the ECG spectrum during VF evolution.

Organization of ventricular fibrillation in the human heart

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

Cardiac defibrillation and the role of mechanoelectric feedback in postshock arrhythmogenesis

Ventricular dilatation increases the defibrillation threshold (DFT). In order to elucidate the mechanisms responsible for this increase, the present article investigates changes in the postshock behavior of the myocardium upon stretch. A two-dimensional electro-mechanical model of cardiac tissue incorporating heterogeneous fiber orientation was used to explore the effect of sustained stretch on postshock behavior via (a) recruitment of mechanosensitive channels (MSC) and (b) tissue deformation and concomitant changes in tissue conductivities. Recruitment of MSC had no influence on vulnerability to electric shocks as compared to control, but increased the complexity of postshock VF patterns. Stretch-induced deformation and changes in tissue conductivities resulted in a decrease in vulnerability to electric shocks.

Mechanistic enquiry into the effect of increased pacing rate on the upper limit of vulnerability

The goal of this study is to investigate the mechanisms responsible for the increase in the upper limit of vulnerability (ULV; highest shock strength that induces arrhythmia) following the increase in pacing rate. To accomplish this goal, the study employs a three-dimensional bidomain finite element model of a slice through the canine ventricles. The preparation was paced eight times at a basic cycle length (BCL) of either 80 or 150ms followed by delivery of shocks of various strengths and timings. Our results demonstrate that the shock strength, which induced an arrhythmia 50% of the time, increased 20% for the faster pacing compared to the slower pacing. Analysis of the mechanisms underlying the increased vulnerability revealed that delayed post-shock activations originating in the tissue depths appear as breakthrough activations on the surfaces of the preparation following an isoelectric window (IW). However, the IW duration was consistently shorter in the faster-paced preparation. Consequently, breakthrough activations appeared on the surfaces of this preparation earlier, when the tissue was less recovered, resulting in higher probability of unidirectional block and reentry. This explains why shocks of the same strength were more likely to result in arrhythmia induction when delivered to a preparation that was rapidly paced.

Defibrillation of the heart: insights into mechanisms from modelling studies

Despite its critical role in restoring cardiac rhythm and thus in saving human life, cardiac defibrillation remains poorly understood. Further mechanistic inquiry is hampered by the inability of presently available experimental techniques to resolve, with sufficient accuracy, electrical behaviour confined to the depth of the ventricles. The objective of this review article is to demonstrate that realistic 3-D simulations of the ventricular defibrillation process in close conjunction with experimental observations are capable of bringing a new level of understanding of the electrical events that ensue from the interaction between fibrillating myocardium and applied shock. The article does this by reviewing the results of two studies, one on vulnerability to electric shocks and another on defibrillation. An overview of the modelling tools used in these studies is also provided.

Defibrillation Depends on Conductivity Fluctuations and the Degree of Disorganization in Reentry Patterns

Defibrillation depends on conductivity and disorganization.

INTRODUCTION

Cardiac fibrillation is the deterioration of the heart's normally well-organized activity into one or more meandering spiral waves, which subsequently break up into many meandering wave fronts. Delivery of an electric shock (defibrillation) is the only effective way of restoring the normal rhythm. This study focuses on examining whether higher degrees of disorganization requires higher shock strengths to defibrillate and whether microscopic conductivity fluctuations favor shock success.

METHODS AND RESULTS

We developed a three-dimensional computer bidomain model of a block of cardiac tissue with straight fibers immersed in a conductive bath. The membrane behavior was described by the Courtemanche human atrial action potential model incorporating electroporation and an acetylcholine- (ACh) dependent potassium current. Intracellular conductivities were varied stochastically around nominal values with variations of up to 50%. A single rotor reentry was initiated and, by adjusting the spatial ACh variation, the level of organization could be controlled. The single rotor could be stabilized or spiral wave breakup could be provoked leading to fibrillatory-like activity. For each level of organization, multiple shock timings and strengths were applied to compute the probability of shock success as a function of shock strength.

CONCLUSIONS

Our results suggest that the level of the small-scale conductivity fluctuations is a very important factor in defibrillation. A higher variation significantly lowers the required shock strength. Further, we demonstrated that success also heavily depends on the level of organization of the fibrillatory episode. In general, higher levels of disorganization require higher shock strengths to defibrillate.

Computer simulations of successful defibrillation in decoupled and non-uniform cardiac tissue

Aim

The aim of the present study is to investigate the origin and effect of virtual electrode polarization in uniform, decoupled and non-uniform cardiac tissue during field stimulation.

Methods

A discrete bidomain model with active membrane behaviour was used to simulate normal cardiac tissue as well as cardiac tissue that is decoupled due to fibrosis and gap junction remodelling. Various uniform and non-uniform electric fields were applied to the external domain of uniform, decoupled and non-uniform resting cardiac tissue as well as cardiac tissue in which spiral waves were induced.

Results

Field stimulation applied on non-uniform tissue results in more virtual electrodes compared with uniform tissue. The spiral waves were terminated in decoupled tissue, but not in uniform, homogeneous tissue. By gradually increasing local differences in intracellular conductivities, the amount and spread of virtual electrodes increased and the spiral waves were terminated.

Conclusion

Fast depolarization of the tissue after field stimulation may be explained by intracellular decoupling and spatial heterogeneity present in normal and pathological cardiac tissue. We demonstrated that termination of spiral waves by means of field stimulation can be achieved when the tissue is modelled as a non-uniform, anisotropic bidomain with active membrane behaviour.

Effect of acute global ischemia on the upper limit of vulnerability: a simulation study

The goal of this modeling research is to provide mechanistic insight into the effect of altered membrane kinetics associated with 5-12 min of acute global ischemia on the upper limit of cardiac vulnerability (ULV) to electric shocks. We simulate electrical activity in a finite-element bidomain model of a 4-mm-thick slice through the canine ventricles that incorporates realistic geometry and fiber architecture. Global acute ischemia is represented by changes in membrane dynamics due to hyperkalemia, acidosis, and hypoxia. Two stages of acute ischemia are simulated corresponding to 5-7 min (stage 1) and 10-12 min (stage 2) after the onset of ischemia. Monophasic shocks are delivered in normoxia and ischemia over a range of coupling intervals, and their outcomes are examined to determine the highest shock strength that resulted in induction of reentrant arrhythmia. Our results demonstrate that acute ischemia stage 1 results in ULV reduction to 0.8A from its normoxic value of 1.4A. In contrast, no arrhythmia is induced regardless of shock strength in acute ischemia stage 2. An investigation of mechanisms underlying this behavior revealed that decreased postshock refractoriness resulting mainly from 1) ischemic electrophysiological substrate and 2) decrease in the extent of areas positively-polarized by the shock is responsible for the change in ULV during stage 1. In contrast, conduction failure is the main cause for the lack of vulnerability in acute ischemia stage 2. The insight provided by this study furthers our understanding of mechanisms by which acute ischemia-induced changes at the ionic level modulate cardiac vulnerability to electric shocks.