Early afterdepolarizations promote transmural reentry in ischemic human ventricles with reduced repolarization reserve

Early afterdepolarizations promote transmural reentry in ischemic human ventricles with reduced repolarization reserve Sara Dutta a , 1 , Ana Minchol e a , 1 , Ernesto Zacur b , T. Alexander Quinn c , Peter Taggart d , Blanca Rodriguez a , * a Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, United Kingdom b Institute of Biomedical Engineering, University of Oxford, Oxford, United Kingdom c Department of Physiology and Biophysics, Dalhousie University, Canada d Institute of Cardiovascular Science, University College London, Bars Heart Centre, United Kingdom article info Article history: Received 1 October 2015 Received in revised form 25 January 2016 Accepted 29 January 2016 Available online 2 February 2016 Keywords: Ischemic heart disease Computer-based model Ventricular arrhythmia Potassium channels Repolarization abstract Aims: Acute ischemia is a major cause of sudden arrhythmic death, further promoted by potassium current blockers. Macro-reentry around the ischemic region and early afterdepolarizations (EADs) caused by electrotonic current have been suggested as potential mechanisms in animal and isolated cell studies. However, ventricular and human-specific arrhythmia mechanisms and their modulation by repolarization reserve remain unclear. The goal of this paper is to unravel multiscale mechanisms underlying the modulation of arrhythmic risk by potassium current (I Kr ) block in human ventricles with acute regional ischemia. Methods and results: A human ventricular biophysically-detailed model, with acute regional ischemia is constructed by integrating experimental knowledge on the electrophysiological ionic alterations caused by coronary occlusion. Arrhythmic risk is evaluated by determining the vulnerable window (VW) for reentry following ectopy at the ischemic border zone. Macro-reentry around the ischemic region is the main reentrant mechanism in the ischemic human ventricle with increased repolarization reserve due to the ATP-sensitive potassium current (I K(ATP) ) activation. Prolongation of refractoriness by 4% caused by 30% I Kr reduction counteracts the establishment of macro-reentry and reduces the VW for reentry (by 23.5%). However, a further decrease in repolarization reserve (50% I Kr reduction) is less anti-arrhythmic despite further prolongation of refractoriness. This is due to the establishment of transmural reentry enabled by electrotonically-triggered EADs in the ischemic border zone. EADs are produced by L-type calcium current (I CaL ) reactivation due to prolonged low amplitude electrotonic current injected during the repolarization phase. Conclusions: Electrotonically-triggered EADs are identified as a potential mechanism facilitating intramural reentry in a regionally-ischemic human ventricles model with reduced repolarization reserve. ©2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 1. Introduction Acute ischemia associated with coronary heart disease is a leading cause of sudden arrhythmic death. Pro-arrhythmic mechanisms of acute ischemia have been extensively investigated, although often in animal models at the single cell or tissue level rather than in human ventricles. Seminal studies by Janse et al. in pig and dog hearts highlight the complexity of the pro-arrhythmic and spatio-temporally dynamic substrate in acute ischemia (Janse and Wit, 1989; Janse et al., 1980, 1979). Heterogeneity in excitability and repolarization properties across the border between the normal and ischemic areas leads to the establishment of reentry around the ischemic region following ectopic excitation. The same studies also showed intramural reentry in certain cases (highlighting the potential variability in the mechanisms), recorded with microelectrodes on and 4 and 8 mm below the epicardial surface in pig hearts affected by left anterior descending coronary artery occlusion. However, the mechanisms that determine reentry formation and intramural patterns in acute ischemia in the three- *Corresponding author. Department of Computer Science, University of Oxford, Parks Road, OX1 3QD, Oxford, United Kingdom. E-mail address: [email protected]x.ac.uk (B. Rodriguez). 1 Equal contribution first authors. Contents lists available at ScienceDirect Progress in Biophysics and Molecular Biology journal homepage: www.elsevier.com/locate/pbiomolbio http://dx.doi.org/10.1016/j.pbiomolbio.2016.01.008 0079-6107/©2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Progress in Biophysics and Molecular Biology 120 (2016) 236e248 dimensional human heart remain unclear, due to low resolution of intramural recordings. Electrotonic current flowing due to gradients in membrane potential in heterogeneous tissue has been hypothesized as a potentially-important pro-arrhythmic factor in acute ischemia (Coronel et al., 1991; Janse et al., 1980). One potential mechanism is that electrotonic current flow could promote the occurrence of early afterdepolarizations (EADs) in acute ischemia (Kumar and Joyner, 1994; Verkerk et al., 2000). Single cell studies showed EADs triggered by electrotonic current in human cardiomyocytes from failing hearts, but only when cells in phase 3 repolarization were coupled to very high potentials corresponding to plateau level (such as 0 mV) (Verkerk et al., 2000). Neighboring cells would therefore have to experience very different membrane potential levels for this to occur. Furthermore, electrotonic interactions in well coupled tissue (as is still the case in the early phase of acute ischemia) are known to decrease the propensity of EAD formation in tissue (Weiss et al., 2010; Xie et al., 2010), even if they can be triggered in single cells as in Verkerk et al. (2000). Therefore, it remains unclear what conditions are needed for EADs to be triggered by electrotonic current, as well as the role of EAD formation in the regionally ischemic human heart. Pro-arrhythmic mechanisms including both reentry and EAD formation strongly depend on the balance of currents during repolarization, and the specific degree of redundancy in repolarizing currents, termed repolarization reserve, available. This is affected by a variety of causes including inter-subject differences in ion channel density, drug effects and mutations. Individuals with reduced repolarization reserve are known to be at higher arrhythmic risk, particularly in the presence of a pro-arrhythmic substrate such as ischemia due to coronary heart disease as shown in clinical studies (Steinberg et al., 2014; Waldo et al., 1996). Furthermore, clinical studies show that interventions that reduce repolarizing currents such as class III drugs result in increased risk of longer-term hospitalization and death in patients with coronary heart disease (Steinberg et al., 2014). Studies, such as the Survival With ORal D-Sotalol (SWORD) trial, also demonstrate that administration of the hERG (human Ether- a-go-go Related Gene) blocker sotalol is associated with increased mortality in patients suffering from recent myocardial infarction, related to ischemic heart disease (Waldo et al., 1996). How changes in the degree of potassium current available during repolarization alters the ischemic human heart to increase arrhythmic risk remains unclear. Investigating the pro-arrhythmic substrate of acute ischemia and its alteration by potassium block experimentally is very challenging. Of particular difficulty is mapping the three-dimensional volume of the heart with enough spatio-temporal resolution to investigate complex and irregular intramural propagation patterns through the ventricular wall. Microelectrode recordings and optical mapping have both been used in animal studies and they both suffer from a different set of limitations. The former exhibit low resolution, whereas the latter provides recordings only of the tissue surface and with photon scattering averaging effects (Bishop et al., 2007, p. 2007). Acute ischemia is further complicated by the fact that the electrophysiological changes induced are dynamic in time and reproducibility of experimental findings is compromised. In this study, we investigate the mechanisms of increased arrhythmic risk in the heterogeneous substrate caused by acute myocardial ischemia for different degrees of reduced repolarization reserve using an anatomically-based biophysically-detailed human ventricles model. We hypothesize that in human ventricles with low repolarization reserve, gradients in membrane potential across the ischemic border zone promote the occurrence of electrotonicallytriggered EADs that contribute to the establishment of transmural micro-reentry. By analyzing high spatio-resolution simulation data, we unravel the three-dimensional mechanisms of EAD formation and propagation dynamics in the acutely-ischemic human ventricles, their modulation by repolarization reserve and their potential contribution to increased arrhythmic risk. 2. Methods 2.1. Human ventricular model in acute regional ischemia An anatomically-based multiscale model of regionally-ischemic human ventricles was developed based on extensive experimental recordings from (Carmeliet, 1999; Coronel et al., 1991, 1988; Durrer et al., 1970; Edvardsson et al., 1980; Glukhov et al., 2010; Lee et al., 1981; Spector et al., 1996; Sutton et al., 2000; Taggart et al., 2001; Wilensky et al., 1986), including biophysically-detailed membrane kinetics, ischemia induced heterogeneities and fiber orientation (based on the Streeter method (Streeter et al., 1969)). Human ventricle membrane kinetics were simulated with the ten Tusscher 2006 action potential (AP) model (TP06) (Tusscher and Panfilov, 2006) with an added ATP-sensitive potassium current (I K(ATP) ) (Michailova et al., 2005). The TP06 has been extensively used and shown to be suitable to simulate reentrant arrhythmias as well as ischemic electrophysiological alterations in human (Dutta et al., 2013; Kazbanov et al., 2014). An expanded version of the methods is available as Supplemental Material. 2.2. Heterogeneous electrophysiological substrate in the ischemic region Size and location of ischemic regions exhibit great variability. We represented the ischemic region in our model as affecting ~40% of the left ventricle, which is within the range of 13e72% reported by Lee et al. (1981).(Fig. 1A) and corresponds to greatest likelihood of ventricular fibrillation (Curtis, 1998). The ischemic region was modeled in agreement with experimental findings by Coronel et al. (1991); Rodríguez et al. (2006); Wilensky et al. (1986) including: (i) the ischemic central zone (ICZ); (ii) the border zone (BZ), (iii) and the endocardial BZ (a layer of endocardial tissue with access to oxygen and nutrients from blood in the ventricular cavity) (Coronel et al., 1988; Wilensky et al., 1986). In the ICZ, ionic alterations caused by acute ischemia were introduced (Coronel et al., 1991; Rodríguez et al., 2006; Wilensky et al., 1986), with a degree of ischemic severity corresponding to highest arrhythmic risk (Carmeliet, 1999; Janse and Wit, 1989; Kazbanov et al., 2014; Rodríguez et al., 2006) during early stages of ischemia (up to 15 min). The BZ included a linear gradient in electrophysiological properties from the ICZ to the normal zone (NZ) tissue as shown in experiments (Coronel et al.,1988; Wilensky et al., 1986). Calcium overload or gap junction changes were not included as they occur at later stages of ischemia (Carmeliet, 1999). The human model also included transmural heterogeneities in AP duration (APD) to reproduce the physiological inverse APD - activation time relationship (Fig. 1B and C) resulting in a positive T wave in the pseudo-ECG, as in data obtained in in vivo human hearts (Franz et al., 1987). More details can be found in the Supplemental Material. 2.3. Stimulation protocol Purkinje-like activation was simulated twice with a cycle length (CL) of 800 ms, by stimulating the endocardium (S1) to mimic the activation sequence in Durrer et al. (1970). This was followed by a premature excitation (S2), applied in a region close to the BZ prone to premature beats that mimics findings by Janse et al. (1980). The coupling interval (CI) of the premature stimulus (e.g., time interval S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248 237 difference between S1 and S2) was varied to quantify the vulnerability window (VW) of reentry. The latter was computed as the range of CIs that resulted in reentry in the human ventricles. As illustrated in Fig. 1, the human ventricles model reproduced the reported electrophysiological heterogeneity caused by regional ischemia, including elevated resting transmembrane potential (V m ) and short APD (Rodríguez et al., 2006) in the ICZ with respect to NZ, in agreement with in vivo human and ex situ animal studies (Carmeliet, 1999; Sutton et al., 2000,p.200;Taggart et al., 2001; Wilensky et al., 1986). Resting potential in the ICZ is 70 mV compared to 86 mV in normal tissue, and APD is 30% shorter in the ICZ than in the NZ, close to the 35% difference reported in human in vivo measurements (Sutton et al., 2000). The simulated pseudo-ECG was computed as in Gima and Rudy (2002), as the extracellular unipolar potential at a position 3.6 cm away from the epicardial surface (Fig. 1A), and yielded a physiological QT interval of ~400 ms and positive T wave, as described in in vivo human hearts (Franz et al., 1987). 2.4. Reduced repolarization reserve In order to investigate mechanistic implications of reduced Fig. 1. Electrophysiological heterogeneity in the anatomically-based model of the human ventricles in acute regional ischemia. A. Resting potential distribution with elevated potentials in ischemic tissue (green solid line), location of the pseudo-ECG probe (black asterisk) and ectopic stimulation site (magenta solid line). B. Action potential duration (APD) distribution highlighting APD shortening in the ischemic central zone (ICZ) (top) and examples of action potentials in normal zone (NZ), border zone (BZ) and ICZ from the locations marked with green dots in the APD map. C. Activation times distribution (top) and pseudo-ECG (bottom) exhibiting first the QRS complex, followed by a positive T wave (due to the inverse relationship between APD and activation times) and ischemia-induced ST elevation. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248238 repolarization reserve in acute ischemia, we considered three scenarios represented by a 0, 30 and 50% decrease of the rapidly activating delayed rectifier potassium current (I Kr ) conductance in the TP06. The simulations yield a 5% and 9% prolongation in QT in the pseudo-ECG, which is similar to the 7% and 12% increase found by Fossa et al. in healthy patients, 2e4 h after being given 160 mg and 320 mg of sotalol respectively (Fossa et al., 2007). One-dimensional simulations were conducted to characterize the effect of reduced repolarization reserve in NZ and ICZ tissue for varying degrees of I Kr reduction. APD was calculated at 90% repolarization at steady state for each CL. Effective refractory period (ERP) was evaluated after reaching steady state at each CL and calculated as the shortest CI to induce sustained propagation in the tissue. Conduction velocity (CV) was calculated at steady state from the time interval separating two activation times at two points in the center of the tissue separated by 2 mm. Fig. 2 illustrates the known electrophysiological effects of acute ischemia, including APD shortening, ERP prolongation and slow CV for all degrees of repolarization reserve. Note that no EADs occurred in any of the single cell and one-dimensional simulations presented in Fig. 2. 3. Results 3.1. Macro-reentrant patterns in the acutely-ischemic human ventricle model support experimental evidence in animal models Simulations show the establishment of macro-reentry around the ischemic region in the human ventricles following premature excitation during the VW for all levels of repolarization reserve. As illustrated in Fig. 3, in all cases, the macro-reentrant pattern was similar, including 3 main phases consistent with experimental patterns (Janse et al., 1980). Firstly, unidirectional conduction block occurs due to prolonged refractoriness in the ICZ, with propagation around this region via the NZ and BZ around the ICZ (first column in Fig. 3). Then, retrograde propagation occurs through the ICZ once tissue distal to premature excitation has recovered (middle column in Fig. 3). Finally, wavefront propagation reenters through the NZ (third column in Fig. 3). In all cases, the macro-reentrant pattern was similar and lasted for up to 3 reentry cycles after which it selfterminated. This may indicate that the establishment of sustained ventricular tachycardia and its transition to ventricular fibrillation may require the contribution of more than one ectopic beat and/or pro-arrhythmic substrate factors. Decreased repolarization reserve by 30% I Kr reduction leads to a decrease of the VW for reentry (65 ms with 30% I Kr reduction vs. 85 ms without I Kr decrease; CI ¼355e420 ms for 30% I Kr reduction vs. CI ¼335e420 ms for no I Kr decrease). The main anti-arrhythmic mechanism of decreasing I Kr is that prolonged ERP favors conduction block both in the NZ and ICZ at short CIs, therefore leading to either bidirectional block just after premature excitation or failure to sustain retrograde conduction in the ICZ. While the patterns of macro-reentry were similar when I Kr reduction was increased to 50%, the anti-arrhythmic activity of I Kr diminished, as shown by an increased width of the VW (80 ms; CI ¼355e435 ms). The causes of this increase in the VW with a further reduction in repolarization reserve are investigated below. 3.2. Electrotonically-triggered EADs facilitate transmural reentry near the ischemic border zone under conditions of reduced repolarization reserve Detailed analysis of the three-dimensional activation patterns in the simulations revealed that the increase in the VW under conditions of reduced repolarization reserve with 50% I Kr block was caused by transmural reentry supported by electrotonicallytriggered EADs near the ischemic BZ, as illustrated in Fig. 4. Firstly, two macro-reentrant cycles are established around the ICZ starting at 1200 ms and 1460 ms (as in Fig. 3). Then, prolonged refractoriness in the ICZ prevents retrograde propagation, and the macro-reentrant propagation dies out (Fig. 4, 1381 ms and 1775 ms panels). However, preserved excitability in the endocardial BZ sustains propagation of a small wavefront (Fig. 4,1775 ms snapshot, circle and arrow), which finds a nonrefractory tunnel for propagation via the lateral BZ, from endocardium to epicardium (Fig. 4, 1831 ms, snapshot, circle and arrow). Following breakthrough on the epicardium, the wavefront reenters back transmurally towards the endocardium and around the ischemic BZ (Fig. 4, 1870 ms, 1920 ms and 1950 ms snapshots). This pattern of transmural propagation is only seen for 50% I Kr block and the cellular-level mechanisms were investigated by analyzing the time course of V m at different locations from endocardium to epicardium in the region sustaining the transmural patterns of reentry, as indicated in Fig. 5A. The AP traces in Fig. 5B show the development of large transmural V m gradients (e.g. at 1831 ms) with fully depolarized epicardial cells (black traces) close to repolarizing endocardial tissue (green traces). The voltage gradients result in electrotonic current flowing towards the endocardial cells, and this facilitates the development of electrotonicallytriggered EADs (Fig. 5B, green and blue traces from 1920 ms). This is a consequence of a weak repolarization reserve due to 50% I Kr block and the electrotonic flow of current that allows for the increase in V m during the plateau of the AP. The electrotonic current that facilitates EAD formation was estimated to be between 1 and Fig. 2. Electrophysiological properties in a one-dimensional homogeneous fiber of human ventricular tissue under control and acute ischemic conditions for varying degrees of the rapidly activating delayed rectifying potassium current (I Kr ) reduction. Results show steady-state values for APD, effective refractory period (ERP) and conduction velocity (CV) for cycle lengths (CL) ranging from 350 to 1500 ms. Normal and ischemic cells were assigned extracellular potassium concentration ([K þ ] o ) of 5.4 and 8.5 mmol/L, ATP-sensitive potassium current (I K(ATP) ) activation of 0 and 5%, and peak conductance of fast sodium current (I Na ) and L-type calcium current (I CaL ) of 100 and 75%, respectively, of their original values in the TP06 model. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248 239 1.5 pA/pF in amplitude and 100e150 ms in duration by computing the diffusive part of the current in the monodomain equation. In the absence of EAD formation triggered by electrotonic current, epicardial and endocardial cells would be refractory at the same time, and this would lead to conduction failure. However, as shown in Fig. 5B, EAD formation prolongs the time during which the endocardial cells are depolarized (green traces), increasing the time available to achieve the full recovery of neighbouring tissue (black and blue traces) and the continuation of the propagation, as was suggested by Kumar and Joyner (1994); Pogwizd and Corr (1990); Yan et al. (2001). This mechanism explains the increased propensity of reentry and VW extension for 50% compared to 30% I Kr reduction, and is further illustrated in Fig. 6 through a comparison of simulations for strong I Kr (0% reduction) and 30% and 50% reduction (top, middle and bottom panels). In all three cases, propagation first proceeds around the ICZ, and then through the ICZ once its tissue has recovered. For 0% and 30% I Kr reduction, propagation completely engulfs the ventricles and reentry fails to be established. However, in the case of 50% I Kr reduction, prolonged refractoriness in the ICZ leads to conduction block close to the BZ, followed by endocardial propagation and EAD formation (1810 ms and Fig. 6A and B), leading to transmural reentry as shown in Figs. 4 and 5. More details can be found in the supplemental movies. 3.3. Prolonged low amplitude current may trigger EADs due to reactivation of the L-type calcium current gates In order to investigate the mechanisms underlying the formation of EADs in the acutely-ischemic human ventricles, we simulated the effect of a current injected during the repolarization phase on the human ventricular ionic dynamics and AP under normal and ischemic conditions with various degrees of repolarization reserve. Based on the conditions shown in Figs. 3, 4 and 6, our hypothesis was that a current of low amplitude and prolonged duration could trigger L-type calcium current (I CaL ) reactivation and EADs if injected during the repolarization phase. Results presented in Figs. 7 and 8 show the following EAD mechanisms. Firstly, that EADs can be triggered when a current of low amplitude (less than 1/20th of the threshold for AP triggering) is injected in the NZ, BZ and ICZ human ventricle cardiomyocytes with normal and reduced repolarization reserve (Fig. 7). This is consistent with results in our ventricular simulations showing EAD formation across the BZ for 50% I Kr block (Figs. 5 and 6). Secondly, analysis of the ionic mechanisms underlying EAD formation shows that reactivation of the I CaL current and its activation gate are the main driving force (Fig. 8, second and third row), while calcium release from the sarcoplasmic reticulum (as speculated by Verkerk et al. (2000)) is not involved (Fig. 8, fourth row). Thirdly, key to EAD formation is that the amplitude of the current injected is above a certain threshold. This is illustrated for example in the BZ and ICZ rows in Fig. 7A: a continuous 1.1 pA/pF current fails to develop EADs for normal and reduced repolarization reserve in the BZ and ICZ respectively, but they are triggered with an increased amplitude of 1.2 pA/pF. Furthermore, reduced repolarization reserve favors EAD occurrence, particularly when the prolonged current is applied during phase 3 of repolarization (i.e. for V m 20 mV) (Fig. 6B). Indeed, the last two columns show that the occurrence of EADs is facilitated by I Kr block. Low amplitude current applied at early stages of repolarization favors the occurrence of EAD formation compared to late repolarization stages (Fig. 8B). Fourthly, the window for EAD formation during the repolarization phase is less in the BZ and ICZ compared to NZ (Fig. 7B). Therefore, EADs occur for low amplitude current applied at membrane potentials from 0 to 35 mV in NZ cardiomyocytes, 0 to 25 mV in BZ cardiomyocytes and 0to20 mV in ICZ cardiomyocytes. Finally overall, ischemia reduces the propensity for EAD development (Fig. 7). This is due both to the stronger repolarization reserve in ischemia caused by I K(ATP) activation and elevated extracellular potassium concentration ([K þ ] o ). The latter decreases the excitability of the cell and therefore the likelihood of regenerative EAD-like behavior. Fig. 3. Macro-reentrant pattern of propagation in the human ventricles around the acute ischemic zone in control (top) and for 30% (middle) and 50% (bottom) I Kr reduction following premature excitation applied with coupling intervals (CIs) of 355 ms, 360 and 361 ms, respectively. The limits of the BZ are marked with green lines and the direction of propagation with white arrows. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248240 4. Discussion High-resolution simulation data obtained with an anatomicallybased model of the regionally-ischemic human ventricles reveal, for the first time, the electrotonic trigger of EADs and their contribution to transmural reentry and increased arrhythmic risk under conditions of reduced repolarization reserve. In our human ventricular model, EADs develop as a consequence of complex patterns of propagation, repolarization and current flow in the heterogeneous ischemic area, whereas they do not develop spontaneously in single cell or in homogeneous tissue with the same membrane kinetics. Simulations show that reduced I Kr causes an anti-arrhythmic prolongation of refractoriness, which counteracts macro-reentry formation. However, it simultaneously causes a proarrhythmic reduction in repolarization reserve, which promotes EAD formation and transmural reentry. Arrhythmic risk in regional ischemia is therefore determined by a fine balance between macroreentrant mechanisms driven by heterogeneity in refractoriness around the ischemic region, and micro-reentrant patterns, electrotonically-triggered EAD formation and transmural reentry. Reactivation of the L-type calcium current is the primary mechanism underlying EAD formation by prolonged low amplitude current injected during repolarization. Our investigations were primarily mechanistic, and thus conducted using a single human ventricles model given the computational demands. We identify key factors that will vary between individuals, and that may explain Fig. 4. Distribution of transmembrane voltage (V m ) throughout the ventricles at different times following ectopic excitation resulting initially in macro-reentry, but failure to support retrograde propagation (1381 ms), followed by propagation through the endocardial BZ (1775 ms) and transmurally (1831 ms), leads to intramural reentry (1831e1950 ms). For each time instant, three different views are shown, from left to right: the epicardium, the endocardium and the depolarized tissue with V m above 20 mV (cells with V m <20 mV are transparent). Isopotential lines are shown in grey linking same potential levels. More details can be found in the supplemental movies. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248 241 inter-subject differences in arrhythmic risk in acute ischemia modulated by repolarization reserve. 4.1. Computational model of the regionally-ischemic human ventricles and varied repolarization reserve For this study, a novel human model was developed, which integrates the biophysical information on membrane kinetics and the pro-arrhythmic alterations induced by acute regional ischemia, as well as fiber orientation and realistic anatomy. Simulation results are in agreement with in vivo and ex situ human experimental studies including: human ventricular activation time (Durrer et al., 1970); inverse APD-activation time relationship (Hanson et al., 2009); ischemia-induced heterogeneity in AP, CV and ERP (Figs. 1 and 2)(Carmeliet, 1999; Rodríguez et al., 2006); size of the ischemic region (Fig. 1)(Lee et al., 1981); pseudo-ECG, and its corresponding QT interval (414 ms), positive T wave, and T peak-toTend value (~100 ms) (Fig. 1)(Browne et al., 1983). The agreement between experiments and simulations lends credibility to our new findings. 4.2. Macro and micro-reentry in the regionally-ischemic human ventricles The macro-reentrant patterns of activation in the regionallyischemic human ventricles model are also consistent with those reported in pig and canine hearts in the seminal work by Janse et al. (1980) and include three necessary conditions (Fig. 3). Firstly, the presence of conduction block is necessary and occurs when heterogeneity in refractoriness is such that at the time of ectopic excitation, NZ tissue is excitable and ICZ tissue refractory (ERP in ICZ >ERP in NZ). Secondly, retrograde conduction occurs when heterogeneity in refractoriness, ischemic area size and CV are such that propagation proceeds around the ICZ, while it remains blocked in the ICZ until ischemic tissue distal to the location of ectopic excitation recovers. A prolonged ERP in ischemic tissue would require increased ischemic region size or decreased CV to allow for retrograde conduction to occur. Finally, the reentry is complete once propagation has retrogradely traversed the ICZ and the tissue in the NZ has recovered. This condition depends on the ERP in the NZ, CV in the ICZ and the size of the ischemic region. A prolonged ERP in the NZ, as with I Kr reduction, would be expected to decrease the likelihood of reentry. Our simulation results show that I Kr reduction slightly lowers the chances of meeting all three conditions for reentry by reducing the heterogeneity in ERP between the ICZ and NZ (from a range of 62e84 ms in control to 57e78 ms for 50% I Kr reduction), while prolonging the ERP in both regions, therefore reducing the chances of retrograde propagation and reentry. This leads to a decrease in the VW for 30% I Kr reduction but the trend becomes different with 50% I Kr decrease. The main reason is the emergence of a newly identified pro-arrhythmic mechanism revealed by our simulations: transmural reentry sustained by electrotonically-triggered EAD development. Even though Janse et al. also reported the occurrence of intramural reentry in some of their acute ischemia experiments in pig and dog, they were unable to resolve the underlying mechanisms due to the low resolution of their microelectrode recordings (Janse et al., 1980). Lukas and Antzelevitch also pointed at the possibility of phase 2 reentry in simulated ischemic tissue using thin Fig. 5. Electrotonically-triggered EADs are recorded near the BZ in the acutely-ischemic human ventricles with 50% I Kr reduction facilitating transmural patterns of reentry. A. Snapshots of a transmural cross-section illustrating propagation of electrical excitation from endocardium to epicardium at the times indicated. B. Time course of V m at the points indicated in panel A from endocardium to epicardium. Legend indicates the transmural location corresponding to each AP trace with % indicating transmural distance from the endocardium. Ectopic excitation is applied at t ¼1181 ms (CI ¼361 ms as in Fig. 3). S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248242 epicardial sheets in acidotic “ischemic”solution and highlighted the need for validation of their hypothesis in “future studies using high-resolution techniques”(Lukas and Antzelevitch, 1996). Our simulations indeed yield high-resolution datasets consistent with a wide range of established knowledge, and reveal for the first time the development of electrotonically-triggered EADs implicated in intramural reentry caused by the heterogeneous disease substrate in the three-dimensional volume of the human ventricles. 4.3. Mechanisms of EAD development in tissue: electrotonic kick or intrinsic instability? EADs have long been identified as a potentially important proarrhythmic player (Levy and Wiseman, 1991; Marban et al., 1986) and they have been demonstrated primarily in single cell studies with highly compromised repolarization reserve (caused by disease, mutations and/or pharmacological challenge) (Vandersickel et al., 2014; Zeng and Rudy, 1995). Despite tissue simulation studies showing their occurrence at the center of spiral wave rotors (Zemlin and Pertsov, 2007; Zemlin et al., 2009), their contribution to pro-arrhythmic tissue dynamics is still controversial, as electrotonic coupling has also been shown to suppress EAD formation and propagation (Huelsing et al., 2000; Mendez et al., 1969; Pueyo et al., 2011). Theoretical studies have demonstrated that EAD formation in tissue requires a huge reduction of repolarization reserve by simultaneously applying a 4-fold increase in I CaL and 80% reduction of I Kr ,asinVandersickel et al. (2014). This may occur in certain extreme conditions, but it does not explain increased arrhythmic risk in disease conditions caused by mildly reduced repolarization reserve as in clinical trials in patients with coronary artery disease, but without genetic disorders (Steinberg et al., 2014). Our simulations identify a more likely pro-arrhythmic mechanism of EAD formation and reentry, sustained by the Fig. 6. Intramural reentry is facilitated by prolonged APD due to EAD formation in the acutely-ischemic human ventricles with reduced repolarization reserve for long CI ¼430 ms. A. Ectopic excitation at t ¼1253 ms (corresponding to CI ¼430 ms) fails to induce reentry for 0 and 30% I Kr reduction (first and second row, respectively), but leads to establishment of transmural reentry for 50% I Kr reduction (third and fourth row). B. Time course of the action potential in the area of transmural reentry for 50% I Kr reduction marked with a blue circle in panel A. Legends indicate the transmural location for each trace (percentage of transmural distance from the endocardium). More details can be found in the supplemental figure. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248 243 heterogeneous ventricular substrate (Figs. 4e6). Reduced repolarization reserve facilitates EAD formation, but is not the primary trigger for repolarization instability, as spontaneous EADs do not occur in single cell or homogeneous tissue preparations with the same membrane kinetics. Our results show, for the first time, that EADs may be triggered in the ischemic human heart by the cell's environment in a heterogeneous substrate, which creates complex electrotonic flow of current across the BZ. Repolarization reserve and electrotonically-triggered EAD formation are therefore identified as potentially important determinants of pro-arrhythmic risk in the heterogeneous substrate of acute ischemia in human. 4.4. Implications The mechanisms highlighted in our study reveal the complex balance of electrophysiological properties that determine proFig. 7. Characterization of EADs induced by low amplitude current injected during the repolarization phase for varying current amplitudes (A) and repolarization levels (B) for human ventricular cardiomyocytes representative of NZ (top), BZ (middle) and ICZ (bottom) for 0, 30 and 50% I Kr reduction (short dash, long dash and continuous lines, respectively) (BZ: [K þ ] o ¼7 mmol/L, I K(ATP) activation ¼3%, I Na and I CaL peak conductance ¼85% and ICZ: [K þ ] o ¼8.5 mmol/L, I K(ATP) activation ¼5%, I Na and I CaL peak conductance ¼75%). In panel A, current is applied at 20 mV with amplitude 0.7, 0.9, 1.1 and 1.2 pA/pF (from left to right) and 170 ms duration. In panel B, current of 1.1 pA/pF amplitude and 170 ms duration is applied at different levels of repolarization from 0 to 40 mV transmembrane potential. The stimulation protocol consisted of a train of 100 stimuli with amplitude 30 pA/pF during 1 ms (1.3*threshold for AP trigger), at a CL of 500 ms. S. Dutta et al. / Progress in Biophysics and Molecular Biology 120 (2016) 236e248244