Incidence of Dispersion of Refractoriness and Cellular Coupling Resistance on Cardiac Reentries and Ventricular Fibrillation

We used computer simulations to study the possible role of the dispersion of cellular coupling, refractoriness or both, in the mechanisms underlying cardiac arrhythmias. Local ischemia was first assumed to induce cell to cell dispersion of the coupling resistance (Case 1), refractory period (Case 2), or both of them (Case 3). Our numerical experiments based on the van Capelle and Durrer model showed that vortices could not be induced by cell to cell variations. With cellular properties dispersed in a patchy way within the ischemic zone, a single activation wave could give rise to abnormal activities. This demonstrates the stability of the wave front under small inhomogeneities. Probabilities of reentry, estimated for the three cases cited above showed that a severe alteration of the coupling resistance may be an important factor in the genesis of reentry. Moreover, use of isochronal maps revealed that vortices were both stable and sustained with an alteration of the coupling alone or combined with a reduction of the action potential duration. Conversely, simulations with reduction of the refractoriness alone, inducing only transient patterns, could exhibit functionally determined reentries.     

Role of the dispersion of refractoriness on cardiac reentries

We used computer simulation to study the possible role of the dispersion of cellular coupling, refractoriness or both, in the mechanisms underlying cardiac arrhythmias. Local ischemia was first assumed to induce cell to cell dispersion of the coupling resistance (case 1), refractory period (case 2), or both (case 3). Our numerical experiments based on the van Capelle and Durrer model showed that vortices could not be induced. On the other hand, with cellular properties dispersed in a patchy way within the ischemic zone, a single activation wave could give rise to abnormal activities. This demonstrates the stability of the wave front under small inhomogeneities. Probabilities of reentry, estimated for the three cases cited above showed that a severe alteration of the coupling resistance may be an important factor in the genesis of reentry. Moreover, use of isochronal maps revealed that vortices were both stable and sustained with an alteration of the coupling alone or along with a reduction of the action potential duration. Conversely, simulations with reduction of the refractoriness alone, inducing only transient patterns, could exhibit functionally determined reentries.     

Interpretation of epicardial mapping by means of computer simulations: Applications to calcium,lidocaine and to BRL 34915

The aim of this work was to compare experimental investigations on effects of lidocaine, calcium and, BRL 34915 on reentries to simulated data obtained by use of a model of propagation based on the Huygens' constriction method already described in previous works. Calcium and lidocaine effects are investigated on anisotropic conduction conditions. In both cases, reduction in conduction velocities are observed. In lidocaine case, a refractory area is located along the longitudinal axis. In agreement with experimental electrical mapping, the simulations show that the stabilization of reentrant excitation is mainly due to the existence of this refractory area around which the reentrant circuit can develop. The experimental study shows that BRL 34915 has both arrhythmogenic and antiarrhythmic effects. A detailed electrophysiological analysis has shown that drug infusion act on normal cardiac cells by decreasing the relative and absolute refractory period. BRL 34915 action is simulated by a decrease in the refractory period showing that the time frequency of the reentrant activity is increased and that the spatial size where the reentry is developing is becoming smaller. These two effects are arrhythmogenic, the simulated data being so in good agreement with the experimental ones.     

Effects of simulated ischemia on spiral wave stability

Regional hyperkalemia during acute myocardial ischemia is a major factor promoting electrophysiological abnormalities leading to ventricular fibrillation (VF). However, steep action potential duration restitution, recently proposed to be a major determinant of VF, is typically decreased rather than increased by hyperkalemia and acute ischemia. To investigate this apparent contradiction, we simulated the effects of regional hyperkalemia and other ischemic components (anoxia and acidosis) on the stability of spiral wave reentry in simulated two-dimensional cardiac tissue by use of the Luo-Rudy ventricular action potential model. We found that the hyperkalemic "ischemic" area promotes wavebreak in the surrounding normal tissue by accelerating the rate of spiral wave reentry, even after the depolarized ischemic area itself has become unexcitable. Furthermore, wavebreak and fibrillation can be prevented if the dynamical instability of the normal tissue is reduced significantly by targeting electrical restitution properties, suggesting a novel therapeutic approach.     

Simulated Interactions between a Class III Antiarrhythmic Drug and a Figure 8 Reentry

Ventricular Fibrillation is responsible for a majority of sudden cardiac death, but little is known about how ventricular tachycardia (VT) degenerates into ventricular fibrillation. Several clinical studies focused only on preventing VT with a class III antiarrhythmic drug resulted in many deaths. Our simulations investigate the interactions between an antiarrhythmic drug likely to suppress a VT and a Figure 8 reentry. A parameter AAR is introduced to increase the action potential duration and therefore simulate various Class III drugs. Simulations are ran under several conditions (phases of the reentry, values of AAR, durations). They show that a VT can be suppressed whatever the phase of the reentry but it strongly depends on the duration of the effect. It confirms that a drug which can suppress a reentry can also worsen it. It also shows a great variety of activation patterns and thus the complexity of antiarrhythmic drugs effects. Simulations also demonstrate that suppressing VT is an increasing function of AAR.     

Theoretical Study of Cardiac Transient Conduction Blocks on Reentries Induction. Applications to Antiarrhythmic Drugs

Limitations of antiarrhythmic drugs on cardiac sudden death prevention appeared since the early 80's. The "Cardiac Arrhythmia Suppression Trial"(CAST) showed more recently that mortality was significantly higher inpatients treated with some particular antiarrhythmic drugs than in non-treated patients. In this field, our group recently demonstrated that a bolus of a Class 1B antiarrhythmic drug was able to trigger a ventricular fibrillation due to transient blocks induction. The aim of the present work was to systematically study, by use of the van Capelle and Durrer (VCD) model which allows to simulate ventricular activation wave propagation, the link between arrhythmogenic effects and the ability of transient blocks to possibly degenerate in severe arrhythmias. A fragment of the ventricular wall is represented by an array of 16384elements electrically coupled. Effects of induction of one or several transient blocks, as the effects of their size and duration on possible induction of reentries have been studied. Results obtained show that various combinations between these different parameters may trigger reentries, ventricular tachycardia and/or more complex patterns assimilable to ventricular fibrillation. These results clearly evidence the fact that possible induction of transient blocks may directly be related to risk factor associated to arrhythmogenic effects of antiarrhythmic drugs. This revised version was published online in June 2006 with corrections to the Cover Date.     

Mechanisms of unpinning and termination of ventricular tachycardia

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

Vulnerable window for conduction block in a one-dimensional cable of cardiac cells, 2: multiple extrasystoles

Unidirectional conduction block of premature extrasystoles can lead to initiation of cardiac reentry, causing lethal arrhythmias including ventricular fibrillation. Multiple extrasystoles are often more effective at inducing unidirectional conduction block and reentry than a single extrasystole. Since the substrate for conduction block is spatial dispersion of refractoriness, in this study we investigate how the first extrasystole modulates this dispersion to influence the "vulnerable window" for conduction block by subsequent extrasystoles, particularly in relation to action potential duration restitution and conduction velocity restitution properties. Using a kinematic model to represent wavefront-waveback interactions and simulations with the Luo-Rudy model in a one-dimensional cable of cardiac cells, we show that in homogeneous tissue, a premature extrasystole can create a large dispersion of refractoriness leading to conduction block of a subsequent extrasystole. In heterogeneous tissue, however, a premature extrasystole can either reduce or enhance the dispersion of refractoriness depending on its propagation direction with respect to the previous beat. With multiple extrasystoles at random coupling intervals, vulnerability to conduction block is proportional to their number. In general, steep action potential duration restitution and broad conduction velocity restitution promote dispersion of refractoriness in response to multiple extrasystoles, and thus enhance vulnerability to conduction block. These restitution properties also promote spatially discordant alternans, a setting which is particularly prone to conduction block. The equivalent dispersion of refractoriness created dynamically in homogeneous tissue by spatially discordant alternans is more likely to cause conduction block than a comparable degree of preexisting dispersion in heterogeneous tissue.     

A physical approach to remove anatomical reentries: a bidomain study

Controlling cardiac chaos is often achieved by applying a large damaging electric shock-defibrillation. It removes all waves, without differentiating reentries and normal waves, anatomical and functional reentries. Anatomical reentries can be removed by anti-tachycardia pacing (ATP) as well. But ATP requires the knowledge of the position of the reentry, and an access to it with an invasive stimulating electrode. We show that the physics of electric field distribution between cardiac cells permits one to deliver an electric pulse exactly to the core of an anatomical reentry, without knowing its position and even to locations where access with a stimulating electrode is not possible. The energy needed is two orders of magnitude less than defibrillation energy. The results are insensitive to both a detailed ionic model and to the geometry of the fibers.     

全心缺血早期阶段室性心律失常仿真研究

为了分析全心缺血早期阶段对心脏电生理活动的影响,以及探讨诱发的室性心律失常机制,本研究考虑了缺血情况下高钾、酸液过多、局部缺氧的情况,结合详细的人类心室细胞生物物理上的动力学特征,开发了一个人体心室细胞和组织全心缺血模型.实验结果表明,全心缺血缩短了动作电位时程(action potential duration,APD),且减缓了兴奋的传导速率(conduction velocity,CV).同时,由于全心缺血降低了动作电位时程曲线(action potential duration restitution,APDR)斜率,且增大了有效不应期(effective refractory period,ERP),因此有利于维持折返波的稳定传导,使得室速不易转化为室颤.另一方面,尽管全心缺血导致了组织易感性的增加,但是由于其需要更长的异位刺激长度来保证折返波的形成,因此也在一定程度上降低了心律失常的发生概率.     

Effects of Na(+) channel and cell coupling abnormalities on vulnerability to reentry: a simulation study

The role of dynamic instabilities in the initiation of reentry in diseased (remodeled) hearts remains poorly explored. Using computer simulations, we studied the effects of altered Na(+) channel and cell coupling properties on the vulnerable window (VW) for reentry in simulated two-dimensional cardiac tissue with and without dynamic instabilities. We related the VW for reentry to effects on conduction velocity, action potential duration (APD), effective refractory period dispersion and restitution, and concordant and discordant APD alternans. We found the following: 1). reduced Na(+) current density and slowed recovery promoted postrepolarization refractoriness and enhanced concordant and discordant APD alternans, which increased the VW for reentry; 2). uniformly reduced cell coupling had little effect on cellular electrophysiological properties and the VW for reentry. However, randomly reduced cell coupling combined with decoupling promoted APD dispersion and alternans, which also increased the VW for reentry; 3). the combination of decreased Na(+) channel conductance, slowed Na(+) channel recovery, and cellular uncoupling synergistically increased the VW for reentry; and 4) the VW for reentry was greater when APD restitution slope was steep than when it was flat. In summary, altered Na(+) channel and cellular coupling properties increase vulnerability to reentrant arrhythmias. In remodeled hearts with altered Na(+) channel properties and cellular uncoupling, dynamic instabilities arising from electrical restitution exert important influences on the VW for reentry.     

Is there a correlation between ventricular fibrillation cycle length and electrophysiological and anatomic properties of the canine left ventricle?

We hypothesized that myocardial infarction-related alterations in ventricular fibrillation (VF) cycle length (VFCL) would correlate with changes in local cardiac electrophysiological and anatomic properties. An electrophysiological study was performed in normal, subacute, and chronic infarction mongrel dogs. VF was induced by programmed electrical stimulation and mean and minimum early and late VFCL was determined and correlated with local electrophysiological and anatomic properties. Effective refractory period (ERP), activation recovery time (ART), ERP/ART ratio, threshold, and ERP and ART dispersion were determined at 112 sites on the anterior left ventricle. Wave front progression was analyzed over a 2-s period. The extent of local tissue necrosis and of myocardial fiber disarray was also evaluated. The early mean VFCL was significantly longer in the subacute infarction (149 +/- 35 ms) and chronic infarction dogs (129 +/- 18 ms) compared with control dogs (102 +/- 15 ms; P < 0.0001 for both comparisons) as was the early minimum VFCL with similar trends seen during late VF. Complete epicardial reentrant circuits were significantly more common in normal dogs (4.3 +/- 2.4, 22.4% of cycles) than in subacute (0.75 +/- 0.96, 5.3% of cycles, P < 0.05 vs. normal) and chronic infarction dogs (1.3 +/- 1.3, 7.5% of cycles, P < 0.05 vs. normal). There was a poor correlation between the mean and minimum early and late VFCL and local electrophysiological and anatomic properties (R(2) < 0.2 for all comparisons) with a much better correlation between average mean and minimum VFCL (over the entire plaque) and global ERP and ART dispersion during early and late VF. In conclusion, VFCL in normal and infarcted myocardium shows a poor correlation with local ventricular electrophysiological and anatomic properties measured in sinus rhythm. However, there was a much better correlation between the average VFCL with global dispersion of repolarization. The lack of correlation between local VFCL and refractoriness and the infrequent occurrence of epicardial reentry suggests that intramural reentry may be the primary mechanism of VF in this model.     

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

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

Electrical refractory period restitution and spiral wave reentry in simulated cardiac tissue

Theoretical and experimental studies have shown that restitution of the cardiac action potential (AP) duration (APD) plays a major role in predisposing ventricular tachycardia to degenerate to ventricular fibrillation, whereas its role in atrial fibrillation is unclear. We used the Courtemanche human atrial cell model and the Luo-Rudy guinea pig ventricular model to compare the roles of electrical restitution in destabilizing spiral wave reentry in simulated two-dimensional homogeneous atrial and ventricular tissue. Because atrial AP morphology is complex, we also validated the usefulness of effective refractory period (ERP) restitution. ERP restitution correlated best with APD restitution at transmembrane potentials greater than or equal to -62 mV, and its steepness was a reliable predictor of spiral wave phenotype (stable, meandering, hypermeandering, and breakup) in both atrial and ventricular tissue. Spiral breakup or single hypermeandering spirals occurred when the slope of ERP restitution exceeded 1 at short diastolic intervals. Thus ERP restitution, which is easier to measure clinically than APD restitution, is a reliable determinant of spiral wave stability in simulated atrial and ventricular tissue.     

Dynamics of intramural scroll waves in three-dimensional continuous myocardium with rotational anisotropy.

It has been suggested that reentrant activity in three-dimensional cardiac muscle may be organized as a scroll wave rotating around a singularity line called the filament. Experimental studies indicate that filaments are often concealed inside the ventricular wall and consequently, scroll waves do not manifest reentrant activity on the surface. Here we analyse how such concealed scroll waves are affected by a twisted anisotropy resulting from rotation of layers of muscle fibers inside the ventricular wall. We used a computer model of a ventricular slab (15x15x15 mm(3)) with a fiber twist of 120 degrees from endocardium to epicardium. The action potential was simulated using FitzHugh-Nagumo equations. Scroll waves with rectilinear filaments were initiated at various depths of the slab and at different angles with respect to fiber orientation. The analysis shows that independent of initial conditions, after a certain transitional period, the filament aligns with the local fiber orientation. The alignment of the filament is determined by the directional variations in cell coupling due to fiber rotation and by boundary conditions. Our findings provide a mechanistic explanation for the prevalence of intramural reentry over transmural reentry during polymorphic ventricular tachycardia and fibrillation.     

Short-term cardiac memory and mother rotor fibrillation

Short-term cardiac memory refers to the effects of pacing history on action potential duration (APD). Although the ionic mechanisms for short-term memory occurring over many heartbeats (also called APD accommodation) are poorly understood, they may have important effects on reentry and fibrillation. To explore this issue, we incorporated a generic memory current into the Phase I Luo and Rudy action potential model, which lacks short-term memory. The properties of this current were matched to simulate quantitatively human ventricular monophasic action potential accommodation. We show that, theoretically, short-term memory can resolve the paradox of how mother rotor fibrillation is initiated in heterogeneous tissue by physiological pacing. In simulated heterogeneous two-dimensional tissue and three-dimensional ventricles containing an inward rectifier K(+) current gradient, short-term memory could spontaneously convert multiple wavelet fibrillation to mother rotor fibrillation or to a mixture of both fibrillation types. This was due to progressive acceleration and stabilization of rotors as accumulation of memory shortened APD and flattened APD restitution slope nonuniformly throughout the tissue.     

Computer simulation of the reentrant cardiac arrhythmias in ischemic myocardium

Computer simulation was performed to determine how reentrant activity could occur due to the spatial heterogeneity in refractoriness induced by the regional ischemia. Two regional ischemic models were developed by decreasing the intracellular ATP concentration, reducing conductance of the inward Na+ current and increasing the extracellular K+ concentration on the two-dimensional sheet. Operator splitting method was used to integrate the models. The vulnerability to reentry was estimated from the timings of premature stimuli on the constructed models, which could result in unidirectionally propagating action potentials. Two kinds of sustained spiral waves and their Pseudo-Electroscardiograms were observed in numerical simulation. The results showed that the dispersion of refractory period increased with ischemic aggravation, and led to augment of the vulnerable window. A permature stimulation within the vulnerable window could easily induce spiral reentry. The Pseudo-Electrocardiograms of the spiral waves exhibited monomorphic tachycardiac waveforms. Thus, the spatial heterogeneity in refractoriness could be a substrate for reentrant ventricular tachyarrhythmias on the regional ischemic tissue.     

Myocardial fibrosis blunts nitric oxide synthase-related preload reserve in human dilated cardiomyopathy

It is clear that ischemia inhibits successful defibrillation by altering regional electro-physiology. However, the exact mechanisms are unclear. This study investigated whether regional gap junction inhibition increases biphasic shock defibrillation thresholds (DFT). Sixteen swine were instrumented with a mid-left anterior descending (LAD) perfusion catheter for regional infusion of 0.5 mM/h heptanol (n = 8) or saline (n = 8). DFT values and effective refractory periods (ERP) at five myocardial sites were determined. Regional conduction velocity (CV) was determined in an LAD drug-perfused and nondrug-perfused region in an additional seven swine. Regional heptanol infusion increased 50% DFT values by 33% (P = 0.01) and slowed CV by 42-59% (P < 0.01) but did not affect ERP. Regional heptanol also increased CV dispersion by approximately 270% (P < 0.05) but did not change ERP dispersion. Regional placebo did not alter any of these parameters. Furthermore, regional heptanol infusion induced spontaneous ventricular fibrillation in eight of eight animals. Increasing spatial conduction velocity dispersion by impairing regional gap junction conductance increased DFT values. Dispersion in conduction velocity slowing during regional ischemia may be an important determinant of defibrillation efficacy.     

基于虚拟心脏的早期后除极导致室颤的仿真研究

为了分析早期后除极(early afterdepolarizations,EADs)诱发室颤的机理,本研究基于精细的浦肯野纤维网络与心室解剖数据,构建了一个三维心室电传导模型.基于该模型,模拟了产生早期后除极的电生理变化,探讨了三种心室细胞的早期后除极的易感性,分析了早期后除极易感细胞对折返波的影响,最后定量比较早期后除极诱发室颤的伪心电图的改变情况.实验结果表明:中间层细胞早期后除极易感性最强,中间层细胞早期后除极的产生能够导致折返波破裂,并且在心电图中表现为紊乱的不规则的颤动心律,这与之前在动物实验观察得到的现象一致,因此中间层细胞可能是一个诱发室颤的重要靶点.     

Oscillation in cycle length induces transient discordant and steady-state concordant alternans in the heart

Alternans is a beat-to-beat alternation of the cardiac action potential duration (APD) or intracellular calcium (Ca(i)) transient. In cardiac tissue, alternans can be spatially concordant or discordant, of which the latter has been shown to increase dispersion of repolarization and promote a substrate for initiation of ventricular fibrillation. Alternans has been studied almost exclusively under constant cycle length pacing conditions. However, heart rate varies greatly on a beat-by-beat basis in normal and pathological conditions. The purpose of this study was to determine if applying a repetitive but non-constant pacing pattern, specifically cycle length oscillation (CLO), promotes or suppresses a proarrhythmic substrate. We performed computational simulations and optical mapping experiments to investigate the potential consequences of CLO. In a single cell computational model, CLO induced APD and Ca(i) alternans, which became "phase-matched" with the applied oscillation. As a consequence of the phase-matching, in one-dimensional cable simulations, neonatal rat ventricular myocyte monolayers, and isolated adult guinea pig hearts CLO could transiently induce spatial and electromechanical discordant alternans followed by a steady-state of concordance. Our results demonstrated that under certain conditions, CLO can initiate ventricular fibrillation in the isolated hearts. On the other hand, CLO can also exert an antiarrhythmic effect by converting an existing state of discordant alternans to concordant alternans.