Spatial correlation of left atrial low voltage substrate in sinus rhythm versus atrial fibrillation: The rhythm specificity of atrial low voltage substrate.

INTRODUCTION: Improved sinus rhythm (SR) maintenance rates have been achieved in patients with persistent atrial fibrillation (AF) undergoing pulmonary vein isolation plus additional ablation of low voltage substrate (LVS) during SR. However, voltage mapping during SR may be hindered in persistent and long-persistent AF patients by immediate AF recurrence after electrical cardioversion. We assess correlations between LVS extent and location during SR and AF, aiming to identify regional voltage thresholds for rhythm-independent delineation/detection of LVS areas. (1) Identification of voltage dissimilarities between mapping in SR and AF. (2) Identification of regional voltage thresholds that improve cross-rhythm substrate detection. (3) Comparison of LVS between SR and native versus induced AF. METHODS: Forty-one ablation-naive persistent AF patients underwent high-definition (1 mm electrodes; >1200 left atrial (LA) mapping sites per rhythm) voltage mapping in SR and AF. Global and regional voltage thresholds in AF were identified which best match LVS < 0.5 mV and <1.0 mV in SR. Additionally, the correlation between SR-LVS with induced versus native AF-LVS was assessed. RESULTS: Substantial voltage differences (median: 0.52, interquartile range: 0.33-0.69, maximum: 1.19 mV) with a predominance of the posterior/inferior LA wall exist between the rhythms. An AF threshold of 0.34 mV for the entire left atrium provides an accuracy, sensitivity and specificity of 69%, 67%, and 69% to identify SR-LVS < 0.5 mV, respectively. Lower thresholds for the posterior wall (0.27 mV) and inferior wall (0.3 mV) result in higher spatial concordance to SR-LVS (4% and 7% increase). Concordance with SR-LVS was higher for induced AF compared to native AF (area under the curve[AUC]: 0.80 vs. 0.73). AF-LVS < 0.5 mV corresponds to SR-LVS < 0.97 mV (AUC: 0.73). CONCLUSION: Although the proposed region-specific voltage thresholds during AF improve the consistency of LVS identification as determined during SR, the concordance in LVS between SR and AF remains moderate, with larger LVS detection during AF. Voltage-based substrate ablation should preferentially be performed during SR to limit the amount of ablated atrial myocardium.

Efficient time splitting schemes for the monodomain equation in cardiac electrophysiology.

Approximating the fast dynamics of depolarization waves in the human heart described by the monodomain model is numerically challenging. Splitting methods for the PDE-ODE coupling enable the computation with very fine space and time discretizations. Here, we compare different splitting approaches regarding convergence, accuracy, and efficiency. Simulations were performed for a benchmark problem with the Beeler-Reuter cell model on a truncated ellipsoid approximating the left ventricle including a localized stimulation. For this configuration, we provide a reference solution for the transmembrane potential. We found a semi-implicit approach with state variable interpolation to be the most efficient scheme. The results are transferred to a more physiological setup using a bi-ventricular domain with a complex external stimulation pattern to evaluate the accuracy of the activation time for different resolutions in space and time.

The reverse mode of the Na+/Ca2+ exchanger contributes to the pacemaker mechanism in rabbit sinus node cells.

Sinus node (SN) pacemaking is based on a coupling between surface membrane ion-channels and intracellular Ca2+-handling. The fundamental role of the inward Na+/Ca2+ exchanger (NCX) is firmly established. However, little is known about the reverse mode exchange. A simulation study attributed important role to reverse NCX activity, however experimental evidence is still missing. Whole-cell and perforated patch-clamp experiments were performed on rabbit SN cells supplemented with fluorescent Ca2+-tracking. We established 2 and 8 mM pipette NaCl groups to suppress and enable reverse NCX. NCX was assessed by specific block with 1 µM ORM-10962. Mechanistic simulations were performed by Maltsev-Lakatta minimal computational SN model. Active reverse NCX resulted in larger Ca2+-transient amplitude with larger SR Ca2+-content. Spontaneous action potential (AP) frequency increased with 8 mM NaCl. When reverse NCX was facilitated by 1 µM strophantin the Ca2+i and spontaneous rate increased. ORM-10962 applied prior to strophantin prevented Ca2+i and AP cycle change. Computational simulations indicated gradually increasing reverse NCX current, Ca2+i and heart rate with increasing Na+i. Our results provide further evidence for the role of reverse NCX in SN pacemaking. The reverse NCX activity may provide additional Ca2+-influx that could increase SR Ca2+-content, which consequently leads to enhanced pacemaking activity.

The openCARP simulation environment for cardiac electrophysiology.

BACKGROUND AND OBJECTIVE: Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community. METHODS AND RESULTS: openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential. CONCLUSION: As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.

Specific Electrogram Characteristics Identify the Extra-Pulmonary Vein Arrhythmogenic Sources of Persistent Atrial Fibrillation - Characterization of the Arrhythmogenic Electrogram Patterns During Atrial Fibrillation and Sinus Rhythm.

Identification of atrial sites that perpetuate atrial fibrillation (AF), and ablation thereof terminates AF, is challenging. We hypothesized that specific electrogram (EGM) characteristics identify AF-termination sites (AFTS). Twenty-one patients in whom low-voltage-guided ablation after pulmonary vein isolation terminated clinical persistent AF were included. Patients were included if short RF-delivery for <8sec at a given atrial site was associated with acute termination of clinical persistent AF. EGM-characteristics at 21 AFTS, 105 targeted sites without termination and 105 non-targeted control sites were analyzed. Alteration of EGM-characteristics by local fibrosis was evaluated in a three-dimensional high resolution (100 µm)-computational AF model. AFTS demonstrated lower EGM-voltage, higher EGM-cycle-length-coverage, shorter AF-cycle-length and higher pattern consistency than control sites (0.49 ± 0.39 mV vs. 0.83 ± 0.76 mV, p < 0.0001; 79 ± 16% vs. 59 ± 22%, p = 0.0022; 173 ± 49 ms vs. 198 ± 34 ms, p = 0.047; 80% vs. 30%, p < 0.01). Among targeted sites, AFTS had higher EGM-cycle-length coverage, shorter local AF-cycle-length and higher pattern consistency than targeted sites without AF-termination (79 ± 16% vs. 63 ± 23%, p = 0.02; 173 ± 49 ms vs. 210 ± 44 ms, p = 0.002; 80% vs. 40%, p = 0.01). Low voltage (0.52 ± 0.3 mV) fractionated EGMs (79 ± 24 ms) with delayed components in sinus rhythm ('atrial late potentials', respectively 'ALP') were observed at 71% of AFTS. EGMs recorded from fibrotic areas in computational models demonstrated comparable EGM-characteristics both in simulated AF and sinus rhythm. AFTS may therefore be identified by locally consistent, fractionated low-voltage EGMs with high cycle-length-coverage and rapid activity in AF, with low-voltage, fractionated EGMs with delayed components/ 'atrial late potentials' (ALP) persisting in sinus rhythm.

Pharmacologic TWIK-Related Acid-Sensitive K+ Channel (TASK-1) Potassium Channel Inhibitor A293 Facilitates Acute Cardioversion of Paroxysmal Atrial Fibrillation in a Porcine Large Animal Model.

Background The tandem of P domains in a weak inward rectifying K+ channel (TWIK)-related acid-sensitive K+ channel (TASK-1; hK2P3.1) two-pore-domain potassium channel was recently shown to regulate the atrial action potential duration. In the human heart, TASK-1 channels are specifically expressed in the atria. Furthermore, upregulation of atrial TASK-1 currents was described in patients suffering from atrial fibrillation (AF). We therefore hypothesized that TASK-1 channels represent an ideal target for antiarrhythmic therapy of AF. In the present study, we tested the antiarrhythmic effects of the high-affinity TASK-1 inhibitor A293 on cardioversion in a porcine model of paroxysmal AF. Methods and Results Heterologously expressed human and porcine TASK-1 channels are blocked by A293 to a similar extent. Patch clamp measurements from isolated human and porcine atrial cardiomyocytes showed comparable TASK-1 currents. Computational modeling was used to investigate the conditions under which A293 would be antiarrhythmic. German landrace pigs underwent electrophysiological studies under general anesthesia. Paroxysmal AF was induced by right atrial burst stimulation. After induction of AF episodes, intravenous administration of A293 restored sinus rhythm within cardioversion times of 177±63 seconds. Intravenous administration of A293 resulted in significant prolongation of the atrial effective refractory period, measured at cycle lengths of 300, 400 and 500 ms, whereas the surface ECG parameters and the ventricular effective refractory period lengths remained unchanged. Conclusions Pharmacological inhibition of atrial TASK-1 currents exerts antiarrhythmic effects in vivo as well as in silico, resulting in acute cardioversion of paroxysmal AF. Taken together, these experiments indicate the therapeutic potential of A293 for AF treatment.

Mapping and Removing the Ventricular Far Field Component in Unipolar Atrial Electrograms.

OBJECTIVE: Unipolar intracardiac electrograms (uEGMs) measured inside the atria during electro-anatomic mapping contain diagnostic information about cardiac excitation and tissue properties. The ventricular far field (VFF) caused by ventricular depolarization compromises these signals. Current signal processing techniques require several seconds of local uEGMs to remove the VFF component and thus prolong the clinical mapping procedure. We developed an approach to remove the VFF component using data obtained during initial anatomy acquisition. METHODS: We developed two models which can approximate the spatio-temporal distribution of the VFF component based on acquired EGM data: Polynomial fit, and dipole fit. Both were benchmarked based on simulated cardiac excitation in two models of the human heart and applied to clinical data. RESULTS: VFF data acquired in one atrium were used to estimate model parameters. Under realistic noise conditions, a dipole model approximated the VFF with a median deviation of 0.029 mV, yielding a median VFF attenuation of 142. In a different setup, only VFF data acquired at distances of more than 5 mm to the atrial endocardium were used to estimate the model parameters. The VFF component was then extrapolated for a layer of 5 mm thickness lining the endocardial tissue. A median deviation of 0.082 mV (median VFF attenuation of 49x) was achieved under realistic noise conditions. CONCLUSION: It is feasible to model the VFF component in a personalized way and effectively remove it from uEGMs. SIGNIFICANCE: Application of our novel, simple and computationally inexpensive methods allows immediate diagnostic assessment of uEGM data without prolonging data acquisition.

Recurrence Quantification Analysis for Investigating Atrial Fibrillation Dynamics in a Heterogeneous Simulation Setup.

The outcomes of ablation targeting either reentry activations or fractionated activity during persistent atrial fibrillation (AF) therapy remain suboptimal due to, among others, the intricate underlying AF dynamics. In the present work, we sought to investigate such AF dynamics in a heterogeneous simulation setup using recurrence quantification analysis (RQA). AF was simulated in a spherical model of the left atrium, from which 412 unipolar atrial electrograms (AEGs) were extracted (2 s duration; 5 mm spacing). The phase was calculated using the Hilbert transform, followed by the identification of points of singularity (PS). Three regions were defined according to the occurrence of PSs: 1) no rotors; 2) transient rotors and; 3) long-standing rotors. Bipolar AEGs (1114) were calculated from pairs of unipolar nodes and bandpass filtered (30-300 Hz). The CARTO criterion (Biosense Webster) was used for AEGs classification (normal vs. fractionated). RQA attributes were calculated from the filtered bipolar AEGs: determinism (DET); recurrence rate (RR); laminarity (LAM). Sample entropy (SampEn) and dominant frequency (DF) were also calculated from the AEGs. Regions with longstanding rotors have shown significantly lower RQA attributes and SampEn when compared to the other regions, suggesting a higher irregular behaviour (P≤0.01 for all cases). Normal and fractionated AEGs were found in all regions (respectively; Region 1: 387 vs. 15; Region 2: 221 vs. 13; Region 3: 415 vs. 63). Region 1 vs. Region 3 have shown significant differences in normal AEGs (P≤0.0001 for all RQA attributes and SampEn), and significant differences in fractionated AEGs for LAM, RR and SampEn (P=0.0071, P=0.0221 and P=0.0086, respectively). Our results suggest the co-existence of normal and fractionated AEGs within long-standing rotors. RQA has unveiled distinct dynamic patterns-irrespective of AEGs classification-related to regularity structures and their nonstationary behaviour in a rigorous deterministic context.

Regional conduction velocity calculation from clinical multichannel electrograms in human atria.

BACKGROUND: During atrial fibrillation, heterogeneities and anisotropies result in a chaotic propagation of the depolarization wavefront. The electrophysiological parameter called conduction velocity (CV) influences the propagation pattern over the atrium. We present a method that determines the regional CV for deformed catheter shapes, which result due to the catheter movement and changing wall contact. METHODS: The algorithm selects stable catheter positions, finds the local activation times (LAT), considers the wall contact and calculates all CV estimates within the area covered by the catheter. The method is evaluated with simulated data and then applied to four clinical data sets. Both sinus rhythm activity as well as depolarization wavefronts initiated by stimulation are analyzed. The regional CV is compared with the fractionation duration (FD) and peak-to-peak (P2P) voltages. A speed of 0.5 m/s was defined to create the simulated LAT. RESULTS: After analyzing the simulated LAT with clinical catheter spatial coordinates, the median CV of 0.5 m/s with an interquartile range of 0.22 and exact CV direction vectors were obtained. For clinical cases, the CV magnitude range of 0.08 m/s to 1.0 m/s was obtained. The P2P amplitude of 0.7 mV to 3.7 mV and the mean FD from 40.79ms to 48.66ms was obtained. The correlation of 0.86 was observed between CV and P2P amplitude, and 0.62 between CV and FD. CONCLUSION: In this paper, a method is presented and validated which calculates the CV for the deformed catheter and changing wall contact. In an exemplary clinical data set correlation between regional CV with FD and the P2P voltage was observed.