Control of electrical alternans in canine cardiac purkinje fibers

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

Role of ephaptic coupling in discordant alternans domain sizes and action potential propagation in the heart

Discordant alternans, the spatially out-of-phase alternation of the durations of propagating action potentials in the heart, has been linked to the onset of fibrillation, a major cardiac rhythm disorder. The sizes of the regions, or domains, within which these alternations are synchronized are critical in this link. However, computer models employing standard gap junction-based coupling between cells have been unable to reproduce simultaneously the small domain sizes and rapid action potential propagation speeds seen in experiments. Here we use computational methods to show that rapid wave speeds and small domain sizes are possible when a more detailed model of intercellular coupling that accounts for so-called ephaptic effects is used. We provide evidence that the smaller domain sizes are possible, because different coupling strengths can exist on the wavefronts, for which both ephaptic and gap-junction coupling are involved, in contrast to the wavebacks, where only gap-junction coupling plays an active role. The differences in coupling strength are due to the high density of fast-inward (sodium) channels known to localize on the ends of cardiac cells, which are only active (and thus engage ephaptic coupling) during wavefront propagation. Thus, our results suggest that this distribution of fast-inward channels, as well as other factors responsible for the critical involvement of ephaptic coupling in wave propagation, including intercellular cleft spacing, play important roles in increasing the vulnerability of the heart to life-threatening tachyarrhythmias. Our results, combined with the absence of short-wavelength discordant alternans domains in standard gap-junction-dominated coupling models, also provide evidence that both gap-junction and ephaptic coupling are critical in wavefront propagation and waveback dynamics.

A machine-learning approach for long-term prediction of experimental cardiac action potential time series using an autoencoder and echo state networks

Computational modeling and experimental/clinical prediction of the complex signals during cardiac arrhythmias have the potential to lead to new approaches for prevention and treatment. Machine-learning (ML) and deep-learning approaches can be used for time-series forecasting and have recently been applied to cardiac electrophysiology. While the high spatiotemporal nonlinearity of cardiac electrical dynamics has hindered application of these approaches, the fact that cardiac voltage time series are not random suggests that reliable and efficient ML methods have the potential to predict future action potentials. This work introduces and evaluates an integrated architecture in which a long short-term memory autoencoder (AE) is integrated into the echo state network (ESN) framework. In this approach, the AE learns a compressed representation of the input nonlinear time series. Then, the trained encoder serves as a feature-extraction component, feeding the learned features into the recurrent ESN reservoir. The proposed AE-ESN approach is evaluated using synthetic and experimental voltage time series from cardiac cells, which exhibit nonlinear and chaotic behavior. Compared to the baseline and physics-informed ESN approaches, the AE-ESN yields mean absolute errors in predicted voltage 6-14 times smaller when forecasting approximately 20 future action potentials for the datasets considered. The AE-ESN also demonstrates less sensitivity to algorithmic parameter settings. Furthermore, the representation provided by the feature-extraction component removes the requirement in previous work for explicitly introducing external stimulus currents, which may not be easily extracted from real-world datasets, as additional time series, thereby making the AE-ESN easier to apply to clinical data.

Pro-Arrhythmic Effects of Discontinuous Conduction at the Purkinje Fiber-Ventricle Junction Arising From Heart Failure-Induced Ionic Remodeling - Insights From Computational Modelling

Heart failure is associated with electrical remodeling of the electrical properties and kinetics of the ion channels and transporters that are responsible for cardiac action potentials. However, it is still unclear whether heart failure-induced ionic remodeling can affect the conduction of excitation waves at the Purkinje fiber-ventricle junction contributing to pro-arrhythmic effects of heart failure, as the complexity of the heart impedes a detailed experimental analysis. The aim of this study was to employ computational models to investigate the pro-arrhythmic effects of heart failure-induced ionic remodeling on the cardiac action potentials and excitation wave conduction at the Purkinje fiber-ventricle junction. Single cell models of canine Purkinje fiber and ventricular myocytes were developed for control and heart failure. These single cell models were then incorporated into one-dimensional strand and three-dimensional wedge models to investigate the effects of heart failure-induced remodeling on propagation of action potentials in Purkinje fiber and ventricular tissue and at the Purkinje fiber-ventricle junction. This revealed that heart failure-induced ionic remodeling of Purkinje fiber and ventricular tissue reduced conduction safety and increased tissue vulnerability to the genesis of the unidirectional conduction block. This was marked at the Purkinje fiber-ventricle junction, forming a potential substrate for the genesis of conduction failure that led to re-entry. This study provides new insights into proarrhythmic consequences of heart failure-induced ionic remodeling.

Controllability and state feedback control of a cardiac ionic cell model

A phenomenon called alternans, which is a beat-to-beat alternation in action potential (AP) duration, sometimes precedes fatal cardiac arrhythmias. Alternans-suppressing electrical stimulus protocols are often represented as perturbations to the dynamics of membrane potential or AP duration variables in nonlinear models of cardiac tissue. Controllability analysis has occasionally been applied to cardiac AP models to determine whether different control or perturbation strategies are capable of suppressing alternans or other unwanted behavior. Since almost all previous cardiac controllability studies have focused on low-dimensional models, we conducted the present study to assess controllability of a higher-dimensional model, specifically the Luo Rudy dynamic (LRd) model of a cardiac ventricular myocyte. Higher-dimensional models are of interest because they provide information on the influence of a wider range of measurable quantities, including ionic concentrations, on controllability. After computing modal controllability measures, we found that larger eigenvalues of a linearized LRd model were on average more strongly controllable through perturbations to calcium-ion concentrations compared with perturbations to other variables. When only membrane potential was adjusted, the best time to apply perturbations (in the sense of maximizing controllability of the largest alternans eigenvalue) was near the AP peak time for shorter cycle lengths. Controllability results were found to be similar for both the default model parameters and for an alternans-promoting parameter set. Additionally, we developed several alternans-suppressing state feedback controllers that were tested in simulations. For the scenarios examined, our controllability measures correctly predicted which strategies and perturbation timings would lead to better feedback controller performance.

Long-Time Prediction of Arrhythmic Cardiac Action Potentials Using Recurrent Neural Networks and Reservoir Computing

The electrical signals triggering the heart's contraction are governed by non-linear processes that can produce complex irregular activity, especially during or preceding the onset of cardiac arrhythmias. Forecasts of cardiac voltage time series in such conditions could allow new opportunities for intervention and control but would require efficient computation of highly accurate predictions. Although machine-learning (ML) approaches hold promise for delivering such results, non-linear time-series forecasting poses significant challenges. In this manuscript, we study the performance of two recurrent neural network (RNN) approaches along with echo state networks (ESNs) from the reservoir computing (RC) paradigm in predicting cardiac voltage data in terms of accuracy, efficiency, and robustness. We show that these ML time-series prediction methods can forecast synthetic and experimental cardiac action potentials for at least 15–20 beats with a high degree of accuracy, with ESNs typically two orders of magnitude faster than RNN approaches for the same network size.

A space-fractional bidomain framework for cardiac electrophysiology: 1D alternans dynamics

Cardiac electrophysiology modeling deals with a complex network of excitable cells forming an intricate syncytium: the heart. The electrical activity of the heart shows recurrent spatial patterns of activation, known as cardiac alternans, featuring multiscale emerging behavior. On these grounds, we propose a novel mathematical formulation for cardiac electrophysiology modeling and simulation incorporating spatially non-local couplings within a physiological reaction-diffusion scenario. In particular, we formulate, a space-fractional electrophysiological framework, extending and generalizing similar works conducted for the monodomain model. We characterize one-dimensional excitation patterns by performing an extended numerical analysis encompassing a broad spectrum of space-fractional derivative powers and various intra- and extracellular conductivity combinations. Our numerical study demonstrates that (i) symmetric properties occur in the conductivity parameters' space following the proposed theoretical framework, (ii) the degree of non-local coupling affects the onset and evolution of discordant alternans dynamics, and (iii) the theoretical framework fully recovers classical formulations and is amenable for parametric tuning relying on experimental conduction velocity and action potential morphology.

Data-driven phase-isostable reduction for optimal nonfeedback stabilization of cardiac alternans

Phase-isostable reduction is an emerging model reduction strategy that can be used to accurately replicate nonlinear behaviors in systems for which standard phase reduction techniques fail. In this work, we derive relationships between the cycle-to-cycle variance of the reduced isostable coordinates for systems subject to both additive white noise and periodic stimulation. Using this information, we propose a data-driven technique for inferring nonlinear terms of the phase-isostable coordinate reduction framework. We apply the proposed model inference strategy to the biologically motivated problem of eliminating cardiac alternans, an arrhythmia that is widely considered to be a precursor to more deadly cardiac arrhythmias. Using this strategy, by simply measuring a series of action potential durations in response to periodic stimulation, we are able to identify energy-optimal, nonfeedback control inputs to stabilize a period-1, alternans-free solution.

Clinical Potential of Beat-to-Beat Diastolic Interval Control in Preventing Cardiac Arrhythmias

Life‐threatening ventricular arrhythmias and sudden cardiac death are often preceded by cardiac alternans, a beat‐to‐beat oscillation in the T‐wave morphology or duration. However, given the spatiotemporal and structural complexity of the human heart, designing algorithms to effectively suppress alternans and prevent fatal rhythms is challenging. Recently, an antiarrhythmic constant diastolic interval pacing protocol was proposed and shown to be effective in suppressing alternans in 0‐, 1‐, and 2‐dimensional in silico studies as well as in ex vivo whole heart experiments. Herein, we provide a systematic review of the electrophysiological conditions and mechanisms that enable constant diastolic interval pacing to be an effective antiarrhythmic pacing strategy. We also demonstrate a successful translation of the constant diastolic interval pacing protocol into an ECG‐based real‐time control system capable of modulating beat‐to‐beat cardiac electrical activity and preventing alternans. Furthermore, we present evidence of the clinical utility of real‐time alternans suppression in reducing arrhythmia susceptibility in vivo. We provide a comprehensive overview of this promising pacing technique, which can potentially be translated into a clinically viable device that could radically improve the quality of life of patients experiencing abnormal cardiac rhythms.

Controllability of voltage- and calcium-driven cardiac alternans in a map model

Certain cardiac arrhythmias are preceded by electrical alternans, a state characterized by beat-to-beat alternation in cellular action potential duration. Cardiac alternans may arise from different mechanisms including instabilities in voltage or intracellular calcium cycling. Although a number of techniques have been proposed to suppress alternans, these methods have mainly been tested using models that do not support calcium-driven alternans. Therefore, it is important to understand how control methods may perform when alternans is driven by instabilities in calcium cycling. In this study, we applied controllability analysis to a discrete map of alternans dynamics in a cardiac cell. We compared two different controllability measures to determine to what extent different control strategies could suppress alternans and tested these predictions using three feedback controllers. We found a modal controllability measure, unlike the minimum singular value of the controllability matrix, consistently indicated the control strategies requiring the least control effort and yielding the smallest closed-loop eigenvalue. In addition, action potential duration was identified as the most effective variable through which control can be applied, regardless of alternans mechanism, although sarcoplasmic reticulum calcium load was also useful for the calcium-driven alternans cases.

Universal mechanisms for self-termination of rapid cardiac rhythm

Excitable media sustain circulating waves. In the heart, sustained circulating waves can lead to serious impairment or even death. To investigate factors affecting the stability of such waves, we have used optogenetic techniques to stimulate a region at the apex of a mouse heart at a fixed delay after the detection of excitation at the base of the heart. For long delays, rapid circulating rhythms can be sustained, whereas for shorter delays, there are paroxysmal bursts of activity that start and stop spontaneously. By considering the dependence of the action potential and conduction velocity on the preceding recovery time using restitution curves, as well as the reduced excitability (fatigue) due to the rapid excitation, we model prominent features of the dynamics including alternation of the duration of the excited phases and conduction times, as well as termination of the bursts for short delays. We propose that this illustrates universal mechanisms that exist in biological systems for the self-termination of such activities.

Stability of spatially discordant repolarization alternans in cardiac tissue

Cardiac alternans, a period-2 behavior of excitation and contraction of the heart, is a precursor of ventricular arrhythmias and sudden cardiac death. One form of alternans is repolarization or action potential duration alternans. In cardiac tissue, repolarization alternans can be spatially in-phase, called spatially concordant alternans, or spatially out-of-phase, called spatially discordant alternans (SDA). In SDA, the border between two out-of-phase regions is called a node in a one-dimensional cable or a nodal line in a two-dimensional tissue. In this study, we investigate the stability and dynamics of the nodes and nodal lines of repolarization alternans driven by voltage instabilities. We use amplitude equation and coupled map lattice models to derive theoretical results, which are compared with simulation results from the ionic model. Both conduction velocity restitution induced SDA and non-conduction velocity restitution induced SDA are investigated. We show that the stability and dynamics of the SDA nodes or nodal lines are determined by the balance of the tensions generated by conduction velocity restitution, convection due to action potential propagation, curvature of the nodal lines, and repolarization and coupling heterogeneities. Our study provides mechanistic insights into the different SDA behaviors observed in experiments.

Effect of constant-DI pacing on single cell pacing dynamics

Cardiac alternans, beat-to-beat alternations in action potential duration, is a precursor to fatal arrhythmias such as ventricular fibrillation. Previous research has shown that voltage driven alternans can be suppressed by application of a constant diastolic interval (DI) pacing protocol. However, the effect of constant-DI pacing on cardiac cell dynamics and its interaction with the intracellular calcium cycle remains to be determined. Therefore, we aimed to examine the effects of constant-DI pacing on the dynamical behavior of a single-cell numerical model of cardiac action potential and the influence of voltage-calcium (V-Ca) coupling on it. Single cell dynamics were analyzed in the vicinity of the bifurcation point using a hybrid pacing protocol, a combination of constant-basic cycle length (BCL) and constant-DI pacing. We demonstrated that in a small region beneath the bifurcation point, constant-DI pacing caused the cardiac cell to remain alternans-free after switching to the constant-BCL pacing, thus introducing a region of bistability (RB). The size of the RB increased with stronger V-Ca coupling and was diminished with weaker V-Ca coupling. Overall, our findings demonstrate that the application of constant-DI pacing on cardiac cells with strong V-Ca coupling may induce permanent changes to cardiac cell dynamics increasing the utility of constant-DI pacing.

Global vs local control of cardiac alternans in a 1D numerical model of human ventricular tissue

Cardiac alternans is a proarrhythmic state in which the action potential duration (APD) of cardiac myocytes alternate between long and short values and often occurs under conditions of rapid pacing of cardiac tissue. In the ventricles, alternans is especially dangerous due to the life-threatening risk of developing arrhythmias, such as ventricular fibrillation. Alternans can be formed in periodically paced tissue as a result of pacing itself. Recently, it has been demonstrated that this pacing-induced alternans can be prevented by performing constant diastolic interval (DI) pacing, in which DI is independent of APD. However, constant DI pacing is difficult to implement in experimental settings since it requires the real-time measurement of APD. A more practical way was proposed based on electrocardiograms (ECGs), which give an indirect measure of the global DI relaxation period through the TR interval assessment. Previously, we demonstrated that constant TR pacing prevented alternans formation in isolated Langendorff-perfused rabbit hearts. However, the efficacy of "local" constant DI pacing vs "global" constant TR pacing in preventing alternans formation has never been investigated. Thus, the purpose of this study was to implement an ECG-based constant TR pacing in a 1D numerical model of human ventricular tissue and to compare the dynamical behavior of cardiac tissue with that resulted from a constant DI pacing. The results showed that both constant TR and constant DI pacing prevented the onset of alternans until lower basic cycle length when compared to periodic pacing. For longer cable lengths, constant TR pacing was shown to exhibit greater control on alternans than constant DI pacing.

Feedback control of calcium driven alternans in cardiac myocytes

Cardiac alternans is a beat-to-beat alternation of the action potential duration (APD), which has been implicated as a possible cause of ventricular fibrillation. Previous studies have shown that alternans can originate via a period doubling bifurcation caused by the nonlinear dependence of the APD on the previous diastolic interval. In this case, it has been demonstrated that alternans can be eliminated by applying feedback control on the pacing cycle length. However, studies have shown that alternans can also originate due to unstable calcium (Ca) cycling in cardiac myocytes. In this study, we explore the effectiveness of APD feedback control to suppress alternans when the underlying instability is due to unstable Ca cycling. In particular, we explore the role of the bi-directional coupling between Ca and voltage and determine the effectiveness of feedback control under a wide range of conditions. We also analyze the applicability of feedback control on a coupled two cell system and show that APD control induces spatially out-of-phase alternans. We analyze the onset and the necessary conditions for the emergence of these out-of-phase patterns and assess the effectiveness of feedback control to suppress Ca driven alternans in a multi-cellular system.

Key aspects for effective mathematical modelling of fractional-diffusion in cardiac electrophysiology: a quantitative study

Microscopic structural features of cardiac tissue play a fundamental role in determining complex spatio-temporal excitation dynamics at the macroscopic level. Recent efforts have been devoted to the development of mathematical models accounting for non-local spatio-temporal coupling able to capture these complex dynamics without the need of resolving tissue heterogeneities down to the micro-scale. In this work, we analyse in detail several important aspects affecting the overall predictive power of these modelling tools and provide some guidelines for an effective use of space-fractional models of cardiac electrophysiology in practical applications. Through an extensive computational study in simplified computational domains, we highlight the robustness of models belonging to different categories, i.e., physiological and phenomenological descriptions, against the introduction of non-locality, and lay down the foundations for future research and model validation against experimental data. A modern genetic algorithm framework is used to investigate proper parameterisations of the considered models, and the crucial role played by the boundary assumptions in the considered settings is discussed. Several numerical results are provided to support our claims.

Spatially Discordant Repolarization Alternans in the Absence of Conduction Velocity Restitution

Spatially discordant alternans (SDA) of action potential duration (APD) has been widely observed in cardiac tissue and is linked to cardiac arrhythmogenesis. Theoretical studies have shown that conduction velocity restitution (CVR) is required for the formation of SDA. However, this theory is not completely supported by experiments, indicating that other mechanisms may exist. In this study, we carried out computer simulations using mathematical models of action potentials to investigate the mechanisms of SDA in cardiac tissue. We show that when CVR is present and engaged, such as fast pacing from one side of the tissue, the spatial pattern of APD in the tissue undergoes either spatially concordant alternans or SDA, independent of initial conditions or tissue heterogeneities. When CVR is not engaged, such as simultaneous pacing of the whole tissue or under normal/slow heart rates, the spatial pattern of APD in the tissue can have multiple solutions, including spatially concordant alternans and different SDA patterns, depending on heterogeneous initial conditions or pre-existing repolarization heterogeneities. In homogeneous tissue, curved nodal lines are not stable, which either evolve into straight lines or disappear. However, in heterogeneous itssue, curved nodal lines can be stable, depending on their initial locations and shapes relative to the structure of the heterogeneity. Therefore, CVR-induced SDA and non-CVR-induced SDA exhibit different dynamical properties, which may be responsible for the different SDA properties observed in experimental studies and arrhythmogenesis in different clinical settings.

Delay differential equation-based models of cardiac tissue: Efficient implementation and effects on spiral-wave dynamics

Delay differential equations (DDEs) recently have been used in models of cardiac electrophysiology, particularly in studies focusing on electrical alternans, instabilities, and chaos. A number of processes within cardiac cells involve delays, and DDEs can potentially represent mechanisms that result in complex dynamics both at the cellular level and at the tissue level, including cardiac arrhythmias. However, DDE-based formulations introduce new computational challenges due to the need for storing and retrieving past values of variables at each spatial location. Cardiac tissue simulations that use DDEs may require over 28 GB of memory if the history of variables is not managed carefully. This paper addresses both computational and dynamical issues. First, we present new methods for the numerical solution of DDEs in tissue to mitigate the memory requirements associated with the history of variables. The new methods exploit the different time scales of an action potential to dynamically optimize history size. We find that the proposed methods decrease memory usage by up to 95% in cardiac tissue simulations compared to straightforward history-management algorithms. Second, we use the optimized methods to analyze for the first time the dynamics of wave propagation in two-dimensional cardiac tissue for models that include DDEs. In particular, we study the effects of DDEs on spiral-wave dynamics, including wave breakup and chaos, using a canine myocyte model. We find that by introducing delays to the gating variables governing the calcium current, DDEs can induce spiral-wave breakup in 2D cardiac tissue domains.

Simulating waves, chaos and synchronization with a microcontroller

The spatiotemporal dynamics of complex systems have been studied traditionally and visualized numerically using high-end computers. However, due to advances in microcontrollers, it is now possible to run what once were considered large-scale simulations using a very small and inexpensive single integrated circuit that can furthermore send and receive information to and from the outside world in real time. In this paper, we show how microcontrollers can be used to perform simulations of nonlinear ordinary differential equations with spatial coupling and to visualize their dynamics using arrays of light-emitting diodes and/or touchscreens. We demonstrate these abilities using three different models: two reaction-diffusion models (one neural and one cardiac) and a generic model of network oscillators. These models are commonly used to simulate various phenomena in biophysical systems, including bifurcations, waves, chaos, and synchronization. We also demonstrate how simple it is to integrate real-time user interaction with the simulations by showing examples with a light sensor, touchscreen, and web browser.

Optimal phase control of biological oscillators using augmented phase reduction

We develop a novel optimal control algorithm to change the phase of an oscillator using a minimum energy input, which also minimizes the oscillator's transversal distance to the uncontrolled periodic orbit. Our algorithm uses a two-dimensional reduction technique based on both isochrons and isostables. We develop a novel method to eliminate cardiac alternans by connecting our control algorithm with the underlying physiological problem. We also describe how the devised algorithm can be used for spike timing control which can potentially help with motor symptoms of essential and parkinsonian tremor, and aid in treating jet lag. To demonstrate the advantages of this algorithm, we compare it with a previously proposed optimal control algorithm based on standard phase reduction for the Hopf bifurcation normal form, and models for cardiac pacemaker cells, thalamic neurons, and circadian gene regulation cycle in the suprachiasmatic nucleus. We show that our control algorithm is effective even when a large phase change is required or when the nontrivial Floquet multiplier is close to unity; in such cases, the previously proposed control algorithm fails.

Phase reduction and phase-based optimal control for biological systems: a tutorial

A powerful technique for the analysis of nonlinear oscillators is the rigorous reduction to phase models, with a single variable describing the phase of the oscillation with respect to some reference state. An analog to phase reduction has recently been proposed for systems with a stable fixed point, and phase reduction for periodic orbits has recently been extended to take into account transverse directions and higher-order terms. This tutorial gives a unified treatment of such phase reduction techniques and illustrates their use through mathematical and biological examples. It also covers the use of phase reduction for designing control algorithms which optimally change properties of the system, such as the phase of the oscillation. The control techniques are illustrated for example neural and cardiac systems.

Isostable reduction of oscillators with piecewise smooth dynamics and complex Floquet multipliers

Phase-amplitude reduction is a widely applied technique in the study of limit cycle oscillators with the ability to represent a complicated and high-dimensional dynamical system in a more analytically tractable set of coordinates. Recent work has focused on the use of isostable coordinates, which characterize the transient decay of solutions toward a periodic orbit, and can ultimately be used to increase the accuracy of these reduced models. The breadth of systems to which this phase-amplitude reduction strategy can be applied, however, is still rather limited. In this work, the theory of phase-amplitude reduction using isostable coordinates is further developed to accommodate a broader set of dynamical systems. In the first part, limit cycles of piecewise smooth dynamical systems are considered and strategies are developed to compute the associated reduced equations. In the second part, the notion of isostable coordinates for complex-valued Floquet multipliers is introduced, resulting in one phaselike coordinate and one amplitudelike coordinate for each pair of complex conjugate Floquet multipliers. Examples are given with relevance to piecewise smooth representations of excitable cardiomyocytes and the relationship between the reduced coordinate system and the emergence of cardiac alternans is discussed. Also, phase-amplitude reduction is implemented for a chaotic, externally forced pendulum with complex Floquet multipliers and a resulting control strategy for the stabilization of its periodic solution is investigated.

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.

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.

Real-time optical manipulation of cardiac conduction in intact hearts

Key points

Although optogenetics has clearly demonstrated the feasibility of cardiac manipulation, current optical stimulation strategies lack the capability to react acutely to ongoing cardiac wave dynamics.

Here, we developed an all‐optical platform to monitor and control electrical activity in real‐time.

The methodology was applied to restore normal electrical activity after atrioventricular block and to manipulate the intraventricular propagation of the electrical wavefront.

The closed‐loop approach was also applied to simulate a re‐entrant circuit across the ventricle.

The development of this innovative optical methodology provides the first proof‐of‐concept that a real‐time all‐optical stimulation can control cardiac rhythm in normal and abnormal conditions.

Abstract

Optogenetics has provided new insights in cardiovascular research, leading to new methods for cardiac pacing, resynchronization therapy and cardioversion. Although these interventions have clearly demonstrated the feasibility of cardiac manipulation, current optical stimulation strategies do not take into account cardiac wave dynamics in real time. Here, we developed an all‐optical platform complemented by integrated, newly developed software to monitor and control electrical activity in intact mouse hearts. The system combined a wide‐field mesoscope with a digital projector for optogenetic activation. Cardiac functionality could be manipulated either in free‐run mode with submillisecond temporal resolution or in a closed‐loop fashion: a tailored hardware and software platform allowed real‐time intervention capable of reacting within 2 ms. The methodology was applied to restore normal electrical activity after atrioventricular block, by triggering the ventricle in response to optically mapped atrial activity with appropriate timing. Real‐time intraventricular manipulation of the propagating electrical wavefront was also demonstrated, opening the prospect for real‐time resynchronization therapy and cardiac defibrillation. Furthermore, the closed‐loop approach was applied to simulate a re‐entrant circuit across the ventricle demonstrating the capability of our system to manipulate heart conduction with high versatility even in arrhythmogenic conditions. The development of this innovative optical methodology provides the first proof‐of‐concept that a real‐time optically based stimulation can control cardiac rhythm in normal and abnormal conditions, promising a new approach for the investigation of the (patho)physiology of the heart.

Although optogenetics has clearly demonstrated the feasibility of cardiac manipulation, current optical stimulation strategies lack the capability to react acutely to ongoing cardiac wave dynamics.

Here, we developed an all‐optical platform to monitor and control electrical activity in real‐time.

The methodology was applied to restore normal electrical activity after atrioventricular block and to manipulate the intraventricular propagation of the electrical wavefront.

The closed‐loop approach was also applied to simulate a re‐entrant circuit across the ventricle.

The development of this innovative optical methodology provides the first proof‐of‐concept that a real‐time all‐optical stimulation can control cardiac rhythm in normal and abnormal conditions.

Mechanisms linking T-wave alternans to spontaneous initiation of ventricular arrhythmias in rabbit models of long QT syndrome

KEY POINTS

T-wave alternans (TWA) and T-wave lability (TWL) are precursors of ventricular arrhythmias in long QT syndrome; however, the mechanistic link remains to be clarified. Computer simulations show that action potential duration (APD) prolongation and slowed heart rates promote APD alternans and chaos, manifesting as TWA and TWL, respectively. Regional APD alternans and chaos can exacerbate pre-existing or induce de novo APD dispersion, which combines with enhanced I_Ca,L to result in premature ventricular complexes (PVCs) originating from the APD gradient region. These PVCs can directly degenerate into re-entrant arrhythmias without the need for an additional tissue substrate or further exacerbate the APD dispersion to cause spontaneous initiation of ventricular arrhythmias. Experiments conducted in transgenic long QT rabbits show that PVC alternans occurs at slow heart rates, preceding spontaneous intuition of ventricular arrhythmias.

ABSTRACT

T-wave alternans (TWA) and irregular beat-to-beat T-wave variability or T-wave lability (TWL), the ECG manifestations of action potential duration (APD) alternans and variability, are precursors of ventricular arrhythmias in long QT syndromes. TWA and TWL in patients tend to occur at normal heart rates and are usually potentiated by bradycardia. Whether or how TWA and TWL at normal or slow heart rates are causally linked to arrhythmogenesis remains unknown. In the present study, we used computer simulations and experiments of a transgenic rabbit model of long QT syndrome to investigate the underlying mechanisms. Computer simulations showed that APD prolongation and slowed heart rates caused early afterdepolarization-mediated APD alternans and chaos, manifesting as TWA and TWL, respectively. Regional APD alternans and chaos exacerbated pre-existing APD dispersion and, in addition, APD chaos could also induce APD dispersion de novo via chaos desynchronization. Increased APD dispersion, combined with substantially enhanced I_Ca,L , resulted in a tissue-scale dynamical instability that gave rise to the spontaneous occurrence of unidirectionally propagating premature ventricular complexes (PVCs) originating from the APD gradient region. These PVCs could directly degenerate into re-entrant arrhythmias without the need for an additional tissue substrate or could block the following sinus beat to result in a longer RR interval, which further exacerbated the APD dispersion giving rise to the spontaneous occurrence of ventricular arrhythmias. Slow heart rate-induced PVC alternans was observed in experiments of transgenic LQT2 rabbits under isoproterenol, which was associated with increased APD dispersion and spontaneous occurrence of ventricular arrhythmias, in agreement with the theoretical predictions.

Real-Time Closed Loop Diastolic Interval Control Prevents Cardiac Alternans in Isolated Whole Rabbit Hearts

Cardiac alternans, a beat-to-beat alternation in action potential duration (APD), can lead to fatal arrhythmias. During periodic pacing, changes in diastolic interval (DI) depend on subsequent changes in APD, thus enhancing cardiac instabilities through a 'feedback' mechanism. Recently, an anti-arrhythmic Constant DI pacing protocol was proposed and shown to be effective in suppressing alternans in 0D and 1D in silico studies. However, previous experimental validation of Constant DI pacing in the heart has been unsuccessful due to the spatio-temporal complexity of 2D cardiac tissue and the technical challenges in its real-time implementation. Here, we developed a novel closed loop system to detect T-waves from real-time ECG data, enabling successful implementation of Constant DI pacing protocol, and performed high-resolution optical mapping experiments on isolated whole rabbit hearts to validate its anti-arrhythmic effects. The results were compared with: (1) Periodic pacing (feedback inherent) and (2) pacing with heart rate variability (HRV) (feedback modulation) introduced by using either Gaussian or Physiological patterns. We observed that Constant DI pacing significantly suppressed alternans in the heart, while maintaining APD spatial dispersion and flattening the slope of the APD restitution curve, compared to traditional Periodic pacing. In addition, introduction of HRV in Periodic pacing failed to prevent cardiac alternans, and was arrhythmogenic.

Stochastic Pacing Inhibits Spatially Discordant Cardiac Alternans

Depressed heart rate variability is a well-established risk factor for sudden cardiac death in survivors of acute myocardial infarction and for those with congestive heart failure. Although measurements of heart rate variability provide a valuable prognostic tool, it is unclear whether reduced heart rate variability itself is proarrhythmic or if it simply correlates with the severity of autonomic nervous system dysfunction. In this work, we investigate a possible mechanism by which heart rate variability could protect against cardiac arrhythmia. Specifically, in numerical simulations, we observe an inverse relationship between the variance of stochastic pacing and the occurrence of spatially discordant alternans, an arrhythmia that is widely believed to facilitate the development of cardiac fibrillation. By analyzing the effects of conduction velocity restitution, cellular dynamics, electrotonic coupling, and stochastic pacing on the nodal dynamics of spatially discordant alternans, we provide intuition for this observed behavior and propose control strategies to inhibit discordant alternans.

Alternans promotion in cardiac electrophysiology models by delay differential equations

Cardiac electrical alternans is a state of alternation between long and short action potentials and is frequently associated with harmful cardiac conditions. Different dynamic mechanisms can give rise to alternans; however, many cardiac models based on ordinary differential equations are not able to reproduce this phenomenon. A previous study showed that alternans can be induced by the introduction of delay differential equations (DDEs) in the formulations of the ion channel gating variables of a canine myocyte model. The present work demonstrates that this technique is not model-specific by successfully promoting alternans using DDEs for five cardiac electrophysiology models that describe different types of myocytes, with varying degrees of complexity. By analyzing results across the different models, we observe two potential requirements for alternans promotion via DDEs for ionic gates: (i) the gate must have a significant influence on the action potential duration and (ii) a delay must significantly impair the gate's recovery between consecutive action potentials.

Theory of the development of alternans in the heart during controlled diastolic interval pacing

The beat-to-beat alternation in action potential durations (APDs) in the heart, called APD alternans, has been linked to the development of serious cardiac rhythm disorders, including ventricular tachycardia and fibrillation. The length of the period between action potentials, called the diastolic interval (DI), is a key dynamical variable in the standard theory of alternans development. Thus, methods that control the DI may be useful in preventing dangerous cardiac rhythms. In this study, we examine the dynamics of alternans during controlled-DI pacing using a series of single-cell and one-dimensional (1D) fiber models of alternans dynamics. We find that a model that combines a so-called memory model with a calcium cycling model can reasonably explain two key experimental results: the possibility of alternans during constant-DI pacing and the phase lag of APDs behind DIs during sinusoidal-DI pacing. We also find that these results can be replicated by incorporating the memory model into an amplitude equation description of a 1D fiber. The 1D fiber result is potentially concerning because it seems to suggest that constant-DI control of alternans can only be effective over only a limited region in space.

Constant DI pacing suppresses cardiac alternans formation in numerical cable models

Cardiac repolarization alternans describe the sequential alternation of the action potential duration (APD) and can develop during rapid pacing. In the ventricles, such alternans may rapidly turn into life risking arrhythmias under conditions of spatial heterogeneity. Thus, suppression of alternans by artificial pacing protocols, or alternans control, has been the subject of numerous theoretical, numerical, and experimental studies. Yet, previous attempts that were inspired by chaos control theories were successful only for a short spatial extent (<2 cm) from the pacing electrode. Previously, we demonstrated in a single cell model that pacing with a constant diastolic interval (DI) can suppress the formation of alternans at high rates of activation. We attributed this effect to the elimination of feedback between the pacing cycle length and the last APD, effectively preventing restitution-dependent alternans from developing. Here, we extend this idea into cable models to study the extent by which constant DI pacing can control alternans during wave propagation conditions. Constant DI pacing was applied to ventricular cable models of up to 5 cm, using human kinetics. Our results show that constant DI pacing significantly shifts the onset of both cardiac alternans and conduction blocks to higher pacing rates in comparison to pacing with constant cycle length. We also demonstrate that constant DI pacing reduces the propensity of spatially discordant alternans, a precursor of wavebreaks. We finally found that the protective effect of constant DI pacing is stronger for increased electrotonic coupling along the fiber in the sense that the onset of alternans is further shifted to higher activation rates. Overall, these results support the potential clinical applicability of such type of pacing in improving protocols of implanted pacemakers, in order to reduce the risk of life-threatening arrhythmias. Future research should be conducted in order to experimentally validate these promising results.

Isostable reduction of periodic orbits

The well-established method of phase reduction neglects information about a limit-cycle oscillator's approach towards its periodic orbit. Consequently, phase reduction suffers in practicality unless the magnitude of the Floquet multipliers of the underlying limit cycle are small in magnitude. By defining isostable coordinates of a periodic orbit, we present an augmentation to classical phase reduction which obviates this restriction on the Floquet multipliers. This framework allows for the study and understanding of periodic dynamics for which standard phase reduction alone is inadequate. Most notably, isostable reduction allows for a convenient and self-contained characterization of the dynamics near unstable periodic orbits.

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.

Effects and underlying mechanisms of refractory period pacing on repolarization dynamics in the human heart

Repolarization alternans is related to the initiation of life threatening cardiac arrhythmias. Experimental and computational studies suggest that the abolishment of alternans using dynamic pacing protocols may prevent abnormal heart rhythms. In a recent animal study, refractory period pacing (RPP) on every other beat has shown promising results in alternans reduction. However, the cellular mechanisms underlying this therapy and its efficiency in human patients remain unclear. In this study, in vivo unipolar electrograms acquired during RPP from 240 epicardial sites from one patient were analysed. Current clamp of 18 channels was performed in silico to elucidate the ionic mechanisms underlying action potential modulation by RPP. Its efficacy with positive and negative polarities was tested on a population of 87 calibrated human ventricular models exhibiting alternans. In vivo electrograms showed significant changes in T-wave alternans when applying RPP. In silico, results showed APD shortening for RPP with positive polarity and APD prolongation with RPP negative. Under current clamp protocols, voltage rectification of L-type Ca2+ (ICaL) and inward rectifier K+ (IK1) currents were identified as the key determinants for the observed changes. RPP pacing successfully reduced alternans on the in silico models using a negative polarity stimulus in the short beat.

Isostable reduction with applications to time-dependent partial differential equations

Isostables and isostable reduction, analogous to isochrons and phase reduction for oscillatory systems, are useful in the study of nonlinear equations which asymptotically approach a stationary solution. In this work, we present a general method for isostable reduction of partial differential equations, with the potential power to reduce the dimensionality of a nonlinear system from infinity to 1. We illustrate the utility of this reduction by applying it to two different models with biological relevance. In the first example, isostable reduction of the Fokker-Planck equation provides the necessary framework to design a simple control strategy to desynchronize a population of pathologically synchronized oscillatory neurons, as might be relevant to Parkinson's disease. Another example analyzes a nonlinear reaction-diffusion equation with relevance to action potential propagation in a cardiac system.

Image-Based Structural Modeling of the Cardiac Purkinje Network

The Purkinje network is a specialized conduction system within the heart that ensures the proper activation of the ventricles to produce effective contraction. Its role during ventricular arrhythmias is less clear, but some experimental studies have suggested that the Purkinje network may significantly affect the genesis and maintenance of ventricular arrhythmias. Despite its importance, few structural models of the Purkinje network have been developed, primarily because current physical limitations prevent examination of the intact Purkinje network. In previous modeling efforts Purkinje-like structures have been developed through either automated or hand-drawn procedures, but these networks have been created according to general principles rather than based on real networks. To allow for greater realism in Purkinje structural models, we present a method for creating three-dimensional Purkinje networks based directly on imaging data. Our approach uses Purkinje network structures extracted from photographs of dissected ventricles and projects these flat networks onto realistic endocardial surfaces. Using this method, we create models for the combined ventricle-Purkinje system that can fully activate the ventricles through a stimulus delivered to the Purkinje network and can produce simulated activation sequences that match experimental observations. The combined models have the potential to help elucidate Purkinje network contributions during ventricular arrhythmias.

Role of sodium and calcium dysregulation in tachyarrhythmias in sudden cardiac death

Despite improvements in the therapy of underlying heart disease, sudden cardiac death is a major cause of death worldwide. Disturbed Na and Ca handling is known to be a major predisposing factor for life-threatening tachyarrhythmias. In cardiomyocytes, many ion channels and transporters, including voltage-gated Na and Ca channels, cardiac ryanodine receptors, Na/Ca-exchanger, and SR Ca-ATPase are involved in this regulation. We have learned a lot about the pathophysiological relevance of disturbed ion channel function from monogenetic disorders. Changes in the gating of a single ion channel and the activity of an ion pump suffice to dramatically increase the propensity for arrhythmias even in structurally normal hearts. Nevertheless, patients with heart failure with acquired dysfunction in many ion channels and transporters exhibit profound dysregulation of Na and Ca handling and Ca/calmodulin-dependent protein kinase and are especially prone to arrhythmias. A deeper understanding of the underlying arrhythmic principles is mandatory if we are to improve their outcome. This review addresses basic tachyarrhythmic mechanisms, the underlying ionic mechanisms and the consequences for ion homeostasis, and the situation in complex diseases like heart failure.

Determining individual phase response curves from aggregate population data

Phase reduction is an invaluable technique for investigating the dynamics of nonlinear limit cycle oscillators. Central to the implementation of phase reduction is the ability to calculate phase response curves (PRCs), which describe an oscillator's response to an external perturbation. Current experimental techniques for inferring PRCs require data from individual oscillators, which can be impractical to obtain when the oscillator is part of a much larger population. Here we present a simple methodology to calculate PRCs of individual oscillators using an aggregate signal from a large homogeneous population. This methodology is shown to be accurate in the presence of interoscillator coupling and noise and can also provide a good estimate of an average PRC of a heterogeneous population. We also find that standard experimental techniques for PRC measurement can produce misleading results when applied to aggregate population data.

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.

Voltage and calcium dynamics both underlie cellular alternans in cardiac myocytes

Cardiac alternans, a putative trigger event for cardiac reentry, is a beat-to-beat alternation in membrane potential and calcium transient. Alternans was originally attributed to instabilities in transmembrane ion channel dynamics (i.e., the voltage mechanism). As of this writing, the predominant view is that instabilities in subcellular calcium handling are the main underlying mechanism. That being said, because the voltage and calcium systems are bidirectionally coupled, theoretical studies have suggested that both mechanisms can contribute. To date, to our knowledge, no experimental evidence of such a dual role within the same cell has been reported. Here, a combined electrophysiological and calcium imaging approach was developed and used to illuminate the contributions of voltage and calcium dynamics to alternans. An experimentally feasible protocol, quantification of subcellular calcium alternans and restitution slope during cycle-length ramping alternans control, was designed and validated. This approach allows simultaneous illumination of the contributions of voltage and calcium-driven instability to total cellular instability as a function of cycle-length. Application of this protocol in in vitro guinea-pig left-ventricular myocytes demonstrated that both voltage- and calcium-driven instabilities underlie alternans, and that the relative contributions of the two systems change as a function of pacing rate.

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.

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.

Continuous-time control of alternans in long Purkinje fibers

Alternans-an arrhythmic response of cardiac tissue to periodic pacing-often serves as a precursor to a more dangerous, and potentially lethal, state of fibrillation. Suppression of alternans using feedback control may be a plausible method to prevent fibrillation. Several approaches based on impulsive control have been proposed previously, where feedback is applied for a brief instance of time during each pacing interval. This paper presents a continuous-time approach, where feedback current is applied at all times, which is capable of suppressing alternans in fibers of significantly greater length (up to at least 4 cm), compared with impulsive control (less than 1 cm), and for a wide range of pacing cycle lengths.

Control of cardiac alternans by mechanical and electrical feedback

A persistent alternation in the cardiac action potential duration has been linked to the onset of ventricular arrhythmia, which may lead to sudden cardiac death. A coupling between these cardiac alternans and the intracellular calcium dynamics has also been identified in previous studies. In this paper, the system of PDEs describing the small amplitude of alternans and the alternation of peak intracellular Ca(2+) are stabilized by optimal boundary and spatially distributed actuation. A simulation study demonstrating the successful annihilation of both alternans on a one-dimensional cable of cardiac cells by utilizing the full-state feedback controller is presented. Complimentary to these studies, a three variable Nash-Panfilov model is used to investigate alternans annihilation via mechanical (or stretch) perturbations. The coupled model includes the active stress which defines the mechanical properties of the tissue and is utilized in the feedback algorithm as an independent input from the pacing based controller realization in alternans annihilation. Simulation studies of both control methods demonstrate that the proposed methods can successfully annihilate alternans in cables that are significantly longer than 1 cm, thus overcoming the limitations of earlier control efforts.

Spatiotemporal dynamics of calcium-driven cardiac alternans

We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We focus on the common case where the bidirectional coupling of intracellular calcium concentration and membrane voltage dynamics produces calcium and voltage alternans that are temporally in phase. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolarization and calcium transient (CaT) alternans, similarly to the well-studied case of voltage-driven alternans. In contrast, further away from the bifurcation, the amplitude of CaT alternans jumps discontinuously at the nodes separating out-of-phase regions, while the amplitude of repolarization alternans remains smooth. We identify universal dynamical features of SDA pattern formation and evolution in the presence of those jumps. We show that node motion of discontinuous SDA patterns is strongly hysteretic even in homogeneous tissue due to the novel phenomenon of "unidirectional pinning": node movement can only be induced towards, but not away from, the pacing site in response to a change of pacing rate or physiological parameter. In addition, we show that the wavelength of discontinuous SDA patterns scales linearly with the conduction velocity restitution length scale, in contrast to the wavelength of smooth patterns that scales sublinearly with this length scale. Those results are also shown to be robust against cell-to-cell fluctuations due to the property that unidirectional node motion collapses multiple jumps accumulating in nodal regions into a single jump. Amplitude equation predictions are in good overall agreement with ionic model simulations. Finally, we briefly discuss physiological implications of our findings. In particular, we suggest that due to the tendency of conduction blocks to form near nodes, the presence of unidirectional pinning makes calcium-driven alternans potentially more arrhythmogenic than voltage-driven alternans.

The past, present, and future of real-time control in cellular electrophysiology

For over 60 years, real-time control has been an important technique in the study of excitable cells. Two such control-based technologies are reviewed here. First, voltage-clamp methods revolutionized the study of excitable cells. In this family of techniques, membrane potential is controlled, allowing one to parameterize a powerful class of models that describe the voltage-current relationship of cell membranes simply, flexibly, and accurately. Second, dynamic-clamp methods allow the addition of new, "virtual" membrane mechanisms to living cells. Dynamic clamp allows researchers unprecedented ways of testing computationally based hypotheses in biological preparations. The review ends with predictions of how control-based technologies will be improved and adapted for new uses in the near future.