Control of cardiac alternans by mechanical and electrical feedback

Unsupervised stochastic learning and reduced order modeling for global sensitivity analysis in cardiac electrophysiology models

BACKGROUND AND OBJECTIVE

Numerical simulations in electrocardiology are often affected by various uncertainties inherited from the lack of precise knowledge regarding input values including those related to the cardiac cell model, domain geometry, and boundary or initial conditions used in the mathematical modeling. Conventional techniques for uncertainty quantification in modeling electrical activities of the heart encounter significant challenges, primarily due to the high computational costs associated with fine temporal and spatial scales. Additionally, the need for numerous model evaluations to quantify ubiquitous uncertainties increases the computational challenges even further.

METHODS

In the present study, we propose a non-intrusive surrogate model to perform uncertainty quantification and global sensitivity analysis in cardiac electrophysiology models. The proposed method combines an unsupervised machine learning technique with the polynomial chaos expansion to reconstruct a surrogate model for the propagation and quantification of uncertainties in the electrical activity of the heart. The proposed methodology not only accurately quantifies uncertainties at a very low computational cost but more importantly, it captures the targeted quantity of interest as either the whole spatial field or the whole temporal period. In order to perform sensitivity analysis, aggregated Sobol indices are estimated directly from the spectral mode of the polynomial chaos expansion.

RESULTS

We conduct Uncertainty Quantification (UQ) and global Sensitivity Analysis (SA) considering both spatial and temporal variations, rather than limiting the analysis to specific Quantities of Interest (QoIs). To assess the comprehensive performance of our methodology in simulating cardiac electrical activity, we utilize the monodomain model. Additionally, sensitivity analysis is performed on the parameters of the Mitchell-Schaeffer cell model.

CONCLUSIONS

Unlike conventional techniques for uncertainty quantification in modeling electrical activities, the proposed methodology performs at a low computational cost the sensitivity analysis on the cardiac electrical activity parameters. The results are fully reproducible and easily accessible, while the proposed reduced-order model represents a significant contribution to enhancing global sensitivity analysis in cardiac electrophysiology.

Thermal modulation of epicardial Ca2+ dynamics uncovers molecular mechanisms of Ca2+ alternans

Ca²⁺ alternans (Ca-Alts) are alternating beat-to-beat changes in the amplitude of Ca²⁺ transients that frequently occur during tachycardia, ischemia, or hypothermia that can lead to sudden cardiac death. Ca-Alts appear to result from a variation in the amount of Ca²⁺ released from the sarcoplasmic reticulum (SR) between two consecutive heartbeats. This variable Ca²⁺ release has been attributed to the alternation of the action potential duration, delay in the recovery from inactivation of RYR Ca²⁺ release channel (RYR2), or an incomplete Ca²⁺ refilling of the SR. In all three cases, the RYR2 mobilizes less Ca²⁺ from the SR in an alternating manner, thereby generating an alternating profile of the Ca²⁺ transients. We used a new experimental approach, fluorescence local field optical mapping (FLOM), to record at the epicardial layer of an intact heart with subcellular resolution. In conjunction with a local cold finger, a series of images were recorded within an area where the local cooling induced a temperature gradient. Ca-Alts were larger in colder regions and occurred without changes in action potential duration. Analysis of the change in the enthalpy and Q₁₀ of several kinetic processes defining intracellular Ca²⁺ dynamics indicated that the effects of temperature change on the relaxation of intracellular Ca²⁺ transients involved both passive and active mechanisms. The steep temperature dependency of Ca-Alts during tachycardia suggests Ca-Alts are generated by insufficient SERCA-mediated Ca²⁺ uptake into the SR. We found that Ca-Alts are heavily dependent on intra-SR Ca²⁺ and can be promoted through partial pharmacologic inhibition of SERCA2a. Finally, the FLOM experimental approach has the potential to help us understand how arrhythmogenesis correlates with the spatial distribution of metabolically impaired myocytes along the myocardium.

A Simulation Study of the Role of Mechanical Stretch in Arrhythmogenesis during Cardiac Alternans

The deformation of the heart tissue due to the contraction can modulate the excitation, a phenomenon referred to as mechanoelectrical feedback (MEF), via stretch-activated channels. The effects of MEF on the electrophysiology at high pacing rates are shown to be proarrhythmic in general. However, more studies need to be done to elucidate the underlying mechanism. In this work, we investigate the effects of MEF on cardiac alternans, which is an alternation in the width of the action potential that typically occurs when the heart is paced at high rates, using a biophysically detailed electromechanical model of cardiac tissue. We observe that the transition from spatially concordant alternans to spatially discordant alternans, which is more arrhythmogenic than concordant alternans, may occur in the presence of MEF and when its strength is sufficiently large. We show that this transition is due to the increase of the dispersion of conduction velocity. In addition, our results also show that the MEF effects, depending on the stretch-activated channels’ conductances and reversal potentials, can result in blocking action potential propagation.

A survey of pathways for mechano-electric coupling in the atria

Mechano-electric coupling (MEC) in atrial tissue has received sparse investigation to date, despite the well-known association between chronic atrial dilation and atrial fibrillation (AF). Of note, no fewer than six different mechanisms pertaining to stretch-activated channels, cellular capacitance and geometric effects have been identified in the literature as potential players. In this mini review, we briefly survey each of these pathways to MEC. We then perform computational simulations using single cell and tissue models in presence of various stretch regimes and MEC pathways. This allows us to assess the relative significance of each pathway in determining action potential duration, conduction velocity and rotor stability. For chronic atrial stretch, we find that stretch-induced alterations in membrane capacitance decrease conduction velocity and increase action potential duration, in agreement with experimental findings. In the presence of time-dependent passive atrial stretch, stretch-activated channels play the largest role, leading to after-depolarizations and rotor hypermeandering. These findings suggest that physiological atrial stretches, such as passive stretch during the atrial reservoir phase, may play an important part in the mechanisms of atrial arrhythmogenesis.

•

Passive strains caused by ventricular contraction need to be considered when incorporating mechano-electro feedback in atrial electrophysiology models.

•

In chronic stretch, stretch-induced capacitance changes dominate.

•

Chronic stretch leads to an increase in action potential duration and a reduction in conduction velocity, consistent with experimental studies.

•

In the presence of passive stretch, stretch-activated channels can induce delayed after-depolarisations and lead to rotor hypermeandering.

•

Mechano-electro feedback is thus likely to have implications for the genesis and maintenance of atrial arrhythmias.

Cardiac Mechano-Electric Coupling: Acute Effects of Mechanical Stimulation on Heart Rate and Rhythm

The heart is vital for biological function in almost all chordates, including humans. It beats continually throughout our life, supplying the body with oxygen and nutrients while removing waste products. If it stops, so does life. The heartbeat involves precise coordination of the activity of billions of individual cells, as well as their swift and well-coordinated adaption to changes in physiological demand. Much of the vital control of cardiac function occurs at the level of individual cardiac muscle cells, including acute beat-by-beat feedback from the local mechanical environment to electrical activity (as opposed to longer term changes in gene expression and functional or structural remodeling). This process is known as mechano-electric coupling (MEC). In the current review, we present evidence for, and implications of, MEC in health and disease in human; summarize our understanding of MEC effects gained from whole animal, organ, tissue, and cell studies; identify potential molecular mediators of MEC responses; and demonstrate the power of computational modeling in developing a more comprehensive understanding of ‟what makes the heart tick.ˮ.

Effects of mechano-electrical feedback on the onset of alternans: A computational study

Cardiac alternans is a heart rhythm instability that is associated with cardiac arrhythmias and may lead to sudden cardiac death. The onset of this instability, which is linked to period-doubling bifurcation and may be a route to chaos, is of particular interest. Mechano-electric feedback depicts the effects of tissue deformation on cardiac excitation. The main effect of mechano-electric feedback is delivered via the so-called stretch-activated ion channels and is caused by stretch-activated currents. Mechano-electric feedback, which is believed to have proarrhythmic and antiarrhythmic effects on cardiac electrophysiology, affects the action potential duration in a manner dependent on cycle length, but the mechanisms by which this occurs remain to be elucidated. In this study, a biophysically detailed electromechanical model of cardiac tissue is employed to show how a stretch-activated current can affect the action potential duration at cellular and tissue levels, illustrating its effects on the onset of alternans. Also, using a two-dimensional iterated map that incorporates stretch-activated current effects, we apply linear stability analysis to study the stability of the bifurcation. We show that alternans bifurcation can be prevented depending on the strength of the stretch-activated current.

Mechanical perturbation control of cardiac alternans

Cardiac alternans is a disturbance in heart rhythm that is linked to the onset of lethal cardiac arrhythmias. Mechanical perturbation control has been recently used to suppress alternans in cardiac tissue of relevant size. In this control strategy, cardiac tissue mechanics are perturbed via active tension generated by the heart's electrical activity, which alters the tissue's electric wave profile through mechanoelectric coupling. We analyze the effects of mechanical perturbation on the dynamics of a map model that couples the membrane voltage and active tension systems at the cellular level. Therefore, a two-dimensional iterative map of the heart beat-to-beat dynamics is introduced, and a stability analysis of the system of coupled maps is performed in the presence of a mechanical perturbation algorithm. To this end, a bidirectional coupling between the membrane voltage and active tension systems in a single cardiac cell is provided, and a discrete form of the proposed control algorithm, that can be incorporated in the coupled maps, is derived. In addition, a realistic electromechanical model of cardiac tissue is employed to explore the feasibility of suppressing alternans at cellular and tissue levels. Electrical activity is represented in two detailed ionic models, the Luo-Rudy 1 and the Fox models, while two active contractile tension models, namely a smooth variant of the Nash-Panfilov model and the Niederer-Hunter-Smith model, are used to represent mechanical activity in the heart. The Mooney-Rivlin passive elasticity model is employed to describe passive mechanical behavior of the myocardium.

Transient and periodic spatiotemporal structures in a reaction-diffusion-mechanics system

We study transient spatiotemporal structures induced by a weak space-time localized stimulus in an excitable contractile fiber within a two-component globally coupled reaction-diffusion model. The model which we develop allows us to analyze various regimes of excitation spreading and determine origin of the induced structures for various contraction types (defined by the fiber fixation) and global coupling strengths. One of the most notable effects we observed is the after-excitation effect. It leads to emergence of multiple excitation pulses excited by a single external stimulus and can result in long-lasting transient activity and appearance of new oscillatory attractor regimes, including the ones with multiple phase clusters.

Control of cardiac alternans in an electromechanical model of cardiac tissue

Electrical alternations in cardiac action potential duration have been shown to be a precursor to arrhythmias and sudden cardiac death. Through the mechanism of excitation-contraction coupling, the presence of electrical alternans induces alternations in the heart muscle contractile activity. Also, contraction of cardiac tissue affects the process of cardiac electric wave propagation through the mechanism of the so-called mechanoelectrical feedback. Electrical excitation and contraction of cardiac tissue can be linked by an electromechanical model such as the Nash-Panfilov model. In this work, we explore the feasibility of suppressing cardiac alternans in the Nash-Panfilov model which is employed for small and large deformations. Several electrical pacing and mechanical perturbation feedback strategies are considered to demonstrate successful suppression of alternans on a one-dimensional cable. This is the first attempt to combine electrophysiologically relevant cardiac models of electrical wave propagation and contractility of cardiac tissue in a synergistic effort to suppress cardiac alternans. Numerical examples are provided to illustrate the feasibility and the effects of the proposed algorithms to suppress cardiac alternans.

Coexisting chaotic and multi-periodic dynamics in a model of cardiac alternans

The spatiotemporal dynamics of cardiac tissue is an active area of research for biologists, physicists, and mathematicians. Of particular interest is the study of period-doubling bifurcations and chaos due to their link with cardiac arrhythmogenesis. In this paper, we study the spatiotemporal dynamics of a recently developed model for calcium-driven alternans in a one dimensional cable of tissue. In particular, we observe in the cable coexistence of regions with chaotic and multi-periodic dynamics over wide ranges of parameters. We study these dynamics using global and local Lyapunov exponents and spatial trajectory correlations. Interestingly, near nodes-or phase reversals-low-periodic dynamics prevail, while away from the nodes, the dynamics tend to be higher-periodic and eventually chaotic. Finally, we show that similar coexisting multi-periodic and chaotic dynamics can also be observed in a detailed ionic model.