Simulation and prediction of functional block in the presence of structural and ionic heterogeneity

Amiodarone prevents wave front-tail interactions in patients with heart failure: an in silico study

Amiodarone (AM) is an antiarrhythmic drug whose chronic use has proved effective in preventing ventricular arrhythmias in a variety of patient populations, including those with heart failure (HF). AM has both class III [i.e., it prolongs the action potential duration (APD) via blocking potassium channels) and class I (i.e., it affects the rapid sodium channel) properties; however, the specific mechanism(s) by which it prevents reentry formation in patients with HF remains unknown. We tested the hypothesis that AM prevents reentry induction in HF during programmed electrical stimulation (PES) via its ability to induce postrepolarization refractoriness (PRR) via its class I effects on sodium channels. Here we extend our previous human action potential model to represent the effects of both HF and AM separately by calibrating to human tissue and clinical PES data, respectively. We then combine these models (HF + AM) to test our hypothesis. Results from simulations in cells and cables suggest that AM acts to increase PRR and decrease the elevation of takeoff potential. The ability of AM to prevent reentry was studied in silico in two-dimensional sheets in which a variety of APD gradients (ΔAPD) were imposed. Reentrant activity was induced in all HF simulations but was prevented in 23 of 24 HF + AM models. Eliminating the AM-induced slowing of the recovery of inactivation of the sodium channel restored the ability to induce reentry. In conclusion, in silico testing suggests that chronic AM treatment prevents reentry induction in patients with HF during PES via its class I effect to induce PRR.NEW & NOTEWORTHY This work presents a new model of the action potential of the human, which reproduces the complex dynamics during premature stimulation in heart failure patients with and without amiodarone. A specific mechanism of the ability of amiodarone to prevent reentrant arrhythmias is presented.

Discordant Alternans as a Mechanism for Initiation of Ventricular Fibrillation In Vitro

Background Ventricular tachyarrhythmias are often preceded by short sequences of premature ventricular complexes. In a previous study, a restitution-based computational model predicted which sequences of stimulated premature complexes were most likely to induce ventricular fibrillation in canines in vivo. However, the underlying mechanism, based on discordant-alternans dynamics, could not be verified in that study. The current study seeks to elucidate the mechanism by determining whether the spatiotemporal evolution of action potentials and initiation of ventricular fibrillation in in vitro experiments are consistent with model predictions. Methods and Results Optical mapping voltage signals from canine right-ventricular tissue (n=9) were obtained simultaneously from the entire epicardium and endocardium during and after premature stimulus sequences. Model predictions of action potential propagation along a 1-dimensional cable were developed using action potential duration versus diastolic interval data. The model predicted sign-change patterns in action potential duration and diastolic interval spatial gradients with posterior probabilities of 91.1%, and 82.1%, respectively. The model predicted conduction block with 64% sensitivity and 100% specificity. A generalized estimating equation logistic-regression approach showed that model-prediction effects were significant for both conduction block ( P<1×10-15, coefficient 44.36) and sustained ventricular fibrillation ( P=0.0046, coefficient, 1.63) events. Conclusions The observed sign-change patterns favored discordant alternans, and the model successfully identified sequences of premature stimuli that induced conduction block. This suggests that the relatively simple discordant-alternans-based process that led to block in the model may often be responsible for ventricular fibrillation onset when preceded by premature beats. These observations may aid in developing improved methods for anticipating block and ventricular fibrillation.

The Major Role of IK1 in Mechanisms of Rotor Drift in the Atria: A Computational Study

Maintenance of paroxysmal atrial fibrillation (AF) by fast rotors in the left atrium (LA) or at the pulmonary veins (PVs) is not fully understood. This review describes the role of the heterogeneous distribution of transmembrane currents in the PVs and LA junction (PV-LAJ) in the localization of rotors in the PVs. Experimentally observed heterogeneities in IK1, IKs, IKr, Ito, and ICaL in the PV-LAJ were incorporated into models of human atrial kinetics to simulate various conditions and investigate rotor drifting mechanisms. Spatial gradients in the currents resulted in shorter action potential duration, less negative minimum diastolic potential, slower upstroke and conduction velocity for rotors in the PV region than in the LA. Rotors under such conditions drifted toward the PV and stabilized at the less excitable region. Our simulations suggest that IK1 heterogeneity is dominant in determining the drift direction through its impact on the excitability gradient. These results provide a novel framework for understanding the complex dynamics of rotors in AF.

Fitting local repolarization parameters in cardiac reaction-diffusion models in the presence of electrotonic coupling

BACKGROUND

Repolarization gradients contribute to arrhythmogenicity. In reaction-diffusion models of cardiac tissue, heterogeneities in action potential duration (APD) can be created by locally modifying an intrinsic membrane kinetics parameter. Electrotonic coupling, however, acts as a confounding factor that modulates APD dispersion.

METHOD

We developed an algorithm based on a quasi-Newton method that iteratively adjusts the spatial distribution of a membrane parameter to reproduce a pre-defined target APD map in a coupled tissue. The method assumes that the relation between the adjustable parameter and APD is bijective in an isolated cell. Each iteration of the algorithm involved simulating the cardiac reaction-diffusion system with the updated parameter profile for one beat and extracting the APD map. The algorithm was extended to simultaneous estimation of two parameter profiles based on two APD maps at different repolarization thresholds.

RESULTS

The method was validated in 1D, 2D and 3D atrial tissues using synthetic target APD maps with controllable total variation and maximum APD gradient. The adjustable parameter was local acetylcholine concentration. The iterations converged provided that APD gradients were not too steep. Convergence was found to be faster 2-5 iterations) when the maximal gradient was less steep, when APD range was smaller and when tissue conductivity was reduced.

CONCLUSION

This algorithm provides a tool to automatically generate arrhythmogenic substrates with controllable repolarization gradients and possibly incorporate experimental APD maps into computer models.

Post-repolarization refractoriness increases vulnerability to block and initiation of reentrant impulses in heterogeneous infarcted myocardium

Myocardial infarction causes remodeling of the tissue structure and the density and kinetics of several ion channels in the cell membrane. Heterogeneities in refractory period (ERP) have been shown to occur in the infarct border zone and have been proposed to lead to initiation of arrhythmias. The purpose of this study is to quantify the window of vulnerability (WV) to block and initiation of reentrant impulses in myocardium with ERP heterogeneities using computer simulations. We found that ERP transitions at the border between normal ventricular cells (NZ) with different ERPs are smooth, whereas ERP transitions between NZ and infarct border zone cells (IZ) are abrupt. The profile of the ERP transitions is a combination of electrotonic interaction between NZ and IZ cells and the characteristic post-repolarization refractoriness (PRR) of IZ cells. ERP heterogeneities between NZ and IZ cells are more vulnerable to block and initiation of reentrant impulses than ERP heterogeneities between NZ cells. The relationship between coupling intervals of premature impulses (V1V2) and coupling intervals between premature and first reentrant impulses (V2T1) at NZ/NZ and NZ/IZ borders is inverse (i.e. the longer the coupling intervals of premature impulses the shorter the coupling interval between the premature and first reentrant impulses); this is in contrast with the reported V1V2/V2T1 relationship measured during initiation of reentrant impulses in canine infarcted hearts which is direct.

IN CONCLUSION

(1) ERP transitions at the NZ-IZ border are abrupt as a consequence of PRR; (2) PRR increases the vulnerability to block and initiation of reentrant impulses in heterogeneous myocardium; (3) V1V2/V2T1 relationships measured at ERP heterogeneities in the computer model and in experimental canine infarcts are not consistent. Therefore, it is likely that other mechanisms like micro and/or macro structural heterogeneities also contribute to initiation of reentrant impulses in infarcted hearts.

Dynamics of propagation of premature impulses in structurally remodeled infarcted myocardium: a computational analysis

Initiation of cardiac arrhythmias typically follows one or more premature impulses either occurring spontaneously or applied externally. In this study, we characterize the dynamics of propagation of single (S2) and double premature impulses (S3), and the mechanisms of block of premature impulses at structural heterogeneities caused by remodeling of gap junctional conductance (Gj) in infarcted myocardium. Using a sub-cellular computer model of infarcted tissue, we found that |I_Na,max|, prematurity (coupling interval with the previous impulse), and conduction velocity (CV) of premature impulses change dynamically as they propagate away from the site of initiation. There are fundamental differences between the dynamics of propagation of S2 and S3 premature impulses: for S2 impulses |I_Na,max| recovers fast, prematurity decreases and CV increases as propagation proceeds; for S3 impulses low values of |I_Na,max| persist, prematurity could increase, and CV could decrease as impulses propagate away from the site of initiation. As a consequence it is more likely that S3 impulses block at sites of structural heterogeneities causing source/sink mismatch than S2 impulses block. Whether premature impulses block at Gj heterogeneities or not is also determined by the values of Gj (and the space constant λ) in the regions proximal and distal to the heterogeneity: when λ in the direction of propagation increases >40%, premature impulses could block. The maximum slope of CV restitution curves for S2 impulses is larger than for S3 impulses. In conclusion: (1) The dynamics of propagation of premature impulses make more likely that S3 impulses block at sites of structural heterogeneities than S2 impulses block; (2) Structural heterogeneities causing an increase in λ (or CV) of >40% could result in block of premature impulses; (3) A decrease in the maximum slope of CV restitution curves of propagating premature impulses is indicative of an increased potential for block at structural heterogeneities.

Attraction of rotors to the pulmonary veins in paroxysmal atrial fibrillation: a modeling study

Maintenance of paroxysmal atrial fibrillation (AF) by fast rotors in the left atrium (LA) or at the pulmonary veins (PVs) is not fully understood. To gain insight into this dynamic and complex process, we studied the role of the heterogeneous distribution of transmembrane currents in the PVs and LA junction (PV-LAJ) in the localization of rotors in the PVs. We also investigated whether simple pacing protocols could be used to predict rotor drift in the PV-LAJ. Experimentally observed heterogeneities in IK1, IKs, IKr, Ito, and ICaL in the PV-LAJ were incorporated into two- and pseudo three-dimensional models of Courtemanche-Ramirez-Nattel-Kneller human atrial kinetics to simulate various conditions and investigate rotor drifting mechanisms. Spatial gradients in the currents resulted in shorter action potential duration, minimum diastolic potential that was less negative, and slower upstroke and conduction velocity for rotors in the PV region than in the LA. Rotors under such conditions drifted toward the PV and stabilized at the shortest action potential duration and less-excitable region, consistent with drift direction under intercellular coupling heterogeneities and regardless of the geometrical constraint in the PVs. Simulations with various IK1 gradient conditions and current-voltage relationships substantiated its major role in the rotor drift. In our 1:1 pacing protocol, we found that among various action potential properties, only the minimum diastolic potential gradient was a rate-independent predictor of rotor drift direction. Consistent with experimental and clinical AF studies, simulations in an electrophysiologically heterogeneous model of the PV-LAJ showed rotor attraction toward the PV. Our simulations suggest that IK1 heterogeneity is dominant compared to other currents in determining the drift direction through its impact on the excitability gradient. These results provide a believed novel framework for understanding the complex dynamics of rotors in AF.

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.

Spatial profiles of electrical mismatch determine vulnerability to conduction failure across a host-donor cell interface

BACKGROUND

Electrophysiological mismatch between host cardiomyocytes and donor cells can directly affect the electrical safety of cardiac cell therapies; however, the ability to study host-donor interactions at the microscopic scale in situ is severely limited. We systematically explored how action potential (AP) differences between cardiomyocytes and other excitable cells modulate vulnerability to conduction failure in vitro.

METHODS AND RESULTS

AP propagation was optically mapped at 75 μm resolution in micropatterned strands (n=152) in which host neonatal rat ventricular myocytes (AP duration=153.2±2.3 ms, conduction velocity=22.3±0.3 cm/s) seamlessly interfaced with genetically engineered excitable donor cells expressing inward rectifier potassium (Kir2.1) and cardiac sodium (Na(v)1.5) channels with either weak (conduction velocity=3.1±0.1 cm/s) or strong (conduction velocity=22.1±0.4 cm/s) electrical coupling. Selective prolongation of engineered donor cell AP duration (31.9-139.1 ms) by low-dose BaCl2 generated a wide range of host-donor repolarization time (RT) profiles with maximum gradients (∇RT(max)) of 5.5 to 257 ms/mm. During programmed stimulation of donor cells, the vulnerable time window for conduction block across the host-donor interface most strongly correlated with ∇RT(max). Compared with well-coupled donor cells, the interface composed of poorly coupled cells significantly shortened the RT profile width by 19.7% and increased ∇RT(max) and vulnerable time window by 22.2% and 19%, respectively. Flattening the RT profile by perfusion of 50 μmol/L BaCl2 eliminated coupling-induced differences in vulnerability to block.

CONCLUSIONS

Our results quantify how the degree of electrical mismatch across a cardiomyocyte-donor cell interface affects vulnerability to conduction block, with important implications for the design of safe cardiac cell and gene therapies.

Mechanisms of ventricular arrhythmias: a dynamical systems-based perspective

Defining the cellular electrophysiological mechanisms for ventricular tachyarrhythmias is difficult, given the wide array of potential mechanisms, ranging from abnormal automaticity to various types of reentry and kk activity. The degree of difficulty is increased further by the fact that any particular mechanism may be influenced by the evolving ionic and anatomic environments associated with many forms of heart disease. Consequently, static measures of a single electrophysiological characteristic are unlikely to be useful in establishing mechanisms. Rather, the dynamics of the electrophysiological triggers and substrates that predispose to arrhythmia development need to be considered. Moreover, the dynamics need to be considered in the context of a system, one that displays certain predictable behaviors, but also one that may contain seemingly stochastic elements. It also is essential to recognize that even the predictable behaviors of this complex nonlinear system are subject to small changes in the state of the system at any given time. Here we briefly review some of the short-, medium-, and long-term alterations of the electrophysiological substrate that accompany myocardial disease and their potential impact on the initiation and maintenance of ventricular arrhythmias. We also provide examples of cases in which small changes in the electrophysiological substrate can result in rather large differences in arrhythmia outcome. These results suggest that an interrogation of cardiac electrical dynamics is required to provide a meaningful assessment of the immediate risk for arrhythmia development and for evaluating the effects of putative antiarrhythmic interventions.

Conduction block in micropatterned cardiomyocyte cultures replicating the structure of ventricular cross-sections

AIMS

Structural and functional heterogeneities in cardiac tissue have been implicated in conduction block and arrhythmogenesis. However, the propensity of specific sites within the heart to initiate conduction block has not been systematically explored. We utilized cardiomyocyte cultures replicating the realistic, magnetic resonance imaging-measured tissue boundaries and fibre directions of ventricular cross-sections to investigate their roles in the development of conduction block.

METHODS AND RESULTS

The Sprague-Dawley neonatal rat cardiomyocytes were micropatterned to obtain cultures with realistic ventricular tissue boundaries and either random or realistic fibre directions. Rapid pacing was applied at multiple sites, with action potential propagation optically mapped. Excitation either failed at the stimulus site or conduction block developed remotely, often initiating reentry. The incidence of conduction block in isotropic monolayers (0% of cultures) increased with the inclusion of realistic tissue boundaries (17%) and further with realistic fibre directions (34%). Conduction block incidence was stimulus site-dependent and highest (77%) with rapid pacing from the right ventricular (RV) free wall. Furthermore, conduction block occurred exclusively at the insertion of the RV free wall into the septum, where structure-mediated current source-load mismatches acutely reduced wavefront and waveback velocity. Tissue boundaries and sharp gradients in fibre direction uniquely determined the evolution, shape, and position of conduction block lines.

CONCLUSION

Our study suggests that specific micro- and macrostructural features of the ventricle determine the incidence and spatiotemporal characteristics of conduction block, independent of spatial heterogeneities in ion channel expression.

Effect of resistive barrier location on the relationship between T-wave alternans and cellular repolarization alternans: a 1-D modeling study

Structural inhomogeneities in cardiac tissue have been associated with increased cellular repolarization alternans in animal experiments and increased T-wave alternans (TWA) in clinical studies. However, the effect of structural inhomogeneities on the relationship between cellular alternans and TWA has not been thoroughly investigated. We created 1-dimensional multicellular fiber models with and without a resistive barrier in various fiber regions and paced each model to induce cellular alternans. The models demonstrate that a resistive barrier in one fiber region substantially alters cellular repolarization alternans throughout the fiber. A midmyocardial or subepicardial barrier increase both TWA amplitude and maximum cellular alternans magnitude, relative to a fiber without a barrier. In addition, a direct relationship exists between TWA amplitude and maximum cellular alternans magnitude, which was highly dependent on barrier location. These results suggest that the position of a structural inhomogeneity within the myocardium may have substantial effects on dynamic repolarization instability and arrhythmogenicity.

Transmural ultrasound-based visualization of patterns of action potential wave propagation in cardiac tissue

The pattern of action potential propagation during various tachyarrhythmias is strongly suspected to be composed of multiple re-entrant waves, but has never been imaged in detail deep within myocardial tissue. An understanding of the nature and dynamics of these waves is important in the development of appropriate electrical or pharmacological treatments for these pathological conditions. We propose a new imaging modality that uses ultrasound to visualize the patterns of propagation of these waves through the mechanical deformations they induce. The new method would have the distinct advantage of being able to visualize these waves deep within cardiac tissue. In this article, we describe one step that would be necessary in this imaging process—the conversion of these deformations into the action potential induced active stresses that produced them. We demonstrate that, because the active stress induced by an action potential is, to a good approximation, only nonzero along the local fiber direction, the problem in our case is actually overdetermined, allowing us to obtain a complete solution. Use of two- rather than three-dimensional displacement data, noise in these displacements, and/or errors in the measurements of the fiber orientations all produce substantial but acceptable errors in the solution. We conclude that the reconstruction of action potential-induced active stress from the deformation it causes appears possible, and that, therefore, the path is open to the development of the new imaging modality.

Cardiac electrical dynamics: maximizing dynamical heterogeneity

The relationships between key features of the cardiac electrical activity, such as electrical restitution, discordant alternans, wavebreak, and reentry, and the onset of ventricular tachyarrhythmias have been characterized extensively under the condition of constant rapid pacing. However, it is unlikely that this scenario applies directly to the clinical situation, where the induction of ventricular tachycardia (VT) typically is associated with the interruption of normal cardiac rhythm by several premature beats. To address this issue, we have developed a general theory to explain why specific patterns of premature stimuli increase dynamic heterogeneity of repolarization and precipitate conduction block. The theory predicts that conduction block is caused by (1) creation of a spatial gradient in diastolic interval (DI) by waves traveling at slightly different velocities (ie, conduction velocity dispersion) and (2) amplification of the spatial gradient in DI over subsequent action potentials, secondary to a strong dependence of action potential duration on the preceding DI (ie, a steep action potential duration restitution function). Tests of this theory have been conducted in computer models of homogeneous tissue, where increased spatial dispersion of repolarization during premature stimulation can be attributed solely to the development of dynamical heterogeneity, and in a canine model exhibiting spontaneously occurring VT and sudden death. Our results thus far indicate that the probability of inducing ventricular fibrillation (VF) in the animal model is highest for those sequences predicted to cause conduction block in the computer model. An understanding of the mechanisms underlying these observations will help to identify key electrical phenomena in the onset of VT and fibrillation. Drug and electrical therapies can then be improved by targeting these specific phenomena.

Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence

Ventricular fibrillation is a lethal arrhythmia characterized by multiple wavelets usually starting from a single or figure-of-eight re-entrant circuit. Understanding the factors regulating vulnerability to the re-entry is essential for developing effective therapeutic strategies to prevent ventricular fibrillation. In this study, we investigated how pre-existing tissue heterogeneities and electrical restitution properties affect the initiation of re-entry by premature extrastimuli in two-dimensional cardiac tissue models. We studied two pacing protocols for inducing re-entry following the "sinus" rhythm (S1) beat: (1) a single premature (S2) extrastimulus in heterogeneous tissue; (2) two premature extrastimuli (S2 and S3) in homogeneous tissue. In the first case, the vulnerable window of re-entry is determined by the spatial dimension and extent of the heterogeneity, and is also affected by electrical restitution properties and the location of the premature stimulus. The vulnerable window first increases as the action potential duration (APD) difference between the inside and outside of the heterogeneous region increases, but then decreases as this difference increases further. Steeper APD restitution reduces the vulnerable window of re-entry. In the second case, electrical restitution plays an essential role. When APD restitution is flat, no re-entry can be induced. When APD restitution is steep, re-entry can be induced by an S3 over a range of S1S2 intervals, which is also affected by conduction velocity restitution. When APD restitution is even steeper, the vulnerable window is reduced due to collision of the spiral tips.

Action potential duration dispersion and alternans in simulated heterogeneous cardiac tissue with a structural barrier

Structural barriers to wave propagation in cardiac tissue are associated with a decreased threshold for repolarization alternans both experimentally and clinically. Using computer simulations, we investigated the effects of a structural barrier on the onset of spatially concordant and discordant alternans. We used two-dimensional tissue geometry with heterogeneity in selected potassium conductances to mimic known apex-base gradients. Although we found that the actual onset of alternans was similar with and without the structural barrier, the increase in alternans magnitude with faster pacing was steeper with the barrier--giving the appearance of an earlier alternans onset in its presence. This is consistent with both experimental structural barrier findings and the clinical observation of T-wave alternans occurring at slower pacing rates in patients with structural heart disease. In ionically homogeneous tissue, discordant alternans induced by the presence of the structural barrier arose at intermediate pacing rates due to a source-sink mismatch behind the barrier. In heterogeneous tissue, discordant alternans occurred during fast pacing due to a barrier-induced decoupling of tissue with different restitution properties. Our results demonstrate a causal relationship between the presence of a structural barrier and increased alternans magnitude and action potential duration dispersion, which may contribute to why patients with structural heart disease are at higher risk for ventricular tachyarrhythmias.

Vulnerable window for conduction block in a one-dimensional cable of cardiac cells, 1: single extrasystoles

Spatial dispersion of refractoriness, which is amplified by genetic diseases, drugs, and electrical and structural remodeling during heart disease, is recognized as a major factor increasing the risk of lethal arrhythmias and sudden cardiac death. Dispersion forms the substrate for unidirectional conduction block, which is required for the initiation of reentry by extrasystoles or rapid pacing. In this study, we examine theoretically and numerically how preexisting gradients in refractoriness control the vulnerable window for unidirectional conduction block by a single premature extrasystole. Using a kinematic model to represent wavefront-waveback interactions, we first analytically derived the relationship (under simplified conditions) between the vulnerable window and various electrophysiological parameters such as action potential duration gradients, refractoriness barriers, conduction velocity restitution, etc. We then compared these findings to numerical simulations using the kinematic model or the Luo-Rudy action potential model in a one-dimensional cable of cardiac cells. The results from all three methods agreed well. We show that a critical gradient in action potential duration for conduction block can be analytically derived, and once this critical gradient is exceeded, the vulnerable window increases proportionately with the refractory barrier and is modulated by conduction velocity restitution and gap junctional conductance. Moreover, the critical gradient for conduction block is higher for an extrasystole traveling in the opposite direction from the sinus beat than for one traveling in the same direction (e.g., an epicardial extrasystole versus an endocardial extrasystole).

Effects of rapid and slow potassium repolarization currents and calcium dynamics on hysteresis in restitution of action potential duration

We used a mathematical model to investigate effects of repolarizing currents I(kr) and I(ks), calcium (Ca) current I(CaL), and Ca dynamics in network sarcoplasmic reticulum and junctional sarcoplasmic reticulum (JSR) on hysteresis in restitution of action potential duration. Enhanced I(kr) increased slope of restitution, hysteresis loop thickness, and delay between peaks of diastolic intervals and action potential duration. Increase in I(ks) decreased loop thickness and peak delay. Decrease in I(CaL) had effects similar to increasing I(kr), except slope of restitution decreased markedly. Uptake of Ca into the network sarcoplasmic reticulum had less effect on hysteresis than transfer of Ca into JSR. Faster transfer of Ca into JSR markedly decreased loop thickness and peak delay. Our results provide insight into mechanisms responsible for this newly identified property of restitution. Such information will be valuable in studies where modification of hysteresis is used to investigate its role in arrhythmogenesis.

Whole heart action potential duration restitution properties in cardiac patients: a combined clinical and modelling study

Steep action potential duration (APD) restitution has been shown to facilitate wavebreak and ventricular fibrillation. The global APD restitution properties in cardiac patients are unknown. We report a combined clinical electrophysiology and computer modelling study to: (1) determine global APD restitution properties in cardiac patients; and (2) examine the interaction of the observed APD restitution with known arrhythmia mechanisms. In 14 patients aged 52-85 years undergoing routine cardiac surgery, 256 electrode epicardial mapping was performed. Activation-recovery intervals (ARI; a surrogate for APD) were recorded over the entire ventricular surface. Mono-exponential restitution curves were constructed for each electrode site using a standard S1-S2 pacing protocol. The median maximum restitution slope was 0.91, with 27% of all electrode sites with slopes<0.5, 29% between 0.5 and 1.0, and 20% between 1.0 and 1.5. Eleven per cent of restitution curves maintained slope>1 over a range of diastolic intervals of at least 30 ms; and 0.3% for at least 50 ms. Activation-recovery interval restitution was spatially heterogeneous, showing regional organization with multiple discrete areas of steep and shallow slope. We used a simplified computer model of 2-D cardiac tissue to investigate how heterogeneous APD restitution can influence vulnerability to, and stability of re-entry. Our model showed that heterogeneity of restitution can act as a potent arrhythmogenic substrate, as well as influencing the stability of re-entrant arrhythmias. Global epicardial mapping in humans showed that APD restitution slopes were organized into regions of shallow and steep slopes. This heterogeneous organization of restitution may provide a substrate for arrhythmia.

An integrative model of mouse cardiac electrophysiology from cell to torso

AIMS

Although the transgenic mouse has become an important new tool in the study of human diseases and the design of new therapies, a complete picture of cardiac electrophysiology in the mouse, from genome to body surface, is lacking. A computational model of the mouse heart is presented, which is used to study the impact of ion-channel and structural manipulations on the distributions of extracellular potentials on the heart and body surface.

METHODS

A model of the mouse heart anatomy, fibre organization and torso geometry was constructed from DTMRI images. An anisotropic bidomain model, with a modified Pandit et al. model for the ionic currents, was used to represent the electrical properties of the tissue. Spatial heterogeneity in the ion currents was introduced by modulating the transient outward current. A sinus beat was simulated in hearts with different tissue and membrane properties and the extracellular potentials were computed at both the heart and body surface.

RESULTS

The simulated transmembrane patterns in the heart, and the timing and morphology of the simulated ECG waveforms were consistent with experimental measurements. In addition, the patterns of activation and recovery and the waveforms of the corresponding ECG were found to be relatively insensitive to changes in cell type distribution and tissue anisotropy.

CONCLUSION

Because of the small size of the heart, an integrative model of mouse electrophysiology can be simulated from cell to torso, enabling a new tool to study how extracellular signals might be used to detect molecular changes underlying an arrhythmogenic substrate.

An ionic model of stretch-activated and stretch-modulated currents in rabbit ventricular myocytes

AIMS

To develop an ionic model of stretch-activated and stretch-modulated currents in rabbit ventricular myocytes consistent with experimental observations, that can be used to investigate the role of these currents in intact myocardium.

METHODS AND RESULTS

A non-specific cation-selective stretch-activated current I(ns), was incorporated into the Puglisi-Bers ionic model of epicardial, endocardial and midmyocardial ventricular myocytes. Using the model, we predict a reduction in action potential duration at 20% repolarization (APD(20)) and action potential amplitude, an elevated resting transmembrane potential and either an increase or decrease in APD(90), depending on the reversal potential of I(ns). A stretch-induced decrease in I(K1) (70%), plus a small I(ns) current (g(ns) = 10 pS), results in a reduction in APD(20) and increase in APD(90), and a reduced safety factor for conduction. Increasing I(K1) (150%) plus a large I(ns) current (g(ns) = 40 pS), also leads to a reduction in APD(20) and increase in APD(90), but with a greater safety factor. Endocardial and midmyocardial cells appear to be the most sensitive to stretch-induced changes in action potential. The addition of the K(+)-specific stretch-activated current (SAC) I(Ko) results in action potential shortening.

CONCLUSION

Transmural heterogeneity of I(Ko) may reduce repolarization gradients in intact myocardium caused by intrinsic ion channel densities, nonuniform strains and electrotonic effects.

Regional differences in APD restitution can initiate wavebreak and re-entry in cardiac tissue: a computational study

Background

Regional differences in action potential duration (APD) restitution in the heart favour arrhythmias, but the mechanism is not well understood.

Methods

We simulated a 150 × 150 mm 2D sheet of cardiac ventricular tissue using a simplified computational model. We investigated wavebreak and re-entry initiated by an S1S2S3 stimulus protocol in tissue sheets with two regions, each with different APD restitution. The two regions had a different APD at short diastolic interval (DI), but similar APD at long DI. Simulations were performed twice; once with both regions having steep (slope > 1), and once with both regions having flat (slope < 1) APD restitution.

Results

Wavebreak and re-entry were readily initiated using the S1S2S3 protocol in tissue sheets with two regions having different APD restitution properties. Initiation occurred irrespective of whether the APD restitution slopes were steep or flat. With steep APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms with S1S2 of 250 ms, to 75 ms (S1S2 180 ms). With flat APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms (S1S2 250 ms), to 21 ms (S1S2 340 ms) and then 11 ms (S1S2 400 ms).

Conclusion

Regional differences in APD restitution are an arrhythmogenic substrate that can be concealed at normal heart rates. A premature stimulus produces regional differences in repolarisation, and a further premature stimulus can then result in wavebreak and initiate re-entry. This mechanism for initiating re-entry is independent of the steepness of the APD restitution curve.

Head-tail interactions in numerical simulations of reentry in a ring of cardiac tissue

BACKGROUND

The relationships between action potential duration (APD) and conduction velocity (CV) to the previous diastolic interval (DI) are known as the APD and CV restitution relationships. There is considerable debate regarding the importance of these relationships in the development and stability of reentry.

OBJECTIVES

The purpose of this study was to increase the understanding of the ionic basis for restitution during reentry.

METHODS

APD and CV were studied numerically during one-dimensional reentry as ring length (L) was shortened. A three-state variable model (u, v, w) was used to analyze the effect of gating variables of the fast (v) and slow (w) currents on the spatial and temporal dynamics of transmembrane potential (u). Three parameter sets were used corresponding to three APD and CV restitution curves.

RESULTS

Sustained spatial oscillations of APD and CV larger than the ring length were observed in two of the parameter sets (cytochalasin-D model [CYTO] and model 3 [M3]) before block occurred at L = 6 cm. The last model (diacetyl monoxime [DAM]) resulted in uniform APD and CV for all L until block occurred at L = 3 cm. Multivalued APD and CV restitution relationships due to "dephasing" of w and v with DI were observed in M3 and CYTO simulations. Overall, these dynamics could be explained by the wavelength-to-ring length ratio and the sensitivity of APD on the value of the gating variables w and v.

CONCLUSION

Propagation stability is mostly controlled by APD sensitivity to w, but the APD restitution slope does not always reflect this sensitivity. The interaction of the dynamic history (i.e., memory) of the fast and slow currents and electrotonic effects resulted in multivalued restitution curves.

Head-tail interactions in numerical simulations of reentry in a ring of cardiac tissue

BACKGROUND

The relationships between action potential duration (APD) and conduction velocity (CV) to the previous diastolic interval (DI) are known as the APD and CV restitution relationships. There is considerable debate regarding the importance of these relationships in the development and stability of reentry.

OBJECTIVES

The purpose of this study was to increase the understanding of the ionic basis for restitution during reentry.

METHODS

APD and CV were studied numerically during one-dimensional reentry as ring length (L) was shortened. A three-state variable model (u, v, w) was used to analyze the effect of gating variables of the fast (v) and slow (w) currents on the spatial and temporal dynamics of transmembrane potential (u). Three parameter sets were used corresponding to three APD and CV restitution curves.

RESULTS

Sustained spatial oscillations of APD and CV larger than the ring length were observed in two of the parameter sets (cytochalasin-D model [CYTO] and model 3 [M3]) before block occurred at L = 6 cm. The last model (diacetyl monoxime [DAM]) resulted in uniform APD and CV for all L until block occurred at L = 3 cm. Multivalued APD and CV restitution relationships due to "dephasing" of w and v with DI were observed in M3 and CYTO simulations. Overall, these dynamics could be explained by the wavelength-to-ring length ratio and the sensitivity of APD on the value of the gating variables w and v.

CONCLUSION

Propagation stability is mostly controlled by APD sensitivity to w, but the APD restitution slope does not always reflect this sensitivity. The interaction of the dynamic history (i.e., memory) of the fast and slow currents and electrotonic effects resulted in multivalued restitution curves.

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.

Dispersion of cardiac action potential duration and the initiation of re-entry: a computational study

BACKGROUND

The initiation of re-entrant cardiac arrhythmias is associated with increased dispersion of repolarisation, but the details are difficult to investigate either experimentally or clinically. We used a computational model of cardiac tissue to study systematically the association between action potential duration (APD) dispersion and susceptibility to re-entry.

METHODS

We simulated a 60 x 60 mm2 D sheet of cardiac ventricular tissue using the Luo-Rudy phase 1 model, with maximal conductance of the K+ channel gKmax set to 0.004 mS mm(-2). Within the central 40 x 40 mm region we introduced square regions with prolonged APD by reducing gKmax to between 0.001 and 0.003 mS mm(-2). We varied (i) the spatial scale of these regions, (ii) the magnitude of gKmax in these regions, and (iii) cell-to-cell coupling.

RESULTS

Changing spatial scale from 5 to 20 mm increased APD dispersion from 49 to 102 ms, and the susceptible window from 31 to 86 ms. Decreasing gKmax in regions with prolonged APD from 0.003 to 0.001 mS mm-2 increased APD dispersion from 22 to 70 ms, and the susceptible window from <1 to 56 ms. Decreasing cell-to-cell coupling by changing the diffusion coefficient from 0.2 to 0.05 mm2 ms(-1) increased APD dispersion from 57 to 88 ms, and increased the susceptible window from 41 to 74 ms.

CONCLUSION

We found a close association between increased APD dispersion and susceptibility to re-entrant arrhythmias, when APD dispersion is increased by larger spatial scale of heterogeneity, greater electrophysiological heterogeneity, and weaker cell-to-cell coupling.

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.

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 Ca2+ handling and transmural heterogeneity in IKs and Ito. 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.

Effects of Na+ channel and cell coupling abnormalities on vulnerability to reentry: a simulation study

The role of dynamic instabilities in the initiation of reentry in diseased (remodeled) hearts remains poorly explored. Using computer simulations, we studied the effects of altered Na+ channel and cell coupling properties on the vulnerable window (VW) for reentry in simulated two-dimensional cardiac tissue with and without dynamic instabilities. We related the VW for reentry to effects on conduction velocity, action potential duration (APD), effective refractory period dispersion and restitution, and concordant and discordant APD alternans. We found the following: 1) reduced Na+ current density and slowed recovery promoted postrepolarization refractoriness and enhanced concordant and discordant APD alternans, which increased the VW for reentry; 2) uniformly reduced cell coupling had little effect on cellular electrophysiological properties and the VW for reentry. However, randomly reduced cell coupling combined with decoupling promoted APD dispersion and alternans, which also increased the VW for reentry; 3) the combination of decreased Na+ channel conductance, slowed Na+ channel recovery, and cellular uncoupling synergistically increased the VW for reentry; and 4) the VW for reentry was greater when APD restitution slope was steep than when it was flat. In summary, altered Na+ channel and cellular coupling properties increase vulnerability to reentrant arrhythmias. In remodeled hearts with altered Na+ channel properties and cellular uncoupling, dynamic instabilities arising from electrical restitution exert important influences on the VW for reentry.

Tissue discontinuities affect conduction velocity restitution: a mechanism by which structural barriers may promote wave break

BACKGROUND

The mechanism by which structural barriers promote wave break and fibrillation is unclear. Conduction velocity (CV) restitution is an important determinant of wave break. Abnormal CV restitution is associated with ventricular fibrillation in patients with heart disease and arises preferentially in fibrotic myocardium. We hypothesize that tissue discontinuities imposed by structural barriers cause abnormal CV restitution.

METHODS AND RESULTS

Tissue discontinuities were simulated in cultures of neonatal rat heart cells grown in 8-armed star patterns. Premature stimulation was applied at the extremity of 1 arm (n=12) while extracellular electrograms were recorded at 24 sites throughout the star. Action potentials were recorded at the following 3 sites: in the stimulated arm and at the discontinuity both proximal to and distal from the star center. Extracellular recordings revealed progressive increases in activation delay (indicative for abnormal CV restitution) only at the discontinuity from arms proximal to the star center. The mean increase in delay was 0.81+/-0.41 ms/10 ms for recording sites proximal to and 3.13+/-0.58 ms/10 ms for sites distal from this discontinuity. Depolarizing currents were determined in single cells during premature stimulation and for voltage configurations similar to those arising at the discontinuity. Both voltage-clamp measurements and computer simulations showed that delay at the discontinuity was associated with biphasic, prolonged activation and delayed inactivation of depolarizing current.

CONCLUSIONS

Tissue discontinuities cause abnormal CV restitution. Rapid increase in activation after an initial slow activation and delayed inactivation at the discontinuity lengthen the duration of depolarizing current and cause the abnormal restitution.

Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity

It has become widely accepted that the most dangerous cardiac arrhythmias are due to reentrant waves, i.e., electrical wave(s) that recirculate repeatedly throughout the tissue at a higher frequency than the waves produced by the heart's natural pacemaker (sinoatrial node). However, the complicated structure of cardiac tissue, as well as the complex ionic currents in the cell, have made it extremely difficult to pinpoint the detailed dynamics of these life-threatening reentrant arrhythmias. A simplified ionic model of the cardiac action potential (AP), which can be fitted to a wide variety of experimentally and numerically obtained mesoscopic characteristics of cardiac tissue such as AP shape and restitution of AP duration and conduction velocity, is used to explain many different mechanisms of spiral wave breakup which in principle can occur in cardiac tissue. Some, but not all, of these mechanisms have been observed before using other models; therefore, the purpose of this paper is to demonstrate them using just one framework model and to explain the different parameter regimes or physiological properties necessary for each mechanism (such as high or low excitability, corresponding to normal or ischemic tissue, spiral tip trajectory types, and tissue structures such as rotational anisotropy and periodic boundary conditions). Each mechanism is compared with data from other ionic models or experiments to illustrate that they are not model-specific phenomena. Movies showing all the breakup mechanisms are available at http://arrhythmia.hofstra.edu/breakup and at ftp://ftp.aip.org/epaps/chaos/E-CHAOEH-12-039203/ INDEX.html. The fact that many different breakup mechanisms exist has important implications for antiarrhythmic drug design and for comparisons of fibrillation experiments using different species, electromechanical uncoupling drugs, and initiation protocols. (c) 2002 American Institute of Physics.