Monophasic action potentials generated by bidomain modeling as a tool for detecting cardiac repolarization times

Understanding repolarization in the intracardiac unipolar electrogram: A long-lasting controversy revisited

Background: The optimal way to determine repolarization time (RT) from the intracardiac unipolar electrogram (UEG) has been a topic of debate for decades. RT is typically determined by either the Wyatt method or the "alternative method," which both consider UEG T-wave slope, but differently. Objective: To determine the optimal method to measure RT on the UEG. Methods: Seven pig hearts surrounded by an epicardial sock with 100 electrodes were Langendorff-perfused with selective cannulation of the left anterior descending (LAD) coronary artery and submersed in a torso-shaped tank containing 256 electrodes on the torso surface. Repolarization was prolonged in the non-LAD-regions by infusing dofetilide and shortened in the LAD-region using pinacidil. RT was determined by the Wyatt (t_Wyatt) and alternative (t_Alt) methods, in both invasive (recorded with epicardial electrodes) and in non-invasive UEGs (reconstructed with electrocardiographic imaging). t_Wyatt and t_Alt were compared to local effective refractory period (ERP). Results: With contact mapping, mean absolute error (MAE) of t_Wyatt and t_Alt vs. ERP were 21 ms and 71 ms, respectively. Positive T-waves typically had an earlier ERP than negative T-waves, in line with theory. t_Wyatt -but not t_Alt-shortened by local infusion of pinacidil. Similar results were found for the non-invasive UEGs (MAE of t_Wyatt and t_Alt vs. ERP were 30 ms and 92 ms, respectively). Conclusion: The Wyatt method is the most accurate to determine RT from (non) invasive UEGs, based on novel and historical analyses. Using it to determine RT could unify and facilitate repolarization assessment and amplify its role in cardiac electrophysiology.

A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms

Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research.

Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study

The aim of this work is to investigate, by means of numerical simulations, the influence of myocardial deformation due to muscle contraction and relaxation on the cardiac repolarization process in presence of transmural intrinsic action potential duration (APD) heterogeneities. The three-dimensional electromechanical model considered consists of the following four coupled components: the quasi-static transversely isotropic finite elasticity equations for the deformation of the cardiac tissue; the active tension model for the intracellular calcium dynamics and cross-bridge binding; the anisotropic Bidomain model for the electrical current flow through the deforming cardiac tissue; the membrane model of ventricular myocytes, including stretch-activated channels. The numerical simulations are based on our finite element parallel solver, which employs Multilevel Additive Schwarz preconditioners for the solution of the discretized Bidomain equations and Newton-Krylov methods for the solution of the discretized non-linear finite elasticity equations. Our findings show that: (i) the presence of intrinsic transmural cellular APD heterogeneities is not fully masked by electrotonic current flow or by the presence of the mechanical deformation; (ii) despite the presence of transmural APD heterogeneities, the recovery process follows the activation sequence and there is no significant transmural repolarization gradient; (iii) with or without transmural APD heterogeneities, epicardial electrograms always display the same wave shape and discordance between the polarity of QRS complex and T-wave; (iv) the main effects of the mechanical deformation are an increase of the dispersion of repolarization time and APD, when computed over the total cardiac domain and over the endo- and epicardial surfaces, while there is a slight decrease along the transmural direction.

Action potential wavelength restitution predicts alternans and arrhythmia in murine Scn5a(+/-) hearts

Reductions in cardiac action potential wavelength, and the consequent wavebreak, have been implicated in arrhythmogenesis. Tachyarrhythmias are more common in the Brugada syndrome, particularly following pharmacological challenge, previously modelled using Scn5a^+/− murine hearts. Propagation latencies and action potential durations (APDs) from monophasic action potential recordings were used to assess wavelength changes with heart rate in Langendorff-perfused wild-type (WT) and Scn5a^+/− hearts. Recordings were obtained from right (RV) and left (LV) ventricular, epicardial and endocardial surfaces during incremental pacing, before and following flecainide or quinidine challenge. Conduction velocities (θ′), action potential wavelengths (λ′= APD ×θ′), and their corresponding alternans depended non-linearly upon diastolic interval (DI). Maximum θ′ was lower in Scn5a^+/− RV epicardium than endocardium. Flecainide further reduced θ′, accentuating this RV conduction block. Quinidine reduced maximum θ′ in WT and caused earlier conduction failure in the RV of both Scn5a^+/− and WT. Use of recovery wavelengths (λ′₀= DI ×θ′) rather than DI, provided novel λ restitution plots of λ′ against λ′₀, which sum to a basic cycle distance permitting feedback analysis. λ′ restitution gradient better correlated with alternans magnitude than either APD or θ restitution gradient. The large differences in θ′ and APD restitution contrasted with minor differences in maximum λ′ between epi- and endocardia of untreated hearts, and quinidine-treated WT hearts. Strikingly, all regions and conditions converged to a common instability point, implying a conserved relationship. Flecainide or quinidine decreased the pacing rates at which this occurred, through reducing basic cycle distance, in the Scn5a^+/− RV epicardium, directly predictive of its arrhythmic phenotype.

Nonlinearity between action potential alternans and restitution, which both predict ventricular arrhythmic properties in Scn5a+/- and wild-type murine hearts

Electrocardiographic QT- and T-wave alternans, presaging ventricular arrhythmia, reflects compromised adaptation of action potential (AP) duration (APD) to altered heart rate, classically attributed to incomplete Na(v)1.5 channel recovery prior to subsequent stimulation. The restitution hypothesis suggests a function whose slope directly relates to APD alternans magnitude, predicting a critical instability condition, potentially generating arrhythmia. The present experiments directly test for such correlations among arrhythmia, APD alternans and restitution. Mice haploinsufficient in the Scn5a, cardiac Na(+) channel gene (Scn5a(+/-)), previously used to replicate Brugada syndrome, were used, owing to their established arrhythmic properties increased by flecainide and decreased by quinidine, particularly in right ventricular (RV) epicardium. Monophasic APs, obtained during pacing with progressively decrementing cycle lengths, were systematically compared at RV and left ventricular epicardial and endocardial recording sites in Langendorff-perfused Scn5a(+/-) and wild-type hearts before and following flecainide (10 μM) or quinidine (5 μM) application. The extent of alternans was assessed using a novel algorithm. Scn5a(+/-) hearts showed greater frequencies of arrhythmic endpoints with increased incidences of ventricular tachycardia, diminished by quinidine, and earlier onsets of ventricular fibrillation, particularly following flecainide challenge. These features correlated directly with increased refractory periods, specifically in the RV, and abnormal restitution and alternans properties in the RV epicardium. The latter variables were related by a unique, continuous higher-order function, rather than a linear relationship with an unstable threshold. These findings demonstrate a specific relationship between alternans and restitution, as well as confirming their capacity to predict arrhythmia, but implicate mechanisms additional to the voltage feedback suggested in the restitution hypothesis.

Validation of a simple model for the morphology of the T wave in unipolar electrograms

Local unipolar electrograms (UEGs) permit assessment of local activation and repolarization times at multiple sites simultaneously. However, UEG-based indexes of local repolarization are still debated, in particular for positive T waves. Previous experimental and computer modeling studies have not been able to terminate the debate. In this study we validate a simple theoretical model of the UEG and use it to explain how repolarization statistics in the UEG relate to those in the action potential. The model reconstructs the UEG by taking the difference between an inverted local action potential and a position-independent remote signal. In normal tissue, this extremely simple model predicts T-wave morphology with surprising accuracy while explaining in a readily understandable way why the instant of repolarization is always related to the steepest upstroke of the UEG, both in positive and negative T waves, and why positive T waves are related to early repolarizing sites, whereas negative T waves are related to late repolarizing sites.

A reliability analysis of cardiac repolarization time markers

Only a limited number of studies have addressed the reliability of extracellular markers of cardiac repolarization time, such as the classical marker RT(eg) defined as the time of maximum upslope of the electrogram T wave. This work presents an extensive three-dimensional simulation study of cardiac repolarization time, extending the previous one-dimensional simulation study of a myocardial strand by Steinhaus [B.M. Steinhaus, Estimating cardiac transmembrane activation and recovery times from unipolar and bipolar extracellular electrograms: a simulation study, Circ. Res. 64 (3) (1989) 449]. The simulations are based on the bidomain - Luo-Rudy phase I system with rotational fiber anisotropy and homogeneous or heterogeneous transmural intrinsic membrane properties. The classical extracellular marker RT(eg) is compared with the gold standard of fastest repolarization time RT(tap), defined as the time of minimum derivative during the downstroke of the transmembrane action potential (TAP). Additionally, a new extracellular marker RT90(eg) is compared with the gold standard of late repolarization time RT90(tap), defined as the time when the TAP reaches 90% of its resting value. The results show a good global match between the extracellular and transmembrane repolarization markers, with small relative mean discrepancy (<or=1.6%) and high correlation coefficients (>or=0.92), ensuring a reasonably good global match between the associated repolarization sequences. However, large local discrepancies of the extracellular versus transmembrane markers may ensue in regions where the curvature of the repolarization front changes abruptly (e.g. near front collisions) or is negligible (e.g. where repolarization proceeds almost uniformly across fiber). As a consequence, the spatial distribution of activation-recovery intervals (ARI) may provide an inaccurate estimate of (and weakly correlated with) the spatial distribution of action potential durations (APD).