Ventricular fibrillation: mechanisms of initiation and maintenance

Noise-induced spatiotemporal patterns in Hodgkin-Huxley neuronal network

The effect of noise on the pattern selection in a regular network of Hodgkin-Huxley neurons is investigated, and the transition of pattern in the network is measured from subexcitable to excitable media. Extensive numerical results confirm that kinds of travelling wave such as spiral wave, circle wave and target wave could be developed and kept alive in the subexcitable network due to the noise. In the case of excitable media under noise, the developed spiral wave and target wave could coexist and new target-like wave is induced near to the border of media. The averaged membrane potentials over all neurons in the network are calculated to detect the periodicity of the time series and the generated traveling wave. Furthermore, the firing probabilities of neurons in networks are also calculated to analyze the collective behavior of networks.

Ventricular fibrillation conduction through an isthmus of preserved myocardium between radiofrequency lesions

BACKGROUND

Selective local acceleration of myocardial activation during ventricular fibrillation (VF) contributes information on the interactions between neighboring zones during the arrhythmia. This study analyzes these interactions, centering the observations on an isthmus of myocardium between two radiofrequency (RF) lesions.

METHODS

In nine isolated rabbit hearts, a gap of preserved myocardium was established between two RF lesions in the anterolateral left ventricle (LV) wall. Before, during, and after increasing the spatial heterogeneity of VF by local myocardial stretching, VF epicardial recordings were obtained.

RESULTS

Local stretch in the anterior LV wall decreased the excitable window (17 ± 7 ms vs 26 ± 7 ms; P < 0.05) and increased the dominant frequency (DFr; 18.9 ± 5.0 Hz vs 15.2 ± 3.6 Hz; P < 0.05) in this zone, without changes in the non-stretched posterolateral zone (25 ± 4 ms vs 27 ± 6 ms, ns and 14.1 ± 2.7 Hz vs 14.3 ± 3.0 Hz, ns). The DFr ratio at both sides of the gap was inversely correlated to the excitable window ratio (R = -0.57; P = 0.002). Before (31% vs 26%), during (29% vs 22%), and after stretch suppression (35% vs 25%), the wavefronts passing through the gap from the posterolateral to the anterior LV wall were seen to predominate. The number of wavefronts that passed from the anterior to the posterolateral LV wall was related to the excitable window in this zone (R = 0.41; P = 0.03).

CONCLUSIONS

The VF acceleration induced in the stretched zone does not increase the flow of wavefronts toward the non-stretched zone in the adjacent gap of preserved myocardium. The absence of significant changes in the electrophysiological parameters of the non-stretched myocardium limits the arrival of wavefronts in this zone.

Spatial gradients in action potential duration created by regional magnetofection of hERG are a substrate for wavebreak and turbulent propagation in cardiomyocyte monolayers

Spatial dispersion of action potential duration (APD) is a substrate for the maintenance of cardiac fibrillation, but the mechanisms are poorly understood. We investigated the role played by spatial APD dispersion in fibrillatory dynamics. We used an in vitro model in which spatial gradients in the expression of ether-à-go-go-related (hERG) protein, and thus rapid delayed rectifying K(+) current (I(Kr)) density, served to generate APD dispersion, high-frequency rotor formation, wavebreak and fibrillatory conduction. A unique adenovirus-mediated magnetofection technique generated well-controlled gradients in hERG and green fluorescent protein (GFP) expression in neonatal rat ventricular myocyte monolayers. Computer simulations using a realistic neonatal rat ventricular myocyte monolayer model provided crucial insight into the underlying mechanisms. Regional hERG overexpression shortened APD and increased rotor incidence in the hERG overexpressing region. An APD profile at 75 percent repolarization with a 16.6 ± 0.72 ms gradient followed the spatial profile of hERG-GFP expression; conduction velocity was not altered. Rotors in the infected region whose maximal dominant frequency was 12.9 Hz resulted in wavebreak at the interface (border zone) between infected and non-infected regions; dominant frequency distribution was uniform when the maximal dominant frequency was <12.9 Hz or the rotors resided in the uninfected region. Regularity at the border zone was lowest when rotors resided in the infected region. In simulations, a fivefold regional increase in I(Kr) abbreviated the APD and hyperpolarized the resting potential. However, the steep APD gradient at the border zone proved to be the primary mechanism of wavebreak and fibrillatory conduction. This study provides insight at the molecular level into the mechanisms by which spatial APD dispersion contributes to wavebreak, rotor stabilization and fibrillatory conduction.

Inhibition of intercellular coupling stabilizes spiral-wave reentry, whereas enhancement of the coupling destabilizes the reentry in favor of early termination

Spiral-wave (SW) reentry is a major organizing principle of ventricular tachycardia/fibrillation (VT/VF). We tested a hypothesis that pharmacological modification of gap junction (GJ) conductance affects the stability of SW reentry in a two-dimensional (2D) epicardial ventricular muscle layer prepared by endocardial cryoablation of Langendorff-perfused rabbit hearts. Action potential signals were recorded and analyzed by high-resolution optical mapping. Carbenoxolone (CBX; 30 μM) and rotigaptide (RG, 0.1 μM) were used to inhibit and enhance GJ coupling, respectively. CBX decreased the space constant (λ) by 36%, whereas RG increased it by 22-24% (n = 5; P < 0.01). During centrifugal propagation, there was a linear relationship between the wavefront curvature (κ) and local conduction velocity (LCV): LCV = LCV(0) - D·κ (D, diffusion coefficient; LCV(0), LCV at κ = 0). CBX decreased LCV(0) and D by 27 ± 3 and 57 ± 3%, respectively (n = 5; P < 0.01). RG increased LCV(0) and D by 18 ± 3 and 54 ± 5%, respectively (n = 5, P < 0.01). The regression lines with and without RG crossed, resulting in a paradoxical decrease of LCV with RG at κ > ~60 cm(-1). SW reentry induced after CBX was stable, and the incidence of sustained VTs (>30 s) increased from 38 ± 4 to 85 ± 4% after CBX (n = 18; P < 0.01). SW reentry induced after RG was characterized by decremental conduction near the rotation center, prominent drift and self-termination by collision with the anatomical boundaries, and the incidence of sustained VTs decreased from 40 ± 5 to 17 ± 6% after RG (n = 13; P < 0.05). These results suggest that decreased intercellular coupling stabilizes SW reentry in 2D cardiac muscle, whereas increased coupling facilitates its early self-termination.

Mechanisms of ventricular arrhythmias: a dynamical systems-based perspective

Defining the cellular electrophysiological mechanisms for ventricular tachyarrhythmias is difficult, given the wide array of potential mechanisms, ranging from abnormal automaticity to various types of reentry and kk activity. The degree of difficulty is increased further by the fact that any particular mechanism may be influenced by the evolving ionic and anatomic environments associated with many forms of heart disease. Consequently, static measures of a single electrophysiological characteristic are unlikely to be useful in establishing mechanisms. Rather, the dynamics of the electrophysiological triggers and substrates that predispose to arrhythmia development need to be considered. Moreover, the dynamics need to be considered in the context of a system, one that displays certain predictable behaviors, but also one that may contain seemingly stochastic elements. It also is essential to recognize that even the predictable behaviors of this complex nonlinear system are subject to small changes in the state of the system at any given time. Here we briefly review some of the short-, medium-, and long-term alterations of the electrophysiological substrate that accompany myocardial disease and their potential impact on the initiation and maintenance of ventricular arrhythmias. We also provide examples of cases in which small changes in the electrophysiological substrate can result in rather large differences in arrhythmia outcome. These results suggest that an interrogation of cardiac electrical dynamics is required to provide a meaningful assessment of the immediate risk for arrhythmia development and for evaluating the effects of putative antiarrhythmic interventions.

Slowed propagation across the compacta-trabeculata interface: a consequence of fiber and sheet anisotropy

Transmural myocardial activation is influenced by myocardial structure, including structural differences between the compacta (Cta) and the trabeculata (Tta), although this has not been fully explained. Hearts from rats were Langendorff perfused, stained with DI-4-ANEPPS, the apex was cut off and fluorescence acquired from the exposed short-axis surface. The hearts were stimulated at 160 ms cycle length at the anterior, lateral, posterior left ventricle (LV) and septal sub-epicardial sites. Conduction velocity perpendicular to the wave front orientation was measured in each pixel using a gradient-based approach. After optical mapping the cut surface was imaged using a light microscope and the extent of the Cta and Tta mapped and validated against 50 u, m isotropic MRI images. We used a 3D rat ventricle computational model, with architecture obtained from 200 u, m isotropic diffusion tensor MRI and kinetics from the modified Pandit model to determine the relative roles of fibers and sheets on propagation. We show in the experimental study that circumferential propagation around the LV cavity is fast in the Cta: 63.2±19.5 and is slower in the Tta: 32.7±11.0(∗) (mean ± s.d cms-1, ∗ p < 0.01 by two sample t test). In the simulation study the pattern and velocity are not replicated in an isotropic model (I), are partially replicated in a simulation study including fiber anisotropy (A) and is more fully replicated in orthotropic (O) ventricles (fiber and sheet anisotropy), where the circumferential propagation velocity is, I: Cta: 54.2±3.9; Tta:54.3±3.9; A: Cta:43.6±3.2; Tta: 40.6±6.6; O: Cta: 63.2±19.5; Tta: 32.7±11.9(∗). We show that sheet orientation is important in understanding activation differences between Cta and Tta.

Rigidity matching between cells and the extracellular matrix leads to the stabilization of cardiac conduction

Biomechanical dynamic interactions between cells and the extracellular environment dynamically regulate physiological tissue behavior in living organisms, such as that seen in tissue maintenance and remodeling. In this study, the substrate-induced modulation of synchronized beating in cultured cardiomyocyte tissue was systematically characterized on elasticity-tunable substrates to elucidate the effect of biomechanical coupling. We found that myocardial conduction is significantly promoted when the rigidity of the cell culture environment matches that of the cardiac cells (4 kiloPascals). The stability of spontaneous target wave activity and calcium transient alternans in high frequency-paced tissue were both enhanced when the cell substrate and cell tissue showed the same rigidity. By adapting a simple theoretical model, we reproduced the experimental trend on the rigidity matching for the synchronized excitation. We conclude that rigidity matching in cell-to-substrate interactions critically improves cardiomyocyte-tissue synchronization, suggesting that mechanical coupling plays an essential role in the dynamic activity of the beating heart.

The statistics of calcium-mediated focal excitations on a one-dimensional cable

It is well known that various cardiac arrhythmias are initiated by an ill-timed excitation that originates from a focal region of the heart. However, up to now, it is not known what governs the timing, location, and morphology of these focal excitations. Recent studies have shown that these excitations can be caused by abnormalities in the calcium (Ca) cycling system. However, the cause-and-effect relationships linking subcellular Ca dynamics and focal activity in cardiac tissue is not completely understood. In this article, we present a minimal model of Ca-mediated focal excitations in cardiac tissue. This model accounts for the stochastic nature of spontaneous Ca release on a one-dimensional cable of cardiac cells. Using this model, we show that the timing of focal excitations is equivalent to a first passage time problem in a spatially extended system. In particular, we find that for a short cable the mean first passage time increases exponentially with the number of cells in tissue, and is critically dependent on the ratio of inward to outward currents near the threshold for an action potential. For long cables excitations occurs due to ectopic foci that occur on a length scale determined by the minimum length of tissue that can induce an action potential. Furthermore, we find that for long cables the mean first passage time decreases as a power law in the number cells. These results provide precise criteria for the occurrence of focal excitations in cardiac tissue, and will serve as a guide to determine the propensity of Ca-mediated triggered arrhythmias in the heart.

Control of turbulence in heterogeneous excitable media

Control of turbulence in two kinds of typical heterogeneous excitable media by applying a combined method is investigated. It is found that local-low-amplitude and high-frequency pacing (LHP) is effective to suppress turbulence if the deviation of the heterogeneity is minor. However, LHP is invalid when the deviation is large. Studies show that an additional radial electric field can greatly increase the efficiency of LHP. The underlying mechanisms of successful control in the two kinds of cases are different and are discussed separately. Since the developed strategy of combining LHP with a radial electric field can terminate turbulence in excitable media with a high degree of inhomogeneity, it has the potential contribution to promote the practical low-amplitude defibrillation approach.

The role of photon scattering in voltage-calcium fluorescent recordings of ventricular fibrillation

Recent optical mapping studies of cardiac tissue suggest that membrane voltage (V(m)) and intracellular calcium concentrations (Ca) become dissociated during ventricular fibrillation (VF), generating a proarrhythmic substrate. However, experimental methods used in these studies may accentuate measured dissociation due to differences in fluorescent emission wavelengths of optical voltage/calcium (V(opt)/Ca(opt)) signals. Here, we simulate dual voltage-calcium optical mapping experiments using a monodomain-Luo-Rudy ventricular-tissue model coupled to a photon-diffusion model. Dissociation of both electrical, V(m)/Ca, and optical, V(opt)/Ca(opt), signals is quantified by calculating mutual information (MI) for VF and rapid pacing protocols. We find that photon scattering decreases MI of V(opt)/Ca(opt) signals by 23% compared to unscattered V(m)/Ca signals during VF. Scattering effects are amplified by increasing wavelength separation between fluorescent voltage/calcium signals and respective measurement-location misalignment. In contrast, photon scattering does not affect MI during rapid pacing, but high calcium dye affinity can decrease MI by attenuating alternans in Ca(opt) but not in V(opt). We conclude that some dissociation exists between voltage and calcium at the cellular level during VF, but MI differences are amplified by current optical mapping methods.

New index for categorising cardiac reentrant wave: in silico evaluation

Based on the similarity between a reentrant wave in cardiac tissue and a vortex in fluid dynamics, the authors hypothesised that a new non-dimensional index, like the Reynolds number in fluid dynamics, may play a critical role in categorising reentrant wave dynamics. Therefore the goal of the present study is to devise a new index to characterise electric wave conduction in cardiac tissue and examined whether this index can be used as a biomarker for categorising the reentrant wave pattern in cardiac tissue. Similar to the procedure used to derive the Reynolds number in fluid dynamics, the authors used a non-dimensionalisation technique to obtain the new index. Its usefulness was verified using a two-dimensional simulation model of electrical wave propagation in cardiac tissue. The simulation results showed that electrical waves in cardiac tissue move into an unstable region when the index exceeds a threshold value.

Regional cooling facilitates termination of spiral-wave reentry through unpinning of rotors in rabbit hearts

BACKGROUND

Moderate global cooling of myocardial tissue was shown to destabilize 2-dimensional (2-D) reentry and facilitate its termination.

OBJECTIVE

This study sought to test the hypothesis that regional cooling destabilizes rotors and facilitates termination of spontaneous and DC shock-induced subepicardial reentry in isolated, endocardially ablated rabbit hearts.

METHODS

Fluorescent action potential signals were recorded from 2-D subepicardial ventricular myocardium of Langendorff-perfused rabbit hearts. Regional cooling (by 5.9°C ± 1.3°C) was applied to the left ventricular anterior wall using a transparent cooling device (10 mm in diameter).

RESULTS

Regional cooling during constant stimulation (2.5 Hz) prolonged the action potential duration (by 36% ± 9%) and slightly reduced conduction velocity (by 4% ± 4%) in the cooled region. Ventricular tachycardias (VTs) induced during regional cooling terminated earlier than those without cooling (control): VTs lasting >30 seconds were reduced from 17 of 39 to 1 of 61. When regional cooling was applied during sustained VTs (>120 seconds), 16 of 33 (48%) sustained VTs self-terminated in 12.5 ± 5.1 seconds. VT termination was the result of rotor destabilization, which was characterized by unpinning, drift toward the periphery of the cooled region, and subsequent collision with boundaries. The DC shock intensity required for cardioversion of the sustained VTs decreased significantly by regional cooling (22.8 ± 4.1 V, n = 16, vs 40.5 ± 17.6 V, n = 21). The major mode of reentry termination by DC shocks was phase resetting in the absence of cooling, whereas it was unpinning in the presence of cooling.

CONCLUSION

Regional cooling facilitates termination of 2-D reentry through unpinning of rotors.

Ventricular fibrillation: dynamics and ion channel determinants

Ventricular fibrillation (VF) is the leading cause of sudden cardiac death. This brief review addresses issues relevant to the dynamics of the rotors responsible for functional reentry and VF. It also makes an attempt to summarize present-day knowledge of the manner in which the dynamic interplay between inward and outward transmembrane currents and the heterogeneous cardiac structure establish a substrate for the initiation and maintenance of rotors and VF. The fragmentary nature of our current understanding of ionic VF mechanisms does not even allow an approach toward a "Theory of VF". Yet some hope is provided by recently obtained insight into the roles played in VF by some of the sarcolemmal ion channels that control the excitation-recovery process. For example, strong evidence supports the idea that the interplay between the rapid-inward sodium current and the inward-rectifier potassium current controls rotor formation, as well as rotor stability and frequency. Solid evidence also exists for an involvement of L-type calcium current in the control of rotor frequency and in determining VF-to-ventricular tachycardia conversion. Less clear, however, is whether or not time dependent outward currents through voltage-gated potassium channels affect the fibrillatory process. Hopefully, taking advantage of currently available approaches of structural, molecular and cellular biology, together with computational and imaging techniques, will afford us the opportunity to further advance knowledge on VF mechanisms.

Quantitative prediction of the arrhythmogenic effects of de novo hERG mutations in computational models of human ventricular tissues

Mutations to hERG which result in changes to the rapid delayed rectifier current I(Kr) can cause long and short QT syndromes and are associated with an increased risk of cardiac arrhythmias. Experimental recordings of I(Kr) reveal the effects of mutations at the channel level, but how these changes translate to the cell and tissue levels remains unclear. We used computational models of human ventricular myocytes and tissues to predict and quantify the effects that de novo hERG mutations would have on cell and tissue electrophysiology. Mutations that decreased I(Kr) maximum conductance resulted in an increased cell and tissue action potential duration (APD) and a long QT interval on the electrocardiogram (ECG), whereas those that caused a positive shift in the inactivation curve resulted in a decreased APD and a short QT. Tissue vulnerability to re-entrant arrhythmias was correlated with transmural dispersion of repolarisation, and any change to this vulnerability could be inferred from the ECG QT interval or T wave peak-to-end time. Faster I(Kr) activation kinetics caused cell APD alternans to appear over a wider range of pacing rates and with a larger magnitude, and spatial heterogeneity in these cellular alternans resulted in discordant alternans at the tissue level. Thus, from channel kinetic data, we can predict the tissue-level electrophysiological effects of any hERG mutations and identify how the mutation would manifest clinically, as either a long or short QT syndrome with or without an increased risk of alternans and re-entrant arrhythmias.

Bepridil facilitates early termination of spiral-wave reentry in two-dimensional cardiac muscle through an increase of intercellular electrical coupling

Bepridil is effective for conversion of atrial fibrillation to sinus rhythm and in the treatment of drug-refractory ventricular tachyarrhythmias. We investigated the effects of bepridil on electrophysiological properties and spiral-wave (SW) reentry in a 2-dimensional ventricular muscle layer of isolated rabbit hearts by optical mapping. Ventricular tachycardia (VT) induced in the presence of bepridil (1 µM) terminated earlier than in the control. Bepridil increased action potential duration (APD) by 5% - 8% under constant pacing and significantly increased the space constant. There was a linear relationship between the wavefront curvature (κ) and local conduction velocity: LCV = LCV₀ - D·κ (D, diffusion coefficient; LCV₀, LCV at κ = 0). Bepridil significantly increased D and LCV₀. The regression lines with and without bepridil crossed at κ = 20 - 40 cm⁻¹, resulting in a paradoxical decrease of LCV at κ > 40 cm⁻¹. Dye transfer assay in cultured rat cardiomyocytes confirmed that bepridil increased intercellular coupling. SW reentry in the presence of bepridil was characterized by decremental conduction near the rotation center, prominent drift, and self-termination by collision with boundaries. These results indicate that bepridil causes an increase of intercellular coupling and a moderate APD prolongation, and this combination compromises wavefront propagation near the rotation center of SW reentry, leading to its drift and early termination.

Construction and validation of anisotropic and orthotropic ventricular geometries for quantitative predictive cardiac electrophysiology

Reaction-diffusion computational models of cardiac electrophysiology require both dynamic excitation models that reconstruct the action potentials of myocytes as well as datasets of cardiac geometry and architecture that provide the electrical diffusion tensor D, which determines how excitation spreads through the tissue. We illustrate an experimental pipeline we have developed in our laboratories for constructing and validating such datasets. The tensor D changes with location in the myocardium, and is determined by tissue architecture. Diffusion tensor magnetic resonance imaging (DT-MRI) provides three eigenvectors e(i) and eigenvalues λ(i) at each voxel throughout the tissue that can be used to reconstruct this architecture. The primary eigenvector e(1) is a histologically validated measure of myocyte orientation (responsible for anisotropic propagation). The secondary and tertiary eigenvectors (e(2) and e(3)) specify the directions of any orthotropic structure if λ(2) is significantly greater than λ(3)-this orthotropy has been identified with sheets or cleavage planes. For simulations, the components of D are scaled in the fibre and cross-fibre directions for anisotropic simulations (or fibre, sheet and sheet normal directions for orthotropic tissues) so that simulated conduction velocities match values from optical imaging or plunge electrode experiments. The simulated pattern of propagation of action potentials in the models is partially validated by optical recordings of spatio-temporal activity on the surfaces of hearts. We also describe several techniques that enhance components of the pipeline, or that allow the pipeline to be applied to different areas of research: Q ball imaging provides evidence for multi-modal orientation distributions within a fraction of voxels, infarcts can be identified by changes in the anisotropic structure-irregularity in myocyte orientation and a decrease in fractional anisotropy, clinical imaging provides human ventricular geometry and can identify ischaemic and infarcted regions, and simulations in human geometries examine the roles of anisotropic and orthotropic architecture in the initiation of arrhythmias.

Alcohol ablation at the posterior papillary muscle prevents ventricular fibrillation in swine without affecting mitral valve function

AIMS

Radiofrequency ablation at the posterior papillary muscle (PM) significantly reduced ventricular fibrillation (VF) inducibility in rabbits and dogs, suggesting that PM may be involved in the generation of VF. However, the effect of ablation at the PM on VF inducibility remains unknown in normal intact swine hearts because in this species radiofrequency energy delivered at PM provoked incessant VF.

METHODS AND RESULTS

Twelve anesthetized swine underwent median sternotomy. Under the ultrasonographic guidance, chemical ablation was performed via injection of dehydrated alcohol into the base of the posterior PM (group PM, n = 6) or anterior wall (control group, n = 6) in the left ventricle. Ventricular fibrillation inducibility and mitral valve function were measured pre- and post-ablation. Hearts were explanted and the ablated myocardium was stained with haematoxylin and eosin. Ventricular fibrillation inducibility was significantly decreased from 100 ± 0% pre-ablation to 11.9 ± 7.8% post-ablation in group PM (P = 0.001), whereas it was not statistically different in the control group (100 ± 0 vs. 92.9 ± 7.1%, pre-ablation vs. post-ablation). Haemorrhage and cellular necrosis was observed in the centre of ablated myocardium and no significant mitral regurgitation was observed following ablation at the posterior PM.

CONCLUSION

Alcohol ablation of the left posterior PM reduced VF inducibility in normal intact swine hearts, with no significant mitral regurgitation. This suggests that the posterior PM may be involved in the generation of VF, and the recurrence of VF may be prevented by chemical ablation at the posterior PM.

Déjà vu in the theories of atrial fibrillation dynamics

This brief review looks back to the major theoretical, experimental, and clinical work on the dynamics and mechanisms of atrial fibrillation (AF). Its goal is to highlight the most important issues, controversies, and advances that have driven the field of investigation into AF mechanisms at any given time during the last ∼100 years. It emphasizes that while the history of AF research has been full of controversies from the start, such controversies have led to new information, and individual scientists have learned from those that have preceded them. However, in the face of the most common sustained cardiac arrhythmia seen in clinical practice, we are yet to fully understand its fundamental mechanisms and learn how to treat it effectively. Future research into AF dynamics and mechanisms should focus on the development and validation of new numerical and animal models. Such models should be relevant to and accurately reproduce the important substrates associated with ageing and with diseases such as hypertension, heart failure, and ischaemic heart disease which cause AF in the vast majority of patients. Knowledge derived from such models may help to greatly advance the field and hopefully lead to more effective prevention and therapy.

A major role for HERG in determining frequency of reentry in neonatal rat ventricular myocyte monolayer

RATIONALE

the rapid delayed rectifier potassium current, I(Kr), which flows through the human ether-a-go-go-related (hERG) channel, is a major determinant of the shape and duration of the human cardiac action potential (APD). However, it is unknown whether the time dependency of I(Kr) enables it to control APD, conduction velocity (CV), and wavelength (WL) at the exceedingly high activation frequencies that are relevant to cardiac reentry and fibrillation.

OBJECTIVE

to test the hypothesis that upregulation of hERG increases functional reentry frequency and contributes to its stability.

METHODS AND RESULTS

using optical mapping, we investigated the effects of I(Kr) upregulation on reentry frequency, APD, CV, and WL in neonatal rat ventricular myocyte (NRVM) monolayers infected with GFP (control), hERG (I(Kr)), or dominant negative mutant hERG G628S. Reentry frequency was higher in the I(Kr)-infected monolayers (21.12 ± 0.8 Hz; n=43 versus 9.21 ± 0.58 Hz; n=16; P<0.001) but slightly reduced in G628S-infected monolayers. APD(80) in the I(Kr)-infected monolayers was shorter (>50%) than control during pacing at 1 to 5 Hz. CV was similar in both groups at low frequency pacing. In contrast, during high-frequency reentry, the CV measured at varying distances from the center of rotation was significantly faster in I(Kr)-infected monolayers than controls. Simulations using a modified NRVM model predicted that rotor acceleration was attributable, in part, to a transient hyperpolarization immediately following the AP. The transient hyperpolarization was confirmed experimentally.

CONCLUSIONS

hERG overexpression dramatically accelerates reentry frequency in NRVM monolayers. Both APD and WL shortening, together with transient hyperpolarization, underlies the increased rotor frequency and stability.

Collision-based spiral acceleration in cardiac media: roles of wavefront curvature and excitable gap

We have previously shown in experimental cardiac cell monolayers that rapid point pacing can convert basic functional reentry (single spiral) into a stable multiwave spiral that activates the tissue at an accelerated rate. Here, our goal is to further elucidate the biophysical mechanisms of this rate acceleration without the potential confounding effects of microscopic tissue heterogeneities inherent to experimental preparations. We use computer simulations to show that, similar to experimental observations, single spirals can be converted by point stimuli into stable multiwave spirals. In multiwave spirals, individual waves collide, yielding regions with negative wavefront curvature. When a sufficient excitable gap is present and the negative-curvature regions are close to spiral tips, an electrotonic spread of excitatory currents from these regions propels each colliding spiral to rotate faster than the single spiral, causing an overall rate acceleration. As observed experimentally, the degree of rate acceleration increases with the number of colliding spiral waves. Conversely, if collision sites are far from spiral tips, excitatory currents have no effect on spiral rotation and multiple spirals rotate independently, without rate acceleration. Understanding the mechanisms of spiral rate acceleration may yield new strategies for preventing the transition from monomorphic tachycardia to polymorphic tachycardia and fibrillation.

Mechanisms underlying the antifibrillatory action of hyperkalemia in Guinea pig hearts

Hyperkalemia increases the organization of ventricular fibrillation (VF) and may also terminate it by mechanisms that remain unclear. We previously showed that the left-to-right heterogeneity of excitation and wave fragmentation present in fibrillating guinea pig hearts is mediated by chamber-specific outward conductance differences in the inward rectifier potassium current (I(K1)). We hypothesized that hyperkalemia-mediated depolarization of the reversal potential of I(K1) (E(K1)) would reduce excitability and thereby reduce VF excitation frequencies and left-to-right heterogeneity. We induced VF in Langendroff-perfused guinea pig hearts and increased the extracellular K(+) concentration ([K(+)](o)) from control (4 mM) to 7 mM (n = 5) or 10 mM (n = 7). Optical mapping enabled spatial characterization of excitation dominant frequencies (DFs) and wavebreaks, and identification of sustained rotors (>4 cycles). During VF, hyperkalemia reduced the maximum DF of the left ventricle (LV) from 31.5 +/- 4.7 Hz (control) to 23.0 +/- 4.7 Hz (7.0 mM) or 19.5 +/- 3.6 Hz (10.0 mM; p < 0.006), the left-to-right DF gradient from 14.7 +/- 3.6 Hz (control) to 4.4 +/- 1.3 Hz (7 mM) and 3.2 +/- 1.4 Hz (10 mM), the number of DF domains, and the incidence of wavebreak in the LV and interventricular regions. During 10 mM [K(+)](o), the rotation period and core area of sustained rotors in the LV increased, and VF often terminated. Two-dimensional computer simulations mimicking experimental VF predicted that clamping E(K1) to normokalemic values during simulated hyperkalemia prevented all of the hyperkalemia-induced VF changes. During hyperkalemia, despite the shortening of the action potential duration, depolarization of E(K1) increased refractoriness, leading to a slowing of VF, which effectively superseded the influence of I(K1) conductance differences on VF organization. This reduced the left-to-right excitation gradients and heterogeneous wavebreak formation. Overall, these results provide, to our knowledge, the first direct mechanistic insight into the organization and/or termination of VF by hyperkalemia.

Acute amiodarone promotes drift and early termination of spiral wave re-entry

Intravenous application of amiodarone is commonly used in the treatment of life-threatening arrhythmias, but the underlying mechanism is not fully understood. The purpose of the present study is to investigate the acute effects of amiodarone on spiral wave (SW) re-entry, the primary organization machinery of ventricular tachycardia/fibrillation (VT/VF), in comparison with lidocaine. A two-dimensional ventricular myocardial layer was obtained from 24 Langendorff-perfused rabbit hearts, and epicardial excitations were analyzed by high-resolution optical mapping. During basic stimulation, amiodarone (5 microM) caused prolongation of action potential duration (APD) by 5.6%-9.1%, whereas lidocaine (15 microM) caused APD shortening by 5.0%-6.4%. Amiodarone and lidocaine reduced conduction velocity similarly. Ventricular tachycardias induced by DC stimulation in the presence of amiodarone were of shorter duration (sustained-VTs >30 s/total VTs: 2/58, amiodarone vs 13/52, control), whereas those with lidocaine were of longer duration (22/73, lidocaine vs 14/58, control). Amiodarone caused prolongation of VT cycle length and destabilization of SW re-entry, which is characterized by marked prolongation of functional block lines, frequent wavefront-tail interactions near the rotation center, and considerable drift, leading to its early annihilation via collision with anatomical boundaries. Spiral wave re-entry in the presence of lidocaine was more stabilized than in control. In the anisotropic ventricular myocardium, amiodarone destabilizes SW re-entry facilitating its early termination. Lidocaine, in contrast, stabilizes SW re-entry resulting in its persistence.

Repolarization of the cardiac action potential. Does an increase in repolarization capacity constitute a new anti-arrhythmic principle?

The cardiac action potential can be divided into five distinct phases designated phases 0-4. The exact shape of the action potential comes about primarily as an orchestrated function of ion channels. The present review will give an overview of ion channels involved in generating the cardiac action potential with special emphasis on potassium channels involved in phase 3 repolarization. In humans, these channels are primarily K(v)11.1 (hERG1), K(v)7.1 (KCNQ1) and K(ir)2.1 (KCNJ2) being the responsible alpha-subunits for conducting I(Kr), I(Ks) and I(K1). An account will be given about molecular components, biophysical properties, regulation, interaction with other proteins and involvement in diseases. Both loss and gain of function of these currents are associated with different arrhythmogenic diseases. The second part of this review will therefore elucidate arrhythmias and subsequently focus on newly developed chemical entities having the ability to increase the activity of I(Kr), I(Ks) and I(K1). An evaluation will be given addressing the possibility that this novel class of compounds have the ability to constitute a new anti-arrhythmic principle. Experimental evidence from in vitro, ex vivo and in vivo settings will be included. Furthermore, conceptual differences between the short QT syndrome and I(Kr) activation will be accounted for.

Atrial septopulmonary bundle of the posterior left atrium provides a substrate for atrial fibrillation initiation in a model of vagally mediated pulmonary vein tachycardia of the structurally normal heart

BACKGROUND

The posterior left atrium (PLA) and pulmonary veins (PVs) have been shown to be critical for atrial fibrillation (AF) initiation. However, the detailed mechanisms of reentry and AF initiation by PV impulses are poorly understood. We hypothesized that PV impulses trigger reentry and AF by undergoing wavebreaks as a result of sink-to-source mismatch at specific PV-PLA transitions along the septopulmonary bundle, where there are changes in thickness and fiber direction.

METHODS AND RESULTS

In 7 Langendorff-perfused sheep hearts AF was initiated by a burst of 6 pulses (CL 80 to 150ms) delivered to the left inferior or right superior PV ostium 100 to 150 ms after the sinus impulse in the presence of 0.5 micromol/L acetylcholine. The exposed septal-PLA endocardial area was mapped with high spatio-temporal resolution (DI-4-ANEPPS, 1000-fr/s) during AF initiation. Isochronal maps for each paced beat preceding AF onset were constructed to localize areas of conduction delay and block. Phase movies allowed the determination of the wavebreak sites at the onset of AF. Thereafter, the PLA myocardial wall thickness was quantified by echocardiography, and the fiber direction in the optical field of view was determined after peeling off the endocardium. Finally, isochrone, phase and conduction velocity maps were superimposed on the corresponding anatomic pictures for each of the 28 episodes of AF initiation. The longest delays of the paced PV impulses, as well as the first wavebreak, occurred at those boundaries along the septopulmonary bundle that showed sharp changes in fiber direction and the largest and most abrupt increase in myocardial thickness.

CONCLUSION

Waves propagating from the PVs into the PLA originating from a simulated PV tachycardia triggered reentry and vagally mediated AF by breaking at boundaries along the septopulmonary bundle where abrupt changes in thickness and fiber direction resulted in sink-to-source mismatch and low safety for propagation.

Inward rectifier potassium channels control rotor frequency in ventricular fibrillation

Ventricular fibrillation (VF) is the most important cause of sudden cardiac death. While traditionally thought to result from random activation of the ventricles by multiple independent wavelets, recent evidence suggests that VF may be determined by the sustained activation of a relatively small number of reentrant sources. In addition, recent experimental data in various species as well as computer simulations have provided important clues about its ionic and molecular mechanisms, particularly in regards to the role of potassium currents in such mechanisms. The results strongly argue that the inward rectifier current, I(K1,) is an important current during functional reentry because it mediates the electrotonic interactions between the unexcited core and its immediate surroundings. In addition, I(K1) is a stabilizer of reentry due to its ability to shorten action potential duration and reduce conduction velocity near the center of rotation. Increased I(K1) prevents wave front-wave tail interactions and thus averts rotor destabilization and breakup. Other studies have shown that while the slow component of the delayed rectifier potassium current I(Ks) does not significantly modify rotor frequency or stability, it plays a major role in postrepolarization refractoriness and wave break formation. Therefore, the interplay between I(K1) and the rapid sodium inward current (I(Na)) is a major factor in the control of cardiac excitability and thus the stability and frequency of reentry, while I(Ks) is an important determinant of fibrillatory conduction.

Arrhythmogenic mechanisms of the Purkinje system during electric shocks: a modeling study

BACKGROUND

The function of the Purkinje system (PS) is to ensure fast and uniform activation of the heart. Although this vital role during sinus rhythm is well understood, this is not the case when shocks are applied to the heart, especially in the case of failed defibrillation. The PS activates differently from the myocardium, has different electrophysiological properties, and provides alternate propagation pathways; thus, there are many ways in which it can contribute to postshock behavior.

OBJECTIVE

The purpose of this study was to elucidate the role of the PS in the initiation and maintenance of postshock arrhythmias.

METHODS

A computer model of the ventricles including the PS was subjected to different reentry induction protocols.

RESULTS

The PS facilitated reentry induction at relatively weaker shocks. Disconnecting the PS from the ventricles during the postshock interval revealed that the PS helps stabilize early-stage reentry by providing focal breakthroughs. During later stages, the PS contributed to reentry by leading to higher frequency rotors. The PS also promoted wave front splitting during reentry due to electrotonic coupling, which prolongs action potential durations at PS-myocyte junctions. The presence of a PS results in the anchoring of reentrant activations that propagate through the pathways provided by the PS.

CONCLUSIONS

The PS is proarrhythmic in that it provides pathways that prolong activity, and it plays a supplementary role in maintaining the later stages of reentry (>800 ms).

Model-based control of cardiac alternans on a ring

Cardiac alternans, a beat-to-beat alternation of cardiac electrical dynamics, and ventricular tachycardia, generally associated with a spiral wave of electrical activity, have been identified as frequent precursors of the life-threatening spatiotemporally chaotic electrical state of ventricular fibrillation (VF). Schemes for the elimination of alternans and the stabilization of spiral waves through the injection of weak external currents have been proposed as methods to prevent VF but have not performed at the level required for clinical implementation. In this paper we propose a control method based on linear-quadratic regulator (LQR) control. Unlike most previously proposed approaches, our method incorporates information from the underlying model to increase efficiency. We use a one-dimensional ringlike geometry, with a single control electrode, to compare the performance of our method with that of two other approaches, quasi-instantaneous suppression of unstable modes (QISUM) and time-delay autosynchronization (TDAS). We find that QISUM fails to suppress alternans due to conduction block. Although both TDAS and LQR succeed in suppressing alternans, LQR is able to suppress the alternans faster and using a much weaker control current. Our results highlight the benefits of a model-based control approach despite its inherent complexity compared with nonmodel-based control such as TDAS.

The relative role of refractoriness and source-sink relationship in reentry generation during simulated acute ischemia

During acute myocardial ischemia, reentrant episodes may lead to ventricular fibrillation (VF), giving rise to potentially mortal arrhythmias. VF has been traditionally related to dispersion of refractoriness and more recently to the source-sink relationship. Our goal is to theoretically investigate the relative role of dispersion of refractoriness and source-sink mismatch in vulnerability to reentry in the specific situation of regional myocardial acute ischemia. The electrical activity of a regionally ischemic tissue was simulated using a modified version of the Luo-Rudy dynamic model. Ischemic conditions were varied to simulate the time-course of acute ischemia. Our results showed that dispersion of refractoriness increased with the severity of ischemia. However, no correlation between dispersion of refractoriness and the width of the vulnerable window was found. Additionally, in approximately 50% of the reentries, unidirectional block (UDB) took place in cells completely recovered from refractoriness. We examined patterns of activation after premature stimulation and they were intimately related to the source-sink relationship, quantified by the safety factor (SF). Moreover, the isoline where the SF dropped below unity matched the area where propagation failed. It was concluded that the mismatch of the source-sink relationship, rather than solely refractoriness, was the ultimate cause of the UDB leading to reentry. The SF represents a very powerful tool to study the mechanisms responsible for reentry.

Suppression of spirals and turbulence in inhomogeneous excitable media

Suppression of spiral and turbulence in inhomogeneous media due to local heterogeneity with higher excitability is investigated numerically. When the inhomogeneity is small, control tactics by boundary periodic forcing (BPF) is effective against the existing spiral and turbulence. When the inhomogeneity of excitability is large, a rotating electric field (REF) is utilized to "smooth" regional heterogeneity based on driven synchronization. Consequently, a control approach combining BPF with REF is proposed to suppress the spiral and turbulence. The underlying mechanism of successful suppression is discussed in terms of dispersion relation.

Synchronization of chaotic early afterdepolarizations in the genesis of cardiac arrhythmias

The synchronization of coupled oscillators plays an important role in many biological systems, including the heart. In heart diseases, cardiac myocytes can exhibit abnormal electrical oscillations, such as early afterdepolarizations (EADs), which are associated with lethal arrhythmias. A key unanswered question is how cellular EADs partially synchronize in tissue, as is required for them to propagate. Here, we present evidence, from computational simulations and experiments in isolated myocytes, that irregular EAD behavior is dynamical chaos. We then show in electrically homogeneous tissue models that chaotic EADs synchronize globally when the tissue is smaller than a critical size. However, when the tissue exceeds the critical size, electrotonic coupling can no longer globally synchronize EADs, resulting in regions of partial synchronization that shift in time and space. These regional partially synchronized EADs then form premature ventricular complexes that propagate into recovered tissue without EADs. This process creates multiple premature ventricular complexes that propagate as [corrected] "shifting" foci resembling polymorphic ventricular tachycardia. Shifting foci encountering shifting repolarization gradients can also develop localized wave breaks leading to reentry and fibrillation. As predicted by the theory, rabbit hearts exposed to oxidative stress (H(2)O(2)) exhibited multiple shifting foci causing polymorphic tachycardia and fibrillation. This mechanism explains how collective cellular behavior integrates at the tissue scale to generate lethal cardiac arrhythmias over a wide range of heart rates.

A computational study of mother rotor VF in the human ventricles

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

Early termination of spiral wave reentry by combined blockade of Na+ and L-type Ca2+ currents in a perfused two-dimensional epicardial layer of rabbit ventricular myocardium

BACKGROUND

Modification of spiral wave (SW) reentry by antiarrhythmic drugs is a central issue to be challenged for better understanding of their benefits and risks.

OBJECTIVE

We investigated the effects of pilsicainide and/or verapamil, which block sodium and L-type calcium currents (I(Na) and I(Ca,L)), respectively, on SW reentry.

METHODS

A two-dimensional epicardial ventricular muscle layer was created in rabbit hearts by cryoablation (n = 32), and action potential signals were analyzed by high-resolution optical mapping.

RESULTS

During constant stimulation, pilsicainide (3-5 microM) caused a frequency-dependent decrease of conduction velocity (CV; by 20%-54% at 5 Hz) without affecting action potential duration (APD). Verapamil (3 microM) caused APD shortening (by 16% at 5 Hz) without affecting CV. Ventricular tachycardias (VTs) that were induced were more sustained in the presence of either pilsicainide or verapamil. The incidence of sustained VTs (>30 s)/all VTs per heart was 58% +/- 9% for 5 microM pilsicainide vs. 22% +/- 9% for controls and 62% +/- 10% for 3 microM verapamil vs. 22% +/- 8% for controls. The SWs with pilsicainide were characterized by slower rotation around longer functional block lines (FBLs), whereas those with verapamil were characterized by faster rotation around shorter FBLs. Combined application of 3 microM pilsicainide and 3 microM verapamil resulted in early termination of VTs (sustained VTs/all VTs per heart: 2% +/- 2% vs. 29% +/- 9% for controls); SWs showed extensive drift and decremental conduction, leading to their spontaneous annihilation.

CONCLUSION

Blockade of either I(Na) or I(Ca,L) stabilizes SWs in a two-dimensional epicardial layer of rabbit ventricular myocardium to help their persistence, whereas blockade of both currents destabilizes SWs to facilitate their termination.

IK1 heterogeneity affects genesis and stability of spiral waves in cardiac myocyte monolayers

Previous studies have postulated an important role for the inwardly rectifying potassium current (I(K1)) in controlling the dynamics of electrophysiological spiral waves responsible for ventricular tachycardia and fibrillation. In this study, we developed a novel tissue model of cultured neonatal rat ventricular myocytes (NRVMs) with uniform or heterogeneous Kir2.1expression achieved by lentiviral transfer to elucidate the role of I(K1) in cardiac arrhythmogenesis. Kir2.1-overexpressed NRVMs showed increased I(K1) density, hyperpolarized resting membrane potential, and increased action potential upstroke velocity compared with green fluorescent protein-transduced NRVMs. Opposite results were observed in Kir2.1-suppressed NRVMs. Optical mapping of uniformly Kir2.1 gene-modified monolayers showed altered conduction velocity and action potential duration compared with nontransduced and empty vector-transduced monolayers, but functional reentrant waves could not be induced. In monolayers with an island of altered Kir2.1 expression, conduction velocity and action potential duration of the locally transduced and nontransduced regions were similar to those of the uniformly transduced and nontransduced monolayers, respectively, and functional reentrant waves could be induced. The waves were anchored to islands of Kir2.1 overexpression and remained stable but dropped in frequency and meandered away from islands of Kir2.1 suppression. In monolayers with an inverse pattern of I(K1) heterogeneity, stable high frequency spiral waves were present with I(K1) overexpression, whereas lower frequency, meandering spiral waves were observed with I(K1) suppression. Our study provides direct evidence for the contribution of I(K1) heterogeneity and level to the genesis and stability of spiral waves and highlights the potential importance of I(K1) as an antiarrhythmia target.

Eliminating spiral waves pinned to an anatomical obstacle in cardiac myocytes by high-frequency stimuli

The unpinning of spiral waves by the application of high-frequency wave trains was studied in cultured cardiac myocytes. Successful unpinning was observed when the frequency of the paced waves exceeded a critical level. The unpinning process was analyzed by a numerical simulation with a model of cardiac tissue. The mechanism of unpinning by high-frequency stimuli is discussed in terms of local entrainment failure, through a reduction of the two-dimensional spatial characteristics into one dimension.

Role of potassium currents in cardiac arrhythmias

Abnormal excitability of myocardial cells may give rise to ectopic beats and initiate re-entry around an anatomical or functional obstacle. As K(+) currents control the repolarization process of the cardiac action potential (AP), the K(+) channel function determines membrane potential and refractoriness of the myocardium. Both gain and loss of the K(+) channel function can lead to arrhythmia. The former because abbreviation of the active potential duration (APD) shortens refractoriness and wave length, and thereby facilitates re-entry and the latter because excessive prolongation of APD may lead to torsades de pointes (TdP) arrhythmia and sudden cardiac death. The pro-arrhythmic consequences of malfunctioning K(+) channels in ventricular and atrial tissue are discussed in the light of three pathophysiologically relevant aspects: genetic background, drug action, and disease-induced remodelling. In the ventricles, loss-of-function mutations in the genes encoding for K(+) channels and many drugs (mainly hERG channel blockers) are related to hereditary and acquired long-QT syndrome, respectively, that put individuals at high risk for developing TdP arrhythmias and life-threatening ventricular fibrillation. Similarly, down-regulation of K(+) channels in heart failure also increases the risk for sudden cardiac death. Mutations and polymorphisms in genes encoding for atrial K(+) channels can be associated with gain-of-function and shortened, or with loss-of-function and prolonged APs. The block of atrial K(+) channels becomes a particular therapeutic challenge when trying to ameliorate atrial fibrillation (AF). This arrhythmia has a strong tendency to cause electrical remodelling, which affects many K(+) channels. Atrial-selective drugs for the treatment of AF without affecting the ventricles could target structures such as I(Kur) or constitutively active I(K,ACh) channels.