Regional differences in APD restitution can initiate wavebreak and re-entry in cardiac tissue: a computational study

Remodelling of potassium currents underlies arrhythmic action potential prolongation under beta-adrenergic stimulation in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) patients often present an enhanced arrhythmogenicity that can lead to lethal arrhythmias, especially during exercise. Recent studies have indicated an abnormal response of HCM cardiomyocytes to β-adrenergic receptor stimulation (β-ARS), with prolongation of their action potential rather than shortening. The mechanisms underlying this aberrant response to sympathetic stimulation and its possible proarrhythmic role remain unknown. The aims of this study are to investigate the key ionic mechanisms underlying the HCM abnormal response to β-ARS and the resultant repolarisation abnormalities using human-based experimental and computational methodologies. We integrated and calibrated the latest models of human ventricular electrophysiology and β-ARS using experimental measurements of human adult cardiomyocytes from control and HCM patients. Our major findings include: (1) the developed in silico models of β-ARS capture the behaviour observed in the experimental data, including the aberrant response of HCM cardiomyocytes to β-ARS; (2) the reduced increase of potassium currents under β-ARS was identified as the main mechanism of action potential prolongation in HCM, rather than a more sustained inward calcium current; (3) action potential duration differences between healthy and HCM cardiomyocytes were increased upon β-ARS, while endocardial to epicardial differences in HCM cardiomyocytes were reduced; (4) models presenting repolarisation abnormalities were characterised by downregulation of the rapid delayed rectifier potassium current and the sodium‑potassium pump, while inward currents were upregulated. In conclusion, our results identify causal relationships between the HCM phenotype and its arrhythmogenic response to β-ARS through the downregulation of potassium currents.

Application of two novel electrical restitution-based ECG markers of ventricular arrhythmia to patients with nonischemic cardiomyopathy

INTRODUCTION

Sudden cardiac death (SCD) risk assessment is limited, particularly in patients with nonischemic cardiomyopathies. This is the first application, in patients with cardiomyopathies, of two novel risk markers, regional restitution instability index (R2I2) and peak electrocardiogram restitution slope (PERS), which have been shown to be predictive of ventricular arrhythmias (VA) or death in ischemic heart disease patients.

METHODS

Blinded retrospective study of 50 patients: 33 dilated cardiomyopathy and 17 other; undergoing electrophysiological study (EPS) for SCD risk stratification, and 29 controls with structurally normal hearts undergoing EPS. R2I2 was calculated from an EPS using electrocardiogram surrogates for action potential duration and diastolic interval. Cut-offs for high and low R2I2/PERS were predefined.

RESULTS

R2I2 was significantly higher in study than control patients (0.99 ± 0.05 vs. 0.63 ± 0.04, p < .001). PERS showed a trend to higher values in the study group (1.18[0.63] vs. 1.09[0.54], p = .07). During median follow up of 5.6 years [interquartile range 1.9 years], nine study patients reached the endpoint of VA/death. Patients who experienced VA/death showed trends to higher mean R2I2 (1.14 ± 0.07 vs.0.95 ± 0.05, p = .12) and PERS (1.46[0.49] vs. 1.13[0.62], p = .22). A Cox proportional hazards model using grouped markers: R2I2 < 1.03 + PERS < 1.21/either R2I2 ≥ 1.03 or PERS ≥ 1.21/R2I2 ≥ 1.03 + PERS ≥ 1.21; significantly predicted VA/death (p = .02) with a hazard ratio per positive component of 3.2 (95% confidence interval 1.2-8.8).

CONCLUSION

R2I2≥ 1.03 + PERS ≥ 1.21 may predict VA/death in patients with cardiomyopathies. R2I2 ≥ 1.03 + PERS ≥ 1.21 have the potential to play an important role in SCD risk stratification in cardiomyopathies but their validity should be confirmed in a larger study.

In-silico study of the cardiac arrhythmogenic potential of biomaterial injection therapy

Biomaterial injection is a novel therapy to treat ischemic heart failure (HF) that has shown to reduce remodeling and restore cardiac function in recent preclinical studies. While the effect of biomaterial injection in reducing mechanical wall stress has been recently demonstrated, the influence of biomaterials on the electrical behavior of treated hearts has not been elucidated. In this work, we developed computational models of swine hearts to study the electrophysiological vulnerability associated with biomaterial injection therapy. The propagation of action potentials on realistic biventricular geometries was simulated by numerically solving the monodomain electrophysiology equations on anatomically-detailed models of normal, HF untreated, and HF treated hearts. Heart geometries were constructed from high-resolution magnetic resonance images (MRI) where the healthy, peri-infarcted, infarcted and gel regions were identified, and the orientation of cardiac fibers was informed from diffusion-tensor MRI. Regional restitution properties in each case were evaluated by constructing a probability density function of the action potential duration (APD) at different cycle lengths. A comparative analysis of the ventricular fibrillation (VF) dynamics for every heart was carried out by measuring the number of filaments formed after wave braking. Our results suggest that biomaterial injection therapy does not affect the regional dispersion of repolarization when comparing untreated and treated failing hearts. Further, we found that the treated failing heart is more prone to sustain VF than the normal heart, and is at least as susceptible to sustained VF as the untreated failing heart. Moreover, we show that the main features of VF dynamics in a treated failing heart are not affected by the level of electrical conductivity of the biogel injectates. This work represents a novel proof-of-concept study demonstrating the feasibility of computer simulations of the heart in understanding the arrhythmic behavior in novel therapies for HF.

Restitution metrics in Brugada syndrome: a systematic review and meta-analysis

BACKGROUND

Brugada syndrome (BrS) is an ion channelopathy that predisposes affected subjects to ventricular tachycardia/fibrillation (VT/VF) and sudden cardiac death. Restitution analysis has been examined in BrS patients but not all studies have reported significant differences between BrS patients and controls. Therefore, we conducted a systematic review and meta-analysis to investigate the different restitution indices used in BrS.

METHODS

PubMed and Embase were searched until April 7, 2019, identifying 20 and 27 studies.

RESULTS

A total of ten studies involving 178 BrS (mean age 38 years old, 63% male) and 102 controls (mean age 31 years old, 42% male) were included in this systematic review. Pacing was carried out at the right ventricular outflow tract (RVOT)/right ventricular apex (RPA) (n = 4), RPA (n = 4), or right atrium (RA) (n = 1). Basic cycle lengths of 400 (n = 4), 500 (n = 2), 600 (n = 6) and 750 ms (n = 1) were used. Recording methods include electrograms (n = 4), monophasic action potentials (n = 5), and electrocardiograms (n = 1). Signals were obtained from the RVOT (n = 8), RVA (n = 3), RA (n = 1), or the body surface (n = 1). The maximum restitution slope for endocardial repolarization at the RVOT was 0.87 for BrS patients (n = 5; 95% confidence interval [CI] 0.68-1.07) compared with 0.74 in control subjects (n = 4; 95% CI 0.42-1.06), with a significant mean difference of 0.40 (n = 4; 95% CI 0.11-0.69; P = 0.007).

CONCLUSIONS

Steeper endocardial repolarization restitution slopes are found in BrS patients compared with controls at baseline. Restitution analysis can provide important information for risk stratification in BrS.

Investigation of the relationship between two novel electrocardiogram-based sudden cardiac death risk markers and autonomic function

BACKGROUND

Regional Restitution Instability Index (R2I2) and Peak ECG Restitution Slope (PERS) are promising sudden cardiac death (SCD) risk markers. R2I2 and PERS use the standard 12‑lead ECG to measure properties of electrical restitution implicated in ventricular arrhythmogenesis. We investigated the relationship between R2I2, PERS and autonomic function to inform future application of these risk markers.

METHODS

Blinded, prospective, observational study of 44 patients with ischaemic cardiomyopathy undergoing risk stratification for an ICD. Patients underwent an electrophysiological study for determination of R2I2 and PERS. 24-hour ambulatory ECG monitoring was carried out for determination of time-domain heart rate variability (HRV).

RESULTS

During median follow up of 22 months, 11 patients experienced ventricular arrhythmia (VA)/SCD. Weak inverse correlation was seen between R2I2 and HRV-i (rho: -0.36, p = 0.02). R2I2 and PERS were significantly higher in patients experiencing VA/SCD than those not (mean ± SEM:1.14 ± 0.11 vs 0.84 ± 0.05, p = 0.01) and (1.73 ± 0.27 vs 1.07 ± 0.08, p = 0.002) respectively. Patients with low HRV-i and high PERS had an incidence rate ratio for VA/SCD 14.5 times that of patients with high HRV-i and low PERS (p = 0.02).

CONCLUSION

This small study suggests that there is minimal correlation between R2I2, PERS and autonomic function as measured by HRV. Combining PERS with HRV identified patients at particularly high risk of ventricular arrhythmia/SCD. A combined PERS+HRV risk marker may improve SCD risk stratification in patients with ischaemic cardiomyopathy.

LifeMap: towards the development of a new technology in sudden cardiac death risk stratification for clinical use

Sudden cardiac death (SCD) is a major cause of mortality presenting a significant unmet clinical need. Patients at risk of SCD are implanted with implantable cardioverter-defibrillators (ICDs) according to international guidelines based on clinical trial evidence. Implantable cardioverter-defibrillators are not inexpensive and not without problem in terms of inappropriate shocks and infection risk. Also, only a minority of patients implanted with the ICD ever use the device during its battery lifetime highlighting the fact that methods used for SCD risk stratification are inadequate. Better ways of predicting who is at risk of SCD are needed. In addition, there is no effective prevention due to the lack of understanding of the electrical mechanisms underlying SCD. Our group has been investigating the electrophysiological basis of ventricular fibrillation and have successfully applied our preclinical findings to translational studies in patients with ischaemic cardiomyopathy. We have developed two ECG markers which have been shown to be strong predictors of ventricular arrhythmias and SCD. Ongoing clinical studies are being carried out including a multicentre UK study to consolidate the evidence base. They are being incorporated into the technology, LifeMap, with the aim to develop a successful clinical tool for the assessment of SCD risk. We hereby present the scientific data leading to the technology and the development to date. The information provided here was presented at the European Heart Rhythm Association (EHRA) Europace/Cardiostim conference at which LifeMap won the EHRA Inventors Award 2016.

Restitution slope is principally determined by steady-state action potential duration

Aims

The steepness of the action potential duration (APD) restitution curve and local tissue refractoriness are both thought to play important roles in arrhythmogenesis. Despite this, there has been little recognition of the apparent association between steady-state APD and the slope of the restitution curve. The objective of this study was to test the hypothesis that restitution slope is determined by APD and to examine the relationship between restitution slope, refractoriness and susceptibility to VF.

Methods and results

Experiments were conducted in isolated hearts and ventricular myocytes from adult guinea pigs and rabbits. Restitution curves were measured under control conditions and following intervention to prolong (clofilium, veratridine, bretylium, low [Ca]e, chronic transverse aortic constriction) or shorten (catecholamines, rapid pacing) ventricular APD. Despite markedly differing mechanisms of action, all interventions that prolonged the action potential led to a steepening of the restitution curve (and vice versa). Normalizing the restitution curve as a % of steady-state APD abolished the difference in restitution curves with all interventions. Effects on restitution were preserved when APD was modulated by current injection in myocytes pre-treated with the calcium chelator BAPTA-AM – to abolish the intracellular calcium transient. The non-linear relation between APD and the rate of repolarization of the action potential is shown to underpin the common influence of APD on the slope of the restitution curve. Susceptibility to VF was found to parallel changes in APD/refractoriness, rather than restitution slope.

Conclusion(s)

Steady-state APD is the principal determinant of the slope of the ventricular electrical restitution curve. In the absence of post-repolarization refractoriness, factors that prolong the action potential would be expected to steepen the restitution curve. However, concomitant changes in tissue refractoriness act to reduce susceptibility to sustained VF. Dependence on steady-state APD may contribute to the failure of restitution slope to predict sudden cardiac death.

Arrhythmogenic drugs can amplify spatial heterogeneities in the electrical restitution in perfused guinea-pig heart: An evidence from assessments of monophasic action potential durations and JT intervals

Non-uniform shortening of the action potential duration (APD90) in different myocardial regions upon heart rate acceleration can set abnormal repolarization gradients and promote arrhythmia. This study examined whether spatial heterogeneities in APD90 restitution can be amplified by drugs with clinically proved proarrhythmic potential (dofetilide, quinidine, procainamide, and flecainide) and, if so, whether these effects can translate to the appropriate changes of the ECG metrics of ventricular repolarization, such as JT intervals. In isolated, perfused guinea-pig heart preparations, monophasic action potentials and volume-conducted ECG were recorded at progressively increased pacing rates. The APD90 measured at distinct ventricular sites, as well as the JTpeak and JTend values were plotted as a function of preceding diastolic interval, and the maximum slopes of the restitution curves were determined at baseline and upon drug administration. Dofetilide, quinidine, and procainamide reverse rate-dependently prolonged APD90 and steepened the restitution curve, with effects being greater at the endocardium than epicardium, and in the right ventricular (RV) vs. the left ventricular (LV) chamber. The restitution slope was increased to a greater extent for the JTend vs. the JTpeak interval. In contrast, flecainide reduced the APD90 restitution slope at LV epicardium without producing effect at LV endocardium and RV epicardium, and reduced the JTpeak restitution slope without changing the JTend restitution. Nevertheless, with all agents, these effects translated to the amplified epicardial-to-endocardial and the LV-to-RV non-uniformities in APD90 restitution, paralleled by the increased JTend vs. JTpeak difference in the restitution slope. In summary, these findings suggest that arrhythmic drug profiles are partly attributable to the accentuated regional heterogeneities in APD90 restitution, which can be indirectly determined through ECG assessments of the JTend vs. JTpeak dynamics at variable pacing rates.

Discrepant effects of heart failure on electrophysiological property in right ventricular outflow tract and left ventricular outflow tract cardiomyocytes

Ventricular arrhythmias commonly arise from the right (RVOT) and left ventricular outflow tracts (LVOT) in patients without structural heart disease. Heart failure (HF) significantly increases the risk of ventricular arrhythmias. The regional differences and how HF affects the electrophysiological characteristics of RVOT and LVOT cardiomyocytes remain unclear. The whole-cell patch-clamp technique was used to investigate the action potentials and ionic currents in isolated single RVOT and LVOT cardiomyocytes from control rabbits and rabbits with HF induced by rapid ventricular pacing. Comparison with control LVOT cardiomyocytes showed that control RVOT cardiomyocytes have a shorter action potential duration (APD), smaller late Na+ currents (INa-late), larger transient outward (Ito) and larger delayed rectifier K+ currents (IKr-tail), but had similar L-type Ca2+ currents (ICa-L) and Na+/Ca2+ exchanger (NCX) current. HF increased APD, INa-late and NCX, but decreased ICa-L and Ito in RVOT cardiomyocytes. In contrast with this, HF decreased APD and ICa-L, but increased Ito and IKr-tail in LVOT cardiomyocytes. In conclusion, RVOT and LVOT cardiomyocytes had distinctive electrophysiological characteristics. HF differentially modulates action potential morphology and ionic currents in RVOT and LVOT cardiomyocytes.

Beat-to-beat ECG restitution: A review and proposal for a new biomarker to assess cardiac stress and ventricular tachyarrhythmia vulnerability

BACKGROUND

Cardiac restitution is the ability of the heart to recover from one beat to the next. Ventricular arrhythmia vulnerability can occur when the heart does not properly adjust to sudden changes in rate or in hemodynamics leading to excessive temporal and/or spatial heterogeneity in conduction or repolarization. Restitution has historically been used to study, by invasive means, the dynamics of the relationship between action potential duration (APD) and diastolic interval (DI) in sedated subjects using various pacing protocols. Even though the analogous measures of APD and DI can be obtained using the surface ECG to acquire the respective QT and TQ intervals for ECG restitution, this methodology has not been widely adopted for a number of reasons.

METHODS

Recent development of more advanced software algorithms enables ECG intervals to be measured accurately, on a continuous beat-to-beat basis, in an automated manner, and under highly dynamic conditions (i.e., ambulatory or exercise) providing information beyond that available in the typical resting state.

RESULTS

Current breakthroughs in ECG technology will allow ECG restitution measures to become a practical approach for providing quantitative measures of the risks for ventricular arrhythmias as well as cardiac stress in general.

CONCLUSIONS

In addition to a review of the underlying principles and caveats of ECG restitution, a new approach toward an advancement of more integrated restitution biomarkers is proposed.

The control of cardiac ventricular excitability by autonomic pathways

Central to the genesis of ventricular cardiac arrhythmia are variations in determinants of excitability. These involve individual ionic channels and transporters in cardiac myocytes but also tissue factors such as variable conduction of the excitation wave, fibrosis and source-sink mismatch. It is also known that in certain diseases and particularly the channelopathies critical events occur with specific stressors. For example, in hereditary long QT syndrome due to mutations in KCNQ1 arrhythmic episodes are provoked by exercise and in particular swimming. Thus not only is the static substrate important but also how this is modified by dynamic signalling events associated with common physiological responses. In this review, we examine the regulation of ventricular excitability by signalling pathways from a cellular and tissue perspective in an effort to identify key processes, effectors and potential therapeutic approaches. We specifically focus on the autonomic nervous system and related signalling pathways.

EFFECTS OF ACUTE GLOBAL ISCHEMIA ON RE-ENTRANT ARRHYTHMOGENESIS: A SIMULATION STUDY

Sudden cardiac death is mainly caused by arrhythmogenesis. For a functional abnormal heart, such as an ischemic heart, the probability of arrhythmia occurring is greatly increased. During myocardial ischemia, re-entry is prone to degenerate into ventricular fibrillation (VF). Therefore it has important meaning to investigate the intricate mechanisms underlying VF under an ischemic condition in order to better facilitate therapeutic interventions. In this paper, to analyze the functional influence of acute global ischemia on cardiac electrical activity and subsequently on re-entrant arrhythmogenesis, we take into account three main pathophysiological consequences of ischemia: hyperkalaemia, acidosis, and anoxia, and develop a 3D human ventricular ischemic model that combines a detailed biophysical description of the excitation kinetics of human ventricular cells with an integrated geometry of human ventricular tissue which incorporates fiber direction anisotropy and the stimulation activation sequence. The results show that under acute global ischemia, the tissue excitability and the slope of ventricular cellular action potential duration restitution (APDR) are greatly decreased. As a result, the complexity of VF activation patterns is reduced. For the three components of ischemia, hyperkalaemia is the dominant contributor to the stability of re-entry under acute global ischemia. Increasing [K+]o acts to prolong the cell refractory period, reduce the tissue excitability and slow the conduction velocity. Our results also show that VF can be eliminated by decreasing cellular excitability, primarily by elevating the concentration value of extracellular K+.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 2D simulation study

BACKGROUND

Heart failure is operationally defined as the inability of the heart to maintain blood flow to meet the needs of the body and it is the final common pathway of various cardiac pathologies. Electrophysiological remodeling, intercellular uncoupling and a pro-fibrotic response have been identified as major arrhythmogenic factors in heart failure.

OBJECTIVE

In this study we investigate vulnerability to reentry under heart failure conditions by incorporating established electrophysiological and anatomical remodeling using computer simulations.

METHODS

The electrical activity of human transmural ventricular tissue (5 cm × 5 cm) was simulated using the human ventricular action potential model Grandi et al. under control and heart failure conditions. The MacCannell et al. model was used to model fibroblast electrical activity, and their electrotonic interactions with myocytes. Selected degrees of diffuse fibrosis and variations in intercellular coupling were considered and the vulnerable window (VW) for reentry was evaluated following cross-field stimulation.

RESULTS

No reentry was observed in normal conditions or in the presence of HF ionic remodeling. However, defined amount of fibrosis and/or cellular uncoupling were sufficient to elicit reentrant activity. Under conditions where reentry was generated, HF electrophysiological remodeling did not alter the width of the VW. However, intermediate fibrosis and cellular uncoupling significantly widened the VW. In addition, biphasic behavior was observed, as very high fibrotic content or very low tissue conductivity hampered the development of reentry. Detailed phase analysis of reentry dynamics revealed an increase of phase singularities with progressive fibrotic components.

CONCLUSION

Structural remodeling is a key factor in the genesis of vulnerability to reentry. A range of intermediate levels of fibrosis and intercellular uncoupling can combine to favor reentrant activity.

In vivo human left-to-right ventricular differences in rate adaptation transiently increase pro-arrhythmic risk following rate acceleration

Left-to-right ventricular (LV/RV) differences in repolarization have been implicated in lethal arrhythmias in animal models. Our goal is to quantify LV/RV differences in action potential duration (APD) and APD rate adaptation and their contribution to arrhythmogenic substrates in the in vivo human heart using combined in vivo and in silico studies. Electrograms were acquired from 10 LV and 10 RV endocardial sites in 15 patients with normal ventricles. APD and APD adaptation were measured during an increase in heart rate. Analysis of in vivo electrograms revealed longer APD in LV than RV (207.8 ± 21.5 vs 196.7 ± 20.1 ms; P<0.05), and slower APD adaptation in LV than RV (time constant τ(s) =47.0 ± 14.3 vs 35.6 ± 6.5 s; P<0.05). Following rate acceleration, LV/RV APD dispersion experienced an increase of up to 91% in 12 patients, showing a strong correlation (r(2) =0.90) with both initial dispersion and LV/RV difference in slow adaptation. Pro-arrhythmic implications of measured LV/RV functional differences were studied using in silico simulations. Results show that LV/RV APD and APD adaptation heterogeneities promote unidirectional block following rate acceleration, albeit being insufficient for establishment of reentry in normal hearts. However, in the presence of an ischemic region at the LV/RV junction, LV/RV heterogeneity in APD and APD rate adaptation promotes reentrant activity and its degeneration into fibrillatory activity. Our results suggest that LV/RV heterogeneities in APD adaptation cause a transient increase in APD dispersion in the human ventricles following rate acceleration, which promotes unidirectional block and wave-break at the LV/RV junction, and may potentiate the arrhythmogenic substrate, particularly in patients with ischemic heart disease.

A novel surface electrocardiogram-based marker of ventricular arrhythmia risk in patients with ischemic cardiomyopathy

Background

Better sudden cardiac death risk markers are needed in ischemic cardiomyopathy (ICM). Increased heterogeneity of electrical restitution is an important mechanism underlying the risk of ventricular arrhythmia (VA). Our aim was to develop and test a novel quantitative surface electrocardiogram–based measure of VA risk in patients with ICM: the Regional Restitution Instability Index (R2I2).

Methods and Results

R2I2, the mean of the standard deviation of residuals from the mean gradient for each ECG lead at a range of diastolic intervals, was measured retrospectively from high-resolution 12-lead ECGs recorded during an electrophysiology study. Patient groups were as follows: Study group, 26 patients with ICM being assessed for implantable defibrillator; Control group, 29 patients with supraventricular tachycardia undergoing electrophysiology study; and Replication group, 40 further patients with ICM. R2I2 was significantly higher in the Study patients than in Controls (mean ± standard error of the mean: 1.09±0.06 versus 0.63±0.04, P<0.001). Over a median follow-up period of 23 months, 6 of 26 Study group patients had VA or death. R2I2 predicted VA or death independently of demographic factors, electrophysiology study result, left ventricular ejection fraction, or QRS duration (Cox model, P=0.029). R2I2 correlated with peri-infarct zone as assessed by cardiac magnetic resonance imaging (r=0.51, P=0.024). The findings were replicated in the Replication group: R2I2 was significantly higher in 11 of 40 Replication patients experiencing VA (1.18±0.10 versus 0.92±0.05, P=0.019). In combined analysis of ICM cohorts, R2I2 ≥1.03 identified subjects with significantly higher risk of VA or death (43%) compared with R2I2 <1.03 (11%) (P=0.004).

Conclusions

In this pilot study, we have developed a novel VA risk marker, R2I2, and have shown that it correlated with a structural measure of arrhythmic risk and predicted risk of VA or death in patients with ICM. R2I2 may improve risk stratification and merits further evaluation. (J Am Heart Assoc. 2012;1:e001552 doi: 10.1161/JAHA.112.001552.)

Mechanisms of ventricular arrhythmias: a dynamical systems-based perspective

Defining the cellular electrophysiological mechanisms for ventricular tachyarrhythmias is difficult, given the wide array of potential mechanisms, ranging from abnormal automaticity to various types of reentry and kk activity. The degree of difficulty is increased further by the fact that any particular mechanism may be influenced by the evolving ionic and anatomic environments associated with many forms of heart disease. Consequently, static measures of a single electrophysiological characteristic are unlikely to be useful in establishing mechanisms. Rather, the dynamics of the electrophysiological triggers and substrates that predispose to arrhythmia development need to be considered. Moreover, the dynamics need to be considered in the context of a system, one that displays certain predictable behaviors, but also one that may contain seemingly stochastic elements. It also is essential to recognize that even the predictable behaviors of this complex nonlinear system are subject to small changes in the state of the system at any given time. Here we briefly review some of the short-, medium-, and long-term alterations of the electrophysiological substrate that accompany myocardial disease and their potential impact on the initiation and maintenance of ventricular arrhythmias. We also provide examples of cases in which small changes in the electrophysiological substrate can result in rather large differences in arrhythmia outcome. These results suggest that an interrogation of cardiac electrical dynamics is required to provide a meaningful assessment of the immediate risk for arrhythmia development and for evaluating the effects of putative antiarrhythmic interventions.

Teaching cardiac electrophysiology modeling to undergraduate students: laboratory exercises and GPU programming for the study of arrhythmias and spiral wave dynamics

As part of a 3-wk intersession workshop funded by a National Science Foundation Expeditions in Computing award, 15 undergraduate students from the City University of New York(1) collaborated on a study aimed at characterizing the voltage dynamics and arrhythmogenic behavior of cardiac cells for a broad range of physiologically relevant conditions using an in silico model. The primary goal of the workshop was to cultivate student interest in computational modeling and analysis of complex systems by introducing them through lectures and laboratory activities to current research in cardiac modeling and by engaging them in a hands-on research experience. The success of the workshop lay in the exposure of the students to active researchers and experts in their fields, the use of hands-on activities to communicate important concepts, active engagement of the students in research, and explanations of the significance of results as the students generated them. The workshop content addressed how spiral waves of electrical activity are initiated in the heart and how different parameter values affect the dynamics of these reentrant waves. Spiral waves are clinically associated with tachycardia, when the waves remain stable, and with fibrillation, when the waves exhibit breakup. All in silico experiments were conducted by simulating a mathematical model of cardiac cells on graphics processing units instead of the standard central processing units of desktop computers. This approach decreased the run time for each simulation to almost real time, thereby allowing the students to quickly analyze and characterize the simulated arrhythmias. Results from these simulations, as well as some of the background and methodology taught during the workshop, is presented in this article along with the programming code and the explanations of simulation results in an effort to allow other teachers and students to perform their own demonstrations, simulations, and studies.

Oscillatory patterns of respiration: consequences for the stability and control of cardiac electrophysiology

Periodic breathing patterns known as Central Sleep Apnea (CSA) are often observed in congestive heart failure. This phenomenon is associated with increased risk of sudden cardiac death, but the mechanism for that outcome has not been exposed. Endocardial electrograms were recorded during spontaneous episodes of CSR and PB in patients in conscious and unconscious states. Analysis exposed a regular bidirectional phase-walk in the relationship between respiration and arterial blood pressure. Recently developed signal processing techniques revealed that respiration also modulates cardiac repolarization properties at multiple simultaneous frequencies, and the effect was heterogeneous across measurement sites in both ventricles. These measurements offer unique evidence of the electro-physiological manifestations of these breathing patterns. Analysis of phase relationships suggested a mechanism by which the behavior may predispose patients to cardiac arrhythmias. Such predisposition would be easily measured to direct treatment priorities and improve risk stratification.

Heterogeneous memory in restitution of action potential duration in pig ventricles

BACKGROUND

Restitution of action potential duration and memory importantly affect electrical stability in ventricles. Studies have reported heterogeneous restitution among different regions of the ventricles. However, existence of heterogeneity in memory is not as well investigated.

METHODS

Transmembrane potentials were recorded in endocardial and epicardial tissues from both ventricles of farm pigs. Pacing protocols with sinusoidally changing diastolic intervals were used to reveal hysteresis in restitution, from which quantitative measures of memory were calculated.

RESULTS

Larger measures of hysteresis were observed in the endocardium than the epicardium (P < .05): loop thickness (in milliseconds), 26.9 vs 16.2; overall tilt, 0.376 vs 0.249; and loop area (in square milliseconds), 7288 vs 4146. Except for overall tilt, no significant differences in these measures were observed between ventricles.

CONCLUSION

Heterogeneity in memory exists in pig ventricles. Because regions with the steepest restitution may also have the largest memory, our results suggest that heterogeneity in memory should also be factored in when predicting electrical stability.

Interaction of Activation–Repolarization Coupling and Restitution Properties in Humans

Background— Dynamic modulation of repolarization is important in arrhythmogenesis. An inverse relation exists in myocardium between activation time (AT) and action potential duration (APD). We hypothesized that resulting gradients of APD and diastolic interval (DI) interact with restitution properties and modulate the timing of repolarization.

Methods and Results— Activation-recovery intervals (ARI) were acquired from reconstructed noncontact unipolar electrograms from the left ventricular endocardium in 9 patients (7 male) with normal ventricles. At a basic paced cycle length (median, 450 ms), ARIs shortened along the path of activation, with a mandatory reciprocal increase of DIs. In the median patient, this range of DIs started at 230 ms at the site of earliest activation and increased to 279 ms at the site of latest activation at a basic cycle length of 450 ms. Four consecutive standard S1 to S2 restitution curves were performed. At sites with a longer ARI (and therefore shorter DI) close to the site of stimulation, premature stimulation produced more shortening of ARIs; therefore, the time course of restitution was steeper than at more distal sites. At normal heart rate, the decrease in ARIs along the conduction pathway compensated for later activation. Thus, dispersion in repolarization time (RT) is smaller than dispersion in ARI in a heart with a steep negative AT–ARI relationship. This protective effect is lost in hearts without such a relationship. In the patients with a steep AT–ARI relationship at basic cycle length, this relation is lost after premature stimulation and is a function of prematurity. Thus, dispersion in RT is larger after shortly coupled extra stimuli in patients with a steep AT–ARI relationship.

Conclusions— A complex interplay exists between activation–repolarization coupling and restitution properties, largely driven by ARI and DI gradients. This plays a significant role in the dynamics of repolarization in humans.

Termination of equine atrial fibrillation by quinidine: an optical mapping study

OBJECTIVE

To perform the first optical mapping studies of equine atrium to assess the spatiotemporal dynamics of atrial fibrillation (AF) and of its termination by quinidine.

ANIMALS

Intact, perfused atrial preparations obtained from four horses with normal cardiovascular examinations.

MATERIALS AND METHODS

AF was induced by a rapid pacing protocol with or without acetylcholine perfusion, and optical mapping was used to determine spatial dominant frequency distributions, electrical activity maps, and single-pixel optical signals. Following induction of AF, quinidine gluconate was perfused into the preparation and these parameters were monitored during quinidine-induced termination of AF.

RESULTS

Equine AF develops in the context of spatial gradients in action potential duration (APD) and diastolic interval (DI) that produce alternans, conduction block, and Wenckebach conduction in different regions at fast pacing rates. Quinidine terminates AF and prevents subsequent reinduction by reducing the maximal frequency and increasing frequency homogeneity.

CONCLUSIONS

Heterogeneity of APD and DI promote alternans and conduction block at fast pacing rates in the equine atrium, predisposing to the development of AF. Quinidine terminates AF by reducing maximum frequency and increasing frequency homogeneity. Our results are consistent with the hypothesis that quinidine increases effective refractory period, thereby decreasing frequency.

Formation of fast spirals on heterogeneities of an excitable medium

We study the process of formation of spiral waves in a heterogeneous excitable medium under external stimulation, using numerical and analytical methods. We show that in an excitable medium with several heterogeneities with respect to refractory period, fast rotating spiral waves can be generated. These fast spirals are formed as a result of a phenomenon of period decrease, which is the generation by a heterogeneity of waves with a period shorter than the period of the external stimulation.

Cardiac mechano-electric feedback and electrical restitution in humans

Electrical restitution in the heart is the property whereby the action potential duration and conduction velocity of a beat of altered cycle length vary according to its immediacy to the preceding basic beat--the coupling interval, usually the diastolic interval. In general, action potential duration (APD) increases with increasing coupling interval, and the relation between action potential duration and the preceding diastolic interval describes the APD restitution curve. The latter has recently been the focus of considerable interest since the steepness of the initial part of the restitution curve plays an important role in electrical stability and arrhythmogenesis. Mechanical stretch has been shown to alter APD and hence refractoriness either through stretch activated channels or by influencing calcium cycling. Such an effect on refractoriness has been proposed as a mechanism of arrhythmogenesis particularly if spatially inhomogeneities manifest within the heart. Here, we review (1) the spatial and temporal characteristics of APD restitution in humans; (2) previously reported work showing that mechanical loading differentially effects APD of interpolated beats of altered cycle length, and hence alters the slope of the APD restitution curve; and (3) evidence that inhomogeneity of APD restitution slope may be an important factor in arrhythmogenesis.

Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence

Ventricular fibrillation is a lethal arrhythmia characterized by multiple wavelets usually starting from a single or figure-of-eight re-entrant circuit. Understanding the factors regulating vulnerability to the re-entry is essential for developing effective therapeutic strategies to prevent ventricular fibrillation. In this study, we investigated how pre-existing tissue heterogeneities and electrical restitution properties affect the initiation of re-entry by premature extrastimuli in two-dimensional cardiac tissue models. We studied two pacing protocols for inducing re-entry following the "sinus" rhythm (S1) beat: (1) a single premature (S2) extrastimulus in heterogeneous tissue; (2) two premature extrastimuli (S2 and S3) in homogeneous tissue. In the first case, the vulnerable window of re-entry is determined by the spatial dimension and extent of the heterogeneity, and is also affected by electrical restitution properties and the location of the premature stimulus. The vulnerable window first increases as the action potential duration (APD) difference between the inside and outside of the heterogeneous region increases, but then decreases as this difference increases further. Steeper APD restitution reduces the vulnerable window of re-entry. In the second case, electrical restitution plays an essential role. When APD restitution is flat, no re-entry can be induced. When APD restitution is steep, re-entry can be induced by an S3 over a range of S1S2 intervals, which is also affected by conduction velocity restitution. When APD restitution is even steeper, the vulnerable window is reduced due to collision of the spiral tips.

A guide to modelling cardiac electrical activity in anatomically detailed ventricles

One of the most recent trends in cardiac electrophysiology is the development of integrative anatomically accurate models of the heart, which include description of cardiac activity from sub-cellular and cellular level to the level of the whole organ. In order to construct this type of model, a researcher needs to collect a wide range of information from books and journal articles on various aspects of biology, physiology, electrophysiology, numerical mathematics and computer programming. The aim of this methodological article is to survey recent developments in integrative modelling of electrical activity in the ventricles of the heart, and to provide a practical guide to the resources and tools that are available for work in this exciting and challenging area.

Steeper restitution slopes across right ventricular endocardium in patients with cardiomyopathy at high risk of ventricular arrhythmias

Steep action potential duration (APD) restitution slopes (>1) and spatial APD restitution heterogeneity provide the substrate for ventricular fibrillation in computational models and experimental studies. Their relationship to ventricular arrhythmia vulnerability in human cardiomyopathy has not been defined. Patients with cardiomyopathy [left ventricular (LV) ejection fraction <40%] and no history of ventricular arrhythmias underwent risk stratification with programmed electrical stimulation or T wave alternans (TWA). Low-risk patients ( n = 10) had no inducible ventricular tachycardia (VT) or negative TWA, while high-risk patients ( n = 8) had inducible VT or positive TWA. Activation recovery interval (ARI) restitution slopes were measured simultaneously from 10 right ventricular (RV) endocardial sites during an S1-S2 pacing protocol. ARI restitution slope heterogeneity was defined as the coefficient of variation of slopes. Mean ARI restitution slope was significantly steeper in the high-risk group compared with the low-risk group [1.16 (SD 0.31) vs. 0.59 (SD 0.19), P = 0.0002]. The proportion of endocardial recording sites with a slope >1 was significantly larger in the high-risk patients [47% (SD 35) vs. 13% (SD 21), P = 0.022]. Spatial heterogeneity of ARI restitution slopes was similar between the two groups [29% (SD 16) vs. 39% (SD 34), P = 0.48]. There was an inverse linear relationship between the ARI restitution slope and the minimum diastolic interval ( P < 0.001). In cardiomyopathic patients at high risk of ventricular arrhythmias, ARI restitution slopes along the RV endocardium are steeper, but restitution slope heterogeneity is similar compared with those at low risk. Steeper ARI restitution slopes may increase the propensity for ventricular arrhythmias in patients with impaired left ventricular function.

Theory of action potential wave block at-a-distance in the heart

Propagation failure of an action potential wave at a finite distance from its source (so-called type-II block) may cause spiral wave formation or wave breakup in the heart, phenomena that are believed to underlie lethal and nonlethal heart rhythm disorders. In this study, we develop a sufficient condition for this type of block in a homogeneous, spatially one-dimensional system. Using a topological argument, we find that type-II block of a wave will always occur when launched within a finite range of times if the velocity of the trailing edge of the preceding wave, as measured at the stimulus site, is smaller than the velocity of a wave launched with the minimum diastolic interval (DI) for which propagation is possible. This "blocking condition" is robust, remaining valid even when memory and waveback electrotonic effects are included. The condition suggests that type-II block is greatly facilitated when waves are initiated at irregular intervals in time such that (1) the velocities of consecutive waves are as different as possible and (2) the DIs preceding each wave fall on the steeply sloped portion of the action potential duration restitution curve as often as possible. The set of timing intervals between stimuli that are predicted by the blocking condition to produce block are found to be consistent with these guidelines, and also to agree well with a coupled-maps computer simulation model, for the case of waves launched by four rapidly and irregularly timed stimuli.

The missing link between cardiovascular rhythm control and myocardial cell modeling

Cardiac arrhythmia is currently investigated from two different points of view. One considers ECG bio-signal analysis and investigates heart rate variability, baroreflex control, heart rate turbulence, alternans phenomena, etc. The other involves building computer models of the heart based on ion channels, bio-domain models and forward calculations to finally reach ECG and body surface potential maps. Both approaches aim to support the cardiologist in better understanding of arrhythmia, improving diagnosis and reliable risk stratification, and optimizing therapy. This article summarizes recent results and aims to trigger new research to bridge the different views.

Whole heart action potential duration restitution properties in cardiac patients: a combined clinical and modelling study

Steep action potential duration (APD) restitution has been shown to facilitate wavebreak and ventricular fibrillation. The global APD restitution properties in cardiac patients are unknown. We report a combined clinical electrophysiology and computer modelling study to: (1) determine global APD restitution properties in cardiac patients; and (2) examine the interaction of the observed APD restitution with known arrhythmia mechanisms. In 14 patients aged 52-85 years undergoing routine cardiac surgery, 256 electrode epicardial mapping was performed. Activation-recovery intervals (ARI; a surrogate for APD) were recorded over the entire ventricular surface. Mono-exponential restitution curves were constructed for each electrode site using a standard S1-S2 pacing protocol. The median maximum restitution slope was 0.91, with 27% of all electrode sites with slopes<0.5, 29% between 0.5 and 1.0, and 20% between 1.0 and 1.5. Eleven per cent of restitution curves maintained slope>1 over a range of diastolic intervals of at least 30 ms; and 0.3% for at least 50 ms. Activation-recovery interval restitution was spatially heterogeneous, showing regional organization with multiple discrete areas of steep and shallow slope. We used a simplified computer model of 2-D cardiac tissue to investigate how heterogeneous APD restitution can influence vulnerability to, and stability of re-entry. Our model showed that heterogeneity of restitution can act as a potent arrhythmogenic substrate, as well as influencing the stability of re-entrant arrhythmias. Global epicardial mapping in humans showed that APD restitution slopes were organized into regions of shallow and steep slopes. This heterogeneous organization of restitution may provide a substrate for arrhythmia.