Dynamic origin of spatially discordant alternans in cardiac tissue

Alternans, a condition in which there is a beat-to-beat alternation in the electromechanical response of a periodically stimulated cardiac cell, has been linked to the genesis of life-threatening ventricular arrhythmias. Optical mapping of membrane voltage (V_m) and intracellular calcium (Ca_i) on the surface of animal hearts reveals complex spatial patterns of alternans. In particular, spatially discordant alternans has been observed in which regions with a large-small-large action potential duration (APD) alternate out-of-phase adjacent to regions of small-large-small APD. However, the underlying mechanisms that lead to the initiation of discordant alternans and govern its spatiotemporal properties are not well understood. Using mathematical modeling, we show that dynamic changes in the spatial distribution of discordant alternans can be used to pinpoint the underlying mechanisms. Optical mapping of V_m and Ca_i in paced rabbit hearts revealed that spatially discordant alternans induced by rapid pacing exhibits properties consistent with a purely dynamical mechanism as shown in theoretical studies. Our results support the viewpoint that spatially discordant alternans in the heart can be formed via a dynamical pattern formation process which does not require tissue heterogeneity.     

Action potential duration and force-frequency relationship in isolated rabbit,guinea pig and rat cardiac muscle

The effect of action potential duration and elevated cytosolic sodium concentration on the forcefrequency relationship in isolated rabbit, guinea pig and rat papillary muscle preparations was studied. Shortening of action potential duration in guinea pig and rabbit from 150–200 ms to values characteristic of rat (20–40 ms), using the K_(ATP) channel activator levkromakalim (15 mol·l^–1), markedly reduced the force of contraction and converted the positive force-frequency relationship into negative one at longer pacing cycle lengths. This conversion was greatly enhanced in the presence of acetylstrophanthidin (0.2–1 mol·l^–1), an inhibitor of the Na-K pump. Acetylstrophanthidin (1 mol·l^–1) alone, however, had no effect on the forcefrequency relationship. Prolongation of action potential duration in rat with inhibitors of cardiac K channels (4-aminopyridine [10 mmol·l^–1] plus tetraethylammonium [2 mmol·l^–1) increased the force of contraction and abolished the negative force-frequency relationship observed in rat at longer pacing-cycle lengths. It is concluded that both action potential duration and cytosolic sodium concentration are major determinants of the force-frequency relationship in mammalian myocardium.Abbreviations AC acetylstrophanthidin - APD action potential duration - APD ₅₀ and APD ₉₀ action potential duration measured at 50% and 90% level of repolarization, respectively - SR sarcoplasmic reticulum     

Reentry in heterogeneous cardiac tissue described by the Luo-Rudy ventricular action potential model

Heterogeneity of cardiac tissue is an important factor determining the initiation and dynamics of cardiac arrhythmias. In this paper, we studied the effects of gradients of electrophysiological heterogeneity on reentrant excitation patterns using computer simulations. We investigated the dynamics of spiral waves in a two-dimensional sheet of cardiac tissue described by the Luo-Rudy phase 1 (LR1) ventricular action potential model. A gradient of action potential duration (APD) was imposed by gradually varying the local current density of K(+) current or inward rectifying K(+) current along one axis of the tissue sheet. We show that a gradient of APD resulted in spiral wave drift. This drift consisted of two components. The longitudinal (along the gradient) component was always directed toward regions of longer spiral wave period. The transverse (perpendicular to the gradient) component had a direction dependent on the direction of rotation of the spiral wave. We estimated the velocity of the drift as a function of the magnitude of the gradient and discuss its implications.     

Action potential propagation in inhomogeneous cardiac tissue: safety factor considerations and ionic mechanism

Heterogeneity of myocardial structure and membrane excitability is accentuated by pathology and remodeling. In this study, a detailed model of the ventricular myocyte in a multicellular fiber was used to compute a location-dependent quantitative measure of conduction (safety factor, SF) and to determine the kinetics and contribution of sodium current (I(Na)) and L-type calcium current [I(Ca(L))] during conduction. We obtained the following results. 1) SF decreases sharply for propagation into regions of increased electrical load (tissue expansion, increased gap junction coupling, reduced excitability, hyperkalemia); it can be <1 locally (a value indicating conduction failure) and can recover beyond the transition region to resume propagation. 2) SF and propagation across inhomogeneities involve major contribution from I(Ca(L)). 3) Modulating I(Na) or I(Ca(L)) (by blocking agents or calcium overload) can cause unidirectional block in the inhomogeneous region. 4) Structural inhomogeneity causes local augmentation of I(Ca(L)) and suppression of I(Na) in a feedback fashion. 5) Propagation across regions of suppressed I(Na) is achieved via a I(Ca(L))-dependent mechanism. 6) Reduced intercellular coupling can effectively compensate for reduced SF caused by tissue expansion but not by reduced membrane excitability.     

Mechanisms of calcium transient and action potential alternans in cardiac cells and tissues

Alternation of cardiac action potential duration (APD) from beat to beat and concurrent alternation of the amplitude of the calcium transient are regarded as important arrhythmia mechanisms. These phenomena are causally interrelated and can be reliably evoked by an increase in beat frequency or by ischemia. The first part of this historical review deals with the physiology of APD alternans. Sections recounting the evolution of knowledge about calcium-activated ion currents and calcium transient alternans are interspersed among sections describing the growth of the so-called "restitution hypothesis," which involves time-dependent recovery of potassium channels (including their passage through pre-open states) as a function of diastolic interval. Major developments are generally in chronological order, but it is necessary to move back and forth between the two theories to respect the overall time line, which runs from about l965 to the present. The concluding two sections deal with the pathophysiology of calcium transient and APD alternans during ischemia, which may be the basis for out-of-hospital cardiac arrest during the initial stages of acute myocardial infarction.     

Propagation of normal beats and re-entry in a computational model of ventricular cardiac tissue with regional differences in action potential shape and duration

There is substantial experimental evidence from studies using both intact tissue and isolated single cells to support the existence of different cell types within the ventricular wall of the heart, each possessing different electrical properties. However other studies have failed to find these differences, and instead support the idea that electrical coupling in vivo between regions with different cell types smoothes out differences in action potential shape and duration. In this study we have used a computational model of electrical activation in heterogenous 2D and 3D cardiac tissue to investigate the propagation of both normal beats and arrhythmias. We used the Luo–Rudy dynamic model for guinea pig ventricular cells, with simplified Ca²⁺ handling and transmural heterogeneity in I_Ks and I_to. With normal cell-to-cell coupling, a layer of M cells was not necessary for the formation of an upright T wave in the simulated electrocardiogram, and the amplitude and configuration of the T wave was not greatly affected by the thickness and configuration of the M cell layer. Transmural gradients in repolarisation pushed re-entrant waves with an intramural filament towards either the base or the apex of the ventricles, and caused transient break up of re-entry with a transmural filament.     

Action potential morphology influences intracellular calcium handling stability and the occurrence of alternans

Instability in the intracellular Ca2+ handling system leading to Ca2+ alternans is hypothesized to be an underlying cause of electrical alternans. The highly coupled nature of membrane voltage and Ca2+ regulation suggests that there should be reciprocal effects of membrane voltage on the stability of the Ca2+ handling system. We investigated such effects using a mathematical model of the cardiac intracellular Ca2+ handling system. We found that the morphology of the action potential has a significant effect on the stability of the Ca2+ handling system at any given pacing rate, with small changes in action potential morphology resulting in a transition from stable nonalternating Ca2+ transients to stable alternating Ca2+ transients. This bifurcation occurs as the alternans eigen value of the system changes from absolute value <1 to absolute value >1. These results suggest that the stability of the intracellular Ca2+ handling system and the occurrence of Ca2+ alternans are not dictated solely by the Ca2+ handling system itself, but are also modulated to a significant degree by membrane voltage (through its influence on sarcolemmal Ca2+ currents) and, therefore, by all ionic currents that affect membrane voltage.     

Electrophysiological heterogeneity and stability of reentry in simulated cardiac tissue

Generation of wave break is a characteristic feature of cardiac fibrillation. In this study, we investigated how dynamic factors and fixed electrophysiological heterogeneity interact to promote wave break in simulated two-dimensional cardiac tissue, by using the Luo-Rudy (LR1) ventricular action potential model. The degree of dynamic instability of the action potential model was controlled by varying the maximal amplitude of the slow inward Ca(2+) current to produce spiral waves in homogeneous tissue that were either nearly stable, meandering, hypermeandering, or in breakup regimes. Fixed electrophysiological heterogeneity was modeled by randomly varying action potential duration over different spatial scales to create dispersion of refractoriness. We found that the degree of dispersion of refractoriness required to induce wave break decreased markedly as dynamic instability of the cardiac model increased. These findings suggest that reducing the dynamic instability of cardiac cells by interventions, such as decreasing the steepness of action potential duration restitution, may still have merit as an antifibrillatory strategy.     

Electrotonic influences on action potential duration dispersion in small hearts: a simulation study

Intrinsic spatial variations in repolarization currents in the heart can produce spatial gradients in action potential duration (APD) that serve as possible sites for conduction block and the initiation of reentrant activity. In well-coupled myocardium, however, electrotonic influences at the stimulus site and wavefront collision sites act to modulate any intrinsic heterogeneity in APD. These effects alter APD gradients over an extent larger than that suggested by the length constant associated with propagation and, thus, are hypothesized to play a greater role in smaller hearts used as experimental models of human disease. This study uses computer simulation to investigate how heart size, tissue properties, and the spatial assignment of cell types affect functional APD dispersion. Simulations were carried out using the murine ventricular myocyte model of Pandit et al. or the Luo-Rudy mammalian model in three-dimensional models of mouse and rabbit ventricular geometries. Results show that the spatial extent of the APD dispersion is related to the dynamic changes in transmembrane resistance during recovery. Also, because of the small dimensions of the mouse heart, electrotonic effects on APD primarily determine the functional dispersion of refractoriness, even in the presence of large intrinsic cellular heterogeneity and reduced coupling. APD dispersion, however, is found to increase significantly when the heart size increases to the size of a rabbit heart, unmasking intrinsic cell types.     

Voltage noise influences action potential duration in cardiac myocytes

Stochastic gating of ion channels introduces noise to membrane currents in cardiac muscle cells (myocytes). Since membrane currents drive membrane potential, noise thereby influences action potential duration (APD) in myocytes. To assess the influence of noise on APD, membrane potential is in this study formulated as a stochastic process known as a diffusion process, which describes both the current-voltage relationship and voltage noise. In this framework, the response of APD voltage noise and the dependence of response on the shape of the current-voltage relationship can be characterized analytically. We find that in response to an increase in noise level, action potential in a canine ventricular myocytes is typically prolonged and that distribution of APDs becomes more skewed towards long APDs, which may lead to an increased frequency of early after-depolarization formation. This is a novel mechanism by which voltage noise may influence APD. The results are in good agreement with those obtained from more biophysically-detailed mathematical models, and increased voltage noise (due to gating noise) may partially underlie an increased incidence of early after-depolarizations in heart failure.     

Nonlinear dynamics of two-dimensional cardiac action potential duration mapping model with memory

The aim of this work is the analysis of the nonlinear dynamics of two-dimensional mapping model of cardiac action potential duration (2D-map APD) with memory derived from one dimensional map (1D-map). Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting cardiac arrhythmias. For a constant rate of stimulation the short action potential during alternans is followed by a longer DI and inversely. It has been suggested that these differences in DI are responsible for the occurrence and maintenance of APD alternans. We focus our attention on the observed bifurcations produced by a change in the stimulation period and a fixed value of a particular parameter in the model. This parameter provides new information about the dynamics of the APD with memory, such as the occurrence of bistabilities not previously described in the literature, as well as the fact that synchronization rhythms occur in different ways and in a new fashion as the stimulation frequency increases. Moreover, we show that this model is flexible enough as to accurately reflect the chaotic dynamics properties of the APD: we have highlighted the fractal structure of the strange attractor of the 2D-map APD, and we have characterized chaos by tools such as the calculation of the Lyapunov exponents, the fractal dimension and the Kolmogorov entropy, with the next objective of refining the study of the nonlinear dynamics of the duration of the action potential and to apply methods of controlling chaos.

     

Effects of transmural electrical heterogeneities and electrotonic interactions on the dispersion of cardiac repolarization and action potential duration: A simulation study

It has been shown in the literature that myocytes isolated from the ventricular walls at various intramural depths have different action potential durations (APDs). When these myocytes are embedded in the ventricular wall, their inhomogeneous properties affect the sequence of repolarization and the actual distribution of the APDs in the entire wall. In this article, we implement a mathematical model to simulate the combined effect of (a) the non-homogeneous intrinsic membrane properties (in particular the non-homogeneous APDs) and (b) the electrotonic currents that modulate the APDs when the myocytes are embedded in the ventricular myocardium. In particular, we study the effect of (a) and (b) on the excitation and repolarization sequences and on the distribution of APDs in the ventricles. We implement a Monodomain tissue representation that includes orthotropic anisotropy, transmural fiber rotation and homogeneous or heterogeneous transmural intrinsic membrane properties, modeled according to the phase I Luo-Rudy membrane ionic model. Three-dimensional simulations are performed in a cartesian slab with a parallel finite element solver employing structured isoparametric trilinear finite elements in space and a semi-implicit adaptive method in time. Simulations of excitation and repolarization sequences elicited by epicardial or endocardial pacing show that in a homogeneous slab the repolarization pathways approximately follow the activation sequence. Conversely, in the heterogeneous cases considered in this study, we observed two repolarization wavefronts that started from the epi and the endocardial faces respectively and collided in the thickness of the wall and in one case an additional repolarization wave starting from an intramural site. Introducing the heterogeneities along the transmural epi-endocardial direction affected both the repolarization sequence and the APD dispersion, but these effects were clearly discernible only in transmural planes. By contrast, in planes parallel to epi- and endocardium the APD distribution remained remarkably similar to that observed in the homogeneous model. Therefore, the patterns of the repolarization sequence and APD dispersion on the epicardial surface (or any other intramural surface parallel to it) do not reveal the uniform transmural heterogeneity.     

Electrical refractory period restitution and spiral wave reentry in simulated cardiac tissue

Theoretical and experimental studies have shown that restitution of the cardiac action potential (AP) duration (APD) plays a major role in predisposing ventricular tachycardia to degenerate to ventricular fibrillation, whereas its role in atrial fibrillation is unclear. We used the Courtemanche human atrial cell model and the Luo-Rudy guinea pig ventricular model to compare the roles of electrical restitution in destabilizing spiral wave reentry in simulated two-dimensional homogeneous atrial and ventricular tissue. Because atrial AP morphology is complex, we also validated the usefulness of effective refractory period (ERP) restitution. ERP restitution correlated best with APD restitution at transmembrane potentials greater than or equal to -62 mV, and its steepness was a reliable predictor of spiral wave phenotype (stable, meandering, hypermeandering, and breakup) in both atrial and ventricular tissue. Spiral breakup or single hypermeandering spirals occurred when the slope of ERP restitution exceeded 1 at short diastolic intervals. Thus ERP restitution, which is easier to measure clinically than APD restitution, is a reliable determinant of spiral wave stability in simulated atrial and ventricular tissue.     

Action potential duration restitution and ventricular fibrillation due to rapid focal excitation

The focal source hypothesis of ventricular fibrillation (VF) posits that rapid activation from a focal source, rather than action potential duration (APD) restitution properties, is responsible for the maintenance of VF. We injected aconitine (100 microg) into normal isolated perfused swine right ventricles (RVs) stained with 4-[beta-[2-(di-n-butylamino)-6-naphthyl]vinyl]pyridinium (di-4-ANEPPS) for optical mapping studies. Within 97 +/- 163 s, aconitine induced ventricular tachycardia (VT) with a mean cycle length 268 +/- 37 ms, which accelerated before converting to VF. Drugs that flatten the APD restitution slope, including diacetyl monoxime (10-20 mM, n = 6), bretylium (10-20 microg/ml, n = 3), and verapamil (2-4 microg/ml, n = 3), reversibly converted VF to VT in all cases. In two RVs, VF persisted despite of the excision of the aconitine site. Simulations in two-dimensional cardiac tissue showed that once VF was initiated, it remained sustained even after the "aconitine" site was eliminated. In this model of focal source VF, the VT-to-VF transition occurred due to a wave break outside the aconitine site, and drugs that flattened the APD restitution slope converted VF to VT despite continuous activation from aconitine site.     

Action potential duration and activation of ATP-sensitive potassium current in isolated guinea-pig ventricular myocytes

It is difficult to associate the ATP-sensitive potassium (K-ATP) channel of cardiac muscle with hypoxia/ischemia induced action potential shortening because this occurs before intracellular ATP falls to levels associated in vitro with channel opening. This leaves the cardiac K-ATP channel without any obvious physiological function. We have quantitatively examined the relationship between action potential duration and K-ATP channel activity in enzymatically isolated ventricular myocytes of the guinea-pig. In whole-cell voltage-clamp recording experiments when the K-ATP channel opener SR 44866 (2-10 microM) stimulated an outward membrane current greater than 50 pA at 0 mV membrane potential (the equivalent of 30 open K-ATP channels or 1% of the cell K-ATP channel population) action potential duration was reduced by more than 50%. In the majority of cell-attached membrane patch recordings metabolic inhibition stimulated K-ATP channel open probability of 1-2% which continued for long periods (7-25 min) before cell contracture and coincident major K-ATP channel activation (open probability 65%). Our quantitative analysis thus shows that physiologically relevant activity of K-ATP channels in cardiac muscle is confined to a very small percentage of the possible cell K-ATP current and thus intracellular ATP would not have to fall very far before the opening of K-ATP channels would influence cardiac excitability.     

Mechanisms and determinants of ultralong action potential duration and slow rate-dependence in cardiac myocytes

In normal cardiac myocytes, the action potential duration (APD) is several hundred milliseconds. However, experimental studies showed that under certain conditions, APD could be excessively long (or ultralong), up to several seconds. Unlike the normal APD, the ultralong APD increases sensitively with pacing cycle length even when the pacing rate is very slow, exhibiting a sensitive slow rate-dependence. In addition, these long action potentials may or may not exhibit early afterdepolarizations (EADs). Although these phenomena are well known, the underlying mechanisms and ionic determinants remain incompletely understood. In this study, computer simulations were performed with a simplified action potential model. Modifications to the L-type calcium current (I(Ca,L)) kinetics and the activation time constant of the delayed rectifier K current were used to investigate their effects on APD. We show that: 1) the ultralong APD and its sensitive slow rate-dependence are determined by the steady-state window and pedestal I(Ca,L) currents and the activation speed and the recovery of the delayed rectifier K current; 2) whether an ultralong action potential exhibits EADs or not depends on the kinetics of I(Ca,L); 3) increasing inward currents elevates the plateau voltage, which in general prolongs APD, however, this can also shorten APD when the APD is already ultralong under certain conditions; and 4) APD alternans occurs at slow pacing rates due to the sensitive slow rate-dependence and the ionic determinants are different from the ones causing APD alternans at fast heart rates.     

Plasmachemical and heterogeneous processes in ozonizers with oxygen activation by a dielectric barrier discharge

Plasmachemical and heterogeneous processes of generation and loss of ozone in the atmosphericpressure dielectric barrier discharge in oxygen are studied theoretically. Plasmachemical and electronic kinetics in the stage of development and decay of a single plasma filament (microdischarge) are calculated numerically with and without allowance for the effects of ozone vibrational excitation and high initial ozone concentration. The developed analytical approach is applied to determine the output ozone concentration taking into account ozone heterogeneous losses on the Al₂O₃ dielectric surface. Using the results of quantummechanical calculations by the method of density functional theory, a multistage catalytic mechanism of heterogeneous ozone loss based on the initial passivation of a pure Al₂O₃ surface by ozone and the subsequent interaction of O₃ molecules with the passivated surface is proposed. It is shown that the conversion reaction 2O₃ → 3O₂ of a gas-phase ozone molecule with a physically adsorbed ozone molecule can result in the saturation of the maximum achievable ozone concentration at high specific energy depositions, the nonstationarity of the output ozone concentration, and its dependence on the prehistory of ozonizer operation.     

Simulating patterns of excitation, repolarization and action potential duration with cardiac Bidomain and Monodomain models

Parallel numerical simulations of excitation and recovery in three-dimensional myocardial domains are presented. The simulations are based on the anisotropic Bidomain and Monodomain models, including intramural fiber rotation and orthotropic or axisymmetric anisotropy of the intra- and extra-cellular conductivity tensors. The Bidomain model consist of a system of two reaction-diffusion equations, while the Monodomain model consists of one reaction-diffusion equation. Both models are coupled with the phase I Luo-Rudy membrane model describing the ionic currents. Simulations of excitation and repolarization sequences on myocardial slabs of different sizes show how the distribution of the action potential durations (APD) is influenced by both the anisotropic electrical conduction and the fiber rotation. This influence occurs in spite of the homogeneous intrinsic properties of the cell membrane. The APD dispersion patterns are closely correlated to the anisotropic curvature of the excitation wavefront.