Effect of heterogeneities in the cellular microstructure on propagation of the cardiac action potential

Microheterogeneity-induced conduction slowing and wavefront collisions govern macroscopic conduction behavior: A computational and experimental study

The incidence of cardiac arrhythmias is known to be associated with tissue heterogeneities including fibrosis. However, the impact of microscopic structural heterogeneities on conduction in excitable tissues remains poorly understood. In this study, we investigated how acellular microheterogeneities affect macroscopic conduction under conditions of normal and reduced excitability by utilizing a novel platform of paired in vitro and in silico studies to examine the mechanisms of conduction. Regular patterns of nonconductive micro-obstacles were created in confluent monolayers of the previously described engineered-excitable Ex293 cell line. Increasing the relative ratio of obstacle size to intra-obstacle strand width resulted in significant conduction slowing up to 23.6% and a significant increase in wavefront curvature anisotropy, a measure of spatial variation in wavefront shape. Changes in bulk electrical conductivity and in path tortuosity were insufficient to explain these observed macroscopic changes. Rather, microscale behaviors including local conduction slowing due to microscale branching, and conduction acceleration due to wavefront merging were shown to contribute to macroscopic phenomena. Conditions of reduced excitability led to further conduction slowing and a reversal of wavefront curvature anisotropy due to spatially non-uniform effects on microscopic slowing and acceleration. This unique experimental and computation platform provided critical mechanistic insights in the impact of microscopic heterogeneities on macroscopic conduction, pertinent to settings of fibrotic heart disease.

It is well known that perturbations in the heart structure are associated with the initiation and maintenance of clinically significant cardiac arrhythmia. While previous studies have examined how single structural perturbations affect local electrical conduction, our understanding of how numerous microscopic heterogeneities act in aggregate to alter macroscopic electrical behavior is limited. In this study, we utilized simplified engineered excitable cells that contain the minimal machinery of excitability and can be directly computationally modeled. By pairing experimental and computational studies, we showed that the microscopic branching and collisions of electrical waves slow and speed conduction, respectively, resulting in macroscopic changes in the speed and pattern of electrical activation. These microscale behaviors are significantly altered under reduced excitability, resulting in exaggerated collision effects. Overall, this study helps improve our understanding of how microscopic structural heterogeneities in excitable tissue lead to abnormal action potential propagation, conducive to arrhythmias.

The roles of mid-myocardial and epicardial cells in T-wave alternans development: a simulation study

Background

The occurrence of T-wave alternans in electrocardiographic signals was recently linked to susceptibility to ventricular arrhythmias and sudden cardiac death. Thus, by detecting and comprehending the origins of T-wave alternans, it might be possible to prevent such events.

Results

Here, we simulated T-wave alternans in a computer-generated human heart model by modulating the action potential duration and amplitude during the first part of the repolarization phase. We hypothesized that changes in the intracardiac alternans patterns of action potential properties would differentially influence T-wave alternans measurements at the body surface. Specifically, changes were simulated globally in the whole left and right ventricles to simulate concordant T-wave alternans, and locally in selected regions to simulate discordant and regional discordant, hereinafter referred to as “regional”, T-wave alternans. Body surface potential maps and 12-lead electrocardiographic signals were then computed. In depth discrimination, the influence of epicardial layers on T-wave alternans development was significantly higher than that of mid-myocardial cells. Meanwhile, spatial discrimination revealed that discordant and regional action potential property changes had a higher influence on T-wave alternans amplitude than concordant changes. Notably, varying T-wave alternans sources yielded distinct body surface potential map patterns for T-wave alternans amplitude, which can be used for location of regions within hearts exhibiting impaired repolarization. The highest ability for T-wave alternans detection was achieved in lead V1. Ultimately, we proposed new parameters Vector Magnitude Alternans and Vector Angle Alternans, with higher ability for T-wave alternans detection when using multi-lead electrocardiographic signals processing than for single leads. Finally, QT alternans was found to be associated with the process of T-wave alternans generation.

Conclusions

The distributions of the body surface T-wave alternans amplitude have been shown to have unique patterns depending on the type of alternans (concordant, discordant or regional) and the location of the disturbance in the heart. The influence of epicardial cells on T-wave alternans development is significantly higher than that of mid-myocardial cells, among which the sub-endocardial layer exerted the highest influence. QT interval alternans is identified as a phenomenon that correlate with T-wave alternans.

Modeling dynamics in diseased cardiac tissue: Impact of model choice

Cardiac arrhythmias have been traditionally simulated using continuous models that assume tissue homogeneity and use a relatively large spatial discretization. However, it is believed that the tissue fibrosis and collagen deposition, which occur on a micron-level, are critical factors in arrhythmogenesis in diseased tissues. Consequently, it remains unclear how well continuous models, which use averaged electrical properties, are able to accurately capture complex conduction behaviors such as re-entry in fibrotic tissues. The objective of this study was to compare re-entrant behavior in discrete microstructural models of fibrosis and in two types of equivalent continuous models, a homogenous continuous model and a hybrid continuous model with distinct heterogeneities. In the discrete model, increasing levels of tissue fibrosis lead to a substantial increase in the re-entrant cycle length which is inadequately reflected in the homogenous continuous models. These cycle length increases appear to be primarily due to increases in the tip path length and to altered restitution behavior, and suggest that it is critical to consider the discrete effects of fibrosis on conduction when studying arrhythmogenesis in fibrotic myocardium. Hybrid models are able to accurately capture some aspects of re-entry and, if carefully tuned, may provide a framework for simulating conduction in diseased tissues with both accuracy and efficiency.

Improved robust T-wave alternans detectors

New statistical and spectral detectors, the modified matched pairs t test, the extended spectral method and the modified spectral method, were proposed for T-wave alternans (TWA) detection gaining robustness according to trend and single-frequency interferences. They were compared to classic detectors such as matched pairs t test, unpaired t test, spectral method, generalized likelihood ratio test and estimated TWA amplitude within a simulation framework and applied to real data. The optimal detection threshold was selected by using a full Monte-Carlo simulation where signals, with and without alternans episodes, were corrupted by Gaussian noise with different power and single-frequency interferences with different tones. All the combinations of noise and frequency were selected and repeated 500 times in order to compute probability of detection ([Formula: see text]) and the false alarm probability ([Formula: see text]), providing ROC curves. The study group consisted of 50 patients with implantable cardioverter-defibrillator (age: [Formula: see text]; LVEF: [Formula: see text]), who were paced (ventricular pacing) at 100 bpm. Two-minute recordings were analyzed. The XYZ orthogonal lead system was used. The best performance was reached by using the modified matched pairs t test (in comparison with the spectral method and other reference methods).

Dynamics of propagation of premature impulses in structurally remodeled infarcted myocardium: a computational analysis

Initiation of cardiac arrhythmias typically follows one or more premature impulses either occurring spontaneously or applied externally. In this study, we characterize the dynamics of propagation of single (S2) and double premature impulses (S3), and the mechanisms of block of premature impulses at structural heterogeneities caused by remodeling of gap junctional conductance (Gj) in infarcted myocardium. Using a sub-cellular computer model of infarcted tissue, we found that |I_Na,max|, prematurity (coupling interval with the previous impulse), and conduction velocity (CV) of premature impulses change dynamically as they propagate away from the site of initiation. There are fundamental differences between the dynamics of propagation of S2 and S3 premature impulses: for S2 impulses |I_Na,max| recovers fast, prematurity decreases and CV increases as propagation proceeds; for S3 impulses low values of |I_Na,max| persist, prematurity could increase, and CV could decrease as impulses propagate away from the site of initiation. As a consequence it is more likely that S3 impulses block at sites of structural heterogeneities causing source/sink mismatch than S2 impulses block. Whether premature impulses block at Gj heterogeneities or not is also determined by the values of Gj (and the space constant λ) in the regions proximal and distal to the heterogeneity: when λ in the direction of propagation increases >40%, premature impulses could block. The maximum slope of CV restitution curves for S2 impulses is larger than for S3 impulses. In conclusion: (1) The dynamics of propagation of premature impulses make more likely that S3 impulses block at sites of structural heterogeneities than S2 impulses block; (2) Structural heterogeneities causing an increase in λ (or CV) of >40% could result in block of premature impulses; (3) A decrease in the maximum slope of CV restitution curves of propagating premature impulses is indicative of an increased potential for block at structural heterogeneities.