Low Energy Defibrillation in Human Cardiac Tissue: A Simulation Study

We aim to assess the effectiveness of feedback-controlled resonant drift pacing as a method for low energy defibrillation. Antitachycardia pacing is the only low energy defibrillation approach to have gained clinical significance, but it is still suboptimal. Low energy defibrillation would avoid adverse side effects associated with high voltage shocks and allow the application of implantable cardioverter defibrillator (ICD) therapy, in cases where such therapy is not tolerated today. We present results of computer simulations of a bidomain model of cardiac tissue with human atrial ionic kinetics. Reentry was initiated and low energy shocks were applied with the same period as the reentry, using feedback to maintain resonance. We demonstrate that such stimulation can move the core of reentrant patterns, in the direction that depends on the location of the electrodes and the time delay in the feedback. Termination of reentry is achieved with shock strength one-order-of-magnitude weaker than in conventional single-shock defibrillation. We conclude that resonant drift pacing can terminate reentry at a fraction of the shock strength currently used for defibrillation and can potentially work where antitachycardia pacing fails, due to the feedback mechanisms. Success depends on a number of details that these numerical simulations have uncovered.     

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.     

Differences between left and right ventricular anatomy determine the types of reentrant circuits induced by an external electric shock. A rabbit heart simulation study

Despite the fact that elucidating the mechanisms of cardiac vulnerability to electric shocks is crucial to understanding why defibrillation shocks fail, important aspects of cardiac vulnerability remain unknown. This research utilizes a novel anatomically based bidomain finite-element model of the rabbit ventricles to investigate the effect of shock polarity reversal on the reentrant activity induced by an external defibrillation-strength shock in the paced ventricles. The specific goal of the study is to examine how differences between left and right ventricular chamber anatomy result in differences in the types of reentrant circuits established by the shock. Truncated exponential monophasic shocks of duration 8 ms were delivered via two external electrodes at various timings. Vulnerability grids were constructed for shocks of reversed polarity (referred to as RV- or LV- when either the RV or the LV electrode is a cathode). Our results demonstrate that reversing electrode polarity from RV- to LV- changes the dominant type of post-shock reentry: it is figure-of-eight for RV- and quatrefoil for LV- shocks. Differences in secondary types of post-shock arrhythmia also occur following shock polarity reversal. These effects of polarity reversal are primarily due to the fact that the LV wall is thicker than the RV, resulting in a post-shock excitable gap that is predominantly within the LV wall for RV- shocks and in the septum for LV- shocks.     

Reentrant waves in a ring of embryonic chick ventricular cells imaged with a Ca2+ sensitive dye

According to the classic model initially formulated by Mines, reentrant cardiac arrhythmias may be associated with waves circulating in a ring geometry. This study was designed to study the dynamics of reentry in a ring geometry of cardiac tissue culture. Reentrant calcium waves in rings of cultured embryonic chick cardiac myocytes were imaged using a macroscope to monitor the fluorescence of intracellular Calcium Green-1 dye. The rings displayed a variety of stable rhythms including pacemaker activity and spontaneous reentry. Waves originating from a localized pacemaker could lead to reentry as a consequence of unidirectional block. In addition, more complex patterns were observed due to the interactions between reentrant and pacemaker rhythms. These rhythms included instances in which pacemakers accelerated the reentrant rhythm, and instances in which the excitation was blocked in the vicinity of pacemakers. During reentrant activity an appropriately timed electrical stimulus could induce resetting of activity or cause complete annihilation of the propagating waves. This experimental preparation reveals many spontaneously occuring complex rhythms. These complex rhythms are hypothesized to reflect interactions between spontaneous pacemakers, wave propagation, refractory period, and overdrive suppression. This preparation may serve as a useful model system to further investigate complex dynamics arising during reentrant rhythms in cardiac 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.     

The effect of spatial scale of resistive inhomogeneity on defibrillation of cardiac tissue

Defibrillation of cardiac tissue can be viewed in the context of dynamical systems theory as the attempt to move a dynamical system from the basin of attraction of one attractor (fibrillation) to another (the uniform rest state) by applying a stimulus whose form is physically constrained. Here we give an introduction to the physical mechanism of cardiac defibrillation from this dynamical perspective and examine the role of resistive inhomogeneity on defibrillation efficacy. Using numerical simulations with rotating waves on a one-dimensional periodic ring, we study the role of the spatial scale of resistive inhomogeneity on defibrillation. For a rotating wave on a periodic ring there are three stable attractors, namely the uniform rest state, a wave traveling clockwise and a wave traveling counterclockwise. As a result, the application of a stimulus has the potential for three different outcomes, namely elimination of the wave, phase resetting of the wave, and reversal of the wave. The results presented here show that with resistive inhomogeneities of large spatial scale, all three of these transitions are possible with large amplitude shocks, so that the probability of defibrillation is bounded well below one, independent of stimulus amplitude. On the other hand, resistive inhomogeneities of small spatial scale produce a defibrillation threshold that is qualitatively consistent with that found experimentally, namely the probability of defibrillation success is an increasing function that approaches one for large enough stimulus amplitude. Extending these results to higher dimensions, we describe conditions for successful defibrillation of functional reentry with large scale spatial inhomogeneity, but find that elimination of anatomical reentry is quite difficult. With small spatial scale inhomogeneity, there are no similar restrictions.     

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.     

Heterogeneous gap junction remodeling stabilizes reentrant circuits in the epicardial border zone of the healing canine infarct: a computational study

The ventricular tachycardias (VTs) that originate in the 5-day epicardial border zone (EBZ) of the healing canine infarcted heart are due to reentrant excitation. In cells surviving in the EBZ, both sarcolemmal ionic channels and gap junction conductance and distribution are remodeled. We previously showed that the heterogeneities in sodium current (I(Na)) and L-type calcium channel current (I(CaL)) of the center and outer pathway cells result in a homogenization of the refractory period that in turn stabilizes reentrant VTs for approximately 10 beats. To understand how heterogeneities in transverse gap junctional conductance remodeling reported experimentally contribute to the stability of these tachycardias, we studied the dynamics of reentering waves in two-dimensional computer models of the EBZ. First we used a computer model with homogeneous ionic channel properties [infarcted border zone cell model (IZ)]. These simulations show that, in the absence of heterogeneities in ionic channel properties, reentrant waves tend to drift to localized regions of uncoupling and stabilize there. Second, we used a computer model with a more realistic representation of the heterogeneous EBZ, including cellular models for both the center (IZ(c)) and outer (IZ(o)) pathway cells. These simulations show that neither a region of uniform uncoupling nor a step transition between two regions with different side-to-side (transverse) cell coupling stabilizes reentry in this substrate. However, an area of localized uncoupling did stabilize reentry in such a model. We propose that in addition to the heterogeneities in I(Na) and I(CaL) properties, heterogeneities in gap junctional conductance in the EBZ causing regions of localized uncoupling stabilize VT in the EBZ. Previous experimental in situ activation maps of the 5-day EBZ show that the lines of block form in regions of slow transverse propagation. This is consistent with our findings that areas of localized uncoupling stabilize reentry.     

Functional reentry and circus movement arrhythmias in the small intestine of normal and diabetic rats

In a few recent studies, the presence of arrhythmias based on reentry and circus movement of the slow wave have been shown to occur in normal and diseased stomachs. To date, however, reentry has not been demonstrated before in any other part of the gastrointestinal system. No animals had to be killed for this study. Use was made of materials obtained during the course of another study in which 11 rats were treated with streptozotocin and housed with age-matched controls. After 3 and 7 mo, segments of duodenum, jejunum, and ileum were isolated and positioned in a tissue bath. Slow wave propagation was recorded with 121 extracellular electrodes. After the experiment, the propagation of the slow waves was reconstructed. In 10 of a total of 66 intestinal segments (15%), a circus movement of the slow wave was detected. These reentries were seen in control (n = 2) as well as in 3-mo (n = 2) and 7-mo (n = 6) diabetic rats. Local conduction velocities and beat-to-beat intervals during the reentries were measured (0.42 ± 0.15 and 3.03 ± 0.67 cm/s, respectively) leading to a wavelength of 1.3 ± 0.5 cm and a circuit diameter of 4.1 ± 1.5 mm. This is the first demonstration of a reentrant arrhythmia in the small intestine of control and diabetic rats. Calculations of the size of the circuits indicate that they are small enough to fit inside the intestinal wall. Extrapolation based on measured velocities and rates indicate that reentrant arrhythmias are also possible in the distal small intestine of larger animals including humans.     

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.     

Electrophysiological and Structural Remodeling in Heart Failure Modulate Arrhythmogenesis. 2D Simulation Study

BackgroundHeart failure is operationally defined as the inability of the heart to maintain blood flow to meet the needs of the body and it is the final common pathway of various cardiac pathologies. Electrophysiological remodeling, intercellular uncoupling and a pro-fibrotic response have been identified as major arrhythmogenic factors in heart failure.ObjectiveIn this study we investigate vulnerability to reentry under heart failure conditions by incorporating established electrophysiological and anatomical remodeling using computer simulations.MethodsThe electrical activity of human transmural ventricular tissue (5 cm×5 cm) was simulated using the human ventricular action potential model Grandi et al. under control and heart failure conditions. The MacCannell et al. model was used to model fibroblast electrical activity, and their electrotonic interactions with myocytes. Selected degrees of diffuse fibrosis and variations in intercellular coupling were considered and the vulnerable window (VW) for reentry was evaluated following cross-field stimulation.ResultsNo reentry was observed in normal conditions or in the presence of HF ionic remodeling. However, defined amount of fibrosis and/or cellular uncoupling were sufficient to elicit reentrant activity. Under conditions where reentry was generated, HF electrophysiological remodeling did not alter the width of the VW. However, intermediate fibrosis and cellular uncoupling significantly widened the VW. In addition, biphasic behavior was observed, as very high fibrotic content or very low tissue conductivity hampered the development of reentry. Detailed phase analysis of reentry dynamics revealed an increase of phase singularities with progressive fibrotic components.ConclusionStructural remodeling is a key factor in the genesis of vulnerability to reentry. A range of intermediate levels of fibrosis and intercellular uncoupling can combine to favor reentrant activity.     

Region [corrected] of slowed conduction acts as core for spiral wave reentry in cardiac cell monolayers

Pathophysiological heterogeneity in cardiac tissue is related to the occurrence of arrhythmias. Of importance are regions of slowed conduction, which have been implicated in the formation of conduction block and reentry. Experimentally, it has been a challenge to produce local heterogeneity in a manner that is both reversible and well controlled. Consequently, we developed a dual-zone superfusion chamber that can dynamically create a small (5 mm) central island of heterogeneity in cultured cardiac cell monolayers. Three different conditions were studied to explore the effect of regionally slowed conduction on wave propagation and reentry: depolarization by elevated extracellular potassium, sodium channel inhibition with lidocaine, and cell-cell decoupling with palmitoleic acid. Using optical mapping of transmembrane voltage, we found that the central region of slowed conduction always served as the core region around which a spiral wave formed and then revolved following a period of rapid pacing. Because of the localized slowing in the core region, we observed experimentally for the first time an S shape of the spiral wave front near its tip. These results indicate that a small region of slowed conduction can play a crucial role in the formation, anchoring, and modulation of reentrant spiral waves.     

Intracellular calcium dynamics at the core of endocardial stationary spiral waves in Langendorff-perfused rabbit hearts

In vitro models of sustained monomorphic ventricular tachycardia (MVT) are rare and do not usually show spiral reentry on the epicardium. We hypothesized that MVT is associated with the spiral wave in the endocardium and that this stable reentrant propagation is supported by a persistently elevated intracellular calcium (Ca(i)) transient at the core of the spiral wave. We performed dual optical mapping of transmembrane potential (V(m)) and Ca(i) dynamics of the right ventricular (RV) endocardium in Langendorff-perfused rabbit hearts (n = 12). Among 64 induced arrhythmias, 55% were sustained MVT (>10 min). Eighty percent of MVT showed stationary spiral waves (>10 cycles, cycle length: 128 +/- 14.6 ms) in the endocardial mapped region, anchoring to the anatomic discontinuities. No reentry activity was observed in the epicardium. During reentry, the amplitudes of V(m) and Ca(i) signals were higher in the periphery and gradually decreased toward the core. At the core, maximal V(m) and Ca(i) amplitudes were 42.95 +/- 5.89% and 43.95 +/- 9.46%, respectively, of the control (P < 0.001). However, the trough of the V(m) and Ca(i) signals at the core were higher than those in the periphery, indicating persistent V(m) and Ca(i) elevations during reentry. BAPTA-AM, a calcium chelator, significantly reduced the maximal Ca(i) transient amplitude and prevented sustained MVT and spiral wave formation in the mapped region. These findings indicate that endocardial spiral waves often anchor to anatomic discontinuities causing stable MVT in normal rabbit ventricles. The spiral core is characterized by diminished V(m) and Ca(i) amplitudes and persistent V(m) and Ca(i) elevations during reentry.     

Mechanisms of superiority of ascending ramp waveforms: new insights into mechanisms of shock-induced vulnerability and defibrillation

Monophasic ascending ramp (AR) and descending ramp (DR) waveforms are known to have significantly different defibrillation thresholds. We hypothesized that this difference arises due to differences in mechanisms of arrhythmia induction for the two waveforms. Rabbit hearts (n = 10) were Langendorff perfused, and AR and DR waveforms (7, 20, and 40 ms) were randomly delivered from two line electrodes placed 10 mm apart on the anterior ventricular epicardium. We optically mapped cellular responses to shocks of various strengths (5, 10, and 20 V/cm) and coupling intervals (CIs; 120, 180, and 300 ms). Optical mapping revealed that maximum virtual electrode polarization (VEP) was reached at significantly different times for AR and DR of the same duration (P < 0.05) for all tested CIs. As a result, VEP for AR were stronger than for DR at the end of the shock. Postshock break excitation resulting from AR generated faster propagation and typically could not form reentry. In contrast, partially dissipated VEP resulting from DR generated slower propagation; the wavefront was able to propagate into deexcited tissue and thus formed a shock-induced reentry circuit. Therefore, for the same delivered energy, AR was less proarrhythmic compared with DR. An active bidomain model was used to confirm the electrophysiological results. The VEP hypothesis explains differences in vulnerability associated with monophasic AR and DR waveforms and, by extension, the superior defibrillation efficacy of the AR waveform compared with the DR waveform.     

Theoretical analysis of the magnetocardiographic pattern for reentry wave propagation in a three-dimensional human heart model

We present a computational study of reentry wave propagation using electrophysiological models of human cardiac cells and the associated magnetic field map of a human heart. We examined the details of magnetic field variation and related physiological parameters for reentry waves in two-dimensional (2-D) human atrial tissue and a three-dimensional (3-D) human ventricle model. A 3-D mesh system representing the human ventricle was reconstructed from the surface geometry of a human heart. We used existing human cardiac cell models to simulate action potential (AP) propagation in atrial tissue and 3-D ventricular geometry, and a finite element method and the Galerkin approximation to discretize the 3-D domain spatially. The reentry wave was generated using an S1-S2 protocol. The calculations of the magnetic field pattern assumed a horizontally layered conductor for reentry wave propagation in the 3-D ventricle. We also compared the AP and magnetocardiograph (MCG) magnitudes during reentry wave propagation to those during normal wave propagation. The temporal changes in the reentry wave motion and magnetic field map patterns were also analyzed using two well-known MCG parameters: the current dipole direction and strength. The current vector in a reentry wave forms a rotating spiral. We delineated the magnetic field using the changes in the vector angle during a reentry wave, demonstrating that the MCG pattern can be helpful for theoretical analysis of reentry waves.     

The role of photon scattering in optical signal distortion during arrhythmia and defibrillation

Optical mapping of arrhythmias and defibrillation provides important insights; however, a limitation of the technique is signal distortion due to photon scattering. The goal of this experimental/simulation study is to investigate the role of three-dimensional photon scattering in optical signal distortion during ventricular tachycardia (VT) and defibrillation. A three-dimensional realistic bidomain rabbit ventricular model was combined with a model of photon transport. Shocks were applied via external electrodes to induce sustained VT, and transmembrane potentials (V(m)) were compared with synthesized optical signals (V(opt)). Fluorescent recordings were conducted in isolated rabbit hearts to validate simulation results. Results demonstrate that shock-induced membrane polarization magnitude is smaller in V(opt) and in experimental signals as compared to V(m). This is due to transduction of potentials from weakly polarized midmyocardium to the epicardium. During shock-induced reentry and in sustained VT, photon scattering, combined with complex wavefront dynamics, results in optical action potentials near a filament exhibiting i), elevated resting potential, ii), reduced amplitude relative to pacing, and iii), dual-humped morphologies. A shift of up to 4 mm in the phase singularity location was observed in V(opt) maps when compared to V(m). This combined experimental/simulation study provides an interpretation of optical recordings during VT and defibrillation.     

High-density mapping of pulmonary veins and left atrium during ibutilide administration in a canine model of sustained atrial fibrillation

Ibutilide can prolong refractory period and terminate reentry. Whether ibutilide has the same effects on pulmonary vein (PV) focal discharge (FD) is unclear. We induced sustained atrial fibrillation (AF) in seven dogs by rapid left atrial (LA) pacing for 74 +/- 46 days. Ibutilide was repeatedly infused until it terminated AF (0.02 +/- 0.01 mg/kg) or when a cumulative dose was reached (0.04 mg/kg). High-resolution computerized epicardial mapping was performed. We found intermittent FD at the PVs and reentry at the PV-LA junction during AF. Ibutilide increased the cycle length of consecutive reentry from 97 +/- 13 to 112 +/- 18 ms and increased FD from 96 +/- 7 to 113 +/- 9 ms. In four dogs with both FD and reentry at the PVs, the incidence of reentry decreased from 3.5 +/- 1.9/s at baseline to 2.2 +/- 1.8/s after ibutilide administration. However, the incidence of FD remained unchanged. The conducted wave fronts between PV and LA were significantly reduced by ibutilide (10.4 +/- 2.0/s vs. 8.0 +/- 1.6/s). The ibutilide dose needed to terminate AF correlated negatively with the baseline effective refractory period of PV and LA. We conclude that ibutilide reduces reentrant wave fronts but not PV FD in a canine model of pacing-induced sustained AF. These findings suggest that the PV FD during AF is due to nonreentrant mechanisms. High doses of ibutilide may completely terminate all reentrant activity, converting AF to PV tachycardia before the resumption of sinus rhythm.     

Multiparametric Optical Mapping of the Langendorff-perfused Rabbit Heart

Optical imaging and fluorescent probes have significantly advanced research methodology in the field of cardiac electrophysiology in ways that could not have been accomplished by other approaches¹. With the use of the calcium- and voltage-sensitive dyes, optical mapping allows measurement of transmembrane action potentials and calcium transients with high spatial resolution without the physical contact with the tissue. This makes measurements of the cardiac electrical activity possible under many conditions where the use of electrodes is inconvenient or impossible¹. For example, optical recordings provide accurate morphological changes of membrane potential during and immediately after stimulation and defibrillation, while conventional electrode techniques suffer from stimulus-induced artifacts during and after stimuli due to electrode polarization¹. The Langendorff-perfused rabbit heart is one of the most studied models of human heart physiology and pathophysiology. Many types of arrhythmias observed clinically could be recapitulated in the rabbit heart model. It was shown that wave patterns in the rabbit heart during ventricular arrhythmias, determined by effective size of the heart and the wavelength of reentry, are very similar to that in the human heart². It was also shown that critical aspects of excitation-contraction (EC) coupling in rabbit myocardium, such as the relative contribution of sarcoplasmic reticulum (SR), is very similar to human EC coupling³. Here we present the basic procedures of optical mapping experiments in Langendorff-perfused rabbit hearts, including the Langendorff perfusion system setup, the optical mapping systems setup, the isolation and cannulation of the heart, perfusion and dye-staining of the heart, excitation-contraction uncoupling, and collection of optical signals. These methods could be also applied to the heart from species other than rabbit with adjustments to flow rates, optics, solutions, etc.Two optical mapping systems are described. The panoramic mapping system is used to map the entire epicardium of the rabbit heart^4-7. This system provides a global view of the evolution of reentrant circuits during arrhythmogenesis and defibrillation, and has been used to study the mechanisms of arrhythmias and antiarrhythmia therapy^8,9. The dual mapping system is used to map the action potential (AP) and calcium transient (CaT) simultaneously from the same field of view^10-13. This approach has enhanced our understanding of the important role of calcium in the electrical alternans and the induction of arrhythmia^14-16.     

Investigating a Novel Activation-Repolarisation Time Metric to Predict Localised Vulnerability to Reentry Using Computational Modelling

Exit sites associated with scar-related reentrant arrhythmias represent important targets for catheter ablation therapy. However, their accurate location in a safe and robust manner remains a significant clinical challenge. We recently proposed a novel quantitative metric (termed the Reentry Vulnerability Index, RVI) to determine the difference between activation and repolarisation intervals measured from pairs of spatial locations during premature stimulation to accurately locate the critical site of reentry formation. In the clinic, the method showed potential to identify regions of low RVI corresponding to areas vulnerable to reentry, subsequently identified as ventricular tachycardia (VT) circuit exit sites. Here, we perform an in silico investigation of the RVI metric in order to aid the acquisition and interpretation of RVI maps and optimise its future usage within the clinic. Within idealised 2D sheet models we show that the RVI produces lower values under correspondingly more arrhythmogenic conditions, with even low resolution (8 mm electrode separation) recordings still able to locate vulnerable regions. When applied to models of infarct scars, the surface RVI maps successfully identified exit sites of the reentrant circuit, even in scenarios where the scar was wholly intramural. Within highly complex infarct scar anatomies with multiple reentrant pathways, the identified exit sites were dependent upon the specific pacing location used to compute the endocardial RVI maps. However, simulated ablation of these sites successfully prevented the reentry re-initiation. We conclude that endocardial surface RVI maps are able to successfully locate regions vulnerable to reentry corresponding to critical exit sites during sustained scar-related VT. The method is robust against highly complex and intramural scar anatomies and low resolution clinical data acquisition. Optimal location of all relevant sites requires RVI maps to be computed from multiple pacing locations.     

Comparison of the duration and power spectral changes of monopolar and bipolar M waves caused by alterations in muscle fibre conduction velocity

The muscle compound action potential (M wave) recorded under monopolar configuration reflects both the propagation of the action potentials along the muscle fibres and their extinction at the tendon. M waves recorded under a bipolar configuration contain less cross talk and noise than monopolar M waves, but they do not contain the entire informative content of the propagating potential. The objective of this study was to compare the effect of changes in muscle fibre conduction velocity (MFCV) on monopolar and bipolar M waves and how this effect depends on the distance between the recording electrodes and tendon. The study was based on a simulation approach and on an experimental investigation of the characteristics of surface M waves evoked in the vastus lateralis during 4-s step-wise isometric contractions in knee extension at 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90% MVC. The peak-to-peak duration (Durpp) and median frequency (Fmedian) of the M waves were calculated. For monopolar M waves, changes in Durpp and Fmedian produced by MFCV depended on the distance from the electrode to the tendon, whereas, for bipolar M waves, changes in Durpp and Fmedian were largely independent of the electrode-to-tendon distance. When the distance between the detection point and tendon lay between approximately 15 and 40 mm, changes in Durpp of bipolar M waves were more pronounced than those of distal monopolar M waves but less marked than those of proximal monopolar M waves, and the opposite occurred for Fmedian. Since, for bipolar M waves, changes in duration and power spectral features produced by alterations in MFCV are not influenced by the electrode-to-tendon distance, the bipolar electrode configuration is a preferable choice over monopolar arrangements to estimate changes in conduction velocity.