Model for the onset of fibrillation following coronary artery occlusion

Paradoxical Onset of Arrhythmic Waves from Depolarized Areas in Cardiac Tissue Due to Curvature-Dependent Instability

The generation of abnormal excitations in pathological regions of the heart is a main trigger for lethal cardiac arrhythmias. Such abnormal excitations, also called ectopic activity, often arise from areas with local tissue heterogeneity or damage accompanied by localized depolarization. Finding the conditions that lead to ectopy is important to understand the basic biophysical principles underlying arrhythmia initiation and might further refine clinical procedures. In this study, we are the first to address the question of how geometry of the abnormal region affects the onset of ectopy using a combination of experimental, in silico, and theoretical approaches. We paradoxically find that, for any studied geometry of the depolarized region in optogenetically modified monolayers of cardiac cells, primary ectopic excitation originates at areas of maximal curvature of the boundary, where the stimulating electrotonic currents are minimal. It contradicts the standard critical nucleation theory applied to nonlinear waves in reaction-diffusion systems, where a higher stimulus is expected to produce excitation more easily. Our in silico studies reveal that the nonconventional ectopic activity is caused by an oscillatory instability at the boundary of the damaged region, the occurrence of which depends on the curvature of that boundary. The onset of this instability is confirmed using the Schrödinger equation methodology proposed by Rinzel and Keener [SIAM J. Appl. Math. 43, 907 (1983)]. Overall, we show distinctively novel insight into how the geometry of a heterogeneous cardiac region determines ectopic activity, which can be used in the future to predict the conditions that can trigger cardiac arrhythmias.

Mechanisms of ventricular arrhythmias: from molecular fluctuations to electrical turbulence

Ventricular arrhythmias have complex causes and mechanisms. Despite extensive investigation involving many clinical, experimental, and computational studies, effective biological therapeutics are still very limited. In this article, we review our current understanding of the mechanisms of ventricular arrhythmias by summarizing the state of knowledge spanning from the molecular scale to electrical wave behavior at the tissue and organ scales and how the complex nonlinear interactions integrate into the dynamics of arrhythmias in the heart. We discuss the challenges that we face in synthesizing these dynamics to develop safe and effective novel therapeutic approaches.

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.

Modeling atrial arrhythmias: impact on clinical diagnosis and therapies

Atrial arrhythmias are the most frequent sustained rhythm disorders in humans and often lead to severe complications such as heart failure and stroke. Despite the important insights provided by animal models into the mechanisms of atrial arrhythmias, direct translation of experimental findings to new therapies in patients has not been straightforward. With the advances in computer technology, large-scale electroanatomical computer models of the atria that integrate information from the molecular to organ scale have reached a level of sophistication that they can be used to interpret the outcome of experimental and clinical studies and aid in the rational design of therapies. This paper reviews the state-of-the-art of computer models of the electrical dynamics of the atria and discusses the evolving role of simulation in assisting the clinical diagnosis and treatment of atrial arrhythmias.

Computing the stability of steady-state solutions of mathematical models of the electrical activity in the heart

Instabilities in the electro-chemical resting state of the heart can generate ectopic waves that in turn can initiate arrhythmias. We derive methods for computing the resting state for mathematical models of the electro-chemical process underpinning a heartbeat, and we estimate the stability of the resting state by invoking the largest real part of the eigenvalues of a linearized model. The implementation of the methods is described and a number of numerical experiments illustrate the feasibility of the methods. In particular, we test the methods for problems where we can compare the solutions with analytical results, and problems where we have solutions computed by independent software. The software is also tested for a fairly realistic 3D model.

Automaticity in acute ischemia: bifurcation analysis of a human ventricular model

Acute ischemia (restriction in blood supply to part of the heart as a result of myocardial infarction) induces major changes in the electrophysiological properties of the ventricular tissue. Extracellular potassium concentration ([K(o)(+)]) increases in the ischemic zone, leading to an elevation of the resting membrane potential that creates an "injury current" (I(S)) between the infarcted and the healthy zone. In addition, the lack of oxygen impairs the metabolic activity of the myocytes and decreases ATP production, thereby affecting ATP-sensitive potassium channels (I(Katp)). Frequent complications of myocardial infarction are tachycardia, fibrillation, and sudden cardiac death, but the mechanisms underlying their initiation are still debated. One hypothesis is that these arrhythmias may be triggered by abnormal automaticity. We investigated the effect of ischemia on myocyte automaticity by performing a comprehensive bifurcation analysis (fixed points, cycles, and their stability) of a human ventricular myocyte model [K. H. W. J. ten Tusscher and A. V. Panfilov, Am. J. Physiol. Heart Circ. Physiol. 291, H1088 (2006)] as a function of three ischemia-relevant parameters [K(o)(+)], I(S), and I(Katp). In this single-cell model, we found that automatic activity was possible only in the presence of an injury current. Changes in [K(o)(+)] and I(Katp) significantly altered the bifurcation structure of I(S), including the occurrence of early-after depolarization. The results provide a sound basis for studying higher-dimensional tissue structures representing an ischemic heart.

A condition for setting off ectopic waves in computational models of excitable cells

The purpose of this paper is to study the stability of steady state solutions of the Monodomain model equipped with Luo-Rudy I kinetics. It is well established that re-entrant arrhythmias can be created in computational models of excitable cells. Such arrhythmias can be initiated by applying an external stimulus that interacts with a partially refractory region, and spawn breaking waves that can eventually generate extremely complex wave patterns commonly referred to as fibrillation. An ectopic wave is one possible stimulus that may initiate fibrillation. Physiologically, it is well known that ectopic waves exist, but the mechanism for initiating ectopic waves in a large collection of cells is poorly understood. In the present paper we consider computational models of collections of excitable cells in one and two spatial dimensions. The cells are modeled by Luo-Rudy I kinetics, and we assume that the spatial dynamics is governed by the Monodomain model. The mathematical analysis is carried out for a reduced model that is known to provide good approximations of the initial phase of solutions of the Luo-Rudy I model. A further simplification is also introduced to motivate and explain the results for the more complicated models. In the analysis the cells are divided into two regions; one region (N) consists of normal cells as model by the standard Luo-Rudy I model, and another region (A) where the cells are automatic in the sense that they would act as pacemaker cells if they where isolated from their surroundings. We let delta denote the spatial diffusion and a denote a characteristic length of the automatic region. It has previously been shown that reducing diffusion or increasing the automatic region enhances ectopic activity. Here we derive a condition for the transition from stable resting state to ectopic wave spread. Under suitable assumptions on the model we provide mathematical and computational arguments indicating that there is a constant eta such that a steady state solution of this system is stable whenever delta approximately > etaa(2), and unstable whenever delta approximately < etaa(2).

Pacemaker activity resulting from the coupling with nonexcitable cells

Fibroblasts are nonexcitable cells that are sometimes coupled with excitable cells (cardiomyocytes). Due to a higher resting potential, these cells may act as a current source or sink and therefore disturb the electrical activity of the surrounding excitable cells. The possible occurrence of spontaneous pacemaker activity resulting from these electrotonic interactions was investigated in a theoretical model of two coupled cells as well as in a multicellular fiber model based on the Courtemanche kinetics. The results indicate that repeated spontaneous activations can be observed after an alteration in the activation and recovery properties of the sodium current (changes in excitability properties), provided that the difference in the resting potential as well as the coupling between the excitable and nonexcitable cells is sufficiently high. This may constitute a mechanism of focal sources triggering arrhythmias such as atrial fibrillation.

Reentry wave formation in excitable media with stochastically generated inhomogeneities

Clinical research shows that the frequency of arrhythmia events depends on the number and area of the border zones of infarct scars. We investigate the possibility that arrhythmia is initiated by reentry waves generated by the inhomogeneity of conduction velocity at the border zone. The interaction of a plane wave with a spatially extended inhomogeneity is simulated in the FitzHugh- Nagumo model. The inhomogeneity is introduced into the model by modifying the spatial dependence of the diffusion coefficient in a stochastic manner. This results in a rich variety of spatial distributions of conductivity. A plane wave propagating through such a system may break up on the regions with low conductivity and produce numerous spiral waves. The frequency of reentry wave formation is studied as a function of the parameters of the inhomogeneity generation algorithm. Three main scenarios of reentry wave formation were found: unidirectional block, main wave-wavelet collision, and wave break up during collision, on a region in which a conduction velocity gradient occurs. These scenarios are likely candidates for the mechanisms of arrhythmia initiation in a damaged tissue, e.g., the border zone of an infarct scar.