Restitution Curve Research Articles

Restitution Slope Affects the Outcome of Dominant Frequency Ablation in Persistent Atrial Fibrillation: CUVIA-AF2 Post-Hoc Analysis Based on Computational Modeling Study.

IntroductionAlthough the dominant frequency (DF) localizes the reentrant drivers and the maximal slope of the action potential duration (APD) restitution curve (Smax) reflects the tendency of the wave-break, their interaction has never been studied. We hypothesized that DF ablation has different effects on atrial fibrillation (AF) depending on Smax.MethodsWe studied the DF and Smax in 25 realistic human persistent AF model samples (68% male, 60 ± 10 years old). Virtual AF was induced by ramp pacing measuring Smax, followed by spatiotemporal DF evaluation for 34 s. We assessed the DF ablation effect depending on Smax in both computational modeling and a previous clinical trial, CUVIA-AF (170 patients with persistent AF, 70.6% male, 60 ± 11 years old).ResultsMean DF had an inverse relationship with Smax regardless of AF acquisition timing (p < 0.001). Virtual DF ablations increased the defragmentation rate compared to pulmonary vein isolation (PVI) alone (p = 0.015), especially at Smax <1 (61.5 vs. 7.7%, p = 0.011). In post-DF ablation defragmentation episodes, DF was significantly higher (p = 0.002), and Smax was lower (p = 0.003) than in episodes without defragmentation. In the post-hoc analysis of CUVIA-AF2, we replicated the inverse relationship between Smax and DF (r = −0.47, p < 0.001), and we observed better rhythm outcomes of clinical DF ablations in addition to a PVI than of empirical PVI at Smax <1 [hazard ratio 0.45, 95% CI (0.22–0.89), p = 0.022; log-rank p = 0.021] but not at ≥ 1 (log-rank p = 0.177).ConclusionWe found an inverse relationship between DF and Smax and the outcome of DF ablation after PVI was superior at the condition with Smax <1 in both in-silico and clinical trials.

Interpretable machine learning of action potential duration restitution properties in single-cell models of atrial cardiomyocyte

Action potential duration (APD) restitution curve and its maximal slope (Smax) reflect single cell-level instability for inducing chaotic heart rhythms. Conventional parameter sensitivity analysis often fails to describe nonlinear relationships between ion channel parameters and electrophysiological phenotypes, such as Smax. We explored the parameter-phenotype mapping in in silico atrial cell models through interpretable machine learning (ML) approaches. We generated a population of 5,000 single-cell atrial cell models by log-uniformly sampling parameter combinations with the Sobol sequence scheme.

Effects of dynamic change of action potential on evolution behavior of spiral wave

It is observed in cardiac patients that the steepnesses of action potential duration (APD) restitution curve of cardiomyocytes in different regions of the ventricle are significantly different from region to region. However, the steep APD restitution curve can either lead the spiral wave to break up and set up the ventricular fibrillation in certain conditions or result in no breakup of spiral wave in other conditions. The relationship between the dynamic behavior of spiral wave and steep APD restitution curve is still not completely clear. Therefore, further research is needed. In this paper, a two-dimensional excitable medium cellular automata model is used to study the influences of the APD restitution curves with different steepnesses on the dynamic behavior of spiral wave. Numerical simulation results show that the steep APD restitution curve can stabilize the meandering spiral wave, causing the stable spiral wave to roam or break, and even to disappear. When the total average slope of APD restitution curve is greater than 1, it is observed that spiral wave may be still stable or meandering. When the total average slope of APD restitution curve is much smaller than 1, the breakup of spiral waves may occur. Three types of spiral wave breakups are observed. They are the Doppler instability, Eckhaus instability, and APD alternation. The Doppler instability and Eckhaus instability are related to the total average slope of APD restitution curve greater than 1, and the spiral wave breakup caused by APD alternans may occur when the total average slope of APD restitution curve is much smaller than 1. When the total average slope of APD restitution curve is greater than 1, the phenomena that spiral waves disappear by meandering out of the system boundary and conduction barriers are observed. In addition, we also find that increasing cellular APD is beneficial to preventing spiral wave from breaking up. The physical mechanisms behind those phenomena are heuristically analyzed.

Clinical Usefulness of Virtual Ablation Guided Catheter Ablation of Atrial Fibrillation Targeting Restitution Parameter-Guided Catheter Ablation: CUVIA-REGAB Prospective Randomized Study.

We investigated whether extra-pulmonary vein (PV) ablation targeting a high maximal slope of the action potential duration restitution curve (Smax) improves the rhythm outcome of persistent atrial fibrillation (PeAF) ablation. In this open-label, multi-center, randomized, and controlled trial, 178 PeAF patients were randomized with 1:1 ratio to computational modeling-guided virtual Smax ablation (V-Smax) or empirical ablation (E-ABL) groups. Smax maps were generated by computational modeling based on atrial substrate maps acquired during clinical procedures in sinus rhythm. Smax maps were generated during the clinical PV isolation (PVI). The V-Smax group underwent an additional extra-PV ablation after PVI targeting the virtual high Smax sites. After a mean follow-up period of 12.3±5.2 months, the clinical recurrence rates (25.6% vs. 23.9% in the V-Smax and the E-ABL group, p=0.880) or recurrence appearing as atrial tachycardia (11.1% vs. 5.7%, p=0.169) did not differ between the 2 groups. The post-ablation cardioversion rate was higher in the V-Smax group than E-ABL group (14.4% vs. 5.7%, p=0.027). Among antiarrhythmic drug-free patients (n=129), the AF freedom rate was 78.7% in the V-Smax group and 80.9% in the E-ABL group (p=0.776). The total procedure time was longer in the V-Smax group (p=0.008), but no significant difference was found in the major complication rates (p=0.497) between the groups. Unlike a dominant frequency ablation, the computational modeling-guided V-Smax ablation did not improve the rhythm outcome of the PeAF ablation and had a longer procedure time. ClinicalTrials.gov Identifier: NCT02558699.

A computationally efficient dynamic model of human epicardial tissue.

We present a new phenomenological model of human ventricular epicardial cells and we test its reentry dynamics. The model is derived from the Rogers-McCulloch formulation of the FitzHugh-Nagumo equations and represents the total ionic current divided into three contributions corresponding to the excitatory, recovery and transient outward currents. Our model reproduces the main characteristics of human epicardial tissue, including action potential amplitude and morphology, upstroke velocity, and action potential duration and conduction velocity restitution curves. The reentry dynamics is stable, and the dominant period is about 270 ms, which is comparable to clinical values. The proposed model is the first phenomenological model able to accurately resemble human experimental data by using only 3 state variables and 17 parameters. Indeed, it is more computationally efficient than existing models (i.e., almost two times faster than the minimal ventricular model). Beyond the computational efficiency, the low number of parameters facilitates the process of fitting the model to the experimental data.

Isolated superfused rats atrial model for the investigation of atrial fibrillation mechanisms and treatment

Abstract Background One of the major barriers to an improved mechanistic understanding of atrial fibrillation (AF), and thus in the pipeline of drug development, has been a lack of appropriate tissue models, especially in small animals. Aim We propose an advanced anatomical ex-vivo model based on rat atria for acute assessment of AF susceptibility. This novel model could yield a better understanding of arrhythmia mechanisms as well as the development of potential therapeutic strategies for the prevention or termination of atrial arrhythmias. Methods Wistar rats atria (N=25) were isolated, flattened and pinned to a custom-made silicon plate. Atria were superfused with an oxygenized Tyrode's solution. Tissues were then loaded with a voltage-sensitive dye and mapped using a high-resolution optical mapping system. AF was induced with 1uM carbamylcholine (N=23) coupled with pacing maneuvers and treated with 30uM Vernakalant (N=10) or 10uM Flecainide (N=10). Finally, the feasibility of a new ablation technique (electroporation) was evaluated. Results Optical mapping results suggested that the superfusion procedure led to a fast atrial recovery. Sinus activity was conserved for all atria for a long period. All the anatomical landmarks were clearly visualized. The acquired optical signals were analyzed during sinus rhythm and pacing, which allowed the creation of detailed activation maps and measurements of action potential duration (APD) and conduction velocity (CV) at different pacing rates. The resulting APD restitution curves revealed electrical excitation at high pacing rates (cycle length between 50ms and 300ms) with a relatively flattened curve. AF was successfully induced and optically mapping confirmed the presence of reentrant activity. AF was successfully treated using Vernacalant and Flecainide. Finally, we demonstrated the feasibility of a new ablation approach (electroporation) for creation of a continuous linear lesion serving as a functional block. Conclusion The isolated superfused atria model, coupled with voltage-sensitive dyes, can be utilized for long-term high-resolution functional imaging of the atria during sinus rhythm, pacing and arrhythmogenic activity. This allows the study of the atrial electrophysiological properties, the mechanisms involved in AF initiation, perpetuation, and termination as well as the study of drug and new ablation modalities. Funding Acknowledgement Type of funding sources: Public grant(s) – EU funding. Main funding source(s): European Research Council (ERC) Spontaneous activation of isolated atria

Neuromodulation With Thoracic Dorsal Root Ganglion Stimulation Reduces Ventricular Arrhythmogenicity.

Introduction: Sympathetic hyperactivity is strongly associated with ventricular arrhythmias and sudden cardiac death. Neuromodulation provides therapeutic options for ventricular arrhythmias by modulating cardiospinal reflexes and reducing sympathetic output at the level of the spinal cord. Dorsal root ganglion stimulation (DRGS) is a recent neuromodulatory approach; however, its role in reducing ventricular arrhythmias has not been evaluated. The aim of this study was to determine if DRGS can reduce cardiac sympathoexcitation and the indices for ventricular arrhythmogenicity induced by programmed ventricular extrastimulation. We evaluated the efficacy of thoracic DRGS at both low (20 Hz) and high (1 kHz) stimulation frequencies.Methods: Cardiac sympathoexcitation was induced in Yorkshire pigs (n = 8) with ventricular extrastimulation (S1/S2 pacing), before and after DRGS. A DRG-stimulating catheter was placed at the left T2 spinal level, and animals were randomized to receive low-frequency (20 Hz and 0.4 ms) or high-frequency (1 kHz and 0.03 ms) DRGS for 30 min. High-fidelity cardiac electrophysiological recordings were performed with an epicardial electrode array measuring the indices of ventricular arrhythmogenicity—activation recovery intervals (ARIs), electrical restitution curve (Smax), and Tpeak–Tend interval (Tp-Te interval).Results: Dorsal root ganglion stimulation, at both 20 Hz and 1 kHz, decreased S1/S2 pacing-induced ARI shortening (20 Hz DRGS −21±7 ms, Control −50±9 ms, P = 0.007; 1 kHz DRGS −13 ± 2 ms, Control −46 ± 8 ms, P = 0.001). DRGS also reduced arrhythmogenicity as measured by a decrease in Smax (20 Hz DRGS 0.5 ± 0.07, Control 0.7 ± 0.04, P = 0.006; 1 kHz DRGS 0.5 ± 0.04, Control 0.7 ± 0.03, P = 0.007), and a decrease in Tp-Te interval/QTc (20 Hz DRGS 2.7 ± 0.13, Control 3.3 ± 0.12, P = 0.001; 1 kHz DRGS 2.8 ± 0.08, Control; 3.1 ± 0.03, P = 0.007).Conclusions: In a porcine model, we show that thoracic DRGS decreased cardiac sympathoexcitation and indices associated with ventricular arrhythmogenicity during programmed ventricular extrastimulation. In addition, we demonstrate that both low-frequency and high-frequency DRGS can be effective neuromodulatory approaches for reducing cardiac excitability during sympathetic hyperactivity.

Spatial Changes in the Atrial Fibrillation Wave-Dynamics After Using Antiarrhythmic Drugs: A Computational Modeling Study.

Background: We previously reported that a computational modeling-guided antiarrhythmic drug (AAD) test was feasible for evaluating multiple AADs in patients with atrial fibrillation (AF). We explored the anti-AF mechanisms of AADs and spatial change in the AF wave-dynamics by a realistic computational model.Methods: We used realistic computational modeling of 25 AF patients (68% male, 59.8 ± 9.8 years old, 32.0% paroxysmal AF) reflecting the anatomy, histology, and electrophysiology of the left atrium (LA) to characterize the effects of five AADs (amiodarone, sotalol, dronedarone, flecainide, and propafenone). We evaluated the spatial change in the AF wave-dynamics by measuring the mean dominant frequency (DF) and its coefficient of variation [dominant frequency-coefficient of variation (DF-COV)] in 10 segments of the LA. The mean DF and DF-COV were compared according to the pulmonary vein (PV) vs. extra-PV, maximal slope of the restitution curves (Smax), and defragmentation of AF.Results: The mean DF decreased after the administration of AADs in the dose dependent manner (p < 0.001). Under AADs, the DF was significantly lower (p < 0.001) and COV-DF higher (p = 0.003) in the PV than extra-PV region. The mean DF was significantly lower at a high Smax (≥1.4) than a lower Smax condition under AADs. During the episodes of AF defragmentation, the mean DF was lower (p < 0.001), but the COV-DF was higher (p < 0.001) than that in those without defragmentation.Conclusions: The DF reduction with AADs is predominant in the PVs and during a high Smax condition and causes AF termination or defragmentation during a lower DF and spatially unstable (higher DF-COV) condition.

The development of L-type Ca2+ current mediated alternans does not depend on the restitution slope in canine ventricular myocardium

Cardiac alternans have crucial importance in the onset of ventricular fibrillation. The early explanation for alternans development was the voltage-driven mechanism, where the action potential (AP) restitution steepness was considered as crucial determining factor. Recent results suggest that restitution slope is an inadequate predictor for alternans development, but several studies still claim the role of membrane potential as underlying mechanism of alternans. These controversial data indicate that the relationship of restitution and alternans development is not completely understood. APs were measured by conventional microelectrode technique from canine right ventricular papillary muscles. Ionic currents combined with fluorescent measurements were recorded by patch-clamp technique. APs combined with fluorescent measurements were monitored by sharp microelectrodes. Rapid pacing evoked restitution-independent AP duration (APD) alternans. When non-alternating AP voltage command was used, Ca2+i-transient (CaT) alternans were not observed. When alternating rectangular voltage pulses were applied, CaT alternans were proportional to ICaL amplitude alternans. Selective ICaL inhibition did not influence the fast phase of APD restitution. In this study we found that ICaL has minor contribution in shaping the fast phase of restitution curve suggesting that ICaL—if it plays important role in the alternans mechanism—could be an additional factor that attenuates the reliability of APD restitution slope to predict alternans.

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators.

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel “restitution curve emulators” as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

Predicting the duration of chaotic transients in excitable media

The spatiotemporal dynamics of excitable media may exhibit chaotic transients. We investigate this transient chaos in the 2D Fenton–Karma model describing the propagation of electrical excitation waves in cardiac tissue and compute the average duration of chaotic transients in dependence on model parameter values. Furthermore, other characteristics like the dominant frequency, the size of the excitable gap, pseudo ECGs, the number of phase singularities and parameters characterizing the action potential duration restitution curve are determined and it is shown that these quantities can be used to predict the average transient time using polynomial regression.

Arrhythmogenic Effects of Genetic Mutations Affecting Potassium Channels in Human Atrial Fibrillation: A Simulation Study.

Genetic mutations in genes encoding for potassium channel protein structures have been recently associated with episodes of atrial fibrillation in asymptomatic patients. The aim of this study is to investigate the potential arrhythmogenicity of three gain-of-function mutations related to atrial fibrillation—namely, KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M—using modeling and simulation of the electrophysiological activity of the heart. A genetic algorithm was used to tune the parameters’ value of the original ionic currents to reproduce the alterations experimentally observed caused by the mutations. The effects on action potentials, ionic currents, and restitution properties were analyzed using versions of the Courtemanche human atrial myocyte model in different tissues: pulmonary vein, right, and left atrium. Atrial susceptibility of the tissues to spiral wave generation was also investigated studying the temporal vulnerability. The presence of the three mutations resulted in an overall more arrhythmogenic substrate. Higher current density, action potential duration shortening, and flattening of the restitution curves were the major effects of the three mutations at the single-cell level. The genetic mutations at the tissue level induced a higher temporal vulnerability to the rotor’s initiation and progression, by sustaining spiral waves that perpetuate until the end of the simulation. The mutation with the highest pro-arrhythmic effects, exhibiting the widest sustained VW and the smallest meandering rotor’s tip areas, was KCNE3-V17M. Moreover, the increased susceptibility to arrhythmias and rotor’s stability was tissue-dependent. Pulmonary vein tissues were more prone to rotor’s initiation, while in left atrium tissues rotors were more easily sustained. Re-entries were also progressively more stable in pulmonary vein tissue, followed by the left atrium, and finally the right atrium. The presence of the genetic mutations increased the susceptibility to arrhythmias by promoting the rotor’s initiation and maintenance. The study provides useful insights into the mechanisms underlying fibrillatory events caused by KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M and might aid the planning of patient-specific targeted therapies.

Repolarization characteristics indicate electrical instability in ventricular arrhythmia originating from papillary muscle

Abstract Funding Acknowledgements Type of funding sources: Public grant(s) – National budget only. Main funding source(s): DZHK (German Center for Cardiovascular Research) Introduction Cardiac arrhythmia originating from the papillary muscle can trigger ventricular fibrillation and cause sudden cardiac death even in the absence of structural heart disease. Yet, no clinical parameters are known to reflect the propensity of arrhythmia to degenerate into ventricular fibrillation. Purpose We aimed at identification of parameters associated with degeneration of ventricular arrhythmia into ventricular fibrillation. Methods Ventricular arrhythmia was induced by aconitine injection into the papillary muscle of healthy sheep (n = 12) in an open-chest model. Endocardial high-density-mapping and epicardial mapping were performed. We determined repolarization time and activation-recovery-interval according to the Wyatt method. Results During focal arrhythmia faster conduction occurred in longitudinal and basal direction than transversal and apical direction. The electrical restitution curve, modelling relation of diastolic interval and activation-recovery-interval, is steeper in aconitine-induced ventricular arrhythmia than in ventricular pacing or sinus rhythm. Steeper restitution curves reflect electrical instability and propensity to degenerate into ventricular fibrillation. The repolarization-related parameters activation-recovery interval and repolarization time exhibit higher heterogeneity per beat in ventricular arrhythmia preceding degeneration into ventricular fibrillation. Repolarization time in relation to cycle length (RT/CL), which can easily be measured during electrophysiological procedures, differentiates self-limiting from degenerating arrhythmia with high specificity and sensitivity. Conclusion In structurally normal ovine hearts, the ratio of repolarization and cycle length (RT/CL) in ventricular arrhythmia and greater dispersion of repolarization are indicators of subsequent degeneration into ventricular fibrillation. While dispersion of repolarization is not easily acquired in clinical routine, a simple index (RT/CL) may be sufficient to differentiate between self-limiting and electrically instable arrhythmia with propensity to degenerate to ventricular fibrillation. Abstract Figure. Repolarization in VA model

CVAR-Seg: An Automated Signal Segmentation Pipeline for Conduction Velocity and Amplitude Restitution.

BackgroundRate-varying S1S2 stimulation protocols can be used for restitution studies to characterize atrial substrate, ionic remodeling, and atrial fibrillation risk. Clinical restitution studies with numerous patients create large amounts of these data. Thus, an automated pipeline to evaluate clinically acquired S1S2 stimulation protocol data necessitates consistent, robust, reproducible, and precise evaluation of local activation times, electrogram amplitude, and conduction velocity. Here, we present the CVAR-Seg pipeline, developed focusing on three challenges: (i) No previous knowledge of the stimulation parameters is available, thus, arbitrary protocols are supported. (ii) The pipeline remains robust under different noise conditions. (iii) The pipeline supports segmentation of atrial activities in close temporal proximity to the stimulation artifact, which is challenging due to larger amplitude and slope of the stimulus compared to the atrial activity.Methods and ResultsThe S1 basic cycle length was estimated by time interval detection. Stimulation time windows were segmented by detecting synchronous peaks in different channels surpassing an amplitude threshold and identifying time intervals between detected stimuli. Elimination of the stimulation artifact by a matched filter allowed detection of local activation times in temporal proximity. A non-linear signal energy operator was used to segment periods of atrial activity. Geodesic and Euclidean inter electrode distances allowed approximation of conduction velocity. The automatic segmentation performance of the CVAR-Seg pipeline was evaluated on 37 synthetic datasets with decreasing signal-to-noise ratios. Noise was modeled by reconstructing the frequency spectrum of clinical noise. The pipeline retained a median local activation time error below a single sample (1 ms) for signal-to-noise ratios as low as 0 dB representing a high clinical noise level. As a proof of concept, the pipeline was tested on a CARTO case of a paroxysmal atrial fibrillation patient and yielded plausible restitution curves for conduction speed and amplitude.ConclusionThe proposed openly available CVAR-Seg pipeline promises fast, fully automated, robust, and accurate evaluations of atrial signals even with low signal-to-noise ratios. This is achieved by solving the proximity problem of stimulation and atrial activity to enable standardized evaluation without introducing human bias for large data sets.

Severe Bradycardia Increases the Incidence and Severity of Torsade de Pointes Arrhythmias by Augmenting Preexistent Spatial Dispersion of Repolarization in the CAVB Dog Model.

IntroductionTorsade de pointes arrhythmias (TdP) in the chronic atrioventricular block (CAVB) dog model result from proarrhythmic factors, which trigger TdP and/or reinforce the arrhythmic substrate. This study investigated electrophysiological and arrhythmogenic consequences of severe bradycardia for TdP.MethodsDofetilide (25 μg/kg per 5 min) was administered to eight anesthetized, idioventricular rhythm (IVR) remodeled CAVB dogs in two serial experiments: once under 60 beats per minute (bpm), right ventricular apex paced (RVA60) conditions, once under more bradycardic IVR conditions. Recordings included surface electrocardiogram and short-term variability (STV) of repolarization from endocardial unipolar electrograms. TdP inducibility (three or more episodes within 10 min after start of dofetilide) and arrhythmic activity scores (AS) were established. Mapping experiments in 10 additional dogs determined the effect of lowering rate on STV and spatial dispersion of repolarization (SDR) in baseline.ResultsIVR-tested animals had longer baseline RR-interval (1,403 ± 271 ms) and repolarization intervals than RVA60 animals. Dofetilide increased STV similarly under both rhythm strategies. Nevertheless, TdP inducibility and AS were higher under IVR conditions (6/8 and 37 ± 27 vs. 1/8 and 8 ± 12 in RVA60, respectively, both p < 0.05). Mapping: Pacing from high (128 ± 10 bpm) to middle (88 ± 10 bpm) to experimental rate (61 ± 3 bpm) increased all electrophysiological parameters, including interventricular dispersion, due to steeper left ventricular restitution curves, and intraventricular SDR: maximal cubic dispersion from 60 ± 14 (high) to 69 ± 17 (middle) to 84 ± 22 ms (p < 0.05 vs. high and middle rate).ConclusionIn CAVB dogs, severe bradycardia increases the probability and severity of arrhythmic events by heterogeneously causing electrophysiological instability, which is mainly reflected in an increased spatial, and to a lesser extent temporal, dispersion of repolarization.

The transient outward potassium current plays a key role in spiral wave breakup in ventricular tissue.

Spiral wave reentry as a mechanism of lethal ventricular arrhythmias has been widely demonstrated in animal experiments and recordings from human hearts. It has been shown that in structurally normal hearts spiral waves are unstable, breaking up into multiple wavelets via dynamical instabilities. However, many of the second-generation action potential models give rise only to stable spiral waves, raising issues regarding the underlying mechanisms of spiral wave breakup. In this study, we carried out computer simulations of two-dimensional homogeneous tissues using five ventricular action potential models. We show that the transient outward potassium current (Ito), although it is not required, plays a key role in promoting spiral wave breakup in all five models. As the maximum conductance of Ito increases, it first promotes spiral wave breakup and then stabilizes the spiral waves. In the absence of Ito, speeding up the L-type calcium kinetics can prevent spiral wave breakup, however, with the same speedup kinetics, spiral wave breakup can be promoted by increasing Ito. Increasing Ito promotes single-cell dynamical instabilities, including action potential duration alternans and chaos, and increasing Ito further suppresses these action potential dynamics. These cellular properties agree with the observation that increasing Ito first promotes spiral wave breakup and then stabilizes spiral waves in tissue. Implications of our observations to spiral wave dynamics in the real hearts and action potential model improvements are discussed.NEW & NOTEWORTHY Spiral wave breakup manifesting as multiple wavelets is a mechanism of ventricular fibrillation. It has been known that spiral wave breakup in cardiac tissue can be caused by a steeply sloped action potential duration restitution curve, a property mainly determined by the recovery of L-type calcium current. Here, we show that the transient outward potassium current (Ito) is another current that plays a key role in spiral wave breakup, that is, spiral waves can be stable for low and high maximum Ito conductance but breakup occurs for intermediate maximum Ito conductance. Since Ito is present in normal hearts of many species and required for Brugada syndrome, it may play an important role in the spiral wave stability and arrhythmogenesis under both normal condition and Brugada syndrome.

Universal mechanisms for self-termination of rapid cardiac rhythm.

Excitable media sustain circulating waves. In the heart, sustained circulating waves can lead to serious impairment or even death. To investigate factors affecting the stability of such waves, we have used optogenetic techniques to stimulate a region at the apex of a mouse heart at a fixed delay after the detection of excitation at the base of the heart. For long delays, rapid circulating rhythms can be sustained, whereas for shorter delays, there are paroxysmal bursts of activity that start and stop spontaneously. By considering the dependence of the action potential and conduction velocity on the preceding recovery time using restitution curves, as well as the reduced excitability (fatigue) due to the rapid excitation, we model prominent features of the dynamics including alternation of the duration of the excited phases and conduction times, as well as termination of the bursts for short delays. We propose that this illustrates universal mechanisms that exist in biological systems for the self-termination of such activities.

Overdrive pacing of spiral waves in a model of human ventricular tissue

High-voltage electrical defibrillation remains the only reliable method of quickly controlling life-threatening cardiac arrhythmias. This paper is devoted to studying an alternative approach, low-voltage cardioversion (LVC), which is based on ideas from non-linear dynamics and aims to remove sources of cardiac arrhythmias by applying high-frequency stimulation to cardiac tissue. We perform a detailed in-silico study of the elimination of arrhythmias caused by rotating spiral waves in a TP06 model of human cardiac tissue. We consider three parameter sets with slopes of the APD restitution curve of 0.7, 1.1 and 1.4, and we study LVC at the baseline and under the blocking of INa and ICaL and under the application of the drugs verapamil and amiodarone. We show that pacing can remove spiral waves; however, its efficiency can be substantially reduced by dynamic instabilities. We classify these instabilities and show that the blocking of INa and the application of amiodarone increase the efficiency of the method, while the blocking of ICaL and the application of verapamil decrease the efficiency. We discuss the mechanisms and the possible clinical applications resulting from our study.

Critical repolarization gradients determine the induction of reentry-based torsades de pointes arrhythmia in models of long QT syndrome

Torsades de pointes arrhythmia is a potentially lethal polymorphic ventricular tachyarrhythmia (pVT) in the setting of long QT syndrome. Arrhythmia susceptibility is influenced by risk factors modifying repolarization. The purpose of this article was to characterize repolarization duration and heterogeneity in relation to pVT inducibility and maintenance. Sotalol was infused regionally or globally in isolated Langendorff blood-perfused pig hearts (N = 7) to create repolarization time (RT) heterogeneities. Programmed stimulation and epicardial activation and repolarization mapping were performed. The role of RT (heterogeneities) was studied in more detail using a computer model of the human heart. pVTs (n = 11) were inducible at a critical combination of RT and RT heterogeneities. The pVT cycle lengths were similar in the short and long RT regions. Short-lasting pVTs were maintained by focal activity while longer-lasting pVTs by reentry wandering along the interface between the 2 regions. Local restitution curves from the long and short RT regions crossed. This was associated with T-wave inversion at coupling intervals at either side of the crossing point. These experimental observations were confirmed by the computer simulations. pVTs are inducible within a critical range of RT and RT heterogeneities and are maintained by reentry wandering along the repolarization gradient. Double potentials localize at the core of the reentrant circuit and reflect phase singularities. RT gradient and T waves invert with short-coupled premature beats in the long RT region as a result of the crossing of the restitution curves allowing reentry initiation.

Proarrhythmic Effects of Electrolyte Imbalance in Virtual Human Atrial and Ventricular Cardiomyocytes.

Dialysis is prescribed to renal failure patients as a long-term chronic treatment. Whereas dialysis therapeutically normalizes serum electrolytes and removes small toxin molecules, it fails to alleviate fibroblast induced structural fibrosis, and unresponsive uremia. The simultaneous presence of altered electrolytes and fibrosis or uremia is thought to be pro-arrhythmogenic. This study explored potential arrhythmogenesis under pre-dialysis (high electrolyte levels) and post-dialysis (low physiological electrolyte levels) in the presence of fibrosis and uremia in human atrial and ventricular model cardiomyocytes.Two validated human cardiomyocyte models were used in this study that permitted simulation of cardiac atrial and ventricular detailed electrophysiology. Pathological conditions simulating active fibrosis and uremia were implemented in both models. Pre- and post-dialysis conditions were simulated using high and low electrolyte levels respectively. Arrythmogenesis was quantified by computing restitution curves that permitted identification of action potential duration and calcium transient alternans instabilities.In comparison to control conditions, fibrosis abbreviated action potential durations while uremia prolonged the same. Under pre-dialysis conditions, an elevation of serum electrolyte levels caused action potential durations to be abbreviated under both fibrosis and uremia. Alternans instability was observed in the ventricular cardiomyocyte model. Under post-dialysis conditions, lower levels of serum electrolytes promoted an abbreviated action potential duration under fibrosis but caused a large increase of the control and uremic action potential durations. Alternans instabilities were observed in the atrial cardiomyocyte model under post-dialysis conditions at physiological heart rates. The calcium transient restitution showed similar alternans instabilities.Co-existing conditions such as fibrosis and uremia in the presence of unphysiological electrolyte levels promote arrhythmogenesis and may require additional treatment to improve dialysis outcomes.Clinical Relevance. Knowledge of model response to clinically relevant conditions permits use of in silico modeling to better understand and dissect underlying arrhythmia mechanisms.