The role of short term memory and conduction velocity restitution in alternans formation

Fast creation of data-driven low-order predictive cardiac tissue excitation models from recorded activation patterns

Excitable systems give rise to important phenomena such as heat waves, epidemics and cardiac arrhythmias. Understanding, forecasting and controlling such systems requires reliable mathematical representations. For cardiac tissue, computational models are commonly generated in a reaction-diffusion framework based on detailed measurements of ionic currents in dedicated single-cell experiments. Here, we show that recorded movies at the tissue-level of stochastic pacing in a single variable are sufficient to generate a mathematical model. Via exponentially weighed moving averages, we create additional state variables, and use simple polynomial regression in the augmented state space to quantify excitation wave dynamics. A spatial gradient-sensing term replaces the classical diffusion as it is more robust to noise. Our pipeline for model creation is demonstrated for an in-silico model and optical voltage mapping recordings of cultured human atrial myocytes and only takes a few minutes. Our findings have the potential for widespread generation, use and on-the-fly refinement of personalised computer models for non-linear phenomena in biology and medicine, such as predictive cardiac digital twins.

Mechanistic insights into spontaneous transition from cellular alternans to ventricular fibrillation

T-wave alternans (TWA) has been used for predicting the risk of malignant cardiac arrhythmias and sudden cardiac death (SCD) in multiple clinical settings; however, possible mechanism(s) underlying the spontaneous transition from cellular alternans reflected by TWA to arrhythmias in impaired repolarization remains unclear. The healthy guinea pig ventricular myocytes under E-4031 blocking IKr (0.1 μM, N = 12; 0.3 μM, N = 10; 1 μM, N = 10) were evaluated using whole-cell patch-clamp. The electrophysiological properties of isolated perfused guinea pig hearts under E-4031 (0.1 μM, N = 5; 0.3 μM, N = 5; 1 μM, N = 5) were evaluated using dual- optical mapping. The amplitude/threshold/restitution curves of action potential duration (APD) alternans and potential mechanism(s) underlying the spontaneous transition of cellular alternans to ventricular fibrillation (VF) were examined. There were longer APD80 and increased amplitude and threshold of APD alternans in E-4031 group compared with baseline group, which was reflected by more pronounced arrhythmogenesis at the tissue level, and were associated with steep restitution curves of the APD and the conduction velocity (CV). Conduction of AP alternans augmented tissue's functional spatiotemporal heterogeneity of regional AP/Ca alternans, as well as the AP/Ca dispersion, leading to localized uni-directional conduction block that spontaneous facilitated the formation of reentrant excitation waves without the need for additional premature stimulus. Our results provide a possible mechanism for the spontaneous transition from cardiac electrical alternans in cellular action potentials and intercellular conduction without the involvement of premature excitations, and explain the increased susceptibility to ventricular arrhythmias in impaired repolarization. In this study, we implemented voltage-clamp and dual-optical mapping approaches to investigate the underlying mechanism(s) for the arrhythmogenesis of cardiac alternans in the guinea pig heart at cellular and tissue levels. Our results demonstrated a spontaneous development of reentry from cellular alternans, arising from a combined actions of restitution properties of action potential duration, conduction velocity of excitation wave and interplay between alternants of action potential and the intracellular Ca handling. We believe this study provides new insights into underlying the mechanism, by which cellular cardiac alternans spontaneously evolves into cardiac arrhythmias.

Electrophysiological Mechanisms Underlying T-Wave Alternans and Their Role in Arrhythmogenesis

T-wave alternans (TWA) reflects every-other-beat alterations in the morphology of the electrocardiogram ST segment or T wave in the setting of a constant heart rate, hence, in the absence of heart rate variability. It is believed to be associated with the dispersion of repolarization and has been used as a non-invasive marker for predicting the risk of malignant cardiac arrhythmias and sudden cardiac death as numerous studies have shown. This review aims to provide up-to-date review on both experimental and simulation studies in elucidating possible mechanisms underlying the genesis of TWA at the cellular level, as well as the genesis of spatially concordant/discordant alternans at the tissue level, and their transition to cardiac arrhythmia. Recent progress and future perspectives in antiarrhythmic therapies associated with TWA are also discussed.

Dual regulation by subcellular calcium heterogeneity and heart rate variability on cardiac electromechanical dynamics

Heart rate constantly varies under physiological conditions, termed heart rate variability (HRV), and in clinical studies, low HRV is associated with a greater risk of cardiac arrhythmias. Prior work has shown that HRV influences the temporal patterns of electrical activity, specifically the formation of pro-arrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca) levels. We previously showed that HRV may be anti-arrhythmic by disrupting APD and Ca alternations in a homogeneous cardiac myocyte. Here, we expand on our previous work, incorporating variation in subcellular Ca handling (also known to influence alternans) into a nonlinear map model of a cardiac myocyte composed of diffusively coupled Ca release units (CRUs). Ca-related parameters and initial conditions of each CRU are varied to mimic subcellular Ca heterogeneity, and a stochastic pacing sequence reproduces HRV. We find that subcellular Ca heterogeneity promotes the formation of spatially discordant subcellular alternans patterns, which decreases whole cell Ca and APD alternation for low and moderate HRV, while high subcellular Ca heterogeneity and HRV both promote electromechanical desynchronization. Finally, we find that for low and moderate HRV, both the specific subcellular Ca-related parameters and the pacing sequences influence measures of electromechanical dynamics, while for high HRV, these measures depend predominantly on the pacing sequence. Our results suggest that pro-arrhythmic subcellular discordant alternans tend to form for low levels of HRV, while high HRV may be anti-arrhythmic due to mitigated influence from subcellular Ca heterogeneity and desynchronization of APD from Ca instabilities.

Control of voltage-driven instabilities in cardiac myocytes with memory

Sudden cardiac death is known to be associated with dynamical instabilities in the heart, and thus control of dynamical instabilities is considered a potential therapeutic strategy. Different control methods were developed previously, including time-delayed feedback pacing control and constant diastolic interval pacing control. Experimental, theoretical, and simulation studies have examined the efficacy of these control methods in stabilizing action potential dynamics. In this study, we apply these control methods to control complex action potential (AP) dynamics under two diseased conditions: early repolarization syndrome and long QT syndrome, in which voltage-driven instabilities occur in the presence of short-term cardiac memory. In addition, we also develop a feedback pacing method to stabilize these instabilities. We perform theoretical analyses using iterated map models and carry out numerical simulations of AP models. We show that under the normal condition where the memory effect is minimal, all three methods can effectively control the action potential duration (APD) dynamics. Under the two diseased conditions where the memory effect is exacerbated, constant diastolic pacing control is least effective, while the feedback pacing control is most effective. Under a very strong memory effect, all three methods fail to stabilize the voltage-driven instabilities. The failure of effective control is due to memory and the all-or-none AP dynamics which results in very steep changes in APD.

Memory-induced nonlinear dynamics of excitation in cardiac diseases

Excitable cells, such as cardiac myocytes, exhibit short-term memory, i.e., the state of the cell depends on its history of excitation. Memory can originate from slow recovery of membrane ion channels or from accumulation of intracellular ion concentrations, such as calcium ion or sodium ion concentration accumulation. Here we examine the effects of memory on excitation dynamics in cardiac myocytes under two diseased conditions, early repolarization and reduced repolarization reserve, each with memory from two different sources: slow recovery of a potassium ion channel and slow accumulation of the intracellular calcium ion concentration. We first carry out computer simulations of action potential models described by differential equations to demonstrate complex excitation dynamics, such as chaos. We then develop iterated map models that incorporate memory, which accurately capture the complex excitation dynamics and bifurcations of the action potential models. Finally, we carry out theoretical analyses of the iterated map models to reveal the underlying mechanisms of memory-induced nonlinear dynamics. Our study demonstrates that the memory effect can be unmasked or greatly exacerbated under certain diseased conditions, which promotes complex excitation dynamics, such as chaos. The iterated map models reveal that memory converts a monotonic iterated map function into a nonmonotonic one to promote the bifurcations leading to high periodicity and chaos.

Memory-Induced Chaos in Cardiac Excitation

Excitable systems display memory, but how memory affects the excitation dynamics of such systems remains to be elucidated. Here we use computer simulation of cardiac action potential models to demonstrate that memory can cause dynamical instabilities that result in complex excitation dynamics and chaos. We develop an iterated map model that correctly describes these dynamics and show that memory converts a monotonic first return map of action potential duration into a nonmonotonic one, resulting in a period-doubling bifurcation route to chaos.

Interactions between Activation and Repolarization Restitution Properties in the Intact Human Heart: In-Vivo Whole-Heart Data and Mathematical Description

BACKGROUND

The restitution of the action potential duration (APDR) and conduction velocity (CVR) are mechanisms whereby cardiac excitation and repolarization adapt to changes in heart rate. They modulate the vulnerability to dangerous arrhythmia, but the mechanistic link between restitution and arrhythmogenesis remains only partially understood.

METHODS

This paper provides an experimental and theoretical study of repolarization and excitation restitution properties and their interactions in the intact human epicardium. The interdependence between excitation and repolarization dynamic is studied in 8 patients (14 restitution protocols, 1722 restitution curves) undergoing global epicardial mapping with multi-electrode socks before open heart surgery. A mathematical description of the contribution of both repolarization and conduction dynamics to the steepness of the APDR slope is proposed.

RESULTS

This study demonstrates that the APDR slope is a function of both activation and repolarization dynamics. At short cycle length, conduction delay significantly increases the APDR slope by interacting with the diastolic interval. As predicted by the proposed mathematical formulation, the APDR slope was more sensitive to activation time prolongation than to the simultaneous shortening of repolarization time. A steep APDR slope was frequently identified, with 61% of all cardiac sites exhibiting an APDR slope > 1, suggesting that a slope > 1 may not necessarily promote electrical instability in the human epicardium. APDR slope did not change for different activation or repolarization times, and it was not a function of local baseline APD. However, it was affected by the spatial organization of electrical excitation, suggesting that in tissue APDR is not a unique function of local electrophysiological properties. Spatial heterogeneity in both activation and repolarization restitution contributed to the increase in the modulated dispersion of repolarization, which for short cycle length was as high as 250 ms. Heterogeneity in conduction velocity restitution can translate into both activation and repolarization dispersion and increase cardiac instability. The proposed mathematical formulation shows an excellent agreement with the experimental data (correlation coefficient r = 0.94) and provides a useful tool for the understanding of the complex interactions between activation and repolarization restitution properties as well as between their measurements.

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that is triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media with applications to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact for cardiac arrhythmias.