Reentrant waves and their elimination in a model of mammalian ventricular tissue

Evaluating pro-arrhythmogenic effects of the T634S-hERG mutation: insights from a simulation study

A mutation to serine of a conserved threonine (T634S) in the hERG K+ channel S6 pore region has been identified as a variant of uncertain significance, showing a loss-of-function effect. However, its potential consequences for ventricular excitation and arrhythmogenesis have not been reported. This study evaluated possible functional effects of the T634S-hERG mutation on ventricular excitation and arrhythmogenesis by using multi-scale computer models of the human ventricle. A Markov chain model of the rapid delayed rectifier potassium current (IKr) was reconstructed for wild-type and T634S-hERG mutant conditions and incorporated into the ten Tusscher et al. models of human ventricles at cell and tissue (1D, 2D and 3D) levels. Possible functional impacts of the T634S-hERG mutation were evaluated by its effects on action potential durations (APDs) and their rate-dependence (APDr) at the cell level; and on the QT interval of pseudo-ECGs, tissue vulnerability to unidirectional conduction block (VW), spiral wave dynamics and repolarization dispersion at the tissue level. It was found that the T634S-hERG mutation prolonged cellular APDs, steepened APDr, prolonged the QT interval, increased VW, destablized re-entry and augmented repolarization dispersion across the ventricle. Collectively, these results imply potential pro-arrhythmic effects of the T634S-hERG mutation, consistent with LQT2.

Investigation of the Role of Myocyte Orientations in Cardiac Arrhythmia Using Image-Based Models

Cardiac electrical excitation-propagation is influenced by myocyte orientations (cellular organization). Quantitatively understanding this relationship presents a significant research challenge, especially during arrhythmias in which excitation patterns become complex. Tissue-scale simulations of cardiac electrophysiology, incorporating both dynamic action potential behavior and image-based myocardial architecture, provide an approach to investigate three-dimensional (3D) propagation of excitation waves in the heart. In this study, we aimed to assess the importance of natural variation in myocyte orientations on cardiac arrhythmogenesis using 3D tissue electrophysiology simulations. Three anatomical models (i.e., describing myocyte orientations) of healthy rat ventricles—obtained using diffusion tensor imaging at 100 μm resolution—were registered to a single biventricular geometry (i.e., a single cardiac shape), in which the myocyte orientations could be represented by each of the diffusion tensor imaging data sets or by an idealized rule-based description. The Fenton-Karma cellular excitation model was modified to reproduce rat ventricular action potential duration restitution to create reaction-diffusion cardiac electrophysiology models. Over 250 3D simulations were performed to investigate the effects of myocyte orientations on the following: 1) ventricular activation, 2) location-dependent arrhythmia induction via rapid pacing, and 3) dynamics of re-entry averaged over multiple episodes. It was shown that 1) myocyte orientation differences manifested themselves in local activation times, but the influence on total activation time was small; 2) differences in myocyte orientations could critically affect the inducibility and persistence of arrhythmias for specific stimulus-location/cycle-length combinations; and 3) myocyte orientations alone could be an important determinant of scroll wave break, although no significant differences were observed in averaged arrhythmia dynamics between the four myocyte orientation scenarios considered. Our results show that myocyte orientations are an important determinant of arrhythmia inducibility, persistence, and scroll wave break. These findings suggest that where specificity is desired (for example, when predicting location-dependent, patient-specific arrhythmia inducibility), subject-specific myocyte orientations may be important.

Arrhythmia mechanisms and spontaneous calcium release: Bi-directional coupling between re-entrant and focal excitation

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation.

In silico Assessment of Pharmacotherapy for Human Atrial Patho-Electrophysiology Associated With hERG-Linked Short QT Syndrome

Short QT syndrome variant 1 (SQT1) arises due to gain-of-function mutations to the human Ether-à-go-go-Related Gene (hERG), which encodes the α subunit of channels carrying rapid delayed rectifier potassium current, I_Kr. In addition to QT interval shortening and ventricular arrhythmias, SQT1 is associated with increased risk of atrial fibrillation (AF), which is often the only clinical presentation. However, the underlying basis of AF and its pharmacological treatment remain incompletely understood in the context of SQT1. In this study, computational modeling was used to investigate mechanisms of human atrial arrhythmogenesis consequent to a SQT1 mutation, as well as pharmacotherapeutic effects of selected class I drugs–disopyramide, quinidine, and propafenone. A Markov chain formulation describing wild type (WT) and N588K-hERG mutant I_Kr was incorporated into a contemporary human atrial action potential (AP) model, which was integrated into one-dimensional (1D) tissue strands, idealized 2D sheets, and a 3D heterogeneous, anatomical human atria model. Multi-channel pharmacological effects of disopyramide, quinidine, and propafenone, including binding kinetics for I_Kr/hERG and sodium current, I_Na, were considered. Heterozygous and homozygous formulations of the N588K-hERG mutation shortened the AP duration (APD) by 53 and 86 ms, respectively, which abbreviated the effective refractory period (ERP) and excitation wavelength in tissue, increasing the lifespan and dominant frequency (DF) of scroll waves in the 3D anatomical human atria. At the concentrations tested in this study, quinidine most effectively prolonged the APD and ERP in the setting of SQT1, followed by disopyramide and propafenone. In 2D simulations, disopyramide and quinidine promoted re-entry termination by increasing the re-entry wavelength, whereas propafenone induced secondary waves which destabilized the re-entrant circuit. In 3D simulations, the DF of re-entry was reduced in a dose-dependent manner for disopyramide and quinidine, and propafenone to a lesser extent. All of the anti-arrhythmic agents promoted pharmacological conversion, most frequently terminating re-entry in the order quinidine > propafenone = disopyramide. Our findings provide further insight into mechanisms of SQT1-related AF and a rational basis for the pursuit of combined I_Kr and I_Na block based pharmacological strategies in the treatment of SQT1-linked AF.

Human Atrial Arrhythmogenesis and Sinus Bradycardia in KCNQ1-Linked Short QT Syndrome: Insights From Computational Modelling

Atrial fibrillation (AF) and sinus bradycardia have been reported in patients with short QT syndrome variant 2 (SQT2), which is underlain by gain-of-function mutations in KCNQ1 encoding the α subunit of channels carrying slow delayed rectifier potassium current, I_Ks. However, the mechanism(s) underlying the increased atrial arrhythmogenesis and impaired cardiac pacemaking activity arising from increased I_Ks remain unclear. Possible pharmacological interventions of AF in the SQT2 condition also remain to be elucidated. Using computational modelling, we assessed the functional impact of SQT2 mutations on human sinoatrial node (SAN) pacemaking, atrial repolarisation and arrhythmogenesis, and efficacy of the anti-arrhythmic drug quinidine. Markov chain formulations of I_Ks describing two KCNQ1 mutations – V141M and V307L – were developed from voltage-clamp experimental data and then incorporated into contemporary action potential (AP) models of human atrial and SAN cells, the former of which were integrated into idealised and anatomically detailed tissue models. Both mutations shortened atrial AP duration (APD) through distinct I_Ks ‘gain-of-function’ mechanisms, whereas SAN pacemaking rate was slowed markedly only by the V141M mutation. Differences in APD restitution steepness influenced re-entry dynamics in tissue – the V141M mutation promoted stationary and stable spiral waves whereas the V307L mutation promoted non-stationary and unstable re-entrant waves. Both mutations shortened tissue excitation wavelength through reduced effective refractory period but not conduction velocity, which served to increase the lifespan of re-entrant excitation in a 3D anatomical human atria model, as well as the dominant frequency (DF), which was higher for the V141M mutation. Quinidine was effective at terminating arrhythmic excitation waves associated with the V307L but not V141M mutation, and reduced the DF in a dose-dependent manner under both mutation conditions. This study provides mechanistic insights into different AF/bradycardia phenotypes in SQT2 and the efficacy of quinidine pharmacotherapy.

Dispersion of Recovery and Vulnerability to Re-entry in a Model of Human Atrial Tissue With Simulated Diffuse and Focal Patterns of Fibrosis

Fibrosis in atrial tissue can act as a substrate for persistent atrial fibrillation, and can be focal or diffuse. Regions of fibrosis are associated with slowed or blocked conduction, and several approaches have been used to model these effects. In this study a computational model of 2D atrial tissue was used to investigate how the spatial scale of regions of simulated fibrosis influenced the dispersion of action potential duration (APD) and vulnerability to re-entry in simulated normal human atrial tissue, and human tissue that has undergone remodeling as a result of persistent atrial fibrillation. Electrical activity was simulated in a 10 × 10 cm square 2D domain, with a spatially varying diffusion coefficient as described below. Cellular electrophysiology was represented by the Courtemanche model for human atrial cells, with the model parameters set for normal and remodeled cells. The effect of fibrosis was modeled with a smoothly varying diffusion coefficient, obtained from sampling a Gaussian random field (GRF) with length scales of between 1.25 and 10.0 mm. Twenty samples were drawn from each field, and used to allocate a value of diffusion coefficient between 0.05 and 0.2 mm²/ms. Dispersion of APD was assessed in each sample by pacing at a cycle length of 1,000 ms, followed by a premature beat with a coupling interval of 400 ms. Vulnerability to re-entry was assessed with an aggressive pacing protocol with pacing cycle lengths decreasing from 450 to 250 ms in 25 ms intervals for normal tissue and 300–150 ms for remodeled tissue. Simulated fibrosis at smaller spatial scales tended to lengthen APD, increase APD dispersion, and increase vulnerability to sustained re-entry relative to fibrosis at larger spatial scales. This study shows that when fibrosis is represented by smoothly varying tissue diffusion, the spatial scale of fibrosis has important effects on both dispersion of recovery and vulnerability to re-entry.

Computational Analysis of the Mode of Action of Disopyramide and Quinidine on hERG-Linked Short QT Syndrome in Human Ventricles

The short QT syndrome (SQTS) is a rare cardiac disorder associated with arrhythmias and sudden death. Gain-of-function mutations to potassium channels mediating the rapid delayed rectifier current, I_Kr, underlie SQTS variant 1 (SQT1), in which treatment with Na⁺ and K⁺ channel blocking class Ia anti-arrhythmic agents has demonstrated some efficacy. This study used computational modeling to gain mechanistic insights into the actions of two such drugs, disopyramide and quinidine, in the setting of SQT1. The O'Hara-Rudy (ORd) human ventricle model was modified to incorporate a Markov chain formulation of I_Kr describing wild type (WT) and SQT1 mutant conditions. Effects of multi-channel block by disopyramide and quinidine, including binding kinetics and altered potency of I_Kr/hERG channel block in SQT1 and state-dependent block of sodium channels, were simulated on action potential and multicellular tissue models. A one-dimensional (1D) transmural ventricular strand model was used to assess prolongation of the QT interval, effective refractory period (ERP), and re-entry wavelength (WL) by both drugs. Dynamics of re-entrant excitation waves were investigated using a 3D human left ventricular wedge model. In the setting of SQT1, disopyramide, and quinidine both produced a dose-dependent prolongation in (i) the QT interval, which was primarily due to I_Kr block, and (ii) the ERP, which was mediated by a synergistic combination of I_Kr and I_Na block. Over the same range of concentrations quinidine was more effective in restoring the QT interval, due to more potent block of I_Kr. Both drugs demonstrated an anti-arrhythmic increase in the WL of re-entrant circuits. In the 3D wedge, disopyramide and quinidine at clinically-relevant concentrations decreased the dominant frequency of re-entrant excitations and exhibited anti-fibrillatory effects; preventing formation of multiple, chaotic wavelets which developed in SQT1, and could terminate arrhythmias. This computational modeling study provides novel insights into the clinical efficacy of disopyramide and quinidine in the setting of SQT1; it also dissects ionic mechanisms underlying QT and ERP prolongation. Our findings show that both drugs demonstrate efficacy in reversing the SQT1 phenotype, and indicate that disopyramide warrants further investigation as an alternative to quinidine in the treatment of SQT1.

Computational assessment of the functional role of sinoatrial node exit pathways in the human heart

AIM

The human right atrium and sinoatrial node (SAN) anatomy is complex. Optical mapping experiments suggest that the SAN is functionally insulated from atrial tissue except at discrete SAN-atrial electrical junctions called SAN exit pathways, SEPs. Additionally, histological imaging suggests the presence of a secondary pacemaker close to the SAN. We hypothesise that a) an insulating border-SEP anatomical configuration is related to SAN arrhythmia; and b) a secondary pacemaker, the paranodal area, is an alternate pacemaker but accentuates tachycardia. A 3D electro-anatomical computational model was used to test these hypotheses.

METHODS

A detailed 3D human SAN electro-anatomical mathematical model was developed based on our previous anatomical reconstruction. Electrical activity was simulated using tissue specific variants of the Fenton-Karma action potential equations. Simulation experiments were designed to deploy this complex electro-anatomical system to assess the roles of border-SEPs and paranodal area by mimicking experimentally observed SAN arrhythmia. Robust and accurate numerical algorithms were implemented for solving the mono domain reaction-diffusion equation implicitly, calculating 3D filament traces, and computing dominant frequency among other quantitative measurements.

RESULTS

A centre to periphery gradient of increasing diffusion was sufficient to permit initiation of pacemaking at the centre of the 3D SAN. Re-entry within the SAN, micro re-entry, was possible by imposing significant SAN fibrosis in the presence of the insulating border. SEPs promoted the micro re-entry to generate more complex SAN-atrial tachycardia. Simulation of macro re-entry, i.e. re-entry around the SAN, was possible by inclusion of atrial fibrosis in the presence of the insulating border. The border shielded the SAN from atrial tachycardia. However, SAN micro-structure intercellular gap junctional coupling and the paranodal area contributed to prolonged atrial fibrillation. Finally, the micro-structure was found to be sufficient to explain shifts of leading pacemaker site location.

CONCLUSIONS

The simulations establish a relationship between anatomy and SAN electrical function. Microstructure, in the form of intercellular gap junction coupling, was found to regulate SAN function and arrhythmia.

In silico assessment of genetic variation in KCNA5 reveals multiple mechanisms of human atrial arrhythmogenesis

A recent experimental study investigating patients with lone atrial fibrillation identified six novel mutations in the KCNA5 gene. The mutants exhibited both gain- and loss-of-function of the atrial specific ultra-rapid delayed rectifier K⁺ current, I_Kur. The aim of this study is to elucidate and quantify the functional impact of these KCNA5 mutations on atrial electrical activity. A multi-scale model of the human atria was updated to incorporate detailed experimental data on I_Kur from both wild-type and mutants. The effects of the mutations on human atrial action potential and rate dependence were investigated at the cellular level. In tissue, we assessed the effects of the mutations on the vulnerability to unidirectional conduction patterns and dynamics of re-entrant excitation waves. Gain-of-function mutations shortened the action potential duration in single cells, and stabilised and accelerated re-entrant excitation in tissue. Loss-of-function mutations had heterogeneous effects on action potential duration and promoted early-after-depolarisations following beta-adrenergic stimulation. In the tissue model, loss-of-function mutations facilitated breakdown of excitation waves at more physiological excitation rates than the wild-type, and the generation of early-after-depolarisations promoted unidirectional patterns of excitation. Gain- and loss-of-function I_Kur mutations produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups of mutations. This study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

In a recent study, six mutations resulting in either gain-of-function or loss-of-function in the ultra-rapid delayed rectifier potassium current I_Kur, were identified to be associated with atrial fibrillation (AF). However, the causative link between the mutant I_Kur (either gain- or loss-of-function) and AF genesis, especially the difference and similarity between the two mutant groups, has not been elucidated. In our study, we used multiscale computational models to investigate the mechanism of arrhythmogenesis mediated by the two groups of mutations. The results suggest that the gain-of-function mutations shortened atrial action potential duration, stabilised and accelerated re-entrant excitation waves in tissue; the loss-of-function mutation promoted early-after-depolarisations following beta-adrenergic stimulation and thus wave breaks in tissue. We show these two groups of mutants carrying I_Kur produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups. Our study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in I_K1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on I_K1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of I_K1 and I_Kr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas I_Kr block alone may be a safer alternative for the E299V mutation. Combined inhibition of I_Kr and I_Kur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, and is characterised by complex and irregular electrical activation of the upper chambers of the heart. One rare, genetic condition associated with increased risk of AF is the short QT syndrome (SQTS), which is caused by mutations in genes involved in normal electrical function of the heart. Underlying mechanisms by which SQTS-related gene mutations facilitate development of arrhythmias in the human atria are not well understood. In this study, sophisticated computer models representing ‘virtual’ human atria, incorporating detailed electrophysiological data at the ‘ion channel’ protein level into both idealised and realistic multi-scale tissue geometries, were used to dissect mechanisms by which two mutations in the KCNJ2 gene responsible for SQTS variant 3 (SQT3) promote initiation and sustenance of arrhythmias. It was found that the D172N and E299V mutations to KCNJ2 accelerated the repolarisation process at the cellular level through distinct mechanisms. This, along with the way the mutations affected heterogeneity in electrical behaviour at the organ level, mediated stability of arrhythmias and response to simulated ion channel block. This study improves understanding of mechanisms underlying increased AF risk associated with D172N and E299V KCNJ2 mutations, and outlines potential therapeutic strategies.

Novel non-invasive algorithm to identify the origins of re-entry and ectopic foci in the atria from 64-lead ECGs: A computational study

Atrial tachy-arrhytmias, such as atrial fibrillation (AF), are characterised by irregular electrical activity in the atria, generally associated with erratic excitation underlain by re-entrant scroll waves, fibrillatory conduction of multiple wavelets or rapid focal activity. Epidemiological studies have shown an increase in AF prevalence in the developed world associated with an ageing society, highlighting the need for effective treatment options. Catheter ablation therapy, commonly used in the treatment of AF, requires spatial information on atrial electrical excitation. The standard 12-lead electrocardiogram (ECG) provides a method for non-invasive identification of the presence of arrhythmia, due to irregularity in the ECG signal associated with atrial activation compared to sinus rhythm, but has limitations in providing specific spatial information. There is therefore a pressing need to develop novel methods to identify and locate the origin of arrhythmic excitation. Invasive methods provide direct information on atrial activity, but may induce clinical complications. Non-invasive methods avoid such complications, but their development presents a greater challenge due to the non-direct nature of monitoring. Algorithms based on the ECG signals in multiple leads (e.g. a 64-lead vest) may provide a viable approach. In this study, we used a biophysically detailed model of the human atria and torso to investigate the correlation between the morphology of the ECG signals from a 64-lead vest and the location of the origin of rapid atrial excitation arising from rapid focal activity and/or re-entrant scroll waves. A focus-location algorithm was then constructed from this correlation. The algorithm had success rates of 93% and 76% for correctly identifying the origin of focal and re-entrant excitation with a spatial resolution of 40 mm, respectively. The general approach allows its application to any multi-lead ECG system. This represents a significant extension to our previously developed algorithms to predict the AF origins in association with focal activities.

Atrial tachy-arrhythmias are associated with irregular excitation waves arising from re-entrant excitation, multiple wavelets or rapid focal activity. Identifying the origin of the irregular activity may be vital for diagnosis and treatment of the disorder. Where invasive and non-invasive methods provide approaches for such identification, both have associated disadvantages. In this study, we used a biophysically detailed model of the human atria and torso to develop an algorithm based on the correlation between the electrocardiogram (ECG) signal from a 64-lead vest and the location of rapid focal and re-entrant excitation. Using the properties of the atrial activation and the ECG signals, we developed a focus-location algorithm which is able to distinguish rapid focal activity from re-entrant scroll waves centred in the same location. Based on simulated data, the algorithm had success rates of 93% and 76% for correctly identifying the origin of focal and re-entrant excitation, respectively, and 88% for distinguish focal and re-entrant excitation, with no false positives. Inherited from our previous algorithm, it is also easily generalizable to any multi-lead ECG system.

A Computer Simulation Study of Anatomy Induced Drift of Spiral Waves in the Human Atrium

The interaction of spiral waves of excitation with atrial anatomy remains unclear. This simulation study isolates the role of atrial anatomical structures on spiral wave spontaneous drift in the human atrium. We implemented realistic and idealised 3D human atria models to investigate the functional impact of anatomical structures on the long-term (∼40 s) behaviour of spiral waves. The drift of a spiral wave was quantified by tracing its tip trajectory, which was correlated to atrial anatomical features. The interaction of spiral waves with the following idealised geometries was investigated: (a) a wedge-like structure with a continuously varying atrial wall thickness; (b) a ridge-like structure with a sudden change in atrial wall thickness; (c) multiple bridge-like structures consisting of a bridge connected to the atrial wall. Spiral waves drifted from thicker to thinner regions and along ridge-like structures. Breakthrough patterns caused by pectinate muscles (PM) bridges were also observed, albeit infrequently. Apparent anchoring close to PM-atrial wall junctions was observed. These observations were similar in both the realistic and the idealised models. We conclude that spatially altering atrial wall thickness is a significant cause of drift of spiral waves. PM bridges cause breakthrough patterns and induce transient anchoring of spiral waves.

Effective dynamics of twisted and curved scroll waves using virtual filaments

Scroll waves are three-dimensional excitation patterns that rotate around a central filament curve; they occur in many physical, biological, and chemical systems. We explicitly derive the equations of motion for scroll wave filaments in reaction-diffusion systems with isotropic diffusion up to third order in the filament's twist and curvature. The net drift components define at every instance of time a virtual filament which lies close to the instantaneous filament. Importantly, virtual filaments obey simpler, time-independent laws of motion which we analytically derive here and illustrate with numerical examples. Stability analysis of scroll waves is performed using virtual filaments, showing that filament curvature and twist add as quadratic terms to the nominal filament tension. Applications to oscillating chemical reactions and cardiac tissue are discussed.

Simulating the role of anisotropy in human atrial cardioversion

This computational study quantifies the effectiveness of feedback controlled low energy cardioversion in the anisotropic human atria. An established biophysical human cell model was adopted to reproduce Control and chronic atrial fibrillation (CAF) action potentials. The cell model was combined with a detailed human atrial geometry to construct a 3D realistic human atrial model. Scroll waves were simulated under Control and CAF conditions and the cardioversion parameters of stimulation strength and pacing duration were evaluated for scroll wave termination. Scroll waves were initiated at two locations in the atria to elicit the effects of scroll wave location. The role of anisotropy was highlighted by comparison to results from the isotropic case. Under Control conditions, scroll wave self-termination was rapid in the anisotropic case. Under CAF conditions, anisotropy caused the initiated scroll wave to degenerate into multiple scrolls with each evolving erratically or pinning to anatomical defects. The cardioversion successfully terminated scroll waves within 10 s, but the stimulus strength had a strong correlation to the location of the scroll wave. The low energy stimulation strength was always lower than the threshold stimulus. Anisotropy plays an important role in atrial electrical properties. Anisotropy aggravates CAF and leads to high frequency atrial pacing. The efficacy of cardioversion is significantly affected by anisotropy.

Resonant drift of spiral waves in the complex ginzburg-landau equation

Weak periodic external perturbations of an autowave medium can cause large-distance directed motion of the spiral wave. This happens when the period of the perturbation coincides with, or is close to the rotation period of a spiral wave, or its multiple. Such motion is called resonant or parametric drift. It may be used for low-voltage defibrillation of heart tissue. Theory of the resonant drift exists, but so far was used only qualitatively. In this paper, we show good quantitative agreement of the theory with direct numerical simulations. This is done for Complex Ginzburg-Landau Equation, one of the simplest autowave models.

Nonequilibrium arrhythmic states and transitions in a mathematical model for diffuse fibrosis in human cardiac tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to fibroblasts via a gap-junctional conductance , reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.

Experiment-model interaction for analysis of epicardial activation during human ventricular fibrillation with global myocardial ischaemia

We describe a combined experiment-modelling framework to investigate the effects of ischaemia on the organisation of ventricular fibrillation in the human heart. In a series of experimental studies epicardial activity was recorded from 10 patients undergoing routine cardiac surgery. Ventricular fibrillation was induced by burst pacing, and recording continued during 2.5 min of global cardiac ischaemia followed by 30 s of coronary reflow. Modelling used a 2D description of human ventricular tissue. Global cardiac ischaemia was simulated by (i) decreased intracellular ATP concentration and subsequent activation of an ATP sensitive K⁺ current, (ii) elevated extracellular K⁺ concentration, and (iii) acidosis resulting in reduced magnitude of the L-type Ca²⁺ current I(Ca,L). Simulated ischaemia acted to shorten action potential duration, reduce conduction velocity, increase effective refractory period, and flatten restitution. In the model, these effects resulted in slower re-entrant activity that was qualitatively consistent with our observations in the human heart. However, the flattening of restitution also resulted in the collapse of many re-entrant waves to several stable re-entrant waves, which was different to the overall trend we observed in the experimental data. These findings highlight a potential role for other factors, such as structural or functional heterogeneity in sustaining wavebreak during human ventricular fibrillation with global myocardial ischaemia.

Hysteresis and bistability in periodically paced cardiac tissue

Hysteresis in periodically paced cardiac tissue is an important issue due to its relevance to cardiac arrhythmias. In the present paper, the mechanism of hysteresis formation and the related properties are interpreted by numerically investigating the phase I Luo-Rudy model. A formula calculating the width of hysteresis is proposed and well confirmed by numerical simulations. We also find that hysteresis in cardiac tissue shows several characteristics due to couplings among cardiac cells which are absent in a single cell. The influences of the physiological parameters are studied in detail. The model dependence of hysteresis is elucidated by considering a number of well-known models of excitable media. Moreover, the influence of bistability on controlling arrhythmias is revealed.

Spiral-wave turbulence and its control in the presence of inhomogeneities in four mathematical models of cardiac tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control schemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

Vortex filament dynamics in computational models of ventricular fibrillation in the heart

In three-dimensional cardiac tissue, the re-entrant waves that sustain ventricular fibrillation rotate around a line of phase singularity or vortex filament. The aim of this study was to investigate how the behavior of these vortex filaments is influenced by membrane kinetics, initial conditions, and tissue geometry in computational models of excitable tissue. A monodomain model of cardiac tissue was used, with kinetics described by a three-variable simplified ionic model (3V-SIM). Two versions of 3V-SIM were used, one with steep action potential duration restitution, and one with reduced excitability. Re-entrant fibrillation was then simulated in three tissue geometries: a cube, a slab, and an anatomically detailed model of rabbit ventricles. Filaments were identified using a phase-based method, and the number, size, origin, and orientation of filaments was tracked throughout each simulation. The main finding of this study is that kinetics, initial conditions, geometry, and anisotropy all affected the number, proliferation, and orientation of vortex filaments in re-entrant fibrillation. An important finding of this study was that the behavior of vortex filaments in simplified slab geometry representing part of the ventricular wall did not necessarily predict behavior in an anatomically detailed model of the rabbit ventricles.

Spatiotemporal frequency analysis of ventricular fibrillation in explanted human hearts

Ventricular fibrillation (VF) is a medical condition that occurs due to rapid and irregular electrical activity of heart. If undiagnosed or untreated, VF leads to sudden cardiac death. VF has been studied by researchers for over 100 years to elucidate the mechanism that maintains VF, and thus to arrive at therapeutic options. VF is a nonstationary process, and it manifests into variations in the waveform morphology, phase, and frequency dynamics of the surface electrograms. Dominant frequency analysis (DF maps) and phase maps are two widely used complementary approaches in assessing the evolution of VF process. These techniques are applied to electrograms or fluorescence signals obtained with voltage-sensitive dyes. In spite of VF being a nonstationary process, most of the existing literature limits frequency analysis to a segmented, time-averaged spectral analysis, where valuable information on the instantaneous temporal evolution of the spectral characteristics is lost. In order to resolve this issue, in this paper, we present a joint time-frequency approach that is suited for VF analysis and demonstrate the application of instantaneous mean frequency (IMF) in interpreting VF episodes. Human VF sources are rarely anatomically stable and are migratory. Traditional DF techniques fail in tracking this migratory behavior. IMF, on the other hand, can deal with these migratory sources and conduction blocks better than DF approaches. Results of the analysis using the electrograms of 204 VF segments obtained from 13 isolated human hearts (explanted during cardiac transplantation) indicate that in 81% of the VF segments, there were significant changes in the spatiotemporal evolution of the frequency, suggesting that IMF provides better mechanistic insight of these signals. The IMF tool presented in this paper demonstrates potential for applications in tracking frequency patterns, conduction blocks, and arriving at newer therapies to modulate VF.

Active and passive control of spiral turbulence in excitable media

The influence of a spatially localized heterogeneity defect, defined by failure of the diffusion effect, on spiral turbulence suppression in two-dimensional excitable media is studied numerically, based on the Bär model. It is shown that in certain parameter regions spiral turbulence without the defect can be suppressed by a boundary periodic forcing (called active control) if the forcing frequency is properly chosen. However, with a sufficiently large defect this active control method no longer works due to the wake turbulence following the defect. We suggest an auxiliary method of enclosing the defect with a thin layer of material of high excitability (called passive control) to screen the interaction between the defect and the turbulence and to restore the global control effect of the periodic forcing. The possible application of the method in cardiac defibrillation is discussed.

The canine virtual ventricular wall: A platform for dissecting pharmacological effects on propagation and arrhythmogenesis

We have constructed computational models of canine ventricular cells and tissues, ultimately combining detailed tissue architecture and heterogeneous transmural electrophysiology. The heterogeneity is introduced by modifying the Hund-Rudy canine cell model in order to reproduce experimentally reported electrophysiological properties of endocardial, midmyocardial (M) and epicardial cells. These models are validated against experimental data for individual ionic current and action potential characteristics, and their rate dependencies. 1D and 3D heterogeneous virtual tissues are constructed, with detailed tissue architecture (anisotropy and orthotropy, due to fibre orientation and sheet structure) of the left ventricular wall wedge extracted from a diffusion tensor imaging data set. The models are used to study the effects of tissue heterogeneity and class III drugs on transmural propagation and tissue vulnerability to re-entry. We have determined relationships between the transmural dispersion of action potential duration (APD) and the vulnerable window in the 1D virtual ventricular wall, and demonstrated how changes in the transmural heterogeneity, and hence tissue vulnerability, can lead to generation of re-entry in the 3D ventricular wedge. Two class III drugs with opposite qualitative effects on transmural APD heterogeneity are considered: d-sotalol that increases transmural APD dispersion, and amiodarone that decreases it. Simulations with the 1D virtual ventricular wall show that under d-sotalol conditions the vulnerable window is substantially wider compared to amiodarone conditions, primarily in the epicardial region where unidirectional conduction block persists until the adjacent M cells are fully repolarised. Further simulations with the 3D ventricular wedge have shown that ectopic stimulation of the epicardial region results in generation of sustained re-entry under d-sotalol conditions, but not under amiodarone conditions or in control. Again, APD increase in M cells was identified as the major contributor to tissue vulnerability--re-entry was initiated primarily due to ectopic excitation propagating around the unidirectional conduction block in the M cell region. This suggests an electrophysiological mechanism for the anti- and proarrhythmic effects of the class III drugs: the relative safety of amiodarone in comparison to d-sotalol can be explained by relatively low transmural APD dispersion, and hence, a narrow vulnerable window and low probability of re-entry in the tissue.

Patterns of phase-dependent spiral wave attenuation in excitable media

We show that introducing periodic planar fronts with long excitation duration can lead to spiral attenuation. The attenuation occurs periodically over cycles of several planar fronts, forming a variety of complex spatiotemporal patterns. We find that these attenuation patterns occur only at specific phases of the descending fronts relative to the rotational phase of the spiral. These patterns fall into two general classes, each defined by a specific expression for the number of attenuated spirals per cycle of planar fronts, and represented by a structured diagram in parameter space. The spiral attenuation patterns we observe remain stable in time and do not change during the evolution of the system.

Virtual cell and tissue dynamics of ectopic activation of the ventricles

Cardiac ventricular cells and tissues are normally excitable, and are activated by propagating waves of excitation that are initiated in the specialized pacemaking region of the heart. However, isolated or repetitive activity can be initiated at abnormal (ectopic) sites in the ventricles. To trigger an endogenous ectopic beat, there must be a compact focus of cells with changed membrane excitation parameters and kinetics, which initiate activity by after-depolarizations triggered by propagating activity, or that have bifurcated into autorhythmicity. This ectopic focus needs to be large enough, and adequately coupled, to drive the surrounding tissue. We investigate the initiation of ectopic excitation in computational models of human ventricular cells triggered by after-depolarizations and by up/down regulation of specific membrane conductance systems, and the propagation and evolution of ectopic activity in homogeneous or heterogeneous and isotropic, anisotropic, or orthotropic tissues.

Patterns of spiral wave attenuation by low-frequency periodic planar fronts

There is evidence that spiral waves and their breakup underlie mechanisms related to a wide spectrum of phenomena ranging from spatially extended chemical reactions to fatal cardiac arrhythmias [A. T. Winfree, The Geometry of Biological Time (Springer-Verlag, New York, 2001); J. Schutze, O. Steinbock, and S. C. Muller, Nature 356, 45 (1992); S. Sawai, P. A. Thomason, and E. C. Cox, Nature 433, 323 (2005); L. Glass and M. C. Mackey, From Clocks to Chaos: The Rhythms of Life (Princeton University Press, Princeton, 1988); R. A. Gray et al., Science 270, 1222 (1995); F. X. Witkowski et al., Nature 392, 78 (1998)]. Once initiated, spiral waves cannot be suppressed by periodic planar fronts, since the domains of the spiral waves grow at the expense of the fronts [A. N. Zaikin and A. M. Zhabotinsky, Nature 225, 535 (1970); A. T. Stamp, G. V. Osipov, and J. J. Collins, Chaos 12, 931 (2002); I. Aranson, H. Levine, and L. Tsimring, Phys. Rev. Lett. 76, 1170 (1996); K. J. Lee, Phys. Rev. Lett. 79, 2907 (1997); F. Xie, Z. Qu, J. N. Weiss, and A. Garfinkel, Phys. Rev. E 59, 2203 (1999)]. Here, we show that introducing periodic planar waves with long excitation duration and a period longer than the rotational period of the spiral can lead to spiral attenuation. The attenuation is not due to spiral drift and occurs periodically over cycles of several fronts, forming a variety of complex spatiotemporal patterns, which fall into two distinct general classes. Further, we find that these attenuation patterns only occur at specific phases of the descending fronts relative to the rotational phase of the spiral. We demonstrate these dynamics of phase-dependent spiral attenuation by performing numerical simulations of wave propagation in the excitable medium of myocardial cells. The effect of phase-dependent spiral attenuation we observe can lead to a general approach to spiral control in physical and biological systems with relevance for medical applications.

Ventricular fibrillation in myopathic human hearts: mechanistic insights from in vivo global endocardial and epicardial mapping

Ventricular fibrillation (VF) is an important cause of sudden cardiac death and cardiovascular mortality in patients with cardiomyopathy. Although it was generally believed that chaotic reentrant wavefronts underlie VF in humans, there is emerging evidence of spatiotemporal organization during early VF. The mechanism of this organization of electrical activity in early VF is unknown in myopathic hearts. We studied early VF in vivo, intraoperatively in five cardiomyopathic patients. Simultaneous electrograms were obtained from the epicardium and endocardium in left ventricular cardiomyopathy and from the endocardium in right ventricular myopathy. The Hilbert transform was used to derive the phase of the electrograms. Rotors were identified by isolating phase singularity points. Rotors were present in all of the myopathic hearts studied during VF and cumulatively lasted a mean of 3.2 +/- 2.0 s of the 7.0 +/- 4.0 s of the VF segments analyzed. For each surface mapped, 3.6 +/- 2.9 rotors were identified for the duration mapped. The average number of cycles completed by these rotors was 4.9 +/- 4.9. The longest rotor lasted 10.2 +/- 6.2 rotations and lasted 2.0 +/- 1.2 s. The rotors on the endocardium had a cycle length of 192 +/- 33 ms compared with 220 +/- 15 ms on the epicardium (P=0.08). There is centrifugal activation of electrical activity from these rotors, and they give rise to domains that activate at faster rates with evidence of conduction block at the border with slower domains. These rotors frequently localized to border regions of myocardium with bipolar electrogram amplitude of <0.5 mV. The organization of electrical activity during early VF in myopathic human hearts is characterized by wavefronts emanating from a few rotors.

Spiral-wave dynamics depend sensitively on inhomogeneities in mathematical models of ventricular tissue

Every sixth death in industrialized countries occurs because of cardiac arrhythmias such as ventricular tachycardia (VT) and ventricular fibrillation (VF). There is growing consensus that VT is associated with an unbroken spiral wave of electrical activation on cardiac tissue but VF with broken waves, spiral turbulence, spatiotemporal chaos and rapid, irregular activation. Thus spiral-wave activity in cardiac tissue has been studied extensively. Nevertheless, many aspects of such spiral dynamics remain elusive because of the intrinsically high-dimensional nature of the cardiac-dynamical system. In particular, the role of tissue heterogeneities in the stability of cardiac spiral waves is still being investigated. Experiments with conduction inhomogeneities in cardiac tissue yield a variety of results: some suggest that conduction inhomogeneities can eliminate VF partially or completely, leading to VT or quiescence, but others show that VF is unaffected by obstacles. We propose theoretically that this variety of results is a natural manifestation of a complex, fractal-like boundary that must separate the basins of the attractors associated, respectively, with spiral breakup and single spiral wave. We substantiate this with extensive numerical studies of Panfilov and Luo-Rudy I models, where we show that the suppression of spiral breakup depends sensitively on the position, size, and nature of the inhomogeneity.

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.

Preventing alternans-induced spiral wave breakup in cardiac tissue: an ion-channel-based approach

The detailed processes involved in spiral wave breakup, believed to be one major mechanism by which tachycardia evolves into fibrillation, are still poorly understood. This has rendered difficult the proper design of an efficient and practical control stimulus protocol to eliminate such events. In order to gain new insights into the underlying electrophysiological and dynamical mechanisms of breakup, we applied linear perturbation theory to a steadily rotating spiral wave in two spatial dimensions. The tissue was composed of cells modeled using the Fenton-Karma equations whose parameters were chosen to emphasize alternans as a primary mechanism for breakup. Along with one meandering mode, not just one but several unstable alternans modes were found with differing growth rates, frequencies, and spatial structures. As the conductance of the fast inward current was increased, the instability of the modes increased, consistent with increased meandering and propensity for spiral breakup in simulations. We also explored a promising new approach, based on the theory, for the design of an energy efficient electrical stimulus protocol to control spiral wave breakup. The novelty lies in addressing the problem directly at the ion channel level and taking advantage of the inherent two dimensional nature of the rotating wave. With the help of the eigenmode method, we were able to calculate the exact timing and amplitude of the stimulus, and locate it optimally to maximize efficiency. The analysis led to a special-case example that demonstrated that a single, properly timed stimulus can have a global effect, suppressing all growing alternans modes over the entire tissue, thus inhibiting spiral wave breakup.

Pinning force in active media

Pinning of vortices by defects plays an important role in various physical (superconductivity, superfluidity, etc.) or biological (propagation in cardiac muscle) situations. Which defects act as pinning centers? We propose a way to study this general problem by using an advection field to quantify the attraction between an obstacle and a vortex. A full solution is obtained for the real Ginzburg-Landau equation (RGLE). Two pinning mechanisms are found in excitable media. Our results suggest strong analogies with the RGLE when the heterogeneity is excitable. Unpinning from an unexcitable obstacle is qualitatively harder, resulting in a stronger pinning force. We discuss the implications of our results to control vortices and propose experiments in a chemical active medium and in cardiac tissue.

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.

Filament behavior in a computational model of ventricular fibrillation in the canine heart

The aim of this paper was to quantify the behavior of filaments in a computational model of re-entrant ventricular fibrillation. We simulated cardiac activation in an anisotropic monodomain with excitation described by the Fenton-Karma model with Beeler-Reuter restitution, and geometry by the Auckland canine ventricle. We initiated re-entry in the left and right ventricular free walls, as well as the septum. The number of filaments increased during the first 1.5 s before reaching a plateau with a mean value of about 36 in each simulation. Most re-entrant filaments were between 10 and 20 mm long. The proportion of filaments touching the epicardial surface was 65%, but most of these were visible for much less than one period of re-entry. This paper shows that useful information about filament dynamics can be gleaned from models of fibrillation in complex geometries, and suggests that the interplay of filament creation and destruction may offer a target for antifibrillatory therapy.

Endocardial detection of repolarization alternans

Repolarization alternans (RPA) is prognostic of sudden cardiac death and is thought to be mechanistically linked to the initiation of ventricular tachyarrhythmias. Thus, implantable cardiac device detection of RPA may be therapeutically valuable. Because alternans detection is currently limited to surface electrocardiograms, we investigated whether RPA could be measured using a single right-ventricular endocardial lead in humans. Such a location was chosen because it is consistent with the requirements for long-term implantable-device implementation. During diagnostic electrophysiological testing, 28 patients (23 male, 5 female; 61 +/- 15 years) were evaluated for surface T-wave alternans (TWA; the current "gold standard" for RPA detection) and endocardial RPA during 5 min of 550-ms right-atrial pacing. Power spectral analysis indicated that 11/28 patients had both surface TWA and endocardial RPA, 9/28 patients had neither, and 8/28 patients had discordant results (71% concordance; p = 0.02). Importantly, unlike surface TWA, endocardial RPA was detectable on a beat-to-beat basis. Given the putative mechanistic link between RPA and ventricular arrhythmias, beat-to-beat endocardial RPA detection might be of diagnostic or therapeutic utility.

Analysing the ventricular fibrillation waveform

The surface electrocardiogram associated with ventricular fibrillation has been of interest to researchers for some time. Over the last few decades, techniques have been developed to analyse this signal in an attempt to obtain more information about the state of the myocardium and the chances of successful defibrillation. This review looks at the implications of analysing the VF waveform and discusses the various techniques that have been used, including fast Fourier transform analysis, wavelet transform analysis and mathematical techniques such as chaos theory.

Mathematical modelling for the new millenium: medicine by numbers

Physicists, engineers and mathematicians are accustomed to the combination of elegance, rigour and utility that characterise mathematical models. They are familiar with the need to dip into their mathematical toolbox to select the technique of choice. However, medicine and biology have not been characterised, in general, by a mathematical formalism. The relative paucity of mathematical models in biology and medicine reflects in part the difficulty in making accurate and appropriate experimental measurements in the field. Signal noise, the lack of appropriate sensors, and uncertainty as to what constitutes the significant measurements are largely to blame for this. The purpose of this paper is to characterise a 'good' model, encourage the development and application of such models to new areas, and outline future developments in the field. It is proposed that a good model will be accurate, predictive, economical, unique and elegant. These principles will be illustrated with reference to four models: radiosensitisation of tumours, modelling solute clearance in haemodialysis, the myogenic response in reactive hyperaemia and cardiac electrical activity. It is suggested that, in the immediate future, the mathematical model will become a useful adjunct to laboratory experiment (and possibly clinical trial), and the provision of 'in silico' models will become routine.

Introduction: Mapping and control of complex cardiac arrhythmias

This paper serves as an introduction to the Focus Issue on mapping and control of complex cardiac arrhythmias. We first introduce basic concepts of cardiac electrophysiology and describe the main clinical methods being used to treat arrhythmia. We then provide a brief summary of the main themes contained in the articles in this Focus Issue. In recent years there have been important advances in the ability to map the spread of excitation in intact hearts and in laboratory settings. This work has been combined with simulations that use increasingly realistic geometry and physiology. Waves of excitation and contraction in the heart do not always propagate with constant velocity but are often subject to instabilities that may lead to fluctuations in velocity and cycle time. Such instabilities are often treated best in the context of simple one- or two-dimensional geometries. An understanding of the mechanisms of propagation and wave stability is leading to the implementation of different stimulation protocols in an effort to modify or eliminate abnormal rhythms. (c) 2002 American Institute of Physics.

Enhanced self-termination of re-entrant arrhythmias as a pharmacological strategy for antiarrhythmic action

Ventricular tachycardia and fibrillation are potentially lethal cardiac arrhythmias generated by high frequency, irregular spatio-temporal electrical activity. Re-entrant propagation has been demonstrated as a mechanism generating these arrhythmias in computational and in vitro animal models of these arrhythmias. Re-entry can be idealised in homogenous isotropic virtual cardiac tissues as spiral and scroll wave solutions of reaction-diffusion equations. A spiral wave in a bounded medium can be terminated if its core reaches a boundary. Ventricular tachyarrhythmias in patients are sometimes observed to spontaneously self-terminate. One possible mechanism for self-termination of a spiral wave is meander of its core to an inexcitable boundary. We have previously proposed the hypothesis that the spatial extent of meander of a re-entrant wave in the heart can be directly related to its probability of self-termination, and so inversely related to its lethality. Meander in two-dimensional virtual ventricular tissues based on the Oxsoft family of cell models, with membrane excitation parameters simulating the inherited long Q-T syndromes has been shown to be consistent with this hypothesis: the largest meander is seen in the syndrome with the lowest probability of death per arrhythmic episode. Here we extend our previous results to virtual tissues based on the Luo-Rudy family of models. Consistent with our hypothesis, for both families of models, whose different ionic mechanisms produce different patterns of meander, the LQT virtual tissue with the larger meander simulates the syndrome with the lower probability of death per episode. Further, we search the parameter space of the repolarizing currents to find their conductance parameter values that give increased meander of spiral waves. These parameters may provide targets for antiarrhythmic drugs designed to act by increasing the likelihood of self-termination of re-entrant arrhythmias. (c) 2002 American Institute of Physics.

A method to quantify the dynamics and complexity of re-entry in computational models of ventricular fibrillation

Ventricular fibrillation is a deadly cardiac arrhythmia. There is evidence that electrical activity in cardiac tissue is sustained during fibrillation by re-entrant waves that rotate around filaments. In this paper we develop a method for identifying and tracking filaments in a computational model of ventricular fibrillation. This method identifies the birth, death, bifurcation and amalgamation of filaments and these events are summarized on a directed graph. The approach described in this study provides ways to quantify the complex patterns of electrical activity seen in computational models of fibrillation, to relate the behaviour of computational models to experimental data and thus to gain insights into the underlying mechanisms of this dangerous arrhythmia.

Characterization of patterned irregularity in locally interacting, spatially extended systems: Ventricular fibrillation

The re-entrant ventricular arrhythmias of monomorphic ventricular tachycardia and fibrillation are produced by abnormal spatio-temporal patterns of propagation in the ventricular myocardium. These behaviors can be described by solutions of reaction-diffusion equation excitable medium models. The direct comparison of such solutions with existing experimental observations is virtually impossible as there are too many factors to be taken into account, including not only the complicated dynamics of the re-entrant waves of excitation in the tissue, but also the way the appearance of these waves on the surface is modified by the inhomogeneity, anisotropy and three-dimensional nature of heart tissue. One way of indirect comparison is to compare characteristics of the complexity of the model and the real data, that are invariant under these modifications of the signal. Karhunen-Loeve decomposition is a standard tool for evaluating the complexity of multidimensional signals. A comparison of the separate and conjoint complexities of the signals on the opposite sides of the preparation can be considered as an indicator how much three-dimensional effects are essential in the preparation behavior. (c) 2001 American Institute of Physics.

Computational models of normal and abnormal action potential propagation in cardiac tissue: linking experimental and clinical cardiology

Computational models have the potential to make a huge impact on our understanding of normal and abnormal cardiac function. The aim of this article is to review tools that have been developed to simulate the electrophysiology of cardiac cells and tissue, and to show how computational models have been used to gain insight into normal and abnormal action potential propagation. Some of the practical problems experienced in the development and application of these models are described, and examples are given.

Nonlinear-dynamical arrhythmia control in humans

Nonlinear-dynamical control techniques, also known as chaos control, have been used with great success to control a wide range of physical systems. Such techniques have been used to control the behavior of in vitro excitable biological tissue, suggesting their potential for clinical utility. However, the feasibility of using such techniques to control physiological processes has not been demonstrated in humans. Here we show that nonlinear-dynamical control can modulate human cardiac electrophysiological dynamics by rapidly stabilizing an unstable target rhythm. Specifically, in 52/54 control attempts in five patients, we successfully terminated pacing-induced period-2 atrioventricular-nodal conduction alternans by stabilizing the underlying unstable steady-state conduction. This proof-of-concept demonstration shows that nonlinear-dynamical control techniques are clinically feasible and provides a foundation for developing such techniques for more complex forms of clinical arrhythmia.

Re-entrant cardiac arrhythmias in computational models of long QT myocardium

The long QT syndrome (LQTS) is an inherited disorder in which repolarization of cardiac ventricular cells is prolonged. Patients with the LQTS are at an increased risk of ventricular cardiac arrhythmias. Two phenotypes of the inherited LQTS are caused by defects in K(+)channels (LQT1 and LQT2) and one by defects in Na(+)channels (LQT3). Patients with LQT1 are more likely to have self-terminating arrhythmias than those with LQT3. The aim of this computational study was to propose an explanation for this finding by comparing the vulnerability of normal and LQT tissue to re-entry, and estimating the likelihood of self-termination by motion of re-entrant waves to an inexcitable boundary in simulated LQT1, LQT2 and LQT3 tissue. We modified a model of mammalian cardiac cells to simulate LQT1 by reducing maximal I(K(s))conductance, LQT2 by reducing maximal I(K(r))conductance, and LQT3 by preventing complete inactivation of I(Na)channels. Each simulated phenotype was incorporated into a computational model of action potential propagation in one- and two-dimensional homogeneous tissue. Simulated LQT tissue was no more vulnerable to re-entry than simulated normal tissue, but the motion of re-entrant waves in simulated LQT1 tissue was between 2 and 5 times greater than the motion of re-entrant waves in simulated LQT2 and LQT3 tissue. These findings suggest that LQT arrhythmias do not result from increased vulnerability to re-entry, and that re-entry once initiated is more likely to self-terminate by moving to an inexcitable tissue boundary in LQT1 than in LQT2 and LQT3. This finding is consistent with clinical observations.

Deterministic nonlinearity in ventricular fibrillation

We provide numerical evidence that the electrocardiogram data collected from pigs during induced ventricular fibrillation cannot be described by a monotonic nonlinear transformation of linearly filtered noise. To establish this we use surrogate techniques and apply two test statistics: (1) the Takens' maximum likelihood estimator of the Grassberger-Procaccia correlation dimension and (2) an improved correlation dimension estimation routine. The improved dimension estimates provide evidence that the correlation dimension of the underlying dynamics during the episode of VF in the first 30 s is slightly less than 6. This result is consistent and reproducible among subjects. (c) 2000 American Institute of Physics.

The biphasic mystery: why a biphasic shock is more effective than a monophasic shock for defibrillation

We demonstrate that a biphasic shock is more effective than a monophasic shock at eliminating reentrant electrical activity in an ionic model of cardiac ventricular electrical activity. This effectiveness results from early hyperpolarization that enhances the recovery of sodium inactivation, thereby enabling earlier activation of recovering cells. The effect can be seen easily in a model of a single cell and also in a cable model with a ring of excitable cells. Finally, we demonstrate the phenomenon in a two-dimensional model of cardiac tissue.

Evolving perspectives during 12 years of electrical turbulence

This Focus issue describes a problem in electrical dynamics which has fascinated generations of physiologists. There are today so many views of fibrillation that only the rarest generalization can embrace all of them. Fifty-two prominent investigators collaborate here to present aspects of the problem in these eighteen articles (including this introduction) tailored for readers whose principal expertise lies elsewhere. In "The High One's Lay" (Norse Runes, ca. 800) Odin remarks, "Much too early I came to many places: the beer was not yet ready, or was already drunk em leader " but to this one we come at very nearly the right time in 1998. This introduction attempts to guide newcomers by noting the changed or multiple meanings of novel technical terms while sorting the key facts and ideas into an order that facilitates comparison and contrast with those of a dozen years ago. This Focus issue is authored by some of the foremost innovators of both theory and experiment in this area. By assimilating their presentations the readers of Chaos can become well poised to appreciate and evaluate the definitive evidence expected in the next few years. (c) 1998 American Institute of Physics.

Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation

Wave propagation in ventricular muscle is rendered highly anisotropic by the intramural rotation of the fiber. This rotational anisotropy is especially important because it can produce a twist of electrical vortices, which measures the rate of rotation (in degree/mm) of activation wavefronts in successive planes perpendicular to a line of phase singularity, or filament. This twist can then significantly alter the dynamics of the filament. This paper explores this dynamics via numerical simulation. After a review of the literature, we present modeling tools that include: (i) a simplified ionic model with three membrane currents that approximates well the restitution properties and spiral wave behavior of more complex ionic models of cardiac action potential (Beeler-Reuter and others), and (ii) a semi-implicit algorithm for the fast solution of monodomain cable equations with rotational anisotropy. We then discuss selected results of a simulation study of vortex dynamics in a parallelepipedal slab of ventricular muscle of varying wall thickness (S) and fiber rotation rate (theta(z)). The main finding is that rotational anisotropy generates a sufficiently large twist to destabilize a single transmural filament and cause a transition to a wave turbulent state characterized by a high density of chaotically moving filaments. This instability is manifested by the propagation of localized disturbances along the filament and has no previously known analog in isotropic excitable media. These disturbances correspond to highly twisted and distorted regions of filament, or "twistons," that create vortex rings when colliding with the natural boundaries of the ventricle. Moreover, when sufficiently twisted, these rings expand and create additional filaments by further colliding with boundaries. This instability mechanism is distinct from the commonly invoked patchy failure or wave breakup that is not observed here during the initial instability. For modified Beeler-Reuter-like kinetics with stable reentry in two dimensions, decay into turbulence occurs in the left ventricle in about one second above a critical wall thickness in the range of 4-6 mm that matches experiment. However this decay is suppressed by uniformly decreasing excitability. Specific experiments to test these results, and a method to characterize the filament density during fibrillation are discussed. Results are contrasted with other mechanisms of fibrillation and future prospects are summarized. (c)1998 American Institute of Physics.

The effect of gap junctional distribution on defibrillation

We summarize a mathematical theory for direct activation and defibrillation of cardiac tissue. We show that the direct stimulus and defibrillation thresholds are likely to be strongly affected by the gap junctional distribution and density, suggesting an indirect experimental test of the theory. (c) 1998 American Institute of Physics.