Suppression of alternans and conduction blocks despite steep APD restitution: electrotonic, memory, and conduction velocity restitution effects

We examine the utility of the action potential (AP) duration (APD) restitution curve slope in predicting the onset of electrical alternans when electrotonic and memory effects are considered. We develop and use two ionic cell models without memory that have the same restitution curve with slope >1 but different AP shapes and, therefore, different electrotonic effects. We also study a third cell model that incorporates short-term memory of previous cycle lengths, so that it has a family of S1-S2 restitution curves as well as a dynamic restitution curve with slope >1. Our results indicate that both electrotonic and memory effects can suppress alternans, even when the APD restitution curve is steep. In the absence of memory, electrotonic currents related to the shape of the AP, as well as conduction velocity restitution, can affect how alternans develops in tissue and, in some cases, can prevent its induction entirely, even when isolated cells exhibit alternans. When short-term memory is included, alternans may not occur in isolated cells, despite a steep APD restitution curve, and may or may not occur in tissue, depending on conduction velocity restitution. We show for the first time that electrotonic and memory effects can prevent conduction blocks and stabilize reentrant waves in two and three dimensions. Thus we find that the slope of the APD restitution curve alone does not always well predict the onset of alternans and that incorporating electrotonic and memory effects may provide a more useful alternans criterion.     

Spatial heterogeneity of the restitution portrait in rabbit epicardium

Spatial heterogeneity of repolarization can provide a substrate for reentry to occur in myocardium. This heterogeneity may result from spatial differences in action potential duration (APD) restitution. The restitution portrait (RP) measures many aspects of rate-dependent restitution: the dynamic restitution curve (RC), S1-S2 RC, and short-term memory response. We used the RP to characterize epicardial patterns of spatial heterogeneity of restitution that were repeatable across animals. New Zealand White rabbit ventricles were paced from the epicardial apex, midventricle, or base, and optical action potentials were recorded from the same three regions. A perturbed downsweep pacing protocol was applied that measured the RP over a range of cycle lengths from 1,000 to 140 ms. The time constant of short-term memory measured close to the stimulus was dependent on location. In the midventricle the mean time constant was 19.1 +/- 1.1 s, but it was 39% longer at the apex (P < 0.01) and 23% longer at the base (P = 0.03). The S1-S2 RC slope was dependent on pacing site (P = 0.015), with steeper slope when pacing from the apex than from the base. There were no significant repeatable spatial patterns in steady-state APD at all cycle lengths or in dynamic RC slope. These results indicate that transient patterns of epicardial heterogeneity of APD may occur after a change in pacing rate. Thus it may affect cardiac electrical stability at the onset of a tachycardia or during a series of ectopic beats. Differences in restitution with respect to pacing site suggest that vulnerability may be affected by the location of reentry or ectopic foci.     

Asymptotics of Conduction Velocity Restitution in Models of Electrical Excitation in the Heart

We extend a non-Tikhonov asymptotic embedding, proposed earlier, for calculation of conduction velocity restitution curves in ionic models of cardiac excitability. Conduction velocity restitution is the simplest non-trivial spatially extended problem in excitable media, and in the case of cardiac tissue it is an important tool for prediction of cardiac arrhythmias and fibrillation. An idealized conduction velocity restitution curve requires solving a non-linear eigenvalue problem with periodic boundary conditions, which in the cardiac case is very stiff and calls for the use of asymptotic methods. We compare asymptotics of restitution curves in four examples, two generic excitable media models, and two ionic cardiac models. The generic models include the classical FitzHugh–Nagumo model and its variation by Barkley. They are treated with standard singular perturbation techniques. The ionic models include a simplified “caricature” of Noble (J. Physiol. Lond. 160:317–352, 1962) model and Beeler and Reuter (J. Physiol. Lond. 268:177–210, 1977) model, which lead to non-Tikhonov problems where known asymptotic results do not apply. The Caricature Noble model is considered with particular care to demonstrate the well-posedness of the corresponding boundary-value problem. The developed method for calculation of conduction velocity restitution is then applied to the Beeler–Reuter model. We discuss new mathematical features appearing in cardiac ionic models and possible applications of the developed method.     

Dynamics of human atrial cell models: restitution, memory, and intracellular calcium dynamics in single cells

Mathematical models of cardiac cells have become important tools for investigating the electrophysiological properties and behavior of the heart. As the number of published models increases, it becomes more difficult to choose a model appropriate for the conditions to be studied, especially when multiple models describing the species and region of the heart of interest are available. In this paper, we will review and compare two detailed ionic models of human atrial myocytes, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models include the same transmembrane currents and are largely based on the same experimental data from human atrial cells, the two models exhibit vastly different properties, especially in their dynamical behavior, including restitution and memory effects. The CM produces pronounced rate adaptation of action potential duration (APD) with limited memory effects, while the NM exhibits strong rate dependence of resting membrane potential (RMP), limited APD restitution, and stronger memory, as well as delayed afterdepolarizations and auto-oscillatory behavior upon cessation of rapid pacing. Channel conductance modifications based on experimentally measured changes during atrial fibrillation modify rate adaptation and memory in both models, but do not change the primary rate-dependent properties of APD and RMP for the CM and NM, respectively. Two sets of proposed changes to the NM that yield a spike-and-dome action potential morphology qualitatively similar to the CM at slow pacing rates similarly do not change the underlying dynamics of the model. Moreover, interchanging the formulations of all transmembrane currents between the two models while leaving calcium handling and ionic concentrations intact indicates that the currents strongly influence memory and the rate adaptation of RMP, while intracellular calcium dynamics primarily determine APD rate adaptation. Our results suggest that differences in intracellular calcium handling between the two human atrial myocyte models are responsible for marked dynamical differences and may prevent reconciliation between the models by straightforward channel conductance modifications.     

A Quantitative Comparison of the Behavior of Human Ventricular Cardiac Electrophysiology Models in Tissue

Numerical integration of mathematical models of heart cell electrophysiology provides an important computational tool for studying cardiac arrhythmias, but the abundance of available models complicates selecting an appropriate model. We study the behavior of two recently published models of human ventricular action potentials, the Grandi-Pasqualini-Bers (GPB) and the O''Hara-Virág-Varró-Rudy (OVVR) models, and compare the results with four previously published models and with available experimental and clinical data. We find the shapes and durations of action potentials and calcium transients differ between the GPB and OVVR models, as do the magnitudes and rate-dependent properties of transmembrane currents and the calcium transient. Differences also occur in the steady-state and S1–S2 action potential duration and conduction velocity restitution curves, including a maximum conduction velocity for the OVVR model roughly half that of the GPB model and well below clinical values. Between single cells and tissue, both models exhibit differences in properties, including maximum upstroke velocity, action potential amplitude, and minimum diastolic interval. Compared to experimental data, action potential durations for the GPB and OVVR models agree fairly well (although OVVR epicardial action potentials are shorter), but maximum slopes of steady-state restitution curves are smaller. Although studies show alternans in normal hearts, it occurs only in the OVVR model, and only for a narrow range of cycle lengths. We find initiated spiral waves do not progress to sustained breakup for either model. The dominant spiral wave period of the GPB model falls within clinically relevant values for ventricular tachycardia (VT), but for the OVVR model, the dominant period is longer than periods associated with VT. Our results should facilitate choosing a model to match properties of interest in human cardiac tissue and to replicate arrhythmia behavior more closely. Furthermore, by indicating areas where existing models disagree, our findings suggest avenues for further experimental work.     

Action potential duration restitution portraits of mammalian ventricular myocytes: role of calcium current

Construction of the action potential duration (APD) restitution portrait allows visualization of multiple aspects of the dynamics of periodically paced myocytes at various basic cycle lengths (BCLs). For the first time, we obtained the restitution portrait of isolated rabbit and guinea pig cardiac ventricular myocytes and analyzed the time constant, tau, of APD accommodation and the slopes of different types of restitution curves, Sdyn and S12, measured at varying BCLs. Our results indicate that both tau and the individual slopes are species and pacing dependent. In contrast, the mutual relationship between slopes Sdyn and S12 does not depend on pacing history, being a generic feature of the species. In addition, the maximum slope S12, measured in the restitution portrait at the lowest BCL, predicts the onset of alternans. Further, we investigated the role of the L-type calcium current, ICa-L, in the restitution portrait. We found that ICa-L dramatically affects APD accommodation, as well as the individual slopes Sdyn and S12 measured in the restitution portrait. However, peak calcium current plays a role only at small values of BCL. In conclusion, the results demonstrate that the restitution portrait is a powerful technique to investigate restitution properties of periodically paced cardiac myocytes and the onset of alternans, in particular. Moreover, the data also show that ICa-L plays a crucial role in multiple aspects of cardiac dynamics measured through the restitution portrait.     

Effects of diacetyl monoxime and cytochalasin D on ventricular fibrillation in swine right ventricles

Whether or not the excitation-contraction (E-C) uncoupler diacetyl monoxime (DAM) and cytochalacin D (Cyto D) alter the ventricular fibrillation (VF) activation patterns is unclear. We recorded single cell action potentials and performed optical mapping in isolated perfused swine right ventricles (RV) at different concentrations of DAM and Cyto D. Increasing the concentration of DAM results in progressively shortened action potential duration (APD) measured to 90% repolarization, reduced the slope of the APD restitition curve, decreased Kolmogorov-Sinai entropy, and reduced the number of VF wave fronts. In all RVs, 15-20 mmol/l DAM converted VF to ventricular tachycardia (VT). The VF could be reinduced after the DAM was washed out. In comparison, Cyto D (10-40 micromol/l) has no effects on APD restitution curve or the dynamics of VF. The effects of DAM on VF are associated with a reduced number of wave fronts and dynamic complexities in VF. These results are compatible with the restitution hypothesis of VF and suggest that DAM may be unsuitable as an E-C uncoupler for optical mapping studies of VF in the swine RVs.     

Dynamical effects of diffusive cell coupling on cardiac excitation and propagation: a simulation study

Cell coupling is considered to be important for cardiac action potential propagation and arrhythmogenesis. We carried out computer simulations to investigate the effects of stimulation strength and cell-to-cell coupling on action potential duration (APD) restitution, APD alternans, and stability of reentry in models of isolated cell, one-dimensional cable, and two-dimensional tissue. Phase I formulation of the Luo and Rudy action potential model was used. We found that stronger stimulation resulted in a shallower APD restitution curve and onset of APD alternans at a faster pacing rate. Reducing diffusive coupling between cells prolonged APD. Weaker diffusive currents along the direction of propagation steepened APD restitution and caused APD alternans to occur at a slower pacing rate in tissue. Diffusive current due to curvature changed APD but had little effect on APD restitution slope and onset of instability. Heterogeneous cell coupling caused APD inhomogeneities in space. Reduction in coupling strength either uniformly or randomly had little effect on the rotation period and stability of a reentry, but random cell decoupling slowed the rotation period and, thus, stabilized the reentry, preventing it from breaking up into multiple waves. Therefore, in addition to its effects on action potential conduction velocity, diffusive cell coupling also affects APD in a rate-dependent manner, causes electrophysiological heterogeneities, and thus modulates the dynamics of cardiac excitation. These effects are brought about by the modulation of ionic current activation and inactivation.     

Effect of beta-adrenergic blockade on dynamic electrical restitution in vivo

The slope of the action potential duration (APD) restitution curve may be a significant determinant of the propensity to develop ventricular fibrillation, with steeper slopes associated with a more arrhythmogenic substrate. We hypothesized that one mechanism by which beta-blockers reduce sudden cardiac death is by flattening the APD restitution curve. Therefore, we investigated whether infusion of esmolol modulates the APD restitution curve in vivo. In 10 Yorkshire pigs, dynamic APD restitution curves were determined from measurements of APD at 90% repolarization with a monophasic action potential catheter positioned against the right ventricular septum during right ventricular apical pacing in the basal state and during infusion of esmolol. APD restitution curves were fitted to the three-parameter (a, b, c) exponential equation, APD = a.[1 - e((-b.DI))] + c, where DI is the diastolic interval. Esmolol decreased the maximal APD slope, 0.68 +/- 0.14 vs. 0.94 +/- 0.24 (baseline), P = 0.002, and flattened the APD restitution curve at shorter DIs, 75 and 100 ms (P < 0.05). To compare the slopes of the APD restitution curves at similar steady states, slopes were also computed at points of intersection between the restitution curve and the lines representing pacing at a fixed cycle length (CL) of 200, 225, 250, 275, and 300 ms using the relationship CL = APD + DI. Esmolol decreased APD restitution slopes at CLs 200-275 ms (P < 0.05). Esmolol flattens the cardiac APD restitution curve in vivo, particularly at shorter CLs and DIs. This may represent a novel mechanism by which beta-blockers prevent sudden cardiac death.     

Properties of two human atrial cell models in tissue: restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.     

A model for human ventricular tissue

The experimental and clinical possibilities for studying cardiac arrhythmias in human ventricular myocardium are very limited. Therefore, the use of alternative methods such as computer simulations is of great importance. In this article we introduce a mathematical model of the action potential of human ventricular cells that, while including a high level of electrophysiological detail, is computationally cost-effective enough to be applied in large-scale spatial simulations for the study of reentrant arrhythmias. The model is based on recent experimental data on most of the major ionic currents: the fast sodium, L-type calcium, transient outward, rapid and slow delayed rectifier, and inward rectifier currents. The model includes a basic calcium dynamics, allowing for the realistic modeling of calcium transients, calcium current inactivation, and the contraction staircase. We are able to reproduce human epicardial, endocardial, and M cell action potentials and show that differences can be explained by differences in the transient outward and slow delayed rectifier currents. Our model reproduces the experimentally observed data on action potential duration restitution, which is an important characteristic for reentrant arrhythmias. The conduction velocity restitution of our model is broader than in other models and agrees better with available data. Finally, we model the dynamics of spiral wave rotation in a two-dimensional sheet of human ventricular tissue and show that the spiral wave follows a complex meandering pattern and has a period of 265 ms. We conclude that the proposed model reproduces a variety of electrophysiological behaviors and provides a basis for studies of reentrant arrhythmias in human ventricular tissue.     

One-Dimensional Mathematical Model of the Atrioventricular Node Including Atrio-Nodal, Nodal, and Nodal-His Cells

Mathematical models are a repository of knowledge as well as research and teaching tools. Although action potential models have been developed for most regions of the heart, there is no model for the atrioventricular node (AVN). We have developed action potential models for single atrio-nodal, nodal, and nodal-His cells. The models have the same action potential shapes and refractoriness as observed in experiments. Using these models, together with models for the sinoatrial node (SAN) and atrial muscle, we have developed a one-dimensional (1D) multicellular model including the SAN and AVN. The multicellular model has slow and fast pathways into the AVN and using it we have analyzed the rich behavior of the AVN. Under normal conditions, action potentials were initiated in the SAN center and then propagated through the atrium and AVN. The relationship between the AVN conduction time and the timing of a premature stimulus (conduction curve) is consistent with experimental data. After premature stimulation, atrioventricular nodal reentry could occur. After slow pathway ablation or block of the L-type Ca²⁺ current, atrioventricular nodal reentry was abolished. During atrial fibrillation, the AVN limited the number of action potentials transmitted to the ventricle. In the absence of SAN pacemaking, the inferior nodal extension acted as the pacemaker. In conclusion, we have developed what we believe is the first detailed mathematical model of the AVN and it shows the typical physiological and pathophysiological characteristics of the tissue. The model can be used as a tool to analyze the complex structure and behavior of the AVN.     

Short-term cardiac memory and mother rotor fibrillation

Short-term cardiac memory refers to the effects of pacing history on action potential duration (APD). Although the ionic mechanisms for short-term memory occurring over many heartbeats (also called APD accommodation) are poorly understood, they may have important effects on reentry and fibrillation. To explore this issue, we incorporated a generic memory current into the Phase I Luo and Rudy action potential model, which lacks short-term memory. The properties of this current were matched to simulate quantitatively human ventricular monophasic action potential accommodation. We show that, theoretically, short-term memory can resolve the paradox of how mother rotor fibrillation is initiated in heterogeneous tissue by physiological pacing. In simulated heterogeneous two-dimensional tissue and three-dimensional ventricles containing an inward rectifier K(+) current gradient, short-term memory could spontaneously convert multiple wavelet fibrillation to mother rotor fibrillation or to a mixture of both fibrillation types. This was due to progressive acceleration and stabilization of rotors as accumulation of memory shortened APD and flattened APD restitution slope nonuniformly throughout the tissue.     

Mechanisms of the acute ischemia-induced arrhythmogenesis--a simulation study

The underlying ionic mechanisms of ischemic-induced arrhythmia were studied by the computer simulation method. To approximate the real situation, ischemic cells were simulated by considering the three major component conditions of acute ischemia (elevated extracellular K(+) concentration, acidosis and anoxia) at the level of ionic currents and ionic concentrations, and a round ischemic zone was introduced into a homogeneous healthy sheet to avoid sharp angle of the ischemic tissue. The constructed models were solved using the operator splitting and adaptive time step methods, and the perturbation finite difference (PFD) scheme was first used to integrate the partial differential equations (PDEs) in the model. The numerical experiments showed that the action potential durations (APDs) of ischemic cells did not exhibited rate adaptation characteristic, resulting in flattening of the APD restitution curve. With reduction of sodium channel availability and long recovery of excitability, refractory period of the ischemic tissue was significantly prolonged, and could no longer be considered as same as APD. Slope of the conduction velocity (CV) restitution curve increased both in normal and ischemic region when pacing cycle length (PCL) was short, and refractory period dispersion increased with shortening of PCL as well. Therefore, dynamic changes of CV and dispersion of refractory period rather than APD were suggested to be the fundamental mechanisms of arrhythmia in regional ischemic myocardium.     

Toward prediction of the local onset of alternans in the heart

A beat-to-beat variation in the cardiac action potential duration is a phenomenon known as alternans. Alternans has been linked to ventricular fibrillation, and thus the ability to predict the onset of alternans could be clinically beneficial. Theoretically, it has been proposed that the slope of a restitution curve, which relates the duration of the action potential to the preceding diastolic interval, can predict the onset of alternans. Experimentally, however, this hypothesis has not been consistently proven, mainly because of the intrinsic complexity of the dynamics of cardiac tissue. It was recently shown that the restitution portrait, which combines several restitution curves simultaneously, is associated with the onset of alternans in isolated myocytes. Our main purpose in this study was to determine whether the restitution portrait is correlated with the onset of alternans in the heart, where the dynamics include a spatial complexity. We performed optical mapping experiments in isolated Langendorff-perfused rabbit hearts in which alternans was induced by periodic pacing at different frequencies, and identified the local onset of alternans, B^onset. We identified two regions of the heart: the area that exhibited alternans at B^onset (1:1_alt) and the area that did not (1:1). We constructed two-dimensional restitution portraits for the epicardial surface of the heart and measured the spatial distribution of three different slopes (the dynamic restitution slope, , and two local S1-S2 slopes, S₁₂ and ) separately for these two regions. We found that the S₁₂ and slopes differed significantly between the 1:1_alt and 1:1 regions just before the onset of alternans, and slopes were statistically similar. In addition, we found that the slopes of the dynamic restitution curve S_dyn were also statistically similar between these two regions. On the other hand, the quantitative values of all slopes were significantly different from the theoretically predicted value of one. These results demonstrate that the slopes measured in the restitution portrait correlate with the onset of alternans in the heart.     

Fitting Membrane Resistance along with Action Potential Shape in Cardiac Myocytes Improves Convergence: Application of a Multi-Objective Parallel Genetic Algorithm

Fitting parameter sets of non-linear equations in cardiac single cell ionic models to reproduce experimental behavior is a time consuming process. The standard procedure is to adjust maximum channel conductances in ionic models to reproduce action potentials (APs) recorded in isolated cells. However, vastly different sets of parameters can produce similar APs. Furthermore, even with an excellent AP match in case of single cell, tissue behaviour may be very different. We hypothesize that this uncertainty can be reduced by additionally fitting membrane resistance (R_m). To investigate the importance of R_m, we developed a genetic algorithm approach which incorporated R_m data calculated at a few points in the cycle, in addition to AP morphology. Performance was compared to a genetic algorithm using only AP morphology data. The optimal parameter sets and goodness of fit as computed by the different methods were compared. First, we fit an ionic model to itself, starting from a random parameter set. Next, we fit the AP of one ionic model to that of another. Finally, we fit an ionic model to experimentally recorded rabbit action potentials. Adding the extra objective (R_m, at a few voltages) to the AP fit, lead to much better convergence. Typically, a smaller MSE (mean square error, defined as the average of the squared error between the target AP and AP that is to be fitted) was achieved in one fifth of the number of generations compared to using only AP data. Importantly, the variability in fit parameters was also greatly reduced, with many parameters showing an order of magnitude decrease in variability. Adding R_m to the objective function improves the robustness of fitting, better preserving tissue level behavior, and should be incorporated.     

Simple experiments to understand the ionic origins and characteristics of the ventricular cardiac action potential

Electrophysiological experiments are helpful for students to understand the role of electrical activity in heart function. Papillary muscle, which belongs to the ventricle, offers the advantage of being easily studied using glass microelectrodes. In addition, there is commercially available software that simulates ventricular electrical activity and can help overcome some difficulties, such as voltage clamp experiments, which need expensive apparatus when used for studies on living preparations. Here, we present a class practical session that is taken by undergraduate students at our University. In the first part of this class, students record action potentials from papillary muscles with the use of glass microelectrodes, and they change extracellular conditions to study the ionic basis of the action potential. In the second part of the class, students simulate action potentials using the Oxsoft Heart model (v. 4.0) and model their previous experiments on papillary muscle to quantify the effects. In particular, the model is very helpful in promoting understanding of the effect that extracellular potassium has on cardiac action potential by simulating voltage clamp experiments. This twin approach of papillary muscle experiments and computer modeling leads to a good understanding of the functioning of the action potential and can help introduce discussion of some abnormal cardiac functioning.     

Cardiac dynamics: a simplified model for action potential propagation

This paper analyzes a new semiphysiological ionic model, used recently to study reexitations and reentry in cardiac tissue [I.R. Cantalapiedra et al, PRE 82 011907 (2010)]. The aim of the model is to reproduce action potencial morphologies and restitution curves obtained, either from experimental data, or from more complex electrophysiological models. The model divides all ion currents into four groups according to their function, thus resulting into fast-slow and inward-outward currents. We show that this simplified model is flexible enough as to accurately capture the electrical properties of cardiac myocytes, having the advantage of being less computational demanding than detailed electrophysiological models. Under some conditions, it has been shown to be amenable to mathematical analysis. The model reproduces the action potential (AP) change with stimulation rate observed both experimentally and in realistic models of healthy human and guinea pig myocytes (TNNP and LRd models, respectively). When simulated in a cable it also gives the right dependence of the conduction velocity (CV) with stimulation rate. Besides reproducing correctly these restitution properties, it also gives a good fit for the morphology of the AP, including the notch typical of phase 1. Finally, we perform simulations in a realistic geometric model of the rabbit’s ventricles, finding a good qualitative agreement in AP propagation and the ECG. Thus, this simplified model represents an alternative to more complex models when studying instabilities in wave propagation.